Genetic Basis of Fiber Improvement and Decreased Stress Tolerance in Cultivated Versus Semi-Domesticated Upland Cotton

Genetic Analysis of an F2 Population Derived from the Cotton Landrace Hopi Identified Novel Loci for Boll Glanding

Landraces are an important reservoir of genetic variation that can expand the narrow genetic base of cultivated cotton. In this study, quantitative trait loci (QTL) analysis was conducted using an F2 population developed from crosses between the landrace Hopi and inbred TM-1. A high-density genetic map spanning 2253.11 and 1932.21 cM for the A and D sub-genomes, respectively, with an average marker interval of 1.14 cM, was generated using the CottonSNP63K array. The linkage map showed a strong co-linearity with the physical map of cotton. A total of 21 QTLs were identified, controlling plant height (1), bract type (1), boll number (1), stem color (2), boll pitting (2), fuzz fiber development (2), boll shape (3), boll point (4), and boll glanding (5). In silico analysis of the novel QTLs for boll glanding identified a total of 13 candidate genes. Analysis of tissue-specific expression of the candidate genes suggests roles for the transcription factors bHLH1, MYB2, and ZF1 in gland formation. Comparative sequencing of open reading frames identified early stop codons in all three transcription factors in Hopi. Functional validation of these genes offers avenues to reduce glanding and, consequently, lower gossypol levels in cottonseeds without compromising the defense mechanisms of the plant against biotic stresses.

Transcriptome Revealed GhPP2C43-A Negatively Regulates Salinity Tolerance in an Introgression Line from a Semi-wild Upland Cotton

Salt damage is one of the major threats to sustainable cotton production owing to the limited arable land in China mainly occupied by the production of staple food crops. Salt-stress tolerant cotton varieties are lacking in production and, the mechanisms underpinning salt-stress tolerance in cotton remain enigmatic. Here, DM37, an intraspecific introgression line from G. hirsutum race yucatanense acc TX-1046 into the G. hirsutum acc TM-1 background, was found to be highly tolerant to salt stress. Its seed germination rate and germination potential were significantly higher than the recipient TM-1 under salt stress. Physiological analysis showed DM37 had higher proline content and Peroxidase activity, as well as lower Na+/K+ ratios at the seedling stage, consistent with higher seedling survival rate after durable salt stress. Furthermore, comparative transcriptome analysis revealed that responsive patterns to salt stress in DM37 were different from TM-1. Weighted Correlation Network Analysis (WGCNA) demonstrated that co-expression modules associated with salt stress in DM37 also differed from TM-1. Out of them, GhPP2C43-A, a phosphatase gene, exhibited negative regulation of salt-stress tolerance verified by VIGS and transgenic Arabidopsis. Gene expression showed GhPP2C43-A in TM-1 was induced by durable salt stress but not in DM37 probably attributing to the variation of cis-element in its promoter, thereby being conferred different salt-stress tolerance. Our result would provide new genes/germplasms from semi-wild cotton in salt-stress tolerant cotton breeding. This study would give us new insights into the mechanisms underpinning the salt-stress tolerance in cotton.

Regulation of inositol 1,2,4,5,6-pentakisphosphate and inositol hexakisphosphate levels in Gossypium hirsutum by IPK1

The IPK1 genes, which code for 2-kinases that can synthesize Ins(1,2,4,5,6)P₅ from Ins(1,4,5,6)P₄, are expressed throughout cotton plants, resulting in the highest Ins(1,2,4,5,6)P₅ concentrations in young leaves and flower buds. Cotton leaves contain large amounts of Ins(1,2,4,5,6)P₅ and InsP₆ compared to plants not in the Malvaceae family. The inositol polyphosphate pathway has been linked to stress tolerance in numerous plant species. Accordingly, we sought to determine why cotton and other Malvaceae have such high levels of these inositol phosphates. We have quantified the levels of InsP₅ and InsP₆ in different tissues of cotton plants and determined the expression of IPK1 (inositol 1,3,4,5,6-pentakisphosphate 2-kinase gene) in vegetative and reproductive tissues. Gossypium hirsutum was found to contain four IPK1 genes that were grouped into two pair (AB, CD) where each pair consists of very similar sequences that were measured together. More IPK1AB is expressed in leaves than in roots, whereas more IPK1CD is expressed in roots than in leaves. Leaves and flower buds have more InsP₅ and InsP₆ than stems and roots. Leaves and roots contain more InsP₅ than InsP₆, whereas flower buds and stems contain more InsP₆ than InsP₅. Dark-grown seedlings contain more InsP₅ and InsP₆ than those grown under lights, and the ratio of InsP₅ to InsP₆ is greater in the light-grown seedlings. During 35 days of the life cycle of the third true leaf, InsP₅ and InsP₆ gradually decreased by more than 50%. Silencing IPK1AB and IPK1CD with Cotton Leaf Crumple Virus-induced gene silencing (VIGS) resulted in plants with an intense viral phenotype, reduced IPK1AB expression and lowered amounts of InsP₅. The results are consistent with Ins(1,2,4,5,6)P₅ synthesis from Ins(1,4,5,6)P₄ by IPK1. This study detailed the central role of IPK1 in cotton inositol polyphosphate metabolism, which has potential to be harnessed to improve the resistance of plants to different kinds of stress.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

The Complete Chloroplast Genome Sequences of Eight Fagopyrum Species: Insights Into Genome Evolution and Phylogenetic Relationships

Buckwheat (Fagopyrum genus, Polygonaceae), is an annual or perennial, herbaceous or semi-shrub dicotyledonous plant. There are mainly three cultivated buckwheat species, common buckwheat (Fagopyrum esculentum) is widely cultivated in Asia, Europe, and America, while Tartary buckwheat (F. tataricum) and F. cymosum (also known as F. dibotrys) are mainly cultivated in China. The genus Fagopyrum is taxonomically confusing due to the complex phenotypes of different Fagopyrum species. In this study, the chloroplast (cp) genomes of three Fagopyrum species, F. longistylum, F. leptopodum, F. urophyllum, were sequenced, and five published cp genomes of Fagopyrum were retrieved for comparative analyses. We determined the sequence differentiation, repeated sequences of the cp genomes, and the phylogeny of Fagopyrum species. The eight cp genomes ranged, gene number, gene order, and GC content were presented. Most of variations of Fagopyrum species cp genomes existed in the LSC and SSC regions. Among eight Fagopyrum chloroplast genomes, six variable regions (ndhF-rpl32, trnS-trnG, trnC, trnE-trnT, psbD, and trnV) were detected as promising DNA barcodes. In addition, a total of 66 different SSR (simple sequence repeats) types were found in the eight Fagopyrum species, ranging from 8 to 16 bp. Interestingly, many SSRs showed significant differences especially in some photosystem genes, which provided valuable information for understanding the differences in light adaptation among different Fagopyrum species. Genus Fagopyrum has shown a typical branch that is distinguished from the Rumex, Rheum, and Reynoutria, which supports the unique taxonomic status in Fagopyrum among the Polygonaceae. In addition, phylogenetic analysis based on the cp genomes strongly supported the division of eight Fagopyrum species into two independent evolutionary directions, suggesting that the separation of cymosum group and urophyllum group may be earlier than the flower type differentiation in Fagopyrum plants. The results of the chloroplast-based phylogenetic tree were further supported by the matK and Internal Transcribed Spacer (ITS) sequences of 17 Fagopyrum species, which may help to further anchor the taxonomic status of other members in the urophyllum group in Fagopyrum. This study provides valuable information and high-quality cp genomes for identifying species and evolutionary analysis for future Fagopyrum research.

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.

Dissecting Source-Sink Relationship of Subtending Leaf for Yield and Fiber Quality Attributes in Upland Cotton (Gossypium hirsutum L.)

Photosynthesis as a source is a significant contributor to the reproductive sink affecting cotton yield and fiber quality. Moreover, carbon assimilation from subtending leaves adds up a significant proportion to the reproductive sink. Therefore, this study aimed to address the source-sink relationship of boll subtending leaf with fiber quality and yield related traits in upland cotton. A core collection of 355 upland cotton accessions was subjected to subtending leaf removal treatment effects across 2 years. The analysis of variance suggested a significant effect range in the source-sink relationship under subtending leaf removal effects at different growth stages. Further insight into the variation was provided by the correlation analysis and principal component analysis. A significant positive correlation between different traits was observed and the multivariate analysis including hierarchical clustering and principal component analysis (PCA) categorised germplasm accessions into three groups on the basis of four subtending leaf removal treatment effects across 2 years. A set of genotypes with the lowest and highest treatment effects has been identified. Selected accessions and the outcome of the current study may provide a basis for a further study to explore the molecular mechanism of source-sink relationship of boll subtending leaf and utilization of breeding programs focused on cotton improvement.

Genome-wide association analysis reveals quantitative trait loci and candidate genes involved in yield components under multiple field environments in cotton (Gossypium hirsutum)

BACKGROUND

Numerous quantitative trait loci (QTLs) and candidate genes associated with yield-related traits have been identified in cotton by genome-wide association study (GWAS) analysis. However, most of the phenotypic data were from a single or few environments, and the stable loci remained to be validated under multiple field environments.

RESULTS

Here, 242 upland cotton accessions collected from different origins were continuously investigated for phenotypic data of four main yield components, including boll weight (BW) and lint percentage (LP) under 13 field environments, and boll number per plant (BN) and seed index (SI) under 11 environments. Correlation analysis revealed a positive correlation between BN and LP, BW and SI, while SI had a negative correlation with LP and BN. Genetic analysis indicated that LP had the highest heritability estimates of 94.97%, followed by 92.08% for SI, 86.09% for BW, and 72.92% for BN, indicating LP and SI were more suitable traits for genetic improvement. Based on 56,010 high-quality single nucleotide polymorphisms (SNPs) and GWAS analysis, a total of 95 non-redundant QTLs were identified, including 12 of BN, 23 of BW, 45 of LP, and 33 of SI, respectively. Of them, 10 pairs of homologous QTLs were detected between A and D sub-genomes. We also found that 15 co-located QTLs with more than two traits and 12 high-confidence QTLs were detected under more than six environments, respectively. Further, two NET genes (GH_A08G0716 and GH_A08G0783), located in a novel QTL hotspot (qtl24, qtl25 and qlt26) were predominately expressed in early fiber development stages, exhibited significant correlation with LP and SI. The GH_A07G1389 in the stable qtl19 region encoded a tetratricopeptide repeat (TPR)-like superfamily protein and was a homologous gene involved in short fiber mutant ligon lintless-y (Li_y), implying important roles in cotton yield.

CONCLUSIONS

The present study provides a foundation for understanding the regulatory mechanisms of yield components and may enhance yield improvement through molecular breeding in cotton.

Evolutionary Origins of Drought Tolerance in Spermatophytes

It is commonly known that drought stress is a major constraint limiting crop production. Drought stress and associated drought tolerance mechanisms are therefore under intense investigation with the view to future production of drought tolerant crops. With an ever-growing population and variable climate, novel approaches need to be considered to sustainably feed future generations. In this context, definitions of drought tolerance are highly variable, which poses a major challenge for the systematic assessment of this trait across the plant kingdom. Furthermore, drought tolerance is a polygenic trait and understanding the evolution of this complex trait may inform us about patterns of gene gain and loss in relation to diverse drought adaptations. We look at the transition of plants from water to land, and the role of drought tolerance in enabling this transition, before discussing the first drought tolerant plant and common drought responses amongst vascular plants. We reviewed the distribution of a combined "drought tolerance" trait in very broad terms to encompass different experimental systems and definitions used in the current literature and assigned a binary trait "tolerance vs. sensitivity" in 178 extant plant species. By simplifying drought responses of plants into this "binary" trait we were able to explore the evolution of drought tolerance across the wider plant kingdom, compared to previous studies. We show how this binary "drought tolerance/sensitivity" trait has evolved and discuss how incorporating this information into an evolutionary genomics framework could provide insights into the molecular mechanisms underlying extreme drought adaptations.

Genome-Wide Association Analysis Reveals Loci and Candidate Genes Involved in Fiber Quality Traits Under Multiple Field Environments in Cotton (Gossypium hirsutum)

Fiber length, fiber strength, and fiber micronaire are the main fiber quality parameters in cotton. Thus, mining the elite and stable loci/alleles related to fiber quality traits and elucidating the relationship between the two may accelerate genetic improvement of fiber quality in cotton. Here, genome-wide association analysis (GWAS) was performed for fiber quality parameters based on phenotypic data, and 56,010 high-quality single nucleotide polymorphisms (SNPs) using 242 upland cotton accessions under 12 field environments were obtained. Phenotypic analysis exhibited that fiber length (FL) had a positive correlation with fiber strength (FS) and had a negative correlation with fiber micronaire (Mic). Genetic analysis also indicated that FL, FS, and Mic had high heritability of more than 80%. A total of 67 stable quantitative trait loci (QTLs) were identified through GWAS analysis, including 31 for FL, 21 for FS, and 22 for Mic. Of them, three pairs homologous QTLs were detected between A and D subgenomes, and seven co-located QTLs with two fiber quality parameters were found. Compared with the reported QTLs, 34 co-located with previous studies, and 33 were newly revealed. Integrated with transcriptome analysis, we selected 256, 244, and 149 candidate genes for FL, FS, and Mic, respectively. Gene Ontology (GO) analysis showed that most of the genes located in QTLs interval of the three fiber quality traits were involved in sugar biosynthesis, sugar metabolism, microtubule, and cytoskeleton organization, which played crucial roles in fiber development. Through correlation analysis between haplotypes and phenotypes, three genes (GH_A05G1494, GH_D11G3097, and GH_A05G1082) predominately expressed in fiber development stages were indicated to be potentially responsible for FL, FS, and Mic, respectively. The GH_A05G1494 encoded a protein containing SGS-domain, which is related to tubulin-binding and ubiquitin-protein ligase binding. The GH_D11G3097 encoded 20S proteasome beta subunit G1, and was involved in the ubiquitin-dependent protein catabolic process. The GH_A05G1082 encoded RAN binding protein 1 with a molecular function of GTPase activator activity. These results provide new insights and candidate loci/genes for the improvement of fiber quality in cotton.