Structure of linkage disequilibrium and phenotypic associations in the maize genome

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Genome-wide association study identifies loci and candidate genes for grain micronutrients and quality traits in wheat (Triticum aestivum L.)

Malnutrition due to micronutrients and protein deficiency is recognized among the major global health issues. Genetic biofortification of wheat is a cost-effective and sustainable strategy to mitigate the global micronutrient and protein malnutrition. Genomic regions governing grain zinc concentration (GZnC), grain iron concentration (GFeC), grain protein content (GPC), test weight (TW), and thousand kernel weight (TKW) were investigated in a set of 184 diverse bread wheat genotypes through genome-wide association study (GWAS). The GWAS panel was genotyped using Breeders' 35 K Axiom Array and phenotyped in three different environments during 2019-2020. A total of 55 marker-trait associations (MTAs) were identified representing all three sub-genomes of wheat. The highest number of MTAs were identified for GPC (23), followed by TKW (15), TW (11), GFeC (4), and GZnC (2). Further, a stable SNP was identified for TKW, and also pleiotropic regions were identified for GPC and TKW. In silico analysis revealed important putative candidate genes underlying the identified genomic regions such as F-box-like domain superfamily, Zinc finger CCCH-type proteins, Serine-threonine/tyrosine-protein kinase, Histone deacetylase domain superfamily, and SANT/Myb domain superfamily proteins, etc. The identified novel MTAs will be validated to estimate their effects in different genetic backgrounds for subsequent use in marker-assisted selection.

Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.)

Powdery mildew (PM), caused by the fungus Oidium lini in flax, can cause defoliation and reduce seed yield and quality. To date, one major dominant gene (Pm1) and three quantitative trait loci (QTL) on chromosomes 1, 7 and 9 have been reported for PM resistance. To fully dissect the genetic architecture of PM resistance and identify QTL, a diverse flax core collection of 372 accessions augmented with an additional 75 breeding lines were sequenced, and PM resistance was evaluated in the field for eight years (2010–2017) in Morden, Manitoba, Canada. Genome-wide association studies (GWAS) were performed using two single-locus and seven multi-locus statistical models with 247,160 single nucleotide polymorphisms (SNPs) and the phenotypes of the 447 individuals for each year separately as well as the means over years. A total of 349 quantitative trait nucleotides (QTNs) were identified, of which 44 large-effect QTNs (R² = 10–30%) were highly stable over years. The total number of favourable alleles per accession was significantly correlated with PM resistance (r = 0.74), and genomic selection (GS) models using all identified QTNs generated significantly higher predictive ability (r = 0.93) than those constructed using the 247,160 genome-wide random SNP (r = 0.69), validating the overall reliability of the QTNs and showing the additivity of PM resistance in flax. The QTNs were clustered on the distal ends of all 15 chromosomes, especially on chromosome 5 (0.4–5.6 Mb and 9.4–16.9 Mb) and 13 (4.7–5.2 Mb). To identify candidate genes, a dataset of 3230 SNPs located in resistance gene analogues (RGAs) was used as input for GWAS, from which an additional 39 RGA-specific QTNs were identified. Overall, 269 QTN loci harboured 445 RGAs within the 200 Kb regions spanning the QTNs, including 45 QTNs located within the RGAs. These RGAs supported by significant QTN/SNP allele effects were mostly nucleotide binding site and leucine-rich repeat receptors (NLRs) belonging to either coiled-coil (CC) NLR (CNL) or toll interleukin-1 (TIR) NLR (TNL), receptor-like kinase (RLK), receptor-like protein kinase (RLP), transmembrane-coiled-coil (TM-CC), WRKY, and mildew locus O (MLO) genes. These results constitute an important genomic tool for resistance breeding and gene cloning for PM in flax.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

Variation in upstream open reading frames contributes to allelic diversity in maize protein abundance

SignificanceProteins are the machinery which execute essential cellular functions. However, measuring their abundance within an organism can be difficult and resource-intensive. Cells use a variety of mechanisms to control protein synthesis from mRNA, including short open reading frames (uORFs) that lie upstream of the main coding sequence. Ribosomes can preferentially translate uORFs instead of the main coding sequence, leading to reduced translation of the main protein. In this study, we show that uORF sequence variation between individuals can lead to different rates of protein translation and thus variable protein abundances. We also demonstrate that natural variation in uORFs occurs frequently and can be linked to whole-plant phenotypes, indicating that uORF sequence variation likely contributes to plant adaptation.

Genome-wide Association Study for Starch Pasting Properties in Chinese Spring Wheat

In order to understand the genetic basis of starch pasting viscosity characteristics of Chinese spring wheat, we assessed the genetic variation of RVA parameters determined by the Rapid Visco Analyser in a panel of 192 Chinese spring wheat accessions grown in Er'shi, Shihezi and Zhaosu during 2012 and 2013 cropping seasons. A genome-wide association study with 47,362 single nucleotide polymorphism (SNP) markers was conducted to detect marker-trait associations using mixed linear model. Phenotypic variations of RVA parameters ranged from 1.6 to 30.7% and broad-sense heritabilities ranged from 0.62 to 0.91. Forty-one SNP markers at 25 loci were significantly associated with seven RVA traits in at least two environments; among these, 20 SNPs were located in coding sequences (CDS) of 18 annotation genes, which can lead to discovering novel genes underpinning starch gelatinization in spring wheat. Haplotype analysis revealed one block for breakdown (BD) on chromosome 3B and two blocks for pasting temperature (T) on chromosome 7B. Cultivars with superior haplotypes at these loci showed better starch pasting viscosity than the average of all cultivars surveyed. The identified loci and associated markers provide valuable sources for future functional characterization and genetic improvement of starch quality in wheat.

Genetic Diversity, Population Structure and Linkage Disequilibrium Analyses in Tropical Maize Using Genotyping by Sequencing

Several maize breeding programs in India have developed numerous inbred lines but the lines have not been characterized using high-density molecular markers. Here, we studied the molecular diversity, population structure, and linkage disequilibrium (LD) patterns in a panel of 314 tropical normal corn, two sweet corn, and six popcorn inbred lines developed by 17 research centers in India, and 62 normal corn from the International Maize and Wheat Improvement Center (CIMMYT). The 384 inbred lines were genotyped with 60,227 polymorphic single nucleotide polymorphisms (SNPs). Most of the pair-wise relative kinship coefficients (58.5%) were equal or close to 0, which suggests the lack of redundancy in the genomic composition in the majority of inbred lines. Genetic distance among most pairs of lines (98.3%) varied from 0.20 to 0.34 as compared with just 1.7% of the pairs of lines that differed by <0.20, which suggests greater genetic variation even among sister lines. The overall average of 17% heterogeneity was observed in the panel indicated the need for further inbreeding in the high heterogeneous genotypes. The mean nucleotide diversity and frequency of polymorphic sites observed in the panel were 0.28 and 0.02, respectively. The model-based population structure, principal component analysis, and phylogenetic analysis revealed three to six groups with no clear patterns of clustering by centers-wise breeding lines, types of corn, kernel characteristics, maturity, plant height, and ear placement. However, genotypes were grouped partially based on their source germplasm from where they derived.

Genome-Wide Association Mapping Indicates Quantitative Genetic Control of Spot Blotch Resistance in Bread Wheat and the Favorable Effects of Some Spot Blotch Loci on Grain Yield

Spot blotch caused by the fungus Bipolaris sorokiniana poses a serious threat to bread wheat production in warm and humid wheat-growing regions of the world. Hence, the major objective of this study was to identify consistent genotyping-by-sequencing (GBS) markers associated with spot blotch resistance using genome-wide association mapping on a large set of 6,736 advanced bread wheat breeding lines from the International Maize and Wheat Improvement Center. These lines were phenotyped as seven panels at Agua Fria, Mexico between the 2013-2014 and 2019-2020 crop cycles. We identified 214 significant spot blotch associated GBS markers in all the panels, among which only 96 were significant in more than one panel, indicating a strong environmental effect on the trait and highlights the need for multiple phenotypic evaluations to identify lines with stable spot blotch resistance. The 96 consistent GBS markers were on chromosomes 1A, 1B, 1D, 2A, 3B, 4A, 5B, 5D, 6B, 7A, 7B, and 7D, including markers possibly linked to the Lr46, Sb1, Sb2 and Sb3 genes. We also report the association of the 2NS translocation from Aegilops ventricosa with spot blotch resistance in some environments. Moreover, the spot blotch favorable alleles at the 2NS translocation and two markers on chromosome 3BS (3B_2280114 and 3B_5601689) were associated with increased grain yield evaluated at several environments in Mexico and India, implying that selection for favorable alleles at these loci could enable simultaneous improvement for high grain yield and spot blotch resistance. Furthermore, a significant relationship between the percentage of favorable alleles in the lines and their spot blotch response was observed, which taken together with the multiple minor effect loci identified to be associated with spot blotch in this study, indicate quantitative genetic control of resistance. Overall, the results presented here have extended our knowledge on the genetic basis of spot blotch resistance in bread wheat and further efforts to improve genetic resistance to the disease are needed for reducing current and future losses under climate change.

Identification of Genomic Regions Associated with Agronomic and Disease Resistance Traits in a Large Set of Multiple DH Populations

Breeding maize lines with the improved level of desired agronomic traits under optimum and drought conditions as well as increased levels of resistance to several diseases such as maize lethal necrosis (MLN) is one of the most sustainable approaches for the sub-Saharan African region. In this study, 879 doubled haploid (DH) lines derived from 26 biparental populations were evaluated under artificial inoculation of MLN, as well as under well-watered (WW) and water-stressed (WS) conditions for grain yield and other agronomic traits. All DH lines were used for analyses of genotypic variability, association studies, and genomic predictions for the grain yield and other yield-related traits. Genome-wide association study (GWAS) using a mixed linear FarmCPU model identified SNPs associated with the studied traits i.e., about seven and eight SNPs for the grain yield; 16 and 12 for anthesis date; seven and eight for anthesis silking interval; 14 and 5 for both ear and plant height; and 15 and 5 for moisture under both WW and WS environments, respectively. Similarly, about 13 and 11 SNPs associated with gray leaf spot and turcicum leaf blight were identified. Eleven SNPs associated with senescence under WS management that had depicted drought-stress-tolerant QTLs were identified. Under MLN artificial inoculation, a total of 12 and 10 SNPs associated with MLN disease severity and AUDPC traits, respectively, were identified. Genomic prediction under WW, WS, and MLN disease artificial inoculation revealed moderate-to-high prediction accuracy. The findings of this study provide useful information on understanding the genetic basis for the MLN resistance, grain yield, and other agronomic traits under MLN artificial inoculation, WW, and WS conditions. Therefore, the obtained information can be used for further validation and developing functional molecular markers for marker-assisted selection and for implementing genomic prediction to develop superior elite lines.

Genomic Analysis of Resistance to Fall Armyworm (Spodoptera frugiperda) in CIMMYT Maize Lines

The recent invasion, rapid spread, and widescale destruction of the maize crop by the fall armyworm (FAW; Spodoptera frugiperda (J.E. Smith)) is likely to worsen the food insecurity situation in Africa. In the present study, a set of 424 maize lines were screened for their responses to FAW under artificial infestation to dissect the genetic basis of resistance. All lines were evaluated for two seasons under screen houses and genotyped with the DArTseq platform. Foliar damage was rated on a scale of 1 (highly resistant) to 9 (highly susceptible) and scored at 7, 14, and 21 days after artificial infestation. Analyses of variance revealed significant genotypic and genotype by environment interaction variances for all traits. Heritability estimates for leaf damage scores were moderately high and ranged from 0.38 to 0.58. Grain yield was negatively correlated with a high magnitude to foliar damage scores, ear rot, and ear damage traits. The genome-wide association study (GWAS) revealed 56 significant marker–trait associations and the predicted functions of the putative candidate genes varied from a defense response to several genes of unknown function. Overall, the study revealed that native genetic resistance to FAW is quantitative in nature and is controlled by many loci with minor effects.

Genetic diversity and population structure of rice (Oryza sativa L.) landraces from Kerala, India analyzed through genotyping-by-sequencing

Researchers stand at the vanguard of advancement and application of next-generation sequencing technology for developing dominant strategies for the sustainable management of genetically diverse crops. We attempt to fill the existing research lacuna in the molecular characterization of potent rice landraces in Kerala. Genotyping-by-sequencing (GBS) was performed on 96 Kerala rice accessions to identify single-nucleotide polymorphisms (SNPs), to examine the genetic diversity, population structure, and to delineate linkage disequilibrium (LD) pattern. GBS identified 5856 high-quality SNPs. The structure analysis indicated three subpopulations with the highest probability for population clustering with significant genetic differentiation, confirmed by principal component analysis. The genome-wide LD decay distance was 772 kb, at which the r² dropped to half its maximum value. The analysis of genetic properties of the identified SNP panel with an average polymorphism information content (PIC) value of 0.22 and a minor allele frequency (MAF) > 0.1 unveiled their efficacy in genome-wide association studies (GWAS). High F_ST (0.266) and low Nm (0.692) portray a strong genetic differentiation among the rice landraces, complementing the genetic structuring observed in the studied population. Slow LD decay in the rice landraces reflects their self-pollinating behavior and the indirect selection of desired traits by domestication. Moreover, the high LD entails only a minimum number of SNP markers for detecting marker-trait association. The diverse germplasm utilized in this study can be further utilized to disclose genetic variants associated with phenotypic traits and define signatures of selection via GWAS and selective sweep, respectively.

Phenotypic and molecular characterization of a set of tropical maize inbred lines from a public breeding program in Brazil

Background

The characterization of genetic diversity and population differentiation for maize inbred lines from breeding programs is of great value in assisting breeders in maintaining and potentially increasing the rate of genetic gain. In our study, we characterized a set of 187 tropical maize inbred lines from the public breeding program of the Universidade Federal de Viçosa (UFV) in Brazil based on 18 agronomic traits and 3,083 single nucleotide polymorphisms (SNP) markers to evaluate whether this set of inbred lines represents a panel of tropical maize inbred lines for association mapping analysis and investigate the population structure and patterns of relationships among the inbred lines from UFV for better exploitation in our maize breeding program.

Results

Our results showed that there was large phenotypic and genotypic variation in the set of tropical maize inbred lines from the UFV maize breeding program. We also found high genetic diversity (GD = 0.34) and low pairwise kinship coefficients among the maize inbred lines (only approximately 4.00 % of the pairwise relative kinship was above 0.50) in the set of inbred lines. The LD decay distance over all ten chromosomes in the entire set of maize lines with r² = 0.1 was 276,237 kb. Concerning the population structure, our results from the model-based STRUCTURE and principal component analysis methods distinguished the inbred lines into three subpopulations, with high consistency maintained between both results. Additionally, the clustering analysis based on phenotypic and molecular data grouped the inbred lines into 14 and 22 genetic divergence clusters, respectively.

Conclusions

Our results indicate that the set of tropical maize inbred lines from UFV maize breeding programs can comprise a panel of tropical maize inbred lines suitable for a genome-wide association study to dissect the variation of complex quantitative traits in maize, mainly in tropical environments. In addition, our results will be very useful for assisting us in the assignment of heterotic groups and the selection of the best parental combinations for new breeding crosses, mapping populations, mapping synthetic populations, guiding crosses that target highly heterotic and yielding hybrids, and predicting untested hybrids in the public breeding program UFV.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08127-7.

Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut (Arachis hypogaea L.)

The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here, we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils. Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect-increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes. With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanuts, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold (α = 1%). We detected 7682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (P-value < 1.0E-16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Development of high-resolution multiple-SNP arrays for genetic analyses and molecular breeding through genotyping by target sequencing and liquid chip

Genotyping platforms, as critical supports for genomics, genetics, and molecular breeding, have been well implemented at national institutions/universities in developed countries and multinational seed companies that possess high-throughput, automatic, large-scale, and shared facilities. In this study, we integrated an improved genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology to develop a multiple single-nucleotide polymorphism (mSNP) approach in which mSNPs can be captured from a single amplicon. From one 40K maize mSNP panel, we developed three types of markers (40K mSNPs, 251K SNPs, and 690K haplotypes), and generated multiple panels with various marker densities (1K-40K mSNPs) by sequencing at different depths. Comparative genetic diversity analysis was performed with genic versus intergenic markers and di-allelic SNPs versus non-typical SNPs. Compared with the one-amplicon-one-SNP system, mSNPs and within-mSNP haplotypes are more powerful for genetic diversity detection, linkage disequilibrium decay analysis, and genome-wide association studies. The technologies, protocols, and application scenarios developed for maize in this study will serve as a model for the development of mSNP arrays and highly efficient GBTS systems in animals, plants, and microorganisms.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Genome-wide association mapping uncovers sex-associated copy number variation markers and female hemizygous regions on the W chromosome in Salix viminalis

Background

Sex chromosomes are in some species largely undifferentiated (homomorphic) with restricted sex determination regions. Homomorphic but different sex chromosomes are found in the closely related genera Populus and Salix indicating flexible sex determination systems, ideal for studies of processes involved in sex chromosome evolution. We have performed genome-wide association studies of sex and analysed sex chromosomes in a population of 265 wild collected Salix viminalis accessions and studied the sex determining locus.

Results

A total of 19,592 markers were used in association analyses using both Fisher’s exact tests and a single-marker mixed linear model, which resulted in 48 and 41 sex-associated (SA) markers respectively. Across all 48 SA markers, females were much more often heterozygous than males, which is expected if females were the heterogametic sex. The majority of the SA markers were, based on positions in the S. purpurea genome, located on chromosome 15, previously demonstrated to be the sex chromosome. Interestingly, when mapping the genotyping-by-sequencing sequence tag harbouring the two SA markers with the highest significance to the S. viminalis genomic scaffolds, five regions of very high similarity were found: three on a scaffold that represents a part of chromosome 15, one on a scaffold that represents a part of chromosome 9 and one on a scaffold not anchored to the genome. Based on segregation differences of the alleles at the two marker positions and on differences in PCR amplification between females and males we conclude that females had multiple copies of this DNA fragment (chromosome 9 and 15), whereas males only had one (chromosome 9). We therefore postulate that the female specific sequences have been copied from chromosome 9 and inserted on chromosome 15, subsequently developing into a hemizygous W chromosome linked region.

Conclusions

Our results support that sex determination in S. viminalis is controlled by one locus on chromosome 15. The segregation patterns observed at the SA markers furthermore confirm that S. viminalis females are the heterogametic sex. We also identified a translocation from chromosome 9 to the W chromosome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08021-2.

Gene-based mapping of trehalose biosynthetic pathway genes reveals association with source- and sink-related yield traits in a spring wheat panel

Trehalose 6‐phosphate (T6P) signalling regulates carbon use and allocation and is a target to improve crop yields. However, the specific contributions of trehalose phosphate synthase (TPS) and trehalose phosphate phosphatase (TPP) genes to source‐ and sink‐related traits remain largely unknown. We used enrichment capture sequencing on TPS and TPP genes to estimate and partition the genetic variation of yield‐related traits in a spring wheat (Triticum aestivum) breeding panel specifically built to capture the diversity across the 75,000 CIMMYT wheat cultivar collection. Twelve phenotypes were correlated to variation in TPS and TPP genes including plant height and biomass (source), spikelets per spike, spike growth and grain filling traits (sink) which showed indications of both positive and negative gene selection. Individual genes explained proportions of heritability for biomass and grain‐related traits. Three TPS1 homologues were particularly significant for trait variation. Epistatic interactions were found within and between the TPS and TPP gene families for both plant height and grain‐related traits. Gene‐based prediction improved predictive ability for grain weight when gene effects were combined with the whole‐genome markers. Our study has generated a wealth of information on natural variation of TPS and TPP genes related to yield potential which confirms the role for T6P in resource allocation and in affecting traits such as grain number and size confirming other studies which now opens up the possibility of harnessing natural genetic variation more widely to better understand the contribution of native genes to yield traits for incorporation into breeding programmes.

Population-level deep sequencing reveals the interplay of clonal and sexual reproduction in the fungal wheat pathogen Zymoseptoria tritici

Pathogens cause significant challenges to global food security. On annual crops, pathogens must re-infect from environmental sources in every growing season. Fungal pathogens have evolved mixed reproductive strategies to cope with the distinct challenges of colonizing growing plants. However, how pathogen diversity evolves during growing seasons remains largely unknown. Here, we performed a deep hierarchical sampling in a single experimental wheat field infected by the major fungal pathogen Zymoseptoria tritici. We analysed whole genome sequences of 177 isolates collected from 12 distinct cultivars replicated in space at three time points of the growing season to maximize capture of genetic diversity. The field population was highly diverse with 37 SNPs per kilobase, a linkage disequilibrium decay within 200-700 bp and a high effective population size. Using experimental infections, we tested a subset of the collected isolates on the dominant cultivar planted in the field. However, we found no significant difference in virulence of isolates collected from the same cultivar compared to isolates collected on other cultivars. About 20 % of the isolate genotypes were grouped into 15 clonal groups. Pairs of clones were disproportionally found at short distances (<5 m), consistent with experimental estimates for per-generation dispersal distances performed in the same field. This confirms predominant leaf-to-leaf transmission during the growing season. Surprisingly, levels of clonality did not increase over time in the field although reproduction is thought to be exclusively asexual during the growing season. Our study shows that the pathogen establishes vast and stable gene pools in single fields. Monitoring short-term evolutionary changes in crop pathogens will inform more durable strategies to contain diseases.

Genetic Dissection of Grain Yield Component Traits Under High Nighttime Temperature Stress in a Rice Diversity Panel

To dissect the genetic complexity of rice grain yield (GY) and quality in response to heat stress at the reproductive stage, a diverse panel of 190 rice accessions in the United States Department of Agriculture (USDA) rice mini-core collection (URMC) diversity panel were treated with high nighttime temperature (HNT) stress at the reproductive stage of panicle initiation. The quantifiable yield component response traits were then measured. The traits, panicle length (PL), and number of spikelets per panicle (NSP) were evaluated in subsets of the panel comprising the rice subspecies Oryza sativa ssp. Indica and ssp. Japonica. Under HNT stress, the Japonica ssp. exhibited lower reductions in PL and NSP and a higher level of genetic variation compared with the other subpopulations. Whole genome sequencing identified 6.5 million single nucleotide polymorphisms (SNPs) that were used for the genome-wide association studies (GWASs) of the PL and NSP traits. The GWAS analysis in the Combined, Indica, and Japonica populations under HNT stress identified 83, 60, and 803 highly significant SNPs associated with PL, compared to the 30, 30, and 11 highly significant SNPs associated with NSP. Among these trait-associated SNPs, 140 were coincident with genomic regions previously reported for major GY component quantitative trait loci (QTLs) under heat stress. Using extents of linkage disequilibrium in the rice populations, Venn diagram analysis showed that the highest number of putative candidate genes were identified in the Japonica population, with 20 putative candidate genes being common in the Combined, Indica and Japonica populations. Network analysis of the genes linked to significant SNPs associated with PL and NSP identified modules that were involved in primary and secondary metabolisms. The findings in this study could be useful to understand the pathways/mechanisms involved in rice GY and its components under HNT stress for the acceleration of rice-breeding programs and further functional analysis by molecular geneticists.

Introgression of Maize Diversity for Drought Tolerance: Subtropical Maize Landraces as Source of New Positive Variants

Current climate change models predict an increased frequency and intensity of drought for much of the developing world within the next 30 years. These events will negatively affect maize yields, potentially leading to economic and social instability in many smallholder farming communities. Knowledge about the genetic resources available for traits related to drought tolerance has great importance in developing breeding program strategies. The aim of this research was to study a maize landrace introgression panel to identify chromosomal regions associated with a drought tolerance index. For that, we performed Genome-Wide Association Study (GWAS) on 1326 landrace progenies developed by the CIMMYT Genetic Resources Program, originating from 20 landraces populations collected in arid regions. Phenotypic data were obtained from early testcross trials conducted in three sites and two contrasting irrigation environments, full irrigation (well-watered) and reduced irrigation (drought). The populations were genotyped using the DArTSeq^® platform, and a final set of 5,695 SNPs markers was used. The genotypic values were estimated using spatial adjustment in a two-stage analysis. First, we performed the individual analysis for each site/irrigation treatment combination. The best linear unbiased estimates (BLUEs) were used to calculate the Harmonic Mean of Relative Performance (HMRP) as a drought tolerance index for each testcross. The second stage was a joint analysis, which was performed using the HMRP to obtain the best linear unbiased predictions (BLUPs) of the index for each genotype. Then, GWAS was performed to determine the marker-index associations and the marker-Grain Yield (GY) associations for the two irrigation treatments. We detected two significant markers associated with the drought-tolerance index, four associated with GY in drought condition, and other four associated with GY in irrigated conditions each. Although each of these markers explained less than 0.1% of the phenotypic variation for the index and GY, we found two genes likely related to the plant response to drought stress. For these markers, alleles from landraces provide a slightly higher yield under drought conditions. Our results indicate that the positive diversity delivered by landraces are still present on the backcrosses and this is a potential breeding strategy for improving maize for drought tolerance and for trait introgression bringing new superior allelic diversity from landraces to breeding populations.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

The genetic basis of adaptation in phenology in an introduced population of Black Cottonwood (Populus trichocarpa, Torr. & Gray)

BACKGROUND

Entering and exiting winter dormancy present important trade-offs between growth and survival at northern latitudes. Many forest trees display local adaptation across latitude in traits associated with these phenology transitions. Transfers of a species outside its native range introduce the species to novel combinations of environmental conditions potentially requiring different combinations of alleles to optimize growth and survival. In this study, we performed genome wide association analyses and a selection scan in a P. trichocarpa mapping population derived from crossings between clones collected across the native range and introduced into Sweden. GWAS analyses were performed using phenotypic data collected across two field seasons and in a controlled phytotron experiment.

RESULTS

We uncovered 584 putative candidate genes associated with spring and autumn phenology traits as well as with growth. Many regions harboring variation significantly associated with the initiation of leaf shed and leaf autumn coloring appeared to have been evolving under positive selection in the native environments of P. trichocarpa. A comparison between the candidate genes identified with results from earlier GWAS analyses performed in the native environment found a smaller overlap for spring phenology traits than for autumn phenology traits, aligning well with earlier observations that spring phenology transitions have a more complex genetic basis than autumn phenology transitions.

CONCLUSIONS

In a small and structured introduced population of P. trichocarpa, we find complex genetic architectures underlying all phenology and growth traits, and identify multiple putative candidate genes despite the limitations of the study population.

Genomic regions associated with heat stress tolerance in tropical maize (Zea mays L.)

With progressive climate change and the associated increase in mean temperature, heat stress tolerance has emerged as one of the key traits in the product profile of the maize breeding pipeline for lowland tropics. The present study aims to identify the genomic regions associated with heat stress tolerance in tropical maize. An association mapping panel, called the heat tolerant association mapping (HTAM) panel, was constituted by involving a total of 543 tropical maize inbred lines from diverse genetic backgrounds, test-crossed and phenotyped across nine locations in South Asia under natural heat stress. The panel was genotyped using a genotyping-by-sequencing (GBS) platform. Considering the large variations in vapor pressure deficit (VPD) at high temperature (Tmax) across different phenotyping locations, genome-wide association study (GWAS) was conducted separately for each location. The individual location GWAS identified a total of 269 novel significant single nucleotide polymorphisms (SNPs) for grain yield under heat stress at a p value of < 10^–5. A total of 175 SNPs were found in 140 unique gene models implicated in various biological pathway responses to different abiotic stresses. Haplotype trend regression (HTR) analysis of the significant SNPs identified 26 haplotype blocks and 96 single SNP variants significant across one to five locations. The genomic regions identified based on GWAS and HTR analysis considering genomic region x environment interactions are useful for breeding efforts aimed at developing heat stress resilient maize cultivars for current and future climatic conditions through marker-assisted introgression into elite genetic backgrounds and/or genome-wide selection.

Spatial pattern of genetic diversity in field populations of Fusarium incarnatum-equiseti species complex

Fusarium is associated with a number of wilt, blight, scab, and rot diseases in a range of economically important staple food crops worldwide. An assessment of the genetic structure and population stratification of Fusarium incarnatum-equiseti species complex (FIESC) pathogen populations is important to understand the evolutionary potential of such populations in adapting to environmental change. Based on intersimple sequence repeat polymerase chain reaction (ISSR-PCR), it was found that the pathogen population was structured into three genetic clusters for which genetic differentiation was higher within than among populations. There was high intrapopulation genetic diversity for population 1 (94.63%) which consisted largely of isolates collected from North Trinidad. Populations 2 and 3 had a low level of admixture among the populations based on overall population differentiation. Population 1 accounted for the highest amount of genetic variation (95.82%) followed by populations 2 and 3. Population stratification was reflected in the dendrogram topology, which consisted of three main genetic clusters and which coincided with the outcome of Bayesian and PCoA analyses. The populations were isolated by distance, and Voronoi tessellations indicated physical or structural barriers to gene flow which contributed to restricted admixture between two of three populations. These findings suggest a high evolutionary potential for this FIESC pathogen population, the implications of which directly affect disease management strategies.

Polyphenol oxidase genes as integral part of the evolutionary history of domesticated tetraploid wheat

Studying and understanding the genetic basis of polyphenol oxidases (PPO)-related traits plays a crucial role in genetic improvement of crops. A tetraploid wheat collection (T. turgidum ssp., TWC) was analyzed using the 90K wheat SNP iSelect assay and phenotyped for PPO activity. A total of 21,347 polymorphic SNPs were used to perform genome-wide association analysis (GWA) in TWC and durum wheat sub-groups, detecting 23 and 85 marker-trait associations (MTA). In addition, candidate genes responsible for PPO activity were predicted. Based on the 23 MTAs detected in TWC, two haplotypes associated with low and high PPO activity were identified. Four SNPs were developed and validated providing one reliable marker (IWB75732) for marker assisted selection. The 23 MTAs were used to evaluate the genetic divergence (F_ST > 0.25) between the T. turgidum subspecies, providing new information important for understanding the domestication process of Triticum turgidum ssp. and in particular of ssp. carthlicum.

Genetic diversity and selection signatures in maize landraces compared across 50 years of in situ and ex situ conservation

Genomics-based, longitudinal comparisons between ex situ and in situ agrobiodiversity conservation strategies can contribute to a better understanding of their underlying effects. However, landrace designations, ambiguous common names, and gaps in sampling information complicate the identification of matching ex situ and in situ seed lots. Here we report a 50-year longitudinal comparison of the genetic diversity of a set of 13 accessions from the state of Morelos, Mexico, conserved ex situ since 1967 and retrieved in situ from the same donor families in 2017. We interviewed farmer families who donated in situ landraces to understand their germplasm selection criteria. Samples were genotyped by sequencing, producing 74,739 SNPs. Comparing the two sample groups, we show that ex situ and in situ genome-wide diversity was similar. In situ samples had 3.1% fewer SNPs and lower pairwise genetic distances (F_st 0.008-0.113) than ex situ samples (F_st 0.031-0.128), but displayed the same heterozygosity. Despite genome-wide similarities across samples, we could identify several loci under selection when comparing in situ and ex situ seed lots, suggesting ongoing evolution in farmer fields. Eight loci in chromosomes 3, 5, 6, and 10 showed evidence of selection in situ that could be related with farmers' selection criteria surveyed with focus groups and interviews at the sampling site in 2017, including wider kernels and larger ear size. Our results have implications for ex situ collection resampling strategies and the in situ conservation of threatened landraces.

Toward the identification of molecular markers associated with phytochemical traits in the Iranian sumac (Rhus coriaria L.) population

Sumac (Rhus coriaria L.) is one of the important forest species dispersed in the northwest of Iran. It is one of the important spice in Iran and the Middle East because of active components containing organic acids, phenolic acids, flavonoids, anthocyanins, tannins and terpenoids. This study aimed to investigate population structure and linkage disequilibrium (LD) extent within Rhus coriaria L. genotypes using ISSR markers and identify molecular markers associated with phytochemical traits using association analysis. In the molecular part of the experiment, the genetic diversity of 75 sumac genotypes from five different areas of northwest Iran was assessed by 18 ISSR primers. In the phenotypic assessment, the fruits of the sumac genotypes were analyzed using HPLC-LC/MSMS for determining phytochemical components including maleic acid, ellagic acid, maleic acid hexoside, gallic acid, coumaric acid, quercetin, caftaric acid, and linoleic acid. The phenotypic data analysis revealed the great phenotypic diversity among and within Iranian sumac populations for the studied phytochemical traits. The studied sumac genotypes were divided into two subpopulations based on molecular marker-based structure analysis. A significant level of LD was observed in 11.64% of the ISSR marker pairs (p < .05). The mixed linear model procedure showed that 12 loci had a significant association with investigated traits. The ISSR loci identified in this study can be potentially used in marker-assisted selection in sumac breeding programs.

Genome-wide association mapping identifies yellow rust resistance loci in Ethiopian durum wheat germplasm

Durum wheat is an important cereal grown in Ethiopia, a country which is also its center for genetic diversity. Yellow (stripe) rust caused by Puccinia striiformis fsp tritici is one of the most devastating diseases threatening Ethiopian wheat production. To identify sources of genetic resistance and combat this pathogen, we conducted a genome wide association study of yellow rust resistance on 300 durum wheat accessions comprising 261 landraces and 39 cultivars. The accessions were evaluated for their field resistance using a modified Cobb scale at Meraro, Kulumsa and Chefe Donsa in the 2015 and 2016 main growing seasons. Analysis of the 35K Axiom Array genotyping data of the panel resulted in a total of 8,797 polymorphic SNPs of which 7,093 were used in subsequent analyses. Population structure analysis suggested two groups in which the cultivars clearly stood out separately from the landraces. Eleven SNPs significantly associated with yellow rust resistance were identified on four chromosomes (1A, 1B, 2B, and 5A) which defined at least five genomic loci. Six of the SNPs were consistently identified on chromosome 1B singly at each and combined overall environments which explained 62.6–64.0% of the phenotypic variation (R²). Resistant allele frequency ranged from 14.0–71.0%; Zooming in to the identified resistance loci revealed the presence of disease resistance related genes involved in the plant defense system such as the ABC transporter gene family, disease resistance protein RPM1 (NBS-LRR class), Receptor kinases and Protein kinases. This study has provided SNPs for tracking the loci associated with yellow rust resistance and a diversity panel which can be used for association study of other agriculturally important traits in durum wheat.

Single nucleotide polymorphisms reveal genetic diversity in New Mexican chile peppers (Capsicum spp.)

BACKGROUND

Chile peppers (Capsicum spp.) are among the most important horticultural crops in the world due to their number of uses. They are considered a major cultural and economic crop in the state of New Mexico in the United States. Evaluating genetic diversity in current New Mexican germplasm would facilitate genetic improvement for different traits. This study assessed genetic diversity, population structure, and linkage disequilibrium (LD) among 165 chile pepper genotypes using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS).

RESULTS

A GBS approach identified 66,750 high-quality SNP markers with known map positions distributed across the 12 chromosomes of Capsicum. Principal components analysis revealed four distinct clusters based on species. Neighbor-joining phylogenetic analysis among New Mexico State University (NMSU) chile pepper cultivars showed two main clusters, where the C. annuum genotypes grouped together based on fruit or pod type. A Bayesian clustering approach for the Capsicum population inferred K = 2 as the optimal number of clusters, where the C. chinense and C. frutescens grouped in a single cluster. Analysis of molecular variance revealed majority of variation to be between the Capsicum species (76.08 %). Extensive LD decay (~ 5.59 Mb) across the whole Capsicum population was observed, demonstrating that a lower number of markers would be required for implementing genome wide association studies for different traits in New Mexican type chile peppers. Tajima's D values demonstrated positive selection, population bottleneck, and balancing selection for the New Mexico Capsicum population. Genetic diversity for the New Mexican chile peppers was relatively low, indicating the need to introduce new alleles in the breeding program to broaden the genetic base of current germplasm.

CONCLUSIONS

Genetic diversity among New Mexican chile peppers was evaluated using GBS-derived SNP markers and genetic relatedness on the species level was observed. Introducing novel alleles from other breeding programs or from wild species could help increase diversity in current germplasm. We present valuable information for future association mapping and genomic selection for different traits for New Mexican chile peppers for genetic improvement through marker-assisted breeding.

Development of an Australian Bread Wheat Nested Association Mapping Population, a New Genetic Diversity Resource for Breeding under Dry and Hot Climates

Genetic diversity, knowledge of the genetic architecture of the traits of interest and efficient means of transferring the desired genetic diversity into the relevant genetic background are prerequisites for plant breeding. Exotic germplasm is a rich source of genetic diversity; however, they harbor undesirable traits that limit their suitability for modern agriculture. Nested association mapping (NAM) populations are valuable genetic resources that enable incorporation of genetic diversity, dissection of complex traits and providing germplasm to breeding programs. We developed the OzNAM by crossing and backcrossing 73 diverse exotic parents to two Australian elite varieties Gladius and Scout. The NAM parents were genotyped using the iSelect wheat 90K Infinium SNP array, and the progeny were genotyped using a custom targeted genotyping-by-sequencing assay based on molecular inversion probes designed to target 12,179 SNPs chosen from the iSelect wheat 90K Infinium SNP array of the parents. In total, 3535 BC1F4:6 RILs from 125 families with 21 to 76 lines per family were genotyped and we found 4964 polymorphic and multi-allelic haplotype markers that spanned the whole genome. A subset of 530 lines from 28 families were evaluated in multi-environment trials over three years. To demonstrate the utility of the population in QTL mapping, we chose to map QTL for maturity and plant height using the RTM-GWAS approach and we identified novel and known QTL for maturity and plant height.

Sex determination through X-Y heterogamety in Salix nigra

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

Association mapping reveals novel genomic regions controlling some root and stolon traits in tetraploid potato (Solanum tuberosum L.)

Tuber crops have measurable biological variation in root and stolon phenotyping and thus may be utilized to identify genomic regions associated with these variations. This is the first comprehensive association mapping study related to potato root and stolon traits. A diverse panel of 192 tetraploid potato (Solanum tuberosum L.) genotypes were grown in aeroponics to reveal a biologically significant variation and detection of genomic regions associated with the root and stolon traits. Phenotyping of root traits was performed by image analysis software "WinRHIZO" (a root scanning method), and stolon traits was measured manually, while SolCAP 25K potato array was used for genotyping. Significant variation was observed between the potato genotypes for root and stolon traits along with high heritabilities (0.80 in TNS to 0.95 in SL). For marker-trait associations, Q + K linear mixed model was implemented and 50 novel genomic regions were detected. Significantly associated SNPs with stolon traits were located on chr 4, chr 6, chr 7, chr 9, chr 11 and chr 12, while those linked to root traits on chr 1, chr 2, chr 3, chr 9, chr 11, and chr 12. Structure and PCA analysis grouped genotypes into four sub-populations disclosing population genetic diversity. LD decay was observed at 2.316 Mbps (r 2 = 0.29) in the population. The identified SNPs were associated with genes performing vital functions such as root signaling and signal transduction in stress environments (GT-2 factors, protein kinases SAPK2-like and protein phosphatases "StPP1"), transcriptional and post-transcriptional gene regulation (RNA-binding proteins), sucrose synthesis and transporter families (UGPase, Sus3, SuSy, and StSUT1) and PVY resistance (Ry sto). The findings of our study can be employed in future breeding programs for improvement in potato production.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02727-6.

Genome-Wide Association of Stem Carbohydrate Accumulation and Remobilization during Grain Growth in Bread Wheat (Triticum aestivum L.) in Mediterranean Environments

Water deficit represents an important challenge for wheat production in many regions of the world. Accumulation and remobilization of water-soluble carbohydrates (WSCs) in stems are part of the physiological responses regulated by plants to cope with water stress and, in turn, determine grain yield (GY). The genetic mechanisms underlying the variation in WSC are only partially understood. In this study, we aimed to identify Single Nucleotide Polymorphism (SNP) markers that account for variation in a suite of WSC and GY, evaluated in 225 cultivars and advanced lines of spring wheat. These genotypes were established in two sites in the Mediterranean region of Central Chile, under water-limited and full irrigation conditions, and assessed in two growing seasons, namely anthesis and maturity growth periods. A genome-wide association study (GWAS) was performed by using 3243 SNP markers. Genetic variance accounted for 5 to 52% of phenotypic variation of the assessed traits. A rapid linkage disequilibrium decay was observed across chromosomes (r² ≤ 0.2 at 2.52 kbp). Marker-trait association tests identified 96 SNPs related to stem weight (SW), WSCs, and GY, among other traits, at the different sites, growing seasons, and growth periods. The percentage of SNPs that were part of the gene-coding regions was 34%. Most of these genes are involved in the defensive response to drought and biotic stress. A complimentary analysis detected significant effects of different haplotypes on WSC and SW, in anthesis and maturity. Our results evidence both genetic and environmental influence on WSC dynamics in spring wheat. At the same time, they provide a series of markers suitable for supporting assisted selection approaches and functional characterization of genes.

Elucidating the genetics of grain yield and stress-resilience in bread wheat using a large-scale genome-wide association mapping study with 55,568 lines

Wheat grain yield (GY) improvement using genomic tools is important for achieving yield breakthroughs. To dissect the genetic architecture of wheat GY potential and stress-resilience, we have designed this large-scale genome-wide association study using 100 datasets, comprising 105,000 GY observations from 55,568 wheat lines evaluated between 2003 and 2019 by the International Maize and Wheat Improvement Center and national partners. We report 801 GY-associated genotyping-by-sequencing markers significant in more than one dataset and the highest number of them were on chromosomes 2A, 6B, 6A, 5B, 1B and 7B. We then used the linkage disequilibrium (LD) between the consistently significant markers to designate 214 GY-associated LD-blocks and observed that 84.5% of the 58 GY-associated LD-blocks in severe-drought, 100% of the 48 GY-associated LD-blocks in early-heat and 85.9% of the 71 GY-associated LD-blocks in late-heat, overlapped with the GY-associated LD-blocks in the irrigated-bed planting environment, substantiating that simultaneous improvement for GY potential and stress-resilience is feasible. Furthermore, we generated the GY-associated marker profiles and analyzed the GY favorable allele frequencies for a large panel of 73,142 wheat lines, resulting in 44.5 million datapoints. Overall, the extensive resources presented in this study provide great opportunities to accelerate breeding for high-yielding and stress-resilient wheat varieties.

Association analysis of physiological traits in spring barley (Hordeum vulgare L.) under water-deficit conditions

In the present study, 148 commercial barley cultivars were assessed by 14 AFLP primer combinations and 32 SSRs primer pairs. Population structure, linkage disequilibrium, and genomic regions associated with physiological traits under drought stress were investigated. The phenotypic results showed a high level of diversity between studied cultivars. The studied barley cultivars were divided into two subgroups. Linkage disequilibrium analysis revealed that r 2 values among all possible marker pairs have an average value of 0.0178. The mixed linear model procedure showed that totally, 207 loci had a significant association with investigated traits. 120 QTLs out of 207 were detected for traits under normal conditions, and 90 QTLs were detected for traits under drought stress conditions. Identified QTLs after validation and transferring to SCAR markers in the case of AFLPs can be used to develop MAS strategies for barley breeding programs. Some common markers were identified for a particular trait or some traits across normal and drought stress conditions. These markers show low interaction with environmental conditions (stable markers); therefore, selection by them for a trait under normal conditions will improve the trait value under stress conditions, too.

Genetic and molecular control of grain yield in maize

Understanding the genetic and molecular basis of grain yield is important for maize improvement. Here, we identified 49 consensus quantitative trait loci (cQTL) controlling maize yield-related traits using QTL meta-analysis. Then, we collected yield-related traits associated SNPs detected by association mapping and identified 17 consensus significant loci. Comparing the physical positions of cQTL with those of significant SNPs revealed that 47 significant SNPs were located within 20 cQTL regions. Furthermore, intensive reviews of 31 genes regulating maize yield-related traits found that the functions of many genes were conservative in maize and other plant species. The functional conservation indicated that some of the 575 maize genes (orthologous to 247 genes controlling yield or seed traits in other plant species) might be functionally related to maize yield-related traits, especially the 49 maize orthologous genes in cQTL regions, and 41 orthologous genes close to the physical positions of significant SNPs. In the end, we prospected on the integration of the public sources for exploring the genetic and molecular mechanisms of maize yield-related traits, and on the utilization of genetic and molecular mechanisms for maize improvement.

Supplementary Information

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

Association mapping for maize stover yield and saccharification efficiency using a multiparent advanced generation intercross (MAGIC) population

Cellulosic ethanol derived from fast growing C4 grasses could become an alternative to finite fossil fuels. With the potential to generate a major source of lignocellulosic biomass, maize has gained importance as an outstanding model plant for studying the complex cell wall network and also to optimize crop breeding strategies in bioenergy grasses. A genome-wide association study (GWAS) was conducted using a subset of 408 Recombinant Inbred Lines (RILs) from a Multi-Parent Advanced Generation Intercross (MAGIC) Population in order to identify single nucleotide polymorphisms (SNPs) associated with yield and saccharification efficiency of maize stover. We identified 13 SNPs significantly associated with increased stover yield that corresponded to 13 QTL, and 2 SNPs significantly associated with improved saccharification efficiency, that could be clustered into 2 QTL. We have pointed out the most interesting SNPs to be implemented in breeding programs based on results from analyses of averaged and yearly data. Association mapping in this MAGIC population highlight genomic regions directly linked to traits that influence the final use of maize. Markers linked to these QTL could be used in genomic or marker-assisted selection programs to improve biomass quality for ethanol production. This study opens a possible optimisation path for improving the viability of second-generation biofuels.

Temporal dynamics of QTL effects on vegetative growth in Arabidopsis thaliana

We assessed early vegetative growth in a population of 382 accessions of Arabidopsis thaliana using automated non-invasive high-throughput phenotyping. All accessions were imaged daily from 7 d to 18 d after sowing in three independent experiments and genotyped using the Affymetrix 250k SNP array. Projected leaf area (PLA) was derived from image analysis and used to calculate relative growth rates (RGRs). In addition, initial seed size was determined. The generated datasets were used jointly for a genome-wide association study that identified 238 marker-trait associations (MTAs) individually explaining up to 8% of the total phenotypic variation. Co-localization of MTAs occurred at 33 genomic positions. At 21 of these positions, sequential co-localization of MTAs for 2-9 consecutive days was observed. The detected MTAs for PLA and RGR could be grouped according to their temporal expression patterns, emphasizing that temporal variation of MTA action can be observed even during the vegetative growth phase, a period of continuous formation and enlargement of seemingly similar rosette leaves. This indicates that causal genes may be differentially expressed in successive periods. Analyses of the temporal dynamics of biological processes are needed to gain important insight into the molecular mechanisms of growth-controlling processes in plants.

Genomic prediction using training population design in interspecific soybean populations

Agronomically important traits generally have complex genetic architecture, where many genes have a small and largely additive effect. Genomic prediction has been demonstrated to increase genetic gain and efficiency in plant breeding programs beyond marker-assisted selection and phenotypic selection. The objective of this study was to evaluate the impact of allelic origin, marker density, training population size, and cross-validation schemes on the accuracy of genomic prediction models in an interspecific soybean nested association mapping (NAM) panel. Three cross-validation schemes were used: (a) Within-Family (WF): training population and predictions are made exclusively within each family; (b) Across All families (AF): all the individuals from the three families were randomly assigned to either the training or validation set; (c) Leave one Family out (LFO): each family is predicted using a training set that contains the other two families. Predictive abilities increased with training population size up to 350 individuals, but no significant gains were noted beyond 250 individuals in the training population. The number of markers had a limited impact on the observed predictive ability across traits; increasing markers used in the model above 1000 revealed no significant increases in prediction accuracy. Predictive abilities for AF were not significantly different from the WF method, and predictive abilities across populations for the WF method had a range of 0.58 to 0.70 for maturity, protein, meal, and oil. Our results also showed encouraging prediction accuracies for grain yield (0.58-0.69) using the WF method. Partitioning genomic prediction between G. max and G. soja alleles revealed useful information to select material with a larger allele contribution from both parents and could accelerate allele introgression from exotic germplasm into the elite soybean gene pool.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01203-6.

Validated MAGIC and GWAS population mapping reveals the link between vitamin E content and natural variation in chorismate metabolism in tomato

Tocochromanols constitute the different forms of vitamin E (VTE), essential components of the human diet, and display a high membrane protectant activity. By combining interval mapping and genome-wide association studies (GWAS), we unveiled the genetic determinants of tocochromanol accumulation in tomato (Solanum lycopersicum) fruits. To enhance the nutritional value of this highly consumed vegetable, we dissected the natural intraspecific variability of tocochromanols in tomato fruits and genetically engineered their biosynthetic pathway. These analyses allowed the identification of a total of 25 quantitative trait loci interspersed across the genome pinpointing the chorismate-tyrosine pathway as a regulatory hub controlling the supply of the aromatic head group for tocochromanol biosynthesis. To validate the link between the chorismate-tyrosine pathway and VTE, we engineered tomato plants to bypass the pathway at the arogenate branch point. Transgenic tomatoes showed moderate increments in tocopherols (up to approximately 20%) and a massive accumulation of tocotrienols (up to approximately 3400%). Gene expression analyses of these plants reveal a trade-off between VTE and natural variation in chorismate metabolism explained by transcriptional reprogramming of specific structural genes of the pathway. By restoring the accumulation of alpha-tocotrienols (α-t3) in fruits, the plants produced here are of high pharmacological and nutritional interest.

Genome-Wide Association Study (GWAS) to Identify Salt-Tolerance QTLs Carrying Novel Candidate Genes in Rice During Early Vegetative Stage

BACKGROUND

Rice is considered as a salt-sensitive plant, particularly at early vegetative stage, and its production is suffered from salinity due to expansion of salt affected land in areas under cultivation. Hence, significant increase of rice productivity on salinized lands is really necessary. Today genome-wide association study (GWAS) is a method of choice for fine mapping of QTLs involved in plant responses to abiotic stresses including salinity stress at early vegetative stage. In this study using > 33,000 SNP markers we identified rice genomic regions associated to early stage salinity tolerance. Eight salinity-related traits including shoot length (SL), root length (RL), root dry weight (RDW), root fresh weight (RFW), shoot fresh weight (SFW), shoot dry weight (SDW), relative water content (RWC) and TW, and 4 derived traits including SL-R, RL-R, RDW-R and RFW-R in a diverse panel of rice were evaluated under salinity (100 mM NaCl) and normal conditions in growth chamber. Genome-wide association study (GWAS) was applied based on MLM(+Q + K) model.

RESULTS

Under stress conditions 151 trait-marker associations were identified that were scattered on 10 chromosomes of rice that arranged in 29 genomic regions. A genomic region on chromosome 1 (11.26 Mbp) was identified which co-located with a known QTL region SalTol1 for salinity tolerance at vegetative stage. A candidate gene (Os01g0304100) was identified in this region which encodes a cation chloride cotransporter. Furthermore, on this chromosome two other candidate genes, Os01g0624700 (24.95 Mbp) and Os01g0812000 (34.51 Mbp), were identified that encode a WRKY transcription factor (WRKY 12) and a transcriptional activator of gibberellin-dependent alpha-amylase expression (GAMyb), respectively. Also, a narrow interval on the same chromosome (40.79-42.98 Mbp) carries 12 candidate genes, some of them were not so far reported for salinity tolerance at seedling stage. Two of more interesting genes are Os01g0966000 and Os01g0963000, encoding a plasma membrane (PM) H⁺-ATPase and a peroxidase BP1 protein. A candidate gene was identified on chromosome 2 (Os02g0730300 at 30.4 Mbp) encoding a high affinity K⁺ transporter (HAK). On chromosome 6 a DnaJ-encoding gene and pseudouridine synthase gene were identified. Two novel genes on chromosome 8 including the ABI/VP1 transcription factor and retinoblastoma-related protein (RBR), and 3 novel genes on chromosome 11 including a Lox, F-box and Na⁺/H⁺ antiporter, were also identified.

CONCLUSION

Known or novel candidate genes in this research were identified that can be used for improvement of salinity tolerance in molecular breeding programmes of rice. Further study and identification of effective genes on salinity tolerance by the use of candidate gene-association analysis can help to precisely uncover the mechanisms of salinity tolerance at molecular level. A time dependent relationship between salt tolerance and expression level of candidate genes could be recognized.

Identifying Anti-Oxidant Biosynthesis Genes in Pearl Millet [Pennisetum glaucum (L.) R. Br.] Using Genome-Wide Association Analysis

Pearl millet [Pennisetum glaucum (L.) R Br.] is an important staple food crop in the semi-arid tropics of Asia and Africa. It is a cereal grain that has the prospect to be used as a substitute for wheat flour for celiac patients. It is an important antioxidant food resource present with a wide range of phenolic compounds that are good sources of natural antioxidants. The present study aimed to identify the total antioxidant content of pearl millet flour and apply it to evaluate the antioxidant activity of its 222 genotypes drawn randomly from the pearl millet inbred germplasm association panel (PMiGAP), a world diversity panel of this crop. The total phenolic content (TPC) significantly correlated with DPPH (1,1-diphenyl-2-picrylhydrazyl) radical scavenging activity (% inhibition), which ranged from 2.32 to 112.45% and ferric-reducing antioxidant power (FRAP) activity ranging from 21.68 to 179.66 (mg ascorbic acid eq./100 g). Genome-wide association studies (GWAS) were conducted using 222 diverse accessions and 67 K SNPs distributed across all the seven pearl millet chromosomes. Approximately, 218 SNPs were found to be strongly associated with DPPH and FRAP activity at high confidence [-log (p) > 3.0-7.4]. Furthermore, flanking regions of significantly associated SNPs were explored for candidate gene harvesting. This identified 18 candidate genes related to antioxidant pathway genes (flavanone 7-O-beta-glycosyltransferase, GDSL esterase/lipase, glutathione S-transferase) residing within or near the association signal that can be selected for further functional characterization. Patterns of genetic variability and the associated genes reported in this study are useful findings, which would need further validation before their utilization in molecular breeding for high antioxidant-containing pearl millet cultivars.

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass

Long-term climate change and periodic environmental extremes threaten food and fuel security¹ and global crop productivity^2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience⁵, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation⁶-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

Genome-wide association studies in tropical maize germplasm reveal novel and known genomic regions for resistance to Northern corn leaf blight

Northern Corn Leaf Blight (NCLB) caused by Setosphaeria turcica, is one of the most important diseases of maize world-wide, and one of the major reasons behind yield losses in maize crop in Asia. In the present investigation, a high-resolution genome wide association study (GWAS) was conducted for NCLB resistance in three association mapping panels, predominantly consisting of tropical lines adapted to different agro-ecologies. These panels were phenotyped for disease severity across three locations with high disease prevalence in India. High density SNPs from Genotyping-by-sequencing were used in GWAS, after controlling for population structure and kinship matrices, based on single locus mixed linear model (MLM). Twenty-two SNPs were identified, that revealed a significant association with NCLB in the three mapping panels. Haplotype regression analysis revealed association of 17 significant haplotypes at FDR ≤ 0.05, with two common haplotypes across three maize panels. Several of the significantly associated SNPs/haplotypes were found to be co-located in chromosomal bins previously reported for major genes like Ht2, Ht3 and Htn1 and QTL for NCLB resistance and multiple foliar disease resistance. Phenotypic variance explained by these significant SNPs/haplotypes ranged from low to moderate, suggesting a breeding strategy of combining multiple resistance alleles towards resistance for NCLB.

Genetic Dissection of Resistance to Gray Leaf Spot by Combining Genome-Wide Association, Linkage Mapping, and Genomic Prediction in Tropical Maize Germplasm

Gray leaf spot (GLS) is one of the major maize foliar diseases in sub-Saharan Africa. Resistance to GLS is controlled by multiple genes with additive effect and is influenced by both genotype and environment. The objectives of the study were to dissect the genetic architecture of GLS resistance through linkage mapping and genome-wide association study (GWAS) and assessing the potential of genomic prediction (GP). We used both biparental populations and an association mapping panel of 410 diverse tropical/subtropical inbred lines that were genotyped using genotype by sequencing. Phenotypic evaluation in two to four environments revealed significant genotypic variation and moderate to high heritability estimates ranging from 0.43 to 0.69. GLS was negatively and significantly correlated with grain yield, anthesis date, and plant height. Linkage mapping in five populations revealed 22 quantitative trait loci (QTLs) for GLS resistance. A QTL on chromosome 7 (qGLS7-105) is a major-effect QTL that explained 28.2% of phenotypic variance. Together, all the detected QTLs explained 10.50, 49.70, 23.67, 18.05, and 28.71% of phenotypic variance in doubled haploid (DH) populations 1, 2, 3, and F3 populations 4 and 5, respectively. Joint linkage association mapping across three DH populations detected 14 QTLs that individually explained 0.10-15.7% of phenotypic variance. GWAS revealed 10 significantly (p < 9.5 × 10-6) associated SNPs distributed on chromosomes 1, 2, 6, 7, and 8, which individually explained 6-8% of phenotypic variance. A set of nine candidate genes co-located or in physical proximity to the significant SNPs with roles in plant defense against pathogens were identified. GP revealed low to moderate prediction correlations of 0.39, 0.37, 0.56, 0.30, 0.29, and 0.38 for within IMAS association panel, DH pop1, DH pop2, DH pop3, F3 pop4, and F3 po5, respectively, and accuracy was increased substantially to 0.84 for prediction across three DH populations. When the diversity panel was used as training set to predict the accuracy of GLS resistance in biparental population, there was 20-50% reduction compared to prediction within populations. Overall, the study revealed that resistance to GLS is quantitative in nature and is controlled by many loci with a few major and many minor effects. The SNPs/QTLs identified by GWAS and linkage mapping can be potential targets in improving GLS resistance in breeding programs, while GP further consolidates the development of high GLS-resistant lines by incorporating most of the major- and minor-effect genes.

Genetic dissection of bread wheat diversity and identification of adaptive loci in response to elevated tropospheric ozone

Rising tropospheric ozone affects the performance of important cereal crops thus threatening global food security. In this study, genetic variation of wheat regarding its physiological and yield responses to ozone was explored by exposing a diversity panel of 150 wheat genotypes to elevated ozone and control conditions throughout the growing season. Differential responses to ozone were observed for foliar symptom formation quantified as leaf bronzing score (LBS), vegetation indices and yield components. Vegetation indices representing the carotenoid to chlorophyll pigment ratio (such as Lic2) were particularly ozone-responsive and were thus considered suitable for the non-invasive diagnosing of ozone stress. Genetic variation in ozone-responsive traits was dissected by a genome-wide association study (GWAS). Significant marker-trait associations were identified for LBS on chromosome 5A and for vegetation indices (NDVI and Lic2) on chromosomes 6B and 6D. Analysis of linkage disequilibrium (LD) in these chromosomal regions revealed distinct LD blocks containing genes with a putative function in plant redox biology such as cytochrome P450 proteins and peroxidases. This study gives novel insight into the natural genetic variation in wheat ozone response, and lays the foundation for the molecular breeding of tolerant wheat varieties.

Population structure, diversifying selection, and local adaptation in Pinus patula

PREMISE

Climate change is predicted to affect natural and plantation forests. The responses of conifers to overcome changing environments will depend on their adaptation to local conditions; however, intraspecific adaptive genetic variation is unknown for most gymnosperms. Studying genetic diversity associated with phenotypic variability along environmental gradients will enhance our understanding of adaptation and may reveal genetic pools important for conservation and management.

METHODS

We used target enrichment and genome skimming to obtain single nucleotide polymorphisms (SNPs) from 61 individuals of Pinus patula, a pine tree native to Mexico widely used in plantation forestry. We investigated the adaptive genetic variation of two varieties with morphological and distributional differences potentially related to genetic and adaptive divergence.

RESULTS

Population structure and haplotype network analyses revealed that genetic diversity between P. patula var. patula and P. patula var. longipedunculata was structured, even within populations of P. patula var. longipedunculata. We observed high genetic diversity, low inbreeding rate, and rapid linkage disequilibrium (LD) decay in the varieties. Based on outlier tests, loci showing signatures of natural selection were detected in geographically distant P. patula var. longipedunculata populations. For both varieties, we found significant correlations between climate-related environmental variation and SNP diversity at loci involved in abiotic stress, cell transport, defense, and cell wall biogenesis, pointing to local adaptation.

CONCLUSIONS

Overall, significant intraspecific adaptive genetic variation in P. patula was detected, highlighting the presence of different genetic pools and signs of local adaptation that should be considered in forestry and conservation.