Transcriptome analysis and codominant markers development in caper, a drought tolerant orphan crop with medicinal value

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The Capparis spinosa var. herbacea genome provides the first genomic instrument for a diversity and evolution study of the Capparaceae family

BACKGROUND

The caper bush Capparis spinosa L., one of the most economically important species of Capparaceae, is a xerophytic shrub that is well adapted to drought and harsh environments. However, genetic studies on this species are limited because of the lack of its reference genome.

FINDINGS

We sequenced and assembled the Capparis spinosa var. herbacea (Willd.) genome using data obtained from the combination of PacBio circular consensus sequencing and high-throughput chromosome conformation capture. The final genome assembly was approximately 274.53 Mb (contig N50 length of 9.36 Mb, scaffold N50 of 15.15 Mb), 99.23% of which was assigned to 21 chromosomes. In the whole-genome sequence, tandem repeats accounted for 19.28%, and transposable element sequences accounted for 43.98%. The proportion of tandem repeats in the C. spinosa var. herbacea genome was much higher than the average of 8.55% in plant genomes. A total of 21,577 protein-coding genes were predicted, with 98.82% being functionally annotated. The result of species divergence times showed that C. spinosa var. herbacea and Tarenaya hassleriana separated from a common ancestor 43.31 million years ago.

CONCLUSIONS

This study reported a high-quality reference genome assembly and genome features for the Capparaceae family. The assembled C. spinosa var. herbacea genome might provide a system for studying the diversity, speciation, and evolution of this family and serve as an important resource for understanding the mechanism of drought and high-temperature resistance.

Population Genetic Structure and Biodiversity Conservation of a Relict and Medicinal Subshrub Capparis spinosa in Arid Central Asia

As a Tertiary Tethyan relict, Capparis spinosa is a typical wind-preventing and sand-fixing deciduous subshrub in arid central Asia. Due to its medicinal and energy value, this species is at risk of potential threat from human overexploitation, habitat destruction and resource depletion. In this study, our purpose was to evaluate the conservation strategies of C. spinosa according to its genetic structure characteristics and genetic diversity pattern among 37 natural distributional populations. Based on genomic SNP data generated from dd-RAD sequencing, genetic diversity analysis, principal component analysis, maximum likelihood phylogenetic trees and ADMIXTURE clustering, the significant population structure and differentiation were explored. The results showed the following: (1) Six distinct lineages were identified corresponding to geographic locations, and various levels of genetic diversity existed among the lineages for the natural habitat heterogeneity or human interferences; (2) The lineage divergences were influenced by isolation by distances, vicariance and restricted gene flow under complex topographic and climatic conditions. Finally, for the preservation of the genetic integrity of C. spinosa, we suggest that conservation units should be established corresponding to different geographic groups, and that attention should be paid to isolated and peripheral populations that are experiencing biodiversity loss. Simultaneously, monitoring and reducing anthropogenic disturbances in addition to rationally and sustainably utilizing wild resources would be beneficial to guarantee population resilience and evolutionary potential of this xerophyte in response to future environmental changes.

Transcriptomic responses to drought stress in Polygonatum kingianum tuber

BACKGROUND

Polygonatum kingianum Coll. et Hemsl. is an important plant in Traditional Chinese Medicine. The extracts from its tubers are rich in polysaccharides and other metabolites such as saponins. It is a well-known concept that growing medicinal plants in semi-arid (or drought stress) increases their natural compounds concentrations. This study was conducted to explore the morpho-physiological responses of P. kingianum plants and transcriptomic signatures of P. kingianum tubers exposed to mild, moderate, and severe drought and rewatering.

RESULTS

The stress effects on the morpho-physiological parameters were dependent on the intensity of the drought stress. The leaf area, relative water content, chlorophyll content, and shoot fresh weight decreased whereas electrolyte leakage increased with increase in drought stress intensity. A total of 53,081 unigenes were obtained; 59% of which were annotated. We observed that 1352 and 350 core genes were differentially expressed in drought and rewatering, respectively. Drought stress driven differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis, flavonoid biosynthesis, starch and sucrose metabolism, and stilbenoid diarylheptanoid and gingerol biosynthesis, and carotenoid biosynthesis pathways. Pathways such as plant-pathogen interaction and galactose metabolism were differentially regulated between severe drought and rewatering. Drought reduced the expression of lignin, gingerol, and flavonoid biosynthesis related genes and rewatering recovered the tubers from stress by increasing the expression of the genes. Increased expression of carotenoid biosynthesis pathway related genes under drought suggested their important role in stress endurance. An increase in starch and sucrose biosynthesis was evident from transcriptomic changes under drought stress. Rewatering recovered the drought affected tubers as evident from the contrasting expression profiles of genes related to these pathways. P. kingianum tuber experiences an increased biosynthesis of sucrose, starch, and carotenoid under drought stress. Drought decreases the flavonoids, phenylpropanoids, gingerol, and lignin biosynthesis. These changes can be reversed by rewatering the P. kingianum plants.

CONCLUSIONS

These results provide a transcriptome resource for P. kingianum and expands the knowledge on the effect of drought and rewatering on important pathways. This study also provides a large number of candidate genes that could be manipulated for drought stress tolerance and managing the polysaccharide and secondary metabolites' contents in P. kingianum.

Down-regulation of aflatoxin biosynthetic genes in Aspergillus parasiticus by Heliopsis longipes roots and affinin for reduction of aflatoxin production

Affinin present in Heliopsis longipes roots has been identified as an anti-aflatoxin molecule. However, its mechanism of action has yet to be clarified. Aflatoxins biosynthesis involves not less than 27 enzymatic reactions. In this work, the genes aflG, aflH, aflI, aflK, aflL, aflM, aflO, aflP, and aflQ of the aflatoxins cluster and the aflS gene encoding an internal regulatory factor involved in aflatoxins biosynthesis in Aspergillus parasiticus, were studied by qRT-PCR. Results demonstrated that ethanolic extract of H. longipes roots and affinin inhibit aflatoxin biosynthesis and fungal growth in a dose-dependent manner. At 300 µg/mL, ethanolic extract and affinin presented the highest inhibition of radial growth (86% and 94%) and aflatoxin production (68% and 80%). The qRT-PCR analysis demonstrated that nine tested genes were down-regulated by affinin and ethanolic extract. The most down-regulated was the aflK, a gene that encodes an enzyme cyclase with double function during the aflatoxin biosynthesis. While no significant down-regulation was obtaining for aflH gene. Exposure to affinin also resulted in decreased transcript levels of the internal regulator factor aflS. Based on our results, a model showing the regulatory mechanism in aflatoxin biosynthesis and its role in gene expression was proposed. In conclusion, affinin modulates the expression of several aflatoxin biosynthetic genes, leading to mycotoxin biosynthesis inhibition. Therefore, H. longipes roots is a suitable candidate to developed control strategies via lowering gene expressions as a future perspective in reducing aflatoxin contamination.

Transcriptomics and Metabolomics Reveal Purine and Phenylpropanoid Metabolism Response to Drought Stress in Dendrobium sinense, an Endemic Orchid Species in Hainan Island

Drought stress is a bottleneck factor for plant growth and development, especially in epiphytic orchids that absorb moisture mainly from the air. Recent studies have suggested that there are complex transcriptional regulatory networks related to drought stress in Dendrobium sinense. In this study, the transcription and metabolite alterations involved in drought stress response in D. sinense were investigated through RNA-seq and metabolomics. A total of 856 metabolites were identified from stressed and control samples, with 391 metabolites showing significant differences. With PacBio and Illumina RNA sequencing, 72,969 genes were obtained with a mean length of 2,486 bp, and 622 differentially expressed genes (DEGs) were identified. Correlation analysis showed 7 differential genes, and 39 differential metabolites were involved in interaction networks. The network analysis of differential genes and metabolites suggested that the pathways of purine metabolism and phenylpropanoid biosynthesis may play an important role in drought response in D. sinense. These results provide new insights and reference data for culturally important medicinal plants and the protection of endangered orchids.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

The genetic legacy of fragmentation and overexploitation in the threatened medicinal African pepper-bark tree, Warburgia salutaris

The pepper-bark tree (Warburgia salutaris) is one of the most highly valued medicinal plant species worldwide. Native to southern Africa, this species has been extensively harvested for the bark, which is widely used in traditional health practices. Illegal harvesting coupled with habitat degradation has contributed to fragmentation of populations and a severe decline in its distribution. Even though the species is included in the IUCN Red List as Endangered, genetic data that would help conservation efforts and future re-introductions are absent. We therefore developed new molecular markers to understand patterns of genetic diversity, structure, and gene flow of W. salutaris in one of its most important areas of occurrence (Mozambique). In this study, we have shown that, despite fragmentation and overexploitation, this species maintains a relatively high level of genetic diversity supporting the existence of random mating. Two genetic groups were found corresponding to the northern and southern locations. Our study suggests that, if local extinctions occurred in Mozambique, the pepper-bark tree persisted in sufficient numbers to retain a large proportion of genetic diversity. Management plans should concentrate on maintaining this high level of genetic variability through both in and ex-situ conservation actions.

An Updated Checklist of the Sicilian Native Edible Plants: Preserving the Traditional Ecological Knowledge of Century-Old Agro-Pastoral Landscapes

The traditional use of native wild food plants (NWFP) may represent a valuable supplementary food source for the present and future generations. In Sicily, the use of wild plants in the human diet dates back to very ancient times and still plays an important role in some rural communities. Moreover, in this regard, the natural and cultural inheritance of this island is wealthy and diversified for several reasons. First, Sicily hosts a rich vascular flora, with 3,000 native and 350 endemic plants. Second, due to its central position in the Mediterranean, the island has acted as a veritable melting pot for the ethnobotanical knowledge of the rural communities of the entire basin. We reviewed all the available literature and, starting from such omnicomprehensive checklist, partially improved thanks to the data issuing from recent field investigations, we critically revised the whole species list, basing our review on field data issuing from interviews and on our expert knowledge. As a result, we provide a substantially updated list of 292 NWFP growing on the island. Further 34 species, reported as NWFP on previous papers were discarded because they are not native to Sicily, while 45 species were listed separately because their identity, occurrence and local use as food is doubtful and needs to be further investigated. Moreover, we tried to shed light on the ecology (growth form and preferential habitat) of the Sicilian NWFP, with special focus on crop wild relatives (CWR). Our preliminary ecological analyses point out that a high percentage of these plants are linked with the so-called 'cultural' landscapes, patchy semi-natural environments rich in ecotones, leading to the conclusion that the maintenance of century-old agro-pastoral practices may represent an effective way to preserve the local heritage of edible plants. Our study allowed to identify as much as 102 taxa of agronomic interest which could be tested as novel crops in order to face ongoing global changes and to comply with sustainable agriculture policies. Among them, 39 taxa show promising traits in terms of tolerance to one or more environmental stress factors, while 55 more are considered CWR and/or can be easily cultivated and/or show high productivity/yield potential.

Characterization of Sicilian rosemary (Rosmarinus officinalis L.) germplasm through a multidisciplinary approach

In Sicily, small differences exist between wild and cultivated rosemary biotypes; VOCs and genetic profiles may be a useful tool to distinguish them. A germplasm collection of Rosmarinus officinalis L. was harvested from 15 locations in Sicily. Eleven wild and four cultivated populations were collected and, due to the surveyed area covered, they can be considered as a representative panel of Sicilian genetic background of the species. Ex situ plant collection was transferred to the field cultivation in homogeneous conditions for characterizing through a multidisciplinary approach. The study included morphological traits observations (growth habitus, flower color, number and size of leaves, length and number of internodes), VOC profiles using HS-SPME, genome size by flow cytometry analysis, and genetic characterization by means of DNA and nuclear microsatellite (nSSR) investigation. To detect any pattern within- and among-populations variability, all morphological and chemical data were submitted to ANOVA, while clustering and structure population analysis were carried out using genetic profiles. The present work allowed us to distinguish rather well between wild and cultivated genotypes and to underline the biodiversity richness among rosemary Sicilian germplasm, never highlighted, useful for future breeding programs addressed to exploit this important resource.