Search

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. [ABSTRACT FROM AUTHOR]

Rapid advancements in long-read sequencing technologies have transformed read lengths from bps to Mbps, which has enabled chromosome-scale genome assemblies. However, read lengths are now becoming limited by the extraction of pure high-molecular weight DNA suitable for long-read sequencing, which is particularly challenging in plants and fungi. To overcome this, we present a protocol collection; high-molecular weight DNA extraction, clean-up and size selection for long-read sequencing. We optimised a gentle magnetic bead based high-molecular weight DNA extraction, which is presented here in detail. The protocol circumvents spin columns and high-centrifugation, to limit DNA fragmentation. The protocol is scalable based on tissue input, which can be used on many species of plants, fungi, reptiles and bacteria. It is also cost effective compared to kit-based protocols and hence applicable at scale in low resource settings. An optional sorbitol wash is listed and is highly recommended for plant and fungal tissues. To further remove any remaining contaminants such as phenols and polysaccharides, optional DNA clean-up and size selection strategies are given. This protocol collection is suitable for all common long-read sequencing platforms, such as technologies offered by PacBio and Oxford Nanopore. Using these protocols, sequencing on the Oxford Nanopore MinION can achieve read length N50 values of 30–50 kb, with reads exceeding 200 kb and outputs ranging from 15–30 Gbp. This has been routinely achieved with various plant, fungi, animal and bacteria samples. [ABSTRACT FROM AUTHOR]

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 (), an assembly-, and alignment-free estimator of genetic similarity. 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 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. is written in C++, licensed under the GNU GPL, and is available from . [ABSTRACT FROM AUTHOR]

Glyoxalase 1 (Glo1) has been implicated in anxiety-like behavior in mice and in multiple psychiatric diseases in humans. We used mouse Affymetrix exon arrays to detect copy number variants (CNV) among inbred mouse strains and thereby identified a ~475 kb tandem duplication on chromosome 17 that includes Glo1 (30,174,390-30,651,226 Mb; mouse genome build 36). We developed a PCR-based strategy and used it to detect this duplication in 23 of 71 inbred strains tested, and in various outbred and wild-caught mice. Presence of the duplication is associated with a cis-acting expression QTL for Glo1 (LOD>30) in BXD recombinant inbred strains. However, evidence for an eQTL for Glo1 was not obtained when we analyzed single SNPs or 3-SNP haplotypes in a panel of 27 inbred strains. We conclude that association analysis in the inbred strain panel failed to detect an eQTL because the duplication was present on multiple highly divergent haplotypes. Furthermore, we suggest that non-allelic homologous recombination has led to multiple reversions to the non-duplicated state among inbred strains. We show associations between multiple duplication-containing haplotypes, Glo1 expression and anxiety-like behavior in both inbred strain panels and outbred CD-1 mice. Our findings provide a molecular basis for differential expression of Glo1 and further implicate Glo1 in anxiety-like behavior. More broadly, these results identify problems with commonly employed tests for association in inbred strains when CNVs are present. Finally, these data provide an example of biologically significant phenotypic variability in model organisms that can be attributed to CNVs. [ABSTRACT FROM AUTHOR]

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 F2 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. [ABSTRACT FROM AUTHOR]

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

The female gametophyte of flowering plants, the embryo sac, develops within the diploid (sporophytic) tissue of the ovule. While embryo sac-expressed genes are known to be required at multiple stages of the fertilization process, the set of embryo sac-expressed genes has remained poorly defined. In particular, the set of genes responsible for mediating intracellular communication between the embryo sac and the male gametophyte, the pollen grain, is unknown. We used high-throughput cDNA sequencing and whole-genome tiling arrays to compare gene expression in wild-type ovules to that in dif1 ovules, which entirely lack embryo sacs, and myb98 ovules, which are impaired in pollen tube attraction. We identified nearly 400 genes that are downregulated in dif1 ovules. Seventy-eight percent of these embryo sac-dependent genes were predicted to encode for secreted proteins, and 60% belonged to multigenic families. Our results define a large number of candidate extracellular signaling molecules that may act during embryo sac development or fertilization; less than half of these are represented on the widely used ATH1 expression array. In particular, we found that 37 out of 40 genes encoding Domain of Unknown Function 784 (DUF784) domains require the synergid-specific transcription factor MYB98 for expression. Several DUF784 genes were transcribed in synergid cells of the embryo sac, implicating the DUF784 gene family in mediating late stages of embryo sac development or interactions with pollen tubes. The coexpression of highly similar proteins suggests a high degree of functional redundancy among embryo sac genes. [ABSTRACT FROM AUTHOR]

Discovering novel genes involved in immune evasion and drug resistance in the human malaria parasite, Plasmodium falciparum, is of critical importance to global health. Such knowledge may assist in the development of new effective vaccines and in the appropriate use of antimalarial drugs. By performing a full-genome scan of allelic variability in 14 field and laboratory strains of P. falciparum, we comprehensively identified ≈500 genes evolving at higher than neutral rates. The majority of the most variable genes have paralogs within the P. falciparum genome and may be subject to a different evolutionary clock than those without. The group of 211 variable genes without paralogs contains most known immunogens and a few drug targets, consistent with the idea that the human immune system and drug use is driving parasite evolution. We also reveal gene-amplification events including one surrounding pfmdr1, the P. falciparum multidrug-resistance gene, and a previously uncharacterized amplification centered around the P. falciparum GTP cyclohydrolase gene, the first enzyme in the folate biosynthesis pathway. Although GTP cyclohydrolase is not the known target of any current drugs, downstream members of the pathway are targeted by several widely used antimalarials. We speculate that an amplification of the GTP cyclohydrolase enzyme in the folate biosynthesis pathway may increase flux through this pathway and facilitate parasite resistance to antifolate drugs. [ABSTRACT FROM AUTHOR]

Corrections to the article "A Systematic Map of Genetic Variation in Plasmodium Falciparum," by Claire Kidgell et al which was published in volume 2, issue 6 of this periodical are presented.

The genetic model plant Arabidopsis thaliana, like many plant species, experiences a range of edaphic conditions across its natural habitat. Such heterogeneity may drive local adaptation, though the molecular genetic basis remains elusive. Here, we describe a study in which we used genome-wide association mapping, genetic complementation, and gene expression studies to identify cis-regulatory expression level polymorphisms at the AtHKT1;1 locus, encoding a known sodium (Na(+)) transporter, as being a major factor controlling natural variation in leaf Na(+) accumulation capacity across the global A. thaliana population. A weak allele of AtHKT1;1 that drives elevated leaf Na(+) in this population has been previously linked to elevated salinity tolerance. Inspection of the geographical distribution of this allele revealed its significant enrichment in populations associated with the coast and saline soils in Europe. The fixation of this weak AtHKT1;1 allele in these populations is genetic evidence supporting local adaptation to these potentially saline impacted environments., Competing Interests: The authors have declared that no competing interests exist.

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., Competing Interests: The authors have declared that no competing interests exist.

Despite their unprecedented density, current SNP genotyping arrays contain large amounts of redundancy, with up to 40 oligonucleotide features used to query each SNP. By using publicly available reference genotype data from the International HapMap, we show that 93.6% sensitivity at <5% false positive rate can be obtained with only four probes per SNP, compared with 98.3% with the full data set. Removal of this redundancy will allow for more comprehensive whole-genome association studies with increased SNP density and larger sample sizes.

Flowering time in Arabidopsis thaliana is controlled by a large number of genes and various environmental factors, such as light and temperature. The objective of this study was to identify flowering time quantitative trait loci (QTL) under growth conditions simulating seasonal conditions from native geographic locations. Our growth chambers were set to simulate the spring conditions in Spain and Sweden, with appropriate changes in light color, intensity and day length, as well as temperature and relative humidity. Thus the Sweden-like spring conditions changed more dramatically compared to Spain-like spring conditions across the duration of our experiment. We have used these conditions to map QTL responsible for flowering time in the Kas-1/Col-gl1 recombinant inbred lines (RILs) across two replicate blocks. A linkage map from 96 RILs was established using 119 markers including 64 new SNPs markers. One major QTL, mapping to the FRIGIDA (FRI) locus, was detected on the top of chromosome 4 that showed significant gene x seasonal environment interactions. Three other minor QTL also were detected. One QTL mapping near FLOWERING LOCUS M (FLM) showed an epistatic interaction with the QTL at FRI. These QTLxenvironment and QTL x QTL interactions suggest that subtle ecologically relevant changes in light, temperature, and relative humidity are differentially felt by alleles controlling flowering time and may be responsible for adaptation to regional environments.