Transcriptome sequencing, annotation and polymorphism detection in the hop bush, Dodonaea viscosa

De novo transcriptome analysis and identification of defensive genes in garlic (Allium sativum L.) using high-throughput sequencing

BACKGROUND

Garlic (Allium sativum L.) is the second most widely cultivated Allium which is mainly grown in temperate regions and used as a flavoring agent in a wide variety of foods. Garlic contains various bioactive compounds whose metabolic pathways, plant-pathogen interactions, defensive genes, identify interaction networks, and functional genomics were not previously predicted in the garlic at the genomic level. To address this issue, we constructed two garlic Illumina 2000 libraries from tissues of garlic clove and leaf.

RESULTS

Approximately 43 million 125 bp paired-end reads were obtained in the two libraries. A total of 239,973 contigs were generated by de novo assembly of both samples and were compared with the sequences in the NCBI non-redundant protein database (Nr). In total, 42% of contigs were matched to known proteins in public databases including Nr, Gene Ontology (GO), and Cluster Orthologous Gene Database (COG), and then, contigs were mapped to 138 via functional annotation against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, a number of regulatory genes including the CCHC (Zn) family, followed by WD40, bromodomain, bZIP, AP2-EREBP, BED-type (Zn) proteins, and defense response proteins related to different conserved domains, such as RGA3, NBS-LRR, TIR-NBS-LRR, LRR, NBS-ARC, and CC-NBS-LRR were discovered based on the transcriptome dataset. We compared the ortholog gene family of the A. sativum transcriptome to A. thaliana, O. sativa, and Z. mays and found that 12,077 orthologous gene families are specific to A. sativum L. Furthermore, we identified genes involved in plant defense mechanisms, their protein-protein interaction network, and plant-pathogen interaction pathways.

CONCLUSIONS

Our study contains an extensive sequencing and functional gene-annotation analysis of A. sativum L. The findings provide insights into the molecular basis of TFs, defensive genes, and a reference for future studies on the genetics and breeding of A. sativum L.

Development and utilization of microsatellite markers to assess genetic variation coupled with modelling range shifts of Dodonaea viscosa (L.) Jacq. in isolated Taita Hills and Mount Kenya forests

BACKGROUND

Understanding genetic variation is critical for the protection and maintenance of fragmented and highly disturbed habitats. The Taita Hills of Kenya are the northernmost part of the Eastern Arc Mountains and have been identified as one of the world's top ten biodiversity hotspots. Over the past century the current forests in the Taita Hills have become highly fragmented. In order to appraise the influence of anthropological disturbance and fragmentation on plant species in these mountains, we studied the genetic variation and population structure of Dodonaea viscosa (L.) Jacq. (Sapindaceae), using newly developed microsatellite (SSR) markers, combined with ecological niche modelling analyses (ENMs).

METHODS AND RESULTS

We utilized the Illumina paired-end technology to sequence D. viscosa's genome and developed its microsatellite markers. In total, 646,428 sequences were analyzed, and 49,836 SSRs were identified from 42,638 sequences. A total of 18 out of 25 randomly selected primer pairs were designed to test polymorphism among 92 individuals across eight populations. The average observed heterozygosity and expected heterozygosity ranged from 0.119 to 0.982 and from 0.227 to 0.691, respectively. Analysis of molecular variance (AMOVA) revealed 78% variance within populations and only 20% among the eight populations. According to ENM results, D. viscosa's suitable habitats have been gradually reducing since the last glacial maximum (LGM), and the situation will worsen under the extreme pessimist scenario of (representative concentration pathway) RCP 8.5. Moreover, genetic diversity was significantly greater in larger fragments.

CONCLUSIONS

In the present study, we successfully developed and tested SSR markers for D. viscosa. Study results indicate that fragmentation would constitute a severe threat to plant forest species. Therefore, urgent conservation management of smaller fragmented patches is necessary to protect this disturbed region and maintain the genetic resources.

LC-MS-Based Metabolomics for the Chemosystematics of Kenyan Dodonaea viscosa Jacq (Sapindaceae) Populations

Dodonaea viscosa Jacq (Sapindaceae) is a medicinal plant with a worldwide distribution. The species has undergone enormous taxonomic changes which caused confusion amongst plant users. In Kenya, for example, two varieties are known to exist based on morphology, i.e., D. viscosa var. viscosa along the coast, and D. viscosa var. angustifolia in the Kenyan inland. These two taxa are recognized as distinct species in some reports. This prompted us to apply metabolomics to understand the relationship among naturally occurring populations of D. viscosa in Kenya, and to identify compounds that can assist in taxonomic delineation of the different varieties of D. viscosa from different parts of Kenya. The phytochemical variability of Kenyan D. viscosa var. angustifolia populations collected from four different geographical regions (Nanyuki, Machakos, Nairobi, and Narok) and one coastal D. viscosa var. viscosa (the Gazi) were analyzed by LC-MS using a metabolomics-driven approach. Four known compounds, two diterpenoids (dodonic acid (1), hautriwaic acid lactone (3), and two flavonoids (5,7,4',5'-tetrahydroxy-3,6,2'-trimethoxyflavone (2) and catechin (4)) were isolated and purified from the Gazi coastal collection. The presence of these compounds and their relative abundance in other populations was determined by LC-MS analyses. Multivariate statistical analyses of LC-MS data was used for the visualization of the patterns of variation and identification of additional compounds. Eleven discriminant compounds responsible for separating chemometric clusters were tentatively identified. In an antimicrobial assay, hautriwaic acid lactone (3) and catechin (4) were the most active compounds followed by the extract from the coastal (Gazi) population. The clustering pattern of the five populations of D. viscosa suggested that the metabolite profiles were influenced by geo-environmental conditions and did not support the current classification of D. viscosa based on morphology. This study disputes the current classification of D. viscosa in Kenya and recommends revision using tools such as molecular phylogenetics.

Transcriptome sequencing and de novo analysis of the northern snakehead, Ophiocephalus argus

Northern snakehead, Ophiocephalus argus Cantor, is an endemic freshwater fish in China. However, wild stocks of O. argus are dwindling sharply. Further, water conservancy projects, environmental pollution and human activities have caused the decrease of wild stocks, which has attracted much attention. Here, we have investigated the genomic information of O. argus using IlluminaHiseq 4000 sequencing. The transcriptomes of O. argus were sequenced by Illumina technology. A total of 67,564 sequences from 79,500,964 paired-end reads were generated, 33,710 unigenes were annotated based on protein databases (NCBI nonredundant (NR) databases). In total, 7182 unigenes had the clusters of orthologous group (COG) classifications, 33,710 unigenes were assigned to 59 gene ontology (GO) terms. Further, a total of 21,464 simple sequence repeats (SSRs) from 67,564 unigenes and 113,518 single nucleotide polymorphism (SNP) sites among 335 Mclean reads were yielded for O. argus based on a transcriptome-wide search. The new transcriptome data which is presented in this study for O. argus will provide valuable information for gene discovery and downstream applications, such as phylogenetic analysis, gene-expression profiling and identification of genetic markers (SSRs andSNP).

Analysis of the first Taraxacum kok-saghyz transcriptome reveals potential rubber yield related SNPs

Taraxacum kok-saghyz (TK) is a potential alternative crop for natural rubber (NR) production, due to its high molecular weight rubber, short breeding cycle, and diverse environmental adaptation. However, improvements in rubber yield and agronomically relevant traits are still required before it can become a commercially-viable crop. An RNA-Seq based transcriptome was developed from a pool of roots from genotypes with high and low rubber yield. A total of 55,532 transcripts with lengths over 200 bp were de novo assembled. As many as 472 transcripts were significantly homologous to 49 out of 50 known plant putative rubber biosynthesis related genes. 158 transcripts were significantly differentially expressed between high rubber and low rubber genotypes. 21,036 SNPs were different in high and low rubber TK genotypes. Among these, 50 SNPs were found within 39 transcripts highly homologous to 49 publically-searched rubber biosynthesis related genes. 117 SNPs were located within 36 of the differentially expressed gene sequences. This comprehensive TK transcriptomic reference, and large set of SNPs including putative exonic markers associated with rubber related gene homologues and differentially expressed genes, provides a solid foundation for further genetic dissection of rubber related traits, comparative genomics and marker-assisted selection for the breeding of TK.

Targeted capture to assess neutral genomic variation in the narrow-leaf hopbush across a continental biodiversity refugium

The Adelaide geosyncline, a mountainous region in central southern Australia, is purported to be an important continental refugium for Mediterranean and semi-arid Australian biota, yet few population genetic studies have been conducted to test this theory. Here, we focus on a plant species distributed widely throughout the region, the narrow-leaf hopbush, Dodonaea viscosa ssp. angustissima, and examine its genetic diversity and population structure. We used a hybrid-capture target enrichment technique to selectively sequence over 700 genes from 89 individuals across 17 sampling locations. We compared 815 single nucleotide polymorphisms among individuals and populations to investigate population genetic structure. Three distinct genetic clusters were identified; a Flinders/Gammon ranges cluster, an Eastern cluster, and a Kangaroo Island cluster. Higher genetic diversity was identified in the Flinders/Gammon Ranges cluster, indicating that this area is likely to have acted as a refugium during past climate oscillations. We discuss these findings and consider the historical range dynamics of these populations. We also provide methodological considerations for population genomics studies that aim to use novel genomic approaches (such as target capture methods) on non-model systems. The application of our findings to restoration of this species across the region are also considered.