Genetic evidence for T-wave area from 12-lead electrocardiograms to monitor cardiovascular diseases in patients taking diabetes medications

Aims Many studies indicated use of diabetes medications can influence the electrocardiogram (ECG), which remains the simplest and fastest tool for assessing cardiac functions. However, few studies have explored the role of genetic factors in determining the relationship between the use of diabetes medications and ECG trace characteristics (ETC). Methods Genome-wide association studies (GWAS) were performed for 168 ETCs extracted from the 12-lead ECGs of 42,340 Europeans in the UK Biobank. The genetic correlations, causal relationships, and phenotypic relationships of these ETCs with medication usage, as well as the risk of cardiovascular diseases (CVDs), were estimated by linkage disequilibrium score regression (LDSC), Mendelian randomization (MR), and regression model, respectively. Results The GWAS identified 124 independent single nucleotide polymorphisms (SNPs) that were study-wise and genome-wide significantly associated with at least one ETC. Regression model and LDSC identified significant phenotypic and genetic correlations of T-wave area in lead aVR (aVR_T-area) with usage of diabetes medications (ATC code: A10 drugs, and metformin), and the risks of ischemic heart disease (IHD) and coronary atherosclerosis (CA). MR analyses support a putative causal effect of the use of diabetes medications on decreasing aVR_T-area, and on increasing risk of IHD and CA. ConclusionPatients taking diabetes medications are prone to have decreased aVR_T-area and an increased risk of IHD and CA. The aVR_T-area is therefore a potential ECG marker for pre-clinical prediction of IHD and CA in patients taking diabetes medications.

Meta-analysis of 46,000 germline de novo mutations linked to human inherited disease

BACKGROUND

De novo mutations (DNMs) are variants that occur anew in the offspring of noncarrier parents. They are not inherited from either parent but rather result from endogenous mutational processes involving errors of DNA repair/replication. These spontaneous errors play a significant role in the causation of genetic disorders, and their importance in the context of molecular diagnostic medicine has become steadily more apparent as more DNMs have been reported in the literature. In this study, we examined 46,489 disease-associated DNMs annotated by the Human Gene Mutation Database (HGMD) to ascertain their distribution across gene and disease categories.

RESULTS

Most disease-associated DNMs reported to date are found to be associated with developmental and psychiatric disorders, a reflection of the focus of sequencing efforts over the last decade. Of the 13,277 human genes in which DNMs have so far been found, the top-10 genes with the highest proportions of DNM relative to gene size were H3-3 A, DDX3X, CSNK2B, PURA, ZC4H2, STXBP1, SCN1A, SATB2, H3-3B and TUBA1A. The distribution of CADD and REVEL scores for both disease-associated DNMs and those mutations not reported to be de novo revealed a trend towards higher deleteriousness for DNMs, consistent with the likely lower selection pressure impacting them. This contrasts with the non-DNMs, which are presumed to have been subject to continuous negative selection over multiple generations.

CONCLUSION

This meta-analysis provides important information on the occurrence and distribution of disease-associated DNMs in association with heritable disease and should make a significant contribution to our understanding of this major type of mutation.

CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

Divergent Evolutionary Rates of Primate Brain Regions as Revealed by Genomics and Transcriptomics

Although the primate brain contains numerous functionally distinct structures that have experienced diverse genetic changes during the course of evolution and development, these changes remain to be explored in detail. Here we utilize two classic metrics from evolutionary biology, the evolutionary rate index (ERI) and the transcriptome age index (TAI), to investigate the evolutionary alterations that have occurred in each area and developmental stage of the primate brain. We observed a higher evolutionary rate for those genes expressed in the non-cortical areas during primate evolution, particularly in human, with the highest rate of evolution being exhibited at brain developmental stages between late infancy and early childhood. Further, the transcriptome age of the non-cortical areas was lower than that of the cerebral cortex, with the youngest age apparent at brain developmental stages between late infancy and early childhood. Our exploration of the evolutionary patterns manifest in each brain area and developmental stage provides important reference points for further research into primate brain evolution.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

Genome-wide prediction of pathogenic gain- and loss-of-function variants from ensemble learning of a diverse feature set

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available at https://itanlab.shinyapps.io/goflof/ .

Genome-wide identification of dominant polyadenylation hexamers for use in variant classification

Polyadenylation is an essential process for the stabilization and export of mRNAs to the cytoplasm and the polyadenylation signal hexamer (herein referred to as hexamer) plays a key role in this process. Yet, only 14 Mendelian disorders have been associated with hexamer variants. This is likely an under-ascertainment as hexamers are not well defined and not routinely examined in molecular analysis. To facilitate the interrogation of putatively pathogenic hexamer variants, we set out to define functionally important hexamers genome-wide as a resource for research and clinical testing interrogation. We identified predominant polyA sites (herein referred to as pPAS) and putative predominant hexamers across protein coding genes (PAS usage >50% per gene). As a measure of the validity of these sites, the population constraint of 4532 predominant hexamers were measured. The predominant hexamers had fewer observed variants compared to non-predominant hexamers and trimer controls, and CADD scores for variants in these hexamers were significantly higher than controls. Exome data for 1477 individuals were interrogated for hexamer variants and transcriptome data were generated for 76 individuals with 65 variants in predominant hexamers. 3' RNA-seq data showed these variants resulted in alternate polyadenylation events (38%) and in elongated mRNA transcripts (12%). Our list of pPAS and predominant hexamers are available in the UCSC genome browser and on GitHub. We suggest this list of predominant hexamers can be used to interrogate exome and genome data. Variants in these predominant hexamers should be considered candidates for pathogenic variation in human disease, and to that end we suggest pathogenicity criteria for classifying hexamer variants.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

Inferring the genetic relationship between brain imaging-derived phenotypes and risk of complex diseases by Mendelian randomization and genome-wide colocalization

Observational studies consistently disclose brain imaging-derived phenotypes (IDPs) as critical markers for early diagnosis of both brain disorders and cardiovascular diseases. However, it remains unclear about the shared genetic landscape between brain IDPs and the risk of brain disorders and cardiovascular diseases, restricting the applications of potential diagnostic techniques through brain IDPs. Here, we reported genetic correlations and putative causal relationships between 921 brain IDPs, 20 brain disorders and six cardiovascular diseases by leveraging their large-scale genome-wide association study (GWAS) summary statistics. Applications of Mendelian randomization (MR) identified significant putative causal effects of multiple region-specific brain IDPs in relation to the increased risks for amyotrophic lateral sclerosis (ALS), major depressive disorder (MDD), autism spectrum disorder (ASD) and schizophrenia (SCZ). We also found brain IDPs specifically from temporal lobe as a putatively causal consequence of hypertension. The genome-wide colocalization analysis identified three genomic regions in which MDD, ASD and SCZ colocalized with the brain IDPs, and two novel SNPs to be associated with ASD, SCZ, and multiple brain IDPs. Furthermore, we identified a list of candidate genes involved in the shared genetics underlying pairs of brain IDPs and MDD, ASD, SCZ, ALS and hypertension. Our results provide novel insights into the genetic relationships between brain disorders and cardiovascular diseases and brain IDP, which may server as clues for using brain IDPs to predict risks of diseases.

Frameshift coding sequence variants in the LPL gene: identification of two novel events and exploration of the genotype-phenotype relationship for variants reported to date

BACKGROUND

Lipoprotein lipase (LPL) is the rate-limiting enzyme for triglyceride hydrolysis. Homozygous or compound heterozygous LPL variants cause autosomal recessive familial chylomicronemia syndrome (FCS), whereas simple heterozygous LPL variants are associated with hypertriglyceridemia (HTG) and HTG-related disorders. LPL frameshift coding sequence variants usually cause complete functional loss of the affected allele, thereby allowing exploration of the impact of different levels of LPL function in human disease.

METHODS

All exons and flanking intronic regions of LPL were Sanger sequenced in patients with HTG-related acute pancreatitis (HTG-AP) or HTG-AP in pregnancy. Previously reported LPL frameshift coding sequence variants were collated from the Human Gene Mutation Database and through PubMed keyword searching. Original reports were manually evaluated for the following information: zygosity status of the variant, plasma LPL activity of the variant carrier, disease referred for genetic analysis, patient's age at genetic analysis, and patient's disease history. SpliceAI was employed to predict the potential impact of collated variants on splicing.

RESULTS

Two novel rare variants were identified, and 53 known LPL frameshift coding sequence variants were collated. Of the 51 variants informative for zygosity, 30 were simple heterozygotes, 12 were homozygotes, and 9 were compound heterozygotes. Careful evaluation of the 55 variants with respect to their clinical and genetic data generated several interesting findings. First, we conclude that 6-7% residual LPL function could significantly delay the age of onset of FCS and reduce the prevalence of FCS-associated syndromes. Second, whereas a large majority of LPL frameshift coding sequence variants completely disrupt gene function through their "frameshift" nature, a small fraction of these variants may act wholly or partly as "in-frame" variants, leading to the generation of protein products with some residual LPL function. Third, we identified two candidate LPL frameshift coding sequence variants that may retain residual function based on genotype-phenotype correlation or SpliceAI-predicted data.

CONCLUSIONS

This study reported two novel LPL variants and yielded new insights into the genotype-phenotype relationship as it pertains to LPL frameshift coding sequence variants.

The East Asian-specific LPL p.Ala288Thr (c.862G > A) missense variant exerts a mild effect on protein function

BACKGROUND

Lipoprotein lipase (LPL) is the key enzyme responsible for the hydrolysis of triglycerides. Loss-of-function variants in the LPL gene are associated with hypertriglyceridemia (HTG) and HTG-related diseases. Unlike nonsense, frameshift and canonical GT-AG splice site variants, a pathogenic role for clinically identified LPL missense variants should generally be confirmed by functional analysis. Herein, we describe the clinical and functional analysis of a rare LPL missense variant.

METHODS

Chinese patients with HTG-associated acute pancreatitis (HTG-AP) were screened for rare nonsense, frameshift, missense or canonical GT-AG splice site variants in LPL and four other lipid metabolism-related genes (APOC2, APOA5, GPIHBP1 and LMF1) by Sanger sequencing. The functional consequences of the LPL missense variant of interest were characterized by in vitro expression in HEK-293T and COS-7 cells followed by Western blot and LPL activity assays.

RESULTS

Five unrelated HTG-AP patients were found to be heterozygous for a rare East Asian-specific LPL missense variant, c.862G > A (p.Ala288Thr). All five patients were adult males, and all were overweight and had a long history of alcohol consumption. Transfection of LPL wild-type and c.862G > A expression vectors into two cell lines followed by Western blot analysis served to exclude the possibility that the p.Ala288Thr missense variant either impaired protein synthesis or increased protein degradation. Contrary to a previous functional study that claimed that p.Ala288Thr had a severe impact on LPL function (reportedly having 36% normal activity), our experiments consistently demonstrated that the variant had a comparatively mild effect on LPL functional activity, which was mediated through its impact upon LPL protein secretion (~ 20% reduced secretion compared to wild-type).

CONCLUSIONS

In this study, we identified the East Asian-specific LPL c.862G > A (p.Ala288Thr) missense variant in five unrelated HTG-AP patients. We demonstrated that this variant exerted only a relatively mild effect on LPL function in two cell lines. Heterozygosity for this LPL variant may have combined with alcohol consumption to trigger HTG-AP in these patients.

Integrative Omics Reveals Rapidly Evolving Regulatory Sequences Driving Primate Brain Evolution

Although the continual expansion of the brain during primate evolution accounts for our enhanced cognitive capabilities, the drivers of brain evolution have scarcely been explored in these ancestral nodes. Here, we performed large-scale comparative genomic, transcriptomic, and epigenomic analyses to investigate the evolutionary alterations acquired by brain genes and provide comprehensive listings of innovatory genetic elements along the evolutionary path from ancestral primates to human. The regulatory sequences associated with brain-expressed genes experienced rapid change, particularly in the ancestor of the Simiiformes. Extensive comparisons of single-cell and bulk transcriptomic data between primate and nonprimate brains revealed that these regulatory sequences may drive the high expression of certain genes in primate brains. Employing in utero electroporation into mouse embryonic cortex, we show that the primate-specific brain-biased gene BMP7 was recruited, probably in the ancestor of the Simiiformes, to regulate neuronal proliferation in the primate ventricular zone. Our study provides a comprehensive listing of genes and regulatory changes along the brain evolution lineage of ancestral primates leading to human. These data should be invaluable for future functional studies that will deepen our understanding not only of the genetic basis of human brain evolution but also of inherited disease.

Identification of protease-sensitive but not misfolding PNLIP variants in familial and hereditary pancreatitis

Mutations in the PNLIP gene have recently been implicated in chronic pancreatitis. Several PNLIP missense variants have been reported to cause protein misfolding and endoplasmic reticulum stress although genetic evidence supporting their association with chronic pancreatitis is currently lacking. Protease-sensitive PNLIP missense variants have also been associated with early-onset chronic pancreatitis although the underlying pathological mechanism remains enigmatic. Herein, we provide new evidence to support the association of protease-sensitive PNLIP variants (but not misfolding PNLIP variants) with pancreatitis. Specifically, we identified protease-sensitive PNLIP variants in 5 of 373 probands (1.3%) with a positive family history of pancreatitis. The protease-sensitive variants, p.F300L and p.I265R, were found to segregate with the disease in three families, including one exhibiting a classical autosomal dominant inheritance pattern. Consistent with previous findings, protease-sensitive variant-positive patients were often characterized by early-onset disease and invariably experienced recurrent acute pancreatitis, although none has so far developed chronic pancreatitis.

Classification of PRSS1 variants responsible for chronic pancreatitis: An expert perspective from the Franco-Chinese GREPAN study group

BACKGROUND

PRSS1 was the first reported chronic pancreatitis (CP) gene. The existence of both gain-of-function (GoF) and gain-of-proteotoxicity (GoP) pathological PRSS1 variants, together with the fact that PRSS1 variants have been identified in CP subtypes spanning the range from monogenic to multifactorial, has made the classification of PRSS1 variants very challenging.

METHODS

All currently reported PRSS1 variants (derived primarily from two databases) were manually reviewed with respect to their clinical genetics, functional analysis and population allele frequency. They were classified by variant type and pathological mechanism within the framework of our recently proposed ACMG/AMP guidelines-based seven-category system.

RESULTS

The total number of distinct germline PRSS1 variants included for analysis was 100, comprising 3 copy number variants (CNVs), 12 5' and 3' variants, 19 intronic variants, 5 nonsense variants, 1 frameshift deletion variant, 6 synonymous variants, 1 in-frame duplication, 3 gene conversions and 50 missense variants. Based upon a combination of clinical genetic and functional analysis, population data and in silico analysis, we classified 26 variants (all 3 CNVs, the in-frame duplication, all 3 gene conversions and 19 missense) as "pathogenic", 3 variants (missense) as "likely pathogenic", 5 variants (four missense and one promoter) as "predisposing", 13 variants (all missense) as "unknown significance", 2 variants (missense) as "likely benign", and all remaining 51 variants as "benign".

CONCLUSIONS

We describe an expert classification of the 100 PRSS1 variants reported to date. The results have immediate implications for reclassifying many ClinVar-registered PRSS1 variants as well as providing optimal guidelines/standards for reporting PRSS1 variants.

Eighty million years of rapid evolution of the primate Y chromosome

The Y chromosome usually plays a critical role in determining male sex and comprises sequence classes that have experienced unique evolutionary trajectories. Here we generated 19 new primate sex chromosome assemblies, analysed them with 10 existing assemblies and report rapid evolution of the Y chromosome across primates. The pseudoautosomal boundary has shifted at least six times during primate evolution, leading to the formation of a Simiiformes-specific evolutionary stratum and to the independent start of young strata in Catarrhini and Platyrrhini. Different primate lineages experienced different rates of gene loss and structural and chromatin change on their Y chromosomes. Selection on several Y-linked genes has contributed to the evolution of male developmental traits across the primates. Additionally, lineage-specific expansions of ampliconic regions have further increased the diversification of the structure and gene composition of the Y chromosome. Overall, our comprehensive analysis has broadened our knowledge of the evolution of the primate Y chromosome.

Lineage-specific accelerated sequences underlying primate evolution

Understanding the mechanisms underlying phenotypic innovation is a key goal of comparative genomic studies. Here, we investigated the evolutionary landscape of lineage-specific accelerated regions (LinARs) across 49 primate species. Genomic comparison with dense taxa sampling of primate species significantly improved LinAR detection accuracy and revealed many novel human LinARs associated with brain development or disease. Our study also yielded detailed maps of LinARs in other primate lineages that may have influenced lineage-specific phenotypic innovation and adaptation. Functional experimentation identified gibbon LinARs, which could have participated in the developmental regulation of their unique limb structures, whereas some LinARs in the Colobinae were associated with metabolite detoxification which may have been adaptive in relation to their leaf-eating diet. Overall, our study broadens knowledge of the functional roles of LinARs in primate evolution.

Comparative genomics reveals the hybrid origin of a macaque group

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Phylogenomic analyses provide insights into primate evolution

Comparative analysis of primate genomes within a phylogenetic context is essential for understanding the evolution of human genetic architecture and primate diversity. We present such a study of 50 primate species spanning 38 genera and 14 families, including 27 genomes first reported here, with many from previously less well represented groups, the New World monkeys and the Strepsirrhini. Our analyses reveal heterogeneous rates of genomic rearrangement and gene evolution across primate lineages. Thousands of genes under positive selection in different lineages play roles in the nervous, skeletal, and digestive systems and may have contributed to primate innovations and adaptations. Our study reveals that many key genomic innovations occurred in the Simiiformes ancestral node and may have had an impact on the adaptive radiation of the Simiiformes and human evolution.

Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

Identifying shared genetic factors underlying epilepsy and congenital heart disease in Europeans

Epilepsy (EP) and congenital heart disease (CHD) are two apparently unrelated diseases that nevertheless display substantial mutual comorbidity. Thus, while congenital heart defects are associated with an elevated risk of developing epilepsy, the incidence of epilepsy in CHD patients correlates with CHD severity. Although genetic determinants have been postulated to underlie the comorbidity of EP and CHD, the precise genetic etiology is unknown. We performed variant and gene association analyses on EP and CHD patients separately, using whole exomes of genetically identified Europeans from the UK Biobank and Mount Sinai BioMe Biobank. We prioritized biologically plausible candidate genes and investigated the enriched pathways and other identified comorbidities by biological proximity calculation, pathway analyses, and gene-level phenome-wide association studies. Our variant- and gene-level results point to the Voltage-Gated Calcium Channels (VGCC) pathway as being a unifying framework for EP and CHD comorbidity. Additionally, pathway-level analyses indicated that the functions of disease-associated genes partially overlap between the two disease entities. Finally, phenome-wide association analyses of prioritized candidate genes revealed that cerebral blood flow and ulcerative colitis constitute the two main traits associated with both EP and CHD.

Profiling human pathogenic repeat expansion regions by synergistic and multi-level impacts on molecular connections

Whilst DNA repeat expansions cause numerous heritable human disorders, their origins and underlying pathological mechanisms are often unclear. We collated a dataset comprising 224 human repeat expansions encompassing 203 different genes, and performed a systematic analysis with respect to key topological features at the DNA, RNA and protein levels. Comparison with controls without known pathogenicity and genomic regions lacking repeats, allowed the construction of the first tool to discriminate repeat regions harboring pathogenic repeat expansions (DPREx). At the DNA level, pathogenic repeat expansions exhibited stronger signals for DNA regulatory factors (e.g. H3K4me3, transcription factor-binding sites) in exons, promoters, 5'UTRs and 5'genes but were not significantly different from controls in introns, 3'UTRs and 3'genes. Additionally, pathogenic repeat expansions were also found to be enriched in non-B DNA structures. At the RNA level, pathogenic repeat expansions were characterized by lower free energy for forming RNA secondary structure and were closer to splice sites in introns, exons, promoters and 5'genes than controls. At the protein level, pathogenic repeat expansions exhibited a preference to form coil rather than other types of secondary structure, and tended to encode surface-located protein domains. Guided by these features, DPREx ( http://biomed.nscc-gz.cn/zhaolab/geneprediction/# ) achieved an Area Under the Curve (AUC) value of 0.88 in a test on an independent dataset. Pathogenic repeat expansions are thus located such that they exert a synergistic influence on the gene expression pathway involving inter-molecular connections at the DNA, RNA and protein levels.

Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis)

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

The PRSS3P2 and TRY7 deletion copy number variant modifies risk for chronic pancreatitis

BACKGROUND

PRSS1 and PRSS2 constitute the only functional copies of a tandemly-arranged five-trypsinogen-gene cluster (i.e., PRSS1, PRSS3P1, PRSS3P2, TRY7 and PRSS2) on chromosome 7q35. Variants in PRSS1 and PRSS2, including missense and copy number variants (CNVs), have been reported to predispose to or protect against chronic pancreatitis (CP). We wondered whether a common trypsinogen pseudogene deletion CNV (that removes two of the three trypsinogen pseudogenes, PRSS3P2 and TRY7) might be associated with CP causation/predisposition.

METHODS

We analyzed the common PRSS3P2 and TRY7 deletion CNV in a total of 1536 CP patients and 3506 controls from France, Germany, India and Japan by means of quantitative fluorescent multiplex polymerase chain reaction.

RESULTS

We demonstrated that the deletion CNV variant was associated with a protective effect against CP in the French, German and Japanese cohorts whilst a trend toward the same association was noted in the Indian cohort. Meta-analysis under a dominant model yielded a pooled odds ratio (OR) of 0.68 (95% confidence interval (CI) 0.52-0.89; p = 0.005) whereas an allele-based meta-analysis yielded a pooled OR of 0.84 (95% CI 0.77-0.92; p = 0.0001). This protective effect is explicable by reference to the recent finding that the still functional PRSS3P2/TRY7 pseudogene enhancers upregulate pancreatic PRSS2 expression.

CONCLUSIONS

The common PRSS3P2 and TRY7 deletion CNV was associated with a reduced risk for CP. This finding provides additional support for the emerging view that dysregulated PRSS2 expression represents a discrete mechanism underlying CP predisposition or protection.

Genetic evidence for a causal relationship between type 2 diabetes and peripheral artery disease in both Europeans and East Asians

BACKGROUND

Observational studies have revealed that type 2 diabetes (T2D) is associated with an increased risk of peripheral artery disease (PAD). However, whether the two diseases share a genetic basis and whether the relationship is causal remain unclear. It is also unclear as to whether these relationships differ between ethnic groups.

METHODS

By leveraging large-scale genome-wide association study (GWAS) summary statistics of T2D (European-based: N_case = 21,926, N_control = 342,747; East Asian-based: N_case = 36,614, N_control = 155,150) and PAD (European-based: N_case = 5673, N_control = 359,551; East Asian-based: N_case = 3593, N_control = 208,860), we explored the genetic correlation and putative causal relationship between T2D and PAD in both Europeans and East Asians using linkage disequilibrium score regression and seven Mendelian randomization (MR) models. We also performed multi-trait analysis of GWAS and two gene-based analyses to reveal candidate variants and risk genes involved in the shared genetic basis between T2D and PAD.

RESULTS

We observed a strong genetic correlation (r_g) between T2D and PAD in both Europeans (r_g = 0.51; p-value = 9.34 × 10^-15) and East Asians (r_g = 0.46; p-value = 1.67 × 10^-12). The MR analyses provided consistent evidence for a causal effect of T2D on PAD in both ethnicities (odds ratio [OR] = 1.05 to 1.28 for Europeans and 1.15 to 1.27 for East Asians) but not PAD on T2D. This putative causal effect was not influenced by total cholesterol, body mass index, systolic blood pressure, or smoking initiation according to multivariable MR analysis, and the genetic overlap between T2D and PAD was further explored employing an independent European sample through polygenic risk score regression. Multi-trait analysis of GWAS revealed two novel European-specific single nucleotide polymorphisms (rs927742 and rs1734409) associated with the shared genetic basis of T2D and PAD. Gene-based analyses consistently identified one gene ANKFY1 and gene-gene interactions (e.g., STARD10 [European-specific] to AP3S2 [East Asian-specific]; KCNJ11 [European-specific] to KCNQ1 [East Asian-specific]) associated with the trans-ethnic genetic overlap between T2D and PAD, reflecting a common genetic basis for the co-occurrence of T2D and PAD in both Europeans and East Asians.

CONCLUSIONS

Our study provides the first evidence for a genetically causal effect of T2D on PAD in both Europeans and East Asians. Several candidate variants and risk genes were identified as being associated with this genetic overlap. Our findings emphasize the importance of monitoring PAD status in T2D patients and suggest new genetic biomarkers for screening PAD risk among patients with T2D.

Expanding ACMG variant classification guidelines into a general framework

BACKGROUND

The American College of Medical Genetics and Genomics (ACMG)-recommended five variant classification categories (pathogenic, likely pathogenic, uncertain significance, likely benign, and benign) have been widely used in medical genetics. However, these guidelines are fundamentally constrained in practice owing to their focus upon Mendelian disease genes and their dichotomous classification of variants as being either causal or not. Herein, we attempt to expand the ACMG guidelines into a general variant classification framework that takes into account not only the continuum of clinical phenotypes, but also the continuum of the variants' genetic effects, and the different pathological roles of the implicated genes.

MAIN BODY

As a disease model, we employed chronic pancreatitis (CP), which manifests clinically as a spectrum from monogenic to multifactorial. Bearing in mind that any general conceptual proposal should be based upon sound data, we focused our analysis on the four most extensively studied CP genes, PRSS1, CFTR, SPINK1 and CTRC. Based upon several cross-gene and cross-variant comparisons, we first assigned the different genes to two distinct categories in terms of disease causation: CP-causing (PRSS1 and SPINK1) and CP-predisposing (CFTR and CTRC). We then employed two new classificatory categories, "predisposing" and "likely predisposing", to replace ACMG's "pathogenic" and "likely pathogenic" categories in the context of CP-predisposing genes, thereby classifying all pathologically relevant variants in these genes as "predisposing". In the case of CP-causing genes, the two new classificatory categories served to extend the five ACMG categories whilst two thresholds (allele frequency and functional) were introduced to discriminate "pathogenic" from "predisposing" variants.

CONCLUSION

Employing CP as a disease model, we expand ACMG guidelines into a five-category classification system (predisposing, likely predisposing, uncertain significance, likely benign, and benign) and a seven-category classification system (pathogenic, likely pathogenic, predisposing, likely predisposing, uncertain significance, likely benign, and benign) in the context of disease-predisposing and disease-causing genes, respectively. Taken together, the two systems constitute a general variant classification framework that, in principle, should span the entire spectrum of variants in any disease-related gene. The maximal compliance of our five-category and seven-category classification systems with the ACMG guidelines ought to facilitate their practical application.

X-CAP improves pathogenicity prediction of stopgain variants

Stopgain substitutions are the third-largest class of monogenic human disease mutations and often examined first in patient exomes. Existing computational stopgain pathogenicity predictors, however, exhibit poor performance at the high sensitivity required for clinical use. Here, we introduce a new classifier, termed X-CAP, which uses a novel training methodology and unique feature set to improve the AUROC by 18% and decrease the false-positive rate 4-fold on large variant databases. In patient exomes, X-CAP prioritizes causal stopgains better than existing methods do, further illustrating its clinical utility. X-CAP is available at https://github.com/bejerano-lab/X-CAP .

The CEL-HYB1 Hybrid Allele Promotes Digestive Enzyme Misfolding and Pancreatitis in Mice

BACKGROUND & AIMS

A hybrid allele that originated from homologous recombination between CEL and its pseudogene (CELP), CEL-HYB1 increases the risk of chronic pancreatitis (CP). Although suggested to cause digestive enzyme misfolding, definitive in vivo evidence for this postulate has been lacking.

METHODS

CRISPR-Cas9 was used to generate humanized mice harboring the CEL-HYB1 allele on a C57BL/6J background. Humanized CEL mice and C57BL/6J mice were used as controls. Pancreata were collected and analyzed by histology, immunohistochemistry, immunoblotting, and transcriptomics. Isolated pancreatic acini were cultured in vitro to measure the secretion and aggregation of CEL-HYB1 protein. Mice were given caerulein injections to induce acute pancreatitis (AP) and CP.

RESULTS

Pancreata from mice expressing CEL-HYB1 developed pathological features characteristic of focal pancreatitis that included acinar atrophy and vacuolization, inflammatory infiltrates, and fibrosis in a time-dependent manner. CEL-HYB1 expression in pancreatic acini led to decreased secretion and increased intracellular aggregation and triggered endoplasmic reticulum stress compared with CEL. The autophagy levels of pancreata from mice expressing CEL-HYB1 changed at different developmental stages; some aged CEL-HYB1 mice exhibited an accumulation of large autophagic vesicles and impaired autophagy in acinar cells. Administration of caerulein increased the severity of AP/CP in mice expressing CEL-HYB1 compared with control mice, accompanied by higher levels of endoplasmic reticulum stress.

CONCLUSIONS

Expression of a humanized form of CEL-HYB1 in mice promotes endoplasmic reticulum stress and pancreatitis through a misfolding-dependent pathway. Impaired autophagy appears to be involved in the pancreatic injury in aged CEL-HYB1 mice. These mice have the potential to be used as a model to identify therapeutic targets for CP.

Trypsinogen (PRSS1 and PRSS2) gene dosage correlates with pancreatitis risk across genetic and transgenic studies: a systematic review and re-analysis

Trypsinogen (PRSS1, PRSS2) copy number gains and regulatory variants have both been proposed to elevate pancreatitis risk through a gene dosage effect (i.e., by increasing the expression of wild-type protein). However, to date, their impact on pancreatitis risk has not been thoroughly evaluated whilst the underlying pathogenic mechanisms remain to be explicitly investigated in mouse models. Genetic studies of the rare trypsinogen duplication and triplication copy number variants (CNVs), and the common rs10273639C variant, were collated from PubMed and/or ClinVar. Mouse studies that analyzed the influence of a transgenically expressed wild-type human PRSS1 or PRSS2 gene on the development of pancreatitis were identified from PubMed. The genetic effects of the different risk genotypes, in terms of odds ratios, were calculated wherever appropriate. The genetic effects of the rare trypsinogen duplication and triplication CNVs were also evaluated by reference to their associated disease subtypes. We demonstrate a positive correlation between increased trypsinogen gene dosage and pancreatitis risk in the context of the rare duplication and triplication CNVs, and between the level of trypsinogen expression and disease risk in the context of the heterozygous and homozygous rs10273639C-tagged genotypes. We retrospectively identify three mouse transgenic studies that are informative in relation to the pathogenic mechanism underlying the trypsinogen gene dosage effect in pancreatitis. Trypsinogen gene dosage correlates with pancreatitis risk across genetic and transgenic studies, highlighting the fundamental role of dysregulated expression of wild-type trypsinogen in the etiology of pancreatitis. Specifically downregulating trypsinogen expression in the pancreas may serve as a potential therapeutic and/or prevention strategy for pancreatitis.

Challenges in the diagnosis of neurofibromatosis type 1 (NF1) in young children facilitated by means of revised diagnostic criteria including genetic testing for pathogenic NF1 gene variants

Neurofibromatosis type 1 (NF1) is the most frequent disorder associated with multiple café-au-lait macules (CALM) which may either be present at birth or appear during the first year of life. Other NF1-associated features such as skin-fold freckling and Lisch nodules occur later during childhood whereas dermal neurofibromas are rare in young children and usually only arise during early adulthood. The NIH clinical diagnostic criteria for NF1, established in 1988, include the most common NF1-associated features. Since many of these features are age-dependent, arriving at a definitive diagnosis of NF1 by employing these criteria may not be possible in infancy if CALM are the only clinical feature evident. Indeed, approximately 46% of patients who are diagnosed with NF1 later in life do not meet the NIH diagnostic criteria by the age of 1 year. Further, the 1988 diagnostic criteria for NF1 are not specific enough to distinguish NF1 from other related disorders such as Legius syndrome. In this review, we outline the challenges faced in diagnosing NF1 in young children, and evaluate the utility of the recently revised (2021) diagnostic criteria for NF1, which include the presence of pathogenic variants in the NF1 gene and choroidal anomalies, for achieving an early and accurate diagnosis.

Functionally deficient TRPV6 variants contribute to hereditary and familial chronic pancreatitis

The recent discovery of TRPV6 as a pancreatitis susceptibility gene served to identify a novel mechanism of chronic pancreatitis (CP) due to Ca²⁺ dysregulation. Herein, we analyzed TRPV6 in 81 probands with hereditary CP (HCP), 204 probands with familial CP (FCP), and 462 patients with idiopathic CP (ICP) by targeted next-generation sequencing. We identified 25 rare nonsynonymous TRPV6 variants, 18 of which had not been previously reported. All 18 variants were characterized by a Ca²⁺ imaging assay, with 8 being identified as functionally deficient. Evaluation of functionally deficient variants in the three CP cohorts revealed two novel findings: (i) functionally deficient TRPV6 variants appear to occur more frequently in HCP/FCP patients than in ICP patients (3.2% vs. 1.5%) and (ii) functionally deficient TRPV6 variants found in HCP and FCP probands appear to be more frequently coinherited with known risk variants in SPINK1, CTRC, and/or CFTR than those found in ICP patients (66.7% vs 28.6%). Additionally, genetic analysis of available HCP and FCP family members revealed complex patterns of inheritance in some families. Our findings confirm that functionally deficient TRPV6 variants represent an important contributor to CP. Importantly, functionally deficient TRPV6 variants account for a significant proportion of cases of HCP/FCP.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

Atypical NF1 Microdeletions: Challenges and Opportunities for Genotype/Phenotype Correlations in Patients with Large NF1 Deletions

Patients with neurofibromatosis type 1 (NF1) and type 1 NF1 deletions often exhibit more severe clinical manifestations than patients with intragenic NF1 gene mutations, including facial dysmorphic features, overgrowth, severe global developmental delay, severe autistic symptoms and considerably reduced cognitive abilities, all of which are detectable from a very young age. Type 1 NF1 deletions encompass 1.4 Mb and are associated with the loss of 14 protein-coding genes, including NF1 and SUZ12. Atypical NF1 deletions, which do not encompass all 14 protein-coding genes located within the type 1 NF1 deletion region, have the potential to contribute to the delineation of the genotype/phenotype relationship in patients with NF1 microdeletions. Here, we review all atypical NF1 deletions reported to date as well as the clinical phenotype observed in the patients concerned. We compare these findings with those of a newly identified atypical NF1 deletion of 698 kb which, in addition to the NF1 gene, includes five genes located centromeric to NF1. The atypical NF1 deletion in this patient does not include the SUZ12 gene but does encompass CRLF3. Comparative analysis of such atypical NF1 deletions suggests that SUZ12 hemizygosity is likely to contribute significantly to the reduced cognitive abilities, severe global developmental delay and facial dysmorphisms observed in patients with type 1 NF1 deletions.

Classification of NF1 microdeletions and its importance for establishing genotype/phenotype correlations in patients with NF1 microdeletions

An estimated 5–11% of patients with neurofibromatosis type-1 (NF1) harbour large deletions encompassing the NF1 gene and flanking regions. These NF1 microdeletions are subclassified into type 1, 2, 3 and atypical deletions which are distinguishable from each other by their extent and by the number of genes included within the deletion regions as well as the frequency of mosaicism with normal cells. Most common are type-1 NF1 deletions which encompass 1.4-Mb and 14 protein-coding genes. Type-1 deletions are frequently associated with overgrowth, global developmental delay, cognitive disability and dysmorphic facial features which are uncommon in patients with intragenic pathogenic NF1 gene variants. Further, patients with type-1 NF1 deletions frequently exhibit high numbers of neurofibromas and have an increased risk of malignant peripheral nerve sheath tumours. Genes located within the type-1 NF1 microdeletion interval and co-deleted with NF1 are likely to act as modifiers responsible for the severe disease phenotype in patients with NF1 microdeletions, thereby causing the NF1 microdeletion syndrome. Genotype/phenotype correlations in patients with NF1 microdeletions of different lengths are important to identify such modifier genes. However, these correlations are critically dependent upon the accurate characterization of the deletions in terms of their extent. In this review, we outline the utility as well as the shortcomings of multiplex ligation-dependent probe amplification (MLPA) to classify the different types of NF1 microdeletion and indicate the importance of high-resolution microarray analysis for correct classification, a necessary precondition to identify those genes responsible for the NF1 microdeletion syndrome.

Splicing Outcomes of 5' Splice Site GT>GC Variants That Generate Wild-Type Transcripts Differ Significantly Between Full-Length and Minigene Splicing Assays

Combining data derived from a meta-analysis of human disease-associated 5' splice site GT>GC (i.e., +2T>C) variants and a cell culture-based full-length gene splicing assay (FLGSA) of forward engineered +2T>C substitutions, we recently estimated that ∼15-18% of +2T>C variants can generate up to 84% wild-type transcripts relative to their wild-type counterparts. Herein, we analyzed the splicing outcomes of 20 +2T>C variants that generate some wild-type transcripts in two minigene assays. We found a high discordance rate in terms of the generation of wild-type transcripts, not only between FLGSA and the minigene assays but also between the different minigene assays. In the pET01 context, all 20 wild-type minigene constructs generated the expected wild-type transcripts; of the 20 corresponding variant minigene constructs, 14 (70%) generated wild-type transcripts. In the pSPL3 context, only 18 of the 20 wild-type minigene constructs generated the expected wild-type transcripts whereas 8 of the 18 (44%) corresponding variant minigene constructs generated wild-type transcripts. Thus, in the context of a particular type of variant, we raise awareness of the limitations of minigene splicing assays and emphasize the importance of sequence context in regulating splicing. Whether or not our findings apply to other types of splice-altering variant remains to be investigated.

MutationTaster2021

Here we present an update to MutationTaster, our DNA variant effect prediction tool. The new version uses a different prediction model and attains higher accuracy than its predecessor, especially for rare benign variants. In addition, we have integrated many sources of data that only became available after the last release (such as gnomAD and ExAC pLI scores) and changed the splice site prediction model. To more easily assess the relevance of detected known disease mutations to the clinical phenotype of the patient, MutationTaster now provides information on the diseases they cause. Further changes represent a major overhaul of the interfaces to increase user-friendliness whilst many changes under the hood have been designed to accelerate the processing of uploaded VCF files. We also offer an API for the rapid automated query of smaller numbers of variants from within other software. MutationTaster2021 integrates our disease mutation search engine, MutationDistiller, to prioritise variants from VCF files using the patient's clinical phenotype. The novel version is available at https://www.genecascade.org/MutationTaster2021/. This website is free and open to all users and there is no login requirement.

Compensatory epistasis explored by molecular dynamics simulations

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position.

AMELIE speeds Mendelian diagnosis by matching patient phenotype and genotype to primary literature

The diagnosis of Mendelian disorders requires labor-intensive literature research. Trained clinicians can spend hours looking for the right publication(s) supporting a single gene that best explains a patient's disease. AMELIE (Automatic Mendelian Literature Evaluation) greatly accelerates this process. AMELIE parses all 29 million PubMed abstracts and downloads and further parses hundreds of thousands of full-text articles in search of information supporting the causality and associated phenotypes of most published genetic variants. AMELIE then prioritizes patient candidate variants for their likelihood of explaining any patient's given set of phenotypes. Diagnosis of singleton patients (without relatives' exomes) is the most time-consuming scenario, and AMELIE ranked the causative gene at the very top for 66% of 215 diagnosed singleton Mendelian patients from the Deciphering Developmental Disorders project. Evaluating only the top 11 AMELIE-scored genes of 127 (median) candidate genes per patient resulted in a rapid diagnosis in more than 90% of cases. AMELIE-based evaluation of all cases was 3 to 19 times more efficient than hand-curated database-based approaches. We replicated these results on a retrospective cohort of clinical cases from Stanford Children's Health and the Manton Center for Orphan Disease Research. An analysis web portal with our most recent update, programmatic interface, and code is available at AMELIE.stanford.edu.

Common polymorphic OTC variants can act as genetic modifiers of enzymatic activity

Understanding the role of common polymorphisms in modulating the clinical phenotype when they co‐occur with a disease‐causing lesion is of critical importance in medical genetics. We explored the impact of apparently neutral common polymorphisms, using the gene encoding the urea cycle enzyme, ornithine transcarbamylase (OTC), as a model system. Distinct combinations of genetic backgrounds embracing two missense polymorphisms were created in cis with the pathogenic p.Arg40His replacement. In vitro enzymatic assays revealed that the polymorphic variants were able to modulate OTC activity both in the presence or absence of the pathogenic lesion. First, we found that the combination of the minor alleles of polymorphisms p.Lys46Arg and p.Gln270Arg significantly enhanced enzymatic activity in the wild‐type protein. Second, enzymatic assays revealed that the minor allele of the p.Gln270Arg polymorphism was capable of ameliorating OTC activity when combined in cis with the pathogenic p.Arg40His replacement. Structural analysis predicted that the minor allele of the p.Gln270Arg polymorphism would serve to stabilize the OTC wild‐type protein, thereby corroborating the results of the experimental assays. Our findings demonstrate the potential importance of cis‐interactions between common polymorphic variants and pathogenic missense mutations and illustrate how standing genetic variation can modulate protein function.

DNA Methylation, Deamination, and Translesion Synthesis Combine to Generate Footprint Mutations in Cancer Driver Genes in B-Cell Derived Lymphomas and Other Cancers

Cancer genomes harbor numerous genomic alterations and many cancers accumulate thousands of nucleotide sequence variations. A prominent fraction of these mutations arises as a consequence of the off-target activity of DNA/RNA editing cytosine deaminases followed by the replication/repair of edited sites by DNA polymerases (pol), as deduced from the analysis of the DNA sequence context of mutations in different tumor tissues. We have used the weight matrix (sequence profile) approach to analyze mutagenesis due to Activation Induced Deaminase (AID) and two error-prone DNA polymerases. Control experiments using shuffled weight matrices and somatic mutations in immunoglobulin genes confirmed the power of this method. Analysis of somatic mutations in various cancers suggested that AID and DNA polymerases η and θ contribute to mutagenesis in contexts that almost universally correlate with the context of mutations in A:T and G:C sites during the affinity maturation of immunoglobulin genes. Previously, we demonstrated that AID contributes to mutagenesis in (de)methylated genomic DNA in various cancers. Our current analysis of methylation data from malignant lymphomas suggests that driver genes are subject to different (de)methylation processes than non-driver genes and, in addition to AID, the activity of pols η and θ contributes to the establishment of methylation-dependent mutation profiles. This may reflect the functional importance of interplay between mutagenesis in cancer and (de)methylation processes in different groups of genes. The resulting changes in CpG methylation levels and chromatin modifications are likely to cause changes in the expression levels of driver genes that may affect cancer initiation and/or progression.

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.