Mutations of conserved non-coding elements of PITX2 in patients with ocular dysgenesis and developmental glaucoma.

Mutations in FOXC1 and PITX2 constitute the most common causes of ocular anterior segment dysgenesis (ASD), and confer a high risk for secondary glaucoma. The genetic causes underlying ASD in approximately half of patients remain unknown, despite many of them being screened by whole exome sequencing. Here, we performed whole genome sequencing on DNA from two affected individuals from a family with dominantly inherited ASD and glaucoma to identify a 748-kb deletion in a gene desert that contains conserved putative PITX2 regulatory elements. We used CRISPR/Cas9 to delete the orthologous region in zebrafish in order to test the pathogenicity of this structural variant. Deletion in zebrafish reduced pitx2 expression during development and resulted in shallow anterior chambers. We screened additional patients for copy number variation of the putative regulatory elements and found an overlapping deletion in a second family and in a potentially-ancestrally-related index patient with ASD and glaucoma. These data suggest that mutations affecting conserved non-coding elements of PITX2 may constitute an important class of mutations in patients with ASD for whom the molecular cause of their disease have not yet been identified. Improved functional annotation of the human genome and transition to sequencing of patient genomes instead of exomes will be required before the magnitude of this class of mutations is fully understood.

Genetic analysis of a population heavy drinking phenotype identifies risk variants in whites.

Genetic association studies thus far have used detailed diagnoses of alcoholism to identify loci associated with risk. This proof-of-concept analysis examined whether population data of lifetime heaviest alcohol consumption may be used to identify genetic loci that modulate risk. We conducted a genetic association study in European Americans between variants in approximately 2100 genes and alcohol consumption as part of the Candidate gene Association Resource project. We defined cases as individuals with a history of drinking 5 or more drinks per day almost every day of the week and controls as current light drinkers (1-5 drinks per week). We cross-validated identified single nucleotide polymorphisms in a meta-analysis of 2 cohorts of unrelated individuals--Atherosclerosis Risk in Communities (ARIC) and Cardiovascular Health Study (CHS)--and in a separate cohort of related individuals--Framingham Heart Study (FHS). The most significant variant in the meta-analysis of ARIC and CHS was rs6933598 in methylenetetrahydrofolate dehydrogenase (P = 7.46 × 10(-05)) with a P value in FHS of 0.042. The top variants in FHS were rs12249562 in cubulin (P = 3.03 × 10(-05)) and rs9839267 near cholecystokinin (P = 3.05 × 10(-05)) with a P value of 0.019 for rs9839267 in CHS. We have here shown feasibility in evaluating lifetime incidence of heavy alcohol drinking from population-based studies for the purpose of conducting genetic association analyses.

Examination of rare missense variants in the CHRNA5-A3-B4 gene cluster to level of response to alcohol in the San Diego Sibling Pair study.

BACKGROUND: Common variants in the CHRNA5-A3-B4 gene cluster have been shown to be associated with nicotine dependence and alcohol use disorders (AUDs) and related traits, including the level of response (LR) to alcohol. Recently, rare variants (MAF < 0.05) in CHRNB4 have been reported to be associated with a decreased risk of developing nicotine dependence. However, the role of rare variants in the CHRNA5-A3-B4 gene cluster to the LR to alcohol has not yet been established. METHODS: To determine whether rare variants in the CHRNA5-A3-B4 gene cluster contribute to the LR to alcohol, the coding regions of these 3 genes were sequenced in 538 subjects from the San Diego Sibling Pair study. RESULTS: The analyses identified 16 rare missense variants, 9 of which were predicted to be damaging using in silico analysis tools. Carriers of these variants were compared to noncarriers using a family-based design for each gene and for the gene cluster as a whole. In these analyses, a CHRNA5 carrier status was significantly associated with the phenotype related to the feeling of intoxication experienced during the alcohol challenge (p = 0.039). CONCLUSIONS: These results indicate that rare genetic variation in the CHRNA5-A3-B4 gene cluster contributes modestly to the LR to alcohol in the San Diego Sibling Pair study and may protect against AUDs. However, replication studies are needed to confirm our findings.

Detection of identity by descent using next-generation whole genome sequencing data.

BACKGROUND: Identity by descent (IBD) has played a fundamental role in the discovery of genetic loci underlying human diseases. Both pedigree-based and population-based linkage analyses rely on estimating recent IBD, and evidence of ancient IBD can be used to detect population structure in genetic association studies. Various methods for detecting IBD, including those implemented in the soft- ware programs fastIBD and GERMLINE, have been developed in the past several years using population genotype data from microarray platforms. Now, next-generation DNA sequencing data is becoming increasingly available, enabling the comprehensive analysis of genomes, in- cluding identifying rare variants. These sequencing data may provide an opportunity to detect IBD with higher resolution than previously possible, potentially enabling the detection of disease causing loci that were previously undetectable with sparser genetic data. RESULTS: Here, we investigate how different levels of variant coverage in sequencing and microarray genotype data influences the resolution at which IBD can be detected. This includes microarray genotype data from the WTCCC study, denser genotype data from the HapMap Project, low coverage sequencing data from the 1000 Genomes Project, and deep coverage complete genome data from our own projects. With high power (78%), we can detect segments of length 0.4 cM or larger using fastIBD and GERMLINE in sequencing data. This compares to similar power to detect segments of length 1.0 cM or higher with microarray genotype data. We find that GERMLINE has slightly higher power than fastIBD for detecting IBD segments using sequencing data, but also has a much higher false positive rate. CONCLUSION: We further quantify the effect of variant density, conditional on genetic map length, on the power to resolve IBD segments. These investigations into IBD resolution may help guide the design of future next generation sequencing studies that utilize IBD, including family-based association studies, association studies in admixed populations, and homozygosity mapping studies.

Contribution of common PCSK1 genetic variants to obesity in 8,359 subjects from multi-ethnic American population.

Common PCSK1 variants (notably rs6232 and rs6235) have been shown to be associated with obesity in European, Asian and Mexican populations. To determine whether common PCSK1 variants contribute to obesity in American population, we conducted association analyses in 8,359 subjects using two multi-ethnic American studies: the Coronary Artery Risk Development in Young Adults (CARDIA) study and the Multi-Ethnic Study of Atherosclerosis (MESA). By evaluating the contribution of rs6232 and rs6235 in each ethnic group, we found that in European-American subjects from CARDIA, only rs6232 was associated with BMI (P = 0.006) and obesity (P = 0.018) but also increased the obesity incidence during the 20 years of follow-up (HR = 1.53 [1.07-2.19], P = 0.019). Alternatively, in African-American subjects from CARDIA, rs6235 was associated with BMI (P = 0.028) and obesity (P = 0.018). Further, by combining the two case-control ethnic groups from the CARDIA study in a meta-analysis, association between rs6235 and obesity risk remained significant (OR = 1.23 [1.05-1.45], P = 9.5×10(-3)). However, neither rs6232 nor rs6235 was associated with BMI or obesity in the MESA study. Interestingly, rs6232 was associated with BMI (P = 4.2×10(-3)) and obesity (P = 3.4×10(-3)) in the younger European-American group combining samples from the both studies [less than median age (53 years)], but not among the older age group (P = 0.756 and P = 0.935 for BMI and obesity, respectively). By combining all the case-control ethnic groups from CARDIA and MESA in a meta-analysis, we found no significant association for the both variants and obesity risk. Finally, by exploring the full PCSK1 locus, we observed that no variant remained significant after correction for multiple testing. These results indicate that common PCSK1 variants (notably rs6232 and rs6235) contribute modestly to obesity in multi-ethnic American population. Further, these results suggest that the association of rs6232 with obesity may be age-dependent in European-Americans. However, multiple replication studies in multi-ethnic American population are needed to confirm our findings.

Gene-centric analysis of serum cotinine levels in African and European American populations.

To date, most genetic association studies of tobacco use have been conducted in European American subjects using the phenotype of smoking quantity (cigarettes per day). However, smoking quantity is a very imprecise measure of exposure to tobacco smoke constituents. Analyses of alternate phenotypes and populations may improve our understanding of tobacco addiction genetics. Cotinine is the major metabolite of nicotine, and measuring serum cotinine levels in smokers provides a more objective measure of nicotine dose than smoking quantity. Previous genetic association studies of serum cotinine have focused on individual genes. We conducted a genetic association study of the biomarker in African American (N=365) and European American (N=315) subjects from the Coronary Artery Risk Development in Young Adults study using a chip containing densely-spaced tag SNPs in ∼2100 genes. We found that rs11187065, located in the non-coding region (intron 1) of insulin-degrading enzyme (IDE), was the most strongly associated SNP (p=8.91 × 10(-6)) in the African American cohort, whereas rs11763963, located on chromosome 7 outside of a gene transcript, was the most strongly associated SNP in European Americans (p=1.53 × 10(-6)). We then evaluated how the top variant association in each population performed in the other group. We found that the association of rs11187065 in IDE was also associated with the phenotype in European Americans (p=0.044). Our top SNP association in European Americans, rs11763963 was non-polymorphic in our African American sample. It has been previously shown that psychostimulant self-administration is reduced in animals with lower insulin because of interference with dopamine transmission in the brain reward centers. Our finding provides a platform for further investigation of this, or additional mechanisms, involving the relationship between insulin and self-administered nicotine dose.