Crop DNA extraction with lab-made magnetic nanoparticles

Molecular breeding methods, such as marker-assisted selection and genomic selection, require high-throughput and cost-effective methods for isolating genomic DNA from plants, specifically from crop tissue or seed with high polysaccharides, lipids, and proteins. A quick and inexpensive high-throughput method for isolating genomic DNA from seed and leaf tissue from multiple crops was tested with a DNA isolation method that combines CTAB extraction buffer and lab-made SA-coated magnetic nanoparticles. This method is capable of isolating quality genomic DNA from leaf tissue and seeds in less than 2 hours with fewer steps than a standard CTAB extraction method. The yield of the genomic DNA was 582-729 ng per 5 leaf discs or 216-1869 ng per seed in soybean, 2.92-62.6 ng per 5 leaf discs or 78.9-219 ng per seed in wheat, and 30.9-35.4 ng per 5 leaf discs in maize. The isolated DNA was tested with multiple molecular breeding methods and was found to be of sufficient quality and quantity for PCR and targeted genotyping by sequencing methods such as molecular inversion probes (MIPs). The combination of SA-coated magnetic nanoparticles and CTAB extraction buffer is a fast, simple, and environmentally friendly, high-throughput method for both leaf tissues and seed(s) DNA preparation at low cost per sample. The DNA obtained from this method can be deployed in applied breeding programs for marker-assisted selection or genomic selection.

Recombination Hotspots in Soybean [Glycine Max (L.) Merr.]

Recombination allows for the exchange of genetic material between two parents which plant breeders exploit to make improved cultivars. This recombination is not distributed evenly across the chromosome. Recombination mostly occurs in euchromatic regions of the genome and even then, recombination is focused into clusters of crossovers termed recombination hotspots. Understanding the distribution of these hotspots along with the sequence motifs associated with them may lead to methods that enable breeders to better exploit recombination in breeding. To map recombination hotspots and identify sequence motifs associated with hotspots in soybean [Glycine max (L.) Merr.], two bi-parental recombinant inbred lines populations were genotyped with the SoySNP50k Illumina Infinium assay. A total of 451 recombination hotspots were identified in the two populations. Despite being half-sib populations, only 18 hotspots were in common between the two populations. While pericentromeric regions did exhibit extreme suppression of recombination, twenty-seven percent of the detected hotspots were located in the pericentromeric regions of the chromosomes. Two genomic motifs associated with hotspots are similar to human, dog, rice, wheat, drosophila, and arabidopsis. These motifs were a CCN repeat motif and a poly-A motif. Genomic regions spanning other hotspots were significantly enriched with the tourist family of mini-inverted-repeat transposable elements that resides in less than 0.34% of the soybean genome. The characterization of recombination hotspots in these two large soybean bi-parental populations demonstrates that hotspots do occur throughout the soybean genome and are enriched for specific motifs but their locations may not be conserved between different populations.

Development of molecular inversion probes for soybean progeny genomic selection genotyping

Increasing rate of genetic gain for key agronomic traits through genomic selection requires the development of new molecular methods to run genome-wide single-nucleotide polymorphisms (SNPs). The main limitation of current methods is the cost is too high to screen breeding populations. Molecular inversion probes (MIPs) are a targeted genotyping-by-sequencing (GBS) method that could be used for soybean [Glycine max (L.) Merr.] that is both cost-effective, high-throughput, and provides high data quality to screen breeder's germplasm for genomic selection. A 1K MIP SNP set was developed for soybean with uniformly distributed markers across the genome. The SNPs were selected to maximize the number of informative markers in germplasm being tested in soybean breeding programs located in the northern-central and middle-southern regions of the United States. The 1K SNP MIP set was tested on diverse germplasm and a recombinant inbred line (RIL) population. Targeted sequencing with MIPs obtained an 85% enrichment for the targeted SNPs. The MIP genotyping accuracy was 93% overall, whereas homozygous call accuracy was 98% with <10% missing data. The accuracy of MIPs combined with its low per-sample cost makes it a powerful tool to enable genomic selection within soybean breeding programs.

Candidate Genes Modulating Reproductive Timing in Elite US Soybean Lines Identified in Soybean Alleles of Arabidopsis Flowering Orthologs With Divergent Latitude Distribution

Adaptation of soybean cultivars to the photoperiod in which they are grown is critical for optimizing plant yield. However, despite its importance, only the major loci conferring variation in flowering time and maturity of US soybean have been isolated. By contrast, over 200 genes contributing to floral induction in the model organism Arabidopsis thaliana have been described. In this work, putative alleles of a library of soybean orthologs of these Arabidopsis flowering genes were tested for their latitudinal distribution among elite US soybean lines developed in the United States. Furthermore, variants comprising the alleles of genes with significant differences in latitudinal distribution were assessed for amino acid conservation across disparate genera to infer their impact on gene function. From these efforts, several candidate genes from various biological pathways were identified that are likely being exploited toward adaptation of US soybean to various maturity groups.

Soybean BARCSoySNP6K: An assay for soybean genetics and breeding research

The limited number of recombinant events in recombinant inbred lines suggests that for a biparental population with a limited number of recombinant inbred lines, it is unnecessary to genotype the lines with many markers. For genomic prediction and selection, previous studies have demonstrated that only 1000-2000 genome-wide common markers across all lines/accessions are needed to reach maximum efficiency of genomic prediction in populations. Evaluation of too many markers will not only increase the cost but also generate redundant information. We developed a soybean (Glycine max) assay, BARCSoySNP6K, containing 6000 markers, which were carefully chosen from the SoySNP50K assay based on their position in the soybean genome and haplotype block, polymorphism among accessions and genotyping quality. The assay includes 5000 single nucleotide polymorphisms (SNPs) from euchromatic and 1000 from heterochromatic regions. The percentage of SNPs with minor allele frequency >0.10 was 95% and 91% in the euchromatic and heterochromatic regions, respectively. Analysis of progeny from two large families genotyped with SoySNP50K versus BARCSoySNP6K showed that the position of the common markers and number of unique bins along linkage maps were consistent based on the SNPs genotyped with the two assays; however, the rate of redundant markers was dramatically reduced with the BARCSoySNP6K. The BARCSoySNP6K assay is proven as an excellent tool for detecting quantitative trait loci, genomic selection and assessing genetic relationships. The assay is commercialized by Illumina Inc. and being used by soybean breeders and geneticists and the list of SNPs in the assay is an ideal resource for SNP genotyping by targeted amplicon sequencing.

Genetic Architecture of Soybean Yield and Agronomic Traits

Soybean is the world's leading source of vegetable protein and demand for its seed continues to grow. Breeders have successfully increased soybean yield, but the genetic architecture of yield and key agronomic traits is poorly understood. We developed a 40-mating soybean nested association mapping (NAM) population of 5,600 inbred lines that were characterized by single nucleotide polymorphism (SNP) markers and six agronomic traits in field trials in 22 environments. Analysis of the yield, agronomic, and SNP data revealed 23 significant marker-trait associations for yield, 19 for maturity, 15 for plant height, 17 for plant lodging, and 29 for seed mass. A higher frequency of estimated positive yield alleles was evident from elite founder parents than from exotic founders, although unique desirable alleles from the exotic group were identified, demonstrating the value of expanding the genetic base of US soybean breeding.

Genetic Characterization of the Soybean Nested Association Mapping Population

A set of nested association mapping (NAM) families was developed by crossing 40 diverse soybean [ (L.) Merr.] genotypes to the common cultivar. The 41 parents were deeply sequenced for SNP discovery. Based on the polymorphism of the single-nucleotide polymorphisms (SNPs) and other selection criteria, a set of SNPs was selected to be included in the SoyNAM6K BeadChip for genotyping the parents and 5600 RILs from the 40 families. Analysis of the SNP profiles of the RILs showed a low average recombination rate. We constructed genetic linkage maps for each family and a composite linkage map based on recombinant inbred lines (RILs) across the families and identified and annotated 525,772 high confidence SNPs that were used to impute the SNP alleles in the RILs. The segregation distortion in most families significantly favored the alleles from the female parent, and there was no significant difference of residual heterozygosity in the euchromatic vs. heterochromatic regions. The genotypic datasets for the RILs and parents are publicly available and are anticipated to be useful to map quantitative trait loci (QTL) controlling important traits in soybean.

Mapping the low palmitate fap1 mutation and validation of its effects in soybean oil and agronomic traits in three soybean populations

KEY MESSAGE

fap 1 mutation is caused by a G174A change in GmKASIIIA that disrupts a donor splice site recognition and creates a GATCTG motif that enhanced its expression. Soybean oil with reduced palmitic acid content is desirable to reduce the health risks associated with consumption of this fatty acid. The objectives of this study were: to identify the genomic location of the reduced palmitate fap1 mutation, determine its molecular basis, estimate the amount of phenotypic variation in fatty acid composition explained by this locus, determine if there are epistatic interactions between the fap1 and fap nc loci and, determine if the fap1 mutation has pleiotropic effects on seed yield, oil and protein content in three soybean populations. This study detected two major QTL for 16:0 content located in chromosome 5 (GmFATB1a, fap nc) and chromosome 9 near BARCSOYSSR_09_1707 that explained, with their interaction, 66-94 % of the variation in 16:0 content in the three populations. Sequencing results of a putative candidate gene, GmKASIIIA, revealed a single unique polymorphism in the germplasm line C1726, which was predicted to disrupt the donor splice site recognition between exon one and intron one and produce a truncated KASIIIA protein. This G to A change also created the GATCTG motif that enhanced gene expression of the mutated GmKASIIIA gene. Lines homozygous for the GmKASIIIA mutation (fap1) had a significant reduction in 16:0, 18:0, and oil content; and an increase in unsaturated fatty acids content. There were significant epistatic interactions between GmKASIIIA (fap1) and fap nc for 16:0 and oil contents, and seed yield in two populations. In conclusion, the fap1 phenotype is caused by a single unique SNP in the GmKASIIIA gene.

Mapping soybean aphid resistance genes in PI 567598B

The soybean aphid (Aphis glycines Matsumura) has been a major pest of soybean [Glycine max (L.) Merr.] in North America since it was first reported in 2000. Our previous study revealed that the strong aphid resistance of plant introduction (PI) 567598B was controlled by two recessive genes. The objective of this study was to locate these two genes on the soybean genetic linkage map using molecular markers. A mapping population of 282 F4:5 lines derived from IA2070 × E06902 was evaluated for aphid resistance in a field trial in 2009 and a greenhouse trial in 2010. Two quantitative trait loci (QTLs) were identified using the composite and multiple interval mapping methods, and were mapped on chromosomes 7 (linkage group M) and 16 (linkage group J), respectively. E06902, a parent derived from PI 567598B, conferred resistance at both loci. In the 2010 greenhouse trial, each of the two QTLs explained over 30 % of the phenotypic variation. Significant epistatic interaction was also found between these two QTLs. However, in the 2009 field trial, only the QTL on chromosome 16 was found and it explained 56.1 % of the phenotypic variation. These two QTLs and their interaction were confirmed with another population consisting of 94 F2:5 lines in the 2008 and 2009 greenhouse trials. For both trials in the alternative population, these two loci explained about 50 and 80.4 % of the total phenotypic variation, respectively. Our study shows that soybean aphid isolate used in the 2009 field trial defeated the QTL found on chromosome 7. Presence of the QTL on chromosome 16 conferred soybean aphid resistance in all trials. The markers linked to the aphid-resistant QTLs in PI 567598B or its derived lines can be used in marker-assisted breeding for aphid resistance.

Conversion of AFLP bands into high-throughput DNA markers

The conversion of AFLP bands into polymorphic sequence-tagged-site (STS) markers is necessary for high-throughput genotype scoring. Technical hurdles that must be overcome arise from genome complexity (particularly sequence duplication), from the low-molecular-weight nature of the AFLP bands and from the location of the polymorphism within the AFLP band. We generated six STS markers from ten AFLP bands (four AFLPs were from co-dominant pairs of bands) in soybean (Glycine max). The markers were all linked to one of two loci, rhg1 on linkage group G and Rhg4 on linkage group A2, that confer resistance to the soybean cyst nematode (Heterodera glycines I.). When the polymorphic AFLP band sequence contained a duplicated sequence or could not be converted to a locus-specific STS marker, direct sequencing of BAC clones anchored to a physical map generated locus-specific flanking sequences at the polymorphic locus. When the polymorphism was adjacent to the restriction site used in the AFLP analysis, single primer extension was performed to reconstruct the polymorphism. The six converted AFLP markers represented 996 bp of sequence from alleles of each of two cultivars and identified eight insertions or deletions, two microsatellites and eight single-nucleotide polymorphisms (SNPs). The polymorphic sequences were used to design a non-electrophoretic, fluorometric assay (based on the TaqMan technology) and/or develop electrophoretic STS markers for high-throughput genotype determination during marker-assisted breeding for resistance to cyst nematode. We conclude that the converted AFLP markers contained polymorphism at a 10- to 20-fold higher frequency than expected for adapted soybean cultivars and that the efficiency of AFLP band conversion to STS can be improved using BAC libraries and physical maps. The method provides an efficient tool for SNP and STS discovery suitable for marker-assisted breeding and genomics.