Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Identification of novel genes associated with herbicide tolerance in Lentil (Lens culinaris ssp. culinaris Medik.)

Weeds pose a major constraint in lentil cultivation, leading to decrease farmers' revenues by reducing the yield and increasing the management costs. The development of herbicide tolerant cultivars is essential to increase lentil yield. Even though herbicide tolerant lines have been identified in lentils, breeding efforts are still limited and lack proper validation. Marker assisted selection (MAS) can increase selection accuracy at early generations. Total 292 lentil accessions were evaluated under different dosages of two herbicides, metribuzin and imazethapyr, during two seasons at Marchouch, Morocco and Terbol, Lebanon. Highly significant differences among accessions were observed for days to flowering (DF) and maturity (DM), plant height (PH), biological yield (BY), seed yield (SY), number of pods per plant (NP), as well as the reduction indices (RI) for PH, BY, SY and NP. A total of 10,271 SNPs markers uniformly distributed along the lentil genome were assayed using Multispecies Pulse SNP chip developed at Agriculture Victoria, Melbourne. Meta-GWAS analysis was used to detect marker-trait associations, which detected 125 SNPs markers associated with different traits and clustered in 85 unique quantitative trait loci. These findings provide valuable insights for initiating MAS programs aiming to enhance herbicide tolerance in lentil crop.

Integration of genome-wide association studies, metabolomics, and transcriptomics reveals phenolic acid- and flavonoid-associated genes and their regulatory elements under drought stress in rapeseed flowers

Introduction

Biochemical and metabolic processes help plants tolerate the adverse effects of drought. In plants accumulating bioactive compounds, understanding the genetic control of the biosynthesis of biochemical pathways helps the discovery of candidate gene (CG)-metabolite relationships.

Methods

The metabolic profile of flowers in 119 rapeseed (Brassica napus) accessions was assessed over two irrigation treatments, one a well-watered (WW) condition and the other a drought stress (DS) regime. We integrated information gained from 52,157 single-nucleotide polymorphism (SNP) markers, metabolites, and transcriptomes to identify linked SNPs and CGs responsible for the genetic control of flower phenolic compounds and regulatory elements.

Results

In a genome-wide association study (GWAS), of the SNPs tested, 29,310 SNPs were qualified to assess the population structure and linkage disequilibrium (LD), of which several SNPs for radical scavenging activity (RSA) and total flavanol content (TFLC) were common between the two irrigation conditions and pleiotropic SNPs were found for chlorogenic and coumaric acids content. The principal component analysis (PCA) and stepwise regression showed that chlorogenic acid and epicatechin in WW and myricetin in DS conditions were the most important components for RSA. The hierarchical cluster analysis (HCA) showed that vanillic acid, myricetin, gallic acid, and catechin were closely associated in both irrigation conditions. Analysis of GWAS showed that 60 CGs were identified, of which 18 were involved in stress-induced pathways, phenylpropanoid pathway, and flavonoid modifications. Of the CGs, PAL1, CHI, UGT89B1, FLS3, CCR1, and CYP75B137 contributed to flavonoid biosynthetic pathways. The results of RNA sequencing (RNA-seq) revealed that the transcript levels of PAL, CHI, and CYP75B137 known as early flavonoid biosynthesis-related genes and FLS3, CCR1, and UGT89B1 related to the later stages were increased during drought conditions. The transcription factors (TFs) NAC035 and ERF119 related to flavonoids and phenolic acids were upregulated under drought conditions.

Discussion

These findings expand our knowledge on the response mechanisms to DS, particularly regarding the regulation of key phenolic biosynthetic genes in rapeseed. Our data also provided specific linked SNPs for marker-assisted selection (MAS) programs and CGs as resources toward realizing metabolomics-associated breeding of rapeseed.

Identification of candidate genes for LepR1 resistance against Leptosphaeria maculans in Brassica napus

Utilising resistance (R) genes, such as LepR1, against Leptosphaeria maculans, the causal agent of blackleg in canola (Brassica napus), could help manage the disease in the field and increase crop yield. Here we present a genome wide association study (GWAS) in B. napus to identify LepR1 candidate genes. Disease phenotyping of 104 B. napus genotypes revealed 30 resistant and 74 susceptible lines. Whole genome re-sequencing of these cultivars yielded over 3 million high quality single nucleotide polymorphisms (SNPs). GWAS in mixed linear model (MLM) revealed a total of 2,166 significant SNPs associated with LepR1 resistance. Of these SNPs, 2108 (97%) were found on chromosome A02 of B. napus cv. Darmor bzh v9 with a delineated LepR1_mlm1 QTL at 15.11-26.08 Mb. In LepR1_mlm1, there are 30 resistance gene analogs (RGAs) (13 nucleotide-binding site-leucine rich repeats (NLRs), 12 receptor-like kinases (RLKs), and 5 transmembrane-coiled-coil (TM-CCs)). Sequence analysis of alleles in resistant and susceptible lines was undertaken to identify candidate genes. This research provides insights into blackleg resistance in B. napus and assists identification of the functional LepR1 blackleg resistance gene.

Genome-wide association studies dissect the G × E interaction for agronomic traits in a worldwide collection of safflowers (Carthamus tinctorius L.)

Genome-wide association studies were conducted using a globally diverse safflower (Carthamus tinctorius L.) Genebank collection for grain yield (YP), days to flowering (DF), plant height (PH), 500 seed weight (SW), seed oil content (OL), and crude protein content (PR) in four environments (sites) that differed in water availability. Phenotypic variation was observed for all traits. YP exhibited low overall genetic correlations (rGoverall) across sites, while SW and OL had high rGoverall and high pairwise genetic correlations (rGij) across all pairwise sites. In total, 92 marker-trait associations (MTAs) were identified using three methods, single locus genome-wide association studies (GWAS) using a mixed linear model (MLM), the Bayesian multi-locus method (BayesR), and meta-GWAS. MTAs with large effects across all sites were detected for OL, SW, and PR, and MTAs specific for the different water stress sites were identified for all traits. Five MTAs were associated with multiple traits; 4 of 5 MTAs were variously associated with the three traits of SW, OL, and PR. This study provided insights into the phenotypic variability and genetic architecture of important safflower agronomic traits under different environments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01295-8.

Data Integration, Imputation, and Meta-analysis for Genome-Wide Association Studies

Growing genomic and phenotypic datasets require different groups around the world to collaborate and integrate these valuable resources to maximize their benefit and increase reference population sizes for genomic prediction and genome-wide association studies (GWAS). However, different studies use different genotyping techniques which requires a synchronizing step for the genotyped variants called "imputation" before combining them. Optimally, different GWAS datasets can be analysed within a meta-analysis, which recruits summary statistics instead of actual data. This chapter describes the general principles for genotypic imputation and meta-GWAS analysis with a description of study designs and command lines required for such analyses.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization

Array-based single nucleotide polymorphism (SNP) genotyping platforms have low genotype error and missing data rates compared to genotyping-by-sequencing technologies. However, design decisions used to create array-based SNP genotyping assays for both research and breeding applications are critical to their success. We describe a novel approach applicable to any animal or plant species for the design of cost-effective imputation-enabled SNP genotyping arrays with broad utility and demonstrate its application through the development of the Illumina Infinium Wheat Barley 40K SNP array Version 1.0. We show that the approach delivers high quality and high resolution data for wheat and barley, including when samples are jointly hybridised. The new array aims to maximally capture haplotypic diversity in globally diverse wheat and barley germplasm while minimizing ascertainment bias. Comprising mostly biallelic markers that were designed to be species-specific and single-copy, the array permits highly accurate imputation in diverse germplasm to improve the statistical power of genome-wide association studies (GWAS) and genomic selection. The SNP content captures tetraploid wheat (A- and B-genome) and Aegilops tauschii Coss. (D-genome) diversity and delineates synthetic and tetraploid wheat from other wheat, as well as tetraploid species and subgroups. The content includes SNP tagging key trait loci in wheat and barley, as well as direct connections to other genotyping platforms and legacy datasets. The utility of the array is enhanced through the web-based tool, Pretzel (https://plantinformatics.io/) which enables the content of the array to be visualized and interrogated interactively in the context of numerous genetic and genomic resources to be connected more seamlessly to research and breeding. The array is available for use by the international wheat and barley community.

Dissection of Quantitative Blackleg Resistance Reveals Novel Variants of Resistance Gene Rlm9 in Elite Brassica napus

Blackleg is one of the major fungal diseases in oilseed rape/canola worldwide. Most commercial cultivars carry R gene-mediated qualitative resistances that confer a high level of race-specific protection against Leptosphaeria maculans, the causal fungus of blackleg disease. However, monogenic resistances of this kind can potentially be rapidly overcome by mutations in the pathogen's avirulence genes. To counteract pathogen adaptation in this evolutionary arms race, there is a tremendous demand for quantitative background resistance to enhance durability and efficacy of blackleg resistance in oilseed rape. In this study, we characterized genomic regions contributing to quantitative L. maculans resistance by genome-wide association studies in a multiparental mapping population derived from six parental elite varieties exhibiting quantitative resistance, which were all crossed to one common susceptible parental elite variety. Resistance was screened using a fungal isolate with no corresponding avirulence (AvrLm) to major R genes present in the parents of the mapping population. Genome-wide association studies revealed eight significantly associated quantitative trait loci (QTL) on chromosomes A07 and A09, with small effects explaining 3-6% of the phenotypic variance. Unexpectedly, the qualitative blackleg resistance gene Rlm9 was found to be located within a resistance-associated haploblock on chromosome A07. Furthermore, long-range sequence data spanning this haploblock revealed high levels of single-nucleotide and structural variants within the Rlm9 coding sequence among the parents of the mapping population. The results suggest that novel variants of Rlm9 could play a previously unknown role in expression of quantitative disease resistance in oilseed rape.

Genome-wide association study and transcriptome analysis dissect the genetic control of silique length in Brassica napus L

BACKGROUND

Rapeseed is the third-largest oilseed crop after soybeans and palm that produces vegetable oil for human consumption and biofuel for industrial production. Silique length (SL) is an important trait that is strongly related to seed yield in rapeseed. Although many studies related to SL have been reported in rapeseed, only a few candidate genes have been found and cloned, and the genetic mechanisms regulating SL in rapeseed remain unclear. Here, we dissected the genetic basis of SL by genome-wide association studies (GWAS) combined with transcriptome analysis.

RESULTS

We identified quantitative trait locus (QTL) for SL using a recombinant inbred line (RIL) population and two independent GWAS populations. Major QTLs on chromosomes A07, A09, and C08 were stably detected in all environments from all populations. Several candidate genes related to starch and sucrose metabolism, plant hormone signal transmission and phenylpropanoid biosynthesis were detected in the main QTL intervals, such as BnaA9.CP12-2, BnaA9.NST2, BnaA7.MYB63, and BnaA7.ARF17. In addition, the results of RNA-seq and weighted gene co-expression network analysis (WGCNA) showed that starch and sucrose metabolism, photosynthesis, and secondary cell wall biosynthesis play an important role in the development of siliques.

CONCLUSIONS

We propose that photosynthesis, sucrose and starch metabolism, plant hormones, and lignin content play important roles in the development of rapeseed siliques.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

QTL Mapping and GWAS Reveal the Genetic Mechanism Controlling Soluble Solids Content in Brassica napus Shoots

Oilseed-vegetable-dual-purpose (OVDP) rapeseed can effectively alleviate the land contradiction between crops and it supplements vegetable supplies in winter or spring. The soluble solids content (SSC) is an important index that is used to evaluate the quality and sugar content of fruits and vegetables. However, the genetic architecture underlying the SSC in Brassica napus shoots is still unclear. Here, quantitative trait loci (QTLs) for the SSC in B. napus shoots were investigated by performing linkage mapping using a recombinant inbred line population containing 189 lines. A germplasm set comprising 302 accessions was also used to conduct a genome-wide association study (GWAS). The QTL mapping revealed six QTLs located on chromosomes A01, A04, A08, and A09 in two experiments. Among them, two major QTLs, qSSC/21GY.A04-1 and qSSC/21NJ.A08-1, accounted for 12.92% and 10.18% of the phenotypic variance, respectively. In addition, eight single-nucleotide polymorphisms with phenotypic variances between 5.62% and 10.18% were identified by the GWAS method. However, no locus was simultaneously identified by QTL mapping and GWAS. We identified AH174 (7.55 °Brix and 7.9 °Brix), L166 (8.9 °Brix and 8.38 °Brix), and L380 (8.9 °Brix and 7.74 °Brix) accessions can be used as superior parents. These results provide valuable information that increases our understanding of the genetic control of SSC and will facilitate the breeding of high-SSC B. napus shoots.

Status and advances in mining for blackleg (Leptosphaeria maculans) quantitative resistance (QR) in oilseed rape (Brassica napus)

KEY MESSAGE

Quantitative resistance (QR) loci discovered through genetic and genomic analyses are abundant in the Brassica napus genome, providing an opportunity for their utilization in enhancing blackleg resistance. Quantitative resistance (QR) has long been utilized to manage blackleg in Brassica napus (canola, oilseed rape), even before major resistance genes (R-genes) were extensively explored in breeding programmes. In contrast to R-gene-mediated qualitative resistance, QR reduces blackleg symptoms rather than completely eliminating the disease. As a polygenic trait, QR is controlled by numerous genes with modest effects, which exerts less pressure on the pathogen to evolve; hence, its effectiveness is more durable compared to R-gene-mediated resistance. Furthermore, combining QR with major R-genes has been shown to enhance resistance against diseases in important crops, including oilseed rape. For these reasons, there has been a renewed interest among breeders in utilizing QR in crop improvement. However, the mechanisms governing QR are largely unknown, limiting its deployment. Advances in genomics are facilitating the dissection of the genetic and molecular underpinnings of QR, resulting in the discovery of several loci and genes that can be potentially deployed to enhance blackleg resistance. Here, we summarize the efforts undertaken to identify blackleg QR loci in oilseed rape using linkage and association analysis. We update the knowledge on the possible mechanisms governing QR and the advances in searching for the underlying genes. Lastly, we lay out strategies to accelerate the genetic improvement of blackleg QR in oilseed rape using improved phenotyping approaches and genomic prediction tools.

Meta-analysis of genome-wide association studies reveal common loci controlling agronomic and quality traits in a wide range of normal and heat stressed environments

Several stable QTL were detected using metaGWAS analysis for different agronomic and quality traits under 26 normal and heat stressed environments. Heat stress, exacerbated by global warming, has a negative influence on wheat production worldwide and climate resilient cultivars can help mitigate these impacts. Selection decisions should therefore depend on multi-environment experiments representing a range of temperatures at critical stages of development. Here, we applied a meta-genome wide association analysis (metaGWAS) approach to detect stable QTL with significant effects across multiple environments. The metaGWAS was applied to 11 traits scored in 26 trials that were sown at optimal or late times of sowing (TOS1 and TOS2, respectively) at five locations. A total of 2571 unique wheat genotypes (13,959 genotypes across all environments) were included and the analysis conducted on TOS1, TOS2 and both times of sowing combined (TOS1&2). The germplasm was genotyped using a 90 k Infinium chip and imputed to exome sequence level, resulting in 341,195 single nucleotide polymorphisms (SNPs). The average accuracy across all imputed SNPs was high (92.4%). The three metaGWAS analyses revealed 107 QTL for the 11 traits, of which 16 were detected in all three analyses and 23 were detected in TOS1&2 only. The remaining QTL were detected in either TOS1 or TOS2 with or without TOS1&2, reflecting the complex interactions between the environments and the detected QTL. Eight QTL were associated with grain yield and seven with multiple traits. The identified QTL provide an important resource for gene enrichment and fine mapping to further understand the mechanisms of gene × environment interaction under both heat stressed and unstressed conditions.

Genome-wide association study of multiple traits linked to heat tolerance in emmer-derived hexaploid wheat genotypes

Heat stress tolerance in plants is a complex trait controlled by multiple genes of minor effect which are influenced by the environment and this makes breeding and selection complicated. Emmer wheat (Triticum dicoccon Schrank) carries valuable diversity that can be used to improve the heat tolerance of modern bread wheat. A diverse set of emmer-based genotypes was developed by crossing emmer wheat with hexaploid wheat. These materials, along with their hexaploid recurrent parents and commercial cultivars, were evaluated at optimum (E1) and heat stressed (E2) sowing times in the field for three consecutive years (2014-2016). The material was genotyped using the Infinium iSelect SNP 90K SNP Assay. The phenotypic data were combined across years within each sowing time and best linear unbiased estimators calculated for each genotype in each environment. These estimates were used for GWAS analysis. Significant phenotypic and genotypic variation was observed for all traits. A total of 125 and 142 marker-trait associations (MTAs) were identified in E1 and E2, respectively. The highest number of MTAs were observed on the A genome (106), followed by the B (105) and D (56) genomes. MTAs with pleiotropic effects within and across the environments were observed. Many of the MTAs found were reported previously for various traits, and a few significant MTAs under heat stress were new and linked to emmer genome. Genomic regions identified on chromosomes 2B and 3A had a significant positive impact on grain yield under stress with a 7% allelic effect. Genomic regions on chromosomes 1A and 4B contributed 11% and 9% of the variation for thousand kernel weight (TKW) under heat stress respectively. Following fine mapping, these regions could be used for marker-assisted selection to improve heat tolerance in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01222-3.

The Use of Genetic and Gene Technologies in Shaping Modern Rapeseed Cultivars (Brassica napus L.)

Since their domestication, Brassica oilseed species have undergone progressive transformation allied with the development of breeding and molecular technologies. The canola (Brassica napus) crop has rapidly expanded globally in the last 30 years with intensive innovations in canola varieties, providing for a wider range of markets apart from the food industry. The breeding efforts of B. napus, the main source of canola oil and canola meal, have been mainly focused on improving seed yield, oil quality, and meal quality along with disease resistance, abiotic stress tolerance, and herbicide resistance. The revolution in genetics and gene technologies, including genetic mapping, molecular markers, genomic tools, and gene technology, especially gene editing tools, has allowed an understanding of the complex genetic makeup and gene functions in the major bioprocesses of the Brassicales, especially Brassica oil crops. Here, we provide an overview on the contributions of these technologies in improving the major traits of B. napus and discuss their potential use to accomplish new improvement targets.