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.

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.

Reference-guided assembly of four diverse Arabidopsis thaliana genomes.

We present whole-genome assemblies of four divergent Arabidopsis thaliana strains that complement the 125-Mb reference genome sequence released a decade ago. Using a newly developed reference-guided approach, we assembled large contigs from 9 to 42 Gb of Illumina short-read data from the Landsberg erecta (Ler-1), C24, Bur-0, and Kro-0 strains, which have been sequenced as part of the 1,001 Genomes Project for this species. Using alignments against the reference sequence, we first reduced the complexity of the de novo assembly and later integrated reads without similarity to the reference sequence. As an example, half of the noncentromeric C24 genome was covered by scaffolds that are longer than 260 kb, with a maximum of 2.2 Mb. Moreover, over 96% of the reference genome was covered by the reference-guided assembly, compared with only 87% with a complete de novo assembly. Comparisons with 2 Mb of dideoxy sequence reveal that the per-base error rate of the reference-guided assemblies was below 1 in 10,000. Our assemblies provide a detailed, genomewide picture of large-scale differences between A. thaliana individuals, most of which are difficult to access with alignment-consensus methods only. We demonstrate their practical relevance in studying the expression differences of polymorphic genes and show how the analysis of sRNA sequencing data can lead to erroneous conclusions if aligned against the reference genome alone. Genome assemblies, raw reads, and further information are accessible through http://1001genomes.org/projects/assemblies.html.

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.

The scale of population structure in Arabidopsis thaliana.

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales.

Cost-conscious generation of multiplexed short-read DNA libraries for whole-genome sequencing.

Massively parallel, second-generation short-read DNA sequencing has become an integral tool in biology for genomic studies. Offering highly accurate base-pair resolution at the most competitive price, the technology has become widespread. However, high-throughput generation of multiplexed DNA libraries can be costly and cumbersome. Here, we present a cost-conscious protocol for generating multiplexed short-read DNA libraries using a bead-linked transposome from Illumina. We prepare libraries in high-throughput with small reaction volumes that use 1/50th the amount of transposome compared to Illumina DNA Prep tagmentation protocols. By reducing transposome usage and optimising the protocol to circumvent magnetic bead-based clean-ups between steps, we reduce costs, labour time and DNA input requirements. Developing our own dual index primers further reduced costs and enables up to nine 96-well microplate combinations. This facilitates efficient usage of large-scale sequencing platforms, such as the Illumina NovaSeq 6000, which offers up to three terabases of sequencing per S4 flow cell. The protocol presented substantially reduces the cost per library by approximately 1/20th compared to conventional Illumina methods.

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.

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.

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.

Global Diversity of the Brachypodium Species Complex as a Resource for Genome-Wide Association Studies Demonstrated for Agronomic Traits in Response to Climate.

The development of model systems requires a detailed assessment of standing genetic variation across natural populations. The Brachypodium species complex has been promoted as a plant model for grass genomics with translation to small grain and biomass crops. To capture the genetic diversity within this species complex, thousands of Brachypodium accessions from around the globe were collected and genotyped by sequencing. Overall, 1897 samples were classified into two diploid or allopolyploid species, and then further grouped into distinct inbred genotypes. A core set of diverse B. distachyon diploid lines was selected for whole genome sequencing and high resolution phenotyping. Genome-wide association studies across simulated seasonal environments was used to identify candidate genes and pathways tied to key life history and agronomic traits under current and future climatic conditions. A total of 8, 22, and 47 QTL were identified for flowering time, early vigor, and energy traits, respectively. The results highlight the genomic structure of the Brachypodium species complex, and the diploid lines provided a resource that allows complex trait dissection within this grass model species.

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.

Analysis of DNAs associated with coconut foliar decay disease implicates a unique single-stranded DNA virus representing a new taxon.

The unique ecology, pathology and undefined taxonomy of coconut foliar decay virus (CFDV), found associated with coconut foliar decay disease (CFD) in 1986, prompted analyses of old virus samples by modern methods. Rolling circle amplification and deep sequencing applied to nucleic acid extracts from virion preparations and CFD-affected palms identified twelve distinct circular DNAs, eleven of which had a size of about 1.3 kb and one of 641 nt. Mass spectrometry-based protein identification proved that a 24 kDa protein encoded by two 1.3 kb DNAs is the virus capsid protein with highest sequence similarity to that of grabloviruses (family Geminiviridae), even though CFDV particles are not geminate. The nine other 1.3 kb DNAs represent alphasatellites coding for replication initiator proteins that differ clearly from those encoded by nanovirid DNA-R. The 641 nt DNA-gamma is unique and may encode a movement protein. Three DNAs, alphasatellite CFDAR, capsid protein encoding CFDV DNA-S.1 and DNA-gamma share sequence motifs near their replication origins and were consistently present in all samples analysed. These DNAs appear to be integral components of a possibly tripartite CFDV genome, different from those of any Geminiviridae or Nanoviridae family member, implicating CFDV as representative of a new genus and family.

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.

DNA metabarcoding of unfractionated water samples relates phyto-, zoo- and bacterioplankton dynamics and reveals a single-taxon bacterial bloom.

Most studies of aquatic plankton focus on either macroscopic or microbial communities, and on either eukaryotes or prokaryotes. This separation is primarily for methodological reasons, but can overlook potential interactions among groups. Here we tested whether DNA metabarcoding of unfractionated water samples with universal primers could be used to qualitatively and quantitatively study the temporal dynamics of the total plankton community in a shallow temperate lake. Significant changes in the relative proportions of normalized sequence reads of eukaryotic and prokaryotic plankton communities over a 3-month period in spring were found. Patterns followed the same trend as plankton estimates measured using traditional microscopic methods. The bloom of a conditionally rare bacterial taxon belonging to Arcicella was characterized, which rapidly came to dominate the whole lake ecosystem and would have remained unnoticed without metabarcoding. The data demonstrate the potential of universal DNA metabarcoding applied to unfractionated samples for providing a more holistic view of plankton communities.

Evolution. Response to Comment on "A promiscuous intermediate underlies the evolution of LEAFY DNA binding specificity".

Brunkard et al. propose that the identification of novel LEAFY sequences contradicts our model of evolution through promiscuous intermediates. Based on the debate surrounding land plant phylogeny and on our analysis of these interesting novel sequences, we explain why there is no solid evidence to disprove our model.

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.

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.