De novo assembly and characterization of transcriptome using Illumina sequencing and development of twenty five microsatellite markers for an endemic tree Juglans hopeiensis Hu in China

Genome assembly and microsatellite marker development using Illumina and PacBio sequencing in Persicaria maackiana (Polygonaceae) from Korea

BACKGROUND

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level.

OBJECTIVE

Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina.

METHODS

Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers.

RESULTS

The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754.

CONCLUSION

Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

Genetic Diversity of Juglans mandshurica Populations in Northeast China Based on SSR Markers

Juglans mandshurica is a native tree species in Northeast China. Due to habitat destruction and human disturbance, its population size has sharply decreased. Currently, information on molecular markers of J. mandshurica is limited and cannot meet the needs of germplasm resource evaluation and molecular marker-assisted breeding of J. mandshurica. Based on transcriptomic data from three tissues (leaves, bark, and fruit pericarp), we developed expressed sequence tag-simple sequence repeats (EST-SSRs) for J. mandshurica, and 15 polymorphic EST-SSR primers were initially selected. The average number of alleles (Na), expected heterozygosity (He), and the polymorphic information content (PIC) at different loci were 18.27, 0.670, and 0.797, respectively. Population genetic diversity analysis revealed that the average Na, He, and Shannon information indices (I) for 15 J. mandshurica populations were 6.993, 0.670, and 1.455, respectively. Among them, population Hunchun exhibited the highest genetic diversity (Na = 7.933, He = 0.723, and I = 1.617), while population Heihe exhibited the lowest genetic diversity (Na = 4.200, He = 0.605, and I = 1.158). STRUCTURE analysis, neighbor-joining method cluster analysis, and principal coordinate analysis showed that the 343 individuals of J. mandshurica from 15 populations were clustered into three categories. Category 1 (green) had 147 individuals from eight populations in Qingyuan, Caohekou, Jian, Ningan, Yongji, Baishishan, Helong, and Maoershan; category 2 (blue) had 81 individuals from three populations in Hulin, Boli, and Sanchazi; and category 3 (red) had 115 individuals from four populations in Heihe, Hunchun, Fangzheng, and Liangshui. Analysis of molecular variance (AMOVA) showed that genetic variations among and within individuals accounted for 16.22% and 21.10% of the total genetic variation, respectively, indicating that genetic variations within populations were greater than genetic variations among populations. The average genetic differentiation coefficient (Fst) and gene flow (Nm) between different populations were 0.109 and 4.063, respectively, implying moderate levels of genetic differentiation and gene flow. Based on the genetic diversity characteristics of different populations, we proposed various genetic conservation strategies for J. mandshurica.

Phenological Differences, Genetic Diversity, and Population Structure of Genotypes Obtained from Seeds of Kaman-1 Walnut Cultivar

Walnut (Juglans regia L.) is a diploid (2n = 32), deciduous, monoecious, and generally open-pollinated tree with nuts of high nutrient content. In this study, the phenological differences, genetic diversity, and population structure of Kaman-1 and its 79 progenies obtained by open pollination were characterized by ISSR primers and some important phenological traits. As a result of the phenological observations, it was determined that the progenies differ significantly from Kaman-1. Besides, using ISSR primers, walnut genotypes were found to have genetic similarities ranging from 0.52 to 0.99. UPGMA cluster analysis showed that accessions from 2 different groups were classified, and population structure analysis confirmed this finding. Based on the results, a significant variation both phenologically and genetically was found within the walnut accessions. Also, this study confirmed that the progenies obtained from the Kaman-1 walnut cultivar have a quite wide variation and that ISSR primers and phenological traits are an important tool in determining genetic diversity.

Development of EST-SSRs from the ark shell (Scapharca broughtonii) transcriptome and their application in genetic analysis of four populations

BACKGROUND

The ark shell (Scapharca broughtonii) is one of the most economically important mollusks in the Bohai Sea and Yellow Sea of China. In recent years, ark shells from the Korean population were introduced to China for seed propagation and culture.

OBJECTIVE

To explore the impact of the introduction of Korean ark shell on the genetic diversity of native population in China.

METHODS

Simple sequence repeat (SSR) is effective and widely used tool for genetic analysis. In this study, 180 EST-SSRs were selected and verified by polymerase chain reaction and polyacrylamide gel electrophoresis. Subsequently, five polymorphic EST-SSRs were screened and their primers were modified by fluorescein for use in the genetic analysis of four populations.

RESULTS

Genetic analysis showed that 361 alleles amplified by five SSR loci were detected in the four populations. The number of alleles for the five SSRs ranged from 8 to 30, with a mean of 18.05 (standard deviation, SD = 6.492). The effective number of alleles varied from 2.253 to 22.222, with a mean of 10.596 (SD = 4.713). Observed heterozygosity and expected heterozygosity were 0.167-0.833 and 0.566-0.971, with average values of 0.520 (SD = 0.177) and 0.891 (SD = 0.062), respectively. Polymorphic information content ranged from 0.521 to 0.953, with a mean of 0.865 (SD = 0.070). The pairwise genetic differentiation coefficient (F_ST) of the four populations ranged from 0.0267 to 0.0477, showing low genetic differentiation. The phylogenetic tree constructed by neighbor-joining method showed that the genetic distance between the Chinese Dalian native population and three Korean populations was relatively more far than that among those Korean populations.

CONCLUSION

The results indicated that the genetic structure of the Dalian wild population was less affected by the introduced Korean wild populations.

Combining independent de novo assemblies to optimize leaf transcriptome of Persian walnut

Transcriptome resources can facilitate to increase yield and quality of walnuts. Finding the best transcriptome assembly has not been the subject of walnuts research as yet. This research generated 240,179,782 reads from 11 walnut leaves according to cDNA libraries. The reads provided a complete de novo transcriptome assembly. Fifteen different transcriptome assemblies were constructed from five different well-known assemblers used in scientific literature with different k-mer lengths (Bridger, BinPacker, SOAPdenovo-Trans, Trinity and SPAdes) as well as two merging approaches (EvidentialGene and Transfuse). Based on the four quality metrics of assembly, the results indicated an efficiency in the process of merging the assemblies after being generated by de novo assemblers. Finally, EvidentialGene was recognized as the best assembler for the de novo assembly of the leaf transcriptome in walnut. Among a total number of 183,191 transcripts which were generated by EvidentialGene, there were 109,413 transcripts capable of protein potential (59.72%) and 104,926 were recognized as ORFs (57.27%). In addition, 79,185 transcripts were predicted to exist with at least one hit to the Pfam database. A number of 3,931 transcription factors were identified by BLAST searching against PlnTFDB. Furthermore, 6,591 of the predicted peptide sequences contained signaling peptides, while 92,704 contained transmembrane domains. Comparison of the assembled transcripts with transcripts of the walnut and published genome assembly for the 'Chandler' cultivar using the BLAST algorithm led to identify a total number of 27,304 and 19,178 homologue transcripts, respectively. De novo transcriptomes in walnut leaves can be developed for the future studies in functional genomics and genetic studies of walnuts.

Population genetic structure and adaptive differentiation of iron walnut Juglans regia subsp. sigillata in southwestern China

Southwestern (SW) China is an area of active tectonism and erosion, yielding a dynamic, deeply eroded landscape that influences the genetic structure of the resident populations of plants and animals. Iron walnut (Juglans regia subsp. sigillata) is a deciduous tree species endemic to this region of China and cultivated there for its edible nuts. We sampled 36 iron walnut populations from locations throughout the species' range in SW China and genotyped a total of 765 individuals at five chloroplast DNA regions and 22 nuclear microsatellite loci. Species distribution models were produced to predict the evolution and historical biogeography of iron walnut and to estimate the impacts of climate oscillations and orographic environments on the species' demography. Our results indicated that J. regia subsp. sigillata had relatively low genetic diversity, high interpopulation genetic differentiation, and asymmetric interpopulation gene flow. Based on DIYABC analysis, we identified two lineages of J. sigillata in southwestern China. The lineages (subpopulations) diverge during the last glacial period (~1.34 Ma). Southwestern China was a glacial refuge during the last glacial period, but increasingly colder and arid climates might have fostered the fragmentation of J. regia subsp. sigillata within this refugium. Finally, we found that recent habitat fragmentation has led to a reduction in population connectivity and increased genetic differentiation by genetic drift in isolated populations. Our results support a conclusion that geological and climatic factors since the Miocene triggered the differentiation, evolutionary origin, and range shifts of J. sigillata in the studied region.

Gene Introgression among Closely Related Species in Sympatric Populations: A Case Study of Three Walnut (Juglans) Species

Gene introgression usually results from natural hybridization occurring among closely related species in sympatric populations. In this study, we discussed two rare and frequent gene flow phenomena between three species of Juglans plants and analyzed the possible causes for the difference. We collected 656 individuals from 40 populations of Persian walnut (Juglans regia L.), Chinese walnut (J. cathayensis Dode), and Iron walnut (J. sigillata Dode) that were genotyped at 17 expressed sequence tag simple sequence repeat (EST-SSR) loci to analyze the introgressions between J. regia and J. cathayensis, and J. regia and J. sigillata. Our study compared the spatial patterns of expected heterozygosity (HE), allelic richness (Rs), and private allele richness (PAR) so as to vividly infer the biogeographic history of related species of Juglans in the two regions. The results of the PCoA, UPGMA, and STRUCTURE analyses showed that all J. regia and J. sigillata populations clustered into one group, and the J. cathayensis populations clustered into the other group. The results of the historical gene flow analysis indicated that J. regia and J. sigillata have no genetic barriers, and the directional gene flow is mainly from J. regia to J. sigillata. For the three species of Juglans, all the above results indicated that gene flow was common among the same group of Juglans, and only rare and low-level gene flow appeared in distinct groups. Therefore, our study revealed multiple phenomena of gene flow and introgression among closely related species in sympatric populations, thereby providing a theoretical basis for the genetic evolution of the genus Juglans.

Antioxidative Effects and Mechanism Study of Bioactive Peptides from Defatted Walnut ( Juglans regia L.) Meal Hydrolysate

The peptide components of defatted walnut ( Juglans regia L.) meal hydrolysate (DWMH) remain unclear, hindering the investigation of biological mechanisms and exploitation of bioactive peptides. The present study aims to identify the peptide composition of DWMH, followed by to evaluate in vitro antioxidant effects of selected peptides and investigate mechanisms of antioxidative effect. First, more than 1 000 peptides were identified by de novo sequencing in DWMH. Subsequently, a scoring method was established to select promising bioactive peptides by structure based screening. Eight brand new peptides were selected due to their highest scores in two different batches of DWMH. All of them showed potent in vitro antioxidant effects on H₂O₂-injured nerve cells. Four of them even possessed significantly stronger effects than DWMH, making the selected bioactive peptides useful for further research as new bioactive entities. Two mechanisms of hydroxyl radical scavenging and ROS reduction were involved in their antioxidative effects at different degrees. The results showed peptides possessing similar capacity of hydroxyl radical scavenging or ROS reduction may have significantly different in vitro antioxidative effects. Therefore, comprehensive consideration of different antioxidative mechanisms were suggested in selecting antioxidative peptides from DWMH.

Population Structure, Genetic Diversity, and Gene Introgression of Two Closely Related Walnuts (Juglans regia and J. sigillata) in Southwestern China Revealed by EST-SSR Markers

The common walnut (Juglans regia L.) and iron walnut (J. sigillata Dode) are well-known economically important species cultivated for their edible nuts, high-quality wood, and medicinal properties and display a sympatric distribution in southwestern China. However, detailed research on the genetic diversity and introgression of these two closely related walnut species, especially in southwestern China, are lacking. In this study, we analyzed a total of 506 individuals from 28 populations of J. regia and J. sigillata using 25 EST-SSR markers to determine if their gene introgression was related to sympatric distribution. In addition, we compared the genetic diversity estimates between them. Our results indicated that all J. regia populations possess slightly higher genetic diversity than J. sigillata populations. The Geostatistical IDW technique (HO, PPL, NA and PrA) revealed that northern Yunnan and Guizhou provinces had high genetic diversity for J. regia while the northwestern Yunnan province had high genetic diversity for J. sigillata. AMOVA analysis revealed that significant genetic variation was mainly distributed within population as 73% in J. regia and 76% in J. sigillata. The genetic differentiation (FST) was 0.307 between the two walnut species (p < 0.0001), which was higher than FST values within populations (J. regia FST = 0.265 and J. sigillata FST = 0.236). However, the STRUCTURE analysis of the J. regia and J. sigillata populations revealed two genetic clusters in which gene introgression exists, therefore, the boundary of separation between these two walnut species is not clear. Moreover, these results were validated by NJ and UPGMA analysis with additional conformation from the PCoA. Based on the SSR data, our results indicate that J. sigillata is an ecotype of J. regia. Taken together, these results reveal novel information on population genetics and provide specific geographical regions containing high genetic diversity of the Juglans species sampled, which will assist in future conservation management.

The Phytogeographic History of Common Walnut in China

Common walnut (Juglans regia L.) is an economically important hardwood tree species cultivated worldwide for its high quality wood and edible nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of common walnut in China is a matter of debate, however. We estimated the genetic diversity and spatial genetic structure of 31 walnut populations sampled across its Chinese range using 22 microsatellite markers (13 neutral and 9 non-neutral). Using historical data and population genetic analysis, including approximate Bayesian analysis (ABC), we reconstructed the demographic history of J. regia in China. The genetic data indicated the likely presence of J. regia in glacial refugia in the Xinjiang province (Northwest China), Northeastern China (Beijing, Shandong, and Changbai Mountains), Central China (Qinling and Baishan Mountains and Xi'an), and Southwestern China (Tibet, Yunnan, Guizhou, and Sichuan provinces). Based on DIY-ABC analysis, we identified three ancient lineages of J. regia in China. Two lineages (subpopulation A and subpopulation B+C) diverged about 2.79 Mya, while Southwestern China, and Qinling and Baishan Mountains lineages diverged during the Quaternary glaciations (about 1.13 Mya). Remnants of these once-distinct genetic clusters of J. regia may warrant ecological management if they are to be retained as in situ resources. A population size expansion in Northeastern China was detected in the last five centuries. The present distribution of walnut in China resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and later human exploitation.

Completion of the Chloroplast Genomes of Five Chinese Juglans and Their Contribution to Chloroplast Phylogeny

Juglans L. (walnuts and butternuts) is an economically and ecologically important genus in the family Juglandaceae. All Juglans are important nut and timber trees. Juglans regia (Common walnut), J. sigillata (Iron walnut), J. cathayensis (Chinese walnut), J. hopeiensis (Ma walnut), and J. mandshurica (Manchurian walnut) are native to or naturalized in China. A strongly supported phylogeny of these five species is not available due to a lack of informative molecular markers. We compared complete chloroplast genomes and determined the phylogenetic relationships among the five Chinese Juglans using IIumina sequencing. The plastid genomes ranged from 159,714 to 160,367 bp encoding 128 functional genes, including 88 protein-coding genes and 40 tRNA genes each. A complete map of the variability across the genomes of the five Juglans species was produced that included single nucleotide variants, indels (insertions and deletions), and large structural variants, as well as differences in simple sequence repeats (SSR) and repeat sequences. Molecular phylogeny strongly supported division of the five walnut species into two previously recognized sections (Juglans/Dioscaryon and Cardiocaryon) with a 100% bootstrap (BS) value using the complete cp genomes, protein coding sequences (CDS), and the introns and spacers (IGS) data. The availability of these genomes will provide genetic information for identifying species and hybrids, taxonomy, phylogeny, and evolution in Juglans, and also provide insight into utilization of Juglans plants.