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

Identification of a global gene expression signature associated with the genetic risk of catastrophic fracture in iPSC-derived osteoblasts from Thoroughbred horses

Bone fractures are a significant problem in Thoroughbred racehorses. The risk of fracture is influenced by both genetic and environmental factors. To determine the biological processes that are affected in genetically susceptible horses, we utilised polygenic risk scoring to establish induced pluripotent stem cells (iPSCs) from horses at high and low genetic risk. RNA-sequencing on iPSC-derived osteoblasts revealed 112 genes that were significantly differentially expressed. Forty-three of these genes have known roles in bone, 27 are not yet annotated in the equine genome and 42 currently have no described role in bone. However, many of the proteins encoded by the known and unknown genes have reported interactions. Functional enrichment analyses revealed that the differentially expressed genes were overrepresented in processes regulating the extracellular matrix and pathways known to be involved in bone remodelling and bone diseases. Gene set enrichment analysis also detected numerous biological processes and pathways involved in glycolysis with the associated genes having a higher expression in the iPSC-osteoblasts from horses with low polygenic risk scores for fracture. Therefore, the differentially expressed genes may be relevant for maintaining bone homeostasis and contribute to fracture risk. A deeper understanding of the consequences of mis-regulation of these genes and the identification of the DNA variants which underpin their differential expression may reveal more about the molecular mechanisms which are involved in equine bone health and fracture risk.

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.

Potential Adaptive Introgression From Dogs in Iberian Grey Wolves (Canis lupus)

Invading species along with increased anthropogenization may lead to hybridization events between wild species and closely related domesticates. As a consequence, wild species may carry introgressed alleles from domestic species, which is generally assumed to yield adverse effects in wild populations. The opposite evolutionary consequence, adaptive introgression, where introgressed genes are positively selected in the wild species, is possible but has rarely been documented. Grey wolves (Canis lupus) are widely distributed across the Holarctic and frequently coexist with their close relative, the domestic dog (C. familiaris). Despite ample opportunity, hybridization rarely occurs in most populations. Here we studied the geographically isolated grey wolves of the Iberian Peninsula, who have coexisted with a large population of loosely controlled dogs for thousands of years in a human-modified landscape. We assessed the extent and impact of dog introgression on the current Iberian grey wolf population by analysing 150 whole genomes of Iberian and other Eurasian grey wolves as well as dogs originating from across Europe and western Siberia. We identified almost no recent introgression and a small (< 5%) overall ancient dog ancestry. Using a combination of single scan statistics and ancestry enrichment estimates, we identified positive selection on six genes (DAPP1, NSMCE4A, MPPED2, PCDH9, MBTPS1, and CDH13) for which wild Iberian wolves carry alleles introgressed from dogs. The genes with introgressed and positively selected alleles include functions in immune response and brain functions, which may explain some of the unique behavioural phenotypes in Iberian wolves such as their reduced dispersal compared to other wolf populations.

Open Targets Platform: facilitating therapeutic hypotheses building in drug discovery

The Open Targets Platform (https://platform.opentargets.org) is a unique, open-source, publicly-available knowledge base providing data and tooling for systematic drug target identification, annotation, and prioritisation. Since our last report, we have expanded the scope of the Platform through a number of significant enhancements and data updates, with the aim to enable our users to formulate more flexible and impactful therapeutic hypotheses. In this context, we have completely revamped our target-disease associations page with more interactive facets and built-in functionalities to empower users with additional control over their experience using the Platform, and added a new Target Prioritisation view. This enables users to prioritise targets based upon clinical precedence, tractability, doability and safety attributes. We have also implemented a direction of effect assessment for eight sources of target-disease association evidence, showing the effect of genetic variation on the function of a target is associated with risk or protection for a trait to inform on potential mechanisms of modulation suitable for disease treatment. These enhancements and the introduction of new back and front-end technologies to support them have increased the impact and usability of our resource within the drug discovery community.

Harmonizome 3.0: integrated knowledge about genes and proteins from diverse multi-omics resources

By processing and abstracting diverse omics datasets into associations between genes and their attributes, the Harmonizome database enables researchers to explore and integrate knowledge about human genes from many central omics resources. Here, we introduce Harmonizome 3.0, a significant upgrade to the original Harmonizome database. The upgrade adds 26 datasets that contribute nearly 12 million associations between genes and various attribute types such as cells and tissues, diseases, and pathways. The upgrade has a dataset crossing feature to identify gene modules that are shared across datasets. To further explain significantly high gene set overlap between dataset pairs, a large language model (LLM) composes a paragraph that speculates about the reasons behind the high overlap. The upgrade also adds more data formats and visualization options. Datasets are downloadable as knowledge graph (KG) assertions and visualized with Uniform Manifold Approximation and Projection (UMAP) plots. The KG assertions can be explored via a user interface that visualizes gene-attribute associations as ball-and-stick diagrams. Overall, Harmonizome 3.0 is a rich resource of processed omics datasets that are provided in several AI-ready formats. Harmonizome 3.0 is available at https://maayanlab.cloud/Harmonizome/.

Targeting the Hippo pathway in cancer: kidney toxicity as a class effect of TEAD inhibitors?

The Hippo pathway has emerged as a critical player in both cancers and targeted therapy resistance. Recent drug discovery efforts have led to the development of TEAD inhibitors, several of which have already progressed to the clinic. To truly leverage their potential as anticancer therapeutics, safety considerations, particularly in regard to the kidney, warrant additional investigation. This review explores the Hippo pathway's role in cancers, its therapeutic potential, role in kidney development, and the need to evaluate the best strategies to translate its clinical application for long-term patient benefit.

Inactivation of the SLC25A1 gene during embryogenesis induces a unique senescence program controlled by p53

Germline inactivating mutations of the SLC25A1 gene contribute to various human disorders, including Velocardiofacial (VCFS), DiGeorge (DGS) syndromes and combined D/L-2-hydroxyglutaric aciduria (D/L-2HGA), a severe systemic disease characterized by the accumulation of 2-hydroxyglutaric acid (2HG). The mechanisms by which SLC25A1 loss leads to these syndromes remain largely unclear. Here, we describe a mouse model of SLC25A1 deficiency that mimics human VCFS/DGS and D/L-2HGA. Surprisingly, inactivation of both Slc25a1 alleles results in alterations in the development of multiple organs, and in a severe proliferation defect by activating two senescence programs, oncogene-induced senescence (OIS) and mitochondrial dysfunction-induced senescence (MiDAS), which converge upon the induction of the p53 tumor suppressor. Mechanistically, cells and tissues with dysfunctional SLC25A1 protein undergo metabolic and transcriptional rewiring leading to the accumulation of 2HG via a non-canonical pathway and to the depletion of nicotinamide adenine dinucleotide, NAD+, which trigger senescence. Replenishing the pool of NAD+ or promoting the clearance of 2HG rescues the proliferation defect of cells with dysfunctional SLC25A1 in a cooperative fashion. Further, removal of p53 activity via RNA interference restores proliferation, indicating that p53 acts as a critical barrier to the expansion of cells lacking functional SLC25A1. These findings reveal unexpected pathogenic roles of senescence and of p53 in D/L-2HGA and identify potential therapeutic strategies to correct salient molecular alterations driving this disease.

Identification of key genes and variants associated with boar sperm freezability using whole genome resequencing

Cryopreservation induces various cryodamages to the structural or functional aspects of boar sperm, resulting in the deterioration of sperm quality. The extent of cryodamages varies significantly among different individual boars. In our study, 50 boars with either good sperm freezability (GSF) or poor sperm freezability (PSF) were selected from a population of 402 boars. These two groups exhibited significant differences in sperm quality, acrosome integrity, mitochondrial membrane potential, and plasma membrane integrity. Subsequent whole-genome resequencing and FST analysis of GSF and PSF boars uncovered genetic differentiation among 5632 genes. Spermatogenesis, sperm structure, and lipid composition were identified as the potential factors, and 35 genes, such as RNF8, PACRG, CADM1, SPAG16, PPP3CA and MFF, were identified as key candidate genes associated with sperm freezability. Furthermore, by using MassARRAY genotyping and a general linear model, nine variants, including seven single nucleotide polymorphisms (SNPs) and two insertion variants, were found to be significantly correlated with sperm freezability. These variants, located within the candidate genes, represent potential molecular markers for sperm freezability. Our research contributes novel insights into the genetic mechanisms underlying freezability differences in boar sperm and identifies candidate genes and key molecular markers that could be applied in boar breeding practices.

Intersectin1 promotes clathrin-mediated endocytosis by organizing and stabilizing endocytic protein interaction networks

During clathrin-mediated endocytosis (CME), dozens of proteins are recruited to nascent CME sites on the plasma membrane, and their spatial and temporal coordination is crucial for efficient CME. Here, we show that the scaffold protein intersectin1 (ITSN1) promotes CME by organizing and stabilizing endocytic protein interaction networks. Live-cell imaging of genome-edited cells revealed that endogenously labeled ITSN1 is recruited during CME site stabilization and growth and that ITSN1 knockdown impairs endocytic protein recruitment during this stage. Targeting ITSN1 to the mitochondrial surface was sufficient to assemble puncta consisting of the EPS15 and FCHO2 initiation proteins, the AP2 and epsin1 (EPN1) adaptor proteins, and the dynamin2 (DNM2) vesicle scission GTPase. ITSN1 can form puncta and recruit DNM2 independent of EPS15/FCHO2 or EPN1. Our findings redefine ITSN1's primary endocytic role as organizing and stabilizing CME protein interaction networks rather than initiation, providing deeper insights into the multi-step and multi-zone organization of CME site assembly.

DNAJC13 localization to endosomes is opposed by its J domain and its disordered C-terminal tail

Endosomes are a central sorting hub for membrane cargos. DNAJC13/RME-8 plays a critical role in endosomal trafficking by regulating the endosomal recycling or degradative pathways. DNAJC13 localizes to endosomes through its N-terminal Plekstrin Homology (PH)-like domain, which directly binds endosomal phosphoinositol-3-phosphate (PI(3)P). However, little is known about how DNAJC13 localization is regulated. Here, we show that two regions within DNAJC13, its J domain and disordered C-terminal tail, act as negative regulators of its PH-like domain. Using a structure-function approach combined with quantitative proteomics, we mapped these control points to a conserved YLT motif in the C-terminal tail as well as the catalytic HPD triad in its J domain. Mutation of either motif enhanced DNAJC13 endosomal localization in cells and increased binding to PI(3)P in vitro. Further, these effects required the N-terminal PH-like domain. We show that, similar to other PI(3)P binding domains, the N-terminal PH-like domain binds PI(3)P weakly in isolation and requires oligomerization for efficient PI(3)P binding and endosomal localization. Together, these results demonstrate that interaction between DNAJC13 and PI(3)P serves as a molecular control point for regulating DNAJC13 localization to endosomes.

Petagraph: A large-scale unifying knowledge graph framework for integrating biomolecular and biomedical data

Over the past decade, there has been substantial growth in both the quantity and complexity of available biomedical data. In order to more efficiently harness this extensive data and alleviate challenges associated with integration of multi-omics data, we developed Petagraph, a biomedical knowledge graph that encompasses over 32 million nodes and 118 million relationships. Petagraph leverages more than 180 ontologies and standards in the Unified Biomedical Knowledge Graph (UBKG) to embed millions of quantitative genomics data points. Petagraph provides a cohesive data environment that enables users to efficiently analyze, annotate, and discern relationships within and across complex multi-omics datasets supported by UBKG's annotation scaffold. We demonstrate how queries on Petagraph can generate meaningful results across various research contexts and use cases.

Large-scale CRISPR/Cas9 deletions within the WFDC gene cluster uncover gene functionality and critical roles in mammalian reproduction

Despite 96 million years of evolution separating humans and rodents, 11 closely related reproductive tract-specific genes in humans-SPINT3, WFDC6, EPPIN, WFDC8, WFDC9, WFDC10A, WFDC11, WFDC10B, WFDC13, SPINT4, and WFDC3-and the 13 reproductive tract-specific orthologous genes in mice, form highly conserved syntenic gene clusters indicative of conserved, combined critical functions. Further, despite significant progress toward a nonhormonal male contraceptive targeting the protein encoded by one of these genes, epididymal peptidase inhibitor (EPPIN), and associations found between mutations in EPPIN and an increased risk of male infertility, neither EPPIN nor any closely related whey acidic protein four-disulfide core (WFDC) gene have been explored functionally. To clarify the involvement of WFDC genes in male fertility, we strategically used CRISPR/Cas9 to generate mice lacking 13, 10, 5, or 4 genes within the cluster and demonstrated that males with deletions of 13, 10, or 4 genes (Wfdc6a, Eppin, Wfdc8, and Wfdc6a) were sterile due to an arrest in spermatogenesis, preventing formation beyond round spermatids. In contrast, the five gene knockout (KO) males (lacking Wfdc16, Wfdc9, Wfdc10, Wfdc11, and Wfdc13), despite normal spermatogenesis and sperm counts, were infertile due to defects in sperm motility and increased sperm death. Similarly to our previously reported Spint3 single gene KO, Wfdc3 single KO mice were fertile with no obvious reproductive phenotype. Our KO mouse studies to explore the entire WFDC locus of closely related genes have clarified the functional requirements of WFDC locus genes in different aspects of male fertility. Our research has implications for improving clinical diagnoses of male infertility and identifying additional targets for nonhormonal male contraception.

THE FAM53C/DYRK1A axis regulates the G1/S transition of the cell cycle

A growing number of therapies are being developed to target the cell cycle machinery for the treatment of cancer and other human diseases. Consequently, a greater understanding of the factors regulating cell cycle progression becomes essential to help enhance the response to these new therapies. Here, using data from the Cancer Dependency Map, we identified the poorly-studied factor FAM53C as a new regulator of cell cycle progression. We found that FAM53C is critical for this cell cycle transition and that it acts upstream of the CyclinD-CDK4/6-RB axis in the regulation of the G1/S transition. By mass spectrometry, biochemical, and cellular assays, we identified and validated DYRK1A as a cell cycle kinase that is inhibited by and directly interacts with FAM53C. DYRK1A kinase inhibition rescues the G1 arrest induced by FAM53C knock-down. Consistent with the role for FAM53C identified in cells in culture, FAM53C knockout human cortical organoids display increased cell cycle arrest and growth defects. In addition, Fam53C knockout mice show defects in body growth and behavioral phenotypes. Because DYRK1A dysregulation contributes to developmental disorders such as Down syndrome as well as tumorigenesis, future strategies aiming at regulating FAM53C activity may benefit a broad range of patients.

Differing genetics of saline and cocaine self-administration in the hybrid mouse diversity panel

To identify genes involved in regulating the behavioral and brain transcriptomic response to the potentially addictive drug cocaine, we performed genome-wide association studies (GWASs) for intravenous self-administration of cocaine or saline (as a control) over 10 days using a panel of inbred and recombinant inbred mice. A linear mixed model increased statistical power for these longitudinal data and identified 145 loci for responding when saline only was delivered, compared to 17 for the corresponding cocaine GWAS. Only one locus overlapped. Transcriptome-wide association studies (TWASs) using RNA-Seq data from the medial frontal cortex and nucleus accumbens identified 5031434O11Rik and Zfp60 as significant for saline self-administration. Two other genes, Myh4 and Npc1, were nominated based on proximity to loci for multiple endpoints or a cis locus regulating expression. All four genes have previously been implicated in locomotor activity. Our results indicate distinct genetic bases for saline and cocaine self-administration, and suggest some common genes for saline self-administration and locomotor activity.

Primary cartilage transcriptional signatures reflect cell-type-specific molecular pathways underpinning osteoarthritis

Translational efforts in osteoarthritis are hampered by a gap in our understanding of disease processes at the molecular level. Here, we present evidence of pronounced transcriptional changes in high- and low-disease-grade cartilage tissue, pointing to embryonic processes involved in disease progression. We identify shared transcriptional programs between osteoarthritis cartilage and cell populations in the human embryonic and fetal limb, pointing to increases in pre-hypertrophic chondrocytes' transcriptional programs in low-grade cartilage and increases in osteoblastic signatures in high-grade disease tissue. We find that osteoarthritis genetic risk signals are enriched in six gene co-expression modules and show that these transcriptional signatures reflect cell-type-specific expression along the endochondral ossification developmental trajectory. Using this network approach in combination with causal inference analysis, we present evidence of a causal effect on osteoarthritis risk for variants associated with the expression of ten genes that have not been previously reported as effector genes in genome-wide association studies in osteoarthritis. Our findings point to key molecular pathways as drivers of cartilage degeneration and identify high-value drug targets and repurposing opportunities.

Coevolution and Adaptation of Transition Nuclear Proteins and Protamines in Naturally Ascrotal Mammals Support the Black Queen Hypothesis

Protamines (PRMs) and transition nuclear proteins (TNPs) are two key classes of sperm nuclear basic proteins that regulate chromatin reorganization and condensation in the spermatozoon head, playing crucial roles in mammalian spermatogenesis. In scrotal mammals, such as humans, cryptorchidism, the failure of the testes to descend into the scrotal sac is generally associated with higher rates of defective spermatozoon quality and function. However, ascrotal mammals, such as cetaceans, with naturally undescended testes, produce normal spermatozoa similar to their scrotal counterparts. This study investigates the evolutionary pattern and functional changes in PRMs and TNPs to explore the potential molecular mechanisms underlying spermatogenesis in naturally ascrotal mammals. Although we found a conserved genomic arrangement for PRM and TNP genes across mammals, the coevolutionary loss of intact PRM2 and TNP2 was observed in several species, correlating significantly with diverse testicular positions. Notably, in cetaceans, which lack intact PRM2 and TNP2, we detected enhanced thermostability and DNA binding in PRM1, along with superior DNA repair capability in TNP1. These findings suggest that gene loss of PRM2 and TNP2, combined with functional enhancements in PRM1 and TNP1 proteins, evolved in response to physiological challenges posed by natural cryptorchidism in most ascrotal lineages. This evolutionary strategy enhances chromatin condensation efficiency and promotes DNA repair during spermatogenesis in natural cryptorchid mammals, supporting the Black Queen Hypothesis.

The USP12/46 deubiquitinases protect integrins from ESCRT-mediated lysosomal degradation

The functions of integrins are tightly regulated via multiple mechanisms including trafficking and degradation. Integrins are repeatedly internalized, routed into the endosomal system and either degraded by the lysosome or recycled back to the plasma membrane. The ubiquitin system dictates whether internalized proteins are degraded or recycled. Here, we use a genetic screen and proximity-dependent biotin identification to identify deubiquitinase(s) that control integrin surface levels. We find that a ternary deubiquitinating complex, comprised of USP12 (or the homologous USP46), WDR48 and WDR20, stabilizes β1 integrin (Itgb1) by preventing ESCRT-mediated lysosomal degradation. Mechanistically, the USP12/46-WDR48-WDR20 complex removes ubiquitin from the cytoplasmic tail of internalized Itgb1 in early endosomes, which in turn prevents ESCRT-mediated sorting and Itgb1 degradation.

Investigating the effects of a cryptic splice site in the En2 splice acceptor sequence used in the IKMC knockout-first alleles

Targeted mouse mutants are a common tool used to investigate gene function. The International Knockout Mouse Consortium undertook a large-scale screen of mouse mutants, making use of the knockout-first allele design that contains the En2 splice acceptor sequence coupled to the lacZ reporter gene. Although the knockout-first allele was designed to interfere with splicing and thus disrupt gene function, the En2 sequence has been reported to be transcribed within the host gene mRNA due to a cryptic splice site within the En2 sequence which allows splicing to the next exon of the host gene. In some circumstances, this has the potential to permit translation of a mutant protein. Here, we describe our computational analysis of all the mouse protein-coding genes with established knockout-first embryonic stem cell lines, and our predictions of their transcription outcome should the En2 sequence be included. As part of the large-scale mutagenesis program, mutant mice underwent a broad phenotyping screen, and their phenotypes are available. No wide-scale effects on mouse phenotypes reported were found as a result of the predicted En2 insertion. However, the En2 insertion was found experimentally in the transcripts of 24 of 35 mutant alleles examined, including the five already described, two with evidence of readthrough. Splicing from the cryptic splice site also has the potential to disrupt expression of the lacZ reporter gene. It is recommended that mutant transcripts be checked for this insertion as well as for leaky transcription in studies involving knockout-first alleles.

RELT Is Upregulated in Breast Cancer and Induces Death in Breast Cancer Cells

BACKGROUND

Receptor Expressed in Lymphoid Tissues (RELT) is a TNFRSF member that has two paralogs, RELL1 and RELL2; the three proteins are collectively referred to as RELT family members (RELTfms).

METHODS

We sought to evaluate RELT expression in cancerous cells by using real-time PCR, western blotting, flow cytometry, and immunohistochemistry (IHC). The mechanism of RELT-induced cell death was assessed by western blotting, flow cytometry, luciferase assays, and morphology staining. RELT localization was detected through immunofluorescence and western blotting, and co-immunoprecipitation was used to test whether a mutated RELT interacts with the OXSR1 kinase.

RESULTS

RELT and RELL1 protein expression was significantly elevated in cell lines representing breast and lung cancer, whereas RELL2 protein expression was relatively consistent across different cell lines. The surface expression of RELT was highest in monocytes. IHC staining revealed increased RELT expression in malignant breast cancer biopsies compared to patient-matched benign tissue. RELTfm overexpression induced death in MDA-MB-231 (231) breast cancer cells, accompanied by increased phosphatidylserine externalization and Caspase-3/7 activation. The co-transfection of plasmids predicted to block the phosphorylation of RELT by the OXSR1 kinase did not abrogate RELT-induced apoptosis, indicating that the activation of p38 by RELT through the OXSR1 kinase is not required for RELT-induced cell death. Interestingly, nuclear localization of RELT was detected in 231 and HEK-293 cells.

CONCLUSIONS

These results demonstrate that RELT induces death in breast cancer cells through an apoptotic pathway that does not require OXSR1 phosphorylation and that RELT possesses the ability to translocate to the nucleus, a novel finding that warrants further investigation.

The extracellular vesicle transcriptome provides tissue-specific functional genomic annotation relevant to disease susceptibility in obesity

We characterized circulating extracellular vesicles (EVs) in obese and lean humans, identifying transcriptional cargo differentially expressed in obesity. Since circulating EVs may have broad origin, we compared this obesity EV transcriptome to expression from human visceral adipose tissue derived EVs from freshly collected and cultured biopsies from the same obese individuals. Using a comprehensive set of adipose-specific epigenomic and chromatin conformation assays, we found that the differentially expressed transcripts from the EVs were those regulated in adipose by BMI-associated SNPs from a large-scale GWAS. Using a phenome-wide association study of the regulatory SNPs for the EV-derived transcripts, we identified a substantial enrichment for inflammatory phenotypes, including type 2 diabetes. Collectively, these findings represent the convergence of the GWAS (genetics), epigenomics (transcript regulation), and EV (liquid biopsy) fields, enabling powerful future genomic studies of complex diseases.

Impaired branched-chain amino acid (BCAA) catabolism during adipocyte differentiation decreases glycolytic flux

Dysregulated branched-chain amino acid (BCAA) metabolism has emerged as a key metabolic feature associated with the obese insulin-resistant state, and adipose BCAA catabolism is decreased in this context. BCAA catabolism is upregulated early in adipogenesis, but the impact of suppressing this pathway on the broader metabolic functions of the resultant adipocyte remains unclear. Here, we use CRISPR/Cas9 to decrease BCKDHA in 3T3-L1 and human pre-adipocytes, and ACAD8 in 3T3-L1 pre-adipocytes to induce a deficiency in BCAA catabolism through differentiation. We characterize the transcriptional and metabolic phenotype of 3T1-L1 cells using RNAseq and 13C metabolic flux analysis within a network spanning glycolysis, tricarboxylic acid (TCA) metabolism, BCAA catabolism, and fatty acid synthesis. While lipid droplet accumulation is maintained in Bckdha-deficient adipocytes, they display a more fibroblast-like transcriptional signature. In contrast, Acad8 deficiency minimally impacts gene expression. Decreased glycolytic flux emerges as the most distinct metabolic feature of 3T3-L1 Bckdha-deficient cells, accompanied by a ∼40% decrease in lactate secretion, yet pyruvate oxidation and utilization for de novo lipogenesis is increased to compensate for the loss of BCAA carbon. Deletion of BCKDHA in human adipocyte progenitors also led to a decrease in glucose uptake and lactate secretion; however, these cells did not upregulate pyruvate utilization, and lipid droplet accumulation and expression of adipocyte differentiation markers was decreased in BCKDH knockout cells. Overall our data suggest that human adipocyte differentiation may be more sensitive to the impact of decreased BCKDH activity than 3T3-L1 cells and that both metabolic and regulatory cross-talk exist between BCAA catabolism and glycolysis in adipocytes. Suppression of BCAA catabolism associated with metabolic syndrome may result in a metabolically compromised adipocyte.

Ocular Mucous Membrane Pemphigoid Demonstrates a Distinct Autoantibody Profile from Those of Other Autoimmune Blistering Diseases: A Preliminary Study

Background: Ocular predominant mucous membrane pemphigoid (oMMP) is a severe subtype of autoimmune blistering disease (AIBD), which can result in scarring and vision loss. The diagnosis of oMMP is challenging as patients often have undetectable levels of circulating autoantibodies by conventional assays. Likewise, the principal autoantigen in oMMP has been an area of debate. Methods: In this preliminary experiment, we performed Phage Immunoprecipitation Sequencing (PhIP-seq) on sera from patients with oMMP, as well as non-ocular MMP, bullous pemphigoid, and mucocutaneous-type pemphigus vulgaris. Results: We identified several autoantigens unique to oMMP relative to other AIBDs. We then cross-referenced these antigens against previously published single-nuclei datasets, as well as the International Mouse Phenotyping Consortium Database. Several protein hits identified in our study demonstrated enriched expression on the anterior surface epithelia, including TNKS1BP1, SEC16B, FNBP4, CASZ1, GOLGB1, DOT1L, PRDM 15, LARP4B, and RPL6. Likewise, a previous study of mouse knockout models of murine analogs CASZ1, HIP1, and ELOA2 reported that these mice showed abnormalities in terms of the ocular surface and development in the eyes. Notably, PhIP-seq failed to identify the canonical markers of AIBDs such as BP180, BP230, desmogleins 1 and 3, or integrin β4, indicating that the patient autoantibodies react with conformational epitopes rather than linear epitopes. Conclusions: oMMP patients demonstrate a unique autoantibody repertoire relative to the other AIBDs. Further validation of the identified autoantibodies will shed light on their potentially pathogenic role.

Generation of human induced pluripotent stem cell lines derived from patients of cystic biliary atresia

Biliary atresia (BA), resulting from abnormal development of the liver's internal or external bile ducts, can lead to liver damage and potentially fatal cirrhosis. Type I cystic biliary atresia is a relatively uncommon, but clinically significant variant of BA. It is critical to develop experimental models of BA to examine the etiology and pathogenesis, which remain elusive, and to develop future therapeutics. Here, we have successfully generated a panel of human induced pluripotent stem cells (hiPSCs) from five Japanese patients carrying type I cystic BA. These hiPSC lines exhibited characteristics of self-renewal and pluripotency. These cells held normal karyotypes mostly, but one of them carried hemizygous deletions, the clinical significance of which is unknown yet. Whole genome sequence analysis indicated that some of the mutations or single nucleotide polymorphisms (SNPs) commonly found in these patients are related to hepatobiliary abnormality. Given the limited understanding of the molecular pathogenesis of cystic BA, attributed to unknown factors of genetic and environmental causes, these cellular resources will be instrumental in replicating disease phenotypes and in advancing novel therapies for this disease.

Osteocytes contribute to sex-specific differences in osteoarthritic pain

Osteoarthritic (OA) pain affects 18% of females and 9.6% of males aged over 60 worldwide, with 62% of all OA patients being women. The molecular drivers of sex-based differences in OA are unknown. Bone is intricately coupled with the sensory nervous system and one of the only joint tissues known to show changes that correlate with patient pain in OA. There are fundamental sex differences in pain sensation and bone biology which may be intrinsic to OA disease progression, however these differences are vastly under researched. We have utilised three data sets to investigate the hypothesis that potential mediators responsible for sex dependent pain mechanisms displayed in OA are derived from mechanically stimulated osteocytes. Our published dataset of the in vitro human osteocyte mechanosome was independently compared with published data from, sex-based gene expression differences in human long bone, the sex-based gene expression differences during the skeletal maturation of the mouse osteocyte transcriptome and sex specific OA risk factors and effector genes in a large human GWAS. 80 of the 377 sex-specific genes identified in the mouse osteocyte transcriptome were mechanically regulated in osteocytes with enrichment associated with neural crest migration and axon extension, and DISEASES analysis enrichment for the rheumatoid arthritis pathway. 3861 mechanically regulated osteocytic genes displayed sex-specific differences in human long bone with enrichment for genes associated with the synapse, sensory perception of pain, axon guidance, immune responses, distal peripheral sensory neuropathy, sensory neuropathy, and poor wound healing. 32 of 77 effector genes and 1 of 3 female specific OA risk factor genes identified in the human GWAS were differentially expressed in the osteocyte mechanosome and male and female bone. This analysis lends support to the hypothesis that mechanically regulated genes in osteocytes could influence sex specific differences in osteoarthritic pain and highlights pain pathways with approved drugs that could potentially treat elevated pain susceptibility in females with OA.

Multivariate genomic analysis of 5 million people elucidates the genetic architecture of shared components of the metabolic syndrome

Metabolic syndrome (MetS) is a complex hereditary condition comprising various metabolic traits as risk factors. Although the genetics of individual MetS components have been investigated actively through large-scale genome-wide association studies, the conjoint genetic architecture has not been fully elucidated. Here, we performed the largest multivariate genome-wide association study of MetS in Europe (nobserved = 4,947,860) by leveraging genetic correlation between MetS components. We identified 1,307 genetic loci associated with MetS that were enriched primarily in brain tissues. Using transcriptomic data, we identified 11 genes associated strongly with MetS. Our phenome-wide association and Mendelian randomization analyses highlighted associations of MetS with diverse diseases beyond cardiometabolic diseases. Polygenic risk score analysis demonstrated better discrimination of MetS and predictive power in European and East Asian populations. Altogether, our findings will guide future studies aimed at elucidating the genetic architecture of MetS.

Examining the liver-pancreas crosstalk reveals a role for the molybdenum cofactor in β-cell regeneration

Regeneration of insulin-producing β-cells is an alternative avenue to manage diabetes, and it is crucial to unravel this process in vivo during physiological responses to the lack of β-cells. Here, we aimed to characterize how hepatocytes can contribute to β-cell regeneration, either directly or indirectly via secreted proteins or metabolites, in a zebrafish model of β-cell loss. Using lineage tracing, we show that hepatocytes do not directly convert into β-cells even under extreme β-cell ablation conditions. A transcriptomic analysis of isolated hepatocytes after β-cell ablation displayed altered lipid- and glucose-related processes. Based on the transcriptomics, we performed a genetic screen that uncovers a potential role of the molybdenum cofactor (Moco) biosynthetic pathway in β-cell regeneration and glucose metabolism in zebrafish. Consistently, molybdenum cofactor synthesis 2 (Mocs2) haploinsufficiency in mice indicated dysregulated glucose metabolism and liver function. Together, our study sheds light on the liver-pancreas crosstalk and suggests that the molybdenum cofactor biosynthesis pathway should be further studied in relation to glucose metabolism and diabetes.

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.

SERPINA3 is a marker of cartilage differentiation and is essential for the expression of extracellular matrix genes during early chondrogenesis

Serine proteinase inhibitors (serpins) are a family of structurally similar proteins which regulate many diverse biological processes from blood coagulation to extracellular matrix (ECM) remodelling. Chondrogenesis involves the condensation and differentiation of mesenchymal stem cells (MSCs) into chondrocytes which occurs during early development. Here, and for the first time, we demonstrate that one serpin, SERPINA3 (gene name SERPINA3, protein also known as alpha-1 antichymotrypsin), plays a critical role in chondrogenic differentiation. We observed that SERPINA3 expression was markedly induced at early time points during in vitro chondrogenesis. We examined the expression of SERPINA3 in human cartilage development, identifying significant enrichment of SERPINA3 in developing cartilage compared to total limb, which correlated with well-described markers of cartilage differentiation. When SERPINA3 was silenced using siRNA, cartilage pellets were smaller and contained lower proteoglycan as determined by dimethyl methylene blue assay (DMMB) and safranin-O staining. Consistent with this, RNA sequencing revealed significant downregulation of genes associated with cartilage ECM formation perturbing chondrogenesis. Conversely, SERPINA3 silencing had a negligible effect on the gene expression profile during osteogenesis suggesting the role of SERPINA3 is specific to chondrocyte differentiation. The global effect on cartilage formation led us to investigate the effect of SERPINA3 silencing on the master transcriptional regulator of chondrogenesis, SOX9. Indeed, we observed that SOX9 protein levels were markedly reduced at early time points suggesting a role for SERPINA3 in regulating SOX9 expression and activity. In summary, our data support a non-redundant role for SERPINA3 in enabling chondrogenesis via regulation of SOX9 levels.

Exome variant prioritization in a large cohort of hearing-impaired individuals indicates IKZF2 to be associated with non-syndromic hearing loss and guides future research of unsolved cases

Although more than 140 genes have been associated with non-syndromic hereditary hearing loss (HL), at least half of the cases remain unexplained in medical genetic testing. One reason is that pathogenic variants are located in 'novel' deafness genes. A variant prioritization approach was used to identify novel (candidate) genes for HL. Exome-wide sequencing data were assessed for subjects with presumed hereditary HL that remained unexplained in medical genetic testing by gene-panel analysis. Cases in group AD had presumed autosomal dominantly inherited HL (n = 124), and in group AR, presumed autosomal recessive HL (n = 337). Variants in known and candidate deafness genes were prioritized based on allele frequencies and predicted effects. Selected variants were tested for their co-segregation with HL. Two cases were solved by variants in recently identified deafness genes (ABHD12, TRRAP). Variant prioritization also revealed potentially causative variants in candidate genes associated with recessive and X-linked HL. Importantly, missense variants in IKZF2 were found to co-segregate with dominantly inherited non-syndromic HL in three families. These variants specifically affected Zn2+-coordinating cysteine or histidine residues of the zinc finger motifs 2 and 3 of the encoded protein Helios. This finding indicates a complex genotype-phenotype correlation for IKZF2 defects, as this gene was previously associated with non-syndromic dysfunction of the immune system and ICHAD syndrome, including HL. The designed strategy for variant prioritization revealed that IKZF2 variants can underlie non-syndromic HL. The large number of candidate genes for HL and variants therein stress the importance of inclusion of family members for variant prioritization.

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.

Coordinating BNIP3/NIX-mediated mitophagy in space and time

Mitochondria maintain organellar homeostasis through multiple quality control pathways, including the clearance of defective or unwanted mitochondria by selective autophagy. This removal of mitochondria, mitophagy, is controlled in large part by the outer mitochondrial membrane mitophagy receptors BNIP3 and NIX. While it has long been appreciated that BNIP3 and NIX mediate mitophagy by controlling the recruitment of autophagic machinery to the mitochondrial surface, the requirement for the carefully controlled spatiotemporal regulation of receptor-mediated mitophagy has only recently come to light. Several new factors that regulate the BNIP3/NIX-mediated mitophagy pathway have emerged, and various loss-of-function cell and animal models have revealed the dire consequences of their dysregulation. In this mini-review, we discuss new insights into the mechanisms and roles of the regulation of BNIP3 and NIX and highlight questions that have emerged from the identification of these new regulators.

Cellular and molecular mechanisms that shape the development and evolution of tail vertebral proportion in mice and jerboas

Despite the functional importance of the vertebral skeleton, little is known about how individual vertebrae elongate or achieve disproportionate lengths as in the giraffe neck. Rodent tails are an abundantly diverse and more tractable system to understand mechanisms of vertebral growth and proportion. In many rodents, disproportionately long mid-tail vertebrae form a 'crescendo-decrescendo' of lengths in the tail series. In bipedal jerboas, these vertebrae grow exceptionally long such that the adult tail is 1.5x the length of a mouse tail, relative to body length, with four fewer vertebrae. How do vertebrae with the same regional identity elongate differently from their neighbors to establish and diversify adult proportion? Here, we find that vertebral lengths are largely determined by differences in growth cartilage height and the number of cells progressing through endochondral ossification. Hypertrophic chondrocyte size, a major contributor to differential elongation in mammal limb bones, differs only in the longest jerboa mid-tail vertebrae where they are exceptionally large. To uncover candidate molecular mechanisms of disproportionate vertebral growth, we performed intersectional RNA-Seq of mouse and jerboa tail vertebrae with similar and disproportionate elongation rates. Many regulators of posterior axial identity and endochondral elongation are disproportionately differentially expressed in jerboa vertebrae. Among these, the inhibitory natriuretic peptide receptor C (NPR3) appears in multiple studies of rodent and human skeletal proportion suggesting it refines local growth rates broadly in the skeleton and broadly in mammals. Consistent with this hypothesis, NPR3 loss of function mice have abnormal tail and limb proportions. Therefore, in addition to genetic components of the complex process of vertebral evolution, these studies reveal fundamental mechanisms of skeletal growth and proportion.

Systematic identification of therapeutic targets for coronary artery calcification: an integrated transcriptomic and proteomic Mendelian randomization

Background

Coronary artery calcification (CAC) is associated with an increased risk of mortality and cardiovascular events. However, none therapeutic drugs have been proven effective for CAC treatment. The objective of this study was to identify potential therapeutic targets for CAC through the utilization of Mendelian randomization (MR) and colocalization analysis.

Methods

The expression quantitative trait loci (eQTLs) of 16,943 genes from the eQTLGen consortium and protein quantitative trait loci (pQTLs) of 4,412 proteins from a plasma proteome were utilized as genetic instruments. Genetic associations with CAC were derived from a GWAS meta-analysis of 26,909 individuals. The MR and colocalization analysis were utilized to identify potential target genes.

Results

A total of 671 genes were found to be significantly associated with the risk of CAC based on transcriptomic MR analysis at a false discovery rate <0.05, while proteomic MR analysis identified 15 genes with significant associations with CAC at the same threshold. With robust evidence from colocalization analysis, we observed positive associations between CWF19L2, JARID2, and MANBA and the risk of CAC, while KLB exhibited an inverse association. In summary, our study identified 23 potential therapeutic targets for CAC. Further downstream analysis revealed IGFBP3, ABCC6, ULK3, DOT1L, KLB and AMH as promising candidates for repurposing in the treatment of CAC.

Conclusion

The integrated MR analysis of transcriptomic and proteomic data identified multiple potential drug targets for the treatment of CAC. ULK3, DOT1L, and AMH were recognized as novel targets for drug repurposing for CAC and deserve further investigation.

A disease-associated PPP2R3C-MAP3K1 phospho-regulatory module controls centrosome function

Centrosomes have critical roles in microtubule organization, ciliogenesis, and cell signaling.1,2,3,4,5,6,7,8 Centrosomal alterations also contribute to diseases, including microcephaly, cancer, and ciliopathies.9,10,11,12,13 To date, over 150 centrosomal proteins have been identified, including several kinases and phosphatases that control centrosome biogenesis, function, and maintenance.2,3,4,5,14,15,16,17,18,19,20,21 However, the regulatory mechanisms that govern centrosome function are not fully defined, and thus how defects in centrosomal regulation contribute to disease is incompletely understood. Using a systems genetics approach, we find here that PPP2R3C, a poorly characterized PP2A phosphatase subunit, is a distal centriole protein and functional partner of centriolar proteins CEP350 and FOP. We further show that a key function of PPP2R3C is to counteract the kinase activity of MAP3K1. In support of this model, MAP3K1 knockout suppresses growth defects caused by PPP2R3C inactivation, and MAP3K1 and PPP2R3C have opposing effects on basal and microtubule stress-induced JNK signaling. Illustrating the importance of balanced MAP3K1 and PPP2R3C activities, acute overexpression of MAP3K1 severely inhibits centrosome function and triggers rapid centriole disintegration. Additionally, inactivating PPP2R3C mutations and activating MAP3K1 mutations both cause congenital syndromes characterized by gonadal dysgenesis.22,23,24,25,26,27,28 As a syndromic PPP2R3C variant is defective in centriolar localization and binding to centriolar protein FOP, we propose that imbalanced activity of this centrosomal kinase-phosphatase pair is the shared cause of these disorders. Thus, our findings reveal a new centrosomal phospho-regulatory module, shed light on disorders of gonadal development, and illustrate the power of systems genetics to identify previously unrecognized gene functions.

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

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

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

Introduction

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

Methods

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

Results

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

Conclusion

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Working with Miraculous Mice: Mus musculus as a Model Organism

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

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

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

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

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

Overabundant endocannabinoids in neurons are detrimental to cognitive function

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

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

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

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

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

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

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