Mammalian Genes Research Articles

Evolutionarily conserved inhibitory uORFs sensitize Hox mRNA translation to start codon selection stringency

Translation start site selection in eukaryotes is influenced by context nucleotides flanking the AUG codon and by levels of the eukaryotic translation initiation factors eIF1 and eIF5. In a search of mammalian genes, we identified five homeobox (Hox) gene paralogs initiated by AUG codons in conserved suboptimal context as well as 13 Hox genes that contain evolutionarily conserved upstream open reading frames (uORFs) that initiate at AUG codons in poor sequence context. An analysis of published cap analysis of gene expression sequencing (CAGE-seq) data and generated CAGE-seq data for messenger RNAs (mRNAs) from mouse somites revealed that the 5' leaders of Hox mRNAs of interest contain conserved uORFs, are generally much shorter than reported, and lack previously proposed internal ribosome entry site elements. We show that the conserved uORFs inhibit Hox reporter expression and that altering the stringency of start codon selection by overexpressing eIF1 or eIF5 modulates the expression of Hox reporters. We also show that modifying ribosome homeostasis by depleting a large ribosomal subunit protein or treating cells with sublethal concentrations of puromycin leads to lower stringency of start codon selection. Thus, altering global translation can confer gene-specific effects through altered start codon selection stringency.

Prediction of transcription factors associated with DNA demethylation during human cellular development

DNA methylation of CpG dinucleotides is an important epigenetic modification involved in the regulation of mammalian gene expression, with each type of cell developing a specific methylation profile during its differentiation. Recently, it has been shown that a small subgroup of transcription factors (TFs) might promote DNA demethylation at their binding sites. We developed a bioinformatics pipeline to predict from genome-wide DNA methylation data TFs that promote DNA demethylation at their binding site. We applied the pipeline to International Human Epigenome Consortium methylome data and selected 393 candidate transcription factor binding motifs and associated 383 TFs that are likely associated with DNA demethylation. Validation of a subset of the candidate TFs using an in vitro assay suggested that 28 of 49 TFs from various TF families had DNA-demethylation-promoting activity; TF families, such as bHLH and ETS, contained both TFs with and without the activity. The identified TFs showed large demethylated/methylated CpG ratios and their demethylated CpGs showed significant bias toward hypermethylation in original cells. Furthermore, the identified TFs promoted demethylation of distinct sets of CpGs, with slight overlap of the targeted CpGs among TF family members, which was consistent with the results of a gene ontology (GO) term analysis of the identified TFs. Gene expression analysis of the identified TFs revealed that multiple TFs from various families are specifically expressed in human cells and tissues. Together, our results suggest that a large number of TFs from various TF families are associated with cell-type-specific DNA demethylation during human cellular development.

Natural selection shaped the evolution of amino acid usage in mammalian toll like receptor genes

AbstractToll-like receptors (TLRs) are important as they are able to sense diverse set of pathogens associated molecular patterns (PAMPs) as ligands. These receptors are involved in functions such as immune response, development of signaling process and cell adhesion. In the present study we are interested to analyze the influence of evolutionary selection pressure on the mutational diversity of mammalian TLR genes. We observed differential patterns of amino acid usage between primate and non-primate mammalian TLR genes. GC-content of TLR genes and hydrophobicity of the encoded proteins are the most influential factors correlated with the differential pattern of amino acid usage.The influence of the subcellular location on the amino acid usage pattern of TLRs is evident in present study. Purifying selection is uniformly present on TLR genes, positively selected sites are mostly located over the ligand binding domain. Our study clearly demonstrates that natural selection has shaped the evolution of primate and non-primate mammalian TLR genes.

Identification of a Schistosoma japonicum MicroRNA That Suppresses Hepatoma Cell Growth and Migration by Targeting Host FZD4 Gene.

Previous studies have demonstrated miRNAs derived from plants and parasites can modulate mammalian gene expression and cell phenotype in a cross-kingdom manner, leading to occurrence of diseases or strengthening resistance of host to diseases such as cancer. In this study, we identified a schistosome miRNA (named Sja-miR-71a) through screening of 57 Schistosoma japonicum miRNAs that exerts antitumor activity in vitro and in vivo models. We demonstrated presence of this parasite miRNA in liver cells during infection. We showed that Sja-miR-71a arrested cell cycle at G0/G1 phase of hepatoma cell lines and inhibited cell proliferation in vitro. The HepG2 transfected with Sja-miR-71a mimics displayed significant reduction of migration and colony formation. Further, growth of the tumor cells transfected with the Sja-miR-71a mimics was obviously suppressed in a xenograft mouse model. Mechanically, we found the antitumor activity of Sja-miR-71a was through targeting a host gene encoding Frizzled Class Receptor 4 (FZD4), as FZD4 small interfering RNAs (siRNAs) generated a similar inhibitory effect on the tumor. These data indicated that Sja-miR-71a is a tumor suppressor miRNA and suggested this parasite-derived miRNA as a potential therapeutic target for cancer.

A generalized moment-based method for estimating parameters of stochastic gene transcription

Gene transcription is a stochastic process manifested by the mRNA distribution data at the single-cell level, which can be fitted using the mass function Pm that quantifies the probability of m mRNA molecules per cell for the gene of interest. Extensive studies have been conducted to infer gene parameter rates in mathematical models using maximum-likelihood method or moment-based method based on mRNA distribution data. However, current methods usually do not emphasize the exact fit to the data of single P0, where P0 plays a critical role in discriminating different distribution modalities. In this study, we highlighted the fit to the data of P0, and developed a generalized moment-based protocol for more reliable parameter estimation in the classical two-state model. Our protocol is easy to follow and allows the estimation of parameters of E.coli and mammalian genes under varying conditions. We showed that our protocol performs much better than the traditional moment-based method, and more likely captures the bimodal mRNA distribution than that obtained using the maximum-likelihood method.

Manipulating Galectin Expression in Zebrafish (Danio rerio).

Techniques for disrupting gene expression are invaluable tools for the analysis of the biological role of a gene product. Because of its genetic tractability and multiple advantages over conventional mammalian models, the zebrafish (Danio rerio) is recognized as a powerful system for gaining new insight into diverse aspects of human health and disease. Among the multiple mammalian gene families for which the zebrafish has shown promise as an invaluable model for functional studies, the galectins have attracted great interest due to their participation in early development, regulation of immune homeostasis, and recognition of microbial pathogens. Galectins are β-galactosyl-binding lectins with a characteristic sequence motif in their carbohydrate recognition domains (CRDs), that constitute an evolutionary conserved family ubiquitous in eukaryotic taxa. Galectins are emerging as key players in the modulation of many important pathological processes, which include acute and chronic inflammatory diseases, autoimmunity and cancer, thus making them potential molecular targets for innovative drug discovery. Here, we provide a review of the current methods available for the manipulation of gene expression in the zebrafish, with a focus on gene knockdown [morpholino (MO)-derived antisense oligonucleotides] and knockout (CRISPR-Cas) technologies.

Finding and Verifying Enhancers for Endothelial-Expressed Genes.

Identification and analysis of enhancers for endothelial-expressed genes can provide crucial information regarding their upstream transcriptional regulators. However, enhancer identification can be challenging, particularly for people with limited access or experience of bioinformatics, and transgenic analysis of enhancer activity patterns can be prohibitively expensive. Here we describe how to use publicly available datasets displayed on the UCSC Genome Browser to identify putative endothelial enhancers for mammalian genes. Furthermore, we detail how to utilize mosaic Tol2-mediated transgenesis in zebrafish to verify whether a putative enhancer is capable of directing endothelial-specific patterns of gene expression.

Principles for assessing the genotoxicity of carbon nanomaterials in vitro (on the example of carbon nanotubes) (literature review)

Introduction. Genotoxicity of nanomaterials (NM) is becoming a major concern when investigating new NM for their safety. Each mutagen is considered to be potentially carcinogenic, therefore a genotoxicity assessment is necessary. However, a clear strategy for assessing the genotoxic effect of NM has not yet been developed. Material and methods. The material for the analysis have included literature sources from the bibliographic databases PubMed, Scopus, RSCI. Results. Physicochemical characterization of NM is carried out using high-resolution microscopic and light scattering methods. Before testing for genotoxicity, it is necessary to know the cytotoxicity of the tested NM in order to select the appropriate concentration range. The most important and significant tests are based on the cell viability. MTT assay is a colorimetric test that evaluates the metabolic activity of cells. In addition, viability can be determined using microscopy, flow cytometry, determination of lactate dehydrogenase. Genotoxicity evaluation can be carried out only after the preliminary steps. The strategy should include genotoxicity endpoints: DNA damage, gene mutations, chromosomal damage. The in vitro mammalian gene mutation test, usually performed using mouse lymphoma cells, detects a wide range of genetic damage, including gene deletions. The most common test for detecting chromosomal damage is an in vitro micronucleus assay. DNA strand breaks are most often assessed using the comet DNA assay. Conclusion. Compulsory stages in the study of the genotoxicity of nanomaterials should be preliminary studies, including physicochemical characterization and assessment of cytotoxicity, as well as the study of the endpoints of genotoxicity and potential mechanisms.

Novel Insights Into the Multifaceted Functions of RNA n6-Methyladenosine Modification in Degenerative Musculoskeletal Diseases.

N6-methyladenosine (m6A) is an important modification of eukaryotic mRNA. Since the first discovery of the corresponding demethylase and the subsequent identification of m6A as a dynamic modification, the function and mechanism of m6A in mammalian gene regulation have been extensively investigated. “Writer”, “eraser” and “reader” proteins are key proteins involved in the dynamic regulation of m6A modifications, through the anchoring, removal, and interpretation of m6A modifications, respectively. Remarkably, such dynamic modifications can regulate the progression of many diseases by affecting RNA splicing, translation, export and degradation. Emerging evidence has identified the relationship between m6A modifications and degenerative musculoskeletal diseases, such as osteoarthritis, osteoporosis, sarcopenia and degenerative spinal disorders. Here, we have comprehensively summarized the evidence of the pathogenesis of m6A modifications in degenerative musculoskeletal diseases. Moreover, the potential molecular mechanisms, regulatory functions and clinical implications of m6A modifications are thoroughly discussed. Our review may provide potential prospects for addressing key issues in further studies.

RNA sequencing-based screen for reactivation of silenced alleles of autosomal genes.

In mammalian cells, maternal and paternal alleles usually have similar transcriptional activity. Epigenetic mechanisms such as X-chromosome inactivation (XCI) and imprinting were historically viewed as rare exceptions to this rule. Discovery of autosomal monoallelic autosomal expression (MAE) a decade ago revealed an additional allele-specific mode regulating thousands of mammalian genes. Despite MAE prevalence, its mechanistic basis remains unknown. Using an RNA sequencing-based screen for reactivation of silenced alleles, we identified DNA methylation as key mechanism of MAE mitotic maintenance. In contrast with the all-or-nothing allelic choice in XCI, allele-specific expression in MAE loci is tunable, with exact allelic imbalance dependent on the extent of DNA methylation. In a subset of MAE genes, allelic imbalance was insensitive to DNA demethylation, implicating additional mechanisms in MAE maintenance in these loci. Our findings identify a key mechanism of MAE maintenance and provide basis for understanding the biological role of MAE.

Inflammation and convergent placenta gene co-option contributed to a novel reproductive tissue

SummaryThe evolution of pregnancy exposes parental tissues to new, potentially stressful conditions, which can trigger inflammation.1 Inflammation is costly2,3 and can induce embryo rejection, which constrains the evolution of pregnancy.1 In contrast, inflammation can also promote morphological innovation at the maternal-embryonic interface as exemplified by co-option of pro-inflammatory signaling for eutherian embryo implantation.1,4,5 Given its dual function, inflammation could be a key process explaining how innovations such as pregnancy and placentation evolved many times convergently. Pelvic brooding ricefishes evolved a novel “plug” tissue,6,7 which forms inside the female gonoduct after spawning, anchors egg-attaching filaments, and enables pelvic brooders to carry eggs externally until hatching.6,8 Compared to pregnancy, i.e., internal bearing of embryos, external bearing should alleviate constraints on inflammation in the reproductive tract. We thus hypothesized that an ancestral inflammation triggered by the retention of attaching filaments gave rise to pathways orchestrating plug formation. In line with our hypothesis, histological sections of the developing plug revealed signs of gonoduct injuries by egg-attaching filaments in the pelvic brooding ricefish Oryzias eversi. Tissue-specific transcriptomes showed that inflammatory signaling dominates the plug transcriptome and inflammation-induced genes controlling vital processes for plug development such as tissue growth and angiogenesis were overexpressed in the plug. Finally, mammalian placenta genes were enriched in the plug transcriptome, indicating convergent gene co-option for building, attaching, and sustaining a transient tissue in the female reproductive tract. This study highlights the role of gene co-option and suggests that recruiting inflammatory signaling into physiological processes provides a fast-track to evolutionary innovation.

Genotoxicity assessment of cellulose nanofibrils using a standard battery of in vitro and in vivo assays

Cellulose nanofibrils (CNFs) are identified as novel nanomaterials with many potential applications. Since CNFs are fibrous manufactured nanomaterials, their potential carcinogenic effects and mesothelial toxicity raise some concerns. In this study, we conducted a standard battery of in vitro and in vivo assays to evaluate the genotoxicity of two CNF types using different manufacturing methods and physicochemical properties. Namely, one was CNF produced via chemical modification by TEMPO (2,2,6,6-tetramethylpiperidine-1-oxyl radical)-mediated oxidation, while the other was CNF produced via mechanical defibrillation using needle bleached kraft pulp. A bacterial reverse mutation test and a mouse lymphoma TK assay revealed that CNFs at 100 μg/mL did not induce bacterial reverse mutations and in vitro mammalian cell gene mutation. Further, in vitro chromosomal aberration tests demonstrated that CNFs at 100 μg/mL did not induce chromosomal aberration in Chinese hamster lung fibroblasts. From the mammalian erythrocyte micronucleus test, no statistically significant increase was observed in the proportion of micronucleated polychromatic erythrocytes in the bone marrow cells of rats intratracheally instilled with any concentration of CNFs (0.25–1.0 mg/kg) compared with values from respective negative control groups. Therefore, this battery of in vitro and in vivo assays illustrated that the CNFs examined in this study did not induce genotoxicity, suggesting our results provide valuable insight on the future use of these materials in various industrial applications.

The Glutamate Receptor Plays a Role in Defense against Botrytis cinerea through Electrical Signaling in Tomato

Plant glutamate-like receptor genes (GLRs) are homologous to mammalian ionotropic glutamate receptors genes (iGluRs). Although GLRs have been implicated in plant defenses to biotic stress, the relationship between GLR-mediated plant immunity against fungal pathogens and electrical signals remains poorly understood. Here, we found that pretreatment with a GLR inhibitor, 6,7-dinitriquinoxaline-2,3-dione (DNQX), increased the susceptibility of tomato plants to the necrotrophic fungal pathogen Botrytis cinerea. Assessment of the glr3.3, glr3.5 and glr3.3/glr3.5 double-mutants upon B. cinerea infection showed that tomato GLR3.3 and GLR3.5 are essential for plant immunity against B. cinerea, wherein GLR3.3 plays the main role. Analysis of the membrane potential changes induced by glutamate (Glu) or glycine (Gly) revealed that amplitude was significantly reduced by knocking out GLR3.3 in tomato. While treatment with Glu or Gly significantly increased immunity against B. cinerea in wild-type plants, this effect was significantly attenuated in glr3.3 mutants. Thus, our data demonstrate that GLR3.3- and GLR3.5-mediated plant immunity against B. cinerea is associated with electrical signals in tomato plants.

Gene expression changes in mammalian hosts during schistosomiasis: a review

Schistosomiasis affects over 250 million people worldwide with an estimated mortality of more than 200,000 deaths per year in sub-Saharan Africa. Efforts to control schistosomiasis in the affected areas have mainly relied on mass administration of praziquantel, which kills adult but not immature worms of all Schistosoma species. Mammalian hosts respond differently to Schistosoma infection with some being more susceptible than others, which is associated with risk factors such as sociodemographic, epidemiological, immunological and/or genetic. Host genetic factors play a major role in influencing molecular processes in response to schistosomiasis as shown in gene expression studies. These studies highlight gene profiles expressed at different time points of infection using model animals. Immune function related genes; cytokines (Th1 and Th17) are upregulated earlier in infection and Th2 upregulated later indicating a mixed Th1/Th2 response. However, Th1 response has been shown to be sustained in S. japonicum infection. Immune mediators such as matrix metalloproteinases (Mmps) and tissue inhibitors of matrix metalloproteinases (Timps) are expressed later in the infection and these are linked to wound healing and fibrosis. Downregulation of metabolic associated genes is recorded in later stages of infection. Most mammalian host gene expression studies have been done using rodent models, with fewer in larger hosts such as bovines and humans. The majority of these studies have focused on S. japonicum infections and less on S. haematobium and S. mansoni infections (the two species that cause most global infections). The few human schistosomiasis gene expression studies so far have focused on S. japonicum and S. haematobium infections and none on S. mansoni, as far as we are aware. This highlights a paucity of gene expression data in humans, specifically with S. mansoni infection. This data is important to understand the disease pathology, identify biomarkers, diagnostics and possible drug targets.

MUC4 is not expressed in cell lines used for live cell imaging.

Background: The ability to visualise specific mammalian gene loci in living cells is important for understanding the dynamic processes linked to transcription. However, some of the tools used to target mammalian genes for live cell imaging, such as dCas9, have been reported to themselves impede processes linked to transcription. The MUC4 gene is a popular target for live cell imaging studies due to the repetitive nature of sequences within some exons of this gene. Methods: We set out to compare the impact of dCas9 and TALE-based imaging tools on MUC4 expression, including in human cell lines previously reported as expressing MUC4. Results: We were unable to detect MUC4 mRNA in these cell lines. Moreover, analysis of publicly available data for histone modifications associated with transcription, and data for transcription itself, indicate that neither MUC4, nor any of the mucin gene family are significantly expressed in the cell lines where dCas9 targeting has been reported to repress MUC4 and MUC1 expression, or in the cell lines where dCas13 has been used to report MUC4 RNA detection in live cells. Conclusions: Methods for visualising specific gene loci and gene transcripts in live human cells are very challenging. Our data suggest that care should be given to the choice of the most appropriate cell lines for these analyses and that orthogonal methods of assaying gene expression be carefully compared.

Systematic Studies of the Circadian Clock Genes Impact on Temperature Compensation and Cell Proliferation Using CRISPR Tools.

Simple SummaryOne of the major characteristics of the circadian clock is temperature compensation, and previous studies suggested a single clock gene may determine the temperature compensation. In this study, we report the first full collection of clock gene knockout cell lines using CRISPR/Cas9 tools. Our full collections indicate that the temperature compensation is a complex gene regulation system instead of being regulated by any single gene. Besides, we systematically compared the proliferation rates and circadian periods using our full collections, and we found that the cell growth rate is not dependent on the circadian period. Therefore, complex interaction between clock genes and their protein products may underlie the mechanism of temperature compensation, which needs further investigations.Mammalian circadian genes are capable of producing a self-sustained, autonomous oscillation whose period is around 24 h. One of the major characteristics of the circadian clock is temperature compensation. However, the mechanism underlying temperature compensation remains elusive. Previous studies indicate that a single clock gene may determine the temperature compensation in several model organisms. In order to understand the influence of each individual clock gene on the temperature compensation, twenty-three well-known mammalian clock genes plus Timeless and Myc genes were knocked out individually, using a powerful gene-editing tool, CRISPR/Cas9. First, Bmal1, Cry1, and Cry2 were knocked out as examples to verify that deleting genes by CRISPR is effective and precise. Cell lines targeting twenty-two genes were successfully edited in mouse fibroblast NIH3T3 cells, and off-target analysis indicated these genes were correctly knocked out. Through measuring the luciferase reporters, the circadian periods of each cell line were recorded under two different temperatures, 32.5 °C and 37 °C. The temperature compensation coefficient Q10 was subsequently calculated for each cell line. Estimations of the Q10 of these cell lines showed that none of the individual cell lines can adversely affect the temperature compensation. Cells with a longer period at lower temperature tend to have a shorter period at higher temperature, while cells with a shorter period at lower temperature tend to be longer at higher temperature. Thus, the temperature compensation is a fundamental property to keep cellular homeostasis. We further conclude that the temperature compensation is a complex gene regulation system instead of being regulated by any single gene. We also estimated the proliferation rates of these cell lines. After systematically comparing the proliferation rates and circadian periods, we found that the cell growth rate is not dependent on the circadian period.

MUC4 is not expressed in cell lines used for live cell imaging

Background: The ability to visualise specific mammalian gene loci in living cells is important for understanding the dynamic processes linked to transcription. However, some of the tools used to target mammalian genes for live cell imaging, such as dCas9, have been reported to themselves impede processes linked to transcription. The MUC4 gene is a popular target for live cell imaging studies due to the repetitive nature of sequences within some exons of this gene. Methods: We set out to compare the impact of dCas9 and TALE-based imaging tools on MUC4 expression, including in human cell lines previously reported as expressing MUC4. Results: We were unable to detect MUC4 mRNA in these cell lines. Moreover, analysis of publicly available data for histone modifications associated with transcription, and data for transcription itself, indicate that neither MUC4, nor any of the mucin gene family are significantly expressed in the cell lines where dCas9 targeting has been reported to repress MUC4 and MUC1 expression, or in the cell lines where dCas13 has been used to report MUC4 RNA detection in live cells. Conclusions: Methods for visualising specific gene loci and gene transcripts in live human cells are very challenging. Our data suggest that care should be given to the choice of the most appropriate cell lines for these analyses and that orthogonal methods of assaying gene expression be carefully compared.

A de novo genome assembly and annotation of the southern flying squirrel (Glaucomys volans).

Northern (Glaucomys sabrinus) and southern (Glaucomys volans) flying squirrels are widespread species distributed across North America. Northern flying squirrels are common inhabitants of the boreal forest, also occurring in coniferous forest remnants farther south, whereas the southern flying squirrel range is centered in eastern temperate woodlands. These two flying squirrel species exhibit a hybrid zone across a latitudinal gradient in an area of recent secondary contact. Glaucomys hybrid offspring are viable and can successfully backcross with either parental species, however, the fitness implications of such events are currently unknown. Some populations of G. sabrinus are endangered, and thus, interspecific hybridization is a key conservation concern in flying squirrels. To provide a resource for future studies to evaluate hybridization and possible introgression, we sequenced and assembled a de novo long-read genome from a G. volans individual sampled in southern Ontario, Canada, while four short-read genomes (two G. sabrinus and two G. volans, all from Ontario) were resequenced on Illumina platforms. The final genome assembly consisted of approximately 2.40 Gb with a scaffold N50 of 455.26 Kb. Benchmarking Universal Single-Copy Orthologs reconstructed 3,742 (91.2%) complete mammalian genes and genome annotation using RNA-Seq identified the locations of 19,124 protein-coding genes. The four short-read individuals were aligned to our reference genome to investigate the demographic history of the two species. A principal component analysis clearly separated resequenced individuals, while inferring population size history using the Pairwise Sequentially Markovian Coalescent model noted an approximate species split 1 million years ago, and a single, possibly recently introgressed individual.

An ancient retroviral RNA element hidden in mammalian genomes and its involvement in co-opted retroviral gene regulation

BackgroundRetroviruses utilize multiple unique RNA elements to control RNA processing and translation. However, it is unclear what functional RNA elements are present in endogenous retroviruses (ERVs). Gene co-option from ERVs sometimes entails the conservation of viral cis-elements required for gene expression, which might reveal the RNA regulation in ERVs.ResultsHere, we characterized an RNA element found in ERVs consisting of three specific sequence motifs, called SPRE. The SPRE-like elements were found in different ERV families but not in any exogenous viral sequences examined. We observed more than a thousand of copies of the SPRE-like elements in several mammalian genomes; in human and marmoset genomes, they overlapped with lineage-specific ERVs. SPRE was originally found in human syncytin-1 and syncytin-2. Indeed, several mammalian syncytin genes: mac-syncytin-3 of macaque, syncytin-Ten1 of tenrec, and syncytin-Car1 of Carnivora, contained the SPRE-like elements. A reporter assay revealed that the enhancement of gene expression by SPRE depended on the reporter genes. Mutation of SPRE impaired the wild-type syncytin-2 expression while the same mutation did not affect codon-optimized syncytin-2, suggesting that SPRE activity depends on the coding sequence.ConclusionsThese results indicate multiple independent invasions of various mammalian genomes by retroviruses harboring SPRE-like elements. Functional SPRE-like elements are found in several syncytin genes derived from these retroviruses. This element may facilitate the expression of viral genes, which were suppressed due to inefficient codon frequency or repressive elements within the coding sequences. These findings provide new insights into the long-term evolution of RNA elements and molecular mechanisms of gene expression in retroviruses.

GEREA: Prediction of Gene Expression Regulators from Transcriptome Profiling Data to Transition Networks

Background: Mammalian genes are regulated at the transcriptional and posttranscriptional levels. These mechanisms may involve the direct promotion or inhibition of transcription via a regulator or post-transcriptional regulation through factors such as micro (mi)RNAs. Objective: Construct gene regulation relationships modulated by causality inference-based miRNA- (transition factor)-(target gene) networks and analysis gene expression data to identify gene expression regulators. Methods: Mouse gene expression regulation relationships were manually curated from literature using a text mining method which were then employed to generate miRNA-(transition factor)- (target gene) networks. An algorithm was then introduced to identify gene expression regulators from transcriptome profiling data by applying enrichment analysis to these networks. Results: A total of 22,271 mouse gene expression regulation relationships were curated for 4,018 genes and 242 miRNAs. GEREA software was developed to perform the integrated analyses. We applied the algorithm to transcriptome data for synthetic miR-155 oligo-treated mouse CD4+ Tcells and confirmed that miR-155 is an important network regulator. The software was also tested on publicly available transcriptional profiling data for Salmonella infection, resulting in the identification of miR-125b as an important regulator. Conclusion: The causality inference-based miRNA-(transition factor)-(target gene) networks serve as a novel resource for gene expression regulation research, and GEREA is an effective and useful adjunct to the currently available methods. The regulatory networks and the algorithm implemented in the GEREA software package are available under a free academic license at http://www.thua45.cn/gerea.