Extraordinary transgressive phenotypes of hybrid tomato are influenced by epigenetics and small silencing RNAs
2011; Springer Nature; Volume: 31; Issue: 2 Linguagem: Inglês
10.1038/emboj.2011.458
ISSN1460-2075
AutoresP. V. Shivaprasad, Ruth M. Dunn, Bruno Santos, Andrew Bassett, David C. Baulcombe,
Tópico(s)Chromosomal and Genetic Variations
ResumoArticle16 December 2011free access Extraordinary transgressive phenotypes of hybrid tomato are influenced by epigenetics and small silencing RNAs Padubidri V Shivaprasad Padubidri V Shivaprasad Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Ruth M Dunn Ruth M Dunn Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Bruno ACM Santos Bruno ACM Santos Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Andrew Bassett Andrew Bassett Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author David C Baulcombe Corresponding Author David C Baulcombe Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Padubidri V Shivaprasad Padubidri V Shivaprasad Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Ruth M Dunn Ruth M Dunn Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Bruno ACM Santos Bruno ACM Santos Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Andrew Bassett Andrew Bassett Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author David C Baulcombe Corresponding Author David C Baulcombe Department of Plant Sciences, University of Cambridge, Cambridge, UK Search for more papers by this author Author Information Padubidri V Shivaprasad1, Ruth M Dunn1, Bruno ACM Santos1, Andrew Bassett1 and David C Baulcombe 1 1Department of Plant Sciences, University of Cambridge, Cambridge, UK *Corresponding author. Department of Plant Sciences, University of Cambridge, Downing Street, Cambridge CB2 3EA, UK. Tel.: +44 1223 339386; Fax: +44 1223 333953; E-mail: [email protected] The EMBO Journal (2012)31:257-266https://doi.org/10.1038/emboj.2011.458 There is a Have you seen? (January 2012) associated with this Article. PDFDownload PDF of article text and main figures. Peer ReviewDownload a summary of the editorial decision process including editorial decision letters, reviewer comments and author responses to feedback. ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InMendeleyWechatReddit Figures & Info Hybrid organisms may fail to develop, be sterile or they may be more vigorous than either of the parents. Examples of hybrid vigour or hybrid necrosis in the F1 are often not inherited stably in subsequent generations if they are associated with overdominance. There can also be transgressive phenotypes that are inherited stably in these later generations, but the underlying mechanisms are not well understood. Here we have investigated the possibility that stable transgressive phenotypes in the progeny of crosses between cultivated tomato (Solanum lycopersicum cv. M82) and a wild relative (Solanum pennellii, accession LA716) are associated with micro or small interfering(si) RNAs. We identified loci from which these small(s)RNAs were more abundant in hybrids than in either parent and we show that accumulation of such transgressive sRNAs correlated with suppression of the corresponding target genes. In one instance this effect was associated with hypermethylation of the corresponding genomic DNA. Our results illustrate a potential role of transgressive sRNAs in plant breeding and in natural evolution with wild plants. Introduction The phenotypes of hybrid plants and animals are primarily a consequence of genetic effects at one or more loci. With a single-gene trait the phenotype in the hybrid reflects the activity of the dominant allele, whereas, with more complex quantitative traits, the phenotype will be affected by the composite action of the multiple loci and will normally be intermediate between the parents. However, there are exceptions in which hybrids have more extreme phenotypes than the parents. The best-characterised exceptions are in the F1 hybrids of plants and animals, in which there is enhanced vigour (heterosis) or hybrid necrosis. These effects may be due to: complementation by dominant genes from one parent of defective or missing alleles in the other; epistasis between unlinked loci or heterozygote advantage (Bomblies and Weigel, 2007; Lippman and Zamir, 2007; Birchler et al, 2010; Chen, 2010). These F1 hybrid phenotypes are lost in later generations because they are caused by alleles or gene combinations that are inherited in only a proportion of the progeny (Baack and Rieseberg, 2007). There are also heritable phenotypes in hybrids and hybrid lineages that are extreme to those of either parent. The heritability of these phenotypes indicates that they are distinct from those involved in heterosis or hybrid necrosis (Rieseberg et al, 1999), and the term transgressive segregation is often used to describe the niche divergence and phenotypic novelty of these hybrid lineages: the phenotypes transgress the parental range. Many eukaryotes exhibit transgressive segregation, but it is more frequent in plants than animals (Rieseberg et al, 2003). The proposed mechanisms of transgressive segregation involve complementation and epistasis as in heterosis. Complementation would occur if defective genes from one parent are compensated by their functional homologues from the other parent. Epistasis would involve genes from one parent that interact in trans with non-homologous loci from the opposite parent (Rieseberg et al, 1999). Transgressive segregation is often associated with chromosome rearrangements (Rieseberg et al, 1996; Shaked et al, 2001; Lai et al, 2005), transposable element mobilization (Liu and Wendel, 2000; Michalak, 2009), DNA methylation (Shaked et al, 2001; Salmon et al, 2005) and changes to gene expression. Presumably, these changes are secondary to the initiating events associated with complementation or epistasis. Any epistatic interactions associated with transgressive phenotypes could be mediated by proteins, DNA or RNA in various combinations affecting protein complexes or gene expression (Caicedo et al, 2004; Werner et al, 2005; Rowe et al, 2008). A previously unexplored possibility is based on RNA silencing at the chromatin or posttranscriptional level (Baulcombe, 2004), in which small(s)RNAs from one parent find novel targets in the genome or transcriptome of the opposite parent. Such interactions would establish patterns of gene expression at either the transcriptional or posttranscriptional level that would be specific to the hybrids. To investigate the possible involvement of RNA silencing in hybrid phenotypes, we used high-throughput sequencing to characterise sRNAs in young seedlings of S. lycopersicum cv. M82 (hereafter abbreviated as M82) and S. pennellii (accession LA716) and in introgression lines (ILs), in which individual defined regions of the S. pennellii genome are present in the background of the M82 genome (Supplementary Figure S1A; Eshed and Zamir, 1995; Lippman et al, 2007). These lines have been studied extensively for their habit, leaf dissection, fruit morphology/yield and metabolic traits (Holtan and Hake, 2003; Schauer et al, 2006; Semel et al, 2006; Lippman and Zamir, 2007; Lippman et al, 2007). The rationale for our approach was based on the ability of sRNAs in plants to trigger secondary small interfering RNAs (siRNAs; Baulcombe, 2004). We hypothesised that primary sRNA from one parent could initiate secondary siRNA on a target RNA from the opposite parent through an RNA-based mechanism. The initiator primary sRNA would associate with an Argonaute protein and stimulate an RNA-dependent RNA polymerase (RDR) to convert a target RNA into a double-stranded (ds) form. The dsRNA would then be processed into 21-nucleotide (nt) secondary siRNAs by a Dicer-like (DCL) protein. In this mechanism, there would be silencing of multiple RNAs including the primary target RNAs and any others that are complementary to the secondary siRNAs (Baulcombe, 2004). At other loci the primary sRNA may direct methylation of a targeted DNA locus in the opposite parent. The methylated DNA would then become the template for secondary siRNA production via a mechanism that involves Polymerase IV, RDR2 and DCL3 (Wierzbicki et al, 2008; Matzke et al, 2009; Law and Jacobsen, 2010). Unlike the other secondary sRNA pathway this process would involve 24-nt siRNAs, and the silencing mechanism could be transcriptional or posttranscriptional. However, as in the posttranscriptional pathway, there would be silencing of both primary and secondary target loci and transgressive expression of sRNA. Here we describe transgressive effects at many different sRNA loci in the hybrid lines of tomato and we investigated three of them in more detail. We describe how the expression of sRNA at these loci in hybrids affects gene expression and how they could be causal of transgressive phenotypes. Two of these examples involve siRNA loci and they could be explained by RNA-mediated epistasis. However, a third example is not easily reconciled with RNA-mediated epistasis because it involves a micro (mi)RNA whose biogenesis does not involve secondary sRNAs. These examples illustrate how different sRNA-based mechanisms could be involved in transgressive segregation. Results To find out whether miRNAs and siRNAs are transgressively expressed, we cloned and sequenced sRNA from 2-week-old seedlings of parental lines (S. lycopersicum cv. M82 and S. pennellii, accession LA716), their F1 and F2 hybrids and ILs. The sRNA sequences in the resulting data sets are predominantly 21- and 24-nt in length (Supplementary Figure S1B and Materials and methods section) with a 5′-nucleotide bias of U (21-nt) or A (24-nt), respectively (Supplementary Figure S1C). This bias is likely due to selective binding by Argonaute proteins (Baulcombe, 2004; Mi et al, 2008). Among all samples, 24-nt sRNAs accumulate to higher levels than those of 21-nt sRNAs. Small RNA-enriched genomic regions (loci) were identified using SiLoCo (Moxon et al, 2008). The number of reads at each locus was normalised to the number of total reads as an indicator of the expression level (Materials and methods section). A transgressive index was calculated based on a comparison of locus representation in the different sRNA data sets. It was calculated as the ratio of the normalised counts for a genomic locus, for an IL or hybrid compared with the mean of the normalised counts of the most similar parent. Transgressive expression of miR395 associated with salt stress Among miRNAs the most striking transgressive effect was with miR395. This miRNA was abundant in S. pennellii but not in M82 seedlings (Figure 1A and B), and it was at a mid-parental level in the F1. However, in two of the three F2 samples and in a few of the ILs including IL8-1-3 (Figure 1A and B and Supplementary Figure S2A) this miR395 was hyper-abundant. This transgressive expression pattern was confirmed by northern blotting (Figure 1B and Supplementary Table SI). Figure 1.Transgressive expression of miR395. (A) Bar chart showing the abundance of miR395, normalised to the sizes of the sRNA data sets. Relative abundance was calculated using miRprof (UEA sRNA tool kit (Moxon et al, 2008)). (B) Accumulation of miR395 in parents and selected ILs. About 15 μg of total RNA from young seedlings were used for analysis of each sample. Total RNA was electrophoresed in 15% polyacrylamide gel, transferred to membrane and probed with corresponding DNA oligonucleotides (Supplementary Table SI) labelled with γ-32P-ATP. U6 serves as loading control. M, Size marker. Download figure Download PowerPoint The miRNA 395 loci in the tomato genome are on chromosomes 2, 3 and 5, and not in the introgressed region of S. pennellii in IL8-1-3. The enhanced expression of miR395 is therefore likely to have been initiated by a component of the S. pennellii genome that acts either directly or indirectly on the M82 genome. Expression of miR395 was induced by salt stress in S. pennellii (Supplementary Figure S3A), which is a naturally salt tolerant wild species (Frary et al, 2010). This miRNA is also induced by salt stress in tobacco (Frazier et al, 2011), maize (Ding et al, 2009) and poplar (Jia et al, 2009). In Arabidopsis, but not M82 or S. pennellii (Supplementary Figure S3B), miR395 is associated with sulphate stress (Jones-Rhoades and Bartel, 2004), and its target transcripts are ATP sulphurylase (APS; Supplementary Figure S3C) (Ding et al, 2009) and other mRNAs involved in both sulfur assimilation and allocation (Chiou, 2007; Kawashima et al, 2009). Similar to S. pennellii, miR395 was salt inducible in IL8-1-3, but not the other ILs tested (Figure 2A). The high-level expression of miR395 in S. pennellii and IL8-1-3 correlated with a low level of APS1 mRNA in the absence of salt. The APS1 mRNA level was further reduced in these two lines, but not in M82 or three other ILs under salt stress conditions (Figure 2B). There was a correspondingly greater tolerance to salt stress in S. pennellii and IL8-1-3 than in wild type, as indicated by reduced loss of chlorophyll in a leaf disc assay (Figure 2C) and in other assays including a measurement of leaf curling (Supplementary Figure S4). It is likely therefore that overexpression of miR395 in IL8-1-3 and S. pennellii contributes, at least in part, to their enhanced salt tolerance in comparison with M82. However, there are likely to be other determinants of salt tolerance because the miR395 target mRNA was at similar levels in IL8-1-3 and S. pennellii, although they differed in their salt tolerance. We also identified a second salt tolerant line (IL2-5) (Supplementary Figure S4A) in which miR395 levels were not elevated. Figure 2.miR395 expression contributes to a transgressive phenotype. (A) miR395 accumulation in salt treated (NaCl in mM) plants from parental lines (M82 and S. pennellii), IL8-1-3 and three other ILs. (B) Transcript analysis for APS1. RNA extraction and qRT–PCR was performed as described. GAPDH served as control. Data are mean±s.d. of three replicates. (C) Chlorophyll accumulation in salt-treated and control plants. Chlorophyll accumulation without salt treatment was normalised to 100%. Chlorophyll (μg/g fresh tissue weight) from three independent experiments (n>30) is shown. Data are mean±s.d. The double asterisks in B and C represent significance (one-way (B) or two-way (C) analysis of variance test (Holm-Sidak method; α=0.05)) at P<0.01. Download figure Download PowerPoint Transgressive siRNA loci When aligned with the available tomato genomic sequence (Materials and methods section), the M82 data set revealed 24 000 siRNA loci in which there were multiple matches to sRNA reads. The S. pennellii data set revealed 18 000 loci but, as the genome of this species is not available, the divergent sequence loci would not be identified and this number is an underestimate. More than 90% of the shared loci are represented at similar levels in the data sets from the two parental species, indicating that they are expressed at the same level in M82 and S. pennellii. To identify differentially expressed siRNA loci, we restricted our analysis to regions of the genome with at least 100 reads in at least one of the data sets. We identified 153 transgressive loci that were expressed in F2 or IL lines at levels that are either higher than the high expressing parent or lower than the low expresser (Figure 3A and Supplementary Table SII). We used a stringent criterion of transgressivity and restricted our analysis to loci that were identified consistently using three different protocols. According to this criterion, there were no sRNA loci with transgressive expression in the F1 samples but, in the ILs and F2 plants, few siRNA loci were transgressively expressed (Supplementary Table SII). At present we cannot rule out that there are sRNA loci that are transgressively expressed in the F1. However, if such loci exist, their transgressive index would be substantially less than that with the loci identified here. Figure 3.Transgressive siRNA locus H06. (A) Heatmap showing changes in the level of sRNA from the most transgressive genomic loci in IL samples. The H06 locus is marked with *. Replicates of M82 (M82-1 and M82-2) and S. pennellii (S. penn-1 and S. penn-2) were used. Details of TI for each siRNA loci presented in the heatmap are given in Supplementary Table SII. (B) sRNAs from M82, S. pennellii, F1, F2 and introgression line samples aligned to the H06 genomic locus in M82 DNA. Arrows indicate direction of sRNA alignment, and the colour of the arrow indicates length: blue, 24-nt; red, 21-nt; and green, 22-nt. (C) Northern blot analysis of H06 siRNAs in parental lines, F1, F2 and in ILs in four different chromosomes. M, size marker. Oligonucleotide probe sequences (H06a and b) are given in Supplementary Table SI. Download figure Download PowerPoint Two types of transgressive loci are present in these data sets: one type was less abundant in the ILs than in the parental lines (lower half, Figure 3A) and the second type was more abundant in ILs than in the parents (upper half, Figure 3A). The transgressive expression, depending on the locus, was evident either in a subset of ILs or in all of them. Fewer than 1% of all loci were transgressive and they were distributed all over the genome similar to the total sRNA loci (Supplementary Figure S5A). Transgressive loci were also of similar length as any other sRNA loci (Supplementary Figure S5B). The transgressive loci overlap transposons, coding sequence and other genomic features (Supplementary Figure S6) but, with only 153 loci involved, we cannot determine whether the pattern is significantly different from total sRNA loci. The transgressive sRNAs were distinct from the total sRNAs in that they had a greater strand bias than the other sRNAs (Supplementary Figure S7A). They produce 21, 22 and 24-nt sRNAs, although the 24-nt sRNAs were more highly represented than in the total sRNA population (Supplementary Figure S7B). The relationship of these properties and mechanisms is not obvious, but they are consistent with the involvement of multiple biogenesis pathways in transgressive sRNA production. A transgressive siRNA locus associated with hypermethylation of its target DNA One of the 153 transgressive loci was overexpressed relative to the parents in all of the ILs (marked as * in Figure 3A). This locus is referred to as H06 and it overlaps a predicted protein coding gene on chromosome 8 (mRNA:Co.8.11-contig28_33.1.1; Figure 3B). There was only a single H06 siRNA read in the data sets from the S. pennellii parent and none from M82 tomato or the F1 hybrid, but there were hundreds of siRNA reads in each of nine IL data sets (Figure 3B and Supplementary Figure S8). Two F2 lines expressed H06 siRNAs, albeit at lower levels than the ILs (Figure 3B and Supplementary Figure S2B). The activation of H06 siRNA could not be correlated with a particular region of S. pennellii DNA because the lines tested have introgressed DNA from four different chromosomes (1, 2, 8 and 4, Figure 3C). The initiation event was therefore likely to have been in a progenitor of the ILs. This event would have acted in trans to affect the H06 locus on chromosome 8 of the M82 genome. The transgressive F2 lines in this study (Figure 3B) would have recapitulated this initiation event. The H06 siRNAs were predominantly 24-nt in length and their overexpression in the ILs relative to the parental or F1 hybrid lines was confirmed by northern analysis (Figure 3B and C). The 24-nt size class of sRNA is associated with methylation of DNA with a similar sequence (Hamilton et al, 2002; Chan et al, 2004) and we therefore tested for methyl cytosine at the H06 locus using an McrBC methylation assay (Figure 4A). The H06 DNA from the parental lines and F1 was largely undigested in this assay indicating hypomethylation whereas, in the F2 and IL samples, the DNA was completely digested by McrBC indicative of hypermethylation. Using a sodium bisulphite analysis we confirmed that genomic DNA from parents and an F1 plant was largely unmethylated at this locus, but it was densely methylated in the F2-3 line and the two ILs that were analysed (Figure 4B and Supplementary Figure S9) at CG, CHG and CHH contexts. Figure 4.DNA methylation and expression analysis of H06 locus. (A) Methylation of H06 genomic DNA using restriction enzyme McrBC. Completely digested DNA indicates hypermethylation. Error bars indicate the 95% confidence intervals of three replicates. The assay was performed as described in Materials and methods section. (B) Sodium bisulphite analysis of H06 locus. Each sample had at least 10 clones analysed for cytosine methylation in replicate experiments. (C) Expression analysis of H06 mRNA using quantitative RT–PCR. Download figure Download PowerPoint The H06 mRNA fragment detected by quantitative RT–PCR was less abundant in the ILs than in the parents or their F1 hybrid (Figure 4C), indicating that epigenetic modification of H06 DNA, siRNA production and silencing of H06 mRNA are interrelated phenomena. Transgressive siRNA from the PAL loci are not associated with DNA methylation In the absence of a fully sequenced and annotated tomato genome, we also searched for transgressive siRNA loci by aligning the siRNA sequence reads with expressed sequence tag cDNAs (ESTs). The transgressive expression of cDNA-specific siRNAs was more pronounced in ILs than in the F1 and F2 loci (Figure 5A and Supplementary Tables SIII and SIV). Figure 5.Transgressive PAL siRNA accumulate in IL8-3. (A) Heatmap showing changes in the level of sRNA from the most transgressive cDNA loci in the IL samples. The major families of transgressive siRNAs correspond to the following genetic loci: PAL, phenylalanine ammonia–lyase; EMB, Embryo defective; HB5,Homeobox protein family; VTC, 'Vitamin C deficient' family; WIP, Wound inducible protein family. Transgressive siRNA loci are listed in Supplementary Table SIII. (B) sRNAs aligning to the PAL locus in M82, S. pennellii, F1, F2 and ILs. Arrows indicate direction of sRNA alignment, and the colour of the arrow indicates length: blue, 24-nt; red, 21-nt; and green, 22-nt. (C) Northern blot analysis of small RNAs hybridising to a PAL-specific probe. Download figure Download PowerPoint The strongest transgressive effect of this type was in IL8-3, in which siRNAs corresponding to phenylalanine ammonia-lyase (PAL) mRNAs were increased by more than 100-fold relative to the parents or other ILs (Figure 5A and B and Supplementary Figure S2C). The PAL siRNAs in IL8-3 were predominantly 21-nt and their overexpression in IL8-3 was confirmed by northern analysis (Figure 5C). PAL siRNAs mostly aligned to the exons and matched predominantly to sense strands in the parents and most of the ILs. However, the transgressive siRNAs in IL8-3 matched to both sense and antisense strands (Supplementary Figure S10), indicating that they originated from a dsRNA molecule. The PAL-derived sRNAs in IL8-3 matched numerous PAL ESTs (Supplementary Table SIII and SIV). A region in chromosome 9 (CO9SLe0130H12.3) has four PAL gene sequences (named as PAL5A, PAL5B, PAL5C and PAL5D) configured as a divergent inverted repeat (Supplementary Figure S10). These genes are 90–99% identical and also very similar (93–99%) to a PAL5 gene, previously identified from tomato that is expressed abundantly (Chang et al, 2008) (Supplementary Figure S11). There were unique siRNA sequence reads corresponding to each of the four PAL5 genes (Supplementary Table SV), indicating that the corresponding mRNAs are all templates for siRNA production. PAL5B had the highest proportion of unique sequence matches and it is likely to be the most abundant source of siRNAs. The alignment of PAL siRNAs with the PAL sequences was phased (Figure 6A) in a manner that is characteristic of trans-acting si (tasi)RNA in Arabidopsis and other plants (Peragine et al, 2004; Vazquez et al, 2004). An initiation event associated with these transgressive PAL siRNAs would have occurred in a progenitor of IL8-3 and it would have acted in trans to affect siRNA production at this chromosome 9 locus. A plausible scenario is that the initiator is an sRNA that is functionally equivalent to the miRNAs in tasiRNA biogenesis. Figure 6.PAL silencing results in reduced lignification. (A) Phasing of PAL-derived siRNAs along the PAL sequences. Two PAL ESTs with significant P-values (PAL5A, log P value=−9.06; PAL5C, log P-value=−8.09) for phasing are shown. Red bars indicate the start position of sRNAs in phase; and blue, the start position for those out of phase. Rectangles indicate expected phased positions and those in red indicate phased sequences. (B) PAL mRNA expression of parental and hybrid lines. See Materials and methods section for details. Data are mean±s.d. of three biological replicates. The double asterisks represent significance determined by one-way ANOVA at P<0.01 (α=0.05). (C) Phloroglucinol staining for total lignin in parental lines, IL8-3 and three other ILs. Violet stain indicates the presence and distribution of lignin. The micrographs are of sections derived from petioles of seedlings with staining concentrated in the vascular bundles. Size bar, 100 μM. Download figure Download PowerPoint The level of PAL mRNA was assessed by RT–PCR using primers that would have detected all transcripts in the PAL5. It corresponded to the level of siRNAs in that it was lower in the IL8-3 line than in the parents, F1 or F2 hybrids or other ILs (Figure 6B). PAL initiates the phenylpropanoid pathway (Jones, 1984) that is associated with lignification of xylem vessels and, in line with the mRNA levels, the total lignin was lower in IL8-3 (and confined to fewer xylem vessels) than in parents or the other ILs tested (Figure 6C). Discussion Our analysis of tomato hybrids was prompted by the hypothesis that transgressive segregation would be affected by epistatic interactions between siRNAs and their targets from the opposite parent. Such interactions would be specific to the hybrids and the effect would persist in lines that inherit both interacting loci. In the most straightforward iteration of this hypothesis, the transgressive expression of sRNA would be fixed in lines in which the loci are both homozygous. The transgressive H06 and PAL siRNA loci are partly consistent with this hypothesis, because they are located outside the introgressed region of S. pennellii DNA, and their overexpression in the hybrids could have been induced by an epistatic interaction, involving primary sRNAs from S. pennellii DNA interacting with targets in the tomato genome. However, simple epistasis can be ruled out because the transgressive expression of these loci was absent in the F1 when all the pairs of putatively interacting loci would have been present. It seems that there must be secondary effects following the initiating interaction. A formal possibility that the primary events involve mutation in male or female gametophytes seems unlikely, at least for H06, because the transgressive effect could be recapitulated in newly created F2 lines (Figures 1 and 3, and Supplementary Figures S2, S8). Genetic rearrangements are presumed to be rare and irregular and would be difficult to recreate. An alternative explanation for H06 invokes epigenetic marks at H06 that are induced in the F1. These marks would then trigger sRNA expression in the F2 and in subsequent generations that inherit the epigenetically modified H06 DNA. To account for the absence of an effect in the F1, we propose that the epigenetic effects were initiated by interactions occurring during gametogenesis of the F1 and that they were subsequently reinforced by RNA-directed DNA methylation. Consistent with that idea, the transgressive effect of H06 was associated with dense methylation of H06 DNA in CHH as well as CG and CHG (Figure 4A and B; Chan et al, 2004; Teixeira and Colot, 2010). The involvement of the 24-nt class of siRNA (Figure 3B and C) that is associated with RNA-directed DNA methylation and histone modification (Chan et al, 2004; Xie et al, 2004; Matzke and Birchler, 2005) is also consistent with RNA-directed mechanism at H06. An epigenetic effect could also explain the progressive increase of H06 over several generations that is implied by the comparison of F1, F2 and IL lines. Epigenetic effects at various transgene and endogenous siRNA loci similarly increase through several generations (Shiu et al, 2001; Volpe et al, 2002; Teixeira et al, 2009). A clue as to the likely mechanism at the PAL locus is from the similarity of PAL siRNAs in IL8-3 to endogenous ta-siRNAs (Peragine et al, 2004; Vazquez et al, 2004). Both are predominantly 21-nt, derived from both strands of a dsRNA precursor, and their alignment with their target nucleotide sequence is phased (Figure 6A). In addition, the target DNA sequences of PAL siRNAs are not methylated in IL8-3, suggesting a posttranscriptional mechanism (Supplementary Figure S12). The production of ta-siRNAs is initiated by miRNAs that target a non-coding ta-siRNA precursor (Allen et al, 2005). By analogy, in IL8-3, it is likely that there is an initiator sRNA of PAL siRNA biogenesis and that the PAL mRNA is equivalent to the ta-siRNA precursor. Presumably there was a genetic or an epigenetic change in one of the precursors of IL8-3 that activated production of this sRNA initiator of PAL silencing. The transgressive
Referência(s)