The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease

An emerging multi-omic understanding of the genetics of opioid addiction

Opioid misuse, addiction, and associated overdose deaths remain global public health crises. Despite the tremendous need for pharmacological treatments, current options are limited in number, use, and effectiveness. Fundamental leaps forward in our understanding of the biology driving opioid addiction are needed to guide development of more effective medication-assisted therapies. This Review focuses on the omics-identified biological features associated with opioid addiction. Recent GWAS have begun to identify robust genetic associations, including variants in OPRM1, FURIN, and the gene cluster SCAI/PPP6C/RABEPK. An increasing number of omics studies of postmortem human brain tissue examining biological features (e.g., histone modification and gene expression) across different brain regions have identified broad gene dysregulation associated with overdose death among opioid misusers. Drawn together by meta-analysis and multi-omic systems biology, and informed by model organism studies, key biological pathways enriched for opioid addiction-associated genes are emerging, which include specific receptors (e.g., GABAB receptors, GPCR, and Trk) linked to signaling pathways (e.g., Trk, ERK/MAPK, orexin) that are associated with synaptic plasticity and neuronal signaling. Studies leveraging the agnostic discovery power of omics and placing it within the context of functional neurobiology will propel us toward much-needed, field-changing breakthroughs, including identification of actionable targets for drug development to treat this devastating brain disease.

IMPC-based screening revealed that ROBO1 can regulate osteoporosis by inhibiting osteogenic differentiation

Introduction

The utilization of denosumab in treating osteoporosis highlights promising prospects for osteoporosis intervention guided by gene targets. While omics-based research into osteoporosis pathogenesis yields a plethora of potential gene targets for clinical transformation, identifying effective gene targets has posed challenges.

Methods

We first queried the omics data of osteoporosis clinical samples on PubMed, used International Mouse Phenotyping Consortium (IMPC) to screen differentially expressed genes, and conducted preliminary functional verification of candidate genes in human Saos2 cells through osteogenic differentiation and mineralization experiments. We then selected the candidate genes with the most significant effects on osteogenic differentiation and further verified the osteogenic differentiation and mineralization functions in mouse 3T3-E1 and bone marrow mesenchymal stem cells (BMSC). Finally, we used RNA-seq to explore the regulation of osteogenesis by the target gene.

Results

We identified PPP2R2A, RRBP1, HSPB6, SLC22A15, ADAMTS4, ATP8B1, CTNNB1, ROBO1, and EFR3B, which may contribute to osteoporosis. ROBO1 was the most significant regulator of osteogenesis in both human and mouse osteoblast. The inhibitory effect of Robo1 knockdown on osteogenic differentiation may be related to the activation of inflammatory signaling pathways.

Conclusion

Our study provides several novel molecular mechanisms involved in the pathogenesis of osteoporosis. ROBO1 is a potential target for osteoporosis intervention.

In vivo dissection of the mouse tyrosine catabolic pathway with CRISPR-Cas9 identifies modifier genes affecting hereditary tyrosinemia type 1

Hereditary tyrosinemia type 1 is an autosomal recessive disorder caused by mutations (pathogenic variants) in fumarylacetoacetate hydrolase, an enzyme involved in tyrosine degradation. Its loss results in the accumulation of toxic metabolites that mainly affect the liver and kidneys and can lead to severe liver disease and liver cancer. Tyrosinemia type 1 has a global prevalence of approximately 1 in 100,000 births but can reach up to 1 in 1,500 births in some regions of Québec, Canada. Mutating functionally related "modifier' genes (i.e. genes that, when mutated, affect the phenotypic impacts of mutations in other genes) is an emerging strategy for treating human genetic diseases. In vivo somatic genome editing in animal models of these diseases is a powerful means to identify modifier genes and fuel treatment development. In this study, we demonstrate that mutating additional enzymes in the tyrosine catabolic pathway through liver-specific genome editing can relieve or worsen the phenotypic severity of a murine model of tyrosinemia type 1. Neonatal gene delivery using recombinant adeno-associated viral vectors expressing Staphylococcus aureus Cas9 under the control of a liver-specific promoter led to efficient gene disruption and metabolic rewiring of the pathway, with systemic effects that were distinct from the phenotypes observed in whole-body knockout models. Our work illustrates the value of using in vivo genome editing in model organisms to study the direct effects of combining pathological mutations with modifier gene mutations in isogenic settings.

Alzheimer's disease risk gene CD2AP is a dose-sensitive determinant of synaptic structure and plasticity

CD2-Associated protein (CD2AP) is a candidate susceptibility gene for Alzheimer's disease, but its role in the mammalian central nervous system remains largely unknown. We show that CD2AP protein is broadly expressed in the adult mouse brain, including within cortical and hippocampal neurons, where it is detected at pre-synaptic terminals. Deletion of Cd2ap altered dendritic branching and spine density, and impaired ubiquitin-proteasome system activity. Moreover, in mice harboring either one or two copies of a germline Cd2ap null allele, we noted increased paired-pulse facilitation at hippocampal Schaffer-collateral synapses, consistent with a haploinsufficient requirement for pre-synaptic release. Whereas conditional Cd2ap knockout in the brain revealed no gross behavioral deficits in either 3.5- or 12-month-old mice, Cd2ap heterozygous mice demonstrated subtle impairments in discrimination learning using a touchscreen task. Based on unbiased proteomics, partial or complete loss of Cd2ap triggered perturbation of proteins with roles in protein folding, lipid metabolism, proteostasis, and synaptic function. Overall, our results reveal conserved, dose-sensitive requirements for CD2AP in the maintenance of neuronal structure and function, including synaptic homeostasis and plasticity, and inform our understanding of possible cell-type specific mechanisms in Alzheimer's Disease.

Expanding the spectrum of phenotypes for MPDZ: Report of four unrelated families and review of the literature

MPDZ, a gene with diverse functions mediating cell-cell junction interactions, receptor signaling, and binding multivalent scaffold proteins, is associated with a spectrum of clinically heterogeneous phenotypes with biallelic perturbation. Despite its clinical relevance, the mechanistic underpinnings of these variants remain elusive, underscoring the need for extensive case series and functional investigations. In this study, we conducted a systematic review of cases in the literature through two electronic databases following the PRISMA guidelines. We selected nine studies, including 18 patients, with homozygous or compound heterozygous variants in MPDZ and added five patients from four unrelated families with novel MPDZ variants. To evaluate the role of Mpdz on hearing, we analyzed available auditory electrophysiology data from a knockout murine model (Mpdzem1(IMPC)J/em1(IMPC)J) generated by the International Mouse Phenotyping Consortium. Using exome and genome sequencing, we identified three families with compound heterozygous variants, and one family with a homozygous frameshift variant. MPDZ-related disease is clinically heterogenous with hydrocephaly, vision impairment, hearing impairment and cardiovascular disease occurring most frequently. Additionally, we describe two unrelated patients with spasticity, expanding the phenotypic spectrum. Our murine analysis of the Mpdzem1(IMPC)J/em1(IMPC)J allele showed severe hearing impairment. Overall, we expand understanding of MPDZ-related phenotypes and highlight hearing impairment and spasticity among the heterogeneous phenotypes.

Novel molecular, structural and clinical findings in an Italian cohort of congenital cataract

The current genetic diagnostic workup of congenital cataract (CC) is mainly based on NGS panels, whereas exome sequencing (ES) has occasionally been employed. In this multicentre study, we investigated by ES the detection yield, mutational spectrum and genotype-phenotype correlations in a CC cohort recruited between 2020 and mid-2022. The cohort consisted of 67 affected individuals from 51 unrelated families and included both non-syndromic (75%) and syndromic (25%) phenotypes, with extra-CC ocular/visual features present in both groups (48% and 76%, respectively). The functional effect of variants was predicted by 3D modelling and hydropathy properties changes. Variant clustering was used for the in-depth assessment of genotype-phenotype correlations. A diagnostic (pathogenic or likely pathogenic) variant was identified in 19 out of 51 probands/families (~37%). In a further 14 probands/families a candidate variant was identified: in 12 families a VUS was detected, of which 9 were considered plausibly pathogenic (i.e., 4 or 5 points according to ACMG criteria), while in 2 probands ES identified a single variant in an autosomal recessive gene associated with CC. Eighteen probands/families, manifesting primarily non-syndromic CC (15/18, 83%), remained unsolved. The identified variants (8 P, 12 LP, 10 VUS-PP, and 5 VUS), half of which were unreported in the literature, affected five functional categories of genes involved in transcription/splicing, lens formation/homeostasis (i.e., crystallin genes), membrane signalling, cell-cell interaction, and immune response. A phenotype-specific variant clustering was observed in four genes (KIF1A, MAF, PAX6, SPTAN1), whereas variable expressivity and potential phenotypic expansion in two (BCOR, NHS) and five genes (CWC27, KIF1A, IFIH1, PAX6, SPTAN1), respectively. Finally, ES allowed to detect variants in six genes not commonly included in commercial CC panels. These findings broaden the genotype-phenotype correlations in one of the largest CC cohorts tested by ES, providing novel insights into the underlying pathogenetic mechanisms and emphasising the power of ES as first-tier test.

Genetic architecture of oral glucose-stimulated insulin release provides biological insights into type 2 diabetes aetiology

The genetics of β-cell function (BCF) offer valuable insights into the aetiology of type 2 diabetes (T2D)1,2. Previous studies have expanded the catalogue of BCF genetic associations through candidate gene studies3-7, large-scale genome-wide association studies (GWAS) of fasting BCF8,9 or functional islet studies on T2D risk variants10-14. Nonetheless, GWAS focused on BCF traits derived from oral glucose tolerance test (OGTT) data have been limited in sample size15,16 and have often overlooked the potential for related traits to capture distinct genetic features of insulin-producing β-cells17,18. We reasoned that investigating the genetic basis of multiple BCF estimates could provide a broader understanding of β-cell physiology. Here, we aggregate GWAS data of eight OGTT-based BCF traits from ~26,000 individuals of European descent, identifying 55 independent genetic associations at 44 loci. By examining the effects of BCF genetic signals on related phenotypes, we uncover diverse disease mechanisms whereby genetic regulation of BCF may influence T2D risk. Integrating BCF-GWAS data with pancreatic islet transcriptomic and epigenomic datasets reveals 92 candidate effector genes. Gene silencing in β-cell models highlights ACSL1 and FAM46C as key regulators of insulin secretion. Overall, our findings yield insights into the biology of insulin release and the molecular processes linking BCF to T2D risk, shedding light on the heterogeneity of T2D pathophysiology.

A prioritization tool for cilia-associated genes and their in vivo resources unveils new avenues for ciliopathy research

Defects in ciliary signaling or mutations in proteins that localize to primary cilia lead to a class of human diseases known as ciliopathies. Approximately 10% of mammalian genes encode cilia-associated proteins, and a major gap in the cilia research field is knowing which genes to prioritize to study and finding the in vivo vertebrate mutant alleles and reagents available for their study. Here, we present a unified resource listing the cilia-associated human genes cross referenced to available mouse and zebrafish mutant alleles, and their associated phenotypes, as well as expression data in the kidney and functional data for vertebrate Hedgehog signaling. This resource empowers researchers to easily sort and filter genes based on their own expertise and priorities, cross reference with newly generated -omics datasets, and quickly find in vivo resources and phenotypes associated with a gene of interest.

Phoenixin knockout mice show no impairment in fertility or differences in metabolic response to a high-fat diet, but exhibit behavioral differences in an open field test

Phoenixin (PNX) is a conserved secreted peptide that was identified 10 years ago with numerous studies published on its pleiotropic functions. PNX is associated with estrous cycle length, protection from a high-fat diet, and reduction of anxiety behavior. However, no study had yet evaluated the impact of deleting PNX in the whole animal. We sought to evaluate a mouse model lacking the PNX parent gene, small integral membrane protein 20 (Smim20), and the resulting effect on reproduction, energy homeostasis, and anxiety. We found that the Smim20 knockout mice had normal fertility and estrous cycle lengths. Consistent with normal fertility, the hypothalamii of the knockout mice showed no changes in the levels of reproduction-related genes, but the male mice had some changes in energy homeostasis-related genes, such as melanocortin receptor 4 (Mc4r). When placed on a high-fat diet, the wildtype and knockout mice responded similarly, but the male heterozygous mice gained slightly less weight. When placed in an open field test box, the female knockout mice traveled less distance in the outer zone, indicating alterations in anxiety or locomotor behavior. In summary, the homozygous knockout of PNX did not alter fertility and modestly alters a few neuroendocrine genes in response to a high-fat diet, especially in the female mice. However, it altered the behavior of mice in an open field test. PNX therefore may not be crucial for reproductive function or weight, however, we cannot rule out possible compensatory mechanisms in the knockout model. Understanding the role of PNX in physiology may ultimately lead to an enhanced understanding of neuroendocrine mechanisms involving this enigmatic peptide.

Working with Miraculous Mice: Mus musculus as a Model Organism

The laboratory mouse has been described as a "miracle" model organism, providing a window by which we may gain an understanding of ourselves. Since the first recorded mouse experiment in 1664, the mouse has become the most used animal model in biomedical research. Mice are ideally suited as a model organism because of their small size, short gestation period, large litter size, and genetic similarity to humans. This article provides a broad overview of the laboratory mouse as a model organism and is intended for undergraduates and those new to working with mice. We delve into the history of the laboratory mouse and outline important terminology to accurately describe research mice. The types of laboratory mice available to researchers are reviewed, including outbred stocks, inbred strains, immunocompromised mice, and genetically engineered mice. The critical role mice have played in advancing knowledge in the areas of oncology, immunology, and pharmacology is highlighted by examining the significant contribution of mice to Nobel Prize winning research. International mouse mutagenesis programs and accurate phenotyping of mouse models are outlined. We also explain important considerations for working with mice, including animal ethics; the welfare principles of replacement, refinement, and reduction; and the choice of mouse model in experimental design. Finally, we present practical advice for maintaining a mouse colony, which involves adequate training of staff, the logistics of mouse housing, monitoring colony health, and breeding strategies. Useful resources for working with mice are also listed. The aim of this overview is to equip the reader with a broad appreciation of the enormous potential and some of the complexities of working with the laboratory mouse in a quest to improve human health. © 2024 The Author(s). Current Protocols published by Wiley Periodicals LLC.

The Unified Phenotype Ontology (uPheno): A framework for cross-species integrative phenomics

Phenotypic data are critical for understanding biological mechanisms and consequences of genomic variation, and are pivotal for clinical use cases such as disease diagnostics and treatment development. For over a century, vast quantities of phenotype data have been collected in many different contexts covering a variety of organisms. The emerging field of phenomics focuses on integrating and interpreting these data to inform biological hypotheses. A major impediment in phenomics is the wide range of distinct and disconnected approaches to recording the observable characteristics of an organism. Phenotype data are collected and curated using free text, single terms or combinations of terms, using multiple vocabularies, terminologies, or ontologies. Integrating these heterogeneous and often siloed data enables the application of biological knowledge both within and across species. Existing integration efforts are typically limited to mappings between pairs of terminologies; a generic knowledge representation that captures the full range of cross-species phenomics data is much needed. We have developed the Unified Phenotype Ontology (uPheno) framework, a community effort to provide an integration layer over domain-specific phenotype ontologies, as a single, unified, logical representation. uPheno comprises (1) a system for consistent computational definition of phenotype terms using ontology design patterns, maintained as a community library; (2) a hierarchical vocabulary of species-neutral phenotype terms under which their species-specific counterparts are grouped; and (3) mapping tables between species-specific ontologies. This harmonized representation supports use cases such as cross-species integration of genotype-phenotype associations from different organisms and cross-species informed variant prioritization.

Clinical-transcriptional prioritization of the circulating proteome in human heart failure

Given expanding studies in epidemiology and disease-oriented human studies offering hundreds of associations between the human "ome" and disease, prioritizing molecules relevant to disease mechanisms among this growing breadth is important. Here, we link the circulating proteome to human heart failure (HF) propensity (via echocardiographic phenotyping and clinical outcomes) across the lifespan, demonstrating key pathways of fibrosis, inflammation, metabolism, and hypertrophy. We observe a broad array of genes encoding proteins linked to HF phenotypes and outcomes in clinical populations dynamically expressed at a transcriptional level in human myocardium during HF and cardiac recovery (several in a cell-specific fashion). Many identified targets do not have wide precedent in large-scale genomic discovery or human studies, highlighting the complementary roles for proteomic and tissue transcriptomic discovery to focus epidemiological targets to those relevant in human myocardium for further interrogation.

Overabundant endocannabinoids in neurons are detrimental to cognitive function

2-Arachidonoylglycerol (2-AG) is the most prevalent endocannabinoid involved in maintaining brain homeostasis. Previous studies have demonstrated that inactivating monoacylglycerol lipase (MAGL), the primary enzyme responsible for degrading 2-AG in the brain, alleviates neuropathology and prevents synaptic and cognitive decline in animal models of neurodegenerative diseases. However, we show that selectively inhibiting 2-AG metabolism in neurons impairs cognitive function in mice. This cognitive impairment appears to result from decreased expression of synaptic proteins and synapse numbers, impaired long-term synaptic plasticity and cortical circuit functional connectivity, and diminished neurogenesis. Interestingly, the synaptic and cognitive deficits induced by neuronal MAGL inactivation can be counterbalanced by inhibiting astrocytic 2-AG metabolism. Transcriptomic analyses reveal that inhibiting neuronal 2-AG degradation leads to widespread changes in expression of genes associated with synaptic function. These findings suggest that crosstalk in 2-AG signaling between astrocytes and neurons is crucial for maintaining synaptic and cognitive functions and that excessive 2-AG in neurons alone is detrimental to cognitive function.

USP50 suppresses alternative RecQ helicase use and deleterious DNA2 activity during replication

Mammalian DNA replication relies on various DNA helicase and nuclease activities to ensure accurate genetic duplication, but how different helicase and nuclease activities are properly directed remains unclear. Here, we identify the ubiquitin-specific protease, USP50, as a chromatin-associated protein required to promote ongoing replication, fork restart, telomere maintenance, cellular survival following hydroxyurea or pyridostatin treatment, and suppression of DNA breaks near GC-rich sequences. We find that USP50 supports proper WRN-FEN1 localisation at or near stalled replication forks. Nascent DNA in cells lacking USP50 shows increased association of the DNA2 nuclease and RECQL4 and RECQL5 helicases and replication defects in cells lacking USP50, or FEN1 are driven by these proteins. Consequently, suppression of DNA2 or RECQL4/5 improves USP50-depleted cell resistance to agents inducing replicative stress and restores telomere stability. These data define an unexpected regulatory protein that promotes the balance of helicase and nuclease use at ongoing and stalled replication forks.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) mediate alternative splicing, but mechanistic interpretation is hindered by the technical limitations of short-read RNA sequencing (RNA-seq), which cannot directly link splicing events to full-length protein isoforms. Long-read RNA-seq represents a powerful tool to characterize transcript isoforms, and recently, infer protein isoform existence. Here, we present an approach that integrates information from GWASs, splicing quantitative trait loci (sQTLs), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes that colocalized with BMD associations (H4PP ≥ 0.75). We generated PacBio Iso-Seq data (N = ∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were unannotated. By casting the sQTLs onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense-mediated decay and 190 that potentially resulted in the expression of unannotated protein isoforms. Finally, we functionally validated colocalizing sQTLs in TPM2, in which siRNA-mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization but exhibited no effect upon knockdown of the entire gene. Our approach should be to generalize across diverse clinical traits and to provide insights into protein isoform activities modulated by GWAS loci.

A Large-Scale Genome-Wide Gene-Sleep Interaction Study in 732,564 Participants Identifies Lipid Loci Explaining Sleep-Associated Lipid Disturbances

We performed large-scale genome-wide gene-sleep interaction analyses of lipid levels to identify novel genetic variants underpinning the biomolecular pathways of sleep-associated lipid disturbances and to suggest possible druggable targets. We collected data from 55 cohorts with a combined sample size of 732,564 participants (87% European ancestry) with data on lipid traits (high-density lipoprotein [HDL-c] and low-density lipoprotein [LDL-c] cholesterol and triglycerides [TG]). Short (STST) and long (LTST) total sleep time were defined by the extreme 20% of the age- and sex-standardized values within each cohort. Based on cohort-level summary statistics data, we performed meta-analyses for the one-degree of freedom tests of interaction and two-degree of freedom joint tests of the main and interaction effect. In the cross-population meta-analyses, the one-degree of freedom variant-sleep interaction test identified 10 loci (P int <5.0e-9) not previously observed for lipids. Of interest, the ASPH locus (TG, LTST) is a target for aspartic and succinic acid metabolism previously shown to improve sleep and cardiovascular risk. The two-degree of freedom analyses identified an additional 7 loci that showed evidence for variant-sleep interaction (P joint <5.0e-9 in combination with P int <6.6e-6). Of these, the SLC8A1 locus (TG, STST) has been considered a potential treatment target for reduction of ischemic damage after acute myocardial infarction. Collectively, the 17 (9 with STST; 8 with LTST) loci identified in this large-scale initiative provides evidence into the biomolecular mechanisms underpinning sleep-duration-associated changes in lipid levels. The identified druggable targets may contribute to the development of novel therapies for dyslipidemia in people with sleep disturbances.

SLC25A48 influences plasma levels of choline and localizes to the inner mitochondrial membrane

Choline contributes to the biogenesis of methyl groups, neurotransmitters, and cell membranes. Our genome-wide association study (GWAS) of circulating choline in 2228 college students found that alleles in SLC25A48 (rs6596270) influence choline concentrations in men (p = 9.6 × 10-8), but not women. Previously, the subcellular location and function of SLC25A48 were unknown. Using super-resolution immunofluorescence microscopy, we localized SLC25A48 to the inner mitochondrial membrane. Our results suggest that SLC25A48 transports choline across the inner mitochondrial membrane.

Ancestral origins and post-admixture adaptive evolution of highland Tajiks

It remains debatable how many genes and how various the mechanisms are behind human adaptation to extreme environments, such as high altitudes. Despite extensive studies on Tibetans, Andeans and Ethiopians, new insights are expected to be provided with careful analysis of underrepresented highlanders living in a different geographical region, such as the Tajiks, who reside on the Pamir Plateau at an average altitude exceeding 4000 meters. Moreover, genetic admixture, as we observed in the current whole-genome deep-sequencing study of Xinjiang Tajiks (XJT), offers a unique opportunity to explore how admixture may facilitate adaptation to high-altitude environments. Compared with other extensively studied highlanders, XJT showed pronounced admixture patterns: most of their ancestry are derived from West Eurasians (34.5%-48.3%) and South Asians (21.4%-40.0%), and some minor ancestry from East Asians and Siberians (3.62%-17.5%). The greater genetic diversity in XJT than in their ancestral source populations provides a genetic basis for their adaptation to high-altitude environments. The admixture gain of functional adaptive components from ancestral populations could facilitate adaptation to high-altitude environments. Specifically, admixture-facilitated adaptation was strongly associated with skin-related candidate genes that respond to UV radiation (e.g. HERC2 and BNC2) and cardiovascular-system-related genes (e.g. MPI and BEST1). Notably, no adaptive variants of genes showing outstanding natural selection signatures in the Tibetan or Andean highlanders were identified in XJT, including EPAS1 and EGLN1, indicating that a different set of genes contributed to XJT's survival on the Pamir Plateau, although some genes underlying natural selection in XJT have been previously reported in other highlanders. Our results highlight the unique genetic adaptations in XJT and propose that admixture may play a vital role in facilitating high-altitude adaptation. By introducing and elevating diversity, admixture likely induces novel genetic factors that contribute to the survival of populations in extreme environments like the highlands.

SMYD5 is a ribosomal methyltransferase that catalyzes RPL40 lysine methylation to enhance translation output and promote hepatocellular carcinoma

While lysine methylation is well-known for regulating gene expression transcriptionally, its implications in translation have been largely uncharted. Trimethylation at lysine 22 (K22me3) on RPL40, a core ribosomal protein located in the GTPase activation center, was first reported 27 years ago. Yet, its methyltransferase and role in translation remain unexplored. Here, we report that SMYD5 has robust in vitro activity toward RPL40 K22 and primarily catalyzes RPL40 K22me3 in cells. The loss of SMYD5 and RPL40 K22me3 leads to reduced translation output and disturbed elongation as evidenced by increased ribosome collisions. SMYD5 and RPL40 K22me3 are upregulated in hepatocellular carcinoma (HCC) and negatively correlated with patient prognosis. Depleting SMYD5 renders HCC cells hypersensitive to mTOR inhibition in both 2D and 3D cultures. Additionally, the loss of SMYD5 markedly inhibits HCC development and growth in both genetically engineered mouse and patient-derived xenograft (PDX) models, with the inhibitory effect in the PDX model further enhanced by concurrent mTOR suppression. Our findings reveal a novel role of the SMYD5 and RPL40 K22me3 axis in translation elongation and highlight the therapeutic potential of targeting SMYD5 in HCC, particularly with concurrent mTOR inhibition. This work also conceptually broadens the understanding of lysine methylation, extending its significance from transcriptional regulation to translational control.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

The phenotypic spectrum of PTCD3 deficiency

The PTCD3 gene product (protein PTCD3 or MRPS39) forms the entry channel of the mitochondrial small ribosomal subunit and binds to single-stranded mRNA. Here, we expand on the clinical manifestations of PTCD3 pathogenic variants by describing an early-onset patient with Leigh-like syndrome and two patients with milder form of disease, with combined oxidative phosphorylation deficiency. A 34-year-old male and his 33-year-old sister both have horizontal nystagmus, pronounced rough tremor, truncal ataxia, dysmetria, spasticity and hyperreflexia. The basal respiration rate decreased significantly for the male patient and his mother (p < 0.0001) compared to the controls. The whole genome sequencing analysis revealed two heterozygous variants in the PTCD3: c.1182T>A, p.(Tyr394Ter) and c.805C>T, p.(His269Tyr). Tyr394Ter variant ablates the C-terminal half of the protein, including a significant portion of the central fold. In silico modelling for the variant His269Tyr shows that the inclusion of the slightly larger tyrosine sidechain is well tolerated, with no significant change in either the position or the movement of the surrounding area. The third case is a 9-year-old boy, who has a global developmental delay, central hypotonia, hyperreflexia and abnormal MRI. PTCD3 pathogenic variant c.538+4A>G was identified by whole exome sequencing. To test the variant's effect on splicing, an RT-PCR experiment was performed, which revealed skipping of an out-of-frame exon 7.

Deep learning-derived splenic radiomics, genomics, and coronary artery disease

Background

Despite advances in managing traditional risk factors, coronary artery disease (CAD) remains the leading cause of mortality. Circulating hematopoietic cells influence risk for CAD, but the role of a key regulating organ, spleen, is unknown. The understudied spleen is a 3-dimensional structure of the hematopoietic system optimally suited for unbiased radiologic investigations toward novel mechanistic insights.

Methods

Deep learning-based image segmentation and radiomics techniques were utilized to extract splenic radiomic features from abdominal MRIs of 42,059 UK Biobank participants. Regression analysis was used to identify splenic radiomics features associated with CAD. Genome-wide association analyses were applied to identify loci associated with these radiomics features. Overlap between loci associated with CAD and the splenic radiomics features was explored to understand the underlying genetic mechanisms of the role of the spleen in CAD.

Results

We extracted 107 splenic radiomics features from abdominal MRIs, and of these, 10 features were associated with CAD. Genome-wide association analysis of CAD-associated features identified 219 loci, including 35 previously reported CAD loci, 7 of which were not associated with conventional CAD risk factors. Notably, variants at 9p21 were associated with splenic features such as run length non-uniformity.

Conclusions

Our study, combining deep learning with genomics, presents a new framework to uncover the splenic axis of CAD. Notably, our study provides evidence for the underlying genetic connection between the spleen as a candidate causal tissue-type and CAD with insight into the mechanisms of 9p21, whose mechanism is still elusive despite its initial discovery in 2007. More broadly, our study provides a unique application of deep learning radiomics to non-invasively find associations between imaging, genetics, and clinical outcomes.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Genome-wide association study of delay discounting in Heterogeneous Stock rats

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders and multiple co-occurring psychopathologies. Human and animal genetic studies have established that delay discounting is heritable, but only a few associated genes have been identified. We aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogeneous Stock (HS) rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female HS rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at various delays. Preference switch points were calculated and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and one indifference point identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family, was the sole gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression might be responsible for the association with behavior. Adgrl3, which encodes a latrophilin subfamily G-protein coupled receptor, was the sole gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

The 'dispanins' and related proteins in physiology and neurological disease

The dispanins are a family of 15 transmembrane proteins that have diverse and often unclear physiological functions. Many dispanins, including synapse differentiation induced gene 1 (SynDIG1), proline-rich transmembrane protein 1 (PRRT1)/SynDIG4, and PRRT2, are expressed in the central nervous system (CNS), where they are involved in the development of synapses, regulation of neurotransmitter release, and interactions with ion channels, including AMPA receptors (AMPARs). Others, including transmembrane protein 233 (TMEM233) and trafficking regulator of GLUT4-1 (TRARG1), are expressed in the peripheral nervous system (PNS); however, the function of these dispanins is less clear. Recently, a family of neurotoxins isolated from the giant Australian stinging tree was shown to target TMEM233 to modulate the function of voltage-gated sodium (NaV) channels, suggesting that the dispanins are inherently druggable. Here, we review current knowledge about the structure and function of the dispanins, in particular TMEM233 and its two most closely related homologs PRRT2 and TRARG1, which may be drug targets involved in neurological disease.

Bayesian estimation of gene constraint from an evolutionary model with gene features

Measures of selective constraint on genes have been used for many applications, including clinical interpretation of rare coding variants, disease gene discovery and studies of genome evolution. However, widely used metrics are severely underpowered at detecting constraints for the shortest ~25% of genes, potentially causing important pathogenic mutations to be overlooked. Here we developed a framework combining a population genetics model with machine learning on gene features to enable accurate inference of an interpretable constraint metric, shet. Our estimates outperform existing metrics for prioritizing genes important for cell essentiality, human disease and other phenotypes, especially for short genes. Our estimates of selective constraint should have wide utility for characterizing genes relevant to human disease. Finally, our inference framework, GeneBayes, provides a flexible platform that can improve the estimation of many gene-level properties, such as rare variant burden or gene expression differences.

Disruption of the ZFP574-THAP12 complex suppresses B cell malignancies in mice

Despite the availability of life-extending treatments for B cell leukemias and lymphomas, many of these cancers remain incurable. Thus, the development of new molecular targets and therapeutics is needed to expand treatment options. To identify new molecular targets, we used a forward genetic screen in mice to identify genes required for development or survival of lymphocytes. Here, we describe Zfp574, an essential gene encoding a zinc finger protein necessary for normal and malignant lymphocyte survival. We show that ZFP574 interacts with zinc finger protein THAP12 and promotes the G1-to-S-phase transition during cell cycle progression. Mutation of ZFP574 impairs nuclear localization of the ZFP574-THAP12 complex. ZFP574 or THAP12 deficiency results in cell cycle arrest and impaired lymphoproliferation. Germline mutation, acute gene deletion, or targeted degradation of ZFP574 suppressed Myc-driven B cell leukemia in mice, but normal B cells were largely spared, permitting long-term survival, whereas complete lethality was observed in control animals. Our findings support the identification of drugs targeting ZFP574-THAP12 as a unique strategy to treat B cell malignancies.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

PURPOSE

This study aimed to identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs).

METHODS

We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single-nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations.

RESULTS

Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43 of 467 probands (9.2%), and prioritized variants of uncertain significance in 70 of 467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (eg, MYH10 [HGNC:7568], KIF21B [HGNC:29442], TGFBR2 [HGNC:11773], and TUBB6 [HGNC:20776]), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (eg, CDK13 [HGNC:1733], TGFB2 [HGNC:11768]), genes with no reported association with oCCDDs or the syndromic phenotypes (eg, TUBA4A [HGNC:12407], KIF5C [HGNC:6325], CTNNA1 [HGNC:2509], KLB [HGNC:15527], FGF21 [HGNC:3678]), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses.

CONCLUSION

This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to 23 genes. Investigating these genes' functional implications shed light on neurological, thyroidal, bone metabolism, and hematopoietic pathways that necessitate future investigation for blood pressure management that caters to sleep health lifestyle. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausible nature of distinct influences of both sleep duration extremes in cardiovascular health. Several of our loci are specific towards a particular population background or sex, emphasizing the importance of addressing heterogeneity entangled in gene-environment interactions, when considering precision medicine design approaches for blood pressure management.

Perturbation of the insomnia WDR90 genome-wide association studies locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome-wide association studies (GWAS) have identified loci for sleep-related traits, they do not directly uncover the underlying causal variants and corresponding effector genes. The majority of such variants reside in non-coding regions and are therefore presumed to impact cis-regulatory elements. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated phosphatidyl inositol glycan (PIG)-Q as a functionally relevant gene at the insomnia "WDR90" GWAS locus. However, importantly that effort did not characterize the corresponding underlying causal variant. Specifically, our previous 3D genomic datasets nominated a shortlist of three neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium within an intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. We sought to investigate the influence of these SNPs collectively and then individually on PIG-Q modulation to pinpoint the causal "regulatory" variant. Starting with gross level perturbation, deletion of the entire region in NPCs via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from individual luciferase reporter assays for each SNP in iPSCs revealed that the region with the rs3752495 risk allele (RA) induced a ~2.5-fold increase in luciferase expression. Importantly, rs3752495 also exhibited an allele-specific effect, with the RA increasing the luciferase expression by ~2-fold versus the non-RA. In conclusion, our variant-to-function approach and in vitro validation implicate rs3752495 as a causal insomnia variant embedded within WDR90 while modulating the expression of the distally located PIG-Q.

Establishment and characterization of an hACE2/hTMPRSS2 knock-in mouse model to study SARS-CoV-2

Despite a substantial body of research, we lack fundamental understanding of the pathophysiology of COVID-19 caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) including pulmonary and cardiovascular outcomes, in part due to limitations of murine models. Most models use transgenic mice (K18) that express the human (h) angiotensin converting enzyme 2 (ACE2), ACE2 knock-in (KI) mice, or mouse-adapted strains of SARS-CoV-2. Further, many SARS-CoV-2 variants produce fatal neurologic disease in K18 mice and most murine studies focus only on acute disease in the first 14 days post inoculation (dpi). To better enable understanding of both acute (<14 dpi) and post-acute (>14 dpi) infection phases, we describe the development and characterization of a novel non-lethal KI mouse that expresses both the ACE2 and transmembrane serine protease 2 (TMPRSS2) genes (hACE2/hTMPRSS2). The human genes were engineered to replace the orthologous mouse gene loci but remain under control of their respective murine promoters, resulting in expression of ACE2 and TMPRSS2 instead of their murine counterparts. After intranasal inoculation with an omicron strain of SARS-CoV-2, hACE2/hTMPRSS2 KI mice transiently lost weight but recovered by 7 dpi. Infectious SARS-CoV-2 was detected in nasopharyngeal swabs 1-2 dpi and in lung tissues 2-6 dpi, peaking 4 dpi. These outcomes were similar to those in K18 mice that were inoculated in parallel. To determine the extent to which hACE2/hTMPRSS2 KI mice are suitable to model pulmonary and cardiovascular outcomes, physiological assessments measuring locomotion, behavior and reflexes, biomonitoring to measure cardiac activity and respiration, and micro computed tomography to assess lung function were conducted frequently to 6 months post inoculation. Male but not female SARS-CoV-2 inoculated hACE2/hTMPRSS2 KI mice showed a transient reduction in locomotion compared to control saline treated mice. No significant changes in respiration, oxygen saturation, heart rate variability, or conductivity were detected in SARS-CoV-2 inoculated mice of either sex. When re-inoculated 6 months after the first inoculation, hACE2/hTMPRSS2 KI became re-infected with disease signs similar to after the first inoculation. Together these data show that a newly generated hACE2/hTMPRSS2 KI mouse can be used to study mild COVID-19.

Characterising the genetic architecture of changes in adiposity during adulthood using electronic health records

Obesity is a heritable disease, characterised by excess adiposity that is measured by body mass index (BMI). While over 1,000 genetic loci are associated with BMI, less is known about the genetic contribution to adiposity trajectories over adulthood. We derive adiposity-change phenotypes from 24.5 million primary-care health records in over 740,000 individuals in the UK Biobank, Million Veteran Program USA, and Estonian Biobank, to discover and validate the genetic architecture of adiposity trajectories. Using multiple BMI measurements over time increases power to identify genetic factors affecting baseline BMI by 14%. In the largest reported genome-wide study of adiposity-change in adulthood, we identify novel associations with BMI-change at six independent loci, including rs429358 (APOE missense variant). The SNP-based heritability of BMI-change (1.98%) is 9-fold lower than that of BMI. The modest genetic correlation between BMI-change and BMI (45.2%) indicates that genetic studies of longitudinal trajectories could uncover novel biology of quantitative traits in adulthood.

Morphological simulation tests the limits on phenotype discovery in 3D image analysis

In the past few decades, advances in 3D imaging have created new opportunities for reverse genetic screens. Rapidly growing datasets of 3D images of genetic knockouts require high-throughput, automated computational approaches for identifying and characterizing new phenotypes. However, exploratory, discovery-oriented image analysis pipelines used to discover these phenotypes can be difficult to validate because, by their nature, the expected outcome is not known a priori . Introducing known morphological variation through simulation can help distinguish between real phenotypic differences and random variation; elucidate the effects of sample size; and test the sensitivity and reproducibility of morphometric analyses. Here we present a novel approach for 3D morphological simulation that uses open-source, open-access tools available in 3D Slicer, SlicerMorph, and Advanced Normalization Tools in R (ANTsR). While we focus on diffusible-iodine contrast-enhanced micro-CT (diceCT) images, this approach can be used on any volumetric image. We then use our simulated datasets to test whether tensor-based morphometry (TBM) can recover our introduced differences; to test how effect size and sample size affect detectability; and to determine the reproducibility of our results. In our approach to morphological simulation, we first generate a simulated deformation based on a reference image and then propagate this deformation to subjects using inverse transforms obtained from the registration of subjects to the reference. This produces a new dataset with a shifted population mean while retaining individual variability because each sample deforms more or less based on how different or similar it is from the reference. TBM is a widely-used technique that statistically compares local volume differences associated with local deformations. Our results showed that TBM recovered our introduced morphological differences, but that detectability was dependent on the effect size, the sample size, and the region of interest (ROI) included in the analysis. Detectability of subtle phenotypes can be improved both by increasing the sample size and by limiting analyses to specific body regions. However, it is not always feasible to increase sample sizes in screens of essential genes. Therefore, methodical use of ROIs is a promising way to increase the power of TBM to detect subtle phenotypes. Generating known morphological variation through simulation has broad applicability in developmental, evolutionary, and biomedical morphometrics and is a useful way to distinguish between a failure to detect morphological difference and a true lack of morphological difference. Morphological simulation can also be applied to AI-based supervised learning to augment datasets and overcome dataset limitations.

Lethal phenotypes in Mendelian disorders

PURPOSE

Existing resources that characterize the essentiality status of genes are based on either proliferation assessment in human cell lines, viability evaluation in mouse knockouts, or constraint metrics derived from human population sequencing studies. Several repositories document phenotypic annotations for rare disorders; however, there is a lack of comprehensive reporting on lethal phenotypes.

METHODS

We queried Online Mendelian Inheritance in Man for terms related to lethality and classified all Mendelian genes according to the earliest age of death recorded for the associated disorders, from prenatal death to no reports of premature death. We characterized the genes across these lethality categories, examined the evidence on viability from mouse models and explored how this information could be used for novel gene discovery.

RESULTS

We developed the Lethal Phenotypes Portal to showcase this curated catalog of human essential genes. Differences in the mode of inheritance, physiological systems affected, and disease class were found for genes in different lethality categories, as well as discrepancies between the lethal phenotypes observed in mouse and human.

CONCLUSION

We anticipate that this resource will aid clinicians in the diagnosis of early lethal conditions and assist researchers in investigating the properties that make these genes essential for human development.

Systems genetics analysis of human body fat distribution genes identifies adipocyte processes

Excess abdominal fat is a sexually dimorphic risk factor for cardio-metabolic disease and is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Whereas this trait is highly heritable, few causal genes are known. We aimed to identify novel drivers of WHRadjBMI using systems genetics. We used two independent cohorts of adipose tissue gene expression and constructed sex- and depot-specific Bayesian networks to model gene-gene interactions from 8,492 genes. Using key driver analysis, we identified genes that, in silico and putatively in vitro, regulate many others. 51-119 key drivers in each network were replicated in both cohorts. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We overexpressed or down-regulated seven key driver genes in human subcutaneous pre-adipocytes. Key driver genes ANAPC2 and RSPO1 inhibited adipogenesis, whereas PSME3 increased adipogenesis. RSPO1 increased Wnt signaling activity. In differentiated adipocytes, MIGA1 and UBR1 down-regulation led to mitochondrial dysfunction. These five genes regulate adipocyte function, and we hypothesize that they regulate fat distribution.

Physiological basis for xenotransplantation from genetically modified pigs to humans

The collective efforts of scientists over multiple decades have led to advancements in molecular and cellular biology-based technologies including genetic engineering and animal cloning that are now being harnessed to enhance the suitability of pig organs for xenotransplantation into humans. Using organs sourced from pigs with multiple gene deletions and human transgene insertions, investigators have overcome formidable immunological and physiological barriers in pig-to-nonhuman primate (NHP) xenotransplantation and achieved prolonged pig xenograft survival. These studies informed the design of Revivicor's (Revivicor Inc, Blacksburg, VA) genetically engineered pigs with 10 genetic modifications (10 GE) (including the inactivation of 4 endogenous porcine genes and insertion of 6 human transgenes), whose hearts and kidneys have now been studied in preclinical human xenotransplantation models with brain-dead recipients. Additionally, the first two clinical cases of pig-to-human heart xenotransplantation were recently performed with hearts from this 10 GE pig at the University of Maryland. Although this review focuses on xenotransplantation of hearts and kidneys, multiple organs, tissues, and cell types from genetically engineered pigs will provide much-needed therapeutic interventions in the future.

Differential Analysis of Cereblon Neosubstrates in Rabbit Embryos Using Targeted Proteomics

Targeted protein degradation is the selective removal of a protein of interest through hijacking intracellular protein cleanup machinery. This rapidly growing field currently relies heavily on the use of the E3 ligase cereblon (CRBN) to target proteins for degradation, including the immunomodulatory drugs (IMiDs) thalidomide, lenalidomide, and pomalidomide which work through a molecular glue mechanism of action with CRBN. While CRBN recruitment can result in degradation of a specific protein of interest (e.g., efficacy), degradation of other proteins (called CRBN neosubstrates) also occurs. Degradation of one or more of these CRBN neosubstrates is believed to play an important role in thalidomide-related developmental toxicity observed in rabbits and primates. We identified a set of 25 proteins of interest associated with CRBN-related protein homeostasis and/or embryo/fetal development. We developed a targeted assay for these proteins combining peptide immunoaffinity enrichment and high-resolution mass spectrometry and successfully applied this assay to rabbit embryo samples from pregnant rabbits dosed with three IMiDs. We confirmed previously reported in vivo decreases in neosubstrates like SALL4, as well as provided evidence of neosubstrate changes for proteins only examined in vitro previously. While there were many proteins that were similarly decreased by all three IMiDs, no compound had the exact same neosubstrate degradation profile as another. We compared our data to previous literature reports of IMiD-induced degradation and known developmental biology associations. Based on our observations, we recommend monitoring at least a major subset of these neosubstrates in a developmental test system to improve CRBN-binding compound-specific risk assessment. A strength of our assay is that it is configurable, and the target list can be readily adapted to focus on only a subset of proteins of interest or expanded to incorporate new findings as additional information about CRBN biology is discovered.

Systems genetics of metabolic health in the BXD mouse genetic reference population

Susceptibility to metabolic syndrome (MetS) is dependent on genetics, environment, and gene-by-environment interactions, rendering the study of underlying mechanisms challenging. The majority of experiments in model organisms do not incorporate genetic variation and lack specific evaluation criteria for MetS. Here, we derived a continuous metric, the metabolic health score (MHS), based on standard clinical parameters and defined its molecular signatures in the liver and circulation. In human UK Biobank, the MHS associated with MetS status and was predictive of future disease incidence, even in individuals without MetS. Using quantitative trait locus analyses in mice, we found two MHS-associated genetic loci and replicated them in unrelated mouse populations. Through a prioritization scheme in mice and human genetic data, we identified TNKS and MCPH1 as candidates mediating differences in the MHS. Our findings provide insights into the molecular mechanisms sustaining metabolic health across species and uncover likely regulators. A record of this paper's transparent peer review process is included in the supplemental information.

Viable mutations of mouse midnolin suppress B cell malignancies

In a genetic screen, we identified two viable missense alleles of the essential gene Midnolin (Midn) that were associated with reductions in peripheral B cells. Causation was confirmed in mice with targeted deletion of four of six MIDN protein isoforms. MIDN was expressed predominantly in lymphocytes where it augmented proteasome activity. We showed that purified MIDN directly stimulated 26S proteasome activity in vitro in a manner dependent on the ubiquitin-like domain and a C-terminal region. MIDN-deficient B cells displayed aberrant activation of the IRE-1/XBP-1 pathway of the unfolded protein response. Partial or complete MIDN deficiency strongly suppressed Eμ-Myc-driven B cell leukemia and the antiapoptotic effects of Eμ-BCL2 on B cells in vivo and induced death of Sp2/0 hybridoma cells in vitro, but only partially impaired normal lymphocyte development. Thus, MIDN is required for proteasome activity in support of normal lymphopoiesis and is essential for malignant B cell proliferation over a broad range of differentiation states.

Computational identification of disease models through cross-species phenotype comparison

The use of standardised phenotyping screens to identify abnormal phenotypes in mouse knockouts, together with the use of ontologies to describe such phenotypic features, allows the implementation of an automated and unbiased pipeline to identify new models of disease by performing phenotype comparisons across species. Using data from the International Mouse Phenotyping Consortium (IMPC), approximately half of mouse mutants are able to mimic, at least partially, the human ortholog disease phenotypes as computed by the PhenoDigm algorithm. We found the number of phenotypic abnormalities in the mouse and the corresponding Mendelian disorder, the pleiotropy and severity of the disease, and the viability and zygosity status of the mouse knockout to be associated with the ability of mouse models to recapitulate the human disorder. An analysis of the IMPC impact on disease gene discovery through a publication-tracking system revealed that the resource has been implicated in at least 109 validated rare disease-gene associations over the last decade.

GPCR Screening Reveals that the Metabolite Receptor HCAR3 Regulates Epithelial Proliferation, Migration, and Cellular Respiration

Epithelial cells in the skin and other tissues rely on signals from their environment to maintain homeostasis and respond to injury, and GPCRs play a critical role in this communication. A better understanding of the GPCRs expressed in epithelial cells will contribute to understanding the relationship between cells and their niche and could lead to developing new therapies to modulate cell fate. This study used human primary keratinocytes as a model to investigate the specific GPCRs regulating epithelial cell proliferation and differentiation. We identified 3 key receptors-HCAR3, LTB4R, and GPR137-and found that knockdown of these receptors led to changes in numerous gene networks that are important for maintaining cell identity and promoting proliferation while inhibiting differentiation. Our study also revealed that the metabolite receptor HCAR3 regulates keratinocyte migration and cellular metabolism. Knockdown of HCAR3 led to reduced keratinocyte migration and respiration, which could be attributed to altered metabolite use and aberrant mitochondrial morphology caused by the absence of the receptor. This study contributes to understanding the complex interplay between GPCR signaling and epithelial cell fate decisions.

Epigenetic inheritance of diet-induced and sperm-borne mitochondrial RNAs

Spermatozoa harbour a complex and environment-sensitive pool of small non-coding RNAs (sncRNAs)1, which influences offspring development and adult phenotypes1-7. Whether spermatozoa in the epididymis are directly susceptible to environmental cues is not fully understood8. Here we used two distinct paradigms of preconception acute high-fat diet to dissect epididymal versus testicular contributions to the sperm sncRNA pool and offspring health. We show that epididymal spermatozoa, but not developing germ cells, are sensitive to the environment and identify mitochondrial tRNAs (mt-tRNAs) and their fragments (mt-tsRNAs) as sperm-borne factors. In humans, mt-tsRNAs in spermatozoa correlate with body mass index, and paternal overweight at conception doubles offspring obesity risk and compromises metabolic health. Sperm sncRNA sequencing of mice mutant for genes involved in mitochondrial function, and metabolic phenotyping of their wild-type offspring, suggest that the upregulation of mt-tsRNAs is downstream of mitochondrial dysfunction. Single-embryo transcriptomics of genetically hybrid two-cell embryos demonstrated sperm-to-oocyte transfer of mt-tRNAs at fertilization and suggested their involvement in the control of early-embryo transcription. Our study supports the importance of paternal health at conception for offspring metabolism, shows that mt-tRNAs are diet-induced and sperm-borne and demonstrates, in a physiological setting, father-to-offspring transfer of sperm mitochondrial RNAs at fertilization.

Loss of the mitochondrial carrier, SLC25A1, during embryogenesis induces a unique senescence program controlled by p53

Germline inactivating mutations of the SLC25A1 gene contribute to various human developmental disorders, including combined D/L-2-hydroxyglutaric aciduria (D/L-2HGA), a severe systemic syndrome characterized by the accumulation of both enantiomers of 2-hydroxyglutaric acid (2HG). The mechanisms by which SLC25A1 deficiency leads to this disease and the role of 2HG are unclear and no therapies exist. We now show that mice lacking both Slc25a1 alleles display a spectrum of alterations that resemble human D/L-2HGA. Mechanistically, SLC25A1 loss results in a proliferation defect and activates two distinct senescence pathways, oncogene-induced senescence (OIS) and mitochondrial dysfunction-induced senescence (MiDAS), both involving the p53 tumor suppressor and driven by two discernible signals: the accumulation of 2HG, inducing OIS, and mitochondrial dysfunction, triggering MiDAS. Inhibiting these senescence programs or blocking p53 activity reverses the growth defect caused by SLC25A1 dysfunction and restores proliferation. These findings reveal novel pathogenic roles of senescence in human disorders and suggest potential strategies to correct the molecular alterations caused by SLC25A1 loss.

Knockout mice with pituitary malformations help identify human cases of hypopituitarism

BACKGROUND

Congenital hypopituitarism (CH) and its associated syndromes, septo-optic dysplasia (SOD) and holoprosencephaly (HPE), are midline defects that cause significant morbidity for affected people. Variants in 67 genes are associated with CH, but a vast majority of CH cases lack a genetic diagnosis. Whole exome and whole genome sequencing of CH patients identifies sequence variants in genes known to cause CH, and in new candidate genes, but many of these are variants of uncertain significance (VUS).

METHODS

The International Mouse Phenotyping Consortium (IMPC) is an effort to establish gene function by knocking-out all genes in the mouse genome and generating corresponding phenotype data. We used mouse embryonic imaging data generated by the Deciphering Mechanisms of Developmental Disorders (DMDD) project to screen 209 embryonic lethal and sub-viable knockout mouse lines for pituitary malformations.

RESULTS

Of the 209 knockout mouse lines, we identified 51 that have embryonic pituitary malformations. These genes not only represent new candidates for CH, but also reveal new molecular pathways not previously associated with pituitary organogenesis. We used this list of candidate genes to mine whole exome sequencing data of a cohort of patients with CH, and we identified variants in two unrelated cases for two genes, MORC2 and SETD5, with CH and other syndromic features.

CONCLUSIONS

The screening and analysis of IMPC phenotyping data provide proof-of-principle that recessive lethal mouse mutants generated by the knockout mouse project are an excellent source of candidate genes for congenital hypopituitarism in children.

Ancient eukaryotic protein interactions illuminate modern genetic traits and disorders

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Causal network perturbation analysis identifies known and novel type-2 diabetes driver genes

The molecular pathogenesis of diabetes is multifactorial, involving genetic predisposition and environmental factors that are not yet fully understood. However, pancreatic β-cell failure remains among the primary reasons underlying the progression of type-2 diabetes (T2D) making targeting β-cell dysfunction an attractive pathway for diabetes treatment. To identify genetic contributors to β-cell dysfunction, we investigated single-cell gene expression changes in β-cells from healthy (C57BL/6J) and diabetic (NZO/HlLtJ) mice fed with normal or high-fat, high-sugar diet (HFHS). Our study presents an innovative integration of the causal network perturbation assessment (ssNPA) framework with meta-cell transcriptome analysis to explore the genetic underpinnings of type-2 diabetes (T2D). By generating a reference causal network and in silico perturbation, we identified novel genes implicated in T2D and validated our candidates using the Knockout Mouse Phenotyping (KOMP) Project database.

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

Methods that analyze single-cell paired RNA-seq and ATAC-seq multiome data have shown great promise in linking regulatory elements to genes. However, existing methods differ in their modeling assumptions and approaches to account for biological and technical noise-leading to low concordance in their linking scores-and do not capture the effects of genomic distance. We propose pgBoost, an integrative modeling framework that trains a non-linear combination of existing linking strategies (including genomic distance) on fine-mapped eQTL data to assign a probabilistic score to each candidate SNP-gene link. We applied pgBoost to single-cell multiome data from 85k cells representing 6 major immune/blood cell types. pgBoost attained higher enrichment for fine-mapped eSNP-eGene pairs (e.g. 21x at distance >10kb) than existing methods (1.2-10x; p-value for difference = 5e-13 vs. distance-based method and < 4e-35 for each other method), with larger improvements at larger distances (e.g. 35x vs. 0.89-6.6x at distance >100kb; p-value for difference < 0.002 vs. each other method). pgBoost also outperformed existing methods in enrichment for CRISPR-validated links (e.g. 4.8x vs. 1.6-4.1x at distance >10kb; p-value for difference = 0.25 vs. distance-based method and < 2e-5 for each other method), with larger improvements at larger distances (e.g. 15x vs. 1.6-2.5x at distance >100kb; p-value for difference < 0.009 for each other method). Similar improvements in enrichment were observed for links derived from Activity-By-Contact (ABC) scores and GWAS data. We further determined that restricting pgBoost to features from a focal cell type improved the identification of SNP-gene links relevant to that cell type. We highlight several examples where pgBoost linked fine-mapped GWAS variants to experimentally validated or biologically plausible target genes that were not implicated by other methods. In conclusion, a non-linear combination of linking strategies, including genomic distance, improves power to identify target genes underlying GWAS associations.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.