De novo assembly and characterization of the leaf, bud, and fruit transcriptome from the vulnerable tree Juglans mandshurica for the development of 20 new microsatellite markers using Illumina sequencing

Genome-wide identification analysis of the 4-Coumarate: CoA ligase (4CL) gene family expression profiles in Juglans regia and its wild relatives J. Mandshurica resistance and salt stress

Persian walnut (Juglans regia) and Manchurian walnut (Juglans mandshurica) belong to Juglandaceae, which are vulnerable, temperate deciduous perennial trees with high economical, ecological, and industrial values. 4-Coumarate: CoA ligase (4CL) plays an essential function in plant development, growth, and stress. Walnut production is challenged by diverse stresses, such as salinity, drought, and diseases. However, the characteristics and expression levels of 4CL gene family in Juglans species resistance and under salt stress are unknown. Here, we identified 36 Jr4CL genes and 31 Jm4CL genes, respectively. Based on phylogenetic relationship analysis, all 4CL genes were divided into three branches. WGD was the major duplication mode for 4CLs in two Juglans species. The phylogenic and collinearity analyses showed that the 4CLs were relatively conserved during evolution, but the gene structures varied widely. 4CLs promoter region contained multiply cis-acting elements related to phytohormones and stress responses. We found that Jr4CLs may be participated in the regulation of resistance to anthracnose. The expression level and some physiological of 4CLs were changed significantly after salt treatment. According to qRT-PCR results, positive regulation was found to be the main mode of regulation of 4CL genes after salt stress. Overall, J. mandshurica outperformed J. regia. Therefore, J. mandshurica can be used as a walnut rootstock to improve salt tolerance. Our results provide new understanding the potential functions of 4CL genes in stress tolerance, offer the theoretical genetic basis of walnut varieties adapted to salt stress, and provide an important reference for breeding cultivated walnuts for stress tolerance.

A Multiplex PCR System of Novel Microsatellite Loci for Population Genetic Application in Walnuts

Multiplex polymerase chain reaction (PCR) of microsatellite loci allows for simultaneous amplification of two or more pairs of primers in a single PCR reaction; hence, it is cost and time effective. However, very few attempts have been reported in non-model species. In this study, by combining a genome-based de novo development and cross-species application approach, a multiplex PCR system comprising 5 PCR reactions of 33 microsatellites consisting of 26 novel genomic and 7 literature-sourced loci was tested for polymorphisms, cross-species transferability, and the ability to assess genetic diversity and population structure of three walnut species (Juglans spp.). We found that the genome-based approach is more efficient than other methods. An allelic ladder was developed for each locus to enhance consistent genotyping among laboratories. The population genetic analysis results showed that all 33 loci were successfully transferred across the three species, showing high polymorphism and a strong genetic structure. Hence, the multiplex PCR system is highly applicable in walnut species. Furthermore, we propose an efficient pipeline to characterize and genotype polymorphic microsatellite loci. The novel toolbox developed here will aid future ecology and evolution studies in walnut and could serve as a model for other plant species.

An integrated approach of homogenization-ultrasound-assisted extraction followed by high-speed countercurrent chromatography for rapidly preparative separation naphthoquinones and diarylheptanes from Juglands mandshurica Maxim. exocarp

This research investigated the effectiveness of an integrated method for the extraction and separation of naphthoquinones and diarylheptanes from exocarp of Juglands mandshurica Maxim. (namely, green walnut husks). The target compounds were obtained by ultra-turrax homogenization (UTH) coupled with ultrasound-assisted extraction (UAE) technology followed by high-speed countercurrent chromatography (HSCCC). The UTH-UAE extraction method achieved higher efficiency with 2.49- and 2.36-fold to those by UAE, and 1.39- and 1.34-fold to those by UTH in a short time. HSCCC was adopted for further separation and purification; six target compounds, namely, regiolone (RE), juglone (JU), myricatomento-genin (MG), galleon (GA), 2-oxatrycyclo[13.2.2.13,7]eicosa-3,5,7(20),15,17,18-hexaen-10-16-diol (OE), and juglanin A (JA), were separated with more than 95.37% purities and more than 84.71% final recovery rates, respectively. In this study, the integrated strategy of extraction and separation could get high purity compounds quickly, which would provide time and solvent saved method for the natural products separation from plants.

De Novo Transcriptome Profiling for the Generation and Validation of Microsatellite Markers, Transcription Factors, and Database Development for Andrographis paniculata

Andrographis paniculata belongs to the family Acanthaceae and is known for its medicinal properties owing to the presence of unique constituents belonging to the lactones, diterpenoids, diterpene glycosides, flavonoids, and flavonoid glycosides groups of chemicals. Andrographolide, a major therapeutic constituent of A. paniculata, is extracted primarily from the leaves of this plant and exhibits antimicrobial and anti-inflammatory activities. Using 454 GS-FLX pyrosequencing, we have generated a whole transcriptome profile of entire leaves of A. paniculata. A total of 22,402 high-quality transcripts were generated, with an average transcript length and N50 of 884 bp and 1007 bp, respectively. Functional annotation revealed that 19,264 (86%) of the total transcripts showed significant similarity with the NCBI-Nr database and were successfully annotated. Out of the 19,264 BLAST hits, 17,623 transcripts were assigned GO terms and distributed into three major functional categories: molecular function (44.62%), biological processes (29.19%), and cellular component (26.18%) based on BLAST2GO. Transcription factor analysis showed 6669 transcripts, belonging to 57 different transcription factor families. Fifteen TF genes that belong to the NAC, MYB, and bHLH TF categories were validated by RT PCR amplification. In silico analysis of gene families involved in the synthesis of biochemical compounds having medicinal values, such as cytochrome p450, protein kinases, heat shock proteins, and transporters, was completed and a total of 102 different transcripts encoding enzymes involved in the biosynthesis of terpenoids were predicted. Out of these, 33 transcripts belonged to terpenoid backbone biosynthesis. This study also identified 4254 EST-SSRs from 3661 transcripts, representing 16.34% of the total transcripts. Fifty-three novel EST-SSR markers generated from our EST dataset were used to assess the genetic diversity among eighteen A. paniculata accessions. The genetic diversity analysis revealed two distinct sub-clusters and all accessions based on the genetic similarity index were distinct from each other. A database based on EST transcripts, EST-SSR markers, and transcription factors has been developed using data generated from the present study combined with available transcriptomic resources from a public database using Meta transcriptome analysis to make genomic resources available in one place to the researchers working on this medicinal plant.

Genome-wide identification, characterization, and expression pattern of the late embryogenesis abundant (LEA) gene family in Juglans regia and its wild relatives J. mandshurica

BACKGROUND

Late Embryogenesis Abundant (LEA) proteins are a class of proteins associated with plant stress resistance. Two Juglans species, Juglans regia and J. mandshurica, are both diploid (2n = 32), monoecious perennial economic tree species with high edible, pharmaceutical, and timber value. The identification, characterization, and expression patterns of LEA proteins in J. regia and its wild relative, J. mandshurica, would not only provide the genetic basis of this gene family, but it would also supply clues for further studies of the evolution and regulating mechanisms of LEA proteins in other tree species.

RESULTS

In this study, we identified 25 and 20 members of the LEA gene family in Juglans regia and its wild relative, Juglans mandshurica, respectively. The results of phylogenetic analysis showed that the LEA members were divided into eight main subgroups. Predictions of their physicochemical properties showed the variable characteristics of LEA proteins, and the subcellular localization analysis indicated that most LEA proteins are localized in the nucleus. Chromosomal localization analysis and gene replication pattern prediction indicated that WGD is the predominant duplication mode of LEA genes. The results of the comparative analysis indicated a high level of collinearity between the two Juglans species. Analysis of cis-acting elements indicated that LEA genes had a relatively wide range of responses to abiotic stresses and phytohormonal processes, particularly in two phytohormones, methyl jasmonate and abscisic acid. Transcriptome profiling and qRT-PCR experiments showed that JrLEAs are commonly expressed in leaves, green husks, and male and female flowers, and most JmLEAs are more highly expressed in male flowers. We also hypothesized that JrLEAs are involved in the process of anthracnose resistance. Anthracnose-resistant varieties of JrLEAs presented relatively high expression levels at later stages.

CONCLUSION

In this study, we provide a theoretical basis for the functional study of LEA genes in J. regia and J. mandshurica. Analysis of cis-acting elements and gene expression indicated that JrLEAs and JmLEAs play important roles in resistance to biotic stresses in these species.

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.

Genetic Diversity and Structure of Persian Walnut (Juglans regia L.) in Pakistan: Implications for Conservation

Persian (Common) walnut (Juglans regia L.) is a famous fruit tree species valued for its nutritious nuts and high-quality wood. Although walnut is widely distributed and plays an important role in the economy and culture of Pakistan, the genetic diversity and structure of its populations in the country remains poorly understood. Therefore, using 31 nuclear microsatellites, we assessed the genetic diversity and population structure of 12 walnut populations sampled across Pakistan. We also implemented the geostatistical IDW technique in ArcGIS to reveal “hotspots” of genetic diversity. Generally, the studied populations registered relatively low indices of genetic diversity (N_A = 3.839, H_O = 0.558, U_HE = 0.580), and eight populations had positive inbreeding coefficient (F_IS) values. Low among-population differentiation was indicated by AMOVA, pairwise F_ST and D_C. STRUCTURE, PCoA and neighbor joining (NJ) analysis revealed a general lack of clear clustering in the populations except that one population in Upper Dir was clearly genetically distinct from the rest. Furthermore, the Mantel test showed no correlation between the geographic and genetic distance (r = 0.14, p = 0.22), while barrier analysis suggested three statistically significant genetic barriers. Finally, the spatial interpolation results indicated that populations in Ziarat, Kashmir, Dir, Swat, Chitral, and upper Dir had high intrapopulation genetic diversity, suggesting the need to conserve populations in those areas. The results from this study will be important for future breeding improvement and conservation of walnuts in Pakistan.

Improved de novo chromosome-level genome assembly of the vulnerable walnut tree Juglans mandshurica reveals gene family evolution and possible genome basis of resistance to lesion nematode

Manchurian walnut (Juglans mandshurica Maxim.) is a synonym of J. cathayensis, a diploid, vulnerable, temperate deciduous tree valued for its wood and nut. It is also valued as a rootstock for Juglans regia because of its reported tolerance of lesion nematode. Reference genomes are available for several Juglans species, our goal was to produce a de novo, chromosome-level assembly of the J. mandshurica genome. Here, we reported an improved assembly of J. mandshurica with a contig N50 size of 6.49 Mb and a scaffold N50 size of 36.1 Mb. The total genome size was 548 Mb encoding 29,032 protein coding genes which were annotated. The collinearity analysis showed that J. mandshurica and J. regia originated from a common ancestor, with both species undergoing two WGD events. A genomic comparison showed that J. mandshurica was missing 1657 genes found in J. regia, and J. mandshurica includes 2827 genes not found in of the J. regia genome. The J. mandshurica contained 1440 unique paralogues that were highly enriched for flavonoid biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction. Four gene families related to disease resistance notable contraction (rapidly evolving; LEA, WAK, PPR, and PR) in J. mandshurica compared to eight species. JmaPR10 and JmaPR8 contained three orthologous gene pairs with J. regia that were highly expressed in root bark. JmaPR10 is a strong candidate gene for lesion nematodes resistance in J. mandshurica. The J. mandshurica genome should be a useful resource for study of the evolution, breeding, and genetic variation in walnuts (Juglans).

Flower Development of Heterodichogamous Juglans mandshurica (Juglandaceae)

Juglans mandshurica is a monoecious heterodichogamous species with protogynous and protandrous mating strategies that occur at a 1:1 ratio and are randomly distributed in the population. The inconsistent male and female flowering periods of the same mating type result in an imbalance of the ratio of male and female flowers, contributing to the low yield of this species. However, little more is known about its floral development. Following three consecutive years of observations, histological analysis, and scanning electron microscopy, we found that the morphological and anatomical development of the male and female flowers were synchronous. The male floral morphological development of J. mandshurica was divided into seven phases, while that of the female flower was nine. Four stages were shared between the male and female flower's anatomical development. Our findings indicate that there was minimal overlap between sexual functions within the same mating type, guaranteeing synchronization, mutual non-interference, outcrossing, and avoidance of self-fertilization. These results provide a theoretical basis for the improvement of fruit yield and quality through the reasonable allocation of protogynous and protandrous individuals in a population, and for artificial pollination control. Further, these findings lay a foundation for further research on the genetic mechanisms and environmental effects on flower development of heterodichogamous J. mandshurica.

Juglans mandshurica Maxim.: A Review of Its Traditional Usages, Phytochemical Constituents, and Pharmacological Properties

Juglans mandshurica Maxim., also known as “Manchurian walnut” (Chinese) and “Onigurumi” (Japanese), is a medicinal plant widely distributed in Western and Central Asia, especially in China. It has been traditionally used to treat cancer, gastric ulcers, diarrhea, dysentery, dermatosis, uterine prolapse, and leukopenia. To date, more than 400 constituents including quinones (e.g. naphthoquinones, anthraquinones, naphthalenones, tetralones), phenolics, flavonoids, triterpenoids, coumarins, lignans, phenylpropanoids, diarylheptanoids, and steroids, were isolated and structurally identified from different plant parts of J. mandshurica. Among them, quinones, phenolics, triterpenoids, and diarylheptanoids, as the major bioactive substances, have been extensively studied and displayed significant bioactivity. Previous studies have demonstrated that J. mandshurica and a few of its active components exhibit a wide range of pharmacologically important properties, such as antitumor, immunomodulatory, anti-inflammatory, neuroprotective, anti-diabetic, antiviral, antimicrobial, and anti-melanogenesis activities. However, many investigations on biological activities were mainly based on crude extracts of this plant, and the major bioactive ingredients responsible for these bioactivities have not been well identified. Further in vitro and in vivo studies on the mechanisms of action of the pure bioactive compounds, and more elaborate toxicity studies as well as clinical studies are needed to ensure safety and effectiveness of the plant for human use. Taken together, the present review will provide some specific useful suggestions guide to further investigations and applications of this plant in the preparation of medicines and functional foods.

De novo assembly and Transcriptome characterization of an endemic species of Vietnam, Panax vietnamensis Ha et Grushv., including the development of EST-SSR markers for population genetics

BACKGROUND

Understanding the genetic diversity in endangered species that occur inforest remnants is necessary to establish efficient strategies for the species conservation, restoration and management. Panax vietnamensis Ha et Grushv. is medicinally important, endemic and endangered species of Vietnam. However, genetic diversity and structure of population are unknown due to lack of efficient molecular markers.

RESULTS

In this study, we employed Illumina HiSeq™ 4000 sequencing to analyze the transcriptomes of P. vietnamensis (roots, leaves and stems). Raw reads total of 23,741,783 was obtained and then assembled, from which the generated unigenes were 89,271 (average length = 598.3191 nt). The 31,686 unigenes were annotated in different databases i.e. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Nucleotide Collection (NR/NT) and Swiss-Prot for functional annotation. Further, 11,343 EST-SSRs were detected. From 7774 primer pairs, 101 were selected for polymorphism validation, in which; 20 primer pairs were successfully amplified to DNA fragments and significant amounts of polymorphism was observed within population. The nine polymorphic microsatellite loci were used for population structure and diversity analyses. The obtained results revealed high levels of genetic diversity in populations, the average observed and expected heterozygosity were HO = 0.422 and HE = 0.479, respectively. During the Bottleneck analysis using TPM and SMM models (p < 0.01) shows that targeted population is significantly heterozygote deficient. This suggests sign of the bottleneck in all populations. Genetic differentiation between populations was moderate (FST = 0.133) and indicating slightly high level of gene flow (Nm = 1.63). Analysis of molecular variance (AMOVA) showed 63.17% of variation within individuals and 12.45% among populations. Our results shows two genetic clusters related to geographical distances.

CONCLUSION

Our study will assist conservators in future conservation management, breeding, production and habitats restoration of the species.

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.

Development and Characterization of EST-SSR Markers From RNA-Seq Data in Phyllostachys violascens

Bamboo are woody grass species containing important economic and ecological values. Lei bamboo (Phyllostachys violascens) is a kind of shoot-producing bamboo species with the highest economic yield per unit area. However, identifying different varieties of Lei bamboo based on morphological characteristics is difficult. Microsatellites play an important role in plant identification and genetic diversity analysis and are superior to other molecular markers. In this study, we identified 18,356 expressed sequence tag-simple sequence repeat (EST-SSR) loci in Lei bamboo transcriptome data. A total of 11,264 primer pairs were successfully designed from unigenes of all EST-SSR loci, and 96 primer pairs were randomly selected and synthesized. A total of 54 primer pairs were used for classifying 16 Lei bamboo varieties and 10 different Phyllostachys species. The number of polymorphism alleles among the 54 primer pairs ranged from 3 to 12 for P. violascens varieties and 3 to 20 for Phyllostachys. The phylogenetic tree based on polymorphism alleles successfully distinguished 16 P. violascens varieties and 10 Phyllostachys species. Our study provides abundant EST-SSR resources that are useful for genetic diversity analysis and molecular verification of bamboo and suggests that SSR markers developed from Lei bamboo are more efficient and reliable than ISSR, SRAP or AFLP markers.

De Novo Transcriptomic Analysis and Development of EST–SSRs for Styrax japonicus

Styrax japonicus sieb. et Zucc. is widely distributed in China with ornamental and medicinal values. However, the transcriptome of S. japonicus has not yet been reported. In this study, we carried out the first transcriptome analysis of S. japonicus and developed a set of expressed sequence tag–simple sequence repeats (EST–SSRs). We obtained 338,570,222 clean reads in total, of which the mean GC content was 41.58%. In total, 136,071 unigenes were obtained having an average length of 611 bp and 71,226 unigenes were favorably annotated in the database. In total, we identified 55,977 potential EST–SSRs from 38,611 unigenes, of which there was 1 SSR per 6.73 kb. The di-nucleotide repeats (40.40%) were the most identified SSRs. One set of 60 primer pairs was randomly selected, and the amplified products in S. japonicus were validated; 28 primer pairs successfully produced clear amplicons. A total of 21 (35%) polymorphic genic SSR markers were identified between two populations. In total, 15 alleles were detected and the average number was 6. The average of observed heterozygosity and expected heterozygosity was 0.614 and 0.552, respectively. The polymorphism information content (PIC) value fluctuated between 0.074 and 0.855, with a mean value of 0.504, which was also the middle level. This study provides useful information for diversity studies and resource assessments of S. japonicus.

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.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

A new coumarin from Juglans mandshurica Maxim induce apoptosis in hepatocarcinoma cells

In this study, a new coumarin, juglansoside C (1) was isolated from the bark of Juglans mandshurica. Its chemical structure was identified by comprehensive spectroscopic analyses. The in vitro cytotoxicity assay showed that 1 exhibited moderate cytotoxicity against human hepatocellular carcinoma Hep3B cells with an IC₅₀ value of 70.9 μM. Furthermore, Annexin V-FITC/PI staining assay indicated that 1 markedly induced apoptosis in Hep3B cells.

Characterization of the Lycium barbarum fruit transcriptome and development of EST-SSR markers

Lycium barbarum, commonly known as goji, is important in Chinese herbal medicine and its fruit is a very important agricultural and biological product. However, the molecular mechanism of formation of its fruit and associated medicinal and nutritional components is unexplored. Moreover, this species lacks SSR markers due to lack of genomic and transcriptomic information. In this study, a total of 139,333 unigenes with average length of 1049 bp and N50 of 1579 bp are obtained by trinity assembly from Illumina sequencing reads. A total of 92,498 (66.38%) unigenes showed similarities in at least one database including Nr (46.15%), Nt (56.56%), KO (15.56%), Swiss-prot (33.34%), Pfam (33.43%), GO (33.62%) and KOG/COG (17.55%). Genes in flavonoid and taurine biosynthesis pathways were found and validated by RT-qPCR. A total of 50,093 EST-SSRs were identified from 38,922 unigenes, and 22,537 EST-SSR primer pairs were designed. Four hundred pairs of SSR markers were randomly selected to validate assembly quality, of which 352 (88%) were successful in PCR amplification of genomic DNA from 11 Lycium accessions and 210 produced polymorphisms. The polymorphic loci showed that the genetic similarity of the 11 Lycium accessions ranged from 0.50 to 0.99 and the accessions could be divided into 4 groups. These results will facilitate investigations of the molecular mechanism of formation of L. barbarum fruit and associated medicinal and nutritional components, and will be of value to novel gene discovery and functional genomic studies. The EST-SSR markers will be useful for genetic diversity evaluation, genetic mapping and marker-assisted breeding.

A rapid and cost-effective approach for the development of polymorphic microsatellites in non-model species using paired-end RAD sequencing

As one of the most informative and versatile DNA-based markers, microsatellites have been widely used in population and conservation genetic studies. However, the development of microsatellites has traditionally been laborious, time-consuming, and expensive. In the present study, a rapid and cost-effective "RAD-seq-Assembly-Microsatellite" approach was developed to identify abundant microsatellite markers in non-model species using the roughskin sculpin Trachidermus fasciatus as a representative. Overlapping paired-end Illumina reads generated by restriction-site-associated DNA sequencing (RAD-seq) were clustered based on the similarity of reads containing the restriction enzyme recognition site and then assembled into contigs, which were used for microsatellite discovery and primer design. A total of 121,750 RAD contigs were generated with a mean length of 522 bp, and 19,782 contigs contained microsatellite motifs. A total of 156,150 primer pairs were successfully designed based on 16,497 contigs containing priming sites. Experimental validation of 52 randomly selected microsatellite loci demonstrated that 45 (86.54%) loci were successfully amplified and polymorphic in two geographically isolated populations of T. fasciatus. Compared with traditional approaches based on DNA cloning and other approaches based on next-generation sequencing, our newly developed approach could yield thousands of microsatellite loci with much higher successful amplification rate and lower costs, especially for non-model species with shallow background of genomic information. The "RAD-seq-Assembly-Microsatellite" approach holds great promise for microsatellite development in future ecological and evolutionary studies of non-model species.

De Novo Assembly of Transcriptome and Development of Novel EST-SSR Markers in Rhododendron rex Lévl. through Illumina Sequencing

Transcriptome sequences generated by next-generation sequencing (NGS) technologies can be utilized to rapidly detect and characterize a large number of gene-based microsatellites from different plants. Rhododendron rex Lévl. is a perennial woody species from the family Ericaceae and an endangered plant with high ornamental value endemic to Southwestern China. Nevertheless, the genetic and genomic information of R. rex remain unknown. In this study, we performed transcriptome sequencing for R. rex leaf samples, and generated large transcript sequences for functional characterization and development gene-associated SSR markers. A total of 164,242 unigenes were assembled and 115,089 (70.07%) unigenes were successfully annotated in public databases. In addition, a total of 15,314 potential EST-SSRs were identified, and the frequency of SSRs in the R. rex unigenes was 9.32%, with an average of one EST-SSR per 5.65 kb. The most abundant type was repeated di-nucleotide (54.63%), followed by mono- (26.03%) and tri-nucleotide (18.51%) repeats. Based on the SSR-containing sequence, 100 primer pairs were randomly selected and synthesized and used for assessment of the polymorphism. Thirty-six primer pairs were polymorphic and revealed polymorphism among 20 individuals from four R. rex populations. A total of 197 alleles were identified, with an average of 5.472 alleles per locus. The Polymorphism Information Content ranged from 0.154 to 0.870, with a mean of 0.482. The newly developed EST-SSR markers exhibited high transferability (58.33-83.33%) among the six subgenera. Thus, these novel EST-SSR markers developed would provide valuable sequence resources for population structure, genetic diversity analysis, and genetic resource assessments of R. rex and its related species.

Transcriptome Analysis to Identify the Putative Biosynthesis and Transport Genes Associated with the Medicinal Components of Achyranthes bidentata Bl

Achyranthes bidentata is a popular perennial medicine herb used for 1000s of years in China to treat various diseases. Although this herb has multiple pharmaceutical purposes in China, no transcriptomic information has been reported for this species. In addition, the understanding of several key pathways and enzymes involved in the biosynthesis of oleanolic acid and ecdysterone, two pharmacologically active classes of metabolites and major chemical constituents of A. bidentata root extracts, is limited. The aim of the present study was to characterize the transcriptome profile of the roots and leaves of A. bidentata to uncover the biosynthetic and transport mechanisms of the active components. In this study, we identified 100,987 transcripts, with an average length of 1146.8 base pairs. A total of 31,634 (31.33%) unigenes were annotated, and 12,762 unigenes were mapped to 303 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. Moreover, we identified a total of 260 oleanolic acid and ecdysterone genes encoding biosynthetic enzymes. Furthermore, the key enzymes involved in the oleanolic acid and ecdysterone synthesis pathways were analyzed using quantitative real-time polymerase chain reaction, revealing that the roots expressed these enzymes to a greater extent than the leaves. In addition, we identified 85 ATP-binding cassette transporters, some of which might be involved in the translocation of secondary metabolites.

Transcriptome Analysis and Identification of Genes Associated with Floral Transition and Flower Development in Sugar Apple (Annona squamosa L.)

Sugar apple (Annona squamosa L.) is a semi-deciduous subtropical tree that progressively sheds its leaves in the spring. However, little information is available on the mechanism involved in flower developmental pattern. To gain a global perspective on the floral transition and flower development of sugar apple, cDNA libraries were prepared independently from inflorescent meristem and three flowering stages. Illumina sequencing generated 107,197,488 high quality reads that were assembled into 71,948 unigenes, with an average sequence length of 825.40 bp. Among the unigenes, various transcription factor families involved in floral transition and flower development were elucidated. Furthermore, a Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that unigenes exhibiting differential expressions were involved in various phytohormone signal transduction events and circadian rhythms. In addition, 147 unigenes exhibiting sequence similarities to known flowering-related genes from other plants were differentially expressed during flower development. The expression patterns of 20 selected genes were validated using quantitative-PCR. The expression data presented in our study is the most comprehensive dataset available for sugar apple so far and will serve as a resource for investigating the genetics of the flowering process in sugar apple and other Annona species.