A copia-like retrotransposon Tgmr closely linked to the Rps1-k allele that confers race-specific resistance of soybean to Phytophthora sojae

The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene

Soybean growers widely use the Resistance to Heterodera glycines 1 (Rhg1) locus to reduce yield losses caused by soybean cyst nematode (SCN). Rhg1 is a tandemly repeated four gene block. Two classes of SCN resistance-conferring Rhg1 haplotypes are recognized: rhg1-a ("Peking-type," low-copy number, three or fewer Rhg1 repeats) and rhg1-b ("PI 88788-type," high-copy number, four or more Rhg1 repeats). The rhg1-a and rhg1-b haplotypes encode α-SNAP (alpha-Soluble NSF Attachment Protein) variants α-SNAP Rhg1 LC and α-SNAP Rhg1 HC, respectively, with differing atypical C-terminal domains, that contribute to SCN resistance. Here we report that rhg1-a soybean accessions harbor a copia retrotransposon within their Rhg1 Glyma.18G022500 (α-SNAP-encoding) gene. We termed this retrotransposon "RAC," for Rhg1 alpha-SNAP copia. Soybean carries multiple RAC-like retrotransposon sequences. The Rhg1 RAC insertion is in the Glyma.18G022500 genes of all true rhg1-a haplotypes we tested and was not detected in any examined rhg1-b or Rhg1WT (single-copy) soybeans. RAC is an intact element residing within intron 1, anti-sense to the rhg1-a α-SNAP open reading frame. RAC has intrinsic promoter activities, but overt impacts of RAC on transgenic α-SNAP Rhg1 LC mRNA and protein abundance were not detected. From the native rhg1-a RAC+ genomic context, elevated α-SNAP Rhg1 LC protein abundance was observed in syncytium cells, as was previously observed for α-SNAP Rhg1 HC (whose rhg1-b does not carry RAC). Using a SoySNP50K SNP corresponding with RAC presence, just ~42% of USDA accessions bearing previously identified rhg1-a SoySNP50K SNP signatures harbor the RAC insertion. Subsequent analysis of several of these putative rhg1-a accessions lacking RAC revealed that none encoded α-SNAPRhg1LC, and thus, they are not rhg1-a. rhg1-a haplotypes are of rising interest, with Rhg4, for combating SCN populations that exhibit increased virulence against the widely used rhg1-b resistance. The present study reveals another unexpected structural feature of many Rhg1 loci, and a selectable feature that is predictive of rhg1-a haplotypes.

A Novel Phytophthora sojae Resistance Rps12 Gene Mapped to a Genomic Region That Contains Several Rps Genes

Phytophthora sojae Kaufmann and Gerdemann, which causes Phytophthora root rot, is a widespread pathogen that limits soybean production worldwide. Development of Phytophthora resistant cultivars carrying Phytophthora resistance Rps genes is a cost-effective approach in controlling this disease. For this mapping study of a novel Rps gene, 290 recombinant inbred lines (RILs) (F7 families) were developed by crossing the P. sojae resistant cultivar PI399036 with the P. sojae susceptible AR2 line, and were phenotyped for responses to a mixture of three P. sojae isolates that overcome most of the known Rps genes. Of these 290 RILs, 130 were homozygous resistant, 12 heterzygous and segregating for Phytophthora resistance, and 148 were recessive homozygous and susceptible. From this population, 59 RILs homozygous for Phytophthora sojae resistance and 61 susceptible to a mixture of P. sojae isolates R17 and Val12-11 or P7074 that overcome resistance encoded by known Rps genes mapped to Chromosome 18 were selected for mapping novel Rps gene. A single gene accounted for the 1:1 segregation of resistance and susceptibility among the RILs. The gene encoding the Phytophthora resistance mapped to a 5.8 cM interval between the SSR markers BARCSOYSSR_18_1840 and Sat_064 located in the lower arm of Chromosome 18. The gene is mapped 2.2 cM proximal to the NBSRps4/6-like sequence that was reported to co-segregate with the Phytophthora resistance genes Rps4 and Rps6. The gene is mapped to a highly recombinogenic, gene-rich genomic region carrying several nucleotide binding site-leucine rich repeat (NBS-LRR)-like genes. We named this novel gene as Rps12, which is expected to be an invaluable resource in breeding soybeans for Phytophthora resistance.

Diversity and evolution of Ty1-copia retroelements within Chalcidoidea by reverse transcriptase domain analysis

Ty1-copia retrotransposons are widespread and diverse in insects. Some features of their hosts, such as mating and genetic systems, are predicted to influence the spread of selfish genetic elements like Ty1-copia. Using part of the reverse transcriptase gene as a reference, we experimentally surveyed Ty1-copia elements in eight species of fig wasps (Hymenoptera: Chalcidoidea), and performed an in silico analysis of six available genomes of chalcid wasps. Contrary to initial expectations that selfish elements such as Ty1-copia would be purged from the genomes of these species because of inbreeding and haplodiploidy, almost all of these wasps harbour an abundance of diverse Ty1-copia elements. Phylogenetic analyses suggest that the families of Ty1-copia elements found in these species have had a long association with their chalcid hosts. These results suggest an evolutionary scenario in which there was ancestral polymorphism followed by some taxa-specific events including stochastic loss and further diversification. Furthermore, estimating natural selection within the internal and terminal portions of the Ty1-copia phylogenies demonstrated that the elements are under strong evolutionary constraints for their long-term survival, but evolve like pseudogenes in the short term, accompanied by the rise and fall of parasitic elements in the history of wasp lineage.

Isolation and characterization of novel Ty1-copia-like retrotransposons from lily

Species of the genus Lilium are well known for their large genomes. Although expansion of noncoding repeated DNA is believed to account for this genome size, retroelement del Ty3-gypsy is the only one described so far in the genus Lilium. We isolated Ty1-copia elements from Lilium longiflorum and named them LIREs (lily retrotransposons). The long terminal repeats, primer binding site, and polypurine tract sequences are highly similar among the LIRE elements, indicating that they are in the same lineage. Although the protein-coding regions were highly decayed, the sequence motifs of the integrase, reverse transcriptase, and RNase H domains were identifiable as belonging to the order of Ty1-copia elements. Phylogenetic analysis and primer binding site sequences revealed that these elements belonged to the Ale lineage among the six lineages of plant Ty1-copia elements. Base substitutions in the long terminal repeats estimated that the integration times of the LIRE Ty1-copia elements were between 0.7 and 5.5 mya. In situ hybridization showed that the LIRE elements were present in all the chromosomes of L. longiflorum and L. lancifolium, but absent in centromeres, telomeres, and 45S rRNA sites in both species. The LIRE elements were present very abundantly in species of the genus Lilium, but absent in other genera of the family Liliaceae, implying that the LIRE elements might have contributed to the expansion of the genome in the genus Lilium.

Adaptive evolution involving gene duplication and insertion of a novel Ty1/copia-like retrotransposon in soybean

Gene duplication is a major force for generating evolutionary novelties that lead to adaptations to environments. We previously identified two paralogs encoding phytochrome A (phyA), GmphyA1 and GmphyA2, in soybean, a paleopolyploid species. GmphyA2 is encoded by the E4 locus responsible for photoperiod sensitivity. In photoperiod insensitive lines, GmphyA2 is inactivated by the insertion of a retrotransposon in exon 1. Here, we describe the detailed characterization of the element and its evolutionary significance inferred from the distribution of the allele that harbors the element. Structural characteristics indicated that the element, designated SORE-1, is a novel Ty1/copia-like retrotransposon in soybean, which was phylogenetically related to the Sto-4, BARE-1, and RIRE1 elements. The element was transcriptionally active, and the transcription was partially repressed by an epigenetic mechanism. Sequences homologous with SORE-1 were detected in a genome sequence database of soybean, most of which appeared silent. GmphyA2 that harbors the SORE-1 insertion was detected only in cultivated soybean lines grown in northern regions of Japan, consistent with the notion that photoperiod insensitivity caused by the dysfunction of GmphyA2 is one of genetic changes that allowed soybean cultivation at high latitudes. Taking into account that genetic redundancy is conferred by the two phyA genes, we propose a novel model for the consequences of gene duplication and transposition of retrotransposons: when the gene is duplicated, retrotransposon insertion that causes the loss of a gene function can lead to adaptive evolution while the organism is sustained by the buffering effect brought about by gene duplication.

Identification and occurrence of the LTR-Copia-like retrotransposon, PSCR and other Copia-like elements in the genome of Phytophthora sojae

Sequence analysis of the genomic region of Phytophthora sojae close to the Avr4/6 locus specifying virulence on soybean identified a Ty1/Copia-like retrotransposon that we have named Phytophthora sojae Copia-like retrotransposon (PSCR). Twelve near-complete homologs of PSCR were found in the published P. sojae genome sequence, none of which encoded a full-length polyprotein characteristic of Copia-like retrotransposons, or appears to exhibit transcriptional activity or show evidence of recent movement, suggesting they are non-functional and unlikely to have caused pathogenic variability. However, reconstructed consensus PSCR sequence encoding a full-length polyprotein resembles a functional, ancestral retroelement within P. sojae. Homologs were also found in sequence databases of other Phytophthora species. Database searches found other families of Copia-like elements in genomes of P. sojae, P. ramorum and P. infestans that were different from members of the PSCR family and from Copia-like elements reported in other organisms. It is possible that the various families of Copia-like retroelements identified in this study represent introgressions into the genome of ancient ancestor(s) of current Phytophthora species, where they have evolved and diverged considerably during the speciation. Some Copia-like families are transcriptionally active with the potential to transpose and contribute to pathogenic variation in current populations of P. sojae.

The marker SCK13(603) associated with resistance to ascochyta blight in chickpea is located in a region of a putative retrotransposon

The sequence characterized amplified region (SCAR) marker SCK13(603), associated with ascochyta blight resistance in a chickpea recombinant inbred line (RIL) population, was used as anchored sequence for genome walking. The PCRs performed in the walking steps to walk in the same direction produced eight bands in 5' direction and five bands in 3' direction with a length ranking from 530 to 2,871 bp. The assembly of the bands sequences along with the sequence of SCK13(603) resulted in 7,815 bp contig. Blastn analyses showed stretches of DNA sequence mainly distributed from the nucleotides 1,500 to 4,500 significantly similar to Medicago truncatula genomic DNA. Three open reading frames (ORFs) were identified and blastp analysis of predicted amino acids sequences revealed that ORF1, ORF2 and ORF3 had significant similarity to a CCHC zinc finger protein, to an integrase, and to a precursor of the glucoamylase s1/s2, respectively, from M. truncatula. The high homology of the putative proteins derived from ORF1 and ORF2 with retrotransposon proteins and the prediction of the existence of conserved domains usually present in retrotransposon proteins indicate that the marker SCK13(603) is located in a region of a putative retrotransposon. The information generated in this study has contributed to increase the knowledge of this important region for blight resistance in chickpea.

The soybean-Phytophthora resistance locus Rps1-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences

BACKGROUND

A series of Rps (resistance to Pytophthora sojae) genes have been protecting soybean from the root and stem rot disease caused by the Oomycete pathogen, Phytophthora sojae. Five Rps genes were mapped to the Rps1 locus located near the 28 cM map position on molecular linkage group N of the composite genetic soybean map. Among these five genes, Rps1-k was introgressed from the cultivar, Kingwa. Rps1-k has been providing stable and broad-spectrum Phytophthora resistance in the major soybean-producing regions of the United States. Rps1-k has been mapped and isolated. More than one functional Rps1-k gene was identified from the Rps1-k locus. The clustering feature at the Rps1-k locus might have facilitated the expansion of Rps1-k gene numbers and the generation of new recognition specificities. The Rps1-k region was sequenced to understand the possible evolutionary steps that shaped the generation of Phytophthora resistance genes in soybean.

RESULTS

Here the analyses of sequences of three overlapping BAC clones containing the 184,111 bp Rps1-k region are reported. A shotgun sequencing strategy was applied in sequencing the BAC contig. Sequence analysis predicted a few full-length genes including two Rps1-k genes, Rps1-k-1 and Rps1-k-2. Previously reported Rps1-k-3 from this genomic region 1 was evolved through intramolecular recombination between Rps1-k-1 and Rps1-k-2 in Escherichia coli. The majority of the predicted genes are truncated and therefore most likely they are nonfunctional. A member of a highly abundant retroelement, SIRE1, was identified from the Rps1-k region. The Rps1-k region is primarily composed of repetitive sequences. Sixteen simple repeat and 63 tandem repeat sequences were identified from the locus.

CONCLUSION

These data indicate that the Rps1 locus is located in a gene-poor region. The abundance of repetitive sequences in the Rps1-k region suggested that the location of this locus is in or near a heterochromatic region. Poor recombination frequencies combined with presence of two functional Rps genes at this locus has been providing stable Phytophthora resistance in soybean.

LTR retrotransposons and flowering plant genome size: emergence of the increase/decrease model

Long Terminal Repeat (LTR) retrotransposons are ubiquitous components of plant genomes. Because of their copy-and-paste mode of transposition, these elements tend to increase their copy number while they are active. In addition, it is now well established that the differences in genome size observed in the plant kingdom are accompanied by variations in LTR retrotransposon content, suggesting that LTR retrotransposons might be important players in the evolution of plant genome size, along with polyploidy. The recent availability of large genomic sequences for many crop species has made it possible to examine in detail how LTR retrotransposons actually drive genomic changes in plants. In the present paper, we provide a review of the recent publications that have contributed to the knowledge of plant LTR retrotransposons, as structural components of the genomes, as well as from an evolutionary genomic perspective. These studies have shown that plant genomes undergo genome size increases through bursts of retrotransposition, while there is a counteracting process that tends to eliminate the transposed copies from the genomes. This process involves recombination mechanisms that occur either between the LTRs of the elements, leading to the formation of solo-LTRs, or between direct repeats anywhere in the sequence of the element, leading to internal deletions. All these studies have led to the emergence of a new model for plant genome evolution that takes into account both genome size increases (through retrotransposition) and decreases (through solo-LTR and deletion formation). In the conclusion, we discuss this new model and present the future prospects in the study of plant genome evolution in relation to the activity of transposable elements.

Soybean phytophthora resistance gene Rps8 maps closely to the Rps3 region

Root and stem rot is one of the major diseases of soybean. It is caused by the oomycete pathogen Phytophthora sojae. A series of resistance genes (Rps) have been providing soybean with reasonable protection against this pathogen. Among these genes, Rps8, which confers resistance to most P. sojae isolates, recently has been mapped. However, the most closely linked molecular marker was mapped at about 10 cM from Rps8. In this investigation, we attempted to develop a high-density genetic map of the Rps8 region and identify closely linked SSR markers for marker-assisted selection of this invaluable gene. Bulk segregant analysis was conducted for the identification of SSR markers that are tightly linked to Rps8. Polymorphic SSR markers selected from the Rps8 region failed to show cosegregation with Phytophthora resistance. Subsequently, bulk segregant analysis of the whole soybean genome and mapping experiments revealed that the Rps8 gene maps closely to the disease resistance gene-rich Rps3 region.

Identification of a large cluster of coiled coil-nucleotide binding site--leucine rich repeat-type genes from the Rps1 region containing Phytophthora resistance genes in soybean

Fifteen Rps genes confer resistance against the oomycete pathogen Phytophthora sojae, which causes root and stem rot disease in soybean. We have isolated a disease resistance gene-like sequence from the genomic region containing Rps1-k. Four classes of cDNA of the sequence were isolated from etiolated hypocotyl tissues that express the Rps1-k-encoded Phytophthora resistance. Sequence analyses of a cDNA clone showed that the sequence is a member of the coiled coil-nucleotide binding site-leucine rich repeat (CC-NBS-LRR)-type of disease resistance genes. It showed 36% identity to the recently cloned soybean resistance gene Rpg1-b, which confers resistance against Pseudomonas syringae pv. glycinea, and 56% and 38% sequence identity to putative resistance gene sequences from lotus and Medicago truncatula, respectively. The soybean genome contains about 38 copies of the sequence. Most of these copies are clustered in approximately 600 kb of contiguous DNA of the Rps1-k region. We have identified a recombinant that carries both rps1-k- and Rps1-k-haplotype-specific allelomorphs of two Rps1-k-linked molecular markers. An unequal crossover event presumably led to duplication of alleles for these two physically linked molecular markers. We hypothesize that the unequal crossing over was one of the mechanisms involved in tandem duplication of CC-NBS-LRR sequences in the Rps1-k region.

Diaspora, a large family of Ty3-gypsy retrotransposons in Glycine max, is an envelope-less member of an endogenous plant retrovirus lineage

BACKGROUND

The chromosomes of higher plants are littered with retrotransposons that, in many cases, constitute as much as 80% of plant genomes. Long terminal repeat retrotransposons have been especially successful colonizers of the chromosomes of higher plants and examinations of their function, evolution, and dispersal are essential to understanding the evolution of eukaryotic genomes. In soybean, several families of retrotransposons have been identified, including at least two that, by virtue of the presence of an envelope-like gene, may constitute endogenous retroviruses. However, most elements are highly degenerate and are often sequestered in regions of the genome that sequencing projects initially shun. In addition, finding potentially functional copies from genomic DNA is rare. This study provides a mechanism to surmount these issues to generate a consensus sequence that can then be functionally and phylogenetically evaluated.

RESULTS

Diaspora is a multicopy member of the Ty3-gypsy-like family of LTR retrotransposons and comprises at least 0.5% of the soybean genome. Although the Diaspora family is highly degenerate, and with the exception of this report, is not represented in the Genbank nr database, a full-length consensus sequence was generated from short overlapping sequences using a combination of experimental and in silico methods. Diaspora is 11,737 bp in length and contains a single 1892-codon ORF that encodes a gag-pol polyprotein. Phylogenetic analysis indicates that it is closely related to Athila and Calypso retroelements from Arabidopsis and soybean, respectively. These in turn form the framework of an endogenous retrovirus lineage whose members possess an envelope-like gene. Diaspora appears to lack any trace of this coding region.

CONCLUSION

A combination of empirical sequencing and retrieval of unannotated Genome Survey Sequence database entries was successfully used to construct a full-length representative of the Diaspora family in Glycine max. Diaspora is presently the only fully characterized member of a lineage of putative plant endogenous retroviruses that contains virtually no trace of an extra coding region. The loss of an envelope-like coding domain suggests that non-infectious retrotransposons could swiftly evolve from infectious retroviruses, possibly by anomalous splicing of genomic RNA.

The matrix metalloproteinase gene GmMMP2 is activated in response to pathogenic infections in soybean

Matrix metalloproteinases (MMPs) play an important role in host defense responses against pathogens in mammals where their activities lead to the production of antimicrobial peptides. We have identified a novel soybean (Glycine max) metalloproteinase gene, GmMMP2, that is transcriptionally up-regulated in infected tissues. The deduced amino acid sequence indicates that this gene belongs to the MMP family. It is a preproprotein containing an N-terminal signal peptide, a cysteine switch, a zinc-binding catalytic motif, and a C-terminal transmembrane domain. The GmMMP2 expressed in and purified from Escherichia coli exhibited an in vitro enzymatic activity in digesting myelin basic protein. All plant metalloproteinases reported so far have no known functions. However, they have been suggested to be involved in extracellular cell matrix degradation during development or senescence. Our investigations demonstrate that the GmMMP2 transcript levels were rapidly increased in compatible and incompatible interactions of soybean tissues with the oomycete pathogen Phytophthora sojae or the bacterial pathogen Pseudomonas syringae pv. glycinea. In agreement with the GmMMP2 activation, a metalloproteinase activity was gradually increased in suspension-cultured cells following the bacterial infection. GmMMP2 was also activated in response to wounding and dehydration. However, GmMMP2 activation did not correlate with the oxidative burst leading to the hypersensitive response cell death or the tissue senescence progress that involves programmed cell death. Our investigations suggest that GmMMP2 may be involved in a novel defense response of soybean against pathogenic infections.

Applications of retrotransposons as genetic tools in plant biology

Retrotransposons are mobile genetic elements that accomplish transposition via an RNA intermediate that is reverse transcribed before integration into a new location within the host genome. They are ubiquitous in eukaryotic organisms and constitute a major portion of the nuclear genome (often more than half of the total DNA) in plants. Furthermore, they are dispersed as interspersed repetitive sequences throughout most of the length of all host chromosomes. These unique properties of retrotransposons have been exploited as genetic tools for plant genome analysis. Major applications are in determining phylogeny and genetic diversity and in the functional analyses of genes in plants. Here, recent advances in molecular markers, gene tagging and functional genomics technologies using plant retrotransposons are described.

Plant retrotransposons

Retrotransposons are mobile genetic elements that transpose through reverse transcription of an RNA intermediate. Retrotransposons are ubiquitous in plants and play a major role in plant gene and genome evolution. In many cases, retrotransposons comprise over 50% of nuclear DNA content, a situation that can arise in just a few million years. Plant retrotransposons are structurally and functionally similar to the retrotransposons and retroviruses that are found in other eukaryotic organisms. However, there are important differences in the genomic organization of retrotransposons in plants compared to some other eukaryotes, including their often-high copy numbers, their extensively heterogeneous populations, and their chromosomal dispersion patterns. Recent studies are providing valuable insights into the mechanisms involved in regulating the expression and transposition of retrotransposons. This review describes the structure, genomic organization, expression, regulation, and evolution of retrotransposons, and discusses both their contributions to plant genome evolution and their use as genetic tools in plant biology.

Mapping of avirulence genes in the rice blast fungus, Magnaporthe grisea, with RFLP and RAPD markers

Three genetically independent avirulence genes, AVR1-Irat7, AVRI-MedNoi; and AVR1-Ku86, were identified in a cross involving isolates Guy11 and 2/0/3 of the rice blast fungus, Magnaporthe grisea. Using 76 random progeny, we constructed a partial genetic map with restriction fragment length polymorphism (RFLP) markers revealed by probes such as the repeated sequences MGL/MGR583 and Pot3/MGR586, cosmids from the M. grisea genetic map, and a telomere sequence oligonucleotide. Avirulence genes AVR1-MedNoi and AVR1-Ku86 were closely linked to telomere RFLPs such as marker TelG (6 cM from AVR1-MedNoi) and TelF (4.5 cM from AVR1-Ku86). Avirulence gene AVR1-Irat7 was linked to a cosmid RFLP located on chromosome 1 and mapped at 20 cM from the avirulence gene AVR1-CO39. Using bulked segregant analysis, we identified 11 random amplified polymorphic DNA (RAPD) markers closely linked (0 to 10 cM) to the avirulence genes segregating in this cross. Most of these RAPD markers corresponded to junction fragments between known or new transposons and a single-copy sequence. Such junctions or the whole sequences of single-copy RAPD markers were frequently absent in one parental isolate. Single-copy sequences from RAPD markers tightly linked to avirulence genes will be used for positional cloning.

SIRE-1, a copia/Ty1-like retroelement from soybean, encodes a retroviral envelope-like protein

The soybean genome hosts a family of several hundred, relatively homogeneous copies of a large, copia/Ty1-like retroelement designated SIRE-1. A copy of this element has been recovered from a Glycine max genomic library. DNA sequence analysis of two SIRE-1 subclones revealed that SIRE-1 contains a long, uninterrupted, ORF between the 3' end of the pol ORF and the 3' long terminal repeat (LTR), a region that harbors the env gene in retroviral genomes. Conceptual translation of this second ORF produces a 70-kDa protein. Computer analyses of the amino acid sequence predicted patterns of transmembrane domains, alpha-helices, and coiled coils strikingly similar to those found in mammalian retroviral envelope proteins. In addition, a 65-residue, proline-rich domain is characterized by a strong amino acid compositional bias virtually identical to that of the 60-amino acid, proline-rich neutralization domain of the feline leukemia virus surface protein. The assignment of SIRE-1 to the copia/Ty1 family was confirmed by comparison of the conceptual translation of its reverse transcriptase-like domain with those of other retroelements. This finding suggests the presence of a proretrovirus in a plant genome and is the strongest evidence to date for the existence of a retrovirus-like genome closely related to copia/Ty1 retrotransposons.