Subgenome Bias and Temporal Postponement of Gene Expression Contributes to the Distinctions of Fiber Quality in Gossypium Species

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

Allelic expression of AhNSP2-B07 due to parent of origin affects peanut nodulation

Legumes are well-known for establishing a symbiotic relationship with rhizobia in root nodules to fix nitrogen from the atmosphere. The nodulation signaling pathway 2 (NSP2) gene plays a critical role in the symbiotic signaling pathway. In cultivated peanut, an allotetraploid (2n = 4x = 40, AABB) legume crop, natural polymorphisms in a pair of NSP2 homoeologs (N_a and N_b) located on chromosomes A08 and B07, respectively, can cause loss of nodulation. Interestingly, some heterozygous (N_Bn_b) progeny produced nodules, while some others do not, suggesting non-Mendelian inheritance in the segregating population at the N_b locus. In this study, we investigated the non-Mendelian inheritance at the N_B locus. Selfing populations were developed to validate the genotypical and phenotypical segregating ratios. Allelic expression was detected in roots, ovaries, and pollens of heterozygous plants. Bisulfite PCR and sequencing of the N_b gene in gametic tissue were performed to detect the DNA methylation variations of this gene in different gametic tissues. The results showed that only one allele at the Nb locus expressed in peanut roots during symbiosis. In the heterozygous (N_bn_b) plants, if dominant allele expressed, the plants produced nodules, if recessive allele expressed, then no nodules were produced. qRT-PCR experiments revealed that the expression of N_b gene in the ovary was extremely low, about seven times lower than that in pollen, regardless of genotypes or phenotypes of the plants at this locus. The results indicated that N_b gene expression in peanut depends on the parent of origin and is imprinted in female gametes. However, no significant differences of DNA methylation level were detected between these two gametic tissues by bisulfite PCR and sequencing. The results suggested that the remarkable low expression of N_b in female gametes may not be caused by DNA methylation. This study provided a unique genetic basis of a key gene involved in peanut symbiosis, which could facilitate understanding the regulation of gene expression in symbiosis in polyploid legumes.

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

Biased gene expression reveals the contribution of subgenome to altitude adaptation in allopolyploid Isoetes sinensis

Allopolyploids are believed to inherit the genetic characteristics of its progenitors and exhibit stronger adaptability and vigor. The allotetraploid Isoetes sinensis was formed by the natural hybridization and polyploidization of two diploid progenitors, Isoetes taiwanensis and Isoetes yunguiensis, and was believed to have the potential to adapt to plateau environments. To explore the expression pattern of homoeologous genes and their contributions to altitude adaptation, we transplanted natural allotetraploid I. sinensis (TnTnYnYn) along the altitude gradient for a long-term, and harvested them in summer and winter, respectively. One year after transplanting, it still lived well, even in the extreme environment of the Qinghai-Tibet Plateau. Then, we performed high-throughput RNA sequencing to measure their gene expression level. A total of 7801 homoeologous genes were expressed, among which 5786 were identified as shared expression in different altitudes and seasons. We further found that altitude variations could change the subgenome bias trend of I. sinensis, but season could not. Moreover, the functions of uniquely expressed genes indicated that temperature might be an important restrictive factor during the adaptation process. Through the analysis of DEGs and uniquely expressed genes, we found that Y subgenome provided more contributions to high altitude adaptation than T subgenome. These adaptive traits to high altitude may be inherited from its plateau progenitor I. yunguiensis. Through weighted gene co-expression network analysis, pentatricopeptide repeats gene family and glycerophospholipid metabolism pathway were considered to play important roles in high-altitude adaptation. Totally, this study will enrich our understanding of allopolyploid in environmental adaptation.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.