Series-temporal transcriptome profiling of cotton reveals the response mechanism of phosphatidylinositol signaling system in the early stage of drought stress.

Plants are sessile organisms suffering severe environmental conditions. Drought stress is one of the major environmental issues that affect plant growth and productivity. Although complex regulatory gene networks of plants under drought stress have been analyzed extensively, the response mechanism in the early stage of drought stress is still rarely mentioned. Here, we performed transcriptome analyses on cotton samples treated for a short time (10 min, 30 min, 60 min, 180 min) using 10% PEG, which is used to simulate drought stress. The analysis of differently expressed genes (DEGs) showed that the number of DEGs in roots was obviously more than that in stems and leaves at the four time points and maintained >2000 FDEGs (DEGs appearing for the first time) from 10 min, indicating that root tissues of plants respond to drought stress quickly and continuously strongly. Gene ontology (GO) analysis showed that DEGs in roots were mainly enriched in protein modification and microtubule-based process. DEGs were found significantly enriched in phosphatidylinositol signaling system at 10 min through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, implying the great importance of phosphatidylinositol signal in the early stage of drought stress. What was more, two co-expression modules, which were significantly positively correlated with drought stress, were found by Weighted Gene Co-expression Network Analysis (WGCNA). From one of the co-expression modules, we identified a hub-gene Gohir.A07G058200, which is annotated as "phosphatidylinositol 3- and 4-kinase" in phosphatidylinositol signaling system, and found this gene may interact with auxin-responsive protein. This result suggested that Gohir.A07G058200 may be involved in the crosstalk of phosphatidylinositol signal and auxin signal in the early stage of drought stress. In summary, through transcriptome sequencing, we found that phosphatidylinositol signaling system is an important signal transduction pathway in early stage in response to drought stress, and it may interact with auxin signal transduction through phosphatidylinositol 3- and 4-kinase.

The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in Upland cotton.

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.

A genome-wide analysis of pentatricopeptide repeat (PPR) protein-encoding genes in four Gossypium species with an emphasis on their expression in floral buds, ovules, and fibers in upland cotton.

Cotton is the most important natural fiber used in textiles. Breeding for "three-lines", i.e., cytoplasmic male sterility (CMS)-based sterile (A), maintainer (B), and restorer (R) line, is a promising approach to harness hybrid vigor in cotton. Pentatricopeptide repeat (PPR) protein-encoding genes play an important role in plant growth and development including restoration of CMS plants to male fertility. However, PPRs, especially those contributing to CMS and fiber development, remain largely unknown in cotton. In this study, a genome-wide identification and characterization of PPR gene family in four Gossypium species with genome sequences (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) were performed, and expressed PPR genes in developing floral buds, ovules, and fibers were compared to identify possible PPRs related to CMS restoration and fiber development. A total of 539, 558, 1032, and 1055 PPRs were predicted in the above four species, respectively, which were further mapped to chromosomes for a synteny analysis. Through an RNA-seq analysis, 86% (882) PPRs were expressed in flowering buds of upland cotton (G. hirsutum); however, only 11 and 6 were differentially expressed (DE) between restorer R and its near-isogenic (NI) B and between R and its NI A line, respectively. Another RNA-seq analysis identified the expression of only 54% (556) PPRs in 0 and 3 day(s) post-anthesis (DPA) ovules and 24% (247) PPRs in 10 DPA fibers; however, only 59, 6, and 27 PPRs were DE in 0 and 3 DPA ovules, and 10 DPA fibers between two backcross inbred lines (BILs) with differing fiber length, respectively. Only 2 PPRs were DE between Xuzhou 142 and its fiberless and fuzzless mutant. Quantitative RT-PCR analysis confirmed the validity of the RNA-seq results for the gene expression pattern. Therefore, only a very small number of PPRs may be associated with fertility restoration of CMS and genetic differences in fiber initiation and elongation. These results lay a foundation for understanding the roles of PPR genes in cotton, and will be useful in the prioritization of candidate PPR gene functional validation for cotton CMS restoration and fiber development.

Expression profiles of a cytoplasmic male sterile line of Gossypium harknessii and its fertility restorer and maintainer lines revealed by RNA-Seq.

The Gossypium harknessii background cytoplasmic male sterility (CMS) system has been used in cotton hybrid breeding in China. However, the mechanism underlying pollen abortion and fertility restoration in CMS remains to be determined. In this study, we used RNA-seq to identify critical genes and pathways associated with CMS in G. harknessii based CMS lines (588A), the near isogenic restorer lines (588R), and maintainer lines (588B). We performed an assembly of 80,811,676 raw reads into 89,939 high-quality unigenes with an average length of 698 bp. Among these, 72.62% unigenes were annotated in public protein databases and were classified into functional clusters. In addition, we investigated the changes in expression of genes between 588A and 588B (588R); the RNA-seq data showed 742 differentially expressed genes (DEGs) between 588A and 588B and 748 DEGs between 588A and 588R. They were mainly down-regulated in 588A and most of them distributed in metabolic and biosynthesis of secondary metabolites pathways. Further analysis revealed 23 pollen development related genes were differentially expressed between 588A and 588B. Numerous genes associated with tapetum development were down-regulated in 588A, implicating tapetum dysplasia may be a key reason for pollen abortion in CMS lines. Also, among DEGs between 588A and 588R, we identified two PPR genes which were highly up-regulated in restorer line. This study may provide assistance for detailed molecular analysis and a better understanding of harknessii based CMS in cotton.

Cloning and Expression Analysis of Eight Upland Cotton Pentatricopeptide Repeat Family Genes.

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in plants. Most PPR genes are localized in mitochondria and chloroplasts functioning in regulation of plant growth and development, fertility restoration for cytoplasmic male sterility (CMS), and stress defense. In this study, using in silico cloning and PCR amplification with degenerate primers based on Arabidopsis PPR genes, we cloned eight new full-length PPR genes encoding protein sequences ranging from 458 to 875 amino acids, with 8 to 16 repetitive PPR elements in upland cotton and all of them lack introns. Expression analysis revealed that eight PPR genes were differently expressed in roots, stems, leaves, and floral buds. As for GhI12, its expression in floral buds at days 3-5 was significantly higher in line 777R (restorer line) than in line 777A (CMS line). Further tests with real-time PCR showed that GhI12 expression peaked at day 3 in 777R, followed by a gradual decline, while its expression fluctuated in 777A, peaking at day 5 and day 13. In addition, Gh155c17 and GhI12 were upregulated under salt stress. This is the first report of upland cotton PPR genes involved in salt stress response.

Overexpression of three TaEXPA1 homoeologous genes with distinct expression divergence in hexaploid wheat exhibit functional retention in Arabidopsis.

Common wheat is a hexaploid species with most of the genes present as triplicate homoeologs. Expression divergences of homoeologs are frequently observed in wheat as well as in other polyploid plants. However, little is known about functional variances among homologous genes arising from polyploidy. Expansins play diverse roles in plant developmental processes related to the action of cell wall loosening. Expression of the three TaEXPA1 homoeologs varied dynamically at different stages and organs, and epigenetic modifications contribute to the expression divergence of three TaEXPA1 homoeologs during wheat development. Nevertheless, their functions remain to be clarified. We found that over expression of TaEXPA1-A, -B and -D produced similar morphological changes in transgenic Arabidopsis plants, including increased germination and growth rate during seedling and adult stages, indicating that the proteins encoded by these three wheat TaEXPA1 homoeologs have similar (or conserved) functions in Arabidopsis. Collectively, our present study provided an example of a set of homoeologous genes expression divergence in different developmental stages and organs in hexaploid wheat but functional retention in transgenic Arabidopsis plants.

Epigenetic modification contributes to the expression divergence of three TaEXPA1 homoeologs in hexaploid wheat (Triticum aestivum).

Common wheat is a hexaploid species with most of the genes present as triplicate homoeologs. Expression divergences of homoeologs are frequently observed in wheat, as well as in other polyploid plants. However, the mechanisms underlying this phenomenon are poorly understood. Expansin genes play important roles in the regulation of cell size, as well as organ size. We found that all three TaEXPA1 homoeologs were silenced in seedling roots. In seedling leaves, TaEXPA1-A and TaEXPA1-D were expressed, but TaEXPA1-B was silenced. Further analysis revealed that silencing of TaEXPA1-B in leaves occurred after the formation of the hexaploid. Chromatin immunoprecipitation assays revealed that the transcriptional silencing of three TaEXPA1 homoeologs in roots was correlated with an increased level of H3K9 dimethylation and decreased levels of H3K4 trimethylation and H3K9 acetylation. Reactivation of TaEXPA1-A and TaEXPA1-D expression in leaves was correlated with increased levels of H3K4 trimethylation and H3K9 acetylation, and decreased levels of H3K9 dimethylation in their promoters, respectively. Moreover, a higher level of cytosine methylation was detected in the promoter region of TaEXPA1-B, which may contribute to its silencing in leaves. We demonstrated that epigenetic modifications contribute to the expression divergence of three TaEXPA1 homoeologs during wheat development.

Ectopic overexpression of wheat adenosine diphosphate-ribosylation factor, TaARF, increases growth rate in Arabidopsis.

Differential gene expression between hybrids and their parents is considered to be associated with heterosis. However, the physiological functions and possible contribution to heterosis of these differentially expressed genes are unknown. We have isolated one hybrid upregulated gene encoding putative wheat ADP-ribosylation factor, designated TaARF. In this study, real-time quantitative reverse transcription-polymerase chain reaction analysis indicated that the TaARF transcript was preferentially expressed in root, node and crown, and the accumulation of TaARF mRNA in hybrid was more than 1.5-fold higher than that in two parents. In order to understand possible roles of the putative wheat ARF gene, TaARF was overexpressed in Arabidopsis, and the transgenic plants were characterized. We show that ectopic overexpression of TaARF in Arabidopsis leads to increased leaf area, increased growth rate and earlier transition to flowering, suggesting that TaARF plays significant roles in growth and development. This study provides evidence demonstrating that TaARF plays important roles in growth and development and we speculate that the upregulated expression of this gene might contribute to the heterosis observed in wheat root and leaf growth.

Identification of differentially expressed genes in leaf and root between wheat hybrid and its parental inbreds using PCR-based cDNA subtraction.

Heterosis was defined as the advantage of hybrid performance over its parents in terms of growth and productivity. Previous studies showed that differential gene expression between hybrids and their parents is responsible for the heterosis; however, information on systematic identification and characterization of the differentially expressed genes are limited. In this study, an interspecific hybrid between common wheat (Triticum aestivum. L., 2n = 6x = 42, AABBDD) line 3338 and spelt (Triticum spelta L. 2n = 6x = 42, AABBDD) line 2463 was found to be highly heterotic in both aerial growth and root related traits, and was then used for expression assay. A modified suppression subtractive hybridization (SSH) was used to generate four subtracted cDNA libraries, and 748 nonreduandant cDNAs were obtained, among which 465 had high sequence similarity to the GenBank entries and represent diverse of functional categories, such as metabolism, cell growth and maintenance, signal transduction, photosynthesis, response to stress, transcription regulation and others. The expression patterns of 68.2% SSH-derived cDNAs were confirmed by reverse Northern blot, and semi-quantitative RT-PCR exhibited the similar results (72.2%). And it was concluded that the genes differentially expressed between hybrids and their parents involved in diverse physiological process pathway, which might be responsible for the observed heterosis.

[Comparative analysis of genetic diversity revealed by genomic-sSR, EST-SSR and pedigree in wheat (Triticum asetivum L.)].

Two SSR molecular markers, genomic-SSR and EST (expressed sequence tagged)-SSR, were used to measure the genetic diversity among 18 accessions of common wheat with known pedigrees, which were collected from winter wheat production region in Northern China. In addtion, the genetic diversity revealed by pedigree, EST-SSR and genomic-SSR was also compared. The results showed that the average number of alleles per genomic-SSR locus is 3.34, which is higher than that of EST-SSR (2.31), indicating that genomic-SSR is more polymorphic than EST-SSR. Genomic-SSR and EST-SSR were used to calculate the genetic distance (GD) in different materials. The mean GD value of EST-SSR for the 18 wheat genotypes is 0.3996,which is lower than that of genomic-SSR (0.5458). At the locus level, the mean GD calculated using pedigree information is higher (0. 9716) than that of genomic-SSR and EST-SSR markers (0.5458 and 0.3996). Therefore, although polymorphism of EST-SSR is low as compared to genomic SSR,it provides a more accurate evaluation of genetic relationship, especially when accessions are very closely related in pedigree. The strategy for improving the genetic diversity of wheat was also discussed.