The mutational constraint spectrum quantified from variation in 141,456 humans

Across ancestries, HLA-B∗08:01∼DRB1∗03:01 (DR3) and HLA-DQA∗01:02 (DR2) increase the risk to develop juvenile-onset systemic lupus erythematosus through low complement C4 levels

Objective

Systemic lupus erythematosus (SLE) is a systemic autoimmune/inflammatory disease with a strong genetic component. Genetic burden is higher in children when compared to patients with adult-onset SLE, contributing to earlier disease expression and more severe phenotypes. The human leukocyte antigen (HLA) cluster on chromosome 6p21.3 is among the most variable genomic regions, representing a major risk-factor for SLE in adults. Its impact on juvenile-onset (j)SLE remains largely unstudied.

Methods

High-resolution sequencing of HLA class I (A, B, C), class II (DRB1, DQA1, DQB1) and class III (complement C2) was undertaken in the multi-ancestral UK JSLE Cohort including participants of Caucasian (n = 151, 48.8 %), Asian (n = 108, 35.0 %) and African/Caribbean (n = 50, 16.2 %) descent. Considering ancestral variation, clinical associations were tested at the level of alleles (2-field resolution), associated HLA protein sequences (antigen binding domains, 4-field resolution), and extended haplotypes (DRh).

Results

Although important ancestral recombination was reported for HLA-DR2 and -DR3 haplotypes, risk associated with jSLE was conserved at related alleles (DR2h: DRB1∗15:01, DQA∗01:02, DQB1∗06:02; DR3h: C∗07:02 [Asian], B∗08:01, C2 rs9332730 [Asian], DRB1∗03:01). HLA-DR7 haplotypes (DRB1∗07:01, OR = 0.44, 95 % CI:0.27-0.72, p = 0.0004; DQA1∗02:01, OR = 0.34, 95 % CI:0.21-0.56, p = 1.8 × 10-6) protect Asians from jSLE development. Among 23 clinical variables recorded, the main association was found between low levels of complement C4 in Caucasian carriers of HLA-DR3h. This was not the case in Asians due to recombination with HLA-C∗07:02 and integration of the C2 rs9332730 minor allele. Low C4 serum levels associated with HLA-DQA1∗01:02 (DR2h) in Caucasians after excluding HLA-DR3h carriers from the analysis. An association between low white blood cell counts and HLA-A∗03:01P was observed across ancestries.

Conclusion

Genetic variation in the HLA cluster associates with organ domain involvement (hematological) and complement levels in jSLE. Lupus-associated HLA haplotypes vary between ancestral groups, underscoring the importance of multi-ancestral approaches to genetic studies in SLE and other autoimmune/inflammatory diseases.

Long-Term Health Outcomes of Individuals With Pseudodeficiency Alleles in IDUA May Inform Newborn Screening Practices for Mucopolysaccharidosis Type I

Mucopolysaccharidosis type I (MPS I), a lysosomal disorder caused by variants in IDUA, was added to the Recommended Uniform Screening Panel for newborn screening in 2016. Positive screening results for MPS I are commonly due to variants known as "pseudodeficiency alleles," which decrease in vitro alpha-L-iduronidase enzyme activity but are thought to provide sufficient in vivo activity. Despite the historic assumption that these variants are biologically benign, the possibility that they could give rise to complex, multigenic, or attenuated phenotypes has not been systemically evaluated in adults. We completed a retrospective matched cohort study using a hospital-based biorepository with data from 65,309 participants, we identified 1803 individuals harboring homozygous IDUA pseudodeficiency alleles. Using electronic medical records (EMR), we compared the prevalence of features of MPS I in participants with homozygous pseudodeficiency alleles to a cohort of matched control participants. We found no clinically relevant significant differences between cases and controls nor genotype-phenotype associations across four alleles. These findings provide empiric support that adults with homozygous IDUA pseudodeficiency alleles are unlikely to develop mild symptoms of disease compared with controls. This study provides a proof-of-concept model for other nonclassical disease variants related to other inherited metabolic disorders, which is necessary as newborn screening expands.

Mutation Screening of ATXN1, ATXN2, and ATXN3 in Amyotrophic Lateral Sclerosis

Emerging evidence suggests potential disease modifying roles of ATXN1, ATXN2, and ATXN3 in amyotrophic lateral sclerosis (ALS). We aimed to provide a comprehensive variants profile of the ATXN1, ATXN2, and ATXN3 genes and examine the association of these variants with the risk and clinical characteristics of ALS. We screened and analyzed the rare variants in a cohort of 2220 ALS patients from Southwest China, using controls from the Genome Aggregation Database (gnomAD) and the China Metabolic Analytics Project (ChinaMAP). The over-representation of rare variants and their association with disease risk in ALS patients were assessed using Fisher's exact test with Bonferroni correction at both allele and gene levels. Kaplan-Meier analysis was employed to explore the relationship between the distribution of variants and survival. A total of 62 eligible rare missense variants were identified, comprising 32 from ATXN1, 21 from ATXN2, and 9 from ATXN3. Allelic association testing revealed a significant enrichment of the ATXN1 (c.2122C > G, p.Leu708Val) variant and the ATXN2 (c.3778C > G, p.Pro1260Ala) variant in ALS. Gene burden analysis indicated that variants in the ATXN1 and ATXN3 genes had a higher burden in ALS. Substantial heterogeneity in survival time was observed among patients carrying different variants within the same gene. However, there were no significant differences in survival between ALS patients grouped by N-terminal or C-terminal distribution. Our results provided a genetic variation profile of ATXN1, ATXN2, and ATXN3 in ALS patients, along with the clinical characteristics of individuals carrying these variations. This information might offer valuable insights for the ongoing ALS disease-modifying treatments.

Biased signaling in drug discovery and precision medicine

Receptors are crucial for converting chemical and environmental signals into cellular responses, making them prime targets in drug discovery, with about 70% of drugs targeting these receptors. Biased signaling, or functional selectivity, has revolutionized drug development by enabling precise modulation of receptor signaling pathways. This concept is more firmly established in G protein-coupled receptor and has now been applied to other receptor types, including ion channels, receptor tyrosine kinases, and nuclear receptors. Advances in structural biology have further refined our understanding of biased signaling. This targeted approach enhances therapeutic efficacy and potentially reduces side effects. Numerous biased drugs have been developed and approved as therapeutics to treat various diseases, demonstrating their significant therapeutic potential. This review provides a comprehensive overview of biased signaling in drug discovery and disease treatment, highlighting recent advancements and exploring the therapeutic potential of these innovative modulators across various diseases.

The recent blooming of therapeutic aptamers

In the dynamic landscape of biomedical research, therapeutic RNA aptamers have recently come to the forefront, showing significant potential in diagnostics and therapeutics. This review aims to raise awareness of aptamer technology within the scientific community by exploring the progress made in the therapeutic field, from the lessons learned in research to the future opportunities and impact that these innovative molecules are increasingly having on society to meet current health needs, i.e. targeted and personalized therapies.

Exploring the connection between RNA splicing and intellectual disability

Intellectual disability (ID) is a broad diagnostic category that encompasses individuals with impaired cognitive ability. While these disorders have heterogeneous causes, recent developments in next-generation sequencing (NGS) are revealing the prevalence of genetic etiologies. In particular, germline mutations in genes that affect RNA splicing are increasingly common causes of ID disorders. Research to elucidate the functional relationship between splicing and neurodevelopment is critical since molecular therapeutics require a nuanced understanding of the pathological mechanism. In this review, we first summarize the trends that have led to the discovery of the RNA splicing-ID relationship, then discuss recent progress and future directions for research surrounding RNA splicing in neurodevelopment. Finally, we speak on how these results may serve as the foundation for burgeoning therapies.

Congenital Titinopathy: Comprehensive Characterization of the Most Severe End of the Disease Spectrum

Congenital titinopathy has recently emerged as one of the most common congenital muscle disorders.

OBJECTIVE

To better understand the presentation and clinical needs of the under-characterized extreme end of the congenital titinopathy severity spectrum.

METHODS

We comprehensively analyzed the clinical, imaging, pathology, autopsy, and genetic findings in 15 severely affected individuals from 11 families.

RESULTS

Prenatal features included hypokinesia or akinesia and growth restriction. Six pregnancies were terminated. Nine infants were born at or near term with severe-to-profound weakness and required resuscitation. Seven died following withdrawal of life support. Two surviving children require ongoing respiratory support. Most cohort members had at least 1 disease-causing variant predicted to result in some near-normal-length titin expression. The exceptions, from 2 unrelated families, had homozygous truncating variants predicted to induce complete nonsense mediated decay. However, subsequent analyses suggested that the causative variant in each family had an additional previously unrecognized impact on splicing likely to result in some near-normal-length titin expression. This impact was confirmed by minigene assay for 1 variant.

INTERPRETATION

This study confirms the clinical variability of congenital titinopathy. Severely affected individuals succumb prenatally/during infancy, whereas others survive into adulthood. It is likely that this variability is because of differences in the amount and/or length of expressed titin. If confirmed, analysis of titin expression could facilitate clinical prediction and increasing expression might be an effective treatment strategy. Our findings also further-support the hypothesis that some near-normal-length titin expression is essential to early prenatal survival. Sometimes expression of normal/near-normal-length titin is due to disease-causing variants having an additional impact on splicing. ANN NEUROL 2025;97:611-628.

Sequence variants in HECTD1 result in a variable neurodevelopmental disorder

Dysregulation of genes encoding the homologous to E6AP C-terminus (HECT) E3 ubiquitin ligases has been linked to cancer and structural birth defects. One member of this family, the HECT-domain-containing protein 1 (HECTD1), mediates developmental pathways, including cell signaling, gene expression, and embryogenesis. Through GeneMatcher, we identified 14 unrelated individuals with 15 different variants in HECTD1 (10 missense, 3 frameshift, 1 nonsense, and 1 splicing variant) with neurodevelopmental disorders (NDDs), including autism, attention-deficit/hyperactivity disorder, and epilepsy. Of these 15 HECTD1 variants, 10 occurred de novo, 3 had unknown inheritance, and 2 were compound heterozygous. While all individuals in this cohort displayed NDDs, no genotype-phenotype correlation was apparent. Conditional knockout of Hectd1 in the neural lineage in mice resulted in microcephaly, severe hippocampal malformations, and complete agenesis of the corpus callosum, supporting a role for Hectd1 in embryonic brain development. Functional studies of select variants in C. elegans revealed dominant effects, including either change-of-function or loss-of-function/haploinsufficient mechanisms, which may explain phenotypic heterogeneity. Significant enrichment of de novo variants in HECTD1 was also shown in an independent cohort of 53,305 published trios with NDDs or congenital heart disease. Thus, our clinical and functional data support a critical requirement of HECTD1 for human brain development.

Isogenic hiPSC models of Turner syndrome development reveal shared roles of inactive X and Y in the human cranial neural crest network

Viable human aneuploidy can be challenging to model in rodents due to syntenic boundaries or primate-specific biology. Human monosomy-X (45,X) causes Turner syndrome (TS), altering craniofacial, skeletal, endocrine, and cardiovascular development, which in contrast remain unaffected in X-monosomic mice. To learn how monosomy-X may impact embryonic development, we turned to 45,X and isogenic euploid human induced pluripotent stem cells (hiPSCs) from male and female mosaic donors. Because the neural crest (NC) is hypothesized to give rise to craniofacial and cardiovascular changes in TS, we assessed differential expression of hiPSC-derived anterior NC cells (NCCs). Across three independent isogenic panels, 45,X NCCs show impaired acquisition of PAX7+SOX10+ markers and disrupted expression of other NCC-specific genes relative to isogenic euploid controls. Additionally, 45,X NCCs increase cholesterol biosynthesis genes while reducing transcripts with 5' terminal oligopyrimidine (TOP) motifs, including those of ribosomal and nuclear-encoded mitochondrial proteins. Such metabolic pathways are also over-represented in weighted co-expression modules that are preserved in monogenic neurocristopathy and reflect 28% of all TS-associated terms of the human phenotype ontology. We demonstrate that 45,X NCCs reduce protein synthesis despite activation of mammalian target of rapamycin (mTOR) but are partially rescued by mild mTOR suppression. Our analysis identifies specific sex-linked genes that are expressed from two copies in euploid males and females alike and qualify as candidate haploinsufficient drivers of TS phenotypes in NC-derived lineages. This study demonstrates that isogenic hiPSC-derived NCC panels representing monosomy-X can serve as powerful models of early NC development in TS and inform new hypotheses toward its etiology.

Modeling SMAD2 Mutations in Induced Pluripotent Stem Cells Provides Insights Into Cardiovascular Disease Pathogenesis

BACKGROUND

SMAD2 is a coregulator that binds a variety of transcription factors in human development. Heterozygous SMAD2 loss-of-function and missense variants are identified in patients with congenital heart disease (CHD) or arterial aneurysms. Mechanisms that cause distinct cardiovascular phenotypes remain unknown. We aimed to define transcriptional and epigenetic effects of SMAD2 variants and their role in CHD. We also assessed the function of SMAD2 missense variants of uncertain significance.

METHODS AND RESULTS

Rare SMAD2 variants (minor allele frequency ≤10-5) were identified in exome sequencing of 11 336 participants with CHD. We constructed isogenic induced pluripotent stem cells with heterozygous or homozygous loss-of-function and missense SMAD2 variants identified in CHD probands. Wild-type and mutant induced pluripotent stem cells were analyzed using bulk RNA sequencing, chromatin accessibility (Assay for Transposase-Accessible Chromatin With Sequencing), and integrated with published SMAD2/3 chromatin immunoprecipitation data. Cardiomyocyte differentiation and contractility were evaluated. Thirty participants with CHD had heterozygous loss-of-function or missense SMAD2 variants. SMAD2 haploinsufficiency altered chromatin accessibility at promoters and dysregulated expression of 385 SMAD regulated genes, including 10 CHD-associated genes. Motifs enriched in differential Assay for Transposase-Accessible Chromatin peaks predicted that SMAD2 haploinsufficiency disrupts interactions with transcription factors NANOG (homeobox protein NANOG), ETS, TEAD3/4 (transcriptional enhanced associate domain 3/4), CREB1 (cAMP response element binding protein 1), and AP1 (activator protein 1). Compared with SMAD2-haploinsufficient cells, induced pluripotent stem cells with R114C or W274C variants exhibited distinct and shared chromatin accessibility and transcription factor binding changes.

CONCLUSIONS

SMAD2 haploinsufficiency disrupts transcription factor binding and chromatin interactions critical for cardiovascular development. Differences between the molecular consequences of loss-of-function and missense variants likely contribute to phenotypic heterogeneity. These findings indicate opportunities for molecular analyses to improve reclassification of SMAD2 variants of uncertain clinical significance.

Population-Scale Studies of Protein S Abnormalities and Thrombosis

Importance

Clinical decision-making in thrombotic disorders is impeded by long-standing uncertainty regarding the magnitude of venous and arterial thrombosis risk associated with low protein S. Population-scale multiomic datasets offer an unprecedented opportunity to answer questions regarding the epidemiology and clinical impacts of protein S deficiency.

Objective

To evaluate the risk associated with protein S deficiency across multiple thrombosis phenotypes.

Design, Setting, and Participants

Cross-sectional study using longitudinal population cohorts derived from the UK Biobank (n = 426 436) and the US National Institutes of Health All of Us (n = 204 006) biorepositories. UK Biobank participants were enrolled in 2006-2010 (last follow-up, May 19, 2020) and underwent whole exome sequencing, with a subset (n = 44 431) having protein S levels measured by high-throughput plasma proteomics. Recruitment for All of Us began in 2017 and is ongoing, with participants receiving germline whole genome sequencing. Both cohorts include individual-level data on demographics, laboratory measurements, and clinical outcomes.

Exposure

Presence of rare germline genetic variants in PROS1, segmented by functional impact score (FIS), an in silico prediction of the probability that a genetic variant will disrupt protein activity.

Main Outcomes and Measures

Firth logistic regression and linear regression modeling were used to evaluate the thrombosis risk associated with low plasma protein S levels and PROS1 variants across a range of FIS ratings.

Results

The UK Biobank cohort was 54.3% female, with a median age of 58.3 (IQR, 50.5-63.7) years at enrollment. Most participants (95.6%) were of European ancestry, and 18 011 had experienced a venous thromboembolism (VTE). In this population cohort, heterozygosity for the highest-risk PROS1 variants with an FIS of 1.0 (nonsense, frameshift, and essential splice site disruptions) was rare (adjusted prevalence, 0.0091% in the UK and 0.0178% in the US) and associated with markedly increased risk of VTE (odds ratio [OR], 14.01; 95% CI, 6.98-27.14; P = 9.09 × 10-11). Plasma proteomics (n = 44 431) demonstrated that carriers of these variants had total protein S levels that were 48.0% of normal (P = .02 compared with noncarriers). In contrast, less damaging missense variants (FIS ≥0.7) occurred more commonly (adjusted prevalence, 0.22% in the UK and 0.20% in the US) and were associated with marginally reduced plasma protein S concentrations and a smaller point estimate for VTE risk (OR, 1.977; 95% CI, 1.552-2.483; P = 1.95 × 10-7). Associations between PROS1 and VTE at both FIS cutoffs were independently validated in the All of Us cohort with similar effect sizes. No association was detected between the presence of coding PROS1 variants and 3 forms of arterial thrombosis: myocardial infarction, peripheral artery disease, and noncardioembolic ischemic stroke. The presence of PROS1 variants correlated poorly with low plasma protein S levels, and protein S deficiency was significantly associated with VTE and peripheral artery disease regardless of PROS1 variant carrier status. The elevated risk of VTE associated with germline loss of function in PROS1 was evident in Kaplan-Meier survival analysis and appeared to persist throughout life (log-rank P = .0005).

Conclusions and Relevance

True inherited loss of function in PROS1 is rare but represents a stronger risk factor for VTE in the general population than previously understood. Acquired, environmental, or trans-acting genetic factors are more likely to cause circulating protein S deficiency than coding variation in PROS1, and low plasma protein S is associated with VTE.

A novel recurrent ARL3 variant c.209G > A p.(Gly70Glu) causes variable non-syndromic dominant retinal dystrophy with defective lipidated protein transport in human retinal stem cell models

Inherited retinal dystrophies (IRDs) are characterized by their high clinical and genetic heterogeneity. Despite significant advances in the identification of genes associated with IRDs, many individuals and families still have not received a definite molecular diagnosis. Here, we performed clinical examinations and conducted genetic testing in five families with IRD. Whole exome sequencing in the five index cases revealed a heterozygous missense variant, c.209G > A, p.(Gly70Glu) in the ARL3 gene (NM_004311.4). A de novo occurrence was demonstrated in one affected individual and autosomal dominant inheritance in nine affected individuals from four families. Their phenotypes displayed variable expressivity, and ranged from rod-cone to cone-rod dystrophy with photophobia. Human induced pluripotent stem cells (hiPSCs) were generated from dermal fibroblasts from the individual with the de novo ARL3 variant and were differentiated to retinal pigment epithelium cells (RPE) and retinal organoids. Immunofluorescence analyses in these models showed decreased INPP5E localization within the cilia of RPE and connecting cilia of retinal organoids, as well as reduced PDE6⍺ in the organoid outer segments, suggesting that the p.(Gly70Glu) variant causes IRD by defective lipidated protein transport in photoreceptors and/or RPE. This is the first study of ARL3 dysfunction in human retinal cells, highlighting its importance for retinal homeostasis, as well as a variability in the clinical presentation of ARL3-associated IRD.

Understanding the genetic epidemiology of hereditary breast cancer in India using whole genome data from 1029 healthy individuals

Breast cancer is the most highly reported cancer in India. Genetic testing could help tackle the increasing cancer burden by enabling carriers obtain early diagnosis through increased surveillance, and help guide treatment. However, accurate interpretation of variant pathogenicity must be established in a population-specific manner to ensure effective use of genetic testing. Here we query IndiGen data obtained from sequencing 1029 Indian individuals, and perform variant classification of all reported BRCA variants using gold-standard ACMG & AMP guidelines to establish disease epidemiology. To address the high number of VUS variants thus obtained, we further utilize the brca-NOVUS ML tool to obtain pathogenicity predictions in a manner close to ACMG guidelines at scale. Through the manual application of ACMG & AMP guidelines, we determined the genetic prevalence to be the following: 1 in 342 carriers of BRCA1, and 1 in 256 carriers of BRCA2 pathogenic/likely pathogenic variants bear a significant lifetime risk of developing breast / ovarian cancer in India. The high population prevalence and unique variant landscape emphasizes a need for population-scale studies of causative variants to enable effective screening. We advise cautious clinical interpretation, given incomplete penetrance and other complex factors that result in cancer disease. To the best of our knowledge, this is the first and most comprehensive population-scale genetic epidemiological study of BRCA-linked breast cancer variants reported from India.

Coincidence of acral peeling skin syndrome and Nagashima-type palmoplantar keratosis in a Japanese pedigree with acral skin peeling

Acral peeling skin syndrome (APSS; MIM 609796) is a rare genodermatosis characterized by painless focal cutaneous exfoliation of the dorsal hands and feet, typically displaying autosomal recessive inheritance. While cases associated with a founder mutation in TGM5 are relatively common in European Caucasian populations, no APSS cases have been reported from Japan or other East Asian countries. In contrast, Nagashima-type palmoplantar keratosis (NPPK; MIM 615598), caused by variants in SERPINB7, is relatively common in East Asia due to founder mutations. We describe a 27-year-old Japanese woman with spontaneous focal cutaneous exfoliation of the dorsal hand following prolonged glove use, indicative of APSS. Histopathological examination revealed a cleft between the stratum corneum and stratum granulosum and within the horny layer of the epidermis, supporting this diagnosis. However, her mother and maternal uncle exhibited similar symptoms, and there was no reported consanguinity in the patient's parents or grandparents, prompting suspicion of an autosomal dominant genodermatosis. Whole-genome sequencing (WGS) revealed compound heterozygous variants in TGM5 (c.1037G>A and c.684 + 1G>A) as suspected causative variants in the patient, leading to an APSS diagnosis, the first reported in East Asia. On the other hand, her mother and maternal uncle were diagnosed with NPPK due to compound heterozygous pathogenic variants in SERPINB7 (c.796C>T and c.455-1G>A). This case highlights the complexity of diagnosing skin disorders when multiple genodermatoses with similar phenotypes exist within a pedigree. Comprehensive genetic analyses, such as whole-exome sequencing and WGS, are invaluable for identifying causative variants in such complex cases.

Multilevel Analysis of MYC and BCL2 Aberrations in Diffuse Large B-Cell Lymphoma: Identifying a High-Risk Patient Subgroup Across Cell-of-Origin Using Targeted Sequencing

INTRODUCTION

Diffuse large B-cell lymphoma (DLBCL) exhibits striking clinical and biological heterogeneity. Recent studies have identified new subgroups within germinal center B-cell like (GCB) DLBCL, associated with inferior prognosis, irrespective of MYC and BCL2 translocations. We explored the existence of such a DLBCL high-risk subgroup, based on multilevel aberrations, especially focusing on MYC and BCL2.

METHODS

Tissue samples from 111 DLBCL patients were sequenced with a 90-gene lymphoma panel, followed by integrative analyses combining sequencing data, immunohistochemistry, fluorescent in situ hybridization, and clinical data.

RESULTS

We identified a high-risk subgroup in DLBCL defined by: dual immunohistochemical MYC and BCL2 expression (DEL), concurrent MYC and BCL2 translocations (DHL-BCL2), mutations in MYC, CXCR4, or both, and/or BCL2 amplification. The high-risk subgroup constituted 41% of the cohort and included DHL-BCL2, DEL, a GCB subgroup likely representing the recently described GCB subgroups, and a subset of non-GCB patients. In multivariate analysis, high-risk features provided independent predictive value from age and IPI. The 5-year overall survival was 36% in high-risk patients, compared to 76% in non-high-risk patients.

CONCLUSION

We identified a distinct high-risk DLBCL subgroup, characterized by MYC and BCL2 aberrations, beyond conventional DHL-BCL2 and DEL, and irrespective of cell-of-origin, thereby expanding the poor-prognosis group.

Toward understanding the role of genomic repeat elements in neurodegenerative diseases

Neurodegenerative diseases cause great medical and economic burdens for both patients and society; however, the complex molecular mechanisms thereof are not yet well understood. With the development of high-coverage sequencing technology, researchers have started to notice that genomic repeat regions, previously neglected in search of disease culprits, are active contributors to multiple neurodegenerative diseases. In this review, we describe the association between repeat element variants and multiple degenerative diseases through genome-wide association studies and targeted sequencing. We discuss the identification of disease-relevant repeat element variants, further powered by the advancement of long-read sequencing technologies and their related tools, and summarize recent findings in the molecular mechanisms of repeat element variants in brain degeneration, such as those causing transcriptional silencing or RNA-mediated gain of toxic function. Furthermore, we describe how in silico predictions using innovative computational models, such as deep learning language models, could enhance and accelerate our understanding of the functional impact of repeat element variants. Finally, we discuss future directions to advance current findings for a better understanding of neurodegenerative diseases and the clinical applications of genomic repeat elements.

The population-specific Thr44Met OCT3 coding variant affects metformin pharmacokinetics with subsequent effects on insulin sensitivity in C57Bl/6J mice

AIMS/HYPOTHESIS

Metformin is an important first-line treatment for type 2 diabetes and acts by increasing the body's ability to dispose of glucose. Metformin's efficacy can be affected by genetic variants in the transporters that regulate its uptake into cells. The SLC22A3 gene (also known as EMT; EMTH; OCT3) codes for organic cation transporter 3 (OCT3), which is a broad-specificity cation transporter that also transports metformin. Most SLC22A3 variants reduce the rate of metformin transport but the rs8187715 variant (p.Thr44Met) is reported to increase uptake of metformin in vitro. However, the impact of this on in vivo metformin transport and efficacy is unknown. Very few carriers of this variant have been reported globally, but, notably, all were of Pacific Island descent. Therefore, this study aims to understand the prevalence of this variant in Polynesian peoples (Māori and Pacific peoples) and to understand its impact on metformin transport and efficacy in vivo.

METHODS

rs8187715 was genotyped in 310 individuals with Māori and Pacific ancestry recruited in Aotearoa New Zealand. To study this variant in a physiological context, an orthologous knockin mouse model with C57BL/6J background was used. Pharmacokinetic analysis compared uptake rate of metformin into tissues. Plasma growth/differentiation factor 15 (GDF-15) was also measured as a marker of metformin efficacy. Glucose and insulin tolerance was assessed after acute or sustained metformin treatment in knockin and wild-type control mice to examine the impact of the variant on metformin's glycaemic control.

RESULTS

The minor allele frequency of this variant in the Māori and Pacific participants was 15.4%. There was no association of the variant with common metabolic parameters including diabetes status, BMI, blood pressure, lipids, or blood glucose and HbA1c. However, in the orthologous knockin mouse model, the rate of metformin uptake into the blood and tissues was increased. Acute metformin dosing increased insulin sensitivity in variant knockin mice but this effect was lost after longer-term metformin treatment. Metformin's effects on GDF-15 levels were also lost in variant knockin mice with longer-term metformin treatment.

CONCLUSIONS/INTERPRETATION

These data provide evidence that the SLC22A3 rs8187715 variant accelerates metformin uptake rate in vivo. While this acutely improves insulin sensitivity, there was no increased effect of metformin with longer-term dosing. Thus, our finding of a high prevalence of this variant specifically in Māori and Pacific peoples identifies it as a potential population-specific pharmacogenetic marker with potential to guide metformin therapy in these peoples.

Management of ultrarare inherited arrhythmia syndromes

Ultrarare inherited arrhythmia syndromes are increasingly diagnosed as a result of increased awareness as well as increased availability and reduced cost of genetic testing. Yet by definition, their rarity and heterogeneous expression make development of evidence-based management strategies more challenging, typically employing strategies garnered from similar genetic cardiac disorders. For the most part, reliance on anecdotal experiences, expert opinion, and small retrospective cohort studies is the only means to diagnose and to treat these patients. Here we review the management of specific ultrarare inherited arrhythmic syndromes together with the genetic and molecular basis, which will become increasingly important with the development of targeted therapies to correct the biologic basis of these disorders.

MixDeR: A SNP mixture deconvolution workflow for forensic genetic genealogy

The generation of forensic DNA profiles consisting of single nucleotide polymorphisms (SNPs) is now being facilitated by wider adoption of next-generation sequencing (NGS) methods in casework laboratories. At the same time, and in part because of this advance, there is an intense focus on the generation of SNP profiles from evidentiary specimens for so-called forensic or investigative genetic genealogy (FGG or IGG) applications. However, FGG methods are constrained by the algorithms for genealogical database searches, which were designed for use with single-source profiles, and the fact that many forensic samples are mixtures. To enable the use of two-person mixtures for FGG, we developed a workflow, MixDeR, for the deconvolution of mixed SNP profiles. MixDeR, a flexible and easy to use R package and Shiny app, processes ForenSeq Kintelligence® (QIAGEN, Inc.) SNP genotyping results and directs deconvolution of the profiles in EuroForMix (EFM). MixDeR then filters the EFM outputs to produce inferred single-source genotypes in reports formatted for use with GEDmatch® PRO. An optional MixDeR output includes metrics that assist with testing and validation of the workflow. As the Shiny app provides a graphical user interface and the software is designed to be run offline, MixDeR should be suitable for use by any laboratory developing FGG capabilities, no matter their bioinformatic resources or expertise.

Genome-wide association testing beyond SNPs

Decades of genetic association testing in human cohorts have provided important insights into the genetic architecture and biological underpinnings of complex traits and diseases. However, for certain traits, genome-wide association studies (GWAS) for common SNPs are approaching signal saturation, which underscores the need to explore other types of genetic variation to understand the genetic basis of traits and diseases. Copy number variation (CNV) is an important source of heritability that is well known to functionally affect human traits. Recent technological and computational advances enable the large-scale, genome-wide evaluation of CNVs, with implications for downstream applications such as polygenic risk scoring and drug target identification. Here, we review the current state of CNV-GWAS, discuss current limitations in resource infrastructure that need to be overcome to enable the wider uptake of CNV-GWAS results, highlight emerging opportunities and suggest guidelines and standards for future GWAS for genetic variation beyond SNPs at scale.

Clinical Characterization, Natural History, and Detailed Phenotyping of NMNAT1-Associated Leber Congenital Amaurosis

PURPOSE

To characterize the clinical phenotype and disease progression in patients with NMNAT1-associated Leber congenital amaurosis (LCA) within the Korean population.

DESIGN

Retrospective, observational case series.

SUBJECTS

Fourteen patients with LCA with biallelic variants of NMNAT1 at a single tertiary referral center.

METHODS

Electronic medical records were reviewed for medical history, ophthalmic examinations, and molecular diagnoses, both cross-sectionally and longitudinally.

MAIN OUTCOME MEASURES

Ophthalmic examination findings were evaluated and retinal phenotypic characteristics were assessed using multimodal imaging.

RESULTS

All patients exhibited early-onset, rapidly progressive bilateral retinal degeneration with pronounced central involvement. The condition was characterized by multiple atrophic lesions that coalesced into a large central retinal scar by age 2. The condition stabilized around 4 years of age. Fluorescein angiography demonstrated central hypofluorescence with visible choroidal vasculature. Optical coherence tomography showed significant retinal thinning, outer retinal layer disruption, and retinal pigment epithelial atrophy. Most patients maintained light perception vision or better, with minimal deterioration of visual acuity after the age of 2. All patients were hyperopic and exhibited undetectable electroretinography and visual-evoked potential responses.

CONCLUSIONS

NMNAT1-associated LCA is characterized by severe, early-onset retinal degeneration with rapid progression, followed by stabilization. This distinct temporal pattern of disease progression suggests a potential therapeutic window in early childhood, emphasizing the importance of early diagnosis and regular monitoring for potential interventions.

Haplotype-Resolved Genotyping and Association Analysis of 1,020 β-Thalassemia Patients by Targeted Long-Read Sequencing

Despite the well-documented mutation spectra of β-thalassemia, the genetic variants and haplotypes of globin gene clusters modulating its clinical heterogeneity remain incompletely illustrated. Here, a targeted long-read sequencing (T-LRS) is demonstrated to capture 20 genes/loci in 1,020 β-thalassemia patients. This panel permits not only identification of thalassemia mutations at 100% of sensitivity and specificity, but also detection of rare structural variants (SVs) and single nucleotide variants (SNVs) in modifier genes/loci. The highly homologous regions of α-/β-globin gene clusters are then phased and 3 novel haplotypes in HBG1/HBG2 region are reported in this population of β-thalassemia patients. Furthermore, one of the haplotypes is associated with ameliorated symptoms of β-thalassemia. Similarly, 5 major haplotypes are identified in HBA1/HBA2 homologous region while one of them is found highly linked with deletional α-thalassemia mutations. Finally, rare mutations in erythroid transcription factors in DNMT1 and KLF1 associated with increased expression of fetal hemoglobin and reduced transfusion dependencies are identified. This study presents the largest T-LRS study for β-thalassemia patients to date, facilitating precise clinical diagnosis and haplotype phasing of globin gene clusters.

A Prevalent TMEM260 Deletion Causes Conotruncal Heart Defects, Including Truncus Arteriosus

Conotruncal heart defects are severe congenital malformations of the outflow tract, including truncus arteriosus (TA) and double-outlet right ventricle (DORV). TA is a severe congenital heart disease (CHD) in which the main arterial outflow tract of the heart fails to separate. We recently reported TMEM260 (NM_017799.4), c.1617del (p.Trp539Cysfs*9), as a major cause of TA in the Japanese population (TMEM260 Keio-Tohoku variant) comparable to the prevalence of the 22q11.2 deletion syndrome, which accounts for 12%-35% of TA. However, no other major causes of TA have not been identified. Here, we report a family that included a TA patient and a DORV patient, harboring the compound heterozygous variants of TMEM260, a 7066-bp deletion encompassing exons 6-7 and c.1393C > T, p.(Gln465*). The allele frequency of the 7066-bp deletion was particularly high in the Japanese population (0.17%). Based on the allele frequency of this deletion and c.1617del (0.36%) in the Japanese population, TMEM260 variants might be associated with more than half of the Japanese patients with TA. This study showed that TMEM260 pathogenic variants might be the most common cause of TA in the Japanese population and could explain the wide spectrum of phenotypes associated with TMEM260-related CHD, including DORV, demonstrating the usefulness of genetic testing in Japanese patients with TA.

Protein family FAM241 in human and mouse

FAM241B was isolated in a genome-wide inactivation screen for generation of enlarged lysosomes. FAM241B and FAM241A comprise protein family FAM241 encoding proteins of 121 and 132 amino acid residues, respectively. The proteins exhibit 25% amino acid sequence identity and contain a domain of unknown function (DUF4605; pfam15378) that is conserved from primitive multicellular eukaryotes through vertebrates. Phylogenetic comparison indicates that duplication of the ancestral FAM241B gene occurred prior to the origin of fish. FAM241B has been deleted from the avian lineage. Fam241a and Fam241b are widely expressed in mouse tissues. Experimental knockout of mouse Fam241a, Fam241b, and the double knockout, did not generate a visible phenotype. Knockout of Fam241A and Fam241B did not exacerbate the phenotype of FIG4 null mice. RNAseq of brain RNA from double knockout mice detected reduced expression of several genes including Arke1e1 and RnaseL. The human variant p.Val115Gly in FAM241B was identified in a patient with developmental delay. Lysosome morphology in patient-derived fibroblasts was normal. In previous studies, FAM241A and FAM241B appeared to co-localize with proteins of the endoplasmic reticulum. The molecular function of this ancient protein family remains to be determined.

Translational strategies to uncover the etiology of congenital anomalies of the kidney and urinary tract

While up to 50% of children requiring kidney replacement therapy have congenital anomalies of the kidney and urinary tract (CAKUT), they represent only a fraction of the total patient population with CAKUT. The extreme variability in clinical outcome underlines the fundamental need to devise personalized clinical management strategies for individuals with CAKUT. Better understanding of the pathophysiology of abnormal kidney and urinary tract development provides a framework for precise diagnoses and prognostication of patients, the identification of biomarkers and disease modifiers, and, thus, the development of personalized strategies for treatment. In this review, we provide a state-of-the-art overview of the currently known genetic causes, including rare variants in kidney and urinary tract development genes, genomic disorders, and common variants that have been attributed to CAKUT. Furthermore, we discuss the impact of environmental factors and their interactions with developmental genes in kidney and urinary tract malformations. Finally, we present multi-angle translational modalities to validate candidate genes and environmental factors and shed light on future strategies to better understand the molecular underpinnings of CAKUT.

Clinical and genetic characteristics and natural history of Finnish families with familial exudative vitreoretinopathy due to pathogenic FZD4 variants

PURPOSE

To report clinical and genetic characteristics of familial exudative vitreoretinopathy (FEVR) in the Finnish population.

METHODS

Detailed clinical and genetic data of 35 individuals with heterozygous pathogenic variants in FZD4 were gathered and analysed.

RESULTS

Thirty-two individuals with FZD4 c.313A>G variant and three individuals with FZD4 c.40_49del were included in the study. The clinical phenotype was variable even among family members with the same FZD4 variant. Only 34% (N = 12/35) of variant-positive individuals had been clinically diagnosed with FEVR. The median age of the onset of symptoms was 2.3 years, ranging between 0 to 25 years. Median visual acuity was 0.1 logMAR (0.8 Snellen decimal), ranging between light perception and -0.1 logMAR (1.25 Snellen decimal). Most (N = 33/35, 94%) were classified as not visually impaired. Despite unilateral visual loss present in some, they did not meet the criteria of visual impairment according to the WHO classification. Two study patients (N = 2/35, 6%) had severe visual impairment. The most common FEVR stage in study patient's eyes (N = 28/70 eyes, 40%) was FEVR stage 1, that is, avascular periphery or abnormal vascularisation. Most of FZD4-variant-positive study patient's eyes (N = 31/50 eyes, 62%) were myopic. Two individuals presented with persistent hyperplastic primary vitreous expanding the phenotypic spectrum of FEVR. Shared haplotypes extending approximately 0.9 Mb around the recurrent FZD4 c.313A>G variant were identified.

CONCLUSION

Most study patients were unaffected or had mild clinical manifestations by FEVR. Myopia seemed to be overly common in FZD4-variant-positive individuals.

Evaluating reproductive carrier screening using biotinidase deficiency as a model: Variants identified, variant rates, and management

PURPOSE

To review biotinidase gene (BTD) variants identified in a large, diverse, reproductive carrier screening (RCS) cohort and outline management of heterozygotes with pathogenic or likely pathogenic (P/LP) variants.

METHODS

This retrospective observational study included samples tested from January 2020 to September 2022 in a 274-gene panel. The study involved females aged 18 to 55 years. Screening was performed using next-generation sequencing covering exons and 10 base-pair flanking introns. The heterozygote frequency was calculated for P/LP variants for the entire population and individual racial/ethnic groups.

RESULTS

Of the 91,637 women tested, 5625 (6.1%) had a P/LP variant in BTD. NM_000060.4:c.1330G>C p.(Asp444His) (referred to as D444H or D424H) alone, or in combination with another variant, accounted for 5193 (92.3%) of the positive tests. P/LP heterozygote rates differed between racial and ethnic groups. We ascertained 7 novel P/LP variants not previously recorded in databases.

CONCLUSION

The BTD P/LP variants identified through RCS were substantially compatible with those found through positive newborn screening. Therefore, RCS provides a potential for earlier diagnosis. We observed significant differences in P/LP heterozygote rates for biotinidase deficiency among different racial and ethnic groups. Most reported variants can be interpreted without requiring determination of serum biotinidase activity.

Genetic variants and phenotypic data curated for the CAGI6 intellectual disability panel challenge

Neurodevelopmental disorders (NDDs) are common conditions including clinically diverse and genetically heterogeneous diseases, such as intellectual disability, autism spectrum disorders, and epilepsy. The intricate genetic underpinnings of NDDs pose a formidable challenge, given their multifaceted genetic architecture and heterogeneous clinical presentations. This work delves into the intricate interplay between genetic variants and phenotypic manifestations in neurodevelopmental disorders, presenting a dataset curated for the Critical Assessment of Genome Interpretation (CAGI6) ID Panel Challenge. The CAGI6 competition serves as a platform for evaluating the efficacy of computational methods in predicting phenotypic outcomes from genetic data. In this study, a targeted gene panel sequencing has been used to investigate the genetic causes of NDDs in a cohort of 415 paediatric patients. We identified 60 pathogenic and 49 likely pathogenic variants in 102 individuals that accounted for 25% of NDD cases in the cohort. The most mutated genes were ANKRD11, MECP2, ARID1B, ASH1L, CHD8, KDM5C, MED12 and PTCHD1 The majority of pathogenic variants were de novo, with some inherited from mildly affected parents. Loss-of-function variants were the most common type of pathogenic variant. In silico analysis tools were used to assess the potential impact of variants on splicing and structural/functional effects of missense variants. The study highlights the challenges in variant interpretation especially in cases with atypical phenotypic manifestations. Overall, this study provides valuable insights into the genetic causes of NDDs and emphasises the importance of understanding the underlying genetic factors for accurate diagnosis, and intervention development in neurodevelopmental conditions.

Discordance between a deep learning model and clinical-grade variant pathogenicity classification in a rare disease cohort

Genetic testing is essential for diagnosing and managing clinical conditions, particularly rare Mendelian diseases. Although efforts to identify rare phenotype-associated variants have focused on protein-truncating variants, interpreting missense variants remains challenging. Deep learning algorithms excel in various biomedical tasks1,2, yet distinguishing pathogenic from benign missense variants remains elusive3-5. Our investigation of AlphaMissense (AM)5, a deep learning tool for predicting the potential functional impact of missense variants and assessing gene essentiality, reveals limitations in identifying pathogenic missense variants over 45 rare diseases, including very early onset inflammatory bowel disease. For the expert-curated pathogenic variants identified in our cohort, AM's precision was 32.9%, and recall was 57.6%. Notably, AM struggles to evaluate pathogenicity in intrinsically disordered regions (IDRs), resulting in unreliable gene-level essentiality scores for genes containing IDRs. This observation underscores ongoing challenges in clinical genetics, highlighting the need for continued refinement of computational methods in variant pathogenicity prediction.

Interpreting mammalian synonymous site conservation in light of the unwanted transcript hypothesis

Mammalian genomes are biased towards GC bases at third codon positions, likely due to a GC-biased ancestral genome and the selectively neutral recombination-related process of GC-biased gene conversion. The unwanted transcript hypothesis posits that this high GC content at synonymous sites may be beneficial for protecting against spurious transcripts, particularly in species with low effective population sizes. Utilising a 240 placental mammal genome alignment and single-base resolution conservation scores, we interpret sequence conservation at mammalian four-fold degenerate sites in this context and find evidence in support of the unwanted transcript hypothesis, including a strong GC bias, high conservation at sites relating to exon splicing, less human genetic variation at conserved four-fold degenerate sites, and conservation of sites important for epigenetic regulation of developmental genes. Additionally, we show that high conservation of four-fold degenerate sites in essential developmental genes, including homeobox genes, likely relates to the low mutation rates experienced by these genes.

Structure of the MUC5AC VWD3 assembly responsible for the formation of net-like mucin polymers

Gel-forming mucins MUC5AC and MUC5B constitute the main structural component of the mucus in the respiratory system. Secreted mucins interact specifically with each other and other molecules giving mucus specific properties. We determined the cryoEM structures of the wild type D3 assembly of the human MUC5AC mucin and the structural single nucleotide polymorphisms (SNP) variants Arg996Gln and Arg1201Trp that affect intermolecular interactions. Our structures explain the MUC5AC N-terminal non-covalent oligomerization after secretion. The D3 assembly forms covalent dimers that can appear in two alternative conformations, open and closed, where the closed conformation dimers interact through an arginine-rich loop in the TIL3 domain to form tetramers. Our study provides a model to explain MUC5AC net-like structures and how the two SNPs will affect mucus organization, something that might affect lung and other diseases.

The Biorepository and Integrative Genomics resource for inclusive genomics: insights from a diverse pediatric and admixed cohort

The Biorepository and Integrative Genomics (BIG) Initiative in Tennessee has developed a pioneering resource to address gaps in genomic research by linking genomic, phenotypic, and environmental data from a diverse Mid-South population, including underrepresented groups. We analyzed 13,152 exomes from BIG and found significant genetic diversity, with 50% of participants inferred to have non-European or several types of admixed ancestry. Ancestry within the BIG cohort is stratified, with distinct geographic and demographic patterns, as African ancestry is more common in urban areas, while European ancestry is more common in suburban regions. We observe ancestry-specific rates of novel genetic variants, which are enriched for functional or clinical relevance. Disease prevalence analysis linked ancestry and environmental factors, showing higher odds ratios for asthma and obesity in minority groups, particularly in the urban area. Finally, we observe discrepancies between self-reported race and genetic ancestry, with related individuals self-identifying in differing racial categories. These findings underscore the limitations of race as a biomedical variable. BIG has proven to be an effective model for community-centered precision medicine. We integrated genomics education, and fostered great trust among the contributing communities. Future goals include cohort expansion, and enhanced genomic analysis, to ensure equitable healthcare outcomes.

Rare disease gene association discovery in the 100,000 Genomes Project

Up to 80% of rare disease patients remain undiagnosed after genomic sequencing1, with many probably involving pathogenic variants in yet to be discovered disease-gene associations. To search for such associations, we developed a rare variant gene burden analytical framework for Mendelian diseases, and applied it to protein-coding variants from whole-genome sequencing of 34,851 cases and their family members recruited to the 100,000 Genomes Project2. A total of 141 new associations were identified, including five for which independent disease-gene evidence was recently published. Following in silico triaging and clinical expert review, 69 associations were prioritized, of which 30 could be linked to existing experimental evidence. The five associations with strongest overall genetic and experimental evidence were monogenic diabetes with the known β cell regulator3,4 UNC13A, schizophrenia with GPR17, epilepsy with RBFOX3, Charcot-Marie-Tooth disease with ARPC3 and anterior segment ocular abnormalities with POMK. Further confirmation of these and other associations could lead to numerous diagnoses, highlighting the clinical impact of large-scale statistical approaches to rare disease-gene association discovery.

Integrative analysis of KCNQ1 variants reveals molecular mechanisms of type 1 long QT syndrome pathogenesis

Loss-of-function (LOF) pathogenic variants in KCNQ1 encoding a cardiac potassium channel predispose to sudden cardiac death in type 1 congenital long QT syndrome (LQT1). To determine the spectrum of molecular mechanisms responsible for this life-threatening condition, we used an integrative approach to determine the biophysical, functional, and trafficking properties of 61 KCNQ1 variants distributed throughout all domains of the channel. Impaired trafficking to the plasma membrane was the most common cause of LOF across all channel domains, often but not always coinciding with protein instability. However, many LOF variants, particularly in transmembrane domains, trafficked normally, but when coexpressed with KCNE1 exhibited impaired conductance, altered voltage dependence, or abnormal gating kinetics, highlighting diverse pathogenic mechanisms. This indicates a need for personalized treatment approaches for LQT1. Use of our data to benchmark variant pathogenicity prediction methods demonstrated that prediction accuracy depends on the exact mechanism of pathogenicity associated with a given variant.

Prevalence of Pathogenic Variants and Eligibility Criteria for Genetic Testing in Patients Who Visit a Memory Clinic

BACKGROUND AND OBJECTIVES

Identifying genetic causes of dementia in patients visiting memory clinics is important for patient care and family planning. Traditional clinical selection criteria for genetic testing may miss carriers of pathogenic variants in dementia-related genes. This study aimed identify how many carriers we are missing and to optimize criteria for selecting patients for genetic counseling in memory clinics.

METHODS

In this clinical cohort study, we retrospectively genetically tested patients during 2.5 years (2010-2012) visiting the Alzheimer Center Amsterdam, a specialized memory clinic. Genetic tests consisted of a 54-gene dementia panel, focusing on Class IV/V variants per American College of Medical Genetics and Genomics guidelines, including APP duplications and the C9ORF72 repeat expansion. We determined the prevalence of pathogenic variants and propose new eligibility criteria for genetic testing in memory clinics. The eligibility criteria were prospectively applied for 1 year (2021-2022), and results were compared with the retrospective cohort.

RESULTS

Genetic tests were retrospectively performed in in 1,022 of 1,138 patients (90%) who consecutively visited the memory clinic. Among these, 1,022 patients analyzed (mean age 62.1 ± 8.9 years; 40.4% were female), 34 pathogenic variant carriers were identified (3.3%), with 24 being symptomatic. Previous clinical criteria would have identified only 15 carriers (44% of all carriers, 65% of symptomatic carriers). The proposed criteria increased identification to 22 carriers (62.5% of all carriers, 91% of symptomatic carriers). In the prospective cohort, 148 (28.7%) of 515 patients were eligible for testing under the new criteria. Of the 90 eligible patients who consented to testing, 13 pathogenic carriers were identified, representing a 73% increase compared with the previous criteria.

DISCUSSION

We found that patients who visit a memory clinic and carry a pathogenic genetic variant are often not eligible for genetic testing. The proposed new criteria improve the identification of patients with a genetic cause for their cognitive complaints. In systems without practical or financial barriers to genetic testing, the new criteria can enhance personalized care. In other countries where the health care systems differs and in other genetic ancestry groups, the performance of the criteria may be different.

Whole-genome sequencing analysis identifies rare, large-effect noncoding variants and regulatory regions associated with circulating protein levels

The contribution of rare noncoding genetic variation to common phenotypes is largely unknown, as a result of a historical lack of population-scale whole-genome sequencing data and the difficulty of categorizing noncoding variants into functionally similar groups. To begin addressing these challenges, we performed a cis association analysis using whole-genome sequencing data, consisting of 1.1 billion variants, 123 million noncoding aggregate-based tests and 2,907 circulating protein levels in ~50,000 UK Biobank participants. We identified 604 independent rare noncoding single-variant associations with circulating protein levels. Unlike protein-coding variation, rare noncoding genetic variation was almost as likely to increase or decrease protein levels. Rare noncoding aggregate testing identified 357 conditionally independent associated regions. Of these, 74 (21%) were not detectable by single-variant testing alone. Our findings have important implications for the identification, and role, of rare noncoding genetic variation associated with common human phenotypes, including the importance of testing aggregates of noncoding variants.

PTEN mutations impair CSF dynamics and cortical networks by dysregulating periventricular neural progenitors

Enlargement of the cerebrospinal fluid (CSF)-filled brain ventricles (ventriculomegaly) is a defining feature of congenital hydrocephalus (CH) and an under-recognized concomitant of autism. Here, we show that de novo mutations in the autism risk gene PTEN are among the most frequent monogenic causes of CH and primary ventriculomegaly. Mouse Pten-mutant ventriculomegaly results from aqueductal stenosis due to hyperproliferation of periventricular Nkx2.1+ neural progenitor cells (NPCs) and increased CSF production from hyperplastic choroid plexus. Pten-mutant ventriculomegalic cortices exhibit network dysfunction from increased activity of Nkx2.1+ NPC-derived inhibitory interneurons. Raptor deletion or postnatal everolimus treatment corrects ventriculomegaly, rescues cortical deficits and increases survival by antagonizing mTORC1-dependent Nkx2.1+ NPC pathology. Thus, PTEN mutations concurrently alter CSF dynamics and cortical networks by dysregulating Nkx2.1+ NPCs. These results implicate a nonsurgical treatment for CH, demonstrate a genetic association of ventriculomegaly and ASD, and help explain neurodevelopmental phenotypes refractory to CSF shunting in select individuals with CH.

Digenic impairments of haploinsufficient genes in patients with craniosynostosis

Craniosynostosis (CRS) is characterized by the development of abnormal cranial suture ossification and premature fusion. Despite the identification of several associated genetic disorders, the genetic determinants of CRS remain poorly understood. In this study, we conducted integrative analyses on 225 exomes, comprising 121 CRS probands and 104 parental exomes (52 trios). These analyses encompassed de novo and pathogenic variants, and digenic combinations within haploinsufficient genes harboring rare variants. Our analysis unveils a shared molecular network between genes associated with CRS and those linked to skeletal and neurodevelopmental disorders, with a notable enrichment of deleterious variants within haploinsufficient genes. Additionally, we identified a unique digenic pair (IL6ST and TRPS1) within haploinsufficient genes that was present in 2 patients with nonsyndromic CRS but absent in parents or 1,048 population controls. In vitro experiments provided evidence that the identified missense variants were hypomorphs, and accelerated bone mineralization could result from the additive effects of diminished IL6ST and TRPS1 activities in osteoblasts. Overall, our study underscores the important role of rare variations in haploinsufficient genes and suggests that in a subset of undiagnosed patients, the CRS phenotype may arise from multiple genetic variations.

Excess of rare noncoding variants in several type 2 diabetes candidate genes among Asian Indian families

BACKGROUND

Type 2 diabetes (T2D) etiology is highly complex due to its multiple roots of origin. Polygenic risk scores (PRS) based on genome-wide association studies (GWAS) can partially explain T2D risk. Asian Indian people have up to six times higher risk of developing T2D than European people, and underlying causes of this disparity are unknown.

METHODS

We have performed targeted sequencing of ten T2D GWAS/candidate regions using endogamous Punjabi Sikh families and replication studies using unrelated Sikh people and families from three other Indian endogamous ethnic groups (EEGs).

RESULTS

We detect rare and ultra-rare variants (RVs) in KCNJ11-ABCC8 and HNF4A (MODY genes) cosegregated with late-onset T2D. We also identify RV enrichment in two new genes, SLC38A11 and ANPEP, associated with T2D. Gene-burden analysis reveals the highest RV burden contributed by HNF4A (p = 0.0003), followed by KCNJ11/ABCC8 (p = 0.0061) and SLC38A11 (p = 0.03). Some RVs detected in Sikh people are also found in Agarwals from Jaipur, both from Northern India, but were monomorphic in other two EEGs from South Indian people. Despite carrying a high burden of T2D and RVs, most families have a significantly lower burden of PRS. Functional studies show that an intronic regulatory variant (RV) in ABCC8 affects the binding of Pax4 and NF-kB transcription factors, influencing downstream gene regulation.

CONCLUSIONS

The high burden of T2D in these families may stem from the enrichment of noncoding RVs in a small number of major known genes (including MODY genes) with oligogenic inheritance alongside RVs from genes associated with polygenic susceptibility. These findings highlight the need to conduct deeper evaluations of families from non-European ancestries to identify potential novel therapeutics and implement preventative strategies.

Rare missense variants of the leukocyte common antigen related receptor (LAR) display reduced activity in transcellular adhesion and synapse formation

The leukocyte common antigen related receptor (LAR) is a member of the LAR receptor protein tyrosine phosphatase (RPTP) family of synaptic adhesion molecules that contribute to the proper alignment and specialization of synaptic connections in the mammalian brain. LAR-RPTP members have been genetically associated with neuropsychiatric disorders, but the molecular consequences of genetic perturbations of LAR remain unstudied. Using exome sequencing data from psychiatric patients and controls, we identify rare missense variants of LAR that render the extracellular domain (ECD) unstable and susceptible to proteolytic cleavage. Using recombinant and cellular systems, we describe three variants that cause disruption of the LAR:NGL-3 interaction, which results in loss of transcellular adhesion and synaptogenic effects. Furthermore, we show that overexpression of two of these variants elicit altered morphological phenotypes in an imaging-based morphological profiling assay compared to wild type LAR, suggesting that destabilization of the LAR ECD has broad effects on LAR function. In conclusion, our study identifies three rare, missense variants in LAR that could provide insights into LAR involvement with psychiatric pathobiology.

Mapping the regulatory effects of common and rare non-coding variants across cellular and developmental contexts in the brain and heart

Whole genome sequencing has identified over a billion non-coding variants in humans, while GWAS has revealed the non-coding genome as a significant contributor to disease. However, prioritizing causal common and rare non-coding variants in human disease, and understanding how selective pressures have shaped the non-coding genome, remains a significant challenge. Here, we predicted the effects of 15 million variants with deep learning models trained on single-cell ATAC-seq across 132 cellular contexts in adult and fetal brain and heart, producing nearly two billion context-specific predictions. Using these predictions, we distinguish candidate causal variants underlying human traits and diseases and their context-specific effects. While common variant effects are more cell-type-specific, rare variants exert more cell-type-shared regulatory effects, with selective pressures particularly targeting variants affecting fetal brain neurons. To prioritize de novo mutations with extreme regulatory effects, we developed FLARE, a context-specific functional genomic model of constraint. FLARE outperformed other methods in prioritizing case mutations from autism-affected families near syndromic autism-associated genes; for example, identifying mutation outliers near CNTNAP2 that would be missed by alternative approaches. Overall, our findings demonstrate the potential of integrating single-cell maps with population genetics and deep learning-based variant effect prediction to elucidate mechanisms of development and disease-ultimately, supporting the notion that genetic contributions to neurodevelopmental disorders are predominantly rare.

Conservation assessment of human splice site annotation based on a 470-genome alignment

Despite many improvements over the years, the annotation of the human genome remains imperfect. The use of evolutionarily conserved sequences provides a strategy for selecting a high-confidence subset of the annotation. Using the latest whole genome alignment, we found that splice sites from protein-coding genes in the high-quality MANE annotation are consistently conserved across more than 350 species. We also studied splice sites from the RefSeq, GENCODE, and CHESS databases not present in MANE. In addition, we analyzed the completeness of the alignment with respect to the human genome annotations and described a method that would allow us to fix up to 50of the protein-coding exons. We trained a logistic regression classifier to distinguish between the conservation exhibited by sites from MANE versus sites chosen randomly from neutrally evolving sequences. We found that splice sites classified by our model as well-supported have lower SNP rates and better transcriptomic support. We then computed a subset of transcripts using only "well-supported" splice sites or ones from MANE. This subset is enriched in high-confidence transcripts of the major gene catalogs that appear to be under purifying selection and are more likely to be correct and functionally relevant.

Graphical abstract

Predicting Resistance to Small Molecule Kinase Inhibitors

Drug resistance is a critical challenge in treating diseases like cancer and infectious disease. This study presents a novel computational workflow for predicting on-target resistance mutations to small molecule inhibitors (SMIs). The approach integrates genetic models with alchemical free energy perturbation (FEP+) calculations to identify likely resistance mutations. Specifically, a genetic model, RECODE, leverages cancer-specific mutation patterns to prioritize probable amino acid changes. Physics-based calculations assess the impact of these mutations on protein stability, endogenous substrate binding, and inhibitor binding. We apply this approach retrospectively to gefitinib and osimertinib, two clinical epidermal growth factor receptor (EGFR) inhibitors used to treat nonsmall cell lung cancer (NSCLC). Among hundreds of possible mutations, the pipeline accurately predicted 4 out of 11 and 7 out of 19 known binding site mutations for gefitinib and osimertinib, respectively, including the clinically relevant T790M and C797S resistance mutations. This study demonstrates the potential of integrating genetic models and physics-based calculations to predict SMI resistance mutations. This approach can be applied to other kinases and target classes, potentially enabling the design of next-generation inhibitors with improved durability of response in patients.

Genome-wide prediction of dominant and recessive neurodevelopmental disorder-associated genes

Despite great progress, thousands of neurodevelopmental disorder (NDD) risk genes remain to be discovered. We present a computational approach that accelerates NDD risk gene identification using machine learning. First, we demonstrate that models trained solely on single-cell RNA sequencing data can robustly predict genes implicated in autism spectrum disorder (ASD), developmental and epileptic encephalopathy (DEE), and developmental delay (DD). Notably, we find differences in gene expression patterns of genes with monoallelic and bi-allelic inheritance patterns in the developing human cortex. We then integrate expression data with 300 orthogonal features, including intolerance metrics, protein-protein interaction data, and others, in a semi-supervised machine learning framework (mantis-ml) to train inheritance-specific models for these disorders. The models have high predictive power (area under the receiver operator curves [AUCs]: 0.84-0.95), and the top-ranked genes were up to 2-fold (monoallelic models) and 6-fold (bi-allelic models) more enriched for high-confidence NDD risk genes compared to genic intolerance metrics alone. Additionally, genes ranking in the top decile were 45 to 180 times more likely to have literature support than those in the bottom decile. Collectively, this work provides robust NDD risk gene predictions that can complement large-scale gene discovery efforts and underscores the importance of considering inheritance in gene risk prediction.

Assessing the contribution of rare protein-coding germline variants to prostate cancer risk and severity in 37,184 cases

To assess the contribution of rare coding germline genetic variants to prostate cancer risk and severity, we perform here a meta-analysis of 37,184 prostate cancer cases and 331,329 male controls from five cohorts with germline whole exome or genome sequencing data, and one cohort with imputed array data. At the gene level, our case-control collapsing analysis confirms associations between rare damaging variants in four genes and increased prostate cancer risk: SAMHD1, BRCA2 and ATM at the study-wide significance level (P < 1×10-8), and CHEK2 at the suggestive threshold (P < 2.6×10-6). Our case-only analysis, reveals that rare damaging variants in AOX1 are associated with more aggressive disease (OR = 2.60 [1.75-3.83], P = 1.35×10-6), as well as confirming the role of BRCA2 in determining disease severity. At the single-variant level, our study reveals that a rare missense variant in TERT is associated with substantially reduced prostate cancer risk (OR = 0.13 [0.07-0.25], P = 4.67×10-10), and confirms rare non-synonymous variants in a further three genes associated with reduced risk (ANO7, SPDL1, AR) and in three with increased risk (HOXB13, CHEK2, BIK). Altogether, this work provides deeper insights into the genetic architecture and biological basis of prostate cancer risk and severity, with potential implications for clinical risk prediction and therapeutic strategies.

Alternatively spliced NFKB1 transcripts enriched in Andean Aymara modulate inflammation, HIF and hemoglobin

The molecular basis of increased hemoglobin in Andean Aymara highlanders is unknown. We conducted an integrative analysis of whole-genome-sequencing and granulocytes transcriptomics from Aymara and Europeans in Bolivia to explore genetic basis of the Aymara high hemoglobin. Differentially expressed and spliced genes in Aymaras were associated with inflammatory and hypoxia-related pathways. We identified transcripts with 4th or 5th exon skipping of NFKB1 (AS-NFKB1), key part of NF-kB complex, and their splicing quantitative trait loci; these were increased in Aymaras. AS-NFKB1 transcripts correlated with both transcripts and protein levels of inflammatory and HIF-regulated genes, including hemoglobin. While overexpression of the AS-NFKB1 variant led to increased expression of inflammatory and HIF-targeted genes; under inflammatory stress, NF-kB protein translocation to the nucleus was attenuated, resulting in reduced expression of these genes. Our study reveals AS-NFKB1 splicing events correlating with increased hemoglobin in Aymara and their possible protective mechanisms against excessive inflammation.

Interstitial Cystitis: a phenotype and rare variant exome sequencing study

Interstitial cystitis/bladder pain syndrome (IC/BPS) is a poorly understood and underdiagnosed syndrome of chronic bladder/pelvic pain with urinary frequency and urgency. Though IC/BPS can be hereditary, little is known of its genetic etiology. Using the eMERGE data, we confirmed known phenotypic associations such as gastroesophageal reflux disease and irritable bowel syndrome and detected new associations, including osteoarthrosis/osteoarthritis and Barrett's esophagus. An exome wide ultra-rare variants analysis in 348 IC/BPS and 11,981 controls extended the previously reported association with ATP2C1 and ATP2A2, implicated in Mendelian desquamating skin disorders, but did not provide evidence for other previously proposed pathogenic pathways such as bladder development, nociception or inflammation. Pathway analysis detected new associations with "anaphase-promoting complex-dependent catabolic process", the "regulation of MAPK cascade" and "integrin binding". These findings suggest perturbations in biological networks for epithelial integrity and cell cycle progression in IC/BPS pathogenesis, and provide a roadmap for its future investigation.