The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum

Identification of miniature inverted-repeat transposable elements (MITEs) and biogenesis of their siRNAs in the Solanaceae: new functional implications for MITEs

Small RNAs regulate the genome by guiding transcriptional and post-transcriptional silencing machinery to specific target sequences, including genes and transposable elements (TEs). Although miniature inverted-repeat transposable elements (MITEs) are closely associated with euchromatic genes, the broader functional impact of these short TE insertions in genes is largely unknown. We identified 22 families of MITEs in the Solanaceae (MiS1-MiS22) and found abundant MiS insertions in Solanaceae genomic DNA and expressed sequence tags (EST). Several Solanaceae MITEs generate genome changes that potentially affect gene function and regulation, most notably, a MiS insertion that provides a functionally indispensable alternative exon in the tobacco mosaic virus N resistance gene. We show that MITEs generate small RNAs that are primarily 24 nt in length, as detected by Northern blot hybridization and by sequencing small RNAs of Solanum demissum, Nicotiana glutinosa, and Nicotiana benthamiana. Additionally, we show that stable RNAi lines silencing DICER-LIKE3 (DCL3) in tobacco and RNA-dependent RNA polymerase 2 (RDR2) in potato cause a reduction in 24-nt MITE siRNAs, suggesting that, as in Arabidopsis, TE-derived siRNA biogenesis is DCL3 and RDR2 dependent. We provide evidence that DICER-LIKE4 (DCL4) may also play a role in MITE siRNA generation in the Solanaceae.

BAC end sequences corresponding to the B4 resistance gene cluster in common bean: a resource for markers and synteny analyses

In common bean, a complex disease resistance (R) gene cluster, harboring many specific R genes against various pathogens, is located at the end of the linkage group B4. A BAC library of the Meso-american bean genotype BAT93 was screened with PRLJ1, a probe previously shown to be specific to the B4 R gene cluster, leading to the identification of 73 positive BAC clones. BAC-end sequencing (BES) of the 73 positive BACs generated 75 kb of sequence. These BACs were organized into 6 contigs, all mapped at the B4 R gene cluster. To evaluate the potential of BES for marker development, BES-derived specific primers were used to check for linkage with two allelic anthracnose R specificities Co-3 and Co-3 ( 2 ), through the analysis of pairs of Near Isogenic Lines (NILs). Out of 32 primer pairs tested, two revealed polymorphisms between the NILs, confirming the suspected location of Co-3 and Co-3 ( 2 ) at the B4 cluster. In order to identify the orthologous region of the B4 R gene cluster in the two model legume genomes, bean BESs were used as queries in TBLASTX searches of Medicago truncatula and Lotus japonicus BAC clones. Putative orthologous regions were identified on chromosome Mt6 and Lj2, in agreement with the colinearity observed between Mt and Lj for these regions.

Haplotype inference from unphased SNP data in heterozygous polyploids based on SAT

Background

Haplotype inference based on unphased SNP markers is an important task in population genetics. Although there are different approaches to the inference of haplotypes in diploid species, the existing software is not suitable for inferring haplotypes from unphased SNP data in polyploid species, such as the cultivated potato (Solanum tuberosum). Potato species are tetraploid and highly heterozygous.

Results

Here we present the software SATlotyper which is able to handle polyploid and polyallelic data. SATlo-typer uses the Boolean satisfiability problem to formulate Haplotype Inference by Pure Parsimony. The software excludes existing haplotype inferences, thus allowing for calculation of alternative inferences. As it is not known which of the multiple haplotype inferences are best supported by the given unphased data set, we use a bootstrapping procedure that allows for scoring of alternative inferences. Finally, by means of the bootstrapping scores, it is possible to optimise the phased genotypes belonging to a given haplotype inference. The program is evaluated with simulated and experimental SNP data generated for heterozygous tetraploid populations of potato. We show that, instead of taking the first haplotype inference reported by the program, we can significantly improve the quality of the final result by applying additional methods that include scoring of the alternative haplotype inferences and genotype optimisation. For a sub-population of nineteen individuals, the predicted results computed by SATlotyper were directly compared with results obtained by experimental haplotype inference via sequencing of cloned amplicons. Prediction and experiment gave similar results regarding the inferred haplotypes and phased genotypes.

Conclusion

Our results suggest that Haplotype Inference by Pure Parsimony can be solved efficiently by the SAT approach, even for data sets of unphased SNP from heterozygous polyploids. SATlotyper is freeware and is distributed as a Java JAR file. The software can be downloaded from the webpage of the GABI Primary Database at . The application of SATlotyper will provide haplotype information, which can be used in haplotype association mapping studies of polyploid plants.

Analysis of 90 Mb of the potato genome reveals conservation of gene structures and order with tomato but divergence in repetitive sequence composition

Background

The Solanaceae family contains a number of important crop species including potato (Solanum tuberosum) which is grown for its underground storage organ known as a tuber. Albeit the 4^th most important food crop in the world, other than a collection of ~220,000 Expressed Sequence Tags, limited genomic sequence information is currently available for potato and advances in potato yield and nutrition content would be greatly assisted through access to a complete genome sequence. While morphologically diverse, Solanaceae species such as potato, tomato, pepper, and eggplant share not only genes but also gene order thereby permitting highly informative comparative genomic analyses.

Results

In this study, we report on analysis 89.9 Mb of potato genomic sequence representing 10.2% of the genome generated through end sequencing of a potato bacterial artificial chromosome (BAC) clone library (87 Mb) and sequencing of 22 potato BAC clones (2.9 Mb). The GC content of potato is very similar to Solanum lycopersicon (tomato) and other dicotyledonous species yet distinct from the monocotyledonous grass species, Oryza sativa. Parallel analyses of repetitive sequences in potato and tomato revealed substantial differences in their abundance, 34.2% in potato versus 46.3% in tomato, which is consistent with the increased genome size per haploid genome of these two Solanum species. Specific classes and types of repetitive sequences were also differentially represented between these two species including a telomeric-related repetitive sequence, ribosomal DNA, and a number of unclassified repetitive sequences. Comparative analyses between tomato and potato at the gene level revealed a high level of conservation of gene content, genic feature, and gene order although discordances in synteny were observed.

Conclusion

Genomic level analyses of potato and tomato confirm that gene sequence and gene order are conserved between these solanaceous species and that this conservation can be leveraged in genomic applications including cross-species annotation and genome sequencing initiatives. While tomato and potato share genic features, they differ in their repetitive sequence content and composition suggesting that repetitive sequences may have a more significant role in shaping speciation than previously reported.

Genomic organization, rapid evolution and meiotic instability of nucleotide-binding-site-encoding genes in a new fruit crop, "chestnut rose"

From chestnut rose, a promising fruit crop of the Rosa genus, powdery mildew disease-resistant and susceptible genotypes and their F(1) progeny were used to isolate nucleotide-binding-site (NBS)-encoding genes using 19 degenerate primer pairs and an additional cloning method called overlapping extension amplification. A total of 126 genes were harvested; of these, 38 were from a resistant parent, 37 from a susceptible parent, and 51 from F(1) progeny. A phylogenetic tree was constructed, which revealed that NBS sequences from parents and F(1) progeny tend to form a mixture and are well distributed among the branches of the tree. Mapping of these NBS genes suggested that their organization in the genome is a "tandem duplicated cluster" and, to a lesser extent, a "heterogeneous cluster." Intraspecific polymorphisms and interspecific divergence were detected by Southern blotting with NBS-encoding genes as probes. Sequencing on the nucleotide level revealed even more intraspecific variation: for the R4 gene, 9.81% of the nucleotides are polymorphic. Amino acid sites under positive selection were detected in the NBS region. Some NBS-encoding genes were meiotically unstable, which may due to recombination and deletion events. Moreover, a transposon-like element was isolated in the flanking region of NBS genes, implying a possible role for transposon in the evolutionary history of resistance genes.

Frequent sequence exchanges between homologs of RPP8 in Arabidopsis are not necessarily associated with genomic proximity

Disease resistance (R) genes are often clustered in plant genomes and may exhibit heterogeneous rates of evolution. Some (type I R genes) have evolved rapidly through frequent sequence exchanges, while others (type II R genes) have evolved independently and tend to be conserved in different genotypes or related species. The RPP8 resistance gene in Arabidopsis thaliana is located at a complex locus that also harbors the sequence-related resistance genes HRT and RCY1 in different ecotypes. We sequenced 98 homologs of RPP8 from A. thaliana, Arabidopsis arenosa and Arabidopsis lyrata. Three lineages of type II and one lineage of type I RPP8 homologs were identified. Two of the three lineages of type II genes are each represented by a single-copy locus on either chromosomes I or V. Chromosome V contains two small clusters of RPP8 paralogs. One cluster contains both type I and type II genes and the other comprises only type I genes. These multi-copy loci have expanded and contracted through unequal crossovers, which have generated chimeric genes as well as variations in copy number. Sequence exchanges, most likely gene conversions, were detected between RPP8 homologs that are spatially separated by 2.2 Mb and 12 cM. The sequence exchanges between type I homologs within a locus have been more frequent than sequence exchanges between homologs from two different loci, indicating the influence of chromosomal position on the evolution of these R genes. However, physical distance was not the only factor determining the frequency of sequence exchange, because some closely linked paralogs exhibited little sequence exchange. At least two distinct lineages of type II RPP8 homologs were identified in different species, with obvious allelic/orthologous relationships within each lineage. Therefore, the differentiation of type I and type II RPP8 homologs seems to have occurred before speciation of A. thaliana, A. arenosa and A. lyrata.

Identification and fine mapping of Pi39(t), a major gene conferring the broad-spectrum resistance to Magnaporthe oryzae

Blast, caused by the ascomycete fungus Magnaporthe oryzae, is one of the most devastating diseases of rice worldwide. The Chinese native cultivar (cv.) Q15 expresses the broad-spectrum resistance to most of the isolates collected from China. To effectively utilize the resistance, three rounds of linkage analysis were performed in an F(2) population derived from a cross of Q15 and a susceptible cv. Tsuyuake, which segregated into 3:1 (resistant/susceptible) ratio. The first round of linkage analysis employing simple sequence repeat (SSR) markers was carried out in the F(2) population through bulked-segregant assay. A total of 180 SSR markers selected from each chromosome equally were surveyed. The results revealed that only two polymorphic markers, RM247 and RM463, located on chromosome 12, were linked to the resistance (R) gene. To further define the chromosomal location of the R gene locus, the second round of linkage analysis was performed using additional five SSR markers, which located in the region anchored by markers RM247 and RM463. The locus was further mapped to a 0.27 cM region bounded by markers RM27933 and RM27940 in the pericentromeric region towards the short arm. For fine mapping of the R locus, seven new markers were developed in the smaller region for the third round of linkage analysis, based on the reference sequences. The R locus was further mapped to a 0.18 cM region flanked by marker clusters 39M11 and 39M22, which is closest to, but away from the Pita/Pita(2) locus by 0.09 cM. To physically map the locus, all the linked markers were landed on the respective bacterial artificial chromosome clones of the reference cv. Nipponbare. Sequence information of these clones was used to construct a physical map of the locus, in silico, by bioinformatics analysis. The locus was physically defined to an interval of approximately 37 kb. To further characterize the R gene, five R genes mapped near the locus, as well as 10 main R genes those might be exploited in the resistance breeding programs, were selected for differential tests with 475 Chinese isolates. The R gene carrier Q15 conveys resistances distinct from those conditioned by the carriers of the 15 R genes. Together, this valuable R gene was, therefore, designated as Pi39(t). The sequence information of the R gene locus could be used for further marker-based selection and cloning.

Comparative sequence analysis of Solanum and Arabidopsis in a hot spot for pathogen resistance on potato chromosome V reveals a patchwork of conserved and rapidly evolving genome segments

Background

Quantitative phenotypic variation of agronomic characters in crop plants is controlled by environmental and genetic factors (quantitative trait loci = QTL). To understand the molecular basis of such QTL, the identification of the underlying genes is of primary interest and DNA sequence analysis of the genomic regions harboring QTL is a prerequisite for that. QTL mapping in potato (Solanum tuberosum) has identified a region on chromosome V tagged by DNA markers GP21 and GP179, which contains a number of important QTL, among others QTL for resistance to late blight caused by the oomycete Phytophthora infestans and to root cyst nematodes.

Results

To obtain genomic sequence for the targeted region on chromosome V, two local BAC (bacterial artificial chromosome) contigs were constructed and sequenced, which corresponded to parts of the homologous chromosomes of the diploid, heterozygous genotype P6/210. Two contiguous sequences of 417,445 and 202,781 base pairs were assembled and annotated. Gene-by-gene co-linearity was disrupted by non-allelic insertions of retrotransposon elements, stretches of diverged intergenic sequences, differences in gene content and gene order. The latter was caused by inversion of a 70 kbp genomic fragment. These features were also found in comparison to orthologous sequence contigs from three homeologous chromosomes of Solanum demissum, a wild tuber bearing species. Functional annotation of the sequence identified 48 putative open reading frames (ORF) in one contig and 22 in the other, with an average of one ORF every 9 kbp. Ten ORFs were classified as resistance-gene-like, 11 as F-box-containing genes, 13 as transposable elements and three as transcription factors. Comparing potato to Arabidopsis thaliana annotated proteins revealed five micro-syntenic blocks of three to seven ORFs with A. thaliana chromosomes 1, 3 and 5.

Conclusion

Comparative sequence analysis revealed highly conserved collinear regions that flank regions showing high variability and tandem duplicated genes. Sequence annotation revealed that the majority of the ORFs were members of multiple gene families. Comparing potato to Arabidopsis thaliana annotated proteins suggested fragmented structural conservation between these distantly related plant species.

Introgressed and endogenous Mi-1 gene clusters in tomato differ by complex rearrangements in flanking sequences and show sequence exchange and diversifying selection among homologues

Many plant disease resistance genes (R-genes) encode proteins characterized by the presence of a nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) region and occur in clusters of related genes in plant genomes. One such gene, Mi-1, confers isolate-specific resistance against root-knot nematodes, aphids and whiteflies in cultivated tomato, Solanum lycopersicon. The DNA region carrying Mi-1 and six closely related sequences was introgressed into tomato from Solanum peruvianum in the 1940s. For both susceptible and resistant tomato, Mi-1 homologues are present in two clusters with 3 and 4 copies each on the short arm of chromosome 6. Two homologues from each source are pseudogenes, and one homologue from each source encodes a truncated product. DNA sequence identity among the homologues including the truncated genes, but excluding the pseudogenes, ranges from 92.9 to 96.7%. All the non-pseudogene homologues are transcribed. Comparison of homologues suggests that extensive sequence exchange has occurred. Regions of diversifying selection are present in the ARC2 domain of the NBS region and dispersed throughout the LRR region, suggesting that these regions are possible locations of specificity determinants. Other sequences in the introgressed region have similarity to the Arabidopsis auxin-receptor protein TIR1, a jumonji-like transcription factor and a Na(+)/H(+) antiporter. Analysis of sequences flanking the Mi-1-homologues reveals blocks of homology, but complex differences in arrangement of these blocks when susceptible and resistant genotypes are compared indicating that the region has undergone considerable rearrangement during evolution, perhaps contributing to evolution of specificity.

Recent insights into R gene evolution

SUMMARY Plants are under strong evolutionary pressure to maintain surveillance against pathogens. Resistance (R) gene-dependent recognition of pathogen avirulence (Avr) determinants plays a major role in plant defence. Here we highlight recent insights into the molecular mechanisms and selective forces that drive the evolution of NB-LRR (nucleotide binding-leucine-rich repeat) resistance genes. New implications for models of R gene evolution have been raised by demonstrations that R proteins can detect cognate Avr proteins indirectly by 'guarding' virulence targets, and by evidence that R protein signalling is regulated by intramolecular interactions between different R functional domains. Comparative genomic surveys of NB-LRR diversity in different species have revealed ancient NB-LRR lineages that are unequally represented among plant taxa, consistent with a Birth and Death Model of evolution. The physical distribution of NB-LRRs in plant genomes indicates that tandem and segmental duplication are important factors in R gene proliferation. The majority of R genes reside in clusters, and the frequency of recombination between clustered genes can vary strikingly, even within a single cluster. Biotic and abiotic factors have been shown to increase the frequency of recombination in reporter transgene-based assays, suggesting that external stressors can affect genome stability. Fitness penalties have been associated with some R genes, and population studies have provided evidence for maintenance of ancient R allelic diversity by balancing selection. The available data suggest that different R genes can follow strikingly distinct evolutionary trajectories, indicating that it will be difficult to formulate universally applicable models of R gene evolution.