De Novo Assembly and Characterization of Bud, Leaf and Flowers Transcriptome from Juglans Regia L. for the Identification and Characterization of New EST-SSRs

Genetic Identification of Hybrid Walnuts (Juglans × intermedia Carr.) in Hungary, the Hidden Potential for Future Breeding

The question of the hybrid walnut (Juglans × intermedia Carr.) is still under debate in the Central European region. There is not simply just an underutilization, rather, even the existence of these hybrid forms is not broadly accepted. On the contrary, there is an intensive cross-breeding activity in the western part of Europe resulting in commercially available hybrid clones. In Hungary, several individuals have been reported with intermediate morphology from different old black walnut plantations. Due to the lacking information, a preliminary study was conducted in order to prove the difference of these selected trees and to test the hybrid state. For this purpose, DNA fingerprinting was used by applying 13 simple sequence repeat (SSR) markers for the identification of 22 hybrid genotypes selected from one study plot. A comparative analysis with a reference sample set of the ‘parental’ species and other known hybrids was performed as well. The genetic analysis resulted distinct, unique genotypes for all of the samples. Based on the genetic pattern, the analyzed hybrid group was clearly distinguishable from the other two walnut groups. The result of this study also highlights the hidden potential in walnut breeding in the Central European region. Future concepts concerning hybrid walnut utilization in plantation forestry, agroforestry or as breeding material are also discussed.

Inflorescence Transcriptome Sequencing and Development of New EST-SSR Markers in Common Buckwheat (Fagopyrum esculentum)

Common buckwheat (Fagopyrum esculentum M.) is known for its adaptability, good nutrition, and medicinal and health care value. However, genetic studies of buckwheat have been hindered by limited genomic resources and genetic markers. In this study, Illumina HiSeq 4000 high-throughput sequencing technology was used to sequence the transcriptome of green-flower common buckwheat (Gr) with coarse pedicels and white-flower Ukrainian daliqiao (UD) with fine pedicels. A total of 118,448 unigenes were obtained, with an average length of 1248 bp and an N50 of 1850 bp. A total of 39,432 differentially expressed genes (DEGs) were identified, and the DEGs of the porphyrins and chlorophyll metabolic pathway had significantly upregulated expression in Gr. Then, a total of 17,579 sequences containing SSR loci were detected, and 20,756 EST-SSR loci were found. The distribution frequency of EST-SSR in the transcriptome was 17.52%, and the average distribution density was 8.21 kb. A total of 224 pairs of primers were randomly selected for synthesis; 35 varieties of common buckwheat and 13 varieties of Tartary buckwheat were verified through these primers. The clustering results well verified the previous conclusion that common buckwheat and Tartary buckwheat had a distant genetic relationship. The EST-SSR markers identified and developed in this study will be helpful to enrich the transcriptome information and marker-assisted selection breeding of buckwheat.

Genetic Diversity among Some Walnut (Juglans regia L.) Genotypes by SSR Markers

The food needs for increasing population, climatic changes, urbanization and industrialization, along with the destruction of forests, are the main challenges of modern life. Therefore, it is very important to evaluate plant genetic resources in order to cope with these problems. Therefore, in this study, a set of ninety-one walnut (Juglans regia L.) accessions from Central Anatolia region, composed of seventy-four accessions and eight commercial cultivars from Turkey, and nine international reference cultivars, was analyzed using 45 SSR (Simple Sequence Repeats) markers to reveal the genetic diversity. SSR analysis identified 390 alleles for 91 accessions. The number of alleles per locus ranged from 3 to 19 alleles with a mean value of 9 alleles per locus. Genetic dissimilarity coefficients ranged from 0.03 to 0.68. The highest number of alleles was obtained from CUJRA212 locus (Na = 19). The values of polymorphism information content (PIC) ranged from 0.42 (JRHR222528) to 0.86 (CUJRA212) with a mean PIC value of 0.68. Genetic distances were estimated according to the UPGMA (Unweighted Pair Group Method with Arithmetic Average), Principal Coordinates (PCoA), and the Structure-based clustering. The UPGMA and Structure clustering of the accessions depicted five major clusters supporting the PCoA results. The dendrogram revealed the similarities and dissimilarities among the accessions by identifying five major clusters. Based on this study, SSR analyses indicate that Yozgat province has an important genetic diversity pool and rich genetic variance of walnuts.

Is "Esterhazy II", an Old Walnut Variety in the Hungarian Gene Bank, the Original Genotype?

The old walnut (Juglans regia L.) genotype called "Esterhazy II" was well-known in the Austro-Hungarian Monarchy before World War II, and it can still be found in the Austrian, German and Swiss backyard gardens today. Unfortunately, nowadays, vegetatively propagated progenies of the original "Esterhazy II" are not available anymore around the world because walnut grafting started later than this genotype had become well-known. Although various accessions with "Esterhazy II"-"blood" are available, it is difficult to determine which one can be considered true or the most similar to the original one. In this paper, phenological and nut morphological characteristics of an "Esterhazy II" specimen planted in a Hungarian gene bank were compared to the varieties "Milotai 10" and "Chandler". Examined characteristics were: budbreak, blossom time, type of dichogamy, ripening time, nut and kernel features. An additional SSR fingerprinting was used to identify identical genotypes and to demonstrate the relatedness of the analyzed "Esterhazy II" genotype to the other Hungarian walnut cultivars. It can be concluded that under the name "Esterhazy II", several different genotypes can be observed. All the checked characteristics except budbreak fitted well with the previous descriptions. Our results confirmed that the examined "Esterhazy II" genotype shows high similarity to the "original" "Esterhazy II" described in the literature.

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).

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.

Genetic Diversity and Population Genetic Structure of Cinnamomum camphora in South China Revealed by EST-SSR Markers

Cinnamomum camphora is a valuable broad-leaf tree indigenous to South China and East Asia and has been widely cultivated and utilized by humans since ancient times. However, owing to its overutilization for essential oil extraction, the Transplanting Big Trees into Cities Program, and over deforestation to make furniture, its wild populations have been detrimentally affected and are declining rapidly. In the present study, the genetic diversity and population structure of 180 trees sampled from 41 populations in South China were investigated with 22 expressed sequence tag-simple sequence repeat (EST-SSR) markers. In total, 61 alleles were harbored across 180 individuals, and medium genetic diversity level was inferred from the observed heterozygosity (Ho), expected heterozygosity (He), and Nei’ gene diversity (GD), which were 0.45, 0.44, and 0.44, respectively. Among the 41 wild populations, C. camphora had an average of 44 alleles, 2.02 effective alleles, and He ranging from 0.30 (SC) to 0.61 (HK). Analysis of molecular variance (AMOVA) showed that 17% of the variation among populations and the average pairwise genetic differentiation coefficient (FST) between populations was 0.162, indicating relatively low genetic population differentiations. Structure analysis suggested two groups for the 180 individuals, which was consistent with the principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA). Populations grouped to cluster I were nearly all distributed in Jiangxi Province (except population XS in Zhejiang Province), and cluster II mainly comprised populations from other regions, indicating a significant geographical distribution. Moreover, the Mantel test showed that this geographical distance was significantly correlated with genetic distance. The findings of this research will assist in future C. camphora conservation management and breeding programs.

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.

Analysis of genetic diversity and structure in a worldwide walnut (Juglans regia L.) germplasm using SSR markers

Persian or English walnut (Juglans regia L.), the walnut species cultivated for nut production, is one of the oldest food sources known and is grown worldwide in temperate areas. France is the 7^th leading producer as of 2016 with 39 kt. Deciphering walnut genetic diversity and structure is important for efficient management and use of genetic resources. In this work, 253 worldwide accessions from the INRA walnut germplasm collection, containing English walnut and several related species, were genotyped using 13 SSR (Single Sequence Repeat) markers selected from the literature to assess diversity and structure. Genetic diversity parameters showed a deficiency of heterozygotes and, for several SSRs, allele-specificities among the accessions tested. Principal Coordinate Analysis (PCoA) showed the 253 accessions clustered in largely in agreement with the existing botanical classification of the genus. Among the 217 J. regia accessions, two main clusters, accessions from Eastern Europe and Asia, and accessions from Western Europe and America, were identified using STRUCTURE software. This was confirmed by Principal Coordinate Analysis and supported by Neighbor-Joining tree construction using DARwin software. Moreover, a substructure was found within the two clusters, mainly according to geographical origin. A core collection containing 50 accessions was selected using the maximum length sub-tree method and prior knowledge about their phenotype. The present study constitutes a preliminary population genetics overview of INRA walnut genetic resources collection using SSR markers. The resulting estimations of genetic diversity and structure are useful for germplasm management and for future walnut breeding programs.

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.

Genomic Variation Among and Within Six Juglans Species

Genomic analysis in Juglans (walnuts) is expected to transform the breeding and agricultural production of both nuts and lumber. To that end, we report here the determination of reference sequences for six additional relatives of Juglans regia: Juglans sigillata (also from section Dioscaryon), Juglans nigra, Juglans microcarpa, Juglans hindsii (from section Rhysocaryon), Juglans cathayensis (from section Cardiocaryon), and the closely related Pterocarya stenoptera. While these are ‘draft’ genomes, ranging in size between 640Mbp and 990Mbp, their contiguities and accuracies can support powerful annotations of genomic variation that are often the foundation of new avenues of research and breeding. We annotated nucleotide divergence and synteny by creating complete pairwise alignments of each reference genome to the remaining six. In addition, we have re-sequenced a sample of accessions from four Juglans species (including regia). The variation discovered in these surveys comprises a critical resource for experimentation and breeding, as well as a solid complementary annotation. To demonstrate the potential of these resources the structural and sequence variation in and around the polyphenol oxidase loci, PPO1 and PPO2 were investigated. As reported for other seed crops variation in this gene is implicated in the domestication of walnuts. The apparently Juglandaceae specific PPO1 duplicate shows accelerated divergence and an excess of amino acid replacement on the lineage leading to accessions of the domesticated nut crop species, Juglans regia and sigillata.