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 F_4: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 F_2: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.     

Identification of positive yield QTL alleles from exotic soybean germplasm in two backcross populations

Increasing seed yield is an important breeding goal of soybean [Glycine max (L.) Merr.] improvement efforts. Due to the small number of ancestors and subsequent breeding and selection, the genetic base of current soybean cultivars in North America is narrow. The objective of this study was to map quantitative trait loci (QTL) in two backcross populations developed using soybean plant introductions as donor parents. The first population included 116 BC(2)F(3)-derived lines developed using "Elgin" as the recurrent parent and PI 436684 as the donor parent (E population). The second population included 93 BC(3)F(3)-derived lines developed with "Williams 82" as the recurrent parent and PI 90566-1 as the donor parent (W population). The two populations were evaluated with 1,536 SNP markers and during 2?years for seed yield and other agronomic traits. Genotypic and phenotypic data were analyzed using the programs MapQTL and QTLNetwork to identify major QTL and epistatic QTL. In the E population, two yield QTL were identified by both MapQTL and QTLNetwork, and the PI 436684 alleles were associated with yield increases. In the W population, a QTL allele from PI 90566-1 accounted for 30?% of the yield variation; however, the PI region was also associated with later maturity and shorter plant height. No epistasis for seed yield was identified in either population. No yield QTL was previously reported at the regions where these QTL map indicating that exotic germplasm can be a source of new alleles that can improve soybean yield.     

Molecular mapping of soybean rust resistance in soybean accession PI 561356 and SNP haplotype analysis of the Rpp1 region in diverse germplasm

Soybean rust (SBR), caused by Phakopsora pachyrhizi Sydow, is one of the most economically important and destructive diseases of soybean [Glycine max (L.) Merr.] and the discovery of novel SBR resistance genes is needed because of virulence diversity in the pathogen. The objectives of this research were to map SBR resistance in plant introduction (PI) 561356 and to identify single nucleotide polymorphism (SNP) haplotypes within the region on soybean chromosome 18 where the SBR resistance gene Rpp1 maps. One-hundred F(2:3) lines derived from a cross between PI 561356 and the susceptible experimental line LD02-4485 were genotyped with genetic markers and phenotyped for resistance to P. pachyrhizi isolate ZM01-1. The segregation ratio of reddish brown versus tan lesion type in the population supported that resistance was controlled by a single dominant gene. The gene was mapped to a 1-cM region on soybean chromosome 18 corresponding to the same interval as Rpp1. A haplotype analysis of diverse germplasm across a 213-kb interval that included Rpp1 revealed 21 distinct haplotypes of which 4 were present among 5 SBR resistance sources that have a resistance gene in the Rpp1 region. Four major North American soybean ancestors belong to the same SNP haplotype as PI 561356 and seven belong to the same haplotype as PI 594538A, the Rpp1-b source. There were no North American soybean ancestors belonging to the SNP haplotypes found in PI 200492, the source of Rpp1, or PI 587886 and PI 587880A, additional sources with SBR resistance mapping to the Rpp1 region.     

’Forrest’ resistance to the soybean cyst nematode is bigenic: saturation mapping of the Rhg1and Rhg4 loci

Field resistance to cyst nematode (SCN) race 3 (Heterodera glycines I.) in soybean [Glycine max (L.) Merr.] cv ’Forrest’ is conditioned by two QTLs: the underlying genes are presumed to include Rhg1 on linkage group G and Rhg4 on linkage group A2. A population of recombinant inbred lines (RILs) and two populations of near-isogenic lines (NILs) derived from a cross of Forrest×Essex were used to map the loci affecting resistance to SCN. Bulked segregant analysis, with 512 AFLP primer combinations and microsatellite markers, produced a high-density genetic map for the intervals carrying Rhg1 and Rhg4. The two QTLs involved in resistance to SCN were strongly associated with the AFLP marker E_ATGM_CGA87 (P=0.0001, R²=24.5%) on linkage group G, and the AFLP marker E_CCGM_AAC405 (P=0.0001, R² =26.2%) on linkage group A2. Two- way analysis of variance showed epistasic interaction (P=0.0001, R² =16%) between the two loci controlling SCN resistance in Essex×Forrest recombinant inbred lines. Considering the two loci as qualitative genes and the resistance as female index FI <5%, jointly the two loci explained over 98% of the resistance. The locations of the two QTLs were confirmed in the NILs populations. Therefore SCN resistance in Forrest×Essex is bigenic. High-efficiency marker-assisted selection can be performed using the markers to develop cultivars with stable resistance to SCN. Received: 5 November 2000 / Accepted: 23 January 2001     

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.

Abstract

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.     

Two large-insert soybean genomic libraries constructed in a binary vector: applications in chromosome walking and genome wide physical mapping

Large DNA insert libraries in binary T-DNA vectors can assist in the isolation of the gene(s) underlying a quantitative trait locus (QTL). Binary vectors facilitate the transfer of large-insert DNA fragments containing a QTL from E. coli to Agrobacterium sp. and then to plants. We constructed two soybean large-insert libraries from cv. Forrest in the pCLD04541 (V41) binary vector after partial digestion of genomic high-molecular-weight DNA with BamHI or HindIII. The libraries contain 76,800 clones with an average insert size of 125 kb, and therefore represent 9.5-fold haploid genome equivalents. Colony hybridization using a chloroplast-specific probe infers that the libraries contain less than 0.5% clones of chloroplast DNA origin. These two libraries have provided clones for physical mapping of the soybean genome and for isolation of a number of disease resistance genes. One microsatellite marker was identified from the clone that hybridized to the Bng122 RFLP probe. The sequence-tagged site was used for genetic mapping and marker-assisted selection for genes underlying resistance to the soybean cyst nematode and sudden death syndrome. Received: 7 May 1999 / Accepted: 18 February 2000     

Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

Structural variants in the soybean genome localize to clusters of biotic stress-response genes

Genome-wide structural and gene content variations are hypothesized to drive important phenotypic variation within a species. Structural and gene content variations were assessed among four soybean (Glycine max) genotypes using array hybridization and targeted resequencing. Many chromosomes exhibited relatively low rates of structural variation (SV) among genotypes. However, several regions exhibited both copy number and presence-absence variation, the most prominent found on chromosomes 3, 6, 7, 16, and 18. Interestingly, the regions most enriched for SV were specifically localized to gene-rich regions that harbor clustered multigene families. The most abundant classes of gene families associated with these regions were the nucleotide-binding and receptor-like protein classes, both of which are important for plant biotic defense. The colocalization of SV with plant defense response signal transduction pathways provides insight into the mechanisms of soybean resistance gene evolution and may inform the development of new approaches to resistance gene cloning.     

Fine mapping the soybean aphid resistance gene Rag1 in soybean

The soybean aphid (Aphis glycines Matsumura) is an important soybean [Glycine max (L.) Merr.] pest in North America. The dominant aphid resistance gene Rag1 was previously mapped from the cultivar ‘Dowling’ to a 12 cM marker interval on soybean chromosome 7 (formerly linkage group M). The development of additional genetic markers mapping closer to Rag1 was needed to accurately position the gene to improve the effectiveness of marker-assisted selection (MAS) and to eventually clone it. The objectives of this study were to identify single nucleotide polymorphisms (SNPs) near Rag1 and to position these SNPs relative to Rag1. To generate a fine map of the Rag1 interval, 824 BC₄F₂ and 1,000 BC₄F₃ plants segregating for the gene were screened with markers flanking Rag1. Plants with recombination events close to the gene were tested with SNPs identified in previous studies along with new SNPs identified from the preliminary Williams 82 draft soybean genome shotgun sequence using direct re-sequencing and gene-scanning melt-curve analysis. Progeny of these recombinant plants were evaluated for aphid resistance. These efforts resulted in the mapping of Rag1 between the two SNP markers 46169.7 and 21A, which corresponds to a physical distance on the Williams 82 8× draft assembly (Glyma1.01) of 115 kilobase pair (kb). Several candidate genes for Rag1 are present within the 115-kb interval. The markers identified in this study that are closely linked to Rag1 will be a useful resource in MAS for this important aphid resistance gene.     

SNP-PHAGE – High throughput SNP discovery pipeline

Background  

Single nucleotide polymorphisms (SNPs) as defined here are single base sequence changes or short insertion/deletions between or within individuals of a given species. As a result of their abundance and the availability of high throughput analysis technologies SNP markers have begun to replace other traditional markers such as restriction fragment length polymorphisms (RFLPs), amplified fragment length polymorphisms (AFLPs) and simple sequence repeats (SSRs or microsatellite) markers for fine mapping and association studies in several species. For SNP discovery from chromatogram data, several bioinformatics programs have to be combined to generate an analysis pipeline. Results have to be stored in a relational database to facilitate interrogation through queries or to generate data for further analyses such as determination of linkage disequilibrium and identification of common haplotypes. Although these tasks are routinely performed by several groups, an integrated open source SNP discovery pipeline that can be easily adapted by new groups interested in SNP marker development is currently unavailable.