A simple and robust TDT-type test against genotyping error with error rates varying across families

The transmission/disequilibrium test (TDT), a family based test of linkage and association, is a popular test for studies of complex inheritance, as it is nonparametric and robust against spurious conclusions induced by hidden genetic structure, such as stratification or admixture. However, the TDT may be biased by genotyping errors. Undetected genotyping errors may be contributing to an inflated type I error rate among reported TDT-derived associations. To adjust for bias, a popular approach is to assume a genotype error model for describing the pattern of errors and propose association tests using likelihood method. However, all model-based approaches tend to perform unsatisfactorily if the related genotyping error rates are not identical across all families. In this paper, we propose a TDT-type association test which is not only simple, robust against population stratification (and hence the assumption of Hardy-Weinberg equilibrium is not required), but also robust against genotyping error with error rates varying across families. Simulation studies confirm that the new test has very reasonable performance.     

Application and interpretation of transmission/disequilibrium tests: transmission of HLA-DQ haplotypes to unaffected siblings in 526 families with type 1 diabetes

It is widely believed that, if a genetic marker shows a transmission distortion in patients by the transmission/disequilibrium test (TDT), then a transmission distortion in healthy siblings would be seen in the opposite direction. This is also the case in a complex disease. Furthermore, it has been suggested that replacing the McNemar statistics of the TDT with a test of heterogeneity between transmissions to affected and unaffected children could increase the power to detect disease association. To test these two hypotheses empirically, we analyzed the transmission of HLA-DQA1-DQB1 haplotypes in 526 Norwegian families with type 1 diabetic children and healthy siblings, since some DQA1-DQB1 haplotypes represent major genetic risk factors for type 1 diabetes. Despite the strong positive and negative disease associations with particular DQ haplotypes, we observed no significant deviation from 50% for transmission to healthy siblings. This could be explained by the low penetrance of susceptibility alleles, together with the fact that IDDM loci also harbor strongly protective alleles that can override the risk contributed by other loci. Our results suggest that, in genetically complex diseases, detectable distortion in transmission to healthy siblings should not be expected. Furthermore, the original TDT seems more powerful than a heterogeneity test.     

Evaluating the evidence for transmission distortion in human pedigrees

Children of a heterozygous parent are expected to carry either allele with equal probability. Exceptions can occur, however, due to meiotic drive, competition among gametes, or viability selection, which we collectively term "transmission distortion" (TD). Although there are several well-characterized examples of these phenomena, their existence in humans remains unknown. We therefore performed a genome-wide scan for TD by applying the transmission disequilibrium test (TDT) genome-wide to three large sets of human pedigrees of European descent: the Framingham Heart Study (FHS), a founder population of European origin (HUTT), and a subset of the Autism Genetic Resource Exchange (AGRE). Genotyping error is an important confounder in this type of analysis. In FHS and HUTT, despite extensive quality control, we did not find sufficient evidence to exclude genotyping error in the strongest signals. In AGRE, however, many signals extended across multiple SNPs, a pattern highly unlikely to arise from genotyping error. We identified several candidate regions in this data set, notably a locus in 10q26.13 displaying a genome-wide significant TDT in combined female and male transmissions and a signature of recent positive selection, as well as a paternal TD signal in 6p21.1, the same region in which a significant TD signal was previously observed in 30 European males. Neither region replicated in FHS, however, and the paternal signal was not visible in sperm competition assays or as allelic imbalance in sperm. In maternal transmissions, we detected no strong signals near centromeres or telomeres, the regions predicted to be most susceptible to female-specific meiotic drive, but we found a significant enrichment of top signals among genes involved in cell junctions. These results illustrate both the potential benefits and the challenges of using the TDT to study transmission distortion and provide candidates for investigation in future studies.     

Association and linkage analysis of RGS4 polymorphisms with schizophrenia and bipolar disorder in Brazil

Linkage and association studies in five independently ascertained samples have suggested that polymorphisms of the regulator of G-protein signaling 4 (RGS4) may confer risk for schizophrenia (SCZ). Suggestive evidence for association with bipolar disorder (BD) has also been presented. However, the associated alleles and haplotypes have differed among the samples. Data from other independent samples may clarify the putative associations. Hence, we investigated an independent, ethnically diverse Brazilian population comprising patients with SCZ (n=271) or BD1 (n=306), who were contrasted with 576 community-based controls. Parents of 49 SCZ cases and 44 BD cases were available for transmission disequilibrium tests (TDTs). Four RGS4 single-nucleotide polymorphisms (SNPs) 1, 4, 7 and 18 putatively associated with SCZ were investigated. In the SCZ samples, significant case-control differences were not observed for individual SNPs or haplotypes, though the TDT suggested transmission distortion similar to that observed in the initial report. For the BD sample, case-control comparisons revealed no significant differences for individual SNPs, but an omnibus test suggested differences in the overall distribution of haplotypes bearing all four SNPs (SNP-EM Omnibus likelihood ratio test; P=0.003). The TDT revealed over-transmission of allele A at SNP7 (P=0.016), as well as haplotypes incorporating this allele. However, global tests incorporating all haplotypes yielded only suggestive trends for association (P=0.19). In conclusion, association with SCZ was not detected in the present analyses. The failure to detect an association may be related to inadequate power or to confounds related to ethnic admixture. Suggestive associations with BD detected here require further investigation in a larger sample.     

The study of CTLA-4 and vitamin D receptor polymorphisms in the Romanian type 1 diabetes population

Several studies suggested that part of the genetic susceptibility for Type 1 diabetes (TIDM) is encoded by some polymorphisms of CTLA-4 gene (2q33) and of Vitamin D Receptor gene (VDR; 12q12-14). Our aim was to assess their contribution to TIDM genetic susceptibility in the Romanian population. We typed CTLA-4 49 A/G and VDR Fok I (F/f), Apa I (A/a) and Taq I (T/t) polymorphisms by Sequence Specific Primer PCR (SSP-PCR) in 204 Romanian diabetic families (756 individuals: 212 TIDM probands and 544 unaffected parents and siblings). We studied alleles transmission using the Transmission Disequilibrium Test (TDT). We found an increased transmission of CTLA-4 49G allele to diabetics (54.8%, p=0.11). The transmission of F (56.1%, p=0.063), a (55.7%, p=0.061) and T (51.8%, p=0.37) alleles of VDR gene to diabetics was increased but did not reach statistical significance. In conclusion we found the same increased transmission of CTLA-4 49 G allele to diabetics as previously reported. VDR Foq I F allele seems to be predisposing while Taq I T allele seems to be protective.     

Combining evidence for association from transmission disequilibrium and case-control studies using single-nucleotide polymorphisms

The aim of the present analysis is to combine evidence for association from the two most commonly used designs in genetic association analysis, the case-control design and the transmission disequilibrium test (TDT) design. The cases here are affected offspring from nuclear families and are used in both the case-control and TDT designs. As a result, inference from these designs is not independent. We applied a simple logistic regression method for combining evidence for association from case-control and TDT designs to single-nucleotide polymorphism data purchased on a region on chromosome 3, replicate 1 of the Aipotu population. Combining the evidence from the case-control and TDT designs yielded a 5-10% reduction in the standard errors of the relative risk estimates. The authors did not know the results before the analyses were conducted.     

Marker selection for the transmission/disequilibrium test, in recently admixed populations.

Recent admixture between genetically differentiated populations can result in high levels of association between alleles at loci that are <=10 cM apart. The transmission/disequilibrium test (TDT) proposed by Spielman et al. (1993) can be a powerful test of linkage between disease and marker loci in the presence of association and therefore could be a useful test of linkage in admixed populations. The degree of association between alleles at two loci depends on the differences in allele frequencies, at the two loci, in the founding populations; therefore, the choice of marker is important. For a multiallelic marker, one strategy that may improve the power of the TDT is to group marker alleles within a locus, on the basis of information about the founding populations and the admixed population, thereby collapsing the marker into one with fewer alleles. We have examined the consequences of collapsing a microsatellite into a two-allele marker, when two founding populations are assumed for the admixed population, and have found that if there is random mating in the admixed population, then typically there is a collapsing for which the power of the TDT is greater than that for the original microsatellite marker. A method is presented for finding the optimal collapsing that has minimal dependence on the disease and that uses estimates either of marker allele frequencies in the two founding populations or of marker allele frequencies in the current, admixed population and in one of the founding populations. Furthermore, this optimal collapsing is not always the collapsing with the largest difference in allele frequencies in the founding populations. To demonstrate this strategy, we considered a recent data set, published previously, that provides frequency estimates for 30 microsatellites in 13 populations.     

Identification, Inheritance, and Variation of Microsatellite Markers in the Black Scallop Mimachlamys varia

Five polymorphic microsatellite loci were identified in the black scallop Mimachlamys varia after construction of a genomic library enriched for (GT)n. To examine the transmission pattern of microsatellite alleles, several families were created and genotypes scored for three loci. The expected Mendelian ratios were found in 12 of 14 segregations examined. Unexpected segregations may be explained by a genotyping error (allelic dropout), given that when a specific allele was treated as dominant, the phenotypic ratios conformed to Mendelian expectations. The five loci were also examined in two samples from the Spanish coast. The two localities displayed similar mean values for the number of alleles per locus (7.2-8.4), allelic richness (7.2-7.9), and observed (0.389-0.484) and expected heterozygosity (0.545-0.618). Significant Hardy-Weinberg deviations were observed at three loci, with heterozygote deficiency occurring in all cases. Global multilocus θ value and allele frequencies at one locus revealed significant differentiation between the two localities.     

HLA and genomewide allele sharing in dizygotic twins

Gametic selection during fertilization or the effects of specific genotypes on the viability of embryos may cause a skewed transmission of chromosomes to surviving offspring. A recent analysis of transmission distortion in humans reported significant excess sharing among full siblings. Dizygotic (DZ) twin pairs are a special case of the simultaneous survival of two genotypes, and there have been reports of DZ pairs with excess allele sharing around the HLA locus, a candidate locus for embryo survival. We performed an allele-sharing study of 1,592 DZ twin pairs from two independent Australian cohorts, of which 1,561 pairs were informative for linkage on chromosome 6. We also analyzed allele sharing in 336 DZ twin pairs from The Netherlands. We found no evidence of excess allele sharing, either at the HLA locus or in the rest of the genome. In contrast, we found evidence of a small but significant (P=.003 for the Australian sample) genomewide deficit in the proportion of two alleles shared identical by descent among DZ twin pairs. We reconciled conflicting evidence in the literature for excess genomewide allele sharing by performing a simulation study that shows how undetected genotyping errors can lead to an apparent deficit or excess of allele sharing among sibling pairs, dependent on whether parental genotypes are known. Our results imply that gene-mapping studies based on affected sibling pairs that include DZ pairs will not suffer from false-positive results due to loci involved in embryo survival.     

Importance of a pilot study for non-invasive genetic sampling: genotyping errors and population size estimation in red deer

Population size information is critical for managing endangered or harvested populations. Population size can now be estimated from non-invasive genetic sampling. However, pitfalls remain such as genotyping errors (allele dropout and false alleles at microsatellite loci). To evaluate the feasibility of non-invasive sampling (e.g., for population size estimation), a pilot study is required. Here, we present a pilot study consisting of (i) a genetic step to test loci amplification and to estimate allele frequencies and genotyping error rates when using faecal DNA, and (ii) a simulation step to quantify and minimise the effects of errors on estimates of population size. The pilot study was conducted on a population of red deer in a fenced natural area of 5440 ha, in France. Twelve microsatellite loci were tested for amplification and genotyping errors. The genotyping error rates for microsatellite loci were 0–0.83 (mean=0.2) for allele dropout rates and 0–0.14 (mean=0.02) for false allele rates, comparable to rates encountered in other non-invasive studies. Simulation results suggest we must conduct 6 PCR amplifications per sample (per locus) to achieve approximately 97% correct genotypes. The 3% error rate appears to have little influence on the accuracy and precision of population size estimation. This paper illustrates the importance of conducting a pilot study (including genotyping and simulations) when using non-invasive sampling to study threatened or managed populations.     

Handling Missing Data in Transmission Disequilibrium Test in Nuclear Families with One Affected Offspring

The Transmission Disequilibrium Test (TDT) compares frequencies of transmission of two alleles from heterozygote parents to an affected offspring. This test requires all genotypes to be known from all members of the nuclear families. However, obtaining all genotypes in a study might not be possible for some families, in which case, a data set results in missing genotypes. There are many techniques of handling missing genotypes in parents but only a few in offspring. The robust TDT (rTDT) is one of the methods that handles missing genotypes for all members of nuclear families [with one affected offspring]. Even though all family members can be imputed, the rTDT is a conservative test with low power. We propose a new method, Mendelian Inheritance TDT (MITDT-ONE), that controls type I error and has high power. The MITDT-ONE uses Mendelian Inheritance properties, and takes population frequencies of the disease allele and marker allele into account in the rTDT method. One of the advantages of using the MITDT-ONE is that the MITDT-ONE can identify additional significant genes that are not found by the rTDT. We demonstrate the performances of both tests along with Sib-TDT (S-TDT) in Monte Carlo simulation studies. Moreover, we apply our method to the type 1 diabetes data from the Warren families in the United Kingdom to identify significant genes that are related to type 1 diabetes.     

Transmission disequilibrium test with discordant sib pairs when parents are available

The transmission disequilibrium test (TDT) has been employed to map disease susceptibility loci (DSL), while being immune to the problem of population admixture. The customary TDT test (TDT(D)) was developed for affected child(ren) and their parents and was most often applied to case-parent trios. Recently, the TDT has been extended to the situations when (1) parents are not available but affected and nonaffected sibs from each family are available, (2) unrelated control-parent trios are available for combined analyses with case-parent trios (TDT(DC)), and (3) large pedigrees. For many diseases, affected children in the case-parent trios enlisted into the TDT(D) have unaffected sibs who can be recruited. We present an extension of the TDT by effectively incorporating one unaffected sib of each of the affected children in the case-parent trios into a single analysis (TDT(DS), where DS denotes discordant sib pairs). We have developed a general analytical method for computing the statistical power of the TDT(DS) under any genetic model, the accuracy of which is validated by computer simulations. We compare the power of the TDT(D), TDT(DC), and TDT(DS) under a range of parameter space and genetic models. We find that the TDT(DS) is generally more powerful than the TDT(DC) and TDT(D), particularly when the disease is prevalent (>30%) in the population. The relative power of the TDT(D) and the TDT(DS) largely depends upon the allele frequencies and genetic effects at the DSL, whereas the recombination rate, the degree of linkage disequilibrium, and the marker allele frequencies have little effect. Importantly, the TDT(DS) not only may be more powerful, it also has the advantage of being able to test for segregation distortion that may yield false linkage/association in the TDT(D).     

Clonal diversity of a lizard malaria parasite, Plasmodium mexicanum, in its vertebrate host, the western fence lizard: role of variation in transmission intensity over time and space

Within the vertebrate host, infections of a malaria parasite (Plasmodium) could include a single genotype of cells (single-clone infections) or two to several genotypes (multiclone infections). Clonal diversity of infection plays an important role in the biology of the parasite, including its life history, virulence, and transmission. We determined the clonal diversity of Plasmodium mexicanum, a lizard malaria parasite at a study region in northern California, using variable microsatellite markers, the first such study for any malaria parasite of lizards or birds (the most common hosts for Plasmodium species). Multiclonal infections are common (50-88% of infections among samples), and measures of genetic diversity for the metapopulation (expected heterozygosity, number of alleles per locus, allele length variation, and effective population size) all indicated a substantial overall genetic diversity. Comparing years with high prevalence (1996-1998 = 25-32% lizards infected), and years with low prevalence (2001-2005 = 6-12%) found fewer alleles in samples taken from the low-prevalence years, but no reduction in overall diversity (H = 0.64-0.90 among loci). In most cases, rare alleles appeared to be lost as prevalence declined. For sites chronically experiencing low transmission intensity (prevalence approximately 1%), overall diversity was also high (H = 0.79-0.91), but there were fewer multiclonal infections. Theory predicts an apparent excess in expected heterozygosity follows a genetic bottleneck. Evidence for such a distortion in genetic diversity was observed after the drop in parasite prevalence under the infinite alleles mutation model but not for the stepwise mutation model. The results are similar to those reported for the human malaria parasite, Plasmodium falciparum, worldwide, and support the conclusion that malaria parasites maintain high genetic diversity in host populations despite the potential for loss in alleles during the transmission cycle or during periods/locations when transmission intensity is low.     

A multistage testing strategy for detection of quantitative trait Loci affecting disease resistance in Atlantic salmon

A multistage testing strategy to detect QTL for resistance to infectious salmon anemia (ISA) in Atlantic salmon is proposed. First, genotyping of amplified fragment length polymorphisms (AFLP) and a transmission disequilibrium test (TDT) were carried out using dead offspring from a disease resistance challenge test. Second, AFLP genotyping among survivors followed by a Mendelian segregation test was performed. Third, within-family survival analyses using all offspring were developed and applied to significant TDT markers with Mendelian inheritance. Maximum-likelihood methodology was developed for TDT with dominant markers to exploit linkage disequilibrium within families. The strategy was tested with two full-sib families of Atlantic salmon sired by the same male and consisting of 79 offspring in total. All dead offspring from the two families were typed for 64 primer combinations, resulting in 340 scored markers. There were 26 significant results out of 401 TDTs using dead offspring. In the second stage, only 17 marker families showed Mendelian segregation and were tested in survival analysis. A permutation test was performed for all survival analyses to compute experimentwise P-values. Two markers, aaccac356 and agccta150, were significant at P < 0.05 when accounting for multiple testing in the survival analyses. The proposed strategy might be more powerful than current mapping strategies because it reduces the number of tests to be performed in the last testing stage.     

Sequence variants of the brain-derived neurotrophic factor (BDNF) gene are strongly associated with obsessive-compulsive disorder

We evaluated a possible association between the brain-derived neurotrophic factor (BDNF) gene and susceptibility to obsessive-compulsive disorder (OCD) by genotyping a number of single-nucleotide polymorphisms (SNPs) and one microsatellite marker from the extended BDNF locus in 164 triads with OCD. Extensive background linkage disequilibrium was observed at this locus. Single-locus transmission-distortion tests revealed significant evidence of association with the disease for all the BDNF gene markers tested, including a Val66Met variation affecting the sequence of the proBDNF protein. Analysis of multi-SNP haplotypes provided similar results. Haplotype transmission comparisons in this and previous studies point to a functionally distinct BDNF haplotype uniquely marked by the rare Met66 allele, which is undertransmitted and likely confers a protective effect in OCD and other psychiatric disorders.     

Power calculations for the transmission/disequilibrium and affected sib pair tests using elementary probability methods

The transmission/disequilibrium test (TDT) and the affected sib pair test (ASP) both test for the association of a marker allele with some conditions. Here, we present methods for calculating the probability of detecting the association (power) for a study examining a fixed number of families for suitability for the study and for calculating the number of such families to be examined. Both calculations use a genetic model for the association. The model considered posits a bi-allelic marker locus that is linked to a bi-allelic disease locus with a possibly nonzero recombination fraction between the loci. The penetrance of the disease is an increasing function of the number of disease alleles. The TDT tests whether the transmission by a heterozygous parent of a particular allele at a marker locus to an affected offspring occurs with probability greater than 0.5. The ASP tests whether transmission of the same allele to two affected sibs occurs with probability greater than 0.5. In either case, evidence that the probability is greater than 0.5 is evidence for association between the marker and the disease. Study inclusion criteria (IC) can greatly affect the necessary sample size of a TDT or ASP study. IC considered by us include a randomly selected parent at least one parent or both parents required to be heterozygous. It also allows a specified minimum number of affected offspring to be required (TDT only). We use elementary probability calculations rather than complex mathematical manipulations or asymptotic methods (large sample size approximations) to compute power and requisite sample size for a proposed study. The advantages of these methods are simplicity and generality.     

Transmission ratio distortion in the human body louse, Pediculus humanus (Insecta: Phthiraptera)

We studied inheritance at three microsatellite loci in eight F, and two F2 families of the body (clothes) louse of humans, Pediculus humanus. The alleles of heterozygous female-parents were always inherited in a Mendelian fashion in these families. Alleles from heterozygous male-parents, however, were inherited in two different ways: (i) in a Mendelian fashion and (ii) in a non-Mendelian fashion, where males passed to their offspring only one of their two alleles, that is, 100% nonrandom transmission. In male body lice, where there was non-Mendelian inheritance, the paternally inherited set of alleles was eliminated. We interpret this pattern of inheritance as evidence for extreme transmission ratio distortion of paternal alleles in this species.     

Estimation of genotyping error rate from repeat genotyping,unintentional recaptures and known parent–offspring comparisons in 16 microsatellite loci for brown rockfish (Sebastes auriculatus)

Genotyping errors are present in almost all genetic data and can affect biological conclusions of a study, particularly for studies based on individual identification and parentage. Many statistical approaches can incorporate genotyping errors, but usually need accurate estimates of error rates. Here, we used a new microsatellite data set developed for brown rockfish (Sebastes auriculatus) to estimate genotyping error using three approaches: (i) repeat genotyping 5% of samples, (ii) comparing unintentionally recaptured individuals and (iii) Mendelian inheritance error checking for known parent–offspring pairs. In each data set, we quantified genotyping error rate per allele due to allele drop‐out and false alleles. Genotyping error rate per locus revealed an average overall genotyping error rate by direct count of 0.3%, 1.5% and 1.7% (0.002, 0.007 and 0.008 per allele error rate) from replicate genotypes, known parent–offspring pairs and unintentionally recaptured individuals, respectively. By direct‐count error estimates, the recapture and known parent–offspring data sets revealed an error rate four times greater than estimated using repeat genotypes. There was no evidence of correlation between error rates and locus variability for all three data sets, and errors appeared to occur randomly over loci in the repeat genotypes, but not in recaptures and parent–offspring comparisons. Furthermore, there was no correlation in locus‐specific error rates between any two of the three data sets. Our data suggest that repeat genotyping may underestimate true error rates and may not estimate locus‐specific error rates accurately. We therefore suggest using methods for error estimation that correspond to the overall aim of the study (e.g. known parent–offspring comparisons in parentage studies).     

A prezygotic transmission distorter acting equally in female and male zebra finches Taeniopygia guttata

The two parental alleles at a specific locus are usually inherited with equal probability to the offspring. However, at least three processes can lead to an apparent departure from fair segregation: early viability selection, biased gene conversion and various kinds of segregation distortion. Here, we conduct a genome‐wide scan for transmission distortion in a captive population of zebra finches (Taeniopygia guttata) using 1302 single‐nucleotide polymorphisms (SNPs) followed by confirmatory analyses on independent samples from the same population. In the initial genome‐wide scan, we found significant distortion at three linked loci on chromosome Tgu2 and we were able to replicate this finding in each of two follow‐up data sets [overall transmission ratio = 0.567 (95% CI = 0.536–0.600), based on 1101 informative meioses]. Although the driving allele was preferentially transmitted by both heterozygous females [ratio = 0.560 (95% CI = 0.519–0.603)] and heterozygous males [ratio = 0.575 (95% CI = 0.531–0.623)], we could rule out postzygotic viability selection and biased gene conversion as possible mechanisms. Early postzygotic viability selection is unlikely, because it would result in eggs with no visible embryo and hence no opportunity for genotyping, and we confirmed that both females and males heterozygous for the driving allele did not produce a larger proportion of such eggs than homozygous birds. Biased gene conversion is expected to be rather localized, while we could trace transmission distortion in haplotypes of several megabases in a recombination desert. Thus, we here report the rare case of a prezygotically active transmission distorter operating equally effectively in female and male meioses.