Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing.

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

On the causes of gene-body methylation variation in Arabidopsis thaliana.

Gene-body methylation (gbM) refers to sparse CG methylation of coding regions, which is especially prominent in evolutionarily conserved house-keeping genes. It is found in both plants and animals, but is directly and stably (epigenetically) inherited over multiple generations in the former. Studies in Arabidopsis thaliana have demonstrated that plants originating from different parts of the world exhibit genome-wide differences in gbM, which could reflect direct selection on gbM, but which could also reflect an epigenetic memory of ancestral genetic and/or environmental factors. Here we look for evidence of such factors in F2 plants resulting from a cross between a southern Swedish line with low gbM and a northern Swedish line with high gbM, grown at two different temperatures. Using bisulfite-sequencing data with nucleotide-level resolution on hundreds of individuals, we confirm that CG sites are either methylated (nearly 100% methylation across sampled cells) or unmethylated (approximately 0% methylation across sampled cells), and show that the higher level of gbM in the northern line is due to more sites being methylated. Furthermore, methylation variants almost always show Mendelian segregation, consistent with their being directly and stably inherited through meiosis. To explore how the differences between the parental lines could have arisen, we focused on somatic deviations from the inherited state, distinguishing between gains (relative to the inherited 0% methylation) and losses (relative to the inherited 100% methylation) at each site in the F2 generation. We demonstrate that deviations predominantly affect sites that differ between the parental lines, consistent with these sites being more mutable. Gains and losses behave very differently in terms of the genomic distribution, and are influenced by the local chromatin state. We find clear evidence for different trans-acting genetic polymorphism affecting gains and losses, with those affecting gains showing strong environmental interactions (G×E). Direct effects of the environment were minimal. In conclusion, we show that genetic and environmental factors can change gbM at a cellular level, and hypothesize that these factors can also lead to transgenerational differences between individuals via the inclusion of such changes in the zygote. If true, this could explain genographic pattern of gbM with selection, and would cast doubt on estimates of epimutation rates from inbred lines in constant environments.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes.

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

Cell cycle status of male and female gametes during Arabidopsis reproduction.

Fertilization in Arabidopsis (Arabidopsis thaliana) is a highly coordinated process that begins with a pollen tube delivering the 2 sperm cells into the embryo sac. Each sperm cell can then fertilize either the egg or the central cell to initiate embryo or endosperm development, respectively. The success of this double fertilization process requires a tight cell cycle synchrony between the male and female gametes to allow karyogamy (nuclei fusion). However, the cell cycle status of the male and female gametes during fertilization remains elusive as DNA quantification and DNA replication assays have given conflicting results. Here, to reconcile these results, we quantified the DNA replication state by DNA sequencing and performed microscopic analyses of fluorescent markers covering all phases of the cell cycle. We show that male and female Arabidopsis gametes are both arrested prior to DNA replication at maturity and initiate their DNA replication only during fertilization.

Germline replications and somatic mutation accumulation are independent of vegetative life span in Arabidopsis.

In plants, gametogenesis occurs late in development, and somatic mutations can therefore be transmitted to the next generation. Longer periods of growth are believed to result in an increase in the number of cell divisions before gametogenesis, with a concomitant increase in mutations arising due to replication errors. However, there is little experimental evidence addressing how many cell divisions occur before gametogenesis. Here, we measured loss of telomeric DNA and accumulation of replication errors in Arabidopsis with short and long life spans to determine the number of replications in lineages leading to gametes. Surprisingly, the number of cell divisions within the gamete lineage is nearly independent of both life span and vegetative growth. One consequence of the relatively stable number of replications per generation is that older plants may not pass along more somatically acquired mutations to their offspring. We confirmed this hypothesis by genomic sequencing of progeny from young and old plants. This independence can be achieved by hierarchical arrangement of cell divisions in plant meristems where vegetative growth is primarily accomplished by expansion of cells in rapidly dividing meristematic zones, which are only rarely refreshed by occasional divisions of more quiescent cells. We support this model by 5-ethynyl-2'-deoxyuridine retention experiments in shoot and root apical meristems. These results suggest that stem-cell organization has independently evolved in plants and animals to minimize mutations by limiting DNA replication.

Genome Sequencing Reveals the Origin of the Allotetraploid Arabidopsis suecica.

Polyploidy is an example of instantaneous speciation when it involves the formation of a new cytotype that is incompatible with the parental species. Because new polyploid individuals are likely to be rare, establishment of a new species is unlikely unless polyploids are able to reproduce through self-fertilization (selfing), or asexually. Conversely, selfing (or asexuality) makes it possible for polyploid species to originate from a single individual-a bona fide speciation event. The extent to which this happens is not known. Here, we consider the origin of Arabidopsis suecica, a selfing allopolyploid between Arabidopsis thaliana and Arabidopsis arenosa, which has hitherto been considered to be an example of a unique origin. Based on whole-genome re-sequencing of 15 natural A. suecica accessions, we identify ubiquitous shared polymorphism with the parental species, and hence conclusively reject a unique origin in favor of multiple founding individuals. We further estimate that the species originated after the last glacial maximum in Eastern Europe or central Eurasia (rather than Sweden, as the name might suggest). Finally, annotation of the self-incompatibility loci in A. suecica revealed that both loci carry non-functional alleles. The locus inherited from the selfing A. thaliana is fixed for an ancestral non-functional allele, whereas the locus inherited from the outcrossing A. arenosa is fixed for a novel loss-of-function allele. Furthermore, the allele inherited from A. thaliana is predicted to transcriptionally silence the allele inherited from A. arenosa, suggesting that loss of self-incompatibility may have been instantaneous.

GMI1, a structural-maintenance-of-chromosomes-hinge domain-containing protein, is involved in somatic homologous recombination in Arabidopsis.

DNA double-strand breaks (DSBs) pose one of the most severe threats to genome integrity, potentially leading to cell death. After detection of a DSB, the DNA damage and repair response is initiated and the DSB is repaired by non-homologous end joining and/or homologous recombination. Many components of these processes are still unknown in Arabidopsis thaliana. In this work, we characterized γ-irradiation and mitomycin C induced 1 (GMI1), a member of the SMC-hinge domain-containing protein family. RT-PCR analysis and promoter-GUS fusion studies showed that γ-irradiation, the radio-mimetic drug bleocin, and the DNA cross-linking agent mitomycin C strongly enhance GMI1 expression particularly in meristematic tissues. The induction of GMI1 by γ-irradiation depends on the signalling kinase Ataxia telangiectasia-mutated (ATM) but not on ATM and Rad3-related (ATR). Epistasis analysis of single and double mutants demonstrated that ATM acts upstream of GMI1 while the atr gmi1-2 double mutant was more sensitive than the respective single mutants. Comet assay revealed a reduced rate of DNA double-strand break repair in gmi1 mutants during the early recovery phase after exposure to bleocin. Moreover, the rate of homologous recombination of a reporter construct was strongly reduced in gmi1 mutant plants upon exposure to bleocin or mitomycin C. GMI1 is the first member of its protein family known to be involved in DNA repair.

Locally adaptive temperature response of vegetative growth in Arabidopsis thaliana.

We investigated early vegetative growth of natural Arabidopsis thaliana accessions in cold, nonfreezing temperatures, similar to temperatures these plants naturally encounter in fall at northern latitudes. We found that accessions from northern latitudes produced larger seedlings than accessions from southern latitudes, partly as a result of larger seed size. However, their subsequent vegetative growth when exposed to colder temperatures was slower. The difference was too large to be explained by random population differentiation, and is thus suggestive of local adaptation, a notion that is further supported by substantial transcriptome and metabolome changes in northern accessions. We hypothesize that the reduced growth of northern accessions is an adaptive response and a consequence of reallocating resources toward cold acclimation and winter survival.

The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history.

The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea 'Goldsmith', complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analyses reveal extensive allele sharing among species and demonstrate that introgression and selection played a role in the Aquilegia radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differs from that of the rest of the genome - a phenomenon that we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

Unstable Inheritance of 45S rRNA Genes in Arabidopsis thaliana.

The considerable genome size variation in Arabidopsis thaliana has been shown largely to be due to copy number variation (CNV) in 45S ribosomal RNA (rRNA) genes. Surprisingly, attempts to map this variation by means of genome-wide association studies (GWAS) failed to identify either of the two likely sources, namely the nucleolus organizer regions (NORs). Instead, GWAS implicated a trans-acting locus, as if rRNA gene CNV was a phenotype rather than a genotype. To explain these results, we investigated the inheritance and stability of rRNA gene copy number using the variety of genetic resources available in A. thaliana - F2 crosses, recombinant inbred lines, the multiparent advanced-generation inter-cross population, and mutation accumulation lines. Our results clearly show that rRNA gene CNV can be mapped to the NORs themselves, with both loci contributing equally to the variation. However, NOR size is unstably inherited, and dramatic copy number changes are visible already within tens of generations, which explains why it is not possible to map the NORs using GWAS. We did not find any evidence of trans-acting loci in crosses, which is also expected since changes due to such loci would take very many generations to manifest themselves. rRNA gene copy number is thus an interesting example of "missing heritability"-a trait that is heritable in pedigrees, but not in the general population.