Feasibility of High-Throughput Genome-Wide Genotyping using DNA from Stored Buccal Cell Samples

It is unclear if buccal cell samples contain sufficient human DNA with adequately sized fragments for high throughput genetic bioassays. Yet buccal cell sample collection is an attractive alternative to gathering blood samples for genetic epidemiologists engaged in large-scale genetic biomarker studies. We assessed the genotyping efficiency (GE) and genotyping concordance (GC) of buccal cell DNA samples compared to corresponding blood DNA samples, from 32 Nurses’ Health Study (NHS) participants using the Illumina Infinium 660W-Quad platform. We also assessed how GE and GC accuracy varied as a function of DNA concentration using serial dilutions of buccal DNA samples. Finally we determined the nature and genomic distribution of discordant genotypes in buccal DNA samples. The mean GE of undiluted buccal cell DNA samples was high (99.32%), as was the GC between the paired buccal and blood samples (99.29%). GC between the dilutions versus the undiluted buccal DNA was also very high (greater than 97%), though both GE and GC notably declined at DNA concentrations less than 5 ng/μl. Most (greater than 95%) genotype determinations in buccal cell samples were of the “missing call” variety (as opposed to the “alternative genotype call” variety) across the spectrum of buccal DNA concentrations studied. Finally, for buccal DNA concentration above 1.7 ng/ul, discordant genotyping calls did not cluster in any particular chromosome. Buccal cell-derived DNA represents a viable alternative to blood DNA for genotyping on a high-density platform

Feasibility of High-Throughput Genome-Wide Genotyping using DNA from Stored Buccal Cell Samples

It is unclear if buccal cell samples contain sufficient human DNA with adequately sized fragments for high throughput genetic bioassays. Yet buccal cell sample collection is an attractive alternative to gathering blood samples for genetic epidemiologists engaged in large-scale genetic biomarker studies. We assessed the genotyping efficiency (GE) and genotyping concordance (GC) of buccal cell DNA samples compared to corresponding blood DNA samples, from 32 Nurses’ Health Study (NHS) participants using the Illumina Infinium 660W-Quad platform. We also assessed how GE and GC accuracy varied as a function of DNA concentration using serial dilutions of buccal DNA samples. Finally we determined the nature and genomic distribution of discordant genotypes in buccal DNA samples. The mean GE of undiluted buccal cell DNA samples was high (99.32%), as was the GC between the paired buccal and blood samples (99.29%). GC between the dilutions versus the undiluted buccal DNA was also very high (>97%), though both GE and GC notably declined at DNA concentrations less than 5 ng/μl. Most (>95%) genotype determinations in buccal cell samples were of the “missing call” variety (as opposed to the “alternative genotype call” variety) across the spectrum of buccal DNA concentrations studied. Finally, for buccal DNA concentration above 1.7 ng/ul, discordant genotyping calls did not cluster in any particular chromosome. Buccal cell-derived DNA represents a viable alternative to blood DNA for genotyping on a high-density platform

Screening for UGT1A1 Genotype in Study A5257 Would Have Markedly Reduced Premature Discontinuation of Atazanavir for Hyperbilirubinemia

Background. Some patients are not prescribed atazanavir because of concern about possible jaundice. Atazanavir-associated hyperbilirubinemia correlates with UGT1A1 rs887829 genotype. We examined bilirubin-related discontinuation of atazanavir in participants from AIDS Clinical Trials Group Study A5257. Methods. Discriminatory properties of UGT1A1 T/T genotype for predicting bilirubin-related atazanavir discontinuation through 96 weeks after antiretroviral initiation were estimated. Results. Genetic analyses involved 1450 participants, including 481 who initiated randomized atazanavir/ritonavir. Positive predictive values of rs887829 T/T for bilirubin-related discontinuation of atazanavir (with 95% confidence intervals [CIs]) were 20% (CI, 9%–36%) in Black, 60% (CI, 32%–84%) in White, and 29% (CI, 8%–58%) in Hispanic participants; negative predictive values were 97% (CI, 93%–99%), 95% (CI, 90%–98%), and 97% (CI, 90%–100%), respectively. Conclusions. Bilirubin-related discontinuation of atazanavir was rare in participants not homozygous for rs887829 T/T, regardless of race or ethnicity. We hypothesize that the higher rate of discontinuation among White participants homozygous for rs887829 T/T may reflect differences in physical manifestations of jaundice by race and ethnicity. Selective avoidance of atazanavir initiation among individuals with T/T genotypes would markedly reduce the likelihood of bilirubin-related discontinuation of atazanavir while allowing atazanavir to be prescribed to the majority of individuals. This genetic association will also affect atazanavir/cobicistat

Mitochondrial DNA Polymorphism A4917G Is Independently Associated with Age-Related Macular Degeneration

The objective of this study was to determine if MTND2*LHON4917G (4917G), a specific non-synonymous polymorphism in the mitochondrial genome previously associated with neurodegenerative phenotypes, is associated with increased risk for age-related macular degeneration (AMD). A preliminary study of 393 individuals (293 cases and 100 controls) ascertained at Vanderbilt revealed an increased occurrence of 4917G in cases compared to controls (15.4% vs.9.0%, p = 0.11). Since there was a significant age difference between cases and controls in this initial analysis, we extended the study by selecting Caucasian pairs matched at the exact age at examination. From the 1547 individuals in the Vanderbilt/Duke AMD population association study (including 157 in the preliminary study), we were able to match 560 (280 cases and 280 unaffected) on exact age at examination. This study population was genotyped for 4917G plus specific AMD-associated nuclear genome polymorphisms in CFH, LOC387715 and ApoE. Following adjustment for the listed nuclear genome polymorphisms, 4917G independently predicts the presence of AMD (OR = 2.16, 95%CI 1.20–3.91, p = 0.01). In conclusion, a specific mitochondrial polymorphism previously implicated in other neurodegenerative phenotypes (4917G) appears to convey risk for AMD independent of recently discovered nuclear DNA polymorphisms

Mapping of the Disease Locus and Identification of ADAMTS10 As a Candidate Gene in a Canine Model of Primary Open Angle Glaucoma

Primary open angle glaucoma (POAG) is a leading cause of blindness worldwide, with elevated intraocular pressure as an important risk factor. Increased resistance to outflow of aqueous humor through the trabecular meshwork causes elevated intraocular pressure, but the specific mechanisms are unknown. In this study, we used genome-wide SNP arrays to map the disease gene in a colony of Beagle dogs with inherited POAG to within a single 4 Mb locus on canine chromosome 20. The Beagle POAG locus is syntenic to a previously mapped human quantitative trait locus for intraocular pressure on human chromosome 19. Sequence capture and next-generation sequencing of the entire canine POAG locus revealed a total of 2,692 SNPs segregating with disease. Of the disease-segregating SNPs, 54 were within exons, 8 of which result in amino acid substitutions. The strongest candidate variant causes a glycine to arginine substitution in a highly conserved region of the metalloproteinase ADAMTS10. Western blotting revealed ADAMTS10 protein is preferentially expressed in the trabecular meshwork, supporting an effect of the variant specific to aqueous humor outflow. The Gly661Arg variant in ADAMTS10 found in the POAG Beagles suggests that altered processing of extracellular matrix and/or defects in microfibril structure or function may be involved in raising intraocular pressure, offering specific biochemical targets for future research and treatment strategies

Identification of a Rare Coding Variant in Complement 3 Associated with Age-related Macular Degeneration

Macular degeneration is a common cause of blindness in the elderly. To identify rare coding variants associated with a large increase in risk of age-related macular degeneration (AMD), we sequenced 2,335 cases and 789 controls in 10 candidate loci (57 genes). To increase power, we augmented our control set with ancestry-matched exome sequenced controls. An analysis of coding variation in 2,268 AMD cases and 2,268 ancestry matched controls revealed two large-effect rare variants; previously described R1210C in the CFH gene (fcase = 0.51%, fcontrol = 0.02%, OR = 23.11), and newly identified K155Q in the C3 gene (fcase = 1.06%, fcontrol = 0.39%, OR = 2.68). The variants suggest decreased inhibition of C3 by Factor H, resulting in increased activation of the alternative complement pathway, as a key component of disease biology

Using genetic variation and environmental risk factor data to identify individuals at high risk for age-related macular degeneration

A major goal of personalized medicine is to pre-symptomatically identify individuals at high risk for disease using knowledge of each individual's particular genetic profile and constellation of environmental risk factors. With the identification of several well-replicated risk factors for age-related macular degeneration (AMD), the leading cause of legal blindness in older adults, this previously unreachable goal is beginning to seem less elusive. However, recently developed algorithms have either been much less accurate than expected, given the strong effects of the identified risk factors, or have not been applied to independent datasets, leaving unknown how well they would perform in the population at large. We sought to increase accuracy by using novel modeling strategies, including multifactor dimensionality reduction (MDR) and grammatical evolution of neural networks (GENN), in addition to the traditional logistic regression approach. Furthermore, we rigorously designed and tested our models in three distinct datasets: a Vanderbilt-Miami (VM) clinic-based case-control dataset, a VM family dataset, and the population-based Age-related Maculopathy Ancillary (ARMA) Study cohort. Using a consensus approach to combine the results from logistic regression and GENN models, our algorithm was successful in differentiating between high- and low-risk groups (sensitivity 77.0%, specificity 74.1%). In the ARMA cohort, the positive and negative predictive values were 63.3% and 70.7%, respectively. We expect that future efforts to refine this algorithm by increasing the sample size available for model building, including novel susceptibility factors as they are discovered, and by calibrating the model for diverse populations will improve accuracy

Genome-wide and Ordered-Subset linkage analyses provide support for autism loci on 17q and 19p with evidence of phenotypic and interlocus genetic correlates

BACKGROUND: Autism is a neurobehavioral spectrum of phenotypes characterized by deficits in the development of language and social relationships and patterns of repetitive, rigid and compulsive behaviors. Twin and family studies point to a significant genetic etiology, and several groups have performed genomic linkage screens to identify susceptibility loci. METHODS: We performed a genome-wide linkage screen in 158 combined Tufts, Vanderbilt and AGRE (Autism Genetics Research Exchange) multiplex autism families using parametric and nonparametric methods with a categorical autism diagnosis to identify loci of main effect. Hypothesizing interdependence of genetic risk factors prompted us to perform exploratory studies applying the Ordered-Subset Analysis (OSA) approach using LOD scores as the trait covariate for ranking families. We employed OSA to test for interlocus correlations between loci with LOD scores ≥1.5, and empirically determined significance of linkage in optimal OSA subsets using permutation testing. Exploring phenotypic correlates as the basis for linkage increases involved comparison of mean scores for quantitative trait-based subsets of autism between optimal subsets and the remaining families. RESULTS: A genome-wide screen for autism loci identified the best evidence for linkage to 17q11.2 and 19p13, with maximum multipoint heterogeneity LOD scores of 2.9 and 2.6, respectively. Suggestive linkage (LOD scores ≥1.5) at other loci included 3p, 6q, 7q, 12p, and 16p. OSA revealed positive correlations of linkage between the 19p locus and 17q, between 19p and 6q, and between 7q and 5p. While potential phenotypic correlates for these findings were not identified for the chromosome 7/5 combination, differences indicating more rapid achievement of "developmental milestones" was apparent in the chromosome 19 OSA-defined subsets for 17q and 6q. OSA was used to test the hypothesis that 19p linkage involved more rapid achievement of these milestones and it revealed significantly increased LOD* scores at 19p13. CONCLUSIONS: Our results further support 19p13 as harboring an autism susceptibility locus, confirm other linkage findings at 17q11.2, and demonstrate the need to analyze more discreet trait-based subsets of complex phenotypes to improve ability to detect genetic effects

Common Variants at 9p21 and 8q22 Are Associated with Increased Susceptibility to Optic Nerve Degeneration in Glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10−18), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10−11). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10−12) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10−10). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma

CDKN2B-AS1 Genotype–Glaucoma Feature Correlations in Primary Open-Angle Glaucoma Patients From the United States

To assess the association between single nucleotide polymorphisms (SNPs) of the gene region containing cyclin dependent kinase inhibitor 2B antisense noncoding RNA (CDKN2B-AS1) and glaucoma features among primary open-angle glaucoma (POAG) patients