The Mouse Genome Database (MGD): facilitating mouse as a model for human biology and disease

SMCHD1 mutations associated with a rare muscular dystrophy can also cause isolated arhinia and Bosma arhinia microphthalmia syndrome

Arhinia, or absence of the nose, is a rare malformation of unknown etiology that is often accompanied by ocular and reproductive defects. Sequencing of 40 people with arhinia revealed that 84% of probands harbor a missense mutation localized to a constrained region of SMCHD1 encompassing the ATPase domain. SMCHD1 mutations cause facioscapulohumeral muscular dystrophy type 2 (FSHD2) via a trans-acting loss-of-function epigenetic mechanism. We discovered shared mutations and comparable DNA hypomethylation patterning between these distinct disorders. CRISPR/Cas9-mediated alteration of smchd1 in zebrafish yielded arhinia-relevant phenotypes. Transcriptome and protein analyses in arhinia probands and controls showed no differences in SMCHD1 mRNA or protein abundance but revealed regulatory changes in genes and pathways associated with craniofacial patterning. Mutations in SMCHD1 thus contribute to distinct phenotypic spectra, from craniofacial malformation and reproductive disorders to muscular dystrophy, which we speculate to be consistent with oligogenic mechanisms resulting in pleiotropic outcomes.

An Asymmetrically Balanced Organization of Kinases versus Phosphatases across Eukaryotes Determines Their Distinct Impacts

Protein phosphorylation underlies cellular response pathways across eukaryotes and is governed by the opposing actions of phosphorylating kinases and de-phosphorylating phosphatases. While kinases and phosphatases have been extensively studied, their organization and the mechanisms by which they balance each other are not well understood. To address these questions we performed quantitative analyses of large-scale 'omics' datasets from yeast, fly, plant, mouse and human. We uncovered an asymmetric balance of a previously-hidden scale: Each organism contained many different kinase genes, and these were balanced by a small set of highly abundant phosphatase proteins. Kinases were much more responsive to perturbations at the gene and protein levels. In addition, kinases had diverse scales of phenotypic impact when manipulated. Phosphatases, in contrast, were stable, highly robust and flatly organized, with rather uniform impact downstream. We validated aspects of this organization experimentally in nematode, and supported additional aspects by theoretic analysis of the dynamics of protein phosphorylation. Our analyses explain the empirical bias in the protein phosphorylation field toward characterization and therapeutic targeting of kinases at the expense of phosphatases. We show quantitatively and broadly that this is not only a historical bias, but stems from wide-ranging differences in their organization and impact. The asymmetric balance between these opposing regulators of protein phosphorylation is also common to opposing regulators of two other post-translational modification systems, suggesting its fundamental value.

Protein phosphorylation is a reversible modification that underlies cellular responses to stimuli across organisms. Historically, the study of protein phosphorylation concentrated on the role of kinases, which introduce the phosphate, at the expense of phosphatases, which remove it. Many kinases have been associated with specific phenotypes and considered attractive drug targets, while phosphatases remained far less characterized. It has been unclear whether this discrepancy is due to historical biases or reflects real systemic differences between these enzymes. By analyzing large-scale ‘omics’ datasets across genes, transcripts, proteins, interactions, and organisms, we uncovered an asymmetric architecture of kinases versus phosphatases that balances between them, determines their distinct impact patterns, and affects their therapeutic potential. This architecture is conserved from yeast to human and is partially shared by two other protein modification systems, suggesting it is a general feature of these systems.

Identification of a Putative Quantitative Trait Gene for Resistance to Obesity in Mice Using Transcriptome Analysis and Causal Inference Tests

It is still challenging to identify causal genes governing obesity. Pbwg1.5, a quantitative trait locus (QTL) for resistance to obesity, was previously discovered from wild Mus musculus castaneus mice and was fine-mapped to a 2.1-Mb genomic region of mouse chromosome 2, where no known gene with an effect on white adipose tissue (WAT) has been reported. The aim of this study was to identify a strong candidate gene for Pbwg1.5 by an integration approach of transcriptome analysis (RNA-sequencing followed by real-time PCR analysis) and the causal inference test (CIT), a statistical method to infer causal relationships between diplotypes, gene expression and trait values. Body weight, body composition and biochemical traits were measured in F2 mice obtained from an intercross between the C57BL/6JJcl strain and a congenic strain carrying Pbwg1.5 on the C57BL/6JJcl background. The F2 mice showed significant diplotype differences in 12 traits including body weight, WAT weight and serum cholesterol/triglyceride levels. The transcriptome analysis revealed that Ly75, Pla2r1, Fap and Gca genes were differentially expressed in the liver and that Fap, Ifih1 and Grb14 were differentially expressed in WAT. However, CITs indicated statistical evidence that only the liver Ly75 gene mediated between genotype and WAT. Ly75 expression was negatively associated with WAT weight. The results suggested that Ly75 is a putative quantitative trait gene for the obesity-resistant Pbwg1.5 QTL discovered from the wild M. m. castaneus mouse. The finding provides a novel insight into a better understanding of the genetic basis for prevention of obesity.

RAIN: RNA-protein Association and Interaction Networks

Protein association networks can be inferred from a range of resources including experimental data, literature mining and computational predictions. These types of evidence are emerging for non-coding RNAs (ncRNAs) as well. However, integration of ncRNAs into protein association networks is challenging due to data heterogeneity. Here, we present a database of ncRNA-RNA and ncRNA-protein interactions and its integration with the STRING database of protein-protein interactions. These ncRNA associations cover four organisms and have been established from curated examples, experimental data, interaction predictions and automatic literature mining. RAIN uses an integrative scoring scheme to assign a confidence score to each interaction. We demonstrate that RAIN outperforms the underlying microRNA-target predictions in inferring ncRNA interactions. RAIN can be operated through an easily accessible web interface and all interaction data can be downloaded.Database URL: http://rth.dk/resources/rain.

Multiplex shRNA Screening of Germ Cell Development by in Vivo Transfection of Mouse Testis

Spermatozoa are one of the few mammalian cell types that cannot be fully derived in vitro, severely limiting the application of modern genomic techniques to study germ cell biology. The current gold standard approach of characterizing single-gene knockout mice is slow as generation of each mutant line can take 6–9 months. Here, we describe an in vivo approach to rapid functional screening of germline genes based on a new nonsurgical, nonviral in vivo transfection method to deliver nucleic acids into testicular germ cells. By coupling multiplex transfection of short hairpin RNA (shRNA) constructs with pooled amplicon sequencing as a readout, we were able to screen many genes for spermatogenesis function in a quick and inexpensive experiment. We transfected nine mouse testes with a pilot pool of RNA interference (RNAi) against well-characterized genes to show that this system is highly reproducible and accurate. With a false negative rate of 18% and a false positive rate of 12%, this method has similar performance as other RNAi screens in the well-described Drosophila model system. In a separate experiment, we screened 26 uncharacterized genes computationally predicted to be essential for spermatogenesis and found numerous candidates for follow-up studies. Finally, as a control experiment, we performed a long-term selection screen in neuronal N2a cells, sampling shRNA frequencies at five sequential time points. By characterizing the effect of both libraries on N2a cells, we show that our screening results from testis are tissue-specific. Our calculations indicate that the current implementation of this approach could be used to screen thousands of protein-coding genes simultaneously in a single mouse testis. The experimental protocols and analysis scripts provided will enable other groups to use this procedure to study diverse aspects of germ cell biology ranging from epigenetics to cell physiology. This approach also has great promise as an applied tool for validating diagnoses made from medical genome sequencing, or designing synthetic biological sequences that can act as potent and highly specific male contraceptives.

Mouse Genome Informatics (MGI): Resources for Mining Mouse Genetic, Genomic, and Biological Data in Support of Primary and Translational Research

The Mouse Genome Informatics (MGI), resource ( www.informatics.jax.org ) has existed for over 25 years, and over this time its data content, informatics infrastructure, and user interfaces and tools have undergone dramatic changes (Eppig et al., Mamm Genome 26:272-284, 2015). Change has been driven by scientific methodological advances, rapid improvements in computational software, growth in computer hardware capacity, and the ongoing collaborative nature of the mouse genomics community in building resources and sharing data. Here we present an overview of the current data content of MGI, describe its general organization, and provide examples using simple and complex searches, and tools for mining and retrieving sets of data.

Heat Shock Proteins and Maternal Contribution to Oogenesis and Early Embryogenesis

Early embryos develop from fertilized eggs using materials that are stored during oocyte growth and which can be defined as maternal contribution (molecules, factors, or determinants). Several heat shock proteins (HSPs) and the heat shock transcriptional factor (HSF) are part of the maternal contribution that is critical for successful embryogenesis and reproduction. A maternal role for heat shock-related genes was mainly demonstrated in genetic experimental organisms (e.g., fly, nematode, mouse). Nowadays, an increasing number of "omics" data are produced from a large panel of organisms implementing a catalog of maternal and/or embryonic HSPs and HSFs. However, for most of them, it remains to better understand their potential roles in this context. Existing and future genome-wide screens mainly set up to create loss-of-function are likely to improve this situation. This chapter will discuss available data from various experimental organisms following the developmental steps from egg production to early embryogenesis.

Ensembl core software resources: storage and programmatic access for DNA sequence and genome annotation

The Ensembl software resources are a stable infrastructure to store, access and manipulate genome assemblies and their functional annotations. The Ensembl 'Core' database and Application Programming Interface (API) was our first major piece of software infrastructure and remains at the centre of all of our genome resources. Since its initial design more than fifteen years ago, the number of publicly available genomic, transcriptomic and proteomic datasets has grown enormously, accelerated by continuous advances in DNA-sequencing technology. Initially intended to provide annotation for the reference human genome, we have extended our framework to support the genomes of all species as well as richer assembly models. Cross-referenced links to other informatics resources facilitate searching our database with a variety of popular identifiers such as UniProt and RefSeq. Our comprehensive and robust framework storing a large diversity of genome annotations in one location serves as a platform for other groups to generate and maintain their own tailored annotation. We welcome reuse and contributions: our databases and APIs are publicly available, all of our source code is released with a permissive Apache v2.0 licence at http://github.com/Ensembl and we have an active developer mailing list ( http://www.ensembl.org/info/about/contact/index.html ).

Database URL

http://www.ensembl.org.

Protein Bioinformatics Databases and Resources

Many publicly available data repositories and resources have been developed to support protein-related information management, data-driven hypothesis generation, and biological knowledge discovery. To help researchers quickly find the appropriate protein-related informatics resources, we present a comprehensive review (with categorization and description) of major protein bioinformatics databases in this chapter. We also discuss the challenges and opportunities for developing next-generation protein bioinformatics databases and resources to support data integration and data analytics in the Big Data era.

Global Analysis of Protein Expression of Inner Ear Hair Cells

The mammalian inner ear (IE) subserves auditory and vestibular sensations via highly specialized cells and proteins. Sensory receptor hair cells (HCs) are necessary for transducing mechanical inputs and stimulating sensory neurons by using a host of known and as yet unknown protein machinery. To understand the protein composition of these unique postmitotic cells, in which irreversible protein degradation or damage can lead to impaired hearing and balance, we analyzed IE samples by tandem mass spectrometry to generate an unbiased, shotgun-proteomics view of protein identities and abundances. By using Pou4f3/eGFP-transgenic mice in which HCs express GFP driven by Pou4f3, we FACS purified a population of HCs to analyze and compare the HC proteome with other IE subproteomes from sensory epithelia and whole IE. We show that the mammalian HC proteome comprises hundreds of uniquely or highly expressed proteins. Our global proteomic analysis of purified HCs extends the existing HC transcriptome, revealing previously undetected gene products and isoform-specific protein expression. Comparison of our proteomic data with mouse and human databases of genetic auditory/vestibular impairments confirms the critical role of the HC proteome for normal IE function, providing a cell-specific pool of candidates for novel, important HC genes. Several proteins identified exclusively in HCs by proteomics and verified by immunohistochemistry map to human genetic deafness loci, potentially representing new deafness genes.

SIGNIFICANCE STATEMENT

Hearing and balance rely on specialized sensory hair cells (HCs) in the inner ear (IE) to convey information about sound, acceleration, and orientation to the brain. Genetically and environmentally induced perturbations to HC proteins can result in deafness and severe imbalance. We used transgenic mice with GFP-expressing HCs, coupled with FACS sorting and tandem mass spectrometry, to define the most complete HC and IE proteome to date. We show that hundreds of proteins are uniquely identified or enriched in HCs, extending previous gene expression analyses to reveal novel HC proteins and isoforms. Importantly, deafness-linked proteins were significantly enriched in HCs, suggesting that this in-depth proteomic analysis of IE sensory cells may hold potential for deafness gene discovery.

Normalization of Hepatic Homeostasis in the Npc1nmf164 Mouse Model of Niemann-Pick Type C Disease Treated with the Histone Deacetylase Inhibitor Vorinostat

Niemann-Pick type C (NP-C) disease is a fatal genetic lipidosis for which there is no Food and Drug Administration (FDA)-approved therapy. Vorinostat, an FDA-approved inhibitor of histone deacetylases, ameliorates lysosomal lipid accumulation in cultured NP-C patient fibroblasts. To assess the therapeutic potential of histone deacetylase inhibition, we pursued these in vitro observations in two murine models of NP-C disease. Npc1^nmf164 mice, which express a missense mutation in the Npc1 gene, were treated intraperitoneally, from weaning, with the maximum tolerated dose of vorinostat (150 mg/kg, 5 days/week). Disease progression was measured via gene expression, liver function and pathology, serum and tissue lipid levels, body weight, and life span. Transcriptome analyses of treated livers indicated multiple changes consistent with reversal of liver dysfunction that typifies NP-C disease. Significant improvements in liver pathology and function were achieved by this treatment regimen; however, NPC1 protein maturation and levels, disease progression, weight loss, and animal morbidity were not detectably altered. Vorinostat concentrations were >200 μm in the plasma compartment of treated animals but were almost 100-fold lower in brain tissue. Apolipoprotein B metabolism and the expression of key components of lipid homeostasis in primary hepatocytes from null (Npc1^-/-) and missense (Npc1^nmf164 ) mutant mice were altered by vorinostat treatment, consistent with a response by these cells independent of the status of the Npc1 locus. These results suggest that HDAC inhibitors have utility to treat visceral NP-C disease. However, it is clear that improved blood-brain barrier penetration will be required to alleviate the neurological symptoms of human NP-C disease.

Accelerating Adverse Outcome Pathway Development Using Publicly Available Data Sources

The adverse outcome pathway (AOP) concept links molecular perturbations with organism and population-level outcomes to support high-throughput toxicity (HTT) testing. International efforts are underway to define AOPs and store the information supporting these AOPs in a central knowledge base; however, this process is currently labor-intensive and time-consuming. Publicly available data sources provide a wealth of information that could be used to define computationally predicted AOPs (cpAOPs), which could serve as a basis for creating expert-derived AOPs in a much more efficient way. Computational tools for mining large datasets provide the means for extracting and organizing the information captured in these public data sources. Using cpAOPs as a starting point for expert-derived AOPs should accelerate AOP development. Coupling this with tools to coordinate and facilitate the expert development efforts will increase the number and quality of AOPs produced, which should play a key role in advancing the adoption of HTT testing, thereby reducing the use of animals in toxicity testing and greatly increasing the number of chemicals that can be tested.

Mobile-phone radiation-induced perturbation of gene-expression profiling, redox equilibrium and sporadic-apoptosis control in the ovary of Drosophila melanogaster

The daily use by people of wireless communication devices has increased exponentially in the last decade, begetting concerns regarding its potential health hazards. Drosophila melanogaster four days-old adult female flies were exposed for 30 min to radiation emitted by a commercial mobile phone at a SAR of 0.15 W/kg and a SAE of 270 J/kg. ROS levels and apoptotic follicles were assayed in parallel with a genome-wide microarrays analysis. ROS cellular contents were found to increase by 1.6-fold (x), immediately after the end of exposure, in follicles of pre-choriogenic stages (germarium - stage 10), while sporadically generated apoptotic follicles (germarium 2b and stages 7-9) presented with an averaged 2x upregulation in their sub-population mass, 4 h after fly's irradiation with mobile device. Microarray analysis revealed 168 genes being differentially expressed, 2 h post-exposure, in response to radiofrequency (RF) electromagnetic field-radiation exposure (≥1.25x, P < 0.05) and associated with multiple and critical biological processes, such as basic metabolism and cellular subroutines related to stress response and apoptotic death. Exposure of adult flies to mobile-phone radiation for 30 min has an immediate impact on ROS production in animal's ovary, which seems to cause a global, systemic and non-targeted transcriptional reprogramming of gene expression, 2 h post-exposure, being finally followed by induction of apoptosis 4 h after the end of exposure. Conclusively, this unique type of pulsed radiation, mainly being derived from daily used mobile phones, seems capable of mobilizing critical cytopathic mechanisms, and altering fundamental genetic programs and networks in D. melanogaster.

Searching the Mouse Genome Informatics (MGI) Resources for Information on Mouse Biology from Genotype to Phenotype

The Mouse Genome Informatics (MGI) resource provides the research community with access to information on the genetics, genomics, and biology of the laboratory mouse. Core data in MGI include gene characterization and function, phenotype and disease model descriptions, DNA and protein sequence data, gene expression data, vertebrate homologies, SNPs, mapping data, and links to other bioinformatics databases. Semantic integration is supported through the use of standardized nomenclature, and through the use of controlled vocabularies such as the mouse Anatomical Dictionary, the Mammalian Phenotype Ontology, and the Gene Ontologies. MGI extracts and organizes data from primary literature. MGI data are shared with and widely displayed from other bioinformatics resources. The database is updated weekly with curated annotations, and regularly adds new datasets and features. This unit provides a guide to using the MGI bioinformatics resource. © 2016 by John Wiley & Sons, Inc.

Concise Review: Getting to the Core of Inherited Bone Marrow Failures

Bone marrow failure syndromes (BMFS) are a group of disorders with complex pathophysiology characterized by a common phenotype of peripheral cytopenia and/or hypoplastic bone marrow. Understanding genetic factors contributing to the pathophysiology of BMFS has enabled the identification of causative genes and development of diagnostic tests. To date more than 40 mutations in genes involved in maintenance of genomic stability, DNA repair, ribosome and telomere biology have been identified. In addition, pathophysiological studies have provided insights into several biological pathways leading to the characterization of genotype/phenotype correlations as well as the development of diagnostic approaches and management strategies. Recent developments in bone marrow transplant techniques and the choice of conditioning regimens have helped improve transplant outcomes. However, current morbidity and mortality remain unacceptable underlining the need for further research in this area. Studies in mice have largely been unable to mimic disease phenotype in humans due to difficulties in fully replicating the human mutations and the differences between mouse and human cells with regard to telomere length regulation, processing of reactive oxygen species and lifespan. Recent advances in induced pluripotency have provided novel insights into disease pathogenesis and have generated excellent platforms for identifying signaling pathways and functional mapping of haplo‐insufficient genes involved in large‐scale chromosomal deletions–associated disorders. In this review, we have summarized the current state of knowledge in the field of BMFS with specific focus on modeling the inherited forms and how to best utilize these models for the development of targeted therapies. Stem Cells2017;35:284–298

Drug repurposing for glioblastoma based on molecular subtypes

A recent multi-platform analysis by The Cancer Genome Atlas identified four distinct molecular subtypes for glioblastoma (GBM) and demonstrated that the subtypes correlate with clinical phenotypes and treatment responses. In this study, we developed a computational drug repurposing approach to predict GBM drugs based on the molecular subtypes. Our approach leverages the genomic signature for each GBM subtype, and integrates the human cancer genomics with mouse phenotype data to identify the opportunity of reusing the FDA-approved agents to treat specific GBM subtypes. Specifically, we first constructed the phenotype profile for each GBM subtype using their genomic signatures. For each approved drug, we also constructed a phenotype profile using the drug target genes. Then we developed an algorithm to match and prioritize drugs based on their phenotypic similarities to the GBM subtypes. Our approach is highly generalizable for other disorders if provided with a list of disorder-specific genes. We first evaluated the approach in predicting drugs for the whole GBM. For a combined set of approved, potential and off-label GBM drugs, we achieved a median rank of 9.3%, which is significantly higher (p Collapse

Key Words

• Drug discovery
• Glioblastoma
• Molecular subtypes
• Precision medicine

Collapse

MESH Headings

• Algorithms
• Animals
• Brain Neoplasms/drug therapy
• Brain Neoplasms/genetics
• Drug Repositioning
• Genomics
• Glioblastoma/drug therapy
• Glioblastoma/genetics
• Humans
• Mice
• Phenotype

Collapse

Grants

• DP2 HD084068 NICHD NIH HHS
• UL1 TR000439 NCATS NIH HHS

Collapse

Affiliation(s)

Yang Chen Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, OH 44106, United States
Rong Xu Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, OH 44106, United States.

Collapse

Loss of the arginine methyltranserase PRMT7 causes syndromic intellectual disability with microcephaly and brachydactyly

Post-translational protein modifications exponentially expand the functional complement of proteins encoded by the human genome. One such modification is the covalent addition of a methyl group to arginine or lysine residues, which is used to regulate a substantial proportion of the proteome. Arginine and lysine methylation are catalyzed by protein arginine methyltransferase (PRMTs) and protein lysine methyltransferase proteins (PKMTs), respectively; each methyltransferase has a specific set of target substrates. Here, we report a male with severe intellectual disability, facial dysmorphism, microcephaly, short stature, brachydactyly, cryptorchidism and seizures who was found to have a homozygous 15,309 bp deletion encompassing the transcription start site of PRMT7, which we confirmed is functionally a null allele. We show that the patient's cells have decreased levels of protein arginine methylation, and that affected proteins include the essential histones, H2B and H4. Finally, we demonstrate that patient cells have altered Wnt signaling, which may have contributed to the skeletal abnormalities. Our findings confirm the recent disease association of PRMT7, expand the phenotypic manifestations of this disorder and provide insight into the molecular pathogenesis of this new condition.

Influence of early life exposure, host genetics and diet on the mouse gut microbiome and metabolome

Although the gut microbiome plays important roles in host physiology, health and disease¹, we lack understanding of the complex interplay between host genetics and early life environment on the microbial and metabolic composition of the gut. We used the genetically diverse Collaborative Cross mouse system² to discover that early life history impacts the microbiome composition, whereas dietary changes have only a moderate effect. By contrast, the gut metabolome was shaped mostly by diet, with specific non-dietary metabolites explained by microbial metabolism. Quantitative trait analysis identified mouse genetic trait loci (QTL) that impact the abundances of specific microbes. Human orthologues of genes in the mouse QTL are implicated in gastrointestinal cancer. Additionally, genes located in mouse QTL for Lactobacillales abundance are implicated in arthritis, rheumatic disease and diabetes. Furthermore, Lactobacillales abundance was predictive of higher host T-helper cell counts, suggesting an important link between Lactobacillales and host adaptive immunity.

Effects of chronic restraint-induced stress on radiation-induced chromosomal aberrations in mouse splenocytes

Both ionizing radiation (IR) and psychological stress (PS) cause detrimental effects on humans. A recent study showed that chronic restraint-induced PS (CRIPS) diminished the functions of Trp53 and enhanced radiocarcinogenesis in Trp53-heterozygous (Trp53^+/-) mice. These findings had a marked impact on the academic field as well as the general public, particularly among residents living in areas radioactively contaminated by nuclear accidents. In an attempt to elucidate the modifying effects of CRIPS on radiation-induced health consequences in Trp53 wild-type (Trp53^+/+) animals, investigations involving multidisciplinary analyses were performed. We herein demonstrated that CRIPS induced changes in the frequency of IR-induced chromosomal aberrations (CAs) in splenocytes. Five-week-old male Trp53^+/+ C57BL/6J mice were restrained for 6h per day for 28 consecutive days, and total body irradiation (TBI) at a dose of 4Gy was performed on the 8th day. Metaphase chromosome spreads prepared from splenocytes at the end of the 28-day restraint regimen were painted with fluorescence in situ hybridization (FISH) probes for chromosomes 1, 2, and 3. The results obtained showed that CRIPS alone did not induce CAs, while TBI caused significant increases in CAs, mostly translocations. Translocations appeared at a lower frequency in mice exposed to TBI plus CRIPS than in those exposed to TBI alone. No significant differences were observed in the frequencies of the other types of CAs (insertions, dicentrics, and fragments) visualized with FISH between these experimental groups (TBI+CRIPS vs. TBI). These results suggest that CRIPS does not appear to synergize with the clastogenicity of IR.

Towards Functional Annotation of the Preimplantation Transcriptome: An RNAi Screen in Mammalian Embryos

With readily available transcriptome-wide data, understanding the role of each expressed gene is an essential next step. Although RNAi technologies allow for genome-wide screens in cell culture, these approaches cannot replace strategies for discovery in the embryo. Here we present, for the first time, a knockdown screen in mouse preimplantation embryos. Early mammalian development encompasses dynamic cellular, molecular and epigenetic events that are largely conserved from mouse to man. We assayed 712 genes for requirements during preimplantation. We identified 59 genes required for successful development or outgrowth and implantation. We have characterized each phenotype and revealed cellular, molecular, and lineage specific defects following knockdown of transcript. Induced network analyses demonstrate this as a valid approach to identify networks of genes that play important roles during preimplantation. Our approach provides a robust and efficient strategy towards identification of novel phenotypes during mouse preimplantation and facilitates functional annotation of the mammalian transcriptome.

Discovery of Age-Related Protein Folding Stability Differences in the Mouse Brain Proteome

Described here is the application of thermodynamic stability measurements to study age-related differences in the folding and stability of proteins in a rodent model of aging. Thermodynamic stability profiles were generated for 809 proteins in brain cell lysates from mice, aged 6 (n = 7) and 18 months (n = 9) using the Stability of Proteins from Rates of Oxidation (SPROX) technique. The biological variability of the protein stability measurements was low and within the experimental error of SPROX. A total of 83 protein hits were detected with age-related stability differences in the brain samples. Remarkably, the large majority of the brain protein hits were destabilized in the old mice, and the hits were enriched in proteins that have slow turnover rates (p < 0.07). Furthermore, 70% of the hits have been previously linked to aging or age-related diseases. These results help validate the use of thermodynamic stability measurements to capture relevant age-related proteomic changes and establish a new biophysical link between these proteins and aging.

GEN3VA: aggregation and analysis of gene expression signatures from related studies

BACKGROUND

Genome-wide gene expression profiling of mammalian cells is becoming a staple of many published biomedical and biological research studies. Such data is deposited into data repositories such as the Gene Expression Omnibus (GEO) for potential reuse. However, these repositories currently do not provide simple interfaces to systematically analyze collections of related studies.

RESULTS

Here we present GENE Expression and Enrichment Vector Analyzer (GEN3VA), a web-based system that enables the integrative analysis of aggregated collections of tagged gene expression signatures identified and extracted from GEO. Each tagged collection of signatures is presented in a report that consists of heatmaps of the differentially expressed genes; principal component analysis of all signatures; enrichment analysis with several gene set libraries across all signatures, which we term enrichment vector analysis; and global mapping of small molecules that are predicted to reverse or mimic each signature in the aggregate. We demonstrate how GEN3VA can be used to identify common molecular mechanisms of aging by analyzing tagged signatures from 244 studies that compared young vs. old tissues in mammalian systems. In a second case study, we collected 86 signatures from treatment of human cells with dexamethasone, a glucocorticoid receptor (GR) agonist. Our analysis confirms consensus GR target genes and predicts potential drug mimickers.

CONCLUSIONS

GEN3VA can be used to identify, aggregate, and analyze themed collections of gene expression signatures from diverse but related studies. Such integrative analyses can be used to address concerns about data reproducibility, confirm results across labs, and discover new collective knowledge by data reuse. GEN3VA is an open-source web-based system that is freely available at: http://amp.pharm.mssm.edu/gen3va .

Genome Sequencing of Chromosome 1 Substitution Lines Derived from Chinese Wild Mice Revealed a Unique Resource for Genetic Studies of Complex Traits

Mouse resources such as Collaborative Cross, outbred stocks, Hybrid Mouse Diversity Panel, and chromosome substitution strains have been instrumental to many progresses in the studies of complex traits genetics. We have established a population of chromosome 1 (Chr 1) substitution lines (C1SLs) in which donor chromosomes were derived from Chinese wild mice. Genome sequencing of 18 lines of this population showed that Chr 1 had been replaced by the donor chromosome. About 4.5 million unique single nucleotide polymorphisms and indels were discovered on Chr 1, of which 1.3 million were novel. Compared with sequenced classical inbred strains, Chr 1 of each C1SL had fivefold more variants, and more loss of function and potentially regulatory variants. Further haplotype analysis showed that the donor chromosome accumulated more historical recombination events, with the largest haplotype block being only 100 kb, and about 57% of the blocks were <1 kb. Subspecies origin analysis showed that these chromosomes had a mosaic genome structure that dominantly originated from Mus musculus musculus and M. m. castaneus subspecies, except for the C57BL/6J-Chr1^KM line from M. m. domesticus. In addition, phenotyping four of these lines on blood biochemistry suggested that there were substantial phenotypic variations among our lines, especially line C57BL/6J-Chr1^HZ and donor strain C57BL/6J. Further gene ontology enrichment revealed that the differentially expressed genes among liver-expressed genes between C57BL/6J and C57BL/6J-Chr1^HZ were enriched in lipid metabolism biological processes. All these characteristics enable C1SLs to be a unique resource for identifying and fine mapping quantitative trait loci on mouse Chr 1, and carrying out systems genetics studies of complex traits.

Regulated Inflammation and Lipid Metabolism in Colon mRNA Expressions of Obese Germfree Mice Responding to Enterobacter cloacae B29 Combined with the High Fat Diet

Increased evidences have demonstrated that gut microbiota targeted diet intervention can alleviate obesity and related metabolic disorders. The underlying mechanism of interactions among diet, microbiota, and host still remains unclear. Enterobacter cloacae B29, an endotoxin-producing strain dominated in the gut of a morbidly obese volunteer (weight 174.8 kg, BMI 58.8 kg m^-2) was isolated and transplanted to germfree mice (inoculated 10¹⁰ cells of B29 per day for 1 week). Using deep mRNA sequencing technology, we compared different gene expression profiles in the colon samples of the germfree mice treated with/without B29 and/or high fat diet (HFD) for 16 weeks and identified 279 differential expressed genes in total, including up-regulated genes Apoa4 (fold change, 2.77), Ido1 (2.66), Cyp4a10 (7.01), and down-regulated genes Cyp2e1 (0.11), Cyp26b1 (0.34), Akr1b7 (0.42), Adipoq (0.36), Cyp1a1 (0.11), Apoa1 (0.44), Npc1l1 (0.37), Tff2 (0.13), Apoc1 (0.30), Ctla2a (0.34), Mttp (0.49), Lpl (0.48). Fifty-nine GO biological processes and five KEGG pathways, particularly the peroxisome proliferator-activated receptors signaling pathway, were significantly enriched in response to HFD+B29, which were mainly relevant to inflammation and the metabolism of lipid, lipoprotein, and sterols. These functional changes were consistent with the developed obesity, insulin-resistance, and aggravated inflammatory conditions of the HFD+B29 mice. This work provides insight into the gene expression changes in response to HFD+B29, helping to understand the mechanism of the interactions among HFD, B29 and the germfree mice.

Mergeomics: multidimensional data integration to identify pathogenic perturbations to biological systems

BACKGROUND

Complex diseases are characterized by multiple subtle perturbations to biological processes. New omics platforms can detect these perturbations, but translating the diverse molecular and statistical information into testable mechanistic hypotheses is challenging. Therefore, we set out to create a public tool that integrates these data across multiple datasets, platforms, study designs and species in order to detect the most promising targets for further mechanistic studies.

RESULTS

We developed Mergeomics, a computational pipeline consisting of independent modules that 1) leverage multi-omics association data to identify biological processes that are perturbed in disease, and 2) overlay the disease-associated processes onto molecular interaction networks to pinpoint hubs as potential key regulators. Unlike existing tools that are mostly dedicated to specific data type or settings, the Mergeomics pipeline accepts and integrates datasets across platforms, data types and species. We optimized and evaluated the performance of Mergeomics using simulation and multiple independent datasets, and benchmarked the results against alternative methods. We also demonstrate the versatility of Mergeomics in two case studies that include genome-wide, epigenome-wide and transcriptome-wide datasets from human and mouse studies of total cholesterol and fasting glucose. In both cases, the Mergeomics pipeline provided statistical and contextual evidence to prioritize further investigations in the wet lab. The software implementation of Mergeomics is freely available as a Bioconductor R package.

CONCLUSION

Mergeomics is a flexible and robust computational pipeline for multidimensional data integration. It outperforms existing tools, and is easily applicable to datasets from different studies, species and omics data types for the study of complex traits.

Underlying genetic factors of the VATER/VACTERL association with special emphasis on the "Renal" phenotype

The acronym VATER/VACTERL association (OMIM #192350) refers to the rare non-random co-occurrence of the following component features (CFs): vertebral defects (V), anorectal malformations (A), cardiac defects (C), tracheoesophageal fistula with or without esophageal atresia (TE), renal malformations (R), and limb defects (L). According to epidemiological studies, the majority of patients with VATER/VACTERL association present with a "Renal" phenotype comprising a large spectrum of congenital renal anomalies. This finding is supported by evidence linking all of the human disease genes for the VATER/VACTERL association identified to date, namely, FGF8, FOXF1, HOXD13, LPP, TRAP1, and ZIC3, with renal malformations. Here we review these genotype-phenotype correlations and suggest that the elucidation of the genetic causes of the VATER/VACTERL association will ultimately provide insights into the genetic causes of the complete spectrum of congenital renal anomalies per se.

Modeling a human hepatocellular carcinoma subset in mice through coexpression of met and point-mutant β-catenin

Hepatocellular cancer (HCC) remains a significant therapeutic challenge due to its poorly understood molecular basis. In the current study, we investigated two independent cohorts of 249 and 194 HCC cases for any combinatorial molecular aberrations. Specifically we assessed for simultaneous HMET expression or hMet activation and catenin β1 gene (CTNNB1) mutations to address any concomitant Met and Wnt signaling. To investigate cooperation in tumorigenesis, we coexpressed hMet and β-catenin point mutants (S33Y or S45Y) in hepatocytes using sleeping beauty transposon/transposase and hydrodynamic tail vein injection and characterized tumors for growth, signaling, gene signatures, and similarity to human HCC. Missense mutations in exon 3 of CTNNB1 were identified in subsets of HCC patients. Irrespective of amino acid affected, all exon 3 mutations induced similar changes in gene expression. Concomitant HMET overexpression or hMet activation and CTNNB1 mutations were evident in 9%-12.5% of HCCs. Coexpression of hMet and mutant-β-catenin led to notable HCC in mice. Tumors showed active Wnt and hMet signaling with evidence of glutamine synthetase and cyclin D1 positivity and mitogen-activated protein kinase/extracellular signal-regulated kinase, AKT/Ras/mammalian target of rapamycin activation. Introduction of dominant-negative T-cell factor 4 prevented tumorigenesis. The gene expression of mouse tumors in hMet-mutant β-catenin showed high correlation, with subsets of human HCC displaying concomitant hMet activation signature and CTNNB1 mutations.

CONCLUSION

We have identified cooperation of hMet and β-catenin activation in a subset of HCC patients and modeled this human disease in mice with a significant transcriptomic intersection; this model will provide novel insight into the biology of this tumor and allow us to evaluate novel therapies as a step toward precision medicine. (Hepatology 2016;64:1587-1605).

Suppressive subtractive hybridization reveals different gene expression between high and low virulence strains of Cladosporium cladosporioides

Cladosporium cladosporioides is a ubiquitous fungus, causing infections in plants, humans, and animals. Suppression subtractive hybridization (SSH) and quantitative real-time PCR (qRT-PCR) were used in this study to identify differences in gene expression between two C. cladosporioides strains, the highly virulent Z20 strain and the lowly virulent Zt strain. A total of 61 unigenes from the forward library and 42 from the reverse library were identified. Gene ontology (GO) analysis showed that these genes were involved in various biological processes, cellular components and molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the unigenes in the forward library corresponded to 5 different pathways and the reverse library unigenes were involved in 3 different pathways. The qRT-PCR results indicated that expressions of APL1, GUD1, CSE1, SPBC3E7.04c and MFS were significantly different between Z20 and Zt strains, while genes encoding the senescence-associated proteins, pse1, nup107, mip1, pex2, icl1 and α/β hydrolase exhibited no significant differences between the two strains. In addition, we found that 5 unigenes encoding mip1, chk1, icl1, α/β hydrolase and β-glucosidase may be associated with pathogenicity. One unigene (MFS) may be related to the resistance to 14 α-demethylase inhibitor fungicides, and 5 unigenes (PEX2, NUP107, PSE1, APL1, and SPBC3E7.04c) may be related to either low spore yield or earlier aging of the Zt strain. Our study may help better understand the molecular mechanism of C. cladosporioides infection, and therefore improve the treatment and prevention of C. cladosporioides induced diseases.

Relevance of mouse models of cardiac fibrosis and hypertrophy in cardiac research

Heart disease causing cardiac cell death due to ischemia-reperfusion injury is a major cause of morbidity and mortality in the United States. Coronary heart disease and cardiomyopathies are the major cause for congestive heart failure, and thrombosis of the coronary arteries is the most common cause of myocardial infarction. Cardiac injury is followed by post-injury cardiac remodeling or fibrosis. Cardiac fibrosis is characterized by net accumulation of extracellular matrix proteins in the cardiac interstitium and results in both systolic and diastolic dysfunctions. It has been suggested by both experimental and clinical evidence that fibrotic changes in the heart are reversible. Hence, it is vital to understand the mechanism involved in the initiation, progression, and resolution of cardiac fibrosis to design anti-fibrotic treatment modalities. Animal models are of great importance for cardiovascular research studies. With the developing research field, the choice of selecting an animal model for the proposed research study is crucial for its outcome and translational purpose. Compared to large animal models for cardiac research, the mouse model is preferred by many investigators because of genetic manipulations and easier handling. This critical review is focused to provide insight to young researchers about the various mouse models, advantages and disadvantages, and their use in research pertaining to cardiac fibrosis and hypertrophy.

DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants

The information about the genetic basis of human diseases lies at the heart of precision medicine and drug discovery. However, to realize its full potential to support these goals, several problems, such as fragmentation, heterogeneity, availability and different conceptualization of the data must be overcome. To provide the community with a resource free of these hurdles, we have developed DisGeNET (http://www.disgenet.org), one of the largest available collections of genes and variants involved in human diseases. DisGeNET integrates data from expert curated repositories, GWAS catalogues, animal models and the scientific literature. DisGeNET data are homogeneously annotated with controlled vocabularies and community-driven ontologies. Additionally, several original metrics are provided to assist the prioritization of genotype-phenotype relationships. The information is accessible through a web interface, a Cytoscape App, an RDF SPARQL endpoint, scripts in several programming languages and an R package. DisGeNET is a versatile platform that can be used for different research purposes including the investigation of the molecular underpinnings of specific human diseases and their comorbidities, the analysis of the properties of disease genes, the generation of hypothesis on drug therapeutic action and drug adverse effects, the validation of computationally predicted disease genes and the evaluation of text-mining methods performance.

Cross-Species Analysis of Gene Expression and Function in Prefrontal Cortex, Hippocampus and Striatum

Background

Mouse has been extensively used as a tool for investigating the onset and development of human neurological disorders. As a first step to construct a transgenic mouse model of human brain lesions, it is of fundamental importance to clarify the similarity and divergence of genetic background between non-diseased human and mouse brain tissues.

Methods

We systematically compared, based on large scale integrated microarray data, the transcriptomes of three anatomically distinct brain regions; prefrontal cortex (PFC), hippocampus (HIP) and striatum (STR), across human and mouse. The widely used DAVID web server was used to decipher the biological functions of the highly expressed genes that were identified using a previously reported approach. Venn analysis was used to depict the overlapping ratios of the notably enriched biological process (BP) terms (one-tailed Fisher’s exact test and Benjamini correction; adjusted p < 0.01) between two brain tissues. GOSemSim, an R package, was selected to perform GO semantic similarity analysis. Next, we adjusted signal intensities of orthologous genes by the total signals in all samples within species, and used one minus Pearson’s correlation coefficient to assess the expression distance. Hierarchical clustering and principal component analysis (PCA) were selected for expression pattern analysis. Lineage specific expressed orthologous genes were identified by comparison of the most extreme sub-datasets across species and further verified using reverse transcription PCR (RT-PCR) and quantitative real-time PCR (qRT-PCR).

Results

We found that the number of the significantly enriched BP terms of the highly expressed genes in human brain regions is larger than that in mouse corresponding brain regions. The mainly involved BP terms in human brain tissues associated with protein-membrane targeting and selenium metabolism are species-specific. The overlapping ratios of all the significantly enriched BP terms between any two brain tissues across species are lower than that within species, but the pairwise semantic similarities are very high between any two brain tissues from either human or mouse. Hierarchical clustering analysis shows the biological functions of the highly expressed genes in brain tissues are more consistent within species than interspecies; whereas it shows the expression patterns of orthologous genes are evidently conserved between human and mouse equivalent brain tissues. In addition, we identified four orthologous genes (COX5B, WIF1, SLC4A10 and PLA2G7) that are species-specific, which have been widely studied and confirmed to be closely linked with neuro- physiological and pathological functions.

Conclusion

Our study highlights the similarities and divergences in gene function and expression between human and mouse corresponding brain regions, including PFC, HIP and STR.

Genome-wide compendium and functional assessment of in vivo heart enhancers

Whole-genome sequencing is identifying growing numbers of non-coding variants in human disease studies, but the lack of accurate functional annotations prevents their interpretation. We describe the genome-wide landscape of distant-acting enhancers active in the developing and adult human heart, an organ whose impairment is a predominant cause of mortality and morbidity. Using integrative analysis of >35 epigenomic data sets from mouse and human pre- and postnatal hearts we created a comprehensive reference of >80,000 putative human heart enhancers. To illustrate the importance of enhancers in the regulation of genes involved in heart disease, we deleted the mouse orthologs of two human enhancers near cardiac myosin genes. In both cases, we observe in vivo expression changes and cardiac phenotypes consistent with human heart disease. Our study provides a comprehensive catalogue of human heart enhancers for use in clinical whole-genome sequencing studies and highlights the importance of enhancers for cardiac function.

Identification of non-coding variants has outstripped our ability to annotate and interpret them. Dickel et al. present a compendium of over 80,000 putative human heart enhancers and demonstrate that two conserved enhancers are required for proper cardiac function in mice.

Deep sequencing of 10,000 human genomes

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

High-throughput discovery of novel developmental phenotypes

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.

A potential founder variant in CARMIL2/RLTPR in three Norwegian families with warts, molluscum contagiosum, and T-cell dysfunction

Background

Four patients from three Norwegian families presented with a common skin phenotype of warts, molluscum contagiosum, and dermatitis since early childhood, and various other immunological features. Warts are a common manifestation of human papilloma virus (HPV), but when they are overwhelming, disseminated and/or persistent, and presenting together with other immunological features, a primary immunodeficiency disease (PIDD) may be suspected.

Methods and results

The four patients were exome sequenced as part of a larger study for detecting genetic causes of primary immunodeficiencies. No disease‐causing variants were identified in known primary immunodeficiency genes or in other disease‐related OMIM genes. However, the same homozygous missense variant in CARMIL2 (also known as RLTPR) was identified in all four patients. In each family, the variant was located within a narrow region of homozygosity, representing a potential region of autozygosity. CARMIL2 is a protein of undetermined function. A role in T‐cell activation has been suggested and the mouse protein homolog (Rltpr) is essential for costimulation of T‐cell activation via CD28, and for the development of regulatory T cells. Immunophenotyping demonstrated reduced regulatory, CD4+ memory, and CD4+ follicular T cells in all four patients. In addition, they all seem to have a deficiency in IFNγ ‐synthesis in CD4+ T cells and NK cells.

Conclusions

We report a novel primary immunodeficiency, and a differential molecular diagnosis to CXCR4‐,DOCK8‐,GATA2‐,MAGT1‐,MCM4‐,STK4‐,RHOH‐,TMC6‐, and TMC8‐related diseases. The specific variant may represent a Norwegian founder variant segregating on a population‐specific haplotype.

De novo large rare copy-number variations contribute to conotruncal heart disease in Chinese patients

Conotruncal heart anomalies (CTDs) are particularly prevalent congenital heart diseases (CHD) in Hong Kong. We surveyed large (>500 kb), rare (<1% frequency in controls) copy-number variations (CNVs) in Chinese patients with CTDs to identify potentially disease-causing variations. Adults who tested negative for 22q11.2 deletions were recruited from the adult CHD clinic in Hong Kong. Using a stringent calling criteria, high-confidence CNV calls were obtained, and a large control set comprising 3,987 Caucasian and 1,945 Singapore Chinese subjects was used to identify rare CNVs. Ten large rare CNVs were identified, and 3 in 108 individuals were confirmed to harbour de novo CNVs. All three patients were syndromic with a more complex phenotype, and each of these CNVs overlapped regions likely to be important in CHD. One was a 611 kb deletion at 17p13.3, telomeric to the Miller-Dieker syndrome (MDS) critical region, overlapping the NXN gene. Another was a 5 Mb deletion at 13q33.3, within a previously described critical region for CHD. A third CNV, previously unreported, was a large duplication at 2q22.3 overlapping the ZEB2 gene. The commonly reported 1q21.1 recurrent duplication was not observed in this Chinese cohort. We provide detailed phenotypic and genotypic descriptions of large rare genic CNVs that may represent CHD loci in the East Asian population. Larger samples of Chinese origin will be required to determine whether the genome-wide distribution differs from that found in predominantly European CHD cohorts.

GenToS: Use of Orthologous Gene Information to Prioritize Signals from Human GWAS

Genome-wide association studies (GWAS) evaluate associations between genetic variants and a trait or disease of interest free of prior biological hypotheses. GWAS require stringent correction for multiple testing, with genome-wide significance typically defined as association p-value <5*10-8. This study presents a new tool that uses external information about genes to prioritize SNP associations (GenToS). For a given list of candidate genes, GenToS calculates an appropriate statistical significance threshold and then searches for trait-associated variants in summary statistics from human GWAS. It thereby allows for identifying trait-associated genetic variants that do not meet genome-wide significance. The program additionally tests for enrichment of significant candidate gene associations in the human GWAS data compared to the number expected by chance. As proof of principle, this report used external information from a comprehensive resource of genetically manipulated and systematically phenotyped mice. Based on selected murine phenotypes for which human GWAS data for corresponding traits were publicly available, several candidate gene input lists were derived. Using GenToS for the investigation of candidate genes underlying murine skeletal phenotypes in data from a large human discovery GWAS meta-analysis of bone mineral density resulted in the identification of significantly associated variants in 29 genes. Index variants in 28 of these loci were subsequently replicated in an independent GWAS replication step, highlighting that they are true positive associations. One signal, COL11A1, has not been discovered through GWAS so far and represents a novel human candidate gene for altered bone mineral density. The number of observed genes that contained significant SNP associations in human GWAS based on murine candidate gene input lists was much greater than the number expected by chance across several complex human traits (enrichment p-value as low as 10-10). GenToS can be used with any candidate gene list, any GWAS summary file, runs on a desktop computer and is freely available.

RNAseq analysis of heart tissue from mice treated with atenolol and isoproterenol reveals a reciprocal transcriptional response

Background

The transcriptional response to many widely used drugs and its modulation by genetic variability is poorly understood. Here we present an analysis of RNAseq profiles from heart tissue of 18 inbred mouse strains treated with the β-blocker atenolol (ATE) and the β-agonist isoproterenol (ISO).

Results

Differential expression analyses revealed a large set of genes responding to ISO (n = 1770 at FDR = 0.0001) and a comparatively small one responding to ATE (n = 23 at FDR = 0.0001). At a less stringent definition of differential expression, the transcriptional responses to these two antagonistic drugs are reciprocal for many genes, with an overall anti-correlation of r = −0.3. This trend is also observed at the level of most individual strains even though the power to detect differential expression is significantly reduced. The inversely expressed gene sets are enriched with genes annotated for heart-related functions. Modular analysis revealed gene sets that exhibit coherent transcription profiles across some strains and/or treatments. Correlations between these modules and a broad spectrum of cardiovascular traits are stronger than expected by chance. This provides evidence for the overall importance of transcriptional regulation for these organismal responses and explicits links between co-expressed genes and the traits they are associated with. Gene set enrichment analysis of differentially expressed groups of genes pointed to pathways related to heart development and functionality.

Conclusions

Our study provides new insights into the transcriptional response of the heart to perturbations of the β-adrenergic system, implicating several new genes that had not been associated to this system previously.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3059-6) contains supplementary material, which is available to authorized users.

Murine Plasma N-Glycosylation Traits Associated with Sex and Strain

Glycosylation is an abundant and important protein modification with large influence on the properties and interactions of glycoconjugates. Human plasma N-glycosylation has been the subject of frequent investigation, revealing strong associations with physiological and pathological conditions. Less well-characterized is the plasma N-glycosylation of the mouse, the most commonly used animal model for studying human diseases, particularly with regard to differences between strains and sexes. For this reason, we used MALDI-TOF(/TOF)-MS(/MS) assisted by linkage-specific derivatization of the sialic acids to comparatively analyze the plasma N-glycosylation of both male and female mice originating from BALB/c, CD57BL/6, CD-1, and Swiss Webster strains. The combined use of this analytical method and the recently developed data processing software named MassyTools allowed the relative quantification of the N-glycan species within plasma, the distinction between α2,3- and α2,6-linked N-glycolylneuraminic acids (due to respective lactonization and ethyl esterification), the detection of sialic acid O-acetylation, as well as the characterization of branching sialylation (Neu5Gcα2,3-Hex-[Neu5Gcα2,6-]HexNAc). When analyzing the glycosylation according to mouse sex, we found that female mice present a considerably higher degree of core fucosylation (2-4-fold depending on the strain), galactosylation, α2,6-linked sialylation, and larger high-mannose type glycan species compared with their male counterparts. Male mice, on the contrary, showed on average higher α2,3-linked sialylation, branching sialylation, and putative bisection. These differences together with sialic acid acetylation proved to be strain-specific as well. Interestingly, the outbred strains CD-1 and Swiss Webster displayed considerably larger interindividual variation than inbred strains BALB/c and CD57BL/6, suggesting a strong hereditable component of the observed plasma N-glycome.

Identification and Characterization of the Dermal Panniculus Carnosus Muscle Stem Cells

The dermal Panniculus carnosus (PC) muscle is important for wound contraction in lower mammals and represents an interesting model of muscle regeneration due to its high cell turnover. The resident satellite cells (the bona fide muscle stem cells) remain poorly characterized. Here we analyzed PC satellite cells with regard to developmental origin and purported function. Lineage tracing shows that they originate in Myf5⁺, Pax3/Pax7⁺ cell populations. Skin and muscle wounding increased PC myofiber turnover, with the satellite cell progeny being involved in muscle regeneration but with no detectable contribution to the wound-bed myofibroblasts. Since hematopoietic stem cells fuse to PC myofibers in the absence of injury, we also studied the contribution of bone marrow-derived cells to the PC satellite cell compartment, demonstrating that cells of donor origin are capable of repopulating the PC muscle stem cell niche after irradiation and bone marrow transplantation but may not fully acquire the relevant myogenic commitment.

•

PC satellite cells originate from Myf5⁺, Pax3/Pax7⁺ cell lineages

•

Skin and muscle wounding increase PC myofiber turnover

•

Donor bone marrow cells repopulate the PC satellite niche after BMT

•

Dermis-derived myogenesis originates from the PC satellite cell population

Polygenic Causes of Congenital Diaphragmatic Hernia Produce Common Lung Pathologies

Congenital diaphragmatic hernia (CDH) is one of the most common and lethal congenital anomalies, and significant evidence is available in support of a genetic contribution to its etiology, including single-gene knockout mice associated with diaphragmatic defects, rare monogenetic disorders in humans, familial aggregation, and association of CDH with chromosomal abnormalities. Structural lung defects in the form of lung hypoplasia are almost invariably seen in patients with CDH and frequently in animal models of this condition. Better understanding of the mechanisms of pulmonary defects in CDH has the potential for creating targeted therapies, particularly in postnatal stages, when therapeutics can have maximum clinical impact on the surviving cohorts. Successful treatment of CDH is dependent on the integration of human genomic and genetic data with developmental expression profiling, mouse knockouts, and gene network and pathway modeling, which have generated a large number of candidate genes and pathways for follow-up studies. In particular, defective alveolarization appears to be a common and potentially actionable phenotype in both patients and animal models.

Training and evaluation corpora for the extraction of causal relationships encoded in biological expression language (BEL)

Success in extracting biological relationships is mainly dependent on the complexity of the task as well as the availability of high-quality training data. Here, we describe the new corpora in the systems biology modeling language BEL for training and testing biological relationship extraction systems that we prepared for the BioCreative V BEL track. BEL was designed to capture relationships not only between proteins or chemicals, but also complex events such as biological processes or disease states. A BEL nanopub is the smallest unit of information and represents a biological relationship with its provenance. In BEL relationships (called BEL statements), the entities are normalized to defined namespaces mainly derived from public repositories, such as sequence databases, MeSH or publicly available ontologies. In the BEL nanopubs, the BEL statements are associated with citation information and supportive evidence such as a text excerpt. To enable the training of extraction tools, we prepared BEL resources and made them available to the community. We selected a subset of these resources focusing on a reduced set of namespaces, namely, human and mouse genes, ChEBI chemicals, MeSH diseases and GO biological processes, as well as relationship types ‘increases’ and ‘decreases’. The published training corpus contains 11 000 BEL statements from over 6000 supportive text excerpts. For method evaluation, we selected and re-annotated two smaller subcorpora containing 100 text excerpts. For this re-annotation, the inter-annotator agreement was measured by the BEL track evaluation environment and resulted in a maximal F-score of 91.18% for full statement agreement. In addition, for a set of 100 BEL statements, we do not only provide the gold standard expert annotations, but also text excerpts pre-selected by two automated systems. Those text excerpts were evaluated and manually annotated as true or false supportive in the course of the BioCreative V BEL track task.

Database URL:http://wiki.openbel.org/display/BIOC/Datasets

Towards precision medicine-based therapies for glioblastoma: interrogating human disease genomics and mouse phenotypes

BACKGROUND

Glioblastoma (GBM) is the most common and aggressive brain tumors. It has poor prognosis even with optimal radio- and chemo-therapies. Since GBM is highly heterogeneous, drugs that target on specific molecular profiles of individual tumors may achieve maximized efficacy. Currently, the Cancer Genome Atlas (TCGA) projects have identified hundreds of GBM-associated genes. We develop a drug repositioning approach combining disease genomics and mouse phenotype data towards predicting targeted therapies for GBM.

METHODS

We first identified disease specific mouse phenotypes using the most recently discovered GBM genes. Then we systematically searched all FDA-approved drugs for candidates that share similar mouse phenotype profiles with GBM. We evaluated the ranks for approved and novel GBM drugs, and compared with an existing approach, which also use the mouse phenotype data but not the disease genomics data.

RESULTS

We achieved significantly higher ranks for the approved and novel GBM drugs than the earlier approach. For all positive examples of GBM drugs, we achieved a median rank of 9.2 45.6 of the top predictions have been demonstrated effective in inhibiting the growth of human GBM cells.

CONCLUSION

We developed a computational drug repositioning approach based on both genomic and phenotypic data. Our approach prioritized existing GBM drugs and outperformed a recent approach. Overall, our approach shows potential in discovering new targeted therapies for GBM.

Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat

The aim of this work is to better understand the genetic mechanisms determining two complex traits affecting porcine meat quality: intramuscular fat (IMF) content and its fatty acid (FA) composition. With this purpose, expression Genome-Wide Association Study (eGWAS) of 45 lipid-related genes associated with meat quality traits in swine muscle (Longissimus dorsi) of 114 Iberian × Landrace backcross animals was performed. The eGWAS identified 241 SNPs associated with 11 genes: ACSM5, CROT, FABP3, FOS, HIF1AN, IGF2, MGLL, NCOA1, PIK3R1, PLA2G12A and PPARA. Three expression Quantitative Trait Loci (eQTLs) for IGF2, ACSM5 and MGLL were identified, showing cis-acting effects, whereas 16 eQTLs had trans regulatory effects. A polymorphism in the ACSM5 promoter region associated with its expression was identified. In addition, strong candidate genes regulating ACSM5, FOS, PPARA, PIK3R1, PLA2G12A and HIF1AN gene expression were also seen. Notably, the analysis highlighted the NR3C1 transcription factor as a strong candidate gene involved in the regulation of the 45 genes analysed. Finally, the IGF2, MGLL, MC2R, ARHGAP6, and NR3C1 genes were identified as potential regulators co-localizing within QTLs for fatness and growth traits in the IBMAP population. The results obtained increase our knowledge in the functional regulatory mechanisms involved in these complex traits.