Landscape drivers of genomic diversity and divergence in woodland Eucalyptus.

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the "eucalypts") are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole-genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome-wide, interspecific differentiation (FST  = 0.15) and intraspecific differentiation between localities (FST  ≈ 0.01-0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.

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.

RNA Polymerase II Read-Through Promotes Expression of Neighboring Genes in SAL1-PAP-XRN Retrograde Signaling.

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.

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.

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.

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

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.

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.

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/.

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.

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.

Dynamic interface pressure distributions of two transtibial prosthetic socket concepts.

In this study, we investigated and compared the dynamic interface pressure distribution of hands-off and hands-on transtibial prosthetic systems by means of pressure mapping. Of the 48 established unilateral amputees recruited, half (n = 24) had been wearing pressure-cast prostheses (IceCast Compact) and the other half (n = 24) had been wearing hand-cast sockets of the patellar tendon bearing design. We measured the dynamic pressure profile of more than 90% of the area within each prosthetic socket by means of four Tekscan F-Scan socket transducer arrays. We compared the interface pressure between socket concepts. We found that the distribution of dynamic pressure at the limb-socket interface was similar for the two intervention (socket prescription) groups. However, a significant difference was found in the magnitude of the interface pressure between the two socket concepts; the interface pressures recorded in the hands-off sockets were higher than those seen in the hands-on concept. Despite the differences in interface pressure, the level of satisfaction with the sockets was similar between subject groups. The sockets instrumented for this study had been in daily use for at least 6 months, with no residual-limb health problems.

The influence of pulmonary staple line reinforcement on air leaks.

OBJECTIVES: Although uncommon, prolonged postoperative air leaks are a troublesome complication of lung surgery. This study was performed to determine if buttressing pulmonary staple lines would reduce air leakage at varying airway pressures, and if there was a difference between buttressing materials. METHODS: Using cadaver lungs, the development of air leak from staple lines was evaluated at incremental airway pressures. Unreinforced staples were compared to staples reinforced with bovine pericardium and staples reinforced with expanded polytetrafluorethylene (ePTFE). RESULTS: Unreinforced staple lines began to leak air at an airway pressure of 20 mm Hg, and > 90% leaked at a pressure of 35 mm Hg. Both bovine pericardium and ePTFE significantly reduced the incidence of air leak at these airway pressures. At higher airway pressures, ePTFE was superior to bovine pericardium. CONCLUSION: Staple line reinforcement with either material protects against air leak. Patients at risk for elevated airway pressures and/or postoperative ventilator support should be considered for utilization of these staple reinforcing materials.