Deciphering the Genetic Architecture of Plant Virus Resistance by GWAS, State of the Art and Potential Advances

Genome-wide association studies to assess genetic factors controlling cucumber resistance to CABYV and CMV in crop fields and the attractiveness for their Aphis gossypii vector

Cucumber crops face high pressure from pathogens, including various viral species. Mapping quantitative trait loci (QTL) for vegetable resistance to viruses has primarily been conducted after mechanical inoculation in controlled environments, but not in crop field conditions. Moreover, viruses that cannot be mechanically inoculated, e.g. the cucurbit aphid-borne yellows virus (CABYV), have been overlooked in resistance studies. Here, we aimed to identify QTLs reducing epidemics of two prevalent cucumber viruses: CABYV and the cucumber mosaic virus (CMV). We evaluated the resistance of 256 elite cucumber lines and landraces in crop field conditions by screening for the presence of both viruses six-times during the season. We mapped twelve QTLs reducing CABYV epidemics and seven QTLs reducing CMV epidemics by combining multiloci genome-wide association studies and local score approach analyses. We also examined the attractiveness of this cucumber panel for Aphis gossypii, a major cucumber virus vector. We identified five QTLs that reduced the attractiveness, including one co-localizing with a QTL reducing CABYV epidemics. Interestingly, some accessions deemed CMV-resistant after mechanical inoculation in controlled environments showed high infection rates in crop field conditions. Only one QTL for CMV resistance was detected in both conditions, indicating that these phenotypes are regulated by independent QTLs. Local linkage disequilibrium study findings suggested that certain QTLs reducing epidemics were introduced early into elite lines through serendipity or selection. QTLs could be pyramided with other low-effect QTLs through genomic selection to obtain cucumber cultivars with enhanced field resistance to viruses.

Genome-wide association study of cassava brown streak disease resistance in cassava germplasm conserved in South America

Cassava (Manihot esculenta Crantz) is a vital carbohydrate source for over 800 million people globally, yet its production in East Africa is severely affected by cassava brown streak disease (CBSD). Genebanks, through ex-situ conservation, play a pivotal role in preserving crop diversity, providing crucial resources for breeding resilient and disease-resistant crops. This study genotyped 234 South American cassava accessions conserved at the CIAT genebank, previously phenotyped for CBSD resistance by an independent group, to perform a genome-wide association analysis (GWAS) to identify genetic variants associated with CBSD resistance. Our GWAS identified 35 single nucleotide polymorphism (SNP) markers distributed across various chromosomes, associated with disease severity or the presence/absence of viral infection. Markers were annotated within or near genes previously identified with functions related to pathogen recognition and immune response activation. Using the SNP candidates, we screened the world's largest cassava collection for accessions with a higher frequency of favorable genotypes, proposing 35 accessions with potential resistance to CBSD. Our results provide insights into the genetics of CBSD resistance and highlight the importance of genetic resources to equip breeders with the raw materials needed to develop new crop varieties resistant to pests and diseases.

Yellow Dwarf Virus Resistance in Barley: Phenotyping, Remote Imagery, and Virus-Vector Characterization

Yellow dwarf viruses (YDVs) spread by aphids are some of the most economically important barley (Hordeum vulgare) virus-vector complexes worldwide. Detection and control of these viruses are critical components in the production of barley, wheat, and numerous other grasses of agricultural importance. Genetic control of plant diseases is often preferable to chemical control to reduce the environmental and economic cost of foliar insecticides. Accordingly, the objectives of this work were to (i) screen a barley population for resistance to YDVs under natural infection using phenotypic assessment of disease symptoms, (ii) implement drone imagery to further assess resistance and test its utility as a disease screening tool, (iii) identify the prevailing virus and vector types in the experimental environment, and (iv) perform a genome-wide association study to identify genomic regions associated with measured traits. Significant genetic differences were found in a population of 192 barley inbred lines regarding their YDV symptom severity, and symptoms were moderately to highly correlated with grain yield. The YDV severity measured with aerial imaging was highly correlated with on-the-ground estimates (r = 0.65). Three aphid species vectoring three YDV species were identified with no apparent genotypic influence on their distribution. A quantitative trait locus impacting YDV resistance was detected on chromosome 2H, albeit undetected using aerial imaging. However, quantitative trait loci for canopy cover and mean normalized difference vegetation index were successfully mapped using the drone. This work provides a framework for utilizing drone imagery in future resistance breeding efforts for YDVs in cereals and grasses, as well as in other virus-vector disease complexes.

The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Cauliflower mosaic virus disease spectrum uncovers novel susceptibility factor NCED9 in Arabidopsis thaliana

Viruses are intimately linked with their hosts and especially dependent on gene-for-gene interactions to establish successful infections. On the host side, defence mechanisms such as tolerance and resistance can occur within the same species, leading to differing virus accumulation in relation to symptomology and plant fitness. The identification of novel resistance genes against viruses and susceptibility factors is an important part of understanding viral patho-genesis and securing food production. The model plant Arabidopsis thaliana displays a wide symptom spectrum in response to RNA virus infections, and unbiased genome-wide association studies have proven a powerful tool to identify novel disease-genes. In this study we infected natural accessions of A. thaliana with the pararetrovirus cauliflower mosaic virus (CaMV) to study the phenotypic variations between accessions and their correlation with virus accumulation. Through genome-wide association mapping of viral accumulation differences, we identified several susceptibility factors for CaMV, the strongest of which was the abscisic acid synthesis gene NCED9. Further experiments confirmed the importance of abscisic acid homeostasis and its disruption for CaMV disease.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

Unravelling cucumber resistance to several viruses via genome-wide association studies highlighted resistance hotspots and new QTLs

The mapping and introduction of sustainable resistance to viruses in crops is a major challenge in modern breeding, especially regarding vegetables. We hence assembled a panel of cucumber elite lines and landraces from different horticultural groups for testing with six virus species. We mapped 18 quantitative trait loci (QTL) with a multiloci genome wide association studies (GWAS), some of which have already been described in the literature. We detected two resistance hotspots, one on chromosome 5 for resistance to the cucumber mosaic virus (CMV), cucumber vein yellowing virus (CVYV), cucumber green mottle mosaic virus (CGMMV) and watermelon mosaic virus (WMV), colocalizing with the RDR1 gene, and another on chromosome 6 for resistance to the zucchini yellowing mosaic virus (ZYMV) and papaya ringspot virus (PRSV) close to the putative VPS4 gene location. We observed clear structuring of resistance among horticultural groups due to plant virus coevolution and modern breeding which have impacted linkage disequilibrium (LD) in resistance QTLs. The inclusion of genetic structure in GWAS models enhanced the GWAS accuracy in this study. The dissection of resistance hotspots by local LD and haplotype construction helped gain insight into the panel's resistance introduction history. ZYMV and CMV resistance were both introduced from different donors in the panel, resulting in multiple resistant haplotypes at same locus for ZYMV, and in multiple resistant QTLs for CMV.