Development of novel EST-SSR markers in the macaúba palm (Acrocomia aculeata) using transcriptome sequencing and cross-species transferability in Arecaceae species

SSR markers development and their application in genetic diversity of burdock (Arctium lappa L.) germplasm

BACKGROUND

Arctium lappa L. is a medicinal edible homologous plant, commonly known as burdock or bardana, which belongs to the Asteraceae family and is abundant all over the world. Genetic diversity assessment is essential for A. lappa germplasm resource conservation and breeding. The assessment techniques include morphological, biochemical, and DNA marker analysis. However, the limited number of available DNA markers is insufficient to conduct related genetic diversity assessment studies.

RESULTS

In this study, we conducted RNA sequencing of the A. lappa cultivar 'Yanagikawa Ideal' and developed SSR markers to characterize the genetic diversity and population structure of 56 A. lappa accessions and 8 wild relative accessions. A total of 4,851 simple sequence repeats (SSRs) loci were identified. The proportions of mono-, di- and tri-nucleotide repeat motifs were 30.40%, 21.50% and 33.10%, respectively. We developed and verified the reliability of 28 SSR core primer pairs through electronic polymerase chain reaction (ePCR) and the PCR amplification process. The polymorphism information content (PIC) values of the 28 SSR core primer pairs ranged from 0.246 to 0.848, with 14 pairs of SSR primers displaying high polymorphism (PIC > 0.5). The 28 SSR core primer pairs showed 100% mobility in Arctium tomentosum Miller and 96.43% mobility in Synurus deltoides (Aiton) Nakai, indicating their high versatility. The average Shannon information index (I) was 1.231, and the average observed heterozygosity (Ho) was 0.132, the average expected heterozygosity (He) was 0.564. The 64 accessions were divided into three clusters at a genetic distance of 0.558. AMOVA analysis shows 83% genetic variation within populations and 17% among populations, highlighting implications for conservation and breeding strategies.

CONCLUSION

Our study provides 28 newly high-quality SSR markers to enhance genetic resource conservation and breeding programs for A. lappa, as well as to support comparative genomics and cross-species breeding strategies for related species.

Genome-wide association insights into the genomic regions controlling vegetative and oil production traits in Acrocomia aculeata

BACKGROUND

Macauba (Acrocomia aculeata) is a non-domesticated neotropical palm that has been attracting attention for economic use due to its great potential for oil production comparable to the commercially used oil palm (Elaeis guineensis). The discovery of associations between quantitative trait loci and economically important traits represents an advance toward understanding its genetic architecture and can contribute to accelerating macauba domestication. Pursuing this advance, this study performs single-trait and multi-trait GWAS models to identify candidate genes associated with vegetative and oil production traits in macauba. Eighteen phenotypic traits were evaluated from 201 palms within a native population. Genotyping was performed with SNP markers, following the protocol of genotyping-by-sequencing. Given that macauba lacks a reference genome, SNP calling was performed using three different strategies: using i) de novo sequencing, ii) the Elaeis guineenses Jacq. reference genome and iii) the macauba transcriptome sequences. After quality control, we identified a total of 27,410 SNPs in 153 individuals for the de novo genotypic dataset, 10,444 SNPs in 158 individuals using the oil palm genotypic dataset, and 4,329 SNPs in 167 individuals using the transcriptome genotypic dataset. The GWAS analysis was then performed on these three genotypic datasets.

RESULTS

Statistical phenotypic analyses revealed significant differences across all studied traits, with heritability values ranging from 63 to 95%. This indicates that the population contains promising genotypes for selection and the initiation of breeding programs. Genetic correlations between the 18 traits ranged from -0.47 to 0.99. The total number of significant SNPs in the single-trait and multi-trait GWAS was 92 and 6 using the de novo genotypic dataset, 19 and 11 using the oil palm genotypic dataset, and 1 and 2 using the transcriptome genotypic dataset, respectively. Gene annotation identified 12 candidate genes in the single-trait GWAS and four in the multi-trait GWAS, across the 18 phenotypic traits studied, in the three genotypic datasets. Gene mapping of the macauba candidate genes revealed similarities with Elaeis guineensis and Phoenix dactylifera. The candidate genes detected are responsible for metal ion binding and transport, protein transportation, DNA repair, and other cell regulation biological processes.

CONCLUSIONS

We provide new insights into genomic regions that map candidate genes associated with vegetative and oil production traits in macauba. These potential candidate genes require confirmation through targeted functional analyses in the future, and multi-trait associations need to be scrutinized to investigate the presence of pleiotropic or linked genes. Markers linked to traits of interest could serve as valuable resources for the development of marker-assisted selection in macauba for its domestication and pre-breeding.

Training set optimization is a feasible alternative for perennial orphan crop domestication and germplasm management: an Acrocomia aculeata example

Orphan perennial native species are gaining importance as sustainability in agriculture becomes crucial to mitigate climate change. Nevertheless, issues related to the undomesticated status and lack of improved germplasm impede the evolution of formal agricultural initiatives. Acrocomia aculeata - a neotropical palm with potential for oil production - is an example. Breeding efforts can aid the species to reach its full potential and increase market competitiveness. Here, we present genomic information and training set optimization as alternatives to boost orphan perennial native species breeding using Acrocomia aculeata as an example. Furthermore, we compared three SNP calling methods and, for the first time, presented the prediction accuracies of three yield-related traits. We collected data for two years from 201 wild individuals. These trees were genotyped, and three references were used for SNP calling: the oil palm genome, de novo sequencing, and the A. aculeata transcriptome. The traits analyzed were fruit dry mass (FDM), pulp dry mass (PDM), and pulp oil content (OC). We compared the predictive ability of GBLUP and BayesB models in cross- and real validation procedures. Afterwards, we tested several optimization criteria regarding consistency and the ability to provide the optimized training set that yielded less risk in both targeted and untargeted scenarios. Using the oil palm genome as a reference and GBLUP models had better results for the genomic prediction of FDM, OC, and PDM (prediction accuracies of 0.46, 0.45, and 0.39, respectively). Using the criteria PEV, r-score and core collection methodology provides risk-averse decisions. Training set optimization is an alternative to improve decision-making while leveraging genomic information as a cost-saving tool to accelerate plant domestication and breeding. The optimized training set can be used as a reference for the characterization of native species populations, aiding in decisions involving germplasm collection and construction of breeding populations.

Establishment of Novel Simple Sequence Repeat (SSR) Markers from Chimonanthus praecox Transcriptome Data and Their Application in the Identification of Varieties

Chimonanthus praecox, a member of the Calycanthaceae family, is a unique, traditional, and famous flowering economic tree species in China. Despite the existence of several varieties, only a few cultivars have been formally named. Currently, expression sequence tag-simple sequence repeat (EST-SSR) markers are extensively used to identify different species and varieties; a large number of microsatellites can be identified from transcriptome databases. A total of 162,638 unigenes were assembled using RNA-seq; 82,778 unigenes were annotated using the Nr, Nt, Swiss-Prot, Pfam, GO, KOG, and KEGG databases. In total, 13,556 SSR loci were detected from 11,691 unigenes, with trinucleotide repeat motifs being the most abundant among the six repeat motifs. To develop the markers, 64,440 pairs of SSR primers with polymorphism potential were designed, and 75 pairs of primers were randomly selected for amplification. Among these markers, seven pairs produced amplified fragments of the expected size with high polymorphism. Using these markers, 12 C. praecox varieties were clustered into two monophyletic clades. Microsatellites in the transcriptome of C. praecox exhibit rich types, strong specificity, and great polymorphism potential. These EST-SSR markers serve as molecular technical methods for identifying different varieties of C. praecox and facilitate the exploration of a large number of candidate genes associated with important traits.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

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.

Anti-androgenic potential of the fruit extracts of certain Egyptian Sabal species and their genetic variability studies: a metabolomic-molecular modeling approach

This work aimed to evaluate the anti-androgenic activity of S. blackburniana Glazebrook, S. causiarum (O. F. Cook) Becc, and S. palmetto (Walter) Lodd. Ex Schult fruit extracts in rats using Hershberger assay. Furthermore, to annotate secondary metabolites using LC-HRMS technique, to investigate underlying mechanisms responsible for 5-α-reductase inhibitory activity in silico and to compare cytotoxic effects in vitro against human prostatic stromal myofibroblast (WPMY-1) and human benign prostatic hyperplasia (BPH-1) cell lines using MTT, 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (spectrophotometrically). The results showed significant anti-androgenic implications with varying degrees, markedly decreased sex organ weights, reduction in testosterone and increase in LH and FSH serum levels. Genetic diversity study ensured the correct genotype and revealed outperformance of SCoT compared with CBDP markers to interpret polymorphism among selected species. S. blackburniana exhibited selective cytotoxic activity against BPH-1 compared to finasteride. Molecular docking of 59 dereplicated metabolites belonging to various chemical classes revealed that helasaoussazine, pinoresinol and tetra-O-caffeoylquinic acid are the top inhibitors of 5-α-reductase-2. Our study provides an insight into the anti-androgenic activity of selected species of Egyptian Sabal supported by docking study for the first time, demonstrates safety toward liver and kidney and highlights a new potential therapeutic candidate for anti-androgenic related disease such as benign prostatic hyperplasia.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Genetic diversity and population structure analysis in a large collection of Vicia amoena in China with newly developed SSR markers

Vicia amoena is a high-nutritional quality forage similar to alfalfa. However, studies on the genetic background of V. amoena are scarce. In the present study, the genetic variation of 24 V. amoena populations was assessed with newly developed simple sequence repeat (SSR) markers. A total of 8799 SSRs were identified in the V. amoena genomic-enriched sequences, and the most abundant repeat number was four. A total of 569 sampled individuals were assayed to evaluate the genetic diversity of the V. amoena populations based on 21 polymorphic SSR primers. The polymorphism information content (PIC) ranged from 0.896 to 0.968, with an average of 0.931, which indicated that the markers were highly informative. Based on analysis of molecular variance, 88% of the variance occurred within populations, and the remaining 12% of the variance occurred among populations. The high degree of gene flow (Nm= 4.958) also showed slight differentiation among the V. amoena populations. The V. amoena populations were mainly clustered by steppe and mountain habitats based on principal coordinate analysis (PCoA) and STRUCTURE analysis. This indicated that the elevation and special habitat of geographical origins may be important factors affecting the clustered pattern of V. amoena populations. Neighbour-joining (NJ) analysis did not separate the populations well by geographical origin, which indicated that the genetic structure of V. amoena was complex and needs further study. Overall, our results showed that the newly developed SSR markers could benefit the V. amoena research community by providing genetic background information to help establish a foundation for breeding improvement and germplasm resource conservation.

Development of Genomic SSR for the Subtropical Hardwood Tree Dalbergia hupeana and Assessment of Their Transferability to Other Related Species

Dalbergia hupeana Hance (D. hupeana) is a precious hardwood tree of the genus Dalbergia. It is one of the few species widely distributed within subtropical areas and is important for timber production and forest restoration. At present, there is little published genetic information on D. hupeana. Therefore, we performed a genome survey using next generation sequencing (NGS) and developed a set of novel genomic SSR (gSSR) markers from the assembled data, and assessed the transferability of these markers to other Dalbergia species in Asia. The results of the genome survey show the genome size of D. hupeana to be about 664 Mb and highly heterozygous. The assembly of sequencing data produced 2,431,997 contigs, and the initial assembly of the NGS data alone resulted in contig N50 of 393 kb with a total of 720 Mb. A total of 127,742 perfect SSR markers were found in the assembled contigs. A total of 37 highly polymorphic and easily genotyped gSSR markers were developed in D. hupeana, while the majority of gSSR markers could be successfully transferred to nine other Dalbergia species in Asia. The transferability rate of gSSR markers was highest in D. balansae, which is more closely related to D. hupeana. Seven gSSR markers were able to be amplified in all tested species. In addition, a preliminary assessment of the genetic diversity of three tree species in the Dalbergia genus suggested a high level of genetic diversity within populations distributed in the subtropical area in China. However, the determination of the global status of their genetic variation still requires further and more comprehensive assessment. Our findings will enable further studies on the genetic diversity, phylogenetics, germplasm characterization, and taxonomy of various Dalbergia species.

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.

Transferability, development of simple sequence repeat (SSR) markers and application to the analysis of genetic diversity and population structure of the African fan palm (Borassus aethiopum Mart.) in Benin

Background

In Sub-Saharan Africa, Borassus aethiopum Mart. (African fan palm) is an important non-timber forest product-providing palm that faces multiple anthropogenic threats to its genetic diversity. However, this species is so far under-studied, which prevents its sustainable development as a resource. The present work is a first attempt at characterizing the genetic diversity and population structure of B. aethiopum across nine collection sites spanning the three climatic regions of Benin, West Africa, through the use of microsatellite markers.

Results

During a first phase we relied on the reported transferability of primers developed in other palm species. We find that, in disagreement with previously published results, only 22.5% of the markers tested enable amplification of B. aethiopum DNA and polymorphism detection is very low.

In a second phase, we generated a B. aethiopum-specific genomic dataset through high-throughput sequencing and used it for the de novo detection of microsatellite loci. Among the primer pairs targeting these, 11 detected polymorphisms and were further used for analyzing genetic diversity. Across the nine sites, expected heterozygosity (He) ranges from 0.263 to 0.451 with an overall average of 0.354, showing a low genetic diversity. Analysis of molecular variance (AMOVA) shows that within-site variation accounts for 53% of the genetic variation. Accordingly, the low number of migrants and positive values of the fixation index (F) in sites from both the Central (Sudano-Guinean) and the Southern (Guinean) climatic regions suggest limited gene flow between sites. The global correlation between genetic and geographic distances is weak; however, our clustering analyses indicate that B. aethiopum palms from Savè (Center) are genetically more similar to those from the North than to samples from other Central sites.

Conclusions

In the light of our results, we discuss the use of inter-species transfer vs. de novo development of microsatellite markers in genetic diversity analyses targeting under-studied species, and suggest future applications for our molecular resources. We propose that, while prominent short-range pollen and seed dispersal in Benin explain most of our results, gene flux between the Central and Northern regions, as a result of animal and/or human migrations, might underlie the Savè discrepancy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00955-y.

Characterization, validation, and cross-species transferability of EST-SSR markers developed from Lycoris aurea and their application in genetic evaluation of Lycoris species

BACKGROUND

The Lycoris genus includes many ornamentally and medicinally important species. Polyploidization and hybridization are considered modes of speciation in this genus, implying great genetic diversity. However, the lack of effective molecular markers has limited the genetic analysis of this genus.

RESULTS

In this study, mining of EST-SSR markers was performed using transcriptome sequences of L. aurea, and 839 primer pairs for non-redundant EST-SSRs were successfully designed. A subset of 60 pairs was randomly selected for validation, of which 44 pairs could amplify products of the expected size. Cross-species transferability of the 60 primer pairs among Lycoris species were assessed in L. radiata Hreb, L. sprengeri Comes ex Baker, L. chinensis Traub and L. anhuiensis, of which between 38 to 77% of the primers were able to amplify products in these Lycoris species. Furthermore, 20 and 10 amplification products were selected for sequencing verification in L. aurea and L. radiata respectively. All products were validated as expected SSRs. In addition, 15 SSRs, including 10 sequence-verified and 5 unverified SSRs were selected and used to evaluate the genetic diversity of seven L. radiata lines. Among these, there were three sterile lines, three fertile lines and one line represented by the offspring of one fertile line. Unweighted pair group method with arithmetic mean analysis (UPGMA) demonstrated that the outgroup, L. aurea was separated from L. radiata lines and that the seven L. radiata lines were clustered into two groups, consistent with their fertility. Interestingly, even a dendrogram with 34 individuals representing the seven L. radiata lines was almost consistent with fertility.

CONCLUSIONS

This study supplies a pool of potential 839 non-redundant SSR markers for genetic analysis of Lycoris genus, that present high amplification rate, transferability and efficiency, which will facilitate genetic analysis and breeding program in Lycoris.

Characterization and Development of Microsatellite Markers in Pseudotaxus chienii (Taxaceae) Based on Transcriptome Sequencing

Pseudotaxus chienii (Taxaceae) is an endangered conifer species endemic to China. However, a lack of suitable molecular markers hinders the genomic and genetic studies on this species. Here, we characterized and developed the microsatellite markers from a newly sequenced P. chienii transcriptome. A total of 21,835 microsatellite loci were detected from 161,131 non-redundant unigene sequences, and the frequency of SSRs was 13.55%, with an average of one SSR loci per 9.18 kb. Mono-nucleotide, di-nucleotide, and tri-nucleotide were the dominant repeat types, accounting for 50.06, 13.49, and 29.39% of the total SSRs, respectively. In terms of distribution location, the coding regions (CDS) with few microsatellites and mainly consisted of tri-nucleotides. There were significant differences in the length of microsatellite among genic regions and motif types. Functional annotation showed that the unigenes containing microsatellites had a wide range of biological functions, most of which were related to basic metabolism, and a few might be involved in expression regulation of gene and response to environmental stress. In addition, 375 primer pairs were randomly selected and synthesized for the amplification and validation of microsatellite markers. Seventy-seven primer pairs were successfully amplified and 40 primer pairs were found to be polymorphic. Finally, 20 pairs of primers with high polymorphism were selected to assess the genetic diversity in four P. chienii populations. In addition, the newly developed microsatellite markers exhibited high transferability (70%) in Amentotaxus argotaenia. Our study could enable further genetic diversity analysis and functional gene mining on Taxaceae.

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

Populus wulianensis is an endangered species endemic to Shandong Province, China. Despite the economic and ornamental value of this species, few genomics and genetic studies have been performed. In this study, we performed a relevant analysis of the full-length transcriptome sequencing data of P. wulianensis and obtained expressed sequence tag (EST)-simple sequence repeat (SSR) markers with polymorphisms that can be used for further genetic research. In total, 8.18 Gb (3,521,665) clean reads with an average GC content of 42.12% were obtained. From the corrected 64,737 high-quality isoforms, 42,323 transcript sequences were obtained after redundancy analysis with CD-HIT. Among these transcript sequences, 41,876 sequences were annotated successfully. A total of 23,539 potential EST-SSRs were identified from 16,057 sequences. Excluding mononucleotides, the most abundant motifs were trinucleotide SSRs (47.80%), followed by di- (46.80%), tetra- (2.98%), hexa- (1.58%) and pentanucleotide SSRs (0.84%). Among the 100 designed EST-SSRs, 18 were polymorphic with high PIC values (0.721 and 0.683) and could be used for analyses of the genetic diversity and population structure of P. wulianensis. These full-length transcriptome sequencing data will facilitate gene discovery and functional genomics research in P. wulianensis, and the novel EST-SSRs developed in our study will promote molecular-assisted breeding, genetic diversity and conservation biology research in this species.

Prebiotic potential of pulp and kernel cake from Jerivá (Syagrus romanzoffiana) and Macaúba palm fruits (Acrocomia aculeata)

The jerivá (Syagrus romanzoffiana) and the macaúba (Acrocomia aculeata) are palm trees of the Arecaceae family, widely distributed in tropical and subtropical areas of Latin America, which have a low production cost and high productivity throughout the year. Due to the high content of lipids, their fruits have been used for oil extraction, which generates byproducts such as the pulps and the kernel cakes, a nutritionally rich byproduct that can be added into human food and, may have prebiotic potential. Therefore, the objective of this work was to characterize and evaluate the prebiotic potential of jerivá pulp (JP), macaúba pulp (MP), jerivá kernel cake (JC) and macaúba kernel cake (MC). For this, the fruits characterization was carried out through proximate composition, phenolic compounds content, and antioxidant activity, besides evaluating the antimicrobial and fermentative capacity of Bifidobacterium lactis, Lactobacillus casei, and Lactobacillus acidophilus against Escherichia coli. Jerivá and macaúba pulps and kernel cakes presented high levels of dietary fiber (20.45% JP, 37.87% JC, 19.95% MP and 35.81% MC) and high antioxidant activity, especially JP, which also showed the high values found for ABTS and DPPH (2498.49 µMTrolox·g^-1 fruit and 96.97 g fruit·g^-1 DPPH, respectively), has a high total phenolic content (850.62 mg GAE·100 g^-1). Also, JP promoted a better growth of probiotic strains and a more relevant pH reduction when compared to the commercial prebiotic FOS. However, MP, JC, and MC were also able to favor the growth of the strains. Probiotic microorganisms were able to use JP, MP, JC, and MC and produced short-chain fatty acids such as lactic, propionic, butyric, and acetic acid, capable of promoting health benefits. Therefore, the byproducts from jerivá and macaúba oil extraction have characteristics that indicate their prebiotic potential, and maybe interesting components to increase the nutritional value of foods.

Genome-wide development of miRNA-based SSR markers in Cleistogenes songorica and analysis of their transferability to Gramineae/non-Gramineae species

Simple sequence repeat (SSR) markers are commonly used for many genetic applications, such as map construction, fingerprinting, and genetic diversity analyses, due to their high reproducibility, polymorphism, and abundance. Endogenous miRNAs play essential roles in plant development and gene expression under diverse biotic and abiotic stress conditions. In the present study, we predicted 110 miRNA-SSR primer pairs from 287 precursor miRNAs. Among 110 primer pairs, 85 were successfully amplified and examined for transferability to other Gramineae and non-Gramineae species. The results showed that all 82 primer pairs yielded unambiguous and strong amplification, and across the 23 studied Cleistogenes accessions, a total of 385 alleles were polymorphic. The number of alleles produced per primer varied from 3 to 11, with an average of 4.69 per locus. The expected heterozygosity (He) ranged from 0.44 to 0.88, with an average of 0.74 per locus, and the PIC (Polymorphism Information Content) values ranged from 0.34 to 0.87, with an average of 0.69 per locus. Furthermore, 1422 miRNA target genes were predicted and analyzed using the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. In conclusion, the results showed that an miRNA-based microsatellite marker system can be applicable for genetic diversity and marker-assisted breeding studies.

Genotyping by RAD Sequencing Analysis Assessed the Genetic Distinctiveness of Experimental Lines and Narrowed Down the Genomic Region Responsible for Leaf Shape in Endive (Cichorium endivia L.)

The characterization of genetic diversity in elite breeding stocks is crucial for the registration and protection of new varieties. Moreover, experimental population structure analysis and information about the genetic distinctiveness of commercial materials are essential for crop breeding programs. The purpose of our research was to assess the genetic relationships of 32 endive (Cichorium endivia L.) breeding lines, 18 from var. latifolium (escarole) and 14 from var. crispum (curly), using heterologous Cichorium intybus-derived simple sequence repeats (SSR) markers and single-nucleotide polymorphisms (SNP) markers. We found that 14 out of 29 SSR markers were successfully amplified, but only 8 of them were related to polymorphic loci. To overcome the limitation of the low number of informative SSR marker loci, an alternative SNP-based approach was employed. The 4621 SNPs produced by a restriction site-associated DNA marker sequencing approach were able to fully discriminate the 32 endive accessions; most importantly, as many as 50 marker loci were found to distinguish the curly group from the escarole group. Interestingly, 24 of the marker loci mapped within a peripheral segment of chromosome 8 of lettuce (Lactuca sativa L.), spanning a chromosomal region of 49.6 Mb. Following Sanger sequencing-based validation, three genes were determined to carry nonsynonymous SNPs, and one of them matched a putative ortholog of AtELP1, subunit 1 of the Elongator complex. Considering that several previously characterized Elongator complex subunit mutants exhibited elongated and/or curly leaf phenotypes, this gene should be taken into consideration for a better understanding of the underlying mechanism controlling leaf shape in endive.

Genetic Divergence and Relationship Among Opisthopappus Species Identified by Development of EST-SSR Markers

Opisthopappus Shih is an endemic and endangered genus restricted to the Taihang Mountains that has important ornamental and economic value. According to the Flora Reipublicae Popularis Sinicae (FRPS, Chinese version), this genus contains two species (Opisthopappus longilobus and Opisthopappus taihangensis), whereas in the Flora of China (English version) only one species O. taihangensis is present. The interspecific phylogenetic relationship remains unclear and undefined, which might primarily be due to the lack of specific molecular markers for phylogenetic analysis. For this study, 2644 expressed sequence tag-simple sequence repeats (EST-SSRs) from 33,974 unigenes using a de novo transcript assembly of Opisthopappus were identified with a distribution frequency of 7.78% total unigenes. Thereinto, mononucleotides (1200, 45.39%) were the dominant repeat motif, followed by trinucleotides (992, 37.52%), and dinucleotides (410, 15.51%). The most dominant trinucleotide repeat motif was ACC/GGT (207, 20.87%). Based on the identified EST-SSRs, 245 among 1444 designed EST-SSR primers were selected for the development of potential molecular markers. Among these markers, 63 pairs of primers (25.71%) generated clear and reproducible bands with expected sizes. Eventually, 11 primer pairs successfully amplified all individuals from the studied populations. Through the EST-SSR markers, a high level of genetic diversity was detected between Opisthopappus populations. A significant genetic differentiation between the O. longilobus and O. taihangensis populations was found. All studied populations were divided into two clusters by UPGMA, NJ, STRUCTURE, and PCoA. These results fully supported the view of the FRPS, namely, that O. longilobus and O. taihangensis should be regarded as two distinct species. Our study demonstrated that transcriptome sequences, as a valuable tool for the quick and cost-effective development of molecular markers, was helpful toward obtaining comprehensive EST-SSR markers that could contribute to an in-depth assessment of the genetic and phylogenetic relationships between Opisthopappus species.