A comparison of DNA pools constructed following whole genome amplification for two-stage SNP genotyping designs

Genotyping in DNA pools reduces the cost and the time required to complete large genotyping projects. The aim of the present study was to evaluate pooling as part of a strategy for fine mapping in regions of significant linkage. Thirty-nine single nucleotide polymorphisms (SNPs) were analyzed in two genomic DNA pools of 384 individuals each and results compared with data after typing all individuals used in the pools. There were no significant differences using data from either 2 or 8 heterozygous individuals to correct frequency estimates for unequal allelic amplification. After correction, the mean difference between estimates from the genomic pool and individual allele frequencies was .033. A major limitation of the use of DNA pools is the time and effort required to carefully adjust the concentration of each individual DNA sample before mixing aliquots. Pools were also constructed by combining DNA after Multiple Displacement Amplification (MDA). The MDA pools gave similar results to pools constructed after careful DNA quantitation (mean difference from individual genotyping .040) and MDA provides a rapid method to generate pools suitable for some applications. Pools provide a rapid and cost-effective screen to eliminate SNPs that are not polymorphic in a test population and can detect minor allele frequencies as low as 1% in the pooled samples. With current levels of accuracy, pooling is best suited to an initial screen in the SNP validation process that can provide high-throughput comparisons between cases and controls to prioritize SNPs for subsequent individual genotyping.

Genomewide linkage study in 1,176 affected sister pair families identifies a significant susceptibility locus for endometriosis on chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite >50 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families (931 from an Australian group and 245 from a U.K. group), each with at least two members--mainly affected sister pairs--with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 (maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 (MLS = 2.09). Minor peaks (with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

Genomewide significant linkage to migrainous headache on chromosome 5q21

Familial typical migraine is a common, complex disorder that shows strong familial aggregation. Using latent-class analysis (LCA), we identified subgroups of people with migraine/severe headache in a community sample of 12,245 Australian twins (60% female), drawn from two cohorts of individuals aged 23-90 years who completed an interview based on International Headache Society criteria. We report results from genomewide linkage analyses involving 756 twin families containing a total of 790 independent sib pairs (130 affected concordant, 324 discordant, and 336 unaffected concordant for LCA-derived migraine). Quantitative-trait linkage analysis produced evidence of significant linkage on chromosome 5q21 and suggestive linkage on chromosomes 8, 10, and 13. In addition, we replicated previously reported typical-migraine susceptibility loci on chromosomes 6p12.2-p21.1 and 1q21-q23, the latter being within 3 cM of the rare autosomal dominant familial hemiplegic migraine gene (ATP1A2), a finding which potentially implicates ATP1A2 in familial typical migraine for the first time. Linkage analyses of individual migraine symptoms for our six most interesting chromosomes provide tantalizing hints of the phenotypic and genetic complexity of migraine. Specifically, the chromosome 1 locus is most associated with phonophobia; the chromosome 5 peak is predominantly associated with pulsating headache; the chromosome 6 locus is associated with activity-prohibiting headache and photophobia; the chromosome 8 locus is associated with nausea/vomiting and moderate/severe headache; the chromosome 10 peak is most associated with phonophobia and photophobia; and the chromosome 13 peak is completely due to association with photophobia. These results will prove to be invaluable in the design and analysis of future linkage and linkage disequilibrium studies of migraine.

A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other

In this report, we describe a simple correction for multiple testing of single-nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) with each other, on the basis of the spectral decomposition (SpD) of matrices of pairwise LD between SNPs. This method provides a useful alternative to more computationally intensive permutation tests. Additionally, output from SNPSpD includes eigenvalues, principal-component coefficients, and factor "loadings" after varimax rotation, enabling the selection of a subset of SNPs that optimize the information in a genomic region.

Zygosity diagnosis in the absence of genotypic data: An approach using latent class analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

A genomewide screen for autism susceptibility loci

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF

Defects in the IFT-B component IFT172 cause Jeune and Mainzer-Saldino syndromes in humans

Intraflagellar transport (IFT) depends on two evolutionarily conserved modules, subcomplexes A (IFT-A) and B (IFT-B), to drive ciliary assembly and maintenance. All six IFT-A components and their motor protein, DYNC2H1, have been linked to human skeletal ciliopathies, including asphyxiating thoracic dystrophy (ATD; also known as Jeune syndrome), Sensenbrenner syndrome, and Mainzer-Saldino syndrome (MZSDS). Conversely, the 14 subunits in the IFT-B module, with the exception of IFT80, have unknown roles in human disease. To identify additional IFT-B components defective in ciliopathies, we independently performed different mutation analyses: candidate-based sequencing of all IFT-B-encoding genes in 1,467 individuals with a nephronophthisis-related ciliopathy or whole-exome resequencing in 63 individuals with ATD. We thereby detected biallelic mutations in the IFT-B-encoding gene IFT172 in 12 families. All affected individuals displayed abnormalities of the thorax and/or long bones, as well as renal, hepatic, or retinal involvement, consistent with the diagnosis of ATD or MZSDS. Additionally, cerebellar aplasia or hypoplasia characteristic of Joubert syndrome was present in 2 out of 12 families. Fibroblasts from affected individuals showed disturbed ciliary composition, suggesting alteration of ciliary transport and signaling. Knockdown of ift172 in zebrafish recapitulated the human phenotype and demonstrated a genetic interaction between ift172 and ift80. In summary, we have identified defects in IFT172 as a cause of complex ATD and MZSDS. Our findings link the group of skeletal ciliopathies to an additional IFT-B component, IFT172, similar to what has been shown for IFT-A.

A homozygous mutation in LTBP2 causes isolated microspherophakia

Microspherophakia is an autosomal-recessive congenital disorder characterized by small spherical lens. It may be isolated or occur as part of a hereditary systemic disorder, such as Marfan syndrome, autosomal dominant and recessive forms of Weill-Marchesani syndrome, autosomal dominant glaucoma–lens ectopia–microspherophakia–stiVness– shortness syndrome, autosomal dominant microspherophakia with hernia, and microspherophakia-metaphyseal dysplasia. The purpose of this study was to map and identify the gene for isolated microspherophakia in two consanguineous Indian families. Using a whole-genome linkage scan in one family, we identiWed a likely locus for microspherophakia (MSP1) on chromosome 14q24.1–q32.12 between markers D14S588 and D14S1050 in a physical distance of 22.76 Mb. The maximum multi-point lod score was 2.91 between markers D14S1020 and D14S606. The MSP1 candidate region harbors 110 reference genes. DNA sequence analysis of one of the genes, LTBP2, detected a homozygous duplication (insertion) mutation, c.5446dupC, in the last exon (exon 36) in aVected family members. This homozygous mutation is predicted to elongate the LTBP2 protein by replacing the last 6 amino acids with 27 novel amino acids. Microspherophakia in the second family did not map to this locus, suggesting genetic heterogeneity. The present study suggests a role for LTBP2 in the structural stability of ciliary zonules, and growth and development of lens.

A homozygous mutation in LTBP2 causes isolated microspherophakia

Microspherophakia is an autosomal-recessive congenital disorder characterized by small spherical lens. It may be isolated or occur as part of a hereditary systemic disorder, such as Marfan syndrome, autosomal dominant and recessive forms of Weill-Marchesani syndrome, autosomal dominant glaucoma-lens ectopia-microspherophakia-stiffness-shortness syndrome, autosomal dominant microspherophakia with hernia, and microspherophakia-metaphyseal dysplasia. The purpose of this study was to map and identify the gene for isolated microspherophakia in two consanguineous Indian families. Using a whole-genome linkage scan in one family, we identified a likely locus for microspherophakia (MSP1) on chromosome 14q24.1-q32.12 between markers D14S588 and D14S1050 in a physical distance of 22.76 Mb. The maximum multi-point lod score was 2.91 between markers D14S1020 and D14S606. The MSP1 candidate region harbors 110 reference genes. DNA sequence analysis of one of the genes, LTBP2, detected a homozygous duplication (insertion) mutation, c.5446dupC, in the last exon (exon 36) in affected family members. This homozygous mutation is predicted to elongate the LTBP2 protein by replacing the last 6 amino acids with 27 novel amino acids. Microspherophakia in the second family did not map to this locus, suggesting genetic heterogeneity. The present study suggests a role for LTBP2 in the structural stability of ciliary zonules, and growth and development of lens.

Demographic history of India and mtDNA-sequence diversity

The demographic history of India was examined by comparing mtDNA sequences obtained from members of three culturally divergent Indian subpopulations (endogamous caste groups). While an inferred tree revealed some clustering according to caste affiliation, there was no clear separation into three genetically distinct groups along caste lines. Comparison of pairwise nucleotide difference distributions, however, did indicate a difference in growth patterns between two of the castes. The Brahmin population appears to have undergone either a rapid expansion or steady growth. The low-ranking Mukri caste, however, may have either maintained a roughly constant population size or undergone multiple bottlenecks during that period. Comparison of the Indian sequences to those obtained from other populations, using a tree, revealed that the Indian sequences, along with ah other non-African samples, form a starlike cluster. This cluster may represent a major expansion, possibly originating in southern Asia, taking place at some point after modern humans initially left Africa.