Characterization of the global transcriptome for cotton (Gossypium hirsutum L.) anther and development of SSR marker

Transcriptome characterization and generation of marker resource for Himalayan vulnerable species, Ulmus wallichiana

Ulmus wallichiana is a traditional medicinal plant listed as a vulnerable in the IUCN red list data. Genomic and transcriptomic resources for this species are lacking, hindering its genetic exploration. Further, no polymorphic marker resource is available for this species, thus limiting the elucidation of its underlying genetic diversity, which is a pre-requisite for its conservation. This study was therefore aimed to generate a functionally annotated transcriptomic resource and screen it for SSR regions. We used paired-end Illumina based RNAseq technology and trinity based de novo assembly approach to generate full length transcripts, which were screened for SSR regions and functionally annotated. Around 6.6 million raw reads were de novo assembled transcripts, which were clustered into 146,083 unigenes. 19,909 transcripts were provided with 3986 unique KEGG ids, 70,519 transcripts with 6621 unique Pfam domains, and 45,125 transcripts with 7302 unique INTERPRO domains. 1456 transcripts were identified as transcriptions factors (TFs). Further, 8868 unique GO terms were obtained for the unigenes. The transcripts mapped to 23,056 known pre-determined orthology clusters in the eggNOG database. A total of 16,570 SSRs were identified from the unigenes. Out of the 90 SSRs selected for characterization on 20 genotypes, 28 were polymorphic. Mean effective alleles (Ne) of 2.53, mean observed heterozygosity (Ho) of 0.77, and average polymorphic information content (PIC) of 0.57 were found. This study may facilitate the genetic exploration of this species. The polymorphic SSRs would prove useful to explore its genetic diversity patterns, required for its conservation.

Screening of abiotic stress-responsive cotton genes using a cotton full-length cDNA overexpressing Arabidopsis library

Cotton (Gossypium hirsutum L.) is a major crop and the main source of natural fiber worldwide. Because various abiotic and biotic stresses strongly influence cotton fiber yield and quality, improved stress resistance of this crop plant is urgently needed. In this study, we used Gateway technology to construct a normalized full-length cDNA overexpressing (FOX) library from upland cotton cultivar ZM12 under various stress conditions. The library was transformed into Arabidopsis to produce a cotton-FOX-Arabidopsis library. Screening of this library yielded 6,830 transgenic Arabidopsis lines, of which 757 were selected for sequencing to ultimately obtain 659 cotton ESTs. GO and KEGG analyses mapped most of the cotton ESTs to plant biological process, cellular component, and molecular function categories. Next, 156 potential stress-responsive cotton genes were identified from the cotton-FOX-Arabidopsis library under drought, salt, ABA, and other stress conditions. Four stress-related genes identified from the library, designated as GhCAS, GhAPX, GhSDH, and GhPOD, were cloned from cotton complementary DNA, and their expression patterns under stress were analyzed. Phenotypic experiments indicated that overexpression of these cotton genes in Arabidopsis affected the response to abiotic stress. The method developed in this study lays a foundation for high-throughput cloning and rapid identification of cotton functional genes.

De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through Illumina sequencing

Curcuma alismatifolia widely used as an ornamental plant in Thailand and Cambodia. This species of herbaceous perennial from the Zingiberaceae family, includes cultivars with a wide range of colours and long postharvest life, and is used as an ornamental cut flower, as a potted plant, and in exterior landscapes. For further genetic improvement, however, little genomic information and no specific molecular markers are available. The present study used Illumina sequencing and de novo transcriptome assembly of two C. alismatifolia cvs, 'Chiang Mai Pink' and 'UB Snow 701', to develop simple sequence repeat markers for genetic diversity studies. After de novo assembly, 62,105 unigenes were generated and 48,813 (78.60%) showed significant similarities versus six functional protein databases. In addition, 9,351 expressed sequence tag-simple sequence repeats (EST-SSRs) were identified with a distribution frequency of 12.5% total unigenes. Out of 8,955 designed EST-SSR primers, 150 primers were selected for the development of potential molecular markers. Among these markers, 17 EST-SSR markers presented a moderate level of genetic diversity among three C. alismatifolia cultivars, one hybrid, three Curcuma, and two Zingiber species. Three different genetic groups within these species were revealed using EST-SSR markers, indicating that the markers developed in this study can be effectively applied to the population genetic analysis of Curcuma and Zingiber species. This report describes the first analysis of transcriptome data of an important ornamental ginger cultivars, also provides a valuable resource for gene discovery and marker development in the genus Curcuma.

Insights into the cotton anther development through association analysis of transcriptomic and small RNA sequencing

BACKGROUND

Plant anther development is a systematic and complex process precisely controlled by genes. Regulation genes and their regulatory mechanisms for this process remain elusive. In contrast to numerous researches on anther development with respect to mRNAs or miRNAs in many crops, the association analysis combining both omics has not been reported on cotton anther.

RESULTS

In this study, the molecular mechanism of cotton anther development was investigated with the employment of association analysis of transcriptome and small RNA sequencing during the predefined four stages of cotton anther development, sporogenuous cell proliferation (SCP), meiotic phase (MP), microspore release period (MRP) and pollen maturity (PM). Analysis revealed that the differentially expressed genes are increasingly recruited along with the developmental progress. Expression of functional genes differed significantly among developmental stages. The genes related with cell cycle, progesterone-mediated oocyte maturation, and meiosis are predominantly expressed at the early stage of anther development (SCP and MP), and the expression of genes involved in energy metabolism, flavonoid biosynthesis, axon guidance and phospholipase D signaling pathways is mainly enriched at the late stage of anther development (MRP and PM). Analysis of expression patterns revealed that there was the largest number of differentially expressed genes in the MP and the expression profiles of differentially expressed genes were significantly increased, which implied the importance of MP in the entire anther development cycle. In addition, prediction and analysis of miRNA targeted genes suggested that miRNAs play important roles in anther development. The miRNAs ghr-miR393, Dt_chr12_6065 and At_chr9_3080 participated in cell cycle, carbohydrate metabolism and auxin anabolism through the target genes, respectively, to achieve the regulation of anther development.

CONCLUSIONS

Through the association analysis of mRNA and miRNA, our work gives a better understanding of the preferentially expressed genes and regulation in different developmental stages of cotton anther and the importance of meiotic phase, and also the involvement of miRNAs in precise regulation for this process, which would be valuable for clarifying the mechanism of plant anther development in response to internal and external environments.

Transcriptome profiling and cataloging differential gene expression in floral buds of fertile and sterile lines of cotton (Gossypium hirsutum L.)

Cytoplasmic Male Sterility is maternally inherited trait in plants, characterized by failure to produce functional pollen during anther development. Anther development is modulated through the interaction of nuclear and mitochondrial genes. In the present study, differential gene expression of floral buds at the sporogenous stage (SS) and microsporocyte stage (MS) between CGMS and its fertile maintainer line of cotton plants was studied. A total of 320 significantly differentially expressed genes, including 20 down-regulated and 37 up-regulated in CGMS comparing with its maintainer line at the SS stage, as well as and 89 down-regulated and 4 up-regulated in CGMS compared to the fertile line at MS stage. Comparing the two stages in the same line, there were 6 down-regulated differentially expressed genes only induced in CGMS and 9 up-regulated differentially expressed gene only induced in its maintainer. GO analysis revealed essential genes responsible for pollen development, and cytoskeleton category show differential expression between the fertile and CGMS lines. Validation studies by qRT-PCR shows concordance with RNA-seq result. A set of novel SSRs identified in this study can be used in evaluating genetic relationships among cultivars, QTL mapping, and marker-assisted breeding. We reported aberrant expression of genes related to pollen exine formation, and synthesis of pectin lyase, myosine heavy chain, tubulin, actin-beta, heat shock protein and myeloblastosis (MYB) protein as targets for CMS in cotton. The results of this study contribute to basic information for future screening of genes and identification of molecular portraits responsible for CMS as well as to elucidate molecular mechanisms that lead to CMS in cotton.

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.

De novo transcriptome assembly of Zanthoxylum bungeanum using Illumina sequencing for evolutionary analysis and simple sequence repeat marker development

Zanthoxylum, an ancient economic crop in Asia, has a satisfying aromatic taste and immense medicinal values. A lack of genomic information and genetic markers has limited the evolutionary analysis and genetic improvement of Zanthoxylum species and their close relatives. To better understand the evolution, domestication, and divergence of Zanthoxylum, we present a de novo transcriptome analysis of an elite cultivar of Z. bungeanum using Illumina sequencing; we then developed simple sequence repeat markers for identification of Zanthoxylum. In total, we predicted 45,057 unigenes and 22,212 protein coding sequences, approximately 90% of which showed significant similarities to known proteins in databases. Phylogenetic analysis indicated that Zanthoxylum is relatively recent and estimated to have diverged from Citrus ca. 36.5–37.7 million years ago. We also detected a whole-genome duplication event in Zanthoxylum that occurred 14 million years ago. We found no protein coding sequences that were significantly under positive selection by Ka/Ks. Simple sequence repeat analysis divided 31 Zanthoxylum cultivars and landraces into three major groups. This Zanthoxylum reference transcriptome provides crucial information for the evolutionary study of the Zanthoxylum genus and the Rutaceae family, and facilitates the establishment of more effective Zanthoxylum breeding programs.

Development of EST-SSR markers in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee) based on de novo transcriptomic assemblies

Flowering Chinese cabbage is one of the most important vegetable crops in southern China. Genetic improvement of various agronomic traits in this crop is underway to meet high market demand in the region, but the progress is hampered by limited number of molecular markers available in this crop. This study aimed to develop EST-SSR markers from transcriptome sequences generated by next-generation sequencing. RNA-seq of eight cabbage samples identified 48,975 unigenes. Of these unigenes, 23,267 were annotated in 56 gene ontology (GO) categories, 6,033 were mapped to 131 KEGG pathways, and 7,825 were assigned to clusters of orthologous groups (COGs). From the unigenes, 8,165 EST-SSR loci were identified and 98.57% of them were 1-3 nucleotide repeats with 14.32%, 41.08% and 43.17% of mono-, di- and tri-nucleotide repeats, respectively. Fifty-eight types of motifs were identified with A/T, AG/CT, AT/AT, AC/GT, AAG/CTT and AGG/CCT the most abundant. The lengths of repeated nucleotide sequences in all SSR loci ranged from 12 to 60 bp, with most (88.51%) under 20 bp. Among 170 primer pairs were randomly selected from a total of 4,912 SSR primers we designed, 48 yielded unambiguously polymorphic bands with high reproducibility. Cluster analysis using 48 SSRs classified 34 flowering Chinese cabbage cultivars into three groups. A large number of EST-SSR markers identified in this study will facilitate marker-assisted selection in the breeding programs of flowering Chinese cabbage.

A newly isolated and identified vitamin B12 producing strain: Sinorhizobium meliloti 320

Vitamin B12 (Cobalamin, VB12) has several physiological functions and is widely used in pharmaceutical and food industries. A new unicellular species was extracted from China farmland, and the strain could produce VB12 which was identified by HPLC and HPLC-MS/MS. 16S rDNA analysis reveals this strain belongs to the species Sinorhizobium meliloti and we named it S. meliloti 320. Its whole genome information indicates that this strain has a complete VB12 synthetic pathway, which paves the way for further metabolic engineering studies. The optimal carbon and nitrogen sources are sucrose and corn steep liquor (CSL) plus peptone. The optimal combination of sucrose and CSL was obtained by response surface methodology as they are the most suitable carbon and nitrogen sources, respectively. This strain could produce 140 ± 4.2 mg L(-1) vitamin B12 after incubating for 7 days in the optimal medium.

Large-scale developing of simple sequence repeat markers and probing its correlation with ramie (Boehmeria nivea L.) fiber quality

Marker-assisted selection is an important component of the discipline of molecular breeding. Using DNA markers to assist in plant breeding, the efficiency and precision could be greatly increased. However, the scarcity number of identified DNA markers has hindered the research and the breeding process of ramie (Boehmeria nivea L.) in many aspects, especially fiber quality, one of the top-priority breeding objectives of ramie. In this study, 4230 SSR loci were identified in 3969 unigenes (6.80 % of 58,369), which were de novo assembled from the transcriptome involving different ramie fiber developmental stages. Among these SSRs, the dinucleotides (1599, 37.80 %) and trinucleotides (772, 18.25 %) were most abundant; the motifs AG/CT (1140, 26.94 %), AT/AT (407, 9.62 %) and AGA/TCT (246, 8.31 %) comprised the three most abundant repeats. A total of 2431 primer pairs were designed flanking the SSRs and 1050 of them were employed in PCR amplification for their usefulness using three ramie cultivars. The results showed that 88.10 % of these primers could generate positive PCR bands in any of the three cultivars. Further phylogenetic analysis that conducted from the PCR amplification of 52 specifically sifted SSR primers within 17 cultivars approved that the possible correlation may exist between the primers and ramie fiber quality. These developed SSR markers could be applied in downstream studies, like genetic and physical maps, quantitative trait loci mapping, genetic diversity studies and cultivar fingerprinting, and breeding processes of ramie with better fiber quality under further confirmation of the correlation with ramie fiber quality.

Genome-wide mining, characterization, and development of microsatellite markers in gossypium species

Although much research has been conducted to characterize microsatellites and develop markers, the distribution of microsatellites remains ambiguous and the use of microsatellite markers in genomic studies and marker-assisted selection is limited. To identify microsatellites for cotton research, we mined 100,290, 83,160, and 56,937 microsatellites with frequencies of 41.2, 49.1, and 74.8 microsatellites per Mb in the recently sequenced Gossypium species: G. hirsutum, G. arboreum, and G. raimondii, respectively. The distributions of microsatellites in their genomes were non-random and were positively and negatively correlated with genes and transposable elements, respectively. Of the 77,996 developed microsatellite markers, 65,498 were physically anchored to the 26 chromosomes of G. hirsutum with an average marker density of 34 markers per Mb. We confirmed 67,880 (87%) universal and 7,705 (9.9%) new genic microsatellite markers. The polymorphism was estimated in above three species by in silico PCR and validated with 505 markers in G. hirsutum. We further predicted 8,825 polymorphic microsatellite markers within G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124. In our study, genome-wide mining and characterization of microsatellites, and marker development were very useful for the saturation of the allotetraploid genetic linkage map, genome evolution studies and comparative genome mapping.