Molecular characterization of the Gossypium Diversity Reference Set of the US National Cotton Germplasm Collection

Fatty acid unsaturation improves germination of upland cotton (Gossypium hirsutum) under cold stress

Introduction

The level of fatty acid unsaturation in seeds is one of the major determinants of cold germination ability, particularly in oilseeds. The presence of cis double bonds in unsaturated fatty acids creates bends that lowers their melting temperatures compared to saturated fatty acids. Unsaturated fatty acids with low melting points mobilize faster at low temperatures providing seeds with sufficient energy for germination.

Methodology

To investigate the effects of fatty acid unsaturation on the ability of cotton seeds to germinate under cold conditions, four recombinant inbred lines (RILs) of cotton with unique fatty acid profiles were evaluated using a set of developmental and biochemical assays at 12°C (critically low temperature), 15°C (cardinal minimum temperature) and 30°C (optimum temperature). Furthermore, whole seed lipidome profiling using liquid chromatography with mass spectrometry was done to compare the lipid compositional changes at 12°C and 30°C after imbibing cotton seeds of all the six genotypes for 0 hours, 3 hours and 6 hours.

Results and discussion

The RILs with higher unsaturation/saturation ratios registered robust germination performance, lower solute leakage, and optimum water uptake rates under cold stress. Imbibition at 30°C for 8 hours before cold exposure significantly improved the germination of cold sensitive genotypes, indicating that the first few hours of water uptake are critical for cold stress. Whole seed lipidome profiling of all the genotypes specifically associated cold germination ability with higher unsaturation levels of phospholipids during early imbibition. The presence of cis double bonds in phospholipids creates kinks that maintain the fluidity of cell membranes under low temperature. Membrane flexibility under cold conditions is essential for facilitating key germination events including membrane organization and respiration. The current results highlight the importance of fatty acid composition in cold germination ability of upland cotton.

Genetic polymorphism detection in brazilian perennial cottons (Gossypium spp.) using an ISSR marker system and its application for molecular interspecific differentiation

BACKGROUND

The semi-domesticated Brazilian perennial cotton (Gossypium spp.) germplasm is considered a source of variability for creating modern upland cotton varieties. Here we used Inter-simple Sequence Repeat (ISSR) markers to detect intra and interspecific genetic polymorphism in Gossypium hirsutum L. r. marie-galante and Gossypium barbadense L. and to use molecular data to assessing genetic diversity and molecular discrimination of these species.

METHODS AND RESULTS

The sets contained 12 G. barbadense genotypes and 16 G. hirsutum genotypes from a Brazilian collection. The 11 ISSR primers were used for genotyping yielded 101 bands (polymorphism = 47.5%) and were classified as moderately informative (PIC = 0.304). The ISSR markers exposed a greater diversity in G. hirsutum (P = 24.72%; H_E =0.071 and I = 0.111) as compared to G. barbadense (P = 17.98%, H_E = 0.043 and I = 0.070). The AMOVA analysis showed that 89.47% of the genetic variation was partitioned within species which is supported by Nei's genetic differentiation (G_st = 0.598) and gene flow (N_m = 0.338), suggesting that strong reproductive barriers between species. The UPGMA Cluster Analysis, Principal Coordinate Analysis and Bayesian Model-Based Structural Analysis divided the 28 genotypes into two main clades consistent with the taxonomical delimitation.

CONCLUSION

The ISSR marker system offers a new approach to determining molecular differences between two cotton species (G. hirsutum L. r. marie-galante and G. barbadense L.). This study can expand the molecular marker resources for the identification and improvement of our knowledge about the genetic diversity and relationships between perennial cotton genotypes.

Genome-Wide Dissection of the Genetic Basis for Drought Tolerance in Gossypium hirsutum L. Races

Drought seriously threats the growth and development of Gossypium hirsutum L. To dissect the genetic basis for drought tolerance in the G. hirsutum L. germplasm, a population, consisting of 188 accessions of G. hirsutum races and a cultivar (TM-1), was genotyped using the Cotton80KSNP biochip, and 51,268 high-quality single-nucleotide polymorphisms (SNPs) were obtained. Based on the phenotypic data of eight drought relative traits from four environments, we carried out association mapping with five models using GAPIT software. In total, thirty-six SNPs were detected significantly associated at least in two environments or two models. Among these SNPs, 8 and 28 (including 24 SNPs in 5 peak regions) were distributed in the A and D subgenome, respectively; eight SNPs were found to be distributed within separate genes. An SNP, TM73079, located on chromosome D10, was simultaneously associated with leaf fresh weight, leaf wilted weight, and leaf dry weight. Another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported quantitative trait loci for drought tolerance. Moreover, 520 genes located 200 kb up- and down-stream of 36 SNPs were obtained and analyzed based on gene annotation and transcriptome sequencing. The results showed that three candidate genes, Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369, may play important roles in drought tolerance. The current GWAS represents the first investigation into mapping QTL for drought tolerance in G. hirsutum races and provides important information for improving cotton cultivars.

Integrated morpho-biochemical and transcriptome analyses reveal multidimensional response of upland cotton (Gossypium hirsutum L.) to low temperature stress during seedling establishment

Cotton is a tropical/subtropical crop and is innately susceptible to cold. Using an approach that integrates morphological, biochemical, and transcriptome analyses, the study aimed to understand the molecular underpinnings of phenotypic adjustments in cotton seedlings under cold stress. Exposure of six cotton accessions to 15°C during the seedling stage significantly reduced chlorophyll content, stomatal conductance, plant height, and biomass, but increased malondialdehyde and proline production. Comparative transcriptome profiling of the cold-sensitive accession SA 3781 grown under low and normal temperatures showed the upregulation of genes related to the production of reactive oxygen species (ROS) under cold stress. Despite a similar upregulation of genes encoding metabolites that can scavenge ROS and provide osmoprotection for the cell, the stressed plants still exhibited oxidative stress in terms of lipid peroxidation. This may be due in part to the upregulation of abscisic acid synthesis genes and downregulation of chlorophyll synthesis genes effecting lower stomatal conductance and chlorophyll contents, respectively. Additionally, stomatal closure which is required to avoid the cooling effect and dehydration under cold conditions may have contributed in reducing the net photosynthetic rates in plants exposed to low temperature. These findings provide an insight into the expression of key genes regulating the phenotypic changes observed in cotton in response to cold stress.

Finding Needles in a Haystack: Using Geo-References to Enhance the Selection and Utilization of Landraces in Breeding for Climate-Resilient Cultivars of Upland Cotton (Gossypium hirsutum L.)

The genetic uniformity of cultivated cotton as a consequence of domestication and modern breeding makes it extremely vulnerable to abiotic challenges brought about by major climate shifts. To sustain productivity amidst worsening agro-environments, future breeding objectives need to seriously consider introducing new genetic variation from diverse resources into the current germplasm base of cotton. Landraces are genetically heterogeneous, population complexes that have been primarily selected for their adaptability to specific localized or regional environments. This makes them an invaluable genetic resource of novel allelic diversity that can be exploited to enhance the resilience of crops to marginal environments. The utilization of cotton landraces in breeding programs are constrained by the phenology of the plant and the lack of phenotypic information that can facilitate efficient selection of potential donor parents for breeding. In this review, the genetic value of cotton landraces and the major challenges in their utilization in breeding are discussed. Two strategies namely Focused Identification of Germplasm Strategy and Environmental Association Analysis that have been developed to effectively screen large germplasm collections for accessions with adaptive traits using geo-reference-based, mathematical modelling are highlighted. The potential applications of both approaches in mining available cotton landrace collections are also presented.

Parallel and Intertwining Threads of Domestication in Allopolyploid Cotton

The two cultivated allopolyploid cottons, Gossypium hirsutum and Gossypium barbadense, represent a remarkable example of parallel independent domestication, both involving dramatic morphological transformations under selection from wild perennial plants to annualized row crops. Deep resequencing of 643 newly sampled accessions spanning the wild-to-domesticated continuum of both species, and their allopolyploid relatives, are combined with existing data to resolve species relationships and elucidate multiple aspects of their parallel domestication. It is confirmed that wild G. hirsutum and G. barbadense were initially domesticated in the Yucatan Peninsula and NW South America, respectively, and subsequently spread under domestication over 4000-8000 years to encompass most of the American tropics. A robust phylogenomic analysis of infraspecific relationships in each species is presented, quantify genetic diversity in both, and describe genetic bottlenecks associated with domestication and subsequent diffusion. As these species became sympatric over the last several millennia, pervasive genome-wide bidirectional introgression occurred, often with striking asymmetries involving the two co-resident genomes of these allopolyploids. Diversity scans revealed genomic regions and genes unknowingly targeted during domestication and additional subgenomic asymmetries. These analyses provide a comprehensive depiction of the origin, divergence, and adaptation of cotton, and serve as a rich resource for cotton improvement.

Use of hydraulic traits for modeling genotype-specific acclimation in cotton under drought

Understanding the genetic and physiological basis of abiotic stress tolerance under field conditions is key to varietal crop improvement in the face of climate variability. Here, we investigate dynamic physiological responses to water stress in silico and their relationships to genotypic variation in hydraulic traits of cotton (Gossypium hirsutum), an economically important species for renewable textile fiber production. In conjunction with an ecophysiological process-based model, heterogeneous data (plant hydraulic traits, spatially-distributed soil texture, soil water content and canopy temperature) were used to examine hydraulic characteristics of cotton, evaluate their consequences on whole plant performance under drought, and explore potential genotype × environment effects. Cotton was found to have R-shaped hydraulic vulnerability curves (VCs), which were consistent under drought stress initiated at flowering. Stem VCs, expressed as percent loss of conductivity, differed across genotypes, whereas root VCs did not. Simulation results demonstrated how plant physiological stress can depend on the interaction between soil properties and irrigation management, which in turn affect genotypic rankings of transpiration in a time-dependent manner. Our study shows how a process-based modeling framework can be used to link genotypic variation in hydraulic traits to differential acclimating behaviors under drought.

Historical Datasets Support Genomic Selection Models for the Prediction of Cotton Fiber Quality Phenotypes Across Multiple Environments

Genomic selection (GS) has successfully been used in plant breeding to improve selection efficiency and reduce breeding time and cost. However, there has not been a study to evaluate GS prediction models that may be used for predicting cotton breeding lines across multiple environments. In this study, we evaluated the performance of Bayes Ridge Regression, BayesA, BayesB, BayesC and Reproducing Kernel Hilbert Spaces regression models. We then extended the single-site GS model to accommodate genotype × environment interaction (G×E) in order to assess the merits of multi- over single-environment models in a practical breeding and selection context in cotton, a crop for which this has not previously been evaluated. Our study was based on a population of 215 upland cotton (Gossypium hirsutum) breeding lines which were evaluated for fiber length and strength at multiple locations in Australia and genotyped with 13,330 single nucleotide polymorphic (SNP) markers. BayesB, which assumes unique variance for each marker and a proportion of markers to have large effects, while most other markers have zero effect, was the preferred model. GS accuracy for fiber length based on a single-site model varied across sites, ranging from 0.27 to 0.77 (mean = 0.38), while that of fiber strength ranged from 0.19 to 0.58 (mean = 0.35) using randomly selected sub-populations as the training population. Prediction accuracies from the M×E model were higher than those for single-site and across-site models, with an average accuracy of 0.71 and 0.59 for fiber length and strength, respectively. The use of the M×E model could therefore identify which breeding lines have effects that are stable across environments and which ones are responsible for G×E and so reduce the amount of phenotypic screening required in cotton breeding programs to identify adaptable genotypes.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Genome-wide divergence, haplotype distribution and population demographic histories for Gossypium hirsutum and Gossypium barbadense as revealed by genome-anchored SNPs

Use of 10,129 singleton SNPs of known genomic location in tetraploid cotton provided unique opportunities to characterize genome-wide diversity among 440 Gossypium hirsutum and 219 G. barbadense cultivars and landrace accessions of widespread origin. Using the SNPs distributed genome-wide, we examined genetic diversity, haplotype distribution and linkage disequilibrium patterns in the G. hirsutum and G. barbadense genomes to clarify population demographic history. Diversity and identity-by-state analyses have revealed little sharing of alleles between the two cultivated allotetraploid genomes, with a few exceptions that indicated sporadic gene flow. We found a high number of new alleles, representing increased nucleotide diversity, on chromosomes 1 and 2 in cultivated G. hirsutum as compared with low nucleotide diversity on these chromosomes in landrace G. hirsutum. In contrast, G. barbadense chromosomes showed negative Tajima's D on several chromosomes for both cultivated and landrace types, which indicate that speciation of G. barbadense itself, might have occurred with relatively narrow genetic diversity. The presence of conserved linkage disequilibrium (LD) blocks and haplotypes between G. hirsutum and G. barbadense provides strong evidence for comparable patterns of evolution in their domestication processes. Our study illustrates the potential use of population genetic techniques to identify genomic regions for domestication.

Independent Domestication of Two Old World Cotton Species

Domesticated cotton species provide raw material for the majority of the world's textile industry. Two independent domestication events have been identified in allopolyploid cotton, one in Upland cotton (Gossypium hirsutum L.) and the other to Egyptian cotton (Gossypium barbadense L.). However, two diploid cotton species, Gossypium arboreum L. and Gossypium herbaceum L., have been cultivated for several millennia, but their status as independent domesticates has long been in question. Using genome resequencing data, we estimated the global abundance of various repetitive DNAs. We demonstrate that, despite negligible divergence in genome size, the two domesticated diploid cotton species contain different, but compensatory, repeat content and have thus experienced cryptic alterations in repeat abundance despite equivalence in genome size. Evidence of independent origin is bolstered by estimates of divergence times based on molecular evolutionary analysis of f7,000 orthologous genes, for which synonymous substitution rates suggest that G. arboreum and G. herbaceum last shared a common ancestor approximately 0.4-2.5 Ma. These data are incompatible with a shared domestication history during the emergence of agriculture and lead to the conclusion that G. arboreum and G. herbaceum were each domesticated independently.

Genetic Diversity of the Two Commercial Tetraploid Cotton Species in the Gossypium Diversity Reference Set

A diversity reference set has been constructed for the Gossypium accessions in the US National Cotton Germplasm Collection to facilitate more extensive evaluation and utilization of accessions held in the Collection. A set of 105 mapped simple sequence repeat markers was used to study the allelic diversity of 1933 tetraploid Gossypium accessions representative of the range of diversity of the improved and wild accessions of G. hirsutum and G. barbadense. The reference set contained 410 G. barbadense accessions and 1523 G. hirsutum accessions. Observed numbers of polymorphic and private bands indicated a greater diversity in G. hirsutum as compared to G. barbadense as well as in wild-type accessions as compared to improved accessions in both species. The markers clearly differentiated the 2 species. Patterns of diversity within species were observed but not clearly delineated, with much overlap occurring between races and regions of origin for wild accessions and between historical and geographic breeding pools for cultivated accessions. Although the percentage of accessions showing introgression was higher among wild accessions than cultivars in both species, the average level of introgression within individual accessions, as indicated by species-specific bands, was much higher in wild accessions of G. hirsutum than in wild accessions of G. barbadense. The average level of introgression within individual accessions was higher in improved G. barbadense cultivars than in G. hirsutum cultivars. This molecular characterization reveals the levels and distributions of genetic diversity that will allow for better exploration and utilization of cotton genetic resources.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.