Genetic Variation and Population Structure of Oryza glaberrima and Development of a Mini-Core Collection Using DArTseq

A comparative study of population structure and genetic diversity of commercial and indigenous chickens from different agro-ecological zones in Ghana using SilicoDArT and SNP markers

Chicken production, both in the local and commercial sectors, contributes significantly to human livelihood and food security. Precise use of diverse genetic resources is primary in breeding programs. The study analyzed the genetic diversity and population structure of commercial chickens and indigenous chicken ecotypes from three different agro-ecological zones (Semi-Deciduous Rainforest Zone, Guinea Savannah, and Coastal Savannah) using SilicoDArT and SNP markers, utilizing whole-genome sequencing and phenotypic data. Phenotypic data were collected from 72 indigenous chicken ecotypes across the three AEZs, and 32 commercial birds kept at the Kwame Nkrumah University of Science and Technology (KNUST). DNA samples used for sequencing were obtained from 88 chickens (62 indigenous chicken ecotypes and 26 commercial chickens). A total of 54,995 SilicoDArT and 85,396 SNPs markers were generated from DArTseq genotyping. After filtering, 44,784 SilicoDArT and 58,353 SNP were used for genetic diversity and population structure analysis. Both markers showed high reproducibility and call rate. Polymorphic information content (PIC) values ranged from 0.00 to 0.50, while ≥ 50 % showed PIC values more than the median. Furthermore, we obtained FST values, Nei's genetic distance, dendrogram analysis, and principal component analysis (PCA) of commercial and indigenous chickens. The FST and Nei's genetic distance showed that there is high genetic diversity between the commercial chickens and the indigenous chicken ecotypes. However, there was low genetic diversity among the indigenous chicken ecotypes. The PCA analysis indicated a clear separation between the commercial and indigenous chicken ecotypes, while no clear separation was observed between the indigenous chicken ecotypes. The phenotypic data and the dendrogram indicated that naked and frizzle genes do not markedly alter the genetics of indigenous and commercial birds, and their influence on economic traits may be solely determined by the prevailing environmental conditions. The results indicate that there is high genetic differentiation between commercial and indigenous chickens based on SilicoDArT and SNP markers. The indigenous chickens from the agro-ecological zones have low genetic diversity and might have a common origin. Naked neck and frizzle genes do not markedly alter the genetic performance of birds in terms of economic traits. Therefore, the superiority of birds carrying these genes in economic traits may be solely due to environmental variation.

Phenotypic diversity and population structure of Pecan (Carya illinoinensis) collections reveals geographic patterns

Pecan (Carya illinoinensis) is an economically important nut crop known for its genetic diversity and adaptability to various climates. Understanding the growth variability, phenological traits, and population structure of pecan populations is crucial for breeding programs and conservation. In this study, plant growth and phenological traits were evaluated over three consecutive seasons (2015-2017) for 550 genotypes from 26 provenances. Significant variations in plant height, stem diameter, and budbreak were observed among provenances, with Southern provenances exhibiting faster growth and earlier budbreak compared to Northern provenances. Population structure analysis using SNP markers revealed eight distinct subpopulations, reflecting genetic differentiation among provenances. Notably, Southern Mexico collections formed two separate clusters, while Western collections, such as 'Allen 3', 'Allen 4', and 'Riverside', were distinguished from others. 'Burkett' and 'Apache' were grouped together due to their shared maternal parentage. Principal component analysis and phylogenetic tree analysis further supported subpopulation differentiation. Genetic differentiation among the 26 populations was evident, with six clusters highly in agreement with the subpopulations identified by STRUCTURE and fastSTRUCTURE. Principal components analysis (PCA) revealed distinct groups, corresponding to subpopulations identified by genetic analysis. Discriminant analysis of PCA (DAPC) based on provenance origin further supported the genetic structure, with clear separation of provenances into distinct clusters. These findings provide valuable insights into the genetic diversity and growth patterns of pecan populations. Understanding the genetic basis of phenological traits and population structure is essential for selecting superior cultivars adapted to diverse environments. The identified subpopulations can guide breeding efforts to develop resilient rootstocks and contribute to the sustainable management of pecan genetic resources. Overall, this study enhances our understanding of pecan genetic diversity and informs conservation and breeding strategies for the long-term viability of pecan cultivation.

Phylogeography of Morella nana: The Wumeng Mountains as a natural geographical isolation boundary on the Yunnan-Guizhou Plateau

The Yunnan-Guizhou Plateau (YGP) is characterized by the distinctive isolated habitat of the limestone Karst Islands and features the Wumeng Mountains, which divide the YGP into the two Plateaus of Yunnan and Guizhou. This study aimed to assess the effects of geographic isolation and past climate fluctuation on the distribution of flora in the YGP. To achieve this, we carried out the phylogeographical pattern and genetic structure based on chloroplast and nuclear ribosomal DNA sequence in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modeling for Morella nana, an important wild plant resource and endemic to the YGP once considered a vulnerable species. The results suggest that the genetic and chlorotype network structures of M. nana are divided into at least two groups: cpDNA chlorotype H2 (or dominant nrDNA haplotypes h1 and h2), distributed primarily to the east of the Wumeng Mountains, and cpDNA chlorotypes H1 and H3-H10 (or dominant nrDNA haplotype h2 and h3), distributed to the west of the Wumeng Mountains. A deep genetic split was noted within the two groups to reach 25 steps, especially for the cpDNA fragment variation. This east-west divergence reveals the existence of a natural geographical isolation boundary in the form of the Wumeng Mountains, and supports the existence of at least two glacial refuges during the Quaternary glacial period, along with two genetic diversity center, and at least two large geographic protection units for the important species of M. nana. This study indicates that the phylogeographical pattern of M. nana can be attributed to geographic/environmental isolation caused by the Wumeng Mountains and climate fluctuation during the last glacial maximum, and proposes an effective strategy to protecting this important plant resource.

Analysis of population structure and genetic diversity of Camellia tachangensis in Guizhou based on SNP markers

BACKGROUND

Camellia tachangensis F. C. Zhang is a five-compartment species in the ovary of tea group plants, which represents the original germline of early differentiation of some tea group plants.

METHODS AND RESULTS

In this study, we analyzed single-nucleotide polymorphisms (SNPs) at the genome level, constructed a phylogenetic tree, analyzed the genetic diversity, and further investigated the population structure of 100 C. tachangensis accessions using the genotyping-by-sequencing (GBS) method. A total of 91,959 high-quality SNPs were obtained. Population structure analysis showed that the 100 C. tachangensis accessions clustered into three groups: YQ-1 (Village Group), YQ-2 (Forest Group) and YQ-3 (Transition Group), which was further consistent with the results of phylogenetic analysis and principal component analyses (PCA). In addition, a comparative analysis of the genetic diversity among the three populations (Forest, Village, and Transition Groups) detected the highest genetic diversity in the Transition Group and the highest differentiation between Forest and Village Groups.

CONCLUSIONS

C. tachangensis plants growing in the forest had different genetic backgrounds from those growing in villages. This study provides a basis for the effective protection and utilization of C. tachangensis populations and lays a foundation for future C. tachangensis breeding.

Genome-wide association study, population structure, and genetic diversity of the tea plant in Guizhou Plateau

BACKGROUND

Guizhou Plateau, as one of the original centers of tea plant, has a profound multi-ethnic cultural heritage and abundant tea germplasm resources. However, the impact of indigenous community factors on the genetic diversity, population structure and geographical distribution of tea plant is still unclear.

RESULTS

Using the genotyping-by-sequencing (GBS) approach, we collected 415 tea plant accessions from the study sites, estimated genetic diversity, developed a core collection, and conducted a genome-wide association study (GWAS) based on 99,363 high-quality single-nucleotide polymorphisms (SNPs). A total of 415 tea accessions were clustered into six populations (GP01, GP02, GP03, GP04, GP05 and GP06), and the results showed that GP04 and GP05 had the highest and lowest genetic diversity (Pi = 0.214 and Pi = 0.145, respectively). Moreover, 136 tea accessions (33%) were selected to construct the core set that can represent the genetic diversity of the whole collection. By analyzing seven significant SNP markers associated with the traits such as the germination period of one bud and two leaves (OTL) and the germination period of one bud and three leaves (OtL), four candidate genes possibly related to OTL and OtL were identified.

CONCLUSIONS

This study revealed the impact of indigenous communities on the population structure of 415 tea accessions, indicating the importance of cultural practices for protection and utilization of tea plant genetic resources. Four potential candidate genes associated with the OTL and OtL of tea plant were also identified, which will facilitate genetic research, germplasm conservation, and breeding.

Genetic Diversity of Phenotypic and Biochemical Traits in VIR Radish (Raphanus sativus L.) Germplasm Collection

Small radish and radish are economically important root crops that represent an integral part of a healthy human diet. The world collection of Raphanus L. root crops, maintained in the VIR genebank, includes 2810 accessions from 75 countries around the world, of which 2800 (1600 small radish, 1200 radish) belong to R. sativus species, three to R. raphanistrum, three to R. landra, and four to R. caudatus. It is necessary to systematically investigate the historical and modern gene pool of root-bearing plants of R. sativus and provide new material for breeding. The material for our research was a set of small radish and radish accessions of various ecological groups and different geographical origin, fully covering the diversity of the species. The small radish subset included 149 accessions from 37 countries, belonging to 13 types of seven varieties of European and Chinese subspecies. The radish subset included 129 accessions from 21 countries, belonging to 18 types of 11 varieties of European, Chinese, and Japanese subspecies. As a result of the evaluation of R. sativus accessions according to phenological, morphological, and biochemical analyses, a wide variation of these characteristics was revealed, which is due to the large genetic diversity of small radish and radish of various ecological and geographical origins. The investigation of the degree of variation regarding phenotypic and biochemical traits revealed adaptive stable and highly variable characteristics of R. sativus accessions. Such insights are crucial for the establishment and further use of trait collections. Trait collections facilitate germplasm use and contribute significantly to the preservation of genetic diversity of the gene pool.

Physical Mapping of QTL in Four Spring Wheat Populations under Conventional and Organic Management Systems. I. Earliness

In previous studies, we reported quantitative trait loci (QTL) associated with the heading, flowering, and maturity time in four hard red spring wheat recombinant inbred line (RIL) populations but the results are scattered in population-specific genetic maps, which is challenging to exploit efficiently in breeding. Here, we mapped and characterized QTL associated with these three earliness traits using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. Our data consisted of (i) 6526 single nucleotide polymorphisms (SNPs) and two traits evaluated at five conventionally managed environments in the 'Cutler' × 'AC Barrie' population; (ii) 3158 SNPs and two traits evaluated across three organic and seven conventional managements in the 'Attila' × 'CDC Go' population; (iii) 5731 SilicoDArT and SNP markers and the three traits evaluated at four conventional and organic management systems in the 'Peace' × 'Carberry' population; and (iv) 1058 SNPs and two traits evaluated across two conventionally and organically managed environments in the 'Peace' × 'CDC Stanley' population. Using composite interval mapping, the phenotypic data across all environments, and the IWGSC RefSeq v2.0 physical maps, we identified a total of 44 QTL associated with days to heading (11), flowering (10), and maturity (23). Fifteen of the 44 QTL were common to both conventional and organic management systems, and the remaining QTL were specific to either the conventional (21) or organic (8) management systems. Some QTL harbor known genes, including the Vrn-A1, Vrn-B1, Rht-A1, and Rht-B1 that regulate photoperiodism, flowering time, and plant height in wheat, which lays a solid basis for cloning and further characterization.

Genetics and Genomics of African Rice (Oryza glaberrima Steud) Domestication

African rice (Oryza glaberrima Steud) is one of the two independently domesticated rice species, the other one being Asian rice (Oryza sativa L.). Despite major progress being made in understanding the evolutionary and domestication history of African rice, key outstanding issues remain controversial. There appears to be an underlying difficulty in identifying the domestication centre and number of times the crop has been domesticated. Advances in genomics have provided unprecedented opportunities for understanding the genetic architecture of domestication related traits. For most of the domestication traits, the underlying genes and mutations have been identified. Comparative analysis of domestication genes between Asian and African rice has revealed that the two species went through an independent but convergent evolution process. The genetic and developmental basis of some of the domestic traits are conserved not only between Asian and African rice but also with other domesticated crop species. Analysis of genome data and its interpretation is emerging as a major challenge facing studies of domestication in African rice as key studies continue giving contradictory findings and conclusions. Insights obtained on the domestication of this species are vital for guiding crop improvement efforts.

The development of SSR markers based on RNA-sequencing and its validation between and within Carex L. species

BACKGROUND

Carex L. is one of the largest genera in the Cyperaceae family and an important vascular plant in the ecosystem. However, the genetic background of Carex is complex and the classification is not clear. In order to investigate the gene function annotation of Carex, RNA-sequencing analysis was performed. Simple sequence repeats (SSRs) were generated based on the Illumina data and then were utilized to investigate the genetic characteristics of the 79 Carex germplasms.

RESULTS

In this study, 36,403 unigenes with a total length of 41,724,615 bp were obtained and annotated based on GO, KOG, KEGG, NR databases. The results provide a theoretical basis for gene function exploration. Out of 8776 SSRs, 96 pairs of primers were randomly selected. One hundred eighty polymorphic bands were amplified with a polymorphism rate of 100% based on 42 pairs of primers with higher polymorphism levels. The average band number was 4.3 per primer, the average distance value was 0.548, and the polymorphic information content was ranged from 0.133 to 0.494. The number of observed alleles (Na), effective alleles (Ne), Nei's (1973) gene diversity (H), and the Shannon information index (I) were 2.000, 1.376, 0.243, and 0.391, respectively. NJ clustering divided into three groups and the accessions from New Zealand showed a similar genetic attribute and clustered into one group. UPGMA and PCoA analysis also revealed the same result. The analysis of molecular variance (AMOVA) revealed a superior genetic diversity within accessions than between accessions based on geographic origin cluster and NJ cluster. What's more, the fingerprints of 79 Carex species are established in this study. Different combinations of primer pairs can be used to identify multiple Carex at one time, which overcomes the difficulties of traditional identification methods.

CONCLUSIONS

The transcriptomic analysis shed new light on the function categories from the annotated genes and will facilitate future gene functional studies. The genetic characteristics analysis indicated that gene flow was extensive among 79 Carex species. These markers can be used to investigate the evolutionary history of Carex and related species, as well as to serve as a guide in future breeding projects.

Germplasm Acquisition and Distribution by CGIAR Genebanks

The international collections of plant genetic resources for food and agriculture (PGRFA) hosted by 11 CGIAR Centers are important components of the United Nations Food and Agriculture Organization's global system of conservation and use of PGRFA. They also play an important supportive role in realizing Target 2.5 of the Sustainable Development Goals. This paper analyzes CGIAR genebanks' trends in acquiring and distributing PGRFA over the last 35 years, with a particular focus on the last decade. The paper highlights a number of factors influencing the Centers' acquisition of new PGRFA to include in the international collections, including increased capacity to analyze gaps in those collections and precisely target new collecting missions, availability of financial resources, and the state of international and national access and benefit-sharing laws and phytosanitary regulations. Factors contributing to Centers' distributions of PGRFA included the extent of accession-level information, users' capacity to identify the materials they want, and policies. The genebanks' rates of both acquisition and distribution increased over the last decade. The paper ends on a cautionary note concerning the potential of unresolved tensions regarding access and benefit sharing and digital genomic sequence information to undermine international cooperation to conserve and use PGRFA.

Genome-wide genetic diversity and population structure of Garcinia kola (Heckel) in Benin using DArT-Seq technology

Garcinia kola (Heckel) is a versatile tree indigenous to West and Central Africa. All parts of the tree have value in traditional medicine. Natural populations of the species have declined over the years due to overexploitation. Assessment of genetic diversity and population structure of G. kola is important for its management and conservation. The present study investigates the genetic diversity and population structure of G. kola populations in Benin using ultra-high-throughput diversity array technology (DArT) single nucleotide polymorphism (SNP) markers. From the 102 accessions sampled, two were excluded from the final dataset owing to poor genotyping coverage. A total of 43,736 SNPs were reported, of which 12,585 were used for analyses after screening with quality control parameters including Minor allele frequency (≥ 0.05), call rate (≥ 80%), reproducibility (≥ 95%), and polymorphic information content (≥ 1%). Analysis revealed low genetic diversity with expected heterozygosity per population ranging from 0.196 to 0.228. Pairwise F-statistics (FST) revealed low levels of genetic differentiation between populations while an Analysis of molecular variance (AMOVA) indicated that the majority of variation (97.86%) was within populations. Population structure analysis through clustering and discriminant analysis on principal component revealed two admixed clusters, implying little genetic structure. However, the model-based maximum likelihood in Admixture indicated only one genetic cluster. The present study indicated low genetic diversity of G. kola, and interventions are needed to be tailored towards its conservation.

Genomic Analysis of Selected Maize Landraces from Sahel and Coastal West Africa Reveals Their Variability and Potential for Genetic Enhancement

Genetic adaptation of maize to the increasingly unpredictable climatic conditions is an essential prerequisite for achievement of food security and sustainable development goals in sub-Saharan Africa. The landraces of maize; which have not served as sources of improved germplasm; are invaluable sources of novel genetic variability crucial for achieving this objective. The overall goal of this study was to assess the genetic diversity and population structure of a maize panel of 208 accessions; comprising landrace gene pools from Burkina Faso (58), Ghana (43), and Togo (89), together with reference populations (18) from the maize improvement program of the International Institute of Tropical Agriculture (IITA). Genotyping the maize panel with 5974 DArTseq-SNP markers revealed immense genetic diversity indicated by average expected heterozygosity (0.36), observed heterozygosity (0.5), and polymorphic information content (0.29). Model-based population structure; neighbor-joining tree; discriminant analysis of principal component; and principal coordinate analyses all separated the maize panel into three major sub-populations; each capable of providing a wide range of allelic variation. Analysis of molecular variance (AMOVA) showed that 86% of the variation was within individuals; while 14% was attributable to differences among gene pools. The Burkinabe gene pool was strongly differentiated from all the others (genetic differentiation values >0.20), with no gene flow (Nm) to the reference populations (Nm = 0.98). Thus; this gene pool could be a target for novel genetic variation for maize improvement. The results of the present study confirmed the potential of this maize panel as an invaluable genetic resource for future design of association mapping studies to speed-up the introgression of this novel variation into the existing breeding pipelines.

Comparisons of sampling methods for assessing intra- and inter-accession genetic diversity in three rice species using genotyping by sequencing

To minimize the cost of sample preparation and genotyping, most genebank genomics studies in self-pollinating species are conducted on a single individual to represent an accession, which may be heterogeneous with larger than expected intra-accession genetic variation. Here, we compared various population genetics parameters among six DNA (leaf) sampling methods on 90 accessions representing a wild species (O. barthii), cultivated and landraces (O. glaberrima, O. sativa), and improved varieties derived through interspecific hybridizations. A total of 1,527 DNA samples were genotyped with 46,818 polymorphic single nucleotide polymorphisms (SNPs) using DArTseq. Various statistical analyses were performed on eleven datasets corresponding to 5 plants per accession individually and in a bulk (two sets), 10 plants individually and in a bulk (two sets), all 15 plants individually (one set), and a randomly sampled individual repeated six times (six sets). Overall, we arrived at broadly similar conclusions across 11 datasets in terms of SNP polymorphism, heterozygosity/heterogeneity, diversity indices, concordance among genetic dissimilarity matrices, population structure, and genetic differentiation; there were, however, a few discrepancies between some pairs of datasets. Detailed results of each sampling method, the concordance in their outputs, and the technical and cost implications of each method were discussed.

Identification and Fine Mapping of Pi69(t), a New Gene Conferring Broad-Spectrum Resistance Against Magnaporthe oryzae From Oryza glaberrima Steud

The discovery and deployment of new broad-spectrum resistance (R) genes from cultivated rice and its wild relatives is a strategy to broaden the genetic basis of modern rice cultivars to combat rice blast disease. Oryza glaberrima possessing many valuable traits for tolerance to biotic and abiotic stresses, is an elite gene pool for improvement of Asian cultivated rice. An introgression line IL106 derived from O. glaberrima (Acc. IRGC100137) confers complete resistance to Magnaporthe oryzae in blast nursery. Genetic analysis using 2185 BC₆F₂ progenies derived from a cross between IL106 and the recurrent parent Dianjingyou 1 showed that IL106 harbors a single dominant resistance gene against M. oryzae strain 09BSH-10-5A. This gene was preliminarily mapped on the long arm of chromosome 6 of rice in a region of ca. 0.9 cM delimited by two SSR markers (RM20650 and RM20701). In order to finely map this gene, 17,100 additional progenies were further analyzed. As a result, this gene was further narrowed down to a region flanked by two molecular markers STS69-15 and STS69-7, and co-segregated with 3 molecular markers, RM20676, STS69-21 and STS69-22 on the long arm of chromosome 6. Based on reference genome sequences, this R gene was mapped in silico in 76.1-Kb and 67.7-Kb physical intervals, and containing 4 and 3 NBS-LRR candidate genes in O. sativa cultivar Nipponbare and O. glaberrima cultivar CG14, respectively. Because no blast resistance gene was finely mapped in this physical interval before, this R gene was considered as not described yet and designated as Pi69(t), which is the first identified and finely mapped blast R gene from O. glaberrima, as far as we know. Evaluation of IL106 with 151 blast strains collected from 6 countries in Asia showed that 148 strains are avirulent on IL106, suggesting that Pi69(t) is a broad-spectrum blast R gene, and a promising resistant resource for improvement of Asian cultivated rice.

Allele mining unlocks the identification of RYMV resistance genes and alleles in African cultivated rice

BACKGROUND

Rice yellow mottle virus (RYMV) is a major rice pathogen in Africa. Three resistance genes, i.e. RYMV1, RYMV2 and RYMV3, have been previously described. RYMV1 encodes the translation initiation factor eIF(iso)4G1 and the best candidate genes for RYMV2 and RYMV3 encode a homolog of an Arabidopsis nucleoporin (CPR5) and a nucleotide-binding domain and leucine-rich repeat containing domain (NLR) protein, respectively. High resistance is very uncommon in Asian cultivated rice (Oryza sativa), with only two highly resistant accessions identified so far, but it is more frequent in African cultivated rice (Oryza glaberrima).

RESULTS

Here we report the findings of a resistance survey in a reference collection of 268 O. glaberrima accessions. A total of 40 resistant accessions were found, thus confirming the high frequency of resistance to RYMV in this species. We analysed the variability of resistance genes or candidate genes in this collection based on high-depth Illumina data or Sanger sequencing. Alleles previously shown to be associated with resistance were observed in 31 resistant accessions but not in any susceptible ones. Five original alleles with a frameshift or untimely stop codon in the candidate gene for RYMV2 were also identified in resistant accessions. A genetic analysis revealed that these alleles, as well as T-DNA insertions in the candidate gene, were responsible of RYMV resistance. All 40 resistant accessions were ultimately linked to a validated or candidate resistance allele at one of the three resistance genes to RYMV.

CONCLUSION

This study demonstrated that the RYMV2 resistance gene is homologous to the Arabidopsis CPR5 gene and revealed five new resistance alleles at this locus. It also confirmed the close association between resistance and an amino-acid substitution in the leucine-rich repeat of the NLR candidate for RYMV3. We also provide an extensive overview of the genetic diversity of resistance to RYMV in the O. glaberrima species, while underlining the contrasted pattern of diversity between O. glaberrima and O. sativa for this trait. The different resistance genes and alleles will be instrumental in breeding varieties with sustainable field resistance to RYMV.

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.

In-depth genome diversity, population structure and linkage disequilibrium analysis of worldwide diverse safflower (Carthamus tinctorius L.) accessions using NGS data generated by DArTseq technology

Safflower (Carthamus tinctorius L.) is one of the most important oilseed crops for its seed oil rich in unsaturated fatty acids. Precise utilization of diverse genetic resources is fundamental in breeding programs to improve high yield genotypes with desirable traits. In this study, for the first time we report successful application of DArTseq technology; an efficient genotyping-by-sequencing (NGS); to analysis genetic diversity and population structure of 89 safflower accessions from worldwide origins. Totally, 19,639 DArTseq markers (10,130 SilicoDArTs and 9509 SNPs) generated through DArTseq genotyping. After filtering the data, 3431 polymorphic DArTseq markers (1136 SilicoDArTs and 2295 SNPs) used for genetic diversity, population structure and linkage disequilibrium analysis in safflower genotypes. All the SilicoDArT and SNP markers showed high reproducibility and call rate. Polymorphism information content (PIC) values ranged from 0.1 to 0.5, while ≥ 0.50% of SilicoDArTs and ≥ 0.64% SNPs showed PIC values more than median. Genotypes grouping using DArTseq markers resulted in three distinct clusters. Results showed weak correlation between safflower diversity pattern and origins. Analysis of molecular variance revealed that the majority of genetic variation was attributed to the differences among varieties within cluster populations and there was no significant molecular variance between origins. However, safflower of accessions belonged to Iran, Turkey, Pakistan and India indeed appear to be genetically similar and grouped close in referred cluster, while the accessions from Near East (Afghanistan, China) being distinct. Our results were in agreement with hypothesis that safflower domesticated in somewhere west of Fertile Crescent and then expanded through Africa and Europe. Present study using a panel of globally diverse safflower accessions and large number of DArTseq markers set the stage for future analysis of safflower domestication using large germplasm from proposed domestication centers. Also, studied germplasm in this study can be used as a valuable source for future genomic studies in safflower for mapping desirable traits through genome-wide association mapping studies.

Whole-genome diversity, population structure and linkage disequilibrium analysis of globally diverse wheat genotypes using genotyping-by-sequencing DArTseq platform

In this study, 129 wheat genotypes from globally diverse origins were genotyped using DArTseq (SilicoDArT and SNP) markers. After filtering markers for quality-filtering, 14,270 SilicoDArTs and 6484 SNPs were retained and used for genetic diversity, population structure and linkage disequilibrium analyses. The highest number of SilicoDArT and SNP markers mapped on genome A and B compared to genome D. In both marker types, polymorphism information content (PIC) values ranged from 0.1 to 0.5, while > 0.80% of SilicoDArTs and > 0.44% SNPs showed PIC value more than median (0.25%). Un-weighted Neighbor Joining cluster analysis and Bayesian-based model population structure grouped wheat genotypes into three and four clusters, respectively. Principal component analysis and discriminant analysis of principal component results showed highly match with cluster and population structure analysis. Linkage disequilibrium (LD) was more extensive in both marker types, while graphical display of LD decay for both marker types showed that LD declined in the region close to 15 kbp, where r 2-values corresponded to r 2 = 0.16. Overall, our genetic diversity analysis showed high level of variation in studied wheat genotypes, even though there was no relationship between wheat grouping and origins. This might be attributed to admixture level that occurred during long-term natural selection of wheat genotypes in different parts of the world. Highly diverse wheat genotypes used in this study may possess unique genes and are useful sources in breeding programs to improve grain yield and quality.

Constructing a Core Collection of the Medicinal Plant Angelica biserrata Using Genetic and Metabolic Data

Angelica biserrata is an important medicinal plant in Chinese traditional medicine. Its roots, which are known as Duhuo in Chinese, are broadly applied to treat inflammation, arthritis, and headache. With increasing market demand, the wild resources of A. biserrata have been overexploited, and conservation, assessment of genetic resources and breeding for this species is needed. Here, we sequenced the transcriptome of A. biserrata and developed simple sequence repeat (SSR) markers from it to construct a core collection based on 208 samples collected from Changyang-related regions. A total of 132 alleles were obtained for 17 SSR loci used with the polymorphic information content (PIC) ranging from 0.44 to 0.83. Abundant genetic diversity was inferred by Shannon's information index (1.51), observed (0.57) and expected heterozygosity (0.72). The clustering analysis resulted into two sample groups and analysis of molecular variance (AMOVA) showed only 6% genetic variation existed among populations. A further metabolic analysis of these samples revealed the main coumarin contents, such as osthole and columbianadin. According to the genetic and metabolic data, we adopted the least distance stepwise sampling strategy to construct seven preliminary core collections, of which the 20CC collection, which possessed 42 A. biserrata individuals accounting for 90.20% of the genetic diversity of the original germplasm, represented the best core collection. This study will contribute to the conservation and management of A. biserrata wild germplasm resources and provide a material basis for future selection and breeding of this medicinal plant.

Identifying sex-linked markers in Litoria aurea: a novel approach to understanding sex chromosome evolution in an amphibian

Few taxa exhibit the variability of sex-determining modes as amphibians. However, due to the presence of homomorphic sex chromosomes in many species, this phenomenon has been difficult to study. The Australian frog, Litoria aurea, has been relatively well studied over the past 20 years due to widespread declines largely attributable to chytrid fungus. However, it has been subject to few molecular studies and its mode of sex determination remained unknown. We applied DArTseq™ to develop sex-linked single nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 44 phenotypically sexed L. aurea individuals from the Molonglo River in NSW, Australia. We conclusively identified a male heterogametic (XX-XY) sex determination mode in this species, identifying 11 perfectly sex-linked SNP and six strongly sex-linked PA markers. We identified a further 47 moderately sex-linked SNP loci, likely serving as evidence indicative of XY recombination. Furthermore, within these 47 loci, a group of nine males were found to have a feminised Y chromosome that significantly differed to all other males. We postulate ancestral sex-reversal as a means for the evolution of this now pseudoautosomal region on the Y chromosome. Our findings present new evidence for the ‘fountain of youth’ hypothesis for the retention of homomorphic sex chromosomes in amphibians and describe a novel approach for the study of sex chromosome evolution in amphibia.

Advances in Molecular Genetics and Genomics of African Rice (Oryza glaberrima Steud)

African rice (Oryza glaberrima) has a pool of genes for resistance to diverse biotic and abiotic stresses, making it an important genetic resource for rice improvement. African rice has potential for breeding for climate resilience and adapting rice cultivation to climate change. Over the last decade, there have been tremendous technological and analytical advances in genomics that have dramatically altered the landscape of rice research. Here we review the remarkable advances in knowledge that have been witnessed in the last few years in the area of genetics and genomics of African rice. Advances in cheap DNA sequencing technologies have fuelled development of numerous genomic and transcriptomic resources. Genomics has been pivotal in elucidating the genetic architecture of important traits thereby providing a basis for unlocking important trait variation. Whole genome re-sequencing studies have provided great insights on the domestication process, though key studies continue giving conflicting conclusions and theories. However, the genomic resources of African rice appear to be under-utilized as there seems to be little evidence that these vast resources are being productively exploited for example in practical rice improvement programmes. Challenges in deploying African rice genetic resources in rice improvement and the genomics efforts made in addressing them are highlighted.

Whole Genome Diversity, Population Structure, and Linkage Disequilibrium Analysis of Chickpea (Cicer arietinum L.) Genotypes Using Genome-Wide DArTseq-Based SNP Markers

Characterization of genetic diversity, population structure, and linkage disequilibrium is a prerequisite for proper management of breeding programs and conservation of genetic resources. In this study, 186 chickpea genotypes, including advanced "Kabuli" breeding lines and Iranian landrace "Desi" chickpea genotypes, were genotyped using DArTseq-Based single nucleotide polymorphism (SNP) markers. Out of 3339 SNPs, 1152 markers with known chromosomal position were selected for genome diversity analysis. The number of mapped SNP markers varied from 52 (LG8) to 378 (LG4), with an average of 144 SNPs per linkage group. The chromosome size that was covered by SNPs varied from 16,236.36 kbp (LG8) to 67,923.99 kbp (LG5), while LG4 showed a higher number of SNPs, with an average of 6.56 SNPs per Mbp. Polymorphism information content (PIC) value of SNP markers ranged from 0.05 to 0.50, with an average of 0.32, while the markers on LG4, LG6, and LG8 showed higher mean PIC value than average. Unweighted neighbor joining cluster analysis and Bayesian-based model population structure grouped chickpea genotypes into four distinct clusters. Principal component analysis (PCoA) and discriminant analysis of principal component (DAPC) results were consistent with that of the cluster and population structure analysis. Linkage disequilibrium (LD) was extensive and LD decay in chickpea germplasm was relatively low. A few markers showed r2 ≥ 0.8, while 2961 pairs of markers showed complete LD (r2 = 1), and a huge LD block was observed on LG4. High genetic diversity and low kinship value between pairs of genotypes suggest the presence of a high genetic diversity among the studied chickpea genotypes. This study also demonstrates the efficiency of DArTseq-based SNP genotyping for large-scale genome analysis in chickpea. The genotypic markers provided in this study are useful for various association mapping studies when combined with phenotypic data of different traits, such as seed yield, abiotic, and biotic stresses, and therefore can be efficiently used in breeding programs to improve chickpea.

Comparisons of molecular diversity indices, selective sweeps and population structure of African rice with its wild progenitor and Asian rice

The extent of molecular diversity parameters across three rice species was compared using large germplasm collection genotyped with genomewide SNPs and SNPs that fell within selective sweep regions. Previous studies conducted on limited number of accessions have reported very low genetic variation in African rice (Oryza glaberrima Steud.) as compared to its wild progenitor (O. barthii A. Chev.) and to Asian rice (O. sativa L.). Here, we characterized a large collection of African rice and compared its molecular diversity indices and population structure with the two other species using genomewide single nucleotide polymorphisms (SNPs) and SNPs that mapped within selective sweeps. A total of 3245 samples representing African rice (2358), Asian rice (772) and O. barthii (115) were genotyped with 26,073 physically mapped SNPs. Using all SNPs, the level of marker polymorphism, average genetic distance and nucleotide diversity in African rice accounted for 59.1%, 63.2% and 37.1% of that of O. barthii, respectively. SNP polymorphism and overall nucleotide diversity of the African rice accounted for 20.1-32.1 and 16.3-37.3% of that of the Asian rice, respectively. We identified 780 SNPs that fell within 37 candidate selective sweeps in African rice, which were distributed across all 12 rice chromosomes. Nucleotide diversity of the African rice estimated from the 780 SNPs was 8.3 × 10^-4, which is not only 20-fold smaller than the value estimated from all genomewide SNPs (π = 1.6 × 10^-2), but also accounted for just 4.1%, 0.9% and 2.1% of that of O. barthii, lowland Asian rice and upland Asian rice, respectively. The genotype data generated for a large collection of rice accessions conserved at the AfricaRice genebank will be highly useful for the global rice community and promote germplasm use.

Origins and geographic diversification of African rice (Oryza glaberrima)

Rice is a staple food for the majority of the world’s population. Whereas Asian rice (Oryza sativa) has been extensively studied, the exact origins of African rice (Oryza glaberrima) are still contested. Previous studies have supported either a centric or a non-centric geographic origin of African rice domestication. Here we review the evidence for both scenarios through a critical reassessment of 206 whole genome sequences of domesticated and wild African rice. While genetic diversity analyses support a severe bottleneck caused by domestication, signatures of recent and strong positive selection do not unequivocally point to candidate domestication genes, suggesting that domestication proceeded differently than in Asian rice–either by selection on different alleles, or different modes of selection. Population structure analysis revealed five genetic clusters localising to different geographic regions. Isolation by distance was identified in the coastal populations, which could account for parallel adaptation in geographically separated demes. Although genome-wide phylogenetic relationships support an origin in the eastern cultivation range followed by diversification along the Atlantic coast, further analysis of domestication genes shows distinct haplotypes in the southwest—suggesting that at least one of several key domestication traits might have originated there. These findings shed new light on an old controversy concerning plant domestication in Africa by highlighting the divergent roots of African rice cultivation, including a separate centre of domestication activity in the Guinea Highlands. We thus suggest that the commonly accepted centric origin of African rice must be reconsidered in favour of a non-centric or polycentric view.

Genetic diversity and population structure of a mini-core subset from the world cowpea (Vigna unguiculata (L.) Walp.) germplasm collection

The International Institute of Tropical Agriculture maintains the world's largest collection of cowpea germplasm of over 15,000 accessions. A sub-set of 298 lines from the loosely composed mini core collection of 370 landraces were genotyped based on genotyping by sequencing (GBS). Ward's minimum variance hierarchical cluster analysis, model-based ancestry analysis and discriminant analysis of principal component (DAPC) were carried out on this sub-set. Three clusters were identified by the different clustering methods. Principal component analysis further supported the three clusters especially when accessions are scattered along the axes of the first two principal components. The first two principal components explained a total of 22.30% of the variation. Cluster one comprises 115 accessions from the largest number of countries and has the highest gene diversity, heterozygosity and polymorphic information content (PIC) values. Cluster two is made up of 102 accessions, 90 percent of which are from West and Central Africa. Analysis of molecular variance shows that the most variation is among accessions and lowest among clusters. No cluster is made exclusively of accessions from a single country. Based on SNP markers, the sub set of cowpea mini core germplasm collection used in this study encompasses the diversity in the crop.

Development of species diagnostic SNP markers for quality control genotyping in four rice (Oryza L.) species

Species misclassification (misidentification) and handling errors have been frequently reported in various plant species conserved at diverse gene banks, which could restrict use of germplasm for correct purpose. The objectives of the present study were to (i) determine the extent of genotyping error (reproducibility) on DArTseq-based single-nucleotide polymorphisms (SNPs); (ii) determine the proportion of misclassified accessions across 3134 samples representing three African rice species complex (Oryza glaberrima, O. barthii, and O. longistaminata) and an Asian rice (O. sativa), which are conserved at the AfricaRice gene bank; and (iii) develop species- and sub-species (ecotype)-specific diagnostic SNP markers for rapid and low-cost quality control (QC) analysis. Genotyping error estimated from 15 accessions, each replicated from 2 to 16 times, varied from 0.2 to 3.1%, with an overall average of 0.8%. Using a total of 3134 accessions genotyped with 31,739 SNPs, the proportion of misclassified samples was 3.1% (97 of the 3134 accessions). Excluding the 97 misclassified accessions, we identified a total of 332 diagnostic SNPs that clearly discriminated the three indigenous African species complex from Asian rice (156 SNPs), O. longistaminata accessions from both O. barthii and O. glaberrima (131 SNPs), and O. sativa spp. indica from O. sativa spp. japonica (45 SNPs). Using chromosomal position, minor allele frequency, and polymorphic information content as selection criteria, we recommended a subset of 24 to 36 of the 332 diagnostic SNPs for routine QC genotyping, which would be highly useful in determining the genetic identity of each species and correct human errors during routine gene bank operations.

Assessment of Genetic Variation and Population Structure of Diverse Rice Genotypes Adapted to Lowland and Upland Ecologies in Africa Using SNPs

Using interspecific crosses involving Oryza glaberrima Steud. as donor and O. sativa L. as recurrent parents, rice breeders at the Africa Rice Center developed several 'New Rice for Africa (NERICA)' improved varieties. A smaller number of interspecific and intraspecific varieties have also been released as 'Advanced Rice for Africa (ARICA)'. The objective of the present study was to investigate the genetic variation, relatedness, and population structure of 330 widely used rice genotypes in Africa using DArTseq-based single nucleotide polymorphisms (SNPs). A sample of 11 ARICAs, 85 NERICAs, 62 O. sativa spp. japonica, and 172 O. sativa spp. indica genotypes were genotyped with 27,560 SNPs using diversity array technology (DArT)-based sequencing (DArTseq) platform. Nearly 66% of the SNPs were polymorphic, of which 15,020 SNPs were mapped to the 12 rice chromosomes. Genetic distance between pairs of genotypes that belong to indica, japonica, ARICA, and NERICA varied from 0.016 to 0.623, from 0.020 to 0.692, from 0.075 to 0.763, and from 0.014 to 0.644, respectively. The proportion of pairs of genotypes with genetic distance > 0.400 was the largest within NERICAs (35.1% of the pairs) followed by ARICAs (18.2%), japonica (17.4%), and indica (5.6%). We found one pair of japonica, 11 pairs of indica, and 35 pairs of NERICA genotypes differing by <2% of the total scored alleles, which was due to 26 pairs of genotypes with identical pedigrees. Cluster analysis, principal component analysis, and the model-based population structure analysis all revealed two distinct groups corresponding to the lowland (primarily indica and lowland NERICAs) and upland (japonica and upland NERICAs) growing ecologies. Most of the interspecific lowland NERICAs formed a sub-group, likely caused by differences in the O. glaberrima genome as compared with the indica genotypes. Analysis of molecular variance revealed very great genetic differentiation (F_ST = 0.688) between the lowland and upland ecologies, and 31.2% of variation attributable to differences within cluster groups. About 8% (1,197 of 15,020) of the 15,020 SNPs were significantly (P < 0.05) different between the lowland and upland ecologies and formed contrasting haplotypes that could clearly discriminate lowland from upland genotypes. This is the first study using high density markers that characterized NERICA and ARICA varieties in comparison with indica and japonica varieties widely used in Africa, which could aid rice breeders on parent selection for developing new improved rice germplasm.