Genotyping-by-sequencing highlights original diversity patterns within a European collection of 1191 maize flint lines, as compared to the maize USDA genebank

Seedling root system adaptation to water availability during maize domestication and global expansion

The maize root system has been reshaped by indirect selection during global adaptation to new agricultural environments. In this study, we characterized the root systems of more than 9,000 global maize accessions and its wild relatives, defining the geographical signature and genomic basis of variation in seminal root number. We demonstrate that seminal root number has increased during maize domestication followed by a decrease in response to limited water availability in locally adapted varieties. By combining environmental and phenotypic association analyses with linkage mapping, we identified genes linking environmental variation and seminal root number. Functional characterization of the transcription factor ZmHb77 and in silico root modeling provides evidence that reshaping root system architecture by reducing the number of seminal roots and promoting lateral root density is beneficial for the resilience of maize seedlings to drought.

Genetic Diversity within a Collection of Italian Maize Inbred Lines: A Resource for Maize Genomics and Breeding

Genetic diversity is fundamental for studying the complex architecture of the traits of agronomic importance, controlled by major and minor loci. Moreover, well-characterized germplasm collections are essential tools for dissecting and analyzing genetic and phenotypic diversity in crops. A panel of 360 entries, a subset of a larger collection maintained within the GenBank at CREA Bergamo, which includes the inbreds derived from traditional Italian maize open-pollinated (OP) varieties and advanced breeding ones (Elite Inbreds), was analyzed to identify SNP markers using the tGBS® genotyping-by-sequencing technology. A total of 797,368 SNPs were found during the initial analysis. Imputation and filtering processes were carried out based on the percentage of missing data, redundant markers, and rarest allele frequencies, resulting in a final dataset of 15,872 SNP markers for which a physical map position was identified. Using this dataset, the inbred panel was characterized for linkage disequilibrium (LD), genetic diversity, population structure, and genetic relationships. LD decay at a genome-wide level indicates that the collection is a suitable resource for association mapping. Population structure analyses, which were carried out with different clustering methods, showed stable grouping statistics for four groups, broadly corresponding to 'Insubria', 'Microsperma', and 'Scagliolino' genotypes, with a fourth group composed prevalently of elite accessions derived from Italian and US breeding programs. Based on these results, the CREA Italian maize collection, genetically characterized in this study, can be considered an important tool for the mapping and characterization of useful traits and associated loci/alleles, to be used in maize breeding programs.

Assessing the potential of genetic resource introduction into elite germplasm: a collaborative multiparental population for flint maize

KEY MESSAGE

Implementing a collaborative pre-breeding multi-parental population efficiently identifies promising donor x elite pairs to enrich the flint maize elite germplasm. Genetic diversity is crucial for maintaining genetic gains and ensuring breeding programs' long-term success. In a closed breeding program, selection inevitably leads to a loss of genetic diversity. While managing diversity can delay this loss, introducing external sources of diversity is necessary to bring back favorable genetic variation. Genetic resources exhibit greater diversity than elite materials, but their lower performance levels hinder their use. This is the case for European flint maize, for which elite germplasm has incorporated only a limited portion of the diversity available in landraces. To enrich the diversity of this elite genetic pool, we established an original cooperative maize bridging population that involves crosses between private elite materials and diversity donors to create improved genotypes that will facilitate the incorporation of original favorable variations. Twenty donor × elite BC1S2 families were created and phenotyped for hybrid value for yield related traits. Crosses showed contrasted means and variances and therefore contrasted potential in terms of selection as measured by their usefulness criterion (UC). Average expected mean performance gain over the initial elite material was 5%. The most promising donor for each elite line was identified. Results also suggest that one more generation, i.e., 3 in total, of crossing to the elite is required to fully exploit the potential of a donor. Altogether, our results support the usefulness of incorporating genetic resources into elite flint maize. They call for further effort to create fixed diversity donors and identify those most suitable for each elite program.

Genetic diversity and selection signatures in a gene bank panel of maize inbred lines from Southeast Europe compared with two West European panels

Southeast Europe (SEE) is a very important maize-growing region, comparable to the Corn belt region of the United States, with similar dent germplasm (dent by dent hybrids). Historically, this region has undergone several genetic material swaps, following the trends in the US, with one of the most significant swaps related to US aid programs after WWII. The imported accessions used to make double-cross hybrids were also mixed with previously adapted germplasm originating from several more distant OPVs, supporting the transition to single cross-breeding. Many of these materials were deposited at the Maize Gene Bank of the Maize Research Institute Zemun Polje (MRIZP) between the 1960s and 1980s. A part of this Gene Bank (572 inbreds) was genotyped with Affymetrix Axiom Maize Genotyping Array with 616,201 polymorphic variants. Data were merged with two other genotyping datasets with mostly European flint (TUM dataset) and dent (DROPS dataset) germplasm. The final pan-European dataset consisted of 974 inbreds and 460,243 markers. Admixture analysis showed seven ancestral populations representing European flint, B73/B14, Lancaster, B37, Wf9/Oh07, A374, and Iodent pools. Subpanel of inbreds with SEE origin showed a lack of Iodent germplasm, marking its historical context. Several signatures of selection were identified at chromosomes 1, 3, 6, 7, 8, 9, and 10. The regions under selection were mined for protein-coding genes and were used for gene ontology (GO) analysis, showing a highly significant overrepresentation of genes involved in response to stress. Our results suggest the accumulation of favorable allelic diversity, especially in the context of changing climate in the genetic resources of SEE.

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.

Construction of a high density genetic linkage map to define the locus conferring seedlessness from Mukaku Kishu mandarin

Mukaku Kishu ('MK'), a small sized mandarin, is an important source of seedlessness in citrus breeding. Identification and mapping the gene(s) governing 'MK' seedlessness will expedite seedless cultivar development. In this study, two 'MK'-derived mapping populations- LB8-9 Sugar Belle^® ('SB') × 'MK' (N=97) and Daisy ('D') × 'MK' (N=68) were genotyped using an Axiom_Citrus56 Array encompassing 58,433 SNP probe sets, and population specific male and female parent linkage maps were constructed. The parental maps of each population were integrated to produce sub-composite maps, which were further merged to develop a consensus linkage map. All the parental maps (except 'MK_D') had nine major linkage groups, and contained 930 ('SB'), 810 ('MK_SB'), 776 ('D') and 707 ('MK_D') SNPs. The linkage maps displayed 96.9 ('MK_D') to 98.5% ('SB') chromosomal synteny with the reference Clementine genome. The consensus map was comprised of 2588 markers including a phenotypic seedless (Fs)-locus and spanned a genetic distance of 1406.84 cM, with an average marker distance of 0.54 cM, which is substantially lower than the reference Clementine map. For the phenotypic Fs-locus, the distribution of seedy and seedless progenies in both 'SB' × 'MK' (55:42, χ² = 1.74) and 'D' × 'MK' populations (33:35, χ² = 0.06) followed a test cross pattern. The Fs-locus mapped on chromosome 5 with SNP marker 'AX-160417325' at 7.4 cM in 'MK_SB' map and between two SNP markers 'AX-160536283' and 'AX-160906995' at a distance of 2.4 and 4.9 cM, respectively in 'MK_D' map. The SNPs 'AX-160417325' and 'AX-160536283' correctly predicted seedlessness of 25-91.9% progenies in this study. Based on the alignment of flanking SNP markers to the Clementine reference genome, the candidate gene for seedlessness hovered in a ~ 6.0 Mb region between 3.97 Mb (AX-160906995) to 10.00 Mb (AX-160536283). This region has 131 genes of which 13 genes (belonging to seven gene families) reportedly express in seed coat or developing embryo. The findings of the study will prove helpful in directing future research for fine mapping this region and eventually underpinning the exact causative gene governing seedlessness in 'MK'.

Genetic diversity, population structure, and selection of breeder germplasm subsets from the USDA sweetpotato (Ipomoea batatas) collection

Sweetpotato (Ipomoea batatas) is the sixth most important food crop and plays a critical role in maintaining food security worldwide. Support for sweetpotato improvement research in breeding and genetics programs, and maintenance of sweetpotato germplasm collections is essential for preserving food security for future generations. Germplasm collections seek to preserve phenotypic and genotypic diversity through accession characterization. However, due to its genetic complexity, high heterogeneity, polyploid genome, phenotypic plasticity, and high flower production variability, sweetpotato genetic characterization is challenging. Here, we characterize the genetic diversity and population structure of 604 accessions from the sweetpotato germplasm collection maintained by the United States Department of Agriculture (USDA), Agricultural Research Service (ARS), Plant Genetic Resources Conservation Unit (PGRCU) in Griffin, Georgia, United States. Using the genotyping-by-sequencing platform (GBSpoly) and bioinformatic pipelines (ngsComposer and GBSapp), a total of 102,870 polymorphic SNPs with hexaploid dosage calls were identified from the 604 accessions. Discriminant analysis of principal components (DAPC) and Bayesian clustering identified six unique genetic groupings across seven broad geographic regions. Genetic diversity analyses using the hexaploid data set revealed ample genetic diversity among the analyzed collection in concordance with previous analyses. Following population structure and diversity analyses, breeder germplasm subsets of 24, 48, 96, and 384 accessions were established using K-means clustering with manual selection to maintain phenotypic and genotypic diversity. The genetic characterization of the PGRCU sweetpotato germplasm collection and breeder germplasm subsets developed in this study provide the foundation for future association studies and serve as precursors toward phenotyping studies aimed at linking genotype with phenotype.

Population Genomics Study and Implications for the Conservation of Zabelia tyaihyonii Based on Genotyping-By-Sequencing

Zabelia tyaihyonii (Nakai) Hisauti and H. Hara is a perennial shrub endemic to Republic of Korea that grows naturally in only a very limited region of the dolomite areas of Gangwon-do and Chungcheongbuk-do Provinces in the Republic of Korea. Given its geographical characteristics, it is more vulnerable than more widely distributed species. Despite the need for comprehensive information to support conservation, population genetic information for this species is very scarce. In this study, we analyzed the genetic diversity and population structure of 94 individuals from six populations of Z. tyaihyonii using a genotyping-by-sequencing (GBS) approach to provide important information for proper conservation and management. Our results, based on 3088 single nucleotide polymorphisms (SNPs), showed a mean expected heterozygosity (He) of 0.233, no sign of within-population inbreeding (G_IS that was close to or even below zero in all populations), and a high level of genetic differentiation (F_ST = 0.170). Analysis of molecular variance (AMOVA) indicated that the principal molecular variance existed within populations (84.5%) rather than among populations (17.0%). We suggested that six management units were proposed for conservation considering Bayesian structure analysis and phylogenetic analysis, and given the various current situations faced by Z. tyaihyonii, it is believed that not only the in situ conservation but also the ex situ conservation should be considered.

Genetic Diversity, Population Structure and Linkage Disequilibrium Analyses in Tropical Maize Using Genotyping by Sequencing

Several maize breeding programs in India have developed numerous inbred lines but the lines have not been characterized using high-density molecular markers. Here, we studied the molecular diversity, population structure, and linkage disequilibrium (LD) patterns in a panel of 314 tropical normal corn, two sweet corn, and six popcorn inbred lines developed by 17 research centers in India, and 62 normal corn from the International Maize and Wheat Improvement Center (CIMMYT). The 384 inbred lines were genotyped with 60,227 polymorphic single nucleotide polymorphisms (SNPs). Most of the pair-wise relative kinship coefficients (58.5%) were equal or close to 0, which suggests the lack of redundancy in the genomic composition in the majority of inbred lines. Genetic distance among most pairs of lines (98.3%) varied from 0.20 to 0.34 as compared with just 1.7% of the pairs of lines that differed by <0.20, which suggests greater genetic variation even among sister lines. The overall average of 17% heterogeneity was observed in the panel indicated the need for further inbreeding in the high heterogeneous genotypes. The mean nucleotide diversity and frequency of polymorphic sites observed in the panel were 0.28 and 0.02, respectively. The model-based population structure, principal component analysis, and phylogenetic analysis revealed three to six groups with no clear patterns of clustering by centers-wise breeding lines, types of corn, kernel characteristics, maturity, plant height, and ear placement. However, genotypes were grouped partially based on their source germplasm from where they derived.

Deep scoping: a breeding strategy to preserve, reintroduce and exploit genetic variation

The deep scoping method incorporates the use of a gene bank together with different population layers to reintroduce genetic variation into the breeding population, thus maximizing the long-term genetic gain without reducing the short-term genetic gain or increasing the total financial cost. Genomic prediction is often combined with truncation selection to identify superior parental individuals that can pass on favorable quantitative trait locus (QTL) alleles to their offspring. However, truncation selection reduces genetic variation within the breeding population, causing a premature convergence to a sub-optimal genetic value. In order to also increase genetic gain in the long term, different methods have been proposed that better preserve genetic variation. However, when the genetic variation of the breeding population has already been reduced as a result of prior intensive selection, even those methods will not be able to avert such premature convergence. Pre-breeding provides a solution for this problem by reintroducing genetic variation into the breeding population. Unfortunately, as pre-breeding often relies on a separate breeding population to increase the genetic value of wild specimens before introducing them in the elite population, it comes with an increased financial cost. In this paper, on the basis of a simulation study, we propose a new method that reintroduces genetic variation in the breeding population on a continuous basis without the need for a separate pre-breeding program or a larger population size. This way, we are able to introduce favorable QTL alleles into an elite population and maximize the genetic gain in the short as well as in the long term without increasing the financial cost.

Genetic Identity, Diversity, and Population Structure of CIP's Sweetpotato (I. batatas) Germplasm Collection

The in trust sweetpotato collection housed by the International Center of Potato (CIP) is one of the largest assemblages of plant material representing the genetic resources of this important staple crop. The collection currently contains almost 6,000 accessions of Ipomoea batatas (cultivated sweetpotato) and over 1,000 accessions of sweetpotato crop wild relatives (CWRs). In this study, the entire cultivated collection (5,979 accessions) was genotyped with a panel of 20 simple sequence repeat (SSR) markers to assess genetic identity, diversity, and population structure. Genotyping and phenotyping of in vitro plantlets and mother plants were conducted simultaneously on 2,711 accessions (45% of the total collection) to identify and correct possible genetic identity errors which could have occurred at any time over the thirty plus years of maintenance in the in vitro collection. Within this group, 533 accessions (19.6%) had errors in identity. Field evaluations of morphological descriptors were carried out to confirm the marker data. A phylogenetic tree was constructed to reveal the intraspecific relationships in the population which uncovered high levels of redundancy in material from Peru and Latin America. These genotypic data were supported by morphological data. Population structure analysis demonstrated support for four ancestral populations with many of the accessions having lower levels of gene flow from the other populations. This was especially true of germplasm derived from Peru, Ecuador, and Africa. The set of 20 SSR markers was subsequently utilized to examine a subset of 189 accessions from the USDA sweetpotato germplasm collection and to identify and reconcile potential errors in the identification of clones shared between these collections. Marker analysis demonstrated that the USDA subset of material had 65 unique accessions that were not found in the larger CIP collection. As far as the authors are aware, this is the first report of genotyping an entire sweetpotato germplasm collection in its entirety.

A worldwide maize panel revealed new genetic variation for cold tolerance

A large association panel of 836 maize inbreds revealed a broader genetic diversity of cold tolerance, as predominantly favorable QTL with small effects were identified, indicating that genomic selection is the most promising option for breeding maize for cold tolerance. Maize (Zea mays L.) has limited cold tolerance, and breeding for cold tolerance is a noteworthy bottleneck for reaching the high potential of maize production in temperate areas. In this study, we evaluate a large panel of 836 maize inbred lines to detect genetic loci and candidate genes for cold tolerance at the germination and seedling stages. Genetic variation for cold tolerance was larger than in previous reports with moderately high heritability for most traits. We identified 187 significant single-nucleotide polymorphisms (SNPs) that were integrated into 159 quantitative trait loci (QTL) for emergence and traits related to early growth. Most of the QTL have small effects and are specific for each environment, with the majority found under control conditions. Favorable alleles are more frequent in 120 inbreds including all germplasm groups, but mainly from Minnesota and Spain. Therefore, there is a large, potentially novel, genetic variability in the germplasm groups represented by these inbred lines. Most of the candidate genes are involved in metabolic processes and intracellular membrane-bounded organelles. We expect that further evaluations of germplasm with broader genetic diversity could identify additional favorable alleles for cold tolerance. However, it is not likely that further studies will find favorable alleles with large effects for improving cold tolerance in maize.

Discovery of beneficial haplotypes for complex traits in maize landraces

Genetic variation is of crucial importance for crop improvement. Landraces are valuable sources of diversity, but for quantitative traits efficient strategies for their targeted utilization are lacking. Here, we map haplotype-trait associations at high resolution in ~1000 doubled-haploid lines derived from three maize landraces to make their native diversity for early development traits accessible for elite germplasm improvement. A comparative genomic analysis of the discovered haplotypes in the landrace-derived lines and a panel of 65 breeding lines, both genotyped with 600k SNPs, points to untapped beneficial variation for target traits in the landraces. The superior phenotypic performance of lines carrying favorable landrace haplotypes as compared to breeding lines with alternative haplotypes confirms these findings. Stability of haplotype effects across populations and environments as well as their limited effects on undesired traits indicate that our strategy has high potential for harnessing beneficial haplotype variation for quantitative traits from genetic resources.

Genetic variations present in landraces are critical for crop genetic improvement. Here, the authors map haplotype-trait associations in ~1000 doubled haploid lines derived from three European maize landraces and identify beneficial haplotypes for quantitative traits that are not present in breeding lines.

Applications of genotyping-by-sequencing (GBS) in maize genetics and breeding

Genotyping-by-Sequencing (GBS) is a low-cost, high-throughput genotyping method that relies on restriction enzymes to reduce genome complexity. GBS is being widely used for various genetic and breeding applications. In the present study, 2240 individuals from eight maize populations, including two association populations (AM), backcross first generation (BC1), BC1F2, F2, double haploid (DH), intermated B73 × Mo17 (IBM), and a recombinant inbred line (RIL) population, were genotyped using GBS. A total of 955,120 of raw data for SNPs was obtained for each individual, with an average genotyping error of 0.70%. The rate of missing genotypic data for these SNPs was related to the level of multiplex sequencing: ~ 25% missing data for 96-plex and ~ 55% for 384-plex. Imputation can greatly reduce the rate of missing genotypes to 12.65% and 3.72% for AM populations and bi-parental populations, respectively, although it increases total genotyping error. For analysis of genetic diversity and linkage mapping, unimputed data with a low rate of genotyping error is beneficial, whereas, for association mapping, imputed data would result in higher marker density and would improve map resolution. Because imputation does not influence the prediction accuracy, both unimputed and imputed data can be used for genomic prediction. In summary, GBS is a versatile and efficient SNP discovery approach for homozygous materials and can be effectively applied for various purposes in maize genetics and breeding.

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

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

Genomic prediction with a maize collaborative panel: identification of genetic resources to enrich elite breeding programs

Collaborative diversity panels and genomic prediction seem relevant to identify and harness genetic resources for polygenic trait-specific enrichment of elite germplasms. In plant breeding, genetic diversity is important to maintain the pace of genetic gain and the ability to respond to new challenges in a context of climatic and social expectation changes. Many genetic resources are accessible to breeders but cannot all be considered for broadening the genetic diversity of elite germplasm. This study presents the use of genomic predictions trained on a collaborative diversity panel, which assembles genetic resources and elite lines, to identify resources to enrich an elite germplasm. A maize collaborative panel (386 lines) was considered to estimate genome-wide marker effects. Relevant predictive abilities (0.40-0.55) were observed on a large population of private elite materials, which supported the interest of such a collaborative panel for diversity management perspectives. Grain-yield estimated marker effects were used to select a donor that best complements an elite recipient at individual loci or haplotype segments, or that is expected to give the best-performing progeny with the elite. Among existing and new criteria that were compared, some gave more weight to the donor-elite complementarity than to the donor value, and appeared more adapted to long-term objective. We extended this approach to the selection of a set of donors complementing an elite population. We defined a crossing plan between identified donors and elite recipients. Our results illustrated how collaborative projects based on diversity panels including both public resources and elite germplasm can contribute to a better characterization of genetic resources in view of their use to enrich elite germplasm.

High throughput genotyping of structural variations in a complex plant genome using an original Affymetrix® axiom® array

Background

Insertions/deletions (InDels) and more specifically presence/absence variations (PAVs) are pervasive in several species and have strong functional and phenotypic effect by removing or drastically modifying genes. Genotyping of such variants on large panels remains poorly addressed, while necessary for approaches such as association mapping or genomic selection.

Results

We have developed, as a proof of concept, a new high-throughput and affordable approach to genotype InDels. We first identified 141,000 InDels by aligning reads from the B73 line against the genome of three temperate maize inbred lines (F2, PH207, and C103) and reciprocally. Next, we designed an Affymetrix® Axiom® array to target these InDels, with a combination of probes selected at breakpoint sites (13%) or within the InDel sequence, either at polymorphic (25%) or non-polymorphic sites (63%) sites. The final array design is composed of 662,772 probes and targets 105,927 InDels, including PAVs ranging from 35 bp to 129kbp. After Affymetrix® quality control, we successfully genotyped 86,648 polymorphic InDels (82% of all InDels interrogated by the array) on 445 maize DNA samples with 422,369 probes. Genotyping InDels using this approach produced a highly reliable dataset, with low genotyping error (~ 3%), high call rate (~ 98%), and high reproducibility (> 95%). This reliability can be further increased by combining genotyping of several probes calling the same InDels (< 0.1% error rate and > 99.9% of call rate for 5 probes). This “proof of concept” tool was used to estimate the kinship matrix between 362 maize lines with 57,824 polymorphic InDels. This InDels kinship matrix was highly correlated with kinship estimated using SNPs from Illumina 50 K SNP arrays.

Conclusions

We efficiently genotyped thousands of small to large InDels on a sizeable number of individuals using a new Affymetrix® Axiom® array. This powerful approach opens the way to studying the contribution of InDels to trait variation and heterosis in maize. The approach is easily extendable to other species and should contribute to decipher the biological impact of InDels at a larger scale.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genotyping-by-sequencing and SNP-arrays are complementary for detecting quantitative trait loci by tagging different haplotypes in association studies

BACKGROUND

Single Nucleotide Polymorphism (SNP) array and re-sequencing technologies have different properties (e.g. calling rate, minor allele frequency profile) and drawbacks (e.g. ascertainment bias). This lead us to study their complementarity and the consequences of using them separately or combined in diversity analyses and Genome-Wide Association Studies (GWAS). We performed GWAS on three traits (grain yield, plant height and male flowering time) measured in 22 environments on a panel of 247 F1 hybrids obtained by crossing 247 diverse dent maize inbred lines with a same flint line. The 247 lines were genotyped using three genotyping technologies (Genotyping-By-Sequencing, Illumina Infinium 50 K and Affymetrix Axiom 600 K arrays).

RESULTS

The effects of ascertainment bias of the 50 K and 600 K arrays were negligible for deciphering global genetic trends of diversity and for estimating relatedness in this panel. We developed an original approach based on linkage disequilibrium (LD) extent in order to determine whether SNPs significantly associated with a trait and that are physically linked should be considered as a single Quantitative Trait Locus (QTL) or several independent QTLs. Using this approach, we showed that the combination of the three technologies, which have different SNP distributions and densities, allowed us to detect more QTLs (gain in power) and potentially refine the localization of the causal polymorphisms (gain in resolution).

CONCLUSIONS

Conceptually different technologies are complementary for detecting QTLs by tagging different haplotypes in association studies. Considering LD, marker density and the combination of different technologies (SNP-arrays and re-sequencing), the genotypic data available were most likely enough to well represent polymorphisms in the centromeric regions, whereas using more markers would be beneficial for telomeric regions.

Testcross performance of doubled haploid lines from European flint maize landraces is promising for broadening the genetic base of elite germplasm

Selected doubled haploid lines averaged similar testcross performance as their original landraces, and the best of them approached the yields of elite inbreds, demonstrating their potential to broaden the narrow genetic diversity of the flint germplasm pool. Maize landraces represent a rich source of genetic diversity that remains largely idle because the high genetic load and performance gap to elite germplasm hamper their use in modern breeding programs. Production of doubled haploid (DH) lines can mitigate problems associated with the use of landraces in pre-breeding. Our objective was to assess in comparison with modern materials the testcross performance (TP) of the best 89 out of 389 DH lines developed from six landraces and evaluated in previous studies for line per se performance (LP). TP with a dent tester was evaluated for the six original landraces, ~ 15 DH lines from each landrace selected for LP, and six elite flint inbreds together with nine commercial hybrids for grain and silage traits. Mean TP of the DH lines rarely differed significantly from TP of their corresponding landrace, which averaged in comparison with the mean TP of the elite flint inbreds ~ 20% lower grain yield and ~ 10% lower dry matter and methane yield. Trait correlations of DH lines closely agreed with the literature; correlation of TP with LP was zero for grain yield, underpinning the need to evaluate TP in addition to LP. For all traits, we observed substantial variation for TP among the DH lines and the best showed similar TP yields as the elite inbreds. Our results demonstrate the high potential of landraces for broadening the narrow genetic base of the flint heterotic pool and the usefulness of the DH technology for exploiting idle genetic resources from gene banks.

Stepwise cis-Regulatory Changes in ZCN8 Contribute to Maize Flowering-Time Adaptation

Maize (Zea mays ssp. mays) was domesticated in southwestern Mexico ∼9,000 years ago from its wild ancestor, teosinte (Zea mays ssp. parviglumis) [1]. From its center of origin, maize experienced a rapid range expansion and spread over 90° of latitude in the Americas [2-4], which required a novel flowering-time adaptation. ZEA CENTRORADIALIS 8 (ZCN8) is the maize florigen gene and has a central role in mediating flowering [5, 6]. Here, we show that ZCN8 underlies a major quantitative trait locus (QTL) (qDTA8) for flowering time that was consistently detected in multiple maize-teosinte experimental populations. Through association analysis in a large diverse panel of maize inbred lines, we identified a SNP (SNP-1245) in the ZCN8 promoter that showed the strongest association with flowering time. SNP-1245 co-segregated with qDTA8 in maize-teosinte mapping populations. We demonstrate that SNP-1245 is associated with differential binding by the flowering activator ZmMADS1. SNP-1245 was a target of selection during early domestication, which drove the pre-existing early flowering allele to near fixation in maize. Interestingly, we detected an independent association block upstream of SNP-1245, wherein the early flowering allele that most likely originated from Zea mays ssp. mexicana introgressed into the early flowering haplotype of SNP-1245 and contributed to maize adaptation to northern high latitudes. Our study demonstrates how independent cis-regulatory variants at a gene can be selected at different evolutionary times for local adaptation, highlighting how complex cis-regulatory control mechanisms evolve. Finally, we propose a polygenic map for the pre-Columbian spread of maize throughout the Americas.