ARID1B is the most frequently mutated gene in Coffin-Siris syndrome (CSS). To date, the vast majority of causative variants reported in ARID1B are truncating, leading to nonsense-mediated mRNA decay. In the absence of experimental data, only few ARID1B amino acid substitutions have been classified as pathogenic, mainly based on clinical data and their de novo occurrence, while most others are currently interpreted as variants of unknown significance. The present study substantiates the pathogenesis of ARID1B non-truncating/NMD-escaping variants located in the SMARCA4-interacting EHD2 and DNA-binding ARID domains. Overexpression assays in cell lines revealed that the majority of EHD2 variants lead to protein misfolding and formation of cytoplasmic aggresomes surrounded by vimentin cage-like structures and co-localizing with the microtubule organisation center. ARID domain variants exhibited not only aggresomes, but also nuclear aggregates, demonstrating robust pathological effects. Protein levels were not compromised, as shown by quantitative western blot analysis. In silico structural analysis predicted the exposure of amylogenic segments in both domains due to the nearby variants, likely causing this aggregation. Genome-wide transcriptome and methylation analysis in affected individuals revealed expression and methylome patterns consistent with those of the pathogenic haploinsufficiency ARID1B alterations in CSS cases. These results further support pathogenicity and indicate two approaches for disambiguation of such variants in everyday practice. The few affected individuals harbouring EHD2 non-truncating variants described to date exhibit mild CSS clinical traits. In summary, this study paves the way for the re-evaluation of previously unclear ARID1B non-truncating variants and opens a new era in CSS genetic diagnosis.

An elevated resting heart rate (RHR) is associated with increased cardiovascular mortality. Genome-wide association studies (GWAS) have identified > 350 loci. Uniquely, in this study we applied genetic fine-mapping leveraging tissue specific chromatin segmentation and colocalization analyses to identify causal variants and candidate effector genes for RHR. We used RHR GWAS summary statistics from 388,237 individuals of European ancestry from UK Biobank and performed fine mapping using publicly available genomic annotation datasets. High-confidence causal variants (accounting for > 75% posterior probability) were identified, and we collated candidate effector genes using a multi-omics approach that combined evidence from colocalisation with molecular quantitative trait loci (QTLs), and long-range chromatin interaction analyses. Finally, we performed druggability analyses to investigate drug repurposing opportunities. The fine mapping pipeline indicated 442 distinct RHR signals. For 90 signals, a single variant was identified as a high-confidence causal variant, of which 22 were annotated as missense. In trait-relevant tissues, 39 signals colocalised with cis-expression QTLs (eQTLs), 3 with cis-protein QTLs (pQTLs), and 75 had promoter interactions via Hi-C. In total, 262 candidate genes were highlighted (79% had promoter interactions, 15% had a colocalised eQTL, 8% had a missense variant and 1% had a colocalised pQTL), and, for the first time, enrichment in nervous system pathways. Druggability analyses highlighted ACHE, CALCRL, MYT1 and TDP1 as potential targets. Our genetic fine-mapping pipeline prioritised 262 candidate genes for RHR that warrant further investigation in functional studies, and we provide potential therapeutic targets to reduce RHR and cardiovascular mortality.

In recent years, there has been increased focus on exploring the role the non-protein-coding genome plays in Mendelian disorders. One class of particular interest is long non-coding RNAs (lncRNAs), which has recently been implicated in the regulation of diverse molecular processes. However, because lncRNAs do not encode protein, there is uncertainty regarding what constitutes a pathogenic lncRNA variant, and thus annotating such elements is challenging. The Developmental Genome Anatomy Project (DGAP) and similar projects recruit individuals with apparently balanced chromosomal abnormalities (BCAs) that disrupt or dysregulate genes in order to annotate the human genome. We hypothesized that rearrangements disrupting lncRNAs could be the underlying genetic etiology for the phenotypes of a subset of these individuals. Thus, we assessed 279 cases with BCAs and selected 191 cases with simple BCAs (breakpoints at only two genomic locations) for further analysis of lncRNA disruptions. From these, we identified 66 cases in which the chromosomal rearrangements directly disrupt lncRNAs. In 30 cases, no genes of any other class aside from lncRNAs are directly disrupted, consistent with the hypothesis that lncRNA disruptions could underly the phenotypes of these individuals. Strikingly, the lncRNAs MEF2C-AS1 and ENSG00000257522 are each disrupted in two unrelated cases. Furthermore, we experimentally tested the lncRNAs TBX2-AS1 and MEF2C-AS1 and found that knockdown of these lncRNAs resulted in decreased expression of the neighboring transcription factors TBX2 and MEF2C, respectively. To showcase the power of this genomic approach for annotating lncRNAs, here we focus on clinical reports and genetic analysis of seven individuals with likely developmental etiologies due to lncRNA disruptions.

Long non-coding RNA (lncRNA) genes represent a large class of transcripts that are widely expressed across species. As most human lncRNAs are non-conserved, we recently employed a unique humanized liver mouse model to study lncRNAs expressed in human livers. We identified a human hepatocyte-specific lncRNA, hLMR1 (human lncRNA metabolic regulator 1), which is induced by feeding and promotes hepatic cholesterol synthesis. Recent genome-wide association studies (GWAS) found that several single-nucleotide polymorphisms (SNPs) from the hLMR1 gene locus are associated with blood lipids and markers of liver damage. These results suggest that dietary and genetic factors may regulate hLMR1 to affect disease progression. In this study, we first screened for nutritional/hormonal factors and found that hLMR1 was robustly induced by insulin/glucose in cultured human hepatocytes, and this induction is dependent on the transcription factor SREBP1. We then tested if GWAS SNPs genetically linked to hLMR1 could regulate hLMR1 expression. We found that DNA sequences flanking rs9653945, a SNP from the last exon of the hLMR1 gene, functions as an enhancer that can be robustly activated by SREBP1c depending on the presence of rs9653945 major allele (G). We further performed CRISPR base editing in human HepG2 cells and found that rs9653945 major (G) to minor (A) allele modification resulted in blunted insulin/glucose-induced expression of hLMR1. Finally, we performed genotyping and gene expression analyses using a published human NAFLD RNA-seq dataset and found that individuals homozygous for rs9653945-G have a higher expression of hLMR1 and risk of NAFLD. Taken together, our data support a model that rs9653945-G predisposes individuals to insulin/glucose-induced hLMR1, contributing to the development of hyperlipidemia and NAFLD.

Biallelic pathogenic variants in MAP3K20, which encodes a mitogen-activated protein kinase, are a rare cause of split-hand foot malformation (SHFM), hearing loss, and nail abnormalities or congenital myopathy. However, heterozygous variants in this gene have not been definitively associated with a phenotype. Here, we describe the phenotypic spectrum associated with heterozygous de novo variants in the linker region between the kinase domain and leucine zipper domain of MAP3K20. We report five individuals with diverse clinical features, including craniosynostosis, limb anomalies, sensorineural hearing loss, and ectodermal dysplasia-like phenotypes who have heterozygous de novo variants in this specific region of the gene. These individuals exhibit both shared and unique clinical manifestations, highlighting the complexity and variability of the disorder. We propose that the involvement of MAP3K20 in endothelial-mesenchymal transition provides a plausible etiology of these features. Together, these findings characterize a disorder that both expands the phenotypic spectrum associated with MAP3K20 and highlights the need for further studies on its role in early human development.

Identification of genes associated with nonsyndromic hearing loss is a crucial endeavor given the substantial number of individuals who remain without a diagnosis after even the most advanced genetic testing. PKHD1L1 was established as necessary for the formation of the cochlear hair-cell stereociliary coat and causes hearing loss in mice and zebrafish when mutated. We sought to determine if biallelic variants in PKHD1L1 also cause hearing loss in humans. Exome sequencing was performed on DNA of four families segregating autosomal recessive nonsyndromic sensorineural hearing loss. Compound heterozygous p.[(Gly129Ser)];p.[(Gly1314Val)] and p.[(Gly605Arg)];p[(Leu2818TyrfsTer5)], homozygous missense p.(His2479Gln) and nonsense p.(Arg3381Ter) variants were identified in PKHD1L1 that were predicted to be damaging using in silico pathogenicity prediction methods. In vitro functional analysis of two missense variants was performed using purified recombinant PKHD1L1 protein fragments. We then evaluated protein thermodynamic stability with and without the missense variants found in one of the families and performed a minigene splicing assay for another variant. In silico molecular modeling using AlphaFold2 and protein sequence alignment analysis were carried out to further explore potential variant effects on structure. In vitro functional assessment indicated that both engineered PKHD1L1 p.(Gly129Ser) and p.(Gly1314Val) mutant constructs significantly reduced the folding and structural stabilities of the expressed protein fragments, providing further evidence to support pathogenicity of these variants. Minigene assay of the c.1813G>A p.(Gly605Arg) variant, located at the boundary of exon 17, revealed exon skipping leading to an in-frame deletion of 48 amino acids. In silico molecular modeling exposed key structural features that might suggest PKHD1L1 protein destabilization. Multiple lines of evidence collectively associate PKHD1L1 with nonsyndromic mild–moderate to severe sensorineural hearing loss. PKHD1L1 testing in individuals with mild–moderate hearing loss may identify further affected families.

Meniere disease is a complex inner ear disorder with significant familial aggregation. A differential prevalence of familial MD (FMD) has been reported, being 9–10% in Europeans compared to 6% in East Asians. A broad genetic heterogeneity in FMD has been described, OTOG being the most common mutated gene, with a compound heterozygous recessive inheritance. We hypothesize that an OTOG-related founder effect may explain the higher prevalence of FMD in the European population. Therefore, the present study aimed to compare the allele frequency (AF) and distribution of OTOG rare variants across different populations. For this purpose, the coding regions with high constraint (low density of rare variants) were retrieved in the OTOG coding sequence in Non-Finnish European (NFE).. Missense variants (AF < 0.01) were selected from a 100 FMD patient cohort, and their population AF was annotated using gnomAD v2.1. A linkage analysis was performed, and odds ratios were calculated to compare AF between NFE and other populations. Thirteen rare missense variants were observed in 13 FMD patients, with 2 variants (rs61978648 and rs61736002) shared by 5 individuals and another variant (rs117315845) shared by two individuals. The results confirm the observed enrichment of OTOG rare missense variants in FMD. Furthermore, eight variants were enriched in the NFE population, and six of them were in constrained regions. Structural modeling predicts five missense variants that could alter the otogelin stability. We conclude that several variants reported in FMD are in constraint regions, and they may have a founder effect and explain the burden of FMD in the European population.

Auditory neuropathy (AN) is a unique type of language developmental disorder, with no precise rate of genetic contribution that has been deciphered in a large cohort. In a retrospective cohort of 311 patients with AN, pathogenic and likely pathogenic variants of 23 genes were identified in 98 patients (31.5% in 311 patients), and 14 genes were mutated in two or more patients. Among subgroups of patients with AN, the prevalence of pathogenic and likely pathogenic variants was 54.4% and 56.2% in trios and families, while 22.9% in the cases with proband-only; 45.7% and 25.6% in the infant and non-infant group; and 33.7% and 0% in the bilateral and unilateral AN cases. Most of the OTOF gene (96.6%, 28/29) could only be identified in the infant group, while the AIFM1 gene could only be identified in the non-infant group; other genes such as ATP1A3 and OPA1 were identified in both infant and non-infant groups. In conclusion, genes distribution of AN, with the most common genes being OTOF and AIFM1, is totally different from other sensorineural hearing loss. The subgroups with different onset ages showed different genetic spectrums, so did bilateral and unilateral groups and sporadic and familial or trio groups.

Ocular coloboma (OC) is a congenital disorder caused by the incomplete closure of the embryonic ocular fissure. OC can present as a simple anomaly or, in more complex forms, be associated with additional ocular abnormalities. It can occur in isolation or as part of a broader syndrome, exhibiting considerable genetic heterogeneity. Diagnostic yield for OC remains below 30%, indicating the need for further genetic exploration. Mutations in the Wnt receptor FZD5, which is expressed throughout eye development, have been linked to both isolated and complex forms of coloboma. These mutations often result in a dominant-negative effect, where the mutated FZD5 protein disrupts WNT signaling by sequestering WNT ligands. Here, we describe a case of syndromic bilateral OC with additional features such as microcornea, bone developmental anomalies, and mild intellectual disability. Whole exome sequencing revealed a homozygous rare missense variant in FZD5. Consistent with a loss-of-function effect, overexpressing of fzd5 mRNA harboring the missense variant in zebrafish embryos does not influence embryonic development, whereas overexpression of wild-type fzd5 mRNA results in body axis duplications. However, in vitro TOPFlash assays revealed that the missense variant only caused partial loss-of-function, behaving as a hypomorphic mutation. We further showed that the mutant protein still localized to the cell membrane and maintained proper conformation when modeled in silico, suggesting that the impairment lies in signal transduction. This hypothesis is further supported by the fact that the variant affects a highly conserved amino acid known to be crucial for protein-protein interactions.Supplementary InformationThe online version contains supplementary material available at 10.1007/s00439-024-02712-y.

Large-scale association analyses using whole-genome sequence data have become feasible, but understanding the functional impacts of these associations remains challenging. Although many tools are available to predict the functional impacts of genetic variants, it is unclear which tool should be used in practice. This work provides a practical guide to assist in selecting appropriate tools for variant annotation. We conducted a MEDLINE search up to November 10, 2023, and included tools that are applicable to a broad range of phenotypes, can be used locally, and have been recently updated. Tools were categorized based on the types of variants they accept and the functional impacts they predict. Sequence Ontology terms were used for standardization. We identified 118 databases and software packages, encompassing 36 variant types and 161 functional impacts. Combining only three tools, namely SnpEff, FAVOR, and SparkINFERNO, allows predicting 99 (61%) distinct functional impacts. Thirty-seven tools predict 89 functional impacts that are not supported by any other tool, while 75 tools predict pathogenicity and can be used within the ACMG/AMP guidelines in a clinical context. We launched a website allowing researchers to select tools based on desired variants and impacts. In summary, more than 100 tools are already available to predict approximately 160 functional impacts. About 60% of the functional impacts can be predicted by the combination of three tools. Unexpectedly, recent tools do not predict more impacts than older ones. Future research should allow predicting the functionality of so far unsupported variant types, such as gene fusions.URL: https://cardio-care.shinyapps.io/VEP_Finder/.Registration: OSF Registries on November 10, 2023, https://osf.io/s2gct.Supplementary InformationThe online version contains supplementary material available at 10.1007/s00439-024-02670-5.