Full-length transcriptome sequencing analysis and development of EST-SSR markers for the endangered species Populus wulianensis

Genetic diversity, population structure, and DNA fingerprinting of Ailanthus altissima var. erythrocarpa based on EST-SSR markers

Ailanthus altissima var. erythrocarpa is an A. altissima variety with high economic, ecological and ornamental value, but there have been no reports on the development of SSR primers for it. According to the SSR primer information provided by the transcriptome of A. altissima var. erythrocarpa, 120 individuals with different redness levels were used to screen polymorphic primers. Transcriptomic analysis revealed 10,681 SSR loci, of which mononucleotide repeats were dominant (58.3%), followed by dinucleotide and trinucleotide repeats (16.6%, 15.1%) and pentanucleotide repeats (0.2%). Among 140 pairs of randomly selected primers, nineteen pairs of core primers with high polymorphism were obtained. The average number of alleles (Na), average number of effective alleles (Ne), average Shannon's diversity index (I), average observed heterozygosity (Ho), average expected heterozygosity (He), fixation index (F) and polymorphic information content (PIC) were 11.623, 4.098, 1.626, 0.516, 0.696, 0.232 and 0.671, respectively. Nineteen EST-SSR markers were used to study the genetic diversity and population structure of A. altissima var. erythrocarpa. The phylogenetic tree, PCoA, and structure analysis all divided the tested resources into two categories, clearly showing the genetic variation between individuals. The population showed high genetic diversity, mainly derived from intraspecific variation. Among nineteen pairs of primers, 4 pairs (p33, p15, p46, p92) could effectively distinguish and be used for fingerprinting of the tested materials. This study is of great significance for genetic diversity analysis and molecular-assisted breeding of A. altissima var. erythrocarpa.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

Transcriptome Analysis Reveals Genes Associated with Flooding Tolerance in Mulberry Plants

Mulberry (Morus alba), a widely distributed economic plant, can withstand long-term flooding stress. However, the regulatory gene network underlying this tolerance is unknown. In the present study, mulberry plants were subjected to submergence stress. Subsequently, mulberry leaves were collected to perform quantitative reverse-transcription PCR (qRT-PCR) and transcriptome analysis. Genes encoding ascorbate peroxidase and glutathione S-transferase were significantly upregulated after submergence stress, indicating that they could protect the mulberry plant from flood damage by mediating ROS homeostasis. Genes that regulate starch and sucrose metabolism; genes encoding pyruvate kinase, alcohol dehydrogenase, and pyruvate decarboxylase (enzymes involved in glycolysis and ethanol fermentation); and genes encoding malate dehydrogenase and ATPase (enzymes involved in the TCA cycle) were also obviously upregulated. Hence, these genes likely played a key role in mitigating energy shortage during flooding stress. In addition, genes associated with ethylene, cytokinin, abscisic acid, and MAPK signaling; genes involved in phenylpropanoid biosynthesis; and transcription factor genes also showed upregulation under flooding stress in mulberry plants. These results provide further insights into the adaptation mechanisms and genetics of submergence tolerance in mulberry plants and could aid in the molecular breeding of these plants.

Optimization and application of non-native Phragmites australis transcriptome assemblies

Phragmites australis (common reed) has a cosmopolitan distribution and has been suggested as a model organism for the study of invasive plant species. In North America, the non-native subspecies (ssp. australis) is widely distributed across the contiguous 48 states in the United States and large parts of Canada. Even though millions of dollars are spent annually on Phragmites management, insufficient knowledge of P. australis impeded the efficiency of management. To solve this problem, transcriptomic information generated from multiple types of tissue could be a valuable resource for future studies. Here, we constructed forty-nine P. australis transcriptomes assemblies via different assembly tools and multiple parameter settings. The optimal transcriptome assembly for functional annotation and downstream analyses was selected among these transcriptome assemblies by comprehensive assessments. For a total of 422,589 transcripts assembled in this transcriptome assembly, 319,046 transcripts (75.5%) have at least one functional annotation. Within the transcriptome assembly, we further identified 1,495 transcripts showing tissue-specific expression pattern, 10,828 putative transcription factors, and 72,165 candidates for simple sequence repeats markers. The identification and analyses of predicted transcripts related to herbicide- and salinity-resistant genes were shown as two applications of the transcriptomic information to facilitate further research on P. australis. Transcriptome assembly and selection would be important for the transcriptome annotation. With this optimal transcriptome assembly and all relative information from downstream analyses, we have helped to establish foundations for future studies on the mechanisms underlying the invasiveness of non-native P. australis subspecies.

Development and validation of SSR markers related to flower color based on full-length transcriptome sequencing in Chrysanthemum

Chrysanthemum (Chrysanthemum moriforlium Ramat.) is one of the most popular flowers worldwide, with very high ornamental and economic values. However, the limitations of available DNA molecular markers and the lack of full genomic sequences hinder the study of genetic diversity and the molecular breeding of chrysanthemum. Here, we developed simple sequence repeat (SSR) from the full-length transcriptome sequences of chrysanthemum cultivar 'Hechengxinghuo'. A total of 11,699 SSRs with mono-, di-, tri-, tetra-, penta- and hexanucleotide repeats were identified, of which eight out of eighteen SSR loci identified based on sixteen transcripts participated in carotenoid metabolism or anthocyanin synthesis were validated as polymorphic SSR markers. These SSRs were used to classify 117 chrysanthemum accessions with different flower colors at the DNA and cDNA levels. The results showed that four SSR markers of carotenoid metabolic pathway divided 117 chrysanthemum accessions into five groups at cDNA level and all purple chrysanthemum accessions were in the group III. Furthermore, the SSR marker CHS-3, LCYE-1 and 3MaT may be related to green color and the PSY-1b marker may be related to yellow color. Overall, our work may be provide a novel method for mining SSR markers associated with specific traits.

Microsatellite analysis and polymorphic marker development based on the full-length transcriptome of Camellia chekiangoleosa

Camellia chekiangoleosa is a popular variety of Oil-camellia that has high oil production and ornamental value. Microsatellite (SSR) markers are the preferred tool for the molecular marker-assisted breeding of C. chekiangoleosa. By focusing on the problems of the low development efficiency of polymorphic SSR markers and the lack of available functional markers in Oil-camellia, we identified 97,510 SSR loci based on the full-length transcriptome sequence of C. chekiangoleosa. An analysis of SSR characteristics showed that mononucleotide (51.29%) and dinucleotide (34.36%) SSRs were the main repeat types. The main SSR distribution areas based on proportion covered were ordered as follows: 5'UTR > 3'UTR > CDS. By comparing our data with those in databases such as GO and KEGG, we obtained functional annotations of unigene sequences containing SSR sites. The data showed that the amplification efficiency of the SSR primers was 51.72%, and the development efficiency of polymorphic SSR primers was 26.72%. Experiments verified that dinucleotide and pentanucleotide SSRs located in UTR regions could produce more polymorphic markers. An investigation into the genetic diversity of several C. chekiangoleosa populations also suggested that the developed SSR markers had higher levels of polymorphism. This study will provide a reference and high-quality markers for the large-scale development of functional SSR markers and genetic research in Oil-camellia.

Full-length fruit transcriptomes of southern highbush (Vaccinium sp.) and rabbiteye (V. virgatum Ait.) blueberry

Background

Blueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries.

Results

We generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the ‘Draper’ genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in ‘Suziblue’ and ‘Powderblue’ respectively by combining the reads mapped to Northern Highbush blueberry ‘Draper’ genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both ‘Suziblue’ and ‘Powderblue’ with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in ‘Suziblue’ and ‘Powderblue’, respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant.

Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in ‘Powderblue’. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics.

Conclusions

A full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality–regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.

De novo transcriptome assembly, gene annotation, and EST-SSR marker development of an important medicinal and edible crop, Amomum tsaoko (Zingiberaceae)

BACKGROUND

Amomum tsaoko is a medicinal and food dual-use crop that belongs to the Zingiberaceae family. However, the lack of transcriptomic and genomic information has limited the understanding of the genetic basis of this species. Here, we performed transcriptome sequencing of samples from different A. tsaoko tissues, and identified and characterized the expressed sequence tag-simple sequence repeat (EST-SSR) markers.

RESULTS

A total of 58,278,226 high-quality clean reads were obtained and de novo assembled to generate 146,911 unigenes with an N50 length of 2002 bp. A total of 128,174 unigenes were successfully annotated by searching seven protein databases, and 496 unigenes were identified as annotated as putative terpenoid biosynthesis-related genes. Furthermore, a total of 55,590 EST-SSR loci were detected, and 42,333 primer pairs were successfully designed. We randomly selected 80 primer pairs to validate their polymorphism in A. tsaoko; 18 of these primer pairs produced distinct, clear, and reproducible polymorphisms. A total of 98 bands and 96 polymorphic bands were amplified by 18 pairs of EST-SSR primers for the 72 A. tsaoko accessions. The Shannon's information index (I) ranged from 0.477 (AM208) to 1.701 (AM242) with an average of 1.183, and the polymorphism information content (PIC) ranged from 0.223 (AM208) to 0.779 (AM247) with an average of 0.580, indicating that these markers had a high level of polymorphism. Analysis of molecular variance (AMOVA) indicated relatively low genetic differentiation among the six A. tsaoko populations. Cross-species amplification showed that 14 of the 18 EST-SSR primer pairs have transferability between 11 Zingiberaceae species.

CONCLUSIONS

Our study is the first to provide transcriptome data of this important medicinal and edible crop, and these newly developed EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity, and molecular marker-assisted selection in A. tsaoko.

Full-length transcriptome sequencing analysis and characterization, development and validation of microsatellite markers in Kengyilia melanthera

As a typical psammophyte of the Triticeae, Kengyilia melanthera possesses high feeding potential and great utilization values in desertification control in the Qinghai-Tibet Plateau. However, few gene function and genetic studies have been performed in K. melanthera. In this study, single-molecule real-time sequencing technology was used to obtain the full-length transcriptome sequence of K. melanthera, following the functional annotation of transcripts and prediction of coding sequences (CDSs), transcription factors (TFs), and long noncoding RNA (lncRNA) sequences. Meanwhile, a total of 42,433 SSR loci were detected, with 5'-UTRs having the most SSR loci and trinucleotide being the most abundant type. In total, 108,399 SSR markers were designed, and 300 SSR markers were randomly selected for diversity verification of K. melanthera. A total of 49 polymorphic SSR markers were used to construct the genetic relationships of 56 K. melanthera accessions, among which 21 SSR markers showed good cross-species transferability among the related species. In conclusion, the full-length transcriptome sequence of the K. melanthera will assist gene prediction and promote molecular biology and genomics research, and the polymorphic SSR markers will promote molecular-assisted breeding and related research of K. melanthera and its relatives.

Pieces of evidence of enhanced cellulose biosynthesis in the low-Cd cultivar and high expression level of transportation genes in the high-Cd cultivar of Lactuca sativa

To investigate the molecular mechanism of Cd-accumulating difference between Lactuca sativa cultivars, full-length transcriptome comparison, as well as biochemical validation, have been conducted between Cd pollution-safe cultivar (Cd-PSC, cv. LYDL) and high-Cd-accumulating cultivar (cv. HXDWQ). The full-length transcriptome of L. sativa cultivars was achieved for the first time. The results showed high Cd compartmentalization in the cell wall of cv. LYDL was ascribed to the enhanced cell wall biosynthesis under Cd stress, which was consistent with the high cellular debris Cd level (32.10-43.58%). The expression levels of transporter genes in cv. HXDWQ were about 1.19 to 1.21-fold higher than those in cv. LYDL, which was in accordance with the high ratio of easy migrative Cd chemical forms (68.59-81.98%), indicating the high Cd accumulation in the shoot of cv. HXDWQ was ascribed to the higher transportation capacity in cv. HXDWQ. Moreover, the Cd-induced endoplasmic reticulum (ER) stress was associated with the higher Cd detoxification and tolerance in cv. HXDWQ rather than in cv. LYDL. The study provides new insights into the Cd-induced transcriptomic difference between L. sativa cultivars and further contributes to the molecular breeding of L. sativa Cd-PSC.

De novo transcriptome assembly and analysis of gene expression in different tissues of moth bean (Vigna aconitifolia) (Jacq.) Marechal

BACKGROUND

The underutilized species Vigna aconitifolia (Moth Bean) is an important legume crop cultivated in semi-arid conditions and is valued for its seeds for their high protein content. It is also a popular green manure cover crop that offers many agronomic benefits including nitrogen fixation and soil nutrients. Despite its economic potential, genomic resources for this crop are scarce and there is limited knowledge on the developmental process of this plant at a molecular level. In the present communication, we have studied the molecular mechanisms that regulate plant development in V. aconitifolia, with a special focus on flower and seed development. We believe that this study will greatly enrich the genomic resources for this plant in form of differentially expressed genes, transcription factors, and genic molecular markers.

RESULTS

We have performed the de novo transcriptome assembly using six types of tissues from various developmental stages of Vigna aconitifolia (var. RMO-435), namely, leaves, roots, flowers, pods, and seed tissue in the early and late stages of development, using the Illumina NextSeq platform. We assembled the transcriptome to get 150938 unigenes with an average length of 937.78 bp. About 79.9% of these unigenes were annotated in public databases and 12839 of those unigenes showed a significant match in the KEGG database. Most of the unigenes displayed significant differential expression in the late stages of seed development as compared with leaves. We annotated 74082 unigenes as transcription factors and identified 12096 simple sequence repeats (SSRs) in the genic regions of V.aconitifolia. Digital expression analysis revealed specific gene activities in different tissues which were validated using Real-time PCR analysis.

CONCLUSIONS

The Vigna aconitifolia transcriptomic resources generated in this study provide foundational resources for gene discovery with respect to various developmental stages. This study provides the first comprehensive analysis revealing the genes involved in molecular as well as metabolic pathways that regulate seed development and may be responsible for the unique nutritive values of moth bean seeds. Hence, this study would serve as a foundation for characterization of candidate genes which would not only provide novel insights into understanding seed development but also provide resources for improved moth bean and related species genetic enhancement.

Full-Length Transcriptome Characterization and Comparative Analysis of Chosenia arbutifolia

As a unique tree species in the Salicaceae family, Chosenia arbutifolia is used primarily for construction materials and landscape planting in China. Compared with other Salicaceae species members, the genomic resources of C. arbutifolia are extremely scarce. Thus, in the present study, the full-length transcriptome of C. arbutifolia was sequenced by single-molecular real-time sequencing (SMRT) technology based on the PacBio platform. Then, it was compared against those of other Salicaceae species. We generated 17,397,064 subreads and 95,940 polished reads with an average length of 1812 bp, which were acquired through calibration, clustering, and polishing. In total, 50,073 genes were reconstructed, of which 48,174 open reading frames, 4281 long non-coding RNAs, and 3121 transcription factors were discovered. Functional annotation revealed that 47,717 genes had a hit in at least one of five reference databases. Moreover, a set of 12,332 putative SSR markers were screened among the reconstructed genes. Single-copy and special orthogroups, and divergent and conserved genes, were identified and analyzed to find divergence among C. arbutifolia and the five Salicaceae species. To reveal genes involved in a specific function and pathway, enrichment analyses for GO and KEGG were also performed. In conclusion, the present study empirically confirmed that SMRT sequencing realistically depicted the C. arbutifolia transcriptome and provided a comprehensive reference for functional genomic research on Salicaceae species.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

Multi-Omics Analysis of Gene and Protein Candidates Possibly Related to Tetrodotoxin Accumulation in the Skin of Takifugu flavidus

Pufferfish is increasingly regarded by many as a delicacy. However, the tetrodotoxin (TTX) that accumulates in its body can be lethal upon consumption by humans. TTX is known to mainly accumulate in pufferfish skin, but the accumulation mechanisms are poorly understood. In this study, we aimed to explore the possible mechanism of TTX accumulation in the skin of the pufferfish Takifugu flavidus following treatment with TTX. Through liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis, we detected 37.3% of toxin accumulated in the skin at the end of the rearing period (168 h). Transcriptome and proteome analyses revealed the mechanism and pathways of TTX accumulation in the skin of T. flavidus in detail. Gene ontology and the Kyoto Encyclopedia of Genes and Genomes analyses strongly suggest that cardiac muscle contraction and adrenergic signaling in cardiomyocyte pathways play an important role in TTX accumulation. Moreover, some upregulated and downregulated genes, which were determined via RNA-Seq, were verified with qPCR analysis. This study is the first to use multi-omics profiling data to identify novel regulatory network mechanisms of TTX accumulation in the skin of pufferfish.

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Genetic Diversity and Population Structure of Capirona (Calycophyllum spruceanum Benth.) from the Peruvian Amazon Revealed by RAPD Markers

Capirona (Calycophyllum spruceanum Benth.) is a tree species of commercial importance widely distributed in South American forests that is traditionally used for its medicinal properties and wood quality. Studies on this tree species have been focused mainly on wood properties, propagation, and growth. However, genetic studies on capirona have been very limited to date. Currently, it is possible to explore genetic diversity and population structure in a fast and reliable manner by using molecular markers. We here used 10 random amplified polymorphic DNA (RAPD) markers to analyze the genetic diversity and population structure of 59 samples of capirona that were sampled from four provinces located in the eastern region of the Peruvian amazon. A total of 186 bands were manually scored, generating a 59 × 186 presence/absence matrix. A dendrogram was generated using the UPGMA clustering algorithm, and, similar to the principal coordinate analysis (PCoA), it showed four groups that correspond to the geographic origin of the capirona samples (LBS, Irazola, Masisea, Iñapari). Similarly, a discriminant analysis of principal components (DAPC) and STRUCTURE analysis confirmed that capirona is grouped into four clusters. However, we also noticed that a few samples were intermingled. Genetic diversity estimation was conducted considering the four groups (populations) identified by STRUCTURE software. AMOVA revealed the greatest variation within populations (71.56%) and indicated that variability among populations is 28.44%. Population divergence (Fst) between clusters 1 and 4 revealed the highest genetic difference (0.269), and the lowest Fst was observed between clusters 3 and 4 (0.123). RAPD markers were successful and effective. However, more studies are needed, employing other molecular tools. To the best of our knowledge, this is the first investigation employing molecular markers in capirona in Peru considering its natural distribution, and as such it is hoped that this helps to pave the way towards its genetic improvement and the urgent sustainable management of forests in Peru.

The Full-Length Transcriptome Sequencing and Identification of Na+/H+ Antiporter Genes in Halophyte Nitraria tangutorum Bobrov

Nitraria tangutorum Bobrov is a halophyte that is resistant to salt and alkali and is widely distributed in northwestern China. However, its genome has not been sequenced, thereby limiting studies on this particular species. For species without a reference genome, the full-length transcriptome is a convenient and rapid way to obtain reference gene information. To better study N. tangutorum, we used PacBio single-molecule real-time technology to perform full-length transcriptome analysis of this halophyte. In this study, a total of 21.83 Gb of data were obtained, and 198,300 transcripts, 51,875 SSRs (simple sequence repeats), 55,574 CDS (coding sequence), and 74,913 lncRNAs (long non-coding RNA) were identified. In addition, using this full-length transcriptome, we identified the key Na⁺/H⁺ antiporter (NHX) genes that maintain ion balance in plants and found that these are induced to express under salt stress. The results indicate that the full-length transcriptome of N. tangutorum can be used as a database and be utilized in elucidating the salt tolerance mechanism of N. tangutorum.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.