Hantao Wang

Identification of favorable SNP alleles and candidate genes for traits related to early maturity via GWAS in upland cotton

BackgroundEarly maturity is one of the most important and complex agronomic traits in upland cotton (Gossypium hirsutum L). To dissect the genetic architecture of this agronomically important trait, a population consisting of 355 upland cotton germplasm accessions was genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) approach, of which a subset of 185 lines representative of the diversity among the accessions was phenotypically characterized for six early maturity traits in four environments. A genome-wide association study (GWAS) was conducted using the generalized linear model (GLM) and mixed linear model (MLM).ResultsA total of 81,675 SNPs in 355 upland cotton accessions were discovered using SLAF-seq and were subsequently used in GWAS. Thirteen significant associations between eight SNP loci and five early maturity traits were successfully identified using the GLM and MLM; two of the 13 associations were common between the models. By computing phenotypic effect values for the associations detected at each locus, 11 highly favorable SNP alleles were identified for five early maturity traits. Moreover, dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between the numbers of highly favorable alleles and the phenotypic values of the target traits were identified. Most importantly, a major locus (rs13562854) on chromosome Dt3 and a potential candidate gene (CotAD_01947) for early maturity were detected.ConclusionsThis study identified highly favorable SNP alleles and candidate genes associated with early maturity traits in upland cotton. The results demonstrate that GWAS is a powerful tool for dissecting complex traits and identifying candidate genes. The highly favorable SNP alleles and candidate genes for early maturity traits identified in this study should be show high potential for improvement of early maturity in future cotton breeding programs.

Detection of Favorable QTL Alleles and Candidate Genes for Lint Percentage by GWAS in Chinese Upland Cotton

Improving cotton yield is a major breeding goal for Chinese upland cotton. Lint percentage is an important yield component and a critical economic index for cotton cultivars, and raising the lint percentage has a close relationship to improving cotton lint yield. To investigate the genetic architecture of lint percentage, a diversity panel consisting of 355 upland cotton accessions was grown, and the lint percentage was measured in four different environments. Genotyping was performed with specific-locus amplified fragment sequencing (SLAF-seq). Twelve single-nucleotide polymorphisms (SNPs) associated with lint percentage were detected via a genome-wide association study (GWAS), in which five SNP loci distributed on chromosomes At3 (A02) and At4 (A08) and contained two major-effect QTLs, which were detected in the best linear unbiased predictions (BLUPs) and in more than three environments simultaneously. Furthermore, favorable haplotypes (FHs) of two major-effect QTLs and 47 putative candidate genes in the two linkage disequilibrium (LD) blocks of these associated loci were identified. The expression levels of these putative candidate genes were estimated using RNA-seq data from ten upland cotton tissues. We found that Gh_A02G1268 was very highly expressed during the early fiber development stage, whereas the gene was poorly expressed in the seed. These results implied that Gh_A02G1268 may determine the lint percentage by regulating seed and fiber development. The favorable QTL alleles and candidate genes for lint percentage identified in this study will have high potential for improving lint yield in future Chinese cotton breeding programs.

Two genomic regions associated with fiber quality traits in Chinese upland cotton under apparent breeding selection.

Fiber quality is one of the most important agronomic traits of cotton, and understanding the genetic basis of its target traits will accelerate improvements to cotton fiber quality. In this study, a panel comprising 355 upland cotton accessions was used to perform genome-wide association studies (GWASs) of five fiber quality traits in four environments. A total of 16, 10 and 7 SNPs were associated with fiber length (FL), fiber strength (FS) and fiber uniformity (FU), respectively, based on the mixed linear model (MLM). Most importantly, two major genomic regions (MGR1 and MGR2) on chromosome Dt7 and four potential candidate genes for FL were identified. Analyzing the geographical distribution of favorable haplotypes (FHs) among these lines revealed that two favorable haplotype frequencies (FHFs) were higher in accessions from low-latitude regions than in accessions from high-latitude regions. However, the genetic diversity of lines from the low-latitude regions was lower than the diversity of lines from the high-latitude regions in China. Furthermore, the FHFs differed among cultivars developed during different breeding periods. These results indicate that FHs have undergone artificial selection during upland cotton breeding in recent decades in China and provide a foundation for the further improvement of fiber quality traits.

Characterization and Functional Analysis of PEBP Family Genes in Upland Cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a naturally occurring photoperiod-sensitive perennial plant species. However, sensitivity to the day length was lost during domestication. The phosphatidylethanolamine-binding protein (PEBP) gene family, of which three subclades have been identified in angiosperms, functions to promote and suppress flowering in photoperiod pathway. Recent evidence indicates that PEBP family genes play an important role in generating mobile flowering signals. We isolated homologues of the PEBP gene family in upland cotton and examined their regulation and function. Nine PEBP-like genes were cloned and phylogenetic analysis indicated the genes belonged to four subclades (FT, MFT, TFL1 and PEBP). Cotton PEBP-like genes showed distinct expression patterns in relation to different cotton genotypes, photoperiod responsive and cultivar maturity. The GhFT gene expression of a semi-wild race of upland cotton were strongly induced under short day condition, whereas the GhPEBP2 gene expression was induced under long days. We also elucidated that GhFT but not GhPEBP2 interacted with FD-like bZIP transcription factor GhFD and promote flowering under both long- and short-day conditions. The present result indicated that GhPEBP-like genes may perform different functions. This work corroborates the involvement of PEBP-like genes in photoperiod response and regulation of flowering time in different cotton genotypes, and contributes to an improved understanding of the function of PEBP-like genes in cotton.

Global analysis of DNA methylation in young (J1) and senescent (J2) Gossypium hirsutum L. cotyledons by MeDIP-Seq

DNA methylation is an important epigenetic modification regulating gene expression, genomic imprinting, transposon silencing and chromatin structure in plants and plays an important role in leaf senescence. However, the DNA methylation pattern during Gossypium hirsutum L. cotyledon senescence is poorly understood. In this study, global DNA methylation patterns were compared between two cotyledon development stages, young (J1) and senescence (J2), using methylated DNA immunoprecipitation (MeDIP-Seq). Methylated cytosine occurred mostly in repeat elements, especially LTR/Gypsy in both J1 and J2. When comparing J1 against J2, there were 1222 down-methylated genes and 623 up-methylated genes. Methylated genes were significantly enriched in carbohydrate metabolism, biosynthesis of other secondary metabolites and amino acid metabolism pathways. The global DNA methylation level decreased from J1 to J2, especially in gene promoters, transcriptional termination regions and regions around CpG islands. We further investigated the expression patterns of 9 DNA methyltransferase-associated genes and 2 DNA demethyltransferase-associated genes from young to senescent cotyledons, which were down-regulated during cotyledon development. In this paper, we first reported that senescent cotton cotyledons exhibited lower DNA methylation levels, primarily due to decreased DNA methyltransferase activity and which also play important role in regulating secondary metabolite process.

An NAM Domain Gene, GhNAC79, Improves Resistance to Drought Stress in Upland Cotton

Plant-specific NAC proteins comprise one of the largest transcription factor families in plants and play important roles in plant development and the stress response. Gossypium hirsutum L. is a major source of fiber, but its growth and productivity are limited by many biotic and abiotic stresses. In this study, the NAC domain gene GhNAC79 was functionally characterized in detail, and according to information about the cotton genome sequences, it was located on scaffold42.1, containing three exons and two introns. Promoter analysis indicated that the GhNAC79 promoter contained both basic and stress-related elements, and it was especially expressed in the cotyledon of Arabidopsis. A transactivation assay in yeast demonstrated that GhNAC79 was a transcription activator, and its activation domain was located at its C-terminus. The results of qRT-PCR proved that GhNAC79 was preferentially expressed at later stages of cotyledon and fiber development, and it showed high sensitivity to ethylene and meJA treatments. Overexpression of GhNAC79 resulted in an early flowering phenotype in Arabidopsis, and it also improved drought tolerance in both Arabidopsis and cotton. Furthermore, VIGS-induced silencing of GhNAC79 in cotton led to a drought-sensitive phenotype. In summary, GhNAC79 positively regulates drought stress, and it also responds to ethylene and meJA treatments, making it a candidate gene for stress studies in cotton.

Identification of the group IIa WRKY subfamily and the functional analysis of GhWRKY17 in upland cotton (Gossypium hirsutum L.)

WRKY transcription factors play important roles in plant defense, stress response, leaf senescence, and plant growth and development. Previous studies have revealed the important roles of the group IIa GhWRKY genes in cotton. To comprehensively analyze the group IIa GhWRKY genes in upland cotton, we identified 15 candidate group IIa GhWRKY genes in the Gossypium hirsutum genome. The phylogenetic tree, intron-exon structure, motif prediction and Ka/Ks analyses indicated that most group IIa GhWRKY genes shared high similarity and conservation and underwent purifying selection during evolution. In addition, we detected the expression patterns of several group IIa GhWRKY genes in individual tissues as well as during leaf senescence using public RNA sequencing data and real-time quantitative PCR. To better understand the functions of group IIa GhWRKYs in cotton, GhWRKY17 (KF669857) was isolated from upland cotton, and its sequence alignment, promoter cis-acting elements and subcellular localization were characterized. Moreover, the over-expression of GhWRKY17 in Arabidopsis up-regulated the senescence-associated genes AtWRKY53, AtSAG12 and AtSAG13, enhancing the plant’s susceptibility to leaf senescence. These findings lay the foundation for further analysis and study of the functions of WRKY genes in cotton.

High-density genetic linkage map construction by F2 populations and QTL analysis of early-maturity traits in upland cotton (Gossypium hirsutum L.)

Due to China’s rapidly increasing population, the total arable land area has dramatically decreased; as a consequence, the competition for farming land allocated for grain and cotton production has become fierce. Therefore, to overcome the existing contradiction between cotton grain and fiber production and the limited farming land, development of early-maturing cultivars is necessary. In this research, a high-density linkage map of upland cotton was constructed using genotyping by sequencing (GBS) to discover single nucleotide polymorphism (SNP) markers associated with early maturity in 170 F2 individuals derived from a cross between LU28 and ZHONG213. The high-density genetic map, which was composed of 3978 SNP markers across the 26 cotton chromosomes, spanned 2480 cM with an average genetic distance of 0.62 cM. Collinearity analysis showed that the genetic map was of high quality and accurate and agreed well with the Gossypium hirsutum reference genome. Based on this high-density linkage map, QTL analysis was performed on cotton early-maturity traits, including FT, FBP, WGP, NFFB, HNFFB and PH. A total 47 QTLs for the six traits were detected; each of these QTLs explained between 2.61% and 32.57% of the observed phenotypic variation. A major region controlling early-maturity traits in Gossypium hirsutum was identified for FT, FBP, WGP, NFFB and HNFFB on chromosome D03. QTL analyses revealed that phenotypic variation explained (PVE) ranged from 10.42% to 32.57%. Two potential candidate genes, Gh_D03G0885 and Gh_D03G0922, were predicted in a stable QTL region and had higher expression levels in the early-maturity variety ZHONG213 than in the late-maturity variety LU28. However, further evidence is required for functional validation. This study could provide useful information for the dissection of early-maturity traits and guide valuable genetic loci for molecular-assisted selection (MAS) in cotton breeding.

Proteomic Analysis of Differences in Fiber Development between Wild and Cultivated Gossypium hirsutum L.

Upland cotton (Gossypium hirsutum L.) is one of the worlds most important fiber crops, accounting for more than 90% of all cotton production. While their wild progenitors have relatively short and coarse, often tan-colored fibers, modern cotton cultivars possess longer, finer, stronger, and whiter fiber. In this study, the wild and cultivated cottons (YU-3 and TM-1) selected show significant differences on fibers at 10 days postanthesis (DPA), 20 DPA, and mature stages at the morphological level. To explore the effects of domestication, reveal molecular mechanisms underlying these phenotypic differences, and better inform our efforts to further enhance cotton fiber quality, isobaric tags for relative and absolute protein quantification-facilitated proteomic methods were performed on developing fibers. There were 6990 proteins identified; among them, 336 were defined as differentially expressed proteins between fibers of wild versus domesticated cotton. The down- or up-regulated proteins in wild cotton were involved in phenylpropanoid biosynthesis, zeatin biosynthesis, fatty acid elongation, and other processes. Association analysis between transcriptome and proteome showed positive correlations between transcripts and proteins at both 10 DPA and 20 DPA. Differences in proteomics have been verified at the mRNA level by quantitative real-time polymerase chain reaction and have been validated at the physiological and biochemical levels by POD (peroxidase) activity assays and ZA (zeatin) content estimates. This work corroborates the major pathways involved in cotton fiber development and demonstrates that POD activity and zeatin content have a great potential related to fiber elongation and thickening.

Characterization and functional analysis of GhWRKY42, a group IId WRKY gene, in upland cotton (Gossypium hirsutum L.)

BackgroundWRKY transcription factors (TFs) participate in various physiological processes of plants. Although WRKY genes have been well studied in model plants, knowledge of the functional roles of these genes is still extremely limited in cotton.ResultsIn this study, a group IId WRKY gene from cotton, GhWRKY42, was isolated and characterized. Our data showed that GhWRKY42 localized to the nucleus. A transactivation assay in yeast demonstrated that GhWRKY42 was not a transcriptional activator. A β-glucuronidase (GUS) activity assay revealed that the promoter of GhWRKY42 showed fragment deletion activity in Nicotiana tabacum and was mainly expressed in the roots, stems and leaves of ProGhWRKY42::GUS transgenic Arabidopsis plants. Quantitative real-time PCR (qRT-PCR) analysis indicated that GhWRKY42 was up-regulated during leaf senescence and was induced after exposure to abiotic stresses. Constitutive expression of GhWRKY42 in Arabidopsis led to a premature aging phenotype, which was correlated with an increased number of senescent leaves, reduced chlorophyll content and elevated expression of senescence-associated genes (SAGs). In addition, virus-induced gene silencing (VIGS) was used to silence the endogenous GhWRKY42 gene in cotton, and this silencing reduced plant height.ConclusionsOur findings indicate that GhWRKY42 is involved in abiotic stress responses, premature leaf senescence and stem development. This work establishes a solid foundation for further functional analysis of the GhWRKY42 gene in cotton.