Identification of novel Drosophila neural precursor genes using a differential embryonic head cDNA screen

During Drosophila neuroblast lineage development, temporally ordered transitions in neuroblast gene expression have been shown to accompany the changing repertoire of functionally diverse cells generated by neuroblasts. To broaden our understanding of the biological significance of these ordered transitions in neuroblast gene expression and the events that regulate them, additional genes have been sought that participate in the timing and execution of these temporally controlled events. To identify dynamically expressed neural precursor genes, we have performed a differential cDNA hybridization screen on a stage specific embryonic head cDNA library, followed by whole-mount embryo in situ hybridizations. Described here are the embryonic expression profiles of 57 developmentally regulated neural precursor genes. Information about 2389 additional genes identified in this screen, including 1614 uncharacterized genes, is available on-line at 'BrainGenes: a search for Drosophila neural precursor genes' (http://sdb.bio.purdue.edu/fly/brain/ahome.htm).     

Gene expression profiling of coronary artery disease and its relation with different severities

Global gene expression profiling is a powerful tool enabling the understanding of pathophysiology and subsequent management of diseases. This study aims to explore functionally annotated differentially expressed genes (DEGs); their biological processes for coronary artery disease (CAD) and its different severities of atherosclerotic lesions. This study also aims to identify the change in expression patterns of DEGs in atherosclerotic lesions of single-vessel disease (SVD) and triple-vessel disease (TVD). The weight of different severities of lesion was estimated using a modified Gensini score. The gene expression profiling was performed using the Affymetrix microarray platform. The functional annotation for CAD was performed using DAVID v6.8. The biological network gene ontology tool (BiNGO) and ClueGO were used to explore the biological processes of functionally annotated genes of CAD. The changes in gene expression from SVD to TVD were determined by evaluating the fold change. Functionally annotated genes were found in an unique set and could be distinguishing two distinct severities of CAD. The biological processes such as cellular migration, locomotion, cell adhesion, cytokine production, positive regulation of cell death etc. enriched the functionally annotated genes in SVD, whereas, wound healing, negative regulation of cell death, blood coagulation, angiogenesis and fibrinolysis were enriched significantly in TVD patients. The genes THBS1 and CAPN10 were functionally annotated for CAD in both SVD and TVD. The 61 DEGs were identified, those have changes their expression with different severities of atherosclerotic lesions, in which 13 genes had more than two-fold change in expression between SVD and TVD. The consistent findings were obtained on validation of microarray gene expression of selected 10 genes in a separate cohort using real-time PCR. This study identified putative candidate genes and their biological processes predisposing toward and affecting the severity of CAD.     

Comparative Study of Regulatory Circuits in Two Sea Urchin Species Reveals Tight Control of Timing and High Conservation of Expression Dynamics

Accurate temporal control of gene expression is essential for normal development and must be robust to natural genetic and environmental variation. Studying gene expression variation within and between related species can delineate the level of expression variability that development can tolerate. Here we exploit the comprehensive model of sea urchin gene regulatory networks and generate high-density expression profiles of key regulatory genes of the Mediterranean sea urchin, Paracentrotus lividus (Pl). The high resolution of our studies reveals highly reproducible gene initiation times that have lower variation than those of maximal mRNA levels between different individuals of the same species. This observation supports a threshold behavior of gene activation that is less sensitive to input concentrations. We then compare Mediterranean sea urchin gene expression profiles to those of its Pacific Ocean relative, Strongylocentrotus purpuratus (Sp). These species shared a common ancestor about 40 million years ago and show highly similar embryonic morphologies. Our comparative analyses of five regulatory circuits operating in different embryonic territories reveal a high conservation of the temporal order of gene activation but also some cases of divergence. A linear ratio of 1.3-fold between gene initiation times in Pl and Sp is partially explained by scaling of the developmental rates with temperature. Scaling the developmental rates according to the estimated Sp-Pl ratio and normalizing the expression levels reveals a striking conservation of relative dynamics of gene expression between the species. Overall, our findings demonstrate the ability of biological developmental systems to tightly control the timing of gene activation and relative dynamics and overcome expression noise induced by genetic variation and growth conditions.     

Identification of novel members reveals the structural and functional divergence of lepidopteran-specific Lipoprotein_11 family

30K proteins (30KPs) are classified into the lepidopteran-specific Lipoprotein_11 family. They are involved in various physiological processes such as energy storage, embryonic development, and immune response in the silkworm. To date, 30KPs were only found in Bombyx mori and Manduca sexta. Moreover, the C-termini of ENF peptide binding proteins (ENF-BPs) show similarity to 30KPs. ENF peptides are multifunctional insect cytokines and involved in growth regulation and defense reaction, whereas ENF-BPs act as active regulators of ENF peptides. In order to get insights into this gene family in Lepidoptera, we performed an extensive survey of lepidopteran-derived genome and EST datasets. We identified 73 30KP homologous genes in 12 lepidopteran species, of which 56 are novel members. The structural and phylogenetic analyses revealed that these genes could be classified into three groups: ENF-BP genes, typical 30KP genes, and serine/threonine-rich 30KP (S/T-rich 30KP) genes. The C-terminal regions are common to all the three subfamilies, but the N-termini are highly variable. We found a novel subfamily of Lipoprotein_11 and named it S/T-rich 30KP according to its exclusive S/T-rich domain in the N terminus. ENF-BP was also found to contain a special domain in the N terminus, which is homologous to Pp-0912 of Pseudomonas putida. Microarray data and semi-quantitative RT-PCR showed that the three groups have their respective temporal–spatial expression patterns. S/T-rich 30KP genes have enriched expression in the mature testis and might be involved in spermiogenesis or fertilization. Typical 30KP genes are expressed mainly in the fat body and integument at the larvae and pupae stages. ENF-BP genes are expressed predominantly in the hemocyte. The differential spatial–temporal expression profiles revealed the functional divergence of three Lipoprotein_11 subfamilies.     

A mRNA landscape of bovine embryos after standard and MAPK-inhibited culture conditions: a comparative analysis

Background

Genes and signalling pathways involved in pluripotency have been studied extensively in mouse and human pre-implantation embryos and embryonic stem (ES) cells. The unsuccessful attempts to generate ES cell lines from other species including cattle suggests that other genes and pathways are involved in maintaining pluripotency in these species. To investigate which genes are involved in bovine pluripotency, expression profiles were generated from morula, blastocyst, trophectoderm and inner cell mass (ICM) samples using microarray analysis. As MAPK inhibition can increase the NANOG/GATA6 ratio in the inner cell mass, additionally blastocysts were cultured in the presence of a MAPK inhibitor and changes in gene expression in the inner cell mass were analysed.

Results

Between morula and blastocyst 3,774 genes were differentially expressed and the largest differences were found in blastocyst up-regulated genes. Gene ontology (GO) analysis shows lipid metabolic process as the term most enriched with genes expressed at higher levels in blastocysts. Genes with higher expression levels in morulae were enriched in the RNA processing GO term. Of the 497 differentially expressed genes comparing ICM and TE, the expression of NANOG, SOX2 and POU5F1 was increased in the ICM confirming their evolutionary preserved role in pluripotency. Several genes implicated to be involved in differentiation or fate determination were also expressed at higher levels in the ICM. Genes expressed at higher levels in the ICM were enriched in the RNA splicing and regulation of gene expression GO term. Although NANOG expression was elevated upon MAPK inhibition, SOX2 and POU5F1 expression showed little increase. Expression of other genes in the MAPK pathway including DUSP4 and SPRY4, or influenced by MAPK inhibition such as IFNT, was down-regulated.

Conclusion

The data obtained from the microarray studies provide further insight in gene expression during bovine embryonic development. They show an expression profile in pluripotent cells that indicates a pluripotent, epiblast-like state. The inability to culture ICM cells as stem cells in the presence of an inhibitor of MAPK activity together with the reported data indicates that MAPK inhibition alone is not sufficient to maintain a pluripotent character in bovine cells.

Electronic supplementary material

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

Robustness under Functional Constraint: The Genetic Network for Temporal Expression in Drosophila Neurogenesis

Precise temporal coordination of gene expression is crucial for many developmental processes. One central question in developmental biology is how such coordinated expression patterns are robustly controlled. During embryonic development of the Drosophila central nervous system, neural stem cells called neuroblasts express a group of genes in a definite order, which leads to the diversity of cell types. We produced all possible regulatory networks of these genes and examined their expression dynamics numerically. From the analysis, we identified requisite regulations and predicted an unknown factor to reproduce known expression profiles caused by loss-of-function or overexpression of the genes in vivo, as well as in the wild type. Following this, we evaluated the stability of the actual Drosophila network for sequential expression. This network shows the highest robustness against parameter variations and gene expression fluctuations among the possible networks that reproduce the expression profiles. We propose a regulatory module composed of three types of regulations that is responsible for precise sequential expression. This study suggests that the Drosophila network for sequential expression has evolved to generate the robust temporal expression for neuronal specification.     

Relocation of genes generates non-conserved chromosomal segments in Fusarium graminearum that show distinct and co-regulated gene expression patterns

Background

Genome comparisons between closely related species often show non-conserved regions across chromosomes. Some of them are located in specific regions of chromosomes and some are even confined to one or more entire chromosomes. The origin and biological relevance of these non-conserved regions are still largely unknown. Here we used the genome of Fusarium graminearum to elucidate the significance of non-conserved regions.

Results

The genome of F. graminearum harbours thirteen non-conserved regions dispersed over all of the four chromosomes. Using RNA-Seq data from the mycelium of F. graminearum, we found weakly expressed regions on all of the four chromosomes that exactly matched with non-conserved regions. Comparison of gene expression between two different developmental stages (conidia and mycelium) showed that the expression of genes in conserved regions is stable, while gene expression in non-conserved regions is much more influenced by developmental stage. In addition, genes involved in the production of secondary metabolites and secreted proteins are enriched in non-conserved regions, suggesting that these regions could also be important for adaptations to new environments, including adaptation to new hosts. Finally, we found evidence that non-conserved regions are generated by sequestration of genes from multiple locations. Gene relocations may lead to clustering of genes with similar expression patterns or similar biological functions, which was clearly exemplified by the PKS2 gene cluster.

Conclusions

Our results showed that chromosomes can be functionally divided into conserved and non-conserved regions, and both could have specific and distinct roles in genome evolution and regulation of gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-191) contains supplementary material, which is available to authorized users.     

Aberrant expression of imprinted genes in post-implantation rat embryos

AimImprinted genes are known regulators of embryo growth. Studies from our laboratory have demonstrated that treatment of adult male rats with tamoxifen increased post-implantation loss at around midgestation. Expression of insulin like growth factor 2 (Igf2), a paternally expressed imprinted gene was down-regulated in the resorbing embryos obtained at embryonic day 13. Hypomethylation of Igf2-H19 imprint control region was observed in the resorbing embryo sires and spermatozoa obtained from tamoxifen-treated rats thereby suggesting that errors in imprint acquisition during spermatogenesis can result in embryo loss. The present study aims at studying the expression of other imprinted genes, besides Igf2 in the embryos sired by tamoxifen-treated males.Main methodsGene expression profiles of resorbing versus normal embryos were assessed by microarrays. Real time quantitative RT-PCR for six imprinted genes and four genes involved in cell cycle was done to validate gene expression data. The affected pathways and functions were identified in the resorbing embryos and effect on cell cycle was confirmed by flow cytometry.Key findingsAberrant expression of a number of imprinted genes was observed in the resorbing embryos when compared to the normal embryos at embryonic days 11 and 13. Down-regulation of Notch signaling, Wnt signaling and cell cycle pathway was observed in the resorbing embryos.SignificanceThe study suggests that exposure of male germ cells to tamoxifen during adulthood results in aberrant expression of imprinted genes and down-regulation of development associated pathways in the F₁ progeny thereby causing embryo loss.     

Reverse engineering and analysis of genome-wide gene regulatory networks from gene expression profiles using high-performance computing

Regulation of gene expression is a carefully regulated phenomenon in the cell. “Reverse-engineering” algorithms try to reconstruct the regulatory interactions among genes from genome-scale measurements of gene expression profiles (microarrays). Mammalian cells express tens of thousands of genes; hence, hundreds of gene expression profiles are necessary in order to have acceptable statistical evidence of interactions between genes. As the number of profiles to be analyzed increases, so do computational costs and memory requirements. In this work, we designed and developed a parallel computing algorithm to reverse-engineer genome-scale gene regulatory networks from thousands of gene expression profiles. The algorithm is based on computing pairwise Mutual Information between each gene-pair. We successfully tested it to reverse engineer the Mus Musculus (mouse) gene regulatory network in liver from gene expression profiles collected from a public repository. A parallel hierarchical clustering algorithm was implemented to discover “communities” within the gene network. Network communities are enriched for genes involved in the same biological functions. The inferred network was used to identify two mitochondrial proteins.     

A Computational Approach towards a Gene Regulatory Network for the Developing Nematostella vectensis Gut

Background

The starlet sea anemone Nematostella vectensis is a diploblastic cnidarian that expresses a set of conserved genes for gut formation during its early development. During the last decade, the spatial distribution of many of these genes has been visualized with RNA hybridization or protein immunolocalization techniques. However, due to N. vectensis'' curved and changing morphology, quantification of these spatial data is problematic. A method is developed for two-dimensional gene expression quantification, which enables a numerical analysis and dynamic modeling of these spatial patterns.

Methods/Result

In this work, first standardized gene expression profiles are generated from publicly available N. vectensis embryo images that display mRNA and/or protein distributions. Then, genes expressed during gut formation are clustered based on their expression profiles, and further grouped based on temporal appearance of their gene products in embryonic development. Representative expression profiles are manually selected from these clusters, and used as input for a simulation-based optimization scheme. This scheme iteratively fits simulated profiles to the selected profiles, leading to an optimized estimation of the model parameters. Finally, a preliminary gene regulatory network is derived from the optimized model parameters.

Outlook

While the focus of this study is N. vectensis, the approach outlined here is suitable for inferring gene regulatory networks in the embryonic development of any animal, thus allowing to comparatively study gene regulation of gut formation in silico across various species.