Any given human individual carries multiple genetic variants that disrupt protein-coding genes, through structural variation, as well as nucleotide variants and indels. Predicting the phenotypic consequences of a gene disruption remains a significant challenge. Current approaches employ information from a range of biological networks to predict which human genes are haploinsufficient (meaning two copies are required for normal function) or essential (meaning at least one copy is required for viability). Using recently available study gene sets, we show that these approaches are strongly biased towards providing accurate predictions for well-studied genes. By contrast, we derive a haploinsufficiency score from a combination of unbiased large-scale high-throughput datasets, including gene co-expression and genetic variation in over 6000 human exomes. Our approach provides a haploinsufficiency prediction for over twice as many genes currently unassociated with papers listed in Pubmed as three commonly-used approaches, and outperforms these approaches for predicting haploinsufficiency for less-studied genes. We also show that fine-tuning the predictor on a set of well-studied 'gold standard' haploinsufficient genes does not improve the prediction for less-studied genes. This new score can readily be used to prioritize gene disruptions resulting from any genetic variant, including copy number variants, indels and single-nucleotide variants.

A variety of single-gene human diseases are caused by haploinsufficiency, a genetic condition by which mutational inactivation of one allele leads to reduced protein levels and functional impairment. Translational enhancement of the spare allele could exert a therapeutic effect. Here we developed BOOST, a novel gene-editing approach to rescue haploinsufficiency loci by the change of specific single nucleotides in the Kozak sequence, which controls translation by regulating start codon recognition. We evaluated for translational strength 230 Kozak sequences of annotated human haploinsufficient genes and 4621 derived variants, which can be installed by base editing, by a high-throughput reporter assay. Of these variants, 149 increased the translation of 47 Kozak sequences, demonstrating that a substantial proportion of haploinsufficient genes are controlled by suboptimal Kozak sequences. Validation of 18 variants for 8 genes produced an average enhancement in an expression window compatible with the rescue of the genetic imbalance. Base editing of the NCF1 gene, whose monoallelic loss causes chronic granulomatous disease, resulted in the desired increase of NCF1 (p47phox) protein levels in a relevant cell model. We propose BOOST as a fine-tuned approach to modulate translation, applicable to the correction of dozens of haploinsufficient monogenic disorders independently of the causing mutation.

Oxidative DNA damage plays a role in disease development and the aging process. A prominent participant in orchestrating the repair of oxidative DNA damage, particularly single-strand breaks, is the scaffold protein XRCC1. A series of chronological and biological aging parameters in XRCC1 heterozygous (HZ) mice were examined. HZ and wild-type (WT) C57BL/6 mice exhibit a similar median lifespan of ~26 months and a nearly identical maximal life expectancy of ~37 months. However, a number of HZ animals (7 of 92) showed a propensity for abdominal organ rupture, which may stem from developmental abnormalities given the prominent role of XRCC1 in endoderm and mesoderm formation. For other end-points evaluated-weight, fat composition, blood chemistries, condition of major organs, tissues and relevant cell types, behavior, brain volume and function, and chromosome and telomere integrity-HZ mice exhibited by-and-large a normal phenotype. Treatment of animals with the alkylating agent azoxymethane resulted in both liver toxicity and an increased incidence of precancerous lesions in the colon of HZ mice. Our study indicates that XRCC1 haploinsufficiency in mammals has little effect on chronological longevity and many key biological markers of aging in the absence of environmental challenges, but may adversely affect normal animal development or increase disease susceptibility to a relevant genotoxic exposure.

The Mus81-Eme1 complex is a structure-specific endonuclease that preferentially cleaves nicked Holliday junctions, 3'-flap structures and aberrant replication fork structures. Mus81-/- mice have been shown to exhibit spontaneous chromosomal aberrations and, in one of two models, a predisposition to cancers. The molecular mechanisms underlying its role in chromosome integrity, however, are largely unknown. To clarify the role of Mus81 in human cells, we deleted the gene in the human colon cancer cell line HCT116 by gene targeting. Here we demonstrate that Mus81 confers resistance to DNA crosslinking agents and slight resistance to other DNA-damaging agents. Mus81 deficiency spontaneously promotes chromosome damage such as breaks and activates the intra-S-phase checkpoint through the ATM-Chk1/Chk2 pathways. Furthermore, Mus81 deficiency activates the G2/M checkpoint through the ATM-Chk2 pathway and promotes DNA rereplication. Increased rereplication is reversed by the ectopic expression of Cdk1. Haploinsufficiency of Mus81 or Eme1 also causes similar phenotypes. These findings suggest that a complex network of the checkpoint pathways that respond to DNA double-strand breaks may participate in some of the phenotypes associated with Mus81 or Eme1 deficiency.

Telomerase, the essential enzyme that maintains telomere length, contains two core components, TERT and TR. Early studies in yeast and mouse showed that loss of telomerase leads to phenotypes only after several generations, due to telomere shortening. However, recent studies have suggested additional roles for telomerase components in transcription and the response to DNA damage. To examine these potential telomere length maintenance-independent roles of telomerase components, we examined first generation mTR(-/-) and mTERT(-/-) mice with long telomeres. We used gene expression profiling and found no genes that were differentially expressed in mTR(-/-) G1 mice and mTERT(-/-) G1 mice compared with wild-type mice. We also compared the response to DNA damage in mTR(-/-)G1 and mTERT(-/-) G1 mouse embryonic fibroblasts, and found no increase in the response to DNA damage in the absence of either telomerase component compared to wild-type. We conclude that, under physiologic conditions, neither mTR nor mTERT acts as a transcription factor or plays a role in the DNA damage response.

Much of the human genetics variant repertoire is composed of single nucleotide variants (SNV) and small insertion/deletions (indel) but structural variants (SV) remain a major part of our modified DNA. SV detection has often been a complex question to answer either because of the necessity to use different technologies (array CGH, SNP array, Karyotype, Optical Genome Mapping…) to detect each category of SV or to get an appropriate resolution (Whole Genome Sequencing). Thanks to the deluge of pangenomic analysis, Human geneticists are accumulating SV and their interpretation remains time consuming and challenging. The AnnotSV webserver (https://www.lbgi.fr/AnnotSV/) aims at being an efficient tool to (i) annotate and interpret SV potential pathogenicity in the context of human diseases, (ii) recognize potential false positive variants from all the SV identified and (iii) visualize the patient variants repertoire. The most recent developments in the AnnotSV webserver are: (i) updated annotations sources and ranking, (ii) three novel output formats to allow diverse utilization (analysis, pipelines), as well as (iii) two novel user interfaces including an interactive circos view.

A tremendous amount of DNA sequencing data is being produced around the world with the ambition to capture in more detail the mechanisms underlying human diseases. While numerous bioinformatics tools exist that allow the discovery of causal variants in Mendelian diseases, little to no support is provided to do the same for variant combinations, an essential task for the discovery of the causes of oligogenic diseases. ORVAL (the Oligogenic Resource for Variant AnaLysis), which is presented here, provides an answer to this problem by focusing on generating networks of candidate pathogenic variant combinations in gene pairs, as opposed to isolated variants in unique genes. This online platform integrates innovative machine learning methods for combinatorial variant pathogenicity prediction with visualization techniques, offering several interactive and exploratory tools, such as pathogenic gene and protein interaction networks, a ranking of pathogenic gene pairs, as well as visual mappings of the cellular location and pathway information. ORVAL is the first web-based exploration platform dedicated to identifying networks of candidate pathogenic variant combinations with the sole ambition to help in uncovering oligogenic causes for patients that cannot rely on the classical disease analysis tools. ORVAL is available at https://orval.ibsquare.be.

With the dramatic increase of pangenomic analysis, Human geneticists have generated large amount of genomic data including millions of small variants (SNV/indel) but also thousands of structural variations (SV) mainly from next-generation sequencing and array-based techniques. While the identification of the complete SV repertoire of a patient is getting possible, the interpretation of each SV remains challenging. To help identifying human pathogenic SV, we have developed a web server dedicated to their annotation and ranking (AnnotSV) as well as their visualization and interpretation (knotAnnotSV) freely available at the following address: https://www.lbgi.fr/AnnotSV/. A large amount of annotations from >20 sources is integrated in our web server including among others genes, haploinsufficiency, triplosensitivity, regulatory elements, known pathogenic or benign genomic regions, phenotypic data. An ACMG/ClinGen compliant prioritization module allows the scoring and the ranking of SV into 5 SV classes from pathogenic to benign. Finally, the visualization interface displays the annotated SV in an interactive way including popups, search fields, filtering options, advanced colouring to highlight pathogenic SV and hyperlinks to the UCSC genome browser or other public databases. This web server is designed for diagnostic and research analysis by providing important resources to the user.

In recent years great progress has been made in identification of structural variants (SV) in the human genome. However, the interpretation of SVs, especially located in non-coding DNA, remains challenging. One of the reasons stems in the lack of tools exclusively designed for clinical SVs evaluation acknowledging the 3D chromatin architecture. Therefore, we present TADeus2 a web server dedicated for a quick investigation of chromatin conformation changes, providing a visual framework for the interpretation of SVs affecting topologically associating domains (TADs). This tool provides a convenient visual inspection of SVs, both in a continuous genome view as well as from a rearrangement's breakpoint perspective. Additionally, TADeus2 allows the user to assess the influence of analyzed SVs within flaking coding/non-coding regions based on the Hi-C matrix. Importantly, the SVs pathogenicity is quantified and ranked using TADA, ClassifyCNV tools and sampling-based P-value. TADeus2 is publicly available at https://tadeus2.mimuw.edu.pl.

Mutations in the gene for telomerase reverse transcriptase (hTERT) are associated with diseases including dyskeratosis congenita, aplastic anemia, pulmonary fibrosis and cancer. Understanding the molecular basis of these telomerase-associated diseases requires dependable quantitative measurements of telomerase enzyme activity. Furthermore, recent findings that the human POT1-TPP1 chromosome end-binding protein complex stimulates telomerase activity and processivity provide incentive for testing variant telomerases in the presence of these factors. In the present work, we compare multiple disease-associated hTERT variants reconstituted with the RNA subunit hTR in two systems (rabbit reticulocyte lysates and human cell lines) with respect to telomerase enzymatic activity, processivity and activation by telomere proteins. Surprisingly, many of the previously reported disease-associated hTERT alleles give near-normal telomerase enzyme activity. It is possible that a small deficit in telomerase activity is sufficient to cause telomere shortening over many years. Alternatively, mutations may perturb functions such as the recruitment of telomerase to telomeres, which are essential in vivo but not revealed by simple enzyme assays.

Haploinsufficiency drives Darwinian evolution. Siblings, while alike in many aspects, differ due to monoallelic differences inherited from each parent. In cancer, solid tumors exhibit aneuploid genetics resulting in hundreds to thousands of monoallelic gene-level copy-number alterations (CNAs) in each tumor. Aneuploidy patterns are heterogeneous, posing a challenge to identify drivers in this high-noise genetic environment. Here, we developed Shifted Weighted Annotation Network (SWAN) analysis to assess biology impacted by cumulative monoallelic changes. SWAN enables an integrated pathway-network analysis of CNAs, RNA expression, and mutations via a simple web platform. SWAN is optimized to best prioritize known and novel tumor suppressors and oncogenes, thereby identifying drivers and potential druggable vulnerabilities within cancer CNAs. Protein homeostasis, phospholipid dephosphorylation, and ion transport pathways are commonly suppressed. An atlas of CNA pathways altered in each cancer type is released. These CNA network shifts highlight new, attractive targets to exploit in solid tumors.

Modifications of ribosomal RNA expand the nucleotide repertoire and thereby contribute to ribosome heterogeneity and translational regulation of gene expression. One particular m5C modification of 25S ribosomal RNA, which is introduced by Rcm1p, was previously shown to modulate stress responses and lifespan in yeast and other small organisms. Here, we report that NSUN5 is the functional orthologue of Rcm1p, introducing m5C3782 into human and m5C3438 into mouse 28S ribosomal RNA. Haploinsufficiency of the NSUN5 gene in fibroblasts from William Beuren syndrome patients causes partial loss of this modification. The N-terminal domain of NSUN5 is required for targeting to nucleoli, while two evolutionary highly conserved cysteines mediate catalysis. Phenotypic consequences of NSUN5 deficiency in mammalian cells include decreased proliferation and size, which can be attributed to a reduction in total protein synthesis by altered ribosomes. Strikingly, Nsun5 knockout in mice causes decreased body weight and lean mass without alterations in food intake, as well as a trend towards reduced protein synthesis in several tissues. Together, our findings emphasize the importance of single RNA modifications for ribosome function and normal cellular and organismal physiology.

Kleefstra syndrome, a disease with intellectual disability, autism spectrum disorders and other developmental defects is caused in humans by haploinsufficiency of EHMT1. Although EHMT1 and its paralog EHMT2 were shown to be histone methyltransferases responsible for deposition of the di-methylated H3K9 (H3K9me2), the exact nature of epigenetic dysfunctions in Kleefstra syndrome remains unknown. Here, we found that the epigenome of Ehmt1+/- adult mouse brain displays a marked increase of H3K9me2/3 which correlates with impaired expression of protocadherins, master regulators of neuronal diversity. Increased H3K9me3 was present already at birth, indicating that aberrant methylation patterns are established during embryogenesis. Interestingly, we found that Ehmt2+/- mice do not present neither the marked increase of H3K9me2/3 nor the cognitive deficits found in Ehmt1+/- mice, indicating an evolutionary diversification of functions. Our finding of increased H3K9me3 in Ehmt1+/- mice is the first one supporting the notion that EHMT1 can quench the deposition of tri-methylation by other Histone methyltransferases, ultimately leading to impaired neurocognitive functioning. Our insights into the epigenetic pathophysiology of Kleefstra syndrome may offer guidance for future developments of therapeutic strategies for this disease.

Haploinsufficiency of EFTUD2 (Elongation Factor Tu GTP Binding Domain Containing 2) is linked to human mandibulofacial dysostosis, Guion-Almeida type (MFDGA), but the underlying cellular and molecular mechanisms remain to be addressed. We report here the isolation, cloning and functional analysis of the mutated eftud2 (snu114) in a novel neuronal mutant fn10a in zebrafish. This mutant displayed abnormal brain development with evident neuronal apoptosis while the development of other organs appeared less affected. Positional cloning revealed a nonsense mutation such that the mutant eftud2 mRNA encoded a truncated Eftud2 protein and was subjected to nonsense-mediated decay. Disruption of eftud2 led to increased apoptosis and mitosis of neural progenitors while it had little effect on differentiated neurons. Further RNA-seq and functional analyses revealed a transcriptome-wide RNA splicing deficiency and a large amount of intron-retaining and exon-skipping transcripts, which resulted in inadequate nonsense-mediated RNA decay and activation of the p53 pathway in fn10a mutants. Therefore, our study has established that eftud2 functions in RNA splicing during neural development and provides a suitable zebrafish model for studying the molecular pathology of the neurological disease MFDGA.

SOX9 encodes a transcription factor that presides over the specification and differentiation of numerous progenitor and differentiated cell types, and although SOX9 haploinsufficiency and overexpression cause severe diseases in humans, including campomelic dysplasia, sex reversal and cancer, the mechanisms underlying SOX9 transcription remain largely unsolved. We identify here an evolutionarily conserved enhancer located 70-kb upstream of mouse Sox9 and call it SOM because it specifically activates a Sox9 promoter reporter in most Sox9-expressing somatic tissues in transgenic mice. Moreover, SOM-null fetuses and pups reduce Sox9 expression by 18-37% in the pancreas, lung, kidney, salivary gland, gut and liver. Weanlings exhibit half-size pancreatic islets and underproduce insulin and glucagon, and adults slowly recover from acute pancreatitis due to a 2-fold impairment in Sox9 upregulation. Molecular and genetic experiments reveal that Sox9 protein dimers bind to multiple recognition sites in the SOM sequence and are thereby both necessary and sufficient for enhancer activity. These findings thus uncover that Sox9 directly enhances its functions in somatic tissue development and adult regeneration through SOM-mediated positive auto-regulation. They provide thereby novel insights on molecular mechanisms controlling developmental and disease processes and suggest new strategies to improve disease treatments.

BORIS, like other members of the 'cancer/testis antigen' family, is normally expressed in testicular germ cells and repressed in somatic cells, but is aberrantly activated in cancers. To understand regulatory mechanisms governing human BORIS expression, we characterized its 5'-flanking region. Using 5' RACE, we identified three promoters, designated A, B and C, corresponding to transcription start sites at -1447, -899 and -658 bp upstream of the first ATG. Alternative promoter usage generated at least five alternatively spliced BORIS mRNAs with different half-lives determined by varying 5'-UTRs. In normal testis, BORIS is transcribed from all three promoters, but 84% of the 30 cancer cell lines tested used only promoter(s) A and/or C while the others utilized primarily promoters B and C. The differences in promoter usage between normal and cancer cells suggested that they were subject to differential regulation. We found that DNA methylation and functional p53 contributes to the negative regulation of each promoter. Moreover, reduction of CTCF in normally BORIS-negative human fibroblasts resulted in derepression of BORIS promoters. These results provide a mechanistic basis for understanding cancer-related associations between haploinsufficiency of CTCF and BORIS derepression, and between the lack of functional p53 and aberrant activation of BORIS.

One third of tumor suppressors are haploinsufficient transcriptional regulators, yet it remains unknown how a 50% reduction of a transcription factor is translated at the cis-regulatory level into a malignant transcriptional program. We studied CUX1, a haploinsufficient transcription factor that is recurrently mutated in hematopoietic and solid tumors. We determined CUX1 DNA-binding and target gene regulation in the wildtype and haploinsufficient states. CUX1 binds with transcriptional activators and cohesin at distal enhancers across three different human cell types. Haploinsufficiency of CUX1 altered the expression of a large number of genes, including cell cycle regulators, with concomitant increased cellular proliferation. Surprisingly, CUX1 occupancy decreased genome-wide in the haploinsufficient state, and binding site affinity did not correlate with differential gene expression. Instead, differentially expressed genes had multiple, low-affinity CUX1 binding sites, features of analog gene regulation. A machine-learning algorithm determined that chromatin accessibility, enhancer activity, and distance to the transcription start site are features of dose-sensitive CUX1 transcriptional regulation. Moreover, CUX1 is enriched at sites of DNA looping, as determined by Hi-C analysis, and these loops connect CUX1 to the promoters of regulated genes. We propose an analog model for haploinsufficient transcriptional deregulation mediated by higher order genome architecture.

DIDA (DIgenic diseases DAtabase) is a novel database that provides for the first time detailed information on genes and associated genetic variants involved in digenic diseases, the simplest form of oligogenic inheritance. The database is accessible via http://dida.ibsquare.be and currently includes 213 digenic combinations involved in 44 different digenic diseases. These combinations are composed of 364 distinct variants, which are distributed over 136 distinct genes. The web interface provides browsing and search functionalities, as well as documentation and help pages, general database statistics and references to the original publications from which the data have been collected. The possibility to submit novel digenic data to DIDA is also provided. Creating this new repository was essential as current databases do not allow one to retrieve detailed records regarding digenic combinations. Genes, variants, diseases and digenic combinations in DIDA are annotated with manually curated information and information mined from other online resources. Next to providing a unique resource for the development of new analysis methods, DIDA gives clinical and molecular geneticists a tool to find the most comprehensive information on the digenic nature of their diseases of interest.