Faculty Bibliography

BACKGROUND:Leber congenital amaurosis (LCA) is a severe form of retinal dystrophy with marked underlying genetic heterogeneity. Until recently, allele-specific assays and Sanger sequencing of targeted segments were the only available approaches for attempted genetic diagnosis in this condition. A broader next-generation sequencing (NGS) strategy, such as whole exome sequencing, provides an improved molecular genetic diagnostic capacity for patients with these conditions. MATERIALS AND METHODS/METHODS:In a child with LCA, an allele-specific assay analyzing 135 known LCA-causing variations, followed by targeted segment sequencing of 61 regions in 14 causative genes was performed. Subsequently, exome sequencing was undertaken in the proband, unaffected consanguineous parents and two unaffected siblings. Bioinformatic analysis used two independent pipelines, BWA-GATK and SOAP, followed by Annovar and SnpEff to annotate the variants. RESULTS:No disease-causing variants were found using the allele-specific or targeted segment Sanger sequencing assays. Analysis of variants in the exome sequence data revealed a novel homozygous nonsense mutation (c.1081C > T, p.Arg361*) in TULP1, a gene with roles in photoreceptor function where mutations were previously shown to cause LCA and retinitis pigmentosa. The identified homozygous variant was the top candidate using both bioinformatic pipelines. CONCLUSIONS:This study highlights the value of the broad sequencing strategy of exome sequencing for disease gene identification in LCA, over other existing methods. NGS is particularly beneficial in LCA where there are a large number of causative disease genes, few distinguishing clinical features for precise candidate disease gene selection, and few mutation hotspots in any of the known disease genes.

This study aimed to identify the genetics underlying dominant forms of inherited retinal dystrophies using whole exome sequencing (WES) in six families extensively screened for known mutations or genes. Thirty-eight individuals were subjected to WES. Causative variants were searched among single nucleotide variants (SNVs) and insertion/deletion variants (indels) and whenever no potential candidate emerged, copy number variant (CNV) analysis was performed. Variants or regions harboring a candidate variant were prioritized and segregation of the variant with the disease was further assessed using Sanger sequencing in case of SNVs and indels, and quantitative PCR (qPCR) for CNVs. SNV and indel analysis led to the identification of a previously reported mutation in PRPH2. Two additional mutations linked to different forms of retinal dystrophies were identified in two families: a known frameshift deletion in RPGR, a gene responsible for X-linked retinitis pigmentosa and p.Ser163Arg in C1QTNF5 associated with Late-Onset Retinal Degeneration. A novel heterozygous deletion spanning the entire region of PRPF31 was also identified in the affected members of a fourth family, which was confirmed with qPCR. This study allowed the identification of the genetic cause of the retinal dystrophy and the establishment of a correct diagnosis in four families, including a large heterozygous deletion in PRPF31, typically considered one of the pitfalls of this method. Since all findings in this study are restricted to known genes, we propose that targeted sequencing using gene-panel is an optimal first approach for the genetic screening and that once known genetic causes are ruled out, WES might be used to uncover new genes involved in inherited retinal dystrophies.

BACKGROUND:Despite heritability estimates of 40-70 % for obesity, less than 2 % of its variation is explained by Body Mass Index (BMI) associated loci that have been identified so far. Epistasis, or gene-gene interactions are a plausible source to explain portions of the missing heritability of BMI. METHODS:Using genotypic data from 18,686 individuals across five study cohorts - ARIC, CARDIA, FHS, CHS, MESA - we filtered SNPs (Single Nucleotide Polymorphisms) using two parallel approaches. SNPs were filtered either on the strength of their main effects of association with BMI, or on the number of knowledge sources supporting a specific SNP-SNP interaction in the context of BMI. Filtered SNPs were specifically analyzed for interactions that are highly associated with BMI using QMDR (Quantitative Multifactor Dimensionality Reduction). QMDR is a nonparametric, genetic model-free method that detects non-linear interactions associated with a quantitative trait. RESULTS:We identified seven novel, epistatic models with a Bonferroni corrected p-value of association < 0.1. Prior experimental evidence helps explain the plausible biological interactions highlighted within our results and their relationship with obesity. We identified interactions between genes involved in mitochondrial dysfunction (POLG2), cholesterol metabolism (SOAT2), lipid metabolism (CYP11B2), cell adhesion (EZR), cell proliferation (MAP2K5), and insulin resistance (IGF1R). Moreover, we found an 8.8 % increase in the variance in BMI explained by these seven SNP-SNP interactions, beyond what is explained by the main effects of an index FTO SNP and the SNPs within these interactions. We also replicated one of these interactions and 58 proxy SNP-SNP models representing it in an independent dataset from the eMERGE study. CONCLUSION/CONCLUSIONS:This study highlights a novel approach for discovering gene-gene interactions by combining methods such as QMDR with traditional statistics.

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22,981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10(-8)). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

White blood cell (WBC) count is a common clinical measure used as a predictor of certain aspects of human health, including immunity and infection status. WBC count is also a complex trait that varies among individuals and ancestry groups. Differences in linkage disequilibrium structure and heterogeneity in allelic effects are expected to play a role in the associations observed between populations. Prior genome-wide association study (GWAS) meta-analyses have identified genomic loci associated with WBC and its subtypes, but much of the heritability of these phenotypes remains unexplained. Using GWAS summary statistics for over 50 000 individuals from three diverse populations (Japanese, African-American and European ancestry), a Bayesian model methodology was employed to account for heterogeneity between ancestry groups. This approach was used to perform a trans-ethnic meta-analysis of total WBC, neutrophil and monocyte counts. Ten previously known associations were replicated and six new loci were identified, including several regions harboring genes related to inflammation and immune cell function. Ninety-five percent credible interval regions were calculated to narrow the association signals and fine-map the putatively causal variants within loci. Finally, a conditional analysis was performed on the most significant SNPs identified by the trans-ethnic meta-analysis (MA), and nine secondary signals within loci previously associated with WBC or its subtypes were identified. This work illustrates the potential of trans-ethnic analysis and ascribes a critical role to multi-ethnic cohorts and consortia in exploring complex phenotypes with respect to variants that lie outside the European-biased GWAS pool.

BACKGROUND:Phosphoribosyl pyrophosphate synthetase (PRS) I deficiency is a rare medical condition caused by missense mutations in PRPS1 that lead to three different phenotypes: Arts Syndrome (MIM 301835), X-linked Charcot-Marie-Tooth (CMTX5, MIM 311070) or X-linked non-syndromic sensorineural deafness (DFN2, MIM 304500). All three are X-linked recessively inherited and males affected display variable degree of central and peripheral neuropathy. We applied whole exome sequencing to a three-generation family with optic atrophy followed by retinitis pigmentosa (RP) in all three cases, and ataxia, progressive peripheral neuropathy and hearing loss with variable presentation. METHODS:Whole exome sequencing was performed in two affecteds and one unaffected member of the family. Sanger sequencing was used to validate and segregate the 12 candidate mutations in the family and to confirm the absence of the novel variant in PRPS1 in 191 controls. The pathogenic role of the novel mutation in PRPS1 was assessed in silico and confirmed by enzymatic determination of PRS activity, mRNA expression and sequencing, and X-chromosome inactivation. RESULTS:A novel missense mutation was identified in PRPS1 in the affected females. Age of onset, presentation and severity of the phenotype are highly variable in the family: both the proband and her mother have neurological and ophthalmological symptoms, whereas the phenotype of the affected sister is milder and currently confined to the eye. Moreover, only the proband displayed a complete lack of expression of the wild type allele in leukocytes that seems to correlate with the degree of PRS deficiency and the severity of the phenotype. Interestingly, optic atrophy and RP are the only common manifestations to all three females and the only phenotype correlating with the degree of enzyme deficiency. CONCLUSIONS:These results are in line with recent evidence of the existence of intermediate phenotypes in PRS-I deficiency syndromes and demonstrate that females can exhibit a disease phenotype as severe and complex as their male counterparts.

BACKGROUND:Mendelian randomization studies have so far restricted attention to linear associations relating the genetic instrument to the exposure, and the exposure to the outcome. In some cases, however, observational data suggest a non-linear association between exposure and outcome. For example, alcohol consumption is consistently reported as having a U-shaped association with cardiovascular events. In principle, Mendelian randomization could address concerns that the apparent protective effect of light-to-moderate drinking might reflect 'sick-quitters' and confounding. METHODS:The Alcohol-ADH1B Consortium was established to study the causal effects of alcohol consumption on cardiovascular events and biomarkers, using the single nucleotide polymorphism rs1229984 in ADH1B as a genetic instrument. To assess non-linear causal effects in this study, we propose a novel method based on estimating local average treatment effects for discrete levels of the exposure range, then testing for a linear trend in those effects. Our method requires an assumption that the instrument has the same effect on exposure in all individuals. We conduct simulations examining the robustness of the method to violations of this assumption, and apply the method to the Alcohol-ADH1B Consortium data. RESULTS:Our method gave a conservative test for non-linearity under realistic violations of the key assumption. We found evidence for a non-linear causal effect of alcohol intake on several cardiovascular traits. CONCLUSIONS:We believe our method is useful for inferring departure from linearity when only a binary instrument is available. We estimated non-linear causal effects of alcohol intake which could not have been estimated through standard instrumental variable approaches.

BACKGROUND:Conflicting evidence exists on whether smoking acts as an effect modifier of the association between APOE genotype and risk of coronary heart disease (CHD). METHODS AND RESULTS/RESULTS:We searched PubMed and EMBASE to June 11, 2013 for published studies reporting APOE genotype, smoking status and CHD events and added unpublished data from population cohorts. We tested for presence of effect modification by smoking status in the relationship between APOE genotype and risk of CHD using likelihood ratio test. In total 13 studies (including unpublished data from eight cohorts) with 10,134 CHD events in 130,004 individuals of European descent were identified. The odds ratio (OR) for CHD risk from APOE genotype (ε4 carriers versus non-carriers) was 1.06 (95% confidence interval (CI): 1.01, 1.12) and for smoking (present vs. past/never smokers) was OR 2.05 (95%CI: 1.95, 2.14). When the association between APOE genotype and CHD was stratified by smoking status, compared to non-ε4 carriers, ε4 carriers had an OR of 1.11 (95%CI: 1.02, 1.21) in 28,789 present smokers and an OR of 1.04 (95%CI 0.98, 1.10) in 101,215 previous/never smokers, with no evidence of effect modification (P-value for heterogeneity = 0.19). Analysis of pack years in individual participant data of >60,000 with adjustment for cardiovascular traits also failed to identify evidence of effect modification. CONCLUSIONS:In the largest analysis to date, we identified no evidence for effect modification by smoking status in the association between APOE genotype and risk of CHD.