Comparison of SSR and SNP markers in estimation of genetic diversity and population structure of Indian rice varieties

Identification of QTLs associated with yield-related traits and superior genotype prediction using recombinant inbred line population in tobacco

Tobacco is an economically significant industrial crop and model plant for genetic research, yet little is known about its genetic architecture. Quantitative trait loci (QTL) analysis was performed for six agronomic traits on an F_7 population of 341 genotypes, parents, and F1 plants using 1974 SSR markers across two environments. 31 QTLs contributing single-locus additive effects on 13 linkage groups (LGs) and 6 QTL pairs contributing epistatic effects on 6 LGs, were detected by the QTLNetwork 2.0 which was developed for the mixed-linear-model-based composite interval mapping (MCIM). Notably, 5 QTLs and 1 epistatic QTL pair were found to have pleiotropic effects on some genetically related traits. Moreover, the Broad sense heritability of the detected QTLs ranged from 1.05% to 43.33%, while genotype-by-environment interaction heritability spanned from 27.09% to 56.25%. Based on the results of QTL mapping, the potential superior lines for all or specific environments were designed and evaluated. Five major QTLs were finely dissected based on the tobacco reference genome of K326, and 31 candidate genes were predicted. This study offered new insights into the complicated genetic architecture and QTL resources for efficient breeding design for genetic improvement of agronomic traits in tobacco.

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

Estimation of Genetic Diversity and Number of Unique Genotypes of Cassava Germplasm from Burkina Faso Using Microsatellite Markers

Genetic diversity is very important in crop improvement. This study was carried out to assess the genetic diversity and the number of unique multilocus genotypes (MLGs) in a cassava collection in Burkina Faso. To achieve this objective, 130 cassava accessions were genotyped using 32 simple sequence repeat (SSR) markers. The results revealed that among these markers, twelve (12) were highly informative, with polymorphic information content (PIC) values greater than 0.50; twelve (12) were moderately informative, with PIC values ranging between 0.25 and 0.50; and eight (8) were not very informative, with PIC values lower than 0.25. A moderate level of genetic diversity was found for the population, indicated by the average expected heterozygosity (0.45) and the observed heterozygosity (0.48). About 83.8% of unique multilocus genotypes were found in the cassava collection, indicating that SSR markers seem to be most appropriate for MLG identification. Population structure analysis based on hierarchical clustering identified two subpopulations and the Bayesian approach suggested five clusters. Additionally, discriminant analysis of principal components (DAPC) separated the cassava accessions into 13 subpopulations. A comparison of these results and those of a previous study using single nucleotide polymorphisms (SNP) suggests that each type of marker can be used to assess the genetic structure of cassava grown in Burkina Faso.

Genetic differentiation of a southern Africa tepary bean (Phaseolus acutifolius A Gray) germplasm collection using high-density DArTseq SNP markers

Genetic resources of tepary bean (Phaseolus acutifolius A. Gray) germplasm collections are not well characterized due to a lack of dedicated genomic resources. There is a need to assemble genomic resources specific to tepary bean for germplasm characterization, heterotic grouping, and breeding. Therefore, the objectives of this study were to deduce the genetic groups in tepary bean germplasm collection using high-density Diversity Array Technology (DArT) based single nucleotide polymorphism (SNP) markers and select contrasting genotypes for breeding. Seventy-eight tepary bean accessions were genotyped using 10527 SNPs markers, and genetic parameters were estimated. Population structure was delineated using principal component and admixture analyses. A mean polymorphic information content (PIC) of 0.27 was recorded, indicating a relatively low genetic resolution of the developed SNPs markers. Low genetic variation (with a genetic distance [GD] = 0.32) existed in the assessed tepary bean germplasm collection. Population structure analysis identified five sub-populations through sparse non-negative matrix factorization (snmf) with high admixtures. Analysis of molecular variance indicated high genetic differentiation within populations (61.88%) and low between populations (38.12%), indicating high gene exchange. The five sub-populations exhibited variable fixation index (FST). The following genetically distant accessions were selected: Cluster 1:Tars-Tep 112, Tars-Tep 10, Tars-Tep 23, Tars-Tep-86, Tars-Tep-83, and Tars-Tep 85; Cluster 3: G40022, Tars-Tep-93, and Tars-Tep-100; Cluster 5: Zimbabwe landrace, G40017, G40143, and G40150. The distantly related and contrasting accessions are useful to initiate crosses to enhance genetic variation and for the selection of economic traits in tepary bean.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian cotton.

Molecular cytological analysis of alien introgressions in common wheat lines created by crossing of Triticum aestivum with T. dicoccoides and T. dicoccum

Wild and domesticated emmer (ВВАА, 2n = 28) are of significant interest for expanding the genetic diversity of common wheat as sources of a high protein and microelement grain content, resistance to many biotic and abiotic factors. Particular interest in these species is also determined by their close relationship with Triticum aestivum L., which facilitates interspecific hybridization. The objective of this work was to analyze the nature of alien introgressions in hybrid lines from crossing common wheat varieties with T. dicoccoides and T. dicoccum, and to assess the effect of their genome fragments on the cytological stability of introgression lines. A C-banding technique and genotyping with SNP and SSR markers were used to determine localization and length of introgression fragments. Assessment of cytological stability was carried out on the basis of chromosome behavior in microsporogenesis. A molecular cytogenetic analysis of introgression wheat lines indicated that the inclusion of the genetic material of wild and domesticated emmer was carried out mainly in the form of whole arms or large fragments in the chromosomes of the B genome and less extended inserts in the A genome. At the same time, the highest frequency of introgressions of the emmer genome was observed in chromosomes 1A, 1B, 2B, and 3B. The analysis of the final stage of meiosis showed a high level of cytological stability in the vast majority of introgression wheat lines (meiotic index was 83.0-99.0 %), which ensures the formation of functional gametes in an amount sufficient for successful reproduction. These lines are of interest for the selection of promising material with agronomically valuable traits and their subsequent inclusion in the breeding process.

Genetic Diversity Analysis of Soybean Collection Using Simple Sequence Repeat Markers

Soybean [Glycine max (L.) Merr.] is a nutrient-rich crop that offers a sustainable source of dietary protein and edible oil. Determining the level of genetic diversity and relationships between various genetic resources involved in breeding programs is very important in crop improvement strategies. This study evaluated 100 soybean accessions with diverse origins for 10 important agronomic traits, including plant height (PH), an important plant adaptation-related trait impacting yield, in conditions in southeastern Kazakhstan for 2 years. The comparison of different groups of PH (tall, middle, and short) using a t-test suggested that the group of plants with the tallest PH provided a higher yield (p < 0.001) in relatively dry field conditions. The genetic diversity of the accessions was estimated using 25 simple sequence repeat (SSR) markers previously known to be associated with plant height. The results showed a significant variation among different groups of origin for all measured agronomic traits, as well as high genetic diversity, with the PIC (polymorphism information content) varying from 0.140 to 0.732, with an average of 0.524. Nei's diversity index ranged between 0.152 and 0.747, with an average of 0.526. The principal coordinate analysis (PCoA) of the studied soybean collection showed that Kazakhstan accessions were genetically distant from European, East Asian, and North American cultivars. Twelve out of twenty-five SSR markers demonstrated significant associations with ten studied agronomic traits, including PH (p < 0.05). Six SSRs with pleiotropic effects for studied traits were selected, and their haplotypes with phenotypic effects were generated for each soybean accession. The obtained results can be used in soybean improvement programs, including molecular-assisted breeding projects.

A Comprehensive Genome-Wide Investigation of the Cytochrome 71 (OsCYP71) Gene Family: Revealing the Impact of Promoter and Gene Variants (Ser33Leu) of OsCYP71P6 on Yield-Related Traits in Indica Rice (Oryza sativa L.)

The cytochrome P450 (CYP450) gene family plays a critical role in plant growth and developmental processes, nutrition, and detoxification of xenobiotics in plants. In the present research, a comprehensive set of 105 OsCYP71 family genes was pinpointed within the genome of indica rice. These genes were categorized into twelve distinct subfamilies, where members within the same subgroup exhibited comparable gene structures and conserved motifs. In addition, 105 OsCYP71 genes were distributed across 11 chromosomes, and 36 pairs of OsCYP71 involved in gene duplication events. Within the promoter region of OsCYP71, there exists an extensive array of cis-elements that are associated with light responsiveness, hormonal regulation, and stress-related signaling. Further, transcriptome profiling revealed that a majority of the genes exhibited responsiveness to hormones and were activated across diverse tissues and developmental stages in rice. The OsCYP71P6 gene is involved in insect resistance, senescence, and yield-related traits in rice. Hence, understanding the association between OsCYP71P6 genetic variants and yield-related traits in rice varieties could provide novel insights for rice improvement. Through the utilization of linear regression models, a total of eight promoters were identified, and a specific gene variant (Ser33Leu) within OsCYP71P6 was found to be linked to spikelet fertility. Additionally, different alleles of the OsCYP71P6 gene identified through in/dels polymorphism in 131 rice varieties were validated for their allelic effects on yield-related traits. Furthermore, the single-plant yield, spikelet number, panicle length, panicle weight, and unfilled grain per panicle for the OsCYP71P6-1 promoter insertion variant were found to contribute 20.19%, 13.65%, 5.637%, 8.79%, and 36.86% more than the deletion variant, respectively. These findings establish a robust groundwork for delving deeper into the functions of OsCYP71-family genes across a range of biological processes. Moreover, these findings provide evidence that allelic variation in the promoter and amino acid substitution of Ser33Leu in the OsCYP71P6 gene could potentially impact traits related to rice yield. Therefore, the identified promoter variants in the OsCYP71P6 gene could be harnessed to amplify rice yields.

Genetic diversity and population structure analysis in cultivated soybean (Glycine max [L.] Merr.) using SSR and EST-SSR markers

Soybean (Glycine max) is an important legume that is used to fulfill the need of protein and oil of large number of population across the world. There are large numbers of soybean germplasm present in the USDA germplasm resources. Finding and understanding genetically diverse germplasm is a top priority for crop improvement programs. The current study used 20 functional EST-SSR and 80 SSR markers to characterize 96 soybean accessions from diverse geographic backgrounds. Ninety-six of the 100 markers were polymorphic, with 262 alleles (average 2.79 per locus). The molecular markers had an average polymorphic information content (PIC) value of 0.44, with 28 markers ≥ 0.50. The average major allele frequency was 0.57. The observed heterozygosity of the population ranged from 0-0.184 (average 0.02), while the expected heterozygosity ranged from 0.20-0.73 (average 0.51). The lower value for observed heterozygosity than expected heterozygosity suggests the likelihood of a population structure among the germplasm. The phylogenetic analysis and principal coordinate analysis (PCoA) divided the total population into two major groups (G1 and G2), with G1 comprising most of the USA lines and the Australian and Brazilian lines. Furthermore, the phylogenetic analysis and PCoA divided the USA lines into three major clusters without any specific differentiation, supported by the model-based STRUCTURE analysis. Analysis of molecular variance (AMOVA) showed 94% variation among individuals in the total population, with 2% among the populations. For the USA lines, 93% of the variation occurred among individuals, with only 2% among lines from different US states. Pairwise population distance indicated more similarity between the lines from continental America and Australia (189.371) than Asia (199.518). Overall, the 96 soybean lines had a high degree of genetic diversity.

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.

Genetic diversity and population structure of Uganda's yam (Dioscorea spp.) genetic resource based on DArTseq

Assessing the genetic diversity of yam germplasm from different geographical origins for cultivation and breeding purposes is an essential step for crop genetic resource conservation and genetic improvement, especially where the crop faces minimal attention. This study aimed to classify the population structure, and assess the extent of genetic diversity in 207 Dioscorea rotundata genotypes sourced from three different geographical origins. A total of 4,957 (16.2%) single nucleotide polymorphism markers were used to assess genetic diversity. The SNP markers were informative, with polymorphic information content ranging from 0.238 to 0.288 and a mean of 0.260 across all the genotypes. The observed and expected heterozygosity was 0.12 and 0.23, respectively while the minor allele frequency ranged from 0.093 to 0.124 with a mean of 0.109. The principal coordinate analysis, model-based structure and discriminant analysis of principal components, and the Euclidean distance matrix method grouped 207 yam genotypes into three main clusters. Genotypes from West Africa (Ghana and Nigeria) had significant similarities with those from Uganda. Analysis of molecular variance revealed that within-population variation across three different geographical origins accounted for 93% of the observed variation. This study, therefore, showed that yam improvement in Uganda is possible, and the outcome will constitute a foundation for the genetic improvement of yams in Uganda.

Comparative analysis of the efficiency of seed protein profiles in assessing genetic variation and population structure among indigenous Manipur black rice cultivars

BACKGROUND

The state of Manipur, North East India has distinct topology of hill and valley regions with vast agroclimatic variability, being considered as one of the centers of rice diversity. The indigenous Manipur black rice cultivars exhibit wide range of diversity in morphology, pericarp color, shape and size of grain, aroma, glutinous or non-glutinous features but remain less characterised. Many of these cultivars, such as those named Chakhao, are endowed with multiple health benefits due to high anthocyanins, and hold special importance for the local people. It is important to analyse the genetic diversity and population structure for this germplasm with unique allelic combinations to utilize in rice breeding programme.

METHODS AND RESULTS

We characterized total soluble seed protein fractions to not only fingerprint the 45 indigenous black rice cultivars but assess their genetic relatedness. Cluster analyses generated mainly two groups, complemented by PCoA scatter plot ascertaining geographical distinction. The hill black rice were more diverse. The population structure analysis revealed seven subpopulations indicating high genetic variability. The 24 polymorphic bands were scored in the range of 127.8 to 10.3 kDa comprised of four protein fractions. Three polypeptide bands each were ascribed to known fractions of glutelins and prolamins, while one band each could be described for albumin and globulin fractions, besides other diagnostic bands.

CONCLUSION

Some diverse cultivars were Amubi, Chedo Anal, Chipi Buh, Athebu, Poireton, BuPu Mui, Kotha Chahao II. These cultivars can be used in future black rice breeding programmes. This can further prevent genetic erosion and protect intellectual property rights.

RP, R. L. R. Population structure and genetic diversity characterization of soybean for seed longevity

Seed longevity is an important trait in the context of germplasm conservation and economics of seed production. The identification of populations with high level of genetic variability for seed longevity and associated traits will become a valuable resource for superior alleles for seed longevity. In this study, Genotyping-by-sequencing (GBS)-single nucleotide polymorphism (SNP) approach, simple sequence repeats (SSR) markers and agro-morphological traits have been explored to investigate the diversity and population structure of assembled 96 genotypes. The GBS technique performed on 96 genotypes of soybean (Glycine max (L.) Merrill) resulted in 37,897 SNPs on sequences aligned to the reference genome sequence. The average genome coverage was 6.81X with a mapping rate of 99.56% covering the entire genome. Totally, 29,955 high quality SNPs were identified after stringent filtering and most of them were detected in non-coding regions. The 96 genotypes were phenotyped for eight quantitative and ten qualitative traits by growing in field by following augmented design. The STRUCTURE (Bayesian-model based algorithm), UPGMA (Un-weighed Pair Group Method with Arithmetic mean) and principal component analysis (PCA) approaches using SSR, SNP as well as quantitative and qualitative traits revealed population structure and diversity in assembled population. The Bayesian-model based STRUCTURE using SNP markers could effectively identify clusters with higher seed longevity associated with seed coat colour and size which were subsequently validated by UPGMA and PCA based on SSR and agro-morphological traits. The results of STRUCTURE, PCA and UPGMA cluster analysis showed high degree of similarity and provided complementary data that helped to identify genotypes with higher longevity. Six black colour genotypes, viz., Local black soybean, Kalitur, ACC Nos. 39, 109, 101 and 37 showed higher seed longevity during accelerated ageing. Higher coefficient of variability observed for plant height, number of pods per plant, seed yield per plant, 100 seed weight and seed longevity confirms the diversity in assembled population and its suitability for quantitative trait loci (QTL) mapping.

DArTSeq SNP-based genetic diversity and population structure studies among taro [(Colocasia esculenta (L.) Schott] accessions sourced from Nigeria and Vanuatu

Taro is a valuable staple food crop among resource-poor rural people in countries such as Nigeria and Ghana, among others. Characterization of genetic diversity is a prerequisite for proper management of breeding programs and conservation of genetic resources. Two hundred seventy one taro accessions obtained from Nigeria and Vanuatu were genotyped using DArTseq-based SNP markers with the objectives of investigating the genetic diversity and population structure. In the analysis, 10,391 SNP markers were filtered from the sequence and used. The analysis revealed higher transition than transversion types of SNPs in the ratio of 1.43:1. The polymorphism ranged from 0.26 to 0.29 for the markers, indicating moderate genetic diversity. A model-based Bayesian clustering analysis of taro accessions yielded five subgroups and revealed the admixture situation in 19.19% of all accessions in the study. Vanuatu taro accessions exhibited more genetic diversity than Nigerian taro accessions. The population diversity estimate (PhiPt) was relatively higher (0.52) for accessions originating from Vanuatu than for Nigerian accessions. Analysis of molecular variance (AMOVA) revealed that most variation existed among individuals within a population at 52%. Nei’s genetic distance showed that relatedness is based on geographical proximity. Collection of taro genetic resources should give more emphasis to within regions to utilize diversity in taro breeding program. This study also demonstrated the efficiency of DArTseq-based SNP genotyping for large-scale genome analysis in taro. The genotypic markers provided in this study are useful for association mapping studies.

Molecular characterization of intergeneric hybrids between Malus and Pyrus

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

Descriptive Genomic Analysis and Sequence Genotyping of the Two Papaya Species (Vasconcellea pubescens and Vasconcellea chilensis) Using GBS Tools

A genotyping by sequencing (GBS) approach was used to analyze the organization of genetic diversity in V. pubescens and V. chilensis. GBS identified 4675 and 4451 SNPs/INDELs in two papaya species. The cultivated orchards of V. pubescens exhibited scarce genetic diversity and low but significant genetic differentiation. The neutrality test yielded a negative and significant result, suggesting that V. pubescens suffered a selective sweep or a rapid expansion after a bottleneck during domestication. In contrast, V. chilensis exhibited a high level of genetic diversity. The genetic differentiation among the populations was slight, but it was possible to distinguish the two genetic groups. The neutrality test indicated no evidence that natural selection and genetic drift affect the natural population of V. chilensis. Using the Carica papaya genome as a reference, we identified critical SNPs/INDELs associated with putative genes. Most of the identified genes are related to stress responses (salt and nematode) and vegetative and reproductive development. These results will be helpful for future breeding and conservation programs of the Caricaceae family.

Consequences of introgression and gene flow on the genetic structure and diversity of Lima bean (Phaseolus lunatus L.) in its Mesoamerican diversity area

We evaluated the role of gene flow and wild-crop introgression on the structure and genetic diversity of Lima bean (Phaseolus lunatus) in the Yucatan Peninsula, an important Mesoamerican diversity area for this crop, using a genotyping-by-sequencing approach (15,168 SNP markers) and two scales. At the local scale, STRUCTURE and NGSEP analyses showed predominantly crop-to-wild introgression, but also evidence of a bidirectional gene flow in the two wild-weedy-crop complexes studied (Itzinté and Dzitnup). The ABBA-BABA tests showed a higher introgression in Itzinté (the older complex) than in Dzitnup (the younger one); at the allelic level, the wild-crop introgression in Itzinté was similar in both directions, in Dzitnup it was higher from crop-to-wild; and at the chromosomal level, introgression in Itzinté was from wild-to-crop, whereas in Dzitnup it occured in the opposite direction. Also, we found H _E values slightly higher in the domesticated accessions than in the wild ones, in both complexes (Itzinté: wild = 0.31, domesticated = 0.34; Dzinup: wild = 0.27, domesticated = 0.36), but %P and π estimators were higher in the wild accessions than in the domesticated ones. At a regional scale, STRUCTURE and MIGRATE showed a low gene flow, predominantly from crop-to-wild; and STRUCTURE, Neighbor-Joining and PCoA analyses indicated the existence of two wild groups and one domesticated group, with a marked genetic structure based in the existence of domesticated MI and wild MII gene pools. Also, at the regional scale, we found a higher genetic diversity in the wild accessions than in the domesticated ones, in all estimators used (e.g., H _E = 0.27 and H _E = 0.17, respectively). Our results indicate that gene flow and introgression are playing an important role at the local scale, but its consequences on the structure and genetic diversity of the Lima bean are not clearly reflected at the regional scale, where diversity patterns between wild and domesticated populations could be reflecting historical events.

RAD-seq as an effective strategy for heterogenous variety identification in plants-a case study in Italian Ryegrass (Lolium multiflorum)

The primary approach for variety distinction in Italian ryegrass is currently the DUS (distinctness, uniformity and stability) test based on phenotypic traits. Considering the diverse genetic background within the population and the complexity of the environment, however, it is challenging to accurately distinguish varieties based on DUS criteria alone. In this study, we proposed the application of high-throughput RAD-seq to distinguish 11 Italian ryegrass varieties with three bulks of 50 individuals per variety. Our findings revealed significant differences among the 11 tested varieties. The PCA, DAPC and STRUCTURE analysis indicated a heterogeneous genetic background for all of them, and the AMOVA analysis also showed large genetic variance among these varieties (Φ_ST = 0.373), which were clearly distinguished based on phylogenetic analysis. Further nucleotide diversity (Pi) analysis showed that the variety 'Changjiang No.2' had the best intra-variety consistency among 11 tested varieties. Our findings suggest that the RAD-seq could be an effectively alternative method for the variety distinction of Italian ryegrass, as well as a potential tool for open-pollinated varieties (OPVs) of other allogamous species.

Two fingerprinting sets for Humulus lupulus based on KASP and microsatellite markers

Verification of clonal identity of hop (Humulus lupulus L.) cultivars within breeding programs and germplasm collections is vital to conserving genetic resources. Accurate and economic DNA-based tools are needed in dioecious hop to confirm identity and parentage, neither of which can be reliably determined from morphological observations. In this study, we developed two fingerprinting sets for hop: a 9-SSR fingerprinting set containing high-core repeats that can be run in a single PCR reaction and a kompetitive allele specific PCR (KASP) assay of 25 single nucleotide polymorphisms (SNPs). The SSR set contains a sex-linked primer pair, HI-AGA7, that was used to genotype 629 hop accessions from the US Department of Agriculture (USDA) National Clonal Germplasm Repository (NCGR), the USDA Forage Seed and Cereal Research (FSCR), and the University of Nebraska-Lincoln (UNL) collections. The SSR set identified unique genotypes except for 89 sets of synonymous samples. These synonyms included: cultivars with different designations, the same cultivars from different sources, heat-treated clones, and clonal variants. Population structure analysis clustered accessions into wild North American (WNA) and cultivated groups. Diversity was slightly higher in the cultivated samples due to larger sample size. Parentage and sib-ship analyses were used to identify true-to-type cultivars. The HI-AGA7 marker generated two male- and nine female-specific alleles among the cultivated and WNA samples. The SSR and KASP fingerprinting sets were compared in 190 samples consisting of cultivated and WNA accession for their ability to confirm identity and assess diversity and population structure. The SSR fingerprinting set distinguished cultivars, selections and WNA accessions while the KASP assays were unable to distinguish the WNA samples and had lower diversity estimates than the SSR set. Both fingerprinting sets are valuable tools for identity confirmation and parentage analysis in hop for different purposes. The 9-SSR assay is cost efficient when genotyping a small number of wild and cultivated hop samples (<96) while the KASP assay is easy to interpret and cost efficient for genotyping a large number of cultivated samples (multiples of 96).

Genetic Diversity Revealed by Microsatellites in Genus Carya

The genus Carya consists of 17 species divided into 3 sections: Carya or the true hickories, Apocarya or the pecan hickories, and Sinocarya or the Asian hickories. Interspecific hybrids exist and have been used in pecan cultivar development. Nuclear and plastid microsatellite or SSR markers have been useful in distinguishing species, sections, and populations. They provide evidence for hybridity between species and can confirm heredity within crosses. As more sophisticated methods of genomic evaluation are cooperatively developed for use in pecan breeding and selection, the use of these methods will be supplemented and informed by the lessons provided by microsatellite markers, as interpreted across broad germplasm collections. In this study, over 400 Carya accessions from diverse diploid and tetraploid taxa and their interspecific hybrids, maintained at the USDA National Collection of Genetic Resources for Carya (NCGR-Carya), were analyzed using 14 nuclear and 3 plastid microsatellite markers. Principal coordinate analysis showed clear taxonomic classifications at multiple taxonomic levels along with patterns of interspecific hybridity. Evidence was also found for genetic differences associated with geographic distribution. The results indicate that this group of markers is useful in examining and characterizing populations and hybrids in the genus Carya and may help delineate the composition of a core collection to help characterize the NCGR-Carya repository collection for use in its pecan breeding program. The SSR fingerprints of the inventories of the USDA NCGR-Carya repository can also be used as a reference for identifying unknown pecan trees for growers.

Population Structure and Relatedness for Genome-Wide Association Studies

The estimation of the population structure and genetic relatedness between individuals within a collection of accessions is important in the formation of core collections for the conservation of genetic resources, uncovering the demographic history of the population under study, as well as for association studies. With the recent development of high-throughput genotyping technologies, several algorithms and methods have been developed and implemented in software to estimate the extent of genetic diversity between individuals. In this chapter, our objective is to describe methods to capture population structure and relatedness in a step-by-step fashion. To exemplify the process, two pruned datasets (14K and 243K SNP markers) were used to investigate population structure and relatedness among a soybean GWAS panel using different approaches and methods.

Grapevine Diversity and Genetic Relationships in Northeast Portugal Old Vineyards

More than 100 grapevine varieties are registered as suitable for wine production in “Douro” and “Trás-os-Montes” Protected Designations of Origin regions; however, only a few are actually used for winemaking. The identification of varieties cultivated in past times can be an important step to take advantage of all the potential of these regions grape biodiversity. The conservation of the vanishing genetic resources boosts greater product diversification, and it can be considered strategic in the valorisation of these wine regions. Hence, one goal of the present study was to prospect and characterise, through molecular markers, 310 plants of 11 old vineyards that constitute a broad representation of the grape genetic patrimony of “Douro” and “Trás-os-Montes” wine regions; 280 samples, grouped into 52 distinct known varieties, were identified through comparison of their genetic profiles generated via 6 nuclear SSR and 43 informative SNP loci amplification; the remaining 30 samples, accounting for 13 different genotypes, did not match with any profile in the consulted databases and were considered as new genotypes. This study also aimed at evaluating the population structure among the 65 non-redundant genotypes identified, which were grouped into two ancestral genetic groups. The mean probability of identity values of 0.072 and 0.510 (for the 6 SSR and 226 SNP sets, respectively) were determined. Minor differences were observed between frequencies of chlorotypes A and D within the non-redundant genotypes studied. Twenty-seven pedigrees were confirmed and nine new trios were established. Ancestors of eight genotypes remain unknown.

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.

Evaluation of gene flow and genetic diversity in rice accessions across Kurdistan region-iraq using SSR markers

Background

In recent years, farmers have complained that the only way to obtain seeds is to select plants that show good performance under local climate conditions in the region. This study aimed to investigate the diversity of rice accessions grown in the region to build a breeding program.

Methods and results

A total of 62 accessions of rice from farmers and research stations were collected from the Kurdistan region, including short-grain and long-grain types, for molecular genetics and diversity analysis. In this study, 37 polymorphic simple sequence repeat (SSR) markers were selected with several molecular genetics software programs. The results show that these SSR markers are very effective for this investigation, generating a total of 152 observed alleles (Na), 75.166 Effective number of alleles (Ne) and an average of 4.1 and 2.03 alleles per locus, respectively. The average polymorphic information content (PIC) per locus was recorded as 0.404. The research presented here confirms two subpopulations, japonica (C1 and C2) and indica (C3), based on molecular genetics data analysis. Analysis of molecular variance revealed that the 72% variance was due to the variation among populations and 28% within the population.

Conclusions

Altogether, these results indicate that there is very low gene flow. These results show the importance of the study of genetic diversity and relationships for starting breeding and improvement programs for rice in the Kurdistan region.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06920-x.

The details are in the genome-wide SNPs: Fine scale evolution of the Malaysian weedy rice

Weedy rice (Oryza spp.) is a major nuisance to rice farmers from all over the world. Although the emergence of weedy rice in East Malaysia on the island of Borneo is very recent, the threat to rice yield has reached an alarming stage. Using 47,027 genotyping-by-sequencing (GBS)-derived SNPs and candidate gene analysis of the plant architecture domestication gene TAC1, we assessed the genetic variations and evolutionary origin of weedy rice in East Malaysia. Our findings revealed two major evolutionary paths for genetically distinct weedy rice types. Whilst the cultivar-like weedy rice are very likely to be the weedy descendant of local coexisting cultivars, the wild-like weedy rice appeared to have arisen through two possible routes: (i) accidental introduction from Peninsular Malaysia weedy rice populations, and (ii) weedy descendants of coexisting cultivars. The outcome of our genetic analyses supports the notion that Sabah cultivars and Peninsular Malaysia weedy rice are the potential progenitors of Sabah weedy rice. Similar TAC1 haplotypes were shared between Malaysian cultivated and weedy rice populations, which further supported the findings of our GBS-SNP analyses. These different strains of weedy rice have convergently evolved shared traits, such as seeds shattering and open tillers. A comparison with our previous simple-sequence repeat-based population genetic analyses highlights the strength of genome-wide SNPs, including detection of admixtures and low-level introgression events. These findings could inform better strategic management for controlling the spread of weedy rice in the region.

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.

SSR marker-based genetic diversity analysis and SNP haplotyping of genes associating abiotic and biotic stress tolerance, rice growth and development and yield across 93 rice landraces

BACKGROUND

As rice is the staple food for more than half of the world's population, enhancing grain yield irrespective of the variable climatic conditions is indispensable. Many traditionally cultivated rice landraces are well adapted to severe environmental conditions and have high genetic diversity that could play an important role in crop improvement.

METHODS AND RESULTS

The present study revealed a high level of genetic diversity among the unexploited rice landraces cultivated by the farmers of Kerala. Twelve polymorphic markers detected a total of seventy- seven alleles with an average of 6.416 alleles per locus. Polymorphic Information Content (PIC) value ranged from 0.459 to 0.809, and to differentiate the rice genotypes, RM 242 was found to be the most appropriate marker with a high value of 0.809. The current study indicated that the rice landraces are highly diverse with higher values of the adequate number of alleles, PIC, and Shannon information index. Utilizing these informative SSR markers for future molecular characterization and population genetic studies in rice landraces are advisable. Haplotypes are sets of genomic regions within a chromosome inherited together, and haplotype-based breeding is a promising strategy for designing next-generation rice varieties. Here, haplotype analysis explored 270 haplotype blocks and 775 haplotypes from all the chromosomes of landraces under study. The number of SNPs in each haplotype block ranged from two to 28. Haplotypes of genes related to biotic and abiotic stress tolerance, yield-enhancing, and growth and development in rice landraces were also elucidated in the current study.

CONCLUSIONS

The present investigation revealed the genetic diversity of rice landraces and the haplotype analysis will open the way for genome-wide association studies, QTL identification, and marker-assisted selection in the unexplored rice landraces collected from Kerala.

Molecular genetic diversity and population structure analyses of rutabaga accessions from Nordic countries as revealed by single nucleotide polymorphism markers

Background

Rutabaga or swede (Brassica napus ssp. napobrassica (L.) Hanelt) varies in root and leaf shape and colour, flesh colour, foliage growth habits, maturity date, seed quality parameters, disease resistance and other traits. Despite these morphological differences, no in-depth molecular analyses of genetic diversity have been conducted in this crop. Understanding this diversity is important for conservation and broadening the use of this resource.

Results

This study investigated the genetic diversity within and among 124 rutabaga accessions from five Nordic countries (Norway, Sweden, Finland, Denmark and Iceland) using a 15 K single nucleotide polymorphism (SNP) Brassica array. After excluding markers that did not amplify genomic DNA, monomorphic and low coverage site markers, the accessions were analyzedwith 6861 SNP markers. Allelic frequency statistics, including polymorphism information content (PIC), minor allele frequency (MAF) and mean expected heterozygosity (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{H} $$\end{document}H¯e) and population differentiation statistics such as Wright’s F-statistics (F_ST) and analysis of molecular variance (AMOVA) indicated that the rutabaga accessions from Norway, Sweden, Finland and Denmark were not genetically different from each other. In contrast, accessions from these countries were significantly different from the accessions from Iceland (P < 0.05). Bayesian analysis with the software STRUCTURE placed 66.9% of the rutabaga accessions into three to four clusters, while the remaining 33.1% constituted admixtures. Three multivariate analyses: principal coordinate analysis (PCoA), the unweighted pair group method with arithmetic mean (UPGMA) and neighbour-joining (NJ) clustering methods grouped the 124 accessions into four to six subgroups.

Conclusion

Overall, the correlation of the accessions with their geographic origin was very low, except for the accessions from Iceland. Thus, Icelandic rutabaga accessions can offer valuable germplasm for crop improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07762-4.

Multiple DNA marker-assisted diversity analysis of Indian mango (Mangifera indica L.) populations

Arbitrary (65 RAPD, 25 ISSR, 23 DAMD), gene-targeted (22 SCoT, 33 CBDP) and co-dominant sequence specific (40 SSR) markers were used individually, or in combinations, to examine the genetic variability within and among 70 selected Indian mango genotypes based on geographic origin (East India, West India, North India, South India) and fruit status (Selection, Hybrid, Landrace). The highest genetic variability was demonstrated by the East Indian populations, followed by those from South India, West India, and North India, when measured in terms of Na, Ne, H, I, PB%, Ht and Hs. Interestingly, the local genotypes of Odisha, which forms a part of East Indian populations, showed the highest diversity compared to hybrid or selection groups, suggesting that the indigenous genotypes hold a greater potential for exploiting the unique and favourable alleles. The maximum genetic variability was detected in geographical/fruit status populations with SSRs (Na-1.76/1.88, Ne-1.48/1.51, H-0.28/0.30, I-0.41/0.45, PB%-76.1/86.9, Ht-0.31/0.32 and Hs-0.28/0.30), followed by CBDPs and SCoTs, reflecting their preeminence for examining the level of genetic polymorphism and diversity. Outcome of AMOVA based analyses as well as low-to-moderate coefficient of genetic differentiation (Gst) and high gene flow (Nm) indicated a greater amount of intra-population genetic variation compared to heterogeneity at inter-population level. Information generated through this investigation could facilitate conservation and further exploitation of mango germplasm including genetic improvement through breeding.

Development of CAPS Markers for Evaluation of Genetic Diversity and Population Structure in the Germplasm of Button Mushroom (Agaricus bisporus)

Agaricus bisporus is a globally cultivated mushroom with high economic value. Despite its widespread cultivation, commercial button mushroom strains have little genetic diversity and discrimination of strains for identification and breeding purposes is challenging. Molecular markers suitable for diversity analyses of germplasms with similar genotypes and discrimination between accessions are needed to support the development of new varieties. To develop cleaved amplified polymorphic sequences (CAPs) markers, single nucleotide polymorphism (SNP) mining was performed based on the A. bisporus genome and resequencing data. A total of 70 sets of CAPs markers were developed and applied to 41 A. bisporus accessions for diversity, multivariate, and population structure analyses. Of the 70 SNPs, 62.85% (44/70) were transitions (G/A or C/T) and 37.15% (26/70) were transversions (A/C, A/T, C/G, or G/T). The number of alleles per locus was 1 or 2 (average = 1.9), and expected heterozygosity and gene diversity were 0.0-0.499 (mean = 0.265) and 0.0-0.9367 (mean = 0.3599), respectively. Multivariate and cluster analyses of accessions produced similar groups, with F-statistic values of 0.134 and 0.153 for distance-based and model-based groups, respectively. A minimum set of 10 markers optimized for accession identification were selected based on high index of genetic diversity (GD, range 0.299-0.499) and major allele frequency (MAF, range 0.524-0.817). The CAPS markers can be used to evaluate genetic diversity and population structure and will facilitate the management of emerging genetic resources.

Population Genetic Analyses of Botrytis cinerea Isolates From Michigan Vineyards Using a High-Throughput Marker System Approach

As sequencing costs continue to decrease, new tools are being developed for assessing pathogen diversity and population structure. Traditional marker types, such as microsatellites, are often more cost effective than single-nucleotide polymorphism (SNP) panels when working with small numbers of individuals, but may not allow for fine scale evaluation of low or moderate structure in populations. Botrytis cinerea is a necrotrophic plant pathogen with high genetic variability that can infect more than 200 plant species worldwide. A panel of 52 amplicons were sequenced for 82 isolates collected from four Michigan vineyards representing 2 years of collection and varying fungicide resistance. A panel of nine microsatellite markers previously described was also tested across 74 isolates from the same population. A microsatellite and SNP marker analysis of B. cinerea populations was performed to assess the genetic diversity and population structure of Michigan vineyards, and the results from both marker types were compared. Both methods were able to detect population structure associated with resistance to the individual fungicides thiabendazole and boscalid, and multiple fungicide resistance (MFR). Microsatellites were also able to differentiate population structure associated with another fungicide, fluopyram, while SNPs were able to additionally differentiate structure based on year. For both methods, AMOVA results were similar, with microsatellite results explaining a smaller portion of the variation compared with the SNP results. The SNP-based markers presented here were able to successfully differentiate population structure similar to microsatellite results. These SNP markers represent new tools to discriminate B. cinerea isolates within closely related populations using multiple targeted sequences.

Genetic variability and biodiversity of Ukrainian Gray cattle by the BoLA-DRB3 gene

At the current stage of genetic studies of cattle, more and more attention is being drawn to autochthonous breeds. Native cattle have a number of prominent phenotypic traits and have preserved unique genes and their combinations lost by modern commercial breeds, which would be valuable to use in selective programs. We surveyed polymorphism of the Ukrainian autochthonous Gray breed according to alleles of exon 2 of the BoLA-DRB3 gene. The uniqueness of the gene lies in the broad variability of its allele variants. Significant informativeness at DNA level is quite important for genetic studies. We surveyed allele polymorphism using the PCR-RLFP method on DNA isolated from 88 samples of blood of cows and 5 samples of sperm. We identified 28 alleles, of which 23 variants were nomenclature ones and 5 (jba, *jab, *jbb, *nad and *nda) were “without established nomenclature”, their share accounting for 8.9%. Four alleles *06, *12, *16 and *jba had a frequency above 5% and occupied 69.9% of the breed’s allele fund overall. The commonest allele was BoLA-DRB3.2*16 (44.1%). In total, we found 40 genotypes. Considering the significant dominance of variant *16, as expected, 5 genotypes with its inclusion occurred: *16/*16, *12/*16, *06/*16, *16/*24 and *jba/*16. It was present in the genotype of two out three studied animals. Parameters of heterozygosity, effective number of alleles, Shannon and Pielou indices indicate that Ukrainian Gray cattle are characterized by lowest level of genetic variability and biodiversity according to the BoLA-DRB3 gene compared with other breeds. Due to significant dominance of allele *16, the breed has no inbred motifs. We noted deviation toward increase in homozygosity without deviations from the norm of the distribution according to Hardy-Weinberg equilibrium. The obtained results will be used for genetic-populational programs with the purpose of improving the genetic potential of cattle breeds in terms of economically beneficial traits and diseases of cattle.

Population genetic structure and association mapping for iron toxicity tolerance in rice

Iron (Fe) toxicity is a major abiotic stress which severely reduces rice yield in many countries of the world. Genetic variation for this stress tolerance exists in rice germplasms. Mapping of gene(s)/QTL controlling the stress tolerance and transfer of the traits into high yielding rice varieties are essential for improvement against the stress. A panel population of 119 genotypes from 352 germplasm lines was constituted for detecting the candidate gene(s)/QTL through association mapping. STRUCTURE, GenAlEx and Darwin softwares were used to classify the population. The marker-trait association was detected by considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) analyses. Wide genetic variation was observed among the genotypes present in the panel population for the stress tolerance. Linkage disequilibrium was detected in the population for iron toxicity tolerance. The population was categorized into three genetic structure groups. Marker-trait association study considering both the Generalized Linear Model (GLM) and Mixed Linear Model (MLM) showed significant association of leaf browning index (LBI) with markers RM471, RM3, RM590 and RM243. Three novel QTL controlling Fe-toxicity tolerance were detected and designated as qFeTox4.3, qFeTox6.1 and qFeTox10.1. A QTL reported earlier in the marker interval of C955-C885 on chromosome 1 is validated using this panel population. The present study showed that QTL controlling Fe-toxicity tolerance to be co-localized with the QTL for Fe-biofortification of rice grain indicating involvement of common pathway for Fe toxicity tolerance and Fe content in rice grain. Fe-toxicity tolerance QTL qFeTox6.1 was co-localized with grain Fe-biofortification QTLs qFe6.1 and qFe6.2 on chromosome 6, whereas qFeTox10.1 was co-localized with qFe10.1 on chromosome 10. The Fe-toxicity tolerance QTL detected from this mapping study will be useful in marker-assisted breeding programs.

Generation of Chloroplast Molecular Markers to Differentiate Sophora toromiro and Its Hybrids as a First Approach to Its Reintroduction in Rapa Nui (Easter Island)

Sophora toromiro is an endemic tree of Rapa Nui with religious and cultural relevance that despite being extinct in the wild, still persists in botanical gardens and private collections around the world. The authenticity of some toromiro trees has been questioned because the similarities among hybrid lines leads to misclassification of the species. The conservation program of toromiro has the objective of its reinsertion into Rapa Nui, but it requires the exact genotyping and certification of the selected plants in order to efficiently reintroduce the species. In this study, we present for the first time the complete chloroplast genome of S. toromiro and four other Sophora specimens, which were sequenced de-novo and assembled after mapping the raw reads to a chloroplast database. The length of the chloroplast genomes ranges from 154,239 to 154,473 bp. A total of 130-143 simple sequence repeats (SSR) loci and 577 single nucleotide polymorphisms (SNPs) were identified.

Genetic Diversity and Genome-Wide Association Study of Seed Aspect Ratio Using a High-Density SNP Array in Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is one of the important oil crops of the world. In this study, we aimed to evaluate the genetic diversity of 384 peanut germplasms including 100 Korean germplasms and 284 core collections from the United States Department of Agriculture (USDA) using an Axiom_Arachis array with 58K single-nucleotide polymorphisms (SNPs). We evaluated the evolutionary relationships among 384 peanut germplasms using a genome-wide association study (GWAS) of seed aspect ratio data processed by ImageJ software. In total, 14,030 filtered polymorphic SNPs were identified from the peanut 58K SNP array. We identified five SNPs with significant associations to seed aspect ratio on chromosomes Aradu.A09, Aradu.A10, Araip.B08, and Araip.B09. AX-177640219 on chromosome Araip.B08 was the most significantly associated marker in GAPIT and Regularization method. Phosphoenolpyruvate carboxylase (PEPC) was found among the eleven genes within a linkage disequilibrium (LD) of the significant SNPs on Araip.B08 and could have a strong causal effect in determining seed aspect ratio. The results of the present study provide information and methods that are useful for further genetic and genomic studies as well as molecular breeding programs in peanuts.

New Transcriptome-Based SNP Markers for Noug (Guizotia abyssinica) and Their Conversion to KASP Markers for Population Genetics Analyses

The development and use of genomic resources are essential for understanding the population genetics of crops for their efficient conservation and enhancement. Noug (Guizotia abyssinica) is an economically important oilseed crop in Ethiopia and India. The present study sought to develop new DNA markers for this crop. Transcriptome sequencing was conducted on two genotypes and 628 transcript sequences containing 959 single nucleotide polymorphisms (SNPs) were developed. A competitive allele-specific PCR (KASP) assay was developed for the SNPs and used for genotyping of 24 accessions. A total of 554 loci were successfully genotyped across the accessions, and 202 polymorphic loci were used for population genetics analyses. Polymorphism information content (PIC) of the loci varied from 0.01 to 0.37 with a mean of 0.24, and about 49% of the loci showed significant deviation from the Hardy-Weinberg equilibrium. The mean expected heterozygosity was 0.27 suggesting moderately high genetic variation within accessions. Low but significant differentiation existed among accessions (FST = 0.045, p < 0.0001). Landrace populations from isolated areas may have useful mutations and should be conserved and used in breeding this crop. The genomic resources developed in this study were shown to be useful for population genetics research and can also be used in, e.g., association genetics.

Application of Genomic Big Data in Plant Breeding:Past, Present, and Future

Plant breeding has a long history of developing new varieties that have ensured the food security of the human population. During this long journey together with humanity, plant breeders have successfully integrated the latest innovations in science and technologies to accelerate the increase in crop production and quality. For the past two decades, since the completion of human genome sequencing, genomic tools and sequencing technologies have advanced remarkably, and adopting these innovations has enabled us to cost down and/or speed up the plant breeding process. Currently, with the growing mass of genomic data and digitalized biological data, interdisciplinary approaches using new technologies could lead to a new paradigm of plant breeding. In this review, we summarize the overall history and advances of plant breeding, which have been aided by plant genomic research. We highlight the key advances in the field of plant genomics that have impacted plant breeding over the past decades and introduce the current status of innovative approaches such as genomic selection, which could overcome limitations of conventional breeding and enhance the rate of genetic gain.

Genetic Architecture and Anthocyanin Profiling of Aromatic Rice From Manipur Reveals Divergence of Chakhao Landraces

Aromatic rice of Manipur popularly known as Chakhao is a speciality glutinous rice, for which protection under geographical indication in India has been granted recently. The agronomic and nutraceutical variability of the Chakhao rice germplasm is yet to be genetically characterized. To address this gap, characterization of ninety-three landraces for agro-morphological traits, grain pigmentation, antioxidant properties, and molecular genetic variation was carried out to unravel their population genetic structure. Two major groups were identified based on pericarp color, namely, purple and non-purple, which showed a significant variation for plant height, panicle length, and grain yield. Molecular marker analysis revealed three subpopulations that could be associated with pericarp pigmentation. Deep purple genotypes formed POP3, japonica genotypes adapted to hill environment formed POP1, while POP2 comprised of both indica and aus types. Liquid chromatography-mass spectrometry (LC-MS) analysis revealed two major anthocyanin compounds in pigmented rices, namely, cyanidin-3-O-glucoside (C3G) and peonidin-3-O-glucoside (P3G). The total anthocyanin content among pigmented genotypes ranged from 29.8 to 275.8 mg.100g-1 DW. Total phenolics ranged from 66.5 to 700.3 mg GAE.100g-1 DW with radical scavenging activity (RSA) varying between 17.7 and 65.7%. Anthocyanins and phenolics showed a direct relationship with RSA implying the nutraceutical benefits of deep pigmented rice such as Manipur black rice. Aromatic rices from Manipur were found to be genetically diverse. Therefore, efforts need to be made for maintaining the geographic identity of these rice and utilization in breeding for region-specific cultivar improvement.

Association of molecular markers with physio-biochemical traits related to seed vigour in rice

Eighteen physio-biochemical traits influencing seed vigour were studied for their association with molecular markers using a mini core set constituted from 120 germplasm lines. High genetic variation was detected in the parameters namely chlrophyll a, Chlrophyll b, total chlorophyll, carotenoids, total anthocyanin content, gamma-oryzanols, total phenolics content, superoxide dismutase, catalase, guaicol peroxidase, total soluble sugar, total protein, seed vigour index -I and seed vigour index -II. Strong positive correlation of seed vigour index II was observed with amylose content, total anthocyanin content, catalase, total phenolic content and total flavonoid content while a negative association was observed for gamma-oryzanol content. High gene diversity (0.7169) and informative markers value (0.6789) were estimated from the investigation. Three genetic structure groups were observed in the panel population and genotypes were grouped in the subpopulations based on the seed vigour trait. Differences in the fixation indices of the three sub populations indicated existence of linkage disequilibrium in the studied panel population. Association of the traits namely total flavonoids, superoxide dismutase, catalase, chlorophyll a, Chlorophyll b, total chlorophyll, carotenoids, starch, amylose, total anthocyanin, gamma-oryzanol, total phenolics with the molecular markers were detected by Generalized Linear Model and Mixed Linear Model showing > 0.10 R² value. Association of the trait, total flavonoids with marker RM7364 located on chromosome 8 reported in earlier study was validated in this investigation. The validated markers and the novel markers detected showing higher R² value will be useful for improvement of seed vigour in rice.

Tuberculosis causes highly conserved metabolic changes in human patients, mycobacteria-infected mice and zebrafish larvae

Tuberculosis is a highly infectious and potentially fatal disease accompanied by wasting symptoms, which cause severe metabolic changes in infected people. In this study we have compared the effect of mycobacteria infection on the level of metabolites in blood of humans and mice and whole zebrafish larvae using one highly standardized mass spectrometry pipeline, ensuring technical comparability of the results. Quantification of a range of circulating small amines showed that the levels of the majority of these compounds were significantly decreased in all three groups of infected organisms. Ten of these metabolites were common between the three different organisms comprising: methionine, asparagine, cysteine, threonine, serine, tryptophan, leucine, citrulline, ethanolamine and phenylalanine. The metabolomic changes of zebrafish larvae after infection were confirmed by nuclear magnetic resonance spectroscopy. Our study identified common biomarkers for tuberculosis disease in humans, mice and zebrafish, showing across species conservation of metabolic reprogramming processes as a result of disease. Apparently, the mechanisms underlying these processes are independent of environmental, developmental and vertebrate evolutionary factors. The zebrafish larval model is highly suited to further investigate the mechanism of metabolic reprogramming and the connection with wasting syndrome due to infection by mycobacteria.

Genetic diversity, SNP-trait associations and genomic selection accuracy in a west African collection of Kersting's groundnut [Macrotyloma geocarpum(Harms) Maréchal & Baudet]

Understanding the mechanisms governing complex traits variation is a requirement for efficient crop improvement. In this study, the molecular characterization, marker-trait associations and the possibility for genomic selection in a collection of 281 Kersting's groundnut accessions were carried out. The diversity panel was phenotyped using an Alpha lattice design with two replicates in two contrasting environments. Accessions were genotyped using genotyping by sequencing technology. Genome-wide association analyses were performed between single nucleotide polymorphism markers and yield-related traits across tested environments. SNP markers were used to calculate the observed (Ho) and expected heterozygosity (He), and the total gene diversity (Ht). Genetic differentiation among accessions across ecological regions of origin was analysed. Our results revealed 493 quality SNPs of which 113 had a minor allele frequency>0.05, a total gene diversity of 0.43 and average Ho and He values of 0.04 and 0.22, respectively. Four clusters, highly differentiated by seed coat colour (Fst = 0.79), were identified. The population structure analysis showed two subpopulations with high differentiation across ecological regions (Fst = 0.37). The GWAS revealed 10 significant marker-trait associations, of which six SNPs were consistent across environments. The genomic selection through cross-validation showed moderate to high prediction accuracies for leaflet length, seed dimension traits, 100 seed weight, days to 50% flowering and days to maturity. This demonstrates the existence of genetic variability within Kersting's groundnut and shows the potential for the improvement of the species. The findings also provide a first insight into the phenotype-to-genotype relationships in Kersting's groundnut, using SNP markers.

Association mapping for yield traits in Paeonia rockii based on SSR markers within transcription factors of comparative transcriptome

BACKGROUND

Allelic variation underlying the quantitative traits in plants is caused by the extremely complex regulation process. Tree peony originated in China is a peculiar ornamental, medicinal and oil woody plant. Paeonia rockii, one of tree peony species, is a precious emerging woody oil crop. However, in this valuable plant, the study of functional loci associated with yield traits has rarely been identified. Therefore, to explore the genetic architecture of 24 yield quantitative traits, the association mapping was first reported in 420 unrelated cultivated P. rockii individuals based on the next-generation sequencing (NGS) and single-molecule long-read sequencing (SMLRS).

RESULTS

The developed 58 pairs of polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers from 959 candidate transcription factors (TFs) associated with yield were used for genotyping the 420 P. rockii accessions. We observed a high level of genetic diversity (polymorphic information content, PIC = 0.514) and low linkage disequilibrium (LD) between EST-SSRs. Moreover, four subpopulations in the association population were revealed by STRUCTURE analyses. Further, single-marker association analysis identified 141 significant associations, involving 17 quantitative traits and 41 EST-SSRs. These loci were mainly from AP2, TCP, MYB, HSF, bHLH, GATA, and B3 gene families and showed a small proportion of the phenotypic variance (3.79 to 37.45%).

CONCLUSIONS

Our results summarize a valuable collection of functional loci associated with yield traits in P. rockii, and provide a precious resource that reveals allelic variation underlying quantitative traits in Paeonia and other woody oil crops.

Allelic sequence variation in the Sub1A, Sub1B and Sub1C genes among diverse rice cultivars and its association with submergence tolerance

Erratic rainfall leading to flash flooding causes huge yield losses in lowland rice. The traditional varieties and landraces of rice possess variable levels of tolerance to submergence stress, but gene discovery and utilization of these resources has been limited to the Sub1A-1 allele from variety FR13A. Therefore, we analysed the allelic sequence variation in three Sub1 genes in a panel of 179 rice genotypes and its association with submergence tolerance. Population structure and diversity analysis based on a 36-plex genome wide genic-SNP assay grouped these genotypes into two major categories representing Indica and Japonica cultivar groups with further sub-groupings into Indica, Aus, Deepwater and Aromatic-Japonica cultivars. Targetted re-sequencing of the Sub1A, Sub1B and Sub1C genes identfied 7, 7 and 38 SNPs making 8, 9 and 67 SNP haplotypes, respectively. Haplotype networks and phylogenic analysis revealed evolution of Sub1B and Sub1A genes by tandem duplication and divergence of the ancestral Sub1C gene in that order. The alleles of Sub1 genes in tolerant reference variety FR13A seem to have evolved most recently. However, no consistent association could be found between the Sub1 allelic variation and submergence tolerance probably due to low minor allele frequencies and presence of exceptions to the known Sub1A-1 association in the genotype panel. We identified 18 cultivars with non-Sub1A-1 source of submergence tolerance which after further mapping and validation in bi-parental populations will be useful for development of superior flood tolerant rice cultivars.

De novo Transcriptomic Resources in the Brain of Vespa velutina for Invasion Control

The invasion of Vespa velutina presents a great threat to the agriculture economy, the ecological environment, and human health. An effective strategy for this hornet control is urgently required, but the limited genome information of Vespa velutina restricts the application of molecular-genomic tools for targeted hornet management. Therefore, we conducted large-scale transcriptome profiling of the hornet brain to obtain functional target genes and molecular markers. Using an Illumina HiSeq platform, more than 41 million clean reads were obtained and de novo assembled into 182,087 meaningful unigenes. A total of 56,400 unigenes were annotated against publicly available protein sequence databases and a set of reliable Simple Sequence Repeats (SSRs) and Single Nucleotide Polymorphisms (SNP) markers were developed. The homologous genes encoding crucial behavior regulation factors, odorant binding proteins (OBPs), and vitellogenin, were also identified from highly expressed transcripts. This study provides abundant molecular targets and markers for invasive hornet control and further promotes the genetic and molecular study of Vespa velutina.

Determining the Polymorphism Information Content of a molecular marker

This review was carried out with the purpose of contributing to the discussion on the equations used in calculating the Polymorphism Information Content (PIC) of molecular markers. PIC measures the ability of a marker to detect polymorphisms, and therefore has enormous importance in selecting markers for genetic studies. We perform a summary analysis of PIC and its difference in relation to heterozygosity, another parameter used to evaluate the quality of a marker, presenting and discussing the several equations registered in the literature for both dominant and codominant markers. Finally, we present a brief direction on estimating PIC for dominant markers.

Functional genomic approaches to improve crop plant heat stress tolerance

Heat stress as a yield limiting issue has become a major threat for food security as global warming progresses. Being sessile, plants cannot avoid heat stress. They respond to heat stress by activating complex molecular networks, such as signal transduction, metabolite production and expressions of heat stress-associated genes. Some plants have developed an intricate signalling network to respond and adapt it. Heat stress tolerance is a polygenic trait, which is regulated by various genes, transcriptional factors, proteins and hormones. Therefore, to improve heat stress tolerance, a sound knowledge of various mechanisms involved in the response to heat stress is required. The classical breeding methods employed to enhance heat stress tolerance has had limited success. In this era of genomics, next generation sequencing techniques, availability of genome sequences and advanced biotechnological tools open several windows of opportunities to improve heat stress tolerance in crop plants. This review discusses the potential of various functional genomic approaches, such as genome wide association studies, microarray, and suppression subtractive hybridization, in the process of discovering novel genes related to heat stress, and their functional validation using both reverse and forward genetic approaches. This review also discusses how these functionally validated genes can be used to improve heat stress tolerance through plant breeding, transgenics and genome editing approaches.

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 assessment of three Fagaceae species in forest restoration trials

Restoring isolated patches of forest ecosystems in degraded landscapes could potentially lead to genetic loss and inbreeding. Therefore, this study determined the occurrence of genetic diversity among the tree species Castanopsis tribuloides, C. calathiformis, and Lithocarpus polystachyus all of which were proven previously to be effective native tree species in the restoration of upland evergreen forests in northern Thailand when using the seed sample collection method. We tested our hypothesis as to whether the genetic diversity of a plant population that had been planted from the seeds of 4–6 adult trees would be lower and whether incidences of fixation index (Fis) would be higher among the second generation seedlings of these three Fagaceae species in isolated forest restoration trial plots. Microsatellite primers were selected from the entire genome sequence of C. tribuloides and the genetic sequences of C. tribuloides, L. polystachyus, and C. calathiformis were analyzed. Our results indicated a high degree of genetic diversity (He) in C. tribuloides (0.736) and C. calathiformis (0.481); however, a low level of genetic diversity was observed in L. polystachyus (0.281) within the restored forest. The fixation index for the second generation of L. polystachyus and C. calathiformis in the restored forest showed evidence of inbreeding. These results imply the efficiency of the seed sample collection method and verify that it does not reduce the level of genetic diversity in C. tribuloides and C. calathiformis. However, it may result in incidences of an inbreeding phenomena, suggesting the need to increase the number of adult trees used at the seed collection stage.