Analysis of protein-coding genetic variation in 60,706 humans

LARP1 haploinsufficiency is associated with an autosomal dominant neurodevelopmental disorder

Autism spectrum disorder (ASD) is a neurodevelopmental disorder (NDD) that affects approximately 4% of males and 1% of females in the United States. While causes of ASD are multi-factorial, single rare genetic variants contribute to around 20% of cases. Here, we report a case series of seven unrelated probands (6 males, 1 female) with ASD or another variable NDD phenotype attributed to de novo heterozygous loss of function or missense variants in the gene LARP1 (La ribonucleoprotein 1). LARP1 encodes an RNA-binding protein that post-transcriptionally regulates the stability and translation of thousands of mRNAs, including those regulating cellular metabolism and metabolic plasticity. Using lymphocytes collected and immortalized from an index proband who carries a truncating variant in one allele of LARP1, we demonstrated that lower cellular levels of LARP1 protein cause reduced rates of aerobic respiration and glycolysis. As expression of LARP1 increases during neurodevelopment, with higher levels in neurons and astrocytes, we propose that LARP1 haploinsufficiency contributes to ASD or related NDDs through attenuated metabolic activity in the developing fetal brain.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

MDVarP: modifier ~ disease-causing variant pairs predictor

BACKGROUND

Modifiers significantly impact disease phenotypes by modulating the effects of disease-causing variants, resulting in varying disease manifestations among individuals. However, identifying genetic interactions between modifier and disease-causing variants is challenging.

RESULTS

We developed MDVarP, an ensemble model comprising 1000 random forest predictors, to identify modifier ~ disease-causing variant combinations. MDVarP achieves high accuracy and precision, as verified using an independent dataset with published evidence of genetic interactions. We identified 25 novel modifier ~ disease-causing variant combinations and obtained supporting evidence for these associations. MDVarP outputs a class label ("Associated-pair" or "Nonrelevant-pair") and two prediction scores indicating the probability of a true association.

CONCLUSIONS

MDVarP prioritizes variant pairs associated with phenotypic modulations, enabling more effective mapping of functional contributions from disease-causing and modifier variants. This framework interprets genetic interactions underlying phenotypic variations in human diseases, with potential applications in personalized medicine and disease prevention.

Importance of copy number variants in childhood apraxia of speech and other speech sound disorders

Childhood apraxia of speech (CAS) is a severe and rare form of speech sound disorder (SSD). CAS is typically sporadic, but may segregate in families with broader speech and language deficits. We hypothesize that genetic changes may be involved in the etiology of CAS. We conduct whole-genome sequencing in 27 families with CAS, 101 individuals in all. We identify 17 genomic regions including 19 unique copy number variants (CNVs). Three variants are shared across families, but the rest are unique; three events are de novo. In four families, siblings with milder phenotypes co-inherited the same CNVs, demonstrating variable expressivity. We independently validate eight CNVs using microarray technology and find many of these CNVs are present in children with milder forms of SSD. Bioinformatic investigation reveal four CNVs with substantial functional consequences (cytobands 2q24.3, 6p12.3-6p12.2, 11q23.2-11q23.3, and 16p11.2). These discoveries show that CNVs are a heterogeneous, but prevalent, cause of CAS.

UBTF haploinsufficiency associated with UBTF-related global developmental delay and distinctive facial features without neuroregression

BACKGROUND

The Upstream Binding Transcription Factor (UBTF) gene encodes two nucleolar proteins, UBTF1 and UBTF2. UBTF1 regulates rRNA transcription by RNA polymerase I, while UBTF2 regulates mRNA transcription by RNA polymerase II. A recurrent de novo dominant mutation c.628G>A (p.Glu210Lys) has been identified as a gain-of-function mutation associated with childhood onset neurodegeneration with brain atrophy (CONDBA). Evidence from large-scale population databases and Ubtf+/- mouse models indicates that UBTF haploinsufficiency is not tolerated.

METHODS

Three unrelated patients with global developmental delay and distinctive facial features were recruited for the study. Whole exome sequencing (WES) was performed to identify potential genetic abnormalities. Additionally, copy number variation analysis was conducted based on the WES data.

RESULTS

All three patients exhibited intellectual disabilities, social challenges and developmental delays in language and gross motor skills. Distinctive facial features included a wide forehead, sparse eyebrows, hypertelorism, narrow palpebral fissures, single-fold eyelids, a flat nasal bridge, anteverted nares, a long philtrum and a thin upper lip. Additionally, patient C presented with more severe language delay, recurrent hepatic dysfunction and an atrial septal defect. Patient A was found to have a nonsense variant, c.1327C>T (p.R443Ter), in the exon 13 of UBTF. Patients B and C both carried a heterozygous deletion encompassing the UBTF gene.

CONCLUSION

In this study, we analysed the detailed phenotypes associated with UBTF haploinsufficiency, which, to our knowledge, have not been previously reported. We propose that UBTF haploinsufficiency-related global developmental delay and distinctive facial features, without neuroregression, constitute a new syndrome distinct from CONDBA.

A Cysteinyl-tRNA Synthetase Mutation Causes Novel Autosomal-Dominant Inheritance of a Parkinsonism/Spinocerebellar-Ataxia Complex

This study aimed to identify possible pathogenic genes in a 90-member family with a rare combination of multiple neurodegenerative disease phenotypes, which has not been depicted by the known neurodegenerative disease. We performed physical and neurological examinations with International Rating Scales to assess signs of ataxia, Parkinsonism, and cognitive function, as well as brain magnetic resonance imaging scans with seven sequences. We searched for co-segregations of abnormal repeat-expansion loci, pathogenic variants in known spinocerebellar ataxia-related genes, and novel rare mutations via whole-genome sequencing and linkage analysis. A rare co-segregating missense mutation in the CARS gene was validated by Sanger sequencing and the aminoacylation activity of mutant CARS was measured by spectrophotometric assay. This pedigree presented novel late-onset core characteristics including cerebellar ataxia, Parkinsonism, and pyramidal signs in all nine affected members. Brain magnetic resonance imaging showed cerebellar/pons atrophy, pontine-midline linear hyperintensity, decreased rCBF in the bilateral basal ganglia and cerebellar dentate nucleus, and hypo-intensities of the cerebellar dentate nuclei, basal ganglia, mesencephalic red nuclei, and substantia nigra, all of which suggested neurodegeneration. Whole-genome sequencing identified a novel pathogenic heterozygous mutation (E795V) in the CARS gene, meanwhile, exhibited none of the known repeat-expansions or point mutations in pathogenic genes. Remarkably, this CARS mutation causes a 20% decrease in aminoacylation activity to charge tRNACys with L-cysteine in protein synthesis compared with that of the wild type. All family members carrying a heterozygous mutation CARS (E795V) had the same clinical manifestations and neuropathological changes of Parkinsonism and spinocerebellar-ataxia. These findings identify novel pathogenesis of Parkinsonism-spinocerebellar ataxia and provide insights into its genetic architecture.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing, the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare diseases. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery, which should, in turn, increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks such as Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, and researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Molecular genetic diagnostics for inherited retinal dystrophies in the clinical setting

OBJECTIVE

To evaluate the success of diagnostic genetic testing in inherited retinal dystrophy (IRD) patients in the clinical setting.

DESIGN

Retrospective cohort analysis.

PARTICIPANTS

A total of 446 consecutive participants from diverse ethnic backgrounds living in western Canada.

METHODS

Clinical information was collected from participants, including family history, and they underwent a full ophthalmic examination with chart review. Those with a suspected IRD were offered panel-based genetic testing of 351 genes between March 1, 2019, and February 28, 2022. The main outcome measure was effect of the genetic testing results on clinical diagnosis.

RESULTS

Genetic testing established a conclusive molecular diagnosis in 249 of 446 cases (55.8%), a clearly negative result in 90 of 446 cases (20.1%), and an inconclusive diagnosis in 108 of 446 cases (24.2%). Conclusive disease-causing variants were identified in 69 genes, and the most commonly affected genes were ABCA4 (31 variants), USH2A (25 variants), and RPGR (19 variants). The inconclusive group included likely novel autosomal dominant variants or a pathogenic variant with a variant of uncertain significance in the same gene for a recessive phenotype. Notably, an inconclusive molecular genetic diagnosis was seen in as many as 47.3% of East Asian participants with an outer retinal dystrophy.

CONCLUSIONS

This study represents the largest review of molecular genetic testing in IRDs in Canada. That negative or inconclusive results obtained in approximately 45% of cases demonstrates that there is an important need for new research into molecular genetic causes of IRDs. This is particularly true in addressing the problem of interpreting a variant of uncertain significance in ethnic minorities.

A systematic comparison of computational methods for expression forecasting

Expression forecasting methods use machine learning models to predict how a cell will alter its transcriptome upon perturbation. Such methods are enticing because they promise to answer pressing questions in fields ranging from developmental genetics to cell fate engineering and because they are a fast, cheap, and accessible complement to the corresponding experiments. However, the absolute and relative accuracy of these methods is poorly characterized, limiting their informed use, their improvement, and the interpretation of their predictions. To address these issues, we created a benchmarking platform that combines a panel of 11 large-scale perturbation datasets with an expression forecasting software engine that encompasses or interfaces to a wide variety of methods. We used our platform to systematically assess methods, parameters, and sources of auxiliary data, finding that performance strongly depends on the choice of metric, and especially for simple metrics like mean squared error, it is uncommon for expression forecasting methods to out-perform simple baselines. Our platform will serve as a resource to improve methods and to identify contexts in which expression forecasting can succeed.

Determination of the frequency and distribution of APC, PIK3CA, and SMAD4 gene mutations in Ugandan patients with colorectal cancer

Uganda is a developing low-income country with a low incidence of colorectal cancer, which is steadily increasing. Ugandan colorectal cancer (CRC) patients are young and present with advanced-stage disease. In our population, there is a scarcity of genetic oncological studies, therefore, we investigated the mutational status of CRC tissues, focusing in particular on the adenomatous polyposis coli (APC), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA), and SMAD4 genes. Our objective was to determine whether there were any differences between other populations and Ugandan patients. We performed next-generation sequencing on the extracted DNA from formalin-fixed paraffin-embedded adenocarcinoma samples from 127 patients (mean (SD) age: 54.9 (16.0) years; male:female sex ratio: 1.2:1). Most tumours were located in the rectum 56 (44.1%), 14 (11%) tumours were high grade, and 96 (75.6%) were moderate grade CRC. Stage III + IV CRC tumours were found in 109 (85.8%) patients. We identified 48 variants of APC, including 9 novel APC mutations that were all pathogenic or deleterious. For PIK3CA, we found 19 variants, of which 9 were deleterious or pathogenic. Four PIK3CA novel pathogenic or deleterious variants were included (c.1397C > G, c.2399_2400insA, c.2621G > C, c.2632C > G). Three SMAD4 variants were reported, including two pathogenic or deleterious variants (c.1268G > T, c.556dupC) and one tolerant (c.563A > C) variant. One novel SMAD4 deleterious mutation (c.1268G > T) was reported. In conclusion, we provide clinicopathological information and new genetic variation data pertinent to CRC in Uganda.

Polygenic scores stratify neurodevelopmental copy number variant carrier cognitive outcomes in the UK Biobank

Rare copy-number variants associated with neurodevelopmental conditions (ND-CNVs) exhibit variable expressivity of clinical, physical, behavioural outcomes. Findings from clinically ascertained cohorts suggest this variability may be partly due to additional genetic variation. Here, we assessed the impact of polygenic scores (PGS) and rare variants on ND-CNV carrier fluid intelligence (FI) scores in the UK Biobank. Greater PGS for cognition (PSCog) and educational attainment (PSEA) is associated with increased FI scores in all ND-CNVs (n = 1317), 15q11.2 del. (n = 543), and 16p13.11 dup. carriers (n = 275). No association of rare variants associated with intellectual disability, autism, or putatively loss-of-function, brain-expressed genes was found. Positive predictive values in the first deciles of PScog and PSEA showed a two- to five-fold increase in the rate of low FI scores compared to baseline rates. These findings demonstrate that PGS can stratify ND-CNV carrier cognitive outcomes in a population-based cohort.

CNVoyant a machine learning framework for accurate and explainable copy number variant classification

The precise classification of copy number variants (CNVs) presents a significant challenge in genomic medicine, primarily due to the complex nature of CNVs and their diverse impact on rare genetic diseases (RGDs). This complexity is compounded by the limitations of existing methods in accurately distinguishing between benign, uncertain, and pathogenic CNVs. Addressing this gap, we introduce CNVoyant, a machine learning-based multi-class framework designed to enhance the clinical significance classification of CNVs. Trained on a comprehensive dataset of 52,176 ClinVar entries across pathogenic, uncertain, and benign classifications, CNVoyant incorporates a broad spectrum of genomic features, including genome position, disease-gene annotations, dosage sensitivity, and conservation scores. Models to predict the clinical significance of copy number gains and losses were trained independently. Final models were selected after testing 29 machine learning architectures and 10,000 hyperparameter combinations each for deletions and duplications via fivefold cross-validation. We validate the performance of CNVoyant by leveraging a comprehensive set of 21,574 CNVs from the DECIPHER database, a highly regarded resource known for its extensive catalog of chromosomal imbalances linked to clinical outcomes. Compared to alternative approaches, CNVoyant shows marked improvements in precision-recall and ROC AUC metrics for binary pathogenic classifications while going one step further, offering multi-classification of clinical significance and corresponding SHAP explainability plots. Additionally, when provided germline CNV calls from real-world RGD cases with diagnostic CNV(s), CNVoyant correctly classified all diagnostic CNVs as having pathogenic significance with high confidence. This large-scale validation demonstrates CNVoyant's superior accuracy and underscores its potential to aid genomic researchers and clinical geneticists in interpreting the clinical implications of real CNVs.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

PMS2 mutation spectra in Norway and risk of cancer for carriers of pathogenic variants

BACKGROUND

In Norway, we have offered testing of PMS2 since 2006, and have a large national cohort of carriers. The aim of this study was to describe all PMS2 variants identified, and to describe frequency, spectrum and penetrance of cancers in carriers of class 4/5 variants.

METHODS

All detected PMS2 variants were collected from the diagnostic laboratories and reclassified according to ACMG criteria and gene specific guidelines. Data on variant, gender, cancer diagnosis, age at diagnosis, and age at last known follow-up was collected on all carriers of class 4/5 variants from electronic patient records. The Kaplan-Meier algorithm was used to calculate cumulative risk of any cancer, colorectal cancer and endometrial cancer.

RESULTS

In total, 220 different PMS2 variants were detected. Twenty nine class 4/5 variants were identified in 482 carriers. The most common pathogenic variant was the founder mutation c.989-1G > T, detected in 204 patients from 58 families. Eighty seven out of 482 (18.0%) had been diagnosed with colorectal cancer, 10 of these (11.8%) before 40 years. Cumulative risk at 70 years in our cohort was 34.7% for colorectal cancer and 26.1% for endometrial cancer.

CONCLUSIONS

After 15 years of genetic testing, 29 different class 4/5 variants have been detected in Norway. Almost half of Norwegian PMS2 carriers have the founder variant 989-1G > T. Penetrance of colorectal cancer in our cohort was moderate but variable, as 11.5% of those diagnosed were younger than 40 years.

Germline functional variants contribute to somatic mutation and outcomes in neuroblastoma

Germline genetic context may play a significant role in the development and evolution of cancer, particularly in childhood cancers such as neuroblastoma. This study investigates the role of putatively functional germline variants in neuroblastoma, even if they do not directly increase disease risk. Our whole-exome sequencing analysis of 125 patients with neuroblastoma reveals a positive correlation between germline variant burden and somatic mutations. Moreover, patients with higher germline variant burden exhibit worse outcomes. Similar findings are observed in the independent neuroblastoma cohort where a higher germline variant burden correlates with a higher somatic mutational burden and a worse overall survival outcome. However, contrasting results emerge in adult-onset cancer, emphasizing the importance of germline genetics in neuroblastoma. The enrichment of putatively functional germline variants in cancer predisposition genes is borderline significant when compared to healthy populations (P = 0.077; Odds Ratio, 1.45; 95% confidence intervals, 0.94-2.21) and significantly more pronounced against adult-onset cancers (P = 0.016; Odds Ratio, 2.13; 95% confidence intervals, 1.10-3.91). Additionally, the presence of these variants proves to have prognostic significance in neuroblastoma (log-rank P < 0.001), and combining germline with clinical risk factors notably improves survival predictions.

Case Report of Friedreich's Ataxia and ALG1-Related Biochemical Abnormalities in a Patient With Progressive Spastic Paraplegia

Frataxin is an evolutionarily conserved mitochondrial protein responsible for iron homeostasis and metabolism. A deficiency of frataxin (encoded by FXN) leads to Friedreich's ataxia (FRDA), a progressive disorder that affects both the central and peripheral nervous systems, most commonly via a pathogenic GAA trinucleotide expansion. In contrast, pathogenic variants in ALG1 in humans cause a form of congenital disorder of glycosylation. Here, we present a 15-year-old boy with a clinical presentation that raised concern for complex hereditary spastic paraplegia (HSP), with motor features including progressive spastic paraparesis, cervical dystonia, cerebellar dysfunction, and diminished lower extremity reflexes. The proband was initially found to have a novel compound heterozygous variant in ALG1 on exome sequencing, along with N-glycan profiling revealing evidence of defective mannosylation and Western blot analysis demonstrating an 84% reduction in ALG1 expression. Although several of his clinical features could be explained by the ALG1 variant specifically or considered as part of the presentation of CDGs in general, there were additional phenotypes that suggested an alternative, or additional, genetic diagnosis. Subsequently, he was found to have biallelic pathogenic GAA repeat expansions in FXN on genome sequencing, leading to a diagnosis of FRDA. Given that FRDA explained all his clinical features, the ALG1 variant may have been a hypomorphic form and/or a biochemical phenotype. Our findings underscore the importance of considering FRDA as a differential diagnosis in cases of complex HSP and demonstrate the utility of unbiased genome sequencing approaches that include detection of trinucleotide repeat expansions for progressive motor disorders.

Gene expansions contributing to human brain evolution

Genomic drivers of human-specific neurological traits remain largely undiscovered. Duplicated genes expanded uniquely in the human lineage likely contributed to brain evolution, including the increased complexity of synaptic connections between neurons and the dramatic expansion of the neocortex. Discovering duplicate genes is challenging because the similarity of paralogs makes them prone to sequence-assembly errors. To mitigate this issue, we analyzed a complete telomere-to-telomere human genome sequence (T2T-CHM13) and identified 213 duplicated gene families likely containing human-specific paralogs (>98% identity). Positing that genes important in universal human brain features should exist with at least one copy in all modern humans and exhibit expression in the brain, we narrowed in on 362 paralogs with at least one copy across thousands of ancestrally diverse genomes and present in human brain transcriptomes. Of these, 38 paralogs co-express in gene modules enriched for autism-associated genes and potentially contribute to human language and cognition. We narrowed in on 13 duplicate gene families with human-specific paralogs that are fixed among modern humans and show convincing brain expression patterns. Using long-read DNA sequencing revealed hidden variation across 200 modern humans of diverse ancestries, uncovering signatures of selection not previously identified, including possible balancing selection of CD8B . To understand the roles of duplicated genes in brain development, we generated zebrafish CRISPR "knockout" models of nine orthologs and transiently introduced mRNA-encoding paralogs, effectively "humanizing" the larvae. Morphometric, behavioral, and single-cell RNA-seq screening highlighted, for the first time, a possible role for GPR89B in dosage-mediated brain expansion and FRMPD2B function in altered synaptic signaling, both hallmark features of the human brain. Our holistic approach provides important insights into human brain evolution as well as a resource to the community for studying additional gene expansion drivers of human brain evolution.

Abstract short

Duplicated genes expanded in the human lineage likely contributed to brain evolution, yet challenges exist in their discovery due to sequence-assembly errors. We used a complete telomere-to-telomere genome sequence to identify 213 human-specific gene families. From these, 362 paralogs were found in all modern human genomes tested and brain transcriptomes, making them top candidates contributing to human-universal brain features. Choosing a subset of paralogs, we used long-read DNA sequencing of hundreds of modern humans to reveal previously hidden signatures of selection. To understand their roles in brain development, we generated zebrafish CRISPR "knockout" models of nine orthologs and introduced mRNA-encoding paralogs, effectively "humanizing" larvae. Our findings implicate two new genes in possibly contributing to hallmark features of the human brain: GPR89B in dosage-mediated brain expansion and FRMPD2B in altered synapse signaling. Our holistic approach provides new insights and a comprehensive resource for studying gene expansion drivers of human brain evolution.

Common variants increase risk for congenital diaphragmatic hernia within the context of de novo variants

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly often accompanied by other structural anomalies and/or neurobehavioral manifestations. Rare de novo protein-coding variants and copy-number variations contribute to CDH in the population. However, most individuals with CDH remain genetically undiagnosed. Here, we perform integrated de novo and common-variant analyses using 1,469 CDH individuals, including 1,064 child-parent trios and 6,133 ancestry-matched, unaffected controls for the genome-wide association study. We identify candidate CDH variants in 15 genes, including eight novel genes, through deleterious de novo variants. We further identify two genomic loci contributing to CDH risk through common variants with similar effect sizes among Europeans and Latinx. Both loci are in putative transcriptional regulatory regions of developmental patterning genes. Estimated heritability in common variants is ∼19%. Strikingly, there is no significant difference in estimated polygenic risk scores between isolated and complex CDH or between individuals harboring deleterious de novo variants and individuals without these variants. The data support a polygenic model as part of the CDH genetic architecture.

Trio-based whole-exome sequencing of 200 Chinese patients with keratoconus

Keratoconus (KC) is a complex corneal disorder with a well-recognized genetic component. In this study, we aimed to expand the genetic spectrum of 200 Chinese patients with keratoconus and their unaffected parents. Trio-based whole-exome sequencing was performed in 200 patients with sporadic keratoconus and their unaffected parents. The variants identified in candidate genes for keratoconus were analyzed using multiple bioinformatics tools. Finally, we identified 7 variants in 5 candidate genes for keratoconus in 5 patients. The c.T464C variant in the IMPDH1 gene was defined as likely pathogenic according to the guidelines of the American College of Medical Genetics and Genomics, and the remaining variants in candidate genes (TRANK1, SLC4A11, CERKL, IFT172) were defined as uncertain significance. Our results expand the genetic spectrum in KC, highlight the genetic heterogeneity of this disease and provide important clues for future functional validation.

Role of Forkhead box F1 in the Pathobiology of Pulmonary Arterial Hypertension

Rationale

Approximately 80% of patients with non-familial pulmonary arterial hypertension (PAH) lack identifiable pathogenic genetic variants. While most genetic studies of PAH have focused on predicted loss-of-function variants, recent approaches have identified ultra-rare missense variants associated with the disease. FOXF1 encodes a highly conserved transcription factor, essential for angiogenesis and vasculogenesis in human and mouse lungs.

Objectives

We identified a rare FOXF1 missense coding variant in two unrelated probands with PAH. FOXF1 is an evolutionarily conserved transcription factor required for lung vascular development and vascular integrity. Our aims were to determine the frequency of FOXF1 variants in larger PAH cohorts compared to the general population, study FOXF1 expression in explanted lung tissue from PAH patients versus control (failed-donor) lungs, and define potential downstream targets linked to PAH development.

Methods

Three independent, international, multicenter cohorts were analyzed to evaluate the frequency of FOXF1 rare variants. Various composite prediction models assessed the deleteriousness of individual variants. Bulk RNA sequencing datasets from human explanted lung tissues were compared to failed-donor controls to determine FOXF1 expression. Bioinformatic tools identified putative FOXF1 binding targets, which were orthogonally validated using mouse ChIP-seq datasets.

Measurements and Main Results

Seven novel or ultra-rare missense coding variants were identified across three patient cohorts in different regions of the FOXF1 gene, including the DNA binding domain. FOXF1 expression was dysregulated in PAH lungs, correlating with disease severity. Histological analysis showed heterogeneous FOXF1 expression, with the lowest levels in phenotypically abnormal endothelial cells within complex vascular lesions in PAH samples. A hybrid bioinformatic approach identified FOXF1 downstream targets potentially involved in PAH pathogenesis, including BMPR2 .

Conclusions

Large genomic and transcriptomic datasets suggest that decreased FOXF1 expression or predicted dysfunction is associated with PAH.

Estimation of demography and mutation rates from one million haploid genomes

As genetic sequencing costs have plummeted, datasets with sizes previously un-thinkable have begun to appear. Such datasets present new opportunities to learn about evolutionary history, particularly via rare alleles that record the very recent past. However, beyond the computational challenges inherent in the analysis of many large-scale datasets, large population-genetic datasets present theoretical problems. In particular, the majority of population-genetic tools require the assumption that each mutant allele in the sample is the result of a single mutation (the "infinite sites" assumption), which is violated in large samples. Here, we present DR EVIL, a method for estimating mutation rates and recent demographic history from very large samples. DR EVIL avoids the infinite-sites assumption by using a diffusion approximation to a branching-process model with recurrent mutation. The branching-process approach limits the method to rare alleles, but, along with recent results, renders tractable likelihoods with recurrent mutation. We show that DR EVIL performs well in simulations and apply it to rare-variant data from a million haploid samples, identifying a signal of mutation-rate heterogeneity within commonly analyzed classes and predicting that in modern sample sizes, most rare variants at sites with high mutation rates represent the descendants of multiple mutation events.

Saturation mutagenesis-reinforced functional assays for disease-related genes

Interpretation of disease-causing genetic variants remains a challenge in human genetics. Current costs and complexity of deep mutational scanning methods are obstacles for achieving genome-wide resolution of variants in disease-related genes. Our framework, saturation mutagenesis-reinforced functional assays (SMuRF), offers simple and cost-effective saturation mutagenesis paired with streamlined functional assays to enhance the interpretation of unresolved variants. Applying SMuRF to neuromuscular disease genes FKRP and LARGE1, we generated functional scores for all possible coding single-nucleotide variants, which aid in resolving clinically reported variants of uncertain significance. SMuRF also demonstrates utility in predicting disease severity, resolving critical structural regions, and providing training datasets for the development of computational predictors. Overall, our approach enables variant-to-function insights for disease genes in a cost-effective manner that can be broadly implemented by standard research laboratories.

Unraveling the mysteries of MGMT: Implications for neuroendocrine tumors

Neuroendocrine tumors (NETs) are a diverse group of tumors that arise from neuroendocrine cells and are commonly found in various organs. A considerable proportion of NET patients were diagnosed at an advanced or metastatic stage. Alkylating agents are the primary treatment for NET, and O6-methylguanine methyltransferase (MGMT) remains the first-line of defense against DNA damage caused by these agents. Clinical trials have indicated that MGMT promoter methylation or its low/lacked expression can predict a favorable outcome with Temozolomide in NETs. Its status could help select NET patients who can benefit from alkylating agents. Therefore, MGMT status serves as a biomarker to guide decisions on the efficacy of Temozolomide as a personalized treatment option. Additionally, delving into the regulatory mechanisms of MGMT status can lead to the development of MGMT-targeted therapies, benefiting individuals with high levels of MGMT expression. This review aims to explore the polymorphism of MGMT regulation and summarize its clinical implications in NETs, which would help establish the role of MGMT as a biomarker and its potential as a therapeutic target in NETs. Additionally, we explore the benefits of combining Temozolomide and immunotherapy in MGMT hypermethylated subgroups. Future studies can focus on optimizing Temozolomide administration to induce specific immunomodulatory changes.

Germline variant affecting p53β isoforms predisposes to familial cancer

Germline and somatic TP53 variants play a crucial role during tumorigenesis. However, genetic variations that solely affect the alternatively spliced p53 isoforms, p53β and p53γ, are not fully considered in the molecular diagnosis of Li-Fraumeni syndrome and cancer. In our search for additional cancer predisposing variants, we identify a heterozygous stop-lost variant affecting the p53β isoforms (p.*342Serext*17) in four families suspected of an autosomal dominant cancer syndrome with colorectal, breast and papillary thyroid cancers. The stop-lost variant leads to the 17 amino-acid extension of the p53β isoforms, which increases oligomerization to canonical p53α and dysregulates the expression of p53's transcriptional targets. Our study reveals the capacity of p53β mutants to influence p53 signalling and contribute to the susceptibility of different cancer types. These findings underscore the significance of p53 isoforms and the necessity of comprehensive investigation into the entire TP53 gene in understanding cancer predisposition.

Genetics and Traumatic Brain Injury: Findings from an Exome-Based Study of a 50-Patient Case Series

Traumatic brain injury (TBI) is the leading cause of global mortality and morbidity. Because TBI is accident-related, the role of genetics in predisposing to TBI has been largely unexplored. However, the likelihood of injury may not be entirely random and may be associated with certain physical and mental characteristics. In this study, we analyzed the exomes of 50 patients undergoing rehabilitation after TBI. Patients were divided into three groups according to rehabilitation outcome: improvement, no change, and deterioration/death. We focused on rare, potentially functional missense and high-impact variants in genes intolerant to these variants. The concordant results from the three independent groups of patients allowed for the suggestion of the existence of a genetic predisposition to TBI, associated with rare functional variations in intolerant genes, with a prevalent dominant mode of inheritance and neurological manifestations in the genetic phenotypes according to the OMIM database. Forty-four of the 50 patients had one or more rare, potentially deleterious variants in one or more neurological genes. Comparison of these results with those of a 50-sampled matched non-TBI cohort revealed significant differences: P = 2.6 × 10-3, OR = 4.89 (1.77-13.47). There were no differences in the distribution of the genes of interest between the TBI patient groups. Our exploratory study provides new insights into the impact of genetics on TBI risk and is the first to address potential genetic susceptibility to TBI.

Whole genome sequencing in adults with clinical hallmarks of hypophosphatasia negative for ALPL variants

BACKGROUND

Hypophosphatasia (HPP) is a rare disease caused by deficient activity of tissue-nonspecific alkaline phosphatase (ALP), encoded by the ALPL gene. The primary objective was to explore novel ALPL variants by whole genome sequencing (WGS) in patients with HPP who previously tested negative by standard methods for ALPL variants. The secondary objective was to search for genes beyond ALPL that may reduce ALP activity or contribute to HPP symptoms.

METHODS AND RESULTS

WGS was performed in 16 patients clinically diagnosed with HPP who had ALP activity below the normal range and tested negative for ALPL variants. Genetic variants in ALPL and genes possibly associated with low ALP activity or phenotypic overlap with HPP were assessed. All 16 patients had ALP activity below the normal range. WGS did not identify any novel disease-causing ALPL variants. Positive findings for other gene variants were identified in 4 patients: 1 patient presented with variants in COL1A1, NLRP12, and SCN9A, coding for collagen, type, I alpha-1 chain, nod-like receptor pyrin domain containing 12, and sodium voltage-gated channel alpha subunit 9, respectively; 1 presented with a heterozygous, likely pathogenic variant in P3H1 coding for prolyl 3 hydroxylase 1; 1 presented with a heterozygous pathogenic variant in SGCE, coding for sarcoglycan epsilon; and 1 presented with a heterozygous variant of uncertain significance in VDR, encoding vitamin D receptor.

CONCLUSION

Genomic analysis did not identify novel ALPL variants or a pattern of disease-causing variants in genes other than ALPL to explain the HPP phenotype in these patients.

REGISTRATION

Clinicaltrials.gov identifier: NCT04925804.

Genetic investigations on singleton school aged children reveal novel variants and new candidate genes for hearing loss

Hearing loss affects around 5% of the global population. Two preliminary studies have described genetic variants in sporadic individuals with hearing loss from Pakistan. Here we extend these studies to determine the spectrum of variants in a cohort of individuals with no previous history of hearing loss. Individuals with hearing loss born to consanguineous couples were identified from special schools. Audiograms were assessed. DNA from participants negative for GJB2 pathogenic variants was subjected to exome sequencing. Data were filtered to include variants with frequencies < 0.01 in the public databases. The effects of the missense variants on respective amino acids were analyzed by using PyMol software. Among the 44 participants, hearing loss was moderate for two individuals; 14 exhibited moderately-severe hearing loss while 25 had a severe degree of hearing loss. Hearing loss was reported to have been progressive in four participants and was currently profound in three participants. Variants were unambiguously identified in 17 genes, of which the majority affected SLC26A4. CDH23, MYO15A and OTOF were other significant contributors. Deleterious variants detected in two genes suggest new associations for hearing loss. Molecular characterization of hearing loss in our cohort revealed high genetic heterogeneity with a 75% diagnostic rate.

A pleiotropic recurrent dominant ITPR3 variant causes a complex multisystemic disease

Inositol 1,4,5-trisphosphate (IP3) receptor type 1 (ITPR1), 2 (ITPR2), and 3 (ITPR3) encode the IP3 receptor (IP3R), a key player in intracellular calcium release. In four unrelated patients, we report that an identical ITPR3 de novo variant-NM_002224.3:c.7570C>T, p.Arg2524Cys-causes, through a dominant-negative effect, a complex multisystemic disorder with immunodeficiency. This leads to defective calcium homeostasis, mitochondrial malfunction, CD4+ lymphopenia, a quasi-absence of naïve CD4+ and CD8+ cells, an increase in memory cells, and a distinct TCR repertoire. The calcium defect was recapitulated in Jurkat knock-in. Site-directed mutagenesis displayed the exquisite sensitivity of Arg2524 to any amino acid change. Despite the fact that all patients had severe immunodeficiency, they also displayed variable multisystemic involvements, including ectodermal dysplasia, Charcot-Marie-Tooth disease, short stature, and bone marrow failure. In conclusion, unlike previously reported ITPR1-3 deficiencies leading to narrow, mainly neurological phenotypes, a recurrent dominant ITPR3 variant leads to a multisystemic disease, defining a unique role for IP3R3 in the tetrameric IP3R complex.

AI-derived comparative assessment of the performance of pathogenicity prediction tools on missense variants of breast cancer genes

Single nucleotide variants (SNVs) can exert substantial and extremely variable impacts on various cellular functions, making accurate predictions of their consequences challenging, albeit crucial especially in clinical settings such as in oncology. Laboratory-based experimental methods for assessing these effects are time-consuming and often impractical, highlighting the importance of in-silico tools for variant impact prediction. However, the performance metrics of currently available tools on breast cancer missense variants from benchmarking databases have not been thoroughly investigated, creating a knowledge gap in the accurate prediction of pathogenicity. In this study, the benchmarking datasets ClinVar and HGMD were used to evaluate 21 Artificial Intelligence (AI)-derived in-silico tools. Missense variants in breast cancer genes were extracted from ClinVar and HGMD professional v2023.1. The HGMD dataset focused on pathogenic variants only, to ensure balance, benign variants for the same genes were included from the ClinVar database. Interestingly, our analysis of both datasets revealed variants across genes with varying penetrance levels like low and moderate in addition to high, reinforcing the value of disease-specific tools. The top-performing tools on ClinVar dataset identified were MutPred (Accuracy = 0.73), Meta-RNN (Accuracy = 0.72), ClinPred (Accuracy = 0.71), Meta-SVM, REVEL, and Fathmm-XF (Accuracy = 0.70). While on HGMD dataset they were ClinPred (Accuracy = 0.72), MetaRNN (Accuracy = 0.71), CADD (Accuracy = 0.69), Fathmm-MKL (Accuracy = 0.68), and Fathmm-XF (Accuracy = 0.67). These findings offer clinicians and researchers valuable insights for selecting, improving, and developing effective in-silico tools for breast cancer pathogenicity prediction. Bridging this knowledge gap contributes to advancing precision medicine and enhancing diagnostic and therapeutic approaches for breast cancer patients with potential implications for other conditions.

Evaluation of plasma phytosterols in sitosterolemia, their kindreds and hyperlipidemia subjects

BACKGROUND

Patients suffering from sitosterolemia with ABCG5/8 mutation typically present with early-onset or rapidly progressive atherosclerosis. Their kindreds with partial genetic deficiencies of ABCG5/8 are often considered healthy. However, discerning sitosterolemia from its familial kindreds and hyperlipidemia subjects has remained challenging.

METHODS

Here we retrospectively recruited seven families including 8 individuals diagnosed with sitosterolemia subjects, and 14 kindreds carrying single gene mutations. Additionally, 17 individuals with hyperlipidemia and 130 healthy controls served as positive and negative controls, respectively. A total of 6 phytosterols combined with cholesterol absorption indices (including sitosterol, campesterol, stigmasterol, and cholestanol) and cholesterol synthesis markers (desmosterol and 7-dehydrocholesterol), was compared across the aforementioned four groups.

RESULTS

As expected, the sitosterolemia subjects with double mutations demonstrated significantly elevated levels of sitosterol and other cholesterol absorption indices. Meanwhile, sitosterolemia kindreds with single gene mutation showed a similar pattern of activated cholesterol-absorption ability to the hyperlipidemia group, but not as high as the double mutation group. Notably, the cholesterol-synthesis enzyme 7-dehydrocholesterol reductase displayed an increase in the hyperlipidemia group but a decrease in the sitosterolemia kindred group, suggesting a potential discriminative role of 7-dehydrocholesterol in distinguishing between these two groups. The combination of phytosterols was more valuable than clinical lipid index for sitosterolemia diagnosis.

CONCLUSION

Our study revealed mild disruptions of cholesterol absorption capacities in sitosterolemia kindreds with single mutations. Furthermore, the combination of 6 phytosterols proved effective in distinguishing between sitosterolemia, its single mutation carriers, and hyperlipidemia patients.

A novel homozygous splicing mutation in AK7 causes multiple morphological abnormalities of sperm flagella in patients from consanguineous Pakistani families

Multiple morphological abnormalities of the flagella (MMAF) represent a severe form of sperm defects leading to asthenozoospermia and male infertility. In this study, we identified a novel homozygous splicing mutation (c.871-4 ACA>A) in the adenylate kinase 7 (AK7) gene by whole-exome sequencing in infertile individuals. Spermatozoa from affected individuals exhibited typical MMAF characteristics, including coiled, bent, short, absent, and irregular flagella. Transmission electron microscopy analysis showed disorganized axonemal structure and abnormal mitochondrial sheets in sperm flagella. Immunofluorescence staining confirmed the absence of AK7 protein from the patients' spermatozoa, validating the pathogenic nature of the mutation. This study provides direct evidence linking the AK7 gene to MMAF-associated asthenozoospermia in humans, expanding the mutational spectrum of AK7 and enhancing our understanding of the genetic basis of male infertility.

Genetic variation and molecular profiling of congenital malformations of the female genital tract based on whole-genome sequencing

BACKGROUND

Congenital malformations of the female genital tract (CM-FGT) are characterized by abnormal development of the fallopian tubes, uterus, and vagina, often accompanied by malformations in the urinary system, bones and hearing. However, no definitive pathogenic genes and molecular genetic causes have been identified.

METHODS

We present the largest whole-genome sequencing study of CM-FGT to date, analyzing 590 individuals in China: 95 patients, 442 case-controls, and 53 familial controls.

RESULTS

Among the patients, 5.3% carried known CM-FGT-related variants. Pedigree and case-control analyses in two dimensions of coding and non-coding regulatory regions revealed seven novel de novo copy number variations, 12 rare single-nucleotide variations, and 10 rare 3' untranslated region (UTR) mutations in genes related to CM-FGT, particularly highlighting ASH1L as a pathogenic gene. Single-cell sequencing data showed that the majority of CM-FGT-related risk genes are spatiotemporally specifically expressed early in uterus development.

CONCLUSIONS

In conclusion, this study identified novel variants related to CM-FGT, particularly highlighting ASH1L as a pathogenic gene. The findings provide insights into the genetic variants underlying CM-FGT, with single-cell sequencing data revealing spatiotemporal specific expression patterns of key risk genes early in uterine development. This study significantly advances the understanding of CM-FGT etiology and genetic landscape, offering new opportunities for prenatal screening.

The FXR1 network acts as a signaling scaffold for actomyosin remodeling

It is currently not known whether mRNAs fulfill structural roles in the cytoplasm. Here, we report the fragile X-related protein 1 (FXR1) network, an mRNA-protein (mRNP) network present throughout the cytoplasm, formed by FXR1-mediated packaging of exceptionally long mRNAs. These mRNAs serve as an underlying condensate scaffold and concentrate FXR1 molecules. The FXR1 network contains multiple protein binding sites and functions as a signaling scaffold for interacting proteins. We show that it is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. Point mutations in FXR1, found in its homolog FMR1, where they cause fragile X syndrome, disrupt the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. Our findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as an organizer of signaling reactions.

WDR44 Loss-of-Function Promoter Deletion in a Male Newborn With a Ciliopathy Phenotype

Gain-of-function variants in the WDR44 gene have recently been associated with an X-linked ciliopathy-related neurodevelopmental phenotype. Here, we report on a WDR44 loss-of-function (LOF) variant identified in the genome sequence from a male fetus enrolled in the Prenatal Genetic Diagnosis by Genomic Sequencing (PrenatalSEQ) multicenter study. The phenotype is consistent with the described X-linked ciliopathy that includes developmental delay, microcephaly, congenital heart defects, kidney abnormalities, cryptorchidism, musculoskeletal abnormalities, craniofacial dysmorphism, and effusions. This is the first report of a WDR44 LOF variant in an affected individual with a prenatal presentation and supports LOF as a mechanism for the X-linked WDR44 ciliopathy-related phenotype.

Identification of Sweetness Preference-Related Single-Nucleotide Polymorphisms for Polygenic Risk Scores Associated with Obesity

Our study aimed to identify sweetness preference-associated single-nucleotide polymorphisms (SNPs), characterize the related genetic loci, and develop SNP-based polygenic risk scores (PRS) to analyze their associations with obesity. For genotyping, we utilized a pooled genome-wide association study (GWAS) dataset of 18,499 females and 10,878 males. We conducted genome-wide association analyses, functional annotation, and employed the weighted method to calculate the levels of PRS from 677 sweetness preference-related SNPs. We used Cox proportional hazards modeling with time-varying covariates to estimate age-adjusted and multivariable hazard ratios (HRs) and 95% confidence intervals (CIs) for obesity incidence. We also tested the correlation between PRS and environmental factors, including smoking and dietary components, on obesity. Our results showed that in males, the TT genotype of rs4861982 significantly increased obesity risk compared to the GG genotype in the Health Professionals Follow-up Study (HPFS) cohort (HR = 1.565; 95% CI, 1.122-2.184; p = 0.008) and in the pooled analysis (HR = 1.259; 95% CI, 1.030-1.540; p = 0.025). Protein tyrosine phosphatase receptor type O (PTPRO) was identified as strongly associated with sweetness preference, indicating a positive correlation between sweetness preference and obesity risk. Moreover, each 10 pack-year increment in smoking was significantly associated with an increased risk of obesity in the HPFS cohort (HR = 1.024; 95% CI, 1.000-1.048) in males but not in females. In conclusion, significant associations between rs4861982, sweetness preference, and obesity were identified, particularly among males, where environmental factors like smoking are also correlated with obesity risk.

Rare RNF213 variants is related to early-onset intracranial atherosclerosis: A Chinese community-based study

BACKGROUND

The relationship between rare variants in Ring finger protein 213 (RNF213) and intracranial atherosclerosis (ICAS) remained unelucidated. Using whole-exome sequencing (WES) and high-resolution magnetic resonance imaging (HR-MRI), this study aimed at investigating the association between rare RNF213 variants and ICAS within a Chinese community-dwelling population.

METHODS

The present study included 821 participants from Shunyi cohort. Genetic data of rare RNF213 variants were acquired by WES and were categorized by functional domains. Intracranial and extracranial atherosclerosis were assessed by brain HR-MRI and carotid ultrasound, respectively. Logistic regression and generalized linear regression were applied to evaluate the effects of rare RNF213 variants on atherosclerosis. Stratification by age were conducted with 50 years old set as the cutoff value.

RESULTS

Ninety-five participants were identified as carriers of rare RNF213 variants. Carotid plaques were observed in 367 (44.7 %) participants, while ICAS was identified in 306 (37.3 %). Rare variants of RNF213 was not associated with ECAS. Employing HR-MRI, both the presence of rare variants (β = 0.150, P = 0.025) and numerical count of variants (β = 0.182, P = 0.003) were significantly correlated with ICAS within the group of age ≤50 years. Both variant existence (β = 0.154, P = 0.014) and variant count (β = 0.188, P = 0.003) were significantly associated with plaques in middle cerebral arteries within younger subgroup, rather than basilar arteries. Furthermore, a significant association was observed between variants that located outside the N-arm domain and ICAS in the younger subgroup (OR = 2.522, P = 0.030). Statistical results remained robust after adjusted for age, gender, and cardiovascular risk factors.

CONCLUSIONS

Rare variants of RNF213 is associated with age-related ICAS in general Chinese population, highlighting the potential role of RNF213 as a genetic contributor to early-onset ICAS.

De Novo Missense Variations of ATP8B2 Impair Its Phosphatidylcholine Flippase Activity

P4-ATPases comprise a family of lipid flippases that translocate lipids from the exoplasmic (or luminal) to the cytoplasmic leaflet of biological membranes. Of the 14 known human P4-ATPases, ATP8B2 is a phosphatidylcholine flippase at the plasma membrane, but its physiological function is not well understood. Although ATP8B2 could interact with both CDC50A and CDC50B, it required only the CDC50A interaction for its exit from the endoplasmic reticulum and subsequent transport to the plasma membrane. Three de novo monoallelic missense variations of ATP8B2 were found in patients with intellectual disability. None of these variations affected the interaction of ATP8B2 with CDC50A or its localization to the plasma membrane. However, variations of either of two amino acid residues, which are conserved in all P4-ATPases, significantly reduced the phosphatidylcholine flippase activity of ATP8B2. Furthermore, mutations in the corresponding residues of ATP8B1 and ATP11C were found to decrease their flippase activities toward phosphatidylcholine and phosphatidylserine, respectively. These results indicate that the conserved amino acid residues are crucial for the enzymatic activities of the P4-ATPases.

Whole exome sequencing reveals a dual diagnosis of BCAP31-related syndrome and glutaric aciduria III

Background

Biochemical testing is a common first-tier approach in the setting of genetic evaluation of patients with unexplained developmental delay. However, results can be unclear, and a plan for second-tier analysis must be determined based on the patient's biochemical results and clinical presentation - in many cases, triggering a diagnostic odyssey.

Case presentation

A male patient from the United States presenting with unexplained developmental delay, microcephaly, hypotonia, and feeding difficulties was referred for clinical genetic evaluation at age 8 months. Biochemical testing revealed an isolated marked elevation of glutaric acid on urine organic acid profile, without elevations of related metabolites. Further testing included GCDH sequencing, a neurometabolic gene panel, chromosomal microarray, Prader Willi/Angelman testing, and lysosomal disease enzyme panel, all of which were non-diagnostic. The patient had persistent developmental delay and hypotonia, dystonia, sensorineural hearing loss, and abnormal brain myelination on magnetic resonance imaging. Whole exome sequencing (WES) was performed and revealed a dual diagnosis of glutaric aciduria III (GA III) and BCAP31-related disorder, an X-linked intellectual disability syndrome, caused by a novel pathogenic variant.

Conclusions

GA III has historically been considered clinically benign, with few reported cases. This patient's presenting symptoms were similar to those commonly seen in GA I and GA II, however the biochemical abnormalities were not consistent with these disorders, prompting additional molecular and biochemical testing. Ultimately, WES confirmed a diagnosis of BCAP31-related syndrome, a rare neurological disorder, which explained the patient's presenting symptoms. WES also identified a secondary diagnosis of GA III. We present a patient with two rare genetic conditions, highlighting the importance of deep phenotyping and the utility of WES in the setting of a patient with dual genetic diagnoses.

Pontocerebellar Hypoplasia Type 9: A New Case with a Novel Mutation and Review of Literature

Pontocerebellar hypoplasia type 9 (PCH-9) is a very rare autosomal recessive neurodegenerative disorder. Affected infants present early with severe developmental delay, spasticity, with the unique magnetic resonance imaging picture of thin corpus callosum, atrophied pons, and cerebellum. It is caused by loss of function mutations in the AMPD2 gene, encoding for the adenosine monophosphate deaminase enzyme-paralog 2. This gene is expressed in different somatic tissues with high level of expression in cerebellum and its encoded enzyme catalyzes a critical step in de novo biosynthesis of purines and its deficiency in the developing neurons severely affects neuronal differentiation and cell viability. We clinically evaluated an Emirati patient presented with severe developmental and growth delay, as well as corpus callosum agenesis and atrophy of brainstem and cerebellum. We performed exome sequencing, Sanger sequencing, and segregation analysis to identify the genetic cause of the phenotype, followed by in silico and in vitro analysis. We identified the novel variant (NM_004037.9:c.1471G > A) in AMPD2 gene leading to a single amino acid substitution (p.Gly491Arg) in adenosine monophosphate deaminase-2 enzyme. This variant is predicted to be pathogenic using several in silico tools, and resulted in a decrease in the enzyme function in the patient's polymorphonuclear cells by 82% (95% confidence interval: 73.3-91.7%, p  = 0.029) compared with the control. This data establishes that the affected child is affected by PCH-9. Furthermore, we review all reported cases in literature to summarize the main clinical features of this rare disease.

The human proteome size as a technological development function

Changes in information on the number of human proteoforms, post-translational modification (PTM) events, alternative splicing (AS), single-amino acid polymorphisms (SAP) associated with protein-coding genes in the neXtProt database have been retrospectively analyzed. In 2016, our group proposed three mathematical models for predicting the number of different proteins (proteoforms) in the human proteome. Eight years later, we compared the original data of the information resources and their contribution to the prediction results, correlating the differences with new approaches to experimental and bioinformatic analysis of protein modifications. The aim of this work is to update information on the status of records in the databases of identified proteoforms since 2016, as well as to identify trends in changes in the quantities of these records. According to various information models, modern experimental methods may identify from 5 to 125 million different proteoforms: the proteins formed due to alternative splicing, the implementation of single nucleotide polymorphisms at the proteomic level, and post-translational modifications in various combinations. This result reflects an increase in the size of the human proteome by 20 or more times over the past 8 years.

Context-adjusted proportion of singletons (CAPS): a novel metric for assessing negative selection in the human genome

Interpretation of genetic variants remains challenging, partly due to the lack of well-established ways of determining the potential pathogenicity of genetic variation, especially for understudied classes of variants. Addressing this, population genetics methods offer a practical solution by evaluating variant effects through human population distributions. Negative selection influences the ratio of singleton variants and can serve as a proxy for deleteriousness, as exemplified by the Mutability-Adjusted Proportion of Singletons (MAPS) metric. However, MAPS is sensitive to the calibration of the singletons-by-mutability linear model, which results in biased estimates for certain variant classes. Building up on the methodology used in MAPS, we introduce the Context-Adjusted Proportion of Singletons (CAPS) metric for assessing negative selection in the human genome. CAPS produces corrected estimates with more accurate confidence intervals by eliminating the mutability layer in the model. Retaining the advantageous features of MAPS, CAPS emerges as a robust and reliable tool. We believe that CAPS has the potential to enhance the identification of new disease-variant associations in clinical and research settings, offering improved accuracy in assessing negative selection for diverse SNV classes.

Rare coding variant analysis for human diseases across biobanks and ancestries

Large-scale sequencing has enabled unparalleled opportunities to investigate the role of rare coding variation in human phenotypic variability. Here, we present a pan-ancestry analysis of sequencing data from three large biobanks, including the All of Us research program. Using mixed-effects models, we performed gene-based rare variant testing for 601 diseases across 748,879 individuals, including 155,236 with ancestry dissimilar to European. We identified 363 significant associations, which highlighted core genes for the human disease phenome and identified potential novel associations, including UBR3 for cardiometabolic disease and YLPM1 for psychiatric disease. Pan-ancestry burden testing represented an inclusive and useful approach for discovery in diverse datasets, although we also highlight the importance of ancestry-specific sensitivity analyses in this setting. Finally, we found that effect sizes for rare protein-disrupting variants were concordant between samples similar to European ancestry and other genetic ancestries (βDeming = 0.7-1.0). Our results have implications for multi-ancestry and cross-biobank approaches in sequencing association studies for human disease.

A broken network of susceptibility genes in the monocytes of Crohn's disease patients

Genome-wide association studies have identified over 200 genetic loci associated with inflammatory bowel disease; however, the mechanism of such a large amount of susceptibility genes remains uncertain. In this study, we integrated bioinformatics analysis and two independent single-cell transcriptome datasets to investigate the expression network of 232 susceptibility genes in Crohn's disease (CD) patients and healthy controls. The study revealed that most of the susceptibility genes are specifically and strictly expressed in the monocytes of the human intestinal tract. The susceptibility genes established a network within the monocytes of health control. The robustness of a gene network may prevent disease onset that is influenced by the genetic and environmental alteration in the expression of susceptibility genes. In contrast, we showed a sparse network in pediatric/adult CD patients, suggesting the broken network contributed to the CD etiology. The network status of susceptibility genes at the single-cell level of monocytes provided novel insight into the etiology.

The Conundrum of Mechanics Versus Genetics in Congenital Hydrocephalus and Its Implications for Fetal Therapy Approaches: A Scoping Review

Recent advances in gene therapy, particularly for single-gene disorders (SGDs), have led to significant progress in developing innovative precision medicine approaches that hold promise for treating conditions such as primary hydrocephalus (CH), which is characterized by increased cerebrospinal fluid (CSF) volumes and cerebral ventricular dilation as a result of impaired brain development, often due to genetic causes. CH is a significant contributor to childhood morbidity and mortality and a driver of healthcare costs. In many cases, prenatal ultrasound can readily identify ventriculomegaly as early as 14-20 weeks of gestation, with severe cases showing poor neurodevelopmental outcomes. Postnatal surgical approaches, such as ventriculoperitoneal shunts, do not address the underlying genetic causes, have high complication rates, and result in a marginal improvement of neurocognitive deficits. Prenatal somatic cell gene therapy (PSCGT) promises a novel approach to conditions such as CH by targeting genetic mutations in utero, potentially improving long-term outcomes. To better understand the pathophysiology, genetic basis, and molecular pathomechanisms of CH, we conducted a scoping review of the literature that identified over 160 published genes linked to CH. Mutations in L1CAM, TRIM71, MPDZ, and CCDC88C play a critical role in neural stem cell development, subventricular zone architecture, and the maintenance of the neural stem cell niche, driving the development of CH. Early prenatal interventions targeting these genes could curb the development of the expected CH phenotype, improve neurodevelopmental outcomes, and possibly limit the need for surgical approaches. However, further research is needed to establish robust genotype-phenotype correlations and develop safe and effective PSCGT strategies for CH.

Prioritizing disease-related rare variants by integrating gene expression data

Rare variants, comprising the vast majority of human genetic variations, are likely to have more deleterious impact in the context of human diseases compared to common variants. Here we present carrier statistic, a statistical framework to prioritize disease-related rare variants by integrating gene expression data. By quantifying the impact of rare variants on gene expression, carrier statistic can prioritize those rare variants that have large functional consequence in the patients. Through simulation studies and analyzing real multi-omics dataset, we demonstrated that carrier statistic is applicable in studies with limited sample size (a few hundreds) and achieves substantially higher sensitivity than existing rare variants association methods. Application to Alzheimer's disease reveals 16 rare variants within 15 genes with extreme carrier statistics. We also found strong excess of rare variants among the top prioritized genes in patients compared to that in healthy individuals. The carrier statistic method can be applied to various rare variant types and is adaptable to other omics data modalities, offering a powerful tool for investigating the molecular mechanisms underlying complex diseases.

Molecular diagnostic yield for Blau syndrome in previously diagnosed juvenile idiopathic arthritis with uveitis or cutaneous lesions

OBJECTIVE

Diagnostic pitfalls often arise in the community because of potentially misleading similarities between juvenile idiopathic arthritis (JIA) and Blau syndrome, an immune-related disorder caused by NOD2 gene mutations. It remains unclear in which population and to what extent next-generation sequencing techniques can aid in diagnosis.

METHODS

We evaluated clinical usefulness of targeted next-generation sequencing in previously diagnosed JIA. Participants were required to have symptoms and signs suspected of Blau syndrome, including at least uveitis or cutaneous lesions in addition to arthritis. Targeted sequencing was conducted on NOD2 gene to detect diagnostic variants classified as pathogenic or likely pathogenic for Blau syndrome. We assessed the molecular diagnostic yield and clinical implications for patient care.

RESULTS

Between 1 May 2008 and 1 June 2021, sequencing data were accrued from 123 previously diagnosed JIA (median age: 5 years; female: 62.6%). Targeted NOD2 sequencing yielded a positive molecular diagnosis of Blau syndrome in 21.1% (95% CI: 14.9%, 29.2%), encompassing six heterozygous missense mutations classified as pathogenic variants. Among those receiving a molecular diagnosis, changes in clinical management and treatment were considered as having occurred in 38.5%. Nine predictors were identified as being associated with a higher diagnostic yield, providing clinical clues to suspect the possibility of Blau syndrome.

CONCLUSION

Among some patients with paediatric-onset arthritis complicated with uveitis or cutaneous lesions, reassessment of the diagnosis of JIA may be warranted. Targeted NOD2 sequencing established the molecular diagnosis of Blau syndrome in nearly one-fifth of these cases and provided clinically relevant information for patient-care decisions.

DAB2IP associates with hereditary angioedema: Insights into the role of VEGF signaling in HAE pathophysiology

BACKGROUND

In the recent years, there was an important improvement in the understanding of the pathogenesis of hereditary angioedema (HAE). Notwithstanding, in a large portion of patients with unknown mutation (HAE-UNK) the genetic cause remains to be identified.

OBJECTIVES

To identify new genetic targets associated with HAE, a large Argentine family with HAE-UNK spanning 3 generations was studied.

METHODS

Whole exome sequencing was performed on affected family members to identify potential genetic variants associated with HAE-UNK. In silico analyses and experimental studies were applied to assess the role of the identified gene variant.

RESULTS

A missense variant (p.D239N) in DAB2IP was identified. The variant occurred in the C2-domain, the region interacting with vascular endothelial growth factor receptor 2 (VEGFR2). It was found to be rare, and predicted to have a detrimental effect on the functionality of DAB2IP. Protein structure modeling predicted changes in the mutant p.D239N protein structure, impacting protein stability. The p.D239N variant affected the subcellular localization of VEGFR2. Cells transfected with the DAB2IP-239N transcript exhibited an intracellular distribution, and VEGFR2 remained associated with the cell membrane. The altered localization pattern indicated reduced colocalization of the mutant protein with VEGFR2, suggesting a diminished ability of VEGFR2 binding.

CONCLUSIONS

The study identified a novel missense variant (p.D239N) in DAB2IP in a family with HAE-UNK and highlighted the role of dysregulated VEGF-mediated signaling in altered endothelial permeability. DAB2IP loss-of-function pathogenic variants lead to the impairment of the endothelial VEGF/VEGFR2 ligand system and represent a new pathophysiologic cause of HAE-UNK.

Identification of whole-genome mutations and structural variations of bile cell-free DNA in cholangiocarcinoma

Bile cell-free DNA (cfDNA) has been reported as a promising liquid biopsy tool for cholangiocarcinoma (CCA), however, the whole-genome mutation landscape and structural variants (SVs) of bile cfDNA remains unknown. Here we performed whole-genome sequencing on bile cfDNA and analyzed the correlation between mutation characteristics of bile cfDNA and clinical prognosis. TP53 and KRAS were the most frequently mutated genes, and the RTK/RAS, homologous recombination (HR), and HIPPO were top three pathways containing most gene mutations. Ten overlapping putative driver genes were found in bile cfDNA and tumor tissue. SVs such as chromothripsis and kataegis were identified. Moreover, the hazard ratio of HR pathway mutations were 15.77 (95% CI: 1.571-158.4), patients with HR pathway mutations in bile cfDNA exhibited poorer overall survival (P = 0.0049). Our study suggests that bile cfDNA contains genome mutations and SVs, and HR pathway mutations in bile cfDNA can predict poor outcomes of CCA patients.