Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

Polygenic architecture of flowering time and its relationship with local environments in the grass Brachypodium distachyon

Synchronizing the timing of reproduction with the environment is crucial in the wild. Among the multiple mechanisms, annual plants evolved to sense their environment, the requirement of cold-mediated vernalization is a major process that prevents individuals from flowering during winter. In many annual plants including crops, both a long and short vernalization requirement can be observed within species, resulting in so-called early-(spring) and late-(winter) flowering genotypes. Here, using the grass model Brachypodium distachyon, we explored the link between flowering-time-related traits (vernalization requirement and flowering time), environmental variation, and diversity at flowering-time genes by combining measurements under greenhouse and outdoor conditions. These experiments confirmed that B. distachyon natural accessions display large differences regarding vernalization requirements and ultimately flowering time. We underline significant, albeit quantitative effects of current environmental conditions on flowering-time-related traits. While disentangling the confounding effects of population structure on flowering-time-related traits remains challenging, population genomics analyses indicate that well-characterized flowering-time genes may contribute significantly to flowering-time variation and display signs of polygenic selection. Flowering-time genes, however, do not colocalize with genome-wide association peaks obtained with outdoor measurements, suggesting that additional genetic factors contribute to flowering-time variation in the wild. Altogether, our study fosters our understanding of the polygenic architecture of flowering time in a natural grass system and opens new avenues of research to investigate the gene-by-environment interaction at play for this trait.

Artificial selection reveals complex genetic architecture of shoot branching and its response to nitrate supply in Arabidopsis

Quantitative traits may be controlled by many loci, many alleles at each locus, and subject to genotype-by-environment interactions, making them difficult to map. One example of such a complex trait is shoot branching in the model plant Arabidopsis, and its plasticity in response to nitrate. Here, we use artificial selection under contrasting nitrate supplies to dissect the genetic architecture of this complex trait, where loci identified by association mapping failed to explain heritability estimates. We found a consistent response to selection for high branching, with correlated responses in other traits such as plasticity and flowering time. Genome-wide scans for selection and simulations suggest that at least tens of loci control this trait, with a distinct genetic architecture between low and high nitrate treatments. While signals of selection could be detected in the populations selected for high branching on low nitrate, there was very little overlap in the regions selected in three independent populations. Thus the regulatory network controlling shoot branching can be tuned in different ways to give similar phenotypes.

Construction of a high-density genetic map and dissection of genetic architecture of six agronomic traits in tobacco (Nicotiana tabacum L.)

Tobacco (Nicotiana tabacum L.) is an economic crop and a model organism for studies on plant biology and genetics. A population of 271 recombinant inbred lines (RIL) derived from K326 and Y3, two elite flue-cured tobacco parents, has been constructed to investigate the genetic basis of agronomic traits in tobacco. Six agronomic traits including natural plant height (nPH), natural leaf number (nLN), stem girth (SG), inter-node length (IL), length of the largest leaf (LL) and width of the largest leaf (LW) were measured in seven environments, spanning the period between 2018 and 2021. We firstly developed an integrated SNP-indel-SSR linkage map with 43,301 SNPs, 2,086 indels and 937 SSRs, which contained 7,107 bin markers mapped on 24 LGs and covered 3334.88 cM with an average genetic distance of 0.469cM. Based on this high-density genetic map, a total of 70 novel QTLs were detected for six agronomic traits by a full QTL model using the software QTLNetwork, of which 32 QTLs showed significant additive effects, 18 QTLs showed significant additive-by-environment interaction effects, 17 pairs showed significant additive-by-additive epistatic effects and 13 pairs showed significant epistasis-by-environment interaction effects. In addition to additive effect as a major contributor to genetic variation, both epistasis effects and genotype-by-environment interaction effects played an important role in explaining phenotypic variation for each trait. In particular, qnLN6-1 was detected with considerably large main effect and high heritability ( h a 2 =34.80%). Finally, four genes including Nt16g00284.1, Nt16g00767.1, Nt16g00853.1, Nt16g00877.1 were predicted as pleiotropic candidate genes for five traits.

Epistatic interaction effect between chromosome 1BL (Yr29) and a novel locus on 2AL facilitating resistance to stripe rust in Chinese wheat Changwu 357-9

Four stable QTL for adult plant resistance were identified in wheat line Changwu 357-9, including a new QTL on 2AL showing significant interaction with Yr29 to reduce stripe rust severity. Stripe rust (yellow rust) is a serious disease of bread wheat (Triticum aestivum L.) worldwide. Genetic resistance is considered the most economical, effective and environmentally friendly method to control the disease and to minimize the use of fungicides. The current study focused on characterizing the components of stripe rust resistance and understanding the interactions in Changwu 357-9 (CW357-9)/Avocet S RIL population. A genetic linkage map constructed using a new GenoBaits Wheat 16K Panel and the 660K SNP array had 5104 polymorphic SNP markers spanning 3533.11 cM. Four stable QTL, consistently identified across five environments, were detected on chromosome arms 1BL, 2AL, 3DS, and 6BS in Changwu357-9. The most effective QTL QYrCW357-1BL was Yr29. The 6BS QTL was identified as Yr78, which has been combined with the 1BL QTL in many wheat cultivars and breeding lines. The novel QTL on 2AL with moderate effect showed a stable and significant epistatic interaction with Yr29. The QTL on 3DL should be same as QYrsn.nwafu-3DL and enriches the overall stripe rust resistance gene pool for breeding. Polymorphisms of flanking AQP markers AX-110020417 (for QYrCW357-1BL), AX-110974948 (for QYrCW357-2AL), AX-109466386 (for QYrCW357-3DL), and AX-109995005 (for QYrCW357-6BS) were evaluated in a diversity panel including 225 wheat cultivars and breeding lines. These results suggested that these high-throughput markers could be used to introduce QYrCW357-1BL, QYrCW357-2AL, QYrCW357-3DL, and QYrCW357-6BS into commercial wheat cultivars. Combinations of these genes with other APR QTL should lead to higher levels of stripe rust resistance along with the beneficial effects of multi-disease resistance gene Yr29 on improving resistance to other diseases.

Reconstruction of the High Stigma Exsertion Rate Trait in Rice by Pyramiding Multiple QTLs

Asian cultivated rice is a self-pollinating crop, which has already lost some traits of natural outcrossing in the process of domestication. However, male sterility lines (MSLs) need to have a strong outcrossing ability to produce hybrid seeds by outcrossing with restorer lines of male parents in hybrid rice seed production. Stigma exsertion rate (SER) is a trait related to outcrossing ability. Reconstruction of the high-SER trait is essential in the MSL breeding of rice. In previous studies, we detected eighteen quantitative trait loci (QTLs) for SER from Oryza sativa, Oryza glaberrima, and Oryza glumaepatula using single-segment substitution lines (SSSLs) in the genetic background of Huajingxian 74 (HJX74). In this study, eleven of the QTLs were used to develop pyramiding lines. A total of 29 pyramiding lines with 2-6 QTLs were developed from 10 SSSLs carrying QTLs for SER in the HJX74 genetic background. The results showed that the SER increased with increasing QTLs in the pyramiding lines. The SER in the lines with 5-6 QTLs was as high as wild rice with strong outcrossing ability. The epistasis of additive by additive interaction between QTLs in the pyramiding lines was less-than-additive or negative effect. One QTL, qSER3a-sat, showed minor-effect epistasis and increased higher SER than other QTLs in pyramiding lines. The detection of epistasis of QTLs on SER uncovered the genetic architecture of SER, which provides a basis for using these QTLs to improve SER levels in MSL breeding. The reconstruction of the high-SER trait will help to develop the MSLs with strong outcrossing ability in rice.

Cold stress and freezing tolerance negatively affect the fitness of Arabidopsis thaliana accessions under field and controlled conditions

Higher acclimated freezing tolerance improved winter survival, but reduced reproductive fitness of Arabidopsis thaliana accessions under field and controlled conditions. Low temperature is one of the most important abiotic factors influencing plant fitness and geographical distribution. In addition, cold stress is known to influence crop yield and is therefore of great economic importance. Increased freezing tolerance can be acquired by the process of cold acclimation, but this may be associated with a fitness cost. To assess the influence of cold stress on the fitness of plants, long-term field trials over 5 years were performed with six natural accessions of Arabidopsis thaliana ranging from very tolerant to very sensitive to freezing. Fitness parameters, as seed yield and 1000 seed mass, were measured and correlation analyses with temperature and freezing tolerance data performed. The results were compared with fitness parameters from controlled chamber experiments over 3 years with application of cold priming and triggering conditions. Winter survival and seed yield per plant were positively correlated with temperature in field experiments. In addition, winter survival and 1000 seed mass were correlated with the cold-acclimated freezing tolerance of the selected Arabidopsis accessions. The results provide strong evidence for a trade-off between higher freezing tolerance and reproductive fitness in A. thaliana, which might have ecological impacts in the context of global warming.

Comparative proteomics of thylakoids from Arabidopsis grown in laboratory and field conditions

Compared to controlled laboratory conditions, plant growth in the field is rarely optimal since it is frequently challenged by large fluctuations in light and temperature which lower the efficiency of photosynthesis and lead to photo-oxidative stress. Plants grown under natural conditions therefore place an increased onus on the regulatory mechanisms that protect and repair the delicate photosynthetic machinery. Yet, the exact changes in thylakoid proteome composition which allow plants to acclimate to the natural environment remain largely unexplored. Here, we use quantitative label-free proteomics to demonstrate that field-grown Arabidopsis plants incorporate aspects of both the low and high light acclimation strategies previously observed in laboratory-grown plants. Field plants showed increases in the relative abundance of ATP synthase, cytochrome b 6 f, ferredoxin-NADP+ reductases (FNR1 and FNR2) and their membrane tethers TIC62 and TROL, thylakoid architecture proteins CURT1A, CURT1B, RIQ1, and RIQ2, the minor monomeric antenna complex CP29.3, rapidly-relaxing non-photochemical quenching (qE)-related proteins PSBS and VDE, the photosystem II (PSII) repair machinery and the cyclic electron transfer complexes NDH, PGRL1B, and PGR5, in addition to decreases in the amounts of LHCII trimers composed of LHCB1.1, LHCB1.2, LHCB1.4, and LHCB2 proteins and CP29.2, all features typical of a laboratory high light acclimation response. Conversely, field plants also showed increases in the abundance of light harvesting proteins LHCB1.3 and CP29.1, zeaxanthin epoxidase (ZEP) and the slowly-relaxing non-photochemical quenching (qI)-related protein LCNP, changes previously associated with a laboratory low light acclimation response. Field plants also showed distinct changes to the proteome including the appearance of stress-related proteins ELIP1 and ELIP2 and changes to proteins that are largely invariant under laboratory conditions such as state transition related proteins STN7 and TAP38. We discuss the significance of these alterations in the thylakoid proteome considering the unique set of challenges faced by plants growing under natural conditions.

Determining dependency and redundancy for identifying gene-gene interaction associated with complex disease

As interactions among genetic variants in different genes can be an important factor for predicting complex diseases, many computational methods have been proposed to detect if a particular set of genes has interaction with a particular complex disease. However, even though many such methods have been shown to be useful, they can be made more effective if the properties of gene-gene interactions can be better understood. Towards this goal, we have attempted to uncover patterns in gene-gene interactions and the patterns reveal an interesting property that can be reflected in an inequality that describes the relationship between two genotype variables and a disease-status variable. We show, in this paper, that this inequality can be generalized to [Formula: see text] genotype variables. Based on this inequality, we establish a conditional independence and redundancy (CIR)-based definition of gene-gene interaction and the concept of an interaction group. From these new definitions, a novel measure of gene-gene interaction is then derived. We discuss the properties of these concepts and explain how they can be used in a novel algorithm to detect high-order gene-gene interactions. Experimental results using both simulated and real datasets show that the proposed method can be very promising.

A Nested Association Mapping Panel in Arabidopsis thaliana for Mapping and Characterizing Genetic Architecture

Linkage and association mapping populations are crucial public resources that facilitate the characterization of trait genetic architecture in natural and agricultural systems. We define a large nested association mapping panel (NAM) from 14 publicly available recombinant inbred line populations (RILs) of Arabidopsis thaliana, which share a common recurrent parent (Col-0). Using a genotype-by-sequencing approach (GBS), we identified single nucleotide polymorphisms (SNPs; range 563-1525 per population) and subsequently built updated linkage maps in each of the 14 RIL sets. Simulations in individual RIL populations indicate that our GBS markers have improved power to detect small effect QTL and enhanced resolution of QTL support intervals in comparison to original linkage maps. Using these robust linkage maps, we imputed a common set of publicly available parental SNPs into each RIL linkage map, generating overlapping markers across all populations. Though ultimately depending on allele frequencies at causal loci, simulations of the NAM panel suggest that surveying between 4 to 7 of the 14 RIL populations provides high resolution of the genetic architecture of complex traits, relative to a single mapping population.

Revealing biomass heterosis in the allodiploid xBrassicoraphanus, a hybrid between Brassica rapa and Raphanus sativus, through integrated transcriptome and metabolites analysis

BACKGROUND

Heterosis is biologically important but the molecular basis of the phenomenon is poorly understood. We characterized intergeneric hybrids between B. rapa cv. Chiifu and R. sativus cv. WK10039 as an extreme example of heterosis. Taking advantage of clear heterosis phenotypes and the genetic distance between parents, we performed transcriptome and metabolite analysis to decipher the molecular basis of heterosis.

RESULTS

The heterosis was expressed as fresh weight in the field and as inflorescence stem length in the glass house. Flowering time, distributed as a normal segregating population, ranged from the early flowering of one parent to the late flowering of the other, in contrast to the homogeneous flowering time in a typical F1 population, indicating unstable allelic interactions. The transcriptome and metabolome both indicated that sugar metabolism was altered, suggesting that the change in metabolism was linked to the heterosis. Because alleles were not shared between the hybridized genomes, classic models only partly explain this heterosis, indicating that other mechanisms are involved.

CONCLUSION

The differential expression of genes for primary and secondary metabolism, along with the altered metabolite profiles, suggests that heterosis could involve a change in balance between primary and secondary metabolism.

Assessing trait contribution and mapping novel QTL for salinity tolerance using the Bangladeshi rice landrace Capsule

BACKGROUND

Salinity is one of the most widespread abiotic stresses affecting rice productivity worldwide. The purpose of this study was to establish the relative importance of different traits associated with salinity tolerance in rice and to identify new quantitative trait loci (QTL) conferring tolerance to salinity at seedling stage. A total of 231 F_2:3 plants derived from a cross between a sensitive variety BRRI dhan29 (BR29 hereafter) and Capsule, a salt tolerant Bangladeshi indica landrace, were evaluated under salt stress in a phytotron.

RESULTS

Out of the 231 F₂ plants, 47 highly tolerant and 47 most sensitive lines were selected, representing the two extreme tails of the phenotypic distribution. These 94 plants were genotyped for 105 simple sequence repeat (SSR) and insertion/deletion (InDel) markers. A genetic linkage map spanning approximately 1442.9 cM of the 12 linkage groups with an average marker distance of 13.7 cM was constructed. QTL were identified on the long arm of chromosome 1 for Na⁺ concentration, K⁺ concentration, Na⁺-K⁺ ratio and survival; chromosome 3 for Na⁺ concentration, survival and overall phenotypic evaluation using the Standard Evaluation system (SES); and chromosome 5 for SES. A total of 6 pairwise epistatic interactions were also detected between QTL-linked and QTL-unlinked regions. Graphical genotyping indicated an association between the phenotypes of the extreme families and their QTL genotypes. Path coefficient analysis revealed that Na⁺ concentration, survival, Na⁺-K⁺ ratio and the overall phenotypic performance (SES score) are the major traits associated with salinity tolerance of Capsule.

CONCLUSIONS

Capsule provides an alternative source of salinity tolerance aside from Pokkali and Nona Bokra, the two Indian salt tolerant landraces traditionally used for breeding salt tolerant rice varieties. Pyramiding the new QTL identified in this study with previously discovered loci, such as Saltol, will facilitate breeding varieties that are highly tolerant of salt stress.

Population climatic history predicts phenotypic responses in novel environments for Arabidopsis thaliana in North America

PREMISE

Determining how species perform in novel climatic environments is essential for understanding (1) responses to climate change and (2) evolutionary consequences of biological invasions. For the vast majority of species, the number of population characteristics that will predict performance and patterns of natural selection in novel locations in the wild remains limited.

METHODS

We evaluated phenological, vegetative, architectural, and fitness-related traits in experimental gardens in contrasting climates (Ontario, Canada, and South Carolina, USA) in the North American non-native distribution of Arabidopsis thaliana. We assessed the effects of climatic distance, geographic distance, and genetic features of history on performance and patterns of natural selection in the novel garden settings.

RESULTS

We found that plants had greater survivorship, flowered earlier, were larger, and produced more fruit in the south, and that genotype-by-environment interactions were significant between gardens. However, our analyses revealed similar patterns of natural selection between gardens in distinct climate zones. After accounting for genetic ancestry, we also detected that population climatic distance best predicted performance within gardens.

CONCLUSIONS

These data suggest that colonization success in novel, non-native environments is determined by a combination of climate and genetic history. When performance at novel sites was assessed with seed sources from geographically and genetically disparate, established non-native populations, proximity to the garden alone was insufficient to predict performance. Our study highlights the need to evaluate seed sources from diverse origins to describe comprehensively phenotypic responses to novel environments, particularly for taxa in which many source populations may contribute to colonization.

A Low Resolution Epistasis Mapping Approach To Identify Chromosome Arm Interactions in Allohexaploid Wheat

Epistasis is an important contributor to genetic variance. In inbred populations, pairwise epistasis is present as additive by additive interactions. Testing for epistasis presents a multiple testing problem as the pairwise search space for modest numbers of markers is large. Single markers do not necessarily track functional units of interacting chromatin as well as haplotype based methods do. To harness the power of multiple markers while minimizing the number of tests conducted, we present a low resolution test for epistatic interactions across whole chromosome arms. Epistasis covariance matrices were constructed from the additive covariances of individual chromosome arms. These covariances were subsequently used to estimate an epistatic variance parameter while correcting for background additive and epistatic effects. We find significant epistasis for 2% of the interactions tested for four agronomic traits in a winter wheat breeding population. Interactions across homeologous chromosome arms were identified, but were less abundant than other chromosome arm pair interactions. The homeologous chromosome arm pair 4BL/4DL showed a strong negative relationship between additive and interaction effects that may be indicative of functional redundancy. Several chromosome arms appeared to act as hubs in an interaction network, suggesting that they may contain important regulatory factors. The differential patterns of epistasis across different traits demonstrate that detection of epistatic interactions is robust when correcting for background additive and epistatic effects in the population. The low resolution epistasis mapping method presented here identifies important epistatic interactions with a limited number of statistical tests at the cost of low precision.

LSD1-, EDS1- and PAD4-dependent conditional correlation among salicylic acid, hydrogen peroxide, water use efficiency and seed yield in Arabidopsis thaliana

In Arabidopsis thaliana, LESION SIMULATING DISEASE 1 (LSD1), ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) and PHYTOALEXIN DEFICIENT 4 (PAD4) proteins are regulators of cell death (CD) in response to abiotic and biotic stresses. Hormones, such as salicylic acid (SA), and reactive oxygen species, such as hydrogen peroxide (H₂ O₂ ), are key signaling molecules involved in plant CD. The proposed mathematical models presented in this study suggest that LSD1, EDS1 and PAD4 together with SA and H₂ O₂ are involved in the control of plant water use efficiency (WUE), vegetative growth and generative development. The analysis of Arabidopsis wild-type and single mutants lsd1, eds1, and pad4, as well as double mutants eds1/lsd1 and pad4/lsd1, demonstrated the strong conditional correlation between SA/H₂ O₂ and WUE that is dependent on LSD1, EDS1 and PAD4 proteins. Moreover, we found a strong correlation between the SA/H₂ O₂ homeostasis of 4-week-old Arabidopsis leaves and a total seed yield of 9-week-old plants. Altogether, our results prove that SA and H₂ O₂ are conditionally regulated by LSD1/EDS/PAD4 to govern WUE, biomass accumulation and seed yield. Conditional correlation and the proposed models presented in this study can be used as the starting points in the creation of a plant breeding algorithm that would allow to estimate the seed yield at the initial stage of plant growth, based on WUE, SA and H₂ O₂ content.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Study of QTLs linked to awn length and their relationships with chloroplasts under control and saline environments in bread wheat

INTRODUCTION

Some studies in wheat showed that awns may have a useful effect on yield, especially under drought stress. Up to this time few researches has identified the awn length QTLs with different effect in salinity stress.

OBJECTIVE

The primary objective of this study was to examine the additive (a) and the epistatic (aa) QTLs involve in wheat awns length in control and saline environments.

METHODS

A F7 RIL population consisting of 319 sister lines, derived from a cross between wheat cultivars Roshan and Falat (seri82), and the two parents were grown in two environments (control and Saline) based on an alpha lattice design with two replications in each environment. At flowering, awn length was measured for each line. For QTL analysis, the linkage map of the ''Roshan × Falat'' population was used, which included 748 markers including 719 DArT, 29 simple sequenced repeats (SSRs). Additive and pleiotropic QTLs were identified. In order to reveal the relationship between the identified QTL for awns length and the role of the gene or genes that it expresses, the awns length locus location and characteristics of its related CDS, gene, UTRs, ORF, exons and Introns were studied using ensemble plant ( http://plants.ensembl.org/Triticum_aestivum ). Furthermore, the promoter analysis has been done using NSITE-PL.

RESULTS

We identified 6 additive QTLs for awn length by QTL Cartographer program using single-environment phenotypical values. Also, we detected three additive and two epistatic QTLs for awn length by the QTLNetwork program using multi-environment phenotypical values. Our results showed that none of the additive and epistatic QTLs had interactions with environment. One of the additive QTLs located on chromosome 4A was co-located with QTLs for number of sterile spikelet per spike in both environment and number of seed per spike in control environment.

COCLUSION

Studies of the locus linked to the awns length QTL revealed the role of awn and its chloroplasts in grain filing during abiotic stress could be enhanced by over expression of some genes like GTP-Binding proteins which are enriched in chloroplasts encoded by genes included wPt-5730 locus.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Plants defend themselves against pathogens by activating immune responses. Unfortunately, these can cause unintended collateral damage to the plant itself. Nevertheless, some wild plants have genetic variants that confer a low threshold for the activation of immunity. While these enable a plant to respond particularly quickly to pathogen attack, such variants might be potentially dangerous. We are investigating one such variant of the immune gene ACCELERATED CELL DEATH 6 (ACD6) in the plant Arabidopsis thaliana. We discovered that there are variants at other genetic loci that can mask the effects of an overly active ACD6 gene. One of these genes, SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), codes for a known immune receptor. The SNC1 variant that attenuates ACD6 activity is rather common in A. thaliana populations, suggesting that new combinations of the hyperactive ACD6 variant and this antagonistic SNC1 variant will often arise by natural crosses. Similarly, because the two genes are unlinked, outcrossing will often lead to the hyperactive ACD6 variants being unmasked again. We propose that allelic diversity at SNC1 contributes to the maintenance of the hyperactive ACD6 variant in natural A. thaliana populations.

ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) affects development, photosynthesis, and hormonal homeostasis in hybrid aspen (Populus tremula L. × P. tremuloides)

ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) was first described as a protein involved in salicylic acid (SA)-, ethylene-, and reactive oxygen species (ROS)-dependent defense and acclimation responses. It is a molecular regulator of biotic and abiotic stress-induced programmed cell death. Its role is relatively well known in annual plants, such as Arabidopsis thaliana or Nicotiana benthamiana. However, little is known about its functions in woody plants. Therefore, in this study, we aimed to characterize the function of EDS1 in the Populus tremula L. × P. tremuloides hybrid grown for several seasons in the natural environment. We used two transgenic lines, eds1-7 and eds1-12, with decreased EDS1 expression levels in this study. The observed changes in physiological and biochemical parameters corresponded with the EDS1 silencing level. Both transgenic lines produced more lateral shoots in comparison to the wild-type (WT) plants, which resulted in the modification of tree morphology. Photosynthetic parameters, such as quantum yield of photosystem II (ϕPSII), photochemical and non-photochemical quenching (qP and NPQ, respectively), as well as chlorophyll content were found to be increased in both transgenic lines, which resulted in changes in photosynthetic efficiency. Our data also revealed lower foliar concentrations of SA and ROS, the latter resulting most probably from more efficient antioxidant system in both transgenic lines. In addition, our data indicated significantly decreased rate of leaf senescence during several autumn seasons. Transcriptomic analysis revealed deregulation of 2215 and 376 genes in eds1-12 and eds1-7, respectively, and also revealed 207 genes that were commonly deregulated in both transgenic lines. The deregulation was primarily observed in the genes involved in photosynthesis, signaling, hormonal metabolism, and development, which was found to agree with the results of biochemical and physiological tests. In general, our data proved that poplar EDS1 affects tree morphology, photosynthetic efficiency, ROS and SA metabolism, as well as leaf senescence.

Genetic background effects in quantitative genetics: gene-by-system interactions

Proper cell function depends on networks of proteins that interact physically and functionally to carry out physiological processes. Thus, it seems logical that the impact of sequence variation in one protein could be significantly influenced by genetic variants at other loci in a genome. Nonetheless, the importance of such genetic interactions, known as epistasis, in explaining phenotypic variation remains a matter of debate in genetics. Recent work from our lab revealed that genes implicated from an association study of toxin tolerance in Saccharomyces cerevisiae show extensive interactions with the genetic background: most implicated genes, regardless of allele, are important for toxin tolerance in only one of two tested strains. The prevalence of background effects in our study adds to other reports of widespread genetic-background interactions in model organisms. We suggest that these effects represent many-way interactions with myriad features of the cellular system that vary across classes of individuals. Such gene-by-system interactions may influence diverse traits and require new modeling approaches to accurately represent genotype-phenotype relationships across individuals.

FIT: statistical modeling tool for transcriptome dynamics under fluctuating field conditions

Motivation

Considerable attention has been given to the quantification of environmental effects on organisms. In natural conditions, environmental factors are continuously changing in a complex manner. To reveal the effects of such environmental variations on organisms, transcriptome data in field environments have been collected and analyzed. Nagano et al. proposed a model that describes the relationship between transcriptomic variation and environmental conditions and demonstrated the capability to predict transcriptome variation in rice plants. However, the computational cost of parameter optimization has prevented its wide application.

Results

We propose a new statistical model and efficient parameter optimization based on the previous study. We developed and released FIT, an R package that offers functions for parameter optimization and transcriptome prediction. The proposed method achieves comparable or better prediction performance within a shorter computational time than the previous method. The package will facilitate the study of the environmental effects on transcriptomic variation in field conditions.

Availability and Implementation

Freely available from CRAN (https://cran.r-project.org/web/packages/FIT/).

Supplementary information

Supplementary data are available at Bioinformatics online

Capturing pair-wise epistatic effects associated with three agronomic traits in barley

Genetic association mapping has been widely applied to determine genetic markers favorably associated with a trait of interest and provide information for marker-assisted selection. Many association mapping studies commonly focus on main effects due to intolerable computing intensity. This study aims to select several sets of DNA markers with potential epistasis to maximize genetic variations of some key agronomic traits in barley. By doing so, we integrated a MDR (multifactor dimensionality reduction) method with a forward variable selection approach. This integrated approach was used to determine single nucleotide polymorphism pairs with epistasis effects associated with three agronomic traits: heading date, plant height, and grain yield in barley from the barley Coordinated Agricultural Project. Our results showed that four, seven, and five SNP pairs accounted for 51.06, 45.66 and 40.42% for heading date, plant height, and grain yield, respectively with epistasis being considered, while corresponding contributions to these three traits were 45.32, 31.39, 31.31%, respectively without epistasis being included. The results suggested that epistasis model was more effective than non-epistasis model in this study and can be more preferred for other applications.

High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition

BACKGROUND

Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition.

RESULTS

Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one.

CONCLUSIONS

Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.

Polygenicity and Epistasis Underlie Fitness-Proximal Traits in the Caenorhabditis elegans Multiparental Experimental Evolution (CeMEE) Panel

Understanding the genetic basis of complex traits remains a major challenge in biology. Polygenicity, phenotypic plasticity, and epistasis contribute to phenotypic variance in ways that are rarely clear. This uncertainty can be problematic for estimating heritability, for predicting individual phenotypes from genomic data, and for parameterizing models of phenotypic evolution. Here, we report an advanced recombinant inbred line (RIL) quantitative trait locus mapping panel for the hermaphroditic nematode Caenorhabditis elegans, the C. elegans multiparental experimental evolution (CeMEE) panel. The CeMEE panel, comprising 507 RILs at present, was created by hybridization of 16 wild isolates, experimental evolution for 140-190 generations, and inbreeding by selfing for 13-16 generations. The panel contains 22% of single-nucleotide polymorphisms known to segregate in natural populations, and complements existing C. elegans mapping resources by providing fine resolution and high nucleotide diversity across > 95% of the genome. We apply it to study the genetic basis of two fitness components, fertility and hermaphrodite body size at time of reproduction, with high broad-sense heritability in the CeMEE. While simulations show that we should detect common alleles with additive effects as small as 5%, at gene-level resolution, the genetic architectures of these traits do not feature such alleles. We instead find that a significant fraction of trait variance, approaching 40% for fertility, can be explained by sign epistasis with main effects below the detection limit. In congruence, phenotype prediction from genomic similarity, while generally poor ([Formula: see text]), requires modeling epistasis for optimal accuracy, with most variance attributed to the rapidly evolving chromosome arms.

Mapping QTLs using a novel source of salinity tolerance from Hasawi and their interaction with environments in rice

BACKGROUND

Salinity is one of the most severe and widespread abiotic stresses that affect rice production. The identification of major-effect quantitative trait loci (QTLs) for traits related to salinity tolerance and understanding of QTL × environment interactions (QEIs) can help in more precise and faster development of salinity-tolerant rice varieties through marker-assisted breeding. Recombinant inbred lines (RILs) derived from IR29/Hasawi (a novel source of salinity) were screened for salinity tolerance in the IRRI phytotron in the Philippines (E1) and in two other diverse environments in Senegal (E2) and Tanzania (E3). QTLs were mapped for traits related to salinity tolerance at the seedling stage.

RESULTS

The RILs were genotyped using 194 polymorphic SNPs (single nucleotide polymorphisms). After removing segregation distortion markers (SDM), a total of 145 and 135 SNPs were used to construct a genetic linkage map with a length of 1655 and 1662 cM, with an average marker density of 11.4 cM in E1 and 12.3 cM in E2 and E3, respectively. A total of 34 QTLs were identified on 10 chromosomes for five traits using ICIM-ADD and segregation distortion locus (SDL) mapping (IM-ADD) under salinity stress across environments. Eight major genomic regions on chromosome 1 between 170 and 175 cM (qSES1.3, qSES1.4, qSL1.2, qSL1.3, qRL1.1, qRL1.2, qFWsht1.2, qDWsht1.2), chromosome 4 at 32 cM (qSES4.1, qFWsht4.2, qDWsht4.2), chromosome 6 at 115 cM (qFWsht6.1, qDWsht6.1), chromosome 8 at 105 cM (qFWsht8.1, qDWsht8.1), and chromosome 12 at 78 cM (qFWsht12.1, qDWsht12.1) have co-localized QTLs for the multiple traits that might be governing seedling stage salinity tolerance through multiple traits in different phenotyping environments, thus suggesting these as hot spots for tolerance of salinity. Forty-nine and 30 significant pair-wise epistatic interactions were detected between QTL-linked and QTL-unlinked regions using single-environment and multi-environment analyses.

CONCLUSIONS

The identification of genomic regions for salinity tolerance in the RILs showed that Hasawi possesses alleles that are novel for salinity tolerance. The common regions for the multiple QTLs across environments as co-localized regions on chromosomes 1, 4, 6, 8, and 12 could be due to linkage or pleiotropic effect, which might be helpful for multiple QTL introgression for marker-assisted breeding programs to improve the salinity tolerance of adaptive and popular but otherwise salinity-sensitive rice varieties.

Epistasis × environment interactions among Arabidopsis thaliana glucosinolate genes impact complex traits and fitness in the field

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

Main and epistatic loci studies in soybean for Sclerotinia sclerotiorum resistance reveal multiple modes of resistance in multi-environments

Genome-wide association (GWAS) and epistatic (GWES) studies along with expression studies in soybean [Glycine max (L.) Merr.] were leveraged to dissect the genetics of Sclerotinia stem rot (SSR) [caused by Sclerotinia sclerotiorum (Lib.) de Bary], a significant fungal disease causing yield and quality losses. A large association panel of 466 diverse plant introduction accessions were phenotyped in multiple field and controlled environments to: (1) discover sources of resistance, (2) identify SNPs associated with resistance, and (3) determine putative candidate genes to elucidate the mode of resistance. We report 58 significant main effect loci and 24 significant epistatic interactions associated with SSR resistance, with candidate genes involved in a wide range of processes including cell wall structure, hormone signaling, and sugar allocation related to plant immunity, revealing the complex nature of SSR resistance. Putative candidate genes [for example, PHYTOALEXIN DEFFICIENT 4 (PAD4), ETHYLENE-INSENSITIVE 3-LIKE 1 (EIL3), and ETHYLENE RESPONSE FACTOR 1 (ERF1)] clustered into salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) pathways suggest the involvement of a complex hormonal network typically activated by both necrotrophic (ET/JA) and biotrophic (SA) pathogens supporting that S. sclerotiorum is a hemibiotrophic plant pathogen.

Fitness effects of mutation: testing genetic redundancy in Arabidopsis thaliana

Screens of organisms with disruptive mutations in a single gene often fail to detect phenotypic consequences for the majority of mutants. One explanation for this phenomenon is that the presence of paralogous loci provides genetic redundancy. However, it is also possible that the assayed traits are affected by few loci, that effects could be subtle or that phenotypic effects are restricted to certain environments. We assayed a set of T-DNA insertion mutant lines of Arabidopsis thaliana to determine the frequency with which mutation affected fitness-related phenotypes. We found that between 8% and 42% of the assayed lines had altered fitness from the wild type. Furthermore, many of these lines exhibited fitness greater than the wild type. In a second experiment, we grew a subset of the lines in multiple environments and found whether a T-DNA insert increased or decreased fitness traits depended on the assay environment. Overall, our evidence contradicts the hypothesis that genetic redundancy is a common phenomenon in A. thaliana for fitness traits. Evidence for redundancy from prior screens of knockout mutants may often be an artefact of the design of the phenotypic assays which have focused on less complex phenotypes than fitness and have used single environments. Finally, our study adds to evidence that beneficial mutations may represent a significant component of the mutational spectrum of A. thaliana.

Promises and Challenges of Eco-Physiological Genomics in the Field: Tests of Drought Responses in Switchgrass

Identifying the physiological and genetic basis of stress tolerance in plants has proven to be critical to understanding adaptation in both agricultural and natural systems. However, many discoveries were initially made in the controlled conditions of greenhouses or laboratories, not in the field. To test the comparability of drought responses across field and greenhouse environments, we undertook three independent experiments using the switchgrass reference genotype Alamo AP13. We analyzed physiological and gene expression variation across four locations, two sampling times, and three years. Relatively similar physiological responses and expression coefficients of variation across experiments masked highly dissimilar gene expression responses to drought. Critically, a drought experiment utilizing small pots in the greenhouse elicited nearly identical physiological changes as an experiment conducted in the field, but an order of magnitude more differentially expressed genes. However, we were able to define a suite of several hundred genes that were differentially expressed across all experiments. This list was strongly enriched in photosynthesis, water status, and reactive oxygen species responsive genes. The strong across-experiment correlations between physiological plasticity-but not differential gene expression-highlight the complex and diverse genetic mechanisms that can produce phenotypically similar responses to various soil water deficits.

Analysis of recombination QTLs, segregation distortion, and epistasis for fitness in maize multiple populations using ultra-high-density markers

Using two nested association mapping populations and high-density markers, some important genomic regions controlling recombination frequency and segregation distortion were detected. Understanding the maize genomic features would be useful for the study of genetic diversity and evolution and for maize breeding. Here, we used two maize nested association mapping (NAM) populations separately derived in China (CN-NAM) and the US (US-NAM) to explore the maize genomic features. The two populations containing 36 families and about 7000 recombinant inbred lines were evaluated with genotyping-by-sequencing. Through the comparison between the two NAMs, we revealed that segregation distortion is little, whereas epistasis for fitness is present in the two maize NAM populations. When conducting quantitative trait loci (QTL) mapping for the total number of recombination events, we detected 14 QTLs controlling recombination. Using high-density markers to identify segregation distortion regions (SDRs), a total of 445 SDRs were detected within the 36 families, among which 15 common SDRs were found in at least ten families. About 80 % of the known maize gametophytic factors (ga) genes controlling segregation distortion were overlapped with highly significant SDRs. In addition, we also found that the regions with high recombination rate and high gene density usually tended to have little segregation distortion. This study will facilitate population genetic studies and gene cloning affecting recombination variation and segregation distortion in maize, which can improve plant breeding progress.

Genotypic variation in biomass allocation in response to field drought has a greater affect on yield than gas exchange or phenology

BACKGROUND

Plant performance in agricultural and natural settings varies with moisture availability, and understanding the range of potential drought responses and the underlying genetic architecture is important for understanding how plants will respond to both natural and artificial selection in various water regimes. Here, we raised genotypes of Brassica rapa under well-watered and drought treatments in the field. Our primary goal was to understand the genetic architecture and yield effects of different drought-escape and dehydration-avoidance strategies.

RESULTS

Drought treatments reduced soil moisture by 62 % of field capacity. Drought decreased biomass accumulation and fruit production by as much as 48 %, whereas instantaneous water-use efficiency and root:shoot ratio increased. Genotypes differed in the mean value of all traits and in the sensitivity of biomass accumulation, root:shoot ratio, and fruit production to drought. Bivariate correlations involving gas-exchange and phenology were largely constant across environments, whereas those involving root:shoot varied across treatments. Although root:shoot was typically unrelated to gas-exchange or yield under well-watered conditions, genotypes with low to moderate increases in root:shoot allocation in response to drought survived the growing season, maintained maximum photosynthesis levels, and produced more fruit than genotypes with the greatest root allocation under drought. QTL for gas-exchange and yield components (total biomass or fruit production) had common effects across environments while those for root:shoot were often environment-specific.

CONCLUSIONS

Increases in root allocation beyond those needed to survive and maintain favorable water relations came at the cost of fruit production. The environment-specific effects of root:shoot ratio on yield and the differential expression of QTL for this trait across water regimes have important implications for efforts to improve crops for drought resistance.

Naturally segregating loci exhibit epistasis for fitness

The extent to which gene interaction or epistasis contributes to fitness variation within populations remains poorly understood, despite its importance to a myriad of evolutionary questions. Here, we report a multi-year field study estimating fitness of Mimulus guttatus genetic lines in which pairs of naturally segregating loci exist in an otherwise uniform background. An allele at QTL x5b-a locus originally mapped for its effect on flower size-positively affects survival if combined with one genotype at quantitative trait locus x10a (aa) but has negative effects when combined with the other genotypes (Aa and AA). The viability differences between genotypes parallel phenotypic differences for the time and node at which a plant flowers. Viability is negatively correlated with fecundity across genotypes, indicating antagonistic pleiotropy for fitness components. This trade-off reduces the genetic variance for total fitness relative to the individual fitness components and thus may serve to maintain variation. Additionally, we find that the effects of each locus and their interaction often vary with the environment.

Identification of Quantitative Trait Loci and Water Environmental Interactions for Developmental Behaviors of Leaf Greenness in Wheat

The maintenance of leaf greenness in wheat, highly responsible for yield potential and resistance to drought stress, has been proved to be quantitatively inherited and susceptible to interact with environments by traditional genetic analysis. In order to further dissect the developmental genetic behaviors of flag leaf greenness under terminal drought, unconditional and conditional QTL mapping strategies were performed with a mixed linear model in 120 F8-derived recombinant inbred lines (RILs) from two Chinese common wheat cultivars (Longjian 19 × Q9086) in different water environments. A total of 65 additive QTLs (A-QTLs) and 42 pairs of epistatic QTLs (AA-QTLs) were identified as distribution on almost all 21 chromosomes except 5A, explaining from 0.24 to 3.29 % of the phenotypic variation. Of these, 22 A-QTLs and 25 pairs of AA-QTLs were common in two sets of mapping methods but the others differed. These putative QTLs were essentially characteristic of time- and environmentally-dependent expression patterns. Indeed some loci were expressed at two or more stages, while no single QTL was continually active through whole measuring duration. More loci were detected in early growth periods but most of QTL × water environment interactions (QEIs) happened in mid-anaphase, where drought stress was more conducted with negative regulation on QTL expressions. Compared to other genetic components, epistatic effects and additive QEIs effects could be predominant in regulating phenotypic variations during the ontogeny of leaf greenness. Several QTL cluster regions were suggestive of tight linkage or expression pleiotropy in the inheritance of these traits. Some reproducibly-expressed QTLs or common loci consistent with previously detected would be useful to the genetic improvement of staygreen types in wheat through MAS, especially in water-deficit environments.

Overexpression of GA20-OXIDASE1 impacts plant height, biomass allocation and saccharification efficiency in maize

Increased biomass yield and quality are of great importance for the improvement of feedstock for the biorefinery. For the production of bioethanol, both stem biomass yield and the conversion efficiency of the polysaccharides in the cell wall to fermentable sugars are of relevance. Increasing the endogenous levels of gibberellic acid (GA) by ectopic expression of GA20-OXIDASE1 (GA20-OX1), the rate-limiting step in GA biosynthesis, is known to affect cell division and cell expansion, resulting in larger plants and organs in several plant species. In this study, we examined biomass yield and quality traits of maize plants overexpressing GA20-OX1 (GA20-OX1). GA20-OX1 plants accumulated more vegetative biomass than control plants in greenhouse experiments, but not consistently over two years of field trials. The stems of these plants were longer but also more slender. Investigation of GA20-OX1 biomass quality using biochemical analyses showed the presence of more cellulose, lignin and cell wall residue. Cell wall analysis as well as expression analysis of lignin biosynthetic genes in developing stems revealed that cellulose and lignin were deposited earlier in development. Pretreatment of GA20-OX1 biomass with NaOH resulted in a higher saccharification efficiency per unit of dry weight, in agreement with the higher cellulose content. On the other hand, the cellulose-to-glucose conversion was slower upon HCl or hot-water pretreatment, presumably due to the higher lignin content. This study showed that biomass yield and quality traits can be interconnected, which is important for the development of future breeding strategies to improve lignocellulosic feedstock for bioethanol production.

Multiple abiotic stimuli are integrated in the regulation of rice gene expression under field conditions

Plants rely on transcriptional dynamics to respond to multiple climatic fluctuations and contexts in nature. We analyzed the genome-wide gene expression patterns of rice (Oryza sativa) growing in rainfed and irrigated fields during two distinct tropical seasons and determined simple linear models that relate transcriptomic variation to climatic fluctuations. These models combine multiple environmental parameters to account for patterns of expression in the field of co-expressed gene clusters. We examined the similarities of our environmental models between tropical and temperate field conditions, using previously published data. We found that field type and macroclimate had broad impacts on transcriptional responses to environmental fluctuations, especially for genes involved in photosynthesis and development. Nevertheless, variation in solar radiation and temperature at the timescale of hours had reproducible effects across environmental contexts. These results provide a basis for broad-based predictive modeling of plant gene expression in the field.

DOI:http://dx.doi.org/10.7554/eLife.08411.001

Plants need to be able to sense and respond to changes in temperature, light levels and other aspects of their environment. One way in which plants can rapidly respond to these changes is to modify how genes involved in growth and other processes are expressed. Therefore, understanding how this happens may help us to improve the ability of crops to grow when exposed to drought or other extreme environmental conditions.

Most previous studies into the effect of the environment on plant gene expression have been carried out under controlled conditions in a laboratory. These findings cannot reflect the full range of gene expression patterns that occur in the natural environment, where multiple factors (e.g. sunlight, water, nutrients) may vary at the same time. Therefore, it is important to also analyze the effect of fluctuations in multiple environmental factors in more complex field experiments.

Plessis et al. developed mathematical models to analyze the gene expression patterns of rice plants grown in the tropical environment of the Philippines using two different farming practices. One field of rice was flooded and constantly supplied with fresh water (referred to as the irrigated field), while the other field was dry and only received water from rainfall (the rainfed field). The experiments show that temperature and levels of sunlight (including UV radiation) have a strong impact on gene expression in the rice plants. Short-term variations in temperature and sunlight levels also have the most consistent effect across the different fields and seasons tested. However, for many genes, the plants grown in the irrigated field responded to the changes in environmental conditions in a different way to the plants grown in the rainfed field.

Further analysis identified groups of genes whose expression combined responses to several environmental factors at the same time. For example, certain genes that responded to increases in sunlight in the absence of drought responded to both sunlight levels and the shortage of water when a drought occurred. The next step is to test more types of environments and climates to be able to predict gene expression responses under future climatic conditions.

DOI:http://dx.doi.org/10.7554/eLife.08411.002

Natural genetic variation in Arabidopsis thaliana defense metabolism genes modulates field fitness

Natural populations persist in complex environments, where biotic stressors, such as pathogen and insect communities, fluctuate temporally and spatially. These shifting biotic pressures generate heterogeneous selective forces that can maintain standing natural variation within a species. To directly test if genes containing causal variation for the Arabidopsis thaliana defensive compounds, glucosinolates (GSL) control field fitness and are therefore subject to natural selection, we conducted a multi-year field trial using lines that vary in only specific causal genes. Interestingly, we found that variation in these naturally polymorphic GSL genes affected fitness in each of our environments but the pattern fluctuated such that highly fit genotypes in one trial displayed lower fitness in another and that no GSL genotype or genotypes consistently out-performed the others. This was true both across locations and within the same location across years. These results indicate that environmental heterogeneity may contribute to the maintenance of GSL variation observed within Arabidopsis thaliana.

DOI:http://dx.doi.org/10.7554/eLife.05604.001

‘Genetic variation’ describes the naturally occurring differences in DNA sequences that are found among individuals of the same species. These genetic differences arise from random mutations and may be passed on to their offspring. Some of these mutations may improve the ability of an individual to survive and reproduce—known as fitness—and are likely to become more common in the population. Other mutations may reduce an individual's fitness and are likely to be lost. However, it is believed that most of the mutations will have no effect on the fitness of individuals.

It is not known why many of these ‘neutral’ genetic differences are maintained in populations. Some researchers have proposed that they are kept by chance and that there is no direct advantage to the population of keeping them unless these neutral mutations later become beneficial. However, other researchers think that the genetic variation itself may improve the fitness of the population by allowing it to quickly adapt to changes in the environment.

Arabidopsis thaliana is a small plant that lives in many different environments and has high levels of genetic variation in many of its physical traits. One of these traits is the production of molecules called glucosinolates, which help the plants to defend against herbivores and infection by microbes. Previous studies have suggested that variation in the genes that make glucosinolates may improve the fitness of A. thaliana populations.

To test this idea, Kerwin et al. carried out a field trial using A. thaliana plants that were genetically identical except for some of the genes involved in the production of glucosinolates. Kerwin et al. grew the plants in several different environments over several years. The field trial shows that variation in these genes affected the fitness of the plants in each of the different environments. However, the fitness benefit depended on the environment, and no single gene variant provided the best fitness across all environments, or over all the years of the trial.

Kerwin et al.'s findings suggest that changes in the environment may contribute to the maintenance of genetic variation in the genes that make glucosinolates. This raises the questions of how many other genes in plants (or other species such as humans) have genetic variation that contributes to fitness across varied environments; and how can this link be tested in natural settings.

DOI:http://dx.doi.org/10.7554/eLife.05604.002

What is crop heterosis: new insights into an old topic

Heterosis (or hybrid vigor) refers to a natural phenomenon whereby hybrid offspring of genetically diverse individuals out-perform their parents in multiple traits including yield, adaptability and resistances to biotic and abiotic stressors. Innovations in technology and research continue to clarify the mechanisms underlying crop heterosis, however the intrinsic relationship between the biological basis of heterosis remain unclear. In this review, we aim to provide insight into the molecular genetic basis of heterosis by presenting recent advances in the 'omics' of heterosis and the role of non-coding regions, particularly in relation to energy-use efficiency. We propose that future research should focus on integrating the expanding datasets from different species and hybrid combinations, to mine key heterotic genes and unravel interactive 'omics' networks associated with heterosis. Improved understanding of heterosis and the biological basis for its manipulation in agriculture should help to streamline its use in enhancing crop productivity.

Endopolyploidy as a potential alternative adaptive strategy for Arabidopsis leaf size variation in response to UV-B

The extent of endoreduplication in leaf growth is group- or even species-specific, and its adaptive role is still unclear. A survey of Arabidopsis accessions for variation at the level of endopolyploidy, cell number, and cell size in leaves revealed extensive genetic variation in endopolyploidy level. High endopolyploidy is associated with increased leaf size, both in natural and in genetically unstructured (mapping) populations. The underlying genes were identified as quantitative trait loci that control endopolyploidy in nature by modulating the progression of successive endocycles during organ development. This complex genetic architecture indicates an adaptive mechanism that allows differential organ growth over a broad geographic range and under stressful environmental conditions. UV-B radiation was identified as a significant positive climatic predictor for high endopolyploidy. Arabidopsis accessions carrying the increasing alleles for endopolyploidy also have enhanced tolerance to UV-B radiation. UV-absorbing secondary metabolites provide an additional protective strategy in accessions that display low endopolyploidy. Taken together, these results demonstrate that high constitutive endopolyploidy is a significant predictor for organ size in natural populations and is likely to contribute to sustaining plant growth under high incident UV radiation. Endopolyploidy may therefore form part of the range of UV-B tolerance mechanisms that exist in natural populations.

The genetic architecture of constitutive and induced trichome density in two new recombinant inbred line populations of Arabidopsis thaliana: phenotypic plasticity, epistasis, and bidirectional leaf damage response

BACKGROUND

Herbivory imposes an important selective pressure on plants. In Arabidopsis thaliana leaf trichomes provide a key defense against insect herbivory; however, trichome production incurs a fitness cost in the absence of herbivory. Previous work on A. thaliana has shown an increase in trichome density in response to leaf damage, suggesting a mechanism by which the cost associated with constitutively high trichome density might be mitigated; however, the genetic basis of trichome density induction has not been studied.

RESULTS

Here, we describe the mapping of quantitative trait loci (QTL) for constitutive and damage induced trichome density in two new recombinant inbred line populations of A. thaliana; mapping for constitutive and induced trichome density also allowed for the investigation of damage response (plasticity) QTL. Both novel and previously identified QTL for constitutive trichome density and the first QTL for induced trichome density and response are identified. Interestingly, two of the four parental accessions and multiple RILs in each population exhibited lower trichome density following leaf damage, a response not previously described in A. thaliana. Importantly, a single QTL was mapped for the response phenotype and allelic variation at this locus appears to determine response trajectory in RILs. The data also show that epistatic interactions are a significant component of the genetic architecture of trichome density.

CONCLUSIONS

Together, our results provide further insights into the genetic architecture of constitutive trichome density and new insights into induced trichome density in A. thaliana specifically and to our understanding of the genetic underpinnings of natural variation generally.

Combining growth-promoting genes leads to positive epistasis in Arabidopsis thaliana

Several genes positively influence final leaf size in Arabidopsis when mutated or overexpressed. The connections between these growth regulators are still poorly understood although such knowledge would further contribute to understand the processes driving leaf growth. In this study, we performed a combinatorial screen with 13 transgenic Arabidopsis lines with an increased leaf size. We found that from 61 analyzed combinations, 39% showed an additional increase in leaf size and most resulted from a positive epistasis on growth. Similar to what is found in other organisms in which such an epistasis assay was performed, only few genes were highly connected in synergistic combinations as we observed a positive epistasis in the majority of the combinations with samba, BRI1(OE) or SAUR19(OE). Furthermore, positive epistasis was found with combinations of genes with a similar mode of action, but also with genes which affect distinct processes, such as cell proliferation and cell expansion.DOI: http://dx.doi.org/10.7554/eLife.02252.001.

QTL architecture of reproductive fitness characters in Brassica rapa

BACKGROUND

Reproductive output is critical to both agronomists seeking to increase seed yield and to evolutionary biologists interested in understanding natural selection. We examine the genetic architecture of diverse reproductive fitness traits in recombinant inbred lines (RILs) developed from a crop (seed oil) × wild-like (rapid cycling) genotype of Brassica rapa in field and greenhouse environments.

RESULTS

Several fitness traits showed strong correlations and QTL-colocalization across environments (days to bolting, fruit length and seed color). Total fruit number was uncorrelated across environments and most QTL affecting this trait were correspondingly environment-specific. Most fitness components were positively correlated, consistent with life-history theory that genotypic variation in resource acquisition masks tradeoffs. Finally, we detected evidence of transgenerational pleiotropy, that is, maternal days to bolting was negatively correlated with days to offspring germination. A QTL for this transgenerational correlation was mapped to a genomic region harboring one copy of FLOWERING LOCUS C, a genetic locus known to affect both days to flowering as well as germination phenotypes.

CONCLUSIONS

This study characterizes the genetic structure of important fitness/yield traits within and between generations in B. rapa. Several identified QTL are suitable candidates for fine-mapping for the improvement of yield in crop Brassicas. Specifically, brFLC1, warrants further investigation as a potential regulator of phenology between generations.

Identification of QTLs for seed and pod traits in soybean and analysis for additive effects and epistatic effects of QTLs among multiple environments

Soybean seed and pod traits are important yield components. Selection for high yield style in seed and pod along with agronomic traits is a goal of many soybean breeders. The intention of this study was to identify quantitative trait loci (QTL) underlying seed and pod traits in soybean among eleven environments in China. 147 recombinant inbred lines were advanced through single-seed-descent method. The population was derived from a cross between Charleston (an American high yield soybean cultivar) and DongNong594 (a Chinese high yield soybean cultivar). A total of 157 polymorphic simple sequence repeat markers were used to construct a genetic linkage map. The phenotypic data of seed and pod traits [number of one-seed pod, number of two-seed pod, number of three-seed pod, number of four-seed pod, number of (two plus three)-seed pod, number of (three plus four)-seed pod, seed weight per plant, number of pod per plant] were recorded in eleven environments. In the analysis of single environment, fourteen main effect QTLs were identified. In the conjoint analysis of multiple environments, twenty-four additive QTLs were identified, and additive QTLs by environments interactions (AE) were evaluated and analyzed at the same time among eleven environments; twenty-three pairs of epistatic QTLs were identified, and epistasis (additive by additive) by environments interactions (AAE) were also analyzed and evaluated among eleven environments. Comparing the results of identification between single environment mapping and multiple environments conjoint mapping, three main effect QTLs with positive additive values and another three main effect QTLs with negative additive values, had no interactions with all environments, supported that these QTLs could be used in molecular assistant breeding in the future. These different effect QTLs could supply a good foundation to the gene clone and molecular asisstant breeding of soybean seed and pod traits.

Genetic mapping of three quantitative trait loci for soybean aphid resistance in PI 567324

Host-plant resistance is an effective method for controlling soybean aphid (Aphis glycines Matsumura), the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America. Recently, resistant soybean lines have been discovered and at least four aphid resistance genes (Rag1, Rag2, Rag3 and rag4) have been mapped on different soybean chromosomes. However, the evolution of new soybean aphid biotypes capable of defeating host-plant resistance conferred by most single genes demonstrates the need for finding germplasm with multigenic resistance to the aphid. This study was conducted to map quantitative trait loci (QTL) for aphid resistance in PI 567324. We identified two major QTL (QTL_13_1 and QTL_13_2) for aphid resistance on soybean chromosome 13 using 184 recombinant inbred lines from a 'Wyandot' × PI 567324 cross. QTL_13_1 was located close to the previously reported Rag2 gene locus, and QTL_13_2 was close to the rag4 locus. A minor QTL (QTL_6_1) was also detected on chromosome 6, where no gene for soybean aphid resistance has been reported so far. These results indicate that PI 567324 possesses oligogenic resistance to the soybean aphid. The molecular markers closely linked to the QTL reported here will be useful for development of cultivars with oligogenic resistance that are expected to provide broader and more durable resistance against soybean aphids compared with cultivars with monogenic resistance.

Multivariable environmental conditions promote photosynthetic adaptation potential in Arabidopsis thaliana

Most of our knowledge on the regulation of photosynthesis originates from studies performed in highly controlled laboratory conditions. However, in their natural habitats plants are simultaneously subjected to a broad range of abiotic and biotic stimuli which influence photosynthetic efficiency; hence there is an emerging need to examine the process of photosynthesis under multivariable field conditions in order to elucidate the mechanisms that underlie its regulation. Such knowledge has potential for providing novel targets that would improve both crop yield and performance. In the present study we compared laboratory- and field-grown Arabidopsis thaliana plants in terms of photosynthetic efficiency in modulated light intensities and CO2 concentrations. We show here that the field-acclimated plants display highly efficient photosynthesis and are more tolerant to variable light intensities and CO2 concentrations than their laboratory-grown counterparts. We also demonstrate that some structural rearrangements of LHCII and PSII, together with altered pigments composition and stomatal density, are responsible for the differences in assimilation and photochemistry. Furthermore, we employ a transcript profiling approach to explain the genetic mechanisms underlying these adaptations and suggest that they are mainly induced by the high and fluctuating light intensities which occur in the natural environment.

Genetic dissection of growth traits in a Chinese indigenous × commercial broiler chicken cross

Background

In China, consumers often prefer indigenous broiler chickens over commercial breeds, as they have characteristic meat qualities requested within traditional culinary customs. However, the growth-rate of these indigenous breeds is slower than that of the commercial broilers, which means they have not yet reached their full economic value. Therefore, combining the valuable meat quality of the native chickens with the efficiency of the commercial broilers is of interest. In this study, we generated an F₂ intercross between the slow growing native broiler breed, Huiyang Beard chicken, and the fast growing commercial broiler breed, High Quality chicken Line A, and used it to map loci explaining the difference in growth rate between these breeds.

Results

A genome scan to identify main-effect loci affecting 24 growth-related traits revealed nine distinct QTL on six chromosomes. Many QTL were pleiotropic and conformed to the correlation patterns observed between phenotypes. Most of the mapped QTL were found in locations where growth QTL have been reported in other populations, although the effects were greater in this population. A genome scan for pairs of interacting loci identified a number of additional QTL in 10 other genomic regions. The epistatic pairs explained 6–8% of the residual phenotypic variance. Seven of the 10 epistatic QTL mapped in regions containing candidate genes in the ubiquitin mediated proteolysis pathway, suggesting the importance of this pathway in the regulation of growth in this chicken population.

Conclusions

The main-effect QTL detected using a standard one-dimensional genome scan accounted for a significant fraction of the observed phenotypic variance in this population. Furthermore, genes in known pathways present interesting candidates for further exploration. This study has thus located several QTL regions as promising candidates for further study, which will increase our understanding of the genetic mechanisms underlying growth-related traits in chickens.

Among- and within-population variation in flowering time of Iberian Arabidopsis thaliana estimated in field and glasshouse conditions

The study of the evolutionary and population genetics of quantitative traits requires the assessment of within- and among-population patterns of variation. We carried out experiments including eight Iberian Arabidopsis thaliana populations (10 individuals per population) in glasshouse and field conditions. We quantified among- and within-population variation for flowering time and for several field life-history traits. Individuals were genotyped with microsatellites, single nucleotide polymorphisms and four well-known flowering genes (FRI, FLC, CRY2 and PHYC). Phenotypic and genotypic data were used to conduct Q(ST)-F(ST) comparisons. Life-history traits varied significantly among- and within-populations. Flowering time also showed substantial within- and among-population variation as well as significant genotype × environment interactions among the various conditions. Individuals bearing FRI truncations exhibited reduced recruitment in field conditions and differential flowering time behavior across experimental conditions, suggesting that FRI contributes to the observed significant genotype × environment interactions. Flowering time estimated in field conditions was the only trait showing significantly higher quantitative genetic differentiation than neutral genetic differentiation values. Overall, our results show that these A. thaliana populations are genetically more differentiated for flowering time than for neutral markers, suggesting that flowering time is likely to be under divergent selection.