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.

Physiological adaptive traits are a potential allele reservoir for maize genetic progress under challenging conditions.

Combined phenomic and genomic approaches are required to evaluate the margin of progress of breeding strategies. Here, we analyze 65 years of genetic progress in maize yield, which was similar (101 kg ha-1 year-1) across most frequent environmental scenarios in the European growing area. Yield gains were linked to physiologically simple traits (plant phenology and architecture) which indirectly affected reproductive development and light interception in all studied environments, marked by significant genomic signatures of selection. Conversely, studied physiological processes involved in stress adaptation remained phenotypically unchanged (e.g. stomatal conductance and growth sensitivity to drought) and showed no signatures of selection. By selecting for yield, breeders indirectly selected traits with stable effects on yield, but not physiological traits whose effects on yield can be positive or negative depending on environmental conditions. Because yield stability under climate change is desirable, novel breeding strategies may be needed for exploiting alleles governing physiological adaptive traits.

The genetic basis of cytoplasmic male sterility and fertility restoration in wheat.

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

Deciphering the genetics of flowering time by an association study on candidate genes in bread wheat (Triticum aestivum L.).

Earliness is very important for the adaptation of wheat to environmental conditions and the achievement of high grain yield. A detailed knowledge of key genetic components of the life cycle would enable an easier control by the breeders. The objective of the study was to investigate the effect of candidate genes on flowering time. Using a collection of hexaploid wheat composed of 235 lines from diverse geographical origins, we conducted an association study for six candidate genes for flowering time and its components (vernalization sensitivity and earliness per se). The effect on the variation of earliness components of polymorphisms within the copies of each gene was tested in ANOVA models accounting for the underlying genetic structure. The collection was structured in five groups that minimized the residual covariance. Vernalization requirement and lateness tend to increase according to the mean latitude of each group. Heading date for an autumnal sowing was mainly determined by the earliness per se. Except for the Constans (CO) gene orthologous of the barley HvCO3, all gene polymorphisms had a significant impact on earliness components. The three traits used to quantify vernalization requirement were primarily associated with polymorphisms at Vrn-1 and then at Vrn-3 and Luminidependens (LD) genes. We found a good correspondence between spring/winter types and genotypes at the three homeologous copies of Vrn-1. Earliness per se was mainly explained by polymorphisms at Vrn-3 and to a lesser extent at Vrn-1, Hd-1 and Gigantea (GI) genes. Vernalization requirement and earliness as a function of geographical origin, as well as the possible role of the breeding practices in the geographical distribution of the alleles and the hypothetical adaptive value of the candidate genes, are discussed.

Cross-genome map based dissection of a nitrogen use efficiency ortho-metaQTL in bread wheat unravels concerted cereal genome evolution.

Monitoring nitrogen use efficiency (NUE) in plants is becoming essential to maintain yield while reducing fertilizer usage. Optimized NUE application in major crops is essential for long-term sustainability of agriculture production. Here, we report the precise identification of 11 major chromosomal regions controlling NUE in wheat that co-localise with key developmental genes such as Ppd (photoperiod sensitivity), Vrn (vernalization requirement), Rht (reduced height) and can be considered as robust markers from a molecular breeding perspective. Physical mapping, sequencing, annotation and candidate gene validation of an NUE metaQTL on wheat chromosome 3B allowed us to propose that a glutamate synthase (GoGAT) gene that is conserved structurally and functionally at orthologous positions in rice, sorghum and maize genomes may contribute to NUE in wheat and other cereals. We propose an evolutionary model for the NUE locus in cereals from a common ancestral region, involving species specific shuffling events such as gene deletion, inversion, transposition and the invasion of repetitive elements.

Genomics in cereals: from genome-wide conserved orthologous set (COS) sequences to candidate genes for trait dissection.

Recent updates in comparative genomics among cereals have provided the opportunity to identify conserved orthologous set (COS) DNA sequences for cross-genome map-based cloning of candidate genes underpinning quantitative traits. New tools are described that are applicable to any cereal genome of interest, namely, alignment criterion for orthologous couples identification, as well as the Intron Spanning Marker software to automatically select intron-spanning primer pairs. In order to test the software, it was applied to the bread wheat genome, and 695 COS markers were assigned to 1,535 wheat loci (on average one marker/2.6 cM) based on 827 robust rice-wheat orthologs. Furthermore, 31 of the 695 COS markers were selected to fine map a pentosan viscosity quantitative trait loci (QTL) on wheat chromosome 7A. Among the 31 COS markers, 14 (45%) were polymorphic between the parental lines and 12 were mapped within the QTL confidence interval with one marker every 0.6 cM defining candidate genes among the rice orthologous region.

Epistatic interactions between Opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize.

Association genetics is a powerful method to track gene polymorphisms responsible for phenotypic variation, since it takes advantage of existing collections and historical recombination to study the correlation between large genetic diversity and phenotypic variation. We used a collection of 375 maize (Zea mays ssp. mays) inbred lines representative of tropical, American, and European diversity, previously characterized for genome-wide neutral markers and population structure, to investigate the roles of two functionally related candidate genes, Opaque2 and CyPPDK1, on kernel quality traits. Opaque2 encodes a basic leucine zipper transcriptional activator specifically expressed during endosperm development that controls the transcription of many target genes, including CyPPDK1, which encodes a cytosolic pyruvate orthophosphate dikinase. Using statistical models that correct for population structure and individual kinship, Opaque2 polymorphism was found to be strongly associated with variation of the essential amino acid lysine. This effect could be due to the direct role of Opaque2 on either zein transcription, zeins being major storage proteins devoid of lysine, or lysine degradation through the activation of lysine ketoglutarate reductase. Moreover, we found that a polymorphism in the Opaque2 coding sequence and several polymorphisms in the CyPPDK1 promoter nonadditively interact to modify both lysine content and the protein-versus-starch balance, thus revealing the role in quantitative variation in plants of epistatic interactions between a transcriptional activator and one of its target genes.

Structure and expression analysis of rice paleo duplications.

Having a well-known history of genome duplication, rice is a good model for studying structural and functional evolution of paleo duplications. Improved sequence alignment criteria were used to characterize 10 major chromosome-to-chromosome duplication relationships associated with 1440 paralogous pairs, covering 47.8% of the rice genome, with 12.6% of genes that are conserved within sister blocks. Using a micro-array experiment, a genome-wide expression map has been produced, in which 2382 genes show significant differences of expression in root, leaf and grain. By integrating both structural (1440 paralogous pairs) and functional information (2382 differentially expressed genes), we identified 115 paralogous gene pairs for which at least one copy is differentially expressed in one of the three tissues. A vast majority of the 115 paralogous gene pairs have been neofunctionalized or subfunctionalized as 88%, 89% and 96% of duplicates, respectively, expressed in grain, leaf and root show distinct expression patterns. On the basis of a Gene Ontology analysis, we have identified and characterized the gene families that have been structurally and functionally preferentially retained in the duplication showing that the vast majority (>85%) of duplicated have been either lost or have been subfunctionalized or neofunctionalized during 50-70 million years of evolution.

A joint transcriptomic, proteomic and metabolic analysis of maize endosperm development and starch filling.

The maize endosperm transcriptome was investigated through cDNA libraries developed at three characteristic stages: (i) lag phase [10 days after pollination (DAP)]; (ii) beginning of storage (14 DAP); and (iii) maximum starch accumulation rate (21 DAP). Expressed sequence tags for 711, 757 and 384 relevant clones, respectively, were obtained and checked manually. The proportion of sequences with no clear function decreased from 35% to 20%, and a large increase in storage protein sequences (i.e. 5% to 38%) was observed from stages (i) to (iii). The remaining major categories included metabolism (11%-13%), transcription-RNA processing-protein synthesis (13%-20%), protein destination (5%-9%), cellular communication (3%-9%) and cell rescue-defence (4%). Good agreement was generally found between category rank in the 10-DAP transcriptome and the recently reported 14-DAP proteome, except that kinases and proteins for RNA processing were not detected in the latter. In the metabolism category, the respiratory pathway transcripts represented the largest proportion (25%-37%), and showed a shift in favour of glycolysis at 21 DAP. At this stage, amino acid metabolism increased to 17%, whereas starch metabolism surprisingly decreased to 7%. A second experiment focused on carbohydrate metabolism by comparing gene expression at three levels (transcripts, proteins and enzyme activities) in relation to substrate or product from 10 to 40 DAP. Here, two distinct patterns were observed: invertases and hexoses were predominant at the beginning, whereas enzyme patterns in the starch pathway, at the three levels, anticipated and paralleled starch accumulation, suggesting that, in most cases, transcriptional control is responsible for the regulation of starch biosynthesis.

MAIZEWALL. Database and developmental gene expression profiling of cell wall biosynthesis and assembly in maize.

An extensive search for maize (Zea mays) genes involved in cell wall biosynthesis and assembly has been performed and 735 sequences have been centralized in a database, MAIZEWALL (http://www.polebio.scsv.ups-tlse.fr/MAIZEWALL). MAIZEWALL contains a bioinformatic analysis for each entry and gene expression data that are accessible via a user-friendly interface. A maize cell wall macroarray composed of a gene-specific tag for each entry was also constructed to monitor global cell wall-related gene expression in different organs and during internode development. By using this macroarray, we identified sets of genes that exhibit organ and internode-stage preferential expression profiles. These data provide a comprehensive fingerprint of cell wall-related gene expression throughout the maize plant. Moreover, an in-depth examination of genes involved in lignin biosynthesis coupled to biochemical and cytological data from different organs and stages of internode development has also been undertaken. These results allow us to trace spatially and developmentally regulated, putative preferential routes of monolignol biosynthesis involving specific gene family members and suggest that, although all of the gene families of the currently accepted monolignol biosynthetic pathway are conserved in maize, there are subtle differences in family size and a high degree of complexity in spatial expression patterns. These differences are in keeping with the diversity of lignified cell types throughout the maize plant.

Two cytosolic glutamine synthetase isoforms of maize are specifically involved in the control of grain production.

The roles of two cytosolic maize glutamine synthetase isoenzymes (GS1), products of the Gln1-3 and Gln1-4 genes, were investigated by examining the impact of knockout mutations on kernel yield. In the gln1-3 and gln1-4 single mutants and the gln1-3 gln1-4 double mutant, GS mRNA expression was impaired, resulting in reduced GS1 protein and activity. The gln1-4 phenotype displayed reduced kernel size and gln1-3 reduced kernel number, with both phenotypes displayed in gln1-3 gln1-4. However, at maturity, shoot biomass production was not modified in either the single mutants or double mutants, suggesting a specific impact on grain production in both mutants. Asn increased in the leaves of the mutants during grain filling, indicating that it probably accumulates to circumvent ammonium buildup resulting from lower GS1 activity. Phloem sap analysis revealed that unlike Gln, Asn is not efficiently transported to developing kernels, apparently causing reduced kernel production. When Gln1-3 was overexpressed constitutively in leaves, kernel number increased by 30%, providing further evidence that GS1-3 plays a major role in kernel yield. Cytoimmunochemistry and in situ hybridization revealed that GS1-3 is present in mesophyll cells, whereas GS1-4 is specifically localized in the bundle sheath cells. The two GS1 isoenzymes play nonredundant roles with respect to their tissue-specific localization.

Single-nucleotide polymorphism frequency in a set of selected lines of bread wheat (Triticum aestivum L.).

Information on single-nucleotide polymorphisms (SNPs) in hexaploid bread wheat is still scarce. The goal of this study was to detect SNPs in wheat and examine their frequency. Twenty-six bread wheat lines from different origins worldwide were used. Specific PCR-products were obtained from 21 genes and directly sequenced. SNPs were discovered from the alignment of these sequences. The overall sequence polymorphism observed in this sample appears to be low; 64 single-base polymorphisms were detected in approximately 21.5 kb (i.e., 1 SNP every 335 bp). The level of polymorphism is highly variable among the different genes studied. Fifty percent of the genes studied contained no sequence polymorphism, whereas most SNPs detected were located in only 2 genes. As expected, taking into account a synthetic line created with a wild Triticum tauschii parent increases the level of polymorphism (101 SNPs; 1 SNP every 212 bp). The detected SNPs are available at http://urgi.versailles.inra.fr/GnpSNP">http://urgi.versailles.inra.fr/GnpSNP. Data on linkage disequilibrium (LD) are still preliminary. They showed a significant level of LD in the 2 most polymorphic genes. To conclude, the genome size of hexaploid wheat and its low level of polymorphism complicate SNP discovery in this species.

Identification of Glu-B1-1 as a candidate gene for the quantity of high-molecular-weight glutenin in bread wheat (Triticum aestivum L.) by means of an association study.

A previous study in wheat (Triticum aestivum L.) identified two candidate genes controlling a quantitative trait locus (QTL) for high-molecular-weight glutenin subunit (HMW-GS) GluBx. These candidates were Glu-B1-1, the structural gene coding for Glu1Bx, and the B homoeologous gene coding for SPA (spa-B), a seed storage protein activator. The goal of this study was to identify the best candidate gene for this QTL. Single nucleotide polymorphisms (SNPs) are an abundant source of DNA polymorphisms that have been successfully used to identify loci associated with particular phenotypes. As no linkage disequilibrium was detected between Glu-B1-1 and spa-B, we performed an association study to identify the individual gene responsible for the QTL. Six SNPs, three located in Glu-B1-1 and three in spa-B, were genotyped by mass spectrometry in a collection of 113 bread wheat lines. These lines were also evaluated for protein content as well as the total quantity of HMW-GSs and of each HMW-GS in seed samples from two harvest years. Significant associations were detected only between Glu-B1-1 polymorphism and most of the traits evaluated. Spa-B was unambiguously discarded as a candidate. To our knowledge, this is the first report on an association study that was successfully used to discriminate between two candidate genes.

UDP-glucose dehydrogenases of maize: a role in cell wall pentose biosynthesis.

UDPGDH (UDP-D-glucose dehydrogenase) oxidizes UDP-Glc (UDP-D-glucose) to UDP-GlcA (UDP-D-glucuronate), the precursor of UDP-D-xylose and UDP-L-arabinose, major cell wall polysaccharide precursors. Maize (Zea mays L.) has at least two putative UDPGDH genes (A and B), according to sequence similarity to a soya bean UDPGDH gene. The predicted maize amino acid sequences have 95% similarity to that of soya bean. Maize mutants with a Mu-element insertion in UDPGDH-A or UDPGDH-B were isolated (udpgdh-A1 and udpgdh-B1 respectively) and studied for changes in wall polysaccharide biosynthesis. The udpgdh-A1 and udpgdh-B1 homozygotes showed no visible phenotype but exhibited 90 and 60-70% less UDPGDH activity respectively than wild-types in a radiochemical assay with 30 microM UDP-glucose. Ethanol dehydrogenase (ADH) activity varied independently of UDPGDH activity, supporting the hypothesis that ADH and UDPGDH activities are due to different enzymes in maize. When extracts from wild-types and udpgdh-A1 homozygotes were assayed with increasing concentrations of UDP-Glc, at least two isoforms of UDPGDH were detected, having K(m) values of approx. 380 and 950 microM for UDP-Glc. Leaf and stem non-cellulosic polysaccharides had lower Ara/Gal and Xyl/Gal ratios in udpgdh-A1 homozygotes than in wild-types, whereas udpgdh-B1 homozygotes exhibited more variability among individual plants, suggesting that UDPGDH-A activity has a more important role than UDPGDH-B in UDP-GlcA synthesis. The fact that mutation of a UDPGDH gene interferes with polysaccharide synthesis suggests a greater importance for the sugar nucleotide oxidation pathway than for the myo-inositol pathway in UDP-GlcA biosynthesis during post-germinative growth of maize.

Linkage mapping of 1454 new maize candidate gene Loci.

Bioinformatic analyses of maize EST sequences have highlighted large numbers of candidate genes putatively involved in agriculturally important traits. To contribute to ongoing efforts toward mapping of these genes, we used two populations of intermated recombinant inbred lines (IRILs), which allow a higher map resolution than nonintermated RILs. The first panel (IBM), derived from B73 x Mo17, is publicly available from the Maize Genetics Cooperation Stock Center. The second panel (LHRF) was developed from F2 x F252 to map loci monomorphic on IBM. We built framework maps of 237 loci from the IBM panel and 271 loci from the LHRF panel. Both maps were used to place 1454 loci (1056 on map IBM_Gnp2004 and 398 on map LHRF_Gnp2004) that corresponded to 954 cDNA probes previously unmapped. RFLP was mostly used, but PCR-based methods were also performed for some cDNAs to map SNPs. Unlike in usual IRIL-based maps published so far, corrected meiotic centimorgan distances were calculated, taking into account the number of intermating generations undergone by the IRILs. The corrected sizes of our framework maps were 1825 cM for IBM_Gnp2004 and 1862 cM for LHRF_Gnp2004. All loci mapped on LHRF_Gnp2004 were also projected on a consensus map IBMconsensus_Gnp2004. cDNA loci formed clusters near the centromeres except for chromosomes 1 and 8.

Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome.

Genetic architecture of flowering time in maize was addressed by synthesizing a total of 313 quantitative trait loci (QTL) available for this trait. These were analyzed first with an overview statistic that highlighted regions of key importance and then with a meta-analysis method that yielded a synthetic genetic model with 62 consensus QTL. Six of these displayed a major effect. Meta-analysis led in this case to a twofold increase in the precision in QTL position estimation, when compared to the most precise initial QTL position within the corresponding region. The 62 consensus QTL were compared first to the positions of the few flowering-time candidate genes that have been mapped in maize. We then projected rice candidate genes onto the maize genome using a synteny conservation approach based on comparative mapping between the maize genetic map and japonica rice physical map. This yielded 19 associations between maize QTL and genes involved in flowering time in rice and in Arabidopsis. Results suggest that the combination of meta-analysis within a species of interest and synteny-based projections from a related model plant can be an efficient strategy for identifying new candidate genes for trait variation.

Nucleotide diversity of the ZmPox3 maize peroxidase gene: relationships between a MITE insertion in exon 2 and variation in forage maize digestibility.

BACKGROUND: Polymorphisms were investigated within the ZmPox3 maize peroxidase gene, possibly involved in lignin biosynthesis because of its colocalization with a cluster of QTL related to lignin content and cell wall digestibility. The purpose of this study was to identify, on the basis of 37 maize lines chosen for their varying degrees of cell wall digestibility and representative of temperate regions germplasm, ZmPox3 haplotypes or individual polymorphisms possibly associated with digestibility. RESULTS: Numerous haplotypes with high diversity were identified. Frequency of nucleotide changes was high with on average one SNP every 57 bp. Nucleotide diversity was not equally distributed among site categories: the estimated pi was on average eight times higher for silent sites than for non-synonymous sites. Numerous sites were in linkage disequilibrium that decayed with increasing physical distance. A zmPox3 mutant allele, carrying an insertion of a transposable element in the second exon, was found in lines derived from the early flint inbred line, F7. This element possesses many structural features of miniature inverted-repeat transposable elements (MITE). The mutant allele encodes a truncated protein lacking important functional sites. An ANOVA performed with a subset of 31 maize lines indicated that the transposable element was significantly associated with cell wall digestibility. This association was confirmed using an additional set of 25 flint lines related to F7. Moreover, RT-PCR experiments revealed a decreased amount of corresponding mRNA in plants with the MITE insertion. CONCLUSION: These results showed that ZmPox3 could possibly be involved in monolignol polymerisation, and that a deficiency in ZmPox3 peroxidase activity seemingly has a negative effect on cell wall digestibility. Also, genetic diversity analyses of ZmPox3 indicated that this peroxidase could be a relevant target for grass digestibility improvement using specific allele introgressions.

Microsatellite-based deletion bin system for the establishment of genetic-physical map relationships in wheat (Triticum aestivum L.).

Because of polyploidy and large genome size, deletion stocks of bread wheat are an ideal material for physically allocating ESTs and genes to small chromosomal regions for targeted mapping. To enhance the utility of deletion stocks for chromosome bin mapping, we characterized a set of 84 deletion lines covering the 21 chromosomes of wheat using 725 microsatellites. We localized these microsatellite loci to 94 breakpoints in a homozygous state (88 distal deletions, 6 interstitial), and 5 in a heterozygous state representing 159 deletion bins. Chromosomes from homoeologous groups 2 and 5 were the best covered (126 and 125 microsatellites, respectively) while the coverage for group 4 was lower (80 microsatellites). We assigned at least one microsatellite in up to 92% of the bins (mean 4.97 SSR/bin). Only a few discrepancies concerning marker order were observed. The cytogenetic maps revealed small genetic distances over large physical regions around the centromeres and large genetic to physical map ratios close to the telomeres. As SSRs are the markers of choice for many genetic and breeding studies, the mapped microsatellite loci will be useful not only for deletion stock verifications but also for allocating associated QTLs to deletion bins where numerous ESTs that could be potential candidate genes are currently assigned.