Genome-wide SNP Genotyping Resolves Signatures of Selection and Tetrasomic Recombination in Peanut.

Peanut (Arachis hypogaea; 2n = 4x = 40) is a nutritious food and a good source of vitamins, minerals, and healthy fats. Expansion of genetic and genomic resources for genetic enhancement of cultivated peanut has gained momentum from the sequenced genomes of the diploid ancestors of cultivated peanut. To facilitate high-throughput genotyping of Arachis species, 20 genotypes were re-sequenced and genome-wide single nucleotide polymorphisms (SNPs) were selected to develop a large-scale SNP genotyping array. For flexibility in genotyping applications, SNPs polymorphic between tetraploid and diploid species were included for use in cultivated and interspecific populations. A set of 384 accessions was used to test the array resulting in 54 564 markers that produced high-quality polymorphic clusters between diploid species, 47 116 polymorphic markers between cultivated and interspecific hybrids, and 15 897 polymorphic markers within A. hypogaea germplasm. An additional 1193 markers were identified that illuminated genomic regions exhibiting tetrasomic recombination. Furthermore, a set of elite cultivars that make up the pedigree of US runner germplasm were genotyped and used to identify genomic regions that have undergone positive selection. These observations provide key insights on the inclusion of new genetic diversity in cultivated peanut and will inform the development of high-resolution mapping populations. Due to its efficiency, scope, and flexibility, the newly developed SNP array will be very useful for further genetic and breeding applications in Arachis.

Genome resources for climate-resilient cowpea, an essential crop for food security.

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Genetic mapping of QTLs controlling fatty acids provided insights into the genetic control of fatty acid synthesis pathway in peanut (Arachis hypogaea L.).

Peanut, a high-oil crop with about 50% oil content, is either crushed for oil or used as edible products. Fatty acid composition determines the oil quality which has high relevance to consumer health, flavor, and shelf life of commercial products. In addition to the major fatty acids, oleic acid (C18:1) and linoleic acid (C18:2) accounting for about 80% of peanut oil, the six other fatty acids namely palmitic acid (C16:0), stearic acid (C18:0), arachidic acid (C20:0), gadoleic acid (C20:1), behenic acid (C22:0), and lignoceric acid (C24:0) are accounted for the rest 20%. To determine the genetic basis and to improve further understanding on effect of FAD2 genes on these fatty acids, two recombinant inbred line (RIL) populations namely S-population (high oleic line 'SunOleic 97R' × low oleic line 'NC94022') and T-population (normal oleic line 'Tifrunner' × low oleic line 'GT-C20') were developed. Genetic maps with 206 and 378 marker loci for the S- and the T-population, respectively were used for quantitative trait locus (QTL) analysis. As a result, a total of 164 main-effect (M-QTLs) and 27 epistatic (E-QTLs) QTLs associated with the minor fatty acids were identified with 0.16% to 40.56% phenotypic variation explained (PVE). Thirty four major QTLs (>10% of PVE) mapped on five linkage groups and 28 clusters containing more than three QTLs were also identified. These results suggest that the major QTLs with large additive effects would play an important role in controlling composition of these minor fatty acids in addition to the oleic and linoleic acids in peanut oil. The interrelationship among these fatty acids should be considered while breeding for improved peanut genotypes with good oil quality and desired fatty acid composition.

Identification of QTLs associated with oil content and mapping FAD2 genes and their relative contribution to oil quality in peanut (Arachis hypogaea L.).

BACKGROUND: Peanut is one of the major source for human consumption worldwide and its seed contain approximately 50% oil. Improvement of oil content and quality traits (high oleic and low linoleic acid) in peanut could be accelerated by exploiting linked markers through molecular breeding. The objective of this study was to identify QTLs associated with oil content, and estimate relative contribution of FAD2 genes (ahFAD2A and ahFAD2B) to oil quality traits in two recombinant inbred line (RIL) populations. RESULTS: Improved genetic linkage maps were developed for S-population (SunOleic 97R × NC94022) with 206 (1780.6 cM) and T-population (Tifrunner × GT-C20) with 378 (2487.4 cM) marker loci. A total of 6 and 9 QTLs controlling oil content were identified in the S- and T-population, respectively. The contribution of each QTL towards oil content variation ranged from 3.07 to 10.23% in the S-population and from 3.93 to 14.07% in the T-population. The mapping positions for ahFAD2A (A sub-genome) and ahFAD2B (B sub-genome) genes were assigned on a09 and b09 linkage groups. The ahFAD2B gene (26.54%, 25.59% and 41.02% PVE) had higher phenotypic effect on oleic acid (C18:1), linoleic acid (C18:2), and oleic/linoleic acid ratio (O/L ratio) than ahFAD2A gene (8.08%, 6.86% and 3.78% PVE). The FAD2 genes had no effect on oil content. This study identified a total of 78 main-effect QTLs (M-QTLs) with up to 42.33% phenotypic variation (PVE) and 10 epistatic QTLs (E-QTLs) up to 3.31% PVE for oil content and quality traits. CONCLUSIONS: A total of 78 main-effect QTLs (M-QTLs) and 10 E-QTLs have been detected for oil content and oil quality traits. One major QTL (more than 10% PVE) was identified in both the populations for oil content with source alleles from NC94022 and GT-C20 parental genotypes. FAD2 genes showed high effect for oleic acid (C18:1), linoleic acid (C18:2), and O/L ratio while no effect on total oil content. The information on phenotypic effect of FAD2 genes for oleic acid, linoleic acid and O/L ratio, and oil content will be applied in breeding selection.

Phylogenetic relationships of species of genus Arachis based on genic sequences.

The genus Arachis (Fabaceae), which originated in South America, consists of 80 species. Based on morphological traits and cross-compatibility among the species, the genus is divided into nine taxonomic sections. Arachis is the largest section including the economically valuable cultivated peanut (A. hypogaea). Seven genic sequences were utilized to better understand the phylogenetic relationships between species of genus Arachis. Our study displayed four clades of species of Arachis. Arachis triseminata was genetically isolated from all other species of Arachis studied, and it formed the basal clade with A. retusa and A. dardani from the most ancient sections Extranervosae and Heteranthae, respectively. Species of section Arachis formed a separated single clade from all other species, within which species having B and D genome clustered in one subgroup and three species characterized with an A genome grouped together in another subgroup. A divergent clade including species from five sections was sister to the clade of section Arachis. Between the sister clades and the basal clade there was a clade containing species from the more advanced sections. Phylogenetic relationships of all the species of Arachis using multiple genic sequences were similar to the phylogenies produced with single-copy genes.

Genotypic effect of ahFAD2 on fatty acid profiles in six segregating peanut (Arachis hypogaea L) populations.

BACKGROUND: Fatty acid composition of oil extracted from peanut (Arachis hypogaea L.) seed is an important quality trait because it may affect the flavor and shelf life of resulting food products. In particular, a high ratio of oleic (C18:1) relative to linoleic (C18:2) fatty acid (O/L ≥ 10) results in a longer shelf life. Previous reports suggest that the high oleic (~80%) trait was controlled by recessive alleles of ahFAD2A and ahFAD2B, the former of which is thought to have a high frequency in US runner- and virginia-type cultivars. Functional mutations, G448A in ahFAD2A and 442insA in ahFAD2B eliminate or knock down desaturase activity and have been demonstrated to produce peanut oil with high O/L ratios. In order to employ marker assisted selection (MAS) to select a high oleic disease resistant peanut and to evaluate genotypic and phenotypic variation, crosses were made between high oleic (~80%) and normal oleic (~50%) peanuts to produce segregating populations. RESULTS: A total of 539 F2 progenies were randomly selected to empirically determine each ahFAD2 genotype and the resulting fatty acid composition. Five of the six crosses segregated for the high oleic trait in a digenic fashion. The remaining cross was consistent with monogenic segregation because both parental genotypes were fixed for the ahFAD2A mutation. Segregation distortion was significant in ahFAD2A in one cross; however, the remaining crosses showed no distortion. Quantitative analyses revealed that dominance was incomplete for the wild type allele of ahFAD2, and both loci showed significant additive effects. Oleic and linoleic acid displayed five unique phenotypes, based on the number of ahFAD2 mutant alleles. Further, the ahFAD2 loci did exhibit pleiotropic interactions with palmitic (C16:0), oleic (C18:1), linoleic (C18:2) acids and the O/L ratio. Fatty acid levels in these progeny were affected by the parental genotype suggesting that other genes also influence fatty acid composition in peanut. As far as the authors are aware, this is the first study in which all of the nine possible ahFAD2 genotypes were quantitatively measured. CONCLUSIONS: The inheritance of the high oleic trait initially was suggested to be controlled by dominant gene action from two homoeologous genes (ahFAD2A and ahFAD2B) exhibiting complete recessivity. Analyzing the ahFAD2 genotypes and fatty acid compositions of these segregating peanut populations clearly demonstrated that the fatty acid contents are quantitative in nature although much of the variability in the predominant fatty acids (oleic, linoleic, and palmitic) is controlled by only two loci.

Oil, fatty acid, flavonoid, and resveratrol content variability and FAD2A functional SNP genotypes in the U.S. peanut mini-core collection.

Peanut seeds contain high amounts of oil and protein as well as some useful bioactive phytochemicals which can contribute to human health. The U.S. peanut mini-core collection is an important genetic resource for improving seed quality and developing new cultivars. Variability of seed chemical composition within the mini-core was evaluated from freshly harvested seeds for two years. Oil, fatty acid composition, and flavonoid/resveratrol content were quantified by NMR, GC, and HPLC, respectively. Significant variability was detected in seed chemical composition among accessions and botanical varieties. Accessions were further genotyped with a functional SNP marker from the FAD2A gene using real-time PCR and classified into three genotypes with significantly different O/L ratios: wild type (G/G with a low O/L ratio <1.7), heterozygote (G/A with O/L ratio >1.4 but <1.7), and mutant (A/A with a high O/L ratio >1.7). The results from real-time PCR genotyping and GC fatty acid analysis were consistent. Accessions with high amounts of oil, quercetin, high seed weight, and O/L ratio were identified. The results from this study may be useful not only for peanut breeders, food processors, and product consumers to select suitable accessions or cultivars but also for curators to potentially expand the mini-core collection.

Always look on both sides: phylogenetic information conveyed by simple sequence repeat allele sequences.

Simple sequence repeat (SSR) markers are widely used tools for inferences about genetic diversity, phylogeography and spatial genetic structure. Their applications assume that variation among alleles is essentially caused by an expansion or contraction of the number of repeats and that, accessorily, mutations in the target sequences follow the stepwise mutation model (SMM). Generally speaking, PCR amplicon sizes are used as direct indicators of the number of SSR repeats composing an allele with the data analysis either ignoring the extent of allele size differences or assuming that there is a direct correlation between differences in amplicon size and evolutionary distance. However, without precisely knowing the kind and distribution of polymorphism within an allele (SSR and the associated flanking region (FR) sequences), it is hard to say what kind of evolutionary message is conveyed by such a synthetic descriptor of polymorphism as DNA amplicon size. In this study, we sequenced several SSR alleles in multiple populations of three divergent tree genera and disentangled the types of polymorphisms contained in each portion of the DNA amplicon containing an SSR. The patterns of diversity provided by amplicon size variation, SSR variation itself, insertions/deletions (indels), and single nucleotide polymorphisms (SNPs) observed in the FRs were compared. Amplicon size variation largely reflected SSR repeat number. The amount of variation was as large in FRs as in the SSR itself. The former contributed significantly to the phylogenetic information and sometimes was the main source of differentiation among individuals and populations contained by FR and SSR regions of SSR markers. The presence of mutations occurring at different rates within a marker's sequence offers the opportunity to analyse evolutionary events occurring on various timescales, but at the same time calls for caution in the interpretation of SSR marker data when the distribution of within-locus polymorphism is not known.

Population structure and marker-trait association analysis of the US peanut (Arachis hypogaea L.) mini-core collection.

Peanut (Arachis hypogaea L.) is one of the most important oilseed and nutritional crops in the world. To efficiently utilize the germplasm collection, a peanut mini-core containing 112 accessions was established in the United States. To determine the population structure and its impact on marker-trait association, this mini-core collection was assessed by genotyping 94 accessions with 81 SSR markers and two functional SNP markers from fatty acid desaturase 2 (FAD2). Seed quality traits (including oil content, fatty acid composition, flavonoids, and resveratrol) were obtained through nuclear magnetic resonance (NMR), gas chromatography (GC), and high-performance liquid chromatography (HPLC) analysis. Genetic diversity and population structure analysis identified four major subpopulations that are related to four botanical varieties. Model comparison with different levels of population structure and kinship control was conducted for each trait and association analyses with the selected models verified that the functional SNP from the FAD2A gene is significantly associated with oleic acid (C18:1), linoleic acid (C18:2), and oleic-to-linoleic (O/L) ratio across this diverse collection. Even though the allele distribution of FAD2A was structured among the four subpopulations, the effect of FAD2A gene remained significant after controlling population structure and had a likelihood-ratio-based R ( 2 ) (R ( LR ) ( 2 ) ) value of 0.05 (oleic acid), 0.09 (linoleic acid), and 0.07 (O/L ratio) because the FAD2A alleles were not completely fixed within subpopulations. Our genetic analysis demonstrated that this peanut mini-core panel is suitable for association mapping. Phenotypic characterization for seed quality traits and association testing of the functional SNP from FAD2A gene provided information for further breeding and genetic research.

FAD2 gene mutations significantly alter fatty acid profiles in cultivated peanuts (Arachis hypogaea).

A panel of 55 peanut lines was analyzed for fatty acid composition with gas chromatography and also genotyped with SNP markers from the FAD2 genes by real-time PCR. Significant variation in fatty acid composition was identified, and the ratio of oleic acid to linoleic acid (O/L) ranged from 1.23 to 56.45. In terms of the FAD2 gene mutation, the assayed lines were classified into four genotypes: wild type (Ol(1)Ol(1)Ol(2)Ol(2)), single functional homozygous mutation on the A genome (ol(1)ol(1)Ol(2)Ol(2)), single functional homozygous mutation on the B genome (Ol(1)Ol(1)ol(2)ol(2)), and a double mutation on both A and B genomes (ol(1)ol(1)ol(2)ol(2)). Each genotype has a significantly different fatty acid profile. Both FAD2A and FAD2B are involved in the conversion of oleic acid to linoleic acid in peanuts. Overall, these results demonstrate the combined power of genetic analysis with biochemical analysis on peanut fatty acid research.

A real-time PCR genotyping assay to detect FAD2A SNPs in peanuts (Arachis hypogaea L)

The high oleic (C18:1) phenotype in peanuts has been previously demonstrated to result from a homozygous recessive genotype (ol1ol1ol2ol2) in two homeologous fatty acid desaturase genes (FAD2A and FAD2B) with two key SNPs. These mutant SNPs, specifically G448A in FAD2A and 442insA in FAD2B, significantly limit the normal function of the desaturase enzyme activity which converts oleic acid into linoleic acid by the addition of a second double bond in the hydrocarbon chain. Previously, a genotyping assay was developed to detect wild type and mutant alleles in FAD2B. A real-time PCR assay has now been developed to detect wild type and mutant alleles (G448A) in FAD2A using either seed or leaf tissue. This assay was demonstrated to be applicable for the detection of homozygous and heterozygous samples. The FAD2A genotyping assay was validated by employing gas chromatography (GC) to determine total fatty acid composition and by genotyping peanut lines that have been well characterized. Overall, development of rapid assays such as real-time PCR which can identify key genotypes associated with important agronomic traits such as oleic acid, will improve breeding efficiency by targeting desirable genotypes at early stages of development.

Genetic diversity and population structure analysis of accessions in the US historic sweet sorghum collection.

Sweet sorghum has the potential to become a versatile feedstock for large-scale bioenergy production given its sugar from stem juice, cellulose/hemicellulose from stalks, and starch from grain. However, for researchers to maximize its feedstock potential a first step includes additional evaluations of the 2,180 accessions with varied origins in the US historic sweet sorghum collection. To assess genetic diversity of this collection for bioenergy breeding and population structure for association mapping, we selected 96 accessions and genotyped them with 95 simple sequence repeat markers. Subsequent genetic diversity and population structure analysis methods identified four subpopulations in this panel, which correlated well with the geographic locations where these accessions originated or were collected. Model comparisons for three quantitative traits revealed different levels of population structure effects on flowering time, plant height, and brix. Our results suggest that diverse germplasm accessions curated from different geographical regions should be considered for plant breeding programs to develop sweet sorghum cultivars or hybrids, and that this sweet sorghum panel can be further explored for association mapping.

Applying genotyping (TILLING) and phenotyping analyses to elucidate gene function in a chemically induced sorghum mutant population.

BACKGROUND: Sorghum [Sorghum bicolor (L.) Moench] is ranked as the fifth most important grain crop and serves as a major food staple and fodder resource for much of the world, especially in arid and semi-arid regions. The recent surge in sorghum research is driven by its tolerance to drought/heat stresses and its strong potential as a bioenergy feedstock. Completion of the sorghum genome sequence has opened new avenues for sorghum functional genomics. However, the availability of genetic resources, specifically mutant lines, is limited. Chemical mutagenesis of sorghum germplasm, followed by screening for mutants altered in important agronomic traits, represents a rapid and effective means of addressing this limitation. Induced mutations in novel genes of interest can be efficiently assessed using the technique known as Targeting Induced Local Lesion IN Genomes (TILLING). RESULTS: A sorghum mutant population consisting of 1,600 lines was generated from the inbred line BTx623 by treatment with the chemical agent ethyl methanesulfonate (EMS). Numerous phenotypes with altered morphological and agronomic traits were observed from M2 and M3 lines in the field. A subset of 768 mutant lines was analyzed by TILLING using four target genes. A total of five mutations were identified resulting in a calculated mutation density of 1/526 kb. Two of the mutations identified by TILLING and verified by sequencing were detected in the gene encoding caffeic acid O-methyltransferase (COMT) in two independent mutant lines. The two mutant lines segregated for the expected brown midrib (bmr) phenotype, a trait associated with altered lignin content and increased digestibility. CONCLUSION: TILLING as a reverse genetic approach has been successfully applied to sorghum. The diversity of the mutant phenotypes observed in the field, and the density of induced mutations calculated from TILLING indicate that this mutant population represents a useful resource for members of the sorghum research community. Moreover, TILLING has been demonstrated to be applicable for sorghum functional genomics by evaluating a small subset of the EMS-induced mutant lines.

Discovering and verifying DNA polymorphisms in a mung bean [V. radiata (L.) R. Wilczek] collection by EcoTILLING and sequencing.

BACKGROUND: Vigna radiata, which is classified in the family Fabaceae, is an important economic crop and a dietary staple in many developing countries. The species radiata can be further subdivided into varieties of which the variety sublobata is currently acknowledged as the putative progenitor of radiata. EcoTILLING was employed to identify single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELS) in a collection of Vigna radiata accessions. FINDINGS: A total of 157 DNA polymorphisms in the collection were produced from ten primer sets when using V. radiata var. sublobata as the reference. The majority of polymorphisms detected were found in putative introns. The banding patterns varied from simple to complex as the number of DNA polymorphisms between two pooled samples increased. Numerous SNPs and INDELS ranging from 4-24 and 1-6, respectively, were detected in all fragments when pooling V. radiata var. sublobata with V. radiata var. radiata. On the other hand, when accessions of V. radiata var. radiata were mixed together and digested with CEL I relatively few SNPs and no INDELS were detected. CONCLUSION: EcoTILLING was utilized to identify polymorphisms in a collection of mung bean, which previously showed limited molecular genetic diversity and limited morphological diversity in the flowers and pod descriptors. Overall, EcoTILLING proved to be a powerful genetic analysis tool providing the rapid identification of naturally occurring variation.

Phylogenetic relationships and genetic diversity of the USDA Vigna germplasm collection revealed by gene-derived markers and sequencing.

Phylogenetic relationships in the USDA Vigna germplasm collection are somewhat unclear and their genetic diversity has not been measured empirically. To reveal interspecific phylogenetic relationships and assess their genetic diversity, 48 accessions representing 12 Vigna species were selected, and 30 gene-derived markers from legumes were employed. Some high-quality amplicons were sequenced. Indels (insertion/deletions) were discovered from the sequence alignments that were specific identifiers for some Vigna species. With regard to revealing polymorphisms, intron-spanning markers were more effective than exon-derived markers. These gene-derived markers were more successful in revealing interspecific polymorphisms than intraspecific polymorphisms at both the DNA fragment and sequence levels. Two different dendrograms were generated from DNA fragment data and sequence data, respectively. The results from these two dendrograms supported each other and showed similar phylogenetic relationships among the Vigna species investigated. The accessions clustered into four main groups and 13 subgroups. Each subgroup represents a subgenus or a species. Phylogenetic analysis revealed that an accession might be misclassified in our collection. The putative misclassified accession was further supported by seed morphology. Limited intraspecific genetic diversity was revealed by these gene-derived markers and/or sequences. The USDA Vigna germplasm collection currently consists of multiple species with many accessions further classified into specific subspecies, but very few subspecies of the total subspecies available exist within the collection. Based on our results, more attention should be paid to the subspecies, wild forms and/or botanical varieties for future curation in order to expand the genetic diversity of Vigna germplasm in the USDA collection.

Genetic diversity of cultivated and wild-type peanuts evaluated with M13-tailed SSR markers and sequencing.

Thirty-one genomic SSR markers with a M13 tail attached were used to assess the genetic diversity of the peanut mini core collection. The M13-tailed method was effective in discriminating almost all the cultivated and wild accessions. A total of 477 alleles were detected with an average of 15.4 alleles per locus. The mean polymorphic information content (PIC) score was 0.687. The cultivated peanut (Arachis hypogaea L.) mini core produced a total of 312 alleles with an average of 10.1 alleles per locus. A neighbour-joining tree was constructed to determine the interspecific and intraspecific relationships in this data set. Almost all the peanut accessions in this data set classified into subspecies and botanical varieties such as subsp. hypogaea var. hypogaea, subsp. fastigiata var. fastigiata, and subsp. fastigiata var. vulgaris clustered with other accessions with the same classification, which lends further support to their current taxonomy. Alleles were sequenced from one of the SSR markers used in this study, which demonstrated that the repeat motif is conserved when transferring the marker across species borders. This study allowed the examination of the diversity and phylogenetic relationships in the peanut mini core which has not been previously reported.

Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs).

Twenty-four simple sequence repeat (SSR) markers were used to detect molecular polymorphisms among 370 mostly sexually derived Citrus accessions from the collection of citrus germplasm maintained at the University of California, Riverside. A total of 275 alleles were detected with an average of 11.5 alleles per locus and an average polymorphism information content of 0.625. Genetic diversity statistics were calculated for each individual SSR marker, the entire population, and for specified Citrus groups. Phylogenetic relationships among all citrus accessions and putative non-hybrid Citrus accessions were determined by constructing neighbor-joining trees. There was strong support for monophyly at the species level when hybrid taxa were removed from the data set. Both of these trees indicate that Fortunella clusters within the genus Citrus but Poncirus is a sister genus to Citrus. Additionally, Citrus accessions were probabilistically assigned to populations or multiple populations if their genotype indicated an admixture by a model-based clustering approach. This approach identified five populations in this data set. These separate analyses (distance and model based) both support the hypothesis that there are only a few naturally occurring species of Citrus and most other types of Citrus arose through various hybridization events between these naturally occurring forms.