Digenic Inheritance and Gene-Environment Interaction in a Patient With Hypertriglyceridemia and Acute Pancreatitis

The etiology of hypertriglyceridemia (HTG) and acute pancreatitis (AP) is complex. Herein, we dissected the underlying etiology in a patient with HTG and AP. The patient had a 20-year history of heavy alcohol consumption and an 8-year history of mild HTG. He was hospitalized for alcohol-triggered AP, with a plasma triglyceride (TG) level up to 21.4 mmol/L. A temporary rise in post-heparin LPL concentration (1.5–2.5 times of controls) was noted during the early days of AP whilst LPL activity was consistently low (50∼70% of controls). His TG level rapidly decreased to normal in response to treatment, and remained normal to borderline high during a ∼3-year follow-up period during which he had abstained completely from alcohol. Sequencing of the five primary HTG genes (i.e., LPL, APOC2, APOA5, GPIHBP1 and LMF1) identified two heterozygous variants. One was the common APOA5 c.553G > T (p.Gly185Cys) variant, which has been previously associated with altered TG levels as well as HTG-induced acute pancreatitis (HTG-AP). The other was a rare variant in the LPL gene, c.756T > G (p.Ile252Met), which was predicted to be likely pathogenic and found experimentally to cause a 40% loss of LPL activity without affecting either protein synthesis or secretion. We provide evidence that both a gene-gene interaction (between the common APOA5 variant and the rare LPL variant) and a gene-environment interaction (between alcohol and digenic inheritance) might have contributed to the development of mild HTG and alcohol-triggered AP in the patient, thereby improving our understanding of the complex etiology of HTG and HTG-AP.

Prioritization of schizophrenia risk genes from GWAS results by integrating multi-omics data

Schizophrenia (SCZ) is a polygenic disease with a heritability approaching 80%. Over 100 SCZ-related loci have so far been identified by genome-wide association studies (GWAS). However, the risk genes associated with these loci often remain unknown. We present a new risk gene predictor, rGAT-omics, that integrates multi-omics data under a Bayesian framework by combining the Hotelling and Box-Cox transformations. The Bayesian framework was constructed using gene ontology, tissue-specific protein-protein networks, and multi-omics data including differentially expressed genes in SCZ and controls, distance from genes to the index single-nucleotide polymorphisms (SNPs), and de novo mutations. The application of rGAT-omics to the 108 loci identified by a recent GWAS study of SCZ predicted 103 high-risk genes (HRGs) that explain a high proportion of SCZ heritability (Enrichment = 43.44 and [Formula: see text]). HRGs were shown to be significantly ([Formula: see text]) enriched in genes associated with neurological activities, and more likely to be expressed in brain tissues and SCZ-associated cell types than background genes. The predicted HRGs included 16 novel genes not present in any existing databases of SCZ-associated genes or previously predicted to be SCZ risk genes by any other method. More importantly, 13 of these 16 genes were not the nearest to the index SNP markers, and them would have been difficult to identify as risk genes by conventional approaches while ten out of the 16 genes are associated with neurological functions that make them prime candidates for pathological involvement in SCZ. Therefore, rGAT-omics has revealed novel insights into the molecular mechanisms underlying SCZ and could provide potential clues to future therapies.

The reversion variant (p.Arg90Leu) at the evolutionarily adaptive p.Arg90 site in CELA3B predisposes to chronic pancreatitis

A gain-of-function missense variant in the CELA3B gene, p.Arg90Cys (c.268C>T), has recently been reported to cause pancreatitis in an extended pedigree. Herein, we sequenced the CELA3B gene in 644 genetically unexplained French chronic pancreatitis (CP) patients (all unrelated) and 566 controls. No obvious loss-of-function variants were identified. None of the six low-frequency or common missense variants detected showed significant association with CP. Nor did the aggregate rare/very rare missense variants (n = 14) show any significant association with CP. However, p.Arg90Leu (c.269G>T), which was found in four patients but no controls, and affects the same amino acid as p.Arg90Cys, serves to revert p.Arg90 to the human elastase ancestral allele. As p.Arg90Leu has previously been shown to exert a similar functional effect to that of p.Arg90Cys, our findings not only confirm the involvement of CELA3B in the etiology of CP but also pinpoint a new evolutionarily adaptive site in the human genome.

A platform for curated products from novel open reading frames prompts reinterpretation of disease variants

Recent evidence from proteomics and deep massively parallel sequencing studies have revealed that eukaryotic genomes contain substantial numbers of as-yet-uncharacterized open reading frames (ORFs). We define these uncharacterized ORFs as novel ORFs (nORFs). nORFs in humans are mostly under 100 codons and are found in diverse regions of the genome, including in long noncoding RNAs, pseudogenes, 3' UTRs, 5' UTRs, and alternative reading frames of canonical protein coding exons. There is therefore a pressing need to evaluate the potential functional importance of these unannotated transcripts and proteins in biological pathways and human disease on a larger scale, rather than one at a time. In this study, we outline the creation of a valuable nORFs data set with experimental evidence of translation for the community, use measures of heritability and selection that reveal signals for functional importance, and show the potential implications for functional interpretation of genetic variants in nORFs. Our results indicate that some variants that were previously classified as being benign or of uncertain significance may have to be reinterpreted.

Heritable pattern of oxidized DNA base repair coincides with pre-targeting of repair complexes to open chromatin

Human genome stability requires efficient repair of oxidized bases, which is initiated via damage recognition and excision by NEIL1 and other base excision repair (BER) pathway DNA glycosylases (DGs). However, the biological mechanisms underlying detection of damaged bases among the million-fold excess of undamaged bases remain enigmatic. Indeed, mutation rates vary greatly within individual genomes, and lesion recognition by purified DGs in the chromatin context is inefficient. Employing super-resolution microscopy and co-immunoprecipitation assays, we find that acetylated NEIL1 (AcNEIL1), but not its non-acetylated form, is predominantly localized in the nucleus in association with epigenetic marks of uncondensed chromatin. Furthermore, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) revealed non-random AcNEIL1 binding near transcription start sites of weakly transcribed genes and along highly transcribed chromatin domains. Bioinformatic analyses revealed a striking correspondence between AcNEIL1 occupancy along the genome and mutation rates, with AcNEIL1-occupied sites exhibiting fewer mutations compared to AcNEIL1-free domains, both in cancer genomes and in population variation. Intriguingly, from the evolutionarily conserved unstructured domain that targets NEIL1 to open chromatin, its damage surveillance of highly oxidation-susceptible sites to preserve essential gene function and to limit instability and cancer likely originated ∼500 million years ago during the buildup of free atmospheric oxygen.

Author Correction: Neuroprotectants attenuate hypobaric hypoxia-induced brain injuries in cynomolgus monkeys

Following the publication of our paper (Zhang et al., 2020), it has come to our attention that we erroneously listed two funding sources unrelated to this study in the “ACKNOWLEDGEMENTS” section. Hereby, we wish to update the “ACKNOWLEDGEMENTS” section as a correction:

EXT1 and EXT2 Variants in 22 Chinese Families With Multiple Osteochondromas: Seven New Variants and Potentiation of Preimplantation Genetic Testing and Prenatal Diagnosis

Multiple osteochondromas (MO), the most common type of benign bone tumor, is an autosomal dominant skeletal disorder characterized by multiple cartilage-capped bony protuberances. In most cases, EXT1 and EXT2, which encode glycosyltransferases involved in the biosynthesis of heparan sulfate, are the genes responsible. Here we describe the clinical, phenotypic and genetic characterization of MO in 22 unrelated Chinese families involving a total of 60 patients. Variant detection was performed by means of a battery of different techniques including Sanger sequencing and whole-exome sequencing (WES). The pathogenicity of the missense and splicing variants was explored by means of in silico prediction algorithms. Sixteen unique pathogenic variants, including 10 in the EXT1 gene and 6 in the EXT2 gene, were identified in 18 (82%) of the 22 families. Fourteen (88%) of the 16 variants were predicted to give rise to truncated proteins whereas the remaining two were missense. Seven variants were newly described here, further expanding the spectrum of MO-causing variants in the EXT1 and EXT2 genes. More importantly, the identification of causative variants allowed us to provide genetic counseling to 8 MO patients in terms either of preimplantation genetic testing (PGT) or prenatal diagnosis, thereby preventing the reoccurrence of MO in the corresponding families. This study is the first to report the successful implementation of PGT in MO families and describes the largest number of subjects undergoing prenatal diagnosis to date.

Verifying nomenclature of DNA variants in submitted manuscripts: Guidance for journals

Documenting variation in our genomes is important for research and clinical care. Accuracy in the description of DNA variants is therefore essential. To address this issue, the Human Variome Project convened a committee to evaluate the feasibility of requiring authors to verify that all variants submitted for publication complied with a widely accepted standard for description. After a pilot study of two journals, the committee agreed that requiring authors to verify that variants complied with Human Genome Variation Society nomenclature is a reasonable step toward standardizing the worldwide inventory of human variation.

Inferring the molecular and phenotypic impact of amino acid variants with MutPred2

Identifying pathogenic variants and underlying functional alterations is challenging. To this end, we introduce MutPred2, a tool that improves the prioritization of pathogenic amino acid substitutions over existing methods, generates molecular mechanisms potentially causative of disease, and returns interpretable pathogenicity score distributions on individual genomes. Whilst its prioritization performance is state-of-the-art, a distinguishing feature of MutPred2 is the probabilistic modeling of variant impact on specific aspects of protein structure and function that can serve to guide experimental studies of phenotype-altering variants. We demonstrate the utility of MutPred2 in the identification of the structural and functional mutational signatures relevant to Mendelian disorders and the prioritization of de novo mutations associated with complex neurodevelopmental disorders. We then experimentally validate the functional impact of several variants identified in patients with such disorders. We argue that mechanism-driven studies of human inherited disease have the potential to significantly accelerate the discovery of clinically actionable variants.

Developmental Gene Expression Differences between Humans and Mammalian Models

Identifying the molecular programs underlying human organ development and how they differ from model species is key for understanding human health and disease. Developmental gene expression profiles provide a window into the genes underlying organ development and a direct means to compare them across species. We use a transcriptomic resource covering the development of seven organs to characterize the temporal profiles of human genes associated with distinct disease classes and to determine, for each human gene, the similarity of its spatiotemporal expression with its orthologs in rhesus macaque, mouse, rat, and rabbit. We find clear associations between spatiotemporal profiles and the phenotypic manifestations of diseases. We also find that half of human genes differ from their mouse orthologs in their temporal trajectories in at least one of the organs. These include more than 200 genes associated with brain, heart, and liver disease for which mouse models should undergo extra scrutiny.

Pathogenic and likely pathogenic variants in at least five genes account for approximately 3% of mild isolated nonsyndromic thrombocytopenia

BACKGROUND

Thrombocytopenia has a variety of different etiologies, both acquired and hereditary. Inherited thrombocytopenia may be associated with other symptoms (syndromic forms) or may be strictly isolated. To date, only about half of all the familial forms of thrombocytopenia have been accounted for in terms of well-defined genetic abnormalities. However, data are limited on the nature and frequency of the underlying causative genetic variants in individuals with mild isolated nonsyndromic thrombocytopenia.

STUDY DESIGN AND METHODS

Thirteen known or candidate genes for isolated thrombocytopenia were included in a gene panel analysis in which targeted next-generation sequencing was performed on 448 French blood donors with mild isolated nonsyndromic thrombocytopenia.

RESULTS

A total of 68 rare variants, including missense, splice site, frameshift, nonsense, and in-frame variants (all heterozygous) were identified in 11 of the 13 genes screened. Twenty-nine percent (N = 20) of the variants detected were absent from both the French Exome Project and gnomAD exome databases. Using stringent criteria and an unbiased approach, we classified seven predicted loss-of-function variants (three in ITGA2B and four in TUBB1) and four missense variants (one in GP1BA, two in ITGB3 and one in ACTN1) as being pathogenic or likely pathogenic. Altogether, they were found in 13 members (approx. 3%) of our studied cohort.

CONCLUSION

We present the results of gene panel sequencing of known and candidate thrombocytopenia genes in mild isolated nonsyndromic thrombocytopenia. Pathogenic and likely pathogenic variants in five known thrombocytopenia genes were identified, accounting for approximately 3% of individuals with the condition.

The Human Gene Mutation Database (HGMD®): optimizing its use in a clinical diagnostic or research setting

The Human Gene Mutation Database (HGMD^®) constitutes a comprehensive collection of published germline mutations in nuclear genes that are thought to underlie, or are closely associated with human inherited disease. At the time of writing (June 2020), the database contains in excess of 289,000 different gene lesions identified in over 11,100 genes manually curated from 72,987 articles published in over 3100 peer-reviewed journals. There are primarily two main groups of users who utilise HGMD on a regular basis; research scientists and clinical diagnosticians. This review aims to highlight how to make the most out of HGMD data in each setting.

Common homozygosity for predicted loss-of-function variants reveals both redundant and advantageous effects of dispensable human genes

Humans homozygous or hemizygous for variants predicted to cause a loss of function (LoF) of the corresponding protein do not necessarily present with overt clinical phenotypes. We report here 190 autosomal genes with 207 predicted LoF variants, for which the frequency of homozygous individuals exceeds 1% in at least one human population from five major ancestry groups. No such genes were identified on the X and Y chromosomes. Manual curation revealed that 28 variants (15%) had been misannotated as LoF. Of the 179 remaining variants in 166 genes, only 11 alleles in 11 genes had previously been confirmed experimentally to be LoF. The set of 166 dispensable genes was enriched in olfactory receptor genes (41 genes). The 41 dispensable olfactory receptor genes displayed a relaxation of selective constraints similar to that observed for other olfactory receptor genes. The 125 dispensable nonolfactory receptor genes also displayed a relaxation of selective constraints consistent with greater redundancy. Sixty-two of these 125 genes were found to be dispensable in at least three human populations, suggesting possible evolution toward pseudogenes. Of the 179 LoF variants, 68 could be tested for two neutrality statistics, and 8 displayed robust signals of positive selection. These latter variants included a known FUT2 variant that confers resistance to intestinal viruses, and an APOL3 variant involved in resistance to parasitic infections. Overall, the identification of 166 genes for which a sizeable proportion of humans are homozygous for predicted LoF alleles reveals both redundancies and advantages of such deficiencies for human survival.

Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

Gene-environment interaction between APOA5 c.553G>T and pregnancy in hypertriglyceridemia-induced acute pancreatitis

BACKGROUND

The etiology of hypertriglyceridemia (HTG) and, consequently, HTG-induced acute pancreatitis (HTG-AP), is complex.

OBJECTIVE

Herein, we explore a possible gene-environment interaction between APOA5 c.553G>T (p.185Gly>Cys, rs2075291), a common variant associated with altered triglyceride levels, and pregnancy in HTG-AP.

METHODS

We enrolled 318 Chinese HTG-AP patients and divided them into 3 distinct groups: Group 1, male patients (n = 183); Group 2, female patients whose disease was unrelated to pregnancy (n = 105); and Group 3, female patients whose disease was related to pregnancy (n = 30). APOA5 rs2075291 genotype status was determined by Sanger sequencing. A total of 362 healthy Han Chinese subjects were used as controls. Data on body mass index, peak triglyceride level, age of disease onset, episode number, and clinical severity of HTG-AP were collected from each patient. Multiple comparisons, between patient groups, between patient groups and controls, or within each patient group, were performed.

RESULTS

A robust association of APOA5 rs2075291 with HTG-AP in general, and HTG-AP during pregnancy in particular, was demonstrated. The minor T allele showed a stronger association with Group 3 patients than with either Group 1 or Group 2 patients. This stronger association was due mainly to the much higher frequency of TT genotype in Group 3 patients (20%) than that (<6%) in Group 1 and Group 2 patients. Moreover, the TT genotype was associated with a significantly higher peak triglyceride level in Group 3 patients compared with the GG genotype.

CONCLUSION

Our findings provide evidence for an interaction between APOA5 rs2075291 and pregnancy in HTG-AP.

SeqTailor: a user-friendly webserver for the extraction of DNA or protein sequences from next-generation sequencing data

Human whole-genome-sequencing reveals about 4 000 000 genomic variants per individual. These data are mostly stored as VCF-format files. Although many variant analysis methods accept VCF as input, many other tools require DNA or protein sequences, particularly for splicing prediction, sequence alignment, phylogenetic analysis, and structure prediction. However, there is no existing webserver capable of extracting DNA/protein sequences for genomic variants from VCF files in a user-friendly and efficient manner. We developed the SeqTailor webserver to bridge this gap, by enabling rapid extraction of (i) DNA sequences around genomic variants, with customizable window sizes and options to annotate the splice sites closest to the variants and to consider the neighboring variants within the window; and (ii) protein sequences encoded by the DNA sequences around genomic variants, with built-in SnpEff annotator and customizable window sizes. SeqTailor supports 11 species, including: human (GRCh37/GRCh38), chimpanzee, mouse, rat, cow, chicken, lizard, zebrafish, fruitfly, Arabidopsis and rice. Standalone programs are provided for command-line-based needs. SeqTailor streamlines the sequence extraction process, and accelerates the analysis of genomic variants with software requiring DNA/protein sequences. It will facilitate the study of genomic variation, by increasing the feasibility of sequence-based analysis and prediction. The SeqTailor webserver is freely available at http://shiva.rockefeller.edu/SeqTailor/.

RegulationSpotter: annotation and interpretation of extratranscriptic DNA variants

RegulationSpotter is a web-based tool for the user-friendly annotation and interpretation of DNA variants located outside of protein-coding transcripts (extratranscriptic variants). It is designed for clinicians and researchers who wish to assess the potential impact of the considerable number of non-coding variants found in Whole Genome Sequencing runs. It annotates individual variants with underlying regulatory features in an intuitive way by assessing over 100 genome-wide annotations. Additionally, it calculates a score, which reflects the regulatory potential of the variant region. Its dichotomous classifications, ‘functional’ or ‘non-functional’, and a human-readable presentation of the underlying evidence allow a biologically meaningful interpretation of the score. The output shows key aspects of every variant and allows rapid access to more detailed information about its possible role in gene regulation. RegulationSpotter can either analyse single variants or complete VCF files. Variants located within protein-coding transcripts are automatically assessed by MutationTaster as well as by RegulationSpotter to account for possible intragenic regulatory effects. RegulationSpotter offers the possibility of using phenotypic data to focus on known disease genes or genomic elements interacting with them. RegulationSpotter is freely available at https://www.regulationspotter.org.

5' splice site GC>GT and GT>GC variants differ markedly in terms of their functionality and pathogenicity

In the human genome, most 5' splice sites (~99%) employ the canonical GT dinucleotide whereas a small minority (~1%) use the noncanonical GC dinucleotide. The functionality and pathogenicity of 5' splice site GT>GC (+2T>C) variants have been extensively studied but we know very little about 5' splice site GC>GT (+2C>T) variants. Herein, we have addressed this deficiency by performing a meta-analysis of reported +2C>T "pathogenic" variants together with a functional analysis of engineered +2C>T substitutions using a cell culture-based full-length gene splicing assay. Our results establish proof of concept that +2C>T variants are qualitatively different from +2T>C variants in terms of their functionality and suggest that, in sharp contrast to +2T>C variants, most if not all +2C>T variants have no pathological relevance. Our findings have important implications for interpreting the clinical relevance of +2C>T variants and understanding the evolutionary switching between GT and GC 5' splice sites in mammalian genomes.

The Experimentally Obtained Functional Impact Assessments of 5' Splice Site GT'GC Variants Differ Markedly from Those Predicted

Introduction:

5' splice site GT>GC or +2T>C variants have been frequently reported to cause human genetic disease and are routinely scored as pathogenic splicing mutations. However, we have recently demonstrated that such variants in human disease genes may not invariably be pathogenic. Moreover, we found that no splicing prediction tools appear to be capable of reliably distinguishing those +2T>C variants that generate wild-type transcripts from those that do not.

Methodology

Herein, we evaluated the performance of a novel deep learning-based tool, SpliceAI, in the context of three datasets of +2T>C variants, all of which had been characterized functionally in terms of their impact on pre-mRNA splicing. The first two datasets refer to our recently described “in vivo” dataset of 45 known disease-causing +2T>C variants and the “in vitro” dataset of 103 +2T>C substitutions subjected to full-length gene splicing assay. The third dataset comprised 12 BRCA1 +2T>C variants that were recently analyzed by saturation genome editing.

Results

Comparison of the SpliceAI-predicted and experimentally obtained functional impact assessments of these variants (and smaller datasets of +2T>A and +2T>G variants) revealed that although SpliceAI performed rather better than other prediction tools, it was still far from perfect. A key issue was that the impact of those +2T>C (and +2T>A) variants that generated wild-type transcripts represents a quantitative change that can vary from barely detectable to an almost full expression of wild-type transcripts, with wild-type transcripts often co-existing with aberrantly spliced transcripts.

Conclusion

Our findings highlight the challenges that we still face in attempting to accurately identify splice-altering variants.

Identification and functional characterization of a novel heterozygous missense variant in the LPL associated with recurrent hypertriglyceridemia-induced acute pancreatitis in pregnancy

Background

Acute pancreatitis in pregnancy (APIP) is a life‐threatening disease for both mother and fetus. To date, only three patients with recurrent hypertriglyceridemia‐induced APIP (HTG‐APIP) have been reported to carry rare variants in the lipoprotein lipase (LPL) gene, which encodes the key enzyme responsible for triglyceride (TG) metabolism. Coincidently, all three patients harbored LPL variants on both alleles and presented with complete or severe LPL deficiency.

Methods

The entire coding regions and splice junctions of LPL and four other TG metabolism genes (APOC2, APOA5, GPIHBP1, and LMF1) were analyzed by Sanger sequencing in a Han Chinese patient who had experienced two episodes of HTG‐APIP. The impact of a novel LPL missense variant on LPL protein expression and activity was analyzed by transient expression in HEK293T cells.

Results

A novel heterozygous LPL missense variant, p.His210Leu (c.629A > T), was identified in our patient. This variant did not affect protein synthesis but significantly impaired LPL secretion and completely abolished the enzymatic activity of the mutant protein.

Conclusion

This report describes the first identification and functional characterization of a heterozygous variant in the LPL that predisposed to recurrent HTG‐APIP. Our findings confirm a major genetic contribution to the etiology of individual predisposition to HTG‐APIP.

Neuroprotectants attenuate hypobaric hypoxia-induced brain injuries in cynomolgus monkeys

Hypobaric hypoxia (HH) exposure can cause serious brain injury as well as life-threatening cerebral edema in severe cases. Previous studies on the mechanisms of HH-induced brain injury have been conducted primarily using non-primate animal models that are genetically distant to humans, thus hindering the development of disease treatment. Here, we report that cynomolgus monkeys (Macacafascicularis) exposed to acute HH developed human-like HH syndrome involving severe brain injury and abnormal behavior. Transcriptome profiling of white blood cells and brain tissue from monkeys exposed to increasing altitude revealed the central role of the HIF-1 and other novel signaling pathways, such as the vitamin D receptor (VDR) signaling pathway, in co-regulating HH-induced inflammation processes. We also observed profound transcriptomic alterations in brains after exposure to acute HH, including the activation of angiogenesis and impairment of aerobic respiration and protein folding processes, which likely underlie the pathological effects of HH-induced brain injury. Administration of progesterone (PROG) and steroid neuroprotectant 5α-androst-3β,5,6β-triol (TRIOL) significantly attenuated brain injuries and rescued the transcriptomic changes induced by acute HH. Functional investigation of the affected genes suggested that these two neuroprotectants protect the brain by targeting different pathways, with PROG enhancing erythropoiesis and TRIOL suppressing glutamate-induced excitotoxicity. Thus, this study advances our understanding of the pathology induced by acute HH and provides potential compounds for the development of neuroprotectant drugs for therapeutic treatment.

RegSNPs-intron: a computational framework for predicting pathogenic impact of intronic single nucleotide variants

Single nucleotide variants (SNVs) in intronic regions have yet to be systematically investigated for their disease-causing potential. Using known pathogenic and neutral intronic SNVs (iSNVs) as training data, we develop the RegSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure, and evolutionary conservation features. RegSNPs-intron showed excellent performance in evaluating the pathogenic impacts of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we evaluate the impact of RegSNPs-intron predictions on splicing outcome. Together, RegSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis.

Compound Heterozygosity for Novel Truncating Variants in the LMOD3 Gene as the Cause of Polyhydramnios in Two Successive Fetuses

Polyhydramnios is sometimes associated with genetic defects. However, establishing an accurate diagnosis and pinpointing the precise genetic cause of polyhydramnios in any given case represents a major challenge because it is known to occur in association with over 200 different conditions. Whole exome sequencing (WES) is now a routine part of the clinical workup, particularly with diseases characterized by atypical manifestations and significant genetic heterogeneity. Here we describe the identification, by means of WES, of novel compound heterozygous truncating variants in the LMOD3 gene [i.e., c.1412delA (p.Lys471Serfs*18) and c.1283dupC (p.Gly429Trpfs*35)] in a Chinese family with two successive fetuses affected with polyhydramnios, thereby potentiating the prenatal diagnosis of nemaline myopathy (NM) in the proband. LMOD3 encodes leiomodin-3, which is localized to the pointed ends of thin filaments and acts as a catalyst of actin nucleation in skeletal and cardiac muscle. This is the first study to describe the prenatal and postnatal manifestations of LMOD3-related NM in the Chinese population. Of all the currently reported NM-causing LMOD3 nonsense and frameshifting variants, c.1412delA generates the shortest truncation at the C-terminal end of the protein, underscoring the critical role of the WH2 domain in LMOD3 structure and function. Survey of the prenatal phenotypes of all known LMOD3-related severe NM cases served to identify fetal edema as a novel presenting feature that may provide an early clue to facilitate prenatal diagnosis of the disease.

Extensive disruption of protein interactions by genetic variants across the allele frequency spectrum in human populations

Each human genome carries tens of thousands of coding variants. The extent to which this variation is functional and the mechanisms by which they exert their influence remains largely unexplored. To address this gap, we leverage the ExAC database of 60,706 human exomes to investigate experimentally the impact of 2009 missense single nucleotide variants (SNVs) across 2185 protein-protein interactions, generating interaction profiles for 4797 SNV-interaction pairs, of which 421 SNVs segregate at > 1% allele frequency in human populations. We find that interaction-disruptive SNVs are prevalent at both rare and common allele frequencies. Furthermore, these results suggest that 10.5% of missense variants carried per individual are disruptive, a higher proportion than previously reported; this indicates that each individual's genetic makeup may be significantly more complex than expected. Finally, we demonstrate that candidate disease-associated mutations can be identified through shared interaction perturbations between variants of interest and known disease mutations.