The mouse Gene Expression Database (GXD): 2019 update

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.

Genome-wide analysis of hepatic DNA methylation reveals impact of epigenetic aging on xenobiotic metabolism and transport genes in an aged mouse model

Hepatic xenobiotic metabolism and transport decline with age, while intact xenobiotic metabolism is associated with longevity. However, few studies have examined the genome-wide impact of epigenetic aging on these processes. We used reduced representation bisulfite sequencing (RRBS) to map DNA methylation changes in liver DNA from mice ages 4 and 24 months. We identified several thousand age-associated differentially methylated sites (a-DMS), many of which overlapped genes encoding Phase I and Phase II drug metabolizing enzymes, in addition to ABC and SLC classes of transporters. Notable genes harboring a-DMS were Cyp1a2, Cyp2d9, and Abcc2 that encode orthologs of the human drug metabolizing enzymes CYP1A2 and CYP2D6, and the multidrug resistance protein 2 (MRP2) transporter. Cyp2d9 hypermethylation with age was significantly associated with reduced gene expression, while Abcc2 expression was unchanged with age. Cyp1a2 lost methylation with age while, counterintuitively, its expression also reduced with age. We hypothesized that age-related dysregulation of the hepatic transcriptional machinery caused down-regulation of genes despite age-related hypomethylation. Bioinformatic analysis of hypomethylated a-DMS in our sample found them to be highly enriched for hepatic nuclear factor 4 alpha (HNF4α) binding sites. HNF4α promotes Cyp1a2 expression and is downregulated with age, which could explain the reduction in Cyp1a2 expression. Overall, our study supports the broad impact of epigenetic aging on xenobiotic metabolism and transport. Future work should evaluate the interplay between hepatic nuclear receptor function and epigenetic aging. These results may have implications for studies of longevity and healthy aging.

Histone H1 facilitates restoration of H3K27me3 during DNA replication by chromatin compaction

During cell renewal, epigenetic information needs to be precisely restored to maintain cell identity and genome integrity following DNA replication. The histone mark H3K27me3 is essential for the formation of facultative heterochromatin and the repression of developmental genes in embryonic stem cells. However, how the restoration of H3K27me3 is precisely achieved following DNA replication is still poorly understood. Here we employ ChOR-seq (Chromatin Occupancy after Replication) to monitor the dynamic re-establishment of H3K27me3 on nascent DNA during DNA replication. We find that the restoration rate of H3K27me3 is highly correlated with dense chromatin states. In addition, we reveal that the linker histone H1 facilitates the rapid post-replication restoration of H3K27me3 on repressed genes and the restoration rate of H3K27me3 on nascent DNA is greatly compromised after partial depletion of H1. Finally, our in vitro biochemical experiments demonstrate that H1 facilitates the propagation of H3K27me3 by PRC2 through compacting chromatin. Collectively, our results indicate that H1-mediated chromatin compaction facilitates the propagation and restoration of H3K27me3 after DNA replication.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Single cell RNA analysis of the left-right organizer transcriptome reveals potential novel heterotaxy genes

The establishment of left-right patterning in mice occurs at a transient structure called the embryonic node or left-right organizer (LRO). Previous analysis of the LRO has proven challenging due to the small cell number and transient nature of this structure. Here, we seek to overcome these difficulties to define the transcriptome of the LRO. Specifically, we used single cell RNA sequencing of 0-1 somite embryos to identify LRO enriched genes which were compared to bulk RNA sequencing of LRO cells isolated by fluorescent activated cell sorting. Gene ontology analysis indicated an enrichment of genes associated with cilia and laterality terms. Furthermore, comparison to previously identified LRO genes identified 127 novel LRO genes, including Ttll3, Syne1 and Sparcl1, for which the expression patterns were validated using whole mount in situ hybridization. This list of novel LRO genes will be a useful resource for further studies on LRO morphogenesis, the establishment of laterality and the genetic causes of heterotaxy.

IL11 activates the placental inflammasome to drive preeclampsia

Introduction

Preeclampsia is a life-threatening disorder of pregnancy unique to humans. Interleukin (IL)11 is elevated in serum from pregnancies that subsequently develop early-onset preeclampsia and pharmacological elevation of IL11 in pregnant mice causes the development of early-onset preeclampsia-like features (hypertension, proteinuria, and fetal growth restriction). However, the mechanism by which IL11 drives preeclampsia is unknown.

Method

Pregnant mice were administered PEGylated (PEG)IL11 or control (PEG) from embryonic day (E)10-16 and the effect on inflammasome activation, systolic blood pressure (during gestation and at 50/90 days post-natal), placental development, and fetal/post-natal pup growth measured. RNAseq analysis was performed on E13 placenta. Human 1^st trimester placental villi were treated with IL11 and the effect on inflammasome activation and pyroptosis identified by immunohistochemistry and ELISA.

Result

PEGIL11 activated the placental inflammasome causing inflammation, fibrosis, and acute and chronic hypertension in wild-type mice. Global and placental-specific loss of the inflammasome adaptor protein Asc and global loss of the Nlrp3 sensor protein prevented PEGIL11-induced fibrosis and hypertension in mice but did not prevent PEGIL11-induced fetal growth restriction or stillbirths. RNA-sequencing and histology identified that PEGIL11 inhibited trophoblast differentiation towards spongiotrophoblast and syncytiotrophoblast lineages in mice and extravillous trophoblast lineages in human placental villi.

Discussion

Inhibition of ASC/NLRP3 inflammasome activity could prevent IL11-induced inflammation and fibrosis in various disease states including preeclampsia.

The HAPSTR2 retrogene buffers stress signaling and resilience in mammals

We recently identified HAPSTR1 (C16orf72) as a key component in a novel pathway which regulates the cellular response to molecular stressors, such as DNA damage, nutrient scarcity, and protein misfolding. Here, we identify a functional paralog to HAPSTR1: HAPSTR2. HAPSTR2 formed early in mammalian evolution, via genomic integration of a reverse transcribed HAPSTR1 transcript, and has since been preserved under purifying selection. HAPSTR2, expressed primarily in neural and germline tissues and a subset of cancers, retains established biochemical features of HAPSTR1 to achieve two functions. In normal physiology, HAPSTR2 directly interacts with HAPSTR1, markedly augmenting HAPSTR1 protein stability in a manner independent from HAPSTR1's canonical E3 ligase, HUWE1. Alternatively, in the context of HAPSTR1 loss, HAPSTR2 expression is sufficient to buffer stress signaling and resilience. Thus, we discover a mammalian retrogene which safeguards fitness.

Roles for the RNA-Binding Protein Caper in Reproductive Output in Drosophila melanogaster

RNA binding proteins (RBPs) play a fundamental role in the post-transcriptional regulation of gene expression within the germline and nervous system. This is underscored by the prevalence of mutations within RBP-encoding genes being implicated in infertility and neurological disease. We previously described roles for the highly conserved RBP Caper in neurite morphogenesis in the Drosophila larval peripheral system and in locomotor behavior. However, caper function has not been investigated outside the nervous system, although it is widely expressed in many different tissue types during embryogenesis. Here, we describe novel roles for Caper in fertility and mating behavior. We find that Caper is expressed in ovarian follicles throughout oogenesis but is dispensable for proper patterning of the egg chamber. Additionally, reduced caper function, through either a genetic lesion or RNA interference-mediated knockdown of caper in the female germline, results in females laying significantly fewer eggs than their control counterparts. Moreover, this phenotype is exacerbated with age. caper dysfunction also results in partial embryonic and larval lethality. Given that caper is highly conserved across metazoa, these findings may also be relevant to vertebrates.

Glial Bmal1 role in mammalian retina daily changes

Visual information processing in the retina requires the rhythmic expression of clock genes. The intrinsic retinal circadian clock is independent of the master clock located in the hypothalamic suprachiasmatic nucleus and emerges from retinal cells, including glia. Less clear is how glial oscillators influence the daily regulation of visual information processing in the mouse retina. Here, we demonstrate that the adult conditional deletion of the gene Bmal1 in GLAST-positive glial cells alters retinal physiology. Specifically, such deletion was sufficient to lower the amplitude of the electroretinogram b-wave recorded under light-adapted conditions. Furthermore, recordings from > 20,000 retinal ganglion cells (RGCs), the retina output, showed a non-uniform effect on RGCs activity in response to light across different cell types and over a 24-h period. Overall, our results suggest a new role of a glial circadian gene in adjusting mammalian retinal output throughout the night-day cycle.

Genetic ablation of the preptin-coding portion of Igf2 impairs pancreatic function in female mice

Preptin is a 34-amino acid peptide derived from the E-peptide of pro-insulin-like growth factor 2 and is co-secreted with insulin from β-cells. Little is understood about the effects of endogenous preptin on whole body glucose metabolism. We developed a novel mouse model in which the preptin portion of Igf2 was genetically ablated in all tissues, hereafter referred to as preptin knockout (KO), and tested the hypothesis that the removal of preptin will lead to a decreased insulin response to a metabolic challenge. Preptin KO and wild-type (WT) mice underwent weekly fasting blood glucose measurements, intraperitoneal insulin tolerance tests (ITT) at 9, 29, and 44 wk of age, and an oral glucose tolerance test (GTT) at 45 wk of age. Preptin KO mice of both sexes had similar Igf2 exon 2-3 mRNA expression in the liver and kidney compared with WT mice, but Igf2 exon 3-4 (preptin) expression was not detectable. Western blot analysis of neonatal serum indicated that processing of pro-IGF2 translated from the KO allele may be altered. Preptin KO mice had similar body weight, body composition, β-cell area, and fasted glucose concentrations compared with WT mice in both sexes up to 47 wk of age. Female KO mice had a diminished ability to mount an insulin response following glucose stimulation in vivo. This effect was absent in male KO mice. Although preptin is not essential for glucose homeostasis, when combined with previous in vitro and ex vivo findings, these data show that preptin positively impacts β-cell function.NEW & NOTEWORTHY This is the first study to describe a model in which the preptin-coding portion of the Igf2 gene has been genetically ablated in mice. The mice do not show reduced size at birth associated with Igf2 knockout suggesting that IGF2 functionality is maintained, yet we demonstrate a change in the processing of mature Igf2. Female knockout mice have diminished glucose-stimulated insulin secretion, whereas the insulin response in males is not different to wild type.

The contribution of sex chromosome conflict to disrupted spermatogenesis in hybrid house mice

Incompatibilities on the sex chromosomes are important in the evolution of hybrid male sterility, but the evolutionary forces underlying this phenomenon are unclear. House mice (Mus musculus) lineages have provided powerful models for understanding the genetic basis of hybrid male sterility. X chromosome-autosome interactions cause strong incompatibilities in M. musculus F1 hybrids, but variation in sterility phenotypes suggests a more complex genetic basis. In addition, XY chromosome conflict has resulted in rapid expansions of ampliconic genes with dosage-dependent expression that is essential to spermatogenesis. Here, we evaluated the contribution of XY lineage mismatch to male fertility and stage-specific gene expression in hybrid mice. We performed backcrosses between two house mouse subspecies to generate reciprocal Y-introgression strains and used these strains to test the effects of XY mismatch in hybrids. Our transcriptome analyses of sorted spermatid cells revealed widespread overexpression of the X chromosome in sterile F1 hybrids independent of Y chromosome subspecies origin. Thus, postmeiotic overexpression of the X chromosome in sterile F1 mouse hybrids is likely a downstream consequence of disrupted meiotic X-inactivation rather than XY gene copy number imbalance. Y chromosome introgression did result in subfertility phenotypes and disrupted expression of several autosomal genes in mice with an otherwise nonhybrid genomic background, suggesting that Y-linked incompatibilities contribute to reproductive barriers, but likely not as a direct consequence of XY conflict. Collectively, these findings suggest that rapid sex chromosome gene family evolution driven by genomic conflict has not resulted in strong male reproductive barriers between these subspecies of house mice.

Homologous gene regulatory networks control development of apical organs and brains in Bilateria

Apical organs are relatively simple larval nervous systems. The extent to which apical organs are evolutionarily related to the more complex nervous systems of other animals remains unclear. To identify common developmental mechanisms, we analyzed the gene regulatory network (GRN) controlling the development of the apical organ in sea urchins. We characterized the developmental expression of 30 transcription factors and identified key regulatory functions for FoxQ2, Hbn, Delta/Notch signaling, and SoxC in the patterning of the apical organ and the specification of neurons. Almost the entire set of apical transcription factors is expressed in the nervous system of worms, flies, zebrafish, frogs, and mice. Furthermore, a regulatory module controlling the axial patterning of the vertebrate brain is expressed in the ectoderm of sea urchin embryos. We conclude that GRNs controlling the formation of bilaterian nervous systems share a common origin and that the apical GRN likely resembles an ancestral regulatory program.

Hepatic FGF21: Its Emerging Role in Inter-Organ Crosstalk and Cancers

Fibroblast growth factor (FGF) 21 is one of the FGF members with special endocrine properties. In the last twenty years, it has attracted intense research and development for its physiological functions that respond to dietary manipulation, pharmacological benefits of improving the macronutrient metabolism, and clinical values as a biomarker of various human diseases. Generally, FGF21 can be produced by major metabolic organs, but only the subgroup from the liver shows canonical endocrine properties, which emphasizes the special value of delineating the unique secretory and functional characteristics of hepatic FGF21. There has been a growth in literature to address the extra-hepatic activities of FGF21, and many striking findings have therefore been published. Yet, they are fragmented and scattered, and controversies are raised from divergent findings. For this reason, there is a need for a systematic and critical evaluation of current research in this aspect. In this review, we focus on the current knowledge about the molecular biology of endocrine FGF21, especially present details on the regulation of circulating levels of FGF21. We also emphasize its emerging roles in inter-organ crosstalk and cancer development.

Versatile Enzymology and Heterogeneous Phenotypes in Cobalamin Complementation Type C Disease

Nutritional deficiency and genetic errors that impair the transport, absorption, and utilization of vitamin B₁₂ (B₁₂) lead to hematological and neurological manifestations. The cblC disease (cobalamin complementation type C) is an autosomal recessive disorder caused by mutations and epi-mutations in the MMACHC gene and the most common inborn error of B₁₂ metabolism. Pathogenic mutations in MMACHC disrupt enzymatic processing of B₁₂, an indispensable step before micronutrient utilization by the two B₁₂-dependent enzymes methionine synthase (MS) and methylmalonyl-CoA mutase (MUT). As a result, patients with cblC disease exhibit plasma elevation of homocysteine (Hcy, substrate of MS) and methylmalonic acid (MMA, degradation product of methylmalonyl-CoA, substrate of MUT). The cblC disorder manifests early in childhood or in late adulthood with heterogeneous multi-organ involvement. This review covers current knowledge on the cblC disease, structure–function relationships of the MMACHC protein, the genotypic and phenotypic spectra in humans, experimental disease models, and promising therapies.

Cytokine Receptor-Like Factor 3 (CRLF3) Contributes to Early Zebrafish Hematopoiesis

Cytokine receptor-like factor 3 (CRLF3) is an ancient protein conserved across metazoans that contains an archetypal cytokine receptor homology domain (CHD). This domain is found in cytokine receptors present in bilateria, including higher vertebrates, that play key roles in a variety of developmental and homeostatic processes, particularly relating to blood and immune cells. However, understanding of CRLF3 itself remains very limited. This study aimed to investigate this evolutionarily significant protein by studying its embryonic expression and function in early development, particularly of blood and immune cells, using zebrafish as a model. Expression of crlf3 was identified in mesoderm-derived tissues in early zebrafish embryos, including the somitic mesoderm and both anterior and posterior lateral plate mesoderm. Later expression was observed in the thymus, brain, retina and exocrine pancreas. Zebrafish crlf3 mutants generated by genome editing technology exhibited a significant reduction in primitive hematopoiesis and early definitive hematopoiesis, with decreased early progenitors impacting on multiple lineages. No other obvious phenotypes were observed in the crlf3 mutants.

Somatostatin receptor ligands suppressed proliferation and lipogenesis in 3T3-L1 preadipocytes

Somatostatin and its analogues, known as somatostatin receptor ligands (SRLs), have been reported to attenuate weight gain in some clinical settings. However, their direct effects on preadipocytes are barely investigated. Therefore, this study aimed to evaluate the influence of SRLs on preadipocytes and to further explore the potential mechanisms. Cell Counting Kit-8 assay, Oil Red O staining, triglyceride contents measurements, qPCR, and western blot were used to investigate the effects of SRLs on preadipocytes. We found that three SRLs (octreotide, TT232, and pasireotide) inhibited cell viability after 8-48 h but not 4 h. Further western blot results showed that they significantly suppressed activation of PI3K/Akt pathway. Besides, lipid accumulation was also significantly inhibited by these SRLs. Moreover, mRNA levels of some critical adipogenic markers, including Pparg, Cebpa, Fasn, Fabp4, Acaca, and Lpl, were downregulated by the treatments of all these SRLs. Consistently, the protein expression of PPARγ, C/EBPα, and FAS were also suppressed by SRLs. SRLs inhibit the proliferation and lipogenesis in preadipocytes. Their inhibitory effects on cell proliferation may be mediated by the downregulated PI3K/Akt pathway, and the suppressive actions on lipogenesis may be related to the decreased PPARγ and C/EBPα expression.

SLC38A10 Knockout Mice Display a Decreased Body Weight and an Increased Risk-Taking Behavior in the Open Field Test

The solute carrier 38 family (SLC38) is a family of 11 members. The most common substrate among these are alanine and glutamine, and members are present in a wide range of tissues with important functions for several biological processes, such as liver and brain function. Some of these transporters are better characterized than others and, in this paper, a behavioral characterization of SLC38A10^−/− mice was carried out. A battery of tests for general activity, emotionality, motor function, and spatial memory was used. Among these tests, the elevated plus maze, Y-maze, marble burying and challenging beam walk have not been tested on the SLC38A10^−/− mice previously, while the open field and the rotarod tests have been performed by the International Mouse Phenotyping Consortium (IMPC). Unlike the results from IMPC, the results from this study showed that SLC38A10^−/− mice spend less time in the wall zone in the open field test than WT mice, implying that SLC38A10-deficient mice have an increased explorative behavior, which suggests an important function of SLC38A10 in brain. The present study also confirmed IMPC's data regarding rotarod performance and weight, showing that SLC38A10^−/− mice do not have an affected motor coordination impairment and have a lower body weight than both SLC38A10^+/− and SLC38A10^+/+ mice. These results imply that a complete deficiency of the SLC38A10 protein might affect body weight homeostasis, but the underlying mechanisms needs to be studied further.

Establishment of H3K9-methylated heterochromatin and its functions in tissue differentiation and maintenance

Heterochromatin is characterized by dimethylated or trimethylated histone H3 Lys9 (H3K9me2 or H3K9me3, respectively) and is found at transposable elements, satellite repeats and genes, where it ensures their transcriptional silencing. The histone methyltransferases (HMTs) that methylate H3K9 — in mammals Suppressor of variegation 3–9 homologue 1 (SUV39H1), SUV39H2, SET domain bifurcated 1 (SETDB1), SETDB2, G9A and G9A-like protein (GLP) — and the ‘readers’ of H3K9me2 or H3K9me3 are highly conserved and show considerable redundancy. Despite their redundancy, genetic ablation or mistargeting of an individual H3K9 methyltransferase can correlate with impaired cell differentiation, loss of tissue identity, premature aging and/or cancer. In this Review, we discuss recent advances in understanding the roles of the known H3K9-specific HMTs in ensuring transcriptional homeostasis during tissue differentiation in mammals. We examine the effects of H3K9-methylation-dependent gene repression in haematopoiesis, muscle differentiation and neurogenesis in mammals, and compare them with mechanistic insights obtained from the study of model organisms, notably Caenorhabditis elegans and Drosophila melanogaster. In all these organisms, H3K9-specific HMTs have both unique and redundant roles that ensure the maintenance of tissue integrity by restricting the binding of transcription factors to lineage-specific promoters and enhancer elements.

Histone H3 Lys9 (H3K9)-methylated heterochromatin ensures transcriptional silencing of repetitive elements and genes, and its deregulation leads to impaired cell and tissue identity, premature aging and cancer. Recent studies in mammals clarified the roles H3K9-specific histone methyltransferases in ensuring transcriptional homeostasis during tissue differentiation.

Incomplete lineage sorting and phenotypic evolution in marsupials

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity

The spatial and temporal domain of a gene's expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate's stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates' stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

Exome sequencing of families from Ghana reveals known and candidate hearing impairment genes

We investigated hearing impairment (HI) in 51 families from Ghana with at least two affected members that were negative for GJB2 pathogenic variants. DNA samples from 184 family members underwent whole-exome sequencing (WES). Variants were found in 14 known non-syndromic HI (NSHI) genes [26/51 (51.0%) families], five genes that can underlie either syndromic HI or NSHI [13/51 (25.5%)], and one syndromic HI gene [1/51 (2.0%)]. Variants in CDH23 and MYO15A contributed the most to HI [31.4% (16/51 families)]. For DSPP, an autosomal recessive mode of inheritance was detected. Post-lingual expression was observed for a family segregating a MARVELD2 variant. To our knowledge, seven novel candidate HI genes were identified (13.7%), with six associated with NSHI (INPP4B, CCDC141, MYO19, DNAH11, POTEI, and SOX9); and one (PAX8) with Waardenburg syndrome. MYO19 and DNAH11 were replicated in unrelated Ghanaian probands. Six of the novel genes were expressed in mouse inner ear. It is known that Pax8-/- mice do not respond to sound, and depletion of Sox9 resulted in defective vestibular structures and abnormal utricle development. Most variants (48/60; 80.0%) have not previously been associated with HI. Identifying seven candidate genes in this study emphasizes the potential of novel HI genes discovery in Africa.

Postnatal expression of IGF2 is the norm in amniote vertebrates

The insulin and insulin-like signalling (IIS) network plays an important role in mediating several life-history traits, including growth, reproduction and senescence. Although insulin-like growth factors (IGFs) 1 and 2 are both key hormones in the vertebrate IIS network, research on IGF2 in juveniles and adults has been largely neglected because early biomedical research on rodents found negligible IGF2 postnatal expression. Here, we challenge this assumption and ask to what degree IGF2 is expressed during postnatal life across amniotes by quantifying the relative gene expression of IGF1 and IGF2 using publicly available RNAseq data for 82 amniote species and quantitative polymerase chain reaction on liver cDNA at embryonic, juvenile and adult stages for two lizard, bird and mouse species. We found that (i) IGF2 is expressed postnatally across amniote species and life stages-often at a higher relative expression than IGF1, contradicting rodent models; (ii) the lack of rodent postnatal IGF2 expression is due to phylogenetic placement, not inbreeding or artificial selection; and (iii) adult IGF2 expression is sex-biased in some species. Our results demonstrate that IGF2 expression is typical for amniotes throughout life, suggesting that a comprehensive understanding of the mechanisms mediating variation in life-history traits will require studies that measure both IGFs.

Whole Genome Sequencing Unravels New Genetic Determinants of Early-Onset Familial Osteoporosis and Low BMD in Malta

Background: Osteoporosis is a skeletal disease with a strong genetic background. The study aimed to identify the genetic determinants of early-onset familial osteoporosis and low bone mineral density (BMD) in a two-generation Maltese family. Methods: Fifteen relatives aged between 28–74 years were recruited. Whole genome sequencing was conducted on 12 relatives and shortlisted variants were genotyped in the Malta Osteoporotic Fracture Study (MOFS) for replication. Results: Sequential variant filtering following a dominant inheritance pattern identified rare missense variants within SELP, TGF-β2 and ADAMTS20, all of which were predicted to be likely pathogenic and participate in osteoimmunology. TGF-β2 c.1136C>T was identified in five individuals from the MOFS in heterozygosity, four of whom had osteopenia/osteoporosis at the lumbar spine and hip, and/or had sustained a low-trauma fracture. Heterozygosity for the ADAMTS20 c.4090A>T was accompanied by lower total hip BMD (p = 0.018) and lower total serum calcium levels in MOFS (p < 0.01), recapitulating the findings from the family. Women carrying at least one copy of the alternative allele (TC/CC) for SELP c.2177T>C exhibited a tendency for lower lumbar spine BMD and/or wrist fracture history relative to women with TT genotype. Conclusions: Our findings suggest that the identified variants, alone or in combination, could be causal factors of familial osteoporosis and low BMD, requiring replication in larger collections.

Sex Differences in Embryonic Gonad Transcriptomes and Benzo[a]pyrene Metabolite Levels After Transplacental Exposure

Polycyclic aromatic hydrocarbons like benzo[a]pyrene (BaP) are generated during incomplete combustion of organic materials. Prior research has demonstrated that BaP is a prenatal ovarian toxicant and carcinogen. However, the metabolic pathways active in the embryo and its developing gonads and the mechanisms by which prenatal exposure to BaP predisposes to ovarian tumors later in life remain to be fully elucidated. To address these data gaps, we orally dosed pregnant female mice with BaP from embryonic day (E) 6.5 to E11.5 (0, 0.2, or 2 mg/kg/day) for metabolite measurement or E9.5 to E11.5 (0 or 3.33 mg/kg/day) for embryonic gonad RNA sequencing. Embryos were harvested at E13.5 for both experiments. The sum of BaP metabolite concentrations increased significantly with dose in the embryos and placentas, and concentrations were significantly higher in female than male embryos and in embryos than placentas. RNA sequencing revealed that enzymes involved in metabolic activation of BaP are expressed at moderate to high levels in embryonic gonads and that greater transcriptomic changes occurred in the ovaries in response to BaP than in the testes. We identified 490 differentially expressed genes (DEGs) with false discovery rate P-values < 0.05 when comparing BaP-exposed to control ovaries but no statistically significant DEGs between BaP-exposed and control testes. Genes related to monocyte/macrophage recruitment and activity, prolactin family genes, and several keratin genes were among the most upregulated genes in the BaP-exposed ovaries. Results show that developing ovaries are more sensitive than testes to prenatal BaP exposure, which may be related to higher concentrations of BaP metabolites in female embryos.

A cis-regulatory-directed pipeline for the identification of genes involved in cardiac development and disease

BACKGROUND

Congenital heart diseases are the major cause of death in newborns, but the genetic etiology of this developmental disorder is not fully known. The conventional approach to identify the disease-causing genes focuses on screening genes that display heart-specific expression during development. However, this approach would have discounted genes that are expressed widely in other tissues but may play critical roles in heart development.

RESULTS

We report an efficient pipeline of genome-wide gene discovery based on the identification of a cardiac-specific cis-regulatory element signature that points to candidate genes involved in heart development and congenital heart disease. With this pipeline, we retrieve 76% of the known cardiac developmental genes and predict 35 novel genes that previously had no known connectivity to heart development. Functional validation of these novel cardiac genes by RNAi-mediated knockdown of the conserved orthologs in Drosophila cardiac tissue reveals that disrupting the activity of 71% of these genes leads to adult mortality. Among these genes, RpL14, RpS24, and Rpn8 are associated with heart phenotypes.

CONCLUSIONS

Our pipeline has enabled the discovery of novel genes with roles in heart development. This workflow, which relies on screening for non-coding cis-regulatory signatures, is amenable for identifying developmental and disease genes for an organ without constraining to genes that are expressed exclusively in the organ of interest.

Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

Expanded phenotype of primary ciliary dyskinesia related to DRC1 pathogenic variant with dysmorphisms and vascular anomalies

We present a case of a female diagnosed with primary ciliary dyskinesia (PCD) type 21 with non-previously reported extrapulmonary symptoms, including facial features and congenital vascular anomalies. Whole genome sequencing in our patient revealed a homozygous pathogenic variant in the DRC1 gene and no other notable structural nor punctual variants. This case demonstrates a unique clinical manifestation of PCD, which is possibly associated with the presence of a homozygous pathogenic DRC1 variant. Therefore, we suggest that analysis of DRC1 be considered with PCD type 21 when such features are present.

The effects of Tbx15 and Pax1 on facial and other physical morphology in mice

DNA variants in or close to the human TBX15 and PAX1 genes have been repeatedly associated with facial morphology in independent genome-wide association studies, while their functional roles in determining facial morphology remain to be understood. We generated Tbx15 knockout (Tbx15 -/-) and Pax1 knockout (Pax1 -/-) mice by applying the one-step CRISPR/Cas9 method. A total of 75 adult mice were used for subsequent phenotype analysis, including 38 Tbx15 mice (10 homozygous Tbx15 -/-, 18 heterozygous Tbx15 +/-, 10 wild-type Tbx15 +/+ WT littermates) and 37 Pax1 mice (12 homozygous Pax1 -/-, 15 heterozygous Pax1 +/-, 10 Pax1 +/+ WT littermates). Facial and other physical morphological phenotypes were obtained from three-dimensional (3D) images acquired with the HandySCAN BLACK scanner. Compared to WT littermates, the Tbx15 -/- mutant mice had significantly shorter faces (p = 1.08E-8, R2 = 0.61) and their ears were in a significantly lower position (p = 3.54E-8, R2 = 0.62) manifesting a "droopy ear" characteristic. Besides these face alternations, Tbx15 -/- mutant mice displayed significantly lower weight as well as shorter body and limb length. Pax1 -/- mutant mice showed significantly longer noses (p = 1.14E-5, R2 = 0.46) relative to WT littermates, but otherwise displayed less obvious morphological alterations than Tbx15 -/- mutant mice did. We provide the first direct functional evidence that two well-known and replicated human face genes, Tbx15 and Pax1, impact facial and other body morphology in mice. The general agreement between our findings in knock-out mice with those from previous GWASs suggests that the functional evidence we established here in mice may also be relevant in humans.

Mouse Genome Informatics (MGI): latest news from MGD and GXD

The Mouse Genome Informatics (MGI) database system combines multiple expertly curated community data resources into a shared knowledge management ecosystem united by common metadata annotation standards. MGI's mission is to facilitate the use of the mouse as an experimental model for understanding the genetic and genomic basis of human health and disease. MGI is the authoritative source for mouse gene, allele, and strain nomenclature and is the primary source of mouse phenotype annotations, functional annotations, developmental gene expression information, and annotations of mouse models with human diseases. MGI maintains mouse anatomy and phenotype ontologies and contributes to the development of the Gene Ontology and Disease Ontology and uses these ontologies as standard terminologies for annotation. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are MGI's two major knowledgebases. Here, we highlight some of the recent changes and enhancements to MGD and GXD that have been implemented in response to changing needs of the biomedical research community and to improve the efficiency of expert curation. MGI can be accessed freely at http://www.informatics.jax.org .

Beyond Neurons: Long Distance Communication in Development and Cancer

Cellular communication is important in all aspects of tissue and organism functioning, from the level of single cells, two discreet populations, and distant tissues of the body. Long distance communication networks integrate individual cells into tissues to maintain a complex organism during development, but when communication between cells goes awry, disease states such as cancer emerge. Herein we discuss the growing body of evidence suggesting that communication methods known to be employed by neurons, also exist in other cell types. We identify three major areas of long-distance communication: bioelectric signaling, tunneling nanotubes (TNTs), and macrophage modulation of networks, and draw comparisons about how these systems operate in the context of development and cancer. Bioelectric signaling occurs between cells through exchange of ions and tissue-level electric fields, leading to changes in biochemical gradients and molecular signaling pathways to control normal development and tumor growth and invasion in cancer. TNTs transport key morphogens and other cargo long distances, mediating electrical coupling, tissue patterning, and malignancy of cancer cells. Lastly macrophages maintain long distance signaling networks through trafficking of vesicles during development, providing communication relays and priming favorable microenvironments for cancer metastasis. By drawing comparisons between non-neural long distance signaling in the context of development and cancer we aim to encourage crosstalk between the two fields to cultivate new hypotheses and potential therapeutic strategies.

Genomic and physiological analyses of the zebrafish atrioventricular canal reveal molecular building blocks of the secondary pacemaker region

The atrioventricular canal (AVC) is the site where key structures responsible for functional division between heart regions are established, most importantly, the atrioventricular (AV) conduction system and cardiac valves. To elucidate the mechanism underlying AVC development and function, we utilized transgenic zebrafish line sqet31Et expressing EGFP in the AVC to isolate this cell population and profile its transcriptome at 48 and 72 hpf. The zebrafish AVC transcriptome exhibits hallmarks of mammalian AV node, including the expression of genes implicated in its development and those encoding connexins forming low conductance gap junctions. Transcriptome analysis uncovered protein-coding and noncoding transcripts enriched in AVC, which have not been previously associated with this structure, as well as dynamic expression of epithelial-to-mesenchymal transition markers and components of TGF-β, Notch, and Wnt signaling pathways likely reflecting ongoing AVC and valve development. Using transgenic line Tg(myl7:mermaid) encoding voltage-sensitive fluorescent protein, we show that abolishing the pacemaker-containing sinoatrial ring (SAR) through Isl1 loss of function resulted in spontaneous activation in the AVC region, suggesting that it possesses inherent automaticity although insufficient to replace the SAR. The SAR and AVC transcriptomes express partially overlapping species of ion channels and gap junction proteins, reflecting their distinct roles. Besides identifying conserved aspects between zebrafish and mammalian conduction systems, our results established molecular hallmarks of the developing AVC which underlies its role in structural and electrophysiological separation between heart chambers. This data constitutes a valuable resource for studying AVC development and function, and identification of novel candidate genes implicated in these processes.

Annotated expression and activity data for murine recombinase alleles and transgenes: the CrePortal resource

Recombinase alleles and transgenes can be used to facilitate spatio-temporal specificity of gene disruption or transgene expression. However, the versatility of this in vivo recombination system relies on having detailed and accurate characterization of recombinase expression and activity to enable selection of the appropriate allele or transgene. The CrePortal (http://www.informatics.jax.org/home/recombinase) leverages the informatics infrastructure of Mouse Genome Informatics to integrate data from the scientific literature, direct data submissions from the scientific community at-large, and from major projects developing new recombinase lines and characterizing recombinase expression and specificity patterns. Searching the CrePortal by recombinase activity or specific recombinase gene driver provides users with a recombinase alleles and transgenes activity tissue summary and matrix comparison of gene expression and recombinase activity with links to generation details, a recombinase activity grid, and associated phenotype annotations. Future improvements will add cell type-based activity annotations. The CrePortal provides a comprehensive presentation of recombinase allele and transgene data to assist researchers in selection of the recombinase allele or transgene based on where and when recombination is desired.

A systems approach using Diversity Outbred mice distinguishes the cardiovascular effects and genetics of circulating GDF11 from those of its homolog, myostatin

Growth differentiation factor 11 (GDF11) is a member of the TGF-β protein family that has been implicated in the development of cardiac hypertrophy. While some studies have suggested that systemic GDF11 protects against cardiomyocyte enlargement and left ventricular wall thickening, there remains uncertainty about the true impact of GDF11 and whether its purported effects are actually attributable to its homolog myostatin. This study was conducted to resolve the statistical and genetic relationships among GDF11, myostatin, and cardiac hypertrophy in a mouse model of human genetics, the Diversity Outbred (DO) stock. In the DO population, serum GDF11 concentrations positively correlated with cardiomyocyte cross-sectional area, while circulating myostatin levels were negatively correlated with body weight, heart weight, and left ventricular wall thickness and mass. Genetic analyses revealed that serum GDF11 concentrations are modestly heritable (0.23) and identified a suggestive peak on murine chromosome 3 in close proximity to the gene Hey1, a transcriptional repressor. Bioinformatic analyses located putative binding sites for the HEY1 protein upstream of the Gdf11 gene in the mouse and human genomes. In contrast, serum myostatin concentrations were more heritable (0.57) than GDF11 concentrations, and mapping identified a significant locus near the gene FoxO1, which has binding motifs within the promoter regions of human and mouse myostatin genes. Together, these findings more precisely define the independent cardiovascular effects of GDF11 and myostatin, as well as their distinct regulatory pathways. Hey1 is a compelling candidate for the regulation of GDF11 and will be further evaluated in future studies.

Coronavirus Disease 2019 and Liver Injury: A Retrospective Analysis of Hospitalized Patients in New York City

BACKGROUND AND AIMS

Coronavirus disease 2019 (COVID-19) is a global threat, affecting more than 100 million people and causing over 2 million deaths. Liver laboratory test abnormalities are an extrapulmonary manifestation of COVID-19, yet characterization of hepatic injury is incomplete. Our objective was to further characterize and identify causes of liver injury in patients with COVID-19.

METHODS

We conducted a retrospective cohort study of 551 patients hospitalized with COVID-19 at NewYork-Presbyterian Hospital/Columbia University Irving Medical Center between March 1, 2020 and May 31, 2020. We analyzed patient demographics, liver laboratory test results, vital signs, other relevant test results, and clinical outcomes (mortality and intensive care unit admission).

RESULTS

Abnormal liver laboratory tests were common on hospital admission for COVID-19 and the incidence increased during hospitalization. Of those with elevated serum alanine aminotransferase and/or alkaline phosphatase activities on admission, 58.2% had a cholestatic injury pattern, 35.2% mixed, and 6.6% hepatocellular. Comorbid liver disease was not associated with outcome; however, abnormal direct bilirubin or albumin on admission were associated with intensive care unit stay and mortality. On average, patients who died had greater magnitudes of abnormalities in all liver laboratory tests than those who survived. Ischemic hepatitis was a mechanism of severe hepatocellular injury in some patients.

CONCLUSIONS

Liver laboratory test abnormalities are common in hospitalized patients with COVID-19, and some are associated with increased odds of intensive care unit stay or death. Severe hepatocellular injury is likely attributable to secondary effects such as systemic inflammatory response syndrome, sepsis, and ischemic hepatitis.

STAT3 maintains skin barrier integrity by modulating SPINK5 and KLK5 expression in keratinocytes

Skin barrier dysfunction induces skin inflammation. Signal transducer and activator of transcription 3 (STAT3) is known to be involved in Th17-mediated immune responses and barrier integrity in the cornea and intestine; however, its role in the skin barrier remains largely unknown. In this study, we elucidated the potential role of STAT3 in the skin barrier and its effect on kallikrein-related peptidase 5 (KLK5) and serine protease inhibitor Kazal-type 5 (SPINK5) expression using a mouse model with keratinocyte-specific ablation of STAT3. Keratinocyte-specific loss of STAT3 induced a cutaneous inflammatory phenotype with pruritus and intense scratching behaviour in mice. Transcriptomic analysis revealed that the genes associated with impaired skin barrier function, including KLK5, were upregulated. The effect of STAT3 on KLK5 expression in keratinocytes was not only substantiated by the increase in KLK5 expression following treatment with STAT3 siRNA but also by its decreased expression following STAT3 overexpression. Overexpression and IL-17A-mediated stimulation of STAT3 increased the expression of SPINK5, which was blocked by STAT3 siRNA. These results suggest that the expression of SPINK5 and KLK5 in keratinocytes could be dependent on STAT3 and that STAT3 might play an essential role in the maintenance of skin barrier homeostasis.

Whole Blood Transcriptome Characterization of 3xTg-AD Mouse and Its Modulation by Transcranial Direct Current Stimulation (tDCS)

The 3xTg-AD mouse is a widely used model in the study of Alzheimer’s Disease (AD). It has been extensively characterized from both the anatomical and behavioral point of view, but poorly studied at the transcriptomic level. For the first time, we characterize the whole blood transcriptome of the 3xTg-AD mouse at three and six months of age and evaluate how its gene expression is modulated by transcranial direct current stimulation (tDCS). RNA-seq analysis revealed 183 differentially expressed genes (DEGs) that represent a direct signature of the genetic background of the mouse. Moreover, in the 6-month-old 3xTg-AD mice, we observed a high number of DEGs that could represent good peripheral biomarkers of AD symptomatology onset. Finally, tDCS was associated with gene expression changes in the 3xTg-AD, but not in the control mice. In conclusion, this study provides an in-depth molecular characterization of the 3xTg-AD mouse and suggests that blood gene expression can be used to identify new biomarkers of AD progression and treatment effects.

Crk adaptor proteins are necessary for the development of the zebrafish retina

BACKGROUND

The development of the central nervous system (CNS) requires critical cell signaling molecules to coordinate cell proliferation and migration in order to structure the adult tissue. Chicken tumor virus #10 Regulator of Kinase (CRK) and CRK-like (CRKL) are adaptor proteins with pre-metazoan ancestry and are known to be required for patterning laminated structures downstream of Reelin (RELN), such as the cerebral cortex, cerebellum, and hippocampus. CRK and CRKL also play crucial roles in a variety of other growth factor and extracellular matrix signaling cascades. The neuronal retina is another highly laminated structure within the CNS that is dependent on migration for proper development, but the cell signaling mechanisms behind neuronal positioning in the retina are only partly understood.

RESULTS

We find that crk and crkl have largely overlapping expression within the developing zebrafish nervous system. We find that their disruption results in smaller eye size and loss of retinal lamination.

CONCLUSIONS

Our data indicate that Crk adaptors are critical for proper development of the zebrafish neural retina in a crk/crkl dose-dependent manner.

Engrailed 2 deficiency and chronic stress alter avoidance and motivation behaviors

Autism spectrum disorder (ASD) is a pervasive neurodevelopmental disorder characterized by impairments in social interaction, cognition, and communication, as well as the presence of repetitive or stereotyped behaviors and interests. ASD is most often studied as a neurodevelopmental disease, but it is a lifelong disorder. Adults with ASD experience more stressful life events and greater perceived stress, and frequently have comorbid mood disorders such as anxiety and depression. It remains unclear whether adult exposure to chronic stress can exacerbate the behavioral and neurodevelopmental phenotypes associated with ASD. To address this issue, we first investigated whether adult male and female Engrailed-2 deficient (En2-KO, En2-/-) mice, which display behavioral disturbances in avoidance tasks and dysregulated monoaminergic neurotransmitter levels, also display impairments in instrumental behaviors associated with motivation, such as the progressive ratio task. We then exposed adult En2-KO mice to chronic environmental stress (CSDS, chronic social defeat stress), to determine if stress exacerbated the behavioral and neuroanatomical effects of En2 deletion. En2-/- mice showed impaired instrumental acquisition and significantly lower breakpoints in a progressive ratio test, demonstrating En2 deficiency decreases motivation to exert effort for reward. Furthermore, adult CSDS exposure increased avoidance behaviors in En2-KO mice. Interestingly, adult CSDS exposure also exacerbated the deleterious effects of En2 deficiency on forebrain-projecting monoaminergic fibers. Our findings thus suggest that adult exposure to stress may exacerbate behavioral and neuroanatomical phenotypes associated with developmental effects of genetic En2 deficiency.

Lifespan-Associated Gene Expression Signatures of Recombinant BXD Mice Implicates Coro7 and Set in Longevity

Although genetic approaches have identified key genes and pathways that promote longevity, systems-level approaches are less utilized. Here, we took advantage of the wealth of omics data characterizing the BXD family of mice. We associated transcript and peptide levels across five tissues from both female and male BXD isogenic lines with their median lifespan. We identified over 5000 genes that showed a longevity correlation in a given tissue. Surprisingly, we found less than 1% overlap among longevity-correlating genes across tissues and sex. These 1% shared genes consist of 51 genes, of which 13 have been shown to alter lifespan. Only two genes -Coro7 and Set- showed a longevity correlation in all tissues and in both sexes. While differential regulation of aging across tissues and sex has been reported, our systems-level analysis reveals two unique genes that may promote healthy aging in unique sex- and tissue-agnostic manner.

RGS14 Regulation of Post-Synaptic Signaling and Spine Plasticity in Brain

The regulator of G-protein signaling 14 (RGS14) is a multifunctional signaling protein that regulates post synaptic plasticity in neurons. RGS14 is expressed in the brain regions essential for learning, memory, emotion, and stimulus-induced behaviors, including the basal ganglia, limbic system, and cortex. Behaviorally, RGS14 regulates spatial and object memory, female-specific responses to cued fear conditioning, and environmental- and psychostimulant-induced locomotion. At the cellular level, RGS14 acts as a scaffolding protein that integrates G protein, Ras/ERK, and calcium/calmodulin signaling pathways essential for spine plasticity and cell signaling, allowing RGS14 to naturally suppress long-term potentiation (LTP) and structural plasticity in hippocampal area CA2 pyramidal cells. Recent proteomics findings indicate that RGS14 also engages the actomyosin system in the brain, perhaps to impact spine morphogenesis. Of note, RGS14 is also a nucleocytoplasmic shuttling protein, where its role in the nucleus remains uncertain. Balanced nuclear import/export and dendritic spine localization are likely essential for RGS14 neuronal functions as a regulator of synaptic plasticity. Supporting this idea, human genetic variants disrupting RGS14 localization also disrupt RGS14’s effects on plasticity. This review will focus on the known and unexplored roles of RGS14 in cell signaling, physiology, disease and behavior.

SOMDE: A scalable method for identifying spatially variable genes with self-organizing map

MOTIVATION

Recent developments of spatial transcriptomic sequencing technologies provide powerful tools for understanding cells in the physical context of tissue microenvironments. A fundamental task in spatial gene expression analysis is to identify genes with spatially variable expression patterns, or spatially variable genes (SVgenes). Several computational methods have been developed for this task. Their high computational complexity limited their scalability to the latest and future large-scale spatial expression data.

RESULTS

We present SOMDE, an efficient method for identifying SVgenes in large-scale spatial expression data. SOMDE uses self-organizing map (SOM) to cluster neighboring cells into nodes, and then uses a Gaussian process to fit the node-level spatial gene expression to identify SVgenes. Experiments show that SOMDE is about 5-50 times faster than existing methods with comparable results. The adjustable resolution of SOMDE makes it the only method that can give results in ∼5 minutes in large datasets of more than 20,000 sequencing sites. SOMDE is available as a python package on PyPI at https://pypi.org/project/somde free for academic use.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In Vitro, In Vivo, and In Silico Methods for Assessment of Muscle Size and Muscle Growth Regulation

Trauma, burn injury, sepsis, and ischemia lead to acute and chronic loss of skeletal muscle mass and function. Healthy muscle is essential for eating, posture, respiration, reproduction, and mobility, as well as for appropriate function of the senses including taste, vision, and hearing. Beyond providing support and contraction, skeletal muscle also exerts essential roles in temperature regulation, metabolism, and overall health. As the primary reservoir for amino acids, skeletal muscle regulates whole-body protein and glucose metabolism by providing substrate for protein synthesis and supporting hepatic gluconeogenesis during illness and starvation. Overall, greater muscle mass is linked to greater insulin sensitivity and glucose disposal, strength, power, and longevity. In contrast, low muscle mass correlates with dysmetabolism, dysmobility, and poor survival. Muscle mass is highly plastic, appropriate to its role as reservoir, and subject to striking genetic control. Defining mechanisms of muscle growth regulation holds significant promise to find interventions that promote health and diminish morbidity and mortality after trauma, sepsis, inflammation, and other systemic insults. In this invited review, we summarize techniques and methods to assess and manipulate muscle size and muscle mass in experimental systems, including cell culture and rodent models. These approaches have utility for studies of myopenia, sarcopenia, cachexia, and acute muscle growth or atrophy in the setting of health or injury.

Assessing evolutionary and developmental transcriptome dynamics in homologous cell types

Background

During development, complex organ patterns emerge through the precise temporal and spatial specification of different cell types. On an evolutionary timescale, these patterns can change, resulting in morphological diversification. It is generally believed that homologous anatomical structures are built—largely—by homologous cell types. However, whether a common evolutionary origin of such cell types is always reflected in the conservation of their intrinsic transcriptional specification programs is less clear.

Results

Here, we developed a user‐friendly bioinformatics workflow to detect gene co‐expression modules and test for their conservation across developmental stages and species boundaries. Using a paradigm of morphological diversification, the tetrapod limb, and single‐cell RNA‐sequencing data from two distantly related species, chicken and mouse, we assessed the transcriptional dynamics of homologous cell types during embryonic patterning. With mouse limb data as reference, we identified 19 gene co‐expression modules with varying tissue or cell type‐restricted activities. Testing for co‐expression conservation revealed modules with high evolutionary turnover, while others seemed maintained—to different degrees, in module make‐up, density or connectivity—over developmental and evolutionary timescales.

Conclusions

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

µ-Crystallin: A thyroid hormone binding protein

µ-Crystallin is a NADPH-regulated thyroid hormone binding protein encoded by the CRYM gene in humans. It is primarily expressed in the brain, muscle, prostate, and kidney, where it binds thyroid hormones, which regulate metabolism and thermogenesis. It also acts as a ketimine reductase in the lysine degradation pathway when it is not bound to thyroid hormone. Mutations in CRYM can result in non-syndromic deafness, while its aberrant expression, predominantly in the brain but also in other tissues, has been associated with psychiatric, neuromuscular, and inflammatory diseases. CRYM expression is highly variable in human skeletal muscle, with 15% of individuals expressing ≥13 fold more CRYM mRNA than the median level. Ablation of the Crym gene in murine models results in the hypertrophy of fast twitch muscle fibers and an increase in fat mass of mice fed a high fat diet. Overexpression of Crym in mice causes a shift in energy utilization away from glycolysis towards an increase in the catabolism of fat via β-oxidation, with commensurate changes of metabolically involved transcripts and proteins. The history, attributes, functions, and diseases associated with CRYM, an important modulator of metabolism, are reviewed.

In Vivo Cell Fate Tracing Provides No Evidence for Mesenchymal to Epithelial Transition in Adult Fallopian Tube and Uterus

The mesenchymal to epithelial transition (MET) is thought to be involved in the maintenance, repair, and carcinogenesis of the fallopian tube (oviduct) and uterine epithelium. However, conclusive evidence for the conversion of mesenchymal cells to epithelial cells in these organs is lacking. Using embryonal cell lineage tracing with reporters driven by mesenchymal cell marker genes of the female reproductive tract (AMHR2, CSPG4, and PDGFRβ), we show that these reporters are also expressed by some oviductal and uterine epithelial cells at birth. These mesenchymal reporter-positive epithelial cells are maintained in adult mice across multiple pregnancies, respond to ovarian hormones, and form organoids. However, no labeled epithelial cells are present in any oviductal or uterine epithelia when mesenchymal cell labeling was induced in adult mice. Organoids developed from mice labeled in adulthood were also negative for mesenchymal reporters. Collectively, our work found no definitive evidence of MET in the adult fallopian tube and uterine epithelium.

Mesenchymal to epithelial transition (MET) is postulated to be involved in the maintenance and regeneration of the epithelium of female reproductive organs. Here, Ghosh et al. report no definitive evidence of MET in the adult epithelium of oviduct and uterus using in vivo cell lineage tracing and organoids.

Prolactin-Releasing Peptide Differentially Regulates Gene Transcriptomic Profiles in Mouse Bone Marrow-Derived Macrophages

Prolactin-releasing Peptide (PrRP) is a neuropeptide whose receptor is GPR10. Recently, the regulatory role of PrRP in the neuroendocrine field has attracted increasing attention. However, the influence of PrRP on macrophages, the critical housekeeper in the neuroendocrine field, has not yet been fully elucidated. Here, we investigated the effect of PrRP on the transcriptome of mouse bone marrow-derived macrophages (BMDMs) with RNA sequencing, bioinformatics, and molecular simulation. BMDMs were exposed to PrRP (18 h) and were subjected to RNA sequencing. Differentially expressed genes (DEGs) were acquired, followed by GO, KEGG, and PPI analysis. Eight qPCR-validated DEGs were chosen as hub genes. Next, the three-dimensional structures of the proteins encoded by these hub genes were modeled by Rosetta and Modeller, followed by molecular dynamics simulation by the Gromacs program. Finally, the binding modes between PrRP and hub proteins were investigated with the Rosetta program. PrRP showed no noticeable effect on the morphology of macrophages. A total of 410 DEGs were acquired, and PrRP regulated multiple BMDM-mediated functional pathways. Besides, the possible docking modes between PrRP and hub proteins were investigated. Moreover, GPR10 was expressed on the cell membrane of BMDMs, which increased after PrRP exposure. Collectively, PrRP significantly changed the transcriptome profile of BMDMs, implying that PrRP may be involved in various physiological activities mastered by macrophages.

Conditional Single Vector CRISPR/SaCas9 Viruses for Efficient Mutagenesis in the Adult Mouse Nervous System

Mice engineered for conditional, cell type-specific gene inactivation have dominated the field of mouse genetics because of the high efficiency of Cre-loxP-mediated recombination. Recent advances in CRISPR/Cas9 technologies have provided alternatives for rapid gene mutagenesis for loss-of-function (LOF) analysis. Whether these strategies can be streamlined for rapid genetic analysis with the efficiencies comparable with those of conventional genetic approaches has yet to be established. We show that a single adeno-associated viral (AAV) vector containing a recombinase-dependent Staphylococcus aureus Cas9 (SaCas9) and a single guide RNA (sgRNA) are as efficient as conventional conditional gene knockout and can be adapted for use in either Cre- or Flp-driver mouse lines. The efficacy of this approach is demonstrated for the analysis of GABAergic, glutamatergic, and monoaminergic neurotransmission. Using this strategy, we reveal insight into the role of GABAergic regulation of midbrain GABA-producing neurons in psychomotor activation.

Hunker et al. generate single adeno-associated viral vectors for conditional gene mutagenesis in the adult mouse nervous system with efficiencies equivalent to those of conventional gene inactivation strategies. On the basis of this efficacy, they provide a resource of Cre- and Flp-dependent constructs for targeting catecholamine, glutamate, and GABA neurotransmission.

Mechanistic Drivers of Müllerian Duct Development and Differentiation Into the Oviduct

The conduits of life; the animal oviducts and human fallopian tubes are of paramount importance for reproduction in amniotes. They connect the ovary with the uterus and are essential for fertility. They provide the appropriate environment for gamete maintenance, fertilization and preimplantation embryonic development. However, serious pathologies, such as ectopic pregnancy, malignancy and severe infections, occur in the oviducts. They can have drastic effects on fertility, and some are life-threatening. Despite the crucial importance of the oviducts in life, relatively little is known about the molecular drivers underpinning the embryonic development of their precursor structures, the Müllerian ducts, and their successive differentiation and maturation. The Müllerian ducts are simple rudimentary tubes comprised of an epithelial lumen surrounded by a mesenchymal layer. They differentiate into most of the adult female reproductive tract (FRT). The earliest sign of Müllerian duct formation is the thickening of the anterior mesonephric coelomic epithelium to form a placode of two distinct progenitor cells. It is proposed that one subset of progenitor cells undergoes partial epithelial-mesenchymal transition (pEMT), differentiating into immature Müllerian luminal cells, and another subset undergoes complete EMT to become Müllerian mesenchymal cells. These cells invaginate and proliferate forming the Müllerian ducts. Subsequently, pEMT would be reversed to generate differentiated epithelial cells lining the fully formed Müllerian lumen. The anterior Müllerian epithelial cells further specialize into the oviduct epithelial subtypes. This review highlights the key established molecular and genetic determinants of the processes involved in Müllerian duct development and the differentiation of its upper segment into oviducts. Furthermore, an extensive genome-wide survey of mouse knockout lines displaying Müllerian or oviduct phenotypes was undertaken. In addition to widely established genetic determinants of Müllerian duct development, our search has identified surprising associations between loss-of-function of several genes and high-penetrance abnormalities in the Müllerian duct and/or oviducts. Remarkably, these associations have not been investigated in any detail. Finally, we discuss future directions for research on Müllerian duct development and oviducts.