Wenfang Gong

Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits

Upland cotton (Gossypium hirsutum) is the most important natural fiber crop in the world. The overall genetic diversity among cultivated species of cotton and the genetic changes that occurred during their improvement are poorly understood. Here we report a comprehensive genomic assessment of modern improved upland cotton based on the genome-wide resequencing of 318 landraces and modern improved cultivars or lines. We detected more associated loci for lint yield than for fiber quality, which suggests that lint yield has stronger selection signatures than other traits. We found that two ethylene-pathway-related genes were associated with increased lint yield in improved cultivars. We evaluated the population frequency of each elite allele in historically released cultivar groups and found that 54.8% of the elite genome-wide association study (GWAS) alleles detected were transferred from three founder landraces: Deltapine 15, Stoneville 2B and Uganda Mian. Our results provide a genomic basis for improving cotton cultivars and for further evolutionary analysis of polyploid crops.

Na(+) compartmentalization related to salinity stress tolerance in upland cotton (Gossypium hirsutum) seedlings.

The capacity for ion compartmentalization among different tissues and cells is the key mechanism regulating salt tolerance in plants. In this study, we investigated the ion compartmentalization capacity of two upland cotton genotypes with different salt tolerances under salt shock at the tissue, cell and molecular levels. We found that the leaf glandular trichome could secrete more salt ions in the salt-tolerant genotype than in the sensitive genotype, demonstrating the excretion of ions from tissue may be a new mechanism to respond to short-term salt shock. Furthermore, an investigation of the ion distribution demonstrated that the ion content was significantly lower in critical tissues and cells of the salt-tolerant genotype, indicating the salt-tolerant genotype had a greater capacity for ion compartmentalization in the shoot. By comparing the membrane H+-ATPase activity and the expression of ion transportation-related genes, we found that the H+-ATPase activity and Na+/H+ antiporter are the key factors determining the capacity for ion compartmentalization in leaves, which might further determine the salt tolerance of cotton. The novel function of the glandular trichome and the comparison of Na+ compartmentalization between two cotton genotypes with contrasting salt tolerances provide a new understanding of the salt tolerance mechanism in cotton.

MicroRNA and mRNA expression profiling analysis revealed the regulation of plant height in Gossypium hirsutum

BackgroundDwarf cottons are more resistant to damage from wind and rain and associated with stable, increased yields, and also desirable source for breeding the machine harvest varieties. In an effort to uncover the transcripts and miRNA networks involved in plant height, the transcriptome and small RNA sequencing were performed based on dwarf mutant Ari1327 (A1), tall-culm mutant Ari3697 (A3) and wild type Ari971 (A9) in Gossypium hirsutum.MethodsThe stem apexes of wild-type upland cotton (Ari971) and its dwarf mutant (Ari1327) and tall-culm mutant (Ari3697) at the fifth true leaf stage were extracted for RNA, respectively. Transcriptome and small RNA libraries were constructed and subjected to next generation sequencing.ResultsThe transcriptome sequencing analysis showed that the enriched pathways of top 3 differentially expressed genes (DEGs) were categorized as carotenoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction in both A1–A9 and A3–A9. The ABA and IAA related factors were differentially expressed in the mutants. Importantly, we found the lower expressed SAUR and elevated expressed GH3, and ABA related genes such as NCED and PP2C maybe relate to reduced growth of the plant height in Ari1327 which was consistent with the higher auxin and ABA content in this mutant. Furthermore, miRNA160 targeted to the auxin response factor (ARF) and miRNA166 (gma-miR166u and gma-miR166h-3p) targeted to ABA responsive element binding factor were related to the mutation in cotton. We have noticed that the cell growth related factors (smg7 targeted by gra-miR482 and 6 novel miRNAs and pectate-lyases targeted by osa-miR159f), the redox reactions related factors (Cytochrome P450 targeted by miR172) and MYB genes targeted by miR828, miR858 and miR159 were also involved in plant height of the cotton mutants. A total of 226 conserved miRNAs representing 32 known miRNA families were obtained, and 38 novel miRNAs corresponding to 23 unique RNA sequences were identified. Total 531 targets for 211 conserved miRNAs were obtained. Using PAREsnip, 27 and 29 miRNA/target conserved interactions were validated in A1–A9 and A3–A9, respectively. Furthermore, miRNA160, miRNA858 and miRNA172 were validated to be up-regulated in A1–A9 but down-regulated in A3–A9, whereas miRNA159 showed the opposite regulation.ConclusionsThis comprehensive interaction of the transcriptome and miRNA at tall-culm and dwarf mutant led to the discovery of regulatory mechanisms in plant height. It also provides the basis for in depth analyses of dwarf mutant genes for further breeding of dwarf cotton.

Comparison of the Transcriptome between Two Cotton Lines of Different Fiber Color and Quality

To understand the mechanism of fiber development and pigmentation formation, the mRNAs of two cotton lines were sequenced: line Z128 (light brown fiber) was a selected mutant from line Z263 (dark brown fiber). The primary walls of the fiber cell in both Z263 and Z128 contain pigments; more pigments were laid in the lumen of the fiber cell in Z263 compared with that in Z128. However, Z263 contained less cellulose than Z128. A total of 71,895 unigenes were generated: 13,278 (20.26%) unigenes were defined as differentially expressed genes (DEGs) by comparing the library of Z128 with that of Z263; 5,345 (8.16%) unigenes were up-regulated and 7,933 (12.10%) unigenes were down-regulated. qRT-PCR and comparative transcriptional analysis demonstrated that the pigmentation formation in brown cotton fiber was possibly the consequence of an interaction between oxidized tannins and glycosylated anthocyanins. Furthermore, our results showed the pigmentation related genes not only regulated the fiber color but also influenced the fiber quality at the fiber elongation stage (10 DPA). The highly expressed flavonoid gene in the fiber elongation stage could be related to the fiber quality. DEGs analyses also revealed that transcript levels of some fiber development genes (Ca2+/CaM, reactive oxygen, ethylene and sucrose phosphate synthase) varied dramatically between these two cotton lines.

Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.The authors resequence a core collection of upland cotton (Gossypium hirsutum) comprising 419 accessions. They analyze genomic variation and conduct a genome-wide association study for 13 fiber quality and yield traits in 12 different environments.

Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits

The ancestors of Gossypium arboreum and Gossypium herbaceum provided the A subgenome for the modern cultivated allotetraploid cotton. Here, we upgraded the G. arboreum genome assembly by integrating different technologies. We resequenced 243 G. arboreum and G. herbaceum accessions to generate a map of genome variations and found that they are equally diverged from Gossypium raimondii. Independent analysis suggested that Chinese G. arboreum originated in South China and was subsequently introduced to the Yangtze and Yellow River regions. Most accessions with domestication-related traits experienced geographic isolation. Genome-wide association study (GWAS) identified 98 significant peak associations for 11 agronomically important traits in G. arboreum. A nonsynonymous substitution (cysteine-to-arginine substitution) of GaKASIII seems to confer substantial fatty acid composition (C16:0 and C16:1) changes in cotton seeds. Resistance to fusarium wilt disease is associated with activation of GaGSTF9 expression. Our work represents a major step toward understanding the evolution of the A genome of cotton.The authors report an improved genome assembly of G. arboretum and resequencing of 243 diploid cotton accessions. GWAS and QTL-seq identify a number of candidate loci that associate with seed oil content, disease resistance and yield traits in cotton.

High-Density Linkage Map Construction and Mapping of Salt-Tolerant QTLs at Seedling Stage in Upland Cotton Using Genotyping by Sequencing (GBS)

Over 6% of agricultural land is affected by salinity. It is becoming obligatory to use saline soils, so growing salt-tolerant plants is a priority. To gain an understanding of the genetic basis of upland cotton tolerance to salinity at seedling stage, an intra-specific cross was developed from CCRI35, tolerant to salinity, as female with Nan Dan (NH), sensitive to salinity, as the male. A genetic map of 5178 SNP markers was developed from 277 F2:3 populations. The map spanned 4768.098 cM, with an average distance of 0.92 cM. A total of 66 QTLs for 10 traits related to salinity were detected in three environments (0, 110, and 150 mM salt treatment). Only 14 QTLs were consistent, accounting for 2.72% to 9.87% of phenotypic variation. Parental contributions were found to be in the ratio of 3:1, 10 QTLs from the sensitive and four QTLs from the resistant parent. Five QTLs were located in At and nine QTLs in the Dt sub-genome. Moreover, eight clusters were identified, in which 12 putative key genes were found to be related to salinity. The GBS-SNPs-based genetic map developed is the first high-density genetic map that has the potential to provide deeper insights into upland cotton salinity tolerance. The 12 key genes found in this study could be used for QTL fine mapping and cloning for further studies.

iTRAQ-Based Comparative Proteomic Analysis of Seedling Leaves of Two Upland Cotton Genotypes Differing in Salt Tolerance

Cotton yields are greatly reduced under high salinity stress conditions, although cotton is considered a moderately salt-tolerant crop. Understanding at the molecular level how cotton responds to salt stress will help in developing salt tolerant varieties. Here, we combined physiological analysis with isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomics of seedling leaves of 2 genotypes differing in salinity tolerance to 200 mM (18.3 dS/m) NaCl stress. Salt stress produced significant stress symptoms in the sensitive genotype Nan Dan Ba Di Da Hua (N), including lower relative water and chlorophyll contents and higher relative electrolyte leakage and Na+/K+ ratio in leaf samples, compared with those in the tolerant genotype Earlistaple 7 (Z). A total of 58 differentially abundant salt-responsive proteins were identified. Asp-Glu-Ala-Asp (DEAD)-box ATP-dependent RNA helicase 3 and protochlorophyllide reductase were markedly suppressed after salt treatment, whereas the phosphate-related differentially abundant proteins (DAPs) phosphoethanolamine N-methyltransferase 1 and 14-3-3-like protein E were induced, and all these proteins may play significant roles in salt stress. Twenty-nine salt-responsive proteins were also genotype specific, and 62.1 and 27.6% of these were related to chloroplast and defense responses, respectively. Based on the Arabidopsis thaliana protein interaction database, orthologs of 25 proteins showed interactions in Arabidopsis, and among these, a calmodulin protein was predicted to have 212 functional partners. In addition, the Golgi apparatus and calcium may be important for salt secretion in cotton. Through integrative proteome and transcriptome analysis, 16 DAPs were matched to differentially expressed genes and verified using qRT-PCR. On the basis of these findings, we proposed that some proteins related to chloroplast, ATP, ribosomal, and phosphate metabolism as well as to the Golgi apparatus and calcium may play key roles in the short-term salt stress response of cotton seedling leaves.

QTL Mapping of Fiber Quality and Yield-Related Traits in an Intra-Specific Upland Cotton Using Genotype by Sequencing (GBS)

Fiber quality and yield improvement are crucial for cotton domestication and breeding. With the transformation in spinning techniques and multiplicity needs, the development of cotton fiber quality and yield is of great importance. A genetic map of 5178 Single Nucleotide Polymorphism (SNP) markers were generated using 277 F2:3 population, from an intra-specific cross between two upland cotton accessions, CCRI35 a high fiber quality as female and Nan Dan Ba Di Da Hua (NH), with good yield properties as male parent. The map spanned 4768.098 cM with an average distance of 0.92 cM. A total of 110 Quantitative Traits Loci (QTLs) were identified for 11 traits, but only 30 QTLs were consistent in at least two environments. The highest percentage of phenotypic variance explained by a single QTL was 15.45%. Two major cluster regions were found, cluster 1 (chromosome17-D03) and cluster 2 (chromosome26-D12). Five candidate genes were identified in the two QTL cluster regions. Based on GO functional annotation, all the genes were highly correlated with fiber development, with functions such as protein kinase and phosphorylation. The five genes were associated with various fiber traits as follows: Gh_D03G0889 linked to qFM-D03_cb, Gh_D12G0093, Gh_D12G0410, Gh_D12G0435 associated with qFS-D12_cb and Gh_D12G0969 linked to qFY-D12_cb. Further structural annotation and fine mapping is needed to determine the specific role played by the five identified genes in fiber quality and yield related pathway.

Comparative transcriptome analysis of TUCPs in Gossypium hirsutum Ligon-lintless-1 mutant and their proposed functions in cotton fiber development.

Transcripts of uncertain coding potential (TUCP) are part of the LncRNAs, which encode some polypeptides. However, the abundance of TUCP transcripts and their roles in Ligon-linless-1 (Li-1) cotton mutant during the early termination of fiber development are still not documented. Li-1 mutant is one of the excellent modules for investigating fiber elongation processes due to its unique fiber developmental stages. To examine the function of TUCP in cotton fiber development, it is important to identify TUCPs and their involvement in fiber development. In this study, we found that 11104 TUCP transcripts were removed by coding potential criteria of Pfam domain scan. Additionally, differential expression levels of TUCP transcripts were detected between Li-1 mutant and the wild-type (WT), which imply their possible functions in cotton fiber development. These results further revealed that a great number of differentially expressed TUCP transcripts in cotton were identified at 8 DPA, followed by 0 DPA and stem. However, these might explain an undesirable function in cotton fiber development. The gene ontology and pathway analysis, based on differential expression patterns of TUCP transcripts on targeted genes, identified the transport process, cytoskeleton structure, membrane permeability and fatty acids. These give new insight into significant involvement in early cessation of cotton fiber development and abnormal stem. The RNA-seq and qRT-PCR expression analyses of TUCP transcripts evidently singled out three possible genes, TUCP_010675, TUCP_001475, TUCP_009444 and other targeted mRNAs. The expression pattern of TUCP transcripts and their mRNA targets provided valuable evidence for further investigations on the biological functions of TUCP in cotton fiber development. The study findings may serve as a useful tool for comparative analysis of TUCP transcripts in cotton species and assist in selection of the applicable candidate genes for further functional analyses, genetic improvement and genetic engineering of cotton fiber development.