Genome-wide association study (GWAS) of leaf cuticular wax components in Camelina sativa identifies genetic loci related to intracellular wax transport.

BACKGROUND: It is important to explore renewable alternatives (e.g. biofuels) that can produce energy sources to help reduce reliance on fossil oils, and reduce greenhouse gases and waste solids resulted from fossil oils consumption. Camelina sativa is an oilseed crop which has received increasing attention due to its short life cycle, broader adaptation regions, high oil content, high level of omega-3 unsaturated fatty acids, and low-input requirements in agriculture practices. To expand its Camelina production areas into arid regions, there is a need to breed for new drought-tolerant cultivars. Leaf cuticular wax is known to facilitate plant development and growth under water-limited conditions. Dissecting the genetic loci underlying leaf cuticular waxes is important to breed for cultivars with improved drought tolerance. RESULTS: Here we combined phenotypic data and single nucleotide polymorphism (SNP) data from a spring C. sativa diversity panel using genotyping-by-sequencing (GBS) technology, to perform a large-scale genome-wide association study (GWAS) on leaf wax compositions. A total of 42 SNP markers were significantly associated with 15 leaf wax traits including major wax components such as total primary alcohols, total alkanes, and total wax esters as well as their constituents. The vast majority of significant SNPs were associated with long-chain carbon monomers (carbon chain length longer than C28), indicating the important effects of long-chain carbon monomers on leaf total wax biosynthesis. These SNP markers are located on genes directly or indirectly related to wax biosynthesis such as maintaining endoplasmic reticulum (ER) morphology and enabling normal wax secretion from ER to plasma membrane or Golgi network-mediated transport. CONCLUSIONS: These loci could potentially serve as candidates for the genetic control involved in intracellular wax transport that might directly or indirectly facilitate leaf wax accumulation in C. sativa and can be used in future marker-assisted selection (MAS) to breed for the cultivars with high wax content to improve drought tolerance.

Transcriptomic changes due to water deficit define a general soybean response and accession-specific pathways for drought avoidance.

BACKGROUND: Among abiotic stresses, drought is the most common reducer of crop yields. The slow-wilting soybean genotype PI 416937 is somewhat robust to water deficit and has been used previously to map the trait in a bi-parental population. Since drought stress response is a complex biological process, whole genome transcriptome analysis was performed to obtain a deeper understanding of the drought response in soybean. RESULTS: Contrasting data from PI 416937 and the cultivar 'Benning', we developed a classification system to identify genes that were either responding to water-deficit in both genotypes or that had a genotype x environment (GxE) response. In spite of very different wilting phenotypes, 90% of classifiable genes had either constant expression in both genotypes (33%) or very similar response profiles (E genes, 57%). By further classifying E genes based on expression profiles, we were able to discern the functional specificity of transcriptional responses at particular stages of water-deficit, noting both the well-known reduction in photosynthesis genes as well as the less understood up-regulation of the protein transport pathway. Two percent of classifiable genes had a well-defined GxE response, many of which are located within slow-wilting QTLs. We consider these strong candidates for possible causal genes underlying PI 416937's unique drought avoidance strategy. CONCLUSIONS: There is a general and functionally significant transcriptional response to water deficit that involves not only known pathways, such as down-regulation of photosynthesis, but also up-regulation of protein transport and chromatin remodeling. Genes that show a genotypic difference are more likely to show an environmental response than genes that are constant between genotypes. In this study, at least five genes that clearly exhibited a genotype x environment response fell within known QTL and are very good candidates for further research into slow-wilting.

Confirmation of delayed canopy wilting QTLs from multiple soybean mapping populations.

KEY MESSAGE: QTLs for delayed canopy wilting from five soybean populations were projected onto the consensus map to identify eight QTL clusters that had QTLs from at least two independent populations. Quantitative trait loci (QTLs) for canopy wilting were identified in five recombinant inbred line (RIL) populations, 93705 KS4895 × Jackson, 08705 KS4895 × Jackson, KS4895 × PI 424140, A5959 × PI 416937, and Benning × PI 416937 in a total of 15 site-years. For most environments, heritability of canopy wilting ranged from 0.65 to 0.85 but was somewhat lower when averaged over environments. Putative QTLs were identified with composite interval mapping and/or multiple interval mapping methods in each population and positioned on the consensus map along with their 95% confidence intervals (CIs). We initially found nine QTL clusters with overlapping CIs on Gm02, Gm05, Gm11, Gm14, Gm17, and Gm19 identified from at least two different populations, but a simulation study indicated that the QTLs on Gm14 could be false positives. A QTL on Gm08 in the 93705 KS4895 × Jackson population co-segregated with a QTL for wilting published previously in a Kefeng1 × Nannong 1138-2 population, indicating that this may be an additional QTL cluster. Excluding the QTL cluster on Gm14, results of the simulation study indicated that the eight remaining QTL clusters and the QTL on Gm08 appeared to be authentic QTLs. QTL × year interactions indicated that QTLs were stable over years except for major QTLs on Gm11 and Gm19. The stability of QTLs located on seven clusters indicates that they are possible candidates for use in marker-assisted selection.

Fine mapping and identification of candidate genes controlling the resistance to southern root-knot nematode in PI 96354.

Meloidogyne incognita (Kofoid and White) Chitwood (Mi) is the most economically damaging species of the root-knot nematode to soybean and other crops in the southern USA. PI 96354 was identified to carry a high level of resistance to galling and Mi egg production. Two Quantitative Trait Locus (QTLs) were found to condition the resistance in PI 96354 including a major QTL and a minor QTL on chromosome 10 and chromosome 18, respectively. To fine map the major QTL on chromosome 10, F5:6 recombinant inbred lines from the cross between PI 96354 and susceptible genotype Bossier were genotyped with Simple Sequence Repeats (SSR) markers to identify recombinational events. Analysis of lines carrying key recombination events placed the Mi-resistant allele on chromosome 10 to a 235-kb region of the 'Williams 82' genome sequence with 30 annotated genes. Candidate gene analysis identified four genes with cell wall modification function that have several mutations in promoter, exon, 5', and 3'UTR regions. qPCR analysis showed significant difference in expression levels of these four genes in Bossier compared to PI 96354 in the presence of Mi. Thirty Mi-resistant soybean lines were found to have same SNPs in these 4 candidate genes as PI 96354 while 12 Mi-susceptible lines possess the 'Bossier' genotype. The mutant SNPs were used to develop KASP assays to detect the resistant allele on chromosome 10. The four candidate genes identified in this study can be used in further studies to investigate the role of cell wall modification genes in conferring Mi resistance in PI 96354.

Phenotypic Diversity in Leaf Cuticular Waxes in Brassica carinata Accessions.

Brassica carinata has received considerable attention as a renewable biofuel crop for semi-arid zones due to its high oil content and polyunsaturated fatty acids contents. It is important to develop new drought-resistant cultivars of B. carinata production to expand its areas into more arid regions. The accumulation of leaf cuticular wax on plant surfaces is one mechanism that reduces non-stomatal water loss, thus increasing drought resistance in plants. To explore phenotypic variations in cuticular wax in B. carinata, leaf waxes were extracted and quantified from a diversity panel consisting of 315 accessions. The results indicate that the accessions have a wide range of total leaf wax content (289-1356 µg dm-2), wax classes, and their components. The C29 and C31 homologues of alkanes, C29 ketone homologue, C29 secondary alcohol, and C30 aldehyde were the most abundant leaf waxes extracted from B. carinata accessions. The high heritability values of these waxes point to the positive selection for high wax content during early generations of future B. carinata breeding programs. Positive correlation coefficients, combined with the effects of these waxes on leaf wax content accumulation, suggest that modifying specific wax content could increase the total wax content and enhance cuticle composition. The identified leaf wax content and compositions in B. carinata will lead to the future discovery of wax biosynthetic pathways, the dissection of its genetic regulatory networks, the identification of candidate genes controlling production of these waxes, and thus, develop and release new B. carinata drought-tolerant cultivars.

Identification of genomic regions associated with the plasticity of carbon 13 ratio in soybean.

Improving water use efficiency (WUE) for soybean [Glycine max (L.) Merr.] through selection for high carbon isotope (C13) ratio may increase drought tolerance, but increased WUE may limit growth in productive environments. An ideal genotype would be plastic for C13 ratio; that is, be able to alter C13 ratio in response to the environment. Our objective was to identify genomic regions associated with C13 ratio plasticity, C13 ratio stability, and overall C13 ratio in two panels of diverse Maturity Group IV soybean accessions. A second objective was to identify accessions that differed in their C13 ratio plasticity. Panel 1 (205 accessions) was evaluated in seven irrigated and four drought environments, and Panel 2 (373 accessions) was evaluated in four environments. Plasticity was quantified as the slope from regressing C13 ratio of individual genotypes against an environmental index calculated based on the mean within and across environments. The regression intercept was considered a measure of C13 ratio over all environments, and the root mean square error was considered a measure of stability. Combined over both panels, genome-wide association mapping (GWAM) identified 19 single nucleotide polymorphisms (SNPs) for plasticity, 39 SNPs for C13 ratio, and 16 SNPs for stability. Among these SNPs, 71 candidate genes had annotations associated with transpiration or water conservation and transport, root development, root hair elongation, and stomatal complex morphogenesis. The genomic regions associated with plasticity and stability identified in the current study will be a useful resource for implementing genomic selection for improving drought tolerance in soybean.

Mapping of quantitative trait loci for canopy-wilting trait in soybean (Glycine max L. Merr).

Drought stress adversely affects [Glycine max (L.) Merr] soybean at most developmental stages, which collectively results in yield reduction. Little information is available on relative contribution and chromosomal locations of quantitative trait loci (QTL) conditioning drought tolerance in soybean. A Japanese germplasm accession, PI 416937, was found to possess drought resistance. Under moisture-deficit conditions, PI 416937 wilted more slowly in the field than elite cultivars and has been used as a parent in breeding programs to improve soybean productivity. A recombinant inbred line (RIL) population was derived from a cross between PI 416937 and Benning, and the population was phenotyped for canopy wilting under rain-fed field conditions in five distinct environments to identify the QTL associated with the canopy-wilting trait. In a combined analysis over environments, seven QTL that explained 75 % of the variation in canopy-wilting trait were identified on different chromosomes, implying the complexity of this trait. Five QTL inherited their positive alleles from PI 416937. Surprisingly, the other two QTL inherited their positive alleles from Benning. These putative QTL were co-localized with other QTL previously identified as related to plant abiotic stresses in soybean, suggesting that canopy-wilting QTL may be associated with additional morpho-physiological traits in soybean. A locus on chromosome 12 (Gm12) from PI 416937 was detected in the combined analysis as well as in each individual environment, and explained 27 % of the variation in canopy-wilting. QTL identified in PI 416937 could provide an efficient means to augment field-oriented development of drought-tolerant soybean cultivars.

Identification of QTL for increased fibrous roots in soybean.

Drought stress adversely affects soybean at various developmental stages, which collectively results in yield reduction. Unpredictable rainfall has been reported to contribute about 36% to variation of yield difference between the rain-fed and irrigated fields. Among the drought resistance mechanisms, drought avoidance in genotypes with fibrous roots was recognized to be associated with drought resistance in soybean. Plant introduction PI416937 was shown to possess fibrous roots and has been used as a parent in breeding programs to improve soybean productivity. Little information is available on relative contribution and chromosomal location of quantitative trait loci (QTL) conditioning fibrous roots in soybean. To identify the genomic locations and genetic bases of this trait, a recombinant inbred line population was derived from a cross between PI416937 and 'Benning'. To detect associated QTLs, phenotypic data were collected and analyzed for 2 years under rain-fed field conditions. The selective genotyping approach was used to reduce the costs and work associated with conducting the QTL analysis. A total of five QTLs were identified on chromosomes Gm01 (Satt383), Gm03 (Satt339), Gm04 (Sct_191), Gm08 (Satt429), and Gm20 (Sat_299), and together explained 51% of the variation in root score. Detected QTLs were co-localized with QTLs related to root morphology, suggesting that fibrous roots QTL may be associated with other morpho-physiological traits and seed yield in soybean. Genetic dissection of the fibrous roots trait at the individual marker loci will allow for marker-assisted selection to develop soybean genotypes with enhanced levels of fibrous roots.

Identification and Confirmation of Loci Associated With Canopy Wilting in Soybean Using Genome-Wide Association Mapping.

Drought causes significant soybean [Glycine max (L.) Merr.] yield losses each year in rain-fed production systems of many regions. Genetic improvement of soybean for drought tolerance is a cost-effective approach to stabilize yield under rain-fed management. The objectives of this study were to confirm previously reported soybean loci and to identify novel loci associated with canopy wilting (CW) using a panel of 200 diverse maturity group (MG) IV accessions. These 200 accessions along with six checks were planted at six site-years using an augmented incomplete block design with three replications under irrigated and rain-fed treatments. Association mapping, using 34,680 single nucleotide polymorphisms (SNPs), identified 188 significant SNPs associated with CW that likely tagged 152 loci. This includes 87 SNPs coincident with previous studies that likely tagged 68 loci and 101 novel SNPs that likely tagged 84 loci. We also determined the ability of genomic estimated breeding values (GEBVs) from previous research studies to predict CW in different genotypes and environments. A positive relationship (P ≤ 0.05;0.37 ≤ r ≤ 0.5) was found between observed CW and GEBVs. In the vicinity of 188 significant SNPs, 183 candidate genes were identified for both coincident SNPs and novel SNPs. Among these 183 candidate genes, 57 SNPs were present within genes coding for proteins with biological functions involved in plant stress responses. These genes may be directly or indirectly associated with transpiration or water conservation. The confirmed genomic regions may be an important resource for pyramiding favorable alleles and, as candidates for genomic selection, enhancing soybean drought tolerance.

Genetic architecture of novel traits in the hopi sunflower.

Following domestication, crop lineages typically undergo diversification either to adapt to disparate habitats or to fill novel agricultural roles. This process has produced the numerous varieties found in modern-day crop germplasm collections. Here, we mapped quantitative trait loci (QTLs) underlying unique traits in the Hopi sunflower, a primitive, Native American domesticate. These traits included a variety of achene (i.e., single-seeded fruit) characters as well as the extremely late flowering time of the Hopi sunflower. Composite interval mapping identified 42 QTLs underlying the 12 traits of interest. Although these QTLs were found on 10 of the 17 sunflower linkage groups, strong genetic correlations were evidenced by the clustering of QTLs across traits in certain genomic regions. The number of QTLs per trait ranged from 2 to 6, and the average QTL explained 14.7% of the variance (range: 2.5-46.9%). The apparent contribution of epistasis was minor, as has previously been observed for domestication-related traits. Unlike typical domestication-related traits in sunflower, the traits under consideration here exhibited a relatively simple genetic basis, with 2 QTL clusters being largely responsible for the unique characteristics of the Hopi sunflower. Based on the rarity of these traits in domesticated sunflower, it would appear that they evolved within the Hopi lineage following domestication. The simple genetic architecture of these traits may be a by-product of genetic constraints imposed by the genetically complex nature of domestication-related traits in sunflower, with the large number of domestication-related QTLs limiting the fraction of the genome that is available for subsequent diversification.

Genome-Wide Association Mapping of Dark Green Color Index using a Diverse Panel of Soybean Accessions.

Nitrogen (N) plays a key role in plants because it is a major component of RuBisCO and chlorophyll. Hence, N is central to both the dark and light reactions of photosynthesis. Genotypic variation in canopy greenness provides insights into the variation of N and chlorophyll concentration, photosynthesis rates, and N2 fixation in legumes. The objective of this study was to identify significant loci associated with the intensity of greenness of the soybean [Glycine max (L.) Merr.] canopy as determined by the Dark Green Color Index (DGCI). A panel of 200 maturity group IV accessions was phenotyped for canopy greenness using DGCI in three environments. Association mapping identified 45 SNPs that were significantly (P ≤ 0.0003) associated with DGCI in three environments, and 16 significant SNPs associated with DGCI averaged across all environments. These SNPs likely tagged 43 putative loci. Out of these 45 SNPs, eight were present in more than one environment. Among the identified loci, 21 were located in regions previously reported for N traits and ureide concentration. Putative loci that were coincident with previously reported genomic regions may be important resources for pyramiding favorable alleles for improved N and chlorophyll concentrations, photosynthesis rates, and N2 fixation in soybean.

A high-throughput quantification of resin and rubber contents in Parthenium argentatum using near-infrared (NIR) spectroscopy.

BACKGROUND: Guayule (Parthenium argentatum A. Gray), a plant native to semi-arid regions of northern Mexico and southern Texas in the United States, is an alternative source for natural rubber (NR). Rapid screening tools are needed to replace the current labor-intensive and cost-inefficient method for quantifying rubber and resin contents. Near-infrared (NIR) spectroscopy is a promising technique that simplifies and speeds up the quantification procedure without losing precision. In this study, two spectral instruments were used to rapidly quantify resin and rubber contents in 315 ground samples harvested from a guayule germplasm collection grown under different irrigation conditions at Maricopa, AZ. The effects of eight different pretreatment approaches on improving prediction models using partial least squares regression (PLSR) were investigated and compared. Important characteristic wavelengths that contribute to prominent absorbance peaks were identified. RESULTS: Using two different NIR devices, ASD FieldSpec®3 performed better than Polychromix Phazir™ in improving R2 and residual predicative deviation (RPD) values of PLSR models. Compared to the models based on full-range spectra (750-2500 nm), using a subset of wavelengths (1100-2400 nm) with high sensitivity to guayule rubber and resin contents could lead to better prediction accuracy. The prediction power of the models for quantifying resin content was better than rubber content. CONCLUSIONS: In summary, the calibrated PLSR models for resin and rubber contents were successfully developed for a diverse guayule germplasm collection and were applied to roughly screen samples in a low-cost and efficient way. This improved efficiency could enable breeders to rapidly screen large guayule populations to identify cultivars that are high in rubber and resin contents.

Genetic Diversity and Population Structure of a Camelina sativa Spring Panel.

There is a need to explore renewable alternatives (e.g., biofuels) that can produce energy sources to help reduce the reliance on fossil oils. In addition, the consumption of fossil oils adversely affects the environment and human health via the generation of waste water, greenhouse gases, and waste solids. Camelina sativa, originated from southeastern Europe and southwestern Asia, is being re-embraced as an industrial oilseed crop due to its high seed oil content (36-47%) and high unsaturated fatty acid composition (>90%), which are suitable for jet fuel, biodiesel, high-value lubricants and animal feed. C. sativa's agronomic advantages include short time to maturation, low water and nutrient requirements, adaptability to adverse environmental conditions and resistance to common pests and pathogens. These characteristics make it an ideal crop for sustainable agricultural systems and regions of marginal land. However, the lack of genetic and genomic resources has slowed the enhancement of this emerging oilseed crop and exploration of its full agronomic and breeding potential. Here, a core of 213 spring C. sativa accessions was collected and genotyped. The genotypic data was used to characterize genetic diversity and population structure to infer how natural selection and plant breeding may have affected the formation and differentiation within the C. sativa natural populations, and how the genetic diversity of this species can be used in future breeding efforts. A total of 6,192 high-quality single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing (GBS) technology. The average polymorphism information content (PIC) value of 0.29 indicate moderate genetic diversity for the C. sativa spring panel evaluated in this report. Population structure and principal coordinates analyses (PCoA) based on SNPs revealed two distinct subpopulations. Sub-population 1 (POP1) contains accessions that mainly originated from Germany while the majority of POP2 accessions (>75%) were collected from Eastern Europe. Analysis of molecular variance (AMOVA) identified 4% variance among and 96% variance within subpopulations, indicating a high gene exchange (or low genetic differentiation) between the two subpopulations. These findings provide important information for future allele/gene identification using genome-wide association studies (GWAS) and marker-assisted selection (MAS) to enhance genetic gain in C. sativa breeding programs.

Fine Mapping and Characterization of Candidate Genes that Control Resistance to Cercospora sojina K. Hara in Two Soybean Germplasm Accessions.

Frogeye leaf spot (FLS), caused by the fungus Cercospora sojina K. Hara, may cause a significant yield loss to soybean growers in regions with a warm and humid climate. Two soybean accessions, PI 594891 and PI 594774, were identified to carry a high level of resistance similar to that conditioned by the Rcs3 gene in 'Davis'. Previously, we reported that the resistance to FLS in these two plant introductions (PIs) was controlled by a novel gene (s) on chromosome 13 that is different from Rcs3. To fine-map the novel FLS resistance gene(s) in these two PIs, F2: 3 seeds from the crosses between PI 594891 and PI 594774, and the FLS susceptible genotype 'Blackhawk' were genotyped with SNP markers that were designed based on the SoySNP50k iSelect BeadChip data to identify recombinant events and locate candidate genes. Analysis of lines possessing key recombination events helped narrow down the FLS-resistance genomic region in PI 594891 from 3.3 Mb to a 72.6 kb region with five annotated genes. The resistance gene in PI 594774 was fine-mapped into a 540 kb region that encompasses the 72.6 kb region found in PI 594891. Sequencing five candidate genes in PI 594891 identified three genes that have several mutations in the promoter, intron, 5', and 3' UTR regions. qPCR analysis showed a difference in expression levels of these genes in both lines compared to Blackhawk in the presence of C. sojina. Based on phenotype, genotype and haplotype analysis results, these two soybean accessions might carry different resistance alleles of the same gene or two different gene(s). The identified SNPs were used to develop Kompetitive Allele Specific PCR (KASP) assays to detect the resistance alleles on chromosome 13 from the two PIs for marker-assisted selection.

An R package for SNP marker-based parent-offspring tests.

BACKGROUND: With the advancement of genotyping technologies, whole genome and high-density SNP markers have been widely used for genotyping of mapping populations and for characterization of germplasm lines in many crops. Before conducting SNP data analysis, it is necessary to check the individuals to ensure the integrity of lines for further data analysis. RESULTS: We have developed an R package to conduct a parent-offspring test of individuals which are genotyped with a fixed set of SNP markers for further genetic studies. The program uses monomorphic SNP loci between parents and their progeny genotypes to calculate the similarity between each offspring and their parents. Based on the similarity of parents and individual offspring, the users can determine the threshold level for the individuals to be included for further data analysis. We used an F5-derived soybean population of '5601T' x PI 157440 that was genotyped with 1,536 SNPs to illustrate the procedure and its application. CONCLUSIONS: The R package 'ParentOffspring' coupled with the available SNP genotyping platforms could be used to detect the possible variants in a specific cross, as well as the potential errors in sample handling and genotyping processes. It can be used in any crop which is genotyped with a fixed set of SNP markers.