Reduced representation approaches produce similar results to whole genome sequencing for some common phylogeographic analyses

When designing phylogeographic investigations researchers can choose to collect many different types of molecular markers, including mitochondrial genes or genomes, SNPs from reduced representation protocols, large sequence capture data sets, and even whole genomes. Given that the statistical power and accuracy of various analyses are expected to differ depending on both the type of marker and the amount of data collected, an exploration of the variance across methodological results as a function of marker type should provide valuable information to researchers. Here we collect mitochondrial Cytochrome b sequences, whole mitochondrial genomes, single nucleotide polymorphisms (SNP)s isolated using a genotype by sequencing (GBS) protocol, sequences from ultraconserved elements, and low-coverage nuclear genomes from the North American water vole (Microtus richardsoni). We estimate genetic distances, population genetic structure, and historical demography using data from each of these datasets and compare the results across markers. As anticipated, the results exhibit differences across marker types, particularly in terms of the resolution offered by different analyses. A cost-benefit analysis indicates that SNPs collected using a GBS protocol are the most cost-effective molecular marker, with inferences that mirror those collected from the whole genome data at a fraction of the cost per sample.

Highly-Selective Harvesting of (6,4) SWCNTs Using the Aqueous Two-Phase Extraction Method and Nonionic Surfactants

Monochiral single-walled carbon nanotubes (SWCNTs) are indispensable for advancing the technology readiness level of nanocarbon-based concepts. In recent times, many separation techniques have been developed to obtain specific SWCNTs from raw unsorted materials to catalyze the development in this area. This work presents how the aqueous two-phase extraction (ATPE) method can be enhanced for the straightforward isolation of (6,4) SWCNTs in one step. Introducing nonionic surfactant into the typically employed mixture of anionic surfactants, which drive the partitioning, is essential to increasing the ATPE system's resolution. A thorough analysis of the parameter space by experiments and modeling reveals the underlying interactions between SWCNTs, surfactants, and phase-forming agents, which drive the partitioning. Based on new insight gained on this front, a separation mechanism is proposed. Notably, the developed method is highly robust, which is proven by isolating (6,4) SWCNTs from several raw SWCNT materials, including SWCNT waste generated over the years in the laboratory.

Assessing model adequacy for Bayesian Skyline plots using posterior predictive simulation

Bayesian skyline plots (BSPs) are a useful tool for making inferences about demographic history. For example, researchers typically apply BSPs to test hypotheses regarding how climate changes have influenced intraspecific genetic diversity over time. Like any method, BSP has assumptions that may be violated in some empirical systems (e.g., the absence of population genetic structure), and the naïve analysis of data collected from these systems may lead to spurious results. To address these issues, we introduce P2C2M.Skyline, an R package designed to assess model adequacy for BSPs using posterior predictive simulation. P2C2M.Skyline uses a phylogenetic tree and the log file output from Bayesian Skyline analyses to simulate posterior predictive datasets and then compares this null distribution to statistics calculated from the empirical data to check for model violations. P2C2M.Skyline was able to correctly identify model violations when simulated datasets were generated assuming genetic structure, which is a clear violation of BSP model assumptions. Conversely, P2C2M.Skyline showed low rates of false positives when models were simulated under the BSP model. We also evaluate the P2C2M.Skyline performance in empirical systems, where we detected model violations when DNA sequences from multiple populations were lumped together. P2C2M.Skyline represents a user-friendly and computationally efficient resource for researchers aiming to make inferences from BSP.

MetaPop: a pipeline for macro- and microdiversity analyses and visualization of microbial and viral metagenome-derived populations

BACKGROUND

Microbes and their viruses are hidden engines driving Earth's ecosystems from the oceans and soils to humans and bioreactors. Though gene marker approaches can now be complemented by genome-resolved studies of inter-(macrodiversity) and intra-(microdiversity) population variation, analytical tools to do so remain scattered or under-developed.

RESULTS

Here, we introduce MetaPop, an open-source bioinformatic pipeline that provides a single interface to analyze and visualize microbial and viral community metagenomes at both the macro- and microdiversity levels. Macrodiversity estimates include population abundances and α- and β-diversity. Microdiversity calculations include identification of single nucleotide polymorphisms, novel codon-constrained linkage of SNPs, nucleotide diversity (π and θ), and selective pressures (pN/pS and Tajima's D) within and fixation indices (FST) between populations. MetaPop will also identify genes with distinct codon usage. Following rigorous validation, we applied MetaPop to the gut viromes of autistic children that underwent fecal microbiota transfers and their neurotypical peers. The macrodiversity results confirmed our prior findings for viral populations (microbial shotgun metagenomes were not available) that diversity did not significantly differ between autistic and neurotypical children. However, by also quantifying microdiversity, MetaPop revealed lower average viral nucleotide diversity (π) in autistic children. Analysis of the percentage of genomes detected under positive selection was also lower among autistic children, suggesting that higher viral π in neurotypical children may be beneficial because it allows populations to better "bet hedge" in changing environments. Further, comparisons of microdiversity pre- and post-FMT in autistic children revealed that the delivery FMT method (oral versus rectal) may influence viral activity and engraftment of microdiverse viral populations, with children who received their FMT rectally having higher microdiversity post-FMT. Overall, these results show that analyses at the macro level alone can miss important biological differences.

CONCLUSIONS

These findings suggest that standardized population and genetic variation analyses will be invaluable for maximizing biological inference, and MetaPop provides a convenient tool package to explore the dual impact of macro- and microdiversity across microbial communities. Video abstract.

Comparing B3LYP and B97 Dispersion-corrected Functionals for Studying Adsorption and Vibrational Spectra in Nitrogen Reduction

Electrochemical ammonia synthesis is being actively studied as a low temperature, low pressure alternative to the Haber-Bosch process. This work studied pure iridium as the catalyst for ammonia synthesis, following promising experimental results of Pt-Ir alloys. The characteristics studied include bond energies, bond lengths, spin densities, and free and adsorbed vibrational frequencies for the molecules N2 , N, NH, NH2 , and NH3 . Overall, these descriptive characteristics explore the use of dispersion-corrected density functional theory methods that can model N2 adsorption - the key reactant for electrochemical ammonia synthesis via transition metal catalysis. Specifically, three methods were tested: hybrid B3LYP, a dispersion-corrected form B3LYP-D3, and semi-empirical B97-D3. The latter semi-empirical method was explored to increase the accuracy obtained in vibrational analysis as well as reduce computational time. Two lattice surfaces, (111) and (100), were compared. The adsorption energies are stronger on (100) and follow the trend EB3LYP >EB3LYP-D3 >EB97-D3 on both surfaces.

FuSpot: a web-based tool for visual evaluation of fusion candidates

BACKGROUND

Gene fusions often occur in cancer cells and in some cases are the main driver of oncogenesis. Correct identification of oncogenic gene fusions thus has implications for targeted cancer therapy. Recognition of this potential has led to the development of a myriad of sequencing-based fusion detection tools. However, given the same input, many of these detectors will find different fusion points or claim different sets of supporting data. Furthermore, the rate at which these tools falsely detect fusion events in data varies greatly. This discrepancy between tools underscores the fact that computation algorithms still cannot perfectly evaluate evidence; especially when provided with small amounts of supporting data as is typical in fusion detection. We assert that when evidence is provided in an easily digestible form, humans are more proficient in identifying true positives from false positives.

RESULTS

We have developed a web tool that, given the genomic coordinates of a candidate fusion breakpoint, will extract fusion and non-fusion reads adjacent to the fusion point from partner transcripts, and color code reads by transcript origin and read orientation for ease of intuitive inspection by the user. Fusion partner transcript read alignments are performed using a novel variant of the Smith-Waterman algorithm.

CONCLUSIONS

Combined with dynamic filtering parameters, the visualization provided by our tool introduces a powerful new investigative step that allows researchers to comprehensively evaluate fusion evidence. Additionally, this allows quick identification of false positives that may deceive most fusion detectors, thus eliminating unnecessary gene fusion validation. We apply our visualization tool to publicly available datasets and provide examples of true as well as false positives reported by open source fusion detection tools.

Novel chloroacetamido compound CWR-J02 is an anti-inflammatory glutaredoxin-1 inhibitor

Glutaredoxin (Grx1) is a ubiquitously expressed thiol-disulfide oxidoreductase that specifically catalyzes reduction of S-glutathionylated substrates. Grx1 is known to be a key regulator of pro-inflammatory signaling, and Grx1 silencing inhibits inflammation in inflammatory disease models. Therefore, we anticipate that inhibition of Grx1 could be an anti-inflammatory therapeutic strategy. We used a rapid screening approach to test 504 novel electrophilic compounds for inhibition of Grx1, which has a highly reactive active-site cysteine residue (pKa 3.5). From this chemical library a chloroacetamido compound, CWR-J02, was identified as a potential lead compound to be characterized. CWR-J02 inhibited isolated Grx1 with an IC₅₀ value of 32 μM in the presence of 1 mM glutathione. Mass spectrometric analysis documented preferential adduction of CWR-J02 to the active site Cys-22 of Grx1, and molecular dynamics simulation identified a potential non-covalent binding site. Treatment of the BV2 microglial cell line with CWR-J02 led to inhibition of intracellular Grx1 activity with an IC₅₀ value (37 μM). CWR-J02 treatment decreased lipopolysaccharide-induced inflammatory gene transcription in the microglial cells in a parallel concentration-dependent manner, documenting the anti-inflammatory potential of CWR-J02. Exploiting the alkyne moiety of CWR-J02, we used click chemistry to link biotin azide to CWR-J02-adducted proteins, isolating them with streptavidin beads. Tandem mass spectrometric analysis identified many CWR-J02-reactive proteins, including Grx1 and several mediators of inflammatory activation. Taken together, these data identify CWR-J02 as an intracellularly effective Grx1 inhibitor that may elicit its anti-inflammatory action in a synergistic manner by also disabling other pro-inflammatory mediators. The CWR-J02 molecule provides a starting point for developing more selective Grx1 inhibitors and anti-inflammatory agents for therapeutic development.

Discovery of riblets in a bird beak (Rynchops) for low fluid drag

Riblet structures found on fast-swimming shark scales, such as those found on a mako shark, have been shown to reduce fluid drag. In previous experimental and modelling studies, riblets have been shown to provide drag reduction by lifting vortices formed in turbulent flow, decreasing overall shear stresses. Skimmer birds (Rynchops) are the only birds to catch fish in flight by flying just above the water surface with a submerged beak to fish for food. Because they need to quickly catch prey, reducing drag on their beak is advantageous. For the first time, riblet structures found on the beak of the skimmer bird have been studied experimentally and computationally for low fluid drag properties. In this study, skimmer replicas were studied for drag reduction through pressure drop in closed-channel, turbulent water flow. Pressure drop measurements are compared for black and yellow skimmer beaks in two configurations, and mako shark skin. In addition, two configurations of skimmer beak were modelled to compare drag properties and vortex structures. Results are discussed, and a conceptual model is presented to explain a possible drag reduction mechanism in skimmers.This article is part of the themed issue 'Bioinspired hierarchically structured surfaces for green science'.