Selection on QTL and complex traits in complex environments

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Genetic Characterization of Multiple Components Contributing to Fusarium Head Blight Resistance of FL62R1, a Canadian Bread Wheat Developed Using Systemic Breeding

Fusarium head blight (FHB) is a devastating fungal disease of small-grain cereals that results in severe yield and quality losses. FHB resistance is controlled by resistance components including incidence, field severity, visual rating index, Fusarium damaged kernels (FDKs), and the accumulation of the mycotoxin deoxynivalenol (DON). Resistance conferred by each of these components is partial and must be combined to achieve resistance sufficient to protect wheat from yield losses. In this study, two biparental mapping populations were analyzed in Canadian FHB nurseries and quantitative trait loci (QTL) mapped for the traits listed above. Nine genomic loci, on 2AS, 2BS, 3BS, 4AS, 4AL, 4BS, 5AS, 5AL, and 5BL, were enriched for the majority of the QTL controlling FHB resistance. The previously validated FHB resistance QTL on 3BS and 5AS affected resistance to severity, FDK, and DON in these populations. The remaining seven genomic loci colocalize with flowering time and/or plant height QTL. The QTL on 4B was a major contributor to all field resistance traits and plant height in the field. QTL on 4AL showed contrasting effects for FHB resistance between Eastern and Western Canada, indicating a local adapted resistance to FHB. In addition, we also found that the 2AS QTL contributed a major effect for DON, and the 2BS for FDK, while the 5AL conferred mainly effect for both FDK/DON. Results presented here provide insight into the genetic architecture underlying these resistant components and insight into how FHB resistance in wheat is controlled by a complex network of interactions between genes controlling flowering time, plant height, local adaption, and FHB resistance components.

Major Novel QTL for Resistance to Cassava Bacterial Blight Identified through a Multi-Environmental Analysis

Cassava, Manihot esculenta Crantz, has been positioned as one of the most promising crops world-wide representing the staple security for more than one billion people mainly in poor countries. Cassava production is constantly threatened by several diseases, including cassava bacterial blight (CBB) caused by Xanthomonas axonopodis pv. manihotis (Xam), it is the most destructive disease causing heavy yield losses. Here, we report the detection and localization on the genetic map of cassava QTL (Quantitative Trait Loci) conferring resistance to CBB. An F1 mapping population of 117 full sibs was tested for resistance to two Xam strains (Xam318 and Xam681) at two locations in Colombia: La Vega, Cundinamarca and Arauca. The evaluation was conducted in rainy and dry seasons and additional tests were carried out under controlled greenhouse conditions. The phenotypic evaluation of the response to Xam revealed continuous variation. Based on composite interval mapping analysis, 5 strain-specific QTL for resistance to Xam explaining between 15.8 and 22.1% of phenotypic variance, were detected and localized on a high resolution SNP-based genetic map of cassava. Four of them show stability among the two evaluated seasons. Genotype by environment analysis detected three QTL by environment interactions and the broad sense heritability for Xam318 and Xam681 were 20 and 53%, respectively. DNA sequence analysis of the QTL intervals revealed 29 candidate defense-related genes (CDRGs), and two of them contain domains related to plant immunity proteins, such as NB-ARC-LRR and WRKY.

Genome-wide association mapping of time-dependent growth responses to moderate drought stress in Arabidopsis

Large areas of arable land are often confronted with irregular rainfall resulting in limited water availability for part(s) of the growing seasons, which demands research for drought tolerance of plants. Natural variation was observed for biomass accumulation upon controlled moderate drought stress in 324 natural accessions of Arabidopsis. Improved performance under drought stress was correlated with early flowering and lack of vernalization requirement, indicating overlap in the regulatory networks of flowering time and drought response or correlated responses of these traits to natural selection. In addition, plant size was negatively correlated with relative water content (RWC) independent of the absolute water content (WC), indicating a prominent role for soluble compounds. Growth in control and drought conditions was determined over time and was modelled by an exponential function. Genome-wide association (GWA) mapping of temporal plant size data and of model parameters resulted in the detection of six time-dependent quantitative trait loci (QTLs) strongly associated with drought. Most QTLs would not have been identified if plant size was determined at a single time point. Analysis of earlier reported gene expression changes upon drought enabled us to identify for each QTL the most likely candidates.