Multi-trait multi-environment quantitative trait loci mapping for a sugarcane commercial cross provides insights on the inheritance of important traits

Genetic loci associated with sorghum drought tolerance in multiple environments and their sensitivity to environmental covariables

KEY MESSAGE

Climate change can limit yields of naturally resilient crops, like sorghum, challenging global food security. Agriculture under an erratic climate requires tapping into a reservoir of flexible adaptive loci that can lead to lasting yield stability under multiple abiotic stress conditions. Domesticated in the hot and dry regions of Africa, sorghum is considered a harsh crop, which is adapted to important stress factors closely related to climate change. To investigate the genetic basis of drought stress adaptation in sorghum, we used a multi-environment multi-locus genome-wide association study (MEML-GWAS) in a subset of a diverse sorghum association panel (SAP) phenotyped for performance both under well-watered and water stress conditions. We selected environments in Brazil that foreshadow agriculture where both drought and temperature stresses coincide as in many tropical agricultural frontiers. Drought reduced average grain yield (Gy) by up to 50% and also affected flowering time (Ft) and plant height (Ph). We found 15 markers associated with Gy on all sorghum chromosomes except for chromosomes 7 and 9, in addition to loci associated with phenology traits. Loci associated with Gy strongly interacted with the environment in a complex way, while loci associated with phenology traits were less affected by G × E. Studying environmental covariables potentially underpinning G × E, increases in relative humidity and evapotranspiration favored and disfavored grain yield, respectively. High temperatures influenced G × E and reduced sorghum yields, with a ~ 100 kg ha-1 average decrease in grain yield for each unit increase in maximum temperature between 29 and 38 °C. Extreme G × E for sorghum stress resilience poses an additional challenge to breed crops for moving, erratic weather conditions.

Validation of candidate gene-based EST-SSR markers for sugar yield in sugarcane

Sugarcane (Saccharum spp.) is a widely cultivated crop that fulfils approximately 75% of the sucrose demand worldwide. Owing to its polyploidy and complex genetic nature, it is difficult to identify and map genes related to complex traits, such as sucrose content. However, association mapping is one of the alternatives for identifying genes or markers for marker-assisted selection. In the present study, EST-SSR primers were obtained from in silico studies. The functionality of each primer was tested using Blast2Go software, and 30 EST-SSR primers related to sugar content were selected. These markers were validated using association analysis. A total of 70 F1 diverse genotypes for sugar content were phenotypes with two check lines. All parameters related to sugar content were recorded. The results showed a significant variation between the genotypes for sugar yield traits such as Brix value, purity, and sucrose content, etc. Correlation studies revealed that the Brix%, sucrose content, and sucrose recovery were significantly correlated. An association analysis was performed using mixed linear model to avoid false positive associations. The association analysis revealed that the SEM 407 marker was significantly associated with Brix% and sucrose content. The SEM 407 primers are putatively related to diphosphate-fructose-6-phosphate 1-phosphotransferase which is associated with Brix% and sucrose content. This functional marker can be used for marker-assisted selection for sugar yield traits in sugarcane that could accelerate the sugarcane breeding program.

Sampling Variation of RAD-Seq Data from Diploid and Tetraploid Potato (Solanum tuberosum L.)

The new sequencing technology enables identification of genome-wide sequence-based variants at a population level and a competitively low cost. The sequence variant-based molecular markers have motivated enormous interest in population and quantitative genetic analyses. Generation of the sequence data involves a sophisticated experimental process embedded with rich non-biological variation. Statistically, the sequencing process indeed involves sampling DNA fragments from an individual sequence. Adequate knowledge of sampling variation of the sequence data generation is one of the key statistical properties for any downstream analysis of the data and for implementing statistically appropriate methods. This paper reports a thorough investigation on modeling the sampling variation of the sequence data from the optimized RAD-seq (Restriction sit associated DNA sequencing) experiments with two parents and their offspring of diploid and autotetraploid potato (Solanum tuberosum L.). The analysis shows significant dispersion in sampling variation of the sequence data over that expected under multinomial distribution as widely assumed in the literature and provides statistical methods for modeling the variation and calculating the model parameters, which may be easily implemented in real sequence datasets. The optimized design of RAD-seq experiments enabled effective control of presentation of undesirable chloroplast DNA and RNA genes in the sequence data generated.

Rhythms of Transcription in Field-Grown Sugarcane Are Highly Organ Specific

Circadian clocks improve plant fitness in a rhythmic environment. As each cell has its own circadian clock, we hypothesized that sets of cells with different functions would have distinct rhythmic behaviour. To test this, we investigated whether different organs in field-grown sugarcane follow the same rhythms in transcription. We assayed the transcriptomes of three organs during a day: leaf, a source organ; internodes 1 and 2, sink organs focused on cell division and elongation; and internode 5, a sink organ focused on sucrose storage. The leaf had twice as many rhythmic transcripts (>68%) as internodes, and the rhythmic transcriptomes of the internodes were more like each other than to those of the leaves. Among the transcripts expressed in all organs, only 7.4% showed the same rhythmic pattern. Surprisingly, the central oscillators of these organs - the networks that generate circadian rhythms - had similar dynamics, albeit with different amplitudes. The differences in rhythmic transcriptomes probably arise from amplitude differences in tissue-specific circadian clocks and different sensitivities to environmental cues, highlighted by the sampling under field conditions. The vast differences suggest that we must study tissue-specific circadian clocks in order to understand how the circadian clock increases the fitness of the whole plant.

Development of an Axiom Sugarcane100K SNP array for genetic map construction and QTL identification

An Axiom Sugarcane100K SNP array has been designed and successfully utilized to construct the sugarcane genetic map and to identify the QTLs associated with SCYLV resistance. To accelerate genetic studies in sugarcane, an Axiom Sugarcane100K single-nucleotide polymorphism (SNP) array was designed and customized in this study. Target enrichment sequencing 300 sugarcane accessions selected from the world collection of sugarcane and related grass species yielded more than four million SNPs, from which a total of 31,449 single-dose (SD) SNPs and 68,648 low-dosage (33,277 SD and 35,371 double dose) SNPs from two datasets, respectively, were selected and tiled on Affymetrix Axiom SNP array. Most of selected SNPs (91.77%) were located within genic regions (12,935 genes), with an average of 7.1 SNPs/gene according to sorghum gene models. This array was used to genotype 469 sugarcane clones, including one F₁ population derived from the cross between Green German and IND81-146, one selfing population derived from CP80-1827, and 11 diverse sugarcane accessions as controls. Results of genotyping revealed a high polymorphic SNP rate (77.04%) among the 469 samples. Three linkage maps were constructed by using SD SNP markers, including a genetic map for Green German with 3482 SD SNP markers spanning 3336 cM, a map for IND81-146 with 1513 SD SNP markers spanning 2615 cM, and a map for CP80-1827 with 536 SD SNP markers spanning 3651 cM. Quantitative trait loci (QTL) analysis identified 18 QTLs controlling Sugarcane yellow leaf virus resistance segregating in the two mapping populations, harboring 27 disease-resistant genes. This study demonstrated the successful development and utilization of a SNP array as an efficient genetic tool for high-throughput genotyping in highly polyploid sugarcane.

A genome-wide association study identified loci for yield component traits in sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection.

Mapping and confirmation of loci for salt tolerance in a novel soybean germplasm, Fiskeby III

KEY MESSAGE

The confirmation of a major locus associated with salt tolerance and mapping of a new locus, which could be beneficial for improving salt tolerance in soybean. Breeding soybean for tolerance to high salt conditions is important in some regions of the USA and world. Soybean cultivar Fiskeby III (PI 438471) in maturity group 000 has been reported to be highly tolerant to multiple abiotic stress conditions, including salinity. In this study, a mapping population of 132 F2 families derived from a cross of cultivar Williams 82 (PI 518671, moderately salt sensitive) and Fiskeby III (salt tolerant) was analyzed to map salt tolerance genes. The evaluation for salt tolerance was performed by analyzing leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC) after treatment with 120 mM NaCl under greenhouse conditions. Genotypic data for the F2 population were obtained using the SoySNP6K Illumina Infinium BeadChip assay. A major allele from Fiskeby III was significantly associated with LSS, CCR, LSC, and LCC on chromosome (Chr.) 03 with LOD scores of 19.1, 11.0, 7.7 and 25.6, respectively. In addition, a second locus associated with salt tolerance for LSC was detected and mapped on Chr. 13 with an LOD score of 4.6 and an R 2 of 0.115. Three gene-based polymorphic molecular markers (Salt-20, Salt14056 and Salt11655) on Chr.03 showed a strong predictive association with phenotypic salt tolerance in the present mapping population. These molecular markers will be useful for marker-assisted selection to improve salt tolerance in soybean.

GBS-based single dosage markers for linkage and QTL mapping allow gene mining for yield-related traits in sugarcane

BACKGROUND

Sugarcane (Saccharum spp.) is predominantly an autopolyploid plant with a variable ploidy level, frequent aneuploidy and a large genome that hampers investigation of its organization. Genetic architecture studies are important for identifying genomic regions associated with traits of interest. However, due to the genetic complexity of sugarcane, the practical applications of genomic tools have been notably delayed in this crop, in contrast to other crops that have already advanced to marker-assisted selection (MAS) and genomic selection. High-throughput next-generation sequencing (NGS) technologies have opened new opportunities for discovering molecular markers, especially single nucleotide polymorphisms (SNPs) and insertion-deletion (indels), at the genome-wide level. The objectives of this study were to (i) establish a pipeline for identifying variants from genotyping-by-sequencing (GBS) data in sugarcane, (ii) construct an integrated genetic map with GBS-based markers plus target region amplification polymorphisms and microsatellites, (iii) detect QTLs related to yield component traits, and (iv) perform annotation of the sequences that originated the associated markers with mapped QTLs to search putative candidate genes.

RESULTS

We used four pseudo-references to align the GBS reads. Depending on the reference, from 3,433 to 15,906 high-quality markers were discovered, and half of them segregated as single-dose markers (SDMs) on average. In addition to 7,049 non-redundant SDMs from GBS, 629 gel-based markers were used in a subsequent linkage analysis. Of 7,678 SDMs, 993 were mapped. These markers were distributed throughout 223 linkage groups, which were clustered in 18 homo(eo)logous groups (HGs), with a cumulative map length of 3,682.04 cM and an average marker density of 3.70 cM. We performed QTL mapping of four traits and found seven QTLs. Our results suggest the presence of a stable QTL across locations. Furthermore, QTLs to soluble solid content (BRIX) and fiber content (FIB) traits had markers linked to putative candidate genes.

CONCLUSIONS

This study is the first to report the use of GBS for large-scale variant discovery and genotyping of a mapping population in sugarcane, providing several insights regarding the use of NGS data in a polyploid, non-model species. The use of GBS generated a large number of markers and still enabled ploidy and allelic dosage estimation. Moreover, we were able to identify seven QTLs, two of which had great potential for validation and future use for molecular breeding in sugarcane.