Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant

Genome-wide identification and expression analysis of the WRKY gene family in Rhododendron henanense subsp. lingbaoense

Background

This work explored the characteristics of the WRKY transcription factor family in Rhododendron henanense subsp. lingbaoense (Rhl) and the expression patterns of these genes under abiotic stress by conducting bioinformatics and expression analyses.

Methods

RhlWRKY genes were identified from a gene library of Rhl. Various aspects of these genes were analyzed, including genetic structures, conserved sequences, physicochemical properties, cis-acting elements, and chromosomal location. RNA-seq was employed to analyze gene expression in five different tissues of Rhl: roots, stems, leaves, flowers, and hypocotyls. Additionally, qRT-PCR was used to detect changes in the expression of five RhlWRKY genes under abiotic stress.

Result

A total of 65 RhlWRKY genes were identified and categorized into three subfamilies based on their structural characteristics: Groups I, II, and III. Group II was further divided into five subtribes, with shared similar genetic structures and conserved motifs among members of the same subtribe. The physicochemical properties of these proteins varied, but the proteins are generally predicted to be hydrophilic. Most proteins are predicted to be in the cell nucleus, and distributed across 12 chromosomes. A total of 84 cis-acting elements were discovered, with many related to responses to biotic stress. Among the identified RhlWRKY genes, there were eight tandem duplicates and 97 segmental duplicates. The majority of duplicate gene pairs exhibited Ka/Ks values <1, indicating purification under environmental pressure. GO annotation analysis indicated that WRKY genes regulate biological processes and participate in a variety of molecular functions. Transcriptome data revealed varying expression levels of 66.15% of WRKY family genes in all five tissue types (roots, stems, leaves, flowers, and hypocotyls). Five RhlWRKY genes were selected for further characterization and there were changes in expression levels for these genes in response to various stresses.

Conclusion

The analysis identified 65 RhlWRKY genes, among which the expression of WRKY_42 and WRKY_17 were mainly modulated by the drought and MeJA, and WRKY_19 was regulated by the low-temperature and high-salinity conditions. This insight into the potential functions of certain genes contributes to understanding the growth regulatory capabilities of Rhl.

REPS2 downregulation facilitates FGF-induced adhesion and migration in human lens epithelial cells through FAK/Cdc42 signaling and contributes to posterior capsule opacification

Posterior capsular opacification (PCO) can cause postoperative visual loss after cataract surgery. Residual human lens epithelial cell (HLEC) proliferation, migration, epithelial-mesenchymal transition (EMT) and synthesis of extracellular matrix (ECM) are the entitative reasons for PCO. Low expression of Ral-binding protein 1-associated Eps domain-containing 2 (REPS2) and high levels of basic fibroblast growth factor (b-FGF) were observed in the lens and postoperative aqueous humor of cataract patients. REPS2 was identified as a negative regulator in growth factor signaling; however, its function in HLECs is unknown. This was first investigated in the present study by evaluating REPS2 expression in anterior lens capsules from cataract patients, a mouse cataract model, and HLE-b3 cells. The biological function of REPS2 in HLE-B3 cells was assessed by REPS2 silencing and Cell Counting Kit 8, wound healing, Transwell migration, F-actin staining, G-protein pulldown and western blot assays. In the present study, REPS2 was significantly downregulated in human and mouse cataract capsules and H₂O₂-treated HLE-B3 cells. REPS2 knockdown increased fibronectin, type I collagen, and α-smooth muscle actin expression levels and stimulated HLECs proliferation and migration; these effects were enhanced by FGF treatment and accompanied with focal adhesion kinase (FAK) phosphorylation, cell division cycle 42 (Cdc42) activation, focal adhesion protein upregulation, and F-actin cytoskeleton reorganization. However, treatment with the FAK inhibitor PF573228 abolished these effects. Thus, REPS2 downregulation in cataract HLECs induces their proliferation and facilitates FGF-induced ECM synthesis, EMT, cell adhesion and migration by activating FAK/Cdc42 signaling, which may underlie PCO pathogenesis.

Molecular Regulation of Cotton Fiber Development: A Review

Cotton (Gossypium spp.) is an economically important natural fiber crop. The quality of cotton fiber has a substantial effect on the quality of cotton textiles. The identification of cotton fiber development-related genes and exploration of their biological functions will not only enhance our understanding of the elongation and developmental mechanisms of cotton fibers but also provide insights that could aid the cultivation of new cotton varieties with improved fiber quality. Cotton fibers are single cells that have been differentiated from the ovule epidermis and serve as a model system for research on single-cell differentiation, growth, and fiber production. Genes and fiber formation mechanisms are examined in this review to shed new light on how important phytohormones, transcription factors, proteins, and genes linked to fiber development work together. Plant hormones, which occur in low quantities, play a critically important role in regulating cotton fiber development. Here, we review recent research that has greatly contributed to our understanding of the roles of different phytohormones in fiber development and regulation. We discuss the mechanisms by which phytohormones regulate the initiation and elongation of fiber cells in cotton, as well as the identification of genes involved in hormone biosynthetic and signaling pathways that regulate the initiation, elongation, and development of cotton fibers.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Comparative Metabolomics Analysis Reveals Sterols and Sphingolipids Play a Role in Cotton Fiber Cell Initiation

Cotton fiber is a seed trichome that protrudes from the outer epidermis of cotton ovule on the day of anthesis (0 day past anthesis, 0 DPA). The initial number and timing of fiber cells are closely related to fiber yield and quality. However, the mechanism underlying fiber initiation is still unclear. Here, we detected and compared the contents and compositions of sphingolipids and sterols in 0 DPA ovules of Xuzhou142 lintless-fuzzless mutants (Xufl) and Xinxiangxiaoji lintless-fuzzless mutants (Xinfl) and upland cotton wild-type Xuzhou142 (XuFL). Nine classes of sphingolipids and sixty-six sphingolipid molecular species were detected in wild-type and mutants. Compared with the wild type, the contents of Sphingosine-1-phosphate (S1P), Sphingosine (Sph), Glucosylceramide (GluCer), and Glycosyl-inositol-phospho-ceramides (GIPC) were decreased in the mutants, while the contents of Ceramide (Cer) were increased. Detail, the contents of two Cer molecular species, d18:1/22:0 and d18:1/24:0, and two Phyto-Cer molecular species, t18:0/22:0 and t18:0/h22:1 were significantly increased, while the contents of all GluCer and GIPC molecular species were decreased. Consistent with this result, the expression levels of seven genes involved in GluCer and GIPC synthesis were decreased in the mutants. Furthermore, exogenous application of a specific inhibitor of GluCer synthase, PDMP (1-phenyl-2-decanoylamino-3-morpholino-1-propanol), in ovule culture system, significantly inhibited the initiation of cotton fiber cells. In addition, five sterols and four sterol esters were detected in wild-type and mutant ovules. Compared with the wild type, the contents of total sterol were not significantly changed. While the contents of stigmasterol and campesterol were significantly increased, the contents of cholesterol were significantly decreased, and the contents of total sterol esters were significantly increased. In particular, the contents of campesterol esters and stigmasterol esters increased significantly in the two mutants. Consistently, the expression levels of some sterol synthase genes and sterol ester synthase genes were also changed in the two mutants. These results suggested that sphingolipids and sterols might have some roles in the initiation of fiber cells. Our results provided a novel insight into the regulatory mechanism of fiber cell initiation.

Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression

Cotton (Gossypium L.) is the most important fiber crop worldwide. Here, transcriptome analysis was conducted on developing fibers of a G. mustelinum introgression line, IL9, and its recurrent parent, PD94042, at 17 and 21 days post-anthesis (dpa). Differentially expressed genes (DEGs) of PD94042 and IL9 were identified. Gene Ontology (GO) enrichment analysis showed that the annotated DEGs were rich in two main biological processes and two main molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis likewise showed that the annotated DEGs were mainly enriched in metabolic pathways and biosynthesis of secondary metabolites. In total, 52 DEGs were selected as candidate genes based on comparison of the DEGs and GO function annotation information. Quantitative real-time PCR (RT-qPCR) analysis results for 12 randomly selected DEGs were consistent with transcriptome analysis. SNP identification based on G. mustelinum chromatin segment introgression showed that 394 SNPs were identified in 268 DEGs, and two genes with known functions were identified within fiber strength quantitative trait loci (QTL) regions or near the confidence intervals. We identified 52 key genes potentially related to high fiber strength in a G. mustelinum introgression line and provided significant insights into the study of cotton fiber quality improvement.

Genome-wide identification and expression analysis of the plant specific LIM genes in Gossypium arboreum under phytohormone, salt and pathogen stress

Asiatic cotton (Gossypium arboreum) cultivated as ‘desi cotton’ in India, is renowned for its climate resilience and robustness against biotic and abiotic stresses. The genome of G. arboreum is therefore, considered as a valued reserve of information for discovering novel genes or gene functions for trait improvements in the present context of cotton cultivation world-wide. In the present study, we carried out genome-wide analysis of LIM gene family in desi cotton and identified twenty LIM domain proteins (GaLIMs) which include sixteen animals CRP-like GaLIMs and four plant specific GaLIMs with presence (GaDA1) or absence (GaDAR) of UIM (Ubiquitin Interacting Motifs). Among the sixteen CRP-like GaLIMs, eleven had two conventional LIM domains while, five had single LIM domain which was not reported in LIM gene family of the plant species studied, except in Brassica rapa. Phylogenetic analysis of these twenty GaLIM proteins in comparison with LIMs of Arabidopsis, chickpea and poplar categorized them into distinct αLIM1, βLIM1, γLIM2, δLIM2 groups in CRP-like LIMs, and GaDA1 and GaDAR in plant specific LIMs group. Domain analysis had revealed consensus [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₂-C-X₁₇-C-X₂-H)] and [(C-X₂-C-X₁₇-H-X₂-C)-X₂-(C-X₄-C-X₁₅-C-X₂-H)] being conserved as first and/or second LIM domains of animal CRP-like GaLIMs, respectively. Interestingly, single LIM domain containing GaLIM15 was found to contain unique consensus with longer inter-zinc-motif spacer but shorter second zinc finger motif. All twenty GaLIMs showed variable spatio-temporal expression patterns and accordingly further categorized into distinct groups of αLIM1, βLIM1, γLIM2 δLIM2 and plant specific LIM (DA1/DAR). For the first time, response of GaDA1/DAR under the influence of biotic and abiotic stresses were studied in cotton, involving treatments with phytohormones (Jasmonic acid and Abscisic acid), salt (NaCl) and wilt causing pathogen (Fusarium oxysporum). Expressions patterns of GaDA1/DAR showed variable response and identified GaDA2 as a probable candidate gene for stress tolerance in G. arboreum.

Identification, characterization and expression profiling of circular RNAs in the early cotton fiber developmental stages

Circular RNA is one of the endogenous non-coding RNAs with a covalently closed loop structure and largely involved in regulating gene expression. However, the abundance of circular RNAs and their regulatory functions during the early stages of fiber development are still not known. In this work, we conducted high-throughput sequencing of the Ligonlintless-1 and its wild-type at 0 DPA, 8 DPA and stem. A total of 2811 circular RNAs were identified and unevenly distributed across cotton chromosomes. We found 34, 142 and 27 circular RNAs were differentially expressed between Ligonlintless-1 and wild-type at 0 DPA, 8 DPA and stem, respectively. Both circular RNA-microRNA-mRNA network and MeJA treatment results in Ligonlintless-1 and wild-type might provide a strong indication of four circular RNAs and ghr_miR169b being important biological molecular associating with fiber development. The results provide new insight into the putative molecular function of circular RNAs in the regulation of fiber development.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.

Salt-Free Dyeing of Modified Cotton through Graft Polymerization with Highly Enhanced Dye Fixation and Good Strength Properties

Modification of cotton fabric with 2-methacryloyloxyethyltrimethyl ammonium chloride (DMC) was achieved through free-radical initiated graft polymerization with K₂S₂O₈/NaHSO₃ as the initiator. Grafting of DMC was confirmed by ATR-IR of the modified cotton. The optimal grafting reaction conditions, including DMC dosage, mole ratio of initiator to DMC, temperature, and time, were determined by cation content and dye fixation results of the modified cotton. The modified fibers were characterized by X-ray diffraction (XRD), scanning electron microscope (SEM), and whiteness measurement. Salt-free dyeing of the modified cotton with commonly used C. I. Reactive Blue 19, C. I. Reactive Yellow 145, and C. I. Reactive Red 195 presented high fixation of 96.8%, 98.7%, and 97.3%, respectively. These results indicated that the modification is effective for changing the surface charge of the fiber and increasing the dye-fiber reactivity. The color fastness and strength property were still very satisfactory. With excellent properties, this dyeing method shows promise in real application for eliminating the usage of salt and reducing environmental pollution.

Genetic Analysis of the Transition from Wild to Domesticated Cotton (Gossypium hirsutum L.)

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F₂ individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Proteomic profiling of cotton fiber developmental transition from cell elongation to secondary wall deposition

Cotton fiber developmental transition from elongation to secondary cell wall biosynthesis is a critical growth shifting phase that affects fiber final length, strength, and other properties. Morphological dynamic analysis indicated an asynchronous fiber developmental pattern between two most important commercial cotton species, Gossypium hirsutum (Gh) and G. barbadense (Gb). Using isobaric tags for relative and absolute quantitation techniques, we examined the temporal changes of protein expression at three representative development periods (15-19, 19-23, and 23-27 dpa) in both species. Strikingly, a large proportion of differentially expressed proteins (DEPs) were identified at 19-23 dpa in Gh and at 23-27 dpa in Gb, corresponding to their fiber developmental transition timing. To better understand fiber transitional development, we comparatively analyzed those DEPs in 19-23 dpa of Gh vs. in 23-27 dpa of Gb, and noted that these cotton species indeed share fundamentally similar fiber developmental features under the biological processes. We also showed that there are limited overlaps in both specific upregulated and downregulated proteins between the two species, suggesting species-specific protein regulations in the development process. Proteomic profiling results revealed dynamic changes of several key proteins and biological processes that are potentially correlated with fiber developmental transition. During the transition, upregulated proteins are mainly involved in carbohydrate/energy metabolism, oxidation-reduction, cytoskeleton, protein turnover, Ca2+ signaling, etc., whereas important downregulated proteins are mostly involved in phenylpropanoid and flavonoid secondary metabolism pathways. The gene expressions of several changed proteins in this key stage were also examined by quantitative reverse transcription polymerase chain reaction. Overall, the present study provides accurate pictures of the regulatory networks of functional proteins during the fiber developmental transition, therefore highlighting candidate genes/proteins and related pathways for the cotton fiber improvement.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.

Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.

Identification of Histone H3 (HH3) Genes in Gossypium hirsutum Revealed Diverse Expression During Ovule Development and Stress Responses

Histone acts as the core for nucleosomes and is a key protein component of chromatin. Among different histone variants, histone H3 (HH3) variants have been reported to play vital roles in plant development. However, biological information and evolutionary relationships of HH3 genes in cotton remain to be elucidated. The current study identified 34 HH3 genes in Gossypium hirsutum. Phylogenetic analysis classified HH3 genes of 19 plant species into eight distinct clades. Sequence logos analysis among Arabidopsis, rice, and G. hirsutum amino acid residues showed higher conservation in amino acids. Using collinearity analysis, we identified 81 orthologous/paralogous gene pairs among the four genomes (A, D, At, and Dt) of cotton. Further, orthologous/paralogous and the Ka/Ks ratio demonstrated that cotton HH3 genes experienced strong purifying selection pressure with restricted functional divergence resulting from segmental and whole genome duplication. Expression pattern analysis indicated that GhHH3 genes were preferentially expressed in cotton ovule tissues. Additionally, GhHH3 gene expression can be regulated by abiotic stresses (cold, heat, sodium chloride (NaCl), and polyethylene glycol (PEG)) and phytohormonal (brassinolide (BL), gibberellic acid (GA), indole-3-acetic acid (IAA), salicylic acid (SA), and methyl jasmonate (MeJA)) treatments, suggesting that GhHH3 genes might play roles in abiotic and hormone stress resistance. Taken together, this work provides important information to decipher complete molecular and physiological functions of HH3 genes in cotton.

Role of RpoN from Labrenzia aggregata LZB033 (Rhodobacteraceae) in Formation of Flagella and Biofilms, Motility, and Environmental Adaptation

Labrenzia aggregata LZB033 (Rhodobacteraceae), which produces dimethylsulfoniopropionate (DMSP) and reduces nitrate to nitrogen, was isolated from seawater of the East China Sea. Its genome encodes a large number of transcriptional regulators which may be important for its adaptation to diverse marine environments. The alternative σ⁵⁴ factor (RpoN) is a central regulator of many bacteria, regulating the transcription of multiple genes and controlling important cellular functions. However, the exact role of RpoN in Labrenzia spp. is unknown. In this study, an in-frame rpoN deletion mutant was constructed in LZB033, and the function of RpoN was determined. To systematically identify RpoN-controlled genes, we performed a detailed analysis of gene expression differences between the wild-type strain and the ΔrpoN mutant using RNA sequencing. The expression of 175 genes was shown to be controlled by RpoN. Subsequent phenotypic assays showed that the ΔrpoN mutant was attenuated in flagellar biosynthesis and swimming motility, utilized up to 13 carbon substrates differently, lacked the ability to assimilate malic acid, and displayed markedly decreased biofilm formation. In addition, stress response assays showed that the ΔrpoN mutant was impaired in the ability to survive under different challenge conditions, including osmotic stress, oxidative stress, temperature changes, and acid stress. Moreover, both the DMSP synthesis and catabolism rates of LZB033 decreased after rpoN was knocked out. Our work provides essential insight into the regulatory function of RpoN, revealing that RpoN is a key determinant for LZB033 flagellar formation, motility, biofilm formation, and environmental fitness, as well as DMSP production and degradation.IMPORTANCE This study established an in-frame gene deletion method in the alphaproteobacterium Labrenzia aggregata LZB033 and generated an rpoN gene mutant. A comparison of the transcriptomes and phenotypic characteristics between the mutant and wild-type strains confirmed the role of RpoN in L. aggregata LZB033 flagellar formation, motility, biofilm formation, and carbon usage. Most importantly, RpoN is a key factor for survival under different environmental challenge conditions. Furthermore, the ability to synthesize and metabolize dimethylsulfoniopropionate (DMSP) was related to RpoN. These features revealed RpoN to be an important regulator of stress resistance and survival for L. aggregata LZB033 in marine environments.

Transcriptomic analysis of cultivated cotton Gossypium hirsutum provides insights into host responses upon whitefly-mediated transmission of cotton leaf curl disease

Cotton is a commercial and economically important crop that generates billions of dollars in annual revenue worldwide. However, cotton yield is affected by a sap-sucking insect Bemisia tabaci (whitefly), and whitefly-borne cotton leaf curl disease (CLCuD). The causative agent of devastating CLCuD is led by the viruses belonging to the genus Begomovirus (family Geminiviridae), collectively called cotton leaf curl viruses. Unfortunately, the extensively cultivated cotton (Gossypium hirsutum) species are highly susceptible and vulnerable to CLCuD. Yet, the concomitant influence of whitefly and CLCuD on the susceptible G. hirsutum transcriptome has not been interpreted. In the present study we have employed an RNA Sequencing (RNA-Seq) transcriptomics approach to explore the differential gene expression in susceptible G. hirsutum variety upon infection with viruliferous whiteflies. Comparative RNA-Seq of control and CLCuD infected plants was done using Illumina HiSeq 2500. This study yielded 468 differentially expressed genes (DEGs). Among them, we identified 220 up and 248 downregulated DEGs involved in disease responses and pathogen defense. We selected ten genes for downstream RT-qPCR analyses on two cultivars, Karishma and MNH 786 that are susceptible to CLCuD. We observed a similar expression pattern of these genes in both susceptible cultivars that was also consistent with our transcriptome data further implying a wider application of our global transcription study on host susceptibility to CLCuD. We next performed weighted gene co-expression network analysis that revealed six modules. This analysis also identified highly co-expressed genes as well as 55 hub genes that co-express with ≥ 50 genes. Intriguingly, most of these hub genes are shown to be downregulated and enriched in cellular processes. Under-expression of such highly co-expressed genes suggests their roles in favoring the virus and enhancing plant susceptibility to CLCuD. We also discuss the potential mechanisms governing the establishment of disease susceptibility. Overall, our study provides a comprehensive differential gene expression analysis of G. hirsutum under whitefly-mediated CLCuD infection. This vital study will advance the understanding of simultaneous effect of whitefly and virus on their host and aid in identifying important G. hirsutum genes which intricate in its susceptibility to CLCuD.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

Roles of TGF β and FGF Signals in the Lens: Tropomyosin Regulation for Posterior Capsule Opacity

Transforming growth factor (TGF) β and fibroblast growth factor (FGF) 2 are related to the development of posterior capsule opacification (PCO) after lens extraction surgery and other processes of epithelial–mesenchymal transition (EMT). Oxidative stress seems to activate TGF β1 largely through reactive oxygen species (ROS) production, which in turn alters the transcription of several survival genes, including lens epithelium-cell derived growth factor (LEDGF). Higher ROS levels attenuate LEDGF function, leading to down-regulation of peroxiredoxin 6 (Prdx6). TGF β is regulated by ROS in Prdx6 knock-out lens epithelial cells (LECs) and induces the up-regulation of tropomyosins (Tpms) 1/2, and EMT of LECs. Mouse and rat PCO are accompanied by elevated expression of Tpm2. Further, the expression of Tpm1/2 is induced by TGF β2 in LECs. Importantly, we previously showed that TGF β2 and FGF2 play regulatory roles in LECs in a contrasting manner. An injury-induced EMT of a mouse lens as a PCO model was attenuated in the absence of Tpm2. In this review, we present findings regarding the roles of TGF β and FGF2 in the differential regulation of EMT in the lens. Tpms may be associated with TGF β2- and FGF2-related EMT and PCO development.

RNA-seq-based monitoring of gene expression changes of viable but non-culturable state of Vibrio cholerae induced by cold seawater

Vibrio cholerae O1 is a natural inhabitant of aquatic environments and causes the acute diarrheal disease cholera. Entry into a viable but non-culturable (VBNC) state is a survival strategy by which V. cholerae withstands natural stresses and is important for the transition between the aquatic and host environments during the V. cholerae life cycle. In this study, the formation of VBNC V. cholerae induced by cold seawater exposure was investigated using RNA sequencing (RNA-seq). The analysis revealed that the expression of 1420 genes was changed on VBNC state formation. In the VBNC cells, genes related to biofilm formation, chitin utilization and stress responses were upregulated, whereas those related to cell division, morphology and ribosomal activity were mainly downregulated. The concurrent acquisition of a carbon source and the arrest of cell division in cells with low metabolic activity help bacteria increase their resistance to unfavourable environments. Moreover, two transcriptional regulators, SlmA and MetJ, were found to play roles in both VBNC formation and intestinal colonization, suggesting that some genes may function in both processes. This acquired knowledge will improve our understanding of the molecular mechanisms of stress tolerance and may help control future cholera infections and outbreaks.

Genome-wide identification of differentially expressed miRNAs induced by ethephon treatment in abscission layer cells of cotton (Gossypium hirsutum)

Cotton is one of the most important economic crops and its production is influenced by various adverse factors. Boll shedding at an inappropriate period causes the severe loss of cotton yield. Ethephon can promote the formation of abscission layer cells, resulting in boll shedding. To genome-widely investigate the differentially expressed miRNAs involved in the formation of abscission layer cells, two sRNA libraries were constructed using abscission tissues of cotton pedicels treated with ethephon or water. A total of 55,264,755 and 59,069,866 clean reads were generated in two libraries. 460 mature miRNAs were identified, including 301 known miRNAs and 159 novel miRNAs. Among them, one known gra-MIR530b and seven novel miRNAs (gar-novel-miR88-3p/gra-novel-miR8-3p, gar-novel-miR16-5p, gar-novel-miR100-3p, gar-novel-miR75-3p/gra-novel-miR92-3p,gar-novel-miR32-5p/gra-novel-miR30-5p, gar-novel-miR9-5p and gra-novel-miR61-3p) were expressed differentially in abscission tissues. Seven-two genes were predicted as potential targets of the eight miRNAs. Gene ontology (GO) analyses revealed that most of these target genes are auxin response factors (ARFs), GTPase activators, and the regulators of oxidative phosphorylation and photosynthesis. Hence, we speculate that these miRNAs might regulate cell separation and aging to affect the formation of bscission layer cells. Our findings will provide a new insight into the regulatory mechanism of miRNAs involved in boll shedding.

Genome-wide comparative and evolutionary analysis of Calmodulin-binding Transcription Activator (CAMTA) family in Gossypium species

The CAMTA gene family is crucial in managing both biotic and abiotic stresses in plants. Our comprehensive analysis of this gene family in cotton resulted in the identification of 6, 7 and 9 CAMTAs in three sequenced cotton species, i.e., Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum, respectively. All cotton CAMTAs were localized in the nucleus and possessed calmodulin-binding domain (CaMBD) as identified computationally. Phylogenetically four significant groups of cotton CAMTAs were identified out of which, Group II CAMTAs experienced higher evolutionary pressure, leading to a faster evolution in diploid cotton. The expansion of cotton CAMTAs in the genome was mainly due to segmental duplication. Purifying selection played a significant role in the evolution of cotton CAMTAs. Expression profiles of GhCAMTAs revealed that GhCAMTA2A.2 and GhCAMTA7A express profoundly in different stages of cotton fiber development. Positive correlation between expression of these two CAMTAs and fiber strength confirmed their functional relevance in fiber development. The promoter region of co-expressing genes network of GhCAMTA2A.2 and GhCAMTA7A showed a higher frequency of occurrence of CAMTA binding motifs. Our present study thus contributes to broad probing into the structure and probable function of CAMTA genes in Gossypium species.

De novo transcriptome assembly and analysis of differential gene expression following lipopolysaccharide challenge in Pelteobagrus fulvidraco

The yellow catfish, Pelteobagrus fulvidraco, has been recognized as an important freshwater aquaculture species in Eastern and Southeast Asia. To gain a better understanding of the immune response in P. fulvidraco, we analyzed its transcriptome following stimulation with lipopolysaccharide (LPS). Phosphate buffer saline (PBS) was used as control. Following assembly and annotation, 72,152 unigenes with an average length of 1090 bp were identified. A total of 370 differentially expressed genes (DEGs) in the P. fulvidraco were observed at 12 h post LPS treatment, including 197 up-regulated genes and 173 down-regulated genes. Clusters of Orthologous Groups of proteins (KOG/COG) annotation demonstrated that a total of 18,819 unigenes classified into 26 categories. Gene ontology (GO) analysis revealed 20 biological process subcategories, 7 cellular component subcategories and 20 molecular function subcategories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified immune responses pathways. Quantitative reverse transcription polymerase chain reaction measured the expression of 18 genes involved in the immune response. CXCL2-like chemokine (CXCL2), goose-type lysozyme (LYZ G), and cathepsin K (CTSK) were significantly up-regulated. This study enriches the P. fulvidraco transcriptome database and provides insight into the immune response of P. fulvidraco against infection.

Transcriptomic profiling of developing fiber in levant cotton (Gossypium herbaceum L.)

Cotton (Gossypium spp.) is an imperative economic crop of the globe due to its natural textile fiber. Molecular mechanisms of fiber development have been greatly revealed in allotetraploid cotton but remained unexplored in Gossypium herbaceum. G. herbaceum can withstand the rigors of nature like drought and pests but produce coarse lint. This undesirable characteristic strongly needs the knowledge of fiber development at molecular basis. The present study reported the transcriptome sequence of the developing fiber of G. herbaceum on pyrosequencing and its analysis. About 1.38 million raw and 1.12 million quality trimmed reads were obtained followed by de novo assembly-generated 20,125 unigenes containing 14,882 coding sequences (CDs). BLASTx-based test of homology indicated that A1-derived transcripts shared a high similarity with Gossypium arboreum (A2). Functional annotation of the CDs using the UniProt categorized them into biological processes, cellular components, and molecular function, COG classification showed that a large number of CDs have significant homology in COG database (6215 CDs), and mapping of CDs with Kyoto Encyclopedia of Genes and Genomes (KEGG) database generated 200 pathways ultimately showing predominant engagement in the fiber development process. Transcription factors were predicted by comparison with Plant Transcription Factor Database, and their differential expression between stages exposed their important regulatory role in fiber development. Differential expression analysis based on reads per kilobase of transcript per million mapped reads (RPKM) value revealed activities of specific gene related to carbohydrate and lipid synthesis, carbon metabolism, energy metabolism, signal transduction, etc., at four stages of fiber development, and was validated by qPCR. Overall, this study will help as a valuable foundation for diploid cotton fiber improvement.

Plant Systems Biology at the Single-Cell Level

Our understanding of plant biology is increasingly being built upon studies using 'omics and system biology approaches performed at the level of the entire plant, organ, or tissue. Although these approaches open new avenues to better understand plant biology, they suffer from the cellular complexity of the analyzed sample. Recent methodological advances now allow plant scientists to overcome this limitation and enable biological analyses of single-cells or single-cell-types. Coupled with the development of bioinformatics and functional genomics resources, these studies provide opportunities for high-resolution systems analyses of plant phenomena. In this review, we describe the recent advances, current challenges, and future directions in exploring the biology of single-cells and single-cell-types to enhance our understanding of plant biology as a system.

Identification of candidate genes involved in the sugar metabolism and accumulation during pear fruit post-harvest ripening of 'Red Clapp's Favorite' (Pyrus communis L.) by transcriptome analysis

Background

Pear (Pyrus spp.) is a popular fruit that is commercially cultivated in most temperate regions. In fruits, sugar metabolism and accumulation are important factors for fruit organoleptic quality. Post-harvest ripening is a special feature of ‘Red Clapp’s Favorite’.

Results

In this study, transcriptome sequencing based on the Illumina platform generated 23.8 - 35.8 million unigenes of nine cDNA libraries constructed using RNAs from the ‘Red Clapp’s Favorite’ pear variety with different treatments, in which 2629 new genes were discovered, and 2121 of them were annotated. A total of 2146 DEGs, 3650 DEGs, 1830 DEGs from each comparison were assembled. Moreover, the gene expression patterns of 8 unigenes related to sugar metabolism revealed by qPCR. The main constituents of soluble sugars were fructose and glucose after pear fruit post-harvest ripening, and five unigenes involved in sugar metabolism were discovered.

Conclusions

Our study not only provides a large-scale assessment of transcriptome resources of ‘Red Clapp’s Favorite’ but also lays the foundation for further research into genes correlated with sugar metabolism.

Electronic supplementary material

The online version of this article (10.1186/s41065-017-0046-0) contains supplementary material, which is available to authorized users.

Researches on Transcriptome Sequencing in the Study of Traditional Chinese Medicine

Due to its incomparable advantages, the application of transcriptome sequencing in the study of traditional Chinese medicine attracts more and more attention of researchers, which greatly promote the development of traditional Chinese medicine. In this paper, the applications of transcriptome sequencing in traditional Chinese medicine were summarized by reviewing recent related papers.

A comparative genomic and transcriptomic analysis at the level of isolated root hair cells reveals new conserved root hair regulatory elements

KEY MESSAGE

A comparative transcriptomic and genomic analysis between Arabidopsis thaliana and Glycine max root hair genes reveals the evolution of the expression of plant genes after speciation and whole genome duplication. Our understanding of the conservation and divergence of the expression patterns of genes between plant species is limited by the quality of the genomic and transcriptomic resources available. Specifically, the transcriptomes generated from plant organs are the reflection of the contribution of the different cell types composing the samples weighted by their relative abundances in the sample. These contributions can vary between plant species leading to the generation of datasets which are difficult to compare. To gain a deeper understanding of the evolution of gene transcription in and between plant species, we performed a comparative transcriptomic and genomic analysis at the level of one single plant cell type, the root hair cell, and between two model plants: Arabidopsis (Arabidopsis thaliana) and soybean (Glycine max). These two species, which diverged 90 million years ago, were selected as models based on the large amount of genomic and root hair transcriptomic information currently available. Our analysis revealed in detail the transcriptional divergence and conservation between soybean paralogs (i.e., the soybean genome is the product of two successive whole genome duplications) and between Arabidopsis and soybean orthologs in this single plant cell type. Taking advantage of this evolutionary study, we combined bioinformatics, molecular, cellular and microscopic tools to characterize plant promoter sequences and the discovery of two root hair regulatory elements (RHE1 and RHE2) consistently and specifically active in plant root hair cells.

Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton

BACKGROUND

Cotton (Gossypium spp.) fibers are single-celled elongated trichomes, the molecular aspects of genetic variation in fiber length (FL) among genotypes are currently unknown. In this study, two backcross inbred lines (BILs), i.e., NMGA-062 ("Long") and NMGA-105 ("Short") with 32.1 vs. 27.2 mm in FL, respectively, were chosen to perform RNA-Seq on developing fibers at 10 days post anthesis (DPA). The two BILs differed in 4 quantitative trait loci (QTL) for FL and were developed from backcrosses between G. hirsutum as the recurrent parent and G. barbadense.

RESULTS

In total, 51.7 and 54.3 million reads were obtained and assembled to 49,508 and 49,448 transcripts in the two genotypes, respectively. Of 1551 differentially expressed genes (DEGs) between the two BILs, 678 were up-regulated and 873 down-regulated in "Long"; and 703 SNPs were identified in 339 DEGs. Further physical mapping showed that 8 DEGs were co-localized with the 4 FL QTL identified in the BIL population containing the two BILs. Four SNP markers in 3 DEGs that showed significant correlations with FL were developed. Among the three candidate genes encoding for proline-rich protein, D-cysteine desulfhydrase, and thaumatin-like protein, a SNP of thaumatin-like protein gene showed consistent correlations with FL across all testing environments.

CONCLUSIONS

This study represents one of the first investigations of positional candidate gene approach of QTL in cotton in integrating transcriptome and SNP identification based on RNA-Seq with linkage and physical mapping of QTL and genes, which will facilitate eventual cloning and identification of genes responsible for FL QTL. The candidate genes may serve as the foundation for further in-depth studies of the molecular mechanism of natural variation in fiber elongation.

MicroRNA expression profiles during cotton (Gossypium hirsutum L) fiber early development

The role of microRNAs (miRNAs) during cotton fiber development remains unclear. Here, a total of 54 miRNAs belonging to 39 families were selected to characterize miRNA regulatory mechanism in eight different fiber development stages in upland cotton cv BM-1. Among 54 miRNAs, 18 miRNAs were involved in cotton fiber initiation and eight miRNAs were related to fiber elongation and secondary wall biosynthesis. Additionally, 3,576 protein-coding genes were candidate target genes of these miRNAs, which are potentially involved in cotton fiber development. We also investigated the regulatory network of miRNAs and corresponding targets in fiber initiation and elongation, and secondary wall formation. Our Gene Ontology-based term classification and KEGG-based pathway enrichment analyses showed that the miRNA targets covered 220 biological processes, 67 molecular functions, 45 cellular components, and 10 KEGG pathways. Three of ten KEGG pathways were involved in lignan synthesis, cell elongation, and fatty acid biosynthesis, all of which have important roles in fiber development. Overall, our study shows the potential regulatory roles of miRNAs in cotton fiber development and the importance of miRNAs in regulating different cell types. This is helpful to design miRNA-based biotechnology for improving fiber quality and yield.

Genome Editing in Cotton with the CRISPR/Cas9 System

Genome editing is an important tool for gene functional studies as well as crop improvement. The recent development of the CRISPR/Cas9 system using single guide RNA molecules (sgRNAs) to direct precise double strand breaks in the genome has the potential to revolutionize agriculture. Unfortunately, not all sgRNAs are equally efficient and it is difficult to predict their efficiency by bioinformatics. In crops such as cotton (Gossypium hirsutum L.), with labor-intensive and lengthy transformation procedures, it is essential to minimize the risk of using an ineffective sgRNA that could result in the production of transgenic plants without the desired CRISPR-induced mutations. In this study, we have developed a fast and efficient method to validate the functionality of sgRNAs in cotton using a transient expression system. We have used this method to validate target sites for three different genes GhPDS, GhCLA1, and GhEF1 and analyzed the nature of the CRISPR/Cas9-induced mutations. In our experiments, the most frequent type of mutations observed in cotton cotyledons were deletions (∼64%). We prove that the CRISPR/Cas9 system can effectively produce mutations in homeologous cotton genes, an important requisite in this allotetraploid crop. We also show that multiple gene targeting can be achieved in cotton with the simultaneous expression of several sgRNAs and have generated mutations in GhPDS and GhEF1 at two target sites. Additionally, we have used the CRISPR/Cas9 system to produce targeted gene fragment deletions in the GhPDS locus. Finally, we obtained transgenic cotton plants containing CRISPR/Cas9-induced gene editing mutations in the GhCLA1 gene. The mutation efficiency was very high, with 80.6% of the transgenic lines containing mutations in the GhCLA1 target site resulting in an intense albino phenotype due to interference with chloroplast biogenesis.

Genome-wide identification and characterization of microRNAs differentially expressed in fibers in a cotton phytochrome A1 RNAi line

Cotton fiber is an important commodity throughout the world. Fiber property determines fiber quality and commercial values. Previous studies showed that silencing phytochrome A1 gene (PHYA1) by RNA interference in Upland cotton (Gossypium hirsutum L. cv. Coker 312) had generated PHYA1 RNAi lines with simultaneous improvements in fiber quality (longer, stronger and finer fiber) and other key agronomic traits. Characterization of the altered molecular processes in these RNAi genotypes and its wild-type controls is a great interest to better understand the PHYA1 RNAi phenotypes. In this study, a total of 77 conserved miRNAs belonging to 61 families were examined in a PHYA1 RNAi line and its parental Coker 312 genotype by using multiplex sequencing. Of these miRNAs, seven (miR7503, miR7514, miR399c, miR399d, miR160, miR169b, and miR2950) were found to be differentially expressed in PHYA1 RNAi cotton. The target genes of these differentially expressed miRNAs were involved in the metabolism and signaling pathways of phytohormones, which included Gibberellin, Auxin and Abscisic Acid. The expression of several MYB transcription factors was also affected by miRNAs in RNAi cotton. In addition, 35 novel miRNAs (novel miR1-novel miR35) were identified in fibers for the first time in this study. Target genes of vast majority of these novel miRNAs were also predicted. Of these, nine novel miRNAs (novel-miR1, 2, 16, 19, 26, 27, 28, 31 and 32) were targeted to cytochrome P450-like TATA box binding protein (TBP). The qRT-PCR confirmed expression levels of several differentially regulated miRNAs. Expression patterns of four miRNAs-targets pairs were also examined via RNA deep sequencing. Together, the results imply that the regulation of miRNA expression might confer to the phenotype of the PHYA1 RNAi line(s) with improved fiber quality.

Genome Editing in Cotton with the CRISPR/Cas9 System

Genome editing is an important tool for gene functional studies as well as crop improvement. The recent development of the CRISPR/Cas9 system using single guide RNA molecules (sgRNAs) to direct precise double strand breaks in the genome has the potential to revolutionize agriculture. Unfortunately, not all sgRNAs are equally efficient and it is difficult to predict their efficiency by bioinformatics. In crops such as cotton (Gossypium hirsutum L.), with labor-intensive and lengthy transformation procedures, it is essential to minimize the risk of using an ineffective sgRNA that could result in the production of transgenic plants without the desired CRISPR-induced mutations. In this study, we have developed a fast and efficient method to validate the functionality of sgRNAs in cotton using a transient expression system. We have used this method to validate target sites for three different genes GhPDS, GhCLA1, and GhEF1 and analyzed the nature of the CRISPR/Cas9-induced mutations. In our experiments, the most frequent type of mutations observed in cotton cotyledons were deletions (∼64%). We prove that the CRISPR/Cas9 system can effectively produce mutations in homeologous cotton genes, an important requisite in this allotetraploid crop. We also show that multiple gene targeting can be achieved in cotton with the simultaneous expression of several sgRNAs and have generated mutations in GhPDS and GhEF1 at two target sites. Additionally, we have used the CRISPR/Cas9 system to produce targeted gene fragment deletions in the GhPDS locus. Finally, we obtained transgenic cotton plants containing CRISPR/Cas9-induced gene editing mutations in the GhCLA1 gene. The mutation efficiency was very high, with 80.6% of the transgenic lines containing mutations in the GhCLA1 target site resulting in an intense albino phenotype due to interference with chloroplast biogenesis.

Effect of oxidative phosphorylation signaling pathway on silkworm midgut following exposure to phoxim

Organophosphate pesticides are applied widely in the world for agricultural purposes, and their exposures often resulted in non-cocooning of Bombyx mori in China. Silkworm midgut is the major organ for digestion and nutrient absorption, importantly it is also a barrier against foreign substances and chemical pesticides. The purpose of this study was to determine the mechanism of oxidative injury in silkworm midgut with phoxim induction. The results showed that the transcription level of oxidative phosphorylation signaling pathway genes of midgut under phoxim stress. Digital gene expression (DGE) analysis revealed that 24 electron transport chain (ETC)-related genes were upregulated. Quantitative real time polymerase chain reaction results indicated that the ETC the genes encoding NADH-CoQ1, Succinic-Q, cyt c reductase-S, cyt c oxidase-S, cytochrome c oxidase polypeptide IV, ATP synthase, and vacuolar H+ ATP synthase were all significantly up-regulated by 1.50-, 1.31-, 1.42-, 1.44-, 1.70-, 2.03- and 1.43-fold, respectively. Phoxim induction enhanced the activity of ETC complex in mitochondria, and induced the accumulation of ROS in midgut. These results indicated that trace phoxim enhanced respiration in midgut, and the imbalance between the activity changes of ETC may led to reactive oxygen species accumulation. The ETC of mitochondria may be potential biomarkers of midgut toxicity in B. mori caused by phoxim exposure. © 2015 Wiley Periodicals, Inc. Environ Toxicol 32: 167-175, 2017.

Transcriptome Analysis of Gene Expression during Chinese Water Chestnut Storage Organ Formation

The product organ (storage organ; corm) of the Chinese water chestnut has become a very popular food in Asian countries because of its unique nutritional value. Corm formation is a complex biological process, and extensive whole genome analysis of transcripts during corm development has not been carried out. In this study, four corm libraries at different developmental stages were constructed, and gene expression was identified using a high-throughput tag sequencing technique. Approximately 4.9 million tags were sequenced, and 4,371,386, 4,372,602, 4,782,494, and 5,276,540 clean tags, including 119,676, 110,701, 100,089, and 101,239 distinct tags, respectively, were obtained after removal of low-quality tags from each library. More than 39% of the distinct tags were unambiguous and could be mapped to reference genes, while 40% were unambiguous tag-mapped genes. After mapping their functions in existing databases, a total of 11,592, 10,949, 10,585, and 7,111 genes were annotated from the B1, B2, B3, and B4 libraries, respectively. Analysis of the differentially expressed genes (DEGs) in B1/B2, B2/B3, and B3/B4 libraries showed that most of the DEGs at the B1/B2 stages were involved in carbohydrate and hormone metabolism, while the majority of DEGs were involved in energy metabolism and carbohydrate metabolism at the B2/B3 and B3/B4 stages. All of the upregulated transcription factors and 9 important genes related to product organ formation in the above four stages were also identified. The expression changes of nine of the identified DEGs were validated using a quantitative PCR approach. This study provides a comprehensive understanding of gene expression during corm formation in the Chinese water chestnut.

Global analysis of the developmental dynamics of Gossypium hirsutum based on strand-specific transcriptome

Cotton is an economically important crop that provides both natural fiber and by-products such as oil and protein. Its global gene expression could provide insight into the biological processes underlying growth and development, which involve suites of genes expressed with temporal and spatial control by regulatory networks. Generally, the goal for cotton breeding is improvement of the fiber; thus, most previous research has focused on identifying genes specific to the fiber. However, seeds may also play an important role in fiber development. In this study, we constructed and systematically analyzed 21 strand-specific RNA-Seq libraries for Gossypium hirsutum, covering different tissues, organs and development stages, from which approximately 970 million reads were generated to provide a global view of gene expression during cotton development. The organ (tissue)-specific gene expression patterns were investigated, providing further insight into the dynamic programming associated with developmental processes and a way to study the coordination of development between fiber cells and ovules. Series of transcription factors and seed-specific genes have been identified as candidate genes that could elucidate key mechanisms and regulatory networks in nutrient accumulation during ovule development and in fiber development. This study reports comprehensive transcriptome dynamics at various stages of cotton development and will serve as a valuable genome-wide transcriptome resource for initial gene discovery and functional characterization of genes in cotton.

New Insights into the Formation of Viable but Nonculturable Escherichia coli O157:H7 Induced by High-Pressure CO2

The formation of viable but nonculturable (VBNC) Escherichia coli O157:H7 induced by high-pressure CO₂ (HPCD) was investigated using RNA sequencing (RNA-Seq) transcriptomics and isobaric tag for relative and absolute quantitation (iTRAQ) proteomic methods. The analyses revealed that 97 genes and 56 proteins were significantly changed upon VBNC state entry. Genes and proteins related to membrane transport, central metabolisms, DNA replication, and cell division were mainly downregulated in the VBNC cells. This caused low metabolic activity concurrently with a division arrest in cells, which may be related to VBNC state formation. Cell division repression and outer membrane overexpression were confirmed to be involved in VBNC state formation by homologous expression of z2046 coding for transcriptional repressor and ompF encoding outer membrane protein F. Upon VBNC state entry, pyruvate catabolism in the cells shifted from the tricarboxylic acid (TCA) cycle toward the fermentative route; this led to a low level of ATP. Combating the low energy supply, ATP production in the VBNC cells was compensated by the degradation of l-serine and l-threonine, the increased AMP generation, and the enhanced electron transfer. Furthermore, tolerance of the cells with respect to HPCD-induced acid, oxidation, and high CO₂ stresses was enhanced by promoting the production of ammonia and NADPH and by reducing CO₂ production during VBNC state formation. Most genes and proteins related to pathogenicity were downregulated in the VBNC cells. This would decrease the cell pathogenicity, which was confirmed by adhesion assays. In conclusion, the decreased metabolic activity, repressed cell division, and enhanced survival ability in E. coli O157:H7 might cause HPCD-induced VBNC state formation.

Escherichia coli O157:H7 has been implicated in large foodborne outbreaks worldwide. It has been reported that the presence of as few as 10 cells in food could cause illness. However, the presence of only 0.73 to 1.5 culturable E. coli O157:H7 cells in salted salmon roe caused infection in Japan. Investigators found that E. coli O157:H7 in the viable but nonculturable (VBNC) state was the source of the outbreak. So far, formation mechanisms of VBNC state are not well known. In a previous study, we demonstrated that high-pressure CO₂ (HPCD) could induce the transition of E. coli O157:H7 into the VBNC state. In this study, we used RNA-Seq transcriptomic analysis combined with the iTRAQ proteomic method to investigate the formation of VBNC E. coli O157:H7 induced by HPCD treatment. Finally, we proposed a putative formation mechanism of the VBNC cells induced by HPCD, which may provide a theoretical foundation for controlling the VBNC state entry induced by HPCD treatment.

Rice (Oryza sativa) Laccases Involved in Modification and Detoxification of Herbicides Atrazine and Isoproturon Residues in Plants

Atrazine (ATR) and isoproturon (IPU) as herbicides have become serious environmental contaminants due to their overuse in crop production. Although ATR and IPU in soils are easily absorbed by many crops, the mechanisms for their degradation or detoxification in plants are poorly understood. This study identified a group of novel genes encoding laccases (EC 1.10.3.2) that are possibly involved in catabolism or detoxification of ATR and IPU residues in rice. Transcriptome profiling shows at least 22 differentially expressed laccase genes in ATR/IPU-exposed rice. Some of the laccase genes were validated by RT-PCR analysis. The biochemical properties of the laccases were analyzed, and their activities in rice were induced under ATR/IPU exposure. To investigate the roles of laccases in degrading or detoxifying ATR/IPU in rice, transgenic yeast cells (Pichia pastoris X-33) expressing two rice laccase genes (LOC_Os01g63180 and LOC_Os12g15680) were generated. Both transformants were found to accumulate less ATR/IPU compared to the control. The ATR/IPU-degraded products in the transformed yeast cells using UPLC-TOF-MS/MS were further characterized. Two metabolites, hydroxy-dehydrogenated atrazine (HDHA) and 2-OH-isopropyl-IPU, catalyzed by laccases were detected in the eukaryotic cells. These results indicate that the laccase-coding genes identified here could confer degradation or detoxification of the herbicides and suggest that the laccases could be one of the important enzymatic pathways responsible for ATR/IPU degradation/detoxification in rice.

Integrated Analysis of the Roles of Long Noncoding RNA and Coding RNA Expression in Sheep (Ovis aries) Skin during Initiation of Secondary Hair Follicle

Initiation of hair follicle (HF) is the first and most important stage of HF morphogenesis. However the precise molecular mechanism of initiation of hair follicle remains elusive. Meanwhile, in previous study, the more attentions had been paid to the function of genes, while the roles of non-coding RNAs (such as long noncoding RNA and microRNA) had not been described. Therefore, the roles of long noncoding RNA(LncRNA) and coding RNA in sheep skin during the initiation of sheep secondary HF were integrated and analyzed, by using strand-specific RNA sequencing (ssRNA-seq).A total of 192 significant differentially expressed genes were detected, including 67 up-regulated genes and 125 down-regulated genes between stage 0 and stage 1 of HF morphogenesis during HF initiation. Only Wnt2, FGF20 were just significant differentially expressed among Wnt, Shh, Notch and BMP signaling pathways. Further expression profile analysis of lncRNAs showed that 884 novel lncRNAs were discovered in sheep skin expression profiles. A total of 15 lncRNAs with significant differential expression were detected, 6 up-regulated and 9 down-regulated. Among of differentially expressed genes and LncRNA, XLOC002437 lncRNA and potential target gene COL6A6 were all significantly down-regulated in stage 1. Furthermore, by using RNAhybrid, XLOC005698 may be as a competing endogenous RNA ''sponges" oar-miR-3955-5p activity. Gene Ontology and KEGG pathway analyses indicated that the significantly enriched pathway was peroxisome proliferator-activated receptors (PPARs) pathway (corrected P-value < 0.05), indicating that PPAR pathway is likely to play significant roles during the initiation of secondary HF.Results suggest that the key differentially expressed genes and LncRNAs may be considered as potential candidate genes for further study on the molecular mechanisms of HF initiation, as well as supplying some potential values for understanding human hair disorders.

Hypobaric Treatment Effects on Chilling Injury, Mitochondrial Dysfunction, and the Ascorbate-Glutathione (AsA-GSH) Cycle in Postharvest Peach Fruit

In this study, hypobaric treatment effects were investigated on chilling injury, mitochondrial dysfunction, and the ascorbate-glutathione (AsA-GSH) cycle in peach fruit stored at 0 °C. Internal browning of peaches was dramatically reduced by applying 10-20 kPa pressure. Hypobaric treatment markedly inhibited membrane fluidity increase, whereas it kept mitochondrial permeability transition pore (MPTP) concentration and cytochrome C oxidase (CCO) and succinic dehydrogenase (SDH) activity relatively high in mitochondria. Similarly, 10-20 kPa pressure treatment reduced the level of decrease observed in AsA and GSH concentrations, while it enhanced ascorbate peroxidase (APX), glutathione reductase (GR), and monodehydroascorbate reductase (MDHAR) activities related to the AsA-GSH cycle. Furthermore, comparative transcriptomic analysis showed that differentially expressed genes (DEGs) associated with the metabolism of glutathione, ascorbate, and aldarate were up-regulated in peaches treated with 10-20 kPa for 30 days at 0 °C. Genes encoding GR, MDHAR, and APX were identified and exhibited higher expression in fruits treated with low pressure than in fruits treated with normal atmospheric pressure. Our findings indicate that the alleviation of chilling injury by hypobaric treatment was associated with preventing mitochondrial dysfunction and triggering the AsA-GSH cycle by the transcriptional up-regulation of related enzymes.

Decipher the Molecular Response of Plant Single Cell Types to Environmental Stresses

The analysis of the molecular response of entire plants or organs to environmental stresses suffers from the cellular complexity of the samples used. Specifically, this cellular complexity masks cell-specific responses to environmental stresses and logically leads to the dilution of the molecular changes occurring in each cell type composing the tissue/organ/plant in response to the stress. Therefore, to generate a more accurate picture of these responses, scientists are focusing on plant single cell type approaches. Several cell types are now considered as models such as the pollen, the trichomes, the cotton fiber, various root cell types including the root hair cell, and the guard cell of stomata. Among them, several have been used to characterize plant response to abiotic and biotic stresses. In this review, we are describing the various -omic studies performed on these different plant single cell type models to better understand plant cell response to biotic and abiotic stresses.

RNA-Seq-Mediated Transcriptome Analysis of a Fiberless Mutant Cotton and Its Possible Origin Based on SNP Markers

As the longest known single-celled trichomes, cotton (Gossypium L.) fibers constitute a classic model system to investigate cell initiation and elongation. In this study, we used a high-throughput transcriptome sequencing technology to identify fiber-initiation-related single nucleotide polymorphism (SNP) markers and differentially expressed genes (DEGs) between the wild-type (WT) Upland cotton (G. hirsutum) Xuzhou 142 and its natural fuzzless-lintless mutant Xuzhou 142 fl. Approximately 700 million high-quality cDNA reads representing over 58 Gb of sequences were obtained, resulting in the identification of 28,610 SNPs--of which 17,479 were novel--from 13,960 expressed genes. Of these SNPs, 50% of SNPs in fl were identical to those of G. barbadense, which suggests the likely origin of the fl mutant from an interspecific hybridization between Xuzhou 142 and an unknown G. barbadense genotype. Of all detected SNPs, 15,555, 12,750, and 305 were classified as non-synonymous, synonymous, and pre-terminated ones, respectively. Moreover, 1,352 insertion/deletion polymorphisms (InDels) were also detected. A total of 865 DEGs were identified between the WT and fl in ovules at -3 and 0 days post-anthesis, with 302 candidate SNPs selected from these DEGs for validation by a high-resolution melting analysis and Sanger sequencing in seven cotton genotypes. The number of genotypic pairwise polymorphisms varied from 43 to 302, indicating that the identified SNPs are reliable. These SNPs should serve as good resources for breeding and genetic studies in cotton.

Comparative Proteomic Analysis of Cotton Fiber Development and Protein Extraction Method Comparison in Late Stage Fibers

The distinct stages of cotton fiber development and maturation serve as a single-celled model for studying the molecular mechanisms of plant cell elongation, cell wall development and cellulose biosynthesis. However, this model system of plant cell development is compromised for proteomic studies due to a lack of an efficient protein extraction method during the later stages of fiber development, because of a recalcitrant cell wall and the presence of abundant phenolic compounds. Here, we compared the quality and quantities of proteins extracted from 25 dpa (days post anthesis) fiber with multiple protein extraction methods and present a comprehensive quantitative proteomic study of fiber development from 10 dpa to 25 dpa. Comparative analysis using a label-free quantification method revealed 287 differentially-expressed proteins in the 10 dpa to 25 dpa fiber developmental period. Proteins involved in cell wall metabolism and regulation, cytoskeleton development and carbohydrate metabolism among other functional categories in four fiber developmental stages were identified. Our studies provide protocols for protein extraction from maturing fiber tissues for mass spectrometry analysis and expand knowledge of the proteomic profile of cotton fiber development.

Identification of low Ca(2+) stress-induced embryo apoptosis response genes in Arachis hypogaea by SSH-associated library lift (SSHaLL)

Calcium is a universal signal in the regulation of wide aspects in biology, but few are known about the function of calcium in the control of early embryo development. Ca(2+) deficiency in soil induces early embryo abortion in peanut, producing empty pods, which is a general problem; however, the underlying mechanism remains unclear. In this study, embryo abortion was characterized to be caused by apoptosis marked with cell wall degradation. Using a method of SSH cDNA libraries associated with library lift (SSHaLL), 62 differentially expressed genes were isolated from young peanut embryos. These genes were classified to be stress responses, catabolic process, carbohydrate and lipid metabolism, embryo morphogenesis, regulation, etc. The cell retardation with cell wall degradation was caused by up-regulated cell wall hydrolases and down-regulated cellular synthases genes. HsfA4a, which was characterized to be important to embryo development, was significantly down-regulated under Ca(2+) -deficient conditions from 15 days after pegging (DAP) to 30 DAP. Two AhCYP707A4 genes, encoding abscisic acid (ABA) 8'-hydroxylases, key enzymes for ABA catabolism, were up-regulated by 21-fold under Ca(2+) -deficient conditions upstream of HsfA4a, reducing the ABA level in early embryos. Over-expression of AhCYP707A4 in Nicotiana benthamiana showed a phenotype of low ABA content with high numbers of aborted embryos, small pods and less seeds, which confirms that AhCYP707A4 is a key player in regulation of Ca(2+) deficiency-induced embryo abortion via ABA-mediated apoptosis. The results elucidated the mechanism of low Ca(2+) -induced embryo abortion and described the method for other fields of study.

A cotton miRNA is involved in regulation of plant response to salt stress

The present study functionally identified a new microRNA (microRNA ovual line 5, miRNVL5) with its target gene GhCHR from cotton (Gossypium hirsutum). The sequence of miRNVL5 precursor is 104 nt long, with a well developed secondary structure. GhCHR contains two DC1 and three PHD Cys/His-rich domains, suggesting that GhCHR encodes a zinc-finger domain-containing transcription factor. miRNVL5 and GhCHR express at various developmental stages of cotton. Under salt stress (50-400 mM NaCl), miRNVL5 expression was repressed, with concomitant high expression of GhCHR in cotton seedlings. Ectopic expression of GhCHR in Arabidopsis conferred salt stress tolerance by reducing Na(+) accumulation in plants and improving primary root growth and biomass. Interestingly, Arabidopsis constitutively expressing miRNVL5 showed hypersensitivity to salt stress. A GhCHR orthorlous gene At2g44380 from Arabidopsis that can be cleaved by miRNVL5 was identified by degradome sequencing, but no confidential miRNVL5 homologs in Arabidopsis have been identified. Microarray analysis of miRNVL5 transgenic Arabidopsis showed six downstream genes (CBF1, CBF2, CBF3, ERF4, AT3G22920, and AT3G49200), which were induced by salt stress in wild-type but repressed in miRNVL5-expressing Arabidopsis. These results indicate that miRNVL5 is involved in regulation of plant response to salt stress.

RNA-sequencing Oryza sativa transcriptome in response to herbicide isoprotruon and characterization of genes involved in IPU detoxification

The comprehensive analysis of transcriptome and UPLC-MS/MS in rice was performed to explore the regulatory mechanism of mRNA level and chemical metabolism in response to herbicide isoproturon.