Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus jujuba Assessed with Genome-wide SNPs Developed by Genotyping-by-sequencing and SSR Markers

Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.)

BACKGROUND

Teak (Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak (Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2).

METHODS AND RESULTS

The CAD1, MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m-3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0-21% SNP pairs) at P < 0.05 and high LD decay (33-934 bp) at R2 = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly (P < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06-7.04 × 10-4) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait.

CONCLUSION

Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

Genetic diversity, population genetic structure and gene flow in the rare and endangered wild plant Cypripedium macranthos revealed by genotyping-by-sequencing

BACKGROUND

Genetic diversity, genetic structure, and gene flow in plant populations and their influencing factors are important in conservation biology. Cypripedium macranthos is one of the few wild orchids with high ornamental value in northern China. However, over the past decade, excessive collection, trading, tourism development, habitat fragmentation, deceptive pollination, and seed germination difficulties have all caused a sharp decline in the number of C. macranthos individuals and its population. In order to propose a scientific and effective conservation strategy, the genetic diversity, genetic structure and gene flow of the current CM population are urgent scientific issues to be clarified.

RESULTS

Here, 99 individuals of C. macranthos from north and northeast China were analyzed to evaluate the genetic diversity, gene flow among populations, and genetic structure by genotyping-by-sequencing. More than 68.44 Gb high-quality clean reads and 41,154 SNPs were obtained. Our data based on bioinformatics methods revealed that C. macranthos has lower genetic diversity, high levels of historical gene flow, and moderate-to-high genetic differentiation between populations. The gene migration model revealed that the direction of gene flow was mainly from northeast populations to north populations in China. The results of genetic structure analysis showed that 11 C. macranthos populations can be considered as two groups, and further divided into four subgroups. Moreover, the Mantel test detected no significant "Isolation by Distance" between populations.

CONCLUSIONS

Our study demonstrates that the present genetic diversity and genetic structure of C. macranthos populations were mainly caused by biological characteristics, human interference, habitat fragmentation, and restricted gene flow. Finally, constructive measures, which can provide a basis for the proposal of conservation strategies, have been suggested.

Exploring genetic diversity and population structure of a large grapevine (Vitis vinifera L.) germplasm collection in Türkiye

Grapevine (Vitis Vinifera L.) has been one of the significant perennial crops in widespread temperate climate regions since its domestication around 6000 years ago. Grapevine and its products, particularly wine, table grapes, and raisins, have significant economic importance not only in grapevine-growing countries but also worldwide. Grapevine cultivation in Türkiye dates back to ancient times, and Anatolia is considered one of the main grapevine migration routes around the Mediterranean basin. Turkish germplasm collection, conserved at the Turkish Viticulture Research Institutes, includes cultivars and wild relatives mainly collected in Türkiye, breeding lines, rootstock varieties, and mutants, but also cultivars of international origin. Genotyping with high-throughput markers enables the investigation of genetic diversity, population structure, and linkage disequilibrium, which are crucial for applying genomic-assisted breeding. Here, we present the results of a high-throughput genotyping-by-sequencing (GBS) study of 341 genotypes from grapevine germplasm collection at Manisa Viticulture Research Institute. A total of 272,962 high-quality single nucleotide polymorphisms (SNP) markers on the nineteen chromosomes were identified using genotyping-by-sequencing (GBS) technology. The high-density coverage of SNPs resulted in an average of 14,366 markers per chromosome, an average polymorphism information content (PIC) value of 0.23 and an expected heterozygosity (He) value of 0.28 indicating the genetic diversity within 341 genotypes. LD decayed very fast when r² was between 0.45 and 0.2 and became flat when r² was 0.05. The average LD decay for the entire genome was 30 kb when r^2 = 0.2. The PCA and structure analysis did not distinguish the grapevine genotypes based on different origins, highlighting the occurrence of gene flow and a high amount of admixture. Analysis of molecular variance (AMOVA) results indicated a high level of genetic differentiation within populations, while variation among populations was extremely low. This study provides comprehensive information on the genetic diversity and population structure of Turkish grapevine genotypes.

Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing

BACKGROUND

Studying the genetic characteristics of tea plant (Camellia spp.) leaf traits is essential for improving yield and quality through breeding and selection. Guizhou Plateau, an important part of the original center of tea plants, has rich genetic resources. However, few studies have explored the associations between tea plant leaf traits and single nucleotide polymorphism (SNP) markers in Guizhou.

RESULTS

In this study, we used the genotyping-by-sequencing (GBS) method to identify 100,829 SNP markers from 338 accessions of tea germplasm in Guizhou Plateau, a region with rich genetic resources. We assessed population structure based on high-quality SNPs, constructed phylogenetic relationships, and performed genome-wide association studies (GWASs). Four inferred pure groups (G-I, G-II, G-III, and G-IV) and one inferred admixture group (G-V), were identified by a population structure analysis, and verified by principal component analyses and phylogenetic analyses. Through GWAS, we identified six candidate genes associated with four leaf traits, including mature leaf size, texture, color and shape. Specifically, two candidate genes, located on chromosomes 1 and 9, were significantly associated with mature leaf size, while two genes, located on chromosomes 8 and 11, were significantly associated with mature leaf texture. Additionally, two candidate genes, located on chromosomes 1 and 2 were identified as being associated with mature leaf color and mature leaf shape, respectively. We verified the expression level of two candidate genes was verified using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and designed a derived cleaved amplified polymorphism (dCAPS) marker that co-segregated with mature leaf size, which could be used for marker-assisted selection (MAS) breeding in Camellia sinensis.

CONCLUSIONS

In the present study, by using GWAS approaches with the 338 tea accessions population in Guizhou, we revealed a list of SNPs markers and candidate genes that were significantly associated with four leaf traits. This work provides theoretical and practical basis for the genetic breeding of related traits in tea plant leaves.

Development of Nuclear DNA Markers for Applications in Genetic Diversity Study of Oil Palm-Pollinating Weevil Populations

The oil palm-pollinating weevil (Elaeidobius kamerunicus Faust) was introduced from Cameroon, West Africa, to Malaysia in 1981, and subsequently, to other oil palm-growing countries as well. This study aims to develop a set of robust E. kamerunicus-specific nuclear DNA markers to directly assess the genetic diversity of the weevil populations. A total of 19,148 SNP and 223,200 SSR were discovered from 48 weevils representing three origins (Peninsular Malaysia, Sabah, and Riau) using RAD tag sequencing. Subsequent filtering steps further reduced these to 1000 SNP and 120 SSR. The selected 220 SNP exhibited a polymorphism information content (PIC) of 0.2387 (±0.1280), and 8 SSR had the PIC of 0.5084 (±0.1928). These markers were found to show sufficient polymorphism, making it possible to assign 180 weevils into three major clusters from Ghana, Cameroon, and Southeast Asia (mainly in Malaysia and Indonesia). These DNA markers successfully confirmed the Cameroon origin of the Southeast Asian cluster. However, the presence of null alleles in the SSR markers, due to limited flexibility of the probe design on the short RAD tags, led to an underestimation of heterozygosity within the populations. Hence, the developed SNP markers turned out to be more efficient than the SSR markers in the genetic diversity assessment of the E. kamerunicus populations. The genetic information provides useful insight into developing guidelines for the genetic monitoring and conservation planning of E. kamerunicus.

Assessment of genetic diversity and population structure in wild Ziziphus species from northwest India using SSR marker technique

BACKGROUND

Ziziphus species particularly Ziziphus mauritiana and Ziziphus nummularia constitute an important part of genetic resources in India. They contribute economically as a fruit crop with lots of morphological and pomological variability. In current study, 48 accessions belonging to two wild Ziziphus species, i.e., Z. mauritiana and Z. nummularia, were characterized using SSR markers. In addition, external features were also examined using stereomicroscope.

RESULTS

Present investigation was done to explore the genetic structure of North Indian jujube. In total, 23 SSR markers detected 57 SSR alleles with an average of 2.47 alleles. Highest number of alleles (4) were detected by three primers, namely BFU1178, BFU479, and ZCMS14, while lowest number of alleles (2) were detected by fifteen primers. Highest Polymorphism Information Content (PIC) was 0.500 and shown by two primers, namely BFU528 and BFU1248, while lowest PIC (0.041) was observed in primers BFU286 with mean value of 0.443. Similarly, highest value of marker index (MI) was detected by primer BFU1178 i.e. 1.969, and lowest value of marker index was observed in primer BFU286 i.e. 0.021. Dendrogram generated using SSR markers data and principal component analysis showed two major groups of the analyzed germplasm with intermixing. STRUCTURE analysis also clustered all the accessions into two groups. We did not found correlation between geographic and genetic distances.

CONCLUSIONS

The preliminary results suggest that there is high level of gene pool mixing in these species which can be attributed to their cross-pollination habit. However, more such studies with large numbers of samples are required in future to gain concrete insights of the genetic structure in these species.

Genetic Diversity and Population Structure of Schima superba From Southern China

The tree Schima superba is important for afforestation and fire prevention in southern China. The wood of this tree can also be used for furniture and buildings. However, the lack of genetic background and genomic information for this species has lowered wood yield speed and quality improvement. Here, we aimed to discover genome-wide single nucleotide polymorphisms (SNPs) in 302 S. superba germplasms collected from southern China and to use these SNPs to investigate the population structure. Using genotyping by sequencing, a total of 785 high-quality SNP markers (minor allele frequency [MAF] ≥ 0.05) were identified from 302 accessions collected from seven geographical locations. Population structure analyses and principal coordinate analyses (PCoAs) indicated that these germplasm resources can be clearly separated into different populations. The S. superba accessions originating from Yunnan (YN) and Guangxi (GX) fell into the same population, separate from the accessions originating from Guangdong (GD), which indicated that these two regions should be regarded as major provenances of this species. In addition, two independent core germplasm sets with abundant genetic polymorphisms were constructed to support the breeding work. The identification of SNP markers, analyses of population genetics, and construction of core germplasm sets will greatly promote the molecular breeding work of S. superba.

Genome-wide core sets of SNP markers and Fluidigm assays for rapid and effective genotypic identification of Korean cultivars of lettuce (Lactuca sativa L.)

Lettuce is one of the economically important leaf vegetables and is cultivated mainly in temperate climate areas. Cultivar identification based on the distinctness, uniformity, and stability (DUS) test is a prerequisite for new cultivar registration. However, DUS testing based on morphological features is time-consuming, labor-intensive, and costly, and can also be influenced by environmental factors. Thus, molecular markers have also been used for the identification of genetic diversity as an effective, accurate, and stable method. Currently, genome-wide single nucleotide polymorphisms (SNPs) using next-generation sequencing technology are commonly applied in genetic research on diverse plant species. This study aimed to establish an effective and high-throughput cultivar identification system for lettuce using core sets of SNP markers developed by genotyping by sequencing (GBS). GBS identified 17 877 high-quality SNPs for 90 commercial lettuce cultivars. Genetic differentiation analyses based on the selected SNPs classified the lettuce cultivars into three main groups. Core sets of 192, 96, 48, and 24 markers were further selected and validated using the Fluidigm platform. Phylogenetic analyses based on all core sets of SNPs successfully discriminated individual cultivars that have been currently recognized. These core sets of SNP markers will support the construction of a DNA database of lettuce that can be useful for cultivar identification and purity testing, as well as DUS testing in the plant variety protection system. Additionally, this work will facilitate genetic research to improve breeding in lettuce.

Identification of a Diverse Core Set Panel of Rice From the East Coast Region of India Using SNP Markers

In India, rice (Oryza sativa L.) is cultivated under a variety of climatic conditions. Due to the fragility of the coastal ecosystem, rice farming in these areas has lagged behind. Salinity coupled with floods has added to this trend. Hence, to prevent genetic erosion, conserving and characterizing the coastal rice, is the need of the hour. This work accessed the genetic variation and population structure among 2,242 rice accessions originating from India’s east coast comprising Andhra Pradesh, Orissa, and Tamil Nadu, using 36 SNP markers, and have generated a core set (247 accessions) as well as a mini-core set (30 accessions) of rice germplasm. All the 36 SNP loci were biallelic and 72 alleles found with average two alleles per locus. The genetic relatedness of the total collection was inferred using the un-rooted neighbor-joining tree, which grouped all the genotypes (2,242) into three major clusters. Two groups were obtained with a core set and three groups obtained with a mini core set. The mean PIC value of total collection was 0.24, and those of the core collection and mini core collection were 0.27 and 0.32, respectively. The mean heterozygosity and gene diversity of the overall collection were 0.07 and 0.29, respectively, and the core set and mini core set revealed 0.12 and 0.34, 0.20 and 0.40 values, respectively, representing 99% of distinctiveness in the core and mini core sets. Population structure analysis showed maximum population at K = 4 for total collection and core collection. Accessions were distributed according to their population structure confirmed by PCoA and AMOVA analysis. The identified small and diverse core set panel will be useful in allele mining for biotic and abiotic traits and managing the genetic diversity of the coastal rice collection. Validation of the 36-plex SNP assay was done by comparing the genetic diversity parameters across two different rice core collections, i.e., east coast and northeast rice collection. The same set of SNP markers was found very effective in deciphering diversity at different genetic parameters in both the collections; hence, these marker sets can be utilized for core development and diversity analysis studies.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

Genome-Wide Scanning Enabled SNP Discovery, Linkage Disequilibrium Patterns and Population Structure in a Panel of Fonio (Digitaria exilis [Kippist] Stapf) Germplasm

White fonio (Digitaria exilis) is a staple food for millions of people in arid and semi-arid areas of West Africa. Knowledge about nutritional and health benefits, insights into morphological diversity, and the recent development of genomic resources call for a better understanding of the genetic structure of the extant germplasm gathered throughout the region in order to set up a robust breeding program. We assessed the genetic diversity and population structure of 259 fonio individuals collected from six countries from West Africa (Nigeria, Benin, Guinea, Mali, Burkina Faso and Niger) in this study using 688 putative out of 21,324 DArTseq-derived SNP markers. Due to the inbreeding and small population size, the results revealed a substantial level of genetic variability. Furthermore, two clusters were found irrespective of the geographic origins of accessions. Moreover, the high level of linkage disequilibrium (LD) between loci observed resulted from the mating system of the crop, which is often associated with a low recombination rate. These findings fill the gaps about the molecular diversity and genetic structure of the white fonio germplasm in West Africa. This was required for the application of genomic tools that can potentially speed up the genetic gain in fonio millet breeding for complex traits such as yield, and other nutrient contents.

Genome-Wide Assessment of Population Structure and Genetic Diversity of the Global Finger Millet Germplasm Panel Conserved at the ICRISAT Genebank

Finger millet [Eleusine coracana (L.) Gaertn.] is an important climate-resilient nutrient-dense crop grown as a staple food grain in Asia and Africa. Utilizing the full potential of the crop mainly depends on an in-depth exploration of the vast diversity in its germplasm. In this study, the global finger millet germplasm diversity panel of 314 accessions was genotyped, using the DArTseq approach to assess genetic diversity and population structure. We obtained 33,884 high-quality single nucleotide polymorphism (SNP) markers on 306 accessions after filtering. Finger millet germplasm showed considerable genetic diversity, and the mean polymorphic information content, gene diversity, and Shannon Index were 0.110, 0.114, and 0.194, respectively. The average genetic distance of the entire set was 0.301 (range 0.040 - 0.450). The accessions of the race elongata (0.326) showed the highest average genetic distance, and the least was in the race plana (0.275); and higher genetic divergence was observed between elongata and vulgaris (0.320), while the least was between compacta and plana (0.281). An average, landrace accessions had higher gene diversity (0.144) and genetic distance (0.299) than the breeding lines (0.117 and 0.267, respectively). A similar average gene diversity was observed in the accessions of Asia (0.132) and Africa (0.129), but Asia had slightly higher genetic distance (0.286) than African accessions (0.276), and the distance between these two regions was 0.327. This was also confirmed by a model-based STRUCTURE analysis, genetic distance-based clustering, and principal coordinate analysis, which revealed two major populations representing Asia and Africa. Analysis of molecular variance suggests that the significant population differentiation was mainly due to within individuals between regions or between populations while races had a negligible impact on population structure. Finger millet diversity is structured based on a geographical region of origin, while the racial structure made negligible contribution to population structure. The information generated from this study can provide greater insights into the population structure and genetic diversity within and among regions and races, and an understanding of genomic-assisted finger millet improvement.

Genetic Diversity and Population Differentiation of Pinus koraiensis in China

Pinus koraiensis is a well-known precious tree species in East Asia with high economic, ornamental and ecological value. More than fifty percent of the P. koraiensis forests in the world are distributed in northeast China, a region with abundant germplasm resources. However, these natural P. koraiensis sources are in danger of genetic erosion caused by continuous climate changes, natural disturbances such as wildfire and frequent human activity. Little work has been conducted on the population genetic structure and genetic differentiation of P. koraiensis in China because of the lack of genetic information. In this study, 480 P. koraiensis individuals from 16 natural populations were sampled and genotyped. Fifteen polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers were used to evaluate genetic diversity, population structure and differentiation in P. koraiensis. Analysis of molecular variance (AMOVA) of the EST-SSR marker data showed that 33% of the total genetic variation was among populations and 67% was within populations. A high level of genetic diversity was found across the P. koraiensis populations, and the highest levels of genetic diversity were found in HH, ZH, LS and TL populations. Moreover, pairwise Fst values revealed significant genetic differentiation among populations (mean Fst = 0.177). According to the results of the STRUCTURE and Neighbor-joining (NJ) tree analyses and principal component analysis (PCA), the studied geographical populations cluster into two genetic clusters: cluster 1 from Xiaoxinganling Mountains and cluster 2 from Changbaishan Mountains. These results are consistent with the geographical distributions of the populations. The results provide new genetic information for future genome-wide association studies (GWAS), marker-assisted selection (MAS) and genomic selection (GS) in natural P. koraiensis breeding programs and can aid the development of conservation and management strategies for this valuable conifer species.

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.

Development of core-collections for Guizhou tea genetic resources and GWAS of leaf size using SNP developed by genotyping-by-sequencing

An accurate depiction of the genetic relationship, the development of core collection, and genome-wide association analysis (GWAS) are key for the effective exploitation and utilization of genetic resources. Here, genotyping-by-sequencing (GBS) was used to characterize 415 tea accessions mostly collected from the Guizhou region in China. A total of 30,282 high-quality SNPs was used to estimate the genetic relationships, develop core collections, and perform GWAS. We suggest 198 and 148 accessions to represent the core set and mini-core set, which consist of 47% and 37% of the whole collection, respectively, and contain 93–95% of the total SNPs. Furthermore, the frequencies of all alleles and genotypes in the whole set were very well retained in the core set and mini-core set. The 415 accessions were clustered into 14 groups and the core and the mini-core collections contain accessions from each group, species, cultivation status and growth habit. By analyzing the significant SNP markers associated with multiple traits, nine SNPs were found to be significantly associated with four leaf size traits, namely MLL, MLW, MLA and MLSI (P < 1.655E−06). This study characterized the genetic distance and relationship of tea collections, suggested the core collections, and established an efficient GWAS analysis of GBS result.

High-Density Genetic Map Construction and QTL Mapping of Leaf and Needling Traits in Ziziphus jujuba Mill

The Chinese jujube (Ziziphus jujuba Mill., 2n = 2x = 24), one of the most popular fruit trees in Asia, is widely cultivated and utilized in China, where it is traditionally consumed as both a fresh and dried food resource. A high-density genetic map can provide the necessary framework for quantitative trait loci (QTL) analyses and map-based gene cloning and molecular breeding. In this study, we constructed a new high-density genetic linkage map via a genotyping-by-sequencing approach. For the consensus linkage map, a total of 3,792 markers spanning 2,167.5 cM were mapped onto 12 linkage groups, with an average marker interval distance of 0.358 cM. The genetic map anchored 301 Mb (85.7%) of scaffolds from the sequenced Z. jujuba "Junzao" genome. Based on this genetic map, 30 potential QTLs were detected, including 27 QTLs for leaf traits and 3 QTLs for needling length. This high-density genetic map and the identified QTLs for relevant agronomic traits lay the groundwork for functional genetic mapping, map-based cloning, and marker-assisted selection in jujube.

Characterization of genetic diversity and population structure in wheat using array based SNP markers

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder's Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~ 50%), A (~ 39%), and D (~ 10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon's information index (I) = 0.648, expected heterozygosity (He) = 0.456 and unbiased expected heterozygosity (uHe) = 0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

Comparative transcriptome and metabolite profiling of four tissues from Alisma orientale (Sam.) Juzep reveals its inflorescence developmental and medicinal characteristics

Alisma orientale (Sam.) Juzep (A. orientale) is an important medicinal plant in traditional Chinese medicine. In this study, de novo RNA-seq of A. orientale was performed based on the cDNA libraries from four different tissues, roots, leaves, scapes and inflorescences. A total of 41,685 unigenes were assembled, 25,024 unigene functional annotations were obtained by searching against the five public sequence databases, and 3,411 simple sequence repeats in A. orientale were reported for the first time. 15,402 differentially expressed genes were analysed. The morphological characteristics showed that compared to the other tissues, the leaves had more chlorophyll, the scapes had more vascular bundles, and the inflorescences contained more starch granules and protein. In addition, the metabolic profiles of eight kinds of alisols metabolite profiling, which were measured by ultra-Performance liquid chromatography-triple quadrupole-mass spectrometry showed that alisol B 23-acetate and alisol B were the major components of the four tissues at amounts of 0.068~0.350 mg/g and 0.046~0.587 mg/g, respectively. In addition, qRT-PCR validated that farnesyl pyrophosphate synthase and 3-hydroxy-3-methylglutaryl-CoA reductase should be considered the critical candidate genes involved in alisol biosynthesis. These transcriptome and metabolic profiles of A. orientale may help clarify the molecular mechanisms underlying the medicinal characteristics of A. orientale.

Genetic diversity, linkage disequilibrium, and population structure analysis of the tea plant (Camellia sinensis) from an origin center, Guizhou plateau, using genome-wide SNPs developed by genotyping-by-sequencing

BACKGROUND

To efficiently protect and exploit germplasm resources for marker development and breeding purposes, we must accurately depict the features of the tea populations. This study focuses on the Camellia sinensis (C. sinensis) population and aims to (i) identify single nucleotide polymorphisms (SNPs) on the genome level, (ii) investigate the genetic diversity and population structure, and (iii) characterize the linkage disequilibrium (LD) pattern to facilitate next genome-wide association mapping and marker-assisted selection.

RESULTS

We collected 415 tea accessions from the Origin Center and analyzed the genetic diversity, population structure and LD pattern using the genotyping-by-sequencing (GBS) approach. A total of 79,016 high-quality SNPs were identified; the polymorphism information content (PIC) and genetic diversity (GD) based on these SNPs showed a higher level of genetic diversity in cultivated type than in wild type. The 415 accessions were clustered into three groups by STRUCTURE software and confirmed using principal component analyses (PCA)-wild type, cultivated type, and admixed wild type. However, unweighted pair group method with arithmetic mean (UPGMA) trees indicated the accessions should be grouped into more clusters. Further analyses identified four groups, the Pure Wild Type, Admixed Wild Type, ancient landraces and modern landraces using STRUCTURE, and the results were confirmed by PCA and UPGMA tree method. A higher level of genetic diversity was detected in ancient landraces and Admixed Wild Type than that in the Pure Wild Type and modern landraces. The highest differentiation was between the Pure Wild Type and modern landraces. A relatively fast LD decay with a short range (kb) was observed, and the LD decays of four inferred populations were different.

CONCLUSIONS

This study is, to our knowledge, the first population genetic analysis of tea germplasm from the Origin Center, Guizhou Plateau, using GBS. The LD pattern, population structure and genetic differentiation of the tea population revealed by our study will benefit further genetic studies, germplasm protection, and breeding.

An Integrated Characterization of Jujube (Ziziphus jujuba Mill.) Grown in the North Adriatic Region

Jujube (Ziziphus jujuba Mill.) has favourable horticultural properties including adaptation to arid conditions, abiotic and biotic stresses, as well as positive impact on human health. The present study describes the characterization of genetic diversity of the germplasm of jujube from the Istrian peninsula, the determination of important chemical compounds, antioxidative properties in relation to antibacterial and antifungal activities of jujube fruit extracts, and the determination of nutritional properties of jujube fruit. The results of the genetic analysis showed that most of the samples from the Istrian peninsula belong to two recently introduced varieties, 'Li' and 'Lang', and the most widespread local variety 'Navadna žižola'. The local variety has smaller fruit than the ‘Li’ and ‘Lang’ varieties, with thick and fleshy mesocarp. Chemical analysis indicated that fruits of the local variety contained a valuable source of dietary fibre ((9.7±0.6) g/100 g) and were rich in minerals such as (in g/100 g dry mass): potassium (829±51), calcium (177±11) and phosphorus (129±19). Aqueous extracts showed slight antibacterial activity, while ethanol extracts had higher mass fractions of phenolic compounds (expressed as gallic acid equivalents (GAE), 5.8-8.7 mg/g) than the aqueous extracts, but did not show antimicrobial activity. Compounds other than phenolic compounds in jujube fruit may be more biologically active. Based on the results of these analyses, the local Istrian jujube variety is a promising candidate for cultivation potential.

Genetic Diversity and Population Structure of a Camelina sativa Spring Panel

There is a need to explore renewable alternatives (e.g., biofuels) that can produce energy sources to help reduce the reliance on fossil oils. In addition, the consumption of fossil oils adversely affects the environment and human health via the generation of waste water, greenhouse gases, and waste solids. Camelina sativa, originated from southeastern Europe and southwestern Asia, is being re-embraced as an industrial oilseed crop due to its high seed oil content (36-47%) and high unsaturated fatty acid composition (>90%), which are suitable for jet fuel, biodiesel, high-value lubricants and animal feed. C. sativa's agronomic advantages include short time to maturation, low water and nutrient requirements, adaptability to adverse environmental conditions and resistance to common pests and pathogens. These characteristics make it an ideal crop for sustainable agricultural systems and regions of marginal land. However, the lack of genetic and genomic resources has slowed the enhancement of this emerging oilseed crop and exploration of its full agronomic and breeding potential. Here, a core of 213 spring C. sativa accessions was collected and genotyped. The genotypic data was used to characterize genetic diversity and population structure to infer how natural selection and plant breeding may have affected the formation and differentiation within the C. sativa natural populations, and how the genetic diversity of this species can be used in future breeding efforts. A total of 6,192 high-quality single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing (GBS) technology. The average polymorphism information content (PIC) value of 0.29 indicate moderate genetic diversity for the C. sativa spring panel evaluated in this report. Population structure and principal coordinates analyses (PCoA) based on SNPs revealed two distinct subpopulations. Sub-population 1 (POP1) contains accessions that mainly originated from Germany while the majority of POP2 accessions (>75%) were collected from Eastern Europe. Analysis of molecular variance (AMOVA) identified 4% variance among and 96% variance within subpopulations, indicating a high gene exchange (or low genetic differentiation) between the two subpopulations. These findings provide important information for future allele/gene identification using genome-wide association studies (GWAS) and marker-assisted selection (MAS) to enhance genetic gain in C. sativa breeding programs.