De novo transcriptome assembly of pummelo and molecular marker development

EST-SSR Markers' Development Based on RNA-Sequencing and Their Application in Population Genetic Structure and Diversity Analysis of Eleusine indica in China

Goosegrass (Eleusine indica) is one of the worst agricultural weeds in China. Molecular markers were developed for genetic diversity and population structure analyses. In this study, we identified 8391 expressed sequence tag-simple sequence repeat (EST-SSR) markers from the de novo assembled unigenes of E. indica. Mononucleotides were the most abundant type of repeats (3591, 42.79%), followed by trinucleotides (3162, 37.68%). The most dominant mononucleotide and trinucleotide repeat motifs were A/T (3406, 40.59%) and AAT/ATT (103, 1.5%), respectively. Fourteen pairs of EST-SSR primers were verified and used to analyze the genetic diversity and population structure of 59 goosegrass populations. A total of 49 alleles were amplified, with the number of alleles (N_a) ranging from two to eleven per locus, and the effective number of alleles (N_e) ranged from 1.07 to 4.53. The average polymorphic information content (PIC) was 0.36. Genetic structure analysis (K = 2) and principal coordinate analysis divided 59 E. indica populations into two groups in a manner similar to the unweighted pair-group method (Dice genetic similarity coefficient = 0.700). This study developed a set of EST-SSR markers in E. indica and successfully analyzed the diversity and population genetic structures of 59 E. indica populations in China.

Comprehensive identification and characterization of simple sequence repeats based on the whole-genome sequences of 14 forest and fruit trees

Simple sequence repeats (SSRs) are popular and important molecular markers that exist widely in plants. Here, we conducted a comprehensive identification and comparative analysis of SSRs in 14 tree species. A total of 16, 298 SSRs were identified from 429, 449 genes, and primers were successfully designed for 99.44% of the identified SSRs. Our analysis indicated that tri-nucleotide SSRs were the most abundant, with an average of ~834 per species. Functional enrichment analysis by combining SSR-containing genes in all species, revealed 50 significantly enriched terms, with most belonging to transcription factor families associated with plant development and abiotic stresses such as Myeloblastosis_DNA-bind_4 (Myb_DNA-bind_4), APETALA2 (AP2), and Fantastic Four meristem regulator (FAF). Further functional enrichment analysis showed that 48 terms related to abiotic stress regulation and floral development were significantly enriched in ten species, whereas no significantly enriched terms were found in four species. Interestingly, the largest number of enriched terms was detected in Citrus sinensis (L.) Osbeck, accounting for 54.17% of all significantly enriched functional terms. Finally, we analyzed AP2 and trihelix gene families (Myb_DNA-bind_4) due to their significant enrichment in SSR-containing genes. The results indicated that whole-genome duplication (WGD) and whole genome triplication (WGT) might have played major roles in the expansion of the AP2 gene family but only slightly affected the expansion of the trihelix gene family during evolution. In conclusion, the identification and comprehensive characterization of SSR markers will greatly facilitate future comparative genomics and functional genomics studies.

Development of SSR markers via de novo transcriptome assembly in Akebia trifoliata (Thunb.) Koidz

Owing to its high nutritive, economic, and medicinal values, Akebia trifoliata has received increased attention, making worthy of being used as a new fruit crop for further domestication and commercialization in China. However, molecular research of A. trifoliata has lagged as investigations of its genomic resources and molecular markers are rare. In this study, a cDNA library of A. trifoliata leaves was sequenced using the Illumina NovaSeq. 6000 sequencing system. In total, 101 417 transcripts, 63 757 unigenes, and 9494 simple sequence repeats were assembled and identified from the transcriptome datasets. The majority of the SSRs were di- and trinucleotide repeats. Length and number of SSR motifs ranged from 15 to 66, and 5 to 48 bp, respectively. Of which, the A/T mononucleotide motif and AG/TC and CT/GA dinucleotide motifs were the most abundant. Furthermore, 100 SSR primers were randomly selected to validate amplification and polymorphism, and 88 A. trifoliata accessions were definitively distinguished by 49 primers. With the Qinling mountains and Huaihe River line as the boundaries, the northern and southern accessions were clustered into different groups, but no clear geographical patterns (city or origin) were observed in the southern accessions. These newly identified molecular markers may provide a foundation for the genetic identification and diversity analysis and marker-assisted selection breeding in species of Akebia.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

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.

Transcriptome analysis and identification of genes associated with flower development in Rhododendron pulchrum Sweet (Ericaceae)

Flowering process is essential for plant development. However, the molecular mechanisms driving flower development of ornamental woody Rhododendron pulchrum Sweet are difficult to elucidate due to the lack of genomic data. In this research, high-throughput sequencing and comparative transcriptome analyses of R. pulchrum flowers collected at three key stages were performed: floral bud stage, early flowering stage, and full-flowering stage. Furthermore, expression of genes involved in flower development was also validated with quantitative real-time PCR (qRT-PCR). RNA-seq yielded 96,350,697 bp of clean reads, which were assembled into 98,610 unigenes with an average length of 717 bp. 58,279 (59.10%) unigenes could be annotated, including 324 major unigenes associated with floral development. In addition, ten modules (20,443 mRNAs) were dissected in the co-expression network. Especially, Flowering Locus (FLC) and Flowering Locus T (FT) were co-expressed. 9493 differentially expressed genes (DEGs) were scanned among three stages, and most DEGs existed between flower bud stage and early flowering stage. In particular, 79 DGEs associated with flowering process were enriched in 28 GO terms. Moreover, the expression levels of MYC2, EIN3, and ARR-B were all lowest at early flowering stage, while transcripts of MYC2, TIR1, CYCD3, COL-1, and EIN3 were all peaked at flower bud stage. Transcriptome profile presented here will benefit deep insights into molecular mechanism underlying R. pulchrum flowering process.

Reproduction in woody perennial Citrus: an update on nucellar embryony and self-incompatibility

Review on citrus reproduction. Citrus is one of the most important and widely grown fruit crops. It possesses several special reproductive characteristics, such as nucellar embryony and self-incompatibility. The special phenomenon of nucellar embryony in citrus, also known as the polyembryony, is a kind of sporophytic apomixis. During the past decade, the emergence of novel technologies and the construction of multiple citrus reference genomes have facilitated rapid advances to our understanding of nucellar embryony. Indeed, several research teams have preliminarily determined the genetic basis of citrus apomixis. On the other hand, the phenomenon of self-incompatibility that promotes genetic diversity by rejecting self-pollen and accepting non-self-pollen is difficult to study in citrus because the long juvenile period of citrus presents challenges to identifying candidate genes that control this phenomenon. In this review, we focus on advances to our understanding of reproduction in citrus from the last decade and discuss priorities for the coming decade.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

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.

Genome-wide identification and functional analysis of S-RNase involved in the self-incompatibility of citrus

S-RNase-based self-incompatibility is found in Solanaceae, Rosaceae, and Scrophulariaceae, and is the most widespread mechanism that prevents self-fertilization in plants. Although 'Shatian' pummelo (Citrus grandis), a traditional cultivated variety, possesses the self-incompatible trait, the role of S-RNases in the self-incompatibility of 'Shatian' pummelo is poorly understood. To identify genes associated with self-incompatibility in citrus, we identified 16 genes encoding homologs of ribonucleases in the genomes of sweet orange (Citrus sinensis) and clementine mandarin (Citrus clementine). We preliminarily distinguished S-RNases from S-like RNases with a phylogenetic analysis that classified these homologs into three groups, which is consistent with the previous reports. Expression analysis provided evidence that CsRNS1 and CsRNS6 are S-like RNase genes. The expression level of CsRNS1 was increased during fruit development. The expression of CsRNS6 was increased during the formation of embryogenic callus. In contrast, we found that CsRNS3 possessed several common characteristics of the pistil determinant of self-incompatibility: it has an alkaline isoelectric point (pI), harbors only one intron, and is specifically expressed in style. We obtained a cDNA encoding CgRNS3 from 'Shatian' pummelo and found that it is homolog to CsRNS3 and that CgRNS3 exhibited the same expression pattern as CsRNS3. In an in vitro culture system, the CgRNS3 protein significantly inhibited the growth of self-pollen tubes from 'Shatian' pummelo, but after a heat treatment, this protein did not significantly inhibit the elongation of self- or non-self-pollen tubes. In conclusion, an S-RNase gene, CgRNS3, was obtained by searching the genomes of sweet orange and clementine for genes exhibiting sequence similarity to ribonucleases followed by expression analyses. Using this approach, we identified a protein that significantly inhibited the growth of self-pollen tubes, which is the defining property of an S-RNase.

Transcriptome Analysis of the Tadpole Shrimp (Triops longicaudatus) by Illumina Paired-End Sequencing: Assembly, Annotation, and Marker Discovery

The tadpole shrimp (Triops longicaudatus) is an aquatic crustacean that helps control pest populations. It inhabits freshwater ponds and pools and has been described as a living fossil. T. longicaudatus was officially declared an endangered species South Korea in 2005; however, through subsequent protection and conservation management, it was removed from the endangered species list in 2012. The limited number of available genetic resources on T. longicaudatus makes it difficult to obtain valuable genetic information for marker-aided selection programs. In this study, whole-transcriptome sequencing of T. longicaudatus generated 39.74 GB of clean data and a total of 269,822 contigs using the Illumina HiSeq 2500 platform. After clustering, a total of 208,813 unigenes with an N₅₀ length of 1089 bp were generated. A total of 95,105 unigenes were successfully annotated against Protostome (PANM), Unigene, Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using BLASTX with a cut-off of 1E−5. A total of 57,731 unigenes were assigned to GO terms, and 7247 unigenes were mapped to 129 KEGG pathways. Furthermore, 1595 simple sequence repeats (SSRs) were detected from the unigenes with 1387 potential SSR markers. This is the first report of high-throughput transcriptome analysis of T. longicaudatus, and it provides valuable insights for genetic research and molecular-assisted breeding of this important species.

De Novo Sequencing and Transcriptome Analysis of Pleurotus eryngii subsp. tuoliensis (Bailinggu) Mycelia in Response to Cold Stimulation

Cold stimulation of Bailinggu's mycelia is the main factor that triggers primordia initiation for successful production of fruiting bodies under commercial cultivation. Yet, the molecular-level mechanisms involved in mycelia response to cold stimulation are still unclear. Here, we performed comparative transcriptomic analysis using RNA-Seq technology to better understand the gene expression regulation during different temporal stages of cold stimulation in Bailinggu. A total of 21,558 Bailinggu mycelia unigenes were de novo assembled and annotated from four libraries (control at 25 °C, plus cold stimulation treatments at -3 °C for a duration of 1-2 days, 5-6 days, and 9-10 days). GO and KEGG pathway analysis indicated that functional groups of differentially expressed unigenes associated with cell wall and membrane stabilization, calcium signaling and mitogen-activated protein kinases (MAPK) pathways, and soluble sugars and protein biosynthesis and metabolism pathways play a vital role in Bailinggu's response to cold stimulation. Six hundred and seven potential EST-based SSRs loci were identified in these unigenes, and 100 EST-SSR primers were randomly selected for validation. The overall polymorphism rate was 92% by using 10 wild strains of Bailinggu. Therefore, these results can serve as a valuable resource for a better understanding of the molecular mechanisms associated with Bailinggu's response to cold stimulation.

De Novo Transcriptome Assembly in Firmiana danxiaensis, a Tree Species Endemic to the Danxia Landform

Many Firmiana species are locally endemic, providing an interesting system for studying adaptation and speciation. Among these species, F. danxiaensis is a tree species endemic to Mount Danxia in Guangdong, China, which is an area known for presenting the Danxia landform. How F. danxiaensis could have adapted to the stressful environment of rocky cliffs covered with barren soils in the Danxia landform is still unknown. In this study, we performed de novo assembly of the transcriptome of F. danxiaensis, obtaining 47,221 unigenes with an N50 value of 987 bp. Homology analysis showed that 32,318 of the unigenes presented hits in the NCBI non-redundant database, and 31,857 exhibited significant matches with the protein database of Theobroma cacao. Gene Ontology (GO) annotation showed that hundreds of unigenes participated in responses to various stresses or nutritional starvation, which may help us to understand the adaptation of F. danxiaensis to Danxia landform. Additionally, we found 263 genes related to responses to Cd, partially explaining the high accumulation of Cd observed in Firmiana species. The EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations revealed many genes playing roles in the biosynthesis of secondary metabolites and environmental adaptation, which may also contribute to the survivor and success of Firmiana species in extreme environments. Based on the obtained transcriptome, we further identified a Firmiana-specific whole-genome duplication event that occurred approximately 20 Mya, which may have provided raw materials for the diversification of Firmiana species.