Plant molecular diversity and applications to genomics

Association Analysis in Rice: From Application to Utilization

Association analysis based on linkage disequilibrium (LD) is an efficient way to dissect complex traits and to identify gene functions in rice. Although association analysis is an effective way to construct fine maps for quantitative traits, there are a few issues which need to be addressed. In this review, we will first summarize type, structure, and LD level of populations used for association analysis of rice, and then discuss the genotyping methods and statistical approaches used for association analysis in rice. Moreover, we will review current shortcomings and benefits of association analysis as well as specific types of future research to overcome these shortcomings. Furthermore, we will analyze the reasons for the underutilization of the results within association analysis in rice breeding.

Gene Discovery and Advances in Finger Millet [Eleusine coracana (L.) Gaertn.] Genomics-An Important Nutri-Cereal of Future

The rapid strides in molecular marker technologies followed by genomics, and next generation sequencing advancements in three major crops (rice, maize and wheat) of the world have given opportunities for their use in the orphan, but highly valuable future crops, including finger millet [Eleusine coracana (L.) Gaertn.]. Finger millet has many special agronomic and nutritional characteristics, which make it an indispensable crop in arid, semi-arid, hilly and tribal areas of India and Africa. The crop has proven its adaptability in harsh conditions and has shown resilience to climate change. The adaptability traits of finger millet have shown the advantage over major cereal grains under stress conditions, revealing it as a storehouse of important genomic resources for crop improvement. Although new technologies for genomic studies are now available, progress in identifying and tapping these important alleles or genes is lacking. RAPDs were the default choice for genetic diversity studies in the crop until the last decade, but the subsequent development of SSRs and comparative genomics paved the way for the marker assisted selection in finger millet. Resistance gene homologs from NBS-LRR region of finger millet for blast and sequence variants for nutritional traits from other cereals have been developed and used invariably. Population structure analysis studies exhibit 2-4 sub-populations in the finger millet gene pool with separate grouping of Indian and exotic genotypes. Recently, the omics technologies have been efficiently applied to understand the nutritional variation, drought tolerance and gene mining. Progress has also occurred with respect to transgenics development. This review presents the current biotechnological advancements along with research gaps and future perspective of genomic research in finger millet.

Association mapping of North American spring wheat breeding germplasm reveals loci conferring resistance to Ug99 and other African stem rust races

BACKGROUND

The recently identified Puccinia graminis f. sp. tritici (Pgt) race TTKSK (Ug99) poses a severe threat to global wheat production because of its broad virulence on several widely deployed resistance genes. Additional virulences have been detected in the Ug99 group of races, and the spread of this race group has been documented across wheat growing regions in Africa, the Middle East (Yemen), and West Asia (Iran). Other broadly virulent Pgt races, such as TRTTF and TKTTF, present further difficulties in maintaining abundant genetic resistance for their effective use in wheat breeding against this destructive fungal disease of wheat. In an effort to identify loci conferring resistance to these races, a genome-wide association study was carried out on a panel of 250 spring wheat breeding lines from the International Maize and Wheat Improvement Center (CIMMYT), six wheat breeding programs in the United States and three wheat breeding programs in Canada.

RESULTS

The lines included in this study were grouped into two major clusters, based on the results of principal component analysis using 23,976 SNP markers. Upon screening for adult plant resistance (APR) to Ug99 during 2013 and 2014 in artificial stem rust screening nurseries at Njoro, Kenya and at Debre Zeit, Ethiopia, several wheat lines were found to exhibit APR. The lines were also screened for resistance at the seedling stage against races TTKSK, TRTTF, and TKTTF at USDA-ARS Cereal Disease Laboratory in St. Paul, Minnesota; and only 9 of the 250 lines displayed seedling resistance to all the races. Using a mixed linear model, 27 SNP markers associated with APR against Ug99 were detected, including markers linked with the known APR gene Sr2. Using the same model, 23, 86, and 111 SNP markers associated with seedling resistance against races TTKSK, TRTTF, and TKTTF were identified, respectively. These included markers linked to the genes Sr8a and Sr11 providing seedling resistance to races TRTTF and TKTTF, respectively. We also identified putatively novel Sr resistance genes on chromosomes 3B, 4D, 5A, 5B, 6A, 7A, and 7B.

CONCLUSION

Our results demonstrate that the North American wheat breeding lines have several resistance loci that provide APR and seedling resistance to highly virulent Pgt races. Using the resistant lines and the SNP markers identified in this study, marker-assisted resistance breeding can assist in development of varieties with elevated levels of resistance to virulent stem rust races including TTKSK.

Ultra-High Density, Transcript-Based Genetic Maps of Pepper Define Recombination in the Genome and Synteny Among Related Species

Our ability to assemble complex genomes and construct ultradense genetic maps now allows the determination of recombination rates, translocations, and the extent of genomic collinearity between populations, species, and genera. We developed two ultradense genetic linkage maps for pepper from single-position polymorphisms (SPPs) identified de novo with a 30,173 unigene pepper genotyping array. The Capsicum frutescens × C. annuum interspecific and the C. annuum intraspecific genetic maps were constructed comprising 16,167 and 3,878 unigene markers in 2108 and 783 genetic bins, respectively. Accuracies of marker groupings and orders are validated by the high degree of collinearity between the two maps. Marker density was sufficient to locate the chromosomal breakpoint resulting in the P1/P8 translocation between C. frutescens and C. annuum to a single bin. The two maps aligned to the pepper genome showed varying marker density along the chromosomes. There were extensive chromosomal regions with suppressed recombination and reduced intraspecific marker density. These regions corresponded to the pronounced nonrecombining pericentromeric regions in tomato, a related Solanaceous species. Similar to tomato, the extent of reduced recombination appears to be more pronounced in pepper than in other plant species. Alignment of maps with the tomato and potato genomes shows the presence of previously known translocations and a translocation event that was not observed in previous genetic maps of pepper.

Association Analysis of Simple Sequence Repeat (SSR) Markers with Agronomic Traits in Tall Fescue (Festuca arundinacea Schreb.)

Tall fescue is widely used in temperate regions throughout the world as a dominant forage grass as well as a turfgrass, in pastoral and turf industry. However, the utilization of tall fescue was limited because of its leaf roughness, poor regeneration ability and poor stress resistance. New cultivars were desirable in modern pastoral industries exceed the potential of existing cultivars. Therefore, well understanding the agronomic traits and describing germplasms would help to overcome these constraints, and morphological evaluation of tall fescue germplasm is the key component in selecting rational parents for hybridization breeding. However, describing the morphological traits of tall fescue germplasm is costly and time-consuming. Fortunately, biotechnology approaches can supplement conventional breeding efforts for tall fescue improvement. Association mapping, as a powerful approach to identify association between agronomic traits and molecular markers has been widely used for enhancing the utilization, conservation and management of the tall fescue germplasms. Therefore, in the present research, 115 tall fescue accessions from different origins (25 accessions are cultivars; 31 accessions from America; 32 accessions from European; 7 accessions from Africa; 20 accessions from Asia), were evaluated for agronomic traits and genetic diversity with 90 simple sequence repeat (SSR) markers. The panel displayed significant variation in spike count per plant (SCP) and spike weight (SW). However, BCS performed the lowest CV among all the observed agronomic traits. Three subpopulations were identified within the collections but no obvious relative kinship (K) was found. The GLM model was used to describe the association between SSR and agronomic traits. Fifty-one SSR markers associated with agronomic traits were observed. Twelve single-associated markers were associated with PH; six single-associated markers were associated with BCS; eight single-associated markers were associated with SW; five single-associated markers were associated with SC; seven single-associated markers were associated with SCP; three single-associated markers were associated with SL. Especially, we observed that the genetic variation of SW was explained 11.6 % by M37 marker. It is interesting to observe that nine markers (M1, M2, M35, M54 marker was associated with both BCS and SC; M3, M4 markers were associated with BCS, SW, and SC; M19 marker was associated with both pH and PD, M40 marker was associated with both SCP and SW; and M193 marker was associated with both PH and SL) were associated with more than two agronomic traits. Notably, Branch count per spike (BCS) was explained by four markers (M1, M2, M3, and M4) exceeding 10 %. These identified marker alleles associated with agronomic traits could provide important information and markers for molecular-assisted breeding that facilitate the breeding process in tall fescue.

Nucleotide diversity analysis of three major bacterial blight resistance genes in rice

Nucleotide sequence polymorphisms among R gene alleles influence the process of co-evolutionary interaction between host and pathogen by shaping the response of host plants towards invading pathogens. Here, we present the DNA sequence polymorphisms and diversities present among natural alleles of three rice bacterial blight resistance genes, Xa21, Xa26 and xa5. The diversity was examined across different wild relatives and cultivars of Oryza species. Functional significance of selected alleles was evaluated through semi-quantitative reverse transcription polymerase chain reaction and real time PCR. The greatest nucleotide diversity and singleton variable sites (SVS) were present in Xa26 (π = 0.01958; SVS = 182) followed by xa5 and Xa21 alleles. The highest frequency of single nucleotide polymorphisms were observed in Xa21 alleles and least in xa5. Transition bias was observed in all the genes and ‘G’ to ‘A’ transitions were more favored than other form of transitions. Neutrality tests failed to show the presence of selection at these loci, though negative Tajima’s D values indicate the presence of a rare form of polymorphisms. At the interspecies level, O. nivara exhibited more diversity than O. sativa. We have also identified two nearly identical resistant alleles of xa5 and two sequentially identical alleles of Xa21. The alleles of xa5 showed basal levels of expression while Xa21 alleles were functionally not expressed.

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1226-9) contains supplementary material, which is available to authorized users.

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

BACKGROUND

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

RESULTS

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10(-7)). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

CONCLUSION

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Developing single nucleotide polymorphism (SNP) markers from transcriptome sequences for identification of longan (Dimocarpus longan) germplasm

Longan (Dimocarpus longan Lour.) is an important tropical fruit tree crop. Accurate varietal identification is essential for germplasm management and breeding. Using longan transcriptome sequences from public databases, we developed single nucleotide polymorphism (SNP) markers; validated 60 SNPs in 50 longan germplasm accessions, including cultivated varieties and wild germplasm; and designated 25 SNP markers that unambiguously identified all tested longan varieties with high statistical rigor (P<0.0001). Multiple trees from the same clone were verified and off-type trees were identified. Diversity analysis revealed genetic relationships among analyzed accessions. Cultivated varieties differed significantly from wild populations (F st=0.300; P<0.001), demonstrating untapped genetic diversity for germplasm conservation and utilization. Within cultivated varieties, apparent differences between varieties from China and those from Thailand and Hawaii indicated geographic patterns of genetic differentiation. These SNP markers provide a powerful tool to manage longan genetic resources and breeding, with accurate and efficient genotype identification.

Plant exomics: concepts, applications and methodologies in crop improvement

Molecular breeding has a crucial role in improvement of crops. Conventional breeding techniques have failed to ameliorate food production. Next generation sequencing has established new concepts of molecular breeding. Exome sequencing has proven to be a significant tool for assessing natural evolution in plants, studying host pathogen interactions and betterment of crop production as exons assist in interpretation of allelic variation with respect to their phenotype. This review covers the platforms for exome sequencing, next generation sequencing technologies that have revolutionized exome sequencing and led toward development of third generation sequencing. Also discussed in this review are the uses of these sequencing technologies to improve wheat, rice and cotton yield and how these technologies are used in exploring the biodiversity of crops, providing better understanding of plant-host pathogen interaction and assessing the process of natural evolution in crops and it also covers how exome sequencing identifies the gene pool involved in symbiotic and other co-existential systems. Furthermore, we conclude how integration of other methodologies including whole genome sequencing, proteomics, transcriptomics and metabolomics with plant exomics covers the areas which are left untouched with exomics alone and in the end how these integration will transform the future of crops.

Characterization of a panel of Vietnamese rice varieties using DArT and SNP markers for association mapping purposes

BACKGROUND

The development of genome-wide association studies (GWAS) in crops has made it possible to mine interesting alleles hidden in gene bank resources. However, only a small fraction of the rice genetic diversity of any given country has been exploited in the studies with worldwide sampling conducted to date. This study presents the development of a panel of rice varieties from Vietnam for GWAS purposes.

RESULTS

The panel, initially composed of 270 accessions, was characterized for simple agronomic traits (maturity class, grain shape and endosperm type) commonly used to classify rice varieties. We first genotyped the panel using Diversity Array Technology (DArT) markers. We analyzed the panel structure, identified two subpanels corresponding to the indica and japonica sub-species and selected 182 non-redundant accessions. However, the number of usable DArT markers (241 for an initial library of 6444 clones) was too small for GWAS purposes. Therefore, we characterized the panel of 182 accessions with 25,971 markers using genotyping by sequencing. The same indica and japonica subpanels were identified. The indica subpanel was further divided into six populations (I1 to I6) using a model-based approach. The japonica subpanel, which was more highly differentiated, was divided into 4 populations (J1 to J4), including a temperate type (J2). Passport data and phenotypic traits were used to characterize these populations. Some populations were exclusively composed of glutinous types (I3 and J2). Some of the upland rice varieties appeared to belong to indica populations, which is uncommon in this region of the world. Linkage disequilibrium decayed faster in the indica subpanel (r2 below 0.2 at 101 kb) than in the japonica subpanel (r2 below 0.2 at 425 kb), likely because of the strongest differentiation of the japonica subpanel. A matrix adapted for GWAS was built by eliminating the markers with a minor allele frequency below 5% and imputing the missing data. This matrix contained 21,814 markers. A GWAS was conducted on time to flowering to prove the utility of this panel.

CONCLUSIONS

This publicly available panel constitutes an important resource giving access to original allelic diversity. It will be used for GWAS on root and panicle traits.

Linkage disequilibrium and association analysis of stripe rust resistance in wild emmer wheat (Triticum turgidum ssp. dicoccoides) population in Israel

Rapid LD decay in wild emmer population from Israel allows high-resolution association mapping. Known and putative new stripe rust resistance genes were found. Genome-wide association mapping (GWAM) is becoming an important tool for the discovery and mapping of loci underlying trait variation in crops, but in the wild relatives of crops the use of GWAM has been limited. Critical factors for the use of GWAM are the levels of linkage disequilibrium (LD) and genetic diversity in mapped populations, particularly in those of self-pollinating species. Here, we report LD estimation in a population of 128 accessions of self-pollinating wild emmer, Triticum turgidum ssp. dicoccoides, the progenitor of cultivated wheat, collected in Israel. LD decayed fast along wild emmer chromosomes and reached the background level within 1 cM. We employed GWAM for the discovery and mapping of genes for resistance to three isolates of Puccinia striiformis, the causative agent of wheat stripe rust. The wild emmer population was genotyped with the wheat iSelect assay including 8643 gene-associated SNP markers (wheat 9K Infinium) of which 2,278 were polymorphic. The significance of association between stripe rust resistance and each of the polymorphic SNP was tested using mixed linear model implemented in EMMA software. The model produced satisfactory results and uncovered four significant associations on chromosome arms 1BS, 1BL and 3AL. The locus on 1BS was located in a region known to contain stripe rust resistance genes. These results show that GWAM is an effective strategy for gene discovery and mapping in wild emmer that will accelerate the utilization of this genetic resource in wheat breeding.

Genetic diversity and population structure analysis of European hexaploid bread wheat (Triticum aestivum L.) varieties

Progress in plant breeding is facilitated by accurate information about genetic structure and diversity. Here, Diversity Array Technology (DArT) was used to characterize a population of 94 bread wheat (Triticum aestivum L.) varieties of mainly European origin. In total, 1,849 of 7,000 tested markers were polymorphic and could be used for population structure analysis. Two major subgroups of wheat varieties, GrI and GrII, were identified using the program STRUCTURE, and confirmed by principal component analysis (PCA). These subgroups were largely separated according to origin; GrI comprised varieties from Southern and Eastern Europe, whereas GrII contained mostly modern varieties from Western and Northern Europe. A large proportion of the markers contributing most to the genetic separation of the subgroups were located on chromosome 2D near the Reduced height 8 (Rht8) locus, and PCR-based genotyping suggested that breeding for the Rht8 allele had a major impact on subgroup separation. Consistently, analysis of linkage disequilibrium (LD) suggested that different selective pressures had acted on chromosome 2D in the two subgroups. Our data provides an overview of the allele composition of bread wheat varieties anchored to DArT markers, which will facilitate targeted combination of alleles following DArT-based QTL studies. In addition, the genetic diversity and distance data combined with specific Rht8 genotypes can now be used by breeders to guide selection of crossing parents.

Analytical method evaluation and discovery of variation within maize varieties in the context of food safety: transcript profiling and metabolomics

Profiling techniques such as microarrays, proteomics, and metabolomics are used widely to assess the overall effects of genetic background, environmental stimuli, growth stage, or transgene expression in plants. To assess the potential regulatory use of these techniques in agricultural biotechnology, we carried out microarray and metabolomic studies of 3 different tissues from 11 conventional maize varieties. We measured technical variations for both microarrays and metabolomics, compared results from individual plants and corresponding pooled samples, and documented variations detected among different varieties with individual plants or pooled samples. Both microarray and metabolomic technologies are reproducible and can be used to detect plant-to-plant and variety-to-variety differences. A pooling strategy lowered sample variations for both microarray and metabolomics while capturing variety-to-variety variation. However, unknown genomic sequences differing between maize varieties might hinder the application of microarrays. High-throughput metabolomics could be useful as a tool for the characterization of transgenic crops. However, researchers will have to take into consideration the impact on the detection and quantitation of a wide range of metabolites on experimental design as well as validation and interpretation of results.

Accurate determination of genetic identity for a single cacao bean, using molecular markers with a nanofluidic system, ensures cocoa authentication

Cacao (Theobroma cacao L.), the source of cocoa, is an economically important tropical crop. One problem with the premium cacao market is contamination with off-types adulterating raw premium material. Accurate determination of the genetic identity of single cacao beans is essential for ensuring cocoa authentication. Using nanofluidic single nucleotide polymorphism (SNP) genotyping with 48 SNP markers, we generated SNP fingerprints for small quantities of DNA extracted from the seed coat of single cacao beans. On the basis of the SNP profiles, we identified an assumed adulterant variety, which was unambiguously distinguished from the authentic beans by multilocus matching. Assignment tests based on both Bayesian clustering analysis and allele frequency clearly separated all 30 authentic samples from the non-authentic samples. Distance-based principle coordinate analysis further supported these results. The nanofluidic SNP protocol, together with forensic statistical tools, is sufficiently robust to establish authentication and to verify gourmet cacao varieties. This method shows significant potential for practical application.

Association analysis of SSR markers with phenology, grain, and Stover-yield related traits in pearl millet (Pennisetum glaucum (L.) R. Br.)

Pearl millet is a staple food crop for millions of people living in the arid and semi-arid tropics. Molecular markers have been used to identify genomic regions linked to traits of interest by conventional QTL mapping and association analysis. Phenotypic recurrent selection is known to increase frequencies of favorable alleles and decrease those unfavorable for the traits under selection. This study was undertaken (i) to quantify the response to recurrent selection for phenotypic traits during breeding of the pearl millet open-pollinated cultivar “CO (Cu) 9” and its four immediate progenitor populations and (ii) to assess the ability of simple sequence repeat (SSR) marker alleles to identify genomic regions linked to grain and stover yield-related traits in these populations by association analysis. A total of 159 SSR alleles were detected across 34 selected single-copy SSR loci. SSR marker data revealed presence of subpopulations. Association analysis identified genomic regions associated with flowering time located on linkage group (LG) 6 and plant height on LG4, LG6, and LG7. Marker alleles on LG6 were associated with stover yield, and those on LG7 were associated with grain yield. Findings of this study would give an opportunity to develop marker-assisted recurrent selection (MARS) or marker-assisted population improvement (MAPI) strategies to increase the rate of gain for pearl millet populations undergoing recurrent selection.

Varietal identification of tea (Camellia sinensis) using nanofluidic array of single nucleotide polymorphism (SNP) markers

Apart from water, tea is the world's most widely consumed beverage. Tea is produced in more than 50 countries with an annual production of approximately 4.7 million tons. The market segment for specialty tea has been expanding rapidly owing to increased demand, resulting in higher revenues and profits for tea growers and the industry. Accurate varietal identification is critically important to ensure traceability and authentication of premium tea products, which in turn contribute to on-farm conservation of tea genetic diversity. Using a set of single nucleotide polymorphism (SNP) markers developed from the expressed sequence tag (EST) database of Camilla senensis, we genotyped deoxyribonucleic acid (DNA) samples extracted from a diverse group of tea varieties, including both fresh and processed commercial loose-leaf teas. The validation led to the designation of 60 SNPs that unambiguously identified all 40 tested tea varieties with high statistical rigor (p<0.0001). Varietal authenticity and genetic relationships among the analyzed cultivars were further characterized by ordination and Bayesian clustering analysis. These SNP markers, in combination with a high-throughput genotyping protocol, effectively established and verified specific DNA fingerprints for all tested tea varieties. This method provides a powerful tool for variety authentication and quality control for the tea industry. It is also highly useful for the management of tea genetic resources and breeding, where accurate and efficient genotype identification is essential.

Population genetics of Braun's Rockcress (Boechera perstellata, Brassicaceae), an endangered plant with a disjunct distribution

Boechera perstellata is an endangered plant found only in middle Tennessee and north central Kentucky. After sampling 4 Tennessee and 3 Kentucky populations, genetic variability and population structure were examined for this species using isozymes, chloroplast DNA, and microsatellites (averaging 35, 29, and 27 individuals per population per class of marker, respectively). The only genetic variability detected for 23 isozymes was a fixed difference between Tennessee and Kentucky populations at 1 locus. Fixed differences between populations of the 2 states were also observed for 3 chloroplast markers. Polymorphism at 19 nuclear microsatellites was 74% at the species level and averaged 21% at the population level. However, observed heterozygosity was extremely low in all populations, ranging from 0.000 to 0.005. High F IS values (0.93) suggest that Boechera perstellata is a primarily selfing species. Tennessee populations have more genetic diversity than Kentucky populations of B. perstellata. Microsatellite markers revealed substantial genetic divergence between the states and genetic differences among populations within each state. Analysis of molecular variance indicates that most variability in this species occurs between the 2 states (49%) and among populations within each state (42%), with relatively little variation found within populations (9%). These data indicate that there is very little gene flow among populations of B. perstellata and that it is important to protect as many populations as possible in order to conserve the genetic diversity of this rare species.

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.

Association mapping combined with linkage analysis for aluminum tolerance among soybean cultivars released in Yellow and Changjiang River Valleys in China

Association mapping (AM) combined with linkage mapping (LM) was executed to identify molecular markers and QTL regions associated with aluminum (Al) tolerance using relative root elongation (RRE) in hydroponics as an indicator. A set of 188 soybean cultivars released in Yellow and Changjiang River Valleys and 184 recombinant inbred lines (RIL) derived from a cross KF No. 1 (tolerant) × NN1138-2 (susceptible) was used in the study. Inheritance analysis of the RIL population suggested four major genes and polygenes controlled Al-tolerance. Further, LM indicated four additive and four epistatic QTL pairs plus a collective unmapped minor QTL were responsible for Al-tolerance and explained 29.39, 18.75 and 43.07 % of the phenotypic variation (PV), respectively. In the set of released cultivars, AM identified 11 markers significant at P < 0.03 that explained 85.2 % of PV with six of which at P < 0.01 accounted for 57.9 % of PV. Ten of these eleven AM marker-QTL were mapped within range of ~2.0 cM to ~43.0 cM outside confidence interval of respective Al-tolerance QTL in previous studies. Five markers, Satt209, Sat_364, Sat_240, Sct_190 and Satt284, were located near Al-tolerance QTL regions in this and previous LM studies. Thus, the two methods confirmed these markers as being the most likely candidate regions for Al-tolerance. Allele effects relative to the population mean for the 11 QTL were estimated, and the allele A210 of Satt209 showed greatest phenotypic effect on Al-tolerance. The two most favorable alleles from each of the 11 marker loci and their carriers were identified, and accordingly the genetic constitution of Al-tolerance for the 188 cultivars was dissected as a QTL-allele matrix. Therefore, marker-assisted pairing of crosses and marker-assisted selection of progenies can be carried out to pyramid favorable alleles of all the 11 loci. This marker-assisted breeding procedure was designated as breeding by design using a QTL-allele matrix.

Microarray-assisted fine-mapping of quantitative trait loci for cold tolerance in rice

Many important agronomic traits, including cold stress resistance, are complex and controlled by quantitative trait loci (QTLs). Isolation of these QTLs will greatly benefit the agricultural industry but it is a challenging task. This study explored an integrated strategy by combining microarray with QTL-mapping in order to identify cold-tolerant QTLs from a cold-tolerant variety IL112 at early-seedling stage. All the early seedlings of IL112 survived normally for 9 d at 4-5°C, while Guichao2 (GC2), an indica cultivar, died after 4 d under the same conditions. Using the F2:3 population derived from the progeny of GC2 and IL112, we identified seven QTLs for cold tolerance. Furthermore, we performed Affymetrix rice whole-genome array hybridization and obtained the expression profiles of IL112 and GC2 under both low-temperature and normal conditions. Four genes were selected as cold QTL-related candidates, based on microarray data mining and QTL-mapping. One candidate gene, LOC_Os07g22494, was shown to be highly associated with cold tolerance in a number of rice varieties and in the F2:3 population, and its overexpression transgenic rice plants displayed strong tolerance to low temperature at early-seedling stage. The results indicated that overexpression of this gene (LOC_Os07g22494) could increase cold tolerance in rice seedlings. Therefore, this study provides a promising strategy for identifying candidate genes in defined QTL regions.

Association mapping and the genomic consequences of selection in sunflower

The combination of large-scale population genomic analyses and trait-based mapping approaches has the potential to provide novel insights into the evolutionary history and genome organization of crop plants. Here, we describe the detailed genotypic and phenotypic analysis of a sunflower (Helianthus annuus L.) association mapping population that captures nearly 90% of the allelic diversity present within the cultivated sunflower germplasm collection. We used these data to characterize overall patterns of genomic diversity and to perform association analyses on plant architecture (i.e., branching) and flowering time, successfully identifying numerous associations underlying these agronomically and evolutionarily important traits. Overall, we found variable levels of linkage disequilibrium (LD) across the genome. In general, islands of elevated LD correspond to genomic regions underlying traits that are known to have been targeted by selection during the evolution of cultivated sunflower. In many cases, these regions also showed significantly elevated levels of differentiation between the two major sunflower breeding groups, consistent with the occurrence of divergence due to strong selection. One of these regions, which harbors a major branching locus, spans a surprisingly long genetic interval (ca. 25 cM), indicating the occurrence of an extended selective sweep in an otherwise recombinogenic interval.

Characterization of genetic diversity and linkage disequilibrium of ZmLOX4 and ZmLOX5 loci in maize

Maize (Zea mays L.) lipoxygenases (ZmLOXs) are well recognized as important players in plant defense against pathogens, especially in cross kingdom lipid communication with pathogenic fungi. This study is among the first to investigate genetic diversity at important gene paralogs ZmLOX4 and ZmLOX5. Sequencing of these genes in 400 diverse maize lines showed little genetic diversity and low linkage disequilibrium in the two genes. Importantly, we identified one inbred line in which ZmLOX5 has a disrupted open reading frame, a line missing ZmLOX5, and five lines with a duplication of ZmLOX5. Tajima's D test suggests that both ZmLOX4 and ZmLOX5 have been under neutral selection. Further investigation of haplotype data revealed that within the ZmLOX family members only ZmLOX12, a monocot specific ZmLOX, showed strong linkage disequilibrium that extends further than expected in maize. Linkage disequilibrium patterns at these loci of interest are crucial for future candidate gene association mapping studies. ZmLOX4 and ZmLOX5 mutations and copy number variants are under further investigation for crop improvement.

A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the 'true citrus fruit trees' group (Citrinae, Rutaceae) and the origin of cultivated species

BACKGROUND AND AIMS

Despite differences in morphology, the genera representing 'true citrus fruit trees' are sexually compatible, and their phylogenetic relationships remain unclear. Most of the important commercial 'species' of Citrus are believed to be of interspecific origin. By studying polymorphisms of 27 nuclear genes, the average molecular differentiation between species was estimated and some phylogenetic relationships between 'true citrus fruit trees' were clarified.

METHODS

Sanger sequencing of PCR-amplified fragments from 18 genes involved in metabolite biosynthesis pathways and nine putative genes for salt tolerance was performed for 45 genotypes of Citrus and relatives of Citrus to mine single nucleotide polymorphisms (SNPs) and indel polymorphisms. Fifty nuclear simple sequence repeats (SSRs) were also analysed.

KEY RESULTS

A total of 16 238 kb of DNA was sequenced for each genotype, and 1097 single nucleotide polymorphisms (SNPs) and 50 indels were identified. These polymorphisms were more valuable than SSRs for inter-taxon differentiation. Nuclear phylogenetic analysis revealed that Citrus reticulata and Fortunella form a cluster that is differentiated from the clade that includes three other basic taxa of cultivated citrus (C. maxima, C. medica and C. micrantha). These results confirm the taxonomic subdivision between the subgenera Metacitrus and Archicitrus. A few genes displayed positive selection patterns within or between species, but most of them displayed neutral patterns. The phylogenetic inheritance patterns of the analysed genes were inferred for commercial Citrus spp.

CONCLUSIONS

Numerous molecular polymorphisms (SNPs and indels), which are potentially useful for the analysis of interspecific genetic structures, have been identified. The nuclear phylogenetic network for Citrus and its sexually compatible relatives was consistent with the geographical origins of these genera. The positive selection observed for a few genes will help further works to analyse the molecular basis of the variability of the associated traits. This study presents new insights into the origin of C. sinensis.

Phenotyping for drought adaptation in wheat using physiological traits

Wheat (Triticum spp) is one of the first domesticated food crops. It represents the first source of calories (after rice) and an important source of proteins in developing countries. As a result of the Green Revolution, wheat yield sharply increased due to the use of improved varieties, irrigation, pesticides, and fertilizers. The rate of increase in world wheat production, however, slowed after 1980, except in China, India, and Pakistan. Being adapted to a wide range of moisture conditions, wheat is grown on more land area worldwide than any other crop, including in drought prone areas. In these marginal rain-fed environments where at least 60 m ha of wheat is grown, amount and distribution of rainfall are the predominant factors influencing yield variability. Intensive work has been carried out in the area of drought adaptation over the last decades. Breeding strategies for drought tolerance improvement include: definition of the target environment, choice and characterization of the testing environment, water stress management and characterization, and use of phenotyping traits with high heritability. The use of integrative traits, facilitated by the development and application of new technologies (thermal imaging, spectral reflectance, stable isotopes) is facilitating high throughput phenotyping and indirect selection, consequently favoring yield improvement in drought prone environments.

Bioenergy grass feedstock: current options and prospects for trait improvement using emerging genetic, genomic, and systems biology toolkits

For lignocellulosic bioenergy to become a viable alternative to traditional energy production methods, rapid increases in conversion efficiency and biomass yield must be achieved. Increased productivity in bioenergy production can be achieved through concomitant gains in processing efficiency as well as genetic improvement of feedstock that have the potential for bioenergy production at an industrial scale. The purpose of this review is to explore the genetic and genomic resource landscape for the improvement of a specific bioenergy feedstock group, the C4 bioenergy grasses. First, bioenergy grass feedstock traits relevant to biochemical conversion are examined. Then we outline genetic resources available bioenergy grasses for mapping bioenergy traits to DNA markers and genes. This is followed by a discussion of genomic tools and how they can be applied to understanding bioenergy grass feedstock trait genetic mechanisms leading to further improvement opportunities.

Identifying loci influencing grain number by microsatellite screening in bread wheat (Triticum aestivum L.)

Grain number (GN) is one of three major yield-related components in wheat. We used the Chinese wheat mini core collection to undertake a genome-wide association analysis of grain number using 531 SSR markers randomly located on all 21 chromosomes. Grain numbers of all accessions were measured in four trials, i.e. two environments in four growing seasons. Association analysis based on a mixed linear model (MLM) revealed that 27 SSR loci were significantly associated with mean GN (MGN) estimated by the best linear unbiased predictor (BLUP) method. These included numerous breeder favorable alleles with strong positive effects at 23 loci. Significant or extremely significant differences were detected on MGN between varieties conveying favored allele and varieties with other alleles. Moreover, statistical simulation showed that the favored alleles have additive genetic effects. Although modern varieties combined larger numbers of favored alleles, the numbers of favored alleles were not significantly different from those in landraces, especially those alleles contributing mostly to the phenotypic variation. These results indicate that there is still considerable genetic potential for use of markers for genome selection of GN for high yield in wheat.

Population structure in a wheat core collection and genomic loci associated with yield under contrasting environments

A set of 96 winter wheat accessions sampled from a variety of geographic origins, including cultivars and breeding lines, were characterized with 46 genome-wide SSR loci for genetic diversity and population structure. The genetic diversity within these accessions was examined using a genetic distance-based and a model-based clustering method. The model-based analysis identified an underlying population structure comprising of four distinct sub-populations which corresponded well with distance-based groupings. Information on the population structure is taken into account in an association mapping study of grain yield from a 3-years field trial incorporating fully irrigated, rainfed and drought stress treatments. A total of 21 marker-grain yield associations (P < 0.01) were identified with nine SSR markers. Most associations were detected only in one to three environments (treatment/year combination), with an average R ( 2 ) value around 13 %. However, marker gwm484 (on chromosome 2D) was associated with yield in six environments, including irrigated, rainfed and drought stress treatments, suggesting it could be used to improve grain yield across a range of environments. Variation in grain yield at this locus was associated with earliness, early vigour, kernels per spikelet and harvest index. Microsatellite locus psp3200 (on chromosome 6D) was associated with yield in dry and hot environments, which was related to earliness, early vigour, productive tillering and total biomass per plant. Partial least squares regression, with nine environmental factors, showed that precipitation from tillering to maturity was the main environmental factor causing marker × environment associations for grain yield.

Genetic diversity, population structure and linkage disequilibrium in elite Chinese winter wheat investigated with SSR markers

To ascertain genetic diversity, population structure and linkage disequilibrium (LD) among a representative collection of Chinese winter wheat cultivars and lines, 90 winter wheat accessions were analyzed with 269 SSR markers distributed throughout the wheat genome. A total of 1,358 alleles were detected, with 2 to 10 alleles per locus and a mean genetic richness of 5.05. The average genetic diversity index was 0.60, with values ranging from 0.05 to 0.86. Of the three genomes of wheat, ANOVA revealed that the B genome had the highest genetic diversity (0.63) and the D genome the lowest (0.56); significant differences were observed between these two genomes (P<0.01). The 90 Chinese winter wheat accessions could be divided into three subgroups based on STRUCTURE, UPGMA cluster and principal coordinate analyses. The population structure derived from STRUCTURE clustering was positively correlated to some extent with geographic eco-type. LD analysis revealed that there was a shorter LD decay distance in Chinese winter wheat compared with other wheat germplasm collections. The maximum LD decay distance, estimated by curvilinear regression, was 17.4 cM (r(2)>0.1), with a whole genome LD decay distance of approximately 2.2 cM (r(2)>0.1, P<0.001). Evidence from genetic diversity analyses suggest that wheat germplasm from other countries should be introduced into Chinese winter wheat and distant hybridization should be adopted to create new wheat germplasm with increased genetic diversity. The results of this study should provide valuable information for future association mapping using this Chinese winter wheat collection.

Patterns of linkage disequilibrium and association mapping in diploid alfalfa (M. sativa L.)

Association mapping enables the detection of marker-trait associations in unstructured populations by taking advantage of historical linkage disequilibrium (LD) that exists between a marker and the true causative polymorphism of the trait phenotype. Our first objective was to understand the pattern of LD decay in the diploid alfalfa genome. We used 89 highly polymorphic SSR loci in 374 unimproved diploid alfalfa (Medicago sativa L.) genotypes from 120 accessions to infer chromosome-wide patterns of LD. We also sequenced four lignin biosynthesis candidate genes (caffeoyl-CoA 3-O-methyltransferase (CCoAoMT), ferulate-5-hydroxylase (F5H), caffeic acid-O-methyltransferase (COMT), and phenylalanine amonialyase (PAL 1)) to identify single nucleotide polymorphisms (SNPs) and infer within gene estimates of LD. As the second objective of this study, we conducted association mapping for cell wall components and agronomic traits using the SSR markers and SNPs from the four candidate genes. We found very little LD among SSR markers implying limited value for genomewide association studies. In contrast, within gene LD decayed within 300 bp below an r (2) of 0.2 in three of four candidate genes. We identified one SSR and two highly significant SNPs associated with biomass yield. Based on our results, focusing association mapping on candidate gene sequences will be necessary until a dense set of genome-wide markers is available for alfalfa.

Genome-wide association mapping in tomato (Solanum lycopersicum) is possible using genome admixture of Solanum lycopersicum var. cerasiforme

Genome-wide association mapping is an efficient way to identify quantitative trait loci controlling the variation of phenotypes, but the approach suffers severe limitations when one is studying inbred crops like cultivated tomato (Solanum lycopersicum). Such crops exhibit low rates of molecular polymorphism and high linkage disequilibrium, which reduces mapping resolution. The cherry type tomato (S. lycopersicum var. cerasiforme) genome has been described as an admixture between the cultivated tomato and its wild ancestor, S. pimpinellifolium. We have thus taken advantage of the properties of this admixture to improve the resolution of association mapping in tomato. As a proof of concept, we sequenced 81 DNA fragments distributed on chromosome 2 at different distances in a core collection of 90 tomato accessions, including mostly cherry type tomato accessions. The 81 Sequence Tag Sites revealed 352 SNPs and indels. Molecular diversity was greatest for S. pimpinellifolium accessions, intermediate for S. l. cerasiforme accessions, and lowest for the cultivated group. We assessed the structure of molecular polymorphism and the extent of linkage disequilibrium over genetic and physical distances. Linkage disequilibrium decreased under r² = 0.3 within 1 cM, and minimal estimated value (r² = 0.13) was reached within 20 kb over the physical regions studied. Associations between polymorphisms and fruit weight, locule number, and soluble solid content were detected. Several candidate genes and quantitative trait loci previously identified were validated and new associations detected. This study shows the advantages of using a collection of S. l. cerasiforme accessions to overcome the low resolution of association mapping in tomato.

An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important food security crops in the tropics and increasingly being adopted for agro-industrial processing. Genetic improvement of cassava can be enhanced through marker-assisted breeding. For this, appropriate genomic tools are required to dissect the genetic architecture of economically important traits. Here, a genome-wide SNP-based genetic map of cassava anchored in SSRs is presented. An outbreeder full-sib (F1) family was genotyped on two independent SNP assay platforms: an array of 1,536 SNPs on Illumina's GoldenGate platform was used to genotype a first batch of 60 F1. Of the 1,358 successfully converted SNPs, 600 which were polymorphic in at least one of the parents and was subsequently converted to KBiosciences' KASPar assay platform for genotyping 70 additional F1. High-precision genotyping of 163 informative SSRs using capillary electrophoresis was also carried out. Linkage analysis resulted in a final linkage map of 1,837 centi-Morgans (cM) containing 568 markers (434 SNPs and 134 SSRs) distributed across 19 linkage groups. The average distance between adjacent markers was 3.4 cM. About 94.2% of the mapped SNPs and SSRs have also been localized on scaffolds of version 4.1 assembly of the cassava draft genome sequence. This more saturated genetic linkage map of cassava that combines SSR and SNP markers should find several applications in the improvement of cassava including aligning scaffolds of the cassava genome sequence, genetic analyses of important agro-morphological traits, studying the linkage disequilibrium landscape and comparative genomics.

Dissecting genetic architecture of grape proanthocyanidin composition through quantitative trait locus mapping

BACKGROUND

Proanthocyanidins (PAs), or condensed tannins, are flavonoid polymers, widespread throughout the plant kingdom, which provide protection against herbivores while conferring organoleptic and nutritive values to plant-derived foods, such as wine. However, the genetic basis of qualitative and quantitative PA composition variation is still poorly understood. To elucidate the genetic architecture of the complex grape PA composition, we first carried out quantitative trait locus (QTL) analysis on a 191-individual pseudo-F1 progeny. Three categories of PA variables were assessed: total content, percentages of constitutive subunits and composite ratio variables. For nine functional candidate genes, among which eight co-located with QTLs, we performed association analyses using a diversity panel of 141 grapevine cultivars in order to identify causal SNPs.

RESULTS

Multiple QTL analysis revealed a total of 103 and 43 QTLs, respectively for seed and skin PA variables. Loci were mainly of additive effect while some loci were primarily of dominant effect. Results also showed a large involvement of pairwise epistatic interactions in shaping PA composition. QTLs for PA variables in skin and seeds differed in number, position, involvement of epistatic interaction and allelic effect, thus revealing different genetic determinisms for grape PA composition in seeds and skin. Association results were consistent with QTL analyses in most cases: four out of nine tested candidate genes (VvLAR1, VvMYBPA2, VvCHI1, VvMYBPA1) showed at least one significant association with PA variables, especially VvLAR1 revealed as of great interest for further functional investigation. Some SNP-phenotype associations were observed only in the diversity panel.

CONCLUSIONS

This study presents the first QTL analysis on grape berry PA composition with a comparison between skin and seeds, together with an association study. Our results suggest a complex genetic control for PA traits and different genetic architectures for grape PA composition between berry skin and seeds. This work also uncovers novel genomic regions for further investigation in order to increase our knowledge of the genetic basis of PA composition.

Association study between the gibberellic acid insensitive gene and leaf length in a Lolium perenne L. synthetic variety

BACKGROUND

Association studies are of great interest to identify genes explaining trait variation since they deal with more than just a few alleles like classical QTL analyses. They are usually performed using collections representing a wide range of variability but which could present a genetic substructure. The aim of this paper is to demonstrate that association studies can be performed using synthetic varieties obtained after several panmictic generations. This demonstration is based on an example of association between the gibberellic acid insensitive gene (GAI) polymorphism and leaf length polymorphism in 'Herbie', a synthetic variety of perennial ryegrass.

METHODS

Leaf growth parameters, consisted of leaf length, maximum leaf elongation rate (LERmax) and leaf elongation duration (LED), were evaluated in spring and autumn on 216 plants of Herbie with three replicates. For each plant, a sequence of 370 bp in GAI was analysed for polymorphism.

RESULTS

Genetic effect was highly significant for all traits. Broad sense heritabilities were higher for leaf length and LERmax with about 0.7 in each period and 0.5 considering both periods than for LED with about 0.4 in each period and 0.3 considering both periods. GAI was highly polymorphic with an average of 12 bp between two consecutive SNPs and 39 haplotypes in which 9 were more frequent. Linkage disequilibrium declined rapidly with distance with r 2 values lower than 0.2 beyond 150 bp. Sequence polymorphism of GAI explained 8-14% of leaf growth parameter variation. A single SNP explained 4% of the phenotypic variance of leaf length in both periods which represents a difference of 33 mm on an average of 300 mm.

CONCLUSIONS

Synthetic varieties in which linkage disequilibrium declines rapidly with distance are suitable for association studies using the "candidate gene" approach. GAI polymorphism was found to be associated with leaf length polymorphism which was more correlated to LERmax than to LED in Herbie. It is a good candidate to explain leaf length variation in other plant material.

Population structure and genetic diversity in a rice core collection (Oryza sativa L.) investigated with SSR markers

The assessment of genetic diversity and population structure of a core collection would benefit to make use of these germplasm as well as applying them in association mapping. The objective of this study were to (1) examine the population structure of a rice core collection; (2) investigate the genetic diversity within and among subgroups of the rice core collection; (3) identify the extent of linkage disequilibrium (LD) of the rice core collection. A rice core collection consisting of 150 varieties which was established from 2260 varieties of Ting's collection of rice germplasm were genotyped with 274 SSR markers and used in this study. Two distinct subgroups (i.e. SG 1 and SG 2) were detected within the entire population by different statistical methods, which is in accordance with the differentiation of indica and japonica rice. MCLUST analysis might be an alternative method to STRUCTURE for population structure analysis. A percentage of 26% of the total markers could detect the population structure as the whole SSR marker set did with similar precision. Gene diversity and MRD between the two subspecies varied considerably across the genome, which might be used to identify candidate genes for the traits under domestication and artificial selection of indica and japonica rice. The percentage of SSR loci pairs in significant (P<0.05) LD is 46.8% in the entire population and the ratio of linked to unlinked loci pairs in LD is 1.06. Across the entire population as well as the subgroups and sub-subgroups, LD decays with genetic distance, indicating that linkage is one main cause of LD. The results of this study would provide valuable information for association mapping using the rice core collection in future.

Association mapping and gene-gene interaction for stem rust resistance in CIMMYT spring wheat germplasm

The recent emergence of wheat stem rust Ug99 and evolution of new races within the lineage threatens global wheat production because they overcome widely deployed stem rust resistance (Sr) genes that had been effective for many years. To identify loci conferring adult plant resistance to races of Ug99 in wheat, we employed an association mapping approach for 276 current spring wheat breeding lines from the International Maize and Wheat Improvement Center (CIMMYT). Breeding lines were genotyped with Diversity Array Technology (DArT) and microsatellite markers. Phenotypic data was collected on these lines for stem rust race Ug99 resistance at the adult plant stage in the stem rust resistance screening nursery in Njoro, Kenya in seasons 2008, 2009 and 2010. Fifteen marker loci were found to be significantly associated with stem rust resistance. Several markers appeared to be linked to known Sr genes, while other significant markers were located in chromosome regions where no Sr genes have been previously reported. Most of these new loci colocalized with QTLs identified recently in different biparental populations. Using the same data and Q + K covariate matrices, we investigated the interactions among marker loci using linear regression models to calculate P values for pairwise marker interactions. Resistance marker loci including the Sr2 locus on 3BS and the wPt1859 locus on 7DL had significant interaction effects with other loci in the same chromosome arm and with markers on chromosome 6B. Other resistance marker loci had significant pairwise interactions with markers on different chromosomes. Based on these results, we propose that a complex network of gene-gene interactions is, in part, responsible for resistance to Ug99. Further investigation may provide insight for understanding mechanisms that contribute to this resistance gene network.

Identification of novel genomic regions associated with resistance to Pyrenophora tritici-repentis races 1 and 5 in spring wheat landraces using association analysis

Tan spot, caused by Pyrenophora tritici-repentis, is a major foliar disease of wheat worldwide. Host plant resistance is the best strategy to manage this disease. Traditionally, bi-parental mapping populations have been used to identify and map quantitative trait loci (QTL) affecting tan spot resistance in wheat. The association mapping (AM) could be an alternative approach to identify QTL based on linkage disequilibrium (LD) within a diverse germplasm set. In this study, we assessed resistance to P. tritici-repentis races 1 and 5 in 567 spring wheat landraces from the USDA-ARS National Small Grains Collection (NSGC). Using 832 diversity array technology (DArT) markers, QTL for resistance to P. tritici-repentis races 1 and 5 were identified. A linear model with principal components suggests that at least seven and three DArT markers were significantly associated with resistance to P. tritici-repentis races 1 and 5, respectively. The DArT markers associated with resistance to race 1 were detected on chromosomes 1D, 2A, 2B, 2D, 4A, 5B, and 7D and explained 1.3-3.1% of the phenotypic variance, while markers associated with resistance to race 5 were distributed on 2D, 6A and 7D, and explained 2.2-5.9% of the phenotypic variance. Some of the genomic regions identified in this study correspond to previously identified loci responsible for resistance to P. tritici-repentis, offering validation for our AM approach. Other regions identified were novel and could possess genes useful for resistance breeding. Some DArT markers associated with resistance to race 1 also were localized in the same regions of wheat chromosomes where QTL for resistance to yellow rust, leaf rust and powdery mildew, have been mapped previously. This study demonstrates that AM can be a useful approach to identify and map novel genomic regions involved in resistance to P. tritici-repentis.

Genetic diversity and population structure in cultivated sunflower and a comparison to its wild progenitor, Helianthus annuus L

Crop germplasm collections are valuable resources for ongoing plant breeding efforts. To fully utilize such collections, however, researchers need detailed information about the amount and distribution of genetic diversity present within collections. Here, we report the results of a population genetic analysis of the primary gene pool of sunflower (Helianthus annuus L.) based on a broad sampling of 433 cultivated accessions from North America and Europe, as well as a range-wide collection of 24 wild sunflower populations. Gene diversity across the cultivars was 0.47, as compared with 0.70 in the wilds, indicating that cultivated sunflower harbors roughly two-thirds of the total genetic diversity present in wild sunflower. Population structure analyses revealed that wild sunflower can be subdivided into four genetically distinct population clusters throughout its North American range, whereas the cultivated sunflower gene pool could be split into two main clusters separating restorer lines from the balance of the gene pool. Use of a maximum likelihood method to estimate the contribution of the wild gene pool to the cultivated sunflower germplasm revealed that the bulk of the cultivar diversity is derived from two wild sunflower population genetic clusters that are primarily composed of individuals from the east-central United States, the same general region in which sunflower domestication is believed to have occurred. We also identified a nested subset of accessions that capture as much of the allelic diversity present within the sampled cultivated sunflower germplasm collection as possible. At the high end, a core set of 288 captured nearly 90% of the alleles present in the full set of 433, whereas a core set of just 12 accessions was sufficient to capture nearly 50% of the total allelic diversity present within this sample of cultivated sunflower.

An analysis of selection on candidate genes for regulation, mobilization, uptake, and transport of iron in maize

Insufficient iron (Fe) availability, which frequently occurs in soils with high pH levels, can lead to leaf chlorosis, a reduced Fe content in harvest products, and yield reduction in maize. The objectives of this study were (i) to describe patterns of sequence variation of 14 candidate genes for mobilization, uptake, and transport of Fe in maize, as well as regulatory function on these processes; (ii) to examine whether Fe-efficiency is an adaptive trait by determining if these genes were targets of selection during domestication; and (iii) to test if the allele distribution at these candidate genes is different for the different subpopulations of maize. The nucleotide diversity of Mtk was reduced by 78% in maize compared with teosinte. The results of our study revealed for the genes Naat1, Nas1, Nramp3, Mtk, and Ys1 a selective sweep, which suggests that these genes might be important for the fast adaptation of maize to new environments with different Fe availabilities.

Population structure in sorghum accessions from West Africa differing in race and maturity class

Accounting for population structure to minimize spurious associations in association analyses is of crucial importance. With sorghum genomic sequence information being available, there is a growing interest in performing such association studies for a number of important agronomic traits using a candidate gene approach. The aims of our study were to conduct a systematic survey of molecular genetic diversity and analyze the population structure in cultivated sorghum [Sorghum bicolor (L.) Moench] accessions from West Africa. Our analysis included 219 West African cultivated sorghum accessions with differing maturity intended for a marker-trait association study. A total of 27 SSRs were used, which resulted in detection of 513 alleles. Genetic diversity estimates for the accessions were found to be high. The accessions were divided into two subgroups using a model-based approach. Our findings partly agree with previous studies in that the guinea race accessions could be distinguished clearly from other accessions included in the analysis. Race and geographical origin of the accessions may be responsible for the structure we observed in our material. The extent of linkage disequilibrium for all combinations of SSRs was in agreement with expectations based on the mating system.

A first insight into population structure and linkage disequilibrium in the US peanut minicore collection

Knowledge of genetic diversity, population structure, and degree of linkage disequilibrium (LD) in target association mapping populations is of great importance and is a prerequisite for LD-based mapping. In the present study, 96 genotypes comprising 92 accessions of the US peanut minicore collection, a component line of the tetraploid variety Florunner, diploid progenitors A. duranensis (AA) and A. ipaënsis (BB), and synthetic amphidiploid accession TxAG-6 were investigated with 392 simple sequence repeat (SSR) marker bands amplified using 32 highly-polymorphic SSR primer pairs. Both distance- and model-based (Bayesian) cluster analysis revealed the presence of structured diversity. In general, the wild-species accessions and the synthetic amphidiploid grouped separately from most minicore accessions except for COC155, and were eliminated from most subsequent analyses. UPGMA analysis divided the population into four subgroups, two major subgroups representing subspecies fastigiata and hypogaea, a third group containing individuals from each subspecies or possibly of mixed ancestry, and a fourth group, either consisting of COC155 alone if wild species were excluded, or of COC155, the diploid species, and the synthetic amphidiploid. Model-based clustering identified four subgroups- one each for fastigiata and hypogaea subspecies, a third consisting of individuals of both subspecies or of mixed ancestry predominantly from Africa or Asia, and a fourth group, consisting of individuals predominantly of var fastigiata, peruviana, and aequatoriana accessions from South America, including COC155. Analysis of molecular variance (AMOVA) revealed statistically-significant (P < 0.0001) genetic variance of 16.87% among subgroups. A total of 4.85% of SSR marker pairs revealed significant LD (at r(2) ≥ 0.1). Of the syntenic marker pairs separated by distances < 10 cM, 11-20 cM, 21-50 cM, and > 50 cM, 19.33, 5.19, 6.25 and 5.29% of marker pairs were found in strong LD (P ≤ 0.01), in accord with LD extending to great distances in self pollinated crops. A threshold value of r(2) > 0.035 was found to distinguish mean r(2) values of linkage distance groups statistically from the mean r(2) values of unlinked markers; LD was found to extend to 10 cM over the entire minicore collection by this criterion. However, there were large differences in r(2) values among marker pairs even among tightly-linked markers. The implications of these findings with regard to the possibility of using association mapping for detection of genome-wide SSR marker-phenotype association are discussed.

Population- and genome-specific patterns of linkage disequilibrium and SNP variation in spring and winter wheat (Triticum aestivum L.)

BACKGROUND

Single nucleotide polymorphisms (SNPs) are ideally suited for the construction of high-resolution genetic maps, studying population evolutionary history and performing genome-wide association mapping experiments. Here, we used a genome-wide set of 1536 SNPs to study linkage disequilibrium (LD) and population structure in a panel of 478 spring and winter wheat cultivars (Triticum aestivum) from 17 populations across the United States and Mexico.

RESULTS

Most of the wheat oligo pool assay (OPA) SNPs that were polymorphic within the complete set of 478 cultivars were also polymorphic in all subpopulations. Higher levels of genetic differentiation were observed among wheat lines within populations than among populations. A total of nine genetically distinct clusters were identified, suggesting that some of the pre-defined populations shared significant proportion of genetic ancestry. Estimates of population structure (F(ST)) at individual loci showed a high level of heterogeneity across the genome. In addition, seven genomic regions with elevated F(ST) were detected between the spring and winter wheat populations. Some of these regions overlapped with previously mapped flowering time QTL. Across all populations, the highest extent of significant LD was observed in the wheat D-genome, followed by lower LD in the A- and B-genomes. The differences in the extent of LD among populations and genomes were mostly driven by differences in long-range LD ( > 10 cM).

CONCLUSIONS

Genome- and population-specific patterns of genetic differentiation and LD were discovered in the populations of wheat cultivars from different geographic regions. Our study demonstrated that the estimates of population structure between spring and winter wheat lines can identify genomic regions harboring candidate genes involved in the regulation of growth habit. Variation in LD suggests that breeding and selection had a different impact on each wheat genome both within and among populations. The higher extent of LD in the wheat D-genome versus the A- and B-genomes likely reflects the episodes of recent introgression and population bottleneck accompanying the origin of hexaploid wheat. The assessment of LD and population structure in this assembled panel of diverse lines provides critical information for the development of genetic resources for genome-wide association mapping of agronomically important traits in wheat.

Population- and genome-specific patterns of linkage disequilibrium and SNP variation in spring and winter wheat (Triticum aestivum L.)

BACKGROUND

Single nucleotide polymorphisms (SNPs) are ideally suited for the construction of high-resolution genetic maps, studying population evolutionary history and performing genome-wide association mapping experiments. Here, we used a genome-wide set of 1536 SNPs to study linkage disequilibrium (LD) and population structure in a panel of 478 spring and winter wheat cultivars (Triticum aestivum) from 17 populations across the United States and Mexico.

RESULTS

Most of the wheat oligo pool assay (OPA) SNPs that were polymorphic within the complete set of 478 cultivars were also polymorphic in all subpopulations. Higher levels of genetic differentiation were observed among wheat lines within populations than among populations. A total of nine genetically distinct clusters were identified, suggesting that some of the pre-defined populations shared significant proportion of genetic ancestry. Estimates of population structure (F(ST)) at individual loci showed a high level of heterogeneity across the genome. In addition, seven genomic regions with elevated F(ST) were detected between the spring and winter wheat populations. Some of these regions overlapped with previously mapped flowering time QTL. Across all populations, the highest extent of significant LD was observed in the wheat D-genome, followed by lower LD in the A- and B-genomes. The differences in the extent of LD among populations and genomes were mostly driven by differences in long-range LD ( > 10 cM).

CONCLUSIONS

Genome- and population-specific patterns of genetic differentiation and LD were discovered in the populations of wheat cultivars from different geographic regions. Our study demonstrated that the estimates of population structure between spring and winter wheat lines can identify genomic regions harboring candidate genes involved in the regulation of growth habit. Variation in LD suggests that breeding and selection had a different impact on each wheat genome both within and among populations. The higher extent of LD in the wheat D-genome versus the A- and B-genomes likely reflects the episodes of recent introgression and population bottleneck accompanying the origin of hexaploid wheat. The assessment of LD and population structure in this assembled panel of diverse lines provides critical information for the development of genetic resources for genome-wide association mapping of agronomically important traits in wheat.