Detection and validation of EST-SSR markers associated with sugar-related traits in sugarcane using linkage and association mapping

Identification and characterization of functionally relevant SSR markers in natural Dalbergia odorifera populations

BACKGROUND

Dalbergia odorifera is a rare and precious rosewood specie, which is valued for its amber tones, abstract figural patterns, and impermeability to water and insects. However, the information on genetic diversity and marker-assisted selection breeding of D. odorifera is still limited. Simple sequence repeat (SSR) markers are an ideal tool for genetic diversity analysis and marker-assisted molecular breeding for complex traits.

RESULTS

Here, we have developed SSR markers within candidate genes and used them to explore the genetic diversity among D. odorifera germplasm resources. A total of 635 SSR loci were identified. The proportions of mono-, di- and tri-nucleotide repeat motifs were 52.28%, 22.99% and 21.42%, respectively. From these, a total of 114 SSR primers were synthesized, of which 24 SSR markers displayed polymorphism (polymorphic information content (PIC) > 0.25). Subsequently, these polymorphic markers were used for the genetic diversity analysis of 106 D. odorifera individuals from 11 natural populations. According to the genetic diversity analysis of D. odorifera natural populations, the average observed heterozygosity (Ho) was 0.500, the average expected heterozygosity (He) was 0.524, and the average Shannon's information index (I) was 0.946. These indicated that the natural populations had moderate genetic diversity. AMOVA analysis showed that 5% of the total variation was within the individuals of a population, whereas 95% of the variation was among the individuals of the populations, indicating a high degree of genetic variation between populations. On the basis of their genetic structures, these populations could be divided into four groups.

CONCLUSIONS

Our study provides important experimental resources for genetic studies and assists in the program of molecular breeding of D. odorifera wood formation.

A SNP variation in the Sucrose synthase (SoSUS) gene associated with sugar-related traits in sugarcane

Background

Sugarcane (Saccharum spp.) is an economically significant crop for both the sugar and biofuel industries. Breeding sugarcane cultivars with high-performance agronomic traits is the most effective approach for meeting the rising demand for sugar and biofuels. Molecular markers associated with relevant agronomic traits could drastically reduce the time and resources required to develop new sugarcane varieties. Previous sugarcane candidate gene association analyses have found single nucleotide polymorphism (SNP) markers associated with sugar-related traits. This study aims to validate these associated SNP markers of six genes, including Lesion simulating disease 1 (LSD), Calreticulin (CALR), Sucrose synthase 1 (SUS1), DEAD-box ATP-dependent RNA helicase (RH), KANADI1 (KAN1), and Sodium/hydrogen exchanger 7 (NHX7), in a diverse population in 2-year and two-location evaluations.

Methods

After genotyping of seven targeted SNP markers was performed by PCR Allelic Competitive Extension (PACE) SNP genotyping, the association with sugar-related traits and important cane yield component traits was determined on a set of 159 sugarcane genotypes. The marker-trait relationships were validated and identified by both t-test analysis and an association analysis based on the general linear model.

Results

The mSoSUS1_SNPCh10.T/C and mSoKAN1_SNPCh7.T/C markers that were designed from the SUS1 and KAN1 genes, respectively, showed significant associations with different amounts of sugar-related traits and yield components. The mSoSUS1_SNPCh10.T/C marker was found to have more significant association with sugar-related traits, including pol, CCS, brix, fiber and sugar yield, with p values of 6.08 × 10-6 to 4.35 × 10-2, as well as some cane yield component traits with p values of 1.61 × 10-4 to 3.35 × 10-2. The significant association is consistent across four environments.

Conclusion

Sucrose synthase (SUS) is considered a crucial enzyme involved in sucrose metabolism. This marker is a high potential functional marker that may be used in sugarcane breeding programs to select superior sugarcane with good fiber and high sugar contents.

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.

Development of simple sequence repeat markers for sugarcane from data mining of expressed sequence tags

Sugarcane (Saccharum spp. hybrids) is a worldwide acclaimed important agricultural crop used primarily for sugar production and biofuel. Sugarcane's genetic complexity, aneuploidy, and extreme heterozygosity make it a challenging crop in developing improved varieties. The molecular breeding programs promise to develop nutritionally improved varieties for both direct consumption and commercial application. Therefore, to address these challenges, the development of simple sequence repeats (SSRs) has been proven to be a powerful molecular tool in sugarcane. This study involved the collection of 285216 expressed sequence tags (ESTs) from sugarcane, resulting in 23666 unigenes, including 4547 contigs. Our analysis identified 4120 unigenes containing a total of 4960 SSRs, with the most abundant repeat types being monomeric (44.33%), dimeric (13.10%), and trimeric (39.68%). We further chose 173 primers to analyze the banding pattern in 10 sugarcane accessions by PAGE analysis. Additionally, functional annotation analysis showed that 71.07%, 53.6%, and 10.3% unigenes were annotated by Uniport, GO, and KEGG, respectively. GO annotations and KEGG pathways were distributed across three functional categories: molecular (46.46%), cellular (33.94%), and biological pathways (19.6%). The cluster analysis indicated the formation of four distinct clusters among selected sugarcane accessions, with maximum genetic distance observed among the varieties. We believe that these EST-SSR markers will serve as valuable references for future genetic characterization, species identification, and breeding efforts in sugarcane.

A Large-Scale Candidate-Gene Association Mapping for Drought Tolerance and Agronomic Traits in Sugarcane

Dissection of the genetic loci controlling drought tolerance traits with a complex genetic inheritance is important for drought-tolerant sugarcane improvement. In this study, we conducted a large-scale candidate gene association study of 649 candidate genes in a sugarcane diversity panel to identify genetic variants underlying agronomic traits and drought tolerance indices evaluated in plant cane and ratoon cane under water-stressed (WS) and non-stressed (NS) environments. We identified 197 significant marker-trait associations (MTAs) in 141 candidate genes associated with 18 evaluated traits with the Bonferroni correction threshold (α = 0.05). Out of the total, 95 MTAs in 78 candidate genes and 62 MTAs in 58 candidate genes were detected under NS and WS conditions, respectively. Most MTAs were found only in specific water regimes and crop seasons. These MTAs explained 7.93-30.52% of phenotypic variation. Association mapping results revealed that 34, 59, and 104 MTAs involved physiological and molecular adaptation, phytohormone metabolism, and drought-inducible genes. They identified 19 pleiotropic genes associated with more than one trait and many genes related to drought tolerance indices. The genetic and genomic resources identified in this study will enable the combining of yield-related traits and sugar-related traits with agronomic value to optimize the yield of sugarcane cultivars grown under drought-stressed and non-stressed environments.

Identification of the Efficacy of Ex Situ Conservation of Ammopiptanthus nanus Based on Its ETS-SSR Markers

Ammopiptanthus possesses ancestral traits and, as a tertiary relict, is one of the surviving remnants of the ancient Mediterranean retreat and climate drought. It is also the only genus of super xerophytic, evergreen, broad-leaved shrubs. Ammopiptanthus nanus, one of the two species in this genus, is predominantly found in extremely arid and frigid environments, and is increasingly threatened with extinction. Study of the species' genetic diversity is thus beneficial for its survival and the efficacy of ex situ conservation efforts. Based on transcriptome data, 15 pairs of effective EST-SSR were screened to evaluate A. nanus genetic diversity. In all, 87 samples from three populations were evaluated, the results of which show that ex situ conservation in the Wuqia region needs to be supplemented. Conservation and breeding of individual A. nanus offspring should be strengthened in the future to ensure their progeny continue to exhibit high genetic diversity and variation.

QTL mapping and identification of candidate genes linked to red rot resistance in sugarcane

Sugarcane (Saccharum species hybrid) is one of the most important commercial crops cultivated worldwide for products like white sugar, bagasse, ethanol, etc. Red rot is a major sugarcane disease caused by a hemi-biotrophic fungus, Colletotrichum falcatum Went., which can potentially cause a reduction in yield up to 100%. Breeding for red rot-resistant sugarcane varieties has become cumbersome due to its complex genome and frequent generation of new pathotypes of red rot fungus. In the present study, a genetic linkage map was developed using a selfed population of a popular sugarcane variety CoS 96268. A QTL linked to red rot resistance (qREDROT) was identified, which explained 26% of the total phenotypic variation for the trait. A genotype-phenotype network analysis performed to account for epistatic interactions, identified the key markers involved in red rot resistance. The differential expression of the genes located in the genomic region between the two flanking markers of the qREDROT as well as in the vicinity of the markers identified through the genotype-phenotype network analysis in a set of contrasting genotypes for red rot infection further confirmed the mapping results. Further, the expression analysis revealed that the plant defense-related gene coding 26S protease regulatory subunit is strongly associated with the red rot resistance. The findings can help in the screening of disease resistant genotypes for developing red rot-resistant varieties of sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03481-7.

A Genome-Wide Association Study and Genomic Prediction for Fiber and Sucrose Contents in a Mapping Population of LCP 85-384 Sugarcane

Sugarcane (Saccharum spp. hybrids) is an economically important crop for both sugar and biofuel industries. Fiber and sucrose contents are the two most critical quantitative traits in sugarcane breeding that require multiple-year and multiple-location evaluations. Marker-assisted selection (MAS) could significantly reduce the time and cost of developing new sugarcane varieties. The objectives of this study were to conduct a genome-wide association study (GWAS) to identify DNA markers associated with fiber and sucrose contents and to perform genomic prediction (GP) for the two traits. Fiber and sucrose data were collected from 237 self-pollinated progenies of LCP 85-384, the most popular Louisiana sugarcane cultivar from 1999 to 2007. The GWAS was performed using 1310 polymorphic DNA marker alleles with three models of TASSEL 5, single marker regression (SMR), general linear model (GLM) and mixed linear model (MLM), and the fixed and random model circulating probability unification (FarmCPU) of R package. The results showed that 13 and 9 markers were associated with fiber and sucrose contents, respectively. The GP was performed by cross-prediction with five models, ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB) and Bayesian least absolute shrinkage and selection operator (BL). The accuracy of GP varied from 55.8% to 58.9% for fiber content and 54.6% to 57.2% for sucrose content. Upon validation, these markers can be applied in MAS and genomic selection (GS) to select superior sugarcane with good fiber and high sucrose contents.

Evaluation of a Sugarcane (Saccharum spp.) Hybrid F1 Population Phenotypic Diversity and Construction of a Rapid Sucrose Yield Estimation Model for Breeding

Sugarcane is the major sugar-producing crop worldwide, and hybrid F₁ populations are the primary populations used in breeding. Challenged by the sugarcane genome's complexity and the sucrose yield's quantitative nature, phenotypic selection is still the most commonly used approach for high-sucrose yield sugarcane breeding. In this study, a hybrid F₁ population containing 135 hybrids was constructed and evaluated for 11 traits (sucrose yield (SY) and its related traits) in a randomized complete-block design during two consecutive growing seasons. The results revealed that all the traits exhibited distinct variation, with the coefficient of variation (CV) ranging from 0.09 to 0.35, the Shannon-Wiener diversity index (H') ranging between 2.64 and 2.98, and the broad-sense heritability ranging from 0.75 to 0.84. Correlation analysis revealed complex correlations between the traits, with 30 trait pairs being significantly correlated. Eight traits, including stalk number (SN), stalk diameter (SD), internode length (IL), stalk height (SH), stalk weight (SW), Brix (B), sucrose content (SC), and yield (Y), were significantly positively correlated with sucrose yield (SY). Cluster analysis based on the 11 traits divided the 135 F₁ hybrids into three groups, with 55 hybrids in Group I, 69 hybrids in Group II, and 11 hybrids in Group III. The principal component analysis indicated that the values of the first four major components' vectors were greater than 1 and the cumulative contribution rate reached 80.93%. Based on the main component values of all samples, 24 F₁ genotypes had greater values than the high-yielding parent 'ROC22' and were selected for the next breeding stage. A rapid sucrose yield estimation equation was established using four easily measured sucrose yield-related traits through multivariable linear stepwise regression. The model was subsequently confirmed using 26 sugarcane cultivars and 24 F₁ hybrids. This study concludes that the sugarcane F₁ population holds great genetic diversity in sucrose yield-related traits. The sucrose yield estimation model, ySY=2.01xSN+8.32xSD+0.79xB+3.44xSH-47.64, can aid to breed sugarcane varieties with high sucrose yield.

Development of EST-SSR markers based on transcriptome sequencing for germplasm evaluation of 65 lilies (Lilium)

BACKGROUND

Lilium genus consists of approximately 100 species and numerous varieties, many of which are interspecific hybrids, which result in a complicated genetic background. The germplasm identification, genetic relationship analysis, and breeding of Lilium rely on exploiting genetic information among different accessions. Hence, an attempt was made to develop new EST-SSR markers and study the molecular divergence among 65 genotypes of Lilium.

METHODS AND RESULTS

A total of 5509 EST-SSRs were identified from the high-throughput sequencing database of L. 'Elodie'. After primer screening, six primer pairs with the most abundant polymorphic bands were selected from 100 primer pairs. Combined with the other 10 reported SSR primers, a total of 16 pairs detected genetic information with an average PIC value of 0.7583. The number of alleles per locus varied from four to 33, the expected heterozygosity varied from 0.3289 to 0.9231, and the observed heterozygosity varied from 0.2857 to 0.5122. Based on the phylogenic results, 22 Asiatic hybrids (A), seven Longiflorum × Asiatic hybrids (LA), as well as two native species were grouped. Eighteen Oriental hybrids (O) and nine Oriental × Trumpet (OT) hybrids, four native species, one LO, and one LL (L. pardalinum × L. longiflorum) variety were grouped.

CONCLUSIONS

Two major clusters were reported and a large number of genotypes were grouped based on UPGMA and STRUCTURE analysis methods. The PIC value as well as other parameters revealed that the EST-SSR markers selected were informative. In addition, the clustering pattern displayed better agreement with the cultivar's pedigree. The newly identified SSRs in this study provides molecular markers for germplasm characterization and genetic diversity for Lilium.

Comparative transcriptome analysis on the mangrove Acanthus ilicifolius and its two terrestrial relatives provides insights into adaptation to intertidal habitats

Acanthus is a unique genus covering both mangroves and terrestrial species, and thus is an ideal system to comparatively analyze the mechanisms of mangrove adaptation to intertidal habitats. We performed RNA sequencing of the mangrove plant Acanthus ilicifolius and its two terrestrial relatives, Acanthus leucostachyus and Acanthus mollis. A total of 91,125, 118,290, and 141,640 unigenes were obtained. Simple sequence repeats (SSR) analysis showed that A. ilicifolius had more SSRs, the highest frequency of distribution, and higher in polymorphism potential compared to the two terrestrial relatives. Phylogenetic analyses suggested a relatively recent split between A. ilicifolius and A. leucostachyus, i.e., about 16.76 million years ago (Mya), after their ancestor divergence with A. mollis (32.11 Mya), indicating that speciation of three Acanthus species occurred in the Early to Middle Miocene. Gene Ontology (GO) enrichment revealed that the unique unigenes in A. ilicifolius are predominantly related to rhythmic process, reproductive process and response to stimuli. The accelerated evolution and positive selection analyses indicated that the genus Acanthus migrated from terrestrial to intertidal habitats, where 311 pairs may be under positive selection. Functional enrichment analysis revealed that these genes associated with essential metabolism and biosynthetic pathways such as oxidative phosphorylation, plant hormone signal transduction, photosynthetic carbon fixation and arginine and proline metabolism, are related to the adaptation of A. ilicifolius to intertidal habitats, which are characterized by high salinity and hypoxia. Our results indicate the evolutionary processes and the mechanisms underlying the adaptability of Acanthus to various harsh environments from the arid terrestrial to intertidal habitats.

Morphological, agronomical, physiological and molecular characterization of a high sugar mutant of sugarcane in comparison to mother variety

Sugarcane is a significant crop plant with the capability of accumulating higher amount of sucrose. In the present study, a high sucrose content sugarcane mutant clone, GXB9, has been studied in comparison to the low sucrose mother clone B9 on morphological, agronomical and physiological level in order to scrutinize the variation because of mutation in GXB9 in field under normal environmental condition. The results showed that GXB9 has less germination, tillering rate, stalk height, leaf length, leaf width, leaf area, number of internodes, internode length and internode diameter than B9. Qualitative traits of leaf and stalk displayed significant variation between GXB9 and B9. Endogenous hormones quantity was also showed variation between the two clones. The relative SPAD reading and chlorophyll a, b concentrations also showed variation between GXB9 and B9. The photosynthetic parameter analysis indicated that the GXB9 has significantly higher net photosynthetic rate (Pn), stomatal conductance (gs) and transpiration rate (Tr) than B9. The qRT-PCR analysis of genes encoding enzymes like SPS, SuSy, CWIN, and CeS showed upregulation in GXB9 and downregulation in B9. However, these genes were significantly differentially expressed between the immature and maturing internodes of GXB9. The cane quality trait analysis showed that GXB9 had higher juice rate, juice gravity purity, brix, juice sucrose content and cane sucrose content than B9. The yield and component investigation results indicated that GXB9 had lower single stalk weight, however higher number of millable stalks per hectare than B9, and GXB9 had lower theoretical cane yield than B9. SSR marker analysis showed genetic variation between GXB9 and B9. This study has shown significant variation in the traits of GXB9 in comparison to B9 which advocates that GXB9 is a high sugar mutant clone of B9 and an elite source for future breeding.

Identification of potential MTAs and candidate genes for juice quality- and yield-related traits in Saccharum clones: a genome-wide association and comparative genomic study

Sugarcane is an economically important commercial crop which provides raw material for the production of sugar, jaggery, bioethanol, biomass and other by-products. Sugarcane breeding till today heavily relies on conventional breeding approaches which is time consuming, laborious and costly. Integration of marker-assisted selection (MAS) in sugarcane genetic improvement programs for difficult to select traits like sucrose content, resistance to pests and diseases and tolerance to abiotic stresses will accelerate varietal development. In the present study, association mapping approach was used to identify QTLs and genes associated with sucrose and other important yield-contributing traits. A mapping panel of 110 diverse sugarcane genotypes and 148 microsatellite primers were used for structured association mapping study. An optimal subpopulation number (ΔK) of 5 was identified by structure analysis. GWAS analysis using TASSEL identified a total of 110 MTAs which were localized into 27 QTLs by GLM and MLM (Q + K, PC + K) approaches. Among the 24 QTLs sequenced, 12 were able to identify potential candidate genes, viz., starch branching enzyme, starch synthase 4, sugar transporters and G3P-DH related to carbohydrate metabolism and hormone pathway-related genes ethylene insensitive 3-like 1, reversion to ethylene sensitive1-like, and auxin response factor associated to juice quality- and yield-related traits. Six markers, NKS 5_185, SCB 270_144, SCB 370_256, NKS 46_176 and UGSM 648_245, associated with juice quality traits and marker SMC31CUQ_304 associated with NMC were validated and identified as significantly associated to the traits by one-way ANOVA analysis. In conclusion, 24 potential QTLs identified in the present study could be used in sugarcane breeding programs after further validation in larger population. The candidate genes from carbohydrate and hormone response pathway presented in this study could be manipulated with genome editing approaches to further improve sugarcane crop.

Sugarcane Ratooning Ability: Research Status, Shortcomings, and Prospects

Sugarcane is an important sugar crop and it can be subjected to ratooning for several years. The advantages of ratooning include quality improvement, efficiency enhancement, and reduced costs and energy use. The genotype, environment, cultivation management, and harvesting technology affect the productivity and longevity of ratoon cane, with the genetic basis being the most critical factor. However, the majority of research has been focused on only limited genotypes, and a few studies have evaluated up to 100 sugarcane germplasm resources. They mainly focus on the comparison among different genotypes or among plant cane, different selection strategies for the first and second ratoon crops, together with screening indicators for the selection of stronger ratooning ability. In this paper, previous studies are reviewed in order to analyze the importance of sugarcane ratooning, the indicative traits used to evaluate ratooning ability, the major factors influencing the productivity and longevity of ratooning, the genetic basis of variation in ratooning ability, and the underlying mechanisms. Furthermore, the shortcomings of the existing research on sugarcane ratooning are highlighted. We then discuss the focus of future ratoon sugarcane research and the technical methods that will shorten the selection cycle and increase the genetic gain of ratooning ability, particularly the development of linked markers. This review is expected to provide a reference for understanding the mechanisms underlying the formation of ratooning ability and for breeding sugarcane varieties with a strong ratooning ability.

Genome-wide association study of individual sugar content in fruit of Japanese pear (Pyrus spp.)

BACKGROUND

Understanding mechanisms of sugar accumulation and composition is essential to determining fruit quality and maintaining a desirable balance of sugars in plant storage organs. The major sugars in mature Rosaceae fruits are sucrose, fructose, glucose, and sorbitol. Among these, sucrose and fructose have high sweetness, whereas glucose and sorbitol have low sweetness. Japanese pear has extensive variation in individual sugar contents in mature fruit. Increasing total sugar content and that of individual high-sweetness sugars is a major target of breeding programs. The objective of this study was to identify quantitative trait loci (QTLs) associated with fruit traits including individual sugar accumulation, to infer the candidate genes underlying the QTLs, and to assess the potential of genomic selection for breeding pear fruit traits.

RESULTS

We evaluated 10 fruit traits and conducted genome-wide association studies (GWAS) for 106 cultivars and 17 breeding populations (1112 F1 individuals) using 3484 tag single-nucleotide polymorphisms (SNPs). By implementing a mixed linear model and a Bayesian multiple-QTL model in GWAS, 56 SNPs associated with fruit traits were identified. In particular, a SNP located close to acid invertase gene PPAIV3 on chromosome 7 and a newly identified SNP on chromosome 11 had quite large effects on accumulation of sucrose and glucose, respectively. We used 'Golden Delicious' doubled haploid 13 (GDDH13), an apple reference genome, to infer the candidate genes for the identified SNPs. In the region flanking the SNP on chromosome 11, there is a tandem repeat of early responsive to dehydration (ERD6)-like sugar transporter genes that might play a role in the phenotypes observed.

CONCLUSIONS

SNPs associated with individual sugar accumulation were newly identified at several loci, and candidate genes underlying QTLs were inferred using advanced apple genome information. The candidate genes for the QTLs are conserved across Pyrinae genomes, which will be useful for further fruit quality studies in Rosaceae. The accuracies of genomic selection for sucrose, fructose, and glucose with genomic best linear unbiased prediction (GBLUP) were relatively high (0.67-0.75), suggesting that it would be possible to select individuals having high-sweetness fruit with high sucrose and fructose contents and low glucose content.

Global transcriptome changes of elongating internode of sugarcane in response to mepiquat chloride

Background

Mepiquat chloride (DPC) is a chemical that is extensively used to control internode growth and create compact canopies in cultured plants. Previous studies have suggested that DPC could also inhibit gibberellin biosynthesis in sugarcane. Unfortunately, the molecular mechanism underlying the suppressive effects of DPC on plant growth is still largely unknown.

Results

In the present study, we first obtained high-quality long transcripts from the internodes of sugarcane using the PacBio Sequel System. A total of 72,671 isoforms, with N50 at 3073, were generated. These long isoforms were used as a reference for the subsequent RNA-seq. Afterwards, short reads generated from the Illumina HiSeq 4000 platform were used to compare the differentially expressed genes in both the DPC and the control groups. Transcriptome profiling showed that most significant gene changes occurred after six days post DPC treatment. These genes were related to plant hormone signal transduction and biosynthesis of several metabolites, indicating that DPC affected multiple pathways, in addition to suppressing gibberellin biosynthesis. The network of DPC on the key stage was illustrated by weighted gene co-expression network analysis (WGCNA). Among the 36 constructed modules, the top positive correlated module, at the stage of six days post spraying DPC, was sienna3. Notably, Stf0 sulfotransferase, cyclin-like F-box, and HOX12 were the hub genes in sienna3 that had high correlation with other genes in this module. Furthermore, the qPCR validated the high accuracy of the RNA-seq results.

Conclusion

Taken together, we have demonstrated the key role of these genes in DPC-induced growth inhibition in sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07352-w.

Candidate gene association of gene expression data in sugarcane contrasting for sucrose content

Association mapping of gene expression data, generated from transcriptome and proteome studies, provides a means of understanding the functional significance and trait association potential of candidate genes. In this study, we applied candidate gene association mapping to validate sugarcane genes, using data from the starch and sucrose metabolism pathway, transcriptome, and proteome. We performed multiplex PCR targeted amplicon sequencing of 109 candidate genes, using NGS technology. A range of statistical models, both single-locus and multi-locus, were compared for minimization of false positives in association mapping of four sugar-related traits with different heritability. The Fixed and random model Circulating Probability Unification model effectively suppressed false positives for both low- and high-heritability traits. We identified favorable alleles of the candidate genes involved in signalling and transcriptional regulation. The results will support genetic improvement of sugarcane and may help clarify the genetic architecture of sugar-related traits.

Linkage Map Development by EST-SSR Markers and QTL Analysis for Inflorescence and Leaf Traits in Chrysanthemum (Chrysanthemum morifolium Ramat.)

Chrysanthemums (Chrysanthemum morifolium Ramat.) are famous ornamental crops with high medicinal and industrial values. The inflorescence and leaf traits are key factors that affect the yield and quality of chrysanthemum. However, the genetic improvement of those traits is slow within chrysanthemum because of its hexaploidy, high heterozygosity and enormous genome. To study the genetic control of the important traits and facilitate marker-assisted selection (MAS) in chrysanthemum, it is desirable to populate the genetic maps with an abundance of transferrable markers such as microsatellites (SSRs). A genetic map was constructed with expressed sequence tag–simple sequence repeat (EST-SSR) markers in an F₁ progeny of 192 offspring. A total of 1000 alleles were generated from 223 EST-SSR primer pairs. The preliminary maternal and paternal maps consisted of 265 marker alleles arranged into 49 and 53 linkage groups (LGs), respectively. The recombined parental maps covered 906.3 and 970.1 cM of the genome, respectively. Finally, 264 polymorphic loci were allocated to nine LGs. The integrated map spanned 954.5 cM in length with an average genetic distance of 3.6 cM between two neighbouring loci. Quantitative trait loci (QTLs) analysis was performed using the integrated map for inflorescence diameter (ID), central disc flower diameter (CDFD), number of whorls of ray florets (NWRF), number of ray florets (NRF), number of disc florets (NDF), number of florets (NF), ray floret length (RFL), ray floret width (RFW), ray floret length/width (RFL/W), leaf length (LL), leaf width (LW) and leaf length/width (LL/W). Overall, 36 (21 major) QTLs were identified. The successful mapping of inflorescence and leaf traits QTL demonstrated the utility of the new integrated linkage map. This study is the first report of a genetic map based on EST-SSR markers in chrysanthemum. The EST-SSR markers, genetic map and QTLs reported here could be valuable resources in implementing MAS for chrysanthemums in breeding programs.

An enriched sugarcane diversity panel for utilization in genetic improvement of sugarcane

Sugarcane crop is important for both sugar and biofuels. A world collection of sugarcane and related grasses (WCSRG) maintained at Miami, FL contains > 1,200 non-redundant clones of different species and genera within the Saccharum complex. However, linkage of undesirable alleles with useful genes in wild species has hindered its efficient utilization in sugarcane breeding. A core collection developed previously with smaller number of clones representing WCSRG did not take into account > 120 wild/exotic clones maintained at the USDA-ARS Sugarcane Research Unit in Houma, Louisiana. Moreover, the genome complexity and sub-tropical to temperate growing climate of Louisiana warrant a region-specific core collection that can be used for base-broadening breeding aimed at efficient introgression of desirable alleles. Genetic diversity of 1,485 clones within WCSRG and Louisiana (commercials, wild/exotic) using 423 SSR alleles showed an average gene diversity (h) at 0.208 among all species groups where Erianthus-like Saccharum species (ELSS), Miscanthus spp., and S. spontaneum each formed a distinct cluster, Saccharum robustum, S. officinarum, hybrid cultivars, and S. edule grouped together in a major cluster, and Saccharum sinense and S. barberi formed distinct grouping. A 309-clone diversity panel (SDP1) was developed that captured the genetic diversity based on the combination of maximum length subtree and manual selection to maximize representation of Louisiana clones and minimize import of clones from Miami. SDP1 shared 324 alleles out of the 423 alleles in the entire population of 1,485 clones and captured the genetic diversity of the entire collection with an average gene diversity (h) at 0.163. The variation within (11–17%) and among (83–89%) the populations in SDP1 were comparable with the entire population of 1,485 clones (9–15% and 85–91%, respectively). The breadth of the genetic variation of SDP1 was exemplified by the intra- and inter-specific diversity of a 190-clone mini-core collection with markers derived from known cold-responsive genes. SDP1 will facilitate genome-wide association studies for identification of trait-specific markers for use in marker-assisted breeding in Louisiana and elsewhere.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Association mapping for yield traits in Paeonia rockii based on SSR markers within transcription factors of comparative transcriptome

BACKGROUND

Allelic variation underlying the quantitative traits in plants is caused by the extremely complex regulation process. Tree peony originated in China is a peculiar ornamental, medicinal and oil woody plant. Paeonia rockii, one of tree peony species, is a precious emerging woody oil crop. However, in this valuable plant, the study of functional loci associated with yield traits has rarely been identified. Therefore, to explore the genetic architecture of 24 yield quantitative traits, the association mapping was first reported in 420 unrelated cultivated P. rockii individuals based on the next-generation sequencing (NGS) and single-molecule long-read sequencing (SMLRS).

RESULTS

The developed 58 pairs of polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers from 959 candidate transcription factors (TFs) associated with yield were used for genotyping the 420 P. rockii accessions. We observed a high level of genetic diversity (polymorphic information content, PIC = 0.514) and low linkage disequilibrium (LD) between EST-SSRs. Moreover, four subpopulations in the association population were revealed by STRUCTURE analyses. Further, single-marker association analysis identified 141 significant associations, involving 17 quantitative traits and 41 EST-SSRs. These loci were mainly from AP2, TCP, MYB, HSF, bHLH, GATA, and B3 gene families and showed a small proportion of the phenotypic variance (3.79 to 37.45%).

CONCLUSIONS

Our results summarize a valuable collection of functional loci associated with yield traits in P. rockii, and provide a precious resource that reveals allelic variation underlying quantitative traits in Paeonia and other woody oil crops.

Development and characterization of 20 novel EST-SSR markers for Pteroceltis tatarinowii, a relict tree in China

PREMISE

Pteroceltis tatarinowii (Ulmaceae), the only species of the genus Pteroceltis, is an endangered tree in China. Here, novel expressed sequence tag-simple sequence repeat (EST-SSR) markers were developed to illuminate its genetic diversity for conservation and assisted breeding.

METHODS AND RESULTS

Based on Illumina transcriptome data from P. tatarinowii, a total of 70 EST-SSR markers were initially designed and tested. Forty-eight of 70 loci (68.6%) were successfully amplified, of which 20 were polymorphic. The number of alleles per locus ranged from two to six, and the levels of observed and expected heterozygosity ranged from 0.018 to 0.781 and from 0.023 to 0.702, respectively. Additionally, cross-amplification was successful for 17 loci in two related species, Ulmus gaussenii and U. chenmoui.

CONCLUSIONS

These new EST-SSR markers are valuable transcriptomic resources for P. tatarinowii and will facilitate population genetics and molecular breeding of this species and its relatives in Ulmaceae.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Transcriptome landscaping for gene mining and SSR marker development in Coriander (Coriandrum sativum L.)

Coriander (Coriandrum sativum L.) is an aromatic herb, widely used as a spice and is of great pharmaceutical interest. Despite high medicinal and economic value, there is a dearth of genomic information about profiling as well as the expressed sequence-based genic markers. In this study, transcriptome was sequenced from seeds, leaves, and flower for gene mining and identification of SSR markers. A total of 9746 SSR containing loci were identified, the most abundant type of SSR identified were the di-nucleotide repeat motifs (45.5%), followed by tri- (34.6%), tetra- (4.5%), penta- (1.5%) and hexanucleotide repeats (1%). A total of 3795 primers were designed, out of which 120 randomly selected were validated in 14 accessions of coriander cultivated in India. The current study provides useful information about preliminary transcriptome sketch and genic markers, which can be useful in breeding and genetic diversity estimation of coriander.

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.