A systems-genetics approach and data mining tool to assist in the discovery of genes underlying complex traits in Oryza sativa

Addressing noise in co-expression network construction

Gene co-expression networks (GCNs) provide multiple benefits to molecular research including hypothesis generation and biomarker discovery. Transcriptome profiles serve as input for GCN construction and are derived from increasingly larger studies with samples across multiple experimental conditions, treatments, time points, genotypes, etc. Such experiments with larger numbers of variables confound discovery of true network edges, exclude edges and inhibit discovery of context (or condition) specific network edges. To demonstrate this problem, a 475-sample dataset is used to show that up to 97% of GCN edges can be misleading because correlations are false or incorrect. False and incorrect correlations can occur when tests are applied without ensuring assumptions are met, and pairwise gene expression may not meet test assumptions if the expression of at least one gene in the pairwise comparison is a function of multiple confounding variables. The ‘one-size-fits-all’ approach to GCN construction is therefore problematic for large, multivariable datasets. Recently, the Knowledge Independent Network Construction toolkit has been used in multiple studies to provide a dynamic approach to GCN construction that ensures statistical tests meet assumptions and confounding variables are addressed. Additionally, it can associate experimental context for each edge of the network resulting in context-specific GCNs (csGCNs). To help researchers recognize such challenges in GCN construction, and the creation of csGCNs, we provide a review of the workflow.

Tripal, a community update after 10 years of supporting open source, standards-based genetic, genomic and breeding databases

Online, open access databases for biological knowledge serve as central repositories for research communities to store, find and analyze integrated, multi-disciplinary datasets. With increasing volumes, complexity and the need to integrate genomic, transcriptomic, metabolomic, proteomic, phenomic and environmental data, community databases face tremendous challenges in ongoing maintenance, expansion and upgrades. A common infrastructure framework using community standards shared by many databases can reduce development burden, provide interoperability, ensure use of common standards and support long-term sustainability. Tripal is a mature, open source platform built to meet this need. With ongoing improvement since its first release in 2009, Tripal provides full functionality for searching, browsing, loading and curating numerous types of data and is a primary technology powering at least 31 publicly available databases spanning plants, animals and human data, primarily storing genomics, genetics and breeding data. Tripal software development is managed by a shared, inclusive governance structure including both project management and advisory teams. Here, we report on the most important and innovative aspects of Tripal after 11 years development, including integration of diverse types of biological data, successful collaborative projects across member databases, and support for implementing FAIR principles.

Agro-Physiologic Responses and Stress-Related Gene Expression of Four Doubled Haploid Wheat Lines under Salinity Stress Conditions

Simple Summary

Productivity of wheat can be enhanced using salt-tolerant genotypes. However, the assessment of salt tolerance potential in wheat through agro-physiological traits and stress-related gene expression analysis could potentially minimize the cost of breeding programs and be a powerful way for the selection of the most salt-tolerant genotype. The study evaluated the salt tolerance potential of four doubled haploid lines of wheat and compared them with the check cultivar Sakha-93 using an extensive set of agro-physiologic parameters and salt-stress-related gene expressions. The results indicated that the five genotypes tested displayed reduction in all traits evaluated except the canopy temperature and electrical conductivity, which had the greatest decline occurring in the check cultivar and the least decline in DHL2. The genotypes DHL21 and DHL5 exhibited increased expression rate of salt-stress-related genes under salt stress conditions. The multiple linear regression model and path coefficient analysis showed a coefficient of determination of 0.93. Concluding, the number of spikelets, and/or number of kernels were identified to be unbiased traits for assessing wheat DHLs under salinity conditions, given their contribution and direct impact on the grain yield. Moreover, the two most salt-tolerant genotypes DHL2 and DHL21 can be useful as genetic resources for future breeding programs.

Abstract

Salinity majorly hinders horizontal and vertical expansion in worldwide wheat production. Productivity can be enhanced using salt-tolerant wheat genotypes. However, the assessment of salt tolerance potential in bread wheat doubled haploid lines (DHL) through agro-physiological traits and stress-related gene expression analysis could potentially minimize the cost of breeding programs and be a powerful way for the selection of the most salt-tolerant genotype. We used an extensive set of agro-physiologic parameters and salt-stress-related gene expressions. Multivariate analysis was used to detect phenotypic and genetic variations of wheat genotypes more closely under salinity stress, and we analyzed how these strategies effectively balance each other. Four doubled haploid lines (DHLs) and the check cultivar (Sakha93) were evaluated in two salinity levels (without and 150 mM NaCl) until harvest. The five genotypes showed reduced growth under 150 mM NaCl; however, the check cultivar (Sakha93) died at the beginning of the flowering stage. Salt stress induced reduction traits, except the canopy temperature and initial electrical conductivity, which was found in each of the five genotypes, with the greatest decline occurring in the check cultivar (Sakha-93) and the least decline in DHL2. The genotypes DHL21 and DHL5 exhibited increased expression rate of salt-stress-related genes (TaNHX1, TaHKT1, and TaCAT1) compared with DHL2 and Sakha93 under salt stress conditions. Principle component analysis detection of the first two components explains 70.78% of the overall variation of all traits (28 out of 32 traits). A multiple linear regression model and path coefficient analysis showed a coefficient of determination (R²) of 0.93. The models identified two interpretive variables, number of spikelets, and/or number of kernels, which can be unbiased traits for assessing wheat DHLs under salinity stress conditions, given their contribution and direct impact on the grain yield.

Exploration into biomarker potential of region-specific brain gene co-expression networks

The human brain is a complex organ that consists of several regions each with a unique gene expression pattern. Our intent in this study was to construct a gene co-expression network (GCN) for the normal brain using RNA expression profiles from the Genotype-Tissue Expression (GTEx) project. The brain GCN contains gene correlation relationships that are broadly present in the brain or specific to thirteen brain regions, which we later combined into six overarching brain mini-GCNs based on the brain's structure. Using the expression profiles of brain region-specific GCN edges, we determined how well the brain region samples could be discriminated from each other, visually with t-SNE plots or quantitatively with the Gene Oracle deep learning classifier. Next, we tested these gene sets on their relevance to human tumors of brain and non-brain origin. Interestingly, we found that genes in the six brain mini-GCNs showed markedly higher mutation rates in tumors relative to matched sets of random genes. Further, we found that cortex genes subdivided Head and Neck Squamous Cell Carcinoma (HNSC) tumors and Pheochromocytoma and Paraganglioma (PCPG) tumors into distinct groups. The brain GCN and mini-GCNs are useful resources for the classification of brain regions and identification of biomarker genes for brain related phenotypes.

The Plant-Transpiration Response to Vapor Pressure Deficit (VPD) in Durum Wheat Is Associated With Differential Yield Performance and Specific Expression of Genes Involved in Primary Metabolism and Water Transport

The regulation of plant transpiration was proposed as a key factor affecting transpiration efficiency and agronomical adaptation of wheat to water-limited Mediterranean environments. However, to date no studies have related this trait to crop performance in the field. In this study, the transpiration response to increasing vapor pressure deficit (VPD) of modern Spanish semi-dwarf durum wheat lines was evaluated under controlled conditions at vegetative stage, and the agronomical performance of the same set of lines was assessed at grain filling as well as grain yield at maturity, in Mediterranean environments ranging from water stressed to good agronomical conditions. A group of linear-transpiration response (LTR) lines exhibited better performance in grain yield and biomass compared to segmented-transpiration response (STR) lines, particularly in the wetter environments, whereas the reverse occurred only in the most stressed trial. LTR lines generally exhibited better water status (stomatal conductance) and larger green biomass (vegetation indices) during the reproductive stage than STR lines. In both groups, the responses to growing conditions were associated with the expression levels of dehydration-responsive transcription factors (DREB) leading to different performances of primary metabolism-related enzymes. Thus, the response of LTR lines under fair to good conditions was associated with higher transcription levels of genes involved in nitrogen (GS1 and GOGAT) and carbon (RCBL) metabolism, as well as water transport (TIP1.1). In conclusion, modern durum wheat lines differed in their response to water loss, the linear transpiration seemed to favor uptake and transport of water and nutrients, and photosynthetic metabolism led to higher grain yield except for very harsh drought conditions. The transpiration response to VPD may be a trait to further explore when selecting adaptation to specific water conditions.

Rice Genome-Scale Network Integration Reveals Transcriptional Regulators of Grass Cell Wall Synthesis

Grasses have evolved distinct cell wall composition and patterning relative to dicotyledonous plants. However, despite the importance of this plant family, transcriptional regulation of its cell wall biosynthesis is poorly understood. To identify grass cell wall-associated transcription factors, we constructed the Rice Combined mutual Ranked Network (RCRN). The RCRN covers >90% of annotated rice (Oryza sativa) genes, is high quality, and includes most grass-specific cell wall genes, such as mixed-linkage glucan synthases and hydroxycinnamoyl acyltransferases. Comparing the RCRN and an equivalent Arabidopsis network suggests that grass orthologs of most genetically verified eudicot cell wall regulators also control this process in grasses, but some transcription factors vary significantly in network connectivity between these divergent species. Reverse genetics, yeast-one-hybrid, and protoplast-based assays reveal that OsMYB61a activates a grass-specific acyltransferase promoter, which confirms network predictions and supports grass-specific cell wall synthesis genes being incorporated into conserved regulatory circuits. In addition, 10 of 15 tested transcription factors, including six novel Wall-Associated regulators (WAP1, WACH1, WAHL1, WADH1, OsMYB13a, and OsMYB13b), alter abundance of cell wall-related transcripts when transiently expressed. The results highlight the quality of the RCRN for examining rice biology, provide insight into the evolution of cell wall regulation, and identify network nodes and edges that are possible leads for improving cell wall composition.

MPIGeneNet: Parallel Calculation of Gene Co-Expression Networks on Multicore Clusters

In this work, we present MPIGeneNet, a parallel tool that applies Pearson's correlation and Random Matrix Theory to construct gene co-expression networks. It is based on the state-of-the-art sequential tool RMTGeneNet, which provides networks with high robustness and sensitivity at the expenses of relatively long runtimes for large scale input datasets. MPIGeneNet returns the same results as RMTGeneNet but improves the memory management, reduces the I/O cost, and accelerates the two most computationally demanding steps of co-expression network construction by exploiting the compute capabilities of common multicore CPU clusters. Our performance evaluation on two different systems using three typical input datasets shows that MPIGeneNet is significantly faster than RMTGeneNet. As an example, our tool is up to 175.41 times faster on a cluster with eight nodes, each one containing two 12-core Intel Haswell processors. The source code of MPIGeneNet, as well as a reference manual, are available at https://sourceforge.net/projects/mpigenenet/.

Sorting Five Human Tumor Types Reveals Specific Biomarkers and Background Classification Genes

We applied two state-of-the-art, knowledge independent data-mining methods - Dynamic Quantum Clustering (DQC) and t-Distributed Stochastic Neighbor Embedding (t-SNE) - to data from The Cancer Genome Atlas (TCGA). We showed that the RNA expression patterns for a mixture of 2,016 samples from five tumor types can sort the tumors into groups enriched for relevant annotations including tumor type, gender, tumor stage, and ethnicity. DQC feature selection analysis discovered 48 core biomarker transcripts that clustered tumors by tumor type. When these transcripts were removed, the geometry of tumor relationships changed, but it was still possible to classify the tumors using the RNA expression profiles of the remaining transcripts. We continued to remove the top biomarkers for several iterations and performed cluster analysis. Even though the most informative transcripts were removed from the cluster analysis, the sorting ability of remaining transcripts remained strong after each iteration. Further, in some iterations we detected a repeating pattern of biological function that wasn't detectable with the core biomarker transcripts present. This suggests the existence of a "background classification" potential in which the pattern of gene expression after continued removal of "biomarker" transcripts could still classify tumors in agreement with the tumor type.

Extension modules for storage, visualization and querying of genomic, genetic and breeding data in Tripal databases

Tripal is an open-source database platform primarily used for development of genomic, genetic and breeding databases. We report here on the release of the Chado Loader, Chado Data Display and Chado Search modules to extend the functionality of the core Tripal modules. These new extension modules provide additional tools for (1) data loading, (2) customized visualization and (3) advanced search functions for supported data types such as organism, marker, QTL/Mendelian Trait Loci, germplasm, map, project, phenotype, genotype and their respective metadata. The Chado Loader module provides data collection templates in Excel with defined metadata and data loaders with front end forms. The Chado Data Display module contains tools to visualize each data type and the metadata which can be used as is or customized as desired. The Chado Search module provides search and download functionality for the supported data types. Also included are the tools to visualize map and species summary. The use of materialized views in the Chado Search module enables better performance as well as flexibility of data modeling in Chado, allowing existing Tripal databases with different metadata types to utilize the module. These Tripal Extension modules are implemented in the Genome Database for Rosaceae (rosaceae.org), CottonGen (cottongen.org), Citrus Genome Database (citrusgenomedb.org), Genome Database for Vaccinium (vaccinium.org) and the Cool Season Food Legume Database (coolseasonfoodlegume.org). Database URL: https://www.citrusgenomedb.org/, https://www.coolseasonfoodlegume.org/, https://www.cottongen.org/, https://www.rosaceae.org/, https://www.vaccinium.org/.

Discovering Condition-Specific Gene Co-Expression Patterns Using Gaussian Mixture Models: A Cancer Case Study

A gene co-expression network (GCN) describes associations between genes and points to genetic coordination of biochemical pathways. However, genetic correlations in a GCN are only detectable if they are present in the sampled conditions. With the increasing quantity of gene expression samples available in public repositories, there is greater potential for discovery of genetic correlations from a variety of biologically interesting conditions. However, even if gene correlations are present, their discovery can be masked by noise. Noise is introduced from natural variation (intrinsic and extrinsic), systematic variation (caused by sample measurement protocols and instruments), and algorithmic and statistical variation created by selection of data processing tools. A variety of published studies, approaches and methods attempt to address each of these contributions of variation to reduce noise. Here we describe an approach using Gaussian Mixture Models (GMMs) to address natural extrinsic (condition-specific) variation during network construction from mixed input conditions. To demonstrate utility, we build and analyze a condition-annotated GCN from a compendium of 2,016 mixed gene expression data sets from five tumor subtypes obtained from The Cancer Genome Atlas. Our results show that GMMs help discover tumor subtype specific gene co-expression patterns (modules) that are significantly enriched for clinical attributes.

ATGC transcriptomics: a web-based application to integrate, explore and analyze de novo transcriptomic data

BACKGROUND

In the last years, applications based on massively parallelized RNA sequencing (RNA-seq) have become valuable approaches for studying non-model species, e.g., without a fully sequenced genome. RNA-seq is a useful tool for detecting novel transcripts and genetic variations and for evaluating differential gene expression by digital measurements. The large and complex datasets resulting from functional genomic experiments represent a challenge in data processing, management, and analysis. This problem is especially significant for small research groups working with non-model species.

RESULTS

We developed a web-based application, called ATGC transcriptomics, with a flexible and adaptable interface that allows users to work with new generation sequencing (NGS) transcriptomic analysis results using an ontology-driven database. This new application simplifies data exploration, visualization, and integration for a better comprehension of the results.

CONCLUSIONS

ATGC transcriptomics provides access to non-expert computer users and small research groups to a scalable storage option and simple data integration, including database administration and management. The software is freely available under the terms of GNU public license at http://atgcinta.sourceforge.net .

blend4php: a PHP API for galaxy

Galaxy is a popular framework for execution of complex analytical pipelines typically for large data sets, and is a commonly used for (but not limited to) genomic, genetic and related biological analysis. It provides a web front-end and integrates with high performance computing resources. Here we report the development of the blend4php library that wraps Galaxy’s RESTful API into a PHP-based library. PHP-based web applications can use blend4php to automate execution, monitoring and management of a remote Galaxy server, including its users, workflows, jobs and more. The blend4php library was specifically developed for the integration of Galaxy with Tripal, the open-source toolkit for the creation of online genomic and genetic web sites. However, it was designed as an independent library for use by any application, and is freely available under version 3 of the GNU Lesser General Public License (LPGL v3.0) at https://github.com/galaxyproject/blend4php.

Database URL:https://github.com/galaxyproject/blend4php

Designing climate-resilient rice with ideal grain quality suited for high-temperature stress

To ensure rice food security, the target outputs of future rice breeding programmes should focus on developing climate-resilient rice varieties with emphasis on increased head rice yield coupled with superior grain quality. This challenge is made greater by a world that is increasingly becoming warmer. Such environmental changes dramatically impact head rice and milling yield as well as increasing chalkiness because of impairment in starch accumulation and other storage biosynthetic pathways in the grain. This review highlights the knowledge gained through gene discovery via quantitative trait locus (QTL) cloning and structural-functional genomic strategies to reduce chalk, increase head rice yield, and develop stable lines with optimum grain quality in challenging environments. The newly discovered genes and the knowledge gained on the influence of specific alleles related to stability of grain quality attributes provide a robust platform for marker-assisted selection in breeding to design heat-tolerant rice varieties with superior grain quality. Using the chalkiness trait in rice as a case study, we demonstrate here that the emerging field of systems genetics can help fast-track the identification of novel alleles and gene targets that can be pyramided for the development of environmentally robust rice varieties that possess improved grain quality.

Systems genetics: a paradigm to improve discovery of candidate genes and mechanisms underlying complex traits

Understanding the control of any trait optimally requires the detection of causal genes, gene interaction, and mechanism of action to discover and model the biochemical pathways underlying the expressed phenotype. Functional genomics techniques, including RNA expression profiling via microarray and high-throughput DNA sequencing, allow for the precise genome localization of biological information. Powerful genetic approaches, including quantitative trait locus (QTL) and genome-wide association study mapping, link phenotype with genome positions, yet genetics is less precise in localizing the relevant mechanistic information encoded in DNA. The coupling of salient functional genomic signals with genetically mapped positions is an appealing approach to discover meaningful gene-phenotype relationships. Techniques used to define this genetic-genomic convergence comprise the field of systems genetics. This short review will address an application of systems genetics where RNA profiles are associated with genetically mapped genome positions of individual genes (eQTL mapping) or as gene sets (co-expression network modules). Both approaches can be applied for knowledge independent selection of candidate genes (and possible control mechanisms) underlying complex traits where multiple, likely unlinked, genomic regions might control specific complex traits.

Tripal v1.1: a standards-based toolkit for construction of online genetic and genomic databases

Tripal is an open-source freely available toolkit for construction of online genomic and genetic databases. It aims to facilitate development of community-driven biological websites by integrating the GMOD Chado database schema with Drupal, a popular website creation and content management software. Tripal provides a suite of tools for interaction with a Chado database and display of content therein. The tools are designed to be generic to support the various ways in which data may be stored in Chado. Previous releases of Tripal have supported organisms, genomic libraries, biological stocks, stock collections and genomic features, their alignments and annotations. Also, Tripal and its extension modules provided loaders for commonly used file formats such as FASTA, GFF, OBO, GAF, BLAST XML, KEGG heir files and InterProScan XML. Default generic templates were provided for common views of biological data, which could be customized using an open Application Programming Interface to change the way data are displayed. Here, we report additional tools and functionality that are part of release v1.1 of Tripal. These include (i) a new bulk loader that allows a site curator to import data stored in a custom tab delimited format; (ii) full support of every Chado table for Drupal Views (a powerful tool allowing site developers to construct novel displays and search pages); (iii) new modules including ‘Feature Map’, ‘Genetic’, ‘Publication’, ‘Project’, ‘Contact’ and the ‘Natural Diversity’ modules. Tutorials, mailing lists, download and set-up instructions, extension modules and other documentation can be found at the Tripal website located at http://tripal.info.

Database URL: http://tripal.info/