Applications of single nucleotide polymorphisms in crop genetics

Genomic prediction in Brassica napus: evaluating the benefit of imputed whole-genome sequencing data

Advances in sequencing technology allow whole plant genomes to be sequenced with high quality. Combining genotypic and phenotypic data in genomic prediction helps breeders to select crossing partners in partially phenotyped populations. In plant breeding programs, the cost of sequencing entire breeding populations still exceeds available genotyping budgets. Hence, the method for genotyping is still mainly single nucleotide polymorphism (SNP) arrays; however, arrays are unable to assess the entire genome- and population-wide diversity. A compromise involves genotyping the entire population using an SNP array and a subset of the population with whole-genome sequencing. Both datasets can then be used to impute markers from whole-genome sequencing onto the entire population. Here, we evaluate whether imputation of whole-genome sequencing data enhances genomic predictions, using data from a nested association mapping population of rapeseed (Brassica napus). Employing two cross-validation schemes that mimic scenarios for the prediction of close and distant relatives, we show that imputed marker data do not significantly improve prediction accuracy, likely due to redundancy in relationship estimates and imputation errors. In simulation studies, only small improvements were observed, further corroborating the findings. We conclude that SNP arrays are already equipped with the information that is added by imputation through relationship and linkage disequilibrium.

Modern Breeding Strategies and Tools for Durable Late Blight Resistance in Potato

Cultivated potato (Solanum tuberosum) is a major crop worldwide. It occupies the second place after cereals (corn, rice, and wheat). This important crop is threatened by the Oomycete Phytophthora infestans, the agent of late blight disease. This pathogen was first encountered during the Irish famine during the 1840s and is a reemerging threat to potatoes. It is mainly controlled chemically by using fungicides, but due to health and environmental concerns, the best alternative is resistance. When there is no disease, no treatment is required. In this study, we present a summary of the ongoing efforts concerning resistance breeding of potato against this devastating pathogen, P. infestans. This work begins with the search for and selection of resistance genes, whether they are from within or from outside the species. The genetic methods developed to date for gene mining, such as effectoromics and GWAS, provide researchers with the ability to identify genes of interest more efficiently. Once identified, these genes are cloned using molecular markers (MAS or QRL) and can then be introduced into different cultivars using somatic hybridization or recombinant DNA technology. More innovative technologies have been developed lately, such as gene editing using the CRISPR system or gene silencing, by exploiting iRNA strategies that have emerged as promising tools for managing Phytophthora infestans, which can be employed. Also, gene pyramiding or gene stacking, which involves the accumulation of two or more R genes on the same individual plant, is an innovative method that has yielded many promising results. All these advances related to the development of molecular techniques for obtaining new potato cultivars resistant to P. infestans can contribute not only to reducing losses in agriculture but especially to ensuring food security and safety.

Exploring Disease Resistance in Pepper (Capsicum spp.) Germplasm Collection Using Fluidigm SNP Genotyping

This study utilized a diverse Capsicum accessions (5658) sourced from various species and geographical regions, deposited at the National Agrobiodiversity Center, Genebank. We employed 19 SNP markers through a Fluidigm genotyping system and screened these accessions against eight prevalent diseases of pepper. This study revealed accessions resistant to individual diseases as well as those exhibiting resistance to multiple diseases, including bacterial spot, anthracnose, powdery mildew, phytophthora root rot, and potyvirus. The C. chacoense accessions were identified as resistant materials against bacterial spot, anthracnose, powdery mildew, and phytophthora root rot, underscoring the robust natural defense mechanisms inherent in the wild Capsicum species and its potential uses as sources of resistance for breeding. C. baccatum species also demonstrated to be a promising source of resistance to major pepper diseases. Generally, disease-resistant germplasm has been identified from various Capsicum species. Originating from diverse locations such as Argentina, Bolivia, and the United Kingdom, these accessions consistently demonstrated resistance, indicating the widespread prevalence of disease-resistant traits across varied environments. Additionally, we selected ten pepper accessions based on their resistance to multiple diseases, including CMV, Phytophthora root rot, potyviruses, and TSWV, sourced from diverse geographical regions like Hungary, Peru, the United States, and the Netherlands. This comprehensive analysis provides valuable insights into disease resistance in Capsicum, crucial for fostering sustainable agricultural practices and advancing crop improvement through breeding strategies.

Development of Novel KASP Markers for Improved Germination in Deep-Sown Direct Seeded Rice

BACKGROUND

The lack of stable-high yielding and direct-seeded adapted varieties with better germination ability from deeper soil depth and availability of molecular markers are major limitation in achieving the maximum yield potential of rice under water and resource limited conditions. Development of high-throughput and trait-linked markers are of great interest in genomics-assisted breeding. The aim of present study was to develop and validate novel KASP (Kompetitive Allele-Specific PCR) markers associated with traits improving germination and seedling vigor of deep sown direct seeded rice (DSR).

RESULTS

Out of 58 designed KASP assays, four KASP assays did not show any polymorphism in any of the eleven genetic backgrounds considered in the present study. The 54 polymorphic KASP assays were then validated for their robustness and reliability on the F1s plants developed from eight different crosses considered in the present study. The third next validation was carried out on 256 F3:F4 and 713 BC3F2:3 progenies. Finally, the reliability of the KASP assays was accessed on a set of random 50 samples from F3:F4 and 80-100 samples from BC3F2:3 progenies using the 10 random markers. From the 54 polymorphic KASP, based on the false positive rate, false negative rate, KASP utility in different genetic backgrounds and significant differences in the phenotypic values of the positive (desirable) and negative (undesirable) traits, a total of 12 KASP assays have been selected. These 12 KASP include 5 KASP on chromosome 3, 1 on chromosome 4, 3 on chromosome 7 and 3 on chromosome 8. The two SNPs lying in the exon regions of LOC_Os04g34290 and LOC_Os08g32100 led to non-synonymous mutations indicating a possible deleterious effect of the SNP variants on the protein structure.

CONCLUSION

The present research work will provide trait-linked KASP assays, improved breeding material possessing favourable alleles and breeding material in form of expected pre-direct-seeded adapted rice varieties. The marker can be utilized in introgression program during pyramiding of valuable QTLs/genes providing adaptation to rice under DSR. The functional studies of the genes LOC_Os04g34290 and LOC_Os08g32100 possessing two validated SNPs may provide valuable information about these genes.

Genome-wide QTL mapping for agronomic traits in the winter wheat cultivar Pindong 34 based on 90K SNP array

Introduction

Agronomic traits are key components of wheat yield. Exploitation of the major underlying quantitative trait loci (QTLs) can improve the yield potential in wheat breeding.

Methods

In this study, we constructed a recombinant inbred line (RIL) population from Mingxian 169 (MX169) and Pindong 34 (PD34) to determine the QTLs for grain length (GL), grain width (GW), grain length-to-width ratio (LWR), plant height (PH), spike length (SL), grain number per spike (GNS), and the thousand grain weight (TGW) across four environments using wheat 90K SNP array.

Results

A QTL associated with TGW, i.e., QTGWpd.swust-6BS, was identified on chromosome 6B, which explained approximately 14.1%-16.2% of the phenotypic variation. In addition, eight QTLs associated with GL were detected across six chromosomes in four different test environments. These were QGLpd.swust-1BL, QGLpd.swust-2BL, QGLpd.swust-3BL.1, QGLpd.swust-3BL.2, QGLpd.swust-5DL, QGLpd.swust-6AL, QGLpd.swust-6DL.1, and QGLpd.swust-6DL.2. They accounted for 9.0%-21.3% of the phenotypic variation. Two QTLs, namely, QGWpd.swust-3BS and QGWpd.swust-6DL, were detected for GW on chromosomes 3B and 6D, respectively. These QTLs explained 12.8%-14.6% and 10.8%-15.2% of the phenotypic variation, respectively. In addition, two QTLs, i.e., QLWRpd.swust-7AS.1 and QLWRpd.swust-7AS.2, were detected on chromosome 7A for the grain LWR, which explained 10.9%-11.6% and 11.6%-11.2% of the phenotypic variation, respectively. Another QTL, named QGNSpd-swust-6DS, was discovered on chromosome 6D, which determines the GNS and which accounted for 11.4%-13.8% of the phenotypic variation. Furthermore, five QTLs associated with PH were mapped on chromosomes 2D, 3A, 5A, 6B, and 7B. These QTLs were QPHpd.swust-2DL, QPHpd.swust-3AL, QPHpd.swust-5AL, QPHpd.swust-6BL, and QPHpd.swust-7BS, which accounted for 11.3%-19.3% of the phenotypic variation. Lastly, a QTL named QSLpd.swust-3AL, conferring SL, was detected on chromosome 3A and explained 16.1%-17.6% of the phenotypic variation. All of these QTLs were defined within the physical interval of the Chinese spring reference genome.

Discussion

The findings of this study have significant implications for the development of fine genetic maps, for genomic breeding, and for marker-assisted selection to enhance wheat grain yield.

Computational Analysis of Single Nucleotide Polymorphisms in Human HIC1 Gene

Background A putative tumor suppressor gene called HIC1 (hypermethylated in cancer) is situated at 17p13.3, a locus where the allelic loss occurs often in human malignancies, including breast cancer. Hypermethylated in cancer 1 protein is a protein that in humans is encoded by the HIC1 gene and it's a Homo sapiens (Human). This gene functions as a growth regulatory and tumor repressor gene. The molecular function of HIC1 gene includes DNA-binding transcription factor activity, sequence-specific DNA binding, DNA binding, histone deacetylase binding, protein binding, metal ion binding, nucleic acid binding, DNA-binding transcription repressor activity, RNA polymerase II-specific, DNA-binding transcription factor activity, RNA polymerase II-specific. The biological process of HIC1 gene includes multicellular organism development, negative regulation of Wnt signaling pathway, positive regulation of DNA damage response, signal transduction by p53 class mediator regulation of transcription, DNA-templated, negative regulation of transcription by RNA polymerase II, Wnt signaling pathway, transcription, DNA-templated, intrinsic apoptotic signaling pathway in response to DNA damage, cellular response to DNA damage stimulus. The study aimed to predict the stability and structure of the protein that will arise from single nucleotide polymorphisms (SNPs) in the human HIC1 gene. Methodology To investigate the possible negative effects associated with these SNPs, bioinformatic analysis is typically essential. The following tools were employed for forecasting harmful SNPs: scale-invariant feature transform (SIFT), Protein Analysis Through Evolutionary Relationships (PANTHER), nonsynonymous SNP by Protein Variation Effect Analyzer (PROVEAN), and nonsynonymous SNP by Single Nucleotide Polymorphism Annotation Platform (SNAP). Results The present study identified a total of 36 SNPs using the SIFT approach, which were shown to have functional significance. Twenty-six were determined to be tolerable, whereas 10 were shown to be detrimental. Out of 20 SNPs, seven (P370A, P646S, R654P, A476T, S400S, D666N, D7V) SNPs were predicted as "Possibly damaging" and seven (L9F, G468R, G490R, L482R, S12W, G489D, S12P) were identified as "probably benign", and six (R725G, G620S, A56V, E463D, D394N, L338V) were identified as "probably damaging" according to the predictions made by PANTHER tools. The majority of the pixels on the strip were red, indicating that the gene changes may have dangerous consequences. These results highlight the need for more research to fully comprehend how these mutations affect the hic1 protein's function, which is essential for the emergence of different types of cancer. Conclusion The current research has provided us with essential information about how SNPs might be used as a diagnostic marker for cancer, given that SNPs may be candidates for cellular changes caused by mutations linked to cancer.

Genome-wide genetic architecture for plant maturity and drought tolerance in diploid potatoes

Cultivated potato (Solanum tuberosum) is known to be highly susceptible to drought. With climate change and its frequent episodes of drought, potato growers will face increased challenges to achieving their yield goals. Currently, a high proportion of untapped potato germplasm remains within the diploid potato relatives, and the genetic architecture of the drought tolerance and maturity traits of diploid potatoes is still unknown. As such, a panel of 384 ethyl methanesulfonate-mutagenized diploid potato clones were evaluated for drought tolerance and plant maturity under field conditions. Genome-wide association studies (GWAS) were conducted to dissect the genetic architecture of the traits. The results obtained from the genetic structure analysis of the panel showed five main groups and seven subgroups. Using the Genome Association and Prediction Integrated Tool-mixed linear model GWAS statistical model, 34 and 17 significant quantitative trait nucleotides (QTNs) were found associated with maturity and drought traits, respectively. Chromosome 5 carried most of the QTNs, some of which were also detected by using the restricted two-stage multi-locus multi-allele-GWAS haploblock-based model, and two QTNs were found to be pleiotropic for both maturity and drought traits. Using the non-parametric U-test, one and three QTNs, with 5.13%-7.4% phenotypic variations explained, showed favorable allelic effects that increase the maturity and drought trait values. The quantitaive trait loci (QTLs)/QTNs associated with maturity and drought trait were found co-located in narrow (0.5-1 kb) genomic regions with 56 candidate genes playing roles in plant development and senescence and in abiotic stress responses. A total of 127 potato clones were found to be late maturing and tolerant to drought, while nine were early to moderate-late maturing and tolerant to drought. Taken together, the data show that the studied germplasm panel and the identified candidate genes are prime genetic resources for breeders and biologists in conventional breeding and targeted gene editing as climate adaptation tools.

Using plants in forensics: State-of-the-art and prospects

The increasing use of plant evidence in forensic investigations gave rise to a powerful new discipline - Forensic Botany - that analyses micro- or macroscopic plant materials, such as the totality or fragments of an organ (i.e., leaves, stems, seeds, fruits, roots) and tissue (i.e., pollen grains, spores, fibers, cork) or its chemical composition (i. e., secondary metabolites, isotopes, DNA, starch grains). Forensic botanists frequently use microscopy, chemical analysis, and botanical expertise to identify and interpret evidence crucial to solving civil and criminal issues, collaborating in enforcing laws or regulations, and ensuring public health safeguards. The present work comprehensively examines the current state and future potential of Forensic Botany. The first section conveys the critical steps of plant evidence collection, documentation, and preservation, emphasizing the importance of these initial steps in maintaining the integrity of the items. It explores the different molecular analyses, covering the identification of plant species and varieties or cultivars, and discusses the limitations and challenges of these techniques in forensics. The subsequent section covers the diversity of Forensic Botany approaches, examining how plant evidence exposes food and pharmaceutical frauds, uncovers insufficient or erroneous labeling, traces illegal drug trafficking routes, and combats the illegal collection or trade of protected species and derivatives. National and global security issues, including the implications of biological warfare, bioterrorism, and biocrime are addressed, and a review of the contributions of plant evidence in crime scene investigations is provided, synthesizing a comprehensive overview of the diverse facets of Forensic Botany.

Accurate prediction of quantitative traits with failed SNP calls in canola and maize

In modern plant breeding, genomic selection is becoming the gold standard to select superior genotypes in large breeding populations that are only partially phenotyped. Many breeding programs commonly rely on single-nucleotide polymorphism (SNP) markers to capture genome-wide data for selection candidates. For this purpose, SNP arrays with moderate to high marker density represent a robust and cost-effective tool to generate reproducible, easy-to-handle, high-throughput genotype data from large-scale breeding populations. However, SNP arrays are prone to technical errors that lead to failed allele calls. To overcome this problem, failed calls are often imputed, based on the assumption that failed SNP calls are purely technical. However, this ignores the biological causes for failed calls-for example: deletions-and there is increasing evidence that gene presence-absence and other kinds of genome structural variants can play a role in phenotypic expression. Because deletions are frequently not in linkage disequilibrium with their flanking SNPs, permutation of missing SNP calls can potentially obscure valuable marker-trait associations. In this study, we analyze published datasets for canola and maize using four parametric and two machine learning models and demonstrate that failed allele calls in genomic prediction are highly predictive for important agronomic traits. We present two statistical pipelines, based on population structure and linkage disequilibrium, that enable the filtering of failed SNP calls that are likely caused by biological reasons. For the population and trait examined, prediction accuracy based on these filtered failed allele calls was competitive to standard SNP-based prediction, underlying the potential value of missing data in genomic prediction approaches. The combination of SNPs with all failed allele calls or the filtered allele calls did not outperform predictions with only SNP-based prediction due to redundancy in genomic relationship estimates.

Genetic diversity and population structure of Polygonatum cyrtonema Hua in China using SSR markers

Polygonatum cyrtonema Hua is a perennial herbaceous plant of the Polygonatum genus, belonging to the Liliaceae family, with significant medicinal and nutritional value. In China, this species is a traditional medicinal and edible herb with a long history of application and is widely appreciated by the people. However, as the demand for medicinal herbs continues to grow, excessive harvesting has led to the depletion of wild resources and the risk of genetic erosion. In addition, the chaotic cultivation of varieties and the lack of high quality germplasm resources have led to inconsistent quality of medical materials. Therefore, it is urgent to conduct genetic diversity evaluation of this species and establish a sound conservation plan. This study assessed the genetic diversity and population structure of 96 samples collected from seven regions in China using the simple sequence repeat (SSR) molecular marker technology. In this study, a total of 60 alleles (Na) were detected across the 10 polymorphic SSR markers used, with an average of 6.0 alleles generated per locus. The values of polymorphic information content (PIC) values ranged from 0.3396 to 0.8794, with an average value of 0.6430. The average value of the effective number of alleles (Ne) was 2.761, and the average value of the Shannon's information index (I) was 1.196. The population structure analysis indicates that the Polygonatum cyrtonema Hua germplasm can be classified into three subpopulations (JZ, QY, JD) at the molecular level, which corresponds to the previous subgroups identified based on individual plant phenotypic traits. Analysis of Molecular Variance (AMOVA) showed that 74% of the genetic variation was between individuals within populations in different regions. The phylogenetic analysis of the 96 germplasm samples divided them into three main populations. The QY and JD subpopulations are largely clustered together, which could be attributed to their mountainous distribution and the local climate environment. The genetic differentiation coefficient (Fst) value was low at 0.065, indicating relatively low population differentiation. The ratio of the genetic differentiation coefficient (Fst) between the JZ population and the other two populations (QY and JD) is much higher than the ratio between the QY and JD populations. Based on the clustering results and the ratio of the genetic differentiation coefficient (Fst), it can be inferred that the genetic relationship between the QY and JD subpopulations is closer, with a certain degree of genetic differentiation from the JZ subpopulation. This study supports the conservation of germplasm resources of Polygonatum cyrtonema Hua in China and provides new parental material for germplasm genetic improvement and breeding programs.

Genome-wide association study of soybean (Glycine max [L.] Merr.) germplasm for dissecting the quantitative trait nucleotides and candidate genes underlying yield-related traits

Soybean (Glycine max [L.] Merr.) is one of the most significant crops in the world in terms of oil and protein. Owing to the rising demand for soybean products, there is an increasing need for improved varieties for more productive farming. However, complex correlation patterns among quantitative traits along with genetic interactions pose a challenge for soybean breeding. Association studies play an important role in the identification of accession with useful alleles by locating genomic sites associated with the phenotype in germplasm collections. In the present study, a genome-wide association study was carried out for seven agronomic and yield-related traits. A field experiment was conducted in 2015/2016 at two locations that include 155 diverse soybean germplasm. These germplasms were genotyped using SoySNP50K Illumina Infinium Bead-Chip. A total of 51 markers were identified for node number, plant height, pods per plant, seeds per plant, seed weight per plant, hundred-grain weight, and total yield using a multi-locus linear mixed model (MLMM) in FarmCPU. Among these significant SNPs, 18 were putative novel QTNs, while 33 co-localized with previously reported QTLs. A total of 2,356 genes were found in 250 kb upstream and downstream of significant SNPs, of which 17 genes were functional and the rest were hypothetical proteins. These 17 candidate genes were located in the region of 14 QTNs, of which ss715580365, ss715608427, ss715632502, and ss715620131 are novel QTNs for PH, PPP, SDPP, and TY respectively. Four candidate genes, Glyma.01g199200, Glyma.10g065700, Glyma.18g297900, and Glyma.14g009900, were identified in the vicinity of these novel QTNs, which encode lsd one like 1, Ergosterol biosynthesis ERG4/ERG24 family, HEAT repeat-containing protein, and RbcX2, respectively. Although further experimental validation of these candidate genes is required, several appear to be involved in growth and developmental processes related to the respective agronomic traits when compared with their homologs in Arabidopsis thaliana. This study supports the usefulness of association studies and provides valuable data for functional markers and investigating candidate genes within a diverse germplasm collection in future breeding programs.

Novel SNP based analysis of genetic diversity in Polygonatum verticillatum Linn. across Indian Himalayas

Polygonatum verticillatum Linn. is an endangered medicinal herb from Himalayas whose rhizomes have recently been used to curate symptoms of COVID-19. During present investigation, a gene bank of P. verticillatum was established at High Altitude Herbal Garden of Forest Research Institute, Dehradun at Chakrata, at 2600 m amsl with germplasm collected from different states and union territory of India including Himachal Pradesh, Sikkim, Uttarakhand and Jammu and Kashmir covering a wide range of geographical locations from an altitude of 1800 to 3600 m amsl. Genotyping by sequencing was applied to a set of 66 accessions of P. verticillatum to identify genome-wide high quality single nucleotide polymorphisms (SNPs) for analysis of genetic diversity. Neighbour-joining tree created from the distance matrix data grouped the genotypes into five distinct clusters. The results of principal coordinate analysis and Cluster analysis overlapped to identify Narkanda, Shimla (Himachal Pradesh) and Sunil village, Chamoli (Uttarakhand) as the regions with undisturbed, highly diverse natural populations of P. verticillatum. The species displayed little congruence in terms of genetic similarities with altitudinal range. This investigation is first of its kind on generation and utilization of SNPs to analyse genetic diversity in P. verticillatum with a very vivid sample collection across the entire Himalayan range in India.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03654-4.

Genotyping of DNA pools identifies untapped landraces and genomic regions to develop next-generation varieties

Landraces, that is, traditional varieties, have a large diversity that is underexploited in modern breeding. A novel DNA pooling strategy was implemented to identify promising landraces and genomic regions to enlarge the genetic diversity of modern varieties. As proof of concept, DNA pools from 156 American and European maize landraces representing 2340 individuals were genotyped with an SNP array to assess their genome-wide diversity. They were compared to elite cultivars produced across the 20th century, represented by 327 inbred lines. Detection of selective footprints between landraces of different geographic origin identified genes involved in environmental adaptation (flowering times, growth) and tolerance to abiotic and biotic stress (drought, cold, salinity). Promising landraces were identified by developing two novel indicators that estimate their contribution to the genome of inbred lines: (i) a modified Roger's distance standardized by gene diversity and (ii) the assignation of lines to landraces using supervised analysis. It showed that most landraces do not have closely related lines and that only 10 landraces, including famous landraces as Reid's Yellow Dent, Lancaster Surecrop and Lacaune, cumulated half of the total contribution to inbred lines. Comparison of ancestral lines directly derived from landraces with lines from more advanced breeding cycles showed a decrease in the number of landraces with a large contribution. New inbred lines derived from landraces with limited contributions enriched more the haplotype diversity of reference inbred lines than those with a high contribution. Our approach opens an avenue for the identification of promising landraces for pre-breeding.

Genomic insights into genetic diversity and local adaptation of a dominant desert steppe feather grass, Stipa breviflora Griseb

Investigating the genetic mechanisms of local adaptation is critical to understanding how species adapt to heterogeneous environments. In the present study, we analyzed restriction site-associated DNA sequencing (RADseq) data in order to explore genetic diversity, genetic structure, genetic differentiation, and local adaptation of Stipa breviflora. In total, 135 individual plants were sequenced and 25,786 polymorphic loci were obtained. We found low genetic diversity (He = 0.1284) within populations of S. breviflora. Four genetic clusters were identified along its distribution range. The Mantel test, partial Mantel test, and multiple matrix regression with randomization (MMRR) indicate that population differentiation was caused by both geographic distance and environmental factors. Through the F_ST outlier test and environmental association analysis (EAA), 113 candidate loci were identified as putatively adaptive loci. RPK2 and CPRF1, which are associated with meristem maintenance and light responsiveness, respectively, were annotated. To explore the effects of climatic factors on genetic differentiation and local adaptation of S. breviflora, gradient forest (GF) analysis was applied to 25,786 single nucleotide polymorphisms (SNPs) and 113 candidate loci, respectively. The results showed that both temperature and precipitation affected the genetic differentiation of S. breviflora, and precipitation was strongly related to local adaptation. Our study provides a theoretical basis for understanding the local adaptation of S. breviflora.

Vitis vinifera genotyping toolbox to highlight diversity and germplasm identification

The contribution of vine cultivation to human welfare as well as the stimulation of basic social and cultural features of civilization has been great. The wide temporal and regional distribution created a wide array of genetic variants that have been used as propagating material to promote cultivation. Information on the origin and relationships among cultivars is of great interest from a phylogenetics and biotechnology perspective. Fingerprinting and exploration of the complicated genetic background of varieties may contribute to future breeding programs. In this review, we present the most frequently used molecular markers, which have been used on Vitis germplasm. We discuss the scientific progress that led to the new strategies being implemented utilizing state-of-the-art next generation sequencing technologies. Additionally, we attempted to delimit the discussion on the algorithms used in phylogenetic analyses and differentiation of grape varieties. Lastly, the contribution of epigenetics is highlighted to tackle future roadmaps for breeding and exploitation of Vitis germplasm. The latter will remain in the top of the edge for future breeding and cultivation and the molecular tools presented herein, will serve as a reference point in the challenging years to come.

Novel QTL for chilling tolerance at germination and early seedling stages in sorghum

Sorghum (Sorghum bicolor L.) a drought tolerant staple crop for half a billion people in Africa and Asia, an important source of animal feed throughout the world and a biofuel feedstock of growing importanceorghum’s originated from tropical regions rendering the crop to be cold sensitive. Low temperature stresses such as chilling and frost greatly affect the agronomic performance of sorghum and limit its geographical distribution, posing a major problem in temperate environments when sorghum is planted early. Understanding the genetic basis of wide adaptability and of sorghum would facilitate molecular breeding programs and studies of other C4 crops. The objective of this study is to conduct quantitative trait loci analysis using genotying by sequencing for early seed germination and seedling cold tolerance in two sorghum recombinant inbred lines populations. To accomplish that, we used two populations of recombinant inbred lines (RIL) developed from crosses between cold-tolerant (CT19, ICSV700) and cold-sensitive (TX430, M81E) parents. The derived RIL populations were evaluated for single nucleotide polymorphism (SNP) using genotype-by-sequencing (GBS) in the field and under controlled environments for their response to chilling stress. Linkage maps were constructed with 464 and 875 SNPs for the CT19 X TX430 (C1) and ICSV700 X M81 E (C2) populations respectively. Using quantitative trait loci (QTL) mapping, we identified QTL conferring tolerance to chilling tolerance at the seedling stage. A total of 16 and 39 total QTL were identified in the C1 and C2 populations, respectively. Two major QTL were identified in the C1 population, and three major QTL were mapped in the C2 population. Comparisons between the two populations and with previously identified QTL show a high degree of similarity in QTL locations. Given the amount of co-localization of QTL across traits and the direction of allelic effect supports that these regions have a pleiotropic effect. These QTL regions were also identified to be highly enriched for genes encoding chilling stress and hormonal response genes. This identified QTL can be useful in developing tools for molecular breeding of sorghums with improved low-temperature germinability.

Marker-Assisted Pyramiding of Blast-Resistance Genes in a japonica Elite Rice Cultivar through Forward and Background Selection

Rice blast, caused by Pyricularia oryzae, is one of the main rice diseases worldwide. The pyramiding of blast-resistance (Pi) genes, coupled to Marker-Assisted BackCrossing (MABC), provides broad-spectrum and potentially durable resistance while limiting the donor genome in the background of an elite cultivar. In this work, MABC coupled to foreground and background selections based on KASP marker assays has been applied to introgress four Pi genes (Piz, Pib, Pita, and Pik) in a renowned japonica Italian rice variety, highly susceptible to blast. Molecular analyses on the backcross (BC) lines highlighted the presence of an additional blast-resistance gene, the Pita-linked Pita2/Ptr gene, therefore increasing the number of blast-resistance introgressed genes to five. The recurrent genome was recovered up to 95.65%. Several lines carrying four (including Pita2) Pi genes with high recovery percentage levels were also obtained. Phenotypic evaluations confirmed the effectiveness of the pyramided lines against multivirulent strains, which also had broad patterns of resistance in comparison to those expected based on the pyramided Pi genes. The developed blast-resistant japonica lines represent useful donors of multiple blast-resistance genes for future rice-breeding programs related to the japonica group.

Identification of QTLs and critical genes related to sugarcane mosaic disease resistance

Mosaic viral diseases affect sugarcane productivity worldwide. Mining disease resistance-associated molecular markers or genes is a key component of disease resistance breeding programs. In the present study, 285 F₁ progeny were produced from a cross between Yuetang 93-159, a moderately resistant variety, and ROC22, a highly susceptible variety. The mosaic disease symptoms of these progenies, with ROC22 as the control, were surveyed by natural infection under 11 different environmental conditions in the field and by artificial infections with a mixed sugarcane mosaic virus (SCMV) and sorghum mosaic virus (SrMV) inoculum. Analysis of consolidated survey data enabled the identification of 29 immune, 55 highly resistant, 70 moderately resistant, 62 susceptible, and 40 highly susceptible progenies. The disease response data and a high-quality SNP genetic map were used in quantitative trait locus (QTL) mapping. The results showed that the correlation coefficients (0.26~0.91) between mosaic disease resistance and test environments were significant (p< 0.001), and that mosaic disease resistance was a highly heritable quantitative trait (H² = 0.85). Seven mosaic resistance QTLs were located to the SNP genetic map, each QTL accounted for 3.57% ~ 17.10% of the phenotypic variation explained (PVE). Furthermore, 110 pathogen response genes and 69 transcription factors were identified in the QTLs interval. The expression levels of nine genes (Soffic.07G0015370-1P, Soffic.09G0015410-2T, Soffic.09G0016460-1T, Soffic.09G0016460-1P, Soffic.09G0017080-3C, Soffic.09G0018730-3P, Soffic.09G0018730-3C, Soffic.09G0019920-3C and Soffic.03G0019710-2C) were significantly different between resistant and susceptible progenies, indicating their key roles in sugarcane resistance to SCMV and SrMV infection. The seven QTLs and nine genes can provide a certain scientific reference to help sugarcane breeders develop varieties resistant to mosaic diseases.

In Silico Mining and Characterization of High-Quality SNP/Indels in Some Agro-Economically Important Species Belonging to the Family Euphorbiaceae

(1) Background: To assess the genetic makeup among the agro-economically important members of Euphorbiaceae, the present study was conducted to identify and characterize high-quality single-nucleotide polymorphism (SNP) markers and their comparative distribution in exonic and intronic regions from the publicly available expressed sequence tags (ESTs). (2) Methods: Quality sequences obtained after pre-processing by an EG assembler were assembled into contigs using the CAP3 program at 95% identity; the mining of SNP was performed by QualitySNP; GENSCAN (standalone) was used for detecting the distribution of SNPs in the exonic and intronic regions. (3) Results: A total of 25,432 potential SNPs (pSNP) and 14,351 high-quality SNPs (qSNP), including 2276 indels, were detected from 260,479 EST sequences. The ratio of quality SNP to potential SNP ranged from 0.22 to 0.75. A higher frequency of transitions and transversions was observed more in the exonic than the intronic region, while indels were present more in the intronic region. C↔T (transition) was the most dominant nucleotide substitution, while in transversion, A↔T was the dominant nucleotide substitution, and in indel, A/- was dominant. (4) Conclusions: Detected SNP markers may be useful for linkage mapping; marker-assisted breeding; studying genetic diversity; mapping important phenotypic traits, such as adaptation or oil production; or disease resistance by targeting and screening mutations in important genes.

Genetic diversity and population structure assessment of Western Canadian barley cooperative trials

Studying the population structure and genetic diversity of historical datasets is a proposed use for association analysis. This is particularly important when the dataset contains traits that are time-consuming or costly to measure. A set of 96 elite barley genotypes, developed from eight breeding programs of the Western Canadian Cooperative Trials were used in the current study. Genetic diversity, allelic variation, and linkage disequilibrium (LD) were investigated using 5063 high-quality SNP markers via the Illumina 9K Barley Infinium iSelect SNP assay. The distribution of SNPs markers across the barley genome ranged from 449 markers on chromosome 1H to 1111 markers on chromosome 5H. The average polymorphism information content (PIC) per locus was 0.275 and ranged from 0.094 to 0.375. Bayesian clustering in STRUCTURE and principal coordinate analysis revealed that the populations are differentiated primarily due to the different breeding program origins and ear-row type into five subpopulations. Analysis of molecular variance based on PhiPT values suggested that high values of genetic diversity were observed within populations and accounted for 90% of the total variance. Subpopulation 5 exhibited the most diversity with the highest values of the diversity indices, which represent the breeding program gene pool of AFC, AAFRD, AU, and BARI. With increasing genetic distance, the LD values, expressed as r², declined to below the critical r^{2 =} 0.18 after 3.91 cM, and the same pattern was observed on each chromosome. Our results identified an important pattern of genetic diversity among the Canadian barley panel that was proposed to be representative of target breeding programs and may have important implications for association mapping in the future. This highlight, that efforts to identify novel variability underlying this diversity may present practical breeding opportunities to develop new barley genotypes.

Characterization of fatty acid compositions in longissimus thoracis muscle and identification of candidate gene and SNPs related to polyunsaturated fatty acid in Hu sheep

Fatty acid (FA) composition contributes greatly to the quality and nutritional value of lamb meat. In the present study, FA was measured in longissimus thoracis (LT) muscles of 1,085 Hu sheep using gas chromatography. Comparative transcriptomic analysis was conducted in LT muscles to identify differentially expressed genes (DEGs) between six individuals with high polyunsaturated fatty acids (H-PUFA, 15.27% ± 0.42%) and six with low PUFA (L-PUFA, 5.22% ± 0.25%). Subsequently, the single nucleotide polymorphisms (SNPs) in a candidate gene PLIN2 were correlated with FA traits. The results showed a total of 29 FA compositions and 8 FA groups were identified, with the highest content of monounsaturated fatty acids (MUFA, 46.54%, mainly C18:1n9c), followed by saturated fatty acids (SFA, 44.32%, mainly C16:0), and PUFA (8.72%, mainly C18:2n6c), and significant correlations were observed among the most of FA traits. Transcriptomic analyses identified 110 upregulated and 302 downregulated DEGs between H-PUFA and L-PUFA groups. The functional enrichment analysis revealed three significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 17 gene ontology (GO) terms, in which regulation of lipolysis in adipocytes, the AMPK signaling pathway, and the PPAR signaling pathway may play important roles in FA metabolism and biosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified 37 module genes associated with PUFA-related traits. In general, PLIN1, LIPE, FABP4, LEP, ACACA, ADIPOQ, SCD, PCK2, FASN, PLIN2, LPL, FABP3, THRSP, and ACADVL may have a great impact on PUFA metabolism and lipid deposition. Four SNPs within PLIN2 were significantly associated with FA. Of those, SNP1 (g.287 G>A) was significantly associated with C18:1n9c and MUFA, and SNP4 (g.7807 T>C) was significantly correlated with PUFA (C18:3n3). In addition, the combined genotype of SNP1 (g.287 G>A), SNP3 (g.7664 T>C), and SNP4 (g.7807 T>C) were significantly correlated with C16:1, C17:0, C18:1C6, PUFA (C18:3n3, C22:6n3), and n-6/n-3 PUFA. These results contribute to the knowledge of the biological mechanisms and genetic markers involved in the composition of FA in Hu sheep.

Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis

Tea plants are widely grown all over the world because they are an important economic crop. The purity and authenticity of tea varieties are frequent problems in the conservation and promotion of germplasm resources in recent years, which has brought considerable inconvenience and uncertainty to the selection of parental lines for breeding and the research and cultivation of superior varieties. However, the development of core SNP markers can quickly and accurately identify the germplasm, which plays an important role in germplasm identification and the genetic relationship analysis of tea plants. In this study, based on 179,970 SNP loci from the whole genome of the tea plant, all of 142 cultivars were clearly divided into three groups: Assam type (CSA), Chinese type (CSS), and transitional type. Most CSA cultivars are from Yunnan Province, which confirms that Yunnan Province is the primary center of CSA origin and domestication. Most CSS cultivars are distributed in east China; therefore, we deduced that east China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. Moreover, 45 core markers were screened using strict criteria to 179,970 SNP loci, and we analyzed 117 well-Known tea cultivars in China with 45 core SNP markers. The results were as follows: (1) In total, 117 tea cultivars were distinguished by eight markers, which were selected to construct the DNA fingerprint, and the remaining markers were used as standby markers for germplasm identification. (2) Ten pairs of parent and offspring relationships were confirmed or identified, and among them, seven pairs were well-established pedigree relationships; the other three pairs were newly identified. In this study, the east of China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. The 45 core SNP markers were developed, which provide a scientific basis at the molecular level to identify the superior tea germplasm, undertake genetic relationship analysis, and benefit subsequent breeding work.

Analysis of Genetic Diversity and Population Structure of Cowpea (Vigna unguiculata (L.) Walp) Genotypes Using Single Nucleotide Polymorphism Markers

Cowpea (Vigna unguiculata (L.) Walp) is an important legume crop with immense potential for nutritional and food security, income generation, and livestock feed in Sub-Saharan Africa. The crop is highly tolerant to heat and drought stresses which makes it an extremely important crop for improving resilience in crop production in the face of climate change. This study was carried out to assess the genetic diversity and population structure of 90 cowpea accessions using single nucleotide polymorphism (SNP) markers. Out of 11,940 SNPs used, 5864 SNPs were polymorphic and maintained for genome diversity analysis. Polymorphic information content (PIC) values ranged from 0.22 to 0.32 with a mean value of 0.27. The model-based Bayesian STRUCTURE analysis classified 90 cowpea accessions into four subpopulations at K = 4, while the distance-based cluster analysis grouped the accessions into three distinct clusters. The analysis of molecular variance (AMOVA) revealed that 59% and 69% of the total molecular variation was attributed to among individual variation for model-based and distance-based populations, respectively, and 18% was attributed to within individual variations. Furthermore, the low heterozygosity among cowpea accessions and the high inbreeding coefficient observed in this study suggests that the accessions reached an acceptable level of homozygosity. This study would serve as a reference for future selection and breeding programs of cowpea with desirable traits and systematic conservation of these plant genetic resources.

Construction of an SNP fingerprinting database and population genetic analysis of 329 cauliflower cultivars

Cauliflower is one of the most important vegetable crops grown worldwide. However, the lack of genetic diversity information and efficient molecular markers hinders efforts to improve cauliflower. This study aims to construct DNA fingerprints for 329 cauliflower cultivars based on SNP markers and the KASP system. After rigorous filtering, a total of 1662 candidate SNPs were obtained from nearly 17.9 million SNP loci. The mean values of PIC, MAF, heterozygosity and gene diversity of these SNPs were 0.389, 0.419, 0.075, and 0.506, respectively. We developed a program for in silico simulations on 153 core germplasm samples to generate ideal SNP marker sets from the candidates. Finally, 41 highly polymorphic KASP markers were selected and applied to identify 329 cauliflower cultivars, mainly collected from the public market. Furthermore, based on the KASP genotyping data, we performed phylogenetic analysis and population structure analysis of the 329 cultivars. As a result, these cultivars could be classified into three major clusters, and the classification patterns were significantly related to their curd solidity and geographical origin. Finally, fingerprints of the 329 cultivars and 2D barcodes with the genetic information of each sample were generated. The fingerprinting database developed in this study provides a practical tool for identifying the authenticity and purity of cauliflower seeds and valuable genetic information about the current cauliflower cultivars.

Detection and Discrimination of Single Nucleotide Polymorphisms by Quantification of CRISPR-Cas Catalytic Efficiency

The specificity of CRISPR-Cas12 assays is attractive for the detection of single nucleotide polymorphisms (SNPs) implicated in, e.g., cancer and SARS-CoV-2 variants. Such assays often employ endpoint measurements of SNP or wild type (WT) activated Cas12 trans-cleavage activity; however, the fundamental kinetic effects of SNP versus WT activation remain unknown. We here show that endpoint-based assays are limited by arbitrary experimental choices (like used reporter concentration and assay duration) and work best for known target concentrations. More importantly, we show that SNP (versus WT) activation results in measurable kinetic shifts in the Cas12 trans-cleavage substrate affinity (K_M) and apparent catalytic efficiency (k_cat^*/K_M). To address endpoint-based assay limitations, we then develop an assay based on the quantification of Michaelis-Menten parameters and apply this assay to a 20 base pair WT target of the SARS-CoV-2 E gene. We find that the k_cat^*/K_M measured for WT is 130-fold greater than the lowest k_cat^*/K_M among all 60 measured SNPs (compared to a 4.8-fold for endpoint fluorescence of the same SNP). K_M also offers a strong ability to distinguish SNPs, varies 27-fold over all the cases, and, importantly, is insensitive to the target concentration. Last, we point out trends among kinetic rates and SNP base and location within the CRISPR-Cas12 targeted region.

Combining QTL mapping and RNA-Seq Unravels candidate genes for Alfalfa (Medicago sativa L.) leaf development

BACKGROUND

Leaf size affects crop canopy morphology and photosynthetic efficiency, which can influence forage yield and quality. It is of great significance to mine the key genes controlling leaf development for breeding new alfalfa varieties. In this study, we mapped leaf length (LL), leaf width (LW), and leaf area (LA) in an F1 mapping population derived from a cultivar named ZhongmuNo.1 with larger leaf area and a landrace named Cangzhou with smaller leaf area.

RESULTS

This study showed that the larger LW was more conducive to increasing LA. A total of 24 significant quantitative trait loci (QTL) associated with leaf size were identified on both the paternal and maternal linkage maps. Among them, nine QTL explained about 11.50-22.45% phenotypic variation. RNA-seq analysis identified 2,443 leaf-specific genes and 3,770 differentially expressed genes. Combining QTL mapping, RNA-seq alalysis, and qRT-PCR, we identified seven candidate genes associated with leaf development in five major QTL regions.

CONCLUSION

Our study will provide a theoretical basis for marker-assisted breeding and lay a foundation for further revealing molecular mechanism of leaf development in alfalfa.

Comparative genomic analysis and phylogeny of NAC25 gene from cultivated and wild Coffea species

Coffee is a high value agricultural commodity grown in about 80 countries. Sustainable coffee cultivation is hampered by multiple biotic and abiotic stress conditions predominantly driven by climate change. The NAC proteins are plants specific transcription factors associated with various physiological functions in plants which include cell division, secondary wall formation, formation of shoot apical meristem, leaf senescence, flowering embryo and seed development. Besides, they are also involved in biotic and abiotic stress regulation. Due to their ubiquitous influence, studies on NAC transcription factors have gained momentum in different crop plant species. In the present study, NAC25 like transcription factor was isolated and characterized from two cultivated coffee species, Coffea arabica and Coffea canephora and five Indian wild coffee species for the first time. The full-length NAC25 gene varied from 2,456 bp in Coffea jenkinsii to 2,493 bp in C. arabica. In all the seven coffee species, sequencing of the NAC25 gene revealed 3 exons and 2 introns. The NAC25 gene is characterized by a highly conserved 377 bp NAM domain (N-terminus) and a highly variable C terminus region. The sequence analysis revealed an average of one SNP per every 40.92 bp in the coding region and 37.7 bp in the intronic region. Further, the non-synonymous SNPs are 8-11 fold higher compared to synonymous SNPs in the non-coding and coding region of the NAC25 gene, respectively. The expression of NAC25 gene was studied in six different tissue types in C. canephora and higher expression levels were observed in leaf and flower tissues. Further, the relative expression of NAC25 in comparison with the GAPDH gene revealed four folds and eight folds increase in expression levels in green fruit and ripen fruit, respectively. The evolutionary relationship revealed the independent evolution of the NAC25 gene in coffee.

Identification of loci controlling timing of stem elongation in red clover using genotyping by sequencing of pooled phenotypic extremes

MAIN CONCLUSION

Through selective genotyping of pooled phenotypic extremes, we identified a number of loci and candidate genes putatively controlling timing of stem elongation in red clover. We have identified candidate genes controlling the timing of stem elongation prior to flowering in red clover (Trifolium pratense L.). This trait is of ecological and agronomic significance, as it affects fitness, competitivity, climate adaptation, forage and seed yield, and forage quality. We genotyped replicate pools of phenotypically extreme individuals (early and late-elongating) within cultivar Lea using genotyping-by-sequencing in pools (pool-GBS). After calling and filtering SNPs and GBS locus haplotype polymorphisms, we estimated allele frequencies and searched for markers with significantly different allele frequencies in the two phenotypic groups using BayeScan, an F_ST-based test utilizing replicate pools, and a test based on error variance of replicate pools. Of the three methods, BayeScan was the least stringent, and the error variance-based test the most stringent. Fifteen significant markers were identified in common by all three tests. The candidate genes flanking the markers include genes with potential roles in the vernalization, autonomous, and photoperiod regulation of floral transition, hormonal regulation of stem elongation, and cell growth. These results provide a first insight into the potential genes and mechanisms controlling transition to stem elongation in a perennial legume, which lays a foundation for further functional studies of the genetic determinants regulating this important trait.

Genetic architecture and QTL selection response for Kernza perennial grain domestication traits

Analysis of multi-year breeding program data revealed that the genetic architecture of an intermediate wheatgrass population was highly polygenic for both domestication and agronomic traits, supporting the use of genomic selection for new crop domestication. Perennial grains have the potential to provide food for humans and decrease the negative impacts of annual agriculture. Intermediate wheatgrass (IWG, Thinopyrum intermedium, Kernza®) is a promising perennial grain candidate that The Land Institute has been breeding since 2003. We evaluated four consecutive breeding cycles of IWG from 2016 to 2020 with each cycle containing approximately 1100 unique genets. Using genotyping-by-sequencing markers, quantitative trait loci (QTL) were mapped for 34 different traits using genome-wide association analysis. Combining data across cycles and years, we found 93 marker-trait associations for 16 different traits, with each association explaining 0.8-5.2% of the observed phenotypic variance. Across the four cycles, only three QTL showed an F_ST differentiation > 0.15 with two corresponding to a decrease in floret shattering. Additionally, one marker associated with brittle rachis was 216 bp from an ortholog of the btr2 gene. Power analysis and quantitative genetic theory were used to estimate the effective number of QTL, which ranged from a minimum of 33 up to 558 QTL for individual traits. This study suggests that key agronomic and domestication traits are under polygenic control and that molecular methods like genomic selection are needed to accelerate domestication and improvement of this new crop.

Fruit and Vegetable Consumption Interacts With HNF1A Variants on the C-Reactive Protein

Epidemiological studies have demonstrated the inverse association between the intake of fruits and vegetables and inflammation. However, the mechanisms by which inflammation-related genes interact with fruit and vegetable intake and the role of these combinations in inflammation remain unclear. Therefore, we assessed the effect of interactions between fruit and vegetable intake and the hepatic nuclear factor 1 alpha (HNF1A) genetic variants on the C-reactive protein (CRP) levels. Baseline data from the Ansan and Ansung Cohort Study of the Korean Genome and Epidemiology Study (KoGES) were used. A total of 7,634 participants (3,700 men and 3,934 women) were included in the analyses. Fruit and vegetable intake was assessed using semi-quantitative food frequency questionnaire data. Genotyping information for HNF1A was extracted from the Affymetrix Genome-Wide Human SNP array 5.0. Inflammation was determined after overnight fasting by measuring CRP levels using automated analyzers. Multivariable logistic regression was used to estimate the adjusted odds ratio (AOR) with a 95% confidence interval (CI). In the fully adjusted model, men and women with the GG genotype of HNF1A rs2393791 and high fruit intake had lower odds of elevated CRP levels compared to those with the AA genotype and low fruit intake (AOR 0.50, 95% CI 0.38–0.67; AOR 0.73, 95% CI 0.55–0.97, respectively). Men and women with the rs2393791 GG genotype and high vegetable intake had lower odds of having elevated CRP levels compared to those with the AA genotype and low fruit intake (AOR 0.57, 95% CI 0.43–0.75; AOR 0.65, 95% CI 0.49–0.86, respectively). Men and women with the GG genotype and high total fruit and vegetable intake had lower odds of having elevated CRP levels. These findings indicate that fruit and vegetable intake interacts with HNF1A genetic polymorphisms, consequently influencing the inflammation levels.

Efficient marker-assisted breeding for clubroot resistance in elite Pol-CMS rapeseed varieties by updating the PbBa8.1 locus

Clubroot disease poses a severe threat to rapeseed (Brassica napus) production worldwide and has recently been spreading across China at an unprecedented pace. Breeding and cultivation of resistant varieties constitute a promising and environment-friendly approach to mitigating this threat. In this study, the clubroot resistance locus PbBa8.1 was successfully transferred into SC4, a shared paternal line of three elite varieties in five generations by marker-assisted backcross breeding. Kompetitive allele specific PCR (KASP) markers of clubroot resistance gene PbBa8.1 and its linked high erucic acid gene (FAE1) were designed and applied for foreground selection, and 1,000 single-nucleotide polymorphisms (SNPs) were selected and used for the background selection. This breeding strategy produced recombinants with the highest recovery ratio of the recurrent parent genome (> 95%) at BC2F2 while breaking the linkage with FAE1 during the selection. An updated version of the paternal line (SC4R) was generated at BC2F3, showing significantly improved clubroot resistance at the seedling stage via artificial inoculation, and was comparable to that of the donor parent. Field trials of the three elite varieties and their updated versions in five environments indicated similar agronomic appearance and final yield. The introduced breeding strategy precisely pyramids the PbBa8.1 and FAE1 loci with the assistance of technical markers in a shorter period and could be applied to other desirable traits for directional improvement in the future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01305-9.

Mutations in PmUFGT3 contribute to color variation of fruit skin in Japanese apricot (Prunus mume Sieb. et Zucc.)

BACKGROUND

Japanese apricot (Prunus mume Sieb. et Zucc.) is popular for both ornamental and processing value, fruit color affects the processing quality, and red pigmentation is the most obvious phenotype associated with fruit color variation in Japanese apricot, mutations in structural genes in the anthocyanin pathway can disrupt the red pigmentation, while the formation mechanism of the red color trait in Japanese apricot is still unclear.  RESULTS: One SNP marker (PmuSNP_27) located within PmUFGT3 gene coding region was found highly polymorphic among 44 different fruit skin color cultivars and relative to anthocyanin biosynthesis in Japanese apricot. Meantime, critical mutations were identified in two alleles of PmUFGT3 in the green-skinned type is inactivated by seven nonsense mutations in the coding region, which leads to seven amino acid substitution, resulting in an inactive UFGT enzyme. Overexpression of the PmUFGT3 allele from red-skinned Japanese apricot in green-skinned fruit lines resulted in greater anthocyanin accumulation in fruit skin. Expression of same allele in an Arabidopsis T-DNA mutant deficient in anthocyanidin activity the accumulation of anthocyanins. In addition, using site-directed mutagenesis, we created a single-base substitution mutation (G to T) of PmUFGT3 isolated from green-skinned cultivar, which caused an E to D amino acid substitution and restored the function of the inactive allele of PmUFGT3 from a green-skinned individual.

CONCLUSION

This study confirms the function of PmUFGT3, and provides insight into the mechanism underlying fruit color determination in Japanese apricot, and possible approaches towards genetic engineering of fruit color.

Prediction and expression analysis of deleterious nonsynonymous SNPs of Arabidopsis ACD11 gene by combining computational algorithms and molecular docking approach

Accelerated cell death 11 (ACD11) is an autoimmune gene that suppresses pathogen infection in plants by preventing plant cells from becoming infected by any pathogen. This gene is widely known for growth inhibition, premature leaf chlorosis, and defense-related programmed cell death (PCD) in seedlings before flowering in Arabidopsis plant. Specific amino acid changes in the ACD11 protein’s highly conserved domains are linked to autoimmune symptoms including constitutive defensive responses and necrosis without pathogen awareness. The molecular aspect of the aberrant activity of the ACD11 protein is difficult to ascertain. The purpose of our study was to find the most deleterious mutation position in the ACD11 protein and correlate them with their abnormal expression pattern. Using several computational methods, we discovered PCD vulnerable single nucleotide polymorphisms (SNPs) in ACD11. We analysed the RNA-Seq data, identified the detrimental nonsynonymous SNPs (nsSNP), built genetically mutated protein structures and used molecular docking to assess the impact of mutation. Our results demonstrated that the A15T and A39D mutations in the GLTP domain were likely to be extremely detrimental mutations that inhibit the expression of the ACD11 protein domain by destabilizing its composition, as well as disrupt its catalytic effectiveness. When compared to the A15T mutant, the A39D mutant was more likely to destabilize the protein structure. In conclusion, these mutants can aid in the better understanding of the vast pool of PCD susceptibilities connected to ACD11 gene GLTP domain activation.

Non synonymous single nucleotide polymorphism (nsSNP) is a process in which amino acid sequence of a protein is altered as a result of single nucleotide alteration in the coding region (mRNA) of any living organism. Therefore, the entire protein structure, interactions and stability are altered, which may have a negative impact on living organisms. Hence, to completely comprehend this biological process, we must first solve the unresolved mutational protein structure and mutated protein interactions. The major goal of our research is to identify the most harmful mutation in our target protein structure and how it interacts within cells. However, it was discovered that only a few alterations in residues had the largest negative impact on the protein’s internal structure and also on the protein-ligand interactions. We show that based on the amino acid sequence of a protein computationally, it is feasible to discover mutational positions in the sequence, generate mutation protein structure and interactions with related ligands. Our findings show that the essential mechanisms underlying protein mutations generated by this process are identical. The capacity to correctly detect mutations from sequence allows the annotation and study of protein-ligand interactions throughout a whole organism, which might aid function prediction and gene expression.

DNA fingerprinting reveals varietal composition of Vietnamese cassava germplasm (Manihot esculenta Crantz) from farmers' field and genebank collections

A molecular analysis using informative SNP markers in 1570 clones of cassava from Vietnam reveals varietal composition from farmers' field and genebank collections Cassava is the most important smallholder cash crops in Southeast Asia and is especially used in industrial products. Yet, systematic genetic studies on molecular markers from Vietnamese germplasm have not been considered for breeding and conservation programs. We conducted a molecular analysis of 1570 clones of cassava germplasm from farms across six agro-ecological zones using informative SNP markers. We unraveled the genetic diversity and population structure and provided insights into the value of breeding and conservation programs. Duplicated genotypes comprised 98% of the total sample of the Central Highlands region. Ninety-six SNPs were amplified Central Highlands and South East provinces had the highest allelic richness, covering up to 83% of alleles. The average observed heterozygosity (Ho = 0.43) was slightly higher than expected (He = 0.40) across SNP markers, suggesting an excess of heterozygotes plants. Diversity indexes indicated that cassava populations from North West and Eastern Vietnam are genetically diverse (mean He = 0.40). Genetic parentage tests identified 85 unique genetic groups within the varieties KM94, KM419, BRA1305, KM101, KM140, PER262, KM60, KM57 and two unidentified varieties, which accounted for 82% of the frequency distribution. KM94 is the most dominant variety in Vietnamese farms surveyed (38%), reflecting its superior quality and productivity. Discriminant analysis of principal components (DAPC) revealed four main subgroups, which were partially corroborated by neighbor joining (NJ) analyses. After removing duplicates, 31 unique genotypes were distributed across five of the agro-ecological zones. These were well distributed in the subgroups revealed via DAPC and NJ analyses. The genetic groups identified herein could be used to select unique accessions that should ideally conform with ex situ germplasm collections and identify areas where on-farm conservation programs should be targeted. Newly identified genotypes may also contribute as genetic breeding resources that could be used to adapt cassava to future changes and farmers' needs.

Linkage disequilibrium and population structure in a core collection of Brassica napus (L.)

Estimation of genetic diversity in rapeseed is important for sustainable breeding program to provide an option for the development of new breeding lines. The objective of this study was to elucidate the patterns of genetic diversity within and among different structural groups, and measure the extent of linkage disequilibrium (LD) of 383 globally distributed rapeseed germplasm using 8,502 single nucleotide polymorphism (SNP) markers. We divided the germplasm collection into five subpopulations (P1 to P5) according to geographic and growth habit-related patterns. All subpopulations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise Fst comparison revealed a great degree of divergence (Fst > 0.24) between most of the combinations. The rutabaga type showed highest divergence with spring and winter types. Higher divergence was also found between winter and spring types. Admixture model based structure analysis, principal component and neighbor-joining tree analysis placed all subpopulations into three distinct clusters. Admixed genotype constituted 29.24% of total genotypes, while remaining 70.76% belongs to identified clusters. Overall, mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 45 kb distance for whole genome. The LD decay was slower in C genome (< 93 kb); relative to the A genome (< 21 kb) which was confirmed by availability of larger haplotype blocks in C genome than A genome. The findings regarding LD pattern and population structure will help to utilize the collection as an important resource for association mapping efforts to identify genes useful in crop improvement as well as for selection of parents for hybrid breeding.

Population Genetics and Development of a Core Collection from Elite Germplasms of Xanthoceras sorbifolium Based on Genome-Wide SNPs

Xanthoceras sorbifolium is one of the most important species of woody oil. In this study, whole genome re-sequencing of 119 X. sorbifolium germplasms was conducted and, after filtering, 105,685,557 high-quality SNPs were identified, which were used to perform population genetics and core collection development analyses. The results from the phylogenetic, population structure, and principal component analyses showed a high level of agreement, with 119 germplasms being classified into three main groups. The germplasms were not completely classified based on their geographical origins and flower colors; furthermore, the genetic backgrounds of these germplasms were complex and diverse. The average polymorphsim information content (PIC) values for the three inferred groups clustered by structure analysis and the six classified color groups were 0.2445 and 0.2628, respectively, indicating a low to medium informative degree of genetic diversity. Moreover, a core collection containing 29.4% (35) out of the 119 X. sorbifolium germplasms was established. Our results revealed the genetic diversity and structure of X. sorbifolium germplasms, and the development of a core collection will be useful for the efficient improvement of breeding programs and genome-wide association studies.

Genome-Wide Association Study of Healthful Flavonoids among Diverse Mandarin Accessions

Mandarins have many unique flavonoids with documented health benefits and that help to prevent chronic human diseases. Flavonoids are difficult to measure and cannot be phenotyped without the use of specialized equipment; consequently, citrus breeders have not used flavonoid contents as selection criteria to develop cultivars with increased benefits for human health or increased tolerance to diseases. In this study, peel, pulp, and seed samples collected from many mandarin accessions and their hybrids were analyzed for the presence of selected flavonoids with documented human health benefits. A genome-wide association study (GWAS) was used to identify SNPs associated with biosynthesis of flavonoids in these mandarin accessions, and there were 420 significant SNPs were found to be associated with 28 compounds in peel, pulp, or seed samples. Four candidate genes involved in flavonoid biosynthesis were identified by enrichment analysis. SNPs that were found to be associated with compounds in pulp samples have the potential to be used as markers to select mandarins with improved phytonutrient content to benefit human health. Mandarin cultivars bred with increased flavonoid content may provide value to growers and consumers.

Utilization of a Sugarcane100K Single Nucleotide Polymorphisms Microarray-Derived High-Density Genetic Map in Quantitative Trait Loci Mapping and Function Role Prediction of Genes Related to Chlorophyll Content in Sugarcane

Chlorophyll is the most important pigment for plant photosynthesis that plays an important role in crop growth and production. In this study, the chlorophyll content trait was explored to improve sugarcane yield. Two hundred and eighty-five F₁ progenies from the cross YT93-159 × ROC22 with significantly different chlorophyll contents were included as test materials. The chlorophyll content of the +1 leaves during elongation phase was measured using a SPAD-502 meter through a three-crop cycle (plant cane, first ratoon, and second ratoon). Linkage analysis was conducted on a high-density genetic map constructed based on the sugarcane 100K SNP chip. In addition, Fv/Fm, plant height, stalk diameter, brix data were collected on plant cane during the elongation and maturation phases. The results showed that the +1 leaf SPAD values, which can be used as an important reference to evaluate the growth potential of sugarcane, were significantly and positively correlated with the Fv/Fm during elongation phase, as well as with plant height, stalk diameter, and brix during maturity phase (P < 0.01). The broad sense heritability (H ²) of the chlorophyll content trait was 0.66 for plant cane crop, 0.67 for first ratoon crop, and 0.73 for second ratoon crop, respectively, indicating that this trait was mainly controlled by genetic factors. Thirty-one quantitative trait loci (QTL) were detected by QTL mapping. Among them, a major QTL, qCC-R1, could account for 12.95% of phenotypic variation explained (PVE), and the other 30 minor QTLs explained 2.37-7.99% PVE. Twenty candidate genes related to chlorophyll content were identified in the QTLs plus a 200-Kb extension region within either sides, of which four were homologous genes involved in the chlorophyll synthesis process and the remaining 16 played a certain role in chlorophyll catabolic pathway, chloroplast organization, or photosynthesis. These results provide a theoretical reference for analyzing the genetic mechanism of chlorophyll synthesis and subsequent improvement of photosynthetic characteristics in sugarcane.

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.

Association analysis for agronomic traits in wheat under terminal heat stress

•

Terminal heat stress leads to irreversible damage in wheat.

•

Marker assisted selection and gene pyramiding for portrayal of heat tolerance.

•

Allelic frequency and polymorphic information showed significant variability.

•

Markers xcfa2147 and xwmc671 could be potentail for heat stress tolerance.

Terminal heat stress causes irreversible damage to wheat crop productivity. It reduces the vegetative growth and flowering period that consequently declines the efficiency to capture available stem reserves (carbohydrates) in grains. Markers associated with thermotolerant traits ease in marker assisted selection (MAS) for crop improvement. It identifies the genomic regions associated with thermotolerant traits in wheat, but the scarcity of markers is the major hindrance in crop improvement. Therefore, 158 wheat genotypes were subjected to genotyping with 165 simple sequence repeat markers dispersed on three genomes (A, B and D). Allelic frequency and polymorphic information content values were highest on genome A (5.34 (14% greater than the lowest value at genome D) and 0.715 (3% greater than the lowest value at genome D)), chromosome 4 (5.40 (16% greater than the lowest value at chromosome 2) and 0.725 (5% greater than the lowest value at chromosome 6)) and marker xgwm44 (13.0 (84% greater than the lowest value at marker xbarc148) and 0.916 (46% greater than the lowest value at marker xbarc148)). Bayesian based population structure discriminated the wheat genotypes into seven groups based on genetic similarity indicating their ancestral origin and geographical ecotype. Linkage disequilibrium pattern had highest significant (P < 0.001) linked loci pairs 732 on genome A at r² > 0.1 whereas, 58 on genome B at r² > 0.5. Linkage disequilibrium decay (P < 0.01 and r² > 0.1) had larger LD block (5–10 cM) on genome A. Highly significant MTAs (P < 0.000061) under heat stress conditions were identified for flag leaf area (xwmc336), spikelet per spike (xwmc553), grains per spike (cxfa2147, xwmc418 and xwmc121), biomass (xbarc7) and grain yield (xcfa2147 and xwmc671). The identified markers in this study could facilitate in MAS and gene pyramiding against heat stress in wheat.

Genome-wide association analysis of chickpea germplasms differing for salinity tolerance based on DArTseq markers

Soil salinity is significant abiotic stress that severely limits global crop production. Chickpea (Cicer arietinum L.) is an important grain legume that plays a substantial role in nutritional food security, especially in the developing world. This study used a chickpea population collected from the International Center for Agricultural Research in the Dry Area (ICARDA) genebank using the focused identification of germplasm strategy. The germplasm included 186 genotypes with broad Asian and African origins and genotyped with 1856 DArTseq markers. We conducted phenotyping for salinity in the field (Arish, Sinai, Egypt) and greenhouse hydroponic experiments at 100 mM NaCl concentration. Based on the performance in both hydroponic and field experiments, we identified seven genotypes from Azerbaijan and Pakistan (IGs: 70782, 70430, 70764, 117703, 6057, 8447, and 70249) as potential sources for high salinity tolerance. Multi-trait genome-wide association analysis (mtGWAS) detected one locus on chromosome Ca4 at 10618070 bp associated with salinity tolerance under hydroponic and field conditions. In addition, we located another locus specific to the hydroponic system on chromosome Ca2 at 30537619 bp. Gene annotation analysis revealed the location of rs5825813 within the Embryogenesis-associated protein (EMB8-like), while the location of rs5825939 is within the Ribosomal Protein Large P0 (RPLP0). Utilizing such markers in practical breeding programs can effectively improve the adaptability of current chickpea cultivars in saline soil. Moreover, researchers can use our markers to facilitate the incorporation of new genes into commercial cultivars.

Genetic Diversity and Population Structure of Myanmar Rice (Oryza sativa L.) Varieties Using DArTseq-Based SNP and SilicoDArT Markers

Myanmar is well known as a primary center of plant genetic resources for rice. A considerable number of genetic diversity studies have been conducted in Myanmar using various DNA markers. However, this is the first report using DArTseq technology for exploring the genetic diversity of Myanmar rice. In our study, two ultra-high-throughput diversity array technology markers were employed to investigate the genetic diversity and population structure of local rice varieties in the Ayeyarwady delta, the major region of rice cultivation. The study was performed using 117 rice genotypes with 7643 SNP and 4064 silicoDArT markers derived from the DArT platform. Genetic variance among the genotypes ranged from 0 to 0.753 in SNPs, and from 0.001 to 0.954 in silicoDArT. Two distinct population groups were identified from SNP data analysis. Cluster analysis with both markers clearly separated traditional Pawsan varieties and modern high-yielding varieties. A significant divergence was found between populations according to the Fst values (0.737) obtained from the analysis of molecular variance, which revealed 74% genetic variation at the population level. These findings support rice researchers in identifying useful DNA polymorphisms in genes and pinpointing specific genes conferring desirable phenotypic traits for further genome-wide association studies and parental selection for recombination breeding to enhance rice varietal development and release.

Insights on Genetic Diversity, Population Structure, and Linkage Disequilibrium in Globally Diverse Coconut Accessions Using Genotyping-by-Sequencing

Genotyping-by-sequencing (GBS) has emerged as a cost-effective approach for genome-wide discovery of single-nucleotide polymorphism (SNP) markers and high-throughput genotyping. In this study, 96 coconut palms, representing 16 accessions from globally diverse origins, were genotyped using the GBS strategy. A total of 10,835 high-quality SNPs, which were identified after stringent filtering, were utilized to assess genetic diversity, population structure, and linkage disequilibrium (LD) analyses. The polymorphism information content (PIC) values of SNPs ranged from 0.1 to 0.4, with a large proportion of SNPs (8633 nos.; 79.7%) having a higher PIC in the range of 0.3-0.4. The genetic diversity analysis revealed the existence of a high level of variation in coconut accessions, with an average expected heterozygosity (H_e) value of 0.43. Unweighted neighbor-joining phylogenetic tree and Bayesian-based model population structure grouped coconut genotypes into four main clusters. The accessions are generally clustered based on their height (tall or dwarf), with a few accession clusterings based on geographical origins. Investigation of LD pattern in coconut indicated a relatively rapid LD decay with a short range (9 kb). The results obtained in this study will contribute to enhancing the capacity of coconut researchers to utilize genetic diversity for further genetic improvement. In addition, it would open up possibilities for performing genomic studies such as genome-wide association studies and genomic selection to accelerate the efficiency and speed of coconut genetic improvement.

Genetic diversity analysis of Korean peanut germplasm using 48 K SNPs 'Axiom_Arachis' Array and its application for cultivar differentiation

Cultivated peanut (Arachis hypogaea) is one of the important legume oilseed crops. Cultivated peanut has a narrow genetic base. Therefore, it is necessary to widen its genetic base and diversity for additional use. The objective of the present study was to assess the genetic diversity and population structure of 96 peanut genotypes with 9478 high-resolution SNPs identified from a 48 K 'Axiom_Arachis' SNP array. Korean set genotypes were also compared with a mini-core of US genotypes. These sets of genotypes were used for genetic diversity analysis. Model-based structure analysis at K = 2 indicated the presence of two subpopulations in both sets of genotypes. Phylogenetic and PCA analysis clustered these genotypes into two major groups. However, clear genotype distribution was not observed for categories of subspecies, botanical variety, or origin. The analysis also revealed that current Korean genetic resources lacked variability compared to US mini-core genotypes. These results suggest that Korean genetic resources need to be expanded by creating new allele combinations and widening the genetic pool to offer new genetic variations for Korean peanut improvement programs. High-quality SNP data generated in this study could be used for identifying varietal contaminant, QTL, and genes associated with desirable traits by performing mapping, genome-wide association studies.

Rapid excavating a FLOWERING LOCUS T-regulator NF-YA using genotyping-by-sequencing

Soybean (Glycine max (L.) Merrill) is one of the most important crop plants in the world as an important source of protein for both human consumption and livestock fodder. As flowering time contributes to yield, finding new QTLs and further identifying candidate genes associated with various flowering time are fundamental to enhancing soybean yield. In this study, a set of 120 recombinant inbred lines (RILs) which was developed from a cross of two soybean cultivars, Suinong4 (SN4) and ZK168, were genotyped by genotyping-by-sequencing (GBS) approach and phenotyped to expand the cognitive of flowering time by quantitative trait loci (QTL) analysis. Eventually, three stable QTLs related to flowering time which were detected separately located on chromosome 14, 18, and 19 under long-day (LD) conditions. We predicted candidate genes for each QTL and carried out association analyses between the putative causal alleles and flowering time. Moreover, a transient transfection assay was performed and showed that NUCLEAR FACTOR YA 1b (GmNF-YA1b) as a strong candidate for the QTL on chromosome 19 might affect flowering time by suppressing the expression of FLOWERING LOCUS T (GmFT) genes in soybean. QTLs detected in this study would provide fundamental resources for finding candidate genes and clarify the mechanisms of flowering which would be helpful for breeding novel high-yielding soybean cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01237-w.

MP3RNA-seq: Massively parallel 3' end RNA sequencing for high-throughput gene expression profiling and genotyping

Transcriptome deep sequencing (RNA-seq) has become a routine method for global gene expression profiling. However, its application to large-scale experiments remains limited by cost and labor constraints. Here we describe a massively parallel 3' end RNA-seq (MP3RNA-seq) method that introduces unique sample barcodes during reverse transcription to permit sample pooling immediately following this initial step. MP3RNA-seq allows for handling of hundreds of samples in a single experiment, at a cost of about $6 per sample for library construction and sequencing. MP3RNA-seq is effective for not only high-throughput gene expression profiling, but also genotyping. To demonstrate its utility, we applied MP3RNA-seq to 477 double haploid lines of maize. We identified 19,429 genes expressed in at least 50% of the lines and 35,836 high-quality single nucleotide polymorphisms for genotyping analysis. Armed with these data, we performed expression and agronomic trait quantitative trait locus (QTL) mapping and identified 25,797 expression QTLs for 15,335 genes and 21 QTLs for plant height, ear height, and relative ear height. We conclude that MP3RNA-seq is highly reproducible, accurate, and sensitive for high-throughput gene expression profiling and genotyping, and should be generally applicable to most eukaryotic species.

The Distribution and Origins of Pyrus hopeiensis-"Wild Plant With Tiny Population" Using Whole Genome Resequencing

Pyrus hopeiensis is a valuable but endangered wild resource in the genus Pyrus. It has been listed as one of the 120 wild species with tiny population in China. The specie has been little studied. A preliminary study of propagation modes in P. hopeiensis was performed through seed propagation, hybridization, self-crossing trials, bud grafting, branch grafting, and investigations of natural growth. The results showed that the population size of P. hopeiensis was very small, the distribution range was limited, and the habitat was extremely degraded. In the wild population, natural hybridization and root tiller production were the major modes of propagation. Whole genome re-sequencing of the 23 wild and cultivated accessions from Pyrus species collected was performed using an Illumina HiSeq sequencing platform. The sequencing depth range was 26.56x-44.85x and the average sequencing depth was 32x. Phylogenetic tree and principal component analyses (PCA) based on SNPs showed that the wild Pyrus species, such as PWH06, PWH07, PWH09, PWH10, PWH13, and PWH17, were closely related to both P. hopeiensis HB-1 and P. hopeiensis HB-2. Using these results in combination with morphological characteristics, it speculated that P. hopeiensis populations may form a natural hybrid group with frequent gene exchanges between and within groups. A selective elimination analysis on the P. hopeiensis population were performed using Fst and π radio and a total of 381 overlapping genes including SAUR72, IAA20, HSFA2, and RKP genes were obtained. These genes were analyzed by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) function enrichment. And four KEGG pathways, including lysine degradation, sphingolipid metabolism, other glycan degradation, and betaine biosynthesis were significantly enriched in the P. hopeiensis population. Our study provided information on genetic variation, evolutionary relationships, and gene enrichment in P. hopeiensis population. These data will help reveal the evolutionary history and origin of P. hopeiensis and provide guidelines for subsequent research on the locations of functional genes.