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

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.

Loss of cped1 does not affect bone and lean mass 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 difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we turned to zebrafish, an emerging model for orthopaedic 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, and discuss how these differences, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies demonstrate that cped1 is not required for normal adult zebrafish bone mass, lean mass, or bone and lean mass morphology, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

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.

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.

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.

Burden re-analysis of neurodevelopmental disorder cohorts for prioritization of candidate genes

This study aimed to uncover novel genes associated with neurodevelopmental disorders (NDD) by leveraging recent large-scale de novo burden analysis studies to enhance a virtual gene panel used in a diagnostic setting. We re-analyzed historical trio-exome sequencing data from 745 individuals with NDD according to the most recent diagnostic standards, resulting in a cohort of 567 unsolved individuals. Next, we designed a virtual gene panel containing candidate genes from three large de novo burden analysis studies in NDD and prioritized candidate genes by stringent filtering for ultra-rare de novo variants with high pathogenicity scores. Our analysis revealed an increased burden of de novo variants in our selected candidate genes within the unsolved NDD cohort and identified qualifying de novo variants in seven candidate genes: RIF1, CAMK2D, RAB11FIP4, AGO3, PCBP2, LEO1, and VCP. Clinical data were collected from six new individuals with de novo or inherited LEO1 variants and three new individuals with de novo PCBP2 variants. Our findings add additional evidence for LEO1 as a risk gene for autism and intellectual disability. Furthermore, we prioritize PCBP2 as a candidate gene for NDD associated with motor and language delay. In summary, by leveraging de novo burden analysis studies, employing a stringent variant filtering pipeline, and engaging in targeted patient recruitment, our study contributes to the identification of novel genes implicated in NDDs.

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.

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.

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.

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.

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.

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.

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.

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.

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.

CRISPR-Cas Systems and Genome Editing: Beginning the Era of CRISPR/Cas Therapies for Humans

Harnessing of CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated genes) systems for detection, chemical modification, and sequence editing of nucleic acids dramatically changed many fields of fundamental science, biotechnology, and biomedicine [...].

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.

Mouse Genome Informatics: an integrated knowledgebase system for the laboratory mouse

Mouse Genome Informatics (MGI) is a federation of expertly curated information resources designed to support experimental and computational investigations into genetic and genomic aspects of human biology and disease using the laboratory mouse as a model system. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are core MGI databases that share data and system architecture. MGI serves as the central community resource of integrated information about mouse genome features, variation, expression, gene function, phenotype, and human disease models acquired from peer-reviewed publications, author submissions, and major bioinformatics resources. To facilitate integration and standardization of data, biocuration scientists annotate using terms from controlled metadata vocabularies and biological ontologies (e.g. Mammalian Phenotype Ontology, Mouse Developmental Anatomy, Disease Ontology, Gene Ontology, etc.), and by applying international community standards for gene, allele, and mouse strain nomenclature. MGI serves basic scientists, translational researchers, and data scientists by providing access to FAIR-compliant data in both human-readable and compute-ready formats. The MGI resource is accessible at https://informatics.jax.org. Here, we present an overview of the core data types represented in MGI and highlight recent enhancements to the resource with a focus on new data and functionality for MGD and GXD.

VARista: a free web platform for streamlined whole-genome variant analysis across T2T, hg38, and hg19

With the increasing importance of genomic data in understanding genetic diseases, there is an essential need for efficient and user-friendly tools that simplify variant analysis. Although multiple tools exist, many present barriers such as steep learning curves, limited reference genome compatibility, or costs. We developed VARista, a free web-based tool, to address these challenges and provide a streamlined solution for researchers, particularly those focusing on rare monogenic diseases. VARista offers a user-centric interface that eliminates much of the technical complexity typically associated with variant analysis. The tool directly supports VCF files generated using reference genomes hg19, hg38, and the emerging T2T, with seamless remapping capabilities between them. Features such as gene summaries and links, tissue and cell-specific gene expression data for both adults and fetuses, as well as automated PCR design and integration with tools such as SpliceAI and AlphaMissense, enable users to focus on the biology and the case itself. As we demonstrate, VARista proved effective in narrowing down potential disease-causing variants, prioritizing them effectively, and providing meaningful biological context, facilitating rapid decision-making. VARista stands out as a freely available and comprehensive tool that consolidates various aspects of variant analysis into a single platform that embraces the forefront of genomic advancements. Its design inherently supports a shift in focus from technicalities to critical thinking, thereby promoting better-informed decisions in genetic disease research. Given its unique capabilities and user-centric design, VARista has the potential to become an essential asset for the genomic research community. https://VARista.link.

Deficiency of FRMD5 results in neurodevelopmental dysfunction and autistic-like behavior in mice

The pathophysiology of autism spectrum disorders (ASDs) is causally linked to postsynaptic scaffolding proteins, as evidenced by numerous large-scale genomic studies [1, 2] and in vitro and in vivo neurobiological studies of mutations in animal models [3, 4]. However, due to the distinct phenotypic and genetic heterogeneity observed in ASD patients, individual mutation genes account for only a small proportion (<2%) of cases [1, 5]. Recently, a human genetic study revealed a correlation between de novo variants in FERM domain-containing-5 (FRMD5) and neurodevelopmental abnormalities [6]. In this study, we demonstrate that deficiency of the scaffolding protein FRMD5 leads to neurodevelopmental dysfunction and ASD-like behavior in mice. FRMD5 deficiency results in morphological abnormalities in neurons and synaptic dysfunction in mice. Frmd5-deficient mice display learning and memory dysfunction, impaired social function, and increased repetitive stereotyped behavior. Mechanistically, tandem mass tag (TMT)-labeled quantitative proteomics revealed that FRMD5 deletion affects the distribution of synaptic proteins involved in the pathological process of ASD. Collectively, our findings delineate the critical role of FRMD5 in neurodevelopment and ASD pathophysiology, suggesting potential therapeutic implications for the treatment of ASD.

Bridge-like lipid transfer protein family member 2 suppresses ciliogenesis

Bridge-like lipid transfer protein family member 2 (BLTP2) is an evolutionary conserved protein with unknown function(s). The absence of BLTP2 in Drosophila melanogaster results in impaired cellular secretion and larval death, while in mice (Mus musculus), it causes preweaning lethality. Structural predictions propose that BLTP2 belongs to the repeating β-groove domain-containing (also called the VPS13) protein family, forming a long tube with a hydrophobic core, suggesting that it operates as a lipid transfer protein (LTP). We establish BLTP2 as a negative regulator of ciliogenesis in RPE-1 cells based on a strong genetic interaction with WDR44, a gene that also suppresses ciliogenesis. Like WDR44, BLTP2 localizes to membrane contact sites involving the endoplasmic reticulum and the tubular endosome network in HeLa cells and that BLTP2 depletion enhanced ciliogenesis in RPE-1 cells grown in serum-containing medium, a condition where ciliogenesis is normally suppressed. This study establishes human BLTP2 as a putative LTP acting between tubular endosomes and ER that regulates primary cilium biogenesis.

The role of prickle proteins in vertebrate development and pathology

Prickle is an evolutionarily conserved family of proteins exclusively associated with planar cell polarity (PCP) signalling. This signalling pathway provides directional and positional cues to eukaryotic cells along the plane of an epithelial sheet, orthogonal to both apicobasal and left-right axes. Through studies in the fruit fly Drosophila, we have learned that PCP signalling is manifested by the spatial segregation of two protein complexes, namely Prickle/Vangl and Frizzled/Dishevelled. While Vangl, Frizzled, and Dishevelled proteins have been extensively studied, Prickle has been largely neglected. This is likely because its role in vertebrate development and pathologies is still being explored and is not yet fully understood. The current review aims to address this gap by summarizing our current knowledge on vertebrate Prickle proteins and to cover their broad versatility. Accumulating evidence suggests that Prickle is involved in many developmental events, contributes to homeostasis, and can cause diseases when its expression and signalling properties are deregulated. This review highlights the importance of Prickle in vertebrate development, discusses the implications of Prickle-dependent signalling in pathology, and points out the blind spots or potential links regarding Prickle, which could be studied further.

Dissecting genetic architecture of rare dystonia: genetic, molecular and clinical insights

BACKGROUND

Dystonia is one of the most common movement disorders. To date, the genetic causes of dystonia in populations of European descent have been extensively studied. However, other populations, particularly those from the Middle East, have not been adequately studied. The purpose of this study is to discover the genetic basis of dystonia in a clinically and genetically well-characterised dystonia cohort from Turkey, which harbours poorly studied populations.

METHODS

Exome sequencing analysis was performed in 42 Turkish dystonia families. Using co-expression network (CEN) analysis, identified candidate genes were interrogated for the networks including known dystonia-associated genes and genes further associated with the protein-protein interaction, animal model-based characteristics and clinical findings.

RESULTS

We identified potentially disease-causing variants in the established dystonia genes (PRKRA, SGCE, KMT2B, SLC2A1, GCH1, THAP1, HPCA, TSPOAP1, AOPEP; n=11 families (26%)), in the uncommon forms of dystonia-associated genes (PCCB, CACNA1A, ALDH5A1, PRKN; n=4 families (10%)) and in the candidate genes prioritised based on the pathogenicity of the variants and CEN-based analyses (n=11 families (21%)). The diagnostic yield was found to be 36%. Several pathways and gene ontologies implicated in immune system, transcription, metabolic pathways, endosomal-lysosomal and neurodevelopmental mechanisms were over-represented in our CEN analysis.

CONCLUSIONS

Here, using a structured approach, we have characterised a clinically and genetically well-defined dystonia cohort from Turkey, where dystonia has not been widely studied, and provided an uncovered genetic basis, which will facilitate diagnostic dystonia research.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Mark4 ablation attenuates pathological phenotypes in a mouse model of tauopathy

Accumulation of abnormally phosphorylated tau proteins is linked to various neurodegenerative diseases, including Alzheimer's disease and frontotemporal dementia. Microtubule affinity-regulating kinase 4 (MARK4) has been genetically and pathologically associated with Alzheimer's disease and reported to enhance tau phosphorylation and toxicity in Drosophila and mouse traumatic brain-injury models but not in mammalian tauopathy models. To investigate the role of MARK4 in tau-mediated neuropathology, we crossed P301S tauopathy model (PS19) and Mark4 knockout mice. We performed behaviour, biochemical and histology analyses to evaluate changes in PS19 pathological phenotype with and without Mark4. Here, we demonstrated that Mark4 deletion ameliorated the tau pathology in a mouse model of tauopathy. In particular, we found that PS19 with Mark4 knockout showed improved mortality and memory compared with those bearing an intact Mark4 gene. These phenotypes were accompanied by reduced neurodegeneration and astrogliosis in response to the reduction of pathological forms of tau, such as those phosphorylated at Ser356, AT8-positive tau and thioflavin S-positive tau. Our data indicate that MARK4 critically contributes to tau-mediated neuropathology, suggesting that MARK4 inhibition may serve as a therapeutic avenue for tauopathies.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Age-Invariant Genes: Multi-Tissue Identification and Characterization of Murine Reference Genes

Studies of the aging transcriptome focus on genes that change with age. But what can we learn from age-invariant genes-those that remain unchanged throughout the aging process? These genes also have a practical application: they serve as reference genes (often called housekeeping genes) in expression studies. Reference genes have mostly been identified and validated in young organisms, and no systematic investigation has been done across the lifespan. Here, we build upon a common pipeline for identifying reference genes in RNA-seq datasets to identify age-invariant genes across seventeen C57BL/6 mouse tissues (brain, lung, bone marrow, muscle, white blood cells, heart, small intestine, kidney, liver, pancreas, skin, brown, gonadal, marrow, and subcutaneous adipose tissue) spanning 1 to 21+ months of age. We identify 9 pan-tissue age-invariant genes and many tissue-specific age-invariant genes. These genes are stable across the lifespan and are validated in independent bulk RNA-seq datasets and RT-qPCR. We find age-invariant genes have shorter transcripts on average and are enriched for CpG islands. Interestingly, pathway enrichment analysis for age-invariant genes identifies an overrepresentation of molecular functions associated with some, but not all, hallmarks of aging. Thus, though hallmarks of aging typically involve changes in cell maintenance mechanisms, select genes associated with these hallmarks resist fluctuations in expression with age. Finally, our analysis concludes no classical reference gene is appropriate for aging studies in all tissues. Instead, we provide tissue-specific and pan-tissue genes for assays utilizing reference gene normalization (i.e., RT-qPCR) that can be applied to animals across the lifespan.

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 constraint for the shortest ∼25% of genes, potentially causing important pathogenic mutations to be overlooked. 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 new 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 estimation of many gene-level properties, such as rare variant burden or gene expression differences.

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.

Effects of Tulp4 deficiency on murine embryonic development and adult phenotype

Genetically engineered mouse models have the potential to unravel fundamental biological processes and provide mechanistic insights into the pathogenesis of human diseases. We have previously observed that germline genetic variation at the TULP4 locus influences clinical characteristics in patients with myeloproliferative neoplasms. To elucidate the role of TULP4 in pathological and physiological processes in vivo, we generated a Tulp4 knockout mouse model. Systemic Tulp4 deficiency exerted a strong impact on embryonic development in both Tulp4 homozygous null (Tulp4-/-) and heterozygous (Tulp4+/-) knockout mice, the former exhibiting perinatal lethality. High-resolution episcopic microscopy (HREM) of day 14.5 embryos allowed for the identification of multiple developmental defects in Tulp4-/- mice, including severe heart defects. Moreover, in Tulp4+/- embryos HREM revealed abnormalities of several organ systems, which per se do not affect prenatal or postnatal survival. In adult Tulp4+/- mice, extensive examinations of hematopoietic and cardiovascular features, involving histopathological surveys of multiple tissues as well as blood counts and immunophenotyping, did not provide evidence for anomalies as observed in corresponding embryos. Finally, evaluating a potential obesity-related phenotype as reported for other TULP family members revealed a trend for increased body weight of Tulp4+/- mice. RESEARCH HIGHLIGHTS: To study the role of the TULP4 gene in vivo, we generated a Tulp4 knockout mouse model. Correlative analyses involving HREM revealed a strong impact of Tulp4 deficiency on murine embryonic development.

Improving prenatal diagnosis through standards and aggregation

Advances in sequencing and imaging technologies enable enhanced assessment in the prenatal space, with a goal to diagnose and predict the natural history of disease, to direct targeted therapies, and to implement clinical management, including transfer of care, election of supportive care, and selection of surgical interventions. The current lack of standardization and aggregation stymies variant interpretation and gene discovery, which hinders the provision of prenatal precision medicine, leaving clinicians and patients without an accurate diagnosis. With large amounts of data generated, it is imperative to establish standards for data collection, processing, and aggregation. Aggregated and homogeneously processed genetic and phenotypic data permits dissection of the genomic architecture of prenatal presentations of disease and provides a dataset on which data analysis algorithms can be tuned to the prenatal space. Here we discuss the importance of generating aggregate data sets and how the prenatal space is driving the development of interoperable standards and phenotype-driven tools.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

Purpose

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 pedigrees 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/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), 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.

Decreased calcium permeability caused by biallelic TRPV5 mutation leads to autosomal recessive renal calcium-wasting hypercalciuria

Hypercalciuria is the most common metabolic risk factor in people with kidney stone disease. Its etiology is mostly multifactorial, although monogenetic causes of hypercalciuria have also been described. Despite the increased availability of genetic diagnostic tests, the vast majority of individuals with familial hypercalciuria remain unsolved. In this study, we investigated a consanguineous pedigree with idiopathic hypercalciuria. The proband additionally exhibited severe skeletal deformities and hyperparathyroidism. Whole-exome sequencing of the proband revealed a homozygous ultra-rare variant in TRPV5 (NM_019841.7:c.1792G>A; p.(Val598Met)), which encodes for a renal Ca2+-selective ion channel. The variant segregates with the three individuals with hypercalciuria. The skeletal phenotype unique to the proband was due to an additional pathogenic somatic mutation in GNAS (NM_000516.7:c.601C>T; p.(Arg201Cys)), which leads to polyostotic fibrous dysplasia. The variant in TRPV5 is located in the TRP helix, a characteristic amphipathic helix that is indispensable for the gating movements of TRP channels. Biochemical characterization of the TRPV5 p.(Val598Met) channel revealed a complete loss of Ca2+ transport capability. This defect is caused by reduced expression of the mutant channel, due to misfolding and preferential targeting to the proteasome for degradation. Based on these findings, we conclude that biallelic loss of TRPV5 function causes a novel form of monogenic autosomal recessive hypercalciuria, which we name renal Ca2+-wasting hypercalciuria (RCWH). The recessive inheritance pattern explains the rarity of RCWH and underscores the potential prevalence of RCWH in highly consanguineous populations, emphasizing the importance of exploration of this disorder within such communities.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

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 genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Evidence for the druggability of aldosterone targets in heart failure: A bioinformatics and data science-driven decision-making approach

BACKGROUND

Aldosterone plays a key role in the neurohormonal drive of heart failure. Systematic prioritization of drug targets using bioinformatics and database-driven decision-making can provide a competitive advantage in therapeutic R&D. This study investigated the evidence on the druggability of these aldosterone targets in heart failure.

METHODS

The target disease predictability of mineralocorticoid receptors (MR) and aldosterone synthase (AS) in cardiac failure was evaluated using Open Targets target-disease association scores. The Open Targets database collections were downloaded to MongoDB and queried according to the desired aggregation level, and the results were retrieved from the Europe PMC (data type: text mining), ChEMBL (data type: drugs), Open Targets Genetics Portal (data type: genetic associations), and IMPC (data type: genetic associations) databases. The target tractability of MR and AS in the cardiovascular system was investigated by computing activity scores in a curated ChEMBL database using supervised machine learning.

RESULTS

The medians of the association scores of the MR and AS groups were similar, indicating a comparable predictability of the target disease. The median of the MR activity scores group was significantly lower than that of AS, indicating that AS has higher target tractability than MR [Hodges-Lehmann difference 0.62 (95%CI 0.53-0.70, p < 0.0001]. The cumulative distributions of the overall multiplatform association scores of cardiac diseases with MR were considerably higher than with AS, indicating more advanced investigations on a wider range of disorders evaluated for MR (Kolmogorov-Smirnov D = 0.36, p = 0.0009). In curated ChEMBL, MR had a higher cumulative distribution of activity scores in experimental cardiovascular assays than AS (Kolmogorov-Smirnov D = 0.23, p < 0.0001). Documented clinical trials for MR in heart failures surfaced in database searches, none for AS.

CONCLUSIONS

Although its clinical development has lagged behind that of MR, our findings indicate that AS is a promising therapeutic target for the treatment of cardiac failure. The multiplatform-integrated identification used in this study allowed us to comprehensively explore the available scientific evidence on MR and AS for heart failure therapy.

A transient protein folding response targets aggregation in the early phase of TDP-43-mediated neurodegeneration

Understanding the mechanisms that drive TDP-43 pathology is integral to combating amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration (FTLD) and other neurodegenerative diseases. Here we generated a longitudinal quantitative proteomic map of the cortex from the cytoplasmic TDP-43 rNLS8 mouse model of ALS and FTLD, and developed a complementary open-access webtool, TDP-map ( https://shiny.rcc.uq.edu.au/TDP-map/ ). We identified distinct protein subsets enriched for diverse biological pathways with temporal alterations in protein abundance, including increases in protein folding factors prior to disease onset. This included increased levels of DnaJ homolog subfamily B member 5, DNAJB5, which also co-localized with TDP-43 pathology in diseased human motor cortex. DNAJB5 over-expression decreased TDP-43 aggregation in cell and cortical neuron cultures, and knockout of Dnajb5 exacerbated motor impairments caused by AAV-mediated cytoplasmic TDP-43 expression in mice. Together, these findings reveal molecular mechanisms at distinct stages of ALS and FTLD progression and suggest that protein folding factors could be protective in neurodegenerative diseases.

A metabolic signature for NADSYN1-dependent congenital NAD deficiency disorder

Nicotinamide adenine dinucleotide (NAD) is essential for embryonic development. To date, biallelic loss-of-function variants in 3 genes encoding nonredundant enzymes of the NAD de novo synthesis pathway - KYNU, HAAO, and NADSYN1 - have been identified in humans with congenital malformations defined as congenital NAD deficiency disorder (CNDD). Here, we identified 13 further individuals with biallelic NADSYN1 variants predicted to be damaging, and phenotypes ranging from multiple severe malformations to the complete absence of malformation. Enzymatic assessment of variant deleteriousness in vitro revealed protein domain-specific perturbation, complemented by protein structure modeling in silico. We reproduced NADSYN1-dependent CNDD in mice and assessed various maternal NAD precursor supplementation strategies to prevent adverse pregnancy outcomes. While for Nadsyn1+/- mothers, any B3 vitamer was suitable to raise NAD, preventing embryo loss and malformation, Nadsyn1-/- mothers required supplementation with amidated NAD precursors (nicotinamide or nicotinamide mononucleotide) bypassing their metabolic block. The circulatory NAD metabolome in mice and humans before and after NAD precursor supplementation revealed a consistent metabolic signature with utility for patient identification. Our data collectively improve clinical diagnostics of NADSYN1-dependent CNDD, provide guidance for the therapeutic prevention of CNDD, and suggest an ongoing need to maintain NAD levels via amidated NAD precursor supplementation after birth.

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures

BACKGROUND

Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science.

RESULTS

Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice.

CONCLUSIONS

C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.