Roxanne Denny

Distribution of Microsatellites in the Genome of Medicago truncatula: A Resource of Genetic Markers That Integrate Genetic and Physical Maps

Microsatellites are tandemly repeated short DNA sequences that are favored as molecular-genetic markers due to their high polymorphism index. Plant genomes characterized to date exhibit taxon-specific differences in frequency, genomic location, and motif structure of microsatellites, indicating that extant microsatellites originated recently and turn over quickly. With the goal of using microsatellite markers to integrate the physical and genetic maps of Medicago truncatula, we surveyed the frequency and distribution of perfect microsatellites in 77 Mbp of gene-rich BAC sequences, 27 Mbp of nonredundant transcript sequences, 20 Mbp of random whole genome shotgun sequences, and 49 Mbp of BAC-end sequences. Microsatellites are predominantly located in gene-rich regions of the genome, with a density of one long (i.e., ≥20 nt) microsatellite every 12 kbp, while the frequency of individual motifs varied according to the genome fraction under analysis. A total of 1,236 microsatellites were analyzed for polymorphism between parents of our reference intraspecific mapping population, revealing that motifs (AT)n, (AG)n, (AC)n, and (AAT)n exhibit the highest allelic diversity. A total of 378 genetic markers could be integrated with sequenced BAC clones, anchoring 274 physical contigs that represent 174 Mbp of the genome and composing an estimated 70% of the euchromatic gene space.

Two simple sequence repeat markers to select for soybean cyst nematode resistance coditioned by the rhg1 locus

Abstract The soybean cyst nematode (SCN) (Heterodera glycines Inchinoe) is the most economically significant soybean pest. The principal strategy to reduce or eliminate damage from this pest is the use of resistant cultivars. Identifying resistant segregants in a breeding program is a difficult and expensive process which is complicated by the oligogenic nature of the resistance and genetic variability in the pathogen. Fortunately, resistance at one SCN-resistance locus, rhg1, is generally accepted as a necessity for the development of resistant genotypes using any source of resistance and when challenged by any SCN race. Thus, the development of SCN resistant cultivars would be expedited if an effective and rapid system were available to identify breeding lines carrying a resistance allele at the rhg1 locus. In this study we report two simple sequence repeat (SSR) or microsatellite loci that cosegregate and map 0.4 cM from rhg1. Allelic variation at the first of these loci, BARC-Satt309, distinguished most, if not all, SCN-susceptible genotypes from those carrying resistance at rhg1 derived from the important SCN-resistance sources ’Peking’, PI 437654, and PI 90763. BARC-Satt309 was also effective in distinguishing SCN resistance sources PI 88788 and PI 209332 from many, but not all, susceptible genotypes. BARC-Satt309 cannot be used in marker-assisted selection in populations developed from typical southern US cultivars crossed with the important resistance sources PI 88788 or PI 209332 because these genotypes all carry the identical allele at the BARC-Satt309 locus. A second SSR locus, BARC-Sat_168, was developed from a bacterial artificial chromosome (BAC) clone that was identified using the primers to BARC-Satt309. BARC-Sat_168 distinguished PI 88788 and PI 209332 from southern US cultivars such as ’Lee’, ’Bragg’ and ’Essex’. Both BARC-Satt309 and BARC-Sat_168 were used to assay lines from SCN-susceptible×SCN-resistant crosses and proved to be highly effective in identifying lines carrying rhg1 resistance from those carrying the allele for SCN susceptibility at the rhg1 locus.

Whole-genome nucleotide diversity, recombination, and linkage disequilibrium in the model legume Medicago truncatula

Medicago truncatula is a model for investigating legume genetics, including the genetics and evolution of legume–rhizobia symbiosis. We used whole-genome sequence data to identify and characterize sequence polymorphisms and linkage disequilibrium (LD) in a diverse collection of 26 M. truncatula accessions. Our analyses reveal that M. truncatula harbors both higher diversity and less LD than soybean (Glycine max) and exhibits patterns of LD and recombination similar to Arabidopsis thaliana. The population-scaled recombination rate is approximately one-third of the mutation rate, consistent with expectations for a species with a high selfing rate. Linkage disequilibrium, however, is not extensive, and therefore, the low recombination rate is likely not a major constraint to adaptation. Nucleotide diversity in 100-kb windows was negatively correlated with gene density, which is expected if diversity is shaped by selection acting against slightly deleterious mutations. Among putative coding regions, members of four gene families harbor significantly higher diversity than the genome-wide average. Three of these families are involved in resistance against pathogens; one of these families, the nodule-specific, cysteine-rich gene family, is specific to the galegoid legumes and is involved in control of rhizobial differentiation. The more than 3 million SNPs that we detected, approximately one-half of which are present in more than one accession, are a valuable resource for genome-wide association mapping of genes responsible for phenotypic diversity in legumes, especially traits associated with symbiosis and nodulation.

Targeted isolation of simple sequence repeat markers through the use of bacterial artificial chromosomes

Abstract Simple sequence repeats (SSRs) are versatile DNA markers that are readily assayed and highly informative. Unfortunately, non-targeted approaches to SSR development often leave large genomic regions without SSR markers. In some cases these same genomic regions are already populated by other types of DNA markers, especially restriction fragment length polymorphisms (RFLPs), random amplified polymorphic DNAs (RAPDs), and amplified fragment length polymorphisms (AFLPs). To identify SSR markers in such regions, bacterial artificial chromosome (BAC) clones can be used as intermediaries. First, one or more BAC clones in a region of interest are identified through the use of an existing DNA marker. BAC clones uncovered in this initial step are then used to create a small insert DNA library that can be screened for the presence of SSR-containing clones. Because BAC inserts are often 100-kb pairs or more in size, most contain one or more SSRs. This strategy was applied to two regions of the soybean genome near genes that condition resistance to the soybean cyst nematode on molecular linkage groups G and A2. This targeted approach to identifying new DNA markers can readily be extended to other types of DNA markers, including single nucleotide polymorphisms.

Differential Accumulation of Retroelements and Diversification of NB-LRR Disease Resistance Genes in Duplicated Regions following Polyploidy in the Ancestor of Soybean

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

A bacterial artificial chromosome library for soybean and identification of clones near a major cyst nematode resistance gene

Abstract We constructed a bacterial artificial chromosome (BAC) library for soybean (Glycine max) consisting of approximately 30 000 clones with an average insert size of 120 kilobase pairs. The library was successfully screened with restriction fragment length polymorphism (RFLP) and microsatellite markers tightly linked to a major resistance gene for the cyst nematode, Heterodera glycines. Since many soybean RFLPs hybridize to duplicate loci, BACs homologous to duplicate RFLP loci were distinguished by digestion with the restriction enzyme originally used to map the RFLP, followed by a comparison of the hybridizing fragments. Linkage mapping of BAC clones identified with markers linked to the cyst nematode resistance gene demonstrated that these clones were located at the expected chromosomal positions and that there were no indications of chimeras within the genomic inserts.

Candidate genes and genetic architecture of symbiotic and agronomic traits revealed by whole-genome, sequence-based association genetics in Medicago truncatula.

Genome-wide association study (GWAS) has revolutionized the search for the genetic basis of complex traits. To date, GWAS have generally relied on relatively sparse sampling of nucleotide diversity, which is likely to bias results by preferentially sampling high-frequency SNPs not in complete linkage disequilibrium (LD) with causative SNPs. To avoid these limitations we conducted GWAS with >6 million SNPs identified by sequencing the genomes of 226 accessions of the model legume Medicago truncatula. We used these data to identify candidate genes and the genetic architecture underlying phenotypic variation in plant height, trichome density, flowering time, and nodulation. The characteristics of candidate SNPs differed among traits, with candidates for flowering time and trichome density in distinct clusters of high linkage disequilibrium (LD) and the minor allele frequencies (MAF) of candidates underlying variation in flowering time and height significantly greater than MAF of candidates underlying variation in other traits. Candidate SNPs tagged several characterized genes including nodulation related genes SERK2, MtnodGRP3, MtMMPL1, NFP, CaML3, MtnodGRP3A and flowering time gene MtFD as well as uncharacterized genes that become candidates for further molecular characterization. By comparing sequence-based candidates to candidates identified by in silico 250K SNP arrays, we provide an empirical example of how reliance on even high-density reduced representation genomic makers can bias GWAS results. Depending on the trait, only 30–70% of the top 20 in silico array candidates were within 1 kb of sequence-based candidates. Moreover, the sequence-based candidates tagged by array candidates were heavily biased towards common variants; these comparisons underscore the need for caution when interpreting results from GWAS conducted with sparsely covered genomes.

Differential expression of two soybean apyrases, one of which is an early nodulin

Two cDNA clones were isolated from soybean (Glycine soja) by polymerase chain reaction with primers designed to conserved motifs found in apyrases (nucleotide phosphohydrolase). The two cDNAs are predicted to encode for two, distinct, apyrase proteins of approximately 50 kDa (i.e., GS50) and 52 kDa (i.e., GS52). Phylogenetic analysis indicated that GS52 is orthologous to a family of apyrases recently suggested to play a role in legume nodulation. GS50 is paralogous to this family and, therefore, likely plays a different physiological role. Consistent with this analysis, GS50 mRNA was detected in root, hypocotyls, flowers, and stems, while GS52 mRNA was found in root and flowers. Neither gene was expressed in leaves or cotyledons. Inoculation of roots with Bradyrhizobium japonicum, nitrogen-fixing symbiont of soybean, resulted in the rapid (<6 h) induction of GS52 mRNA expression. The level of GS50 mRNA expression was not affected by bacterial inoculation. Western blot (immunoblot) analysis of GS50 expression mirrored the results obtained by mRNA analysis. However, in contrast to the mRNA results, GS52 protein was found in stems. Interestingly, anti-GS52 antibody recognized a 50-kDa protein found only in nodule extracts. Treatment of roots with anti-GS52 antibody, but not anti-GS50 antibody or preimmune serum, blocked nodulation by B. japonicum. Fractionation of cellular membranes in sucrose density gradients and subsequent Western analysis of the fractions revealed that GS50 colocalized with marker enzymes for the Golgi, while GS52 colocalized with marker enzymes for the plasma membrane. Restriction fragment length polymorphism (RFLP)-based mapping placed the gs52 gene on major linkage group J of the integrated genetic map of soybean. These data suggest that GS50 is likely an endo-apyrase involved in Golgi function, while GS52 is localized on the root surface and appears to play an important role in nodulation.

Comparative genomics of the core and accessory genomes of 48 Sinorhizobium strains comprising five genospecies

BackgroundThe sinorhizobia are amongst the most well studied members of nitrogen-fixing root nodule bacteria and contribute substantial amounts of fixed nitrogen to the biosphere. While the alfalfa symbiont Sinorhizobium meliloti RM 1021 was one of the first rhizobial strains to be completely sequenced, little information is available about the genomes of this large and diverse species group.ResultsHere we report the draft assembly and annotation of 48 strains of Sinorhizobium comprising five genospecies. While S. meliloti and S. medicae are taxonomically related, they displayed different nodulation patterns on diverse Medicago host plants, and have differences in gene content, including those involved in conjugation and organic sulfur utilization. Genes involved in Nod factor and polysaccharide biosynthesis, denitrification and type III, IV, and VI secretion systems also vary within and between species. Symbiotic phenotyping and mutational analyses indicated that some type IV secretion genes are symbiosis-related and involved in nitrogen fixation efficiency. Moreover, there is a correlation between the presence of type IV secretion systems, heme biosynthesis and microaerobic denitrification genes, and symbiotic efficiency.ConclusionsOur results suggest that each Sinorhizobium strain uses a slightly different strategy to obtain maximum compatibility with a host plant. This large genome data set provides useful information to better understand the functional features of five Sinorhizobium species, especially compatibility in legume-Sinorhizobium interactions. The diversity of genes present in the accessory genomes of members of this genus indicates that each bacterium has adopted slightly different strategies to interact with diverse plant genera and soil environments.

Targeted comparative genome analysis and qualitative mapping of a major partial-resistance gene to the soybean cyst nematode.

A major partial-resistance locus to the soybean cyst nematode (Heterodera glycines Ichinohe; SCN) was identified on linkage group ‘G’ of soybean [Glycine max (L.) Merr.] using restriction fragment length polymorphisms (RFLPs). This locus explained 51.4% (LOD=10.35) of the total phenotypic variation in disease response in soybean Plant Introduction (PI) 209332, 52.7% (LOD=15.58) in PI 90763, 40.0% (LOD=10.50) in PI 88788, and 28.1% (LOD=6.94) in ‘Peking’. Initially, the region around this major resistance locus was poorly populated with DNA markers. To increase marker density in this genomic region, first random, and later targeted, comparative mapping with RFLPs from mungbean [Vigna radiata (L.) R. Wilcz.] and common bean (Phaseolus vulgaris L.) was performed, eventually leading to one RFLP marker every 2.6 centimorgans (cM). Even with this marker density, the inability to resolve SCN disease response into discrete Mendelian categories posed a major limitation to mapping. Thus, qualitative scoring of SCN disease response was carried out in an F5∶6 recombinant inbred population derived from ‘Evans’xPI 209332 using a 30% disease index cut-off for resistance. Using the computer program JoinMap, an integrated map of the region of interest was created, placing the SCN resistance locus 4.6 cM from RFLP marker B53 and 2.8 cM from Bng30. This study demonstrates how a combination of molecularmapping strategies, including comparative genome analysis, join mapping, and qualitative scoring of a quantitative trait, potentially provide the necessary tools for high-resolution mapping around a quantitative-trait locus.