A skellam model to identify differential patterns of gene expression induced by environmental signals

Integrated metabolite profiling and transcriptome analysis reveal candidate genes involved in the formation of yellow Nelumbo nucifera

As a worldwide major ornamental flower and a edible plant, lotus (Nelumbo nucifera) is also used as medicine and tea beverage. Here, transcriptome and metabolites of yellow (MLQS) and white (YGB) lotus cultivars during five key flower coloration stages were profiled. 2014 differentially expressed genes were detected with 11 carotenoids in lotus were identified for the first time. Then, regulatory networks between and within functional modules was reconstructed, and the correlation between module-metabolites and gene-metabolites was conducted within 3 core modules. 18 candidate genes related to the formation of yellow flower were screened out and a gene regulatory model for the flower color difference between MLQS and YGB were speculated as follows: The substrate competition between F3'H and F3'5'H and substrate specificity of FLS, together with differential expression of CCD4a and CCD4b were contribute to the differences in flavonoids and carotenoids. Besides, UGT73C2, UGT91C1-2 and SGTase, and regulation of UGTs by transcription factors PLATZ, MADS, NAC031, and MYB308 may also play a role in the upstream regulation. The following verification results indicated that functional differences existed in the coding sequences of NnCCD4b and promoters of NnCCD4a of MLQS and YGB. In all, this study preliminarily reveals the mechanism of yellow flower coloration in lotus and provides new ideas for the study of complex ornamental characters of other plants.

Transcriptional differentiation driving Cucumis sativus-Botrytis cinerea interactions based on the Skellam model and Bayesian networks

Robust statistical tools such as the Skellam model and Bayesian networks can capture the count properties of transcriptome sequencing data and clusters of genes among treatments, thereby improving our knowledge of gene functions and networks. In this study, we successfully implemented a model to analyze a transcriptome dataset of Cucumis sativus and Botrytis cinerea before and after their interaction. First, 4200 differentially expressed genes (DEGs) from C. sativus were clustered into 17 distinct groups, and 670 DEGs from B. cinerea were clustered into 12 groups. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied on these DEGs to assess the interactions between C. sativus and B. cinerea. In C. sativus, more DEGs were divided into terms in the molecular function and biological process domains than into cellular components, and 277 DEGs were allocated to 19 KEGG pathways. In B. cinerea, more DEGs were divided into terms in the biological process and cellular component domains than into molecular functions, and 150 DEGs were allocated to 26 KEGG pathways. In this study, we constructed networks of genes that interact with each other to screen hub genes based on a directed graphical model known as Bayesian networks. Through a detailed GO analysis, we excavated hub genes which were biologically meaningful. These results verify that availability of Skellam model and Bayesian networks in clustering gene expression data and sorting out hub genes. These models are instrumental in increasing our knowledge of gene functions and networks in plant–pathogen interaction.

Rejuvenation remodels transcriptional network to improve rhizogenesis in mature Juglans tree

Adventitious rooting of walnut species (Juglans L.) is known to be rather difficult, especially for mature trees. The adventitious root formation (ARF) capacities of mature trees can be significantly improved by rejuvenation. However, the underlying gene regulatory networks (GRNs) of rejuvenation remain largely unknown. To characterize such regulatory networks, we carried out the transcriptomic study using RNA samples of the cambia and peripheral tissues on the bottom of rejuvenated and mature walnut (Juglans hindsii × J. regia) cuttings during the ARF. The RNA sequencing data suggested that zeatin biosynthesis, energy metabolism and substance metabolism were activated by rejuvenation, whereas photosynthesis, fatty acid biosynthesis and the synthesis pathways for secondary metabolites were inhibited. The inter- and intra-module GRNs were constructed using differentially expressed genes. We identified 35 hub genes involved in five modules associated with ARF. Among these hub genes, particularly, beta-glucosidase-like (BGLs) family members involved in auxin metabolism were overexpressed at the early stage of the ARF. Furthermore, BGL12 from the cuttings of Juglans was overexpressed in Populus alba × P. glandulosa. Accelerated ARF and increased number of ARs were observed in the transgenic poplars. These results provide a high-resolution atlas of gene activity during ARF and help to uncover the regulatory modules associated with the ARF promoted by rejuvenation.

SEGN: Inferring real-time gene networks mediating phenotypic plasticity

The capacity of an organism to alter its phenotype in response to environmental perturbations changes over developmental time and is a process determined by multiple genes that are co-expressed in intricate but organized networks. Characterizing the spatiotemporal change of such gene networks can offer insight into the genomic signatures underlying organismic adaptation, but it represents a major methodological challenge. Here, we integrate the holistic view of systems biology and the interactive notion of evolutionary game theory to reconstruct so-called systems evolutionary game networks (SEGN) that can autonomously detect, track, and visualize environment-induced gene networks along the time axis. The SEGN overcomes the limitations of traditional approaches by inferring context-specific networks, encapsulating bidirectional, signed, and weighted gene-gene interactions into fully informative networks, and monitoring the process of how networks topologically alter across environmental and developmental cues. Based on the design principle of SEGN, we perform a transcriptional plasticity study by culturing Euphrates poplar, a tree that can grow in the saline desert, in saline-free and saline-stress conditions. SEGN characterize previously unknown gene co-regulation that modulates the time trajectories of the trees' response to salt stress. As a marriage of multiple disciplines, SEGN shows its potential to interpret gene interdependence, predict how transcriptional co-regulation responds to various regimes, and provides a hint for exploring the mass, energetic, or signal basis that drives various types of gene interactions.

A block mixture model to map eQTLs for gene clustering and networking

To study how genes function in a cellular and physiological process, a general procedure is to classify gene expression profiles into categories based on their similarity and reconstruct a regulatory network for functional elements. However, this procedure has not been implemented with the genetic mechanisms that underlie the organization of gene clusters and networks, despite much effort made to map expression quantitative trait loci (eQTLs) that affect the expression of individual genes. Here we address this issue by developing a computational approach that integrates gene clustering and network reconstruction with genetic mapping into a unifying framework. The approach can not only identify specific eQTLs that control how genes are clustered and organized toward biological functions, but also enable the investigation of the biological mechanisms that individual eQTLs perturb in a signaling pathway. We applied the new approach to characterize the effects of eQTLs on the structure and organization of gene clusters in Caenorhabditis elegans. This study provides the first characterization, to our knowledge, of the effects of genetic variants on the regulatory network of gene expression. The approach developed can also facilitate the genetic dissection of other dynamic processes, including development, physiology and disease progression in any organisms.

Analyzing allele specific RNA expression using mixture models

Background

Measuring allele-specific RNA expression provides valuable insights into cis-acting genetic and epigenetic regulation of gene expression. Widespread adoption of high-throughput sequencing technologies for studying RNA expression (RNA-Seq) permits measurement of allelic RNA expression imbalance (AEI) at heterozygous single nucleotide polymorphisms (SNPs) across the entire transcriptome, and this approach has become especially popular with the emergence of large databases, such as GTEx. However, the existing binomial-type methods used to model allelic expression from RNA-seq assume a strong negative correlation between reference and variant allele reads, which may not be reasonable biologically.

Results

Here we propose a new strategy for AEI analysis using RNA-seq data. Under the null hypothesis of no AEI, a group of SNPs (possibly across multiple genes) is considered comparable if their respective total sums of the allelic reads are of similar magnitude. Within each group of “comparable” SNPs, we identify SNPs with AEI signal by fitting a mixture of folded Skellam distributions to the absolute values of read differences. By applying this methodology to RNA-Seq data from human autopsy brain tissues, we identified numerous instances of moderate to strong imbalanced allelic RNA expression at heterozygous SNPs. Findings with SLC1A3 mRNA exhibiting known expression differences are discussed as examples.

Conclusion

The folded Skellam mixture model searches for SNPs with significant difference between reference and variant allele reads (adjusted for different library sizes), using information from a group of “comparable” SNPs across multiple genes. This model is particularly suitable for performing AEI analysis on genes with few heterozygous SNPs available from RNA-seq, and it can fit over-dispersed read counts without specifying the direction of the correlation between reference and variant alleles.

Electronic supplementary material

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

Modeling Expression Plasticity of Genes that Differentiate Drug-sensitive from Drug-resistant Cells to Chemotherapeutic Treatment

By measuring gene expression at an unprecedented resolution and throughput, RNA-seq has played a pivotal role in studying biological functions. Its typical application in clinical medicine is to identify the discrepancies of gene expression between two different types of cancer cells, sensitive and resistant to chemotherapeutic treatment, in a hope to predict drug response. Here we modified and used a mechanistic model to identify distinct patterns of gene expression in response of different types of breast cancer cell lines to chemotherapeutic treatment. This model was founded on a mixture likelihood of Poisson-distributed transcript read data, with each mixture component specified by the Skellam function. By estimating and comparing the amount of gene expression in each environment, the model can test how genes alter their expression in response to environment and how different genes interact with each other in the responsive process. Using the modified model, we identified the alternations of gene expression between two cell lines of breast cancer, resistant and sensitive to tamoxifen, which allows us to interpret the expression mechanism of how genes respond to metabolic differences between the two cell types. The model can have a general implication for studying the plastic pattern of gene expression across different environments measured by RNA-seq.