Robust quantification of polymerase chain reactions using global fitting

Variable Region Sequences Influence 16S rRNA Performance

16S rRNA gene sequences are commonly analyzed for taxonomic and phylogenetic studies because they contain variable regions that can help distinguish different genera. However, intra-genus distinction using variable region homology is often impossible due to the high overall sequence identities among closely related species, even though some residues may be conserved within respective species. Using a computational method that included the allelic diversity within individual genomes, we discovered that certain Escherichia and Shigella species can be distinguished by a multi-allelic 16S rRNA variable region single nucleotide polymorphism (SNP). To evaluate the performance of 16S rRNAs with altered variable regions, we developed an in vivo system that measures the acceptance and distribution of variant 16S rRNAs into a large pool of natural versions supporting normal translation and growth. We found that 16S rRNAs containing evolutionarily disparate variable regions were underpopulated both in ribosomes and in active translation pools, even for an SNP. Overall, this study revealed that variable region sequences can substantially influence the performance of 16S rRNAs and that this biological constraint can be leveraged to justify refining taxonomic assignments of variable region sequence data. IMPORTANCE This study reevaluates the notion that 16S rRNA gene variable region sequences are uninformative for intra-genus classification and that single nucleotide variations within them have no consequence to strains that bear them. We demonstrated that the performance of 16S rRNAs in Escherichia coli can be negatively impacted by sequence changes in variable regions, even for single nucleotide changes that are native to closely related Escherichia and Shigella species; thus, biological performance is likely constraining the evolution of variable regions in bacteria. Further, the native nucleotide variations we tested occur in all strains of their respective species and across their multiple 16S rRNA gene copies, suggesting that these species evolved beyond what would be discerned from a consensus sequence comparison. Therefore, this work also reveals that the multiple 16S rRNA gene alleles found in most bacteria can provide more informative phylogenetic and taxonomic detail than a single reference allele.

Robust and Precise Quantitative Real-Time Polymerase Chain Reaction with an Amplification Efficiency-Aware Reaction Kinetics Model

Quantitative real-time PCR (qPCR) is a method extensively used in nucleic acid testing for plants and animals. During the coronavirus disease 2019 (COVID-19) pandemic, high-precision qPCR analysis was urgently needed since quantitative results obtained from conventional qPCR methods were not accurate and precise, causing misdiagnoses and high rates of false-negative. To achieve more accurate results, we propose a new qPCR data analysis method with an amplification efficiency-aware reaction kinetics model (AERKM). Our reaction kinetics model (RKM) mathematically describes the tendency of the amplification efficiency during the whole qPCR process inferred by biochemical reaction dynamics. Amplification efficiency (AE) was introduced to rectify the fitted data so as to match the real reaction process for individual tests, thus reducing errors. The 5-point 10-fold gradient qPCR tests of 63 genes have been verified. The results of a 0.9% slope bias and an 8.2% ratio bias using AERKM exceed 4.1 and 39.4%, respectively, of the best performance of existing models, which demonstrates higher precision, less fluctuation, and better robustness among different nucleic acids. AERKM also provides a better understanding of the real qPCR process and gives insights into the detection, treatment, and prevention of severe diseases.

Current Advances in Nanotechnology for the Next Generation of Sequencing (NGS)

This communication aims at discussing strategies based on developments from nanotechnology focused on the next generation of sequencing (NGS). In this regard, it should be noted that even in the advanced current situation of many techniques and methods accompanied with developments of technology, there are still existing challenges and needs focused on real samples and low concentrations of genomic materials. The approaches discussed/described adopt spectroscopical techniques and new optical setups. PCR bases are introduced to understand the role of non-covalent interactions by discussing about Nobel prizes related to genomic material detection. The review also discusses colorimetric methods, polymeric transducers, fluorescence detection methods, enhanced plasmonic techniques such as metal-enhanced fluorescence (MEF), semiconductors, and developments in metamaterials. In addition, nano-optics, challenges linked to signal transductions, and how the limitations reported in each technique could be overcome are considered in real samples. Accordingly, this study shows developments where optical active nanoplatforms generate signal detection and transduction with enhanced performances and, in many cases, enhanced signaling from single double-stranded deoxyribonucleic acid (DNA) interactions. Future perspectives on miniaturized instrumentation, chips, and devices aimed at detecting genomic material are analyzed. However, the main concept in this report derives from gained insights into nanochemistry and nano-optics. Such concepts could be incorporated into other higher-sized substrates and experimental and optical setups.

Overexpression of ABI5 Binding Proteins Suppresses Inhibition of Germination Due to Overaccumulation of DELLA Proteins

Abscisic acid (ABA) and gibberellic acid (GA) antagonistically regulate many aspects of plant growth, including seed dormancy and germination. The effects of these hormones are mediated by a complex network of positive and negative regulators of transcription. The DELLA family of proteins repress GA response, and can promote an ABA response via interactions with numerous regulators, including the ABA-insensitive (ABI) transcription factors. The AFP family of ABI5 binding proteins are repressors of the ABA response. This study tested the hypothesis that the AFPs also interact antagonistically with DELLA proteins. Members of these protein families interacted weakly in yeast two-hybrid and bimolecular fluorescence complementation studies. Overexpression of AFPs in sleepy1, a mutant that over-accumulates DELLA proteins, suppressed DELLA-induced overaccumulation of storage proteins, hyperdormancy and hypersensitivity to ABA, but did not alter the dwarf phenotype of the mutant. The interaction appeared to reflect additive effects of the AFPs and DELLAs, consistent with action in convergent pathways.

Size-dependent resistance to amoebic gill disease in naïve Atlantic salmon (Salmo salar)

Amoebic gill disease, caused by the protozoan ectoparasite Neoparamoeba perurans, remains a significant threat to commercial Atlantic salmon aquaculture operations worldwide, despite partial control afforded by selective breeding and therapeutic intervention. Anecdotal reports from commercial producers suggest that historically, smaller Atlantic salmon smolts are more susceptible to AGD than larger smolts. Here, large (>350 g) and small (<200 g) commercially sourced, AGD-naïve Atlantic salmon cohorts were experimentally exposed to 50 N. perurans trophozoites L^-1 without intervention. Progression and severity of AGD in challenged cohorts was evaluated through gill pathology, using gill score and histological examination, and quantification of gill-associated amoebae burden using qPCR. To determine the potential basis for differences in AGD susceptibility between cohorts, transcriptome analysis was conducted using RNA extracted from whole gill arches. Overall, the large Atlantic salmon cohort had significantly lower gill parasite burdens and reduced AGD-related gross pathology compared to the small cohort. Relative gill load of N. perurans appeared to be proportional to gill score in both size classes, with larger smolts typically observed to have comparatively reduced parasite burdens at a given gill score. Moreover, comparison between gene expression profiles of large and small smolts highlighted upregulation of genes consistent with elevated immune activity in large smolts. Combined, the results presented here provide strong evidence of size-dependent resistance to AGD in AGD-naïve Atlantic salmon.

Estimating Real-Time qPCR Amplification Efficiency from Single-Reaction Data

Methods for estimating the qPCR amplification efficiency E from data for single reactions are tested on six multireplicate datasets, with emphasis on their performance as a function of the range of cycles n₁–n₂ included in the analysis. The two-parameter exponential growth (EG) model that has been relied upon almost exclusively does not allow for the decline of E(n) with increasing cycle number n through the growth region and accordingly gives low-biased estimates. Further, the standard procedure of “baselining”—separately estimating and subtracting a baseline before analysis—leads to reduced precision. The three-parameter logistic model (LRE) does allow for such decline and includes a parameter E₀ that represents E through the baseline region. Several four-parameter extensions of this model that accommodate some asymmetry in the growth profiles but still retain the significance of E₀ are tested against the LRE and EG models. The recursion method of Carr and Moore also describes a declining E(n) but tacitly assumes E₀ = 2 in the baseline region. Two modifications that permit varying E₀ are tested, as well as a recursion method that directly fits E(n) to a sigmoidal function. All but the last of these can give E₀ estimates that agree fairly well with calibration-based estimates but perform best when the calculations are extended to only about one cycle below the first-derivative maximum (FDM). The LRE model performs as well as any of the four-parameter forms and is easier to use. Its proper implementation requires fitting to it plus a suitable baseline function, which typically requires four–six adjustable parameters in a nonlinear least-squares fit.

Experimental exposure to low concentrations of Neoparamoeba perurans induces amoebic gill disease in Atlantic salmon

Amoebic gill disease (AGD) is a significant issue in Atlantic salmon mariculture. Research on the development of treatments or vaccines uses experimental challenges where salmon is exposed to amoebae concentrations ranging from 500 to 5,000/L. However, the water concentrations of N. perurans on affected salmon farms are much lower. The lowest concentration of N. perurans previously reported to cause AGD was 10/L. Here, we report that concentrations as low as 0.1/L of N. perurans can cause AGD. We propose that concentrations of N. perurans that reflect those measured on salmon farms should be used for future experimental challenges.

Exopolysaccharide anchoring creates an extreme resistance to sedimentation

By evolving strains of E. coli that hyper-resist sedimentation, we discovered an uncharacterized mechanism that bacteria can use to remain in suspension indefinitely without expending energy. This unusual phenotype was traced to the anchoring of long colanic acid polymers (CAP) that project from the cell surface. Although each characterized mutant activated this same mechanism, the genes responsible and the strengths of the phenotypes varied. Mutations in rcsC, lpp, igaA, or the yjbEFGH operon were sufficient to stimulate sedimentation resistance, while mutations altering the cps promoter, cdgI, or yjbF provided phenotypic enhancements. The sedimentation resistances changed in response to temperature, growth phase, and carbon source and each mutant exhibited significantly reduced biofilm formation. We discovered that the degree of colony mucoidy exhibited by these mutants was not related to the degree of Rcs pathways activation or to the amount of CAP that was produced; rather, it was related to the fraction of CAP that was shed as a true exopolysaccharide. Therefore, these and other mutations that activate this phenotype are likely to be absent from genetic screens that relied on centrifugation to harvest bacteria. We also found that this anchored CAP form is not linked to LPS cores and may not be attached to the outer membrane.IMPORTANCEBacteria can partition in aqueous environments between surface-dwelling, planktonic, sedimentary, and biofilm forms. Residence in each location provides an advantage depending on nutritional and environmental stresses and a community of a single species is often observed to be distributed throughout two or more of these niches. Another adaptive strategy is to produce an extracellular capsule, which provides an environmental shield for the microbe and can allow escape from predators and immune systems. We discovered that bacteria can either shed or stably anchor capsules to dramatically alter their propensity to sediment. The degree to which the bacteria anchor their capsule is controlled by a stress sensing system, suggesting that anchoring may be used as an adaptive response to severe environmental challenges.

Post-reclamation microbial diversity and functions in hexachlorocyclohexane (HCH) contaminated soil in relation to spontaneous HCH tolerant vegetation

The toxicity, volatility and persistence of the obsolete organochlorine pesticide hexachlorocyclohexane (HCH), makes reclamation of contaminated areas a priority for the health and welfare of neighboring human communities. Microbial diversity and functions and their relation to spontaneous vegetation in post-excavation situations, are essential indicators to consider in bioaugmentation or microbe-assisted phytoremediation strategies at field scale. Our study aimed to evaluate the effects of long-term HCH contamination on soil and plant-associated microbial communities, and whether contaminated soil has the potential to act as a bacterial inoculum in post-excavation bioremediation strategies. To scrutinize the role of vegetation, the potential nitrogen fixation of free-living and symbiotic diazotrophs of the legume Lotus tenuis was assessed as a measure of nutrient cycling functions in soil under HCH contamination. Potential nitrogen fixation was generally not affected by HCH, with the exception of lower nifH gene counts in excavated contaminated rhizospheres, most probably a short-term HCH effect on early bacterial succession in this compartment. HCH shaped microbial communities in long-term contaminated bulk soil, where we identified possible HCH tolerants such as Sphingomonas and Altererythrobacter. In L. tenuis rhizosphere, microbial community composition was additionally influenced by plant growth stage. Sphingobium and Massilia were the bacterial genera characteristic for HCH contaminated rhizospheres. Long-term HCH contamination negatively affected L. tenuis growth and development. However, root-associated bacterial community composition was driven solely by plant age, with negligible HCH effect. Results showed that L. tenuis acquired possible HCH tolerant bacteria such as the Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium clade, Sphingomonas, Massilia or Pantoea which could simultaneously offer plant growth promoting (PGP) benefits for the host. Finally, we identified an inoculum with possibly HCH tolerant, PGP bacteria transferred from the contaminated bulk soil to L. tenuis roots through the rhizosphere compartment, consisting of Mesorhizobium loti, Neorhizobium galegae, Novosphingobium lindaniclasticum, Pantoea agglomerans and Lysobacter bugurensis.

Feedback control of the Gpr161-Gαs-PKA axis contributes to basal Hedgehog repression in zebrafish

Hedgehog (Hh) ligands act as morphogens to direct patterning and proliferation during embryonic development. Protein kinase A (PKA) is a central negative regulator of Hh signalling, and in the absence of Hh ligands, PKA activity prevents inappropriate expression of Hh target genes. The orphan G-protein-coupled receptor Gpr161 contributes to the basal Hh repression machinery by activating PKA. Gpr161 acts as an A-kinase-anchoring protein, and is itself phosphorylated by PKA, but the functional significance of PKA phosphorylation of Gpr161 in the context of Hh signalling remains unknown. Here, we show that loss of Gpr161 in zebrafish leads to constitutive activation of medium and low, but not maximal, levels of Hh target gene expression. Furthermore, we find that PKA phosphorylation-deficient forms of Gpr161, which we show directly couple to G_αs, display an increased sensitivity to Shh, resulting in excess high-level Hh signalling. Our results suggest that PKA feedback-mediated phosphorylation of Gpr161 may provide a mechanism for fine-tuning Gpr161 ciliary localisation and PKA activity.

52 years of ecological restoration following a major disturbance by opencast lignite mining does not reassemble microbiome structures of the original arable soils

Opencast mining for lignite continuously creates areas of land that require restoration. Here we applied a chronosequence approach to investigate the development of soil bacterial communities during 52 years as influenced by the restoration process and subsequent changes in soil physico-chemical conditions starting from the initial reclamation of the sites. By comparison with the unaffected soils near the mine, we were able to address the question if soil bacterial communities have reached a steady state within 52 years, which is comparable to the original soil. Our study revealed three distinct phases of the restoration process, each with a specific bacterial community composition. The effect size of these changes was similar to the one observed for seasonal dynamics at our sites. At the beginning of the restoration process Flavobacteriaceae, Cytophagaceae and Sphingobacteriaceae were found as typical members of the bacterial community as well as Rhizobiales as a result of the cultivation of alfalfa on the restored plots. At later stage the families Peptostreptococcaceae, Desulfurellaceae as well as Streptomycetaceae increased in relative abundance and became dominant members of the bacterial community. Even though overall bacterial abundance and richness exhibited values comparable to the original soil already 5 years after the start of the restoration process, main responder analyses reveal differences in the bacterial community structure even 52 years after the start of the restoration process. Mostly Nitrospirae were reduced in abundance in the soils restored for 52 years compared to the original soils. To broaden the significance of our study, we compared our data bioinformatically with published results from other restored areas, which were previously affected by opencast mining. Despite different durations of the different restoration phase, we could observe a large degree of conformity when bacterial patterns of succession were compared indicating common modes of action of ecological restoration tools for bacterial communities.

Ribosome collisions alter frameshifting at translational reprogramming motifs in bacterial mRNAs

Ribosomes move along mRNAs in 3-nucleotide steps as they interpret codons that specify which amino acid is required at each position in the protein. There are multiple examples of genes with DNA sequences that do not match the produced proteins because ribosomes move to a new reading frame in the message before finishing translation (so-called frameshifting). This report shows that, when ribosomes stall at mRNA regions prone to cause frameshifting events, trailing ribosomes that collide with them can significantly change the outcome and potentially regulate protein production. This work highlights the principle that biological macromolecules do not function in isolation, and it provides an example of how physical interactions between neighboring complexes can be used to augment their performance.

Translational frameshifting involves the repositioning of ribosomes on their messages into decoding frames that differ from those dictated during initiation. Some messenger RNAs (mRNAs) contain motifs that promote deliberate frameshifting to regulate production of the encoded proteins. The mechanisms of frameshifting have been investigated in many systems, and the resulting models generally involve single ribosomes responding to stimulator sequences in their engaged mRNAs. We discovered that the abundance of ribosomes on messages containing the IS3, dnaX, and prfB frameshift motifs significantly influences the levels of frameshifting. We show that this phenomenon results from ribosome collisions that occur during translational stalling, which can alter frameshifting in both the stalled and trailing ribosomes. Bacteria missing ribosomal protein bL9 are known to exhibit a reduction in reading frame maintenance and to have a strong dependence on elongation factor P (EFP). We discovered that ribosomes lacking bL9 become compacted closer together during collisions and that the E-sites of the stalled ribosomes appear to become blocked, which suggests subsequent transpeptidation in transiently stalled ribosomes may become compromised in the absence of bL9. In addition, we determined that bL9 can suppress frameshifting of its host ribosome, likely by regulating E-site dynamics. These findings provide mechanistic insight into the behavior of colliding ribosomes during translation and suggest naturally occurring frameshift elements may be regulated by the abundance of ribosomes relative to an mRNA pool.

qPCR data analysis: Better results through iconoclasm

The standard approach for quantitative estimation of genetic materials with qPCR is calibration with known concentrations for the target substance, in which estimates of the quantification cycle (C_q ) are fitted to a straight-line function of log(N ₀), where N ₀ is the initial number of target molecules. The location of C_q for the unknown on this line then yields its N ₀. The most widely used definition for C_q is an absolute threshold that falls in the early growth cycles. This usage is flawed as commonly implemented: threshold set very close to the baseline level, which is estimated separately, from designated "baseline cycles." The absolute threshold is especially poor for dealing with the scale variability often observed for growth profiles. Scale-independent markers, like the first derivative maximum (FDM) and a relative threshold (C_r ) avoid this problem. We describe improved methods for estimating these and other C_q markers and their standard errors, from a nonlinear algorithm that fits growth profiles to a 4-parameter log-logistic function plus a baseline function. Further, by examining six multidilution, multireplicate qPCR data sets, we find that nonlinear expressions are often preferred statistically for the dependence of C_q on log(N ₀). This means that the amplification efficiency E depends on N ₀, in violation of another tenet of qPCR analysis. Neglect of calibration nonlinearity leads to biased estimates of the unknown. By logic, E estimates from calibration fitting pertain to the earliest baseline cycles, not the early growth cycles used to estimate E from growth profiles for single reactions. This raises concern about the use of the latter in lengthy extrapolations to estimate N ₀. Finally, we observe that replicate ensemble standard deviations greatly exceed predictions, implying that much better results can be achieved from qPCR through better experimental procedures, which likely include reducing pipette volume uncertainty.

Synergistic osmoregulatory dysfunction during salmon lice (Lepeophtheirus salmonis) and infectious hematopoietic necrosis virus co-infection in sockeye salmon (Oncorhynchus nerka) smolts

While co-infections are common in both wild and cultured fish, knowledge of the interactive effects of multiple pathogens on host physiology, gene expression and immune response is limited. To evaluate the impact of co-infection on host survival, physiology and gene expression, sockeye salmon Oncorhynchus nerka smolts were infected with the salmon louse Lepeophtheirus salmonis (V-/SL+), infectious hematopoietic necrosis virus (IHNV; V+/SL-), both (V+/SL+), or neither (V-/SL-). Survival in the V+/SL+ group was significantly lower than the V-/SL- and V-/SL+ groups (p = 0.024). Co-infected salmon had elevated osmoregulatory indicators and lowered haematocrit values as compared to the uninfected control. Expression of 12 genes associated with the host immune response was analysed in anterior kidney and skin. The only evidence of L. salmonis-induced modulation of the host antiviral response was down-regulation of mhc I although the possibility of modulation cannot be ruled out for mx-1 and rsad2. Co-infection did not influence the expression of genes associated with the host response to L. salmonis. Therefore, we conclude that the reduced survival in co-infected sockeye salmon resulted from the osmoregulatory consequences of the sea lice infections which were amplified due to infection with IHNV.

Pairwise efficiency: a new mathematical approach to qPCR data analysis increases the precision of the calibration curve assay

Background

The real-time quantitative polymerase chain reaction (qPCR) is routinely used for quantification of nucleic acids and is considered the gold standard in the field of relative nucleic acid measurements. The efficiency of the qPCR reaction is one of the most important parameters in data analysis in qPCR experiments. The Minimum Information for publication of Quantitative real-time PCR Experiments (MIQE) guidelines recommends the calibration curve as the method of choice for estimation of qPCR efficiency. The precision of this method has been reported to be between SD = 0.007 (three replicates) and SD = 0.022 (no replicates).

Results

In this article, we present a novel approach to the analysis of qPCR data which has been obtained by running a dilution series. Unlike previously developed methods, our method, Pairwise Efficiency, involves a new formula that describes pairwise relationships between data points on separate amplification curves and thus enables extensive statistics. The comparison of Pairwise Efficiency with the calibration curve by Monte Carlo simulation shows the two-folds improvement in the precision of estimations of efficiency and gene expression ratios on the same dataset.

Conclusions

The Pairwise Efficiency nearly doubles the precision in qPCR efficiency determinations compared to standard calibration curve method. This paper demonstrates that applications of combinatorial treatment of data provide the improvement of the determination.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2911-5) contains supplementary material, which is available to authorized users.

An annual cycle of gene regulation in the red-legged salamander mental gland: from hypertrophy to expression of rapidly evolving pheromones

Background

Cell differentiation is mediated by synchronized waves of coordinated expression for hundreds to thousands of genes, and must be regulated to produce complex tissues and phenotypes. For many animal species, sexual selection has driven the development of elaborate male ornaments, requiring sex-specific differentiation pathways. One such male ornament is the pheromone-producing mental gland of the red-legged salamander (Plethodon shermani). Mental gland development follows an annual cycle of extreme hypertrophy, production of pheromones for the ~ 2 month mating season, and then complete resorption before repeating the process in the following year. At the peak of the mating season, the transcriptional and translational machinery of the mental gland are almost exclusively redirected to the synthesis of rapidly evolving pheromones. Of these pheromones, Plethodontid Modulating Factor (PMF) has experienced an unusual history: following gene duplication, the protein coding sequence diversified from positive sexual selection while the untranslated regions have been conserved by purifying selection. The molecular underpinnings that bridge the processes of gland hypertrophy, pheromone synthesis, and conservation of the untranslated regions remain to be determined.

Results

Using Illumina sequencing, we prepared a de novo transcriptome of the mental gland at six stages of development. Differential expression analysis and immunohistochemistry revealed that the mental gland initially adopts a highly proliferative, almost tumor-like phenotype, followed by a rapid increase in pheromone mRNA and protein. One likely player in this transition is Cold Inducible RNA Binding Protein (CIRBP), which selectively and cooperatively binds the highly conserved PMF 3′ UTR. CIRBP, along with other proteins associated with stress response, have seemingly been co-opted to aid in mental gland development by helping to regulate pheromone synthesis.

Conclusions

The P. shermani mental gland utilizes a complex system of transcriptional and post-transcriptional gene regulation to facilitate its hypertrophication and pheromone synthesis. The data support the evolutionary interplay of coding and noncoding segments in rapid gene evolution, and necessitate the study of co-evolution between pheromone gene products and their transcriptional/translational regulators. Additionally, the mental gland could be a powerful emerging model of regulated tissue proliferation and subsequent resorption within the dermis and share molecular links to skin cancer biology.

Electronic supplementary material

The online version of this article (10.1186/s12861-019-0190-z) contains supplementary material, which is available to authorized users.

A highly sensitive, non-invasive qPCR-based strategy for direct quantification of Yersinia ruckeri in fish faeces

Finfish with asymptomatic Yersinia ruckeri infections pose a major risk as they can transmit the pathogen and cause clinical outbreaks in stock populations. Current tools have insufficient quantitative ability for accurately detecting the trace levels of Y. ruckeri typically associated with asymptomatic infection, necessitate invasive or lethal sampling, or require long processing times. This study presents a highly sensitive qPCR-based method, targeting part of the Y. ruckeri 16S rRNA sequence, that is capable of detecting extremely low levels of Y. ruckeri in noninvasively collected faecal samples. Quantitative precision and accuracy of faecal sample analysis was consistent, despite the complexity of the faecal matrix. The assay demonstrated linearity over a six log-wide dynamic range. Its limit of detection (LOD) and limit of quantification (LOQ) were 4 and 10 copies of the target sequence, respectively. Sensitivity of the assay was comparable to other qPCR-based methods without requiring invasive or lethal sampling. Applicability as a screening strategy was tested using passively collected faecal samples. Asymptomatic Y. ruckeri infection was detected in all samples, although none of the fish exhibited overt infection. This method will be beneficial for finfish disease management if developed further as a noninvasive, screening tool against asymptomatic Y. ruckeri infection.

Development of a portable electrochemical loop mediated isothermal amplification (LAMP) device for detection of hepatitis B virus

The objective of this study was to develop a simple, inexpensive prototype device for rapid detection of hepatitis B virus (HBV). The device was able to simultaneously amplify, detect and quantify the target HBV DNA. The system was fabricated from a custom-made electrochemical set-up of which the temperature was thermostatically controlled by a water bath. Real-time monitoring of HBV DNA was accomplished by measuring the response of redox indicator in the reaction mixture. Concentration of HBV DNA in the samples was determined from the peak high ratio (PHR) and threshold time relationship. The signal was processed by sigmoidal model fitting to enhance the accuracy of the results. Key parameters including concentrations of redox indicator and reaction temperatures were optimized. Sensitivity and specificity of the method toward HBV DNA were evaluated. The prototype was capable of real-time amplification and detection of HBV DNA with concentration as low as 6.18 fg μl⁻¹. The test showed high specificity against HBV DNA. The system was also able to detect HBV positive serum directly with simple thermal pretreatment instead of tedious DNA extraction. The electrochemical set-up was compatible with microfluidic platforms and can be readily adapted for efficient and high throughput point-of-care (POC) diagnosis of HBV.

A novel prototype device using LAMP and electrochemical drop cell set-up for rapid detection of hepatitis B virus.

Savitzky-Golay smoothing and differentiation for polymerase chain reaction quantification

In quantitative PCR (qPCR), replicates can minimize the impact of intra-assay variation; however, inter-assay variations must be minimized to obtain a robust quantification method. The method proposed in this study uses Savitzky-Golay smoothing and differentiation (SGSD) to identify a derivative-maximum-based cycle of quantification. It does not rely on curve modeling, as is the case with many existing techniques. PCR fluorescence data sets challenged for inter-assay variations (different thermocycler units, different reagents batches, different operators, different standard curves, and different labs) were used for the evaluation. The algorithm was compared with a four-parameter logistic model (4PLM) method, the C_y0 method, and the threshold method. The SGSD method compared favourably with all methods in terms of inter-assay variation. SGSD was statistically different from the 4PLM (P = 0.03), C_y0 (P = 0.05), and threshold (P = 0.004) methods on relative error comparison basis. For intra-assay variations, SGSD outperformed the threshold method (P = 0.005) and equalled the 4PLM and C_y0 methods (P > 0.05) on relative error basis. Our results demonstrate that the SGSD method could potentially be an alternative to sigmoid modeling based methods (4PLM and C_y0) when PCR data are challenged for inter-assay variations.

Gene expression signature predicts human islet integrity and transplant functionality in diabetic mice

There is growing evidence that transplantation of cadaveric human islets is an effective therapy for type 1 diabetes. However, gauging the suitability of islet samples for clinical use remains a challenge. We hypothesized that islet quality is reflected in the expression of specific genes. Therefore, gene expression in 59 human islet preparations was analyzed and correlated with diabetes reversal after transplantation in diabetic mice. Analysis yielded 262 differentially expressed probesets, which together predict islet quality with 83% accuracy. Pathway analysis revealed that failing islet preparations activated inflammatory pathways, while functional islets showed increased regeneration pathway gene expression. Gene expression associated with apoptosis and oxygen consumption showed little overlap with each other or with the 262 probeset classifier, indicating that the three tests are measuring different aspects of islet cell biology. A subset of 36 probesets surpassed the predictive accuracy of the entire set for reversal of diabetes, and was further reduced by logistic regression to sets of 14 and 5 without losing accuracy. These genes were further validated with an independent cohort of 16 samples. We believe this limited number of gene classifiers in combination with other tests may provide complementary verification of islet quality prior to their clinical use.

Tetramethylpyrazine attenuates periorbital allodynia and neuroinflammation in a model of traumatic brain injury

Background

Traumatic brain injury (TBI) is a public health issue. As the major complaint in 51% of TBI patients, chronic pain is an important aspect in TBI treatment. Tetramethylpyrazine (TMP) is an important compound in Ligustrazine, an analgesic drug in traditional Chinese medicine, but its potential in relieving pain symptom in TBI has not been tested. We established a TBI mouse model with controlled cortical impact (CCI), and measured periorbital hypersensitivity with von Frey monofilaments. We examined activated microglia and astrocytes and the levels of substance P (SP) and inducible isoform of nitric oxide synthase (iNOS) with immunohistochemistry, measured mRNA and protein levels of proinflammatory cytokines with qPCR and enzyme-linked immunosorbent assay, respectively. Western blot was employed to detect molecules in NF-κB signaling pathway.

Results

TMP significantly attenuated periorbital hypersensitivity in TBI mice. Within 3 days after CCI, TMP attenuated activation of microglia and astrocytes, levels of SP, iNOS, and CGRP in trigeminal pathway, and levels of proinflammatory cytokines (including IL-6, TNF-α, IL-12). In isolated microglia, TMP attenuated the effects of lipopolysaccharide on the phosphorylation of cytoplasmic IKKα/β and IKB-α, and levels of nucleic p65.

Conclusion

TMP reversed periorbital hypersensitivity by limiting neuroinflammation at the primary stage of TBI, and could be a promising drug for pain treatment in TBI.

The effect of engineered nanotopography of electrospun microfibers on fiber rigidity and macrophage cytokine production

Currently, it is unknown how the mechanical properties of electrospun fibers, and the presentation of surface nanotopography influence macrophage gene expression and protein production. By further elucidating how specific fiber properties (mechanical properties or surface properties) alter macrophage behavior, it may be possible to create electrospun fiber scaffolds capable of initiating unique cellular and tissue responses. In this study, we determined the elastic modulus and rigidity of fibers with varying topographies created by finely controlling humidity and including a non-solvent during electrospinning. In total,five fiber scaffold types were produced. Analysis of fiber physical properties demonstrated no change in fiber diameter amongst the five different fiber groups. However, the four different fibrous scaffolds with nanopits or divots each possessed different numbers of pits with different nanoscale dimensions. Unpolarized bone marrow derived murine macrophages (M0), macrophages polarized towards a pro-inflammatory phenotype (M1), or macrophages polarized towards anti-inflammatory phenotype (M2b) were placed onto each of the scaffolds and cytokine RNA expression and protein production were analyzed. Specific nanotopographies did not appreciably alter cytokine production from undifferentiated macrophages (M0) or anti-inflammatory macrophages (M2b), but a specific fiber (with many small pits) did increase IL-12 transcript and IL-12 protein production compared to fibers with small divots. When analyzing the mechanical properties between fibers with divots or with many small pits,divoted fibers possessed similar elastic moduli but different stiffness values. In total,we present techniques capable of creating unique electrospun fibers. These unique fibers have varying fiber mechanical characteristics and modestly modulate macrophage cytokine expression.

Tramadol differentially regulates M1 and M2 macrophages from human umbilical cord blood

Tramadol is an analgesic drug and relieves pain through activating μ-opioid receptors and inhibiting serotonin and noradrenaline reuptake. Emerging evidence shows that it also stimulates immune cells, including NK cells, splenocytes, and lymphocytes, and elevates IL-2 production. However, it remains unknown whether and how tramadol directly affects macrophages. To answer these questions, we collected human umbilical cord blood, isolated macrophages, and examined their responses to tramadol. Although tramadol did not alter resting macrophages and the antigen-presenting function in lipopolysaccharide-activated macrophages, it regulated M1 and M2 macrophages, which are, respectively, transformed by IFN-γ and IL-4. Interestingly, tramadol inhibits production and secretion of cytokines in M1 macrophages, but facilitates the production of inflammation-responding molecules, synthesized in M2 macrophages. We also found that STAT6 cascade pathway in M2 macrophages was significantly enhanced by tramadol. Therefore, this study reveals that tramadol regulates inflammation by inhibiting M1 macrophages (killing process), but promoting the function of M2 macrophages (healing process).

ABI5-binding proteins (AFPs) alter transcription of ABA-induced genes via a variety of interactions with chromatin modifiers

Overexpression of ABI5/ABF binding proteins (AFPs) results in extreme ABA resistance of seeds via multiple mechanisms repressing ABA response, including interactions with histone deacetylases and the co-repressor TOPLESS. Several ABI5/ABF binding proteins (AFPs) inhibit ABA response, resulting in extreme ABA resistance in transgenic Arabidopsis overexpression lines, but their mechanism of action has remained obscure. By analogy to the related Novel Interactor of JAZ (NINJA) protein, it was suggested that the AFPs interact with the co-repressor TOPLESS to inhibit ABA-regulated gene expression. This study shows that the AFPs that inhibit ABA response have intrinsic repressor activity in a heterologous system, which does not depend on the domain involved in the interaction with TOPLESS. This domain is also not essential for repressing ABA response in transgenic plants, but does contribute to stronger ABA resistance. Additional interactions between some AFPs and histone deacetylase subunits were observed in yeast two-hybrid and bimolecular fluorescence assays, consistent with a more direct mechanism of AFP-mediated repression of gene expression. Chemical inhibition of histone deacetylase activity by trichostatin A suppressed AFP effects on a small fraction of the ABI5-regulated genes tested. Collectively, these results suggest that the AFPs participate in multiple mechanisms modulating ABA response, including both TOPLESS-dependent and -independent chromatin modification.

System-specific periodicity in quantitative real-time polymerase chain reaction data questions threshold-based quantitation

Real-time quantitative polymerase chain reaction (qPCR) data are found to display periodic patterns in the fluorescence intensity as a function of sample number for fixed cycle number. This behavior is seen for technical replicate datasets recorded on several different commercial instruments; it occurs in the baseline region and typically increases with increasing cycle number in the growth and plateau regions. Autocorrelation analysis reveals periodicities of 12 for 96-well systems and 24 for a 384-well system, indicating a correlation with block architecture. Passive dye experiments show that the effect may be from optical detector bias. Importantly, the signal periodicity manifests as periodicity in quantification cycle (C_q) values when these are estimated by the widely applied fixed threshold approach, but not when scale-insensitive markers like first- and second-derivative maxima are used. Accordingly, any scale variability in the growth curves will lead to bias in constant-threshold-based C_qs, making it mandatory that workers should either use scale-insensitive C_qs or normalize their growth curves to constant amplitude before applying the constant threshold method.

De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection

Background

Piscine reovirus (PRV) has been associated with the serious disease known as Heart and Skeletal Muscle Inflammation (HSMI) in cultured Atlantic salmon Salmo salar in Norway. PRV is also prevalent in wild and farmed salmon without overt disease manifestations, suggesting multifactorial triggers or PRV variant-specific factors are required to initiate disease. In this study, we explore the head kidney transcriptome of Sockeye salmon Oncorhynchus nerka during early PRV infection to identify host responses in the absence of disease in hopes of elucidating mechanisms by which PRV may directly alter host functions and contribute to the development of a disease state. We further investigate the role of PRV as a coinfecting agent following superinfection with infectious hematopoietic necrosis virus (IHNV) – a highly pathogenic rhabdovirus endemic to the west coast of North America.

Results

Challenge of Sockeye salmon with PRV resulted in high quantities of viral transcripts to become present in the blood and kidney of infected fish without manifestations of disease. De novo transcriptome assembly of over 2.3 billion paired RNA-seq reads from the head kidneys of 36 fish identified more than 320,000 putative unigenes, of which less than 20 were suggested to be differentially expressed in response to PRV at either 2 or 3 weeks post challenge by DESeq2 and edgeR analysis. Of these, only one, Ependymin, was confirmed to be differentially expressed by qPCR in an expanded sample set. In contrast, IHNV induced substantial transcriptional changes (differential expression of > 20,000 unigenes) which included transcripts involved in antiviral and inflammatory response pathways. Prior infection with PRV had no significant effect on host responses to superinfecting IHNV, nor did host responses initiated by IHNV exposure influence increasing PRV loads.

Conclusions

PRV does not substantially alter the head kidney transcriptome of Sockeye salmon during early (2 to 3 week) infection and dissemination in a period of significant increasing viral load, nor does the presence of PRV change the host transcriptional response to an IHNV superinfection. Further, concurrent infections of PRV and IHNV do not appear to significantly influence the infectivity or severity of IHNV associated disease, or conversely, PRV load.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3196-y) contains supplementary material, which is available to authorized users.

A Tale of Tails: Dissecting the Enhancing Effect of Tailed Primers in Real-Time PCR

Non-specific tail sequences are often added to the 5’-terminus of primers to improve the robustness and overall performance of diagnostic assays. Despite the widespread use of tailed primers, the underlying working mechanism is not well understood. To address this problem, we conducted a detailed in vitro and in silico analysis of the enhancing effect of primer tailing on 2 well-established foot-and-mouth disease virus (FMDV) RT-qPCR assays using an FMDV reference panel. Tailing of the panFMDV-5UTR primers mainly affected the shape of the amplification curves. Modelling of the raw fluorescence data suggested a reduction of the amplification efficiency due to the accumulation of inhibitors. In depth analysis of PCR products indeed revealed the rapid accumulation of forward-primer derived artefacts. More importantly, tailing of the forward primer delayed artefacts formation and concomitantly restored the sigmoidal shape of the amplification curves. Our analysis also showed that primer tailing can alter utilisation patterns of degenerate primers and increase the number of primer variants that are able to participate in the reaction. The impact of tailed primers was less pronounced in the panFMDV-3D assay with only 5 out of 50 isolates showing a clear shift in Cq values. Sequence analysis of the target region of these 5 isolates revealed several mutations in the inter-primer region that extend an existing hairpin structure immediately downstream of the forward primer binding site. Stabilisation of the forward primer with either a tail sequence or cationic spermine units restored the sensitivity of the assay, which suggests that the enhancing effect in the panFMDV-3D assay is due to a more efficient extension of the forward primer. ur results show that primer tailing can alter amplification through various mechanisms that are determined by both the assay and target region. These findings expand our understanding of primer tailing and should enable a more targeted and efficient use of tailed primers.

The large ribosomal subunit protein L9 enables the growth of EF-P deficient cells and enhances small subunit maturation

The loss of the large ribosomal protein L9 causes a reduction in translation fidelity by an unknown mechanism. To identify pathways affected by L9, we identified mutants of E. coli that require L9 for fitness. In a prior study, we characterized L9-dependent mutations in the essential GTPase Der (EngA). Here, we describe a second class of L9-dependent mutations that either compromise or inactivate elongation factor P (EF-P, eIF5A in eukaryotes). Without L9, Δefp cells are practically inviable. Cell fractionation studies revealed that, in both the Der and EF-P mutant cases, L9's activity reduces immature 16S rRNA in 30S particles and partially restores the abundance of monosomes. Inspired by these findings, we discovered that L9 also enhances 16S maturation in wild-type cells. Surprisingly, although the amount of immature 16S in 30S particles was found to be elevated in ΔrplI cells, the amount in polysomes was low and inversely correlated with the immature 16S abundance. These findings provide an explanation for the observed fitness increases afforded by L9 in these mutants and reveal particular physiological conditions in which L9 becomes critical. Additionally, L9 may affect the partitioning of small subunits containing immature 16S rRNA.

A fast-and-robust profiler for improving polymerase chain reaction diagnostics

Polymerase chain reaction (PCR) is an in vitro technology in molecular genetics that progressively amplifies minimal copies of short DNA sequences in a fast and inexpensive manner. However, PCR performance is sensitive to suboptimal processing conditions. Compromised PCR conditions lead to artifacts and bias that downgrade the discriminatory power and reproducibility of the results. Promising attempts to resolve the PCR performance optimization issue have been guided by quality improvement tactics adopted in the past for industrial trials. Thus, orthogonal arrays (OAs) have been employed to program quick-and-easy structured experiments. Profiling of influences facilitates the quantification of effects that may counteract the detectability of amplified DNA fragments. Nevertheless, the attractive feature of reducing greatly the amount of work and expenditures by planning trials with saturated-unreplicated OA schemes is known to be relinquished in the subsequent analysis phase. This is because of an inherent incompatibility of ordinary multi-factorial comparison techniques to convert small yet dense datasets. Treating unreplicated-saturated data with either the analysis of variance (ANOVA) or regression models destroys the information extraction process. Both of those mentioned approaches are rendered blind to error since the examined effects absorb all available degrees of freedom. Therefore, in lack of approximating an experimental uncertainty, any outcome interpretation is rendered subjective. We propose a profiling method that permits the non-linear maximization of amplicon resolution by eliminating the necessity for direct error estimation. Our approach is distribution-free, calibration-free, simulation-free and sparsity-free with well-known power properties. It is also user-friendly by promoting rudimentary analytics. Testing our method on published amplicon count data, we found that the preponderant effect is the concentration of MgCl2 (p<0.05) followed by the primer content (p<0.1) whilst the effects due to either the content of the deoxynucleotide (dNTP) or DNA remained dormant (p>0.1). Comparison of the proposed method with other stochastic approaches is also discussed. Our technique is expected to have extensive applications in genetics and biotechnology where there is a demand for cheap, expedient, and robust information.

Group-effort applied research: expanding opportunities for undergraduate research through original, class-based research projects

Undergraduate research clearly enriches the educational development of participating students, but these experiences are limited by the inherent inefficiency of the standard one student-one mentor model for undergraduate research. Group-effort applied research (GEAR) was developed as a strategy to provide substantial numbers of undergraduates with meaningful research experiences. The GEAR curriculum delivers concept-driven lecture material and provides hands-on training in the context of an active research project from the instructor's laboratory. Because GEAR is structured as a class, participating students benefit from intensive, supervised research training that involves a built-in network of peer support and abundant contact with faculty mentors. The class format also ensures a relatively standardized and consistent research experience. Furthermore, meaningful progress toward a research objective can be achieved more readily with GEAR than with the traditional one student-one mentor model of undergraduate research because sporadic mistakes by individuals in the class are overshadowed by the successes of the group as a whole. Three separate GEAR classes involving three distinct research projects have been offered to date. In this article, we provide an overview of the GEAR format and review some of the recurring themes for GEAR instruction. We propose GEAR can serve as a template to expand student opportunities for life science research without sacrificing the quality of the mentored research experience.

Statistical uncertainty and its propagation in the analysis of quantitative polymerase chain reaction data: comparison of methods

Most methods for analyzing real-time quantitative polymerase chain reaction (qPCR) data for single experiments estimate the hypothetical cycle 0 signal y0 by first estimating the quantification cycle (Cq) and amplification efficiency (E) from least-squares fits of fluorescence intensity data for cycles near the onset of the growth phase. The resulting y0 values are statistically equivalent to the corresponding Cq if and only if E is taken to be error free. But uncertainty in E usually dominates the total uncertainty in y0, making the latter much degraded in precision compared with Cq. Bias in E can be an even greater source of error in y0. So-called mechanistic models achieve higher precision in estimating y0 by tacitly assuming E=2 in the baseline region and so are subject to this bias error. When used in calibration, the mechanistic y0 is statistically comparable to Cq from the other methods. When a signal threshold yq is used to define Cq, best estimation precision is obtained by setting yq near the maximum signal in the range of fitted cycles, in conflict with common practice in the y0 estimation algorithms.

Estimation of the relative sensitivity of qPCR analysis using pooled samples

The high sensitivity of qPCR makes it a desirable diagnostic method in epidemiological surveillance programs. However, due to high costs, the use of pooling has been suggested. In this paper, an algorithm based on the Montecarlo method has been designed and implemented. The algorithm had been tested in many different situations, and finally it was validated with a real dataset. Moreover, based on the results obtained and depending on pooling conditions, a drastic decrease of sensitivity is observed.

Quantification of the pirimicarb resistance allele frequency in pooled cotton aphid (Aphis gossypii Glover) samples by TaqMan SNP genotyping assay

Background

Pesticide resistance monitoring is a crucial part to achieving sustainable integrated pest management (IPM) in agricultural production systems. Monitoring of resistance in arthropod populations is initially performed by bioassay, a method that detects a phenotypic response to pesticides. Molecular diagnostic assays, offering speed and cost improvements, can be developed when the causative mutation for resistance has been identified. However, improvements to throughput are limited as genotyping methods cannot be accurately applied to pooled DNA. Quantifying an allele frequency from pooled DNA would allow faster and cheaper monitoring of pesticide resistance.

Methodology/Principal Findings

We demonstrate a new method to quantify a resistance allele frequency (RAF) from pooled insects via TaqMan assay by using raw fluorescence data to calculate the transformed fluorescence ratio k’ at the inflexion point based on a four parameter sigmoid curve. Our results show that k’ is reproducible and highly correlated with RAF (r >0.99). We also demonstrate that k’ has a non-linear relationship with RAF and that five standard points are sufficient to build a prediction model. Additionally, we identified a non-linear relationship between runs for k’, allowing the combination of samples across multiple runs in a single analysis.

Conclusions/Significance

The transformed fluorescence ratio (k') method can be used to monitor pesticide resistance in IPM and to accurately quantify allele frequency from pooled samples. We have determined that five standards (0.0, 0.2, 0.5, 0.8, and 1.0) are sufficient for accurate prediction and are statistically-equivalent to the 13 standard points used experimentally

Evidence of immune and inflammatory processes in the gills of AGD-affected Atlantic salmon, Salmo salar L

Amoebic gill disease (AGD) is a disease caused by the ectoparasite Neoparamoeba perurans which affects several cultured marine fish worldwide. The characterisation of pro-inflammatory and immune related genes at the mRNA level in AGD-affected Atlantic salmon gills was performed at 10 days post-inoculation using 2D quantitative RT-PCR, a method of mapping transcriptional responses in tissues. The genes of interest were IL-1β, TNF-α, TCR-α chain, CD8, CD4, MHC-IIα, MHC-I, IgM and IgT. A significant increase in expression of the mRNA of all the genes was observed in the gills of AGD-affected fish. Contrary to previous studies, our data suggest that the parasite, N. perurans, elicits a classical inflammatory response in the gills of AGD-affected fish and indicates that the mRNA expression of immune genes within gill lesions misrepresents the cellular immune response in the gills during AGD.

Comparing real-time quantitative polymerase chain reaction analysis methods for precision, linearity, and accuracy of estimating amplification efficiency

New methods are used to compare seven qPCR analysis methods for their performance in estimating the quantification cycle (Cq) and amplification efficiency (E) for a large test data set (94 samples for each of 4 dilutions) from a recent study. Precision and linearity are assessed using chi-square (χ(2)), which is the minimized quantity in least-squares (LS) fitting, equivalent to the variance in unweighted LS, and commonly used to define statistical efficiency. All methods yield Cqs that vary strongly in precision with the starting concentration N0, requiring weighted LS for proper calibration fitting of Cq vs log(N0). Then χ(2) for cubic calibration fits compares the inherent precision of the Cqs, while increases in χ(2) for quadratic and linear fits show the significance of nonlinearity. Nonlinearity is further manifested in unphysical estimates of E from the same Cq data, results which also challenge a tenet of all qPCR analysis methods - that E is constant throughout the baseline region. Constant-threshold (Ct) methods underperform the other methods when the data vary considerably in scale, as these data do.

New concepts in diagnostics for infectious diarrhea

Conventional approaches to the diagnosis of infectious diarrhea must include several modalities to detect an array of potential viruses, bacteria, and parasites. We will provide a general overview of the wide range of diagnostic modalities available for enteropathogens, briefly discuss some of the limitations of conventional methods, and then focus on new molecular methods, including real-time PCR and next-generation sequencing. In particular, we will discuss quantitation of pathogen load with these techniques. We will then describe examples whereby novel diagnostics may help illuminate the etiology of infectious diarrhea, where they may not, and how they may benefit studies of immunity to enteric infections.