Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants

BACKGROUND

Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated.

RESULTS

De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups.

CONCLUSIONS

We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Predicting genomic selection efficiency to optimize calibration set and to assess prediction accuracy in highly structured populations

KEY MESSAGE

We propose a criterion to predict genomic selection efficiency for structured populations. This criterion is useful to define optimal calibration set and to estimate prediction reliability for multiparental populations. Genomic selection refers to the use of genotypic information for predicting the performance of selection candidates. It has been shown that prediction accuracy depends on various parameters including the composition of the calibration set (CS). Assessing the level of accuracy of a given prediction scenario is of highest importance because it can be used to optimize CS sampling before collecting phenotypes, and once the breeding values are predicted it informs the breeders about the reliability of these predictions. Different criteria were proposed to optimize CS sampling in highly diverse panels, which can be useful to screen collections of genotypes. But plant breeders often work on structured material such as biparental or multiparental populations, for which these criteria are less adapted. We derived from the generalized coefficient of determination (CD) theory different criteria to optimize CS sampling and to assess the reliability associated to predictions in structured populations. These criteria were evaluated on two nested association mapping (NAM) populations and two highly diverse panels of maize. They were efficient to sample optimized CS in most situations. They could also estimate at least partly the reliability associated to predictions between NAM families, but they could not estimate differences in the reliability associated to the predictions of NAM families using the highly diverse panels as calibration sets. We illustrated that the CD criteria could be adapted to various prediction scenarios including inter and intra-family predictions, resulting in higher prediction accuracies.

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

Genotyping by sequencing is suitable for analysis of global diversity in maize. We showed the distinctiveness of flint maize inbred lines of interest to enrich the diversity of breeding programs. Genotyping-by-sequencing (GBS) is a highly cost-effective procedure that permits the analysis of large collections of inbred lines. We used it to characterize diversity in 1191 maize flint inbred lines from the INRA collection, the European Cornfed association panel, and lines recently derived from landraces. We analyzed the properties of GBS data obtained with different imputation methods, through comparison with a 50 K SNP array. We identified seven ancestral groups within the Flint collection (dent, Northern flint, Italy, Pyrenees-Galicia, Argentina, Lacaune, Popcorn) in agreement with breeding knowledge. Analysis highlighted many crosses between different origins and the improvement of flint germplasm with dent germplasm. We performed association studies on different agronomic traits, revealing SNPs associated with cob color, kernel color, and male flowering time variation. We compared the diversity of both our collection and the USDA collection which has been previously analyzed by GBS. The population structure of the 4001 inbred lines confirmed the influence of the historical inbred lines (B73, A632, Oh43, Mo17, W182E, PH207, and Wf9) within the dent group. It showed distinctly different tropical and popcorn groups, a sweet-Northern flint group and a flint group sub-structured in Italian and European flint (Pyrenees-Galicia and Lacaune) groups. Interestingly, we identified several selective sweeps between dent, flint, and tropical inbred lines that co-localized with SNPs associated with flowering time variation. The joint analysis of collections by GBS offers opportunities for a global diversity analysis of maize inbred lines.

Genetic diversity and population structure of native maize populations in Latin America and the Caribbean

This study describes the genetic diversity and population structure of 194 native maize populations from 23 countries of Latin America and the Caribbean. The germplasm, representing 131 distinct landraces, was genetically characterized as population bulks using 28 SSR markers. Three main groups of maize germplasm were identified. The first, the Mexico and Southern Andes group, highlights the Pre-Columbian and modern exchange of germplasm between North and South America. The second group, Mesoamerica lowland, supports the hypothesis that two separate human migration events could have contributed to Caribbean maize germplasm. The third, the Andean group, displayed early introduction of maize into the Andes, with little mixing since then, other than a regional interchange zone active in the past. Events and activities in the pre- and post-Columbian Americas including the development and expansion of pre-Columbian cultures and the arrival of Europeans to the Americas are discussed in relation to the history of maize migration from its point of domestication in Mesoamerica to South America and the Caribbean through sea and land routes.

Independent introductions and admixtures have contributed to adaptation of European maize and its American counterparts

Through the local selection of landraces, humans have guided the adaptation of crops to a vast range of climatic and ecological conditions. This is particularly true of maize, which was domesticated in a restricted area of Mexico but now displays one of the broadest cultivated ranges worldwide. Here, we sequenced 67 genomes with an average sequencing depth of 18x to document routes of introduction, admixture and selective history of European maize and its American counterparts. To avoid the confounding effects of recent breeding, we targeted germplasm (lines) directly derived from landraces. Among our lines, we discovered 22,294,769 SNPs and between 0.9% to 4.1% residual heterozygosity. Using a segmentation method, we identified 6,978 segments of unexpectedly high rate of heterozygosity. These segments point to genes potentially involved in inbreeding depression, and to a lesser extent to the presence of structural variants. Genetic structuring and inferences of historical splits revealed 5 genetic groups and two independent European introductions, with modest bottleneck signatures. Our results further revealed admixtures between distinct sources that have contributed to the establishment of 3 groups at intermediate latitudes in North America and Europe. We combined differentiation- and diversity-based statistics to identify both genes and gene networks displaying strong signals of selection. These include genes/gene networks involved in flowering time, drought and cold tolerance, plant defense and starch properties. Overall, our results provide novel insights into the evolutionary history of European maize and highlight a major role of admixture in environmental adaptation, paralleling recent findings in humans.

General and specific combining abilities in a maize (Zea mays L.) test-cross hybrid panel: relative importance of population structure and genetic divergence between parents

General and specific combining abilities of maize hybrids between 288 inbred lines and three tester lines were highly related to population structure and genetic distance inferred from SNP data. Many studies have attempted to provide reliable and quick methods to identify promising parental lines and combinations in hybrid breeding programs. Since the 1950s, maize germplasm has been organized into heterotic groups to facilitate the exploitation of heterosis. Molecular markers have proven efficient tools to address the organization of genetic diversity and the relationship between lines or populations. The aim of the present work was to investigate to what extent marker-based evaluations of population structure and genetic distance may account for general (GCA) and specific (SCA) combining ability components in a population composed of 800 inter and intra-heterotic group hybrids obtained by crossing 288 inbred lines and three testers. Our results illustrate a strong effect of groups identified by population structure analysis on both GCA and SCA components. Including genetic distance between parental lines of hybrids in the model leads to a significant decrease of SCA variance component and an increase in GCA variance component for all the traits. The latter suggests that this approach can be efficient to better estimate the potential combining ability of inbred lines when crossed with unrelated lines, and limits the consequences of tester choice. Significant residual GCA and SCA variance components of models taking into account structure and/or genetic distance highlight the variation available for breeding programs within structure groups.

Estimation of the relatedness coefficients from biallelic markers, application in plant mating designs

The problem of inferring the relatedness distribution between two individuals from biallelic marker data is considered. This problem can be cast as an estimation task in a mixture model: at each marker the latent variable is the relatedness state, and the observed variable is the genotype of the two individuals. In this model, only the prior proportions are unknown, and can be obtained via ML estimation using the EM algorithm. When the markers are biallelic and the data unphased, the identifiability of the model is known not to be guaranteed. In this article, model identifiability is investigated in the case of phased data generated from a crossing design, a classical situation in plant genetics. It is shown that identifiability can be guaranteed under some conditions on the crossing design. The adapted ML estimator is implemented in an R package called Relatedness. The performance of the ML estimator is evaluated and compared to that of the benchmark moment estimator, both on simulated and real data. Compared to its competitor, the ML estimator is shown to be more robust and to provide more realistic estimates.

Association mapping for phenology and plant architecture in maize shows higher power for developmental traits compared with growth influenced traits

Plant architecture, phenology and yield components of cultivated plants have repeatedly been shaped by selection to meet human needs and adaptation to different environments. Here we assessed the genetic architecture of 24 correlated maize traits that interact during plant cycle. Overall, 336 lines were phenotyped in a network of 9 trials and genotyped with 50K single-nucleotide polymorphisms. Phenology was the main factor of differentiation between genetic groups. Then yield components distinguished dents from lower yielding genetic groups. However, most of trait variation occurred within group and we observed similar overall and within group correlations, suggesting a major effect of pleiotropy and/or linkage. We found 34 quantitative trait loci (QTLs) for individual traits and six for trait combinations corresponding to PCA coordinates. Among them, only five were pleiotropic. We found a cluster of QTLs in a 5 Mb region around Tb1 associated with tiller number, ear row number and the first PCA axis, the latter being positively correlated to flowering time and negatively correlated to yield. Kn1 and ZmNIP1 were candidate genes for tillering, ZCN8 for leaf number and Rubisco Activase 1 for kernel weight. Experimental repeatabilities, numbers of QTLs and proportion of explained variation were higher for traits related to plant development such as tillering, leaf number and flowering time, than for traits affected by growth such as yield components. This suggests a simpler genetic determinism with larger individual QTL effects for the first category.

Genome-Wide Analysis of Yield in Europe: Allelic Effects Vary with Drought and Heat Scenarios

Assessing the genetic variability of plant performance under heat and drought scenarios can contribute to reduce the negative effects of climate change. We propose here an approach that consisted of (1) clustering time courses of environmental variables simulated by a crop model in current (35 years × 55 sites) and future conditions into six scenarios of temperature and water deficit as experienced by maize (Zea mays L.) plants; (2) performing 29 field experiments in contrasting conditions across Europe with 244 maize hybrids; (3) assigning individual experiments to scenarios based on environmental conditions as measured in each field experiment; frequencies of temperature scenarios in our experiments corresponded to future heat scenarios (+5°C); (4) analyzing the genetic variation of plant performance for each environmental scenario. Forty-eight quantitative trait loci (QTLs) of yield were identified by association genetics using a multi-environment multi-locus model. Eight and twelve QTLs were associated to tolerances to heat and drought stresses because they were specific to hot and dry scenarios, respectively, with low or even negative allelic effects in favorable scenarios. Twenty-four QTLs improved yield in favorable conditions but showed nonsignificant effects under stress; they were therefore associated with higher sensitivity. Our approach showed a pattern of QTL effects expressed as functions of environmental variables and scenarios, allowing us to suggest hypotheses for mechanisms and candidate genes underlying each QTL. It can be used for assessing the performance of genotypes and the contribution of genomic regions under current and future stress situations and to accelerate breeding for drought-prone environments.

Fortune telling: metabolic markers of plant performance

BACKGROUND

In the last decade, metabolomics has emerged as a powerful diagnostic and predictive tool in many branches of science. Researchers in microbes, animal, food, medical and plant science have generated a large number of targeted or non-targeted metabolic profiles by using a vast array of analytical methods (GC-MS, LC-MS, 1H-NMR….). Comprehensive analysis of such profiles using adapted statistical methods and modeling has opened up the possibility of using single or combinations of metabolites as markers. Metabolic markers have been proposed as proxy, diagnostic or predictors of key traits in a range of model species and accurate predictions of disease outbreak frequency, developmental stages, food sensory evaluation and crop yield have been obtained.

AIM OF REVIEW

(i) To provide a definition of plant performance and metabolic markers, (ii) to highlight recent key applications involving metabolic markers as tools for monitoring or predicting plant performance, and (iii) to propose a workable and cost-efficient pipeline to generate and use metabolic markers with a special focus on plant breeding.

KEY MESSAGE

Using examples in other models and domains, the review proposes that metabolic markers are tending to complement and possibly replace traditional molecular markers in plant science as efficient estimators of performance.

Association mapping for cold tolerance in two large maize inbred panels

BACKGROUND

Breeding for cold tolerance in maize promises to allow increasing growth area and production in temperate zones. The objective of this research was to conduct genome-wide association analyses (GWAS) in temperate maize inbred lines and to find strategies for pyramiding genes for cold tolerance. Two panels of 306 dent and 292 European flint maize inbred lines were evaluated per se and in testcrosses under cold and control conditions in a growth chamber. We recorded indirect measures for cold tolerance as the traits number of days from sowing to emergence, relative leaf chlorophyll content or quantum efficiency of photosystem II. Association mapping for identifying genes associated to cold tolerance in both panels was based on genotyping with 49,585 genome-wide single nucleotide polymorphism (SNP) markers.

RESULTS

We found 275 significant associations, most of them in the inbreds evaluated per se, in the flint panel, and under control conditions. A few candidate genes coincided between the current research and previous reports. A total of 47 flint inbreds harbored the favorable alleles for six significant quantitative trait loci (QTL) detected for inbreds per se evaluated under cold conditions, four of them had also the favorable alleles for the main QTL detected from the testcrosses. Only four dent inbreds (EZ47, F924, NK807 and PHJ40) harbored the favorable alleles for three main QTL detected from the evaluation of the dent inbreds per se under cold conditions. There were more QTL in the flint panel and most of the QTL were associated with days to emergence and ΦPSII.

CONCLUSIONS

These results open new possibilities to genetically improve cold tolerance either with genome-wide selection or with marker assisted selection.

Dent and Flint maize diversity panels reveal important genetic potential for increasing biomass production

Genetic and phenotypic analysis of two complementary maize panels revealed an important variation for biomass yield. Flowering and biomass QTL were discovered by association mapping in both panels. The high whole plant biomass productivity of maize makes it a potential source of energy in animal feeding and biofuel production. The variability and the genetic determinism of traits related to biomass are poorly known. We analyzed two highly diverse panels of Dent and Flint lines representing complementary heterotic groups for Northern Europe. They were genotyped with the 50 k SNP-array and phenotyped as hybrids (crossed to a tester of the complementary pool) in a western European field trial network for traits related to flowering time, plant height, and biomass. The molecular information revealed to be a powerful tool for discovering different levels of structure and relatedness in both panels. This study revealed important variation and potential genetic progress for biomass production, even at constant precocity. Association mapping was run by combining genotypes and phenotypes in a mixed model with a random polygenic effect. This permitted the detection of significant associations, confirming height and flowering time quantitative trait loci (QTL) found in literature. Biomass yield QTL were detected in both panels but were unstable across the environments. Alternative kinship estimator only based on markers unlinked to the tested SNP increased the number of significant associations by around 40% with a satisfying control of the false positive rate. This study gave insights into the variability and the genetic architectures of biomass-related traits in Flint and Dent lines and suggests important potential of these two pools for breeding high biomass yielding hybrid varieties.

Mesenchymal stem cell therapy induces glucocorticoid synthesis in colonic mucosa and suppresses radiation-activated T cells: new insights into MSC immunomodulation

Non-neoplastic tissues around an abdomino-pelvic tumor can be damaged by the radiotherapy protocol, leading to chronic gastrointestinal complications that affect the quality of life with substantial mortality. Stem cell-based approaches using immunosuppressive bone marrow mesenchymal stem cells (MSCs) are promising cell therapy tools. In a rat model of radiation proctitis, we evidenced that a single MSC injection reduces colonic mucosa damages induced by ionizing radiation with improvement of the re-epithelization process for up to 21 days. Immune cell infiltrate and inflammatory molecule expressions in the colonic mucosa were investigated. We report that MSC therapy specifically reduces T-cell infiltration and proliferation, and increases apoptosis of radiation-activated T cells. We assessed the underlying molecular mechanisms and found that interleukin-10 and regulatory T lymphocytes are not involved in the immunosuppressive process in this model. However, an increased level of corticosterone secretion and HSD11b1 (11β-hydroxysteroid dehydrogenase type 1)-steroidogenic enzyme expression was detected in colonic mucosa 21 days after MSC treatment. Moreover, blocking the glucocorticoid (GC) receptor using the RU486 molecule statistically enhances the allogenic lymphocyte proliferation inhibited by MSCs in vitro and abrogates the mucosal protection induced by MSC treatment in vivo. Using the irradiation model, we found evidence for a new MSC immunosuppressive mechanism involving GCs.

Recovering power in association mapping panels with variable levels of linkage disequilibrium

Association mapping has permitted the discovery of major QTL in many species. It can be applied to existing populations and, as a consequence, it is generally necessary to take into account structure and relatedness among individuals in the statistical model to control false positives. We analytically studied power in association studies by computing noncentrality parameter of the tests and its relationship with parameters characterizing diversity (genetic differentiation between groups and allele frequencies) and kinship between individuals. Investigation of three different maize diversity panels genotyped with the 50k SNPs array highlighted contrasted average power among panels and revealed gaps of power of classical mixed models in regions with high linkage disequilibrium (LD). These gaps could be related to the fact that markers are used for both testing association and estimating relatedness. We thus considered two alternative approaches to estimating the kinship matrix to recover power in regions of high LD. In the first one, we estimated the kinship with all the markers that are not located on the same chromosome than the tested SNP. In the second one, correlation between markers was taken into account to weight the contribution of each marker to the kinship. Simulations revealed that these two approaches were efficient to control false positives and were more powerful than classical models.

Out of America: tracing the genetic footprints of the global diffusion of maize

Maize was first domesticated in a restricted valley in south-central Mexico. It was diffused throughout the Americas over thousands of years, and following the discovery of the New World by Columbus, was introduced into Europe. Trade and colonization introduced it further into all parts of the world to which it could adapt. Repeated introductions, local selection and adaptation, a highly diverse gene pool and outcrossing nature, and global trade in maize led to difficulty understanding exactly where the diversity of many of the local maize landraces originated. This is particularly true in Africa and Asia, where historical accounts are scarce or contradictory. Knowledge of post-domestication movements of maize around the world would assist in germplasm conservation and plant breeding efforts. To this end, we used SSR markers to genotype multiple individuals from hundreds of representative landraces from around the world. Applying a multidisciplinary approach combining genetic, linguistic, and historical data, we reconstructed possible patterns of maize diffusion throughout the world from American "contribution" centers, which we propose reflect the origins of maize worldwide. These results shed new light on introductions of maize into Africa and Asia. By providing a first globally comprehensive genetic characterization of landraces using markers appropriate to this evolutionary time frame, we explore the post-domestication evolutionary history of maize and highlight original diversity sources that may be tapped for plant improvement in different regions of the world.

Combined linkage and linkage disequilibrium QTL mapping in multiple families of maize (Zea mays L.) line crosses highlights complementarities between models based on parental haplotype and single locus polymorphism

Advancements in genotyping are rapidly decreasing marker costs and increasing marker density. This opens new possibilities for mapping quantitative trait loci (QTL), in particular by combining linkage disequilibrium information and linkage analysis (LDLA). In this study, we compared different approaches to detect QTL for four traits of agronomical importance in two large multi-parental datasets of maize (Zea mays L.) of 895 and 928 testcross progenies composed of 7 and 21 biparental families, respectively, and genotyped with 491 markers. We compared to traditional linkage-based methods two LDLA models relying on the dense genotyping of parental lines with 17,728 SNP: one based on a clustering approach of parental line segments into ancestral alleles and one based on single marker information. The two LDLA models generally identified more QTL (60 and 52 QTL in total) than classical linkage models (49 and 44 QTL in total). However, they performed inconsistently over datasets and traits suggesting that a compromise must be found between the reduction of allele number for increasing statistical power and the adequacy of the model to potentially complex allelic variation. For some QTL, the model exclusively based on linkage analysis, which assumed that each parental line carried a different QTL allele, was able to capture remaining variation not explained by LDLA models. These complementarities between models clearly suggest that the different QTL mapping approaches must be considered to capture the different levels of allelic variation at QTL involved in complex traits.

Adaptation of maize to temperate climates: mid-density genome-wide association genetics and diversity patterns reveal key genomic regions, with a major contribution of the Vgt2 (ZCN8) locus

The migration of maize from tropical to temperate climates was accompanied by a dramatic evolution in flowering time. To gain insight into the genetic architecture of this adaptive trait, we conducted a 50K SNP-based genome-wide association and diversity investigation on a panel of tropical and temperate American and European representatives. Eighteen genomic regions were associated with flowering time. The number of early alleles cumulated along these regions was highly correlated with flowering time. Polymorphism in the vicinity of the ZCN8 gene, which is the closest maize homologue to Arabidopsis major flowering time (FT) gene, had the strongest effect. This polymorphism is in the vicinity of the causal factor of Vgt2 QTL. Diversity was lower, whereas differentiation and LD were higher for associated loci compared to the rest of the genome, which is consistent with selection acting on flowering time during maize migration. Selection tests also revealed supplementary loci that were highly differentiated among groups and not associated with flowering time in our panel, whereas they were in other linkage-based studies. This suggests that allele fixation led to a lack of statistical power when structure and relatedness were taken into account in a linear mixed model. Complementary designs and analysis methods are necessary to unravel the architecture of complex traits. Based on linkage disequilibrium (LD) estimates corrected for population structure, we concluded that the number of SNPs genotyped should be at least doubled to capture all QTLs contributing to the genetic architecture of polygenic traits in this panel. These results show that maize flowering time is controlled by numerous QTLs of small additive effect and that strong polygenic selection occurred under cool climatic conditions. They should contribute to more efficient genomic predictions of flowering time and facilitate the dissemination of diverse maize genetic resources under a wide range of environments.

OptiMAS: a decision support tool for marker-assisted assembly of diverse alleles

Current advances in plant genotyping lead to major progress in the knowledge of genetic architecture of traits of interest. It is increasingly important to develop decision support tools to help breeders and geneticists to conduct marker-assisted selection methods to assemble favorable alleles that are discovered. Algorithms have been implemented, within an interactive graphical interface, to 1) trace parental alleles throughout generations, 2) propose strategies to select the best plants based on estimated molecular scores, and 3) efficiently intermate them depending on the expected value of their progenies. With the possibility to consider a multi-allelic context, OptiMAS opens new prospects to assemble favorable alleles issued from diverse parents and further accelerate genetic gain.

Intraspecific variation of recombination rate in maize

BACKGROUND

In sexually reproducing organisms, meiotic crossovers ensure the proper segregation of chromosomes and contribute to genetic diversity by shuffling allelic combinations. Such genetic reassortment is exploited in breeding to combine favorable alleles, and in genetic research to identify genetic factors underlying traits of interest via linkage or association-based approaches. Crossover numbers and distributions along chromosomes vary between species, but little is known about their intraspecies variation.

RESULTS

Here, we report on the variation of recombination rates between 22 European maize inbred lines that belong to the Dent and Flint gene pools. We genotype 23 doubled-haploid populations derived from crosses between these lines with a 50 k-SNP array and construct high-density genetic maps, showing good correspondence with the maize B73 genome sequence assembly. By aligning each genetic map to the B73 sequence, we obtain the recombination rates along chromosomes specific to each population. We identify significant differences in recombination rates at the genome-wide, chromosome, and intrachromosomal levels between populations, as well as significant variation for genome-wide recombination rates among maize lines. Crossover interference analysis using a two-pathway modeling framework reveals a negative association between re combination rate and interference strength.

CONCLUSIONS

To our knowledge, the present work provides the most comprehensive study on intraspecific variation of recombination rates and crossover interference strength in eukaryotes. Differences found in recombination rates will allow for selection of high or low recombining lines in crossing programs. Our methodology should pave the way for precise identification of genes controlling recombination rates in maize and other organisms.

A composite map of expressed sequences in maize

A maize genetic map based mainly on expressed sequences has been constructed. The map incorporates data from four segregating populations. Three recombinant inbred line populations were derived from the nonreciprocal crosses between three inbred lines. A map derived from an independent F2 progeny from one of the crosses was also used. With a total of 521 genotyped individuals, accuracy in gene order is expected. Five sources of markers were used: (i) 109 loci corresponding to 69 genes of known function, (ii) 39 loci controlling protein position shifts revealed by two-dimensional electrophoresis, (iii) 8 isozyme loci, (iv) 17 loci corresponding to 14 sequenced cDNAs for which no homology was found in gene banks, and (v) 102 loci corresponding to 81 anonymous probes. As many loci were common to all maps, we tested heterogeneity between recombination fractions. The comparison of recombination fractions revealed: (i) a good correspondence between the maps derived from the same cross, (ii) few significant differences in interval distances, and (iii) global differences, which can reach 20% of the total map length. A composite map of 275 loci covering 1765 cM has been constructed. Key words : Zea mays L., RFLP, genetic map, molecular markers, proteins.

Diversity and linkage disequilibrium features in a composite public/private dent maize panel: consequences for association genetics as evaluated from a case study using flowering time

Recent progress in genotyping and resequencing techniques have opened new opportunities for deciphering quantitative trait variation by looking for associations between traits of interest and polymorphisms in panels of diverse inbred lines. Association mapping raises specific issues related to the choice of appropriate (i) panels and marker-densities and (ii) statistical methods to capture associations. In this study, we used a panel of 314 maize inbred lines from the dent pool, composed of inbred material from public institutes (113 inbred lines) and a private company (201 inbred lines). We showed that local LD was higher and genetic diversity lower in the material of private origin than in the public material. We compared the results obtained by different software for identifying population structure and computing relatedness among lines, and ran association tests for earliness related traits. Our results confirmed the importance of the mite polymorphism of Vgt1 on flowering time, but also showed that its effect can be captured by zmRap2.7 polymorphisms located 70 kb apart. We also highlighted associations with polymorphisms within genes putatively involved in lignin biosynthesis pathway, which deserve further investigations.

BioMercator V3: an upgrade of genetic map compilation and quantitative trait loci meta-analysis algorithms

Summary: Compilation of genetic maps combined to quantitative trait loci (QTL) meta-analysis has proven to be a powerful approach contributing to the identification of candidate genes underlying quantitative traits. BioMercator was the first software offering a complete set of algorithms and visualization tool covering all steps required to perform QTL meta-analysis. Despite several limitations, the software is still widely used. We developed a new version proposing additional up to date methods and improving graphical representation and exploration of large datasets.

Availability and implementation: BioMercator V3 is implemented in JAVA and freely available (http://moulon.inra.fr/biomercator)

Contact:joets@moulon.inra.fr

Maize genetic diversity and association mapping using transposable element insertion polymorphisms

Transposable elements are the major component of the maize genome and presumably highly polymorphic yet they have not been used in population genetics and association analyses. Using the Transposon Display method, we isolated and converted into PCR-based markers 33 Miniature Inverted Repeat Transposable Elements (MITE) polymorphic insertions. These polymorphisms were genotyped on a population-based sample of 26 American landraces for a total of 322 plants. Genetic diversity was high and partitioned within and among landraces. The genetic groups identified using Bayesian clustering were in agreement with published data based on SNPs and SSRs, indicating that MITE polymorphisms reflect maize genetic history. To explore the contribution of MITEs to phenotypic variation, we undertook an association mapping approach in a panel of 367 maize lines phenotyped for 26 traits. We found a highly significant association between the marker ZmV1-9, on chromosome 1, and male flowering time. The variance explained by this association is consistent with a flowering delay of +123 degree-days. This MITE insertion is located at only 289 nucleotides from the 3' end of a Cytochrome P450-like gene, a region that was never identified in previous association mapping or QTL surveys. Interestingly, we found (i) a non-synonymous mutation located in the exon 2 of the gene in strong linkage disequilibrium with the MITE polymorphism, and (ii) a perfect sequence homology between the MITE sequence and a maize siRNA that could therefore potentially interfere with the expression of the Cytochrome P450-like gene. Those two observations among others offer exciting perspectives to validate functionally the role of this region on phenotypic variation.

A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome

SNP genotyping arrays have been useful for many applications that require a large number of molecular markers such as high-density genetic mapping, genome-wide association studies (GWAS), and genomic selection. We report the establishment of a large maize SNP array and its use for diversity analysis and high density linkage mapping. The markers, taken from more than 800,000 SNPs, were selected to be preferentially located in genes and evenly distributed across the genome. The array was tested with a set of maize germplasm including North American and European inbred lines, parent/F1 combinations, and distantly related teosinte material. A total of 49,585 markers, including 33,417 within 17,520 different genes and 16,168 outside genes, were of good quality for genotyping, with an average failure rate of 4% and rates up to 8% in specific germplasm. To demonstrate this array's use in genetic mapping and for the independent validation of the B73 sequence assembly, two intermated maize recombinant inbred line populations - IBM (B73×Mo17) and LHRF (F2×F252) - were genotyped to establish two high density linkage maps with 20,913 and 14,524 markers respectively. 172 mapped markers were absent in the current B73 assembly and their placement can be used for future improvements of the B73 reference sequence. Colinearity of the genetic and physical maps was mostly conserved with some exceptions that suggest errors in the B73 assembly. Five major regions containing non-colinearities were identified on chromosomes 2, 3, 6, 7 and 9, and are supported by both independent genetic maps. Four additional non-colinear regions were found on the LHRF map only; they may be due to a lower density of IBM markers in those regions or to true structural rearrangements between lines. Given the array's high quality, it will be a valuable resource for maize genetics and many aspects of maize breeding.

A common genetic determinism for sensitivities to soil water deficit and evaporative demand: meta-analysis of quantitative trait Loci and introgression lines of maize

Evaporative demand and soil water deficit equally contribute to water stress and to its effect on plant growth. We have compared the genetic architectures of the sensitivities of maize (Zea mays) leaf elongation rate with evaporative demand and soil water deficit. The former was measured via the response to leaf-to-air vapor pressure deficit in well-watered plants, the latter via the response to soil water potential in the absence of evaporative demand. Genetic analyses of each sensitivity were performed over 21 independent experiments with (1) three mapping populations, with temperate or tropical materials, (2) one population resulting from the introgression of a tropical drought-tolerant line in a temperate line, and (3) two introgression libraries genetically independent from mapping populations. A very large genetic variability was observed for both sensitivities. Some lines maintained leaf elongation at very high evaporative demand or water deficit, while others stopped elongation in mild conditions. A complex architecture arose from analyses of mapping populations, with 19 major meta-quantitative trait loci involving strong effects and/or more than one mapping population. A total of 68% of those quantitative trait loci affected sensitivities to both evaporative demand and soil water deficit. In introgressed lines, 73% of the tested genomic regions affected both sensitivities. To our knowledge, this study is the first genetic demonstration that hydraulic processes, which drive the response to evaporative demand, also have a large contribution to the genetic variability of plant growth under water deficit in a large range of genetic material.