Annotation of a 95-kb Populus deltoides genomic sequence reveals a disease resistance gene cluster and novel class I and class II transposable elements

Identification of Genes Underlying the Resistance to Melampsora larici-populina in an R Gene Supercluster of the Populus deltoides Genome

Identification of the particular genes in an R genes supercluster underlying resistance to the rust fungus Melampsora larici-populina in poplar genome remains challenging. Based on the de novo assembly of the Populus deltoides genome, all of the detected major genetic loci conferring resistance to M. larici-populina were confined to a 3.5-Mb region on chromosome 19. The transcriptomes of the resistant and susceptible genotypes were sequenced for a timespan from 0 to 168 hours postinoculation. By mapping the differentially expressed genes to the target genomic region, we identified two constitutive expression R genes and one inducible expression R gene that might confer resistance to M. larici-populina. Nucleotide variations were predicted based on the reconstructed haplotypes for each allele of the candidate genes. We also confirmed that salicylic acid was the phytohormone mediating signal transduction pathways, and PR-1 was identified as a key gene inhibiting rust reproduction. Finally, quantitative reverse transcription PCR assay revealed consistent expressions with the RNA-sequencing data for the detected key genes. This study presents an efficient approach for the identification of particular genes underlying phenotype of interest by the combination of genetic mapping, transcriptome profiling, and candidate gene sequences dissection. The identified key genes would be useful for host resistance diagnosis and for molecular breeding of elite poplar cultivars exhibiting resistance to M. larici-populina infection. The detected R genes are also valuable for testing whether the combination of individual R genes can induce durable quantitative resistance.

QTLs for Woolly Poplar Aphid (Phloeomyzus passerinii L.) Resistance Detected in an Inter-Specific Populus deltoides x P. nigra Mapping Population

The genus Populus represents one of the most economically important groups of forest trees. It is composed by approximately 30 species used for wood and non-wood products, phytoremediation and biomass. Poplar is subjected to several biological and environmental threats although, compared to annual crops, we know far less about the genetic bases of biotic stress resistance. Woolly poplar aphid (Phloeomyzus passerinii) is considered a main pest of cultivated poplars in European and American countries. In this work we present two high density linkage maps in poplar obtained by a genotyping by sequencing (GBS) approach and the identification of QTLs involved in Ph. passerinii resistance. A total of 5,667 polymorphic markers (5,606 SNPs and 61 SSRs) identified on expressed sequences have been used to genotype 131 plants of an F1 population P ×canadensis obtained by an interspecific mate between Populus deltoides (resistant to woolly poplar aphid) and Populus nigra (susceptible to woolly poplar aphid). The two linkage maps, obtained following the two-way pseudo-testcross mapping strategy, have been used to investigate the genetic bases of woolly poplar aphid resistance. One major QTL and two QTLs with minor effects (mapped on LGV, LGXVI and LG XIX) explaining the 65.8% of the genetic variance observed in the progeny in response to Ph. passerinii attack were found. The high density coverage of functional markers allowed the identification of three genes belonging to disease resistance pathway as putative candidates for P. deltoides resistance to woolly poplar aphid. This work is the first report on genetic of woolly poplar aphid genetic resistance and the resistant loci associated markers identified represent a valuable tool in resistance poplar breeding programs.

Association analysis identifies Melampsora ×columbiana poplar leaf rust resistance SNPs

Populus species are currently being domesticated through intensive time- and resource-dependent programs for utilization in phytoremediation, wood and paper products, and conversion to biofuels. Poplar leaf rust disease can greatly reduce wood volume. Genetic resistance is effective in reducing economic losses but major resistance loci have been race-specific and can be readily defeated by the pathogen. Developing durable disease resistance requires the identification of non-race-specific loci. In the presented study, area under the disease progress curve was calculated from natural infection of Melampsora ×columbiana in three consecutive years. Association analysis was performed using 412 P. trichocarpa clones genotyped with 29,355 SNPs covering 3,543 genes. We found 40 SNPs within 26 unique genes significantly associated (permutated P<0.05) with poplar rust severity. Moreover, two SNPs were repeated in all three years suggesting non-race-specificity and three additional SNPs were differentially expressed in other poplar rust interactions. These five SNPs were found in genes that have orthologs in Arabidopsis with functionality in pathogen induced transcriptome reprogramming, Ca²⁺/calmodulin and salicylic acid signaling, and tolerance to reactive oxygen species. The additive effect of non-R gene functional variants may constitute high levels of durable poplar leaf rust resistance. Therefore, these findings are of significance for speeding the genetic improvement of this long-lived, economically important organism.

Qualitative and quantitative resistances to leaf rust finely mapped within two nucleotide-binding site leucine-rich repeat (NBS-LRR)-rich genomic regions of chromosome 19 in poplar

• R(US) is a major dominant gene controlling quantitative resistance, inherited from Populus trichocarpa, whereas R(1) is a gene governing qualitative resistance, inherited from P. deltoides. • Here, we report a reiterative process of concomitant fine-scale genetic and physical mapping guided by the P. trichocarpa genome sequence. The high-resolution linkage maps were developed using a P. deltoides × P. trichocarpa progeny of 1415 individuals. R(US) and R(1) were mapped in a peritelomeric region of chromosome 19. Markers closely linked to R(US) were used to screen a bacterial artificial chromosome (BAC) library constructed from the P. trichocarpa parent, heterozygous at the locus R(US) . • Two local physical maps were developed, one encompassing the R(US) allele and the other spanning r(US) . The alignment of the two haplophysical maps showed structural differences between haplotypes. The genetic and physical maps were anchored to the genome sequence, revealing genome sequence misassembly. Finally, the R(US) locus was localized within a 0.8-cM interval, whereas R(1) was localized upstream of R(US) within a 1.1-cM interval. • The alignment of the genetic and physical maps with the local reorder of the chromosome 19 sequence indicated that R(US) and R(1) belonged to a genomic region rich in nucleotide-binding site leucine-rich repeat (NBS-LRR) and serine threonine kinase (STK) genes.

The Ma gene for complete-spectrum resistance to Meloidogyne species in Prunus is a TNL with a huge repeated C-terminal post-LRR region

Root-knot nematode (RKN) Meloidogyne species are major polyphagous pests of most crops worldwide, and cultivars with durable resistance are urgently needed because of nematicide bans. The Ma gene from the Myrobalan plum (Prunus cerasifera) confers complete-spectrum, heat-stable, and high-level resistance to RKN, which is remarkable in comparison with the Mi-1 gene from tomato (Solanum lycopersicum), the sole RKN resistance gene cloned. We report here the positional cloning and the functional validation of the Ma locus present at the heterozygous state in the P.2175 accession. High-resolution mapping totaling over 3,000 segregants reduced the Ma locus interval to a 32-kb cluster of three Toll/Interleukin1 Receptor-Nucleotide Binding Site-Leucine-Rich Repeat (LRR) genes (TNL1-TNL3), including a pseudogene (TNL2) and a truncated gene (TNL3). The sole complete gene in this interval (TNL1) was validated as Ma, as it conferred the same complete-spectrum and high-level resistance (as in P.2175) using its genomic sequence and native promoter region in Agrobacterium rhizogenes-transformed hairy roots and composite plants. The full-length cDNA (2,048 amino acids) of Ma is the longest of all Resistance genes cloned to date. Its TNL structure is completed by a huge post-LRR (PL) sequence (1,088 amino acids) comprising five repeated carboxyl-terminal PL exons with two conserved motifs. The amino-terminal region (213 amino acids) of the LRR exon is conserved between alleles and contrasts with the high interallelic polymorphisms of its distal region (111 amino acids) and of PL domains. The Ma gene highlights the importance of these uncharacterized PL domains, which may be involved in pathogen recognition through the decoy hypothesis or in nuclear signaling.

Over-expression of the triploid white poplar PtDrl01 gene in tobacco enhances resistance to tobacco mosaic virus

A full-length cDNA, designated as the Populus tomentosa disease resistance-like 01 (PtDrl01) gene, was isolated from triploid white poplar [(Populus tomentosa × P. bolleana) × P. tomentosa]. The protein thought to be produced by the PtDrl01 gene contains a nuclear localisation sequence (NLS), a toll/interleukin-1 receptor (TIR) homologue region, a nucleotide binding site (NBS) and a leucine-rich repeat (LRR) domain. The protein also exhibits a considerable degree of homology to N-like resistance proteins. Real-time quantitative RT-PCR analysis revealed that expression of the PtDrl01 gene in triploid white poplar leaves could be induced by two defence signalling molecules: methyl jasmonate (MeJA) and salicylic acid (SA). Over-expression of the PtDrl01 gene in transgenic tobacco induced enhanced resistance to tobacco mosaic virus (TMV). Long-term resistance from the PtDrl01 gene to TMV infection was also observed in transgenic tobacco plants. Additionally, over-expression of the PtDrl01 gene resulted in transcriptional changes in genes expressing pathogenesis-related proteins in transgenic tobacco under non-stress conditions. These data strongly suggest that the PtDrl01 gene is involved in plant defence responses to pathogen infection.

Mutations in retrotransposon AtCOPIA4 compromises resistance to Hyaloperonospora parasitica in Arabidopsis thaliana

Retrotransposons (RTEs) are a principal component of most eukaryotic genomes, representing 50%-80% of some grass genomes. RTE sequences have been shown to be preferentially present in disease resistance gene clusters in plants. Arabidopsis thaliana has over 1,600 annotated RTE sequences and 56 of these appear to be expressed because of the exact expressed sequence tag (EST) matches and the presence of intact open reading frames. Of the 22 represented in the Affymetrix ATH1 array, AtCOPIA4 was found to be expressed at a higher level than all other RTEs across different developmental stages. Since AtCOPIA4 is located in the RPP5 gene cluster and is adjacent to RPP4 which confers resistance to the downy mildew oomycete Hyaloperonospora parasitica isolate EMWA1, we evaluated AtCOPIA4 mutants for resistance to this pathogen. T-DNA insertional and antisense knockout of AtCOPIA4 was found to reduce the resistance of wild type plants by 2-4 folds. Our results suggest that retrotransposon can be exapted to participate in plant defense response.

Functional analysis of 5' untranslated region of a TIR-NBS-encoding gene from triploid white poplar

Genome-wide analyses have identified a set of TIR-NBS-encoding genes in plants. However, the molecular mechanism underlying the expression of these genes is still unknown. In this study, we presented a TIR-NBS-encoding gene, PtDrl02, that displayed a low level of tissue-specific expression in a triploid white poplar [(Populus tomentosa x P. bolleana) x P. tomentosa], and analyzed the effects of the 5' untranslated region (UTR) on gene expression. The 5' UTR sequence repressed the reporter activity of beta-glucuronidase (GUS) gene under PtDrl02 promoter by 113.5-fold with a staining ratio of 2.97% in the transgenic tobacco plants. Quantitative RT-PCR assays revealed that the 5' UTR sequence decreased the transcript level of the GUS reporter gene by 13.3-fold, implying a regulatory role of 5' UTR in transcription and/or mRNA destabilization. The comparison of GUS activity with the transcript abundance indicated that the 5' UTR sequence decreased the translation efficiency of target gene by 88.3%. Additionally, the analysis of the transgenic P-985/UTRDelta/GUS plants showed that both the exon1 sequence and the leading intron within the 5' UTR region were responsible for the regulation of gene expression. Our results suggested a negative effect of the 5' UTR of PtDrl02 gene on gene expression.

Characterization of resistance gene analogs with a nucleotide binding site isolated from a triploid white poplar

The majority of cloned plant disease resistance genes (R genes) encode a nucleotide binding site (NBS) and a leucine-rich repeat (LRR) domain. In this study, to better understand the R genes in white poplar, 59 resistance gene analogues (RGAs) were identified from a triploid white poplar [(Populus tomentosa x Populus bolleana) x P. tomentosa], based on conserved NBS regions. The 59 RGAs were phylogenetically classified into 10 subfamilies, and 54 RGAs with open-reading frames (ORFs) were further grouped into two classes, toll and interleukin-1 receptor (TIR) and non-TIR. BLAST searches with reference to the genomic sequence of Populus trichocarpa found 96 highly homologous regions distributed in 37 loci, suggesting the abundance and divergence of NBS-encoding genes in the triploid poplar genome. Within subfamilies 1-3, the average non-synonymous/synonymous substitution (omega) rates were < 1, indicating purifying selection on these RGAs, but some sites were clearly under diversifying selection with omega > 1. Many intergenic exchanges were also detected among these RGAs, indicating a probable role in homogenising NBS domains. Quantitative real-time PCR analysis revealed dramatic variations in the transcript level of 18 RGAs in the mature leaves, bark and roots of the triploid poplar, and identified two RGAs that had significantly higher level of transcripts in bark, four RGAs in mature leaves, and 14 in the above-ground portion of poplars, suggesting their probable roles in resistance against diseases attacking the organs. Our results shed light on genetic resources of poplar resistance and will be useful for further resistance gene isolation and exploitation.

Genome-wide identification of NBS resistance genes in Populus trichocarpa

As the largest class of disease resistance R genes, the genes encoding nucleotide binding site and leucine-rich repeat proteins ("NBS-LRR genes") play a critical role in defending plants from a multitude of pathogens and pests. The diversity of NBS-LRR genes was examined in the Populus trichocarpa draft genome sequence. The NBS class of genes in this perennial tree is large and diverse, comprised of approximately 400 genes, at least twice the complement of Arabidopsis. The NBS family can be divided into multiple subfamilies with distinct domain organizations. It includes 119 Coiled-Coil-NBS-LRR genes, 64 TIR-NBS-LRR genes, 34 BED-finger-NBS-LRR, and both truncated and unusual NBS- and NBS-LRR-containing genes. The transcripts of only 34 NBS-LRR genes were detected in rust-infected and non-infected leaves using a whole-genome oligoarray. None showed an altered expression two days post inoculation.

Poplar leaf rusts: model pathogens for a model treeThis minireview is one of a selection of papers published in the Special Issue on Poplar Research in Canada

With the availability of the entire genome of the model tree Populus trichocarpa Torr. & A. Gray and the current genome sequencing project of its rust pathogen Melampsora larici-populina Kleb., rust–poplar interaction research has entered the genomic era. Recent genomics research on poplars has attempted to connect the genetic localizations of loci for qualitative and quantitative disease resistance with putative genes encoding resistance or signalling proteins. The interactions between these putative resistance genes and rust effectors remain unknown. Genomic resources developed for Melampsora spp. promise to contribute to our understanding of the molecular basis of pathogenicity by facilitating the isolation of pathogenicity genes. A multifaceted approach for the identification of such genes that relies largely on trimming and sequence data analysis has been developed. The strategy takes advantage of the resources available and combines EST libraries, bioinformatics data mining for extracellularly expressed secreted proteins, intra- and inter-specific comparative genomics, and testing for the presence of positive selection. It has resulted in the discovery of several putative candidate genes. In silico evidence for candidate genes will be further validated by robust experimental evidence through functional analyses.

Transcript profiling of poplar leaves upon infection with compatible and incompatible strains of the foliar rust Melampsora larici-populina

To understand key processes governing defense mechanisms in poplar (Populus spp.) upon infection with the rust fungus Melampsora larici-populina, we used combined histological and molecular techniques to describe the infection of Populus trichocarpa x Populus deltoides 'Beaupré' leaves by compatible and incompatible fungal strains. Striking differences in host-tissue infection were observed after 48-h postinoculation (hpi) between compatible and incompatible interactions. No reactive oxygen species production could be detected at infection sites, while a strong accumulation of monolignols occurred in the incompatible interaction after 48 hpi, indicating a late plant response once the fungus already penetrated host cells to form haustorial infection structures. P. trichocarpa whole-genome expression oligoarrays and sequencing of cDNAs were used to determine changes in gene expression in both interactions at 48 hpi. Temporal expression profiling of infection-regulated transcripts was further compared by cDNA arrays and reverse transcription-quantitative polymerase chain reaction. Among 1,730 significantly differentially expressed transcripts in the incompatible interaction, 150 showed an increase in concentration > or =3-fold, whereas 62 were decreased by > or =3-fold. Regulated transcripts corresponded to known genes targeted by R genes in plant pathosystems, such as inositol-3-P synthase, glutathione S-transferases, and pathogenesis-related proteins. However, the transcript showing the highest rust-induced up-regulation encodes a putative secreted protein with no known function. In contrast, only a few transcripts showed an altered expression in the compatible interaction, suggesting a delay in defense response between incompatible and compatible interactions in poplar. This comprehensive analysis of early molecular responses of poplar to M. larici-populina infection identified key genes that likely contain the fungus proliferation in planta. Sequences from the SSH library described in this article can be retrieved in GenBank under accession numbers CT 027996 to CT 029994 and CT 033829.

Experimental validation of novel genes predicted in the un-annotated regions of the Arabidopsis genome

Background

Several lines of evidence support the existence of novel genes and other transcribed units which have not yet been annotated in the Arabidopsis genome. Two gene prediction programs which make use of comparative genomic analysis, Twinscan and EuGene, have recently been deployed on the Arabidopsis genome. The ability of these programs to make use of sequence data from other species has allowed both Twinscan and EuGene to predict over 1000 genes that are intergenic with respect to the most recent annotation release. A high throughput RACE pipeline was utilized in an attempt to verify the structure and expression of these novel genes.

Results

1,071 un-annotated loci were targeted by RACE, and full length sequence coverage was obtained for 35% of the targeted genes. We have verified the structure and expression of 378 genes that were not present within the most recent release of the Arabidopsis genome annotation. These 378 genes represent a structurally diverse set of transcripts and encode a functionally diverse set of proteins.

Conclusion

We have investigated the accuracy of the Twinscan and EuGene gene prediction programs and found them to be reliable predictors of gene structure in Arabidopsis. Several hundred previously un-annotated genes were validated by this work. Based upon this information derived from these efforts it is likely that the Arabidopsis genome annotation continues to overlook several hundred protein coding genes.

Populus: a model system for plant biology

With the completion of the Populus trichocarpa genome sequence and the development of various genetic, genomic, and biochemical tools, Populus now offers many possibilities to study questions that cannot be as easily addressed in Arabidopsis and rice, the two prime model systems of plant biology and genomics. Tree-specific traits such as wood formation, long-term perennial growth, and seasonality are obvious areas of research, but research in other areas such as control of flowering, biotic interactions, and evolution of adaptive traits is enriched by adding a tree to the suite of model systems. Furthermore, the reproductive biology of Populus (a dioeceous wind-pollinated long-lived tree) offers both new possibilities and challenges in the study and analysis of natural genetic and phenotypic variation. The relatively close phylogenetic relationship of Populus to Arabidopsis in the Eurosid clade of Eudicotyledonous plants aids in comparative functional studies and comparative genomics, and has the potential to greatly facilitate studies on genome and gene family evolution in eudicots.

Identification of transposons, retroelements, and a gene family predominantly expressed in floral tissues in chromosome 3DS of the hexaploid wheat progenitor Aegilops tauschii

A multigene family expressed during early floral development was identified on the short arm of wheat chromosome 3D in the region of the Ph2 locus, a locus controlling homoeologous chromosome pairing in allohexaploid wheat. Physical, genetic and molecular characterisation of the Wheat Meiosis 1 (WM1) gene family identified seven members that localised within a region of 173-kb. WM1 gene family members were sequenced and they encode mainly type Ia plasma membrane-anchored leucine rich repeat-like receptor proteins. In situ expression profiling suggests the gene family is predominantly expressed in floral tissue. In addition to the WM1 gene family, a number of other genes, gene fragments and pseudogenes were identified. It has been predicted that there is approximately one gene every 19-kb and that this region of the wheat genome contains 23 repetitive elements including BARE-1 and Wis2-1 like sequences. Nearly 50% of the repetitive elements identified were similar to known transposons from the CACTA superfamily. Ty1-copia, Ty3-gypsy and Athila LTR retroelements were also prevalent within the region. The WM1 gene cluster is present on 3DS and on barley 3HS but missing from the A and B genomes of hexaploid wheat. This suggests either recent generation of the cluster or specific deletion of the cluster during wheat polyploidisation. The evolutionary significance of the cluster, its possible roles in disease response or floral and early meiotic development and its location at or near the Ph2 locus are discussed.

Biotechnology and the domestication of forest trees

Wood is one of the major renewable materials. To compensate for the ever-increasing demand for wood and to reduce pressure on native forests, more wood of higher quality will need to be produced on less land by planting highly productive trees. Biotechnology has shown great promise for forest tree improvement and over the past 10 years this field has flourished. Not only has the potential of transgenic trees with optimized yield and quality traits been demonstrated in field trials, but progress in genetical genomics and association genetics promise quantum leaps forward for tree improvement.

Genetic architecture of qualitative and quantitative Melampsora larici-populina leaf rust resistance in hybrid poplar: genetic mapping and QTL detection

In order to elucidate the genetic control of resistance to Melampsora larici-populina leaf rust in hybrid poplars, a Populus deltoides x P. trichocarpa F(1) progeny was analysed for qualitative and quantitative rust resistances. This progeny was evaluated for three components of quantitative resistance (latent period, uredinia number and uredinia size) to seven M. larici-populina strains in controlled conditions, and for one component of field susceptibility (rust colonization on the most infected leaf). One qualitative resistance locus inherited from P. deltoides, R(1), was localized on the genetic map. It segregates 1 : 1 in the F(1) progeny and is effective against four of the studied strains. QTL analysis was performed separately on R(1) and r(1) genotype subsets. An additional detection was conducted on the entire F(1) progeny for the three strains able to overcome R(1) and for MAX2. A total of nine QTLs were detected. Two had large, broad-spectrum effects. One (R(US)) is inherited from the P. trichocarpa parent; the other is inherited from P. deltoides and colocalized with R(1). Seven QTLs had only limited and specific effects. Significant interaction effects were detected mainly between the two major QTLs. Implications of these results for durable resistance breeding strategies, and possible benefits from the Populus genome sequence, are discussed.

Characterization and mapping of NBS-LRR resistance gene analogs in apricot (Prunus armeniaca L.)

Genomic DNA sequences sharing homology with the NBS-LRR (nucleotide binding site-leucine-rich repeat) resistance genes were isolated and cloned from apricot (Prunus armeniaca L.) using a PCR approach with degenerate primers designed from conserved regions of the NBS domain. Restriction digestion and sequence analyses of the amplified fragments led to the identification of 43 unique amino acid sequences grouped into six families of resistance gene analogs (RGAs). All of the RGAs identified belong to the Toll-Interleukin receptor (TIR) group of the plant disease resistance genes (R-genes). RGA-specific primers based on non-conserved regions of the NBS domain were developed from the consensus sequences of each RGA family. These primers were used to develop amplified fragment length polymorphism (AFLP)-RGA markers by means of an AFLP-modified procedure where one standard primer is substituted by an RGA-specific primer. Using this method, 27 polymorphic markers, six of which shared homology with the TIR class of the NBS-LRR R-genes, were obtained from 17 different primer combinations. Of these 27 markers, 16 mapped in an apricot genetic map previously constructed from the self-pollination of the cultivar Lito. The development of AFLP-RGA markers may prove to be useful for marker-assisted selection and map-based cloning of R-genes in apricot.

Genetic and physical mapping of Melampsora rust resistance genes in Populus and characterization of linkage disequilibrium and flanking genomic sequence

•  In an attempt to elucidate the molecular mechanisms of Melampsora rust resistance in Populus trichocarpa, we have mapped two resistance loci, MXC3 and MER, and intensively characterized the flanking genomic sequence for the MXC3 locus and the level of linkage disequilibrium (LD) in natural populations. •  We used an interspecific backcross pedigree and a genetic map that was highly saturated with AFLP and SSR markers, and assembled shotgun-sequence data in the region containing markers linked to MXC3. •  The two loci were mapped to different linkage groups. Linkage disequilibrium for MXC3 was confined to two closely linked regions spanning 34 and 16 kb, respectively. The MXC3 region also contained six disease-resistance candidate genes. •  The MER and MXC3 loci are clearly distinct, and may have different mechanisms of resistance, as different classes of putative resistance genes were present near each locus. The suppressed recombination previously observed in the MXC3 region was possibly caused by extensive hemizygous rearrangements confined to the original parent tree. The relatively low observed LD may facilitate association studies using candidate genes for rust resistance, but will probably inhibit marker-aided selection.