Genome survey of pistachio (Pistacia vera L.) by next generation sequencing: Development of novel SSR markers and genetic diversity in Pistacia species

Identification of Sex-Associated Genetic Markers in Pistacia lentiscus var. chia for Early Male Detection

Pistacia lentiscus var. chia is a valuable crop for its high-added-value mastic, a resin with proven pharmaceutical and cosmeceutical properties harvested from the male tree trunk. To achieve the maximum economic benefits from the cultivation of male mastic trees, it is important to develop early sex diagnosis molecular tools for distinguishing the sex type. Thus far, the work on sex identification has focused on Pistacia vera with promising results; however, the low transferability rates of these markers in P. lentiscus necessitates the development of species-specific sex-linked markers for P. lentiscus var. chia. To our knowledge, this is the first report regarding: (i) the development of species-specific novel transcriptome-based markers for P. lentiscus var. chia and their assessment on male, female and monoecious individuals using PCR-HRM analysis, thus, introducing a cost-effective method for sex identification with high accuracy that can be applied with minimum infrastructure, (ii) the effective sex identification in mastic tree using a combination of different sex-linked ISSR and SCAR markers with 100% accuracy, and (iii) the impact evaluation of sex type on the genetic diversity of different P. lentiscus var. chia cultivars. The results of this study are expected to provide species-specific markers for accurate sex identification that could contribute to the selection process of male mastic trees at an early stage for mass propagation systems and to facilitate future breeding efforts related to sex-linked productivity and quality of mastic resin.

Genome-wide survey and characterization of microsatellites in cashew and design of a web-based microsatellite database: CMDB

The cashew is an edible tree nut crop having a wide range of food and industrial applications. Despite great economic importance, the genome-wide characterization of microsatellites [simple sequence repeats (SSRs)] in cashew is lacking. In this study, we carried out the first comprehensive genome-wide microsatellites/SSRs characterization in cashew and developed polymorphic markers and a web-based microsatellite database. A total of 54526 SSRs were discovered in the cashew genome, with a mean frequency of 153 SSRs/Mb. Among the mined genome-wide SSRs (2-6 bp size motifs), the dinucleotide repeat motifs were dominant (68.98%) followed by the trinucleotides (24.56%). The Class I type of SSRs (≥20 bp) were 45.10%, while Class II repeat motifs (≥12-<20 bp) were 54.89% of the total genomic SSRs discovered here. Further, the AT-rich SSRs occurred more frequently in the cashew genome (84%) compared to the GC-rich SSRs. The validation of the in silico-mined genome-wide SSRs by PCR screening in cashew genotypes resulted in the development of 59 polymorphic SSR markers, and the polymorphism information content (PIC) of the polymorphic SSR markers ranged from 0.19 to 0.84. Further, a web-based database, "Cashew Microsatellite Database (CMDB)," was constructed to provide access to the genome-wide SSRs mined in this study as well as transcriptome-based SSRs from our previous study to the research community through a user-friendly searchable interface. Besides, CMDB provides information on experimentally validated SSRs. CMDB permits the retrieval of SSR markers information with the customized search options. Altogether, the genome-wide SSRs characterization, the polymorphic markers and CMDB database developed in this study would serve as valuable marker resources for DNA fingerprinting, germplasm characterization, genetic studies, and molecular breeding in cashew and related Anacardium species.

Genetic Diversity and Relationship of Shanlan Upland Rice Were Revealed Based on 214 Upland Rice SSR Markers

Shanlan upland rice (Oryza sativa L.) is a unique upland rice variety cultivated by the Li nationality for a long time, which has good drought resistance and high utilization value in drought resistance breeding. To explore the origin of Shanlan upland rice and its genetic relationship with upland rice from other geographical sources, 214 upland rice cultivars from Southeast Asia and five provinces (regions) in southern China were used to study genetic diversity by using SSR markers. Twelve SSR primers were screened and 164 alleles (Na) were detected, with the minimum number of alleles being 8 and the maximum number of alleles being 23, with an average of 13.667. The analysis of genetic diversity and analysis of molecular variance (AMOVA) showed that the differences among the materials mainly came from the individuals of upland rice. The results of gene flow and genetic differentiation revealed the relationship between the upland rice populations, and Hainan Shanlan upland rice presumably originated from upland rice in Guangdong province, and some of them were genetically differentiated from Hunan upland rice. It can be indirectly proved that the Li nationality in Hainan is a descendant of the ancient Baiyue ethnic group, which provides circumstantial evidence for the migration history of the Li nationality in Hainan, and also provides basic data for the advanced protection of Shanlan upland rice, and the innovative utilization of germplasm resources.

Development of EST-SSRs based on the transcriptome of Castanopsis carlesii and cross-species transferability in other Castanopsis species

Castanopsis carlesii (Hemsl.) Hay. is a widely distributed and dominant tree species native to subtropical China with significant ecological and economic value. Due to serious human-related disturbance, its wild resources have been increasingly reduced, and whether may result in the loss of genetic diversity. However, no population genetics studies of natural C. carlesii have been reported to date. Microsatellite markers have been a useful tool in population genetics. Therefore, we developed EST-SSR markers based on the transcriptome sequencing of C. carlesii leaves. A total of 149,380,224 clean reads were obtained, and 63,012 nonredundant unigenes with a mean length of 1,034 bp were assembled and annotated based on sequence similarity searches in the Nr, Nt, KO, SwissProt, PFAM, KOG, and GO databases. The results showed that only 5,559 (8.82%) unigenes were annotated in all seven databases, but 46,338 (73.53%) could be annotated in at least one database. A total of 31,459 potential EST-SSRs were identified in 18,690 unigenes, with an average frequency of one SSR approximately 2 kb. Among the 100 EST-SSR primer pairs designed, 49 primer pairs successfully produced the expected product by amplification, with a success rate of 49%, but only 20 primer pairs showed abundant polymorphisms. Polymorphisms were verified using 25 samples from C. carlesii in Qimen, Anhui. A total of 119 alleles were detected, with a mean number of alleles (Na) of 5.95 per locus and a mean polymorphism information content (PIC) of 0.6125. All the 20 newly developed EST-SSR markers were verified in other Castanopsis species (C. sclerophylla, C. lamontii, C. fargesii, C. eyrei and C. jucunda). Sixteen primer pairs showed successful amplification in all five Castanopsis species (80%), and the transferability ratios ranged from 90% to 100%. These developed EST-SSR markers can be applied to population genetic and germplasm evaluations of C. carlesii and related species.

A first insight into the genomic background of Ilex pubescens (Aquifoliaceae) by flow cytometry and genome survey sequencing

BACKGROUND

Ilex pubescens is an important traditional Chinese medicinal plant with many naturally occurring compounds and multiple pharmacological effects. However, the lack of reference genomic information has led to tardiness in molecular biology research and breeding programs of this plant.

RESULTS

To obtain knowledge on the genomic information of I. pubescens, a genome survey was performed for the first time by next generation sequencing (NGS) together with genome size estimation using flow cytometry. The whole genome survey of I. pubescens generated 46.472 Gb of sequence data with approximately 82.2 × coverage. K-mer analysis indicated that I. pubescens has a small genome of approximately 553 Mb with 1.93% heterozygosity rate and 39.1% repeat rate. Meanwhile, the genome size was estimated to be 722 Mb using flow cytometry, which was possibly more precise for assessment of genome size than k-mer analysis. A total of 45.842 Gb clean reads were assembled into 808,938 scaffolds with a relatively short N50 of 760 bp. The average guanine and cytosine (GC) content was 37.52%. In total, 197,429 microsatellite motifs were detected with a frequency of 2.8 kb, among which mononucleotide motifs were the most abundant (up to 62.47% of the total microsatellite motifs), followed by dinucleotide and trinucleotide motifs.

CONCLUSION

In summary, the genome of I. pubescens is small but complex with a high level of heterozygosity. Even though not successfully applied for estimation of genome size due to its complex genome, the survey sequences will help to design whole genome sequencing strategies and provide genetic information support for resource protection, genetic diversity analysis, genetic improvement and artificial breeding of I. pubescens.

History, Phylogeny, Biodiversity, and New Computer-Based Tools for Efficient Micropropagation and Conservation of Pistachio (Pistacia spp.) Germplasm

The word "pstk" [pistag], used in the ancient Persian language, is the linguistic root from which the current name "pistachio", used worldwide, derives. The word pistachio is generally used to designate the plants and fruits of a single species: Pistacia vera L. Both the plant and its fruits have been used by mankind for thousands of years, specifically the consumption of its fruits by Neanderthals has been dated to about 300,000 years ago. Native to southern Central Asia (including northern Afghanistan and northeastern Iran), its domestication and cultivation occurred about 3000 years ago in this region, spreading to the rest of the Mediterranean basin during the Middle Ages and finally being exported to America and Australia at the end of the 19th century. The edible pistachio is an excellent source of unsaturated fatty acids, carbohydrates, proteins, dietary fiber, vitamins, minerals and bioactive phenolic compounds that help promote human health through their antioxidant capacity and biological activities. The distribution and genetic diversity of wild and domesticated pistachios have been declining due to increasing population pressure and climatic changes, which have destroyed natural pistachio habitats, and the monoculture of selected cultivars. As a result, the current world pistachio industry relies mainly on a very small number of commercial cultivars and rootstocks. In this review we discuss and summarize the current status of: etymology, origin, domestication, taxonomy and phylogeny by molecular analysis (RAPID, RFLP, AFLP, SSR, ISSR, IRAP, eSSR), main characteristics and world production, germplasm biodiversity, main cultivars and rootstocks, current conservation strategies of both conventional propagation (seeds, cutting, and grafting), and non-conventional propagation methods (cryopreservation, slow growth storage, synthetic seed techniques and micropropagation) and the application of computational tools (Design of Experiments (DoE) and Machine Learning: Artificial Neural Networks, Fuzzy logic and Genetic Algorithms) to design efficient micropropagation protocols for the genus Pistacia.

Genome assembly and association tests identify interacting loci associated with vigor, precocity, and sex in interspecific pistachio rootstocks

Understanding the basis of hybrid vigor remains a key question in crop breeding and improvement, especially for rootstock development where F1 hybrids are extensively utilized. Full-sibling UCB-1 F1 seedling rootstocks are widely planted in commercial pistachio orchards that are generated by crossing 2 highly heterozygous outbreeding parental trees of Pistacia atlantica (female) and P. integerrima (male). This results in extensive phenotypic variability, prompting costly removal of low-yielding small trees. To identify the genetic basis of this variability, we assembled chromosome-scale genome assemblies of the parental trees of UCB-1. We genotyped 960 UCB-1 trees in an experimental orchard for which we also collected multiyear phenotypes. We genotyped an additional 1,358 rootstocks in 6 commercial pistachio orchards and collected single-year tree-size data. Genome-wide single marker association tests identified loci associated with tree size and shape, sex, and precocity. In the experimental orchard, we identified multiple trait-associated loci and a strong candidate for ZZ/ZW sex chromosomes. We found significant marker associations unique to different traits and to early vs late phenotypic measures of the same trait. We detected 2 loci strongly associated with rootstock size in commercial orchards. Pseudo-testcross classification of markers demonstrated that the trait-associated alleles for each locus were segregating in the gametes of opposite parents. These 2 loci interact epistatically to generate the bimodal distribution of tree size with undesirable small trees observed by growers. We identified candidate genes within these regions. These findings provide a foundational resource for marker development and genetic selection of vigorous pistachio UCB-1 rootstock.

Genome Survey Sequencing and Genetic Background Characterization of Ilex chinensis Sims (Aquifoliaceae) Based on Next-Generation Sequencing

Ilex chinensis Sims. is an evergreen arbor species with high ornamental and medicinal value that is widely distributed in China. However, there is a lack of molecular and genomic data for this plant, which severely restricts the development of its relevant research. To obtain the whole reference genome, we first conducted a genome survey of I. chinensis by next-generation sequencing (NGS) to perform de novo whole-genome sequencing. As a result, our estimates using k-mer and flow cytometric analysis suggested the genome size of I. chinensis to be around 618-655 Mb, with the GC content, heterozygous rate, and repeat sequence rate of 37.52%, 1.1%, and 38%, respectively. A total of 334,649 microsatellite motifs were detected from the I. chinensis genome data, which will provide basic molecular markers for germplasm characterization, genetic diversity, and QTL mapping studies for I. chinensis. In summary, the I. chinensis genome is complex with high heterozygosity and few repeated sequences. Overall, this is the first report on the genome features of I. chinensis, and the information may lay a strong groundwork for future whole-genome sequencing and molecular breeding studies of this species.

Whole-Genome Survey Analyses Provide a New Perspective for the Evolutionary Biology of Shimofuri Goby, Tridentiger bifasciatus

The shimofuri goby (Tridentiger bifasciatus) is a small and highly adaptable goby, distributed along the coasts of China, the Sea of Japan, and the west coastal and estuarine areas of the Northwest Pacific. Next-generation sequencing was used to generate genome-wide survey data to provide essential characterization of the shimofuri goby genome and for the further mining of genomic information. The genome size of the shimofuri goby was estimated to be approximately 887.60 Mb through K-mer analysis, with a heterozygosity ratio and repeat sequence ratio of 0.47% and 32.60%, respectively. The assembled genome was used to identify microsatellite motifs (Simple Sequence Repeats, SSRs), extract single-copy homologous genes and assemble the mitochondrial genome. A total of 288,730 SSRs were identified. The most frequent SSRs were dinucleotide repeats (with a frequency of 61.15%), followed by trinucleotide (29.87%), tetranucleotide (6.19%), pentanucleotide (1.13%), and hexanucleotide repeats (1.66%). The results of the phylogenetic analysis based on single-copy homologous genes showed that the shimofuri goby and Rhinogobius similis can be clustered into one branch. The shimofuri goby was originally thought to be the same as the chameleon goby (Tridentiger trigonocephalus) due to their close morphological resemblance. However, a complete mitochondrial genome was assembled and the results of the phylogenetic analysis support the inclusion of the shimofuri goby as a separate species. PSMC analysis indicated that the shimofuri goby experienced a bottleneck event during the Pleistocene Glacial Epoch, in which its population size decreased massively, and then it began to recover gradually after the Last Glacial Maximum. This study provides a reference for the further assembly of the complete genome map of the shimofuri goby, and is a valuable genomic resource for the study of its evolutionary biology.

Genome survey sequencing of common vetch (Vicia sativa L.) and genetic diversity analysis of Chinese germplasm with genomic SSR markers

Background

Common vetch (Vicia sativa L.) is an annual legume with excellent suitability in cold and dry regions. Despite its great applied potential, the genomic information regarding common vetch currently remains unavailable.

Methods and results

In the present study, the whole genome survey of common vetch was performed using the next-generation sequencing (NGS). A total of 79.84 Gbp high quality sequence data were obtained and assembled into 3,754,145 scaffolds with an N50 length of 3556 bp. According to the K-mer analyses, the genome size, heterozygosity rate and GC content of common vetch genome were estimated to be 1568 Mbp, 0.4345 and 35%, respectively. In addition, a total of 76,810 putative simple sequence repeats (SSRs) were identified. Among them, dinucleotide was the most abundant SSR type (44.94%), followed by Tri- (35.82%), Tetra- (13.22%), Penta- (4.47%) and Hexanucleotide (1.54%). Furthermore, a total of 58,175 SSR primer pairs were designed and ten of them were validated in Chinese common vetch. Further analysis showed that Chinese common vetch harbored high genetic diversity and could be clustered into two main subgroups.

Conclusion

This is the first report about the genome features of common vetch, and the information will help to design whole genome sequencing strategies. The newly identified SSRs in this study provide basic molecular markers for germplasm characterization, genetic diversity and QTL mapping studies for common vetch.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06875-z.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

Genome-wide development and application of miRNA-SSR markers in Melilotus genus

Genetic diversity of plants is the brace of biodiversity and diversity within species, between species, and of ecosystems. SSR markers are the most preferable molecular marker tool that has been successfully used to study the genetic diversity of plant species. Development of miRNA-SSR markers has been deed in animals but is still limited in plants. In this study, 365 precursors miRNA were extracted from Melilotus albus (Ma) genome and used to design Ma miRNA-SSR primers. 137 Ma primer pairs (79 from known and 58 from novel pre-miRNAs) were obtained. 66 pairs of Ma miRNA-SSR primers were selected with polymorphisms and expected fragment size. The polymorphisms of primers were evaluated in 60 individuals of 15 Ma accessions. A total of 66 primer pairs showed high polymorphism, with average polymorphic information content of 0.49 among 15 Ma accessions and 0.63 among 18 Melilotus species, indicating that these primers have high polymorphisms. The number of alleles produced per primer ranged from 2 to 6 with an average of 3.6 alleles per locus in Ma accessions, and 2 to 10 numbers of alleles with a mean of 5.24 alleles per locus in Melilotus spp. For further studies, the genetic relationship was examined and the cluster analysis showed that 15 Ma accessions were grouped in three groups, on the other hand, 18 Melilotus species clustered into two groups. The analysis of molecular variance (AMOVA) revealed that 64.82% of the variation was found within the species and 35.18% between the species. The population structure analysis showed similar results with PCA analysis in that 18 species were grouped in two groups. In addition, 16,450 miRNA target genes were identified and used for GO and KEGG analysis. This is the first study to develop miRNA-SSR molecular markers in Melilotus spp., which has a great potential for marker-assisted, genetic improvement, genotyping applications, QTL analysis, and molecular-assisted selection studies for plant breeders and other researchers.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01086-z.

Whole genome resequencing data sets of different species from Pistacia genus

OBJECTIVES

Pistacia genus belongs to the flowering plants in the cashew family and contains at least 11 species. The whole-genome resequencing data of different species from Pistacia genus are described herein. The data reported here will be useful for better understand the adaptive evolution, demographic history, genetic diversity, population structure, and domestication of pistachio.

DATA DESCRIPTION

Genomic DNA was isolated from fresh leaves and used to construct libraries with insert size of 350 bp. Sequence libraries were made and sequenced on the Illumina Hiseq 4000 platform to produce 150 bp paired-end reads. A total number of 4,851,118,730 billion reads (ranging from 33,305,900 to 34,990,618 reads per sample) were created across all samples. We produced a total of 727.67 Gbp data which have been deposited in the Genome Sequence Archive (GSA) database with the Accession of CRA000978. All of the data are also available as the sequence read archive (SRA) format in the National Center for Biotechnology Information (NCBI) with identifier of SRP189222, mirroring our deposited data in GSA.

The karyotype, genome survey, and assembly of Mud artemisia (Artemisia selengensis)

BACKGROUND

Artemisia selengensis is traditional Chinese medicine and phytochemical analysis indicated that A. selengensis contains essential oils, fatty acids and phenolic acids. The lack of reference genomic information may lead to tardiness in molecular biology research of A. selengensis.

METHOD AND RESULTS

Karyotype analysis, genome survey, and genome assembly was employed to acquire information on the genome structure of A. selengensis. The chromosome number is 2n = 2x = 36, karyotype formula is 28 m + 8Sm, karyotype asymmetry coefficient is 58.8%, and karyotypes were symmetric to Stebbins' type 2A. Besides, the flow cytometry findings reported that the mean peak value of fluorescent intensity is 1,170,677, 2C DNA content is 12 pg and the genome size was estimated to be approximately 5.87 Gb. Furthermore, the genome survey generates 341,478,078 clean reads, unfortunately, after K-mer analysis, no significant peak can be observed, the heterozygosity, repetitive rate and genome size was unable to estimated. It is speculated that this phenomenon might be due to the complexity of genome structure. 37,266 contigs are preliminary assembled with Oxford Nanopore Technology (ONT) sequencing, totaling 804 Mb and GC content was 34.08%. The total length is 804,475,881 bp, N50 is 29,624 bp, and the largest contig length is 239,792 bp.

CONCLUSION

This study reveals the preliminary information of genome size of A. selengensis. These findings may provide supportive information for sequencing and assembly of whole-genome sequencing and encourage the progress of functional gene discovery, genetic improvement, evolutionary study, and structural studies of A. selengensis.

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.

Genomic characteristics and profile of microsatellite primers for Acanthogobius ommaturus by genome survey sequencing

Acanthogobius ommaturus is one of the suitable species to study the genetic mechanism of adaptive evolution, but there are few reports on its genetics. In the present study, the genomic survey sequencing method was used to analyze the genome characters of A. ommaturus. A total of 50.50 G high-quality sequence data were obtained in the present study. From the 19-mer distribution frequency, the estimated genome size was 928.01 Mb. The calculated sequence repeat rate was about 38.31%, the heterozygosity was approximately 0.17%, and the GC% content was approximately 40.88%. Moreover, 475,724 simple sequence repeats (SSRs) were identified. Among them, dinucleotide repeats were the most (53.70% of the total SSRs), followed by tri- (35.36%), hexa- (4.59%), tetra- (4.57%) and penta- (1.77%) nucleotide repeats type. This is the first genome-wide feature of this species to be reported.

Complete chloroplast genome of Myracrodruon urundeuva and its phylogenetics relationships in Anacardiaceae family

Continuous exploratory use of tree species is threatening the existence of several plants in South America. One of these threatened species is Myracroduron urundeuva, highly exploited due to the high quality and durability of its wood. The chloroplast (cp) has been used for several evolutionary studies as well traceability of timber origin, based on its gene sequences and simple sequence repeats (SSR) variability. Cp genome organization is usually consisting of a large single copy and a small single copy region separated by two inverted repeats regions. We sequenced the complete cp genome from M. urundeuva based on Illumina next-generation sequencing. Our results show that the cp genome is 159,883 bp in size. The 36 SSR identified ranging from mono- to hexanucleotides. Positive selection analysis revealed nine genes related to photosystem, protein synthesis, and DNA replication, and protease are under positive selection. Genome comparison a other Anacardiaceae chloroplast genomes showed great variability in the family. The phylogenetic analysis using complete chloroplast genome sequences of other Anacardiaceae family members showed a close relationship with two other economically important genera, Pistacia and Rhus. These results will help future investigations of timber monitoring and population and evolutionary studies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00989-1.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Whole genome survey analysis and microsatellite motif identification of Sebastiscus marmoratus

The marbled rockfish Sebastiscus marmoratus is an ecologically and economically important marine fish species distributed along the northwestern Pacific coast from Japan to the Philippines. Here, next-generation sequencing was used to generate a whole genome survey dataset to provide fundamental information of its genome and develop genome-wide microsatellite markers for S. marmoratus. The genome size of S. marmoratus was estimated as approximate 800 Mb by using K-mer analyses, and its heterozygosity ratio and repeat sequence ratio were 0.17% and 39.65%, respectively. The preliminary assembled genome was nearly 609 Mb with GC content of 41.3%, and the data were used to develop microsatellite markers. A total of 191,592 microsatellite motifs were identified. The most frequent repeat motif was dinucleotide with a frequency of 76.10%, followed by 19.63% trinucleotide, 3.91% tetranucleotide, and 0.36% pentanucleotide motifs. The AC, GAG, and ATAG repeats were the most abundant motifs of dinucleotide, trinucleotide, and tetranucleotide motifs, respectively. In summary, a wide range of candidate microsatellite markers were identified and characterized in the present study using genome survey analysis. High-quality whole genome sequence based on the “Illumina+PacBio+Hi-C” strategy is warranted for further comparative genomics and evolutionary biology studies in this species.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

Genome survey sequencing and genetic diversity of cultivated Akebia trifoliata assessed via phenotypes and SSR markers

Akebia trifoliata (Lardizabalaceae) is an important medicinal plant with multiple pharmacological effects. However, the lack of genomic information had limited the further excavation and utilization of this plant. An initial survey of the genome A. trifoliata was performed by next-generation sequencing, and then the genome size was inferred by flow cytometry. The whole genome survey of A. trifoliata generated 61.90 Gb of sequence data with approximately 95.51 × coverage. The genome size, heterozygosity and GC content obtained by k-mer analysis were almost 648.07 Mb, 0.72% and 36.11%, respectively. The genome size calculated by flow cytometry was 685.77 Mb, which was consistent with the results of genome survey. A total of 851,957 simple sequence repeats (SSR) were identified in the A. trifoliata genome. Twenty-eight phenotypic traits and thirty pairs of SSR primers were selected for the analysis of the genetic diversity of 43 accessions of cultivated A. trifoliata. The results showed that 216 bands were generated by 30 pairs of SSR primers, of which 189 (87.5%) were polymorphic. In addition, the phenotypes and SSR markers were used for cluster analysis of 43 cultivated accessions. The results of the two clustering methods were partially consistent. The genome survey of A. trifoliata demonstrated that the genome size of this plant was about 648.07 Mb. In the present study, the size and characteristics of the genome of A. trifoliata were reported for the first time, which greatly enriched the genomic resources of A. trifoliata for the further research and utilization.

Genome Survey Sequencing of In Vivo Mother Plant and In Vitro Plantlets of Mikania cordata

Mikania cordata, the only native congener of the invasive weed Mikania micrantha in China, is an ideal species for comparative study to reveal the invasion mechanism. However, its genome resources are lagging far behind its congener, which limits the comparative genomic analysis. Our goal is to characterize the genome of M. cordata by next-generation sequencing and propose a scheme for long-read genome sequencing. Previous studies have shown that the genomic resources of the host plant would be affected by the endophytic microbial DNA. An aseptic sample of M. cordata will ensure the proper genome in downstream analysis. Because endophytes are ubiquitous in the greenhouse-grown M. cordata, the in vitro culture with cefotaxime or timentin treatment was undertaken to obtain the aseptic plantlets. The in vivo mother plant and in vitro plantlets were used to survey the genome. The microbial contamination in M. cordata was recognized by blast search and eliminated from the raw reads. The decontaminated sequencing reads were used to predict the genome size, heterozygosity, and repetitive rate. The in vivo plant was so contaminated that microbes occupied substantial sequencing resources and misled the scaffold assembly. Compared with cefotaxime, treatment with timentin performed better in cultivating robust in vitro plantlets. The survey result from the in vitro plantlets was more accurate due to low levels of contamination. The genome size was estimated to be 1.80 Gb with 0.50% heterozygosity and 78.35% repetitive rate. Additionally, 289,831 SSRs were identified in the genome. The genome is heavily contaminated and repetitive; therefore, the in vitro culture technique and long-read sequencing technology are recommended to generate a high-quality and highly contiguous genome.

Genetic Variation in Damaged Populations of Pistacia atlantica Desf

The Atlas Pistachio tree, Pistacia atlantica Desf., has great importance in the ecological landscape of North Africa, due to its adaptive plasticity, as well as its use as a rootstock in the cultivation of the economically important species, Pistacia vera L. The conservation and valuation of this species require sampling and an assessment of its genetic variability. For the first time in North Africa, the inter-simple sequence repeats (ISSR) molecular marker has been used in genetic-diversity assessment and in the population relationships of P. atlantica subsp. atlantica. The ISSR markers tested showed 74.1% polymorphism, while molecular variance (AMOVA) analysis revealed a high percentage of the total genetic diversity of 55.7% among the four populations studied. Cluster analysis with neighbor-joining (NJ) and principal coordinate analysis (PCO) divided the study sites into four distinct groups according to their geographical locations (Tiaret, Batna, Djelfa, and Bechar). Isolation by distance or Mantel test gave a positive correlation of r = 0.86 between geographical and genetic distances. The results in this study indicate an absence of gene flow, implying that conservation efforts should be taken separately for each population.

Genome survey sequencing of Atractylodes lancea and identification of its SSR markers

Atractylodes lancea (Thunb.) DC. is a traditional Chinese medicine rich in sesquiterpenes that has been widely used in China and Japan for the treatment of viral infections. Despite its important pharmacological value, genomic information regarding A. lancea is currently unavailable. In the present study, the whole genome sequence of A. lancea was obtained using an Illumina sequencing platform. The results revealed an estimated genome size for A. lancea of 4,159.24 Mb, with 2.28% heterozygosity, and a repeat rate of 89.2%, all of which indicate a highly heterozygous genome. Based on the genomic data of A. lancea, 27,582 simple sequence repeat (SSR) markers were identified. The differences in representation among nucleotide repeat types were large, e.g., the mononucleotide repeat type was the most abundant (54.74%) while the pentanucleotide repeats were the least abundant (0.10%), and sequence motifs GA/TC (31.17%) and TTC/GAA (7.23%) were the most abundant among the dinucleotide and trinucleotide repeat motifs, respectively. A total of 93,434 genes matched known genes in common databases including 48,493 genes in the Gene Ontology (GO) database and 34,929 genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. This is the first report to sequence and characterize the whole genome of A. lancea and will provide a theoretical basis and reference for further genome-wide deep sequencing and SSR molecular marker development of A. lancea.

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

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

Comprehensive whole genome survey analyses of male and female brown-spotted flathead fish Platycephalus sp.1

The flathead fish Platycephalus sp.1 is an ecologically and commercially important marine fish in the northwestern Pacific with notable sexual differences in growth and development. Yet the genomic data of this species is lacking. In the present study, whole genome sequencing of two individuals (one male and one female) of Platycephalus sp.1 were conducted to provide fundamental genomic information. The genome sizes were estimated to be 674.96 Mb (male) and 684.15 Mb (female) by using k-mer analyses. The heterozygosity and repeat ratios suggested possible male heterogamety of this species. The draft genome sequences were initially assembled and genome-wide microsatellite motifs were identified. Besides, the complete mitochondrial genome sequences were assembled and the phylogenetic analyses genetically supported the validation of Platycephalus sp.1. The reported genomic data and genetic markers in this study could be useful in future comparative genomics and evolutionary biology studies.

Development of novel g-SSR markers in guava (Psidium guajava L.) cv. Allahabad Safeda and their application in genetic diversity, population structure and cross species transferability studies

Dearth of genomic resources particularly, microsatellite markers in nutritionally and commercially important fruit crop, guava necessitate the development of the novel genomic SSR markers through the library enrichment techniques. Three types of 3' -biotinylated oligonucleotide probes [(CT)14, (GT)12, and (AAC)8] were used to develop microsatellite enriched libraries. A total of 153 transformed colonies were screened of which 111 positive colonies were subjected for Sanger sequencing. The clones having more than five motif repeats were selected for primer designing and a total of 38 novel genomic simple sequence repeats could be identified. The g-SSRs had the motif groups ranging from monomer to pentamer out of which dimer group occurred the most (89.47%). Out of 38 g-SSRs markers developed, 26 were found polymorphic, which showed substantial genetic diversity among the guava genotypes including wild species. The average number of alleles per locus, major allele frequency, gene diversity, expected heterozygosity and polymorphic information content of 26 SSRs were 3.46, 0.56, 0.53, 0.29 and 0.46, respectively. The rate of cross-species transferability of the developed g-SSR loci varied from 38.46 to 80.77% among the studied wild Psidium species. Generation of N-J tree based on 26 SSRs grouped the 40 guava genotypes into six clades with two out-groups, the wild guava species showed genetic distinctness from cultivated genotypes. Furthermore, population structure analysis grouped the guava genotypes into three genetic groups, which were partly supported by PCoA and N-J tree. Further, AMOVA and PCoA deciphered high genetic diversity among the present set of guava genotypes including wild species. Thus, the developed novel g-SSRs were found efficient and informative for diversity and population structure analyses of the guava genotypes. These developed novel g-SSR loci would add to the new genomic resource in guava, which may be utilized in genomic-assisted guava breeding.

Intercontinental long-distance seed dispersal across the Mediterranean Basin explains population genetic structure of a bird-dispersed shrub

Long-distance dispersal (LDD) is a pivotal process for plants determining their range of distribution and promoting gene flow among distant populations. Most fleshy-fruited species rely on frugivorous vertebrates to disperse their seeds across the landscape. While LDD events are difficult to record, a few ecological studies have shown that birds move a sizeable number of ingested seeds across geographic barriers, such as sea straits. The foraging movements of migrant frugivores across distant populations, including those separated by geographic barriers, creates a constant flow of propagules that in turn shapes the spatial distributions of the genetic variation in populations. Here, we have analysed the genetic diversity and structure of 74 populations of Pistacia lentiscus, a fleshy-fruited shrub widely distributed in the Mediterranean Basin, to elucidate whether the Mediterranean Sea acts as a geographic barrier or alternatively whether migratory frugivorous birds promote gene flow among populations located on both sides of the sea. Our results show reduced genetic distances among populations, including intercontinental populations, and they show a significant genetic structure across an eastern-western axis. These findings are consistent with known bird migratory routes that connect the European and African continents following a north-southwards direction during the fruiting season of many fleshy-fruited plants. Further, approximate Bayesian analysis failed to explain the observed patterns as a result of historical population migrations at the end of Last Glacial Maximum. Therefore, anthropic and/or climatic changes that would disrupt the migratory routes of frugivorous birds might have genetic consequences for the plant species they feed upon.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20–30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

Whole genomes and transcriptomes reveal adaptation and domestication of pistachio

BACKGROUND

Pistachio (Pistacia vera), one of the most important commercial nut crops worldwide, is highly adaptable to abiotic stresses and is tolerant to drought and salt stresses.

RESULTS

Here, we provide a draft de novo genome of pistachio as well as large-scale genome resequencing. Comparative genomic analyses reveal stress adaptation of pistachio is likely attributable to the expanded cytochrome P450 and chitinase gene families. Particularly, a comparative transcriptomic analysis shows that the jasmonic acid (JA) biosynthetic pathway plays an important role in salt tolerance in pistachio. Moreover, we resequence 93 cultivars and 14 wild P. vera genomes and 35 closely related wild Pistacia genomes, to provide insights into population structure, genetic diversity, and domestication. We find that frequent genetic admixture occurred among the different wild Pistacia species. Comparative population genomic analyses reveal that pistachio was domesticated about 8000 years ago and suggest that key genes for domestication related to tree and seed size experienced artificial selection.

CONCLUSIONS

Our study provides insight into genetic underpinning of local adaptation and domestication of pistachio. The Pistacia genome sequences should facilitate future studies to understand the genetic basis of agronomically and environmentally related traits of desert crops.

Comprehensive Study of the Phenolic Compound Profile and Antioxidant Activity of Eight Pistachio Cultivars and Their Residual Cakes and Virgin Oils

Phenolic compounds of eight pistachio ( Pistacia vera L.) cultivars and their residual cakes and virgin oils (screw pressing) were studied using high-performance liquid chromatography-diode array detection-electrospray ionization-tandem mass spectrometry. A total of 25 compounds were identified and quantified for pistachio nuts and residual cakes, with the presence of five flavonols, six flavanols, and one gallotannin being reported for the first time. Total phenolics in pistachio nuts showed a concentration from 1359 mg/kg (Kastel) to 4507 mg/kg (Larnaka). Flavanols were the most abundant phenolics, at about 90%, with resulting procyanidin B1 and gallocatechin being the main phenolics, depending upon the cultivar. Other phenolic groups, such as anthocyanins (from 54 to 218 mg/kg), flavonols (from 76 to 130 mg/kg), flavanones (from 12 to 71 mg/kg), and gallotannins (from 4 to 46 mg/kg), were also identified. Residual cakes presented the same phenolic profile but with a concentration almost double because of the concentration effect caused by the oil separation. Virgin pistachio oils showed a very low phenolic content, with eriodyctiol being the only compound identified.

Genome survey and development of 18 microsatellite markers to assess genetic diversity in Spondias tuberosa Arruda Câmara (Anacardiaceae) and cross-amplification in congeneric species

Spondias tuberosa Arr. Cam belongs to the Anacardiaceae family, an economically important family of plants whose fruits are consumed by humans and animals. The aim of this study was to develop microsatellite markers using sequences from high-throughput sequencing and a magnetic bead enrichment method. The sequences were used to obtain contigs with a minimum of 500 nucleotides using Ray software and the mining of the simple sequence repeats (SSR) was performed with Phobos software, while the primers were designed by Primer3. We developed 18 polymorphic nuclear microsatellite markers and successfully cross-amplified them to three Spondias species. In S. tuberosa, the alleles ranged from 2 to 5 for each locus and Hardy-Weinberg equilibrium was found for 16 loci, with an expected and observed heterozygosity at 0.095-0.755 and 0.1-0.75, respectively. Cross-transferability was obtained for all loci in S. bahiensis, S. dulcis and S. purpurea. We concluded that the microsatellite markers developed in this study are useful in genetic population and conservation studies, as well as for investigating the hybrid origins of Spondias species.

Genome Survey Sequencing of Acer truncatum Bunge to Identify Genomic Information, Simple Sequence Repeat (SSR) Markers and Complete Chloroplast Genome

The Acer truncatum Bunge is a particular forest tree species found in the north of China. Due to the recent discovery that its seeds contain a considerable amount of nervonic acid, this species has received more and more attention. However, there have been no reports of the genome in this species. In this study, we report on the Acer truncatum genome sequence produced by genome survey sequencing. In total, we obtained 61.90 Gbp of high-quality data, representing approximately 116x coverage of the Acer truncatum genome. The genomic characteristics of Acer truncatum include a genome size of 529.88 Mbp, a heterozygosis rate of 1.06% and a repeat rate of 48.8%. A total of 392,961 high-quality genomic SSR markers were developed and a graphical map of the annotated circular chloroplast genome was generated. Thus far, this is the first report of de novo whole genome sequencing and assembly of Acer truncatum. We believe that this genome sequence dataset may provide a new resource for future genomic analysis and molecular breeding studies of Acer truncatum.

Genetic Evaluation of Natural Populations of the Endangered Conifer Thuja koraiensis Using Microsatellite Markers by Restriction-Associated DNA Sequencing

Thuja koraiensis Nakai is an endangered conifer of high economic and ecological value in Jilin Province, China. However, studies on its population structure and conservation genetics have been limited by the lack of genomic data. Here, 37,761 microsatellites (simple sequence repeat, SSR) were detected based on 875,792 de novo-assembled contigs using a restriction-associated DNA (RAD) approach. Among these SSRs, 300 were randomly selected to test for polymorphisms and 96 obtained loci were able to amplify a fragment of expected size. Twelve polymorphic SSR markers were developed to analyze the genetic diversity and population structure of three natural populations. High genetic diversity (mean NA = 5.481, HE = 0.548) and moderate population differentiation (pairwise Fst = 0.048–0.078, Nm = 2.940–4.958) were found in this species. Molecular variance analysis suggested that most of the variation (83%) existed within populations. Combining the results of STRUCTURE, principal coordinate, and neighbor-joining analysis, the 232 individuals were divided into three genetic clusters that generally correlated with their geographical distributions. Finally, appropriate conservation strategies were proposed to protect this species. This study provides genetic information for the natural resource conservation and utilization of T. koraiensis and will facilitate further studies of the evolution and phylogeography of the species.

Cross-species transferability of EST-SSR markers developed from the transcriptome of Melilotus and their application to population genetics research

Melilotus is one of the most important legume forages, but the lack of molecular markers has limited the development and utilization of Melilotus germplasm resources. In the present study, 151 M clean reads were generated from various genotypes of Melilotus albus using Illumina sequencing. A total of 19,263 potential EST-SSRs were identified from 104,358 unigene sequences. Moreover, 18,182 primer pairs were successfully designed, and 550 primer pairs were selected using criteria of base repeat type, fragment length and annealing temperature. In addition, 550 primer pairs were screened by using PCR amplification products and used to assess polymorphisms in 15 M. albus accessions. A total of 114 primer pairs were detected as being highly polymorphic, and the average polymorphism information content (PIC) value was 0.79. Furthermore, those 114 polymorphic primer pairs were used to evaluate the transferability to 18 species of the genus Melilotus, and 70 EST-SSR markers were found to be transferable among the 18 Melilotus species. According to the UPGMA dendrogram and STRUCTURE analysis, the 18 Melilotus species were classified into three clusters. This study offers a valuable resource for the genetic diversity and molecular assisted breeding of germplasm resources in the genus Melilotus.

A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery

Background

Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress.

Results

RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars.

Conclusion

This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3989-7) contains supplementary material, which is available to authorized users.