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

Characterisation of the Cinnamomumparthenoxylon (Jack) Meisn (Lauraceae) transcriptome using Illumina paired-end sequencing and EST-SSR markers development for population genetics

Cinnamomumparthenoxylon is an endemic and endangered species with significant economic and ecological value in Vietnam. A better understanding of the genetic architecture of the species will be useful when planning management and conservation. We aimed to characterize the transcriptome of C.parthenoxylon, develop novel molecular markers, and assess the genetic variability of the species. First, transcriptome sequencing of five trees (C.parthenoxylon) based on root, leaf, and stem tissues was performed for functional annotation analysis and development of novel molecular markers. The transcriptomes of C.parthenoxylon were analyzed via an Illumina HiSeqTM 4000 sequencing system. A total of 27,363,199 bases were generated for C.parthenoxylon. De novo assembly indicated that a total of 160,435 unigenes were generated (average length = 548.954 bp). The 51,691 unigenes were compared against different databases, i.e. COG, GO, KEGG, KOG, Pfam, Swiss-Prot, and NR for functional annotation. Furthermore, a total of 12,849 EST-SSRs were identified. Of the 134 primer pairs, 54 were randomly selected for testing, with 15 successfully amplified across nine populations of C.parthenoxylon. We uncovered medium levels of genetic diversity (PIC = 0.52, Na = 3.29, Ne = 2.18, P = 94.07%, Ho = 0.56 and He = 0.47) within the studied populations. The molecular variance was 10% among populations and low genetic differentiation (Fst = 0.06) indicated low gene flow (Nm = 2.16). A reduction in the population size of C.parthenoxylon was detected using BOTTLENECK (VP population). The structure analysis suggested two optimal genetic clusters related to gene flow among the populations. Analysis of molecular variance (AMOVA) revealed higher genetic variation within populations (90%) than among populations (10%). The UPGMA approach and DAPC divided the nine populations into three main clusters. Our findings revealed a significant fraction of the transcriptome sequences and these newlydeveloped novel EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity and molecular marker-assisted selection in C.parthenoxylon. This study provides comprehensive genetic resources for the breeding and conservation of different varieties of C.parthenoxylon.

Development of simple sequence repeat markers for sugarcane from data mining of expressed sequence tags

Sugarcane (Saccharum spp. hybrids) is a worldwide acclaimed important agricultural crop used primarily for sugar production and biofuel. Sugarcane's genetic complexity, aneuploidy, and extreme heterozygosity make it a challenging crop in developing improved varieties. The molecular breeding programs promise to develop nutritionally improved varieties for both direct consumption and commercial application. Therefore, to address these challenges, the development of simple sequence repeats (SSRs) has been proven to be a powerful molecular tool in sugarcane. This study involved the collection of 285216 expressed sequence tags (ESTs) from sugarcane, resulting in 23666 unigenes, including 4547 contigs. Our analysis identified 4120 unigenes containing a total of 4960 SSRs, with the most abundant repeat types being monomeric (44.33%), dimeric (13.10%), and trimeric (39.68%). We further chose 173 primers to analyze the banding pattern in 10 sugarcane accessions by PAGE analysis. Additionally, functional annotation analysis showed that 71.07%, 53.6%, and 10.3% unigenes were annotated by Uniport, GO, and KEGG, respectively. GO annotations and KEGG pathways were distributed across three functional categories: molecular (46.46%), cellular (33.94%), and biological pathways (19.6%). The cluster analysis indicated the formation of four distinct clusters among selected sugarcane accessions, with maximum genetic distance observed among the varieties. We believe that these EST-SSR markers will serve as valuable references for future genetic characterization, species identification, and breeding efforts in sugarcane.

SSR marker based analysis for identification and of genetic diversity of non-heading Chinese cabbage varieties

As a widely cultivated vegetable in China and Southeast Asia, the breeding of non-heading Chinese cabbage (Brassica campestris ssp. chinensis Makino) is widespread; more than 400 varieties have been granted new plant variety rights (PVRs) in China. Distinctness is one of the key requirements for the granting of PVRs, and molecular markers are widely used as a robust supplementary method for similar variety selection in the distinctness test. Although many genome-wide molecular markers have been developed, they have not all been well used in variety identification and tests of distinctness of non-heading Chinese cabbage. In this study, by using 423 non-heading Chinese cabbage varieties collected from different regions of China, 287 simple sequence repeat (SSR) markers were screened for polymorphisms, and 23 core markers were finally selected. The polymorphic information content (PIC) values of the 23 SSR markers ranged from 0.555 to 0.911, with an average of 0.693, and the average number of alleles per marker was 13.65. Using these 23 SSR markers, 418 out of 423 varieties could be distinguished, with a discrimination rate of 99.994%. Field tests indicated that those undistinguished varieties were very similar and could be further distinguished by a few morphological characteristics. According to the clustering results, the 423 varieties could be divided into three groups: pak-choi, caitai, and tacai. The similarity coefficient between the SSR markers and morphological characteristics was moderate (0.53), and the efficiency of variety identification was significantly improved by using a combination of SSR markers and morphological characteristics.

Comprehensive transcriptome analysis reveals heat-responsive genes in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis) using RNA sequencing

Flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee, 2n=20, AA) is a vegetable species in southern parts of China that faces high temperatures in the summer and winter seasons. While heat stress adversely impacts plant productivity and survival, the underlying molecular and biochemical causes are poorly understood. This study investigated the gene expression profiles of heat-sensitive (HS) '3T-6' and heat-tolerant (HT) 'Youlu-501' varieties of flowering Chinese cabbage in response to heat stress using RNA sequencing. Among the 37,958 genes expressed in leaves, 20,680 were differentially expressed genes (DEGs) at 1, 6, and 12 h, with 1,078 simultaneously expressed at all time points in both varieties. Hierarchical clustering analysis identified three clusters comprising 1,958, 556, and 591 down-regulated, up-regulated, and up- and/or down-regulated DEGs (3205 DEGs; 8.44%), which were significantly enriched in MAPK signaling, plant-pathogen interactions, plant hormone signal transduction, and brassinosteroid biosynthesis pathways and involved in stimulus, stress, growth, reproductive, and defense responses. Transcription factors, including MYB (12), NAC (13), WRKY (11), ERF (31), HSF (17), bHLH (16), and regulatory proteins such as PAL, CYP450, and photosystem II, played an essential role as effectors of homeostasis, kinases/phosphatases, and photosynthesis. Among 3205 DEGs, many previously reported genes underlying heat stress were also identified, e.g., BraWRKY25, BraHSP70, BraHSPB27, BraCYP71A23, BraPYL9, and BraA05g032350.3C. The genome-wide comparison of HS and HT provides a solid foundation for understanding the molecular mechanisms of heat tolerance in flowering Chinese cabbage.

De novo transcriptome assembly using Illumina sequencing and development of EST-SSR markers in a monoecious herb Sagittaria trifolia Linn

Background

Sagittaria trifolia Linn. is a widespread macrophyte in Asia and southeast Europe and cultivated in parts of Asia. Although a few genomic studies have been conducted for S. trifolia var. sinensis, a crop breed, there is limited genomic information on the wild species of S. trifolia. Effective microsatellite markers are also lacking.

Objective

To assemble transcriptome sequence and develop effective EST-SSR markers for S. trifolia.

Methods

Here we developed microsatellite markers based on tri-, tetra-, penta-, and hexa-nucleotide repeat sequences by comparatively screening multiple transcriptome sequences of eleven individuals from ten natural populations of S. trifolia.

Results

A total of 107,022 unigenes were de novo assembled, with a mean length of 730 bp and an N50 length of 1,378 bp. The main repeat types were mononucleotide, trinucleotide, and dinucleotide, accounting for 55.83%, 23.51%, and 17.56% of the total repeats, respectively. A total of 86 microsatellite loci were identified with repeats of tri-, tetra-, penta-, and hexa-nucleotide. For SSR verification, 28 polymorphic loci from 41 randomly picked markers were found to produce stable and polymorphic bands, with the number of alleles per locus ranging from 2 to 11 and a mean of 5.2. The range of polymorphic information content (PIC) of each SSR locus varied from 0.25 to 0.80, with an average of 0.58. The expected heterozygosity ranged from 0.29 to 0.82, whereas the observed heterozygosity ranged from 0.25 to 0.90.

Conclusion

The assembled transcriptome and annotated unigenes of S. trifolia provide a basis for future studies on gene functions, pathways, and molecular mechanisms associated with this species and other related. The newly developed EST-SSR markers could be effective in examining population genetic structure, differentiation, and parentage analyses in ecological and evolutionary studies of S. trifolia.

Genetic Diversity and DNA Fingerprinting in Broccoli Carrying Multiple Clubroot Resistance Genes Based on SSR Markers

To identify cultivars quickly and accurately, DNA fingerprinting of 10 broccoli varieties was performed by using simple sequence repeat (SSR) marker technology. Highly informative and polymorphic SSR markers were screened using broccoli and rapeseed. Out of the 93 SSR marker pairs, 21 pairs were selected and found to have good polymorphism. Each marker pair generated 1 to 10 polymorphic bands with an average of 4.29. The average polymorphism information content (PIC) was 0.41 with a range from 0.16 to 0.95. Six selected marker pairs established the fingerprinting of the 10 accessions and their unique fingerprints. Cluster analysis of 10 accessions showed that the genetic similarity coefficient was between 0.57 and 0.91. They can be divided into 3 groups at the genetic similarity coefficient (GSC) of 0.73. The above results indicated that DNA fingerprinting could provide a scientific basis for the identification of broccoli polymerized multiple clubroot resistance genes. Research shows that SSR marker-based DNA fingerprinting further ensures plant seed purity.

Genome-wide identification and development of InDel markers in tobacco (Nicotiana tabacum L.) using RAD-seq

Insertions and deletions (InDels) can be used as molecular markers in genetic studies and marker-assisted selection breeding. However, genetic improvement in tobacco has been hindered by limited genetic diversity information and relatedness within available germplasm. A Chinese tobacco variety, Yueyan-98, was resequenced using restriction-site associated DNA (RAD-seq) approach to develop InDel markers. In total, 32,884 InDel loci were detected between Yueyan-98 and the K326 reference sequence [18,598 (56.55%) deletions and 14,288 (43.45%) insertions], ranging from 1 to 62 bp in length. Of the 6,733 InDels (> 4 bp) that were suitable for polyacrylamide gel electrophoresis, 150 were randomly selected. These 150 InDels were unevenly distributed on 23 chromosomes, and the highest numbers of InDels were observed on chromosomes Nt05, Nt13, and Nt23. The average density of adjacent InDels was 19.36 Mb. Thirty-seven InDels were located in genic regions. Polymerase chain reaction (PCR)-based markers were developed to validate polymorphism; 113 (79.80%) of the 150 InDel markers showed polymorphism and were further used for genetic diversity analysis of 50 tobacco accessions (13 from China, 1 from Mexico, and 36 from the USA). The average expected heterozygosity (He) and polymorphism information content (PIC) values were 0.28 ± 0.16 and 0.38 ± 0.10, respectively. The average Shannon diversity index (I) was 0.34 ± 0.18, with genetic diversity ranging from 0.13-0.57. The 50 accessions were classified into two groups with a genetic similarity coefficient of 0.68. Principal coordinate analysis (PCoA) and population structure analysis showed similar results and divided the population into two groups unrelated to their geographical origins. AMOVA showed 4% variance among the population and the remaining 96% within the population, suggesting low genetic differentiation between two subpopulations. Furthermore, 10 InDels (19 alleles) were significantly identified for tobacco plant height using GLM+Q model at P < 0.005. Among these, three markers (Nt-I-26, Nt-I-41, and Nt-I-44) were detected in at least two environments, with phenotypic variance explained (PVE) ranging from 14.03 to 32.68%. The polymorphic InDel markers developed can be used for hybrid identification, genetic diversity, genetic linkage map construction, gene mapping, and MAS breeding programs of tobacco.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01187-3.

Identification and Analysis of Small Interfering RNAs Associated With Heat Stress in Flowering Chinese Cabbage Using High-Throughput Sequencing

Endogenous small interfering RNAs (siRNAs) are substantial gene regulators in eukaryotes and play key functions in plant development and stress tolerance. Among environmental factors, heat is serious abiotic stress that severely influences the productivity and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, how siRNAs are involved in regulating gene expression during heat stress is not fully understood in flowering Chinese cabbage. Combining bioinformatical and next-generation sequencing approaches, we identified heat-responsive siRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0, 1, 6, and 12 h after a 38°C heat-stress treatment; 536, 816, and 829 siRNAs exhibited substantial differential expression at 1, 6, and 12 h, respectively. Seventy-five upregulated and 69 downregulated differentially expressed siRNAs (DE-siRNAs) were common for the three time points of heat stress. We identified 795 target genes of DE-siRNAs, including serine/threonine-protein kinase SRK2I, CTR1-like, disease resistance protein RML1A-like, and RPP1, which may play a role in regulating heat tolerance. Gene ontology showed that predictive targets of DE-siRNAs may have key roles in the positive regulation of biological processes, organismal processes, responses to temperature stimulus, signaling, and growth and development. These novel results contribute to further understanding how siRNAs modulate the expression of their target genes to control heat tolerance in flowering Chinese cabbage.

Illuminating the Plant Rhabdovirus Landscape through Metatranscriptomics Data

Rhabdoviruses infect a large number of plant species and cause significant crop diseases. They have a negative-sense, single-stranded unsegmented or bisegmented RNA genome. The number of plant-associated rhabdovirid sequences has grown in the last few years in concert with the extensive use of high-throughput sequencing platforms. Here, we report the discovery of 27 novel rhabdovirus genomes associated with 25 different host plant species and one insect, which were hidden in public databases. These viral sequences were identified through homology searches in more than 3000 plant and insect transcriptomes from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) using known plant rhabdovirus sequences as the query. The identification, assembly and curation of raw SRA reads resulted in sixteen viral genome sequences with full-length coding regions and ten partial genomes. Highlights of the obtained sequences include viruses with unique and novel genome organizations among known plant rhabdoviruses. Phylogenetic analysis showed that thirteen of the novel viruses were related to cytorhabdoviruses, one to alphanucleorhabdoviruses, five to betanucleorhabdoviruses, one to dichorhaviruses and seven to varicosaviruses. These findings resulted in the most complete phylogeny of plant rhabdoviruses to date and shed new light on the phylogenetic relationships and evolutionary landscape of this group of plant viruses. Furthermore, this study provided additional evidence for the complexity and diversity of plant rhabdovirus genomes and demonstrated that analyzing SRA public data provides an invaluable tool to accelerate virus discovery, gain evolutionary insights and refine virus taxonomy.

Functional Annotation of a Full-Length Transcriptome and Identification of Genes Associated with Flower Development in Rhododendronsimsii (Ericaceae)

Rhododendronsimsii is one of the top ten famous flowers in China. Due to its historical value and high aesthetic, it is widely popular among Chinese people. Various colors are important breeding objectives in Rhododendron L. The understanding of the molecular mechanism of flower color formation can provide a theoretical basis for the improvement of flower color in Rhododendron L. To generate the R.simsii transcriptome, PacBio sequencing technology has been used. A total of 833,137 full-length non-chimeric reads were obtained and 726,846 high-quality full-length transcripts were found. Moreover, 40,556 total open reading frames were obtained; of which 36,018 were complete. In gene annotation analyses, 39,411, 18,565, 16,102 and 17,450 transcriptions were allocated to GO, Nr, KEGG and COG databases, correspondingly. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with Protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT) and Coding Non Coding Index (CNCI) databases and observed 6170, 2265, 4084 and 1240 lncRNAs, respectively. Based on the results, most genes were enriched in the flavonoid biosynthetic pathway. The eight key genes on the anthocyanin biosynthetic pathway were further selected and analyzed by qRT-PCR. The F3'H and ANS showed an upward trend in the developmental stages of R. simsii. The highest expression of F3'5'H and FLS in the petal color formation of R. simsii was observed. This research provided a huge number of full-length transcripts, which will help to proceed genetic analyses of R.simsii. native, which is a semi-deciduous shrub.

Developing a new genic SSR primer database in faba bean (Vicia faba L.)

The goal of this research was to develop a new genetic database of simple sequence repetition (SSR) primers for faba and classify them according to their target genes and respective biological processes. Approximately 75,605 and 148,196 previously published genomic and transcriptomic faba sequences, respectively, have been used to detect possible SSRs. The numbers of identified SSRs from each dataset were 25,502 and 12,319, respectively. The distribution of different repeat classes indicated that trinucleotides represent the largest number of repeat counts, followed by dinucleotides. The extracted genic SSR sequences were used to design 1091 polymerase chain reaction (PCR) primers, of which only 238 (21.8%) primers target genomic sequences and the other 853 PCR primers targeted transcriptomic sequences. The annotation of gene-targeted SSRs showed that approximately 897 genes were targeted by our SSR primers. Approximately 1890 gene ontology (GO) identification codes have been obtained. The GO keywords were distributed among distinct molecular cell features. The highest redundancies involved 554 technical words, 196 domains, and 160 molecular feature phrases. These GO codes belonged to the general level of GO and included molecular function, cellular component, and biological process (544, 670, and 676 GOs, respectively). Twenty-seven SSR PCR primers were synthesized to 12 Egyptian faba bean genotypes. Approximately 11 SSR provided one to two PCR bands, whereas other SSRs provided only one sharp band with polymorphic band size. There were 13 polymorphic primers. The polymorphism information content was 0.3, which implied moderate informativeness.

Transcriptome Sequencing and Development of Novel Genic SSR Markers From Pistacia vera L

In this study, we aimed to develop novel genic simple sequence repeat (eSSR) markers and to study phylogenetic relationship among Pistacia species. Transcriptome sequencing was performed in different tissues of Siirt and Atl cultivars of pistachio (Pistacia vera). A total of 37.5-Gb data were used in the assembly. The number of total contigs and unigenes was calculated as 98,831, and the length of N50 was 1,333 bp after assembly. A total of 14,308 dinucleotide, trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide SSR motifs (4–17) were detected, and the most abundant SSR repeat types were trinucleotide (29.54%), dinucleotide (24.06%), hexanucleotide (20.67%), pentanucleotide (18.88%), and tetranucleotide (6.85%), respectively. Overall 250 primer pairs were designed randomly and tested in eight Pistacia species for amplification. Of them, 233 were generated polymerase chain reaction products in at least one of the Pistacia species. A total of 55 primer pairs that had amplifications in all tested Pistacia species were used to characterize 11 P. vera cultivars and 78 wild Pistacia genotypes belonging to nine Pistacia species (P. khinjuk, P. eurycarpa, P. atlantica, P. mutica, P. integerrima, P. chinensis, P. terebinthus, P. palaestina, and P. lentiscus). A total of 434 alleles were generated from 55 polymorphic eSSR loci with an average of 7.89 alleles per locus. The mean number of effective allele was 3.40 per locus. Polymorphism information content was 0.61, whereas observed (Ho) and expected heterozygosity (He) values were 0.39 and 0.65, respectively. UPGMA (unweighted pair-group method with arithmetic averages) and STRUCTURE analysis divided 89 Pistacia genotypes into seven populations. The closest species to P. vera was P. khinjuk. P. eurycarpa was closer P. atlantica than P. khinjuk. P. atlantica–P. mutica and P. terebinthus–P. palaestina pairs of species were not clearly separated from each other, and they were suggested as the same species. The present study demonstrated that eSSR markers can be used in the characterization and phylogenetic analysis of Pistacia species and cultivars, as well as genetic linkage mapping and QTL (quantitative trait locus) analysis.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L

Rhododendron lapponicum L. is a familiar ornamental plant worldwide with important ornamental and economic value. However, a full-length R. lapponicum transcriptome is still lacking. In the present study, we used the Pacific Biosciences single-molecule real-time sequencing technology to generate the R. lapponicum transcriptome. A total of 346,270 full-length non-chimeric reads were generated, from which we obtained 75,002 high-quality full-length transcripts. We identified 55,255 complete open reading frames, 7,140 alternative splicing events and 2,011 long non-coding RNAs. In gene annotation analyses, 71,155, 33,653, 30,359 and 31,749 transcripts were assigned to the Nr, GO, COG and KEGG databases, respectively. Additionally, 3,150 transcription factors were detected. KEGG pathway analysis showed that 96 transcripts were identified coding for the enzymes associated with anthocyanin synthesis. Furthermore, we identified 64,327 simple sequence repeats from 45,319 sequences, and 150 pairs of primers were randomly selected to develop SSR markers. This study provides a large number of full-length transcripts, which will facilitate the further study of the genetics of R. lapponicum.

Comparative Analysis of miRNA Expression Profiles Between Heat-Tolerant and Heat-Sensitive Genotypes of Flowering Chinese Cabbage Under Heat Stress Using High-Throughput Sequencing

Heat stress disturbs cellular homeostasis, thus usually impairs yield of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). MicroRNAs (miRNAs) play a significant role in plant responses to different stresses by modulating gene expression at the post-transcriptional level. However, the roles that miRNAs and their target genes may play in heat tolerance of flowering Chinese cabbage remain poorly characterized. The current study sequenced six small RNA libraries generated from leaf tissues of flowering Chinese cabbage collected at 0, 6, and 12 h after 38 °C heat treatment, and identified 49 putative novel miRNAs and 43 known miRNAs that differentially expressed between heat-tolerant and heat-sensitive flowering Chinese cabbage. Among them, 14 novel and nine known miRNAs differentially expressed only in the heat-tolerant genotype under heat-stress, therefore, their target genes including disease resistance protein TAO1-like, RPS6, reticuline oxidase-like protein, etc. might play important roles in enhancing heat-tolerance. Gene Ontology (GO) analysis revealed that targets of these differentially expressed miRNAs may play key roles in responses to temperature stimulus, cell part, cellular process, cell, membrane, biological regulation, binding, and catalytic activities. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified their important functions in signal transduction, environmental adaptation, global and overview maps, as well as in stress adaptation and in MAPK signaling pathways such as cell death. These findings provide insight into the functions of the miRNAs in heat stress tolerance of flowering Chinese cabbage.

Development of SSR markers in Paeonia based on De Novo transcriptomic assemblies

Peony is a famous ornamental and medicinal plant in China, and peony hybrid breeding is an important means of germplasm innovation. However, research on the genome of this species is limited, thereby hindering the genetic and breeding research on peony. In the present study, simple sequence repeat (SSR) locus analysis was performed on expressed sequence tags obtained by the transcriptome sequencing of Paeonia using Microsatellite software. Primers with polymorphism were obtained via polymerase chain reaction amplification and electrophoresis. As a result, a total of 86,195 unigenes were obtained by assembling the transcriptome data of Paeonia. Functional annotations were obtained in seven functional databases including 49,172 (Non-Redundant Protein Sequence Database: 57.05%), 38,352 (Nucleotide Sequence Database: 44.49%), 36,477 (Swiss Prot: 42.32%), 38,905 (Clusters of Orthologous Groups for Eukaryotic Complete Genomes: 45.14%), 37,993 (Kyoto Encyclopedia of Genes and Genomes: 44.08%), 26,832 (Gene Ontology: 31.13%) and 37,758 (Pfam: 43.81%) unigenes. Meanwhile, 21,998 SSR loci were distributed in 17,567 unigenes containing SSR sequences, and the SSR distribution frequency was 25.52%, with an average of one SSR sequence per 4.66 kb. Mononucleotide, dinucleotide, and trinucleotide were the main repeat types, accounting for 55.74%, 25.58%, and 13.21% of the total repeat times, respectively. Forty-five pairs of the 100 pairs of primers selected randomly could amplify clear polymorphic bands. The polymorphic primers of these 45 pairs were used to cluster and analyze 16 species of peony. The new SSR molecular markers can be useful for the study of genetic diversity and marker-assisted breeding of peony.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

Characterization and expression profiles of miRNAs in the triploid hybrids of Brassica napus and Brassica rapa

BACKGROUND

Polyploidy provides a means of interspecific genome transfer to incorporate preferable traits from progenitor to progeny. However, few studies on miRNA expression profiles of interspecific hybrids of B. napus (AnAnCnCn) and B. rapa (ArAr) have been reported.

RESULTS

Here, we apply small RNA sequencing to explore miRNA expression patterns between B. napus, B. rapa and their F1 hybrid. Bioinformatics analysis identified 376, 378, 383 conserved miRNAs and 82, 76, 82 novel miRNAs in B. napus, B. rapa and the F1 hybrid, respectively. Moreover, 213 miRNAs were found to be differentially expressed between B. napus, B. rapa and the F1 hybrid. The present study also shows 211 miRNAs, including 77 upregulated and 134 downregulated miRNAs, to be nonadditively expressed in the F1 hybrid. Furthermore, miRNA synteny analysis revealed high genomic conservation between the genomes of B. napus, B. rapa and their F1 hybrid, with some miRNA loss and gain events in the F1 hybrid.

CONCLUSIONS

This study not only provides useful resources for exploring global miRNA expression patterns and genome structure but also facilitates genetic research on the roles of miRNAs in genomic interactions of Brassica allopolyploids.

Bruceine D Isolated from Brucea Javanica (L.) Merr. as a Systemic Feeding Deterrent for Three Major Lepidopteran Pests

Systemicity is a desirable property for insecticides. Many phytochemicals show good systemic properties and thus are natural sources of novel systemic insecticidal ingredients. Bruceine D, a quassinoid, was identified in Brucea javanica (L.) Merr. and displayed outstanding systemic properties and excellent antifeedant activity against the diamondback moth (DBM, Plutella xylostella L.), beet armyworm ( Spodoptera exigua Hübner), and cotton leafworm ( Spodoptera litura Fabricius). Its antifeedant effect on third instar larvae of DBM was approximately 6.2-fold stronger than that of azadirachtin. When bruceine D was applied to roots at a concentration of 100 μg/mL for 24 and 48 h, its concentration in flowering Chinese cabbage ( Brassica campestris L. ssp. chinensis var. utiliz Tsen et Lee) leaves was 38.69 μg/g (fresh weight, FW) and 108.45 μg/g (FW), respectively. These concentrations could achieve 93.80% and 96.83% antifeedant effects, which were significantly greater than those of azadirachtin. Similar to azadirachtin, bruceine D also posed a potent growth inhibition effect on insect larvae. After feeding with 20 μg/g bruceine D, no pupae were observed. The results demonstrated that bruceine D is an effective botanical insect antifeedant with outstanding systemic properties, causing potent pest growth inhibitory activity.

The transcriptome analysis of Protaetia brevitarsis Lewis larvae

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential.

Development and Characterization of EST-SSR Markers From RNA-Seq Data in Phyllostachys violascens

Bamboo are woody grass species containing important economic and ecological values. Lei bamboo (Phyllostachys violascens) is a kind of shoot-producing bamboo species with the highest economic yield per unit area. However, identifying different varieties of Lei bamboo based on morphological characteristics is difficult. Microsatellites play an important role in plant identification and genetic diversity analysis and are superior to other molecular markers. In this study, we identified 18,356 expressed sequence tag-simple sequence repeat (EST-SSR) loci in Lei bamboo transcriptome data. A total of 11,264 primer pairs were successfully designed from unigenes of all EST-SSR loci, and 96 primer pairs were randomly selected and synthesized. A total of 54 primer pairs were used for classifying 16 Lei bamboo varieties and 10 different Phyllostachys species. The number of polymorphism alleles among the 54 primer pairs ranged from 3 to 12 for P. violascens varieties and 3 to 20 for Phyllostachys. The phylogenetic tree based on polymorphism alleles successfully distinguished 16 P. violascens varieties and 10 Phyllostachys species. Our study provides abundant EST-SSR resources that are useful for genetic diversity analysis and molecular verification of bamboo and suggests that SSR markers developed from Lei bamboo are more efficient and reliable than ISSR, SRAP or AFLP markers.