Small RNA-directed epigenetic natural variation in Arabidopsis thaliana

Progress in epigenetics has revealed mechanisms that can heritably regulate gene function independent of genetic alterations. Nevertheless, little is known about the role of epigenetics in evolution. This is due in part to scant data on epigenetic variation among natural populations. In plants, small interfering RNA (siRNA) is involved in both the initiation and maintenance of gene silencing by directing DNA methylation and/or histone methylation. Here, we report that, in the model plant Arabidopsis thaliana, a cluster of approximately 24 nt siRNAs found at high levels in the ecotype Landsberg erecta (Ler) could direct DNA methylation and heterochromatinization at a hAT element adjacent to the promoter of FLOWERING LOCUS C (FLC), a major repressor of flowering, whereas the same hAT element in ecotype Columbia (Col) with almost identical DNA sequence, generates a set of low abundance siRNAs that do not direct these activities. We have called this hAT element MPF for Methylated region near Promoter of FLC, although de novo methylation triggered by an inverted repeat transgene at this region in Col does not alter its FLC expression. DNA methylation of the Ler allele MPF is dependent on genes in known silencing pathways, and such methylation is transmissible to Col by genetic crosses, although with varying degrees of penetrance. A genome-wide comparison of Ler and Col small RNAs identified at least 68 loci matched by a significant level of approximately 24 nt siRNAs present specifically in Ler but not Col, where nearly half of the loci are related to repeat or TE sequences. Methylation analysis revealed that 88% of the examined loci (37 out of 42) were specifically methylated in Ler but not Col, suggesting that small RNA can direct epigenetic differences between two closely related Arabidopsis ecotypes.     

Natural allelic variation identifies new genes in the Arabidopsis circadian system

We have analysed the circadian rhythm of Arabidopsis thaliana leaf movements in the accession Cvi from the Cape Verde Islands, and in the commonly used laboratory strains Columbia (Col) and Landsberg (erecta) (Ler), which originated in Northern Europe. The parental lines have similar rhythmic periods, but the progeny of crosses among them reveal extensive variation for this trait. An analysis of 48 Ler/Cvi recombinant inbred lines (RILs) and a further 30 Ler/Col RILs allowed us to locate four putative quantitative trait loci (QTLs) that control the period of the circadian clock. Near-isogenic lines (NILs) that contain a QTL in a small, defined chromo- somal region allowed us to confirm the phenotypic effect and to map the positions of three period QTLs, designated ESPRESSO, NON TROPPO and RALENTANDO. Quantitative trait loci at the locations of RALENTANDO and of a fourth QTL, ANDANTE, were identified in both Ler/Cvi and Ler/Col RIL populations. Some QTLs for circadian period are closely linked to loci that control flowering time, including FLC. We show that flc mutations shorten the circadian period such that the known allelic variation in the MADS-box gene FLC can account for the ANDANTE QTL. The QTLs ESPRESSO and RALENTANDO identify new genes that regulate the Arabidopsis circadian system in nature, one of which may be the flowering-time gene GIGANTEA.     

Global analysis of genetic, epigenetic and transcriptional polymorphisms in Arabidopsis thaliana using whole genome tiling arrays

Whole genome tiling arrays provide a high resolution platform for profiling of genetic, epigenetic, and gene expression polymorphisms. In this study we surveyed natural genomic variation in cytosine methylation among Arabidopsis thaliana wild accessions Columbia (Col) and Vancouver (Van) by comparing hybridization intensity difference between genomic DNA digested with either methylation-sensitive (HpaII) or -insensitive (MspI) restriction enzyme. Single Feature Polymorphisms (SFPs) were assayed on a full set of 1,683,620 unique features of Arabidopsis Tiling Array 1.0F (Affymetrix), while constitutive and polymorphic CG methylation were assayed on a subset of 54,519 features, which contain a 5'CCGG3' restriction site. 138,552 SFPs (1% FDR) were identified across enzyme treatments, which preferentially accumulated in pericentromeric regions. Our study also demonstrates that at least 8% of all analyzed CCGG sites were constitutively methylated across the two strains, while about 10% of all analyzed CCGG sites were differentially methylated between the two strains. Within euchromatin arms, both constitutive and polymorphic CG methylation accumulated in central regions of genes but under-represented toward the 5' and 3' ends of the coding sequences. Nevertheless, polymorphic methylation occurred much more frequently in gene ends than constitutive methylation. Inheritance of methylation polymorphisms in reciprocal F1 hybrids was predominantly additive, with F1 plants generally showing levels of methylation intermediate between the parents. By comparing gene expression profiles, using matched tissue samples, we found that magnitude of methylation polymorphism immediately upstream or downstream of the gene was inversely correlated with the degree of expression variation for that gene. In contrast, methylation polymorphism within genic region showed weak positive correlation with expression variation. Our results demonstrated extensive genetic and epigenetic polymorphisms between Arabidopsis accessions and suggested a possible relationship between natural CG methylation variation and gene expression variation.     

FRIGIDA LIKE 2 is a functional allele in Landsberg erecta and compensates for a nonsense allele of FRIGIDA LIKE 1

The Landsberg erecta (Ler) accession of Arabidopsis (Arabidopsis thaliana) has a weak allele of the floral inhibitor FLOWERING LOCUS C (FLC). FLC-Ler is weakly up-regulated by the active San Feliu-2 (Sf2) allele of FRIGIDA (FRI-Sf2), resulting in a moderately late-flowering phenotype. By contrast, the Columbia (Col) allele of FLC is strongly up-regulated by FRI-Sf2, resulting in a very late-flowering phenotype. In Col, the FRI-related gene FRI LIKE 1 (FRL1) is required for FRI-mediated up-regulation of FLC. It is shown here that in Ler, the FRL1-related gene FRI LIKE 2 (FRL2), but not FRL1, is required for FRI-mediated up-regulation of FLC. FRL1-Ler is shown to be a nonsense allele of FRL1 due to a naturally occurring premature stop codon in the middle of the conceptual protein sequence, suggesting that FRL1-Ler is nonfunctional. Compared to FRL2-Col, FRL2-Ler has two amino acid changes in the conceptual protein sequence. Plants homozygous for FRI-Sf2, FLC-Ler, FRL1-Ler, and FRL2-Col have no detectable FLC expression, resulting in an extremely early flowering phenotype. Transformation of a genomic fragment of FRL2-Ler, but not of FRL2-Col, into a recombinant inbred line derived from these plants restores both FRI-mediated up-regulation of FLC expression and a late-flowering phenotype, indicating that FRL2-Ler is the functional allele of FRL2. Taken together, these results suggest that in the two different Arabidopsis accessions Col and Ler, either FRL1 or FRL2, but not both, is functional and required for FRI-mediated up-regulation of FLC.     

Quantitative trait loci for inflorescence development in Arabidopsis thaliana

Variation in inflorescence development patterns is a central factor in the evolutionary ecology of plants. The genetic architectures of 13 traits associated with inflorescence developmental timing, architecture, rosette morphology, and fitness were investigated in Arabidopsis thaliana, a model plant system. There is substantial naturally occurring genetic variation for inflorescence development traits, with broad sense heritabilities computed from 21 Arabidopsis ecotypes ranging from 0.134 to 0.772. Genetic correlations are significant for most (64/78) pairs of traits, suggesting either pleiotropy or tight linkage among loci. Quantitative trait locus (QTL) mapping indicates 47 and 63 QTL for inflorescence developmental traits in Ler x Col and Cvi x Ler recombinant inbred mapping populations, respectively. Several QTL associated with different developmental traits map to the same Arabidopsis chromosomal regions, in agreement with the strong genetic correlations observed. Epistasis among QTL was observed only in the Cvi x Ler population, and only between regions on chromosomes 1 and 5. Examination of the completed Arabidopsis genome sequence in three QTL regions revealed between 375 and 783 genes per region. Previously identified flowering time, inflorescence architecture, floral meristem identity, and hormone signaling genes represent some of the many candidate genes in these regions.     

Quantitative trait loci analysis of nitrate storage in Arabidopsis leading to an investigation of the contribution of the anion channel gene, AtCLC-c, to variation in nitrate levels

Storage of excess nitrate in the vacuole and its subsequent remobilization is an important aspect of a plant's nitrogen economy, but the genes controlling the underlying processes have not all been identified and characterized. Cape Verdi Island (Cvi)/Landsberg erecta (Ler) and Columbia (Col)/Landsberg erecta recombinant inbred line (RIL) populations of Arabidopsis thaliana were used to identify quantitative trait loci (QTL) controlling natural variation in nitrate concentrations. One major and two minor QTLs were found for the Cvi/Ler population and one minor QTL for the Col/Ler RIL. These were designated NA1 to NA4. The major Cvi/Ler QTL (NA3) was located at the bottom of chromosome 5. No interaction among the QTLs was found by two-way ANOVA. By comparing in silico the locations of the QTLs with a physical map of the Arabidopsis genome, candidate genes for each QTL were identified. Several of these were anion channels of the AtCLC family. One of these, AtCLC-c, coincided with NA3 and its role was investigated using a mutant with a transposon insertion in AtCLC-c. Mutant plants homozygous for the insertion (designated clcc-1) had less than 5% of AtCLC-c mRNA compared with wild-type (WT) shoots. They also had significantly lower nitrate concentrations when grown at a range of external nitrate concentrations. The concentrations of chloride, malate, and citrate were also affected in the mutant. In wild-type plants, expression of AtCLC-c was down-regulated in the presence of nitrate, but ammonium had a much smaller effect while chloride and sulphate did not affect expression. These and published results suggest that multiple genes affect nitrate concentrations in plants and that AtCLC-c and other members of the AtCLC gene family play some role in this.     

The Cape Verde Islands allele of cryptochrome 2 enhances cotyledon unfolding in the absence of blue light in Arabidopsis

We analyzed the natural genetic variation between Landsburg erecta (Ler) and Cape Verde Islands (Cvi) accessions by studying 105 recombinant inbred lines to search for players in the regulation of sensitivity to light signals perceived by phytochromes in etiolated seedlings of Arabidopsis. In seedlings grown under hourly pulses of far-red (FR) light, we identified three quantitative trait loci (QTLs; VLF3, VLF4, and VLF5) for hypocotyl growth inhibition and three different QTLs (VLF6, VLF7, and VLF1) for cotyledon unfolding. This indicates that different physiological outputs have selective regulation of sensitivity during de-etiolation. Ler alleles, compared with Cvi alleles, of VLF3, VLF4, VLF5, VLF7, and VLF1 enhanced, whereas the Ler allele of VLF6 reduced, the response to pulses of FR. We confirmed and narrowed down the position of some QTLs by using near-isogenic lines. VLF6 mapped close to the CRY2 (cryptochrome 2) gene. Transgenic Ler seedlings expressing the Cvi allele of CRY2 showed enhanced cotyledon unfolding under hourly pulses of FR compared with the wild type or transgenics expressing the CRY2-Ler allele. This response required phytochrome A. The cry1 cry2 double mutant lacking both cryptochromes showed reduced cotyledon unfolding under FR pulses. Because the CRY2-Cvi is a gain-of-function allele compared with CRY2-Ler, cryptochrome activity correlates positively with cotyledon unfolding under FR pulses. We conclude that the blue light photoreceptor cryptochrome 2 can modulate seedling photomorphogenesis in the absence of blue light. In addition to the nuclear loci, we identified cytoplasmic effects on seedling de-etiolation.     

Molecular marker genotypes, heterozygosity and genetic interactions explain heterosis in Arabidopsis thaliana

The underlying mechanisms for hybrid vigor or heterosis are elusive. Here we report a population of recombinant inbred lines (RILs), derived from the two ecotypes, Col and Ler, which can serve as a permanent resource for studying the molecular basis of hybrid vigor in Arabidopsis. Using a North Carolina mating design III (NCIII), we determined the additive and dominant nature of gene action in this population. We detected heterosis among crosses of RILs with one of the two parents (Col and Ler) and analyzed genotypes and heterozygosities for RILs and test cross families (RILs crossed to Col and Ler) using a total of 446 published molecular markers. The performance of RILs and additive and dominant components in the test cross families were used to analyze QTLs for 16 traits, using QTL cartographer and composite interval mapping with 1,000 permutations for each trait. Our data suggest that locus-specific and/or genome-wide differential heterozygosity, including epistasis, plays an important role in the generation of the observed heterosis. Furthermore, the hybrid vigor occurred between two closely related ecotypes, and provides a general mechanism for novel variation generated between genetically similar materials.     

Pleiotropic effects of the Arabidopsis cryptochrome 2 allelic variation underlie fruit trait-related QTL

The previous molecular identification of a flowering time QTL segregating in the Arabidopsis L er x Cvi cross, demonstrated that natural allelic variation at the blue light photoreceptor CRY2 gene affects flowering time (El-Assal et al., 2001). In addition, previous works on the same cross have mapped several QTL affecting other unrelated life history traits in the CRY2 genomic region. In the present report, we have used a set of Arabidopsis L er transgenic plants carrying four different functional CRY2 transgenes for phenotypic analyses, with the aim of exploring the extent of pleiotropy of CRY2 allelic variation. It is concluded that previously identified QTL affecting fruit length, ovule number per fruit, and percentage of unfertilized ovules are caused by this same Ler/Cvi CRY2 allelic variation. In addition, dose effects of the CRY2-L er allele are detected for fruit length. A seed weight QTL at the map position of CRY2 could not be confirmed and also no effect on seed dormancy was observed. Thus, it is shown that transgenic plants carrying different alleles can be a useful tool to attribute QTL for different complex traits to a specific locus, even when the relationship among the traits has not been previously suggested.     

Identification of quantitative trait loci controlling symptom development during viral infection in Arabidopsis thaliana

In compatible interactions between plants and viruses that result in systemic infection, symptom development is a major phenotypic trait. However, host determinants governing this trait are mostly unknown, and the mechanisms underlying it are still poorly understood. In a previous study on the Arabidopsis thaliana-Plum pox virus (PPV) pathosystem, we showed a large degree of variation in symptom development among susceptible accessions. In particular, Cvi-1 (Cape Verde islands) accumulates viral particules but remains symptomless, Col-0 (Columbia) sometimes shows weak symptoms compared with Ler (Landsberg erecta), which always shows severe symptoms. Genetic analyses of Col x Ler and Cvi x Ler F2 and recombinant inbred line (RIL) populations suggested that symptom development as well as viral accumulation traits are polygenic and quantitative. Three of the symptom quantitative trait loci (QTL) identified could be confirmed in near-isogenic lines, including PSI1 (PPV symptom induction 1), which was identified on the distal part of chromosome 1 in both RIL populations. With respect to viral accumulation, several factors have been detected and, interestingly, in the Col x Ler population, two out of three viral accumulation QTL colocalized with loci controlling symptom development, although correlation analysis showed weak linearity between symptom severity and virus accumulation. In addition, in the Cvi x Ler RIL population, a digenic recessive determinant controlling PPV infection was identified.     

利用拟南芥杂交组合研究siRNA与等位基因DNA甲基化调控中的联系

近年来,越来越多的实验结果表明,表观遗传因子,如DNA甲基化、小RNA、组蛋白修饰等在杂种优势形成中起到重要作用,然而对于这些表观遗传因子在F。中遗传调控方式的认识仍很有限．本实验室先前工作曾以拟南芥C24和Ler两种生态型及其正反交子一代为材料,运用新一代测序方法获得该杂交组合中DNA甲基化及小RNA单碱基分辨率的全基因组图谱．本文进一步对这批数据中的等位基因DNA甲基化水平进行分析,区分DNA甲基化遗传过程中的顺式与反式调控方式,并发现这两种调控方式均有重要的贡献．研究发现,siRNA与DNA甲基化的两种调控方式有密切联系,尤其在DNA甲基化的反式调控中,Fl中DNA甲基化变化程度越大,该区域内siRNA富集程度越强,二者可能存在某种调控机制．通过等位基因表观遗传组的分析研究杂交过程中DNA甲基化和小RNA遗传调控的规律,为更好地理解杂种优势机制提供了帮助．