Modifier effects between regulatory and protein-coding variation

Genome-wide associations have shown a lot of promise in dissecting the genetics of complex traits in humans with single variants, yet a large fraction of the genetic effects is still unaccounted for. Analyzing genetic interactions between variants (epistasis) is one of the potential ways forward. We investigated the abundance and functional impact of a specific type of epistasis, namely the interaction between regulatory and protein-coding variants. Using genotype and gene expression data from the 210 unrelated individuals of the original four HapMap populations, we have explored the combined effects of regulatory and protein-coding single nucleotide polymorphisms (SNPs). We predict that about 18% (1,502 out of 8,233 nsSNPs) of protein-coding variants are differentially expressed among individuals and demonstrate that regulatory variants can modify the functional effect of a coding variant in cis. Furthermore, we show that such interactions in cis can affect the expression of downstream targets of the gene containing the protein-coding SNP. In this way, a cis interaction between regulatory and protein-coding variants has a trans impact on gene expression. Given the abundance of both types of variants in human populations, we propose that joint consideration of regulatory and protein-coding variants may reveal additional genetic effects underlying complex traits and disease and may shed light on causes of differential penetrance of known disease variants.     

Rare variants in APP, PSEN1 and PSEN2 increase risk for AD in late-onset Alzheimer's disease families

Pathogenic mutations in APP, PSEN1, PSEN2, MAPT and GRN have previously been linked to familial early onset forms of dementia. Mutation screening in these genes has been performed in either very small series or in single families with late onset AD (LOAD). Similarly, studies in single families have reported mutations in MAPT and GRN associated with clinical AD but no systematic screen of a large dataset has been performed to determine how frequently this occurs. We report sequence data for 439 probands from late-onset AD families with a history of four or more affected individuals. Sixty sequenced individuals (13.7%) carried a novel or pathogenic mutation. Eight pathogenic variants, (one each in APP and MAPT, two in PSEN1 and four in GRN) three of which are novel, were found in 14 samples. Thirteen additional variants, present in 23 families, did not segregate with disease, but the frequency of these variants is higher in AD cases than controls, indicating that these variants may also modify risk for disease. The frequency of rare variants in these genes in this series is significantly higher than in the 1,000 genome project (p = 5.09×10⁻⁵; OR = 2.21; 95%CI = 1.49–3.28) or an unselected population of 12,481 samples (p = 6.82×10⁻⁵; OR = 2.19; 95%CI = 1.347–3.26). Rare coding variants in APP, PSEN1 and PSEN2, increase risk for or cause late onset AD. The presence of variants in these genes in LOAD and early-onset AD demonstrates that factors other than the mutation can impact the age at onset and penetrance of at least some variants associated with AD. MAPT and GRN mutations can be found in clinical series of AD most likely due to misdiagnosis. This study clearly demonstrates that rare variants in these genes could explain an important proportion of genetic heritability of AD, which is not detected by GWAS.     

Estimating the Efficacy and Efficiency of Cascade Genetic Screening

Screening for genetic variants that predispose individuals or their offspring to disease may be performed at the general population level or may instead be targeted at the relatives of previously identified carriers. The latter strategy has come to be known as "cascade genetic screening." Since the carrier risk of close relatives of known carriers is generally higher than the population risk, cascade screening is more efficient than population screening, in the sense that fewer individuals have to be genotyped per detected carrier. The efficacy of cascade screening, as measured by the overall proportion of carriers detected in a given population, is, however, lower than that of population-wide screening, and the respective inclusion rates vary according to the population frequency and mode of inheritance of the predisposing variants. For dominant mutations, we have developed equations that allow the inclusion rates of cascade screening to be calculated in an iterative fashion, depending upon screening depth and penetrance. For recessive mutations, we derived only equations for the screening of siblings and the children of patients. Owing to their mathematical complexity, it was necessary to study more extended screening strategies by simulation. Cascade screening turned out to result in low inclusion rates (<1%) when aimed at the identification of heterozygous carriers of rare recessive variants. Considerably higher rates are achievable, however, when screening is performed to detect covert homozygotes for frequent recessive mutations with reduced penetrance. This situation is exemplified by hereditary hemochromatosis, for which up to 40% of at-risk individuals may be identifiable through screening of first- to third-degree relatives of overt carriers (i.e., patients); the efficiency of this screening strategy was found to be approximately 50 times higher than that of population-wide screening. For dominant mutations, inclusion rates of cascade screening were estimated to be higher than for recessive variants. Thus, some 80% of all carriers of the factor V Leiden mutation would be detected if screening were to be targeted specifically at first- to third-degree relatives of patients with venous thrombosis. The relative cost efficiency of cascade as compared with population-wide screening (i.e., the overall savings in the extra managerial cost of the condition) is also likely to be higher for dominant than for recessive mutations. This notwithstanding, once screening has become cost-effective at the population level, it can be expected that cascade screening would only transiently represent an economically viable option.     

Reduced Penetrance and Variable Expression of SCN5A Mutations and the Importance of Co-inherited Genetic Variants: Case Report and Review of the Literature

Mutations in the SCN5A gene are responsible for multiple phenotypical presentations including Brugada syndrome, long QT syndrome, progressive familial heart block, sick sinus syndrome, dilated cardiomyopathy, lone atrial fibrillation and multiple overlap syndromes. These different phenotypic expressions of a mutation in a single gene can be explained by variable expression and reduced penetrance. One of the possible explanations of these phenomena is the co-inheritance of genetic variants. We describe a family where the individuals exhibit a compound heterozygosity in the SCN5A gene including a mutation (R1632H) and a new variant (M858L). Individuals with both the mutation and new variant present with a more severe phenotype including spontaneous atrial tachyarrhythmia at young age. We give an overview of the different phenotypes of "SCN5A disease" and discuss the importance of co-inherited genetic variants in the expression of SCN5A disease.     

The Characteristics of Heterozygous Protein Truncating Variants in the Human Genome

Sequencing projects have identified large numbers of rare stop-gain and frameshift variants in the human genome. As most of these are observed in the heterozygous state, they test a gene’s tolerance to haploinsufficiency and dominant loss of function. We analyzed the distribution of truncating variants across 16,260 autosomal protein coding genes in 11,546 individuals. We observed 39,893 truncating variants affecting 12,062 genes, which significantly differed from an expectation of 12,916 genes under a model of neutral de novo mutation (p<10⁻⁴). Extrapolating this to increasing numbers of sequenced individuals, we estimate that 10.8% of human genes do not tolerate heterozygous truncating variants. An additional 10 to 15% of truncated genes may be rescued by incomplete penetrance or compensatory mutations, or because the truncating variants are of limited functional impact. The study of protein truncating variants delineates the essential genome and, more generally, identifies rare heterozygous variants as an unexplored source of diversity of phenotypic traits and diseases.     

Analysis of two Arab families reveals additional support for a DFNB2 nonsyndromic phenotype of MYO7A

Variants in the head and tail domains of the MYO7A gene, encoding myosin VIIA, cause Usher syndrome type 1B (USH1B) and nonsyndromic deafness (DFNB2, DFNA11). In order to identify the genetic defect(s) underling profound deafness in two consanguineous Arab families living in UAE, we have sequenced a panel of 19 genes involved in Usher syndrome and nonsyndromic deafness in the index cases of the two families. This analysis revealed a novel homozygous insertion of AG (c.1952_1953insAG/p.C652fsX11) in exon 17 of the MYO7A gene in an Iraqi family, and a homozygous point mutation (c.5660C>T/p.P1887L) in exon 41 affecting the same gene in a large Palestinian family. Moreover, some individuals from the Palestinian family also harbored a novel heterozygous truncating variant (c.1267C>T/p.R423X) in the DFNB31 gene, which is involved in autosomal recessive nonsyndromic deafness type DFNB31 and Usher syndrome type II. Assuming an autosomal recessive mode of inheritance in the two inbred families, we conclude that the homozygous variants in the MYO7A gene are the disease-causing mutations in these families. Furthermore, given the absence of retinal disease in all affected patients examined, particularly a 28 year old patient, suggests that at least one family may segregate a DFNB2 presentation rather than USH1B. This finding further supports the premise that the MYO7A gene is responsible for two distinct diseases and gives evidence that the p.P1887L mutation in a homozygous state may be responsible for nonsyndromic hearing loss.     

Structural Modeling of a Novel CAPN5 Mutation that Causes Uveitis and Neovascular Retinal Detachment

CAPN5 mutations have been linked to autosomal dominant neovascular inflammatory vitreoretinopathy (ADNIV), a blinding autoimmune eye disease. Here, we link a new CAPN5 mutation to ADNIV and model the three-dimensional structure of the resulting mutant protein. In our study, a kindred with inflammatory vitreoretinopathy was evaluated by clinical eye examinations, DNA sequencing, and protein structural modeling to investigate the disease-causing mutation. Two daughters of an affected mother demonstrated symptoms of stage III ADNIV, with posterior uveitis, cystoid macular edema, intraocular fibrosis, retinal neovascularization, retinal degeneration, and cataract. The women also harbored a novel guanine to thymine (c.750G>T, p.Lys250Asn) missense mutation in exon 6 of CAPN5, a gene that encodes a calcium-activated cysteine protease, calpain-5. Modeling based on the structures of all known calpains revealed the mutation falls within a calcium-sensitive flexible gating loop that controls access to the catalytic groove. Three-dimensional modeling placed the new mutation in a region adjacent to two previously identified disease-causing mutations, all three of which likely disrupt hydrogen bonding within the gating loop, yielding a CAPN5 with altered enzymatic activity. This is the third case of a CAPN5 mutation leading to inherited uveitis and neovascular vitreoretinopathy, suggesting patients with ADNIV features should be tested for CAPN5 mutations. Structural modeling of novel variants can be used to support mechanistic consequences of the disease-causing variants.     

Mutational Analysis of Patients with Neurofibromatosis 2

Neurofibromatosis 2 (NF2) is a genetic disorder characterized by the development of multiple nervous-system tumors in young adulthood. The NF2 gene has recently been isolated and found to encode a new member of the protein 4.1 family of cytoskeletal associated proteins, which we have named merlin. To define the molecular basis of NF2 in affected individuals, we have used SSCP analysis to scan the exons of the NF2 gene from 33 unrelated patients with NF2. Twenty unique SSCP variants were seen in 21 patients; 10 of these individuals were known to be the only affected person in their kindred, while 7 had at least one other known affected relative. In all cases in which family members were available, the SSCP variant segregated with the disease; comparison of sporadic cases with their parents confirmed the de novo variants. DNA sequence analysis revealed that 19 of the 20 variants observed are predicted to lead to a truncated protein due to frameshift, creation of a stop codon, or interference with normal RNA splicing. A single patient carried a 3-bp deletion removing a phenylalanine residue. We conclude that the majority of NF2 patients carry an inactivating mutation of the NF2 gene and that neutral polymorphism in the gene is rare.     

Targeted next-generation sequencing identifies ABCA4 mutations in Chinese families with childhood-onset and adult-onset Stargardt disease

Background: Stargardt disease (STGD) is the most common form of juvenile macular dystrophy associated with progressive central vision loss, and is agenetically and clinically heterogeneous disease. Molecular diagnosis is of great significance in aiding the clinical diagnosis, helping to determine the phenotypic severity and visual prognosis. In the present study, we determined the clinical and genetic features of seven childhood-onset and three adult-onset Chinese STGD families. We performed capture next-generation sequencing (NGS) of the probands and searched for potentially disease-causing genetic variants in previously identified retinal or macular dystrophy genes.Methods: In all, ten unrelated Chinese families were enrolled. Panel-based NGS was performed to identify potentially disease-causing genetic variants in previously identified retinal or macular dystrophy genes, including the five known STGD genes (ABCA4, PROM1, PRPH2, VMD2, and ELOVL4). Variant analysis, Sanger validation, and segregation tests were utilized to validate the disease-causing mutations in these families.Results: Using systematic data analysis with an established bioinformatics pipeline and segregation analysis, 17 pathogenic mutations in ABCA4 were identified in the 10 STGD families. Four of these mutations were novel: c.371delG, c.681T > G, c.5509C > T, and EX37del. Childhood-onset STGD was associated with severe visual loss, generalized retinal dysfunction and was due to more severe variants in ABCA4 than those found in adult-onset disease.Conclusions: We expand the existing spectrum of STGD and reveal the genotype–phenotype relationships of the ABCA4 mutations in Chinese patients. Childhood-onset STGD lies at the severe end of the spectrum of ABCA4-associated retinal phenotypes.     

SNP Linkage Analysis and Whole Exome Sequencing Identify a Novel POU4F3 Mutation in Autosomal Dominant Late-Onset Nonsyndromic Hearing Loss (DFNA15)

Autosomal dominant non-syndromic hearing loss (AD-NSHL) is one of the most common genetic diseases in human and is well-known for the considerable genetic heterogeneity. In this study, we utilized whole exome sequencing (WES) and linkage analysis for direct genetic diagnosis in AD-NSHL. The Korean family had typical AD-NSHL running over 6 generations. Linkage analysis was performed by using genome-wide single nucleotide polymorphism (SNP) chip and pinpointed a genomic region on 5q31 with a significant linkage signal. Sequential filtering of variants obtained from WES, application of the linkage region, bioinformatic analyses, and Sanger sequencing validation identified a novel missense mutation Arg326Lys (c.977G>A) in the POU homeodomain of the POU4F3 gene as the candidate disease-causing mutation in the family. POU4F3 is a known disease gene causing AD-HSLH (DFNA15) described in 5 unrelated families until now each with a unique mutation. Arg326Lys was the first missense mutation affecting the 3^rd alpha helix of the POU homeodomain harboring a bipartite nuclear localization signal sequence. The phenotype findings in our family further supported previously noted intrafamilial and interfamilial variability of DFNA15. This study demonstrated that WES in combination with linkage analysis utilizing bi-allelic SNP markers successfully identified the disease locus and causative mutation in AD-NSHL.     

Exome-Based Mapping and Variant Prioritization for Inherited Mendelian Disorders

Exome sequencing in families affected by rare genetic disorders has the potential to rapidly identify new disease genes (genes in which mutations cause disease), but the identification of a single causal mutation among thousands of variants remains a significant challenge. We developed a scoring algorithm to prioritize potential causal variants within a family according to segregation with the phenotype, population frequency, predicted effect, and gene expression in the tissue(s) of interest. To narrow the search space in families with multiple affected individuals, we also developed two complementary approaches to exome-based mapping of autosomal-dominant disorders. One approach identifies segments of maximum identity by descent among affected individuals; the other nominates regions on the basis of shared rare variants and the absence of homozygous differences between affected individuals. We showcase our methods by using exome sequence data from families affected by autosomal-dominant retinitis pigmentosa (adRP), a rare disorder characterized by night blindness and progressive vision loss. We performed exome capture and sequencing on 91 samples representing 24 families affected by probable adRP but lacking common disease-causing mutations. Eight of 24 families (33%) were revealed to harbor high-scoring, most likely pathogenic (by clinical assessment) mutations affecting known RP genes. Analysis of the remaining 17 families identified candidate variants in a number of interesting genes, some of which have withstood further segregation testing in extended pedigrees. To empower the search for Mendelian-disease genes in family-based sequencing studies, we implemented them in a cross-platform-compatible software package, MendelScan, which is freely available to the research community.     

Novel A14841G mutation is associated with high penetrance of LHON/C4171A family

We report the clinical and genetic characterization of a Chinese LHON family carrying an ND1/C4171A mutation. This family has high penetrance of visual impairment and extremely low frequency of vision recovery, which is in marked contrast to previously reported results for Korean LHON families with the same mutation. Sequence analysis of the complete mtDNA in the partially defined East Asian haplogroup N9a1 revealed the presence of 29 other variants. A novel heteroplasmic A14841G mutation, one of the variants with a serine substituted for a highly conserved asparagine at amino acid 32 of Cytochrome b (Cytb), may play a synergistic role with the C4171A mutation, leading to significantly different clinical manifestations of LHON among these families. The study further confirmed that C4171A was a rare primary LHON mutation, and the mtDNA background could also contribute to the clinical manifestation of the LHON/C4171A mutation.     

TMEM43 mutations associated with arrhythmogenic right ventricular cardiomyopathy in non-Newfoundland populations

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a myocardial disease characterized by fibro-fatty replacement of right ventricular free wall myocardium and life-threatening ventricular arrhythmias. A missense mutation, c.1073C>T (p.S358L) in the transmembrane protein 43 (TMEM43) gene, has been genetically identified to cause ARVC type 5 in a founder population from Newfoundland. It is unclear whether this mutation occurs in other populations outside of this founder population or if other variants of TMEM43 are associated with ARVC disease. We sought to identify non-Newfoundland individuals with TMEM43 variants among patient samples sent for genetic assessment for possible ARVC. Of 195 unrelated individuals with suspected ARVC, mutation of desmosomal proteins was seen in 28 and the p.S358L TMEM43 mutation in six. We identified a de novo p.S358L mutation in a non-Newfoundland patient and five separate rare TMEM43 (four novel) sequence variants in non-Newfoundland patients, each occurring in an evolutionarily conserved amino acid. TMEM43 mutations occur outside of the founder population of the island of Newfoundland where it was originally described. TMEM43 sequencing should be incorporated into clinical genetic testing for ARVC patients.     

Focus: Rare Disease: Pulmonary Alveolar Microlithiasis – A Review

Pulmonary Alveolar Microlithiasis (PAM) is a rare genetic disorder causing widespread deposition of calcium-phosphate crystals in the alveolar space. A hallmark of the disease is the discrepancy between perceived symptoms upon diagnosis compared with the extensive, sandstorm-like appearance of the microliths on chest X-ray or HRCT. Caused by a defective sodium-dependent phosphate transport protein due to loss-of-function variants of the SLC34A2 gene, PAM is an autosomal recessive transmitted disorder, and as such has a high correlation to consanguinity. The most common variants of the SLC34A2 gene are single nucleotide biallelic changes, but larger deletions are described. Initial suspicion of PAM on radiological examination should be followed by genetic testing to verify the diagnosis and identify the disease-causing variant. When not available, the diagnosis can be made by means of invasive techniques, such as transbronchial forceps or cryobiopsy, or a surgical lung biopsy. In families with a history of PAM, genetic counseling should be offered, as well as preimplantation/prenatal testing if necessary. As of writing this review, no definitive treatment exists, and PAM may in some cases progress to severe pulmonary disease with respiratory failure and potential death. Patients with PAM should be offered preventative and symptomatic treatments such as vaccinations and oxygen therapy when needed. In some cases, lung transplantation may be required.     

Missense mutations in disease genes: a Bayesian approach to evaluate causality.

The problem of interpreting missense mutations of disease-causing genes is an increasingly important one. Because these point mutations result in alteration of only a single amino acid of the protein product, it is often unclear whether this change alone is sufficient to cause disease. We propose a Bayesian approach that utilizes genetic information on affected relatives in families ascertained through known missense-mutation carriers. This method is useful in evaluating known disease genes for common disease phenotypes, such as breast cancer or colorectal cancer. The posterior probability that a missense mutation is disease causing is conditioned on the relationship of the relatives to the proband, the population frequency of the mutation, and the phenocopy rate of the disease. The approach is demonstrated in two cancer data sets: BRCA1 R841W and APC I1307K. In both examples, this method helps establish that these mutations are likely to be disease causing, with Bayes factors in favor of causality of 5.09 and 66.97, respectively, and posterior probabilities of .836 and .985. We also develop a simple approximation for rare alleles and consider the case of unknown penetrance and allele frequency.     

A Novel Missense SNRNP200 Mutation Associated with Autosomal Dominant Retinitis Pigmentosa in a Chinese Family

The SNRNP200 gene encodes hBrr2, a helicase essential for pre-mRNA splicing. Six mutations in SNRNP200 have recently been discovered to be associated with autosomal dominant retinitis pigmentosa (adRP). In this work, we analyzed a Chinese family with adRP and identified a novel missense mutation in SNRNP200. To identify the genetic defect in this family, exome of the proband was captured and sequencing analysis was performed to exclude known genetic defects and find possible pathogenic mutations. Subsequently, candidate mutations were validated in affected family members using Sanger sequencing. A novel missense mutation, c.2653C>G transition (p.Q885E), in exon 20 of SNRNP200 was identified. The mutation co-segregated with the disease phenotype over four generations and was absent in 100 normal unaffected individuals. This mutation occurs at highly conserved position in hBrr2 and is predicted to have a functional impact, suggesting that hBrr2-dependent small nuclear riboproteins (snRNPs) unwinding and spliceosome activation is important in the pathogenesis of some variants of RP.     

Evolution and disease converge in the mitochondrion

Mitochondrial DNA (mtDNA) mutations are long known to cause diseases but also underlie tremendous population divergence in humans. It was assumed that the two types of mutations differ in one major trait: functionality. However, evidence from disease association studies, cell culture and animal models support the functionality of common mtDNA genetic variants, leading to the hypothesis that disease-causing mutations and mtDNA genetic variants share considerable common features. Here we provide evidence showing that the two types of mutations obey the rules of evolution, including random genetic drift and natural selection. This similarity does not only converge at the principle level; rather, disease-causing mutations could recapitulate the ancestral DNA sequence state. Thus, the very same mutations could either mark ancient evolutionary changes or cause disease.     

Exome Analysis of Two Limb-Girdle Muscular Dystrophy Families: Mutations Identified and Challenges Encountered

The molecular diagnosis of muscle disorders is challenging: genetic heterogeneity (>100 causal genes for skeletal and cardiac muscle disease) precludes exhaustive clinical testing, prioritizing sequencing of specific genes is difficult due to the similarity of clinical presentation, and the number of variants returned through exome sequencing can make the identification of the disease-causing variant difficult. We have filtered variants found through exome sequencing by prioritizing variants in genes known to be involved in muscle disease while examining the quality and depth of coverage of those genes. We ascertained two families with autosomal dominant limb-girdle muscular dystrophy of unknown etiology. To identify the causal mutations in these families, we performed exome sequencing on five affected individuals using the Agilent SureSelect Human All Exon 50 Mb kit and the Illumina HiSeq 2000 (2×100 bp). We identified causative mutations in desmin (IVS3+3A>G) and filamin C (p.W2710X), and augmented the phenotype data for individuals with muscular dystrophy due to these mutations. We also discuss challenges encountered due to depth of coverage variability at specific sites and the annotation of a functionally proven splice site variant as an intronic variant.     

Mitochondrial ND6 T14502C variant may modulate the phenotypic expression of LHON-associated G11778A mutation in four Chinese families

We report here the clinical, genetic, and molecular evaluations of four Han Chinese families with Leber’s hereditary optic neuropathy. Thirty-one (20 males/11 females) of 83 matrilineal relatives in these families exhibited the variable severity and age-at-onset in visual impairment. The average age-of-onset of vision loss was 22 years old. Strikingly, these penetrances of visual impairment in these Chinese families were higher than those in other 11 Chinese pedigrees carrying the only ND4 G11778A mutation. Molecular analysis identified the known G11778A mutation and distinct sets of variants belonging to the Asian haplogroups M10a and M7c2. Of these, the T14502C mutation caused the substitution of a highly conserved isoleucine for valine at position 58 in ND6. This mutation has been associated with LHON in other Chinese families with very low penetrance of LHON. Thus, the deficient activities of complex I, caused by G11778A mutation, would be worsened by the T14502C mutation in these four Chinese families. As a result, mitochondrial dysfunctions would lead to the high penetrance and expressivity of visual loss in these Chinese families carrying both G11778A and T14502C mutations than other 11 Chinese families carrying only G11778A mutation. These data suggested that the T14502C variant may modulate the phenotypic manifestation of the G11778A mutation in these Chinese pedigrees.