The mouse Gene Expression Database (GXD): 2014 update

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

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.

Imprinted Gene Expression and Function of the Dopa Decarboxylase Gene in the Developing Heart

Dopa decarboxylase (DDC) synthesizes serotonin in the developing mouse heart where it is encoded by Ddc_exon1a, a tissue-specific paternally expressed imprinted gene. Ddc_exon1a shares an imprinting control region (ICR) with the imprinted, maternally expressed (outside of the central nervous system) Grb10 gene on mouse chromosome 11, but little else is known about the tissue-specific imprinted expression of Ddc_exon1a. Fluorescent immunostaining localizes DDC to the developing myocardium in the pre-natal mouse heart, in a region susceptible to abnormal development and implicated in congenital heart defects in human. Ddc_exon1a and Grb10 are not co-expressed in heart nor in brain where Grb10 is also paternally expressed, despite sharing an ICR, indicating they are mechanistically linked by their shared ICR but not by Grb10 gene expression. Evidence from a Ddc_exon1a gene knockout mouse model suggests that it mediates the growth of the developing myocardium and a thinning of the myocardium is observed in a small number of mutant mice examined, with changes in gene expression detected by microarray analysis. Comparative studies in the human developing heart reveal a paternal expression bias with polymorphic imprinting patterns between individual human hearts at DDC_EXON1a, a finding consistent with other imprinted genes in human.

Utilizing Organoid and Air-Liquid Interface Models as a Screening Method in the Development of New Host Defense Peptides

Host defense peptides (HDPs), also known as antimicrobial peptides, are naturally occurring polypeptides (~12–50 residues) composed of cationic and hydrophobic amino acids that adopt an amphipathic conformation upon folding usually after contact with membranes. HDPs have a variety of biological activities including immunomodulatory, anti-inflammatory, anti-bacterial, and anti-biofilm functions. Although HDPs have the potential to address the global threat of antibiotic resistance and to treat immune and inflammatory disorders, they have yet to achieve this promise. Indeed, there are several challenges associated with bringing peptide-based drug candidates from the lab bench to clinical practice, including identifying appropriate indications, stability, toxicity, and cost. These challenges can be addressed in part by the development of innate defense regulator (IDR) peptides and peptidomimetics, which are synthetic derivatives of HDPs with similar or better efficacy, increased stability, and reduced toxicity and cost of the original HDP. However, one of the largest gaps between basic research and clinical application is the validity and translatability of conventional model systems, such as cell lines and animal models, for screening HDPs and their derivatives as potential drug therapies. Indeed, such translation has often relied on animal models, which have only limited validity. Here we discuss the recent development of human organoids for disease modeling and drug screening, assisted by the use of omics analyses. Organoids, developed from primary cells, cell lines, or human pluripotent stem cells, are three-dimensional, self-organizing structures that closely resemble their corresponding in vivo organs with regards to immune responses, tissue organization, and physiological properties; thus, organoids represent a reliable method for studying efficacy, formulation, toxicity and to some extent drug stability and pharmacodynamics. The use of patient-derived organoids enables the study of patient-specific efficacy, toxicogenomics and drug response predictions. We outline how organoids and omics data analysis can be leveraged to aid in the clinical translation of IDR peptides.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment

BACKGROUND

The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually.

RESULTS

Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis.

CONCLUSIONS

Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Duplication 2p16 is associated with perisylvian polymicrogyria

Polymicrogyria (PMG) is a heterogeneous brain malformation that may result from prenatal vascular disruption or infection, or from numerous genetic causes that still remain difficult to identify. We identified three unrelated patients with polymicrogyria and duplications of chromosome 2p, defined the smallest region of overlap, and performed gene pathway analysis using Cytoscape. The smallest region of overlap in all three children involved 2p16.1-p16.3. All three children have bilateral perisylvian polymicrogyria (BPP), intrauterine and postnatal growth deficiency, similar dysmorphic features, and poor feeding. Two of the three children had documented intellectual disability. Gene pathway analysis suggested a number of developmentally relevant genes and gene clusters that were over-represented in the critical region. We narrowed a rare locus for polymicrogyria to a region of 2p16.1-p16.3 that contains 33-34 genes, 23 of which are expressed in cerebral cortex during human fetal development. Using pathway analysis, we showed that several of the duplicated genes contribute to neurodevelopmental pathways including morphogen, cytokine, hormonal and growth factor signaling, regulation of cell cycle progression, cell morphogenesis, axonal guidance, and neuronal migration. These findings strengthen the evidence for a novel locus associated with polymicrogyria on 2p16.1-p16.3, and comprise the first step in defining the underlying genetic etiology.

High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies

We describe a 7-year-old male with high functioning autism spectrum disorder (ASD) and maternally-inherited rare missense variant of Synaptotagmin-like protein 4 (SYTL4) gene (Xq22.1; c.835C>T; p.Arg279Cys) and an unknown missense variant of Transmembrane protein 187 (TMEM187) gene (Xq28; c.708G>T; p. Gln236His). Multiple in-silico predictions described in our study indicate a potentially damaging status for both X-linked genes. Analysis of predicted atomic threading models of the mutant and the native SYTL4 proteins suggest a potential structural change induced by the R279C variant which eliminates the stabilizing Arg279-Asp60 salt bridge in the N-terminal half of the SYTL4, affecting the functionality of the protein's critical RAB-Binding Domain. In the European (Non-Finnish) population, the allele frequency for this variant is 0.00042. The SYTL4 gene is known to directly interact with several members of the RAB family of genes, such as, RAB27A, RAB27B, RAB8A, and RAB3A which are known autism spectrum disorder genes. The SYTL4 gene also directly interacts with three known autism genes: STX1A, SNAP25 and STXBP1. Through a literature-based analytical approach, we identified three of five (60%) autism-associated serum microRNAs (miRs) with high predictive power among the total of 298 mouse Sytl4 associated/predicted microRNA interactions. Five of 13 (38%) miRs were differentially expressed in serum from ASD individuals which were predicted to interact with the mouse equivalent Sytl4 gene. TMEM187 gene, like SYTL4, is a protein-coding gene that belongs to a group of genes which host microRNA genes in their introns or exons. The novel Q236H amino acid variant in the TMEM187 in our patient is near the terminal end region of the protein which is represented by multiple sequence alignments and hidden Markov models, preventing comparative structural analysis of the variant harboring region. Like SYTL4, the TMEM187 gene is expressed in the brain and interacts with four known ASD genes, namely, HCFC1; TMLHE; MECP2; and GPHN. TMM187 is in linkage with MECP2, which is a well-known determinant of brain structure and size and is a well-known autism gene. Other members of the TMEM gene family, TMEM132E and TMEM132D genes are associated with bipolar and panic disorders, respectively, while TMEM231 is a known syndromic autism gene. Together, TMEM187 and SYTL4 genes directly interact with recognized important ASD genes, and their mRNAs are found in extracellular vesicles in the nervous system and stimulate target cells to translate into active protein. Our evidence shows that both these genes should be considered as candidate genes for autism. Additional biological testing is warranted to further determine the pathogenicity of these gene variants in the causation of autism.

Common and distinct transcriptional signatures of mammalian embryonic lethality

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4-36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders.

Identification of TGFβ-related genes regulated in murine osteoarthritis and chondrocyte hypertrophy by comparison of multiple microarray datasets

OBJECTIVE

Osteoarthritis (OA) is a joint disease characterized by progressive degeneration of articular cartilage. Some features of OA, including chondrocyte hypertrophy and focal calcification of articular cartilage, resemble the endochondral ossification processes. Alterations in transforming growth factor β (TGFβ) signaling have been associated with OA as well as with chondrocyte hypertrophy. Our aim was to identify novel candidate genes implicated in chondrocyte hypertrophy during OA pathogenesis by determining which TGFβ-related genes are regulated during murine OA and endochondral ossification.

METHODS

A list of 580 TGFβ-related genes, including TGFβ signaling pathway components and TGFβ-target genes, was generated. Regulation of these TGFβ-related genes was assessed in a microarray of murine OA cartilage: 1, 2 and 6 weeks after destabilization of the medial meniscus (DMM). Subsequently, genes regulated in the DMM model were studied in two independent murine microarray datasets on endochondral ossification: the growth plate and transient embryonic cartilage (joint development).

RESULTS

A total of 106 TGFβ-related genes were differentially expressed in articular cartilage of DMM-operated mice compared to sham-control. From these genes, 43 were similarly regulated during chondrocyte hypertrophy in the growth plate or embryonic joint development. Among these 43 genes, 18 genes have already been associated with OA. The remaining 25 genes were considered as novel candidate genes involved in OA pathogenesis and endochondral ossification. In supplementary data of published human OA microarrays we found indications that 15 of the 25 novel genes are indeed regulated in articular cartilage of human OA patients.

CONCLUSION

By focusing on TGFβ-related genes during OA and chondrocyte hypertrophy in mice, we identified 18 known and 25 new candidate genes potentially implicated in phenotypical changes in chondrocytes leading to OA. We propose that 15 of these candidates warrant further investigation as therapeutic target for OA as they are also regulated in articular cartilage of OA patients.

Hematological findings in non-treated and gamma-irradiated mice deficient for MIC-1/GDF15

Several members of the TGF-beta family are known to effectively regulate the fate of hematopoietic progenitor cells in a complex and context-dependent manner. Growth differentiation factor-15 (GDF15) is a divergent member of the TGF-beta family. This stress-induced cytokine has been proposed to possess immunomodulatory functions and its high expression is often associated with progression of a variety of pathological conditions. GDF15 is also induced by chemotherapy and irradiation. Very few fundamental studies have been published regarding the effect of GDF15 in hematopoiesis. In this study, we analyzed the hematological status of untreated and gamma-irradiated mice deficient for GDF15 as a result of genetic knock-out (KO), in order to clarify the regulatory role of GDF15 in hematopoiesis. Significant differences between GDF15 KO mice and their pertinent WT controls were found in the parameters of blood monocyte numbers, blood platelet size, and distribution width, as well as in the values of bone marrow granulocyte/macrophage progenitor cells. Different tendencies of some hematological parameters in the GDF15 KO mice in normal conditions and those under exposure of the mice to ionizing radiation were registered. These findings are discussed in the context of the GDF15 gene function and its lack under conditions of radiation-induced damage.

MicroRNAs Induce a Permissive Chromatin Environment that Enables Neuronal Subtype-Specific Reprogramming of Adult Human Fibroblasts

Directed reprogramming of human fibroblasts into fully differentiated neurons requires massive changes in epigenetic and transcriptional states. Induction of a chromatin environment permissive for acquiring neuronal subtype identity is therefore a major barrier to fate conversion. Here we show that the brain-enriched miRNAs miR-9/9^∗ and miR-124 (miR-9/9^∗-124) trigger reconfiguration of chromatin accessibility, DNA methylation, and mRNA expression to induce a default neuronal state. miR-9/9^∗-124-induced neurons (miNs) are functionally excitable and uncommitted toward specific subtypes but possess open chromatin at neuronal subtype-specific loci, suggesting that such identity can be imparted by additional lineage-specific transcription factors. Consistently, we show that ISL1 and LHX3 selectively drive conversion to a highly homogeneous population of human spinal cord motor neurons. This study shows that modular synergism between miRNAs and neuronal subtype-specific transcription factors can drive lineage-specific neuronal reprogramming, providing a general platform for high-efficiency generation of distinct subtypes of human neurons.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

Emerging roles for RNA-binding proteins as effectors and regulators of cardiovascular disease

The cardiovascular system comprises multiple cell types that possess the capacity to modulate their phenotype in response to acute or chronic injury. Transcriptional and post-transcriptional mechanisms play a key role in the regulation of remodelling and regenerative responses to damaged cardiovascular tissues. Simultaneously, insufficient regulation of cellular phenotype is tightly coupled with the persistence and exacerbation of cardiovascular disease. Recently, RNA-binding proteins such as Quaking, HuR, Muscleblind, and SRSF1 have emerged as pivotal regulators of these functional adaptations in the cardiovascular system by guiding a wide-ranging number of post-transcriptional events that dramatically impact RNA fate, including alternative splicing, stability, localization and translation. Moreover, homozygous disruption of RNA-binding protein genes is commonly associated with cardiac- and/or vascular complications. Here, we summarize the current knowledge on the versatile role of RNA-binding proteins in regulating the transcriptome during phenotype switching in cardiovascular health and disease. We also detail existing and potential DNA- and RNA-based therapeutic approaches that could impact the treatment of cardiovascular disease in the future.

Susceptibility of different mouse strains to peri-implantitis

BACKGROUND AND OBJECTIVE

Peri-implantitis (PI) is an inflammatory condition that affects the tissues surrounding dental implants. Although the pathogenesis of PI is not fully understood, evidence suggests that the etiology is multifactorial and may include a genetic component. The aim of this study was to investigate the role of genetics in the development of peri-implantitis.

MATERIAL AND METHODS

Four-week-old C57BL/6J, C3H/HeJ and A/J male mice had their left maxillary molars extracted. Implants were placed in the healed extraction sockets. Upon osseointegration, ligatures were placed around the implant head for 1 or 4 weeks to induce PI. Micro-computed tomography scanning was used to measure volumetric bone loss. Histological analyses were also performed to evaluate collagen organization and the presence of neutrophils and osteoclasts.

RESULTS

Radiographically, comparing the ligature-treated mice, C57BL/6J displayed the greatest amount of bone loss, followed by C3H/HeJ and A/J mice at 1 and 4 weeks. Histologically, at 1 week, C57BL/6J mice presented with the highest numbers of neutrophils and osteoclasts. At 4 weeks, C57BL/6J mice presented with the most active bone remodeling compared with the other two strains.

CONCLUSION

There were significant differences in the severity of peri-implantitis among the different mouse strains, suggesting that the genetic framework can affect implant survival and success. Future work is needed to dissect the genetic contribution to the development of peri-implantitis.

The 24th annual Nucleic Acids Research database issue: a look back and upcoming changes

This year's Database Issue of Nucleic Acids Research contains 152 papers that include descriptions of 54 new databases and update papers on 98 databases, of which 16 have not been previously featured in NAR As always, these databases cover a broad range of molecular biology subjects, including genome structure, gene expression and its regulation, proteins, protein domains, and protein-protein interactions. Following the recent trend, an increasing number of new and established databases deal with the issues of human health, from cancer-causing mutations to drugs and drug targets. In accordance with this trend, three recently compiled databases that have been selected by NAR reviewers and editors as 'breakthrough' contributions, denovo-db, the Monarch Initiative, and Open Targets, cover human de novo gene variants, disease-related phenotypes in model organisms, and a bioinformatics platform for therapeutic target identification and validation, respectively. We expect these databases to attract the attention of numerous researchers working in various areas of genetics and genomics. Looking back at the past 12 years, we present here the 'golden set' of databases that have consistently served as authoritative, comprehensive, and convenient data resources widely used by the entire community and offer some lessons on what makes a successful database. The Database Issue is freely available online at the https://academic.oup.com/nar web site. An updated version of the NAR Molecular Biology Database Collection is available at http://www.oxfordjournals.org/nar/database/a/.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.

Comparative studies of vertebrate iduronate 2-sulfatase (IDS) genes and proteins: evolution of A mammalian X-linked gene

IDS is responsible for the lysosomal degradation of heparan sulfate and dermatan sulfate and linked to an X-linked lysosomal storage disease, mucopolysaccharidosis 2 (MPS2), resulting in neurological damage and early death. Comparative IDS amino acid sequences and structures and IDS gene locations were examined using data from several vertebrate genome projects. Vertebrate IDS sequences shared 60–99% identities with each other. Human IDS showed 47% sequence identity with fruit fly (Drosophila melanogaster) IDS. Sequence alignments, key amino acid residues, N-glycosylation sites and conserved predicted secondary and tertiary structures were also studied, including signal peptide, propeptide and active site residues. Mammalian IDS genes usually contained 9 coding exons. The human IDS gene promoter contained a large CpG island (CpG46) and 5 transcription factor binding sites, whereas the 3′-UTR region contained 5 miRNA target sites. These may contribute to IDS gene regulation of expression in the brain and other neural tissues of the body. An IDS pseudogene (IDSP1) was located proximally to the IDS gene on the X-chromosome in primate genomes. Phylogenetic analyses examined the relationships and potential evolutionary origins of the vertebrate IDS gene. These suggested that IDS has originated in an invertebrate ancestral genome and retained throughout vertebrate evolution and conserved on marsupial and eutherian X-chromosomes, with the exception of rat Ids on chromosome 8.

Role of renal urothelium in the development and progression of kidney disease

The clinical and financial impact of chronic kidney disease (CKD) is significant, while its progression and prognosis is variable and often poor. Studies using the megabladder (mgb ^-/- ) model of CKD show that renal urothelium plays a key role in modulating early injury responses following the development of congenital obstruction. The aim of this review is to examine the role that urothelium has in normal urinary tract development and pathogenesis. We discuss normal morphology of renal urothelium and then examine the role that uroplakins (Upks) play in its development. Histologic, biochemical, and molecular characterization of Upk1b ^RFP/RFP mice indicated Upk1b expression is essential for normal urinary tract development, apical plaque/asymmetric membrane unit (AUM) formation, and differentiation and functional integrity of the renal urothelium. Our studies provide the first evidence that Upk1b is directly associated with the development of congenital anomalies of the urinary tract (CAKUT), spontaneous age-dependent hydronephrosis, and dysplastic urothelia. These observations demonstrate the importance of proper urothelial differentiation in normal development and pathogenesis of the urinary tract and provide a unique working model to test the hypothesis that the complex etiology associated with CKD is dependent upon predetermined genetic susceptibilities that establish pathogenic thresholds for disease initiation and progression.

eMouseAtlas: An atlas-based resource for understanding mammalian embryogenesis

The eMouseAtlas resource is an online database of 3D digital models of mouse development, an ontology of mouse embryo anatomy and a gene-expression database with about 30K spatially mapped gene-expression patterns. It is closely linked with the MGI/GXD database at the Jackson Laboratory and holds links to almost all available image-based gene-expression data for the mouse embryo. In this resource article we describe the novel web-based tools we have developed for 3D visualisation of embryo anatomy and gene expression. We show how mapping of gene expression data onto spatial models delivers a framework for capturing gene expression that enhances our understanding of development, and we review the exploratory tools utilised by the EMAGE gene expression database as a means of defining co-expression of in situ hybridisation, immunohistochemistry, and lacZ-omic expression patterns. We report on recent developments of the eHistology atlas and our use of web-services to support embedding of the online 'The Atlas of Mouse Development' in the context of other resources such as the DMDD mouse phenotype database. In addition, we discuss new developments including a cellular-resolution placental atlas, third-party atlas models, clonal analysis data and a new interactive eLearning resource for developmental processes.

Transcriptomic dynamics of breast cancer progression in the MMTV-PyMT mouse model

BACKGROUND

Malignant breast cancer with complex molecular mechanisms of progression and metastasis remains a leading cause of death in women. To improve diagnosis and drug development, it is critical to identify panels of genes and molecular pathways involved in tumor progression and malignant transition. Using the PyMT mouse, a genetically engineered mouse model that has been widely used to study human breast cancer, we profiled and analyzed gene expression from four distinct stages of tumor progression (hyperplasia, adenoma/MIN, early carcinoma and late carcinoma) during which malignant transition occurs.

RESULTS

We found remarkable expression similarity among the four stages, meaning genes altered in the later stages showed trace in the beginning of tumor progression. We identified a large number of differentially expressed genes in PyMT samples of all stages compared with normal mammary glands, enriched in cancer-related pathways. Using co-expression networks, we found panels of genes as signature modules with some hub genes that predict metastatic risk. Time-course analysis revealed genes with expression transition when shifting to malignant stages. These may provide additional insight into the molecular mechanisms beyond pathways.

CONCLUSIONS

Thus, in this study, our various analyses with the PyMT mouse model shed new light on transcriptomic dynamics during breast cancer malignant progression.

ARID3A is required for mammalian placenta development

Previous studies in the mouse indicated that ARID3A plays a critical role in the first cell fate decision required for generation of trophectoderm (TE). Here, we demonstrate that ARID3A is widely expressed during mouse and human placentation and essential for early embryonic viability. ARID3A localizes to trophoblast giant cells and other trophoblast-derived cell subtypes in the junctional and labyrinth zones of the placenta. Conventional Arid3a knockout embryos suffer restricted intrauterine growth with severe defects in placental structural organization. Arid3a null placentas show aberrant expression of subtype-specific markers as well as significant alteration in cytokines, chemokines and inflammatory response-related genes, including previously established markers of human placentation disorders. BMP4-mediated induction of trophoblast stem (TS)-like cells from human induced pluripotent stem cells results in ARID3A up-regulation and cytoplasmic to nuclear translocation. Overexpression of ARID3A in BMP4-mediated TS-like cells up-regulates TE markers, whereas pluripotency markers are down-regulated. Our results reveal an essential, conserved function for ARID3A in mammalian placental development through regulation of both intrinsic and extrinsic developmental programs.

Retinoic acid receptor regulation of epimorphic and homeostatic regeneration in the axolotl

Salamanders are capable of regenerating amputated limbs by generating a mass of lineage-restricted cells called a blastema. Blastemas only generate structures distal to their origin unless treated with retinoic acid (RA), which results in proximodistal (PD) limb duplications. Little is known about the transcriptional network that regulates PD duplication. In this study, we target specific retinoic acid receptors (RARs) to either PD duplicate (RA treatment or RARγ agonist) or truncate (RARβ antagonist) regenerating limbs. RARE-EGFP reporter axolotls showed divergent reporter activity in limbs undergoing PD duplication versus truncation, suggesting differences in patterning and skeletal regeneration. Transcriptomics identified expression patterns that explain PD duplication, including upregulation of proximal homeobox gene expression and silencing of distal-associated genes, whereas limb truncation was associated with disrupted skeletal differentiation. RARβ antagonism in uninjured limbs induced a loss of skeletal integrity leading to long bone regression and loss of skeletal turnover. Overall, mechanisms were identified that regulate the multifaceted roles of RARs in the salamander limb including regulation of skeletal patterning during epimorphic regeneration, skeletal tissue differentiation during regeneration, and homeostatic regeneration of intact limbs.

Search for Master Regulators in Walking Cancer Pathways

In this chapter, we present an approach that allows a causal analysis of multiple "-omics" data with the help of an "upstream analysis" strategy. The goal of this approach is to identify master regulators in gene regulatory networks as potential drug targets for a pathological process. The data analysis strategy includes a state-of-the-art promoter analysis for potential transcription factor (TF)-binding sites using the TRANSFAC^® database combined with an analysis of the upstream signal transduction pathways that control the activity of these TFs. When applied to genes that are associated with a switch to a pathological process, the approach identifies potential key molecules (master regulators) that may exert major control over and maintenance of transient stability of the pathological state. We demonstrate this approach on examples of analysis of multi-omics data sets that contain transcriptomics and epigenomics data in cancer. The results of this analysis helped us to better understand the molecular mechanisms of cancer development and cancer drug resistance. Such an approach promises to be very effective for rapid and accurate identification of cancer drug targets with true potential. The upstream analysis approach is implemented as an automatic workflow in the geneXplain platform ( www.genexplain.com ) using the open-source BioUML framework ( www.biouml.org ).

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.

GeneAnalytics: An Integrative Gene Set Analysis Tool for Next Generation Sequencing, RNAseq and Microarray Data

Postgenomics data are produced in large volumes by life sciences and clinical applications of novel omics diagnostics and therapeutics for precision medicine. To move from “data-to-knowledge-to-innovation,” a crucial missing step in the current era is, however, our limited understanding of biological and clinical contexts associated with data. Prominent among the emerging remedies to this challenge are the gene set enrichment tools. This study reports on GeneAnalytics™ (geneanalytics.genecards.org), a comprehensive and easy-to-apply gene set analysis tool for rapid contextualization of expression patterns and functional signatures embedded in the postgenomics Big Data domains, such as Next Generation Sequencing (NGS), RNAseq, and microarray experiments. GeneAnalytics' differentiating features include in-depth evidence-based scoring algorithms, an intuitive user interface and proprietary unified data. GeneAnalytics employs the LifeMap Science's GeneCards suite, including the GeneCards®—the human gene database; the MalaCards—the human diseases database; and the PathCards—the biological pathways database. Expression-based analysis in GeneAnalytics relies on the LifeMap Discovery®—the embryonic development and stem cells database, which includes manually curated expression data for normal and diseased tissues, enabling advanced matching algorithm for gene–tissue association. This assists in evaluating differentiation protocols and discovering biomarkers for tissues and cells. Results are directly linked to gene, disease, or cell “cards” in the GeneCards suite. Future developments aim to enhance the GeneAnalytics algorithm as well as visualizations, employing varied graphical display items. Such attributes make GeneAnalytics a broadly applicable postgenomics data analyses and interpretation tool for translation of data to knowledge-based innovation in various Big Data fields such as precision medicine, ecogenomics, nutrigenomics, pharmacogenomics, vaccinomics, and others yet to emerge on the postgenomics horizon.

Genetics of Combined Pituitary Hormone Deficiency: Roadmap into the Genome Era

The genetic basis for combined pituitary hormone deficiency (CPHD) is complex, involving 30 genes in a variety of syndromic and nonsyndromic presentations. Molecular diagnosis of this disorder is valuable for predicting disease progression, avoiding unnecessary surgery, and family planning. We expect that the application of high throughput sequencing will uncover additional contributing genes and eventually become a valuable tool for molecular diagnosis. For example, in the last 3 years, six new genes have been implicated in CPHD using whole-exome sequencing. In this review, we present a historical perspective on gene discovery for CPHD and predict approaches that may facilitate future gene identification projects conducted by clinicians and basic scientists. Guidelines for systematic reporting of genetic variants and assigning causality are emerging. We apply these guidelines retrospectively to reports of the genetic basis of CPHD and summarize modes of inheritance and penetrance for each of the known genes. In recent years, there have been great improvements in databases of genetic information for diverse populations. Some issues remain that make molecular diagnosis challenging in some cases. These include the inherent genetic complexity of this disorder, technical challenges like uneven coverage, differing results from variant calling and interpretation pipelines, the number of tolerated genetic alterations, and imperfect methods for predicting pathogenicity. We discuss approaches for future research in the genetics of CPHD.

The mouse Gene Expression Database (GXD): 2017 update

The Gene Expression Database (GXD; www.informatics.jax.org/expression.shtml) is an extensive and well-curated community resource of mouse developmental expression information. Through curation of the scientific literature and by collaborations with large-scale expression projects, GXD collects and integrates data from RNA in situ hybridization, immunohistochemistry, RT-PCR, northern blot and western blot experiments. Expression data from both wild-type and mutant mice are included. The expression data are combined with genetic and phenotypic data in Mouse Genome Informatics (MGI) and made readily accessible to many types of database searches. At present, GXD includes over 1.5 million expression results and more than 300 000 images, all annotated with detailed and standardized metadata. Since our last report in 2014, we have added a large amount of data, we have enhanced data and database infrastructure, and we have implemented many new search and display features. Interface enhancements include: a new Mouse Developmental Anatomy Browser; interactive tissue-by-developmental stage and tissue-by-gene matrix views; capabilities to filter and sort expression data summaries; a batch search utility; gene-based expression overviews; and links to expression data from other species.

Single-Cell Resolution of Temporal Gene Expression during Heart Development

Activation of complex molecular programs in specific cell lineages governs mammalian heart development, from a primordial linear tube to a four-chamber organ. To characterize lineage-specific, spatiotemporal developmental programs, we performed single-cell RNA sequencing of >1,200 murine cells isolated at seven time points spanning embryonic day 9.5 (primordial heart tube) to postnatal day 21 (mature heart). Using unbiased transcriptional data, we classified cardiomyocytes, endothelial cells, and fibroblast-enriched cells, thus identifying markers for temporal and chamber-specific developmental programs. By harnessing these datasets, we defined developmental ages of human and mouse pluripotent stem-cell-derived cardiomyocytes and characterized lineage-specific maturation defects in hearts of mice with heterozygous mutations in Nkx2.5 that cause human heart malformations. This spatiotemporal transcriptome analysis of heart development reveals lineage-specific gene programs underlying normal cardiac development and congenital heart disease.

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.

Multi-omics "upstream analysis" of regulatory genomic regions helps identifying targets against methotrexate resistance of colon cancer

•

Upstream analysis strategy for multi-omics data is proposed.

•

Drug targets are predicted by search for TFBS and analysis of signaling network.

•

Methotrexate resistance data include transcriptomics, proteomics and epigenomics.

•

Predicted targets are: TGFalpha, IGFBP7, alpha9-integrin.

•

Predicted drugs are: zardaverine, divalproex and human metabolite nicotinamide N-oxide.

We present an “upstream analysis” strategy for causal analysis of multiple “-omics” data. It analyzes promoters using the TRANSFAC database, combines it with an analysis of the upstream signal transduction pathways and identifies master regulators as potential drug targets for a pathological process. We applied this approach to a complex multi-omics data set that contains transcriptomics, proteomics and epigenomics data. We identified the following potential drug targets against induced resistance of cancer cells towards chemotherapy by methotrexate (MTX): TGFalpha, IGFBP7, alpha9-integrin, and the following chemical compounds: zardaverine and divalproex as well as human metabolites such as nicotinamide N-oxide.

Transcriptome and proteome characterization of surface ectoderm cells differentiated from human iPSCs

Surface ectoderm (SE) cells give rise to structures including the epidermis and ectodermal associated appendages such as hair, eye, and the mammary gland. In this study, we validate a protocol that utilizes BMP4 and the γ-secretase inhibitor DAPT to induce SE differentiation from human induced pluripotent stem cells (hiPSCs). hiPSC-differentiated SE cells expressed markers suggesting their commitment to the SE lineage. Computational analyses using integrated quantitative transcriptomic and proteomic profiling reveal that TGFβ superfamily signaling pathways are preferentially activated in SE cells compared with hiPSCs. SE differentiation can be enhanced by selectively blocking TGFβ-RI signaling. We also show that SE cells and neural ectoderm cells possess distinct gene expression patterns and signaling networks as indicated by functional Ingenuity Pathway Analysis. Our findings advance current understanding of early human SE cell development and pave the way for modeling of SE-derived tissue development, studying disease pathogenesis, and development of regenerative medicine approaches.

Exome sequencing discloses KALRN homozygous variant as likely cause of intellectual disability and short stature in a consanguineous pedigree

BACKGROUND

The recent availability of whole-exome sequencing has opened new possibilities for the evaluation of individuals with genetically undiagnosed intellectual disability.

RESULTS

We report two affected siblings, offspring of first-cousin parents, with intellectual disability, hypotonia, short stature, growth hormone deficiency, and delayed bone age. All members of the nuclear family were genotyped, and exome sequencing was performed in one of the affected individuals. We used an in-house algorithm (CATCH v1.1) that combines homozygosity mapping with exome sequencing results and provides a list of candidate variants. One identified novel homozygous missense variant in KALRN (NM_003947.4:c.3644C>A: p.(Thr1215Lys)) was predicted to be pathogenic by all pathogenicity prediction software used (SIFT, PolyPhen, Mutation Taster). KALRN encodes the protein kalirin, which is a GTP-exchange factor protein with a reported role in cytoskeletal remodeling and dendritic spine formation in neurons. It is known that mice with ablation of Kalrn exhibit age-dependent functional deficits and behavioral phenotypes.

CONCLUSION

Exome sequencing provided initial evidence linking KALRN to monogenic intellectual disability in man, and we propose that KALRN is the causative gene for the autosomal recessive phenotype in this family.

Using Whole Exome Sequencing to Identify Candidate Genes With Rare Variants In Nonsyndromic Cleft Lip and Palate

Studies suggest that nonsyndromic cleft lip and palate (NSCLP) is polygenic with variable penetrance, presenting a challenge in identifying all causal genetic variants. Despite relatively high prevalence of NSCLP among Amerindian populations, no large whole exome sequencing (WES) studies have been completed in this population. Our goal was to identify candidate genes with rare genetic variants for NSCLP in a Honduran population using WES. WES was performed on two to four members of 27 multiplex Honduran families. Genetic variants with a minor allele frequency > 1% in reference databases were removed. Heterozygous variants consistent with dominant disease with incomplete penetrance were ascertained, and variants with predicted functional consequence were prioritized for analysis. Pedigree-specific P-values were calculated as the probability of all affected members in the pedigree being carriers, given that at least one is a carrier. Preliminary results identified 3,727 heterozygous rare variants; 1,282 were predicted to be functionally consequential. Twenty-three genes had variants of interest in ≥3 families, where some genes had different variants in each family, giving a total of 50 variants. Variant validation via Sanger sequencing of the families and unrelated unaffected controls excluded variants that were sequencing errors or common variants not in databases, leaving four genes with candidate variants in ≥3 families. Of these, candidate variants in two genes consistently segregate with NSCLP as a dominant variant with incomplete penetrance: ACSS2 and PHYH. Rare variants found at the same gene in all affected individuals in several families are likely to be directly related to NSCLP.

Defining the identity of mouse embryonic dermal fibroblasts

Embryonic dermal fibroblasts in the skin have the exceptional ability to initiate hair follicle morphogenesis and contribute to scarless wound healing. Activation of the Wnt signaling pathway is critical for dermal fibroblast fate selection and hair follicle induction. In humans, mutations in Wnt pathway components and target genes lead to congenital focal dermal hypoplasias with diminished hair. The gene expression signature of embryonic dermal fibroblasts during differentiation and its dependence on Wnt signaling is unknown. Here we applied Shannon entropy analysis to identify the gene expression signature of mouse embryonic dermal fibroblasts. We used available human DNase-seq and histone modification ChiP-seq data on various cell-types to demonstrate that genes in the fibroblast cell identity signature can be epigenetically repressed in other cell-types. We found a subset of the signature genes whose expression is dependent on Wnt/β-catenin activity in vivo. With our approach, we have defined and validated a statistically derived gene expression signature that may mediate dermal fibroblast identity and function in development and disease. genesis 54:415-430, 2016. © 2016 Wiley Periodicals, Inc.

Parallel evolution of male germline epigenetic poising and somatic development in animals

Changes in gene regulation frequently underlie changes in morphology during evolution, and differences in chromatin state have been linked with changes in anatomical structure and gene expression across evolutionary time. Here we assess the relationship between evolution of chromatin state in germ cells and evolution of gene regulatory programs governing somatic development. We examined the poised (H3K4me3/H3K27me3 bivalent) epigenetic state in male germ cells from five mammalian and one avian species. We find that core genes poised in germ cells from multiple amniote species are ancient regulators of morphogenesis that sit at the top of transcriptional hierarchies controlling somatic tissue development, whereas genes that gain poising in germ cells from individual species act downstream of core poised genes during development in a species-specific fashion. We propose that critical regulators of animal development gained an epigenetically privileged state in germ cells, manifested in amniotes by H3K4me3/H3K27me3 poising, early in metazoan evolution.

Representation of anatomy in online atlases and databases: a survey and collection of patterns for interface design

BACKGROUND

A large number of online atlases and databases have been developed to mange the rapidly growing amount of data describing embryogenesis. As these community resources continue to evolve, it is important to understand how representations of anatomy can facilitate the sharing and integration of data. In addition, attention to the design of the interfaces is critical to make online resources useful and usable.

RESULTS

I first present a survey of online atlases and gene expression resources for model organisms, with a focus on methods of semantic and spatial representation of anatomy. A total of 14 anatomical atlases and 21 gene expression resources are included. This survey demonstrates how choices in semantic representation, in the form of ontologies, can enhance interface search functions and provide links between relevant information. This survey also reviews methods for spatially representing anatomy in online resources. I then provide a collection of patterns for interface design based on the atlases and databases surveyed. These patterns include methods for displaying graphics, integrating semantic and spatial representations, organizing information, and querying databases to find genes expressed in anatomical structures.

CONCLUSIONS

This collection of patterns for interface design will assist biologists and software developers in planning the interfaces of new atlases and databases or enhancing existing ones. They also show the benefits of standardizing semantic and spatial representations of anatomy by demonstrating how interfaces can use standardization to provide enhanced functionality.

Data Sources for Trait Databases: Comparing the Phenomic Content of Monographs and Evolutionary Matrices

Databases of organismal traits that aggregate information from one or multiple sources can be leveraged for large-scale analyses in biology. Yet the differences among these data streams and how well they capture trait diversity have never been explored. We present the first analysis of the differences between phenotypes captured in free text of descriptive publications ('monographs') and those used in phylogenetic analyses ('matrices'). We focus our analysis on osteological phenotypes of the limbs of four extinct vertebrate taxa critical to our understanding of the fin-to-limb transition. We find that there is low overlap between the anatomical entities used in these two sources of phenotype data, indicating that phenotypes represented in matrices are not simply a subset of those found in monographic descriptions. Perhaps as expected, compared to characters found in matrices, phenotypes in monographs tend to emphasize descriptive and positional morphology, be somewhat more complex, and relate to fewer additional taxa. While based on a small set of focal taxa, these qualitative and quantitative data suggest that either source of phenotypes alone will result in incomplete knowledge of variation for a given taxon. As a broader community develops to use and expand databases characterizing organismal trait diversity, it is important to recognize the limitations of the data sources and develop strategies to more fully characterize variation both within species and across the tree of life.

Correction of the Pathogenic Alternative Splicing, Caused by the Common GNB3 c.825C>T Allele, Using a Novel, Antisense Morpholino

The very common GNB3 c.825C>T polymorphism (rs5443) is present in approximately half of all human chromosomes. Significantly, the presence of the GNB3 825T allele has been strongly associated with predisposition to essential hypertension. Paradoxically the presence of the GNB3 825T allele, in exon 10, introduces a pathogenic alternative RNA splice site into the middle of exon 9. To attempt to correct this pathogenic aberrant splicing, we, therefore, bioinformatically designed, using a Gene Tools(®) algorithm, a GNB3-specific, antisense morpholino. It was hoped that this morpholino would behave in vitro as either a potential splice blocker and/or exon skipper, to both bind and inhibit/reduce the aberrant splicing of the GNB3 825T allele. On transfecting a human lymphoblast cell line homozygous for the 825T allele, with this antisense morpholino, we encouragingly observed both a significant reduction (from ∼58% to ∼5%) in the production of the aberrant smaller GNB3 transcript, and a subsequent increase in the normal GNB3 transcript (from ∼42% to ∼95%). Our results demonstrate the potential use of a GNB3-specific antisense morpholino, as a pharmacogenetic therapy for essential hypertension.

Gene profiling of embryonic skeletal muscle lacking type I ryanodine receptor Ca(2+) release channel

In mature skeletal muscle, the intracellular Ca²⁺ concentration rises dramatically upon membrane depolarization, constituting the link between excitation and contraction. This process requires Ca²⁺ release from the sarcoplasmic reticulum via the type 1 ryanodine receptor (RYR1). However, RYR1’s potential roles in muscle development remain obscure. We used an established RyR1- null mouse model, dyspedic, to investigate the effects of the absence of a functional RYR1 and, consequently, the lack of RyR1-mediated Ca²⁺ signaling, during embryogenesis. Homozygous dyspedic mice die after birth and display small limbs and abnormal skeletal muscle organization. Skeletal muscles from front and hind limbs of dyspedic fetuses (day E18.5) were subjected to microarray analyses, revealing 318 differentially expressed genes. We observed altered expression of multiple transcription factors and members of key signaling pathways. Differential regulation was also observed for genes encoding contractile as well as muscle-specific structural proteins. Additional qRT-PCR analysis revealed altered mRNA levels of the canonical muscle regulatory factors Six1, Six4, Pax7, MyoD, MyoG and MRF4 in mutant muscle, which is in line with the severe developmental retardation seen in dyspedic muscle histology analyses. Taken together, these findings suggest an important non-contractile role of RyR1 or RYR1-mediated Ca²⁺ signaling during muscle organ development.

Prioritizing hypothesis tests for high throughput data

MOTIVATION

The advent of high throughput data has led to a massive increase in the number of hypothesis tests conducted in many types of biological studies and a concomitant increase in stringency of significance thresholds. Filtering methods, which use independent information to eliminate less promising tests and thus reduce multiple testing, have been widely and successfully applied. However, key questions remain about how to best apply them: When is filtering beneficial and when is it detrimental? How good does the independent information need to be in order for filtering to be effective? How should one choose the filter cutoff that separates tests that pass the filter from those that don't?

RESULT

We quantify the effect of the quality of the filter information, the filter cutoff and other factors on the effectiveness of the filter and show a number of results: If the filter has a high probability (e.g. 70%) of ranking true positive features highly (e.g. top 10%), then filtering can lead to dramatic increase (e.g. 10-fold) in discovery probability when there is high redundancy in information between hypothesis tests. Filtering is less effective when there is low redundancy between hypothesis tests and its benefit decreases rapidly as the quality of the filter information decreases. Furthermore, the outcome is highly dependent on the choice of filter cutoff. Choosing the cutoff without reference to the data will often lead to a large loss in discovery probability. However, naïve optimization of the cutoff using the data will lead to inflated type I error. We introduce a data-based method for choosing the cutoff that maintains control of the family-wise error rate via a correction factor to the significance threshold. Application of this approach offers as much as a several-fold advantage in discovery probability relative to no filtering, while maintaining type I error control. We also introduce a closely related method of P-value weighting that further improves performance.

AVAILABILITY AND IMPLEMENTATION

R code for calculating the correction factor is available at http://www.stat.uga.edu/people/faculty/paul-schliekelman

CONTACT

pdschlie@stat.uga.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Animal models to study acute and chronic intestinal inflammation in mammals

Acute and chronic inflammatory diseases of the intestine impart a significant and negative impact on the health and well-being of human and non-human mammalian animals. Understanding the underlying mechanisms of inflammatory disease is mandatory to develop effective treatment and prevention strategies. As inflammatory disease etiologies are multifactorial, the use of appropriate animal models and associated metrics of disease are essential. In this regard, animal models used alone or in combination to study acute and chronic inflammatory disease of the mammalian intestine paired with commonly used inflammation-inducing agents are reviewed. This includes both chemical and biological incitants of inflammation, and both non-mammalian (i.e. nematodes, insects, and fish) and mammalian (i.e. rodents, rabbits, pigs, ruminants, dogs, and non-human primates) models of intestinal inflammation including germ-free, gnotobiotic, as well as surgical, and genetically modified animals. Importantly, chemical and biological incitants induce inflammation via a multitude of mechanisms, and intestinal inflammation and injury can vary greatly according to the incitant and animal model used, allowing studies to ascertain both long-term and short-term effects of inflammation. Thus, researchers and clinicians should be aware of the relative strengths and limitations of the various animal models used to study acute and chronic inflammatory diseases of the mammalian intestine, and the scope and relevance of outcomes achievable based on this knowledge. The ability to induce inflammation to mimic common human diseases is an important factor of a successful animal model, however other mechanisms of disease such as the amount of infective agent to induce disease, invasion mechanisms, and the effect various physiologic changes can have on inducing damage are also important features. In many cases, the use of multiple animal models in combination with both chemical and biological incitants is necessary to answer the specific question being addressed regarding intestinal disease. Some incitants can induce acute responses in certain animal models while others can be used to induce chronic responses; this review aims to illustrate the strengths and weaknesses in each animal model and to guide the choice of an appropriate acute or chronic incitant to facilitate intestinal disease.

Developing the eHistology Atlas

The eMouseAtlas project has undertaken to generate a new resource providing access to high-resolution colour images of the slides used in the renowned textbook 'The Atlas of Mouse Development' by Matthew H. Kaufman. The original histology slides were digitized, and the associated anatomy annotations captured for display in the new resource. These annotations were assigned to objects in the standard reference anatomy ontology, allowing the eHistology resource to be linked to other data resources including the Edinburgh Mouse Atlas Gene-Expression database (EMAGE) an the Mouse Genome Informatics (MGI) gene-expression database (GXD). The provision of the eHistology Atlas resource was assisted greatly by the expertise of the eMouseAtlas project in delivering large image datasets within a web environment, using IIP3D technology. This technology also permits future extensions to the resource through the addition of further layers of data and annotations to the resource. Database URL: www.emouseatlas.org/emap/eHistology/index.php.

Breast Cancer Tumor Suppressors: A Special Emphasis on Novel Protein Nischarin

Tumor suppressor genes regulate cell growth and prevent spontaneous proliferation that could lead to aberrant tissue function. Deletions and mutations of these genes typically lead to progression through the cell-cycle checkpoints, as well as increased cell migration. Studies of these proteins are important as they may provide potential treatments for breast cancers. In this review, we discuss a comprehensive overview on Nischarin, a novel protein discovered by our laboratory. Nischarin, or imidazoline receptor antisera-selected protein, is a protein involved in a vast number of cellular processes, including neuronal protection and hypotension. The NISCH promoter experiences hypermethylation in several cancers, whereas some highly aggressive breast cancer cells exhibit genomic loss of the NISCH locus. Furthermore, we discuss data illustrating a novel role of Nischarin as a tumor suppressor in breast cancer. Analysis of this new paradigm may shed light on various clinical questions. Finally, the therapeutic potential of Nischarin is discussed.

A Novel de novo Mutation in CEACAM16 Associated with Postlingual Hearing Impairment

Mutations in CEACAM16 cause autosomal dominant nonsyndromic hearing loss (DFNA4B). So far, 2 families have been reported with segregating missense mutations, both in the immunoglobulin constant domain A of the CEACAM16 protein. In this study, we used the TruSight One panel to investigate a parent-child trio without familial history of hearing loss and one affected child. When filtering for recessive inheritance and de novo events, we discovered a de novo CEACAM16 mutation (c.1094T>G, p.Leu365Arg) as the sole likely pathogenic variant. The de novo mutation was confirmed by Sanger sequencing and STR analysis. The proband's hearing loss closely matches the described onset and severity for DFNA4B. We present the third CEACAM16 variant and the first de novo mutation in CEACAM16. This de novo mutation is robustly described as a pathogenic mutation according to in silico mutation prediction tools and affects a highly conserved amino acid in the most strongly conserved CEACAM16 N2 domain. Our strategy of screening family trios enhances de novo mutation discovery and the exclusion of other variants of potential interest through pedigree filtering.

A mouse informatics platform for phenotypic and translational discovery

The International Mouse Phenotyping Consortium (IMPC) is providing the world's first functional catalogue of a mammalian genome by characterising a knockout mouse strain for every gene. A robust and highly structured informatics platform has been developed to systematically collate, analyse and disseminate the data produced by the IMPC. As the first phase of the project, in which 5000 new knockout strains are being broadly phenotyped, nears completion, the informatics platform is extending and adapting to support the increasing volume and complexity of the data produced as well as addressing a large volume of users and emerging user groups. An intuitive interface helps researchers explore IMPC data by giving overviews and the ability to find and visualise data that support a phenotype assertion. Dedicated disease pages allow researchers to find new mouse models of human diseases, and novel viewers provide high-resolution images of embryonic and adult dysmorphologies. With each monthly release, the informatics platform will continue to evolve to support the increased data volume and to maintain its position as the primary route of access to IMPC data and as an invaluable resource for clinical and non-clinical researchers.

Using the phenoscape knowledgebase to relate genetic perturbations to phenotypic evolution

The abundance of phenotypic diversity among species can enrich our knowledge of development and genetics beyond the limits of variation that can be observed in model organisms. The Phenoscape Knowledgebase (KB) is designed to enable exploration and discovery of phenotypic variation among species. Because phenotypes in the KB are annotated using standard ontologies, evolutionary phenotypes can be compared with phenotypes from genetic perturbations in model organisms. To illustrate the power of this approach, we review the use of the KB to find taxa showing evolutionary variation similar to that of a query gene. Matches are made between the full set of phenotypes described for a gene and an evolutionary profile, the latter of which is defined as the set of phenotypes that are variable among the daughters of any node on the taxonomic tree. Phenoscape's semantic similarity interface allows the user to assess the statistical significance of each match and flags matches that may only result from differences in annotation coverage between genetic and evolutionary studies. Tools such as this will help meet the challenge of relating the growing volume of genetic knowledge in model organisms to the diversity of phenotypes in nature. The Phenoscape KB is available at http://kb.phenoscape.org.

The atlas of mouse development eHistology resource

The Atlas of Mouse Development by Professor Mathew Kaufman is an essential text for understanding mouse developmental anatomy. This definitive and authoritative atlas is still in production and is essential for any biologist working with the mouse embryo, although the last revision dates back to 1994. Here, we announce the eHistology online resource that provides free access to high-resolution colour images digitized from the original histological sections (www.emouseatlas.org/emap/eHistology/index.php) used by Kaufman for the Atlas. The images are provided with the original annotations and plate numbering of the paper atlas and enable viewing the material to cellular resolution.