Kevin D. Murray

The khmer software package: enabling efficient nucleotide sequence analysis

The khmer package is a freely available software library for working efficiently with fixed length DNA words, or k-mers. khmer provides implementations of a probabilistic k-mer counting data structure, a compressible De Bruijn graph representation, De Bruijn graph partitioning, and digital normalization. khmer is implemented in C++ and Python, and is freely available under the BSD license at https://github.com/dib-lab/khmer/.

TraitCapture: genomic and environment modelling of plant phenomic data

Agriculture requires a second green revolution to provide increased food, fodder, fiber, fuel and soil fertility for a growing population while being more resilient to extreme weather on finite land, water, and nutrient resources. Advances in phenomics, genomics and environmental control/sensing can now be used to directly select yield and resilience traits from large collections of germplasm if software can integrate among the technologies. Traits could be Captured throughout development and across environments from multi-dimensional phenotypes, by applying Genome Wide Association Studies (GWAS) to identify causal genes and background variation and functional structural plant models (FSPMs) to predict plant growth and reproduction in target environments. TraitCapture should be applicable to both controlled and field environments and would allow breeders to simulate regional variety trials to pre-select for increased productivity under challenging environments.

A chromatin modifying enzyme, SDG8, is involved in morphological, gene expression, and epigenetic responses to mechanical stimulation.

Thigmomorphogenesis is viewed as being a response process of acclimation to short repetitive bursts of mechanical stimulation or touch. The underlying molecular mechanisms that coordinate changes in how touch signals lead to long-term morphological changes are enigmatic. Touch responsive gene expression is rapid and transient, and no transcription factor or DNA regulatory motif has been reported that could confer a genome wide mechanical stimulus. We report here on a chromatin modifying enzyme, SDG8/ASHH2, which can regulate the expression of many touch responsive genes identified in Arabidopsis. SDG8 is required for the permissive expression of touch induced genes; and the loss of function of sdg8 perturbs the maximum levels of induction on selected touch gene targets. SDG8 is required to maintain permissive H3K4 trimethylation marks surrounding the Arabidopsis touch-inducible gene TOUCH 3 (TCH3), which encodes a calmodulin-like protein (CML12). The gene neighboring was also slightly down regulated, revealing a new target for SDG8 mediated chromatin modification. Finally, sdg8 mutants show perturbed morphological response to wind-agitated mechanical stimuli, implicating an epigenetic memory-forming process in the acclimation response of thigmomorphogenesis.

Identification and qualification of 500 nuclear, single-copy, orthologous genes for the Eupulmonata (Gastropoda) using transcriptome sequencing and exon capture.

The qualification of orthology is a significant challenge when developing large, multiloci phylogenetic data sets from assembled transcripts. Transcriptome assemblies have various attributes, such as fragmentation, frameshifts and mis‐indexing, which pose problems to automated methods of orthology assessment. Here, we identify a set of orthologous single‐copy genes from transcriptome assemblies for the land snails and slugs (Eupulmonata) using a thorough approach to orthology determination involving manual alignment curation, gene tree assessment and sequencing from genomic DNA. We qualified the orthology of 500 nuclear, protein‐coding genes from the transcriptome assemblies of 21 eupulmonate species to produce the most complete phylogenetic data matrix for a major molluscan lineage to date, both in terms of taxon and character completeness. Exon capture targeting 490 of the 500 genes (those with at least one exon >120 bp) from 22 species of Australian Camaenidae successfully captured sequences of 2825 exons (representing all targeted genes), with only a 3.7% reduction in the data matrix due to the presence of putative paralogs or pseudogenes. The automated pipeline Agalma retrieved the majority of the manually qualified 500 single‐copy gene set and identified a further 375 putative single‐copy genes, although it failed to account for fragmented transcripts resulting in lower data matrix completeness when considering the original 500 genes. This could potentially explain the minor inconsistencies we observed in the supported topologies for the 21 eupulmonate species between the manually curated and ‘Agalma‐equivalent’ data set (sharing 458 genes). Overall, our study confirms the utility of the 500 gene set to resolve phylogenetic relationships at a range of evolutionary depths and highlights the importance of addressing fragmentation at the homolog alignment stage for probe design.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis

Abiotic stress and recovery transcriptomes reveal extremely short RNA half-lives, changes to cotranslational decay, and recovery-specific networks, consistent with active recovery and cellular memory. Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5ʹ-3ʹ RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

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.

Axe: rapid, competitive sequence read demultiplexing using a trie

Summary We describe a rapid algorithm for demultiplexing DNA sequence reads with in‐read indices. Axe selects the optimal index present in a sequence read, even in the presence of sequencing errors. The algorithm is able to handle combinatorial indexing, indices of differing length and several mismatches per index sequence. Availability and implementation Axe is implemented in C, and is used as a command‐line program on Unix‐like systems. Axe is available online at https://github.com/kdmurray91/axe, and is available in Debian/Ubuntu distributions of GNU/Linux as the package axe‐demultiplexer. Supplementary information Supplementary data are available at Bioinformatics online

Population structure of the Brachypodium species complex and genome wide association of 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 translational 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 sequenced using genotyping by sequencing (GBS). Overall, 1,897 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 were 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 QTLs were identified for flowering time, early vigour and energy traits, respectively. Overall, the results highlight the genomic structure of the Brachypodium species complex and allow powerful complex trait dissection within this new grass model species.

RNA Polymerase II read-through promotes expression of neighboring genes in SAL1-PAP-XRN retrograde signaling.

Consequences of transcription out of bounds: a retrograde signal can trigger RNA Polymerase II read-through, upregulating the expression of downstream genes. In plants, the molecular function(s) of the nucleus-localized 5′-3′ EXORIBONUCLEASES (XRNs) are unclear; however, their activity is reported to have a significant effect on gene expression and SAL1-mediated retrograde signaling. Using parallel analysis of RNA ends, we documented a dramatic increase in uncapped RNA substrates of the XRNs in both sal1 and xrn2xrn3 mutants. We found that a major consequence of reducing SAL1 or XRN activity was RNA Polymerase II 3′ read-through. This occurred at 72% of expressed genes, demonstrating a major genome-wide role for the XRN-torpedo model of transcription termination in Arabidopsis (Arabidopsis thaliana). Read-through is speculated to have a negative effect on transcript abundance; however, we did not observe this. Rather, we identified a strong association between read-through and increased transcript abundance of tandemly orientated downstream genes, strongly correlated with the proximity (less than 1,000 bp) and expression of the upstream gene. We observed read-through in the proximity of 903 genes up-regulated in the sal1-8 retrograde signaling mutant; thus, this phenomenon may account directly for up to 23% of genes up-regulated in sal1-8. Using APX2 and AT5G43770 as exemplars, we genetically uncoupled read-through loci from downstream genes to validate the principle of read-through-mediated mRNA regulation, providing one mechanism by which an ostensibly posttranscriptional exoribonuclease that targets uncapped RNAs could modulate gene expression.