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.

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.

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.

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.

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.

Artificial microRNAs for specific gene silencing in rice.

Artificial microRNAs (amiRNAs) have been shown to facilitate efficient gene silencing with high specificity to the intended target gene(s). For the plant breeder, gene silencing by artificial miRNAs will certainly accelerate gene discovery, because it allows targeting of all genes in a mapping interval, independent of the genetic background. In addition, beneficial knockout phenotypes can easily be transferred between varieties and across incompatibility barriers. This chapter describes the generation and application of amiRNAs as a gene silencing tool in rice.

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.

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.

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.

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.

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.

Simultaneous alignment of short reads against multiple genomes.

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 http://1001genomes.org.

Cis-regulatory changes at FLOWERING LOCUS T mediate natural variation in flowering responses of Arabidopsis thaliana.

Flowering time, a critical adaptive trait, is modulated by several environmental cues. These external signals converge on a small set of genes that in turn mediate the flowering response. Mutant analysis and subsequent molecular studies have revealed that one of these integrator genes, FLOWERING LOCUS T (FT), responds to photoperiod and temperature cues, two environmental parameters that greatly influence flowering time. As the central player in the transition to flowering, the protein coding sequence of FT and its function are highly conserved across species. Using QTL mapping with a new advanced intercross-recombinant inbred line (AI-RIL) population, we show that a QTL tightly linked to FT contributes to natural variation in the flowering response to the combined effects of photoperiod and ambient temperature. Using heterogeneous inbred families (HIF) and introgression lines, we fine map the QTL to a 6.7 kb fragment in the FT promoter. We confirm by quantitative complementation that FT has differential activity in the two parental strains. Further support for FT underlying the QTL comes from a new approach, quantitative knockdown with artificial microRNAs (amiRNAs). Consistent with the causal sequence polymorphism being in the promoter, we find that the QTL affects FT expression. Taken together, these results indicate that allelic variation at pathway integrator genes such as FT can underlie phenotypic variability and that this may be achieved through cis-regulatory changes.

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(2) 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.