Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Systems genomics of salinity stress response in rice

Populations can adapt to stressful environments through changes in gene expression. However, the role of gene regulation in mediating stress response and adaptation remains largely unexplored. Here, we use an integrative field dataset obtained from 780 plants of Oryza sativa ssp. indica (rice) grown in a field experiment under normal or moderate salt stress conditions to examine selection and evolution of gene expression variation under salinity stress conditions. We find that salinity stress induces increased selective pressure on gene expression. Further, we show that trans-eQTLs rather than cis-eQTLs are primarily associated with rice's gene expression under salinity stress, potentially via a few master-regulators. Importantly, and contrary to the expectations, we find that cis-trans reinforcement is more common than cis-trans compensation which may be reflective of rice diversification subsequent to domestication. We further identify genetic fixation as the likely mechanism underlying this compensation/reinforcement. Additionally, we show that cis- and trans-eQTLs are under different selection regimes, giving us insights into the evolutionary dynamics of gene expression variation. By examining genomic, transcriptomic, and phenotypic variation across a rice population, we gain insights into the molecular and genetic landscape underlying adaptive salinity stress responses, which is relevant for other crops and other stresses.

Molecular mechanisms underpinning quantitative resistance to Phytophthora sojae in Glycine max using a systems genomics approach

Expression of quantitative disease resistance in many host-pathogen systems is controlled by genes at multiple loci, each contributing a small effect to the overall response. We used a systems genomics approach to study the molecular underpinnings of quantitative disease resistance in the soybean-Phytophthora sojae pathosystem, incorporating expression quantitative trait loci (eQTL) mapping and gene co-expression network analysis to identify the genes putatively regulating transcriptional changes in response to inoculation. These findings were compared to previously mapped phenotypic (phQTL) to identify the molecular mechanisms contributing to the expression of this resistance. A subset of 93 recombinant inbred lines (RILs) from a Conrad × Sloan population were inoculated with P. sojae isolate 1.S.1.1 using the tray-test method; RNA was extracted, sequenced, and the normalized read counts were genetically mapped from tissue collected at the inoculation site 24 h after inoculation from both mock and inoculated samples. In total, more than 100,000 eQTLs were mapped. There was a switch from predominantly cis-eQTLs in the mock treatment to an almost entirely nonoverlapping set of predominantly trans-eQTLs in the inoculated treatment, where greater than 100-fold more eQTLs were mapped relative to mock, indicating vast transcriptional reprogramming due to P. sojae infection occurred. The eQTLs were organized into 36 hotspots, with the four largest hotspots from the inoculated treatment corresponding to more than 70% of the eQTLs, each enriched for genes within plant-pathogen interaction pathways. Genetic regulation of trans-eQTLs in response to the pathogen was predicted to occur through transcription factors and signaling molecules involved in plant-pathogen interactions, plant hormone signal transduction, and MAPK pathways. Network analysis identified three co-expression modules that were correlated with susceptibility to P. sojae and associated with three eQTL hotspots. Among the eQTLs co-localized with phQTLs, two cis-eQTLs with putative functions in the regulation of root architecture or jasmonic acid, as well as the putative master regulators of an eQTL hotspot nearby a phQTL, represent candidates potentially underpinning the molecular control of these phQTLs for resistance.

Integrating Phenotypic and Gene Expression Linkage Mapping to Dissect Rust Resistance in Chickling Pea

Rusts are among the most important foliar biotrophic fungal diseases in legumes. Lathyrus cicera crop can be severely damaged by Uromyces pisi, to which partial resistance has been identified. Nevertheless, the underlying genetic basis and molecular mechanisms of this resistance are poorly understood in L. cicera. To prioritise the causative variants controlling partial resistance to rust in L. cicera, a recombinant inbred line (RIL) population, segregating for response to this pathogen, was used to combine the detection of related phenotypic- and expression-quantitative trait loci (pQTLs and eQTLs, respectively). RILs' U. pisi disease severity (DS) was recorded in three independent screenings at seedling (growth chamber) and in one season of exploratory screening at adult plant stage (semi-controlled field conditions). A continuous DS range was observed in both conditions and used for pQTL mapping. Different pQTLs were identified under the growth chamber and semi-controlled field conditions, indicating a distinct genetic basis depending on the plant developmental stage and/or the environment. Additionally, the expression of nine genes related to U. pisi resistance in L. cicera was quantified for each RIL individual and used for eQTL mapping. One cis-eQTL and one trans-eQTL were identified controlling the expression variation of one gene related to rust resistance - a member of glycosyl hydrolase family 17. Integrating phenotyping, gene expression and linkage mapping allowed prioritising four candidate genes relevant for disease-resistance precision breeding involved in adaptation to biotic stress, cellular, and organelle homeostasis, and proteins directly involved in plant defence.

Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus

Photosynthesis and wood formation underlie the ability of trees to provide renewable resources and perform ecological functions; however, the genetic basis and regulatory pathways coordinating these two linked processes remain unclear. Here, we used a systems genetics strategy, integrating genome-wide association studies, transcriptomic analyses, and transgenic experiments, to investigate the genetic architecture of photosynthesis and wood properties among 435 unrelated individuals of Populus tomentosa, and unravel the coordinated regulatory networks resulting in two trait categories. We detected 222 significant single-nucleotide polymorphisms, annotated to 177 candidate genes, for 10 traits of photosynthesis and wood properties. Epistasis uncovered 74 epistatic interactions for phenotypes. Strikingly, we deciphered the coordinated regulation patterns of pleiotropic genes underlying phenotypic variations for two trait categories. Furthermore, expression quantitative trait nucleotide mapping and coexpression analysis were integrated to unravel the potential transcriptional regulatory networks of candidate genes coordinating photosynthesis and wood properties. Finally, heterologous expression of two pleiotropic genes, PtoMYB62 and PtoMYB80, in Arabidopsis thaliana demonstrated that they control regulatory networks balancing photosynthesis and stem secondary cell wall components, respectively. Our study provides insights into the regulatory mechanisms coordinating photosynthesis and wood formation in poplar, and should facilitate genetic breeding in trees via molecular design.

Chromosomal characteristics of salt stress heritable gene expression in the rice genome

BACKGROUND

Gene expression is potentially an important heritable quantitative trait that mediates between genetic variation and higher-level complex phenotypes through time and condition-dependent regulatory interactions. Therefore, we sought to explore both the genomic and condition-specific characteristics of gene expression heritability within the context of chromosomal structure.

RESULTS

Heritability was estimated for biological gene expression using a diverse, 84-line, Oryza sativa (rice) population under optimal and salt-stressed conditions. Overall, 5936 genes were found to have heritable expression regardless of condition and 1377 genes were found to have heritable expression only during salt stress. These genes with salt-specific heritable expression are enriched for functional terms associated with response to stimulus and transcription factor activity. Additionally, we discovered that highly and lowly expressed genes, and genes with heritable expression are distributed differently along the chromosomes in patterns that follow previously identified high-throughput chromosomal conformation capture (Hi-C) A/B chromatin compartments. Furthermore, multiple genomic hot-spots enriched for genes with salt-specific heritability were identified on chromosomes 1, 4, 6, and 8. These hotspots were found to contain genes functionally enriched for transcriptional regulation and overlaps with a previously identified major QTL for salt-tolerance in rice.

CONCLUSIONS

Investigating the heritability of traits, and in-particular gene expression traits, is important towards developing a basic understanding of how regulatory networks behave across a population. This work provides insights into spatial patterns of heritable gene expression at the chromosomal level.

Detecting selection with a genetic cross

Distinguishing which traits have evolved under natural selection, as opposed to neutral evolution, is a major goal of evolutionary biology. Several tests have been proposed to accomplish this, but these either rely on false assumptions or suffer from low power. Here, I introduce an approach to detecting selection that makes minimal assumptions and only requires phenotypic data from ∼10 individuals. The test compares the phenotypic difference between two populations to what would be expected by chance under neutral evolution, which can be estimated from the phenotypic distribution of an F2 cross between those populations. Simulations show that the test is robust to variation in the number of loci affecting the trait, the distribution of locus effect sizes, heritability, dominance, and epistasis. Comparing its performance to the QTL sign test-an existing test of selection that requires both genotype and phenotype data-the new test achieves comparable power with 50- to 100-fold fewer individuals (and no genotype data). Applying the test to empirical data spanning over a century shows strong directional selection in many crops, as well as on naturally selected traits such as head shape in Hawaiian Drosophila and skin color in humans. Applied to gene expression data, the test reveals that the strength of stabilizing selection acting on mRNA levels in a species is strongly associated with that species' effective population size. In sum, this test is applicable to phenotypic data from almost any genetic cross, allowing selection to be detected more easily and powerfully than previously possible.

Early Responses of Brassica oleracea Roots to Zinc Supply Under Sufficient and Sub-Optimal Phosphorus Supply

Shoot zinc (Zn) concentration in Brassica oleracea is affected by soil Zn and phosphorus (P) supply. Most problematic is the negative impact of P fertilizers on Zn concentrations in crops, which makes balancing yield and mineral quality challenging. To evaluate early molecular mechanisms involved in the accumulation of large shoot Zn concentrations regardless of the P supply, two B. oleracea accessions differing in root architecture and root exudates were grown hydroponically for two weeks with different combinations of P and Zn supply. Ionome profiling and deep RNA sequencing of roots revealed interactions of P and Zn in planta, without apparent phenotypic effects. In addition, increasing P supply did not reduce tissue Zn concentration. Substantial changes in gene expression in response to different P and/or Zn supplies in roots of both accessions ensured nutritionally sufficient P and Zn uptake. Numerous genes were differentially expressed after changing Zn or P supply and most of them were unique to only one accession, highlighting their different strategies in achieving nutrient sufficiency. Thus, different gene networks responded to the changing P and Zn supply in the two accessions. Additionally, enrichment analysis of gene ontology classes revealed that genes involved in lipid metabolism, response to starvation, and anion transport mechanisms were most responsive to differences in P and Zn supply in both accessions. The results agreed with previously studies demonstrating alterations in P and Zn transport and phospholipid metabolism in response to reduced P and Zn supply. It is anticipated that improved knowledge of genes responsive to P or Zn supply will help illuminate the roles in uptake and accumulation of P and Zn and might identify candidate genes for breeding high-yield-high-Zn brassicas.

QTL and candidate genes associated with leaf anion concentrations in response to phosphate supply in Arabidopsis thaliana

BACKGROUND

Phosphorus is often present naturally in the soil as inorganic phosphate, Pi, which bio-availability is limited in many ecosystems due to low soil solubility and mobility. Plants respond to low Pi with a Pi Starvation Response, involving Pi sensing and long-distance signalling. There is extensive cross-talk between Pi homeostasis mechanisms and the homeostasis mechanism for other anions in response to Pi availability.

RESULTS

Recombinant Inbred Line (RIL) and Genome Wide Association (GWA) mapping populations, derived from or composed of natural accessions of Arabidopsis thaliana, were grown under sufficient and deficient Pi supply. Significant treatment effects were found for all traits and significant genotype x treatment interactions for the leaf Pi and sulphate concentrations. Using the RIL/QTL population, we identified 24 QTLs for leaf concentrations of Pi and other anions, including a major QTL for leaf sulphate concentration (SUL2) mapped to the bottom of chromosome (Chr) 1. GWA mapping found 188 SNPs to be associated with the measured traits, corresponding to 152 genes. One of these SNPs, associated with leaf Pi concentration, mapped to PP2A-1, a gene encoding an isoform of the catalytic subunit of a protein phosphatase 2A. Of two additional SNPs, associated with phosphate use efficiency (PUE), one mapped to AT5G49780, encoding a leucine-rich repeat protein kinase involved in signal transduction, and the other to SIZ1, a gene encoding a SUMO E3 ligase, and a known regulator of P starvation-dependent responses. One SNP associated with leaf sulphate concentration was found in SULTR2;1, encoding a sulphate transporter, known to enhance sulphate translocation from root to shoot under P deficiency. Finally, one SNP was mapped to FMO GS-OX4, a gene encoding glucosinolate S-oxygenase involved in glucosinolate biosynthesis, which located within the confidence interval of the SUL2 locus.

CONCLUSION

We identified several candidate genes with known functions related to anion homeostasis in response to Pi availability. Further molecular studies are needed to confirm and validate these candidate genes and understand their roles in examined traits. Such knowledge will contribute to future breeding for improved crop PUE .

Integrating transcriptomic network reconstruction and eQTL analyses reveals mechanistic connections between genomic architecture and Brassica rapa development

Plant developmental dynamics can be heritable, genetically correlated with fitness and yield, and undergo selection. Therefore, characterizing the mechanistic connections between the genetic architecture governing plant development and the resulting ontogenetic dynamics of plants in field settings is critically important for agricultural production and evolutionary ecology. We use hierarchical Bayesian Function-Valued Trait (FVT) models to estimate Brassica rapa growth curves throughout ontogeny, across two treatments, and in two growing seasons. We find genetic variation for plasticity of growth rates and final sizes, but not the inflection point (transition from accelerating to decelerating growth) of growth curves. There are trade-offs between growth rate and duration, indicating that selection for maximum yields at early harvest dates may come at the expense of late harvest yields and vice versa. We generate eigengene modules and determine which are co-expressed with FVT traits using a Weighted Gene Co-expression Analysis. Independently, we seed a Mutual Rank co-expression network model with FVT traits to identify specific genes and gene networks related to FVT. GO-analyses of eigengene modules indicate roles for actin/cytoskeletal genes, herbivore resistance/wounding responses, and cell division, while MR networks demonstrate a close association between metabolic regulation and plant growth. We determine that combining FVT Quantitative Trait Loci (QTL) and MR genes/WGCNA eigengene expression profiles better characterizes phenotypic variation than any single data type (i.e. QTL, gene, or eigengene alone). Our network analysis allows us to employ a targeted eQTL analysis, which we use to identify regulatory hotspots for FVT. We examine cis vs. trans eQTL that mechanistically link FVT QTL with structural trait variation. Colocalization of FVT, gene, and eigengene eQTL provide strong evidence for candidate genes influencing plant height. The study is the first to explore eQTL for FVT, and specifically do so in agroecologically relevant field settings.

We estimate the developmental dynamics of plant growth using mathematical functions to fit continuous functions to discrete plant height data collected throughout growth, and we use the parameters defining these mathematical functions as data. We identify genomic regions controlling plant growth and filter a novel transcriptomic data set using network reconstruction models to identify the genes and eigengenes associated with plant height. We combine these genomic and transcriptomic data to predict variation in plant height, and we use quantitative genetics to mechanistically connect plant genetics, transcriptomics, and development. Our approach demonstrates two powerful methods for the type of data reduction (FVT modeling and gene expression network reconstruction for targeted eQTL analyses) and data integration that will be necessary for driving forward the field of genetics in the post-genomic era. To the best of our knowledge, we are the first to apply these techniques to continuous models of plant development, and the first to do so in agroecologically relevant field settings.

Kernel size-related genes revealed by an integrated eQTL analysis during early maize kernel development

In maize, kernel traits strongly impact overall grain yields, and it is known that sophisticated spatiotemporal programs of gene expression coordinate kernel development, so advancing our knowledge of kernel development can help efforts to improve grain yields. Here, using phenotype, genotype and transcriptomics data of maize kernels at 5 and 15 days after pollination (DAP) for a large association mapping panel, we employed multiple quantitative genetics approaches-genome-wide association studies (GWAS) as well as expression quantitative trait loci (eQTL) and quantitative trait transcript (QTT) analyses-to gain insights about molecular genetic basis of kernel development in maize. This resulted in the identification of 137 putative kernel length-related genes at 5 DAP, of which 43 are located in previously reported QTL regions. Strikingly, we identified an eQTL that overlaps the locus encoding a maize homolog of the recently described m⁶ A methylation reader protein ECT2 from Arabidopsis; this putative ^epi eQTL is associated with 53 genes and may represent a master epi-transcriptomic regulator of kernel development. Notably, among the genes associated with this ^epi eQTL, 10 are for the main storage proteins in the maize endosperm (zeins) and two are known regulators of zein expression or endosperm development (Opaque2 and ZmICE1). Collectively, beyond cataloging and characterizing genomic attributes of a large number of eQTL associated with kernel development in maize, our study highlights how an eQTL approach can bolster the impact of both GWAS and QTT studies and can drive insights about the basic biology of plants.

Effect of phosphate nutrition on growth, physiology and phosphate transporter expression of cucumber seedlings

Although abundantly present in soils, inorganic phosphate (Pi) acquisition by plants is highly dependent on the transmembrane phosphate transporter (PT) gene family. Cucumber (Cucumis sativus) requires a large amount of phosphorus (P). The purpose of this study was to isolate the CsPT2-1 from cucumber roots, and to determine the influence of Pi nutrition on cucumber growth, metabolism and transcript levels of CsPT2-1 in tissues. Full length CsPT2-1 was cloned and phylogenetically identified. In two greenhouse experiments, P-deficient seedlings provided with low or high P concentrations were sampled at 10 and 21 days post treatment, respectively. Addition of P dramatically reduced growth of roots but not shoots. Supplying plants with high P resulted in increased total protein in leaves. Acid phosphatase activity increased significantly in leaves at any rate higher than 4 mM P. Increasing P concentration had a notable decrease in glucose concentrations in leaves of plants supplied with >0.5 mM P. In roots, glucose and starch concentrations increased with increasing P supply. Steady-state transcript levels of CsPT2-1 were high in P-deprived roots, but declined when plants were provided >10 mM P. To our knowledge, this is the first report focusing on a PT and its expression levels in cucumber.

Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum

Fusarium head blight is a devastating disease of small grain cereals such as bread wheat (Triticum aestivum). The pathogen switches from a biotrophic to a nectrotrophic lifestyle in course of disease development forcing its host to adapt its defence strategies. Using a genetical genomics approach, we illustrate genome-wide reconfigurations of genetic control over transcript abundances between two decisive time points after inoculation with the causative pathogen Fusarium graminearum. Whole transcriptome measurements have been recorded for 163 lines of a wheat doubled haploid population segregating for several resistance genes yielding 15 552 at 30 h and 15 888 eQTL at 50 h after inoculation. The genetic map saturated with transcript abundance-derived markers identified of a novel QTL on chromosome 6A, besides the previously reported QTL Fhb1 and Qfhs.ifa-5A. We find a highly different distribution of eQTL between time points with about 40% of eQTL being unique for the respective assessed time points. But also for more than 20% of genes governed by eQTL at either time point, genetic control changes in time. These changes are reflected in the dynamic compositions of three major regulatory hotspots on chromosomes 2B, 4A and 5A. In particular, control of defence-related biological mechanisms concentrated in the hotspot at 4A shift to hotspot 2B as the disease progresses. Hotspots do not colocalize with phenotypic QTL, and within their intervals no higher than expected number of eQTL was detected. Thus, resistance conferred by either QTL is mediated by few or single genes.

Using RNA-Seq for Genomic Scaffold Placement, Correcting Assemblies, and Genetic Map Creation in a Common Brassica rapa Mapping Population

Brassica rapa is a model species for agronomic, ecological, evolutionary, and translational studies. Here, we describe high-density SNP discovery and genetic map construction for a B. rapa recombinant inbred line (RIL) population derived from field collected RNA sequencing (RNA-Seq) data. This high-density genotype data enables the detection and correction of putative genome misassemblies and accurate assignment of scaffold sequences to their likely genomic locations. These assembly improvements represent 7.1-8.0% of the annotated B. rapa genome. We demonstrate how using this new resource leads to a significant improvement for QTL analysis over the current low-density genetic map. Improvements are achieved by the increased mapping resolution and by having known genomic coordinates to anchor the markers for candidate gene discovery. These new molecular resources and improvements in the genome annotation will benefit the Brassicaceae genomics community and may help guide other communities in fine-tuning genome annotations.

A Cross-Species Gene Expression Marker-Based Genetic Map and QTL Analysis in Bambara Groundnut

Bambara groundnut (Vigna subterranea (L.) Verdc.) is an underutilised legume crop, which has long been recognised as a protein-rich and drought-tolerant crop, used extensively in Sub-Saharan Africa. The aim of the study was to identify quantitative trait loci (QTL) involved in agronomic and drought-related traits using an expression marker-based genetic map based on major crop resources developed in soybean. The gene expression markers (GEMs) were generated at the (unmasked) probe-pair level after cross-hybridisation of bambara groundnut leaf RNA to the Affymetrix Soybean Genome GeneChip. A total of 753 markers grouped at an LOD (Logarithm of odds) of three, with 527 markers mapped into linkage groups. From this initial map, a spaced expression marker-based genetic map consisting of 13 linkage groups containing 218 GEMs, spanning 982.7 cM (centimorgan) of the bambara groundnut genome, was developed. Of the QTL detected, 46% were detected in both control and drought treatment populations, suggesting that they are the result of intrinsic trait differences between the parental lines used to construct the cross, with 31% detected in only one of the conditions. The present GEM map in bambara groundnut provides one technically feasible route for the translation of information and resources from major and model plant species to underutilised and resource-poor crops.

QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems

A high-density SNP-based genetic linkage map was constructed and integrated with a previous map in the Tapidor x Ningyou7 (TNDH) Brassica napus population, giving a new map with a total of 2041 molecular markers and an average marker density which increased from 0.39 to 0.97 (0.82 SNP bin) per cM. Root and shoot traits were screened under low and 'normal' phosphate (Pi) supply using a 'pouch and wick' system, and had been screened previously in an agar based system. The P-efficient parent Ningyou7 had a shorter primary root length (PRL), greater lateral root density (LRD) and a greater shoot biomass than the P-inefficient parent Tapidor under both treatments and growth systems. Quantitative trait loci (QTL) analysis identified a total of 131 QTL, and QTL meta-analysis found four integrated QTL across the growth systems. Integration reduced the confidence interval by ~41%. QTL for root and shoot biomass were co-located on chromosome A3 and for lateral root emergence were co-located on chromosomes A4/C4 and C8/C9. There was a major QTL for LRD on chromosome C9 explaining ~18% of the phenotypic variation. QTL underlying an increased LRD may be a useful breeding target for P uptake efficiency in Brassica.

eQTL Regulating Transcript Levels Associated with Diverse Biological Processes in Tomato

Variation in gene expression, in addition to sequence polymorphisms, is known to influence developmental, physiological, and metabolic traits in plants. Genetic mapping populations have facilitated identification of expression quantitative trait loci (eQTL), the genetic determinants of variation in gene expression patterns. We used an introgression population developed from the wild desert-adapted Solanum pennellii and domesticated tomato (Solanum lycopersicum) to identify the genetic basis of transcript level variation. We established the effect of each introgression on the transcriptome and identified approximately 7,200 eQTL regulating the steady-state transcript levels of 5,300 genes. Barnes-Hut t-distributed stochastic neighbor embedding clustering identified 42 modules revealing novel associations between transcript level patterns and biological processes. The results showed a complex genetic architecture of global transcript abundance pattern in tomato. Several genetic hot spots regulating a large number of transcript level patterns relating to diverse biological processes such as plant defense and photosynthesis were identified. Important eQTL regulating transcript level patterns were related to leaf number and complexity as well as hypocotyl length. Genes associated with leaf development showed an inverse correlation with photosynthetic gene expression, but eQTL regulating genes associated with leaf development and photosynthesis were dispersed across the genome. This comprehensive eQTL analysis details the influence of these loci on plant phenotypes and will be a valuable community resource for investigations on the genetic effects of eQTL on phenotypic traits in tomato.

Envirotyping for deciphering environmental impacts on crop plants

Global climate change imposes increasing impacts on our environments and crop production. To decipher environmental impacts on crop plants, the concept "envirotyping" is proposed, as a third "typing" technology, complementing with genotyping and phenotyping. Environmental factors can be collected through multiple environmental trials, geographic and soil information systems, measurement of soil and canopy properties, and evaluation of companion organisms. Envirotyping contributes to crop modeling and phenotype prediction through its functional components, including genotype-by-environment interaction (GEI), genes responsive to environmental signals, biotic and abiotic stresses, and integrative phenotyping. Envirotyping, driven by information and support systems, has a wide range of applications, including environmental characterization, GEI analysis, phenotype prediction, near-iso-environment construction, agronomic genomics, precision agriculture and breeding, and development of a four-dimensional profile of crop science involving genotype (G), phenotype (P), envirotype (E) and time (T) (developmental stage). In the future, envirotyping needs to zoom into specific experimental plots and individual plants, along with the development of high-throughput and precision envirotyping platforms, to integrate genotypic, phenotypic and envirotypic information for establishing a high-efficient precision breeding and sustainable crop production system based on deciphered environmental impacts.

A Systems Genetics Approach Identifies Gene Regulatory Networks Associated with Fatty Acid Composition in Brassica rapa Seed

Fatty acids in seeds affect seed germination and seedling vigor, and fatty acid composition determines the quality of seed oil. In this study, quantitative trait locus (QTL) mapping of fatty acid and transcript abundance was integrated with gene network analysis to unravel the genetic regulation of seed fatty acid composition in a Brassica rapa doubled haploid population from a cross between a yellow sarson oil type and a black-seeded pak choi. The distribution of major QTLs for fatty acids showed a relationship with the fatty acid types: linkage group A03 for monounsaturated fatty acids, A04 for saturated fatty acids, and A05 for polyunsaturated fatty acids. Using a genetical genomics approach, expression quantitative trait locus (eQTL) hotspots were found at major fatty acid QTLs on linkage groups A03, A04, A05, and A09. An eQTL-guided gene coexpression network of lipid metabolism-related genes showed major hubs at the genes BrPLA2-ALPHA, BrWD-40, a number of seed storage protein genes, and the transcription factor BrMD-2, suggesting essential roles for these genes in lipid metabolism. Three subnetworks were extracted for the economically important and most abundant fatty acids erucic, oleic, linoleic, and linolenic acids. Network analysis, combined with comparison of the genome positions of cis- or trans-eQTLs with fatty acid QTLs, allowed the identification of candidate genes for genetic regulation of these fatty acids. The generated insights in the genetic architecture of fatty acid composition and the underlying complex gene regulatory networks in B. rapa seeds are discussed.

Crop epigenetics and the molecular hardware of genotype × environment interactions

Crop plants encounter thermal environments which fluctuate on a diurnal and seasonal basis. Future climate resilient cultivars will need to respond to thermal profiles reflecting more variable conditions, and harness plasticity that involves regulation of epigenetic processes and complex genomic regulatory networks. Compartmentalization within plant cells insulates the genomic central processing unit within the interphase nucleus. This review addresses the properties of the chromatin hardware in which the genome is embedded, focusing on the biophysical and thermodynamic properties of DNA, histones and nucleosomes. It explores the consequences of thermal and ionic variation on the biophysical behavior of epigenetic marks such as DNA cytosine methylation (5mC), and histone variants such as H2A.Z, and how these contribute to maintenance of chromatin integrity in the nucleus, while enabling specific subsets of genes to be regulated. Information is drawn from theoretical molecular in vitro studies as well as model and crop plants and incorporates recent insights into the role epigenetic processes play in mediating between environmental signals and genomic regulation. A preliminary speculative framework is outlined, based on the evidence of what appears to be a cohesive set of interactions at molecular, biophysical and electrostatic level between the various components contributing to chromatin conformation and dynamics. It proposes that within plant nuclei, general and localized ionic homeostasis plays an important role in maintaining chromatin conformation, whilst maintaining complex genomic regulation that involves specific patterns of epigenetic marks. More generally, reversible changes in DNA methylation appear to be consistent with the ability of nuclear chromatin to manage variation in external ionic and temperature environment. Whilst tentative, this framework provides scope to develop experimental approaches to understand in greater detail the internal environment of plant nuclei. It is hoped that this will generate a deeper understanding of the molecular mechanisms underlying genotype × environment interactions that may be beneficial for long-term improvement of crop performance in less predictable climates.

RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

BACKGROUND

The evolutionary basis of reproductive success in different environments is of major interest in the study of plant adaptation. Since the reproductive stage is particularly sensitive to drought, genes affecting reproductive success during this stage are key players in the evolution of adaptive mechanisms. We used an ecological genomics approach to investigate the reproductive response of drought-tolerant and sensitive wild barley accessions originating from different habitats in the Levant.

RESULTS

We sequenced mRNA extracted from spikelets at the flowering stage in drought-treated and control plants. The barley genome was used for a reference-guided assembly and differential expression analysis. Our approach enabled to detect biological processes affecting grain production under drought stress. We detected novel candidate genes and differentially expressed alleles associated with drought tolerance. Drought associated genes were shown to be more conserved than non-associated genes, and drought-tolerance genes were found to evolve more rapidly than other drought associated genes.

CONCLUSIONS

We show that reproductive success under drought stress is not a habitat-specific trait but a shared physiological adaptation that appeared to evolve recently in the evolutionary history of wild barley. Exploring the genomic basis of reproductive success under stress in crop wild progenitors is expected to have considerable ecological and economical applications.

The temporal transcriptomic response of Pinus massoniana seedlings to phosphorus deficiency

BACKGROUND

Phosphorus (P) is an essential macronutrient for plant growth and development. Several genes involved in phosphorus deficiency stress have been identified in various plant species. However, a whole genome understanding of the molecular mechanisms involved in plant adaptations to low P remains elusive, and there is particularly little information on the genetic basis of these acclimations in coniferous trees. Masson pine (Pinus massoniana) is grown mainly in the tropical and subtropical regions in China, many of which are severely lacking in inorganic phosphate (Pi). In previous work, we described an elite P. massoniana genotype demonstrating a high tolerance to Pi-deficiency.

METHODOLOGY/PRINCIPAL FINDINGS

To further investigate the mechanism of tolerance to low P, RNA-seq was performed to give an idea of extent of expression from the two mixed libraries, and microarray whose probes were designed based on the unigenes obtained from RNA-seq was used to elucidate the global gene expression profiles for the long-term phosphorus starvation. A total of 70,896 unigenes with lengths ranging from 201 to 20,490 bp were assembled from 112,108,862 high quality reads derived from RNA-Seq libraries. We identified 1,396 and 943 transcripts that were differentially regulated (P<0.05) under P1 (0.01 mM P) and P2 (0.06 mM P) Pi-deficiency conditions, respectively. Numerous transcripts were consistently differentially regulated under Pi deficiency stress, many of which were also up- or down-regulated in other species under the corresponding conditions, and are therefore ideal candidates for monitoring the P status of plants. The results also demonstrated the impact of different Pi starvation levels on global gene expression in Masson pine.

CONCLUSIONS/SIGNIFICANCE

To our knowledge, this work provides the first insight into the molecular mechanisms involved in acclimation to long-term Pi starvation and different Pi availability levels in coniferous trees.

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. FST outlier analysis revealed that promoter, 5'-UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within-population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Sustainable harvest: managing plasticity for resilient crops

Maintaining crop production to feed a growing world population is a major challenge for this period of rapid global climate change. No consistent conceptual or experimental framework for crop plants integrates information at the levels of genome regulation, metabolism, physiology and response to growing environment. An important role for plasticity in plants is assisting in homeostasis in response to variable environmental conditions. Here, we outline how plant plasticity is facilitated by epigenetic processes that modulate chromatin through dynamic changes in DNA methylation, histone variants, small RNAs and transposable elements. We present examples of plant plasticity in the context of epigenetic regulation of developmental phases and transitions and map these onto the key stages of crop establishment, growth, floral initiation, pollination, seed set and maturation of harvestable product. In particular, we consider how feedback loops of environmental signals and plant nutrition affect plant ontogeny. Recent advances in understanding epigenetic processes enable us to take a fresh look at the crosstalk between regulatory systems that confer plasticity in the context of crop development. We propose that these insights into genotype × environment (G × E) interaction should underpin development of new crop management strategies, both in terms of information-led agronomy and in recognizing the role of epigenetic variation in crop breeding.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities

BACKGROUND AND AIMS

Phosphate (Pi) deficiency in soils is a major limiting factor for crop growth worldwide. Plant growth under low Pi conditions correlates with root architectural traits and it may therefore be possible to select these traits for crop improvement. The aim of this study was to characterize root architectural traits, and to test quantitative trait loci (QTL) associated with these traits, under low Pi (LP) and high Pi (HP) availability in Brassica napus.

METHODS

Root architectural traits were characterized in seedlings of a double haploid (DH) mapping population (n = 190) of B. napus ['Tapidor' × 'Ningyou 7' (TNDH)] using high-throughput phenotyping methods. Primary root length (PRL), lateral root length (LRL), lateral root number (LRN), lateral root density (LRD) and biomass traits were measured 12 d post-germination in agar at LP and HP.

KEY RESULTS

In general, root and biomass traits were highly correlated under LP and HP conditions. 'Ningyou 7' had greater LRL, LRN and LRD than 'Tapidor', at both LP and HP availability, but smaller PRL. A cluster of highly significant QTL for LRN, LRD and biomass traits at LP availability were identified on chromosome A03; QTL for PRL were identified on chromosomes A07 and C06.

CONCLUSIONS

High-throughput phenotyping of Brassica can be used to identify root architectural traits which correlate with shoot biomass. It is feasible that these traits could be used in crop improvement strategies. The identification of QTL linked to root traits under LP and HP conditions provides further insights on the genetic basis of plant tolerance to P deficiency, and these QTL warrant further dissection.

Natural variation and genetic constraints on drought tolerance

Drought is a central abiotic stress for both natural plant populations and agricultural crops. Substantial natural genetic variation in drought resistance traits has been identified in plant populations, crop species, and laboratory model systems. In particular, studies in Arabidopsis thaliana have discovered variation in a number of key physiological traits involved in plant-water relations that may underlie evolved drought stress responses among accessions. Despite this abundant variation, we still know little about the complex genetic architecture of drought tolerance or its role in constraining evolution. Unfortunately, few natural allelic variants have been cloned for drought related traits--progress cloning QTL, the use of RNA-sequencing methods for evaluating gene expression responses to soil drying, and improved methodology for exploring complex multivariate data all hold promise for moving the field forward. In particular, a better understanding of the molecular nature of pleiotropic gene action and the genetics of phenotypic plasticity will give insight into local adaptation in plants and provide new avenues for improving crops.

QTL for yield traits and their association with functional genes in response to phosphorus deficiency in Brassica napus

BACKGROUND

Oilseed rape (Brassica napus L.) is one of the most important oil crops. A primary limitation to the cultivation of this crop is the lack of available phosphorus (P) in soils. To elucidate the genetic control of P deficiency tolerance in Brassica napus, quantitative trait locus (QTL) for seed yield and yield related-traits in response to P deficiency were identified using a double haploid mapping population (TN DH) derived from a cross between a P-efficient cultivar, Ningyou 7 and a P-inefficient cultivar, Tapidor.

RESULTS

Three field trials were conducted to determine seed yield (SY), plant height (PH), number of primary branches (BN), height to the first primary branch (FBH), relative first primary branch height (RBH), pod number per plant (PN), seed number per pod (SN) and seed weight of 1,000 seeds (SW) in 188 lines of TN DH population exposed to low P (LP) and optimal P (OP) conditions. P deficiency decreased PH, BN, SN, PN and SY, and increased FBH and RBH with no effect on SW. Three reproducible LP-specific QTL regions were identified on chromosomes A2, A3 and A5 that controlled SN, PN and SW respectively. In addition, six reproducible constitutive regions were also mapped with two each for SY-LP on A2, and FBH-LP on C6 and one each for PH-LP and SW-LP on A3. About 30 markers derived from 19 orthologous genes involved in Arabidopsis P homeostasis were mapped on 24 QTL regions by comparative mapping between Arabidopsis and Brassica napus. Among these genes, GPT1, MGD2 and SIZ1 were associated with two major loci regulating SY-LP and other yield-related traits on A2 between 77.1 and 95.0 cM.

CONCLUSION

The stable QTLs detected under LP conditions and their candidate genes may provide useful information for marker-assisted selection in breeding high-yield B. napus varieties with improved P efficiency.

Opportunities for improving phosphorus-use efficiency in crop plants

Limitation of grain crop productivity by phosphorus (P) is widespread and will probably increase in the future. Enhanced P efficiency can be achieved by improved uptake of phosphate from soil (P-acquisition efficiency) and by improved productivity per unit P taken up (P-use efficiency). This review focuses on improved P-use efficiency, which can be achieved by plants that have overall lower P concentrations, and by optimal distribution and redistribution of P in the plant allowing maximum growth and biomass allocation to harvestable plant parts. Significant decreases in plant P pools may be possible, for example, through reductions of superfluous ribosomal RNA and replacement of phospholipids by sulfolipids and galactolipids. Improvements in P distribution within the plant may be possible by increased remobilization from tissues that no longer need it (e.g. senescing leaves) and reduced partitioning of P to developing grains. Such changes would prolong and enhance the productive use of P in photosynthesis and have nutritional and environmental benefits. Research considering physiological, metabolic, molecular biological, genetic and phylogenetic aspects of P-use efficiency is urgently needed to allow significant progress to be made in our understanding of this complex trait.

InterStoreDB: a generic integration resource for genetic and genomic data

Associating phenotypic traits and quantitative trait loci (QTL) to causative regions of the underlying genome is a key goal in agricultural research. InterStoreDB is a suite of integrated databases designed to assist in this process. The individual databases are species independent and generic in design, providing access to curated datasets relating to plant populations, phenotypic traits, genetic maps, marker loci and QTL, with links to functional gene annotation and genomic sequence data. Each component database provides access to associated metadata, including data provenance and parameters used in analyses, thus providing users with information to evaluate the relative worth of any associations identified. The databases include CropStoreDB, for management of population, genetic map, QTL and trait measurement data, SeqStoreDB for sequence-related data and AlignStoreDB, which stores sequence alignment information, and allows navigation between genetic and genomic datasets. Genetic maps are visualized and compared using the CMAP tool, and functional annotation from sequenced genomes is provided via an EnsEMBL-based genome browser. This framework facilitates navigation of the multiple biological domains involved in genetics and genomics research in a transparent manner within a single portal. We demonstrate the value of InterStoreDB as a tool for Brassica research. InterStoreDB is available from: http://www.interstoredb.org.

Diversity in expression of phosphorus (P) responsive genes in Cucumis melo L

BACKGROUND

Phosphorus (P) is a major limiting nutrient for plant growth in many soils. Studies in model species have identified genes involved in plant adaptations to low soil P availability. However, little information is available on the genetic bases of these adaptations in vegetable crops. In this respect, sequence data for melon now makes it possible to identify melon orthologues of candidate P responsive genes, and the expression of these genes can be used to explain the diversity in the root system adaptation to low P availability, recently observed in this species.

METHODOLOGY AND FINDINGS

Transcriptional responses to P starvation were studied in nine diverse melon accessions by comparing the expression of eight candidate genes (Cm-PAP10.1, Cm-PAP10.2, Cm-RNS1, Cm-PPCK1, Cm-transferase, Cm-SQD1, Cm-DGD1 and Cm-SPX2) under P replete and P starved conditions. Differences among melon accessions were observed in response to P starvation, including differences in plant morphology, P uptake, P use efficiency (PUE) and gene expression. All studied genes were up regulated under P starvation conditions. Differences in the expression of genes involved in P mobilization and remobilization (Cm-PAP10.1, Cm-PAP10.2 and Cm-RNS1) under P starvation conditions explained part of the differences in P uptake and PUE among melon accessions. The levels of expression of the other studied genes were diverse among melon accessions, but contributed less to the phenotypical response of the accessions.

CONCLUSIONS

This is the first time that these genes have been described in the context of P starvation responses in melon. There exists significant diversity in gene expression levels and P use efficiency among melon accessions as well as significant correlations between gene expression levels and phenotypical measurements.

Lessons from eQTL mapping studies: non-coding regions and their role behind natural phenotypic variation in plants

Even if considerable progress has been achieved towards the understanding of natural variation in plant systems, the contribution of transcript abundance variation to phenotypic diversity remains unappreciated. Over the last decade, efforts to characterise the genome-wide expression variation in natural accessions, structured populations and hybrids have improved our knowledge of the contribution of non-coding polymorphisms to gene expression regulation. Moreover, new studies are helping to unravel the role of expression polymorphisms and their orchestrated performance. Recent advances involving classical linkage analysis, GWAS and improved eQTL mapping strategies will provide a greater resolution to determine the genetic variants shaping the broad diversity in plant systems.

Phosphate Import in Plants: Focus on the PHT1 Transporters

The main source of phosphorus for plants is inorganic phosphate (Pi), which is characterized by its poor availability and low mobility. Uptake of this element from the soil relies heavily upon the PHT1 transporters, a specific family of plant plasma membrane proteins that were identified by homology with the yeast PHO84 Pi transporter. Since the discovery of PHT1 transporters in 1996, various studies have revealed that their function is controlled by a highly complex network of regulation. This review will summarize the current state of research on plant PHT1 multigenic families, including physiological, biochemical, molecular, cellular, and genetics studies.

Gene expression changes in phosphorus deficient potato (Solanum tuberosum L.) leaves and the potential for diagnostic gene expression markers

BACKGROUND

There are compelling economic and environmental reasons to reduce our reliance on inorganic phosphate (Pi) fertilisers. Better management of Pi fertiliser applications is one option to improve the efficiency of Pi fertiliser use, whilst maintaining crop yields. Application rates of Pi fertilisers are traditionally determined from analyses of soil or plant tissues. Alternatively, diagnostic genes with altered expression under Pi limiting conditions that suggest a physiological requirement for Pi fertilisation, could be used to manage Pifertiliser applications, and might be more precise than indirect measurements of soil or tissue samples.

RESULTS

We grew potato (Solanum tuberosum L.) plants hydroponically, under glasshouse conditions, to control their nutrient status accurately. Samples of total leaf RNA taken periodically after Pi was removed from the nutrient solution were labelled and hybridised to potato oligonucleotide arrays. A total of 1,659 genes were significantly differentially expressed following Pi withdrawal. These included genes that encode proteins involved in lipid, protein, and carbohydrate metabolism, characteristic of Pi deficient leaves and included potential novel roles for genes encoding patatin like proteins in potatoes. The array data were analysed using a support vector machine algorithm to identify groups of genes that could predict the Pi status of the crop. These groups of diagnostic genes were tested using field grown potatoes that had either been fertilised or unfertilised. A group of 200 genes could correctly predict the Pi status of field grown potatoes.

CONCLUSIONS

This paper provides a proof-of-concept demonstration for using microarrays and class prediction tools to predict the Pi status of a field grown potato crop. There is potential to develop this technology for other biotic and abiotic stresses in field grown crops. Ultimately, a better understanding of crop stresses may improve our management of the crop, improving the sustainability of agriculture.

Expression level polymorphisms: heritable traits shaping natural variation

Natural accessions of many species harbor a wealth of genetic variation visible in a large array of phenotypes. Although expression level polymorphisms (ELPs) in several genes have been shown to contribute to variation in diverse traits, their general impact on adaptive variation has likely been underestimated. At present, ELPs have predominantly been correlated to quantitative trait loci (eQTLs) that occupy central hubs in signaling networks, which pleiotropically affect numerous traits. To increase the sensitivity of detecting minor effect eQTLs or those that act in a trait-specific manner, we emphasize the need for more systematic approaches. This requires, but is not limited to, refining experimental designs such as reduction of tissue complexity and combinatorial methods including a priori defined networks.