CottonGen: a genomics, genetics and breeding database for cotton research

GhN/AINV13 positively regulates cotton stress tolerance by interacting with the 14-3-3 protein

Neutral/alkaline invertases (N/AINVs) are sucrose hydrolases with important roles in plants. In this study, 15, 15, 15, 29, and 30 N/AINVs were identified in the Gossypium species, G. raimondii, G. herbaceum, G. arboreum, G. hirsutum, and G. barbadense, respectively. Along with two previously discovered branches, α and β, a new clade γ was first discovered in our study. Investigation of gene collinearity showed that whole-genome duplication (WGD) and polyploidization were responsible for the expansion of the N/AINV gene family in allopolyploid Gossypium. Moreover, expression patterns revealed that GhN/AINV3/13/17/23/24/28 from the β clade is highly expressed during the period of fiber initiation. The invertase activity of GhN/AINV13 and GhN/AINV23 were confirmed by restoring defects of invertase-deficient yeast mutant SEY2102. Treatments of abiotic stress showed that most GhN/AINVs were induced in response to polyethylene glycol (PEG) or salt stress. A virus-induced gene-silencing (VIGS) experiment and yeast two-hybrid assay demonstrated that GhN/AINV13 may interact with their positive regulators Gh14-3-3 proteins and participate in the fiber initiation or stress tolerance of cotton. Our results provided fundamental information regarding N/AINVs and highlight their potential functions in cotton stress tolerance.

The NanDeSyn database for Nannochloropsis systems and synthetic biology

Nannochloropsis species, unicellular industrial oleaginous microalgae, are model organisms for microalgal systems and synthetic biology. To facilitate community-based annotation and mining of the rapidly accumulating functional genomics resources, we have initiated an international consortium and present a comprehensive multi-omics resource database named Nannochloropsis Design and Synthesis (NanDeSyn; http://nandesyn.single-cell.cn). Via the Tripal toolkit, it features user-friendly interfaces hosting genomic resources with gene annotations and transcriptomic and proteomic data for six Nannochloropsis species, including two updated genomes of Nannochloropsis oceanica IMET1 and Nannochloropsis salina CCMP1776. Toolboxes for search, Blast, synteny view, enrichment analysis, metabolic pathway analysis, a genome browser, etc. are also included. In addition, functional validation of genes is indicated based on phenotypes of mutants and relevant bibliography. Furthermore, epigenomic resources are also incorporated, especially for sequencing of small RNAs including microRNAs and circular RNAs. Such comprehensive and integrated landscapes of Nannochloropsis genomics and epigenomics will promote and accelerate community efforts in systems and synthetic biology of these industrially important microalgae.

Elucidation of sequence polymorphism in fuzzless-seed cotton lines

Most commercially produced cotton cultivars have two types of fibers on the seed coat, short fuzz and long lint. Lint fiber is used in the textile industry, while fuzz is considered an undesirable trait. Both types of fibers are believed to be controlled by the same regulators; however, their mechanisms of actions are still obscure. Cotton fiber mutants provide an excellent system to study the genes that regulate fiber development. Here we described four uncharacterized and three previously reported cotton mutants with fuzzless seed phenotypes. To evaluate whether or not the genes previously associated with fuzzless seed phenotypes have mutations we sequenced whole genomic DNA of seven mutants and wild type varieties. We identified multiple polymorphic changes among the tested genes. Non-synonymous SNPs in the coding region of the MML3-A gene was common in the six mutant lines tested in this study, showing both dominant and recessive fuzzless phenotypes. We have mapped the locus of the causative mutation for one of the uncharacterized fuzzless lines using an F₂ population that originated from a cross between the dominant fuzzless mutant and a wild type. Further, we have clarified the current knowledge about the causative n₂ mutations by analyzing the sequence data and previously reported mapping data. The key genes and possible mechanisms of fiber differentiation are discussed in this study.

Sequencing Multiple Cotton Genomes Reveals Complex Structures and Lays Foundation for Breeding

Cotton is a major fiber plant, which provides raw materials for clothing, protecting humans from the harsh environment of cold or hot weathers, enriching the culture and custom of human societies. Due to its importance, the diploid and tetraploid genomes of different cotton plants have been repeatedly sequenced to obtain their complete and fine genome sequences. These valuable genome data sets revealed the evolutionary past of the cotton plants, which were recursively affected by polyploidization, with a decaploidization contributing to the formation of the genus Gossypium, and a neo-tetraploidization contributing to the formation of nowadays widely cultivated cotton plants. Post-polyploidization genome instability resulted in numerous structural changes of the genomes, such as gene loss, DNA inversion and translocation, illegitimate recombination, and accumulation of repetitive sequences, and functional innovation accompanied by elevated evolutionary rates of genes. Many these changes have been asymmetric between subgnomes of the tetraploid cottons, rendering their divergent profiles of biological regulation and function. The availability of whole-genome sequences has now paved the way to identify and clone functional genes, e.g., those relating to fiber development, and to enhance breeding efforts to cultivate cottons to produce high-yield and high-quality fibers, and to resist environmental and biological stress.

Application of Pedimap: a pedigree visualization tool to facilitate the decisioning of rice breeding in Sri Lanka

The development of rice cultivars with desirable traits is essential. The decision-making is a crucial step in rice breeding programs. Breeders can make efficient and pragmatic decisions if an organized pedigree visualization platform is available for the accessions and cultivars in rice breeding germplasm. In the present study, the available data of all the rice varieties released by Rice Research and Development Institute, Sri Lanka, and the related landraces and genotypes were arranged in Pedimap, a pedigree visualization tool. Pedimap can showcase pedigree relationships, phenotypic, and molecular data. The identity by descent probabilities were calculated using FlexQTL software and included in the Pedimap database. The parentage selection based on the variations of phenotypic traits, selection of marker alleles for molecular breeding, and detection of the founders of genetic effects can be swiftly conducted using Pedimap. However, the power of harnessing the value of Pedimap for making breeding decisions relies on the availability of data for the traits, markers, and genomic sequences. Thus, it is imperative to characterize the breeding germplasms using standard phenomic and genomic characterization procedures such as the assessment of before organized into Pedimap. Thereby, the worldwide breeding programs can benefit from each other to produce improved varieties to meet global challenges.

Genome-wide association study reveals the genetic basis of fiber quality traits in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Fiber quality is an important economic trait of cotton, and its improvement is a major goal of cotton breeding. To better understand the genetic mechanisms responsible for fiber quality traits, we conducted a genome-wide association study to identify and mine fiber-quality-related quantitative trait loci (QTLs) and genes.

RESULTS

In total, 42 single nucleotide polymorphisms (SNPs) and 31 QTLs were identified as being significantly associated with five fiber quality traits. Twenty-five QTLs were identified in previous studies, and six novel QTLs were firstly identified in this study. In the QTL regions, 822 genes were identified and divided into four clusters based on their expression profiles. We also identified two pleiotropic SNPs. The SNP locus i52359Gb was associated with fiber elongation, strength, length and uniformity, while i11316Gh was associated with fiber strength and length. Moreover, these two SNPs were nonsynonymous and located in genes Gh_D09G2376 and Gh_D06G1908, respectively. RT-qPCR analysis revealed that these two genes were preferentially expressed at one or more stages of cotton fiber development, which was consistent with the RNA-seq data. Thus, Gh_D09G2376 and Gh_D06G1908 may be involved in fiber developmental processes.

CONCLUSIONS

The findings of this study provide insights into the genetic bases of fiber quality traits, and the identified QTLs or genes may be applicable in cotton breeding to improve fiber quality.

Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium

BACKGROUND

Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress.

RESULTS

We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium.

CONCLUSION

In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.

Conservation and Divergence in Duplicated Fiber Coexpression Networks Accompanying Domestication of the Polyploid Gossypium hirsutum L

Gossypium hirsutum L. (Upland cotton) has an evolutionary history involving inter-genomic hybridization, polyploidization, and subsequent domestication. We analyzed the developmental dynamics of the cotton fiber transcriptome accompanying domestication using gene coexpression networks for both joint and homoeologous networks. Remarkably, most genes exhibited expression for at least one homoeolog, confirming previous reports of widespread gene usage in cotton fibers. Most coexpression modules comprising the joint network are preserved in each subgenomic network and are enriched for similar biological processes, showing a general preservation of network modular structure for the two co-resident genomes in the polyploid. Interestingly, only one fifth of homoeologs co-occur in the same module when separated, despite similar modular structures between the joint and homoeologous networks. These results suggest that the genome-wide divergence between homoeologous genes is sufficient to separate their co-expression profiles at the intermodular level, despite conservation of intramodular relationships within each subgenome. Most modules exhibit D-homoeolog expression bias, although specific modules do exhibit A-homoeolog bias. Comparisons between wild and domesticated coexpression networks revealed a much tighter and denser network structure in domesticated fiber, as evidenced by its fewer modules, 13-fold increase in the number of development-related module member genes, and the poor preservation of the wild network topology. These results demonstrate the amazing complexity that underlies the domestication of cotton fiber.

VdNPS, a Nonribosomal Peptide Synthetase, Is Involved in Regulating Virulence in Verticillium dahliae

Nonribosomal peptide synthetases (NPS) are known for the biosynthesis of antibiotics, toxins, and siderophore production. They are also a virulence determinant in different phytopathogens. However, until now, the functional characterization of NPS in Verticillium dahliae has not been reported. Deletion of the NPS gene in V. dahliae led to the decrease of conidia, microsclerotia, and pathogenicity. ΔVdNPS strains were tolerant to H₂O₂, and the genes involved in H₂O₂ detoxification, iron/copper transport, and cytoskeleton were differentially expressed in ΔVdNPS. Interestingly, ΔVdNPS strains exhibited hypersensitivity to salicylic acid (SA), and the genes involved in SA hydroxylation were up-regulated in ΔVdNPS compared with wild-type V. dahliae under SA stress. Additionally, during infection, ΔVdNPS induced more pathogenesis-related gene expression, higher reactive oxygen species production, and stronger SA-mediated signaling transduction in host to overcome pathogen. Uncovering the function of VdNPS in pathogenicity could provide a reliable theoretical basis for the development of cultivars with durable resistance against V. dahliae-associated diseases.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

Genome-wide identification and characterization of multiple C2 domains and transmembrane region proteins in Gossypium hirsutum

Background

Multiple C2 domains and transmembrane region proteins (MCTPs) may act as transport mediators of other regulators. Although increased number of MCTPs in higher plants implies their diverse and specific functions in plant growth and development, only a few plant MCTPs have been studied and no study on the MCTPs in cotton has been reported.

Results

In this study, we identified 31 MCTPs in G. hirsutum, which were classified into five subfamilies according to the phylogenetic analysis. GhMCTPs from subfamily V exhibited isoelectric points (pIs) less than 7, whereas GhMCTPs from subfamily I, II, III and IV exhibited pIs more than 7.5, implying their distinct biological functions. In addition, GhMCTPs within subfamily III, IV and V exhibited more diverse physicochemical properties, domain architectures and expression patterns than GhMCTPs within subfamily I and II, suggesting that GhMCTPs within subfamily III, IV and V diverged to perform more diverse and specific functions. Analyses of conserved motifs and pIs indicated that the N-terminus was more divergent than the C-terminus and GhMCTPs’ functional divergence might be mainly contributed by the N-terminus. Furthermore, yeast two-hybrid assay indicated that the N-terminus was responsible to interact with target proteins. Phylogenetic analysis classified multiple N-terminal C2 domains into four subclades, suggesting that these C2 domains performed different molecular functions in mediating the transport of target proteins.

Conclusions

Our systematic characterization of MCTPs in G. hirsutum will provide helpful information to further research GhMCTPs’ molecular roles in mediating other regulators’ transport to coordinate growth and development of various cotton tissues.

PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants

The Plant Promoter Analysis Navigator (PlantPAN; http://PlantPAN.itps.ncku.edu.tw/) is an effective resource for predicting regulatory elements and reconstructing transcriptional regulatory networks for plant genes. In this release (PlantPAN 3.0), 17 230 TFs were collected from 78 plant species. To explore regulatory landscapes, genomic locations of TFBSs have been captured from 662 public ChIP-seq samples using standard data processing. A total of 1 233 999 regulatory linkages were identified from 99 regulatory factors (TFs, histones and other DNA-binding proteins) and their target genes across seven species. Additionally, this new version added 2449 matrices extracted from ChIP-seq peaks for cis-regulatory element prediction. In addition to integrated ChIP-seq data, four major improvements were provided for more comprehensive information of TF binding events, including (i) 1107 experimentally verified TF matrices from the literature, (ii) gene regulation network comparison between two species, (iii) 3D structures of TFs and TF-DNA complexes and (iv) condition-specific co-expression networks of TFs and their target genes extended to four species. The PlantPAN 3.0 can not only be efficiently used to investigate critical cis- and trans-regulatory elements in plant promoters, but also to reconstruct high-confidence relationships among TF–targets under specific conditions.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a ‘SyntenyViewer’ to view genome synteny between different cucurbit species and an ‘RNA-Seq’ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

MaGenDB: a functional genomics hub for Malvaceae plants

Malvaceae is a family of flowering plants containing many economically important plant species including cotton, cacao and durian. Recently, the genomes of several Malvaceae species have been decoded, and many omics data were generated for individual species. However, no integrative database of multiple species, enabling users to jointly compare and analyse relevant data, is available for Malvaceae. Thus, we developed a user-friendly database named MaGenDB (http://magen.whu.edu.cn) as a functional genomics hub for the plant community. We collected the genomes of 13 Malvaceae species, and comprehensively annotated genes from different perspectives including functional RNA/protein element, gene ontology, KEGG orthology, and gene family. We processed 374 sets of diverse omics data with the ENCODE pipelines and integrated them into a customised genome browser, and designed multiple dynamic charts to present gene/RNA/protein-level knowledge such as dynamic expression profiles and functional elements. We also implemented a smart search system for efficiently mining genes. In addition, we constructed a functional comparison system to help comparative analysis between genes on multiple features in one species or across closely related species. This database and associated tools will allow users to quickly retrieve large-scale functional information for biological discovery.

Cotton Duplicated Genes Produced by Polyploidy Show Significantly Elevated and Unbalanced Evolutionary Rates, Overwhelmingly Perturbing Gene Tree Topology

A phylogenetic tree can be used to illustrate the evolutionary relationship between a group of genes, especially duplicated genes, which are sources of genetic innovation and are often a hotspot of research. However, duplicated genes may have complex phylogenetic topologies due to changes in their evolutionary rates. Here, by constructing phylogenetic trees using different methods, we evaluated the phylogenetic relationships of duplicated genes produced by polyploidization in cotton. We found that at least 83.2% of phylogenetic trees did not conform the expected topology. Moreover, cotton homologous gene copy number has little effect on the topology of duplicated genes. Compared with their cacao orthologs, elevated evolutionary rates of cotton genes are responsible for distorted tree topology. Furthermore, as to both branch and site models, we inferred that positive natural selection during the divergence of fiber-development-related MYB genes was likely, and found that the reconstructed tree topology may often overestimate natural selection, as compared to the inference with the expected trees. Therefore, we emphasize the importance of borrowing precious information from gene collinearity in tree construction and evaluation, and have evidence to alert the citation of thousands of previous reports of adaptivity and functional innovation of duplicated genes.

The Gossypium longicalyx Genome as a Resource for Cotton Breeding and Evolution

Cotton is an important crop that has made significant gains in production over the last century. Emerging pests such as the reniform nematode have threatened cotton production. The rare African diploid species Gossypium longicalyx is a wild species that has been used as an important source of reniform nematode immunity. While mapping and breeding efforts have made some strides in transferring this immunity to the cultivated polyploid species, the complexities of interploidal transfer combined with substantial linkage drag have inhibited progress in this area. Moreover, this species shares its most recent common ancestor with the cultivated A-genome diploid cottons, thereby providing insight into the evolution of long, spinnable fiber. Here we report a newly generated de novo genome assembly of G. longicalyx. This high-quality genome leveraged a combination of PacBio long-read technology, Hi-C chromatin conformation capture, and BioNano optical mapping to achieve a chromosome level assembly. The utility of the G. longicalyx genome for understanding reniform immunity and fiber evolution is discussed.

Identification of NHXs in Gossypium species and the positive role of GhNHX1 in salt tolerance

BACKGROUND

Plant Na⁺/H⁺ antiporters (NHXs) are membrane-localized proteins that maintain cellular Na⁺/K⁺ and pH homeostasis. Considerable evidence highlighted the critical roles of NHX family in plant development and salt response; however, NHXs in cotton are rarely studied.

RESULTS

The comprehensive and systematic comparative study of NHXs in three Gossypium species was performed. We identified 12, 12, and 23 putative NHX proteins from G. arboreum, G. raimondii, and G. hirsutum, respectively. Phylogenetic study revealed that repeated polyploidization of Gossypium spp. contributed to the expansion of NHX family. Gene structure analysis showed that cotton NHXs contain many introns, which will lead to alternative splicing and help plants to adapt to high salt concentrations in soil. The expression changes of NHXs indicate the possible differences in the roles of distinct NHXs in salt response. GhNHX1 was proved to be located in the vacuolar system and intensively induced by salt stress in cotton. Silencing of GhNHX1 resulted in enhanced sensitivity of cotton seedlings to high salt concentrations, which suggests that GhNHX1 positively regulates cotton tolerance to salt stress.

CONCLUSION

We characterized the gene structure, phylogenetic relationship, chromosomal location, and expression pattern of NHX genes from G. arboreum, G. raimondii, and G. hirsutum. Our findings indicated that the cotton NHX genes are regulated meticulously and differently at the transcription level with possible alternative splicing. The tolerance of plants to salt stress may rely on the expression level of a particular NHX, rather than the number of NHXs in the genome. This study could provide significant insights into the function of plant NHXs, as well as propose promising candidate genes for breeding salt-resistant cotton cultivars.

Coriander Genomics Database: a genomic, transcriptomic, and metabolic database for coriander

Coriander (Coriandrum sativum L.), also known as cilantro, is a globally important vegetable and spice crop. Its genome and that of carrot are models for studying the evolution of the Apiaceae family. Here, we developed the Coriander Genomics Database (CGDB, http://cgdb.bio2db.com/) to collect, store, and integrate the genomic, transcriptomic, metabolic, functional annotation, and repeat sequence data of coriander and carrot to serve as a central online platform for Apiaceae and other related plants. Using these data sets in the CGDB, we intriguingly found that seven transcription factor (TF) families showed significantly greater numbers of members in the coriander genome than in the carrot genome. The highest ratio of the numbers of MADS TFs between coriander and carrot reached 3.15, followed by those for tubby protein (TUB) and heat shock factors. As a demonstration of CGDB applications, we identified 17 TUB family genes and conducted systematic comparative and evolutionary analyses. RNA-seq data deposited in the CGDB also suggest dose compensation effects of gene expression in coriander. CGDB allows bulk downloading, significance searches, genome browser analyses, and BLAST searches for comparisons between coriander and other plants regarding genomics, gene families, gene collinearity, gene expression, and the metabolome. A detailed user manual and contact information are also available to provide support to the scientific research community and address scientific questions. CGDB will be continuously updated, and new data will be integrated for comparative and functional genomic analysis in Apiaceae and other related plants.

Comprehensive characterization and gene expression patterns of LBD gene family in Gossypium

A comprehensive account of the LBD gene family of Gossypium was provided in this work. Expression analysis and functional characterization revealed that LBD genes might play different roles in G. hirsutum and G. barbadense. The Lateral Organ Boundaries Domain (LBD) proteins comprise a plant-specific transcription factor family, which plays crucial roles in physiological processes of plant growth, development, and stress tolerance. In the present work, a systematical analysis of LBD gene family from two allotetraploid cotton species, G. hirsutum and G. barbadense, together with their genomic donor species, G. arboreum and G. raimondii, was conducted. There were 131, 128, 62, and 68 LBDs identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. The LBD proteins could be classified into two main classes, class I and class II, based on the structure of their lateral organ boundaries domain and traits of phylogenetic tree, and class I was further divided into five subgroups. The gene structure and motif composition analyses conducted in both G. hirsutum and G. barbadense revealed that LBD genes kept relatively conserved within the subfamilies. Synteny analysis suggested that segmental duplication acted as an important mechanism in expansion of the cotton LBD gene family. Cis-element analysis predicated the possible functions of LBD genes. Public RNA-seq data were investigated to analyze the expression patterns of cotton LBD genes in various tissues as well as gene expression under abiotic stress treatments. Furthermore, RT-qPCR results found that GhLBDs had various expression regulation under MeJA treatments. Expression analysis indicated the differential functions of cotton LBD genes in response to abiotic stress and hormones.

Genome-wide expression analysis suggests glutaredoxin genes response to various stresses in cotton

Oxidative stress reflects an imbalance between the systemic manifestation of reactive oxygen species (ROS) and a biological system's ability to readily detoxify the reactive intermediates or to repair the resulting damage. Glutaredoxins (GRXs) are ubiquitous oxidoreductase enzymes involved in diverse cellular processes and play a key role in oxidative stress responsive mechanisms. This study was aimed to explore the structure-function relationship and to provide a framework for functional validation and biochemical characterization of various GRX members. In this study, our analysis revealed the presence of 127 genes encoding GRX proteins in G. hirsutum. A total of 758 genes from two typical monocot and nine dicot species were naturally divided into four classes based on phylogenetic analysis. The classification was supported with organization of conserved protein motifs and sequence logos comparison between cotton, rice and Arabidopsis. Cotton GRX gene family has underwent strong purifying selection with limited functional divergence. A good collinearity was observed in the synteny analysis of four Gossypium species. Majority of cotton GRXs were influenced by various phytohormones and abiotic stress conditions during expression analysis, suggesting an important role of GRX proteins in response to oxidative stress. Cis-regulatory elements, gene enrichments and co-expression network analysis also support their predicted role against various abiotic stresses. Whole genome and segmental duplication were determined to be the two major impetuses for the expansion of gene numbers during the evolution. The identification of GRX genes showing differential expression in specific tissues or in response to environmental stimuli provides a new avenue for in-depth characterization of selected genes of importance. This study will further broaden our insights into the evolution and functional elucidation of GRX gene family in cotton.

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.

Comprehensive analysis of the Gossypium hirsutum L. respiratory burst oxidase homolog (Ghrboh) gene family

BACKGROUND

Plant NADPH oxidase (NOX), also known as respiratory burst oxidase homolog (rboh), encoded by the rboh gene, is a key enzyme in the reactive oxygen species (ROS) metabolic network. It catalyzes the formation of the superoxide anion (O2•-), a type of ROS. In recent years, various studies had shown that members of the plant rboh gene family were involved in plant growth and developmental processes as well as in biotic and abiotic stress responses, but little is known about its functional role in upland cotton.

RESULTS

In the present study, 26 putative Ghrboh genes were identified and characterized. They were phylogenetically classified into six subfamilies and distributed at different densities across 18 of the 26 chromosomes or scaffolds. Their exon-intron structures, conserved domains, synteny and collinearity, gene family evolution, regulation mediated by cis-acting elements and microRNAs (miRNAs) were predicted and analyzed. Additionally, expression profiles of Ghrboh gene family were analyzed in different tissues/organs and at different developmental stages and under different abiotic stresses, using RNA-Seq data and real-time PCR. These profiling studies indicated that the Ghrboh genes exhibited temporal and spatial specificity with respect to expression, and might play important roles in cotton development and in stress tolerance through modulating NOX-dependent ROS induction and other signaling pathways.

CONCLUSIONS

This comprehensive analysis of the characteristics of the Ghrboh gene family determined features such as sequence, synteny and collinearity, phylogenetic and evolutionary relationship, expression patterns, and cis-element- and miRNA-mediated regulation of gene expression. Our results will provide valuable information to help with further gene cloning, evolutionary analysis, and biological function analysis of cotton rbohs.

Strengthening of banana breeding through data digitalization

Improvement of edible bananas (a triploid and sterile crop) through conventional breeding is a challenging task owing to its recalcitrant nature for seed set, prolonged crop duration. In addition, the need of huge man power at different stages of progeny development and evaluation often leads to mislabeling, poor data management and loss of vital data. All this can be overcome by the application of advanced information technology source. This ensured secure and efficient data management such as storage, retrieval and data analysis and further could assist in tracking the breeding status in real time. Thus, a user-friendly web-based banana breeding tracker (BBT) has been developed using MySQL database with Hypertext Preprocessor (PHP). This BBT works on all operating systems with access to multiple users from anywhere at any time. Quick responsive (QR) code labels can be generated by the tracker, which can be decoded using QR scanner. Also for each and every updated progress in breeding stages, a new QR code can be generated, which in turn reduce labeling errors. Moreover, the tracker has additional tools to search, sort and filter the data from the data sets for efficient retrieval and analysis. This tracker is being upgraded with phenotypic and genotypic data that will be made available in the public domain for hastening the banana improvement program.

Identification of chlorprothixene as a potential drug that induces apoptosis and autophagic cell death in acute myeloid leukemia cells

Although acute myeloid leukemia (AML) is a highly heterogeneous malignance, the common molecular mechanisms shared by different AML subtypes play critical roles in AML development. It is possible to identify new drugs that are effective for various AML subtypes based on the common molecular mechanisms. Therefore, we developed a hypothesis-driven bioinformatic drug screening framework by integrating multiple omics data. In this study, we identified that chlorprothixene, a dopamine receptor antagonist, could effectively inhibit growth of AML cells from different subtypes. RNA-seq analysis suggested that chlorprothixene perturbed a series of crucial biological processes such as cell cycle, apoptosis, and autophagy in AML cells. Further investigations indicated that chlorprothixene could induce both apoptosis and autophagy in AML cells, and apoptosis and autophagy could act as partners to induce cell death cooperatively. Remarkably, chlorprothixene was found to inhibit tumor growth and induce in situ leukemic cell apoptosis in the murine xenograft model. Furthermore, chlorprothixene treatment could reduce the level of oncofusion proteins PML-RARα and AML1-ETO, thus elevate the expression of apoptosis-related genes, and lead to AML cell death. Our results provided new insights for drug repositioning of AML therapy and confirmed that chlorprothixene might be a potential candidate for treatment of different subtypes of AML by reducing expression of oncofusion proteins. DATABASE: RNA-seq data are available in GEO database under the accession number GSE124316.

Multi-strategic RNA-seq analysis reveals a high-resolution transcriptional landscape in cotton

Cotton is an important natural fiber crop, however, its comprehensive and high-resolution gene map is lacking. Here we integrate four complementary high-throughput techniques, including Pacbio long read Iso-seq, strand-specific RNA-seq, CAGE-seq, and PolyA-seq, to systematically explore the transcription landscape across 16 tissues or different organ types in Gossypium arboreum. We devise a computational pipeline, named IGIA, to reconstruct accurate gene structures from the integrated data. Our results reveal a dynamic and diverse transcriptional map in cotton: tissue-specific gene expression, alternative usage of TSSs and polyadenylation sites, hotspot of alternative splicing, and transcriptional read-through. These regulated events affect many genes in various aspects such as gain or loss of functional RNA motifs and protein domains, fine-tuning of DNA binding activity, and co-regulation for genes in the same complex or pathway. The methods and findings provide valuable resources for further functional genomic studies such as understanding natural SNP variations for plant community.

Systematic Comparisons of Positively Selected Genes between Gossypium arboreum and Gossypium raimondii Genomes

Background:

Studies of Positively Selected Genes (PSGs) in microorganisms and mammals have provided insights into the dynamics of genome evolution and the genetic basis of differences between species by using whole genome-wide scans. Systematic investigations and comparisons of PSGs in plants, however, are still limited.

Objective:

A systematic comparison of PSGs between the genomes of two cotton species, Gossypium arboreum (G. arboreum) and G. raimondii, will give the key answer for revealing molecular evolutionary differences in plants.

Methods:

Genome sequences of G. arboreum and G. raimondii were compared, including Whole Genome Duplication (WGD) events and genomic features such as gene number, gene length, codon bias index, evolutionary rate, number of expressed genes, and retention of duplicated copies.

Results:

Unlike the PSGs in G. raimondii, G. arboreum comprised more PSGs, smaller gene size and fewer expressed gene. In addition, the PSGs evolved at a higher rate of synonymous substitutions, but were subjected to lower selection pressure. The PSGs in G. arboreum were also retained with a lower number of duplicate gene copies than G. raimondii after a single WGD event involving Gossypium.

Conclusion:

These data indicate that PSGs in G. arboreum and G. raimondii differ not only in Ka/Ks, but also in their evolutionary, structural, and expression properties, indicating that divergence of G. arboreum and G. raimondii was associated with differences in PSGs in terms of evolutionary rates, gene length, expression patterns, and WGD retention in Gossypium.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Gossypium hirsutum

Background

Many BURP domain-containing proteins, which are unique to plants, have been identified. They performed diverse functions in plant development and the stress response. To date, only a few BURP domain-containing genes have been studied, and no comprehensive analysis of the gene family in cotton has been reported.

Results

In this study, 18, 17 and 30 putative BURP genes were identified in G. raimondii (D₅), G. arboreum (A₂) and G. hirsutum (AD₁), respectively. These BURP genes were phylogenetically classified into eight subfamilies, which were confirmed by analyses of gene structures, motifs and protein domains. The uneven distribution of BURPs in chromosomes and gene duplication analysis indicated that segmental duplication might be the main driving force of the GhBURP family expansion. Promoter regions of all GhBURPs contained at least one putative stress-related cis-elements. Analysis of transcriptomic data and qRT-PCR showed that GhBURPs showed different expression patterns in different organs, and all of them, especially the members of the RD22-like subfamily, could be induced by different stresses, such as abscisic acid (ABA) and salicylic acid (SA), which indicated that the GhBURPs may performed important functions in cotton’s responses to various abiotic stresses.

Conclusions

Our study comprehensively analyzed BURP genes in G. hirsutum, providing insight into the functions of GhBURPs in cotton development and adaptation to stresses.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5948-y) contains supplementary material, which is available to authorized users.

Combining genome-wide and transcriptome-wide analyses reveal the evolutionary conservation and functional diversity of aquaporins in cotton

BACKGROUND

Aquaporins (AQPs) are integral membrane proteins from a larger family of major intrinsic proteins (MIPs) and function in a huge variety of processes such as water transport, plant growth and stress response. The availability of the whole-genome data of different cotton species allows us to study systematic evolution and function of cotton AQPs on a genome-wide level.

RESULTS

Here, a total of 53, 58, 113 and 111 AQP genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. A comprehensive analysis of cotton AQPs, involved in exon/intron structure, functional domains, phylogenetic relationships and gene duplications, divided these AQPs into five subfamilies (PIP, NIP, SIP, TIP and XIP). Comparative genome analysis among 30 species from algae to angiosperm as well as common tandem duplication events in 24 well-studied plants further revealed the evolutionary conservation of AQP family in the organism kingdom. Combining transcriptome analysis and Quantitative Real-time PCR (qRT-PCR) verification, most AQPs exhibited tissue-specific expression patterns both in G. raimondii and G. hirsutum. Meanwhile, a bias of time to peak expression of several AQPs was also detected after treating G. davidsonii and G. hirsutum with 200 mM NaCl. It is interesting that both PIP1;4 h/i/j and PIP2;2a/e showed the highly conserved tandem structure, but differentially contributed to tissue development and stress response in different cotton species.

CONCLUSIONS

These results demonstrated that cotton AQPs were structural conservation while experienced the functional differentiation during the process of evolution and domestication. This study will further broaden our insights into the evolution and functional elucidation of AQP gene family in cotton.

A genomic resource derived from the integration of genome sequences, expressed transcripts and genetic markers in ramie

Background

The redundancy of genomic resources, including transcript and molecular markers, and their uncertain position in the genome have dramatically hindered the study of traits in ramie, an important natural fiber crop.

Results

We obtained a high-quality transcriptome consisting of 30,591 non-redundant transcripts using single-molecule long-read sequencing and proposed it as a universal ramie transcriptome. Additionally, 55,882 single nucleotide polymorphisms (SNPs) were identified and a high-density genetic map was developed. Based on this genetic map, 181.7 Mb ramie genome sequences were assembled into 14 chromosomes. For the convenient use of these resources, 29,286 (~ 95.7%) of the transcripts and all 55,882 SNPs, along with 1827 previously reported sequence repeat markers (SSRs), were mapped into the ramie genome, and 22,343 (~ 73.0%) transcripts, 50,154 (~ 89.7%) SNPs, and 1466 (~ 80.3%) SSRs were assigned to a specific location in the corresponding chromosome.

Conclusion

This is the first study to characterize the ramie transcriptome by long-read sequencing, and the substantial number of transcripts of significant length obtained will accelerate our understanding of ramie growth and development. This integration of genome sequences, expressed transcripts, and genetic markers will provide an extremely useful resource for genetic, molecular, and breeding studies of ramie.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5878-8) contains supplementary material, which is available to authorized users.

Genes regulating gland development in the cotton plant

In seeds and other parts of cultivated, tetraploid cotton (Gossypium hirsutum L.), multicellular groups of cells lysigenously form dark glands containing toxic terpenoids such as gossypol that defend the plant against pests and pathogens. Using RNA-seq analysis of embryos from near-isogenic glanded (Gl2 Gl2 Gl3 Gl3 ) versus glandless (gl2 gl2 gl3 gl3 ) plants, we identified 33 genes that expressed exclusively or at higher levels in embryos just prior to gland formation in glanded plants. Virus-induced gene silencing against three gene pairs led to significant reductions in the number of glands in the leaves, and significantly lower levels of gossypol and related terpenoids. These genes encode transcription factors and have been designated the 'Cotton Gland Formation' (CGF) genes. No sequence differences were found between glanded and glandless cotton for CGF1 and CGF2 gene pairs. The glandless cotton has a transposon insertion within the coding sequence of the GoPGF (synonym CGF3) gene of the A subgenome and extensive mutations in the promoter of D subgenome homeolog. Overexpression of GoPGF (synonym CGF3) led to a dramatic increase in gossypol and related terpenoids in cultured cells, whereas CRISPR/Cas9 knockout of GoPGF (synonym CGF3) genes resulted in glandless phenotype. Taken collectively, the results show that the GoPGF (synonym CGF3) gene plays a critical role in the formation of glands in the cotton plant. Seed-specific silencing of CGF genes, either individually or in combination, could eliminate glands, thus gossypol, from the cottonseed to render it safe as food or feed for monogastrics.

Cumulative and different genetic effects contributed to yield heterosis using maternal and paternal backcross populations in Upland cotton

Heterosis has been utilized in commercial production, but the heterosis mechanism has remained vague. Hybrid cotton is suitable to dissect the heterosis mechanism. In order to explore the genetic basis of heterosis in Upland cotton, we generated paternal and maternal backcross (BC/P and BC/M) populations. Data for yield and yield-component traits were collected over 2 years in three replicated BC/P field trials and four replicated BC/M field trials. At single-locus level, 26 and 27 QTLs were identified in BC/P and BC/M populations, respectively. Six QTLs shared in both BC populations. A total of 27 heterotic loci were detected. Partial dominant and over-dominant QTLs mainly determined yield heterosis in the BC/P and BC/M populations. QTLs for different traits displayed varied genetic effects in two BC populations. Eleven heterotic loci overlapped with QTLs but no common heterotic locus was detected in both BC populations. We resolved the 333 kb (48 genes) and 516 kb (25 genes) physical intervals based on 16 QTL clusters and 35 common QTLs, respectively, in more than one environment or population. We also identified 189 epistatic QTLs and a number of QTL × environment interactions in two BC populations and the corresponding MPH datasets. The results indicated that cumulative effects contributed to yield heterosis in Upland cotton, including epistasis, QTL × environment interaction, additive, partial dominance and over-dominance.

Evolution of PEPC gene family in Gossypium reveals functional diversification and GhPEPC genes responding to abiotic stresses

Phosphoenolpyruvate carboxylase (PEPC) family genes play important roles in regulating plant growth and abiotic stress response. Based on the sequenced Gossypium genomes, we performed comprehensive analysis of PEPC homolog genes in cotton, which six, six, eleven and ten PEPC genes were identified in Gossypium arboreum (A₂), G. raimondii (D₅), G. hirsutum (AD₁) and G. barbadense (AD₂), respectively. These genes were divided into six subgroups: PEPC-i, PEPC-ii, PEPC-iii, PEPC-iv, PEPC-v and PEPC-vi; PEPC genes in each subgroup displayed conserved gene structure and motifs. Segmental duplication and whole genome duplication (WGD) events yielded the expansion of PEPC genes. Expression assays showed that the duplicated PEPC genes displayed diverse expression patterns, indicating that they experienced functional divergence. Of which, genes in PEPC-iv subgroup played crucial role for substrate distribution in cottonseed. Cis-elements, putative miRNAs and expression analyses showed that GhPEPC homologs might respond to abiotic stresses, expression levels of GhPEPC1 and GhPEPC2/GhPEPC2D genes were larger induced than other GhPEPC genes under cold, heat, salt, and drought stresses, indicating the crucial roles in abiotic stresses response. Present study serves new information to decipher the evolution and function of PEPC genes in Gossypium.

Genome-Wide Identification, Sequence Variation, and Expression of the Glycerol-3-Phosphate Acyltransferase (GPAT) Gene Family in Gossypium

Cotton is an economically important crop grown for natural fiber and seed oil production. Cottonseed oil ranks third after soybean oil and colza oil in terms of edible oilseed tonnage worldwide. Glycerol-3-phosphate acyltransferase (GPAT) genes encode enzymes involved in triacylglycerol biosynthesis in plants. In the present study, 85 predicted GPAT genes were identified from the published genome data in Gossypium. Among them, 14, 16, 28, and 27 GPAT homologs were identified in G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. Phylogenetic analysis revealed that a total of 108 GPAT genes from cotton, Arabidopsis and cacao could be classified into three groups. Furthermore, through comparison, the gene structure analyses indicated that GPAT genes from the same group were highly conserved between Arabidopsis and cotton. Segmental duplication could be the major driver for GPAT gene family expansion in the four cotton species above. Expression patterns of GhGPAT genes were diverse in different tissues. Most GhGPAT genes were induced or suppressed after salt or cold stress in Upland cotton. Eight GhGPAT genes were co-localized with oil and protein quantitative trait locus (QTL) regions. Thirty-two single nucleotide polymorphisms (SNPs) were detected from 12 GhGPAT genes, sixteen of which in nine GhGPAT genes were classified as synonymous, and sixteen SNPs in ten GhGPAT genes non-synonymous. Two SNP markers of the GhGPAT16 and GhGPAT26 genes were significantly correlated with cotton oil content in one of the three field tests. This study shed lights on the molecular evolutionary properties of GPAT genes in cotton, and provided reference for improvement of cotton response to abiotic stress and the genetic improvement of cotton oil content.

Co-Expression Network Analysis and Hub Gene Selection for High-Quality Fiber in Upland Cotton (Gossypium hirsutum) Using RNA Sequencing Analysis

Upland cotton (Gossypium hirsutum) is grown for its elite fiber. Understanding differential gene expression patterns during fiber development will help to identify genes associated with fiber quality. In this study, we used two recombinant inbred lines (RILs) differing in fiber quality derived from an intra-hirsutum population to explore expression profiling differences and identify genes associated with high-quality fiber or specific fiber-development stages using RNA sequencing. Overall, 72/27, 1137/1584, 437/393, 1019/184, and 2555/1479 differentially expressed genes were up-/down-regulated in an elite fiber line (L1) relative to a poor-quality fiber line (L2) at 10, 15, 20, 25, and 30 days post-anthesis, respectively. Three-hundred sixty-three differentially expressed genes (DEGs) between two lines were colocalized in fiber strength (FS) quantitative trait loci (QTL). Short Time-series Expression Miner (STEM) analysis discriminated seven expression profiles; gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation were performed to identify difference in function between genes unique to L1 and L2. Co-expression network analysis detected five modules highly associated with specific fiber-development stages, especially for high-quality fiber tissues. The hub genes in each module were identified by weighted gene co-expression network analysis. Hub genes encoding actin 1, Rho GTPase-activating protein with PAK-box, TPX2 protein, bHLH transcription factor, and leucine-rich repeat receptor-like protein kinase were identified. Correlation networks revealed considerable interaction among the hub genes, transcription factors, and other genes.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.

PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants

The Plant Promoter Analysis Navigator (PlantPAN; http://PlantPAN.itps.ncku.edu.tw/) is an effective resource for predicting regulatory elements and reconstructing transcriptional regulatory networks for plant genes. In this release (PlantPAN 3.0), 17 230 TFs were collected from 78 plant species. To explore regulatory landscapes, genomic locations of TFBSs have been captured from 662 public ChIP-seq samples using standard data processing. A total of 1 233 999 regulatory linkages were identified from 99 regulatory factors (TFs, histones and other DNA-binding proteins) and their target genes across seven species. Additionally, this new version added 2449 matrices extracted from ChIP-seq peaks for cis-regulatory element prediction. In addition to integrated ChIP-seq data, four major improvements were provided for more comprehensive information of TF binding events, including (i) 1107 experimentally verified TF matrices from the literature, (ii) gene regulation network comparison between two species, (iii) 3D structures of TFs and TF-DNA complexes and (iv) condition-specific co-expression networks of TFs and their target genes extended to four species. The PlantPAN 3.0 can not only be efficiently used to investigate critical cis- and trans-regulatory elements in plant promoters, but also to reconstruct high-confidence relationships among TF-targets under specific conditions.

Target of Rapamycin (TOR) Regulates the Expression of lncRNAs in Response to Abiotic Stresses in Cotton

TOR (Target of Rapamycin) kinase is an evolutionarily conserved protein kinase, which integrates stress-related cues with growth and metabolic outputs. Long non-coding RNAs (lncRNAs) play a vital role in the regulation of eukaryotic genes. However, little is known about TOR's function in regulating the expression of lncRNAs in plants. In this study, four putative homologous genes encoding the TOR protein were identified by utilizing the recently completed cotton genome. Pharmacological experiments with TOR inhibitor AZD8055 and on silencing GhTOR genes resulted in obvious cotton growth retardation, indicating the conserved role of TOR in plant growth. The expression pattern analyses in different tissues reveal that TOR may play a role in root development, and the transcript levels of TOR genes were changed under different stress conditions. Importantly, we found TOR may be a key player in regulating the expression of long non-coding RNAs (lncRNAs). A total of 10,315 lncRNAs were discovered in cotton seedlings, 90.7% of which were long intergenic ncRNAs. Moreover, we identified the differentially expressed lncRNAs, of which 296 were significantly upregulated and 105 were downregulated in TOR inactivated plants. GO and KEGG analyses of differentially expressed lncRNA neighboring genes reveal that these differentially expressed lncRNA-targeted genes are involved in many life processes, including stress response, glutathione, and ribosomes in cotton. A series of differentially expressed lncRNAs potentially involved in plant stress response was identified under TOR inhibition. Collectively, these results suggest that cotton TOR proteins may directly modulate the expression of putative stress-related lncRNAs and eventually play a potential role in the cotton stress response.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Fiber Quality Improvement in Upland Cotton (Gossypium hirsutum L.): Quantitative Trait Loci Mapping and Marker Assisted Selection Application

Genetic improvement in fiber quality is one of the main challenges for cotton breeders. Fiber quality traits are controlled by multiple genes and are classified as complex quantitative traits, with a negative relationship with yield potential, so the genetic gain is low in traditional genetic improvement by phenotypic selection. The availability of Gossypium genomic sequences facilitates the development of high-throughput molecular markers, quantitative trait loci (QTL) fine mapping and gene identification, which helps us to validate candidate genes and to use marker assisted selection (MAS) on fiber quality in breeding programs. Based on developments of high density linkage maps, QTLs fine mapping, marker selection and omics, we have performed trait dissection on fiber quality traits in diverse populations of upland cotton. QTL mapping combined with multi-omics approaches such as, RNA sequencing datasets to identify differentially expressed genes have benefited the improvement of fiber quality. In this review, we discuss the application of molecular markers, QTL mapping and MAS for fiber quality improvement in upland cotton.

Tripal MapViewer: A tool for interactive visualization and comparison of genetic maps

Tripal is an open-source, resource-efficient toolkit for construction of genomic, genetic and breeding databases. It facilitates development of biological websites by providing tools to integrate and display biological data using the generic database schema, Chado, together with Drupal, a popular website creation and content management system. Tripal MapViewer is a new interactive tool for visualizing genetic map data. Developed as a Tripal replacement for Comparative Map Viewer (CMap), it enables visualization of entire maps or linkage groups and features such as molecular markers, quantitative trait loci (QTLs) and heritable phenotypic markers. It also provides graphical comparison of maps sharing the same markers as well as dot plot and correspondence matrices. MapViewer integrates directly with the Tripal application programming interface framework, improving data searching capability and providing a more seamless experience for site visitors. The Tripal MapViewer interface can be integrated in any Tripal map page and linked from any Tripal page for markers, QTLs, heritable morphological markers or genes. Configuration of the display is available through a control panel and the administration interface. The administration interface also allows configuration of the custom database query for building materialized views, providing better performance and flexibility in the way data is stored in the Chado database schema. MapViewer is implemented with the D3.js technology and is currently being used at the Genome Database for Rosaceae (https://www.rosaceae.org), CottonGen (https://www.cottongen.org), Citrus Genome Database (https://citrusgenomedb.org), Vaccinium Genome Database (https://www.vaccinium.org) and Cool Season Food Legume Database (https://www.coolseasonfoodlegume.org). It is also currently in development on the Hardwood Genomics Web (https://hardwoodgenomics.org) and TreeGenes (https://treegenesdb.org). Database URL: https://gitlab.com/mainlabwsu/tripal_map.

Genome-wide characterization and phylogenetic analysis of GSK gene family in three species of cotton: evidence for a role of some GSKs in fiber development and responses to stress

BACKGROUND

The glycogen synthase kinase 3/shaggy kinase (GSK3) is a serine/threonine kinase with important roles in animals. Although GSK3 genes have been studied for more than 30 years, plant GSK genes have been studied only since the last decade. Previous research has confirmed that plant GSK genes are involved in diverse processes, including floral development, brassinosteroid signaling, and responses to abiotic stresses.

RESULT

In this study, 20, 15 (including 5 different transcripts) and 10 GSK genes were identified in G. hirsutum, G. raimondii and G. arboreum, respectively. A total of 65 genes from Arabidopsis, rice, and cotton were classified into 4 clades. High similarities were found in GSK3 protein sequences, conserved motifs, and gene structures, as well as good concordance in gene pairwise comparisons (G. hirsutum vs. G. arboreum, G. hirsutum vs. G. raimondii, and G. arboreum vs. G. raimondii) were observed. Whole genome duplication (WGD) within At and Dt sub-genomes has been central to the expansion of the GSK gene family. Furthermore, GhSK genes showed diverse expression patterns in various tissues. Additionally, the expression profiles of GhSKs under different stress treatments demonstrated that many are stress-responsive genes. However, none were induced by brassinolide treatment. Finally, nine co-expression sub-networks were observed for GhSKs and the functional annotations of these genes suggested that some GhSKs might be involved in cotton fiber development.

CONCLUSION

In this present work, we identified 45 GSK genes from three cotton species, which were divided into four clades. The gene features, muti-alignment, conversed motifs, and syntenic blocks indicate that they have been highly conserved during evolution. Whole genome duplication was determined to be the dominant factor for GSK gene family expansion. The analysis of co-expressed sub-networks and tissue-specific expression profiles suggested functions of GhSKs during fiber development. Moreover, their different responses to various abiotic stresses indicated great functional diversity amongst the GhSKs. Briefly, data presented herein may serve as the basis for future functional studies of GhSKs.

Zinc Finger-Homeodomain Transcriptional Factors (ZHDs) in Upland Cotton (Gossypium hirsutum): Genome-Wide Identification and Expression Analysis in Fiber Development

Zinc finger-homeodomain (ZHD) genes encode a family of plant-specific transcription factors that not only participate in the regulation of plant growth and development but also play an important role in the response to abiotic stress. The ZHD gene family has been studied in several model plants, including Solanum lycopersicum, Zea mays, Oryza sativa, and Arabidopsis thaliana. However, a comprehensive study of the genes of the ZHD family and their roles in fiber development and pigmentation in upland cotton has not been completed. To address this gap, we selected a brown fiber cultivar for our study; brown color in cotton is one of the most desired colors in the textile industry. The natural colored fibers require less processing and little dying, thereby eliminating dye costs and chemical residues. Using bioinformatics approaches, we identified 37 GhZHD genes from Gossypium hirsutum and then divided these genes into seven groups based on their phylogeny. The GhZHD genes were mostly conserved in each subfamily with minor variations in motif distribution and gene structure. These genes were largely distributed on 19 of the 26 upland cotton chromosomes. Among the Gossypium genomes, the paralogs and orthologs of the GhZHD genes were identified and further characterized. Furthermore, among the paralogs, we observed that the ZHD family duplications in Gossypium genomes (G. hirsutum, G. arboreum, and G. raimondii) were probably derived from segmental duplication or genome-wide duplication (GWD) events. Through a combination of qRT-PCR and proanthocyanidins (PA) accumulation analyses in brown cotton fibers, we concluded that the candidate genes involved in early fiber development and fiber pigment synthesis include the following: GhZHD29, GhZHD35, GhZHD30, GhZHD31, GhZHD11, GhZHD27, GhZHD18, GhZHD15, GhZHD16, GhZHD22, GhZHD6, GhZHD33, GhZHD13, GhZHD5, and GhZHD23. This study delivers insights into the evolution of the GhZHD genes in brown cotton, serves as a valuable resource for further studies, and identifies the conditions necessary for improving the quality of brown cotton fiber.

Genome-wide identification and expression analyses of the pectate lyase (PEL) gene family in cotton (Gossypium hirsutum L.)

Background

Pectin is a major component and structural polysaccharide of the primary cell walls and middle lamella of higher plants. Pectate lyase (PEL, EC 4.2.2.2), a cell wall modification enzyme, degrades de-esterified pectin for cell wall loosening, remodeling and rearrangement. Nevertheless, there have been few studies on PEL genes and no comprehensive analysis of the PEL gene family in cotton.

Results

We identified 53, 42 and 83 putative PEL genes in Gossypium raimondii (D5), Gossypium arboreum (A2), and Gossypium hirsutum (AD1), respectively. These PEL genes were classified into five subfamilies (I-V). Members from the same subfamilies showed relatively conserved gene structures, motifs and protein domains. An analysis of gene chromosomal locations and gene duplication revealed that segmental duplication likely contributed to the expansion of the GhPELs. The 2000 bp upstream sequences of all the GhPELs contained auxin response elements. A transcriptomic data analysis showed that 62 GhPELs were expressed in various tissues. Notably, most (29/32) GhPELs of subfamily IV were preferentially expressed in the stamen, and five GhPELs of subfamily V were prominently expressed at the fiber elongation stage. In addition, qRT-PCR analysis revealed the expression characteristics of 24 GhPELs in four pollen developmental stages and significantly different expression of some GhPELs between long- and short-fiber cultivars. Moreover, some members were responsive to IAA treatment. The results indicate that GhPELs play significant and functionally diverse roles in the development of different tissues.

Conclusions

In this study, we comprehensively analyzed PELs in G. hirsutum, providing a foundation to better understand the functions of GhPELs in different tissues and pathways, especially in pollen, fiber and the auxin signaling pathway.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5047-5) contains supplementary material, which is available to authorized users.

Genome-Wide Analysis Elucidates the Role of CONSTANS-like Genes in Stress Responses of Cotton

The CONSTANS (CO)-like gene family has been well studied for its role in the regulation of plant flowering time. However, their role remains poorly understood in cotton. To better understand the possible roles of CO-like in cotton, we performed a comprehensive genome-wide analysis of CO-like genes in cotton. Phylogenetic tree analysis showed that CO-like genes naturally clustered into three groups. Segmental duplication and whole genome duplication (WGD), which occurred before polyploidy, were important contributors to its expansion within the At (“t” indicates tetraploid) and Dt subgenomes, particularly in Group III. Long-terminal repeat retroelements were identified as the main transposable elements accompanying 18 genes. The genotype of GhCOL12_Dt displayed low diversity; it was a candidate involved in domestication. Selection pressure analyses indicated that relaxed purifying selection might have provided the main impetus during the evolution of CO-like genes in upland cotton. In addition, the high expression in the torus and calycle indicated that CO-like genes might affect flowering. The genes from Group II, and those from Group III involved in segmental duplication or WGD, might play important roles in response to drought and salt stress. Overall, this comprehensive genome-wide study of the CO-like gene family would facilitate further detailed studies in cotton.

Genome-wide identification and characterization of LRR-RLKs reveal functional conservation of the SIF subfamily in cotton (Gossypium hirsutum)

BACKGROUND

As one of the largest subfamilies of the receptor-like protein kinases (RLKs) in plants, Leucine Rich Repeats-RLKs (LRR-RLKs) are involved in many critical biological processes including growth, development and stress responses in addition to various physiological roles. Arabidopsis contains 234 LRR-RLKs, and four members of Stress Induced Factor (SIF) subfamily (AtSIF1-AtSIF4) which are involved in abiotic and biotic stress responses. Herein, we aimed at identification and functional characterization of SIF subfamily in cultivated tetraploid cotton Gossypium hirsutum.

RESULTS

Genome-wide analysis of cotton LRR-RLK gene family identified 543 members and phylogenetic analysis led to the identification of 6 cotton LRR-RLKs with high homology to Arabidopsis SIFs. Of the six SIF homologs, GhSIF1 is highly conserved exhibiting 46-47% of homology with AtSIF subfamily in amino acid sequence. The GhSIF1 was transiently silenced using Virus-Induced Gene Silencing system specifically targeting the 3' Untranslated Region. The transiently silenced cotton seedlings showed enhanced salt tolerance compared to the control plants. Further, the transiently silenced plants showed better growth, lower electrolyte leakage, and higher chlorophyll and biomass contents.

CONCLUSIONS

Overall, 543 LRR-RLK genes were identified using genome-wide analysis in cultivated tetraploid cotton G. hirsutum. The present investigation also demonstrated the conserved salt tolerance function of SIF family member in cotton. The GhSIF1 gene can be knocked out using genome editing technologies to improve salt tolerance in cotton.

Identification of Introgressed Alleles Conferring High Fiber Quality Derived From Gossypium barbadense L. in Secondary Mapping Populations of G. hirsutum L

The improvement of fiber quality is an essential goal in cotton breeding. In our previous studies, several quantitative trait loci (QTLs) contributing to improved fiber quality were identified in different introgressed chromosomal regions from Sea Island cotton (Gossypium barbadense L.) in a primary introgression population (Pop. A) of upland cotton (G. hirsutum L.). In the present study, to finely map introgressed major QTLs and accurately dissect the genetic contribution of the target introgressed chromosomal segments, we backcrossed two selected recombinant inbred lines (RILs) that presented desirable high fiber quality with their high lint-yielding recurrent parent to ultimately develop two secondary mapping populations (Pop. B and Pop. C). Totals of 20 and 27 QTLs for fiber quality were detected in Pop. B and Pop. C, respectively, including four and five for fiber length, four and eight for fiber micronaire, two and four for fiber uniformity, five and four for fiber elongation, and six and four for fiber strength, respectively. Two QTLs for lint percentage were detected only in Pop. C. In addition, seven stable QTLs were identified, including two for both fiber length and fiber strength and three for fiber elongation. Five QTL clusters for fiber quality were identified in the introgressed chromosomal regions, and negative effects of these chromosomal regions on lint percentage (a major lint yield parameter) were not observed. Candidate genes with a QTL-cluster associated with fiber strength and fiber length in the introgressed region of Chr.7 were further identified. The results may be helpful for revealing the genetic basis of superior fiber quality contributed by introgressed alleles from G. barbadense. Possible strategies involving marker-assisted selection (MAS) for simultaneously improving upland cotton fiber quality and lint yield in breeding programs was also discussed.

RNA interference of GhPEPC2 enhanced seed oil accumulation and salt tolerance in Upland cotton

Phosphoenolpyruvate carboxylase (PEPCase) mainly produces oxaloacetic acid for tricarboxylic acid (TCA) cycle. Here we reported that GhPEPC2 silencing with PEPC2-RNAi vector could regulate oil and protein accumulation in cottonseeds. In GhPEPC2 transgenic plants, PEPCase activities in immature embryos were significantly reduced, and the oil content in seed kernel was increased 7.3 percentages, whereas total proteins decreased 5.65 percentages. Compared to wild type, agronomical traits of transgenic plant were obviously unaffected. Furthermore, gene expression profile of GhPEPC2 transgenic seeds were investigated using RNA-seq, most lipid synthesis related genes were up-regulated, but amino acid metabolic related genes were down-regulated. In addition, the GhPEPC2 transgenic cotton seedlings were stressed using sodium salts at seedling stage, and the salt tolerance was significantly enhanced. Our observations of GhPEPC2 in cotton would shade light on understanding the regulation of oil content, protein accumulation and salt tolerance enhancement in other plants.

Genome-Wide Study of YABBY Genes in Upland Cotton and Their Expression Patterns under Different Stresses

Members of the YABBY gene family, a small plant-specific family of genes, have been proposed to function in specifying abaxial cell fate. Although to date little has been learned about cotton YABBY genes, completion of the cotton genome enables a comprehensive genome-wide analysis of YABBY genes in cotton. Here, a total of 12, 12, and 23 YABBY genes were identified in Gossypium arboreum (2n = 26, A₂), G. raimondii (2n = 26, D₅), and G. hirsutum (2n = 4x = 52, [AD]_t), respectively. Sequence analysis showed that the N-terminal zinc-finger and C-terminal YABBY domains in YABBY proteins are highly conserved among cotton, Arabidopsis, and rice. Eighty-five genes from eight sequenced species naturally clustered into five groups, and the YAB2-like group could be divided into three sub-groups, indicating that YABBYs are highly conserved among the examined species. Orthologs from the At and Dt sub-genomes (where “t” indicates tetraploid) showed good collinearity, indicating that YABBY loci are highly conserved between these two sub-genomes. Whole-genome duplication was the primary cause of upland cotton YABBY gene expansion, segmental duplication played important roles in YABBY gene expansion within the At and Dt sub-genomes, and the YAB5-like group was mainly generated by segmental duplication. The long-terminal repeat retroelements Copia and Gypsy were identified as major transposable elements accompanying the appearance of duplicated YABBY genes, suggesting that transposable element expansion might be involved in gene duplication. Selection pressure analyses using PAML revealed that relaxed purifying selection might be the main impetus during evolution of YABBY genes in the examined species. Furthermore, exon/intron pattern and motif analyses indicated that genes within the same group were significantly conserved between Arabidopsis and cotton. In addition, the expression patterns in different tissues suggest that YABBY proteins may play roles in ovule development because YABBYs are highly expressed in ovules. The expression pattern of YABBY genes showed that approximately half of the YABBYs were down-regulated under different stress treatments. Collectively, our results represent a comprehensive genome-wide study of the YABBY gene family, which should be helpful in further detailed studies on the gene function and evolution of YABBY genes in cotton.