Retention of agronomically important variation in germplasm core collections: implications for allele mining

Uncovering natural allelic and structural variants of OsCENH3 gene by targeted resequencing and in silico mining in genus Oryza

Plant breeding efforts to boost rice productivity have focused on developing a haploid development pipeline. CENH3 gene has emerged as a leading player that can be manipulated to engineer haploid induction system. Currently, allele mining for the OsCENH3 gene was done by PCR-based resequencing of 33 wild species accessions of genus Oryza and in silico mining of alleles from pre-existing data. We have identified and characterized CENH3 variants in genus Oryza. Our results indicated that the majority CENH3 alleles present in the Oryza gene pool carry synonymous substitutions. A few non-synonymous substitutions occur in the N-terminal Tail domain (NTT). SNP A/G at position 69 was found in accessions of AA genome and non-AA genome species. Phylogenetic analysis revealed that non-synonymous substitutions carrying alleles follow pre-determined evolutionary patterns. O. longistaminata accessions carry SNPs in four codons along with indels in introns 3 and 6. Fifteen haplotypes were mined from our panel; representative mutant alleles exhibited structural variations upon modeling. Structural analysis indicated that more than one structural variant may be exhibited by different accessions of single species (Oryza barthii). NTT allelic mutants, though not directly implicated in HI, may show variable interactions. HI and interactive behavior could be ascertained in future investigations.

A pan-genome data structure induced by pooled sequencing facilitates variant mining in heterogeneous germplasm

Valuable genetic variation lies unused in gene banks due to the difficulty of exploiting heterogeneous germplasm accessions. Advances in molecular breeding, including transgenics and genome editing, present the opportunity to exploit hidden sequence variation directly. Here we describe the pan-genome data structure induced by whole-genome sequencing of pooled individuals from wild populations of Patellifolia spp., a source of disease resistance genes for the related crop species sugar beet (Beta vulgaris). We represent the pan-genome as a map of reads from pooled sequencing of a heterogeneous population sample to a reference genome, plus a BLAST data base of the mapped reads. We show that this basic data structure can be queried by reference genome position or homology to identify sequence variants present in the wild relative, at genes of agronomic interest in the crop, a process known as allele or variant mining. Further we demonstrate the possibility of cataloging variants in all Patellifolia genomic regions that have corresponding single copy orthologous regions in sugar beet. The data structure, termed a "pooled read archive," can be produced, altered, and queried using standard tools to facilitate discovery of agronomically-important sequence variation.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01308-6.

Genomics Enabled Breeding Strategies for Major Biotic Stresses in Pea (Pisum sativum L.)

Pea (Pisum sativum L.) is one of the most important and productive cool season pulse crops grown throughout the world. Biotic stresses are the crucial constraints in harnessing the potential productivity of pea and warrant dedicated research and developmental efforts to utilize omics resources and advanced breeding techniques to assist rapid and timely development of high-yielding multiple stress-tolerant-resistant varieties. Recently, the pea researcher's community has made notable achievements in conventional and molecular breeding to accelerate its genetic gain. Several quantitative trait loci (QTLs) or markers associated with genes controlling resistance for fusarium wilt, fusarium root rot, powdery mildew, ascochyta blight, rust, common root rot, broomrape, pea enation, and pea seed borne mosaic virus are available for the marker-assisted breeding. The advanced genomic tools such as the availability of comprehensive genetic maps and linked reliable DNA markers hold great promise toward the introgression of resistance genes from different sources to speed up the genetic gain in pea. This review provides a brief account of the achievements made in the recent past regarding genetic and genomic resources' development, inheritance of genes controlling various biotic stress responses and genes controlling pathogenesis in disease causing organisms, genes/QTLs mapping, and transcriptomic and proteomic advances. Moreover, the emerging new breeding approaches such as transgenics, genome editing, genomic selection, epigenetic breeding, and speed breeding hold great promise to transform pea breeding. Overall, the judicious amalgamation of conventional and modern omics-enabled breeding strategies will augment the genetic gain and could hasten the development of biotic stress-resistant cultivars to sustain pea production under changing climate. The present review encompasses at one platform the research accomplishment made so far in pea improvement with respect to major biotic stresses and the way forward to enhance pea productivity through advanced genomic tools and technologies.

Whole genome re-sequencing reveals the genetic diversity and evolutionary patterns of Eucommia ulmoides

Eucommia ulmoides (E. ulmoides) is a deciduous perennial tree belonging to the order Garryales, and is known as "living fossil" plant, along with ginkgo (Ginkgo biloba), metaspaca (Metasequoia glyptostroboides) and dove tree (Davidia involucrata Baill). However, the genetic diversity and population structure of E. ulmoides are still  ambiguous nowdays. In this study, we re-sequenced the genomes of 12 E. ulmoides accessions from different major climatic geography regions in China to elucidate the genetic diversity, population structure and evolutionary pattern. By integration of phylogenetic analysis, principal component analysis and population structure analysis based on a number of high-quality SNPs, a total of 12 E. ulmoides accessions were clustered into four different groups. This result is consistent with their geographical location except for group samples from Shanghai and Hunan province. E. ulmoides accessions from Hunan province exhibited a closer genetic relationship with E. ulmoides accessions from Shanghai in China compared with other regions, which is also supported by the result of population structure analyses. Genetic diversity analysis further revealed that E. ulmoides samples in Shanghai and Hunan province were with higher genetic diversity than those in other regions in this study. In addition, we treated the E. ulmoides materials from Shanghai and Hunan province as group A, and the other materials from other places as group B, and then analyzed the evolutionary pattern of E. ulmoides. The result showed the significant differentiation (Fst = 0.1545) between group A and group B. Some candidate highly divergent genome regions were identified in group A by selective sweep analyses, and the function analysis of candidate genes in these regions showed that biological regulation processes could be correlated with the Eu-rubber biosynthesis. Notably, nine genes were identified from selective sweep regions. They were involved in the Eu-rubber biosynthesis and expressed in rubber containing tissues. The genetic diversity research and evolution model of E. ulmoides were preliminarily explored in this study, which laid the foundation for the protection of germplasm resources and the development and utilization of multipurpose germplasm resources in the future.

Molecular and genetic dissection of the USDA rice mini-core collection using high-density SNP markers

Molecular tools and knowledge of crop germplasm are vital for their effective utilization. In this study, we developed 40,866 high-quality and well distributed SNPs for a rice mini-core collection (RMC) developed by the United States Department of Agriculture (USDA). The high-quality SNPs clustered the USDA-RMC into five subpopulations (Ind, indica; Aus, aus; Afr, African rice; TeJ, temperate japonica; TrJ, tropical japonica) and one admixture (Adm). This classification was further confirmed by phylogenetic and principal component analyses. The rice ARO (aromatic) subpopulation of previous studies was re-assigned with Adm and the WD (wild-type) subpopulation was re-defined to the Afr subpopulation because most of its accessions are African cultivated rice. The Aus and Ind subpopulations had a substantially wider genetic variation than the TrJ and TeJ subpopulations. The genetic diversities were much larger between the Ind or Aus subpopulation and the TrJ or TeJ subpopulation than between the Afr subpopulation and the Ind, Aus, TrJ or TeJ subpopulation. Comparative agronomic trait analysis between the subpopulations also supported the genetic structure and variation of the RMC, and suggested the existence of extensive variation in the genes controlling agronomic traits among them. Furthermore, analysis of ancestral membership of the RMC accessions revealed that reproductive barrier or wide incompatibility existed between the Indica and Japonica groups, while gene flow occurred between them. These results provide high-quality SNPs and knowledge of genetic structure and diversity of the USDA-RMC necessary for enhanced rice research and breeding.

Genetic diversity and evolutionary patterns of Taraxacum kok-saghyz Rodin

Taraxacum kok-saghyz Rodin (TKS) is an important potential alternative source of natural inulin and rubber production, which has great significance for the production of industrial products. In this study, we sequenced 58 wild TKS individuals collected from four different geography regions worldwide to elucidate the population structure, genetic diversity, and the patterns of evolution. Also, the first flowering time, crown diameter, morphological characteristics of leaf, and scape of all TKS individuals were measured and evaluated statistically. Phylogenetic analysis based on SNPs and cluster analysis based on agronomic traits showed that all 58 TKS individuals could be roughly divided into three distinct groups: (a) Zhaosu County in Xinjiang (population AB, including a few individuals from population C and D); (b) Tekes County in Xinjiang (population C); and (c) Tuzkol lake in Kazakhstan (population D). Population D exhibited a closer genetic relationship with population C compared with population AB. Genetic diversity analysis further revealed that population expansion from C and D to AB occurred, as well as gene flow between them. Additionally, some natural selection regions were identified in AB population. Function annotation of candidate genes identified in these regions revealed that they mainly participated in biological regulation processes, such as transporter activity, structural molecule activity, and molecular function regulator. We speculated that the genes identified in selective sweep regions may contribute to TKS adaptation to the Yili River Valley of Xinjiang. In general, this study provides new insights in clarifying population structure and genetic diversity analysis of TKS using SNP molecular markers and agronomic traits.

Molecular diversity assessment of a world collection of safflower genotypes by SRAP and SCoT molecular markers

Safflower (Carthamus tinctorius L.) is considered as an oil crop that is rich in medicinal and industrial properties. In this study, the genetic diversity of safflower was assessed using 12 polymorphic sequence-related amplified polymorphism (SRAPs) and 11 polymorphic start codon targeted (SCoT) markers in 100 genotypes of safflower gathered from different geographical regions of the world. The 23 primers generated a total of 227 polymorphism fragments with a mean of 68.2% within the range of 3 (SCoT 31 and SCoT 35) to 13 (SCoT 35) bands per primer. Polymorphism per primer ranged between 100% (in Me4-Em1) and 18.1% (in SCoT19), with an average of 36.76%. The polymorphism information contents of the SRAP and SCoT markers were 0.35 and 0.30, respectively, indicating that SRAP markers were more effective than SCoT markers for assessing the degree of genetic diversity of the safflower. The results of the analysis of molecular variance showed a significant difference across cultivated safflower genotypes possessing a high intra-population variation. The examined accessions were categorized into five clusters based on similarity centers: the Middle East containing Iran, Iraq, Turkey, and Tajikistan; the Far East, including India, Pakistan, and Korea; Europe; the American continent; and Africa, including Egypt, Sudan and Libya. The present study shows the effectiveness of employing the mixture of SRAP and SCoT markers in the identification of safflower genetic diversity that would be useful for conservation and population genetics of safflower improvement in further studies.

Towards new sources of resistance to the currant-lettuce aphid (Nasonovia ribisnigri)

Domesticated lettuce varieties encompass much morphological variation across a range of crop type groups, with large collections of cultivars and landrace accessions maintained in genebanks. Additional variation not captured during domestication, present in ancestral wild relatives, represents a potentially rich source of alleles that can deliver to sustainable crop production. However, these large collections are difficult and costly to screen for many agronomically important traits. In this paper, we describe the generation of a diversity collection of 96 lettuce and wild species accessions that are amenable to routine phenotypic analysis and their genotypic characterization with a panel of 682 newly developed expressed sequence tag (EST)-linked KASP™ single nucleotide polymorphism (SNP) markers that are anchored to the draft Lactuca sativa genome assembly. To exemplify the utility of these resources, we screened the collection for putative sources of resistance to currant-lettuce aphid (Nasonovia ribisnigri) and carried out association analyses to look for potential SNPs linked to resistance.

A Review of Living Collections with Special Emphasis on Sustainability and Its Impact on Research Across Multiple Disciplines

Formal living collections have unique characteristics that distinguish them from other types of biorepositories. Comprising diverse resources, microbe culture collections, crop and biodiversity plant germplasm collections, and animal germplasm repositories are commonly allied with specific research communities or stakeholder groups. Among living collections, microbial culture collections have very long and unique life histories, with some being older than 100 years. Regulatory, financial, and technical developments have impacted living collections in many ways. International treaty obligations and restrictions on release of genetically modified organisms complicate the activities of living collections. Funding for living collections is a continuing challenge and threatens to create a two-tier system where medically relevant collections are well funded and all other collections are underfunded and hence understaffed. Molecular, genetic, and whole genome sequence analysis of contents of microbes and other living resource collections bring additional value to living collections.

The U.S. Culture Collection Network Lays the Foundation for Progress in Preservation of Valuable Microbial Resources

The U.S. Culture Collection Network was formed in 2012 by a group of culture collection scientists and stakeholders in order to continue the progress established previously through efforts of an ad hoc group. The network is supported by a Research Coordination Network grant from the U.S. National Science Foundation (NSF) and has the goals of promoting interaction among collections, encouraging the adoption of best practices, and protecting endangered or orphaned collections. After prior meetings to discuss best practices, shared data, and synergy with genome programs, the network held a meeting at the U.S. Department of Agriculture (USDA)-Agricultural Research Service (ARS) National Center for Genetic Resources Preservation (NCGRP) in Fort Collins, Colorado in October 2015 specifically to discuss collections that are vulnerable because of changes in funding programs, or are at risk of loss because of retirement or lack of funding. The meeting allowed collection curators who had already backed up their resources at the USDA NCGRP to visit the site, and brought collection owners, managers, and stakeholders together. Eight formal collections have established off-site backups with the USDA-ARS, ensuring that key material will be preserved for future research. All of the collections with backup at the NCGRP are public distributing collections including U.S. NSF-supported genetic stock centers, USDA-ARS collections, and university-supported collections. Facing the retirement of several pioneering researchers, the community discussed the value of preserving personal research collections and agreed that a mechanism to preserve these valuable collections was essential to any future national culture collection system. Additional input from curators of plant and animal collections emphasized that collections of every kind face similar challenges in developing long-range plans for sustainability.

Principal components analysis--K-means transposon element based foxtail millet core collection selection method

Background

Core collections are important tools in genetic resources research and administration. At present, most core collection selection criteria are based on one of the following item characteristics: passport data, genetic markers, or morphological traits, which may lead to inadequate representations of variability in the complete collection. The development of a comprehensive methodology that includes as much element data as possible has been explored poorly. Using a collection of (Setaria italica sbsp. italica (L.) P. Beauv.) as a model, we developed a method for core collection construction based on genotype data and numerical representations of agromorphological traits, thereby improving the selection process.

Results

Principal component analysis allows the selection of the most informative discriminators among the various elements evaluated, regardless of whether they are genetic or morphological, thereby providing an adequate criterion for further K-mean clustering. Overall, the core collections of S. italica constructed using only genotype data demonstrated overall better validation scores than other core collections that we generated. However, core collection based on both genotype and agromorphological characteristics represented the overall diversity adequately.

Conclusions

The inclusion of both genotype and agromorphological characteristics as a comprehensive dataset in this methodology ensures that agricultural traits are considered in the core collection construction. This approach will be beneficial for genetic resources management and research activities for S. italica as well as other genetic resources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0343-z) contains supplementary material, which is available to authorized users.

Variations in DREB1A and VP1.1 Genes Show Association with Salt Tolerance Traits in Wild Tomato (Solanum pimpinellifolium)

Association analysis was conducted in a core collection of 94 genotypes of Solanum pimpinellifolium to identify variations linked to salt tolerance traits (physiological and yield traits under salt stress) in four candidate genes viz., DREB1A, VP1.1, NHX1, and TIP. The candidate gene analysis covered a concatenated length of 4594 bp per individual and identified five SNP/Indels in DREB1A and VP1.1 genes explaining 17.0% to 25.8% phenotypic variation for various salt tolerance traits. Out of these five alleles, one at 297 bp in DREB1A had in-frame deletion of 6 bp (CTGCAT) or 12 bp (CTGCATCTGCAT), resulting in two alleles, viz., SpDREB1A_297_6 and SpDREB1A_297_12. These alleles individually or as haplotypes accounted for maximum phenotypic variance of about 25% for various salt tolerance traits. Design of markers for selection of the favorable alleles/haplotypes will hasten marker-assisted introgression of salt tolerance from S. pimpinellifolium into cultivated tomato.

Analysis of genetic diversity and population structure of rice germplasm from north-eastern region of India and development of a core germplasm set

The North-Eastern region (NER) of India, comprising of Arunachal Pradesh, Assam, Manipur, Meghalaya, Mizoram, Nagaland and Tripura, is a hot spot for genetic diversity and the most probable origin of rice. North-east rice collections are known to possess various agronomically important traits like biotic and abiotic stress tolerance, unique grain and cooking quality. The genetic diversity and associated population structure of 6,984 rice accessions, originating from NER, were assessed using 36 genome wide unlinked single nucleotide polymorphism (SNP) markers distributed across the 12 rice chromosomes. All of the 36 SNP loci were polymorphic and bi-allelic, contained five types of base substitutions and together produced nine types of alleles. The polymorphic information content (PIC) ranged from 0.004 for Tripura to 0.375 for Manipur and major allele frequency ranged from 0.50 for Assam to 0.99 for Tripura. Heterozygosity ranged from 0.002 in Nagaland to 0.42 in Mizoram and gene diversity ranged from 0.006 in Arunachal Pradesh to 0.50 in Manipur. The genetic relatedness among the rice accessions was evaluated using an unrooted phylogenetic tree analysis, which grouped all accessions into three major clusters. For determining population structure, populations K = 1 to K = 20 were tested and population K = 3 was present in all the states, with the exception of Meghalaya and Manipur where, K = 5 and K = 4 populations were present, respectively. Principal Coordinate Analysis (PCoA) showed that accessions were distributed according to their population structure. AMOVA analysis showed that, maximum diversity was partitioned at the individual accession level (73% for Nagaland, 58% for Arunachal Pradesh and 57% for Tripura). Using POWERCORE software, a core set of 701 accessions was obtained, which accounted for approximately 10% of the total NE India collections, representing 99.9% of the allelic diversity. The rice core set developed will be a valuable resource for future genomic studies and crop improvement strategies.

Geographical gradient of the eIF4E alleles conferring resistance to potyviruses in pea (Pisum) germplasm

BACKGROUND

The eukaryotic translation initiation factor 4E was shown to be involved in resistance against several potyviruses in plants, including pea. We combined our knowledge of pea germplasm diversity with that of the eIF4E gene to identify novel genetic diversity.

METHODOLOGY/PRINCIPAL FINDINGS

Germplasm of 2803 pea accessions was screened for eIF4E intron 3 length polymorphism, resulting in the detection of four eIF4E(A-B-C-S) variants, whose distribution was geographically structured. The eIF4E(A) variant conferring resistance to the P1 PSbMV pathotype was found in 53 accessions (1.9%), of which 15 were landraces from India, Afghanistan, Nepal, and 7 were from Ethiopia. A newly discovered variant, eIF4E(B), was present in 328 accessions (11.7%) from Ethiopia (29%), Afghanistan (23%), India (20%), Israel (25%) and China (39%). The eIF4E(C) variant was detected in 91 accessions (3.2% of total) from India (20%), Afghanistan (33%), the Iberian Peninsula (22%) and the Balkans (9.3%). The eIF4E(S) variant for susceptibility predominated as the wild type. Sequencing of 73 samples, identified 34 alleles at the whole gene, 26 at cDNA and 19 protein variants, respectively. Fifteen alleles were virologically tested and 9 alleles (eIF4E(A-1-2-3-4-5-6-7), eIF4E(B-1), eIF4E(C-2)) conferred resistance to the P1 PSbMV pathotype.

CONCLUSIONS/SIGNIFICANCE

This work identified novel eIF4E alleles within geographically structured pea germplasm and indicated their independent evolution from the susceptible eIF4E(S1) allele. Despite high variation present in wild Pisum accessions, none of them possessed resistance alleles, supporting a hypothesis of distinct mode of evolution of resistance in wild as opposed to crop species. The Highlands of Central Asia, the northern regions of the Indian subcontinent, Eastern Africa and China were identified as important centers of pea diversity that correspond with the diversity of the pathogen. The series of alleles identified in this study provides the basis to study the co-evolution of potyviruses and the pea host.

A strategy for finding the optimal scale of plant core collection based on Monte Carlo simulation

Core collection is an ideal resource for genome-wide association studies (GWAS). A subcore collection is a subset of a core collection. A strategy was proposed for finding the optimal sampling percentage on plant subcore collection based on Monte Carlo simulation. A cotton germplasm group of 168 accessions with 20 quantitative traits was used to construct subcore collections. Mixed linear model approach was used to eliminate environment effect and GE (genotype × environment) effect. Least distance stepwise sampling (LDSS) method combining 6 commonly used genetic distances and unweighted pair-group average (UPGMA) cluster method was adopted to construct subcore collections. Homogeneous population assessing method was adopted to assess the validity of 7 evaluating parameters of subcore collection. Monte Carlo simulation was conducted on the sampling percentage, the number of traits, and the evaluating parameters. A new method for “distilling free-form natural laws from experimental data” was adopted to find the best formula to determine the optimal sampling percentages. The results showed that coincidence rate of range (CR) was the most valid evaluating parameter and was suitable to serve as a threshold to find the optimal sampling percentage. The principal component analysis showed that subcore collections constructed by the optimal sampling percentages calculated by present strategy were well representative.

Carotenoids gene markers for sweetpotato (Ipomoea batatas L. Lam): applications in genetic mapping, diversity evaluation and cross-species transference

Carotenoids play essential biological roles in plants, and genes involved in the carotenoid biosynthesis pathway are evolutionarily conserved. Orange sweetpotato is an important source of β-carotene, a precursor of vitamin A. In spite of this, only a few research studies have focussed on the molecular aspects of carotenoid genes regarding their specific sequence and structure. In this study, we used published carotenoid gene sequences from Ipomoea and other species for "exon-primed intron-crossing" approaches. Fifteen pairs of primers representing six carotenoid genes were designed for different introns, eleven of which amplified scorable and reproducible alleles. The sequence of PCR products showed high homology to the original ones. Moreover, the structure and sequence of the introns and exons from five carotenoid structural genes were partially defined. Intron length polymorphism and intron single nucleotide polymorphisms were detected in amplified sequences. Marker dosages and allelic segregations were analysed in a mapping population. The developed markers were evaluated in a set of Ipomoeas batatas accessions so as to analyse genetic diversity and conservation applicability. Using CG strategy combined with EPIC-PCR technique, we developed carotenoid gene markers in sweetpotato. We reported the first set of polymorphic Candidate Gene markers for I. batatas, and demonstrated transferability in seven wild Ipomoea species. We described the sequence and structure of carotenoid genes and introduced new information about genomic constitution and allele dosage.

Genotype imputation reference panel selection using maximal phylogenetic diversity

The recent dramatic cost reduction of next-generation sequencing technology enables investigators to assess most variants in the human genome to identify risk variants for complex diseases. However, sequencing large samples remains very expensive. For a study sample with existing genotype data, such as array data from genome-wide association studies, a cost-effective approach is to sequence a subset of the study sample and then to impute the rest of the study sample, using the sequenced subset as a reference panel. The use of such an internal reference panel identifies population-specific variants and avoids the problem of a substantial mismatch in ancestry background between the study population and the reference population. To efficiently select an internal panel, we introduce an idea of phylogenetic diversity from mathematical phylogenetics and comparative genomics. We propose the “most diverse reference panel”, defined as the subset with the maximal “phylogenetic diversity”, thereby incorporating individuals that span a diverse range of genotypes within the sample. Using data both from simulations and from the 1000 Genomes Project, we show that the most diverse reference panel can substantially improve the imputation accuracy compared to randomly selected reference panels, especially for the imputation of rare variants. The improvement in imputation accuracy holds across different marker densities, reference panel sizes, and lengths for the imputed segments. We thus propose a novel strategy for planning sequencing studies on samples with existing genotype data.

Establishing an efficient way to utilize the drought resistance germplasm population in wheat

Drought resistance breeding provides a hopeful way to improve yield and quality of wheat in arid and semiarid regions. Constructing core collection is an efficient way to evaluate and utilize drought-resistant germplasm resources in wheat. In the present research, 1,683 wheat varieties were divided into five germplasm groups (high resistant, HR; resistant, R; moderate resistant, MR; susceptible, S; and high susceptible, HS). The least distance stepwise sampling (LDSS) method was adopted to select core accessions. Six commonly used genetic distances (Euclidean distance, Euclid; Standardized Euclidean distance, Seuclid; Mahalanobis distance, Mahal; Manhattan distance, Manhat; Cosine distance, Cosine; and Correlation distance, Correlation) were used to assess genetic distances among accessions. Unweighted pair-group average (UPGMA) method was used to perform hierarchical cluster analysis. Coincidence rate of range (CR) and variable rate of coefficient of variation (VR) were adopted to evaluate the representativeness of the core collection. A method for selecting the ideal constructing strategy was suggested in the present research. A wheat core collection for the drought resistance breeding programs was constructed by the strategy selected in the present research. The principal component analysis showed that the genetic diversity was well preserved in that core collection.