Highly specific gene silencing by artificial microRNAs in Arabidopsis

Plant microRNAs (miRNAs) affect only a small number of targets with high sequence complementarity, while animal miRNAs usually have hundreds of targets with limited complementarity. We used artificial miRNAs (amiRNAs) to determine whether the narrow action spectrum of natural plant miRNAs reflects only intrinsic properties of the plant miRNA machinery or whether it is also due to past selection against natural miRNAs with broader specificity. amiRNAs were designed to target individual genes or groups of endogenous genes. Like natural miRNAs, they had varying numbers of target mismatches. Previously determined parameters of target selection for natural miRNAs could accurately predict direct targets of amiRNAs. The specificity of amiRNAs, as deduced from genome-wide expression profiling, was as high as that of natural plant miRNAs, supporting the notion that extensive base pairing with targets is required for plant miRNA function. amiRNAs make an effective tool for specific gene silencing in plants, especially when several related, but not identical, target genes need to be downregulated. We demonstrate that amiRNAs are also active when expressed under tissue-specific or inducible promoters, with limited nonautonomous effects. The design principles for amiRNAs have been generalized and integrated into a Web-based tool (http://wmd.weigelworld.org).

The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana

To take complete advantage of information on within-species polymorphism and divergence from close relatives, one needs to know the rate and the molecular spectrum of spontaneous mutations. To this end, we have searched for de novo spontaneous mutations in the complete nuclear genomes of five Arabidopsis thaliana mutation accumulation lines that had been maintained by single-seed descent for 30 generations. We identified and validated 99 base substitutions and 17 small and large insertions and deletions. Our results imply a spontaneous mutation rate of 7 x 10(-9) base substitutions per site per generation, the majority of which are G:C-->A:T transitions. We explain this very biased spectrum of base substitution mutations as a result of two main processes: deamination of methylated cytosines and ultraviolet light-induced mutagenesis.

Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana

The genomes of individuals from the same species vary in sequence as a result of different evolutionary processes. To examine the patterns of, and the forces shaping, sequence variation in Arabidopsis thaliana, we performed high-density array resequencing of 20 diverse strains (accessions). More than 1 million nonredundant single-nucleotide polymorphisms (SNPs) were identified at moderate false discovery rates (FDRs), and approximately 4% of the genome was identified as being highly dissimilar or deleted relative to the reference genome sequence. Patterns of polymorphism are highly nonrandom among gene families, with genes mediating interaction with the biotic environment having exceptional polymorphism levels. At the chromosomal scale, regional variation in polymorphism was readily apparent. A scan for recent selective sweeps revealed several candidate regions, including a notable example in which almost all variation was removed in a 500-kilobase window. Analyzing the polymorphisms we describe in larger sets of accessions will enable a detailed understanding of forces shaping population-wide sequence variation in A. thaliana.

Export of FT protein from phloem companion cells is sufficient for floral induction in Arabidopsis

Several endogenous and environmental factors need to be integrated to time the onset of flowering. Genetic and molecular analyses, primarily in Arabidopsis thaliana and rice, have shown that CONSTANS (CO) and FLOWERING LOCUS T (FT) play central roles in photoperiod-dependent flowering. The overall picture is that CO acts in the phloem companion cells of leaves and that its main effect is to induce FT mRNA in these cells. Surprisingly, FT, a small globular protein of 20 kDa, interacts at the shoot apex with the bZIP transcription factor FLOWERING LOCUS D (FD) to induce downstream targets. Given that green fluorescent protein (GFP), which as a monomer is 27 kDa, can be easily exported to sink tissue including flowers when expressed in phloem companion cells, the latter finding strongly implied that FT protein is the mobile floral-inductive signal. In agreement with this hypothesis, an FT-GFP fusion, just like GFP, can be exported from the phloem of both rice and Arabidopsis. It has been unknown, however, whether mobile FT protein is sufficient for transmitting the flowering signal. Here we show that FT mRNA is required in phloem companion cells where it acts partially redundant with its paralog TWIN SISTER OF FT (TSF) to induce flowering. Furthermore, we have devised a method that uncouples FT mRNA and protein effects in vivo. We demonstrate that export of FT protein from phloem companion cells is sufficient to induce flowering.

Autoimmune response as a mechanism for a Dobzhansky-Muller-type incompatibility syndrome in plants

Epistatic interactions between genes are a major factor in evolution. Hybrid necrosis is an example of a deleterious phenotype caused by epistatic interactions that is observed in many intra- and interspecific plant hybrids. A large number of hybrid necrosis cases share phenotypic similarities, suggesting a common underlying mechanism across a wide range of plant species. Here, we report that approximately 2% of intraspecific crosses in Arabidopsis thaliana yield F1 progeny that express necrosis when grown under conditions typical of their natural habitats. We show that several independent cases result from epistatic interactions that trigger autoimmune-like responses. In at least one case, an allele of an NB-LRR disease resistance gene homolog is both necessary and sufficient for the induction of hybrid necrosis, when combined with a specific allele at a second locus. The A. thaliana cases provide insights into the molecular causes of hybrid necrosis, and serve as a model for further investigation of intra- and interspecific incompatibilities caused by a simple epistatic interaction. Moreover, our finding that plant immune-system genes are involved in hybrid necrosis suggests that selective pressures related to host-pathogen conflict might cause the evolution of gene flow barriers in plants.

Sequencing of natural strains of Arabidopsis thaliana with short reads

Whole-genome hybridization studies have suggested that the nuclear genomes of accessions (natural strains) of Arabidopsis thaliana can differ by several percent of their sequence. To examine this variation, and as a first step in the 1001 Genomes Project for this species, we produced 15- to 25-fold coverage in Illumina sequencing-by-synthesis (SBS) reads for the reference accession, Col-0, and two divergent strains, Bur-0 and Tsu-1. We aligned reads to the reference genome sequence to assess data quality metrics and to detect polymorphisms. Alignments revealed 823,325 unique single nucleotide polymorphisms (SNPs) and 79,961 unique 1- to 3-bp indels in the divergent accessions at a specificity of >99%, and over 2000 potential errors in the reference genome sequence. We also identified >3.4 Mb of the Bur-0 and Tsu-1 genomes as being either extremely dissimilar, deleted, or duplicated relative to the reference genome. To obtain sequences for these regions, we incorporated the Velvet assembler into a targeted de novo assembly method. This approach yielded 10,921 high-confidence contigs that were anchored to flanking sequences and harbored indels as large as 641 bp. Our methods are broadly applicable for polymorphism discovery in moderate to large genomes even at highly diverged loci, and we established by subsampling the Illumina SBS coverage depth required to inform a broad range of functional and evolutionary studies. Our pipeline for aligning reads and predicting SNPs and indels, SHORE, is available for download at http://1001genomes.org.

Sequence and expression differences underlie functional specialization of Arabidopsis microRNAs miR159 and miR319

Many microRNAs (miRNAs) are encoded by small gene families. In a third of all conserved Arabidopsis miRNA families, members vary at two or more nucleotide positions. We have focused on the related miR159 and miR319 families, which share sequence identity at 17 of 21 nucleotides, yet affect different developmental processes through distinct targets. MiR159 regulates MYB mRNAs, while miR319 predominantly acts on TCP mRNAs. In the case of miR319, MYB targeting plays at most a minor role because miR319 expression levels and domain limit its ability to affect MYB mRNAs. In contrast, in the case of miR159, the miRNA sequence prevents effective TCP targeting. We complement these observations by identifying nucleotide positions relevant for miRNA activity with mutants recovered from a suppressor screen. Together, our findings reveal that functional specialization of miR159 and miR319 is achieved through both expression and sequence differences.

Highly specific gene silencing by artificial miRNAs in rice

Background

Endogenous microRNAs (miRNAs) are potent negative regulators of gene expression in plants and animals. Artificial miRNAs (amiRNAs)–designed to target one or several genes of interest–provide a new and highly specific approach for effective post-transcriptional gene silencing (PTGS) in plants.

Methodology

We devised an amiRNA-based strategy for both japonica and indica type strains of cultivated rice, Oryza sativa. Using an endogenous rice miRNA precursor and customized 21mers, we designed amiRNA constructs targeting three different genes (Pds, Spl11, and Eui1/CYP714D1). Upon constitutive expression of these amiRNAs in the varieties Nipponbare (japonica) and IR64 (indica), the targeted genes are down-regulated by amiRNA-guided cleavage of the transcripts, resulting in the expected mutant phenotypes. The effects are highly specific to the target gene, the transgenes are stably inherited and they remain effective in the progeny.

Conclusion/Significance

Our results not only show that amiRNAs can efficiently trigger gene silencing in a monocot crop, but also that amiRNAs can effectively modulate agronomically important traits in varieties used in modern breeding programs. We provide all software tools and a protocol for the design of rice amiRNA constructs, which can be easily adapted to other crops. The approach is suited for candidate gene validation, comparative functional genomics between different varieties, and for improvement of agronomic performance and nutritional value.

Quantitative trait locus mapping and DNA array hybridization identify an FLM deletion as a cause for natural flowering-time variation

Much of the flowering time variation in wild strains of Arabidopsis thaliana is due to allelic variation at two epistatically acting loci, FRIGIDA (FRI) and FLOWERING LOCUS C (FLC). FLC encodes a MADS (MCM1/AGAMOUS/DEFICIENS/SRF1) domain transcription factor that directly represses a series of flowering-promoting genes. FRI and FLC, however, do not explain all of the observed variation, especially when plants are grown in short days. To identify loci that act in addition to FRI and FLC in controlling flowering of natural accessions, we have analyzed a recombinant inbred line population derived from crosses of accession Niederzenz (Nd) to Columbia, both of which contain natural FRI lesions. Quantitative trait locus mapping and genomic DNA analysis by microarray hybridization were used to identify candidate genes affecting variation in flowering behavior. In both long and short days, the quantitative trait locus of largest effect, termed FLOWERING 1 (FLW1), was found to be associated with a Nd-specific deletion of FLOWERING LOCUS M (FLM), which encodes a floral repressor closely related to FLC. Analysis of near isogenic lines and quantitative transgenic complementation experiments confirmed that the FLM deletion is, in large part, responsible for the early flowering of the Nd accession.

The PHYTOCHROME C photoreceptor gene mediates natural variation in flowering and growth responses of Arabidopsis thaliana

Light has an important role in modulating seedling growth and flowering time. We show that allelic variation at the PHYTOCHROME C (PHYC) photoreceptor locus affects both traits in natural populations of A. thaliana. Two functionally distinct PHYC haplotype groups are distributed in a latitudinal cline dependent on FRIGIDA, a locus that together with FLOWERING LOCUS C explains a large portion of the variation in A. thaliana flowering time. In a genome-wide scan for association of 65 loci with latitude, there was an excess of significant P values, indicative of population structure. Nevertheless, PHYC was the most strongly associated locus across 163 strains, suggesting that PHYC alleles are under diversifying selection in A. thaliana. Our work, together with previous findings, suggests that photoreceptor genes are major agents of natural variation in plant flowering and growth response.

Local-scale patterns of genetic variability, outcrossing, and spatial structure in natural stands of Arabidopsis thaliana

As Arabidopsis thaliana is increasingly employed in evolutionary and ecological studies, it is essential to understand patterns of natural genetic variation and the forces that shape them. Previous work focusing mostly on global and regional scales has demonstrated the importance of historical events such as long-distance migration and colonization. Far less is known about the role of contemporary factors or environmental heterogeneity in generating diversity patterns at local scales. We sampled 1,005 individuals from 77 closely spaced stands in diverse settings around Tübingen, Germany. A set of 436 SNP markers was used to characterize genome-wide patterns of relatedness and recombination. Neighboring genotypes often shared mosaic blocks of alternating marker identity and divergence. We detected recent outcrossing as well as stretches of residual heterozygosity in largely homozygous recombinants. As has been observed for several other selfing species, there was considerable heterogeneity among sites in diversity and outcrossing, with rural stands exhibiting greater diversity and heterozygosity than urban stands. Fine-scale spatial structure was evident as well. Within stands, spatial structure correlated negatively with observed heterozygosity, suggesting that the high homozygosity of natural A. thaliana may be partially attributable to nearest-neighbor mating of related individuals. The large number of markers and extensive local sampling employed here afforded unusual power to characterize local genetic patterns. Contemporary processes such as ongoing outcrossing play an important role in determining distribution of genetic diversity at this scale. Local "outcrossing hotspots" appear to reshuffle genetic information at surprising rates, while other stands contribute comparatively little. Our findings have important implications for sampling and interpreting diversity among A. thaliana accessions.

Simultaneous alignment of short reads against multiple genomes

New software for the alignment of short-read sequence data to multiple genomes allows identification of polymorphisms that cannot be identified by alignment to a single reference genome.

Genome resequencing with short reads generally relies on alignments against a single reference. GenomeMapper supports simultaneous mapping of short reads against multiple genomes by integrating related genomes (e.g., individuals of the same species) into a single graph structure. It constitutes the first approach for handling multiple references and introduces representations for alignments against complex structures. Demonstrated benefits include access to polymorphisms that cannot be identified by alignments against the reference alone. Download GenomeMapper at .

Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation

Despite evolutionary conserved mechanisms to silence transposable element activity, there are drastic differences in the abundance of transposable elements even among closely related plant species. We conducted a de novo assembly for the 375 Mb genome of the perennial model plant, Arabis alpina. Analysing this genome revealed long-lasting and recent transposable element activity predominately driven by Gypsy long terminal repeat retrotransposons, which extended the low-recombining pericentromeres and transformed large formerly euchromatic regions into repeat-rich pericentromeric regions. This reduced capacity for long terminal repeat retrotransposon silencing and removal in A. alpina co-occurs with unexpectedly low levels of DNA methylation. Most remarkably, the striking reduction of symmetrical CG and CHG methylation suggests weakened DNA methylation maintenance in A. alpina compared with Arabidopsis thaliana. Phylogenetic analyses indicate a highly dynamic evolution of some components of methylation maintenance machinery that might be related to the unique methylation in A. alpina.

The recombination landscape in Arabidopsis thaliana F2 populations

Recombination during meiosis shapes the complement of alleles segregating in the progeny of hybrids, and has important consequences for phenotypic variation. We examined allele frequencies, as well as crossover (XO) locations and frequencies in over 7000 plants from 17 F(2) populations derived from crosses between 18 Arabidopsis thaliana accessions. We observed segregation distortion between parental alleles in over half of our populations. The potential causes of distortion include variation in seed dormancy and lethal epistatic interactions. Such a high occurrence of distortion was only detected here because of the large sample size of each population and the number of populations characterized. Most plants carry only one or two XOs per chromosome pair, and therefore inherit very large, non-recombined genomic fragments from each parent. Recombination frequencies vary between populations but consistently increase adjacent to the centromeres. Importantly, recombination rates do not correlate with whole-genome sequence differences between parental accessions, suggesting that sequence diversity within A. thaliana does not normally reach levels that are high enough to exert a major influence on the formation of XOs. A global knowledge of the patterns of recombination in F(2) populations is crucial to better understand the segregation of phenotypic traits in hybrids, in the laboratory or in the wild.

A promiscuous intermediate underlies the evolution of LEAFY DNA binding specificity

Transcription factors (TFs) are key players in evolution. Changes affecting their function can yield novel life forms but may also have deleterious effects. Consequently, gene duplication events that release one gene copy from selective pressure are thought to be the common mechanism by which TFs acquire new activities. Here, we show that LEAFY, a major regulator of flower development and cell division in land plants, underwent changes to its DNA binding specificity, even though plant genomes generally contain a single copy of the LEAFY gene. We examined how these changes occurred at the structural level and identify an intermediate LEAFY form in hornworts that appears to adopt all different specificities. This promiscuous intermediate could have smoothed the evolutionary transitions, thereby allowing LEAFY to evolve new binding specificities while remaining a single-copy gene.

QTL mapping in new Arabidopsis thaliana advanced intercross-recombinant inbred lines

Background

Even when phenotypic differences are large between natural or domesticated strains, the underlying genetic basis is often complex, and causal genomic regions need to be identified by quantitative trait locus (QTL) mapping. Unfortunately, QTL positions typically have large confidence intervals, which can, for example, lead to one QTL being masked by another, when two closely linked loci are detected as a single QTL. One strategy to increase the power of precisely localizing small effect QTL, is the use of an intercross approach before inbreeding to produce Advanced Intercross RILs (AI-RILs).

Methodology/Principal Findings

We present two new AI-RIL populations of Arabidopsis thaliana genotyped with an average intermarker distance of 600 kb. The advanced intercrossing design led to expansion of the genetic map in the two populations, which contain recombination events corresponding to 50 kb/cM in an F₂ population. We used the AI-RILs to map QTL for light response and flowering time, and to identify segregation distortion in one of the AI-RIL populations due to a negative epistatic interaction between two genomic regions.

Conclusions/Significance

The two new AI-RIL populations, EstC and KendC, derived from crosses of Columbia (Col) to Estland (Est-1) and Kendallville (Kend-L) provide an excellent resource for high precision QTL mapping. Moreover, because they have been genotyped with over 100 common markers, they are also excellent material for comparative QTL mapping.

Major-effect alleles at relatively few loci underlie distinct vernalization and flowering variation in Arabidopsis accessions

We have explored the genetic basis of variation in vernalization requirement and response in Arabidopsis accessions, selected on the basis of their phenotypic distinctiveness. Phenotyping of F2 populations in different environments, plus fine mapping, indicated possible causative genes. Our data support the identification of FRI and FLC as candidates for the major-effect QTL underlying variation in vernalization response, and identify a weak FLC allele, caused by a Mutator-like transposon, contributing to flowering time variation in two N. American accessions. They also reveal a number of additional QTL that contribute to flowering time variation after saturating vernalization. One of these was the result of expression variation at the FT locus. Overall, our data suggest that distinct phenotypic variation in the vernalization and flowering response of Arabidopsis accessions is accounted for by variation that has arisen independently at relatively few major-effect loci.

Directed gene silencing with artificial microRNAs

The characterization of gene function typically includes a detailed analysis of loss-of-function alleles. In model plants, such as Arabidopsis thaliana and rice, sequence-indexed insertion collections provide a large resource of potential null alleles that can often be easily accessed through convenient Web sites (e.g., http://signal.salk.edu ). They are, however, not available for nonmodel species, require stacking for knockout of redundant homologs, and do not easily allow for partial or regulated loss of gene function, which is particularly useful when null alleles are lethal. Transgene approaches that employ directed gene silencing can substitute for null alleles and also enable refined studies of gene function, e.g., by tissue-specific and inducible gene-silencing. This chapter describes the generation and application of artificial microRNAs (amiRNAs) as a gene silencing tool in a wide variety of different plant species.

Comparative analysis of the MIR319a microRNA locus in Arabidopsis and related Brassicaceae

MicroRNAs (miRNAs) are important regulators of gene expression in multicellular organisms. Yet, little is known about their molecular evolution. The 20- to 22-nt long miRNAs are processed in plants from foldbacks that are a few hundred base pairs in size. Often, these foldbacks are embedded in much larger precursor transcripts. To investigate functional constraints on sequence evolution of miRNA precursor genes, we have studied sequence variation in the precursor of miR319a, MIR319a, between species from the Brassicaceae. We compared the genomic context in Arabidopsis thaliana, Arabidopsis halleri, and Capsella rubella, using bacterial artificial chromosome clones, and analyzed precursor sequences obtained by polymerase chain reaction from 13 additional species. Phylogenetic shadowing identifies a conserved motif around the transcription start site, which we demonstrate to be functionally important. We further assessed the functionality of MIR319a orthologs from several Brassicaceae species in A. thaliana. The ortholog from kale (Brassica oleracea var. acephala) was found to be largely inactive, at least partially due to mutations in the miRNA itself, but experimental evidence suggests that loss of miR319a function is compensated by other members of the miR319 family. More broadly, we find that the foldback diverges less rapidly than the remainder of the primary transcript. To understand the molecular evolution of miRNA genes, investigations at different levels of phylogenetic divergence are required.

MSQT for choosing SNP assays from multiple DNA alignments

MOTIVATION

One challenging aspect of genotyping and association mapping projects is often the identification of markers that are informative between groups of individuals and to convert these into genotyping assays.

RESULTS

The Multiple SNP Query Tool (MSQT) extracts SNP information from multiple sequence alignments, stores it in a database, provides a web interface to query the database and outputs SNP information in a format directly applicable for SNP-assay design. MSQT was applied to Arabidopsis thaliana sequence data to develop SNP genotyping assays that distinguish a recurrent parent (Col-0) from five other strains. SNPs with intermediate allele frequencies were also identified and developed into markers suitable for efficient genetic mapping among random pairs of wild strains.

AVAILABILITY

The source code for MSQT is available at http://msqt.weigelworld.org, together with an online instance of MSQT containing data on 1214 sequenced fragments from 96 ecotypes (wild inbred strains) of the reference plant A. thaliana. All SNP genotyping assays are available in several formats for broad community use.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Do longer root hairs improve phosphorus uptake? Testing the hypothesis with transgenic Brachypodium distachyon lines overexpressing endogenous RSL genes

Mutants without root hairs show reduced inorganic orthophosphate (Pi) uptake and compromised growth on soils when Pi availability is restricted. What is less clear is whether root hairs that are longer than wild-type provide an additional benefit to phosphorus (P) nutrition. This was tested using transgenic Brachypodium lines with longer root hairs. The lines were transformed with the endogenous BdRSL2 and BdRSL3 genes using either a constitutive promoter or a root hair-specific promoter. Plants were grown for 32 d in soil amended with various Pi concentrations. Plant biomass and P uptake were measured and genotypes were compared on the basis of critical Pi values and P uptake per unit root length. Ectopic expression of RSL2 and RSL3 increased root hair length three-fold but decreased plant biomass. Constitutive expression of BdRSL2, but not expression of BdRSL3, consistently improved P nutrition as measured by lowering the critical Pi values and increasing Pi uptake per unit root length. Increasing root hair length through breeding or biotechnology can improve P uptake efficiency if the pleotropic effects on plant biomass are avoided. Long root hairs, alone, appear to be insufficient to improve Pi uptake and need to be combined with other traits to benefit P nutrition.

kWIP: The k-mer weighted inner product, a de novo estimator of genetic similarity

Modern genomics techniques generate overwhelming quantities of data. Extracting population genetic variation demands computationally efficient methods to determine genetic relatedness between individuals (or "samples") in an unbiased manner, preferably de novo. Rapid estimation of genetic relatedness directly from sequencing data has the potential to overcome reference genome bias, and to verify that individuals belong to the correct genetic lineage before conclusions are drawn using mislabelled, or misidentified samples. We present the k-mer Weighted Inner Product (kWIP), an assembly-, and alignment-free estimator of genetic similarity. kWIP combines a probabilistic data structure with a novel metric, the weighted inner product (WIP), to efficiently calculate pairwise similarity between sequencing runs from their k-mer counts. It produces a distance matrix, which can then be further analysed and visualised. Our method does not require prior knowledge of the underlying genomes and applications include establishing sample identity and detecting mix-up, non-obvious genomic variation, and population structure. We show that kWIP can reconstruct the true relatedness between samples from simulated populations. By re-analysing several published datasets we show that our results are consistent with marker-based analyses. kWIP is written in C++, licensed under the GNU GPL, and is available from https://github.com/kdmurray91/kwip.

ESTs and EST-linked polymorphisms for genetic mapping and phylogenetic reconstruction in the guppy, Poecilia reticulata

BACKGROUND

The guppy, Poecilia reticulata, is a well-known model organism for studying inheritance and variation of male ornamental traits as well as adaptation to different river habitats. However, genomic resources for studying this important model were not previously widely available.

RESULTS

With the aim of generating molecular markers for genetic mapping of the guppy, cDNA libraries were constructed from embryos and different adult organs to generate expressed sequence tags (ESTs). About 18,000 ESTs were annotated according to BLASTN and BLASTX results and the sequence information from the 3' UTRs was exploited to generate PCR primers for re-sequencing of genomic DNA from different wild type strains. By comparison of EST-linked genomic sequences from at least four different ecotypes, about 1,700 polymorphisms were identified, representing about 400 distinct genes. Two interconnected MySQL databases were built to organize the ESTs and markers, respectively. A robust phylogeny of the guppy was reconstructed, based on 10 different nuclear genes.

CONCLUSION

Our EST and marker databases provide useful tools for genetic mapping and phylogenetic studies of the guppy.