DeGNServer: deciphering genome-scale gene networks through high performance reverse engineering analysis

Acid phosphatase involved in phosphate homeostasis in Brassica napus and the functional analysis of BnaPAP10s

Purple acid phosphatases (PAPs) are involved in activating the rhizosphere's organic phosphorus (P) and promoting P recycling during plant development, especially under the long-term P deficiency conditions in acid soil. However, the function of BnaPAPs in response to P deficiency stress in Brassica napus has rarely been explored. In this study, we found that the acid phosphatase activities (APA) of rapeseed shoot and root increased under P deficienct conditions. Genome-wide identification found that 82 PAP genes were unevenly distributed on 19 chromosomes in B. napus, which could be divided into eight subfamilies. The segmental duplication events were the main driving force for expansion during evolution, and the gene structures and conserved motifs of most members within the same subfamily were highly conservative. Moreover, the expression levels of 37 and 23 different expressed genes were induced by low P in leaf and root, respectively. BnaA09.PAP10a and BnaC09.PAP10a were identified as candidate genes via interaction networks. Significantly, both BnaPAP10a overexpression lines significantly increased root-related APA and total phosphate concentration under P deficiency and ATP supply conditions, thereby improving plant growth and root length. In summary, our results provided a valuable foundation for further study of BnaPAP functions.

Physiological Responses of Cigar Tobacco Crop to Nitrogen Deficiency and Genome-Wide Characterization of the NtNPF Family Genes

Tobacco prefers nitrate as a nitrogen (N) source. However, little is known about the molecular components responsible for nitrate uptake and the physiological responses of cigar tobacco to N deficiency. In this study, a total of 117 nitrate transporter 1 (NRT1) and peptide transporter (PTR) family (NPF) genes were comprehensively identified and systematically characterized in the whole tobacco genome. The NtNPF members showed significant genetic diversity within and across subfamilies but showed conservation between subfamilies. The NtNPF genes are dispersed unevenly across the chromosomes. The phylogenetic analysis revealed that eight subfamilies of NtNPF genes are tightly grouped with their orthologues in Arabidopsis. The promoter regions of the NtNPF genes had extensive cis-regulatory elements. Twelve core NtNPF genes, which were strongly induced by N limitation, were identified based on the RNA-seq data. Furthermore, N deprivation severely impaired plant growth of two cigar tobaccos, and CX26 may be more sensitive to N deficiency than CX14. Moreover, 12 hub genes respond differently to N deficiency between the two cultivars, indicating the vital roles in regulating N uptake and transport in cigar tobacco. The findings here contribute towards a better knowledge of the NtNPF genes and lay the foundation for further functional analysis of cigar tobacco.

Nitrogen Assimilation Related Genes in Brassicanapus: Systematic Characterization and Expression Analysis Identified Hub Genes in Multiple Nutrient Stress Responses

Nitrogen (N) is an essential macronutrient for plants. However, little is known about the molecular regulation of N assimilation in Brassica napus, one of the most important oil crops worldwide. Here, we carried out a comprehensive genome-wide analysis of the N assimilation related genes (NAGs) in B. napus. A total of 67 NAGs were identified encoding major enzymes involved in N assimilation, including asparagine synthetase (AS), glutamate dehydrogenase (GDH), glutamine oxoglutarate aminotransferase (GOGAT), glutamine synthetase (GS), nitrite reductase (NiR), nitrate reductase (NR). The syntenic analysis revealed that segmental duplication and whole-genome duplication were the main expansion pattern during gene evolution. Each NAG family showed different degrees of differentiation in characterization, gene structure, conserved motifs and cis-elements. Furthermore, diverse responses of NAG to multiple nutrient stresses were observed. Among them, more NAGs were regulated by N deficiency and ammonium toxicity than by phosphorus and potassium deprivations. Moreover, 12 hub genes responding to N starvation were identified, which may play vital roles in N utilization. Taken together, our results provide a basis for further functional research of NAGs in rapeseed N assimilation and also put forward new points in their responses to contrasting nutrient stresses.

Construction of a hierarchical gene regulatory network centered around a transcription factor

We have modified a multitude of transcription factors (TFs) in numerous plant species and some animal species, and obtained transgenic lines that exhibit phenotypic alterations. Whenever we observe phenotypic changes in a TF's transgenic lines, we are always eager to identify its target genes, collaborative regulators and even upstream high hierarchical regulators. This issue can be addressed by establishing a multilayered hierarchical gene regulatory network (ML-hGRN) centered around a given TF. In this article, a practical approach for constructing an ML-hGRN centered on a TF using a combined approach of top-down and bottom-up network construction methods is described. Strategies for constructing ML-hGRNs are vitally important, as these networks provide key information to advance our understanding of how biological processes are regulated.

The NIN-like protein 5 (ZmNLP5) transcription factor is involved in modulating the nitrogen response in maize

Maize exhibits marked growth and yield response to supplemental nitrogen (N). Here, we report the functional characterization of a maize NIN-like protein ZmNLP5 as a central hub in a molecular network associated with N metabolism. Predominantly expressed and accumulated in roots and vascular tissues, ZmNLP5 was shown to rapidly respond to nitrate treatment. Under limited N supply, compared with that of wild-type (WT) seedlings, the zmnlp5 mutant seedlings accumulated less nitrate and nitrite in the root tissues and ammonium in the shoot tissues. The zmnlp5 mutant plants accumulated less nitrogen than the WT plants in the ear leaves and seed kernels. Furthermore, the mutants carrying the transgenic ZmNLP5 cDNA fragment significantly increased the nitrate content in the root tissues compared with that of the zmnlp5 mutants. In the zmnlp5 mutant plants, loss of the ZmNLP5 function led to changes in expression for a significant number of genes involved in N signalling and metabolism. We further show that ZmNLP5 directly regulates the expression of nitrite reductase 1.1 (ZmNIR1.1) by binding to the nitrate-responsive cis-element at the 5' UTR of the gene. Interestingly, a natural loss-of-function allele of ZmNLP5 in Mo17 conferred less N accumulation in the ear leaves and seed kernels resembling that of the zmnlp5 mutant plants. Our findings show that ZmNLP5 is involved in mediating the plant response to N in maize.

Meta-analysis of mammary RNA seq datasets reveals the molecular understanding of bovine lactation biology

A better understanding of the biology of lactation, both in terms of gene expression and the identification of candidate genes for the production of milk and its components, is made possible by recent advances in RNA seq technology. The purpose of this study was to understand the synthesis of milk components and the molecular pathways involved, as well as to identify candidate genes for milk production traits within whole mammary transcriptomic datasets. We performed a meta-analysis of publically available RNA seq transcriptome datasets of mammary tissue/milk somatic cells. In total, 11 562 genes were commonly identified from all RNA seq based mammary gland transcriptomes. Functional annotation of commonly expressed genes revealed the molecular processes that contribute to the synthesis of fats, proteins, and lactose in mammary secretory cells and the molecular pathways responsible for milk synthesis. In addition, we identified several candidate genes responsible for milk production traits and constructed a gene regulatory network for RNA seq data. In conclusion, this study provides a basic understanding of the lactation biology of cows at the gene expression level.

Analysis of Gene Regulatory Networks of Maize in Response to Nitrogen

Nitrogen (N) fertilizer has a major influence on the yield and quality. Understanding and optimising the response of crop plants to nitrogen fertilizer usage is of central importance in enhancing food security and agricultural sustainability. In this study, the analysis of gene regulatory networks reveals multiple genes and biological processes in response to N. Two microarray studies have been used to infer components of the nitrogen-response network. Since they used different array technologies, a map linking the two probe sets to the maize B73 reference genome has been generated to allow comparison. Putative Arabidopsis homologues of maize genes were used to query the Biological General Repository for Interaction Datasets (BioGRID) network, which yielded the potential involvement of three transcription factors (TFs) (GLK5, MADS64 and bZIP108) and a Calcium-dependent protein kinase. An Artificial Neural Network was used to identify influential genes and retrieved bZIP108 and WRKY36 as significant TFs in both microarray studies, along with genes for Asparagine Synthetase, a dual-specific protein kinase and a protein phosphatase. The output from one study also suggested roles for microRNA (miRNA) 399b and Nin-like Protein 15 (NLP15). Co-expression-network analysis of TFs with closely related profiles to known Nitrate-responsive genes identified GLK5, GLK8 and NLP15 as candidate regulators of genes repressed under low Nitrogen conditions, while bZIP108 might play a role in gene activation.

Construction and Optimization of a Large Gene Coexpression Network in Maize Using RNA-Seq Data

With the emergence of massively parallel sequencing, genomewide expression data production has reached an unprecedented level. This abundance of data has greatly facilitated maize research, but may not be amenable to traditional analysis techniques that were optimized for other data types. Using publicly available data, a gene coexpression network (GCN) can be constructed and used for gene function prediction, candidate gene selection, and improving understanding of regulatory pathways. Several GCN studies have been done in maize (Zea mays), mostly using microarray datasets. To build an optimal GCN from plant materials RNA-Seq data, parameters for expression data normalization and network inference were evaluated. A comprehensive evaluation of these two parameters and a ranked aggregation strategy on network performance, using libraries from 1266 maize samples, were conducted. Three normalization methods and 10 inference methods, including six correlation and four mutual information methods, were tested. The three normalization methods had very similar performance. For network inference, correlation methods performed better than mutual information methods at some genes. Increasing sample size also had a positive effect on GCN. Aggregating single networks together resulted in improved performance compared to single networks.

A Transcriptomic Comparison of Two Bambara Groundnut Landraces under Dehydration Stress

The ability to grow crops under low-water conditions is a significant advantage in relation to global food security. Bambara groundnut is an underutilised crop grown by subsistence farmers in Africa and is known to survive in regions of water deficit. This study focuses on the analysis of the transcriptomic changes in two bambara groundnut landraces in response to dehydration stress. A cross-species hybridisation approach based on the Soybean Affymetrix GeneChip array has been employed. The differential gene expression analysis of a water-limited treatment, however, showed that the two landraces responded with almost completely different sets of genes. Hence, both landraces with very similar genotypes (as assessed by the hybridisation of genomic DNA onto the Soybean Affymetrix GeneChip) showed contrasting transcriptional behaviour in response to dehydration stress. In addition, both genotypes showed a high expression of dehydration-associated genes, even under water-sufficient conditions. Several gene regulators were identified as potentially important. Some are already known, such as WRKY40, but others may also be considered, namely PRR7, ATAUX2-11, CONSTANS-like 1, MYB60, AGL-83, and a Zinc-finger protein. These data provide a basis for drought trait research in the bambara groundnut, which will facilitate functional genomics studies. An analysis of this dataset has identified that both genotypes appear to be in a dehydration-ready state, even in the absence of dehydration stress, and may have adapted in different ways to achieve drought resistance. This will help in understanding the mechanisms underlying the ability of crops to produce viable yields under drought conditions. In addition, cross-species hybridisation to the soybean microarray has been shown to be informative for investigating the bambara groundnut transcriptome.

Bioinformatics analysis of the proteins interacting with LASP-1 and their association with HBV-related hepatocellular carcinoma

LIM and SH3 domain protein (LASP-1) is responsible for the development of several types of human cancers via the interaction with other proteins; however, the precise biological functions of proteins interacting with LASP-1 are not fully clarified. Although the role of LASP-1 in hepatocarcinogenesis has been reported, the implication of LASP-1 interactors in HBV-related hepatocellular carcinoma (HCC) is not clearly evaluated. We obtained information regarding LASP-1 interactors from public databases and published studies. Via bioinformatics analysis, we found that LASP-1 interactors were related to distinct molecular functions and associated with various biological processes. Through an integrated network analysis of the interaction and pathways of LASP-1 interactors, cross-talk between different proteins and associated pathways was found. In addition, LASP-1 and several its interactors are significantly altered in HBV-related HCC through microarray analysis and could form a complex co-expression network. In the disease, LASP-1 and its interactors were further predicted to be regulated by a complex interaction network composed of different transcription factors. Besides, numerous LASP-1 interactors were associated with various clinical factors and related to the survival and recurrence of HBV-related HCC. Taken together, these results could help enrich our understanding of LASP-1 interactors and their relationships with HBV-related HCC.

Identification of DEP domain-containing proteins by a machine learning method and experimental analysis of their expression in human HCC tissues

The Dishevelled/EGL-10/Pleckstrin (DEP) domain-containing (DEPDC) proteins have seven members. However, whether this superfamily can be distinguished from other proteins based only on the amino acid sequences, remains unknown. Here, we describe a computational method to segregate DEPDCs and non-DEPDCs. First, we examined the Pfam numbers of the known DEPDCs and used the longest sequences for each Pfam to construct a phylogenetic tree. Subsequently, we extracted 188-dimensional (188D) and 20D features of DEPDCs and non-DEPDCs and classified them with random forest classifier. We also mined the motifs of human DEPDCs to find the related domains. Finally, we designed experimental verification methods of human DEPDC expression at the mRNA level in hepatocellular carcinoma (HCC) and adjacent normal tissues. The phylogenetic analysis showed that the DEPDCs superfamily can be divided into three clusters. Moreover, the 188D and 20D features can both be used to effectively distinguish the two protein types. Motif analysis revealed that the DEP and RhoGAP domain was common in human DEPDCs, human HCC and the adjacent tissues that widely expressed DEPDCs. However, their regulation was not identical. In conclusion, we successfully constructed a binary classifier for DEPDCs and experimentally verified their expression in human HCC tissues.

Computing and Applying Atomic Regulons to Understand Gene Expression and Regulation

Understanding gene function and regulation is essential for the interpretation, prediction, and ultimate design of cell responses to changes in the environment. An important step toward meeting the challenge of understanding gene function and regulation is the identification of sets of genes that are always co-expressed. These gene sets, Atomic Regulons (ARs), represent fundamental units of function within a cell and could be used to associate genes of unknown function with cellular processes and to enable rational genetic engineering of cellular systems. Here, we describe an approach for inferring ARs that leverages large-scale expression data sets, gene context, and functional relationships among genes. We computed ARs for Escherichia coli based on 907 gene expression experiments and compared our results with gene clusters produced by two prevalent data-driven methods: Hierarchical clustering and k-means clustering. We compared ARs and purely data-driven gene clusters to the curated set of regulatory interactions for E. coli found in RegulonDB, showing that ARs are more consistent with gold standard regulons than are data-driven gene clusters. We further examined the consistency of ARs and data-driven gene clusters in the context of gene interactions predicted by Context Likelihood of Relatedness (CLR) analysis, finding that the ARs show better agreement with CLR predicted interactions. We determined the impact of increasing amounts of expression data on AR construction and find that while more data improve ARs, it is not necessary to use the full set of gene expression experiments available for E. coli to produce high quality ARs. In order to explore the conservation of co-regulated gene sets across different organisms, we computed ARs for Shewanella oneidensis, Pseudomonas aeruginosa, Thermus thermophilus, and Staphylococcus aureus, each of which represents increasing degrees of phylogenetic distance from E. coli. Comparison of the organism-specific ARs showed that the consistency of AR gene membership correlates with phylogenetic distance, but there is clear variability in the regulatory networks of closely related organisms. As large scale expression data sets become increasingly common for model and non-model organisms, comparative analyses of atomic regulons will provide valuable insights into fundamental regulatory modules used across the bacterial domain.

Mining Functional Modules in Heterogeneous Biological Networks Using Multiplex PageRank Approach

Identification of functional modules/sub-networks in large-scale biological networks is one of the important research challenges in current bioinformatics and systems biology. Approaches have been developed to identify functional modules in single-class biological networks; however, methods for systematically and interactively mining multiple classes of heterogeneous biological networks are lacking. In this paper, we present a novel algorithm (called mPageRank) that utilizes the Multiplex PageRank approach to mine functional modules from two classes of biological networks. We demonstrate the capabilities of our approach by successfully mining functional biological modules through integrating expression-based gene-gene association networks and protein-protein interaction networks. We first compared the performance of our method with that of other methods using simulated data. We then applied our method to identify the cell division cycle related functional module and plant signaling defense-related functional module in the model plant Arabidopsis thaliana. Our results demonstrated that the mPageRank method is effective for mining sub-networks in both expression-based gene-gene association networks and protein-protein interaction networks, and has the potential to be adapted for the discovery of functional modules/sub-networks in other heterogeneous biological networks. The mPageRank executable program, source code, the datasets and results of the presented two case studies are publicly and freely available at http://plantgrn.noble.org/MPageRank/.

DTW-MIC Coexpression Networks from Time-Course Data

When modeling coexpression networks from high-throughput time course data, Pearson Correlation Coefficient (PCC) is one of the most effective and popular similarity functions. However, its reliability is limited since it cannot capture non-linear interactions and time shifts. Here we propose to overcome these two issues by employing a novel similarity function, Dynamic Time Warping Maximal Information Coefficient (DTW-MIC), combining a measure taking care of functional interactions of signals (MIC) and a measure identifying time lag (DTW). By using the Hamming-Ipsen-Mikhailov (HIM) metric to quantify network differences, the effectiveness of the DTW-MIC approach is demonstrated on a set of four synthetic and one transcriptomic datasets, also in comparison to TimeDelay ARACNE and Transfer Entropy.

Methanol regulated yeast promoters: production vehicles and toolbox for synthetic biology

Promoters are indispensable elements of a standardized parts collection for synthetic biology. Regulated promoters of a wide variety of well-defined induction ratios and expression strengths are highly interesting for many applications. Exemplarily, we discuss the application of published genome scale transcriptomics data for the primary selection of methanol inducible promoters of the yeast Pichia pastoris (Komagataella sp.). Such a promoter collection can serve as an excellent toolbox for cell and metabolic engineering, and for gene expression to produce heterologous proteins.

LegumeIP 2.0--a platform for the study of gene function and genome evolution in legumes

The LegumeIP 2.0 database hosts large-scale genomics and transcriptomics data and provides integrative bioinformatics tools for the study of gene function and evolution in legumes. Our recent updates in LegumeIP 2.0 include gene and protein sequences, gene models and annotations, syntenic regions, protein families and phylogenetic trees for six legume species: Medicago truncatula, Glycine max (soybean), Lotus japonicus, Phaseolus vulgaris (common bean), Cicer arietinum (chickpea) and Cajanus cajan (pigeon pea) and two outgroup reference species: Arabidopsis thaliana and Poplar trichocarpa. Moreover, the LegumeIP 2.0 features the following new data resources and bioinformatics tools: (i) an integrative gene expression atlas for four model legumes that include 550 array hybridizations from M. truncatula, 962 gene expression profiles of G. max, 276 array hybridizations from L. japonicas and 56 RNA-Seq-based gene expression profiles for C. arietinum. These datasets were manually curated and hierarchically organized based on Experimental Ontology and Plant Ontology so that users can browse, search, and retrieve data for their selected experiments. (ii) New functions/analytical tools to query, mine and visualize large-scale gene sequences, annotations and transcriptome profiles. Users may select a subset of expression experiments and visualize and compare expression profiles for multiple genes. The LegumeIP 2.0 database is freely available to the public at http://plantgrn.noble.org/LegumeIP/.

Functional characterization of drought-responsive modules and genes in Oryza sativa: a network-based approach

Drought is one of the major environmental stress conditions affecting the yield of rice across the globe. Unraveling the functional roles of the drought-responsive genes and their underlying molecular mechanisms will provide important leads to improve the yield of rice. Co-expression relationships derived from condition-dependent gene expression data is an effective way to identify the functional associations between genes that are part of the same biological process and may be under similar transcriptional control. For this purpose, vast amount of freely available transcriptomic data may be used. In this study, we consider gene expression data for different tissues and developmental stages in response to drought stress. We analyze the network of co-expressed genes to identify drought-responsive genes modules in a tissue and stage-specific manner based on differential expression and gene enrichment analysis. Taking cues from the systems-level behavior of these modules, we propose two approaches to identify clusters of tightly co-expressed/co-regulated genes. Using graph-centrality measures and differential gene expression, we identify biologically informative genes that lack any functional annotation. We show that using orthologous information from other plant species, the conserved co-expression patterns of the uncharacterized genes can be identified. Presence of a conserved neighborhood enables us to extrapolate functional annotation. Alternatively, we show that single 'guide-gene' approach can help in understanding tissue-specific transcriptional regulation of uncharacterized genes. Finally, we confirm the predicted roles of uncharacterized genes by the analysis of conserved cis-elements and explain the possible roles of these genes toward drought tolerance.

The personal touch: strategies toward personalized vaccines and predicting immune responses to them

The impact of vaccines on public health and wellbeing has been profound. Smallpox has been eradicated, polio is nearing eradication, and multiple diseases have been eliminated from certain areas of the world. Unfortunately, we now face diseases such as hepatitis C, malaria or tuberculosis, as well as new and re-emerging pathogens for which we lack effective vaccines. Empirical approaches to vaccine development have been successful in the past, but may not be up to the current infectious disease challenges facing us. New, directed approaches to vaccine design, development, and testing need to be developed. Ideally these approaches will capitalize on cutting-edge technologies, advanced analytical and modeling strategies, and up-to-date knowledge of both pathogen and host. These approaches will pay particular attention to the causes of inter-individual variation in vaccine response in order to develop new vaccines tailored to the unique needs of individuals and communities within the population.