Search

Background Dinotefuran (CAS No. 165252–70-0), a neonicotinoid insecticide, has been used to protect various crops against invertebrate pests and has been associated with numerous negative sublethal effects on honey bees. Long noncoding RNAs (lncRNAs) play important roles in mediating various biological and pathological processes, involving transcriptional and gene regulation. The effects of dinotefuran on lncRNA expression and lncRNA function in the honey bee brain are still obscure. Results Through RNA sequencing, a comprehensive analysis of lncRNAs and mRNAs was performed following exposure to 0.01 mg/L dinotefuran for 1, 5, and 10 d. In total, 312 lncRNAs and 1341 mRNAs, 347 lncRNAs and 1458 mRNAs, and 345 lncRNAs and 1155 mRNAs were found to be differentially expressed (DE) on days 1, 5 and 10, respectively. Gene set enrichment analysis (GSEA) indicated that the dinotefuran-treated group showed enrichment in carbohydrate and protein metabolism and immune-inflammatory responses such as glycine, serine and threonine metabolism, pentose and glucuronate interconversion, and Hippo and transforming growth factor-β (TGF-β) signaling pathways. Moreover, the DE lncRNA TCONS_00086519 was shown by fluorescence in situ hybridization (FISH) to be distributed mainly in the cytoplasm, suggesting that it may serve as a competing endogenous RNA and a regulatory factor in the immune response to dinotefuran. Conclusion This study characterized the expression profile of lncRNAs upon exposure to neonicotinoid insecticides in young adult honey bees and provided a framework for further study of the role of lncRNAs in honey bee growth and the immune response.

Background Long noncoding RNAs (lncRNAs) have emerged as an important class of transcriptional regulators in cellular processes. The past decades have witnessed great progress in lncRNA studies in a variety of organisms. The codling moth (Cydia pomonella L.) is an important invasive insect in China. However, the functional impact of lncRNAs in this insect remains unclear. In this study, an atlas of codling moth lncRNAs was constructed based on publicly available RNA-seq datasets. Results In total, 9875 lncRNA transcripts encoded by 9161 loci were identified in the codling moth. As expected, the lncRNAs exhibited shorter transcript lengths, lower GC contents, and lower expression levels than protein-coding genes (PCGs). Additionally, the lncRNAs were more likely to show tissue-specific expression patterns than PCGs. Interestingly, a substantial fraction of the lncRNAs showed a testis-biased expression pattern. Additionally, conservation analysis indicated that lncRNA sequences were weakly conserved across insect species, though additional lncRNAs with homologous relationships could be identified based on synteny, suggesting that synteny could be a more reliable approach for the cross-species comparison of lncRNAs. Furthermore, the correlation analysis of lncRNAs with neighbouring PCGs indicated a stronger correlation between them, suggesting potential cis-acting roles of these lncRNAs in the regulation of gene expression. Conclusions Taken together, our work provides a valuable resource for the comparative and functional study of lncRNAs, which will facilitate the understanding of their mechanistic roles in transcriptional regulation.

Background The red sea urchin Mesocentrotus franciscanus is an ecologically important kelp forest herbivore and an economically valuable wild fishery species. To examine how M. franciscanus responds to its environment on a molecular level, differences in gene expression patterns were observed in embryos raised under combinations of two temperatures (13 °C or 17 °C) and two pCO2 levels (475 μatm or 1050 μatm). These combinations mimic various present-day conditions measured during and between upwelling events in the highly dynamic California Current System with the exception of the 17 °C and 1050 μatm combination, which does not currently occur. However, as ocean warming and acidification continues, warmer temperatures and higher pCO2 conditions are expected to increase in frequency and to occur simultaneously. The transcriptomic responses of the embryos were assessed at two developmental stages (gastrula and prism) in light of previously described plasticity in body size and thermotolerance under these temperature and pCO2 treatments. Results Although transcriptomic patterns primarily varied by developmental stage, there were pronounced differences in gene expression as a result of the treatment conditions. Temperature and pCO2 treatments led to the differential expression of genes related to the cellular stress response, transmembrane transport, metabolic processes, and the regulation of gene expression. At each developmental stage, temperature contributed significantly to the observed variance in gene expression, which was also correlated to the phenotypic attributes of the embryos. On the other hand, the transcriptomic response to pCO2 was relatively muted, particularly at the prism stage. Conclusions M. franciscanus exhibited transcriptomic plasticity under different temperatures, indicating their capacity for a molecular-level response that may facilitate red sea urchins facing ocean warming as climate change continues. In contrast, the lack of a robust transcriptomic response, in combination with observations of decreased body size, under elevated pCO2 levels suggest that this species may be negatively affected by ocean acidification. High present-day pCO2 conditions that occur due to coastal upwelling may already be influencing populations of M. franciscanus.

Background Women account for 60% of all stroke deaths and are more often permanently disabled than men, despite their higher observed stroke incidence. Considering the clinical population affected by stroke, an obvious drawback is that many pre-clinical and clinical studies only investigate young males. To improve therapeutic translation from bench to bedside, we believe that it is advantageous to include both sexes in experimental models of stroke. The aims of this study were to identify early cerebral vascular responses to ischemic stroke in females, compare the differential gene expression patterns with those seen in males, and identify potential new therapeutic targets. Results Transient middle cerebral artery occlusion (tMCAO) was used to induce stroke in both female and male rats, the middle cerebral arteries (MCAs) were isolated 3 h post reperfusion and RNA was extracted. Affymetrix whole transcriptome expression profiling was performed on female (n = 12) MCAs to reveal differentially expressed genes. In total, 1076 genes had an increased expression and 879 genes a decreased expression in the occluded MCAs as compared with the control MCAs from female rats. An enrichment of genes related to apoptosis, regulation of transcription, protein autophosphorylation, inflammation, oxidative stress, and tissue repair and recovery were seen in the occluded MCA. The high expression genes chosen for qPCR verification (Adamts4, Olr1, JunB, Fosl1, Serpine1, S1pr3, Ccl2 and Socs3) were all shown to be upregulated in the same manner in both females and males after tMCAO (p n = 23). When comparing the differentially expressed genes in female MCAs (occluded and non-occluded) with our previous findings in males after tMCAO, a total of 297 genes overlapped (all groups had 32 genes in common). Conclusions The cascades of processes initiated in the vasculature following reperfusion are complex. Dynamic gene expression alterations were observed in the occluded MCAs, and to a less pronounced degree in the non-occluded MCAs. Dysregulation of inflammation and blood-brain barrier breakdown are possible pharmacological targets. The sample of genes (

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs). Specifically, we computationally inferred mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3,942 genes and 186,132 eQTLs for 4,751 genes from 15,122,700 genetic variants for 13,476 genes, respectively. Interesting, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.Author SummaryIn the past decade, many studies have identified genetic variants associated with gene expression level (eQTLs) in different phenotypes, including tissues and diseases. Gene expression is the result of cooperation between transcriptional regulation, such as transcriptional activity, and post-transcriptional regulation, such as mRNA stability. Here, we present a computational framework that take advantage of recently developed methods to estimate mRNA stability from RNA-Seq, which is widely used to estimate gene expression, and then to identify genetic variants associated with mRNA stability (stQTLs) in lung tissue. Compared to eQTLs, we found that genetic variants that affects mRNA stability are more significantly located in the CDS and 3’UTR regions, which are known to interact with RNA-binding proteins (RBPs) or microRNAs to regulate stability. In addition, stQTLs are significantly more likely to overlap the binding sites of RBPs. We show that the six RBPs that most significantly bind to stQTLs are all known to regulate mRNA stability. This pipeline of simultaneously identifying eQTLs and stQTLs using only RNA-Seq data can provide higher resolution than traditional eQTLs study to better understand the molecular mechanisms of genetic variants on the regulation of gene expression.

Chimpanzees are the closest living relatives of humans. The divergence between human and chimpanzee ancestors dates to approximately 6,5–7,5 million years ago. Genetic features distinguishing us from chimpanzees and making us humans are still of a great interest. After divergence of their ancestor lineages, human and chimpanzee genomes underwent multiple changes including single nucleotide substitutions, deletions and duplications of DNA fragments of different size, insertion of transposable elements and chromosomal rearrangements. Human-specific single nucleotide alterations constituted 1.23% of human DNA, whereas more extended deletions and insertions cover ~ 3% of our genome. Moreover, much higher proportion is made by differential chromosomal inversions and translocations comprising several megabase-long regions or even whole chromosomes. However, despite of extensive knowledge of structural genomic changes accompanying human evolution we still cannot identify with certainty the causative genes of human identity. Most structural gene-influential changes happened at the level of expression regulation, which in turn provoked larger alterations of interactome gene regulation networks. In this review, we summarized the available information about genetic differences between humans and chimpanzees and their potential functional impacts on differential molecular, anatomical, physiological and cognitive peculiarities of these species.

Background Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses. Results By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing. Conclusions Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Background Studying plasticity in gene expression in natural systems is crucial, for predicting and managing the effects of climate change on plant species. To understand the contribution of gene expression level variations to abiotic stress compensation in a Himalaya plant (Primula sikkimensis), we carried out a transplant experiment within (Ambient), and beyond (Below Ambient and Above Ambient) the altitudinal range limit of species. We sequenced nine transcriptomes (three each from each altitudinal range condition) using Illumina sequencing technology. We compared the fitness variation of transplants among three transplant conditions. Results A large number of significantly differentially expressed genes (DEGs) between below ambient versus ambient (109) and above ambient versus ambient (85) were identified. Transcripts involved in plant growth and development were mostly up-regulated in below ambient conditions. Transcripts involved in signalling, defence, and membrane transport were mostly up-regulated in above ambient condition. Pathway analysis revealed that most of the genes involved in metabolic processes, secondary metabolism, and flavonoid biosynthesis were differentially expressed in below ambient conditions, whereas most of the genes involved in photosynthesis and plant hormone signalling were differentially expressed in above ambient conditions. In addition, we observed higher reproductive fitness in transplant individuals at below ambient condition compared to above ambient conditions; contrary to what we expect from the cold adaptive P. sikkimensis plants. Conclusions We reveal P. sikkimensis’s capacity for rapid adaptation to climate change through transcriptome variation, which may facilitate the phenotypic plasticity observed in morphological and life history traits. The genes and pathways identified provide a genetic resource for understanding the temperature stress (both the hot and cold stress) tolerance mechanism of P. sikkimensis in their natural environment.

Background The three-dimensional (3D) structure of chromatins plays significant roles during cell differentiation and development. Hi-C and other 3C-based technologies allow us to look deep into the chromatin architectures. Many studies have suggested that topologically associating domains (TAD), as the structure and functional unit, are conserved across different organs. However, our understanding about the underlying mechanism of the TAD boundary formation is still limited. Results We developed a computational method, TAD–Lactuca, to infer this structure by taking the contextual information of the epigenetic modification signals and the primary DNA sequence information on the genome. TAD–Lactuca is found stable in the case of multi-resolutions and different datasets. It could achieve high accuracy and even outperforms the state-of-art methods when the sequence patterns were incorporated. Moreover, several transcript factor binding motifs, besides the well-known CCCTC-binding factor (CTCF) motif, were found significantly enriched on the boundaries. Conclusions We provided a low cost, effective method to predict TAD boundaries. Above results suggested the incorporation of sequence features could significantly improve the performance. The sequence motif enrichment analysis indicates several gene regulation motifs around the boundaries, which is consistent with TADs may serve as the functional units of gene regulation and implies the sequence patterns would be important in chromatin folding.

Background The Metarhizium genus harbors important entomopathogenic fungi. These species have been widely explored as biological control agents, and strategies to improve the fungal virulence are under investigation. Thus, the interaction between Metarhizium species and susceptible hosts have been explored employing different methods in order to characterize putative virulence determinants. However, the impact of epigenetic modulation on the infection cycle of Metarhizium is still an open topic. Among the different epigenetic modifications, DNA methylation of cytosine bases is an important mechanism to control gene expression in several organisms. To better understand if DNA methylation can govern Metarhizium-host interactions, the genome-wide DNA methylation profile of Metarhizium anisopliae was explored in two conditions: tick mimicked infection and a saprophytic-like control. Results Using a genome wide DNA methylation profile based on bisulfite sequencing (BS-Seq), approximately 0.60% of the total cytosines were methylated in saprophytic-like condition, which was lower than the DNA methylation level (0.89%) in tick mimicked infection condition. A total of 670 mRNA genes were found to be putatively methylated, with 390 mRNA genes uniquely methylated in the tick mimicked infection condition. GO terms linked to response to stimuli, cell wall morphogenesis, cytoskeleton morphogenesis and secondary metabolism biosynthesis were over-represented in the tick mimicked infection condition, suggesting that energy metabolism is directed towards the regulation of genes associated with infection. However, recognized virulence determinants known to be expressed at distinct infection steps, such as the destruxin backbone gene and the collagen-like protein gene Mcl1, were found methylated, suggesting that a dynamic pattern of methylation could be found during the infectious process. These results were further endorsed employing RT-qPCR from cultures treated or not with the DNA methyltransferase inhibitor 5-Azacytidine. Conclusions The set of genes here analyzed focused on secondary metabolites associated genes, known to be involved in several processes, including virulence. The BS-Seq pipeline and RT-qPCR analysis employing 5-Azacytidine led to identification of methylated virulence genes in M. anisopliae. The results provided evidences that DNA methylation in M. anisopliae comprises another layer of gene expression regulation, suggesting a main role of DNA methylation regulating putative virulence determinants during M. anisopliae infection cycle.

Background The Myxococcales are well known for their predatory and developmental social processes, and for the molecular complexity of regulation of these processes. Many species within this order have unusually large genomes compared to other bacteria, and their genomes have many genes that are unique to one specific sequenced species or strain. Here, we describe RNAseq based transcriptome analysis of the FruA regulon of Myxococcus xanthus and a comparative RNAseq analysis of two Myxococcus species, M. xanthus and Myxococcus stipitatus, as they respond to starvation and begin forming fruiting bodies. Results We show that both species have large numbers of genes that are developmentally regulated, with over half the genome showing statistically significant changes in expression during development in each species. We also included a non-fruiting mutant of M. xanthus that is missing the transcriptional regulator FruA to identify the direct and indirect FruA regulon and to identify transcriptional changes that are specific to fruiting and not just the starvation response. We then identified Interpro gene ontologies and COG annotations that are significantly up- or down-regulated during development in each species. Our analyses support previous data for M. xanthus showing developmental upregulation of signal transduction genes, and downregulation of genes related to cell-cycle, translation, metabolism, and in some cases, DNA replication. Gene expression in M. stipitatus follows similar trends. Although not all specific genes show similar regulation patterns in both species, many critical developmental genes in M. xanthus have conserved expression patterns in M. stipitatus, and some groups of otherwise unstudied orthologous genes share expression patterns. Conclusions By identifying the FruA regulon and identifying genes that are similarly and uniquely regulated in two different species, this work provides a more complete picture of transcription during Myxococcus development. We also provide an R script to allow other scientists to mine our data for genes whose expression patterns match a user-selected gene of interest.

Background Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a cold-tolerant plant that is vulnerable to high temperature. This study explored the response mechanism of wucai to low temperature. In this study, wucai seedlings were treated with different temperatures, including low temperature (LT), high temperature (HT), and a control. Results According to transcriptomics analysis, the number of differentially expressed genes (DEGs) in HT and LT was 10,702 and 7267, respectively, compared with the control. The key genes associated with the physiological response of wucai to the treatments were analyzed. The Kyoto Encyclopedia of Genes and Genomes and Gene Ontology annotations indicated the importance of the photosynthesis and photosynthetic-antenna protein pathways. We found that a high-temperature environment greatly inhibited the expression of important genes in the photosynthetic pathway (BrLhc superfamily members, PsaD, PsaE, PsaD, PsaD, PsbO, PsbP, PsbQ, PsbR, PsbS, PsbW, PsbY, Psb27, and Psb28), whereas low temperature resulted in the expression of certain key genes (BrLhc superfamily members, Psa F, Psa H, Psb S, Psb H, Psb 28). In addition, the wucai seedlings exhibited better photosynthetic performance under low-temperature conditions than high-temperature conditions. Conclusions Based on the above results, we speculate that upon exposure to low temperature, the plants developed higher cold tolerance by upregulating the expression of genes related to photosynthesis. Conversely, high-temperature stress inhibited the expression of pivotal genes and weakened the self-regulating ability of the plants.

BackgroundThe thermal-manipulation (TM) during egg incubation is a cyclic exposure to hot or cold temperatures during embryogenesis that is associated to long-lasting effects on growth performance, physiology, metabolism and temperature tolerance in birds. An increase of the incubation temperature of Japanese quail eggs affected the embryonic and post-hatch survival, growth, surface temperatures and blood characteristics potentially related to thermoregulation capacities. To gain new insights in the molecular basis of TM in quails, we investigated by RNA-seq the hypothalamus transcriptome of 35 days-old male and female quails that were treated by TM or not (C, control) during embryogenesis and that were exposed (HC) or not (RT) to a 36 °C heat challenge for 7 h before sampling.ResultsFor males, 76, 27, 47 and 0 genes were differentially expressed in the CHC vs. CRT, CRT vs. TMRT, TMHC vs. TMRT and CHC vs. TMHC comparisons, respectively. For females, 17, 0, 342 and 1 genes were differentially expressed within the same respective comparisons. Inter-individual variability of gene expression response was observed particularly when comparing RT and HC female animals. The differential expression of several genes was corroborated by RT-qPCR analysis. Gene Ontology functional analysis of the differentially expressed genes showed a prevalent enrichment of terms related to cellular responses to stimuli and gene expression regulation in both sexes. Gene Ontology terms related to the membrane transport, the endoplasmic reticulum and mitochondrial functions as well as DNA metabolism and repair were also identified in specific comparisons and sexes.ConclusionsTM had little to no effect on the regulation of gene expression in the hypothalamus of 35 days-old Japanese quails. However, the consequences of TM on gene expression were revealed by the HC, with sex-specific and common functions altered. The effects of the HC on gene expression were most prominent in TM females with a ~ 20-fold increase of the number of differentially expressed genes, suggesting that TM may enhance the gene response during challenging conditions in female quail hypothalamus. TM may also promote new cellular strategies in females to help coping to the adverse conditions as illustrated by the identification of differentially expressed genes related to the mitochondrial and heat-response functions.

Background In mammals, the ovary is the essential system of female reproduction for the onset of puberty, and the abnormal puberty has negative outcomes on health. CircRNA is a non-coding RNA produced by non-canonical alternative splicing (AS). Several studies have reported that circRNA is involved in the gene regulation and plays an important role in some human diseases. However, the contribution of circRNA has received little known within the onset of puberty in ovary. Results Here, the profiles of ovarian circRNAs across pre-, in- and post-pubertal stages were established by RNA-sEq. In total, 972 circRNAs were identified, including 631 stage-specific circRNAs and 8 tissue-specific circRNAs. The biological functions of parental genes of circRNAs were enriched in steroid biosynthesis, autophagy-animal, MAPK signaling pathway, progesterone-mediated oocyte maturation and ras signaling pathway. Moreover, 5 circRNAs derived from 4 puberty-related genes (ESR1, JAK2, NF1 and ARNT) were found in this study. The A3SS events were the most alternative splicing, but IR events were likely to be arose in post-pubertal ovaries. Besides, the circRNA-miRNA-gene networks were explored for 10 differentially expressed circRNAs. Furthermore, the head-to-tail exon as well as the expressions of 10 circRNAs were validated by the divergent RT-qPCR and sanger sequencing. Conclusions In summary, the profiles of ovarian circRNAs were provided during pubertal transition in gilts, and these results provided useful information for the investigation on the onset of puberty at the ovarian-circRNAs-level in mammals.

Background Our understanding of genome regulation is ever-evolving with the continuous discovery of new modes of gene regulation, and transcriptomic studies of mammalian genomes have revealed the presence of a considerable population of non-coding RNA molecules among the transcripts expressed. One such non-coding RNA molecule is long non-coding RNA (lncRNA). However, the function of lncRNAs in gene regulation is not well understood; moreover, finding conserved lncRNA across species is a challenging task. Therefore, we propose a novel approach to identify conserved lncRNAs and functionally annotate these molecules. Results In this study, we exploited existing myogenic transcriptome data and identified conserved lncRNAs in mice and humans. We identified the lncRNAs expressing differentially between the early and later stages of muscle development. Differential expression of these lncRNAs was confirmed experimentally in cultured mouse muscle C2C12 cells. We utilized the three-dimensional architecture of the genome and identified topologically associated domains for these lncRNAs. Additionally, we correlated the expression of genes in domains for functional annotation of these trans-lncRNAs in myogenesis. Using this approach, we identified conserved lncRNAs in myogenesis and functionally annotated them. Conclusions With this novel approach, we identified the conserved lncRNAs in myogenesis in humans and mice and functionally annotated them. The method identified a large number of lncRNAs are involved in myogenesis. Further studies are required to investigate the reason for the conservation of the lncRNAs in human and mouse while their sequences are dissimilar. Our approach can be used to identify novel lncRNAs conserved in different species and functionally annotated them.

Background Hibernation is a physiological state exploited by many animals exposed to prolonged adverse environmental conditions associated with winter. Large changes in metabolism and cellular function occur, with many stress response pathways modulated to tolerate physiological challenges that might otherwise be lethal. Many studies have sought to elucidate the molecular mechanisms of mammalian hibernation, but detailed analyses are lacking in reptiles. Here we examine gene expression in the Australian central bearded dragon (Pogona vitticeps) using mRNA-seq and label-free quantitative mass spectrometry in matched brain, heart and skeletal muscle samples from animals at late hibernation, 2 days post-arousal and 2 months post-arousal. Results We identified differentially expressed genes in all tissues between hibernation and post-arousal time points; with 4264 differentially expressed genes in brain, 5340 differentially expressed genes in heart, and 5587 differentially expressed genes in skeletal muscle. Furthermore, we identified 2482 differentially expressed genes across all tissues. Proteomic analysis identified 743 proteins (58 differentially expressed) in brain, 535 (57 differentially expressed) in heart, and 337 (36 differentially expressed) in skeletal muscle. Tissue-specific analyses revealed enrichment of protective mechanisms in all tissues, including neuroprotective pathways in brain, cardiac hypertrophic processes in heart, and atrophy protective pathways in skeletal muscle. In all tissues stress response pathways were induced during hibernation, as well as evidence for gene expression regulation at transcription, translation and post-translation. Conclusions These results reveal critical stress response pathways and protective mechanisms that allow for maintenance of both tissue-specific function, and survival during hibernation in the central bearded dragon. Furthermore, we provide evidence for multiple levels of gene expression regulation during hibernation, particularly enrichment of miRNA-mediated translational repression machinery; a process that would allow for rapid and energy efficient reactivation of translation from mature mRNA molecules at arousal. This study is the first molecular investigation of its kind in a hibernating reptile, and identifies strategies not yet observed in other hibernators to cope stress associated with this remarkable state of metabolic depression. Electronic supplementary material The online version of this article (10.1186/s12864-019-5750-x) contains supplementary material, which is available to authorized users.

Background Cadmium (Cd)-containing chemicals can cause serious damage to biological systems. In animals and plants, Cd exposure can lead to metabolic disorders or death. However, for the most part the effects of Cd on specific biological processes are not known. DNA methylation is an important mechanism for the regulation of gene expression. In this study we examined the effects of Cd exposure on global DNA methylation in a living organism by whole-genome bisulfite sequencing (WGBS) using Drosophila melanogaster as model. Results A total of 71 differentially methylated regions and 63 differentially methylated genes (DMGs) were identified by WGBS. A total of 39 genes were demethylated in the Cd treatment group but not in the control group, whereas 24 showed increased methylation in the former relative to the latter. In most cases, demethylation activated gene expression: genes such as Cdc42 and Mekk1 were upregulated as a result of demethylation. There were 37 DMGs that overlapped with differentially expressed genes from the digital expression library including baz, Act5C, and ss, which are associated with development, reproduction, and energy metabolism. Conclusions DNA methylation actively regulates the physiological response to heavy metal stress in Drosophila in part via activation of apoptosis. Electronic supplementary material The online version of this article (10.1186/s12864-019-5688-z) contains supplementary material, which is available to authorized users.

Background Trait loss is a pervasive phenomenon in evolution, yet the underlying molecular causes have been identified in only a handful of cases. Most of these cases involve loss-of-function mutations in one or more trait-specific genes. However, in parasitoid insects the evolutionary loss of a metabolic trait is not associated with gene decay. Parasitoids have lost the ability to convert dietary sugars into fatty acids. Earlier research suggests that lack of lipogenesis in the parasitoid wasp Nasonia vitripennis is caused by changes in gene regulation. Results We compared transcriptomic responses to sugar-feeding in the non-lipogenic parasitoid species Nasonia vitripennis and the lipogenic Drosophila melanogaster. Both species adjusted their metabolism within 4 hours after sugar-feeding, but there were sharp differences between the expression profiles of the two species, especially in the carbohydrate and lipid metabolic pathways. Several genes coding for key enzymes in acetyl-CoA metabolism, such as malonyl-CoA decarboxylase (mcd) and HMG-CoA synthase (hmgs) differed in expression between the two species. Their combined action likely blocks lipogenesis in the parasitoid species. Network-based analysis showed connectivity of genes to be negatively correlated to the fold change of gene expression. Furthermore, genes involved in the fatty acid metabolic pathway were more connected than the set of genes of all metabolic pathways combined. Conclusion High connectivity of lipogenesis genes is indicative of pleiotropic effects and could explain the absence of gene degradation. We conclude that modification of expression levels of only a few little-connected genes, such as mcd, is sufficient to enable complete loss of lipogenesis in N. vitripennis. Electronic supplementary material The online version of this article (10.1186/s12864-019-5673-6) contains supplementary material, which is available to authorized users.

Background Paulownia withes’-broom (PaWB) disease caused by phytoplasma is a serious infectious disease for Paulownia. However, the underlying molecular pathogenesis is not fully understood. Recent studies have demonstrated that histone modifications could play a role in plant defense responses to pathogens. But there is still no available genome-wide histone modification data in non-model ligneous species infected with phytoplasma. Results Here, we provided the first genome-wide profiles of three histone marks (H3K4me3, H3K36me3 and H3K9ac) in Paulownia fortunei under phytoplasma stress by using chromatin immunoprecipitation sequencing (ChIP-Seq). We found that H3K4me3, H3K36me3 and H3K9ac were mainly enriched in the genic regions in P. fortunei with (PFI) and without (PF) phytoplasma infection. ChIP-Seq analysis revealed 1738, 986, and 2577 genes were differentially modified by H3K4me3, H3K36me3 and H3K9ac marks in PFI under phytoplasma infection, respectively. The functional analysis of these genes suggested that most of them were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction. In addition, the combinational analysis of ChIP-Seq and RNA-Seq showed that differential histone methylation and acetylation only affected a small subset of phytoplasma-responsive genes. Conclusions Taken together, this is the first report of integrated analysis of histone modifications and gene expression involved in Paulownia-phytoplasma interaction. Our results will provide the valuable resources for the mechanism studies of gene regulation in non-model plants upon pathogens attack. Electronic supplementary material The online version of this article (10.1186/s12864-019-5609-1) contains supplementary material, which is available to authorized users.

Background Multicellular organisms require precise gene regulation during ontogeny, and epigenetic modifications, such as DNA methylation and histone modification, facilitate this precise regulation. The conservative reprogramming patterns of DNA methylation in vertebrates have been well described. However, knowledge of how histone modifications are passed on from gametes to early embryos is limited, and whether histone modification reprogramming is conserved is not clear. Results We profiled H3K4me3/H3K27me3 modifications in gametes and early embryos in zebrafish and found that the patterns in gene promoter regions have been largely set to either co-occupied or active states in gametes and then passed on to early embryos. Co-occupied states are partially maintained, while active states are largely restored to nearly match the sperm’s pattern prior to zygotic genome activation (ZGA). However, repressive H3K27me3 modifications in promoter regions are largely discarded in early embryos. Prior to ZGA, patterns of genes that initialize ZGA are converted to nonrepressive states to coordinate gene expression. Moreover, promoter peaks that mark stage-specific genes are hypermethylated, and histone modifications in these regions are erased independently of DNA methylation reprogramming. Furthermore, comparative analysis revealed that the functions of co-occupied and active genes passed on from gametes are conserved in vertebrates. Gene age preferences by co-occupied and active histone modifications are also confirmed in vertebrates. Conclusions Our data provide fundamental resources for understanding H3K4me3/H3K27me3 modifications in early zebrafish embryos. The data also reveal that the reprogramming progress of histone modifications is conserved in vertebrates and coordinates with gene expression during ZGA. Electronic supplementary material The online version of this article (10.1186/s12864-019-5611-7) contains supplementary material, which is available to authorized users.

Background Cardiac cell fate specification occurs through progressive steps, and its gene expression regulation features are still being defined. There has been an increasing interest in understanding the coordination between transcription and post-transcriptional regulation during the differentiation processes. Here, we took advantage of the polysome profiling technique to isolate and high-throughput sequence ribosome-free and polysome-bound RNAs during cardiomyogenesis. Results We showed that polysome-bound RNAs exhibit the cardiomyogenic commitment gene expression and that mesoderm-to-cardiac progenitor stages are strongly regulated. Additionally, we compared ribosome-free and polysome-bound RNAs and found that the post-transcriptional regulation vastly contributes to cardiac phenotype determination, including RNA recruitment to and dissociation from ribosomes. Moreover, we found that protein synthesis is decreased in cardiomyocytes compared to human embryonic stem-cells (hESCs), possibly due to the down-regulation of translation-related genes. Conclusions Our data provided a powerful tool to investigate genes potentially controlled by post-transcriptional mechanisms during the cardiac differentiation of hESC. This work could prospect fundamental tools to develop new therapy and research approaches. Electronic supplementary material The online version of this article (10.1186/s12864-019-5550-3) contains supplementary material, which is available to authorized users.

Background Fruit peel colour is an important agronomic trait for fruit quality. Cytosine methylation plays an important role in gene regulation. Although the DNA methylation level of a single gene is important to affect the phenotype of mutation, there are large unknown of difference of the DNA methylation in plant and its mutants. Results Using bisulfite sequencing (BS-Seq) and RNA-sequencing (RNA-Seq), we analysed three deep-red-skinned apple (Malus × domestica) mutants (Yanfu 3, YF3; Yanfu 8, YF8; Shannonghong, SNH) and their lighter-skinned parents (Nagafu 2, NF2; Yanfu 3, YF3; Ralls, RL) to explore the different changes in methylation patterns associated with anthocyanin concentrations. We identified 13,405, 13,384, and 10,925 differentially methylated regions (DMRs) and 1987, 956, and 1180 differentially expressed genes (DEGs) in the NF2/YF3, YF3/YF8, and RL/SNH comparisons, respectively. And we found two DMR-associated DEGs involved in the anthocyanin pathway: ANS (MD06G1071600) and F3H (MD05G1074200). These genes exhibited upregulated expression in apple mutants, and differences were observed in the methylation patterns of their promoters. These results suggested that both the regulatory and structural genes may be modified by DNA methylation in the anthocyanin pathway. However, the methylation of structural genes was not the primary reason for expression-level changes. The expression of structural genes may be synergistically regulated by transcription factors and methylation changes. Additionally, the expression of the transcription factor gene MYB114 (MD17G1261100) was upregulated in the deep-red-skinned apple. Conclusion Through the analysis of global methylation and transcription, we did not find the correlation between gene expression and the DNA methylation. However, we observed that the upregulated expression of ANS (MD06G1071600) and F3H (MD05G1074200) in apple mutants results in increased anthocyanin contents. Moreover, MYB114 (MD17G1261100) is likely another regulatory gene involved in apple coloration. Our data provided a new understanding about the differences in formation of apple colour mutants. Electronic supplementary material The online version of this article (10.1186/s12864-019-5499-2) contains supplementary material, which is available to authorized users.

Background Eggshell formation takes place in the shell gland of the oviduct of laying hens. The eggshell is rich in calcium and various glycoproteins synthesised in the shell gland. Although studies have identified genes involved in eggshell formation, little is known about the regulation of genes in the shell gland particularly in a temporal manner. The current study investigated the global gene expression profile of the shell gland of laying hens at different time-points of eggshell formation using RNA-Sequencing (RNA-Seq) analysis. Results Gene expression profiles of the shell gland tissue at 5 and 15 h time-points were clearly distinct from each other. Out of the 14,334 genes assessed for differential expression in the shell gland tissue, 278 genes were significantly down-regulated (log2 fold change > 1.5; FDR

Background Alfalfa (Medicago sativa L.) is a forage legume with significant agricultural value worldwide. MicroRNAs (miRNAs) are key components of post-transcriptional gene regulation and essentially regulate many aspects of plant growth and development. Although miRNAs were reported in alfalfa, their expression profiles in different tissues and the discovery of novel miRNAs as well as their targets have not been described in this plant species. Results To identify tissue-specific miRNA profiles in whole plants, shoots and roots of three different alfalfa genotypes (Altet-4, NECS-141and NF08ALF06) were used. Small RNA libraries were generated and sequenced using a high-throughput sequencing platform. Analysis of these libraries enabled identification of100 miRNA families; 21 of them belong to the highly conserved families while the remaining 79 families are conserved at the minimum between M. sativa and the model legume and close relative, M. truncatula. The profiles of the six abundantly expressed miRNA families (miR156, miR159, miR166, miR319, miR396 and miR398) were relatively similar between the whole plants, roots and shoots of these three alfalfa genotypes. In contrast, robust differences between shoots and roots for miR160 and miR408 levels were evident, and their expression was more abundant in the shoots. Additionally, 17 novel miRNAs were identified and the relative abundance of some of these differed between tissue types. Further, the generation and analysis of degradome libraries from the three alfalfa genotypes enabled confirmation of 69 genes as targets for 31 miRNA families in alfalfa. Conclusions The miRNA profiles revealed both similarities and differences in the expression profiles between tissues within a genotype as well as between the genotypes. Among the highly conserved miRNA families, miR166 was the most abundantly expressed in almost all tissues from the three genotypes. The identification of conserved and novel miRNAs as well as their targets in different tissues of multiple genotypes increased our understanding of miRNA-mediated gene regulation in alfalfa and could provide valuable insights for practical research and plant improvement applications in alfalfa and related legume species.

Background The regulatory landscape of a gene locus often consists of several functionally redundant enhancers establishing phenotypic robustness and evolutionary stability of its regulatory program. However, it is unclear what mechanisms are employed by redundant enhancers to cooperatively orchestrate gene expression. Results By comparing redundant enhancers to single enhancers (enhancers present in a single copy in a gene locus), we observed that the DNA sequence encryption differs between these two classes of enhancers, suggesting a difference in their regulatory mechanisms. Initiator enhancers, which are a subset of redundant enhancers and show similar sequence encryption to single enhancers, differ from the rest of redundant enhancers in their sequence encryption, evolutionary conservation and proximity to target genes. Genes hosting initiator enhancers in their loci feature elevated levels of expression. Initiator enhancers show a high level of 3D chromatin contacts with both transcription start sites and regular enhancers, suggesting their roles as primary activators and intermediate catalysts of gene expression, through which the regulatory signals of redundant enhancers are propagated to the target genes. In addition, GWAS and eQTLs variants are significantly enriched in initiator enhancers compared to redundant enhancers, suggesting a key functional role these sequences play in gene regulation. Conclusions The specific characteristics and widespread abundance of initiator enhancers advocate for a possible universal hierarchical mechanism of tissue-specific gene regulation involving multiple redundant enhancers acting through initiator enhancers. Electronic supplementary material The online version of this article (10.1186/s12864-018-5335-0) contains supplementary material, which is available to authorized users.

Background The epithelium is the first line of defense against pathogens. Notably the epithelial cells lining the respiratory track are crucial in sensing airborne microbes and mounting an effective immune response via the expression of target genes such as cytokines and chemokines. Gene expression regulation following microbial recognition is partly regulated by chromatin re-organization and has been described in immune cells but data from epithelial cells is not as detailed. Here, we report genome-wide changes of the H3K27ac mark, characteristic of activated enhancers and promoters, after stimulation of nasopharyngeal epithelial cells with the bacterial endotoxin Lipopolysaccharide (LPS). Results In this study, we have identified 626 regions where the H3K27ac mark showed reproducible increase following LPS induction in epithelial cells. This indicated that sensing of LPS led to opening of the chromatin in our system. Moreover, this phenomenon seemed to happen extensively at enhancers regions and we could observe instances of Super-enhancer formation. As expected, LPS-increased H3K27ac regions were found in the vicinity of genes relevant for LPS response and these changes correlated with up-regulation of their expression. In addition, we found the induction of H3K27ac mark to overlap with the binding of one of the NF-kB members and key regulator of the innate immune response, RELA, following LPS sensing. Indeed, inhibiting the NF-kB pathway abolished the deposition of H3K27ac at the TNF locus, a target of RELA, suggesting that these two phenomena are associated. Conclusions Enhancers’ selection and activation following microbial or inflammatory stimuli has been described previously and shown to be mediated via the NF-kB pathway. Here, we demonstrate that this is also likely to occur in the case of LPS-sensing by nasopharyngeal epithelial cells as well. In addition to validating previous findings, we generated a valuable data set relevant to the host immune response to epithelial cell colonizing or infecting pathogens. Electronic supplementary material The online version of this article (10.1186/s12864-018-5295-4) contains supplementary material, which is available to authorized users.

Backgrounds Macaca fascicularis (M. fascicularis) is a primate model organism that played important role in studying human health. It is vital to better understand the similarity and differences of gene regulation between M. fascicularis and human. Current comparative study of gene regulation between the two species are limited by low quality of gene annotation and lack of regulatory element data on M. fascicularis genome. Results In this study, we improved the M. fascicularis gene annotation with 57 gene expression data from multiple tissues and, more importantly, a manual curation procedure. The new annotation enabled us to map gene expression and identify gene location more accurately. Conclusions Comparing with human gene expression data from the same cell types, we characterized the evolution of expression patterns of homologous genes. Electronic supplementary material The online version of this article (10.1186/s12864-018-5183-y) contains supplementary material, which is available to authorized users.

Background It is still challenging to predict interacting enhancer-promoter pairs (IEPs), partially because of our limited understanding of their characteristics. To understand IEPs better, here we studied the IEPs in nine cell lines and nine primary cell types. Results By measuring the bipartite clustering coefficient of the graphs constructed from these experimentally supported IEPs, we observed that one enhancer is likely to interact with either none or all of the target genes of another enhancer. This observation implies that enhancers form clusters, and every enhancer in the same cluster synchronously interact with almost every member of a set of genes and only this set of genes. We perceived that an enhancer can be up to two megabase pairs away from other enhancers in the same cluster. We also noticed that although a fraction of these clusters of enhancers do overlap with super-enhancers, the majority of the enhancer clusters are different from the known super-enhancers. Conclusions Our study showed a new characteristic of IEPs, which may shed new light on distal gene regulation and the identification of IEPs.

BackgroundBovine Respiratory Syncytial Virus (BRSV) is a cause of Bovine Respiratory Disease (BRD). DNA-based biomarkers contributing to BRD resistance are potentially present in non-protein-coding regulatory regions of the genome, which can be determined using ATAC-Seq. The objectives of this study were to: (i) identify regions of open chromatin in DNA extracted from bronchial lymph nodes (BLN) of healthy dairy calves experimentally challenged with BRSV and compare them with those from non-challenged healthy control calves, (ii) elucidate the chromatin regions that were differentially or uniquely open in the BRSV challenged relative to control calves, and (iii) compare the genes found in regions proximal to the differentially open regions to the genes previously found to be differentially expressed in the BLN in response to BRSV and to previously identified BRD susceptibility loci. This was achieved by challenging clinically healthy Holstein-Friesian calves (mean age 143 ± 14 days) with either BRSV inoculum (n = 12) or with sterile phosphate buffered saline (PBS) (n = 6) and preparing and sequencing ATAC-Seq libraries from fresh BLN tissues.ResultsUsing Diffbind, 9,144 and 5,096 differentially accessible regions (P ConclusionsThe identified open chromatin regions are likely to be involved in the regulatory response of gene transcription induced by infection with BRSV. Consequently, they may contain variants which impact resistance to BRD that could be used in breeding programmes to select healthier, more robust cattle.

Background Yaks that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from cattle and have adapted well to the extreme conditions on the plateau. However, the mechanisms used by these animals for the regulation of gene expression at high altitude are not fully understood. Results Here, we sequenced nine lung transcriptomes of yaks at altitudes of 3400, 4200 and 5000 m, and low-altitude Zaosheng cattle, which is a closely related species, served as controls. The analysis identified 21,764 mRNAs, 1377 circRNAs and 1209 miRNAs. By comparing yaks and cattle, 4975 mRNAs, 252 circRNAs and 75 miRNAs were identified differentially expressed. By comparing yaks at different altitudes, we identified 756 mRNAs, 64 circRNAs and 83 miRNAs that were differentially expressed (fold change ≥2 and P-value Conclusions This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in yaks and might contribute to improvements in the understanding and prevention of hypoxia-related diseases.

Background The ongoing COVID-19 outbreak has caused devastating mortality and posed a significant threat to public health worldwide. Despite the severity of this illness and 2.3 million worldwide deaths, the disease mechanism is mostly unknown. Previous studies that characterized differential gene expression due to SARS-CoV-2 infection lacked robust validation. Although vaccines are now available, effective treatment options are still out of reach. Results To characterize the transcriptional activity of SARS-CoV-2 infection, a gene signature consisting of 25 genes was generated using a publicly available RNA-Sequencing (RNA-Seq) dataset of cultured cells infected with SARS-CoV-2. The signature estimated infection level accurately in bronchoalveolar lavage fluid (BALF) cells and peripheral blood mononuclear cells (PBMCs) from healthy and infected patients (mean 0.001 vs. 0.958; P TNFAIP3, PPP1R15A, NFKBIA, and IFIT2 had shown bimodal gene expression in various immune cells from severely infected patients compared to healthy or moderate infection cases. Finally, this signature was assessed using the publicly available ConnectivityMap database to identify potential disease mechanisms and drug repurposing candidates. Pharmacological classes of tricyclic antidepressants, SRC-inhibitors, HDAC inhibitors, MEK inhibitors, and drugs such as atorvastatin, ibuprofen, and ketoconazole showed strong negative associations (connectivity score Conclusions Thus, using the 25-gene SARS-CoV-2 infection signature, the SARS-CoV-2 infection status was captured in BALF cells, PBMCs and postmortem lung biopsies. In addition, candidate SARS-CoV-2 therapies with known safety profiles were identified. The signature genes could potentially also be used to characterize the COVID-19 disease severity in patients’ expression profiles of BALF cells.

Background Long non-coding RNAs (lncRNAs) play crucial roles in gene regulation at the transcriptional and post-transcriptional levels. LncRNAs are belonging to a large class of transcripts with ≥200 nt in length which do not code for proteins, have been widely investigated in various physiological and pathological contexts by high-throughput sequencing techniques and bioinformatics analysis. However, little is known about the regulatory mechanisms by which lncRNAs regulate genes that are associated with Enterotoxigenic Escherichia coli F4 fimbriae (ETEC-F4ac) adhesion phenotype in small intestine epithelial cells of Large White piglets. To address this, we used RNA sequencing to profile lncRNAs and mRNAs of small intestine epithelial cells in Large White piglets differing in their ETEC-F4 adhesion phenotypes and ITGB5 genotypes. Eight male piglets were used in this study and were divided into two groups on the basis of their adhesion phenotype and ITGB5 genotypes, a candidate gene for F4ac receptor. Non-adhesive group (n = 4) with CC genotype and adhesive group (n = 4) with TT genotype. Results In total, 78 differentially expressed lncRNAs (DE-lncRNA) and 223 differentially expressed mRNAs (log2 |FC| > 1, P cis- and trans-regulatory target genes of DE-lncRNAs were identified, then interaction networks of lncRNAs and their cis- and trans-target differentially expressed genes (DEGs) were constructed separately. A total of 194 cis-targets were involved in the lncRNAs-cis genes interaction network and 61 trans-targets, were involved in lncRNA-trans gene interaction network that we constructed. We determined that cis-target genes were involved in alcoholism, systemic lupus erythematosus, viral carcinogenesis and malaria. Whereas trans-target DEGs were engaged in three important pathways related to the ETEC-F4 adhesion phenotype namely cGMP-PKG signaling pathway, focal adhesion, and adherens junction. The trans-target DEGs which directly involved in these pathways are KCNMB1 in cGMP-PKG signaling pathway, GRB2 in focal adhesion pathway and ACTN4 in focal adhesion and adherens junction pathways. Conclusion The findings of the current study provides an insight into biological functions and epigenetic regulatory mechanism of lncRNAs on porcine small intestine epithelial cells adhesion to ETEC-F4-ac and piglets’ diarrhea susceptibility/resistance.

Background Associations of tree roots with diverse symbiotic mycorrhizal fungi have distinct effects on whole plant functioning. An untested explanation might be that such effect variability is associated with distinct impacts of different fungi on gene expression in local and distant plant organs. Using a large scale transcriptome sequencing approach, we compared the impact of three ectomycorrhizal (EMF) and one orchid mycorrhizal fungi (OMF) on gene regulation in colonized roots (local), non-colonized roots (short distance) and leaves (long distance) of the Quercus robur clone DF159 with reference to the recently published oak genome. Since different mycorrhizal fungi form symbiosis in a different time span and variable extents of apposition structure development, we sampled inoculated but non-mycorrhizal plants, for which however markedly symbiotic effects have been reported. Local root colonization by the fungi was assessed by fungal transcript analysis. Results The EMF induced marked and species specific effects on plant development in the analysed association stage, but the OMF did not. At local level, a common set of plant differentially expressed genes (DEG) was identified with similar patterns of responses to the three EMF, but not to the OMF. Most of these core DEG were down-regulated and correspond to already described but also new functions related to establishment of EMF symbiosis. Analysis of the fungal transcripts of two EMF in highly colonized roots also revealed onset of a symbiosis establishment. In contrast, in the OMF, the DEG were mainly related to plant defence. Already at short distances, high specificities in transcriptomic responses to the four fungi were detected, which were further enhanced at long distance in leaves, where almost no common DEG were found between the treatments. Notably, no correlation between phylogeny of the EMF and gene expression patterns was observed. Conclusions Use of clonal oaks allowed us to identify a core transcriptional program in roots colonized by three different EMF, supporting the existence of a common EMF symbiotic pathway. Conversely, the specific responses in non-colonized organs were more closely related to the specific impacts of the different of EMF on plant performance.

BackgroundDNA methylation is a major epigenetic modification involved in regulating gene expression. The effects of DNA methylation on gene expression differ by genomic location and vary across kingdoms, species and environmental conditions. To identify the functional regulatory roles of DNA methylation, the correlation between DNA methylation changes and alterations in gene expression is crucial. With the advance of next-generation sequencing, genome-wide methylation and gene expression profiling have become feasible. Current bioinformatics tools for investigating such correlation are designed to the assessment of DNA methylation at CG sites. The correlation of non-CG methylation and gene expression is very limited. Some bioinformatics databases allow correlation analysis, but they are limited to specific genomes such as that of humans and do not allow user-provided data.ResultsHere, we developed a bioinformatics web tool, MethGET (Methylation and Gene Expression Teller), that is specialized to analyse the association between genome-wide DNA methylation and gene expression. MethGET is the first web tool to which users can supply their own data from any genome. It is also the tool that correlates gene expression with CG, CHG, and CHH methylation based on whole-genome bisulfite sequencing data. MethGET not only reveals the correlation within an individual sample (single-methylome) but also performs comparisons between two groups of samples (multiple-methylomes). For single-methylome analyses, MethGET provides Pearson correlations and ordinal associations to investigate the relationship between DNA methylation and gene expression. It also groups genes with different gene expression levels to view the methylation distribution at specific genomic regions. Multiple-methylome analyses include comparative analyses and heatmap representations between two groups. These functions enable the detailed investigation of the role of DNA methylation in gene regulation. Additionally, we applied MethGET to rice regeneration data and discovered that CHH methylation in the gene body region may play a role in the tissue culture process, which demonstrates the capability of MethGET for use in epigenomic research.ConclusionsMethGET is a Python software that correlates DNA methylation and gene expression. Its web interface is publicly available athttps://paoyang.ipmb.sinica.edu.tw/Software.html. The stand-alone version and source codes are available on GitHub athttps://github.com/Jason-Teng/MethGET.

Background Recent sequencing projects on early-diverging metazoans such as cnidarians, have unveiled a rich innate immunity gene repertoire; however, little is known about immunity gene regulation in the host’s early response against marine bacterial pathogens over time. Here, we used RNA-seq on the sea anemone Exaiptasia pallida (Ep) strain CC7 as a model to depict the innate immune response during the onset of infection with the marine pathogenic bacteria Vibrio parahaemolyticus (Vp) clinical strain O3:K6, and lipopolysaccharides (LPS) exposure. Pairwise and time series analyses identified the genes responsive to infection as well as the kinetics of innate immune genes over time. Comparisons between the responses to live Vp and purified LPS was then performed. Results Gene expression and functional analyses detected hundreds to thousands of genes responsive to the Vp infection after 1, 3, 6 and 12 h, including a few shared with the response to LPS. Our results bring to light the first indications that non-canonical cytoplasmic pattern recognition receptors (PRRs) such as NOD-like and RIG-I-like receptor homologs take part in the immune response of Ep. Over-expression of several members of the lectin-complement pathways in parallel with novel transmembrane and Ig containing ficolins (CniFLs) suggest an active defense against the pathogen. Although lacking typical Toll-like receptors (TLRs), Ep activates a TLR-like pathway including the up-regulation of MyD88, TRAF6, NF-κB and AP-1 genes, which are not induced under LPS treatment and therefore suggest an alternative ligand-to-PRR trigger. Two cytokine-dependent pathways involving Tumor necrosis factor receptors (TNFRs) and several other potential downstream signaling genes likely lead to inflammation and/or apoptosis. Finally, both the extrinsic and intrinsic apoptotic pathways were strongly supported by over-expression of effector and executioner genes. Conclusions To our knowledge, this pioneering study is first to follow the kinetics of the innate immune response in a cnidarian during the onset of infection with a bacterial pathogen. Overall, our findings reveal the involvement of both novel immune gene candidates such as NLRs, RLRs and CniFLs, and previously identified TLR-like and apoptotic pathways in anthozoan innate immunity with a large amount of transcript-level evidence.

Background The circadian clock drives endogenous 24-h rhythms that allow organisms to adapt and prepare for predictable and repeated changes in their environment throughout the day-night (diurnal) cycle. Many components of the circadian clock in Arabidopsis thaliana have been functionally characterized, but comparatively little is known about circadian clocks in grass species including major crops like maize and sorghum. Results Comparative research based on protein homology and diurnal gene expression patterns suggests the function of some predicted clock components in grasses is conserved with their Arabidopsis counterparts, while others have diverged in function. Our analysis of diurnal gene expression in three panicoid grasses sorghum, maize, and foxtail millet revealed conserved and divergent evolution of expression for core circadian clock genes and for the overall transcriptome. We find that several classes of core circadian clock genes in these grasses differ in copy number compared to Arabidopsis, but mostly exhibit conservation of both protein sequence and diurnal expression pattern with the notable exception of maize paralogous genes. We predict conserved cis-regulatory motifs shared between maize, sorghum, and foxtail millet through identification of diurnal co-expression clusters for a subset of 27,196 orthologous syntenic genes. In this analysis, a Cochran–Mantel–Haenszel based method to control for background variation identified significant enrichment for both expected and novel 6–8 nucleotide motifs in the promoter regions of genes with shared diurnal regulation predicted to function in common physiological activities. Conclusions This study illustrates the divergence and conservation of circadian clocks and diurnal regulatory networks across syntenic orthologous genes in panacoid grass species. Further, conserved local regulatory sequences contribute to the architecture of these diurnal regulatory networks that produce conserved patterns of diurnal gene expression.

Background Circular RNA is a type of non-coding RNA, which has a circular structure. Many circular RNAs are stable and contain exons, but are not translated into proteins. Circular RNA has important functions in gene regulation and plays an important role in some human diseases. Several biological methods, such as RNase R treatment, have been developed to identify circular RNA. Multiple bioinformatics tools have also been developed for circular RNA detection with high-throughput sequence data. Results In this paper, we present circDBG, a new method for circular RNA detection with de Bruijn graph. We conduct various experiments to evaluate the performance of CircDBG based on both simulated and real data. Our results show that CircDBG finds more reliable circRNA with low bias, has more efficiency in running time, and performs better in balancing accuracy and sensitivity than existing methods. As a byproduct, we also introduce a new method to classify circular RNAs based on reads alignment. Finally, we report a potential chimeric circular RNA that is found by CircDBG based on real sequence data. CircDBG can be downloaded from https://github.com/lxwgcool/CircDBG. Conclusions We develop a new method called CircDBG for circular RNA detection, which is based on de Bruijn graph. We conduct extensive experiments and demonstrate CircDBG outperforms existing tools, especially in saving running time, reducing bias and improving capability of balancing accuracy and sensitivity. We also introduce a new method to classify circular RNAs and report a potential case of chimeric circular RNA.

BackgroundGene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.ResultsWe employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression inSaccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1andRAP1-HAP1-MSN4).ConclusionOur findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

Background Intron retention (IR) has been traditionally overlooked as ‘noise’ and received negligible attention in the field of gene expression analysis. In recent years, IR has become an emerging field for interrogating transcriptomes because it has been recognized to carry out important biological functions such as gene expression regulation and it has been found to be associated with complex diseases such as cancers. However, methods for detecting IR today are limited. Thus, there is a need to develop novel methods to improve IR detection. Results Here we present iREAD (intron REtention Analysis and Detector), a tool to detect IR events genome-wide from high-throughput RNA-seq data. The command line interface for iREAD is implemented in Python. iREAD takes as input a BAM file, representing the transcriptome, and a text file containing the intron coordinates of a genome. It then 1) counts all reads that overlap intron regions, 2) detects IR events by analyzing the features of reads such as depth and distribution patterns, and 3) outputs a list of retained introns into a tab-delimited text file. iREAD provides significant added value in detecting IR compared with output from IRFinder with a higher AUC on all datasets tested. Both methods showed low false positive rates and high false negative rates in different regimes, indicating that use together is generally beneficial. The output from iREAD can be directly used for further exploratory analysis such as differential intron expression and functional enrichment. The software is freely available at https://github.com/genemine/iread. Conclusion Being complementary to existing tools, iREAD provides a new and generic tool to interrogate poly-A enriched transcriptomic data of intron regions. Intron retention analysis provides a complementary approach for understanding transcriptome.

Background Tuberculosis remains one of the leading causes of morbidity and mortality worldwide. Therefore, understanding the pathophysiology of Mycobacterium tuberculosis is imperative for developing new drugs. Post-transcriptional regulation plays a significant role in microbial adaptation to different growth conditions. While the proteins associated with gene expression regulation have been extensively studied in the pathogenic strain M. tuberculosis H37Rv, post-transcriptional regulation involving small RNAs (sRNAs) remains poorly understood. Results We developed a novel moving-window based approach to detect sRNA expression using RNA-Seq data. Overlaying ChIP-seq data of RNAP (RNA Polymerase) and NusA suggest that these putative sRNA coding regions are significantly bound by the transcription machinery. Besides capturing many experimentally validated sRNAs, we observe the context-dependent expression of novel sRNAs in the intergenic regions of M. tuberculosis genome. For example, ncRv11806 shows expression only in the stationary phase, suggesting its role in mycobacterial latency which is a key attribute to long term pathogenicity. Also, ncRv11875C showed expression in the iron-limited condition, which is prevalent inside the macrophages of the host cells. Conclusion The systems level analysis of sRNA highlights the condition-specific expression of sRNAs which might enable the pathogen survival by rewiring regulatory circuits.

Background Luman is a member of CREB3 (cAMP responsive element-binding) subfamily of the basic leucine-zipper (bZIP) transcription factors. It may play an important regulatory role during the decidualization process since Luman was highly expressed in the decidual cells. However, the exact molecular mechanisms of how Luman regulating decidualization is unknown. Results Using an in vitro model, we prove that Luman knockdown significantly affects the decidualization process of mice endometrial stromal cells (ESCs) as the expression of two decidual markers PRL8a2 and PRL3c1 were repressed. We employed massively parallel RNA sequencing (RNA-Seq) to understand the changes in the transcriptional landscape associated with knockdown of Luman in ESCs during in vitro decidualization. We found significant dysregulation of genes related to protein processing in the endoplasmic reticulum (ER). Several genes involved in decidualization including bone morphogenetic proteins (e.g. BMP1, BMP4, BMP8A, BMP2, and BMP8B), growth factor-related genes (e.g. VEGFB, FGF10, and FGFR2), and transcription factors (IF4E, IF4A2, WNT4, WNT9A, ETS1, NOTCH1, IRX1, IDB1, IDB2, and IDB3), show altered expression. We also found that the knockdown of Luman is associated with increased expression of cell cycle-related genes including cycA1, cycB1, cycB2, CDK1, CDK2, and PLPK1, which resulted in an increased proportion of ESCs in the G1 phase. Differentially expressed genes (DEGs) were highly enriched on ECM-receptor interaction signaling, endoplasmic reticulum protein processing, focal adhesion, and PI3K-Akt signaling pathways. Conclusions Luman knockdown results in widespread gene dysregulation during decidualization of ESCs. Genes involved in protein processing in ER, bone morphogenetic protein, growth factor, and cell cycle progression were identified as particularly important for explaining the decidual deficiency observed in this in vitro model. Therefore, this study provides clues as to the underlying mechanisms that may expand our understanding of gene regulation during decidualization.

Background Large difference in cerebrum size exist between avian species and populations of the same species and is believed to reflect differences in processing power, i.e. in the speed and efficiency of processing information in this brain region. During domestication chickens developed a larger cerebrum compared to their wild progenitor, the Red jungle fowl. The underlying mechanisms that control cerebrum size and the extent to which genetic regulation is similar across brain regions is not well understood. In this study, we combine measurement of cerebrum size with genome-wide genetical genomics analysis to identify the genetic architecture of the cerebrum, as well as compare the regulation of gene expression in this brain region with gene expression in other regions of the brain (the hypothalamus) and somatic tissue (liver). Results We identify one candidate gene that putatively regulates cerebrum size (MTF2) as well as a large number of eQTL that regulate the transcriptome in cerebrum tissue, with the majority of these eQTL being trans-acting. The overall regulation of gene expression variation in the cerebrum was markedly different to the hypothalamus, with relatively few eQTL in common. In comparison, the cerebrum tissue shared more eQTL with a distant tissue (liver) than with a neighboring tissue (hypothalamus). Conclusion The candidate gene for cerebrum size (MTF2) has previously been linked to brain development making it a good candidate for further investigation as a regulator of inter-population variation in cerebrum size. The lack of shared eQTL between the two brain regions implies that genetic regulation of gene expression appears to be relatively independent between the two brain regions and suggest that coevolution between these two brain regions might be more functionally driven than developmental. These findings have relevance for current brain size evolution theories.

Background Plant NADPH oxidase (NOX), also known as respiratory burst oxidase homolog (rboh), encoded by the rboh gene, is a key enzyme in the reactive oxygen species (ROS) metabolic network. It catalyzes the formation of the superoxide anion (O2•−), a type of ROS. In recent years, various studies had shown that members of the plant rboh gene family were involved in plant growth and developmental processes as well as in biotic and abiotic stress responses, but little is known about its functional role in upland cotton. Results In the present study, 26 putative Ghrboh genes were identified and characterized. They were phylogenetically classified into six subfamilies and distributed at different densities across 18 of the 26 chromosomes or scaffolds. Their exon-intron structures, conserved domains, synteny and collinearity, gene family evolution, regulation mediated by cis-acting elements and microRNAs (miRNAs) were predicted and analyzed. Additionally, expression profiles of Ghrboh gene family were analyzed in different tissues/organs and at different developmental stages and under different abiotic stresses, using RNA-Seq data and real-time PCR. These profiling studies indicated that the Ghrboh genes exhibited temporal and spatial specificity with respect to expression, and might play important roles in cotton development and in stress tolerance through modulating NOX-dependent ROS induction and other signaling pathways. Conclusions This comprehensive analysis of the characteristics of the Ghrboh gene family determined features such as sequence, synteny and collinearity, phylogenetic and evolutionary relationship, expression patterns, and cis-element- and miRNA-mediated regulation of gene expression. Our results will provide valuable information to help with further gene cloning, evolutionary analysis, and biological function analysis of cotton rbohs.

Background Small RNAs (sRNAs) are noncoding molecules that regulate different cellular activities in several bacteria. The role of sRNAs in gene expression regulation is poorly characterized in the etiological agent of porcine enzootic pneumonia Mycoplasma hyopneumoniae. We performed a global analysis of the sRNAs, sRNA target genes and regulatory elements previously identified in their genome and analyzed the expression of some sRNAs and their target genes by quantitative RT-PCR (qPCR) in three different culture conditions. Results Seven of the 145 sRNA target genes are organized as monocistronic genes (mCs) while the other 138 sRNA target genes are organized into transcriptional units (TU). The identification of transcriptional regulatory elements (promoter motif, DNA repeat sequence or intrinsic terminator) was verified in 116 of the 145 sRNA target genes. Moreover, the 29 sRNA target genes without regulatory elements revealed the presence of at least one regulatory element in the boundaries of the TU or in other internal genes of the TU. We verified that 16 sRNAs showed differential expression, seven in heat shock condition and 14 in oxidative stress condition. Analysis of the differential expression of the sRNA target genes showed that the tested sRNAs possibly regulate gene expression. The sRNA target genes were up- or down-regulated possibly in response to sRNA only under oxidative stress condition. Moreover, the sRNA target genes are involved in diverse processes of the cell, some of which could be linked to transcription processes and cell homeostasis. Conclusion Our results indicate that bacterial sRNAs could regulate a number of targets with various outcomes, and different correlations between the levels of sRNA transcripts and their target gene mRNAs were found, which suggest that the regulation of gene expression via sRNAs may play an important role in mycoplasma. Electronic supplementary material The online version of this article (10.1186/s12864-018-5136-5) contains supplementary material, which is available to authorized users.