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.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

High Clinical and Genetic Similarity Between Chronic Pancreatitis Associated With Light-to-Moderate Alcohol Consumption and Classical Alcoholic Chronic Pancreatitis

Background and Aims

Heavy alcohol consumption and genetic factors represent the 2 major etiologies of chronic pancreatitis (CP). However, little is so far known about the clinical features and genetic basis of light-to-moderate alcohol consumption-related CP (LMA-CP).

Methods

A cross-sectional analysis was performed on 1061 Chinese CP patients between 2010 and 2015. CP was classified as classical alcoholic CP (ACP; n = 206), LMA-CP (n = 154), and idiopathic CP (ICP; n = 701). Clinical features and genetic characteristics (PRSS1, SPINK1, CTRC, CFTR variant status) were compared between the different groups. Odds ratios (ORs) with 95% confidence intervals were calculated to ascertain the combinatorial effect of alcohol consumption and gene mutation.

Results

Compared with ICP, the clinical features of LMA-CP were characterized by higher rates of developing pancreatic stones, pseudocyst, diabetes, and steatorrhea, which were similar to those associated with ACP. The prevalence of CP-related gene variants in LMA-CP was 38.3%, similar to ACP (39.8%), although significantly lower than ICP (56.2%). Alcohol consumption enhanced the risk of a poor clinical outcome, whereas genetic factors amplified alcohol's effects. Compared with ICP, LMA-CP and ACP were associated with a high risk of pancreatic stones (patients without variants, OR = 2.01 and 2.54; patients with variants, OR = 2.17 and 1.07), pseudocyst (patients without variants, OR = 1.03 and 1.43; patients with variants, OR = 1.67 and 2.14), diabetes mellitus (patients without variants, OR = 0.86 and 1.31; patients with variants, OR = 2.05 and 1.55), and steatorrhea (patients without variants, OR = 1.56 and 2.10; patients with variants, OR = 2.11 and 1.60).

Conclusion

Evidence was presented to show that LMA-CP was clinically and genetically similar to ACP but significantly different from ICP. Our findings provide support to the growing view that there is no safe level of alcohol consumption.

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.