Attack of the PCR clones: Rates of clonality have little effect on RAD‐seq genotype calls

Interpretation of high‐throughput sequence data requires an understanding of how decisions made during bioinformatic data processing can influence results. One source of bias that is often cited is PCR clones (or PCR duplicates). PCR clones are common in restriction site‐associated sequencing (RAD‐seq) data sets, which are increasingly being used for molecular ecology. To determine the influence PCR clones and the bioinformatic handling of clones have on genotyping, we evaluate four RAD‐seq data sets. Data sets were compared before and after clones were removed to estimate the number of clones present in RAD‐seq data, quantify how often the presence of clones in a data set causes genotype calls to change compared to when clones were removed, investigate the mechanisms that lead to genotype call changes and test whether clones bias heterozygosity estimates. Our RAD‐seq data sets contained 30%–60% PCR clones, but 95% of RAD‐tags had five or fewer clones. Relatively few genotypes changed once clones were removed (5%–10%), and the vast majority of these changes (98%) were associated with genotypes switching from a called to no‐call state or vice versa. PCR clones had a larger influence on genotype calls in individuals with low read depth but appeared to influence genotype calls at all loci similarly. Removal of PCR clones reduced the number of called genotypes by 2% but had almost no influence on estimates of heterozygosity. As such, while steps should be taken to limit PCR clones during library preparation, PCR clones are likely not a substantial source of bias for most RAD‐seq studies.     

Using density surface models to estimate spatio-temporal changes in population densities and trend

Precise measures of population abundance and trend are needed for species conservation; these are most difficult to obtain for rare and rapidly changing populations. We compare uncertainty in densities estimated from spatio–temporal models with that from standard design-based methods. Spatio–temporal models allow us to target priority areas where, and at times when, a population may most benefit. Generalised additive models were fitted to a 31-year time series of point-transect surveys of an endangered Hawaiian forest bird, the Hawai‘i ‘ākepa Loxops coccineus. This allowed us to estimate bird densities over space and time. We used two methods to quantify uncertainty in density estimates from the spatio–temporal model: the delta method (which assumes independence between detection and distribution parameters) and a variance propagation method. With the delta method we observed a 52% decrease in the width of the design-based 95% confidence interval (CI), while we observed a 37% decrease in CI width when propagating the variance. We mapped bird densities as they changed across space and time, allowing managers to evaluate management actions. Integrating detection function modelling with spatio–temporal modelling exploits survey data more efficiently by producing finer-grained abundance estimates than are possible with design-based methods as well as producing more precise abundance estimates. Model-based approaches require switching from making assumptions about the survey design to assumptions about bird distribution. Such a switch warrants consideration. In this case the model-based approach benefits conservation planning through improved management efficiency and reduced costs by taking into account both spatial shifts and temporal changes in population abundance and distribution.     

Effects of thermal history on the responses to thermal stress of a large benthic foraminifera,Calcarina gaudichaudii

Marine ecosystems, particularly coastal environments, are rapidly changing due to anthropogenic impacts resulting in increased global climate change (ocean warming), ocean acidification, hypoxia, and eutrophication. On coral reefs, symbiont-bearing large benthic foraminifera (LBFs) can play a key role as reef constituents and carbonate producers, contributing up to 5% of reef-scale carbonate budgets. However, projected climate change, particularly ocean warming, has the potential to significantly alter the conditions in which marine organisms persist. While the response of LBFs to elevated thermal stress is well documented in laboratory studies, the potential influence of adaptation or acclimatization through prior environmental thermal history on this response remains largely unknown. In this study, specimens of Calcarina gaudichaudii, an LBF from the Penghu Islands, Taiwan, were collected from thermally variable intertidal and thermally stable subtidal (~ 6 m depth) environments representing thermal history. LBFs were then acclimated to laboratory conditions at ambient (25 °C) and elevated (28 °C) temperatures for three weeks, and subsequently exposed to control and heat stress treatments (25 °C, 28 °C, 30 °C, 33 °C) for an additional one week. Photosynthetic rates (determined through oxygen flux measurements) of C. gaudichaudii significantly decreased in specimens collected at subtidal depths acclimated at 25 °C when compared to those acclimated at 28 °C, whereas there was no effect of thermal history on respiration, indicating that symbiont and holobiont responses may differ in LBFs. Additionally, maximum photochemical efficiency (Fv/Fm) significantly decreased as a result of heat stress, although bleaching was not visually observed after one week. These results highlight the plastic responses of the algal microbiome and indicate that thermal history, acclimatization temperature, and heat stress interact to affect the physiological status of C. gaudichaudii. This study adds to the growing literature which highlights the larger implications of understanding thermal history as an important factor to consider to better understand how ecosystem processes (e.g., carbonate production) are altered on modern coral reefs.

     

Potassium channel blocker crafted by α-hairpinin scaffold engineering

The α-Hairpinins are a family of plant defense peptides with a common fold presenting two short α-helices stabilized by two invariant S–S-bridges. We have shown previously that substitution of just two amino acid residues in a wheat α-hairpinin Tk-AMP-X2 leads to Tk-hefu-2 that features specific affinity to voltage-gated potassium channels K_V1.3. Here, we utilize a combined molecular modeling approach based on molecular dynamics simulations and protein surface topography technique to improve the affinity of Tk-hefu-2 to K_V1.3 while preserving its specificity. An important advance of this work compared with our previous studies is transition from the analysis of various physicochemical properties of an isolated toxin molecule to its consideration in complex with its target, a membrane-bound ion channel. As a result, a panel of computationally designed Tk-hefu-2 derivatives was synthesized and tested against K_V1.3. The most active mutant Tk-hefu-10 showed a half-maximal inhibitory concentration of ～150 nM being >10 times more active than Tk-hefu-2 and >200 times more active than the original Tk-hefu. We conclude that α-hairpinins provide an attractive disulfide-stabilized scaffold for the rational design of ion channel inhibitors. Furthermore, the success rate can be considerably increased by the proposed “target-based” iterative strategy of molecular design.