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

Shieldin and CST co-orchestrate DNA polymerase-dependent tailed-end joining reactions independently of 53BP1-governed repair pathway choice

Tumor suppressor p53-binding protein 1 (53BP1) regulates DNA end joining in lymphocytes, diversifying immune antigen receptors. This involves nucleosome-bound 53BP1 at DNA double-stranded breaks (DSBs) recruiting Rap1-interacting factor 1 homolog (RIF1) and shieldin, a poorly understood DNA-binding complex. The 53BP1-RIF1-shieldin axis is pathological in BRCA1-mutated cancers, blocking homologous recombination (HR) and driving illegitimate nonhomologous end joining (NHEJ). However, how this axis regulates DNA end joining and HR suppression remains unresolved. We investigated shieldin and its interplay with the Ctc1-Stn1-Ten1 (CST) complex, which was recently implicated downstream of 53BP1. Immunophenotypically, mice lacking shieldin or CST are equivalent, with class-switch recombination coreliant on both complexes. Ataxia-telangiectasia mutated kinase-dependent DNA damage signaling underpins this cooperation, inducing physical interactions between these complexes that reveal shieldin as a DSB-responsive CST adaptor. Furthermore, DNA polymerase ζ functions downstream of shieldin, establishing DNA fill-in synthesis as the physiological function of shieldin-CST. Lastly, we demonstrate that 53BP1 suppresses HR and promotes NHEJ in BRCA1-deficient mice and cells independently of shieldin. These findings showcase the versatility of the 53BP1 pathway, achieved through the collaboration of chromatin-bound 53BP1 complexes and DNA end-processing effector proteins.

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.

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

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/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.

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.

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 [...].

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.

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.

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.