Effect of linkage disequilibrium between markers in linkage and association analyses

Contributions to Group 17 of the Genetic Analysis Workshop 15 considered dense markers in linkage disequilibrium (LD) in the context of either linkage or association analysis. Three contributions reported on methods for modeling LD or selecting a subset of markers in linkage equilibrium to perform linkage analysis. When all markers were used without modeling LD, inflated evidence for linkage was observed when parental genotypes were missing. All methods for handling LD led to some decreased linkage evidence. Two groups performed a genome-wide association scan using either mixed models to account for known or unknown relatedness between individuals, trend tests or combination statistics. All methods failed to detect four of the eight simulated loci because of low LD in some regions. Three groups performed association analysis using simulated dense markers on chromosome 6, where a simulated HLA-DRB1 locus played a major role in disease susceptibility along with two additional loci of smaller effect. The overall conditional genotype method correctly identified both additional loci while a novel transmission disequilibrium test-statistic to combine studies with non-overlapping markers identified one HLA locus after stratifying on the parental HLA-DRB1 genotypes; LD mapping using the Malécot model mapped two loci in this region, even when using greatly reduced marker density. While LD between markers appears to be a nuisance that may cause spurious linkage results with missing parental genotypes in linkage analysis, association analysis thrives on LD, and disease genes fail to be detected in regions of low LD.     

Genetic analysis of natural recombinant Brazilian Toxoplasma gondii strains by multilocus PCR-RFLP.

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) at eight independent loci was used to determine the types (I, II or III) lineage of 20 Toxoplasma gondii strains isolated from humans and animals in Brazil. RH (type I, highly virulent), ME49 (type II, avirulent) and VEG (type III, avirulent) were used as reference strains. Differently from expected frequencies, all Brazilian strains showed to have recombinant genotypes, with typical alleles of types I, II or III at almost all loci assessed. The cB21-4 locus, a microsatellite marker, showed a higher allelic polymorphism with seven alleles among strains under analysis. Data have also shown that many Brazilian T. gondii strains presented a new haplotype at the L363 locus. When results of the eight loci were combined, 14 schizodemes were characterized out of the 20 T. gondii strains isolated in Brazil. The phenogram representing PCR-RFLP data separated Brazilian strains into two distinct genetic groups associated with murine virulence phenotype, termed groups I-A and I-B. Strains from group I-A (AS28, BV and N) that were highly virulent in BALB/c mice, were clustered with RH reference strain. Only those strains presented the haplotype I at the L363 locus, suggesting that this could be a possible marker of highly virulent strains. Strains from group I-B (cystogenic strains) showed a more heterogeneous behavior regarding virulence: a few of them (EGS, RAR, SAF, D5 and D6) were virulent, others (C4, P and D8) avirulent and most of them (D1, D2, D3, D4, D7, EFP, CH1, CH2 and CH3) intermediate virulent in mice. A significant linkage disequilibrium was observed in the population surveyed. However, the role of sexual recombination in the population structure of T. gondii in Brazil seems to be more central than in Europe and North America, where most studies have been performed.     

Generalization of the extended transmission disequilibrium test to two unlinked disease loci

The extended transmission disequilibrium test (ETDT) of Sham and Curtis [1995] is a powerful test of the null hypothesis of no linkage between a multi-allelic marker locus and a disease susceptibility locus of unknown location in the presence of association between alleles at the two loci. We propose a generalization of the ETDT to test simultaneously for linkage of markers to two unlinked disease loci. Analysis of the simulated data for Problem 2 would suggest improvements in power to detect linkage to pairs of disease loci over single locus tests. This is as a result of the two locus model taking account of epistasis between the disease loci.     

Evidence of substantial recombination among Trypanosoma cruzi II strains from Minas Gerais

Due to the scarcity of evidence of sexuality in Trypanosoma cruzi, the causative agent of Chagas disease, it has been general accepted that the parasite reproduction is essentially clonal with infrequent genetic recombination. This assumption is mainly supported by indirect evidence, such as Hardy–Weinberg imbalances, linkage disequilibrium and a strong correlation between independent sets of genetic markers of T. cruzi populations. However, because the analyzed populations are usually isolated from different geographic regions, the possibility of population substructuring as generating these genetic marker imbalances cannot be eliminated. To investigate this possibility, we firstly compared the allele frequencies and haplotype networks using seven different polymorphic loci (two from mitochondrial and five from different nuclear chromosomes) in two groups of TcII strains: one including isolates obtained from different regions in Latin America and the other including isolates obtained only from patients of the Minas Gerais State in Brazil. Our hypothesis was that if the population structure is essentially clonal, Hardy–Weinberg disequilibrium and a sharp association between the clusters generated by analyzing independent markers should be observed in both strain groups, independent of the geographic origin of the samples. The results demonstrated that the number of microsatellite loci in linkage disequilibrium decreased from 4 to 1 when only strains from Minas Gerais were analyzed. Moreover, we did not observed any correlation between the clusters when analyzing the nuclear and mitochondrial loci, suggesting independent inheritance of these markers among the Minas Gerais strains. Besides, using a second subset of five physically linked microsatellite loci and the Minas Gerais strains, we could also demonstrate evidence of homologous recombination roughly proportional to the relative distance among them. Taken together, our results do not support a clonal population structure for T. cruzi, particularly in TcII, which coexists in the same geographical area, suggesting that genetic exchanges among these strains may occur more frequently than initially expected.     

Genetic divergence of Toxoplasma gondii strains associated with ocular toxoplasmosis, Brazil

Previous studies have shown a high prevalence of toxoplasmosis and the frequent occurrence of ocular disease in Brazil. To identify the genotypes of parasite strains associated with ocular disease, we compared 25 clinical and animal isolates of Toxoplasma gondii from Brazil to previously characterized clonal lineages from North America and Europe. Multilocus nested polymerase chain reaction analysis was combined with direct sequencing of a polymorphic intron to classify strains by phylogenetic methods. The genotypes of T. gondii strains isolated from Brazil were highly divergent when compared to the previously described clonal lineages. Several new predominant genotypes were identified from different regions of Brazil, including 2 small outbreaks attributable to foodborne or waterborne infection. These findings show that the genetic makeup of T. gondii is more complex than previously recognized and suggest that unique or divergent genotypes may contribute to different clinical outcomes of toxoplasmosis in different localities.     

Molecular genotyping of Toxoplasma gondii from Central and South America revealed high diversity within and between populations

Recent population studies revealed that a few major clonal lineages of Toxoplasma gondii dominate in different geographical regions. The Type II and III lineages are widespread in all continents and dominate in Europe, Africa and North America. In addition, the type 12 lineage is the most common type in wildlife in North America, the Africa 1 and 3 are among the major types in Africa, and ToxoDB PCR-RFLP #9 is the major type in China. Overall the T. gondii strains are more diverse in South America than any other regions. Here, we analyzed 164 T. gondii isolates from three countries in Central America (Guatemala, Nicaragua, Costa Rica), from one country in Caribbean (Grenada) and five countries from South America (Venezuela, Colombia, Peru, Chile, and Argentina). The multilocous polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) based genotyping of 11 polymorphic markers (SAG1, SAG2, alt.SAG2, SAG3, BTUB, GRA6, L358, PK1, C22-8, C29-2 and Apico) were applied to 148 free-range chicken (Gallus domesticus) isolates and 16 isolates from domestic cats (Felis catus) in Colombia; 42 genotypes were identified. Linkage disequilibrium analysis indicated more frequent genetic recombination in populations of Nicaragua and Colombia, and to a lesser degree in populations of Costa Rica and Argentina. Bayesian structural analysis identified at least three genetic clusters, and phylogenetic network analysis identified four major groups. The ToxoDB PCR-RFLP #7, Type III and II were major lineages identified from Central and South America, with high frequencies of the closely related ToxoDB PCR-RFLP #7 and Type III lineages. Taken together, this study revealed high diversity within and between T. gondii populations in Central and South America, and the dominance of Type III and its closely related ToxoDB PCR-RFLP #7 lineages.     

Properties of case/pseudocontrol analysis for genetic association studies: Effects of recombination, ascertainment, and multiple affected offspring

The case/pseudocontrol approach is a general framework for family-based association analysis, incorporating several previously proposed methods such as the transmission/disequilibrium test and log-linear modelling of parent-of-origin effects. In this report, I examine the properties of methods based on a case/pseudocontrol approach when applied to a linked marker rather than (or in addition to) the true disease locus or loci, and when applied to sibships that have been ascertained on, or that may simply contain, multiple affected sibs. Through simulations and analytical calculations, I show that the expected values of the observed relative risk parameters (estimating quantities such as effects due to a child's own genotype, maternal genotype, and parent-of-origin) depend crucially on the ascertainment scheme used, as well as on whether there is non-negligible recombination between the true disease locus and the locus under study. In the presence of either recombination or ascertainment on multiple affected offspring, methods based on conditioning on parental genotypes are shown to give unbiased genotype relative risk estimates at the true disease locus (or loci) but biased estimates of population genotype relative risks at a linked marker, suggesting that the resulting estimates may be misleading when used to predict the power of future studies. Methods that allow for exchangeability of parental genotypes are shown (in the presence of either recombination or ascertainment on multiple affected offspring) to produce false-positive evidence of maternal genotype effects when there are true parent-of-origin or mother-child interaction effects, even when analyzing the true locus. These results suggest that care should be taken in both the interpretation and application of parameter estimates obtained from family-based genetic association studies.     

Wide distribution of a high-virulence Borrelia burgdorferi clone in Europe and North America

The A and B clones of Borrelia burgdorferi sensu stricto, distinguished by outer surface protein C (ospC) gene sequences, are commonly associated with disseminated Lyme disease. To resolve phylogenetic relationships among isolates, we sequenced 68 isolates from Europe and North America at 1 chromosomal locus (16S-23S ribosomal RNA spacer) and 3 plasmid loci (ospC,dbpA, and BBD14). The ospC-A clone appeared to be highly prevalent on both continents, and isolates of this clone were uniform in DNA sequences, which suggests a recent trans-oceanic migration. The genetic homogeneity of ospC-A isolates was confirmed by sequences at 6 additional chromosomal housekeeping loci (gap, alr, glpA, xylB, ackA, and tgt). In contrast, the ospC-B group consists of genotypes distinct to each continent, indicating geographic isolation. We conclude that the ospC-A clone has dispersed rapidly and widely in the recent past. The spread of the ospC-A clone may have contributed, and likely continues to contribute, to the rise of Lyme disease incidence.     

Transmission dynamics of Toxoplasma gondii on a pig farm.

Transmission of Toxoplasma gondii infection on a pig farm in New England was investigated using genetic and ecological methods to (i) determine if infection of pigs was a result of a single source, such as in an epizootic situation (e.g. outbreak) or of multiple sources, such as in an enzootic situation, (ii) identify the main source species of infection to pigs and (iii) evaluate the role of the environment surrounding the farm as the source of infection on the farm. Genetic characterization of 25 T. gondii isolates from market pigs revealed three distinct genotypes with no evidence of recombinants. These data imply that at least three distinct exposure events occurred during the 7-month lifespan of these pigs. This genotype diversity is consistent with enzootic transmission of T. gondii on the farm. Cats were suspected as the main source of pig infection based on the high seroprevalence (>95%) in pigs. The presence of the two most common T. gondii genotypes in eight isolates from free ranging chickens on this farm corroborated the role of cats because chickens were probably infected through ingestion of oocysts in the soil. The seroprevalence of toxoplasmosis in 163 wild mammals and birds captured around the pig sties (overall 13.1%) increased with proximity to the pig sties. Thus, transmission of T. gondii was higher near the pig sties than in the surrounding environment probably because of increased density of oocysts there. We propose that the farm does not simply reflect its surroundings in terms of strain composition and risk of infection, but that it acts as a reservoir of strains from which the outflow of new infections into its surrounding environment is higher than the inflow.     

应用VNTR技术对绵阳地区结核分枝杆菌进行基因分型

目的应用数目可变串联重复序列(variable number of tandem repeats,VNTRs)分析技术探讨绵阳地区结核分枝杆菌的基因分型,为结核病的防治提供科学依据。方法选取绵阳地区结核分枝杆菌临床分离菌株,采用聚合酶链反应(PCR)分别对VNTR位点进行检测分析,应用BioNumerics 5.0软件进行聚类分析。结果共对79株结核分枝杆菌的15个VNTR位点进行了检测,结果显示明显的基因多态性。各个位点的分辨能力不同,QUB-11b位点的分辨指数(Hunter-Gaston index,HGI)最高,为0.880;ETR C的HGI最低,为0.234;总HGI达到1。经聚类分析,79株菌可分为7个基因群、79个基因型,以Ⅲ群和Ⅴ群所占比例较大,Ⅲ群含37株菌,占46.8%;Ⅴ群含33株菌,占41.8%,其他菌株呈散在分布。结论绵阳地区结核分枝杆菌VNTRs基因存在明显的多态性,主要流行菌群为Ⅲ群和Ⅴ群,应加强此两群菌株的监控。     

Joint genotype‐ and ancestry‐based genome‐wide association studies in admixed populations

In genome‐wide association studies (GWAS) genetic loci that influence complex traits are localized by inspecting associations between genotypes of genetic markers and the values of the trait of interest. On the other hand, admixture mapping, which is performed in case of populations consisting of a recent mix of two ancestral groups, relies on the ancestry information at each locus (locus‐specific ancestry). Recently it has been proposed to jointly model genotype and locus‐specific ancestry within the framework of single marker tests. Here, we extend this approach for population‐based GWAS in the direction of multimarker models. A modified version of the Bayesian information criterion is developed for building a multilocus model that accounts for the differential correlation structure due to linkage disequilibrium (LD) and admixture LD. Simulation studies and a real data example illustrate the advantages of this new approach compared to single‐marker analysis or modern model selection strategies based on separately analyzing genotype and ancestry data, as well as to single‐marker analysis combining genotypic and ancestry information. Depending on the signal strength, our procedure automatically chooses whether genotypic or locus‐specific ancestry markers are added to the model. This results in a good compromise between the power to detect causal mutations and the precision of their localization. The proposed method has been implemented in R and is available at http://www.math.uni.wroc.pl/~mbogdan/admixtures/ .     

Comparison of admixture and association mapping in admixed families

The family-based admixture mapping test (AMT) identifies disease-related genes using family data from admixed individuals with the disease of interest (cases). The cases' genotypes at a set of markers are used to infer their DNA ancestry as it varies in blocks along the chromosomes. The test compares the cases' inferred ancestries to those expected from their family histories. Deviation between observed and expected ancestries in a region suggests the presence of a disease gene. We use a likelihood-based development of the AMT to compare it with the transmission disequilibrium test (TDT) as applied to admixed populations. The two tests have a common framework but differ significantly when the disease locus is untyped. The TDT infers disease-locus genotypes using the markers with which it is in linkage disequilibrium (LD). In contrast, the AMT infers disease locus ancestries using those of its linked markers. Thus, TDT power depends on LD between disease and marker loci, while AMT power depends on the lengths of the ancestry blocks containing the disease locus. We compare the power of the two tests when applied to cases with descent from two ancestral populations. The AMT outperforms the TDT when case marker ancestries are correctly specified and LD between disease and marker loci is less than one-third its maximal value (Delta' < 1/3). However, the TDT performs better in the presence of uncertain marker ancestries, even for weak LD between disease and marker loci (Delta' = 0.1). These findings have implications for the design of studies using admixed populations.     

Discernible but limited introgression has occurred where Trichinella nativa and the T6 genotype occur in sympatry

The genetic diversity within and among parasite populations provides clues to their evolutionary history. Here, we sought to determine whether mitochondrial and microsatellite DNA variation could be used to evaluate the extent of differentiation, gene flow and historical reproductive isolation among the freeze resistant parasites Trichinella nativa and the Trichinella T6 genotype infecting wolverines (Gulo gulo) in Nunavut, Canada. To this end, we genotyped Trichinella isolates derived from the diaphragms of 39 wolverines from this locale to reference strains of T. nativa and the Trichinella T6 genotype. Results showed that among a subset of 13 isolates examined, individuals resembled T. nativa in their mitochondrial DNA, but resembled the Trichinella T6 genotype when assayed at expansion segment V and the internal transcribed spacer of the nuclear rDNA. To adjudicate among these conflicting diagnoses, we further characterized each isolate at several nuclear microsatellite loci and again compared these to data from reference strains. Statistical assignment established that the nuclear genomes of most Nunavut isolates corresponded to those of the Trichinella T6 genotype; however, two isolates corresponded to T. nativa, and one isolate exhibited equal similarity to both reference strains. Taken as a whole, the evidence suggests that these isolates derive from the T. nativa matrilineage, but that their nuclear genomes resemble individuals previously designated as Trichinella T6. Assuming distinct lineages, this argues for cross-hybridization among these genotypes. Although introgression has occurred, recognizable genetic distinctions persist. One possibility is that selection disfavors the survival of hybrid offspring in most instances. Alternatively, the recent disappearance of glacial barriers might have increased contact, and therefore introgression. Broader geographic sampling will be required to determine the extent to which hybridization occurs beyond this particular geographic focus.     

The population structure of drug-resistant Mycobacterium tuberculosis clinical isolates from Sichuan in China

China ranks second next to India among 22 high-burden countries despite decades' effort on tuberculosis (TB) control. The Sichuan province today contains the second-largest number of TB cases among Chinese provinces, where the prevalence of drug-resistant TB, especially MDR-TB, is much higher than the average level in eastern China. In this study, the population structure and the transmission characteristics of drug-resistant TB in Sichuan province were studied by spoligotyping and 24-locus Mycobacterial interspersed repetitive units-variable number tandem DNA repeats (MIRU-VNTR), applied to a total of 306 clinical isolates. Spoligotyping-based analysis showed that Beijing family represented 69.28% of all isolates and constituted the largest group (66.24%) of MDR-TB in Sichuan. The remaining isolates, accounting for 33.76% of MDR isolates, belonged to the ill-defined T family, Manu2, H3, LAM9, and other minor unassigned clades. The discriminatory power evaluated for spoligotyping was poor (HGI=0.595), but high for 24-locus MIRU-VNTRs (HGI=0.999). The number of the most discriminatory loci (h>0.6) was 12, including locus 424, 802, 960, 1644, 1955, 2163b, 2996, 3007, 3192, 3690, 4348 and 4052. It was concluded that 24-locus MIRU-VNTRs could be a more discriminatory tool for differentiating clinical isolates from Sichuan region. The small clustering size obtained from the current population structure analysis suggested that the high prevalence of drug-resistant TB in this region might be attributed partially to the acquired resistance due to inappropriate drug use rather than active transmission of drug-resistant TB (primary resistance).     

Molecular characterization of Cryptosporidium spp. from wild rats and mice from rural communities in the Philippines

In order to examine the prevalence of Cryptosporidium in wild rodents in the Philippines and understand the role wild rodents play in the transmission of this parasite to humans and livestock, 194 fecal samples from wild rats and mice from Luzon and Mindoro islands were examined. Molecular screening at the 18S and actin gene loci identified an overall prevalence of 25.8% (95%CI: 19.8, 32.5). Sequence and phylogenetic analysis of both loci identified C. parvum, C. muris, C. scrofarum, rat genotypes I-IV and a C. suis-like genotype in the rat-derived isolates and is the first report of C. suis-like and C. scrofarum in rats. Mixed infections were identified in 24% of the Cryptosporidium positive isolates. Rat genotypes II, III and IV showed high intragenotypic variation at the 18S gene locus compared to the actin locus.     

梅州地区TLR4基因连锁不平衡与缺血性脑卒中相关性研究

目的 检测梅州地区缺血性脑卒中患者和健康对照人群TLR4基因的（rs10759932、rs11536879、rs11536891、rs1927914）的基因多态性，经连锁不平衡分析其与缺血性脑卒中的相关性。方法 收集2018年1月1日 - 2018年7月31日住院治疗的突发缺血性脑卒中患者作为病例组，同期在体检中心收集健康人群作为对照组。应用Massarray SNP 分型技术检测两组患者TLR4基因的4个位点基因型，进行Hardy - Weinberg（H - W）平衡检测。采用不同模型分析上述位点不同基因型与脑梗死发病风险的相关性，并通过连锁不平衡分析其与缺血性脑卒中的相关性。结果 病例组纳入病例186名，对照组纳入健康人194名；4种SNP位点均符合H - W 平衡。rs1927914位点G/G基因型在对照组出现频率远远高于病例组（χ2 = 9.267，P＜0.05）。rs10759932位点T/T基因型在女性对照组中出现的频率显著高于男性[OR = 0.38 (0.18 - 0.81)，P＜0.05]。4个SNP位点之间均存在连锁不平衡，TGCG基因型组合在缺血性卒中男性患者出现的频率显著高于女性[OR = 3.54 (1.17 - 10.69)，P＜0.05]。结论 梅州地区rs1927914位点A＞G为保护性基因突变，可以降低缺血性脑卒中发生。4个位点的连锁不平衡与缺血性脑卒中发生存在部分性别差异，TGCG组合为男性人群的危险基因，脑卒中发生率显著升高。     

Multidrug-resistant tuberculosis, People's Republic of China, 2007-2009

We conducted a case-control study to investigate risk factors for multidrug-resistant tuberculosis (MDR TB) in the People's Republic of China. Genotyping analysis was used to estimate the percentage of cases from recent transmission among 100 MDR TB case-patients hospitalized during April 2007-July 2009. Molecular subtyping of isolates showed that 41% of MDR TB strains clustered. Beijing genotype was found in 94% of the MDR TB isolates and 79% of the pan-susceptible isolates. In multivariate analysis, MDR TB was independently associated with Beijing genotype, retreatment for TB, symptoms lasting >3 months before first evaluation at the hospital, lack of health insurance, and being a farmer (vs. being a student). MDR TB was associated with Beijing genotype and lower socioeconomic status. A large percentage of MDR TB cases seemed to result from recent transmission. Early detection, effective treatment, and infection control measures for MDR TB are needed to reduce transmission.     

A spatial probit model for fine-scale mapping of disease genes

We present a novel statistical method for linkage disequilibrium (LD) mapping of disease susceptibility loci in case-control studies. Such studies exploit the statistical correlation or LD that exist between variants physically close along the genome to identify those that correlate with disease status and might thus be close to a causative mutation, generally assumed unobserved. LD structure, however, varies markedly over short distances because of variation in local recombination rates, mutation and genetic drift among other factors. We propose a Bayesian multivariate probit model that flexibly accounts for the local spatial correlation between markers. In a case-control setting, we use a retrospective model that properly reflects the sampling scheme and identify regions where single- or multi-locus marker frequencies differ across cases and controls. We formally quantify these differences using information-theoretic distance measures while the fully Bayesian approach naturally accommodates unphased or missing genotype data. We demonstrate our approach on simulated data and on real data from the CYP2D6 region that has a confirmed role in drug metabolism.     

Case/pseudocontrol analysis in genetic association studies: A unified framework for detection of genotype and haplotype associations, gene-gene and gene-environment interactions, and parent-of-origin effects

Estimation and testing of genetic effects (genotype relative risks) are often performed conditionally on parental genotypes, using data from case-parent trios. This strategy avoids having to estimate nuisance parameters such as parental mating type frequencies, and also avoids generating spurious results due to confounding causes of association such as population stratification. For effects at a single locus, the resulting analysis is equivalent to matched case/control analysis via conditional logistic regression, using the case and three "pseudocontrols" derived from the untransmitted parental alleles. We previously showed that a similar approach can be used for analyzing genotype and haplotype effects at a set of closely linked loci, but with a required adjustment to the conditioning argument that results in varying numbers of pseudocontrols, depending on the disease model that is to be fitted. Here we extend this method to include the analysis of epistatic effects (gene-gene interactions) at unlinked loci, to include parent-of-origin effects at one or more loci, and to allow additional incorporation of gene-environment interactions. The conditional logistic approach provides a natural and flexible framework for incorporating these additional effects. By relaxing the conditioning on parental genotypes to allow exchangeability of parental genotypes, we show how the power of this approach can be increased when studying parent-of-origin effects. Simulations suggest that there is limited power to distinguish between parent-of-origin effects and effects due to interaction between genotypes of mother and child.     

First genetic characterization of Toxoplasma gondii infection in poultry meat intended for human consumption in eastern China

Toxoplasma gondii can infect nearly all warm-blooded animals including birds. However, limited information on the molecular epidemiology and genotypes of T. gondii infecting poultry is available in China. Therefore, the present study characterized T. gondii genotypes in poultry meat in eastern China. During August 2015 and September 2016, muscle tissue samples collecting from 414 poultries (257 chickens, 115 ducks and 42 geese) in Shandong provinces were used to detect the T. gondii B1 gene by a semi-nested PCR, and the positive samples were genotyped at 10 nuclear loci (i.e., SAG1, alternative SAG2, 5′- and 3′-SAG2, SAG3, L358, BTUB, c22-8, GRA6, c29-2, PK1) and an apicoplast locus Apico by multi-locus PCR-RFLP technology. Thirty-two (7.37%) samples from 414 poultry meat were T. gondii B1 gene positive. Chicken had the highest T. gondii prevalence (8.17%), followed by ducks (7.83%) and geese (4.76%). Furthermore, only one genotype (ToxoDB#9) was identified. This is the first genetic characterization of T. gondii isolates from poultry meat in Shandong province, eastern China and also the first report of genotype ToxoDB#9 was found in poultry in China, which provide basic data for the surveillance and control of T. gondii infection in poultry, other animals and humans.