Efficient computation of spaced seed hashing with block indexing

Background

Spaced-seeds, i.e. patterns in which some fixed positions are allowed to be wild-cards, play a crucial role in several bioinformatics applications involving substrings counting and indexing, by often providing better sensitivity with respect to k-mers based approaches. K-mers based approaches are usually fast, being based on efficient hashing and indexing that exploits the large overlap between consecutive k-mers. Spaced-seeds hashing is not as straightforward, and it is usually computed from scratch for each position in the input sequence. Recently, the FSH (Fast Spaced seed Hashing) approach was proposed to improve the time required for computation of the spaced seed hashing of DNA sequences with a speed-up of about 1.5 with respect to standard hashing computation.

Results

In this work we propose a novel algorithm, Fast Indexing for Spaced seed Hashing (FISH), based on the indexing of small blocks that can be combined to obtain the hashing of spaced-seeds of any length. The method exploits the fast computation of the hashing of runs of consecutive 1 in the spaced seeds, that basically correspond to k-mer of the length of the run.

Conclusions

We run several experiments, on NGS data from simulated and synthetic metagenomic experiments, to assess the time required for the computation of the hashing for each position in each read with respect to several spaced seeds. In our experiments, FISH can compute the hashing values of spaced seeds with a speedup, with respect to the traditional approach, between 1.9x to 6.03x, depending on the structure of the spaced seeds.

     

<Emphasis Type="Italic">massPix</Emphasis>: an R package for annotation and interpretation of mass spectrometry imaging data for lipidomics

Introduction

Mass spectrometry imaging (MSI) experiments result in complex multi-dimensional datasets, which require specialist data analysis tools.

Objectives

We have developed massPix—an R package for analysing and interpreting data from MSI of lipids in tissue.

Methods

massPix produces single ion images, performs multivariate statistics and provides putative lipid annotations based on accurate mass matching against generated lipid libraries.

Results

Classification of tissue regions with high spectral similarly can be carried out by principal components analysis (PCA) or k-means clustering.

Conclusion

massPix is an open-source tool for the analysis and statistical interpretation of MSI data, and is particularly useful for lipidomics applications.

     

Optimal choice of word length when comparing two Markov sequences using a <Emphasis Type="Italic">χ</Emphasis><Superscript>2</Superscript>-statistic

Background

Alignment-free sequence comparison using counts of word patterns (grams, k-tuples) has become an active research topic due to the large amount of sequence data from the new sequencing technologies. Genome sequences are frequently modelled by Markov chains and the likelihood ratio test or the corresponding approximate χ ²-statistic has been suggested to compare two sequences. However, it is not known how to best choose the word length k in such studies.

Results

We develop an optimal strategy to choose k by maximizing the statistical power of detecting differences between two sequences. Let the orders of the Markov chains for the two sequences be r ₁ and r ₂, respectively. We show through both simulations and theoretical studies that the optimal k= max(r ₁,r ₂)+1 for both long sequences and next generation sequencing (NGS) read data. The orders of the Markov chains may be unknown and several methods have been developed to estimate the orders of Markov chains based on both long sequences and NGS reads. We study the power loss of the statistics when the estimated orders are used. It is shown that the power loss is minimal for some of the estimators of the orders of Markov chains.

Conclusion

Our studies provide guidelines on choosing the optimal word length for the comparison of Markov sequences.

     

Iridoid and phenylethanoid/phenylpropanoid metabolite profiles of <Emphasis Type="Italic">Scrophularia</Emphasis> and <Emphasis Type="Italic">Verbascum</Emphasis> species used medicinally in North America

Introduction

Botanicals containing iridoid and phenylethanoid/phenylpropanoid glycosides are used worldwide for the treatment of inflammatory musculoskeletal conditions that are primary causes of human years lived with disability, such as arthritis and lower back pain.

Objectives

We report the analysis of candidate anti-inflammatory metabolites of several endemic Scrophularia species and Verbascum thapsus used medicinally by peoples of North America.

Methods

Leaves, stems, and roots were analyzed by ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) and partial least squares-discriminant analysis (PLS-DA) was performed in MetaboAnalyst 3.0 after processing the datasets in Progenesis QI.

Results

Comparison of the datasets revealed significant and differential accumulation of iridoid and phenylethanoid/phenylpropanoid glycosides in the tissues of the endemic Scrophularia species and Verbascum thapsus.

Conclusions

Our investigation identified several species of pharmacological interest as good sources for harpagoside and other important anti-inflammatory metabolites.

     

Metabolic response of porcine colon explants to in vitro infection by <Emphasis Type="Italic">Brachyspira hyodysenteriae</Emphasis>: a leap into disease pathophysiology

Introduction

Swine dysentery caused by Brachyspira hyodysenteriae is a production limiting disease in pig farming. Currently antimicrobial therapy is the only treatment and control method available.

Objective

The aim of this study was to characterize the metabolic response of porcine colon explants to infection by B. hyodysenteriae.

Methods

Porcine colon explants exposed to B. hyodysenteriae were analyzed for histopathological, metabolic and pro-inflammatory gene expression changes.

Results

Significant epithelial necrosis, increased levels of l-citrulline and IL-1α were observed on explants infected with B. hyodysenteriae.

Conclusions

The spirochete induces necrosis in vitro likely through an inflammatory process mediated by IL-1α and NO.

     

A comparative study of <Emphasis Type="Italic">k</Emphasis>-spectrum-based error correction methods for next-generation sequencing data analysis

Background

Innumerable opportunities for new genomic research have been stimulated by advancement in high-throughput next-generation sequencing (NGS). However, the pitfall of NGS data abundance is the complication of distinction between true biological variants and sequence error alterations during downstream analysis. Many error correction methods have been developed to correct erroneous NGS reads before further analysis, but independent evaluation of the impact of such dataset features as read length, genome size, and coverage depth on their performance is lacking. This comparative study aims to investigate the strength and weakness as well as limitations of some newest k-spectrum-based methods and to provide recommendations for users in selecting suitable methods with respect to specific NGS datasets.

Methods

Six k-spectrum-based methods, i.e., Reptile, Musket, Bless, Bloocoo, Lighter, and Trowel, were compared using six simulated sets of paired-end Illumina sequencing data. These NGS datasets varied in coverage depth (10× to 120×), read length (36 to 100 bp), and genome size (4.6 to 143 MB). Error Correction Evaluation Toolkit (ECET) was employed to derive a suite of metrics (i.e., true positives, false positive, false negative, recall, precision, gain, and F-score) for assessing the correction quality of each method.

Results

Results from computational experiments indicate that Musket had the best overall performance across the spectra of examined variants reflected in the six datasets. The lowest accuracy of Musket (F-score?=?0.81) occurred to a dataset with a medium read length (56 bp), a medium coverage (50×), and a small-sized genome (5.4 MB). The other five methods underperformed (F-score?<?0.80) and/or failed to process one or more datasets.

Conclusions

This study demonstrates that various factors such as coverage depth, read length, and genome size may influence performance of individual k-spectrum-based error correction methods. Thus, efforts have to be paid in choosing appropriate methods for error correction of specific NGS datasets. Based on our comparative study, we recommend Musket as the top choice because of its consistently superior performance across all six testing datasets. Further extensive studies are warranted to assess these methods using experimental datasets generated by NGS platforms (e.g., 454, SOLiD, and Ion Torrent) under more diversified parameter settings (k-mer values and edit distances) and to compare them against other non-k-spectrum-based classes of error correction methods.

     

Intra-specific variation of sperm length in the malaria vector <Emphasis Type="Italic">Anopheles gambiae</Emphasis>: males with shorter sperm have higher reproductive success

Background

Intra-specific variation in sperm length influences male reproductive success in several species of insects. In males of the malaria vector Anopheles gambiae, sperm length is highly variable but the significance of this variation is unknown. Understanding what determines the reproductive success of male mosquitoes is critical for controlling malaria, and in particular for replacing natural populations with transgenic, malaria-resistant mosquitoes.

Methods

A laboratory population of A. gambiae males was tested for intra-specific variation in sperm length. A full-sib quantitative genetic design was used to test for a genetic component of sperm length in A. gambiae males and estimate its heritability. This study also tested for a relationship between sperm length and male reproductive success in A. gambiae. Male reproductive success was measured as the proportions of inseminated and ovipositing females.

Results

There was intra-specific variation of sperm length in A. gambiae. There was no significant genetic variation in sperm length and its heritability was low (h² = 0.18) compared to other insects. Sperm length was correlated with male body size (measured as wing length). Males with short sperm had significantly higher reproductive success than males with long sperm and this was independent of body size.

Conclusion

This is the first study to demonstrate intra-specific variation in sperm length in A. gambiae and that males with short sperm have higher reproductive success. That sperm length influences female oviposition is important for any strategy considering the release of transgenic males.

     

<Emphasis Type="Italic">Rhizobium</Emphasis> strains in the biological control of the phytopathogenic fungi <Emphasis Type="Italic">Sclerotium</Emphasis> (<Emphasis Type="Italic">Athelia</Emphasis>) <Emphasis Type="Italic">rolfsii</Emphasis> on the common bean

Aims

To identify Rhizobium strains’ ability to biocontrol Sclerotium rolfsii, a fungus that causes serious damage to the common bean and other important crops, 78 previously isolated rhizobia from common bean were assessed.

Methods

Dual cultures, volatiles, indole-acetic acid (IAA), siderophore production and 16S rRNA sequencing were employed to select strains for pot and field experiments.

Results

Thirty-three antagonistic strains were detected in dual cultures, 16 of which were able to inhibit ≥84% fungus mycelial growth. Antagonistic strains produced up to 36.5 μg mL^?1 of IAA, and a direct correlation was verified between IAA production and mycelium inhibition. SEMIA 460 inhibited 45% of mycelial growth through volatile compounds. 16S rRNA sequences confirmed strains as Rhizobium species. In pot condition, common bean plants grown on S. rolfsii-infested soil and inoculated with SEMIA 4032, 4077, 4088, 4080, 4085, or 439 presented less or no disease symptoms. The most efficient strains under field conditions, SEMIA 439 and 4088, decreased disease incidence by 18.3 and 14.5% of the S. rolfsii-infested control.

Conclusions

Rhizobium strains could be strong antagonists towards S. rolfsii growth. SEMIA 4032, 4077, 4088, 4080, 4085, and 439 are effective in the biological control of the collar rot of the common bean.

     

Screening,tandem expression and immune activity appraisal of <Emphasis Type="Italic">Deinagkistrodon acutus</Emphasis> (pit viper) venom mimotopes from a phage display 12-mer peptide library

Objective

To design a specific polyclonal antibody against Deinagkistrodon acutus venom (DA-pAb) by immunizating New Zealand white rabbits.

Results

The IgG fraction was purified by affinity chromatography, and specific antibodies were purified by immunoaffinity chromatography. Polyclonal antibodies were subjected to ELISA and western blotting to evaluate their immune reactivity. We identified the mimotopes by screening a phage display 12-mer peptide library against D. acutus venom. After three rounds of biopanning with DA-pAb, 30 positive clones were identified. Eighteen phage clones were sequenced, and their corresponding amino acid sequences were deduced. Additional immunoassays with the peptides and DA-pAb identified five sequences as possible epitopes. Recombinant antigens synthesized with the five epitopes were used for the immunization of BALB/c mice.

Conclusion

The antibodies induced by these peptides recognized the recombinant antigen and D. acutus venom and protected mice against the hemorrhagic effects of the venom.

     

NS-kNN: a modified <Emphasis Type="Italic">k</Emphasis>-nearest neighbors approach for imputing metabolomics data

Introduction

A common problem in metabolomics data analysis is the existence of a substantial number of missing values, which can complicate, bias, or even prevent certain downstream analyses. One of the most widely-used solutions to this problem is imputation of missing values using a k-nearest neighbors (kNN) algorithm to estimate missing metabolite abundances. kNN implicitly assumes that missing values are uniformly distributed at random in the dataset, but this is typically not true in metabolomics, where many values are missing because they are below the limit of detection of the analytical instrumentation.

Objectives

Here, we explore the impact of nonuniformly distributed missing values (missing not at random, or MNAR) on imputation performance. We present a new model for generating synthetic missing data and a new algorithm, No-Skip kNN (NS-kNN), that accounts for MNAR values to provide more accurate imputations.

Methods

We compare the imputation errors of the original kNN algorithm using two distance metrics, NS-kNN, and a recently developed algorithm KNN-TN, when applied to multiple experimental datasets with different types and levels of missing data.

Results

Our results show that NS-kNN typically outperforms kNN when at least 20–30% of missing values in a dataset are MNAR. NS-kNN also has lower imputation errors than KNN-TN on realistic datasets when at least 50% of missing values are MNAR.

Conclusion

Accounting for the nonuniform distribution of missing values in metabolomics data can significantly improve the results of imputation algorithms. The NS-kNN method imputes missing metabolomics data more accurately than existing kNN-based approaches when used on realistic datasets.

     

Root morphology and its contribution to a large root system for phosphorus uptake by <Emphasis Type="Italic">Rytidosperma</Emphasis> species (wallaby grass)

Aims

Rytidosperma species are native Australian grasses which have different growth rates and phosphorus (P) requirements. This study examined the role of root morphology traits in response to P supply.

Methods

Nine Rytidosperma species ranging from slow- to fast-growth were examined along with Lolium perenne and Bromus hordeaceus. Plants were grown in a glasshouse for 47 days in soil supplied with six levels of P between 0 and 60 mg P per pot. Root mass, length and diameter, root hair length and density, and extent of mycorrhizal colonisation were measured.

Results

Across all species there was a positive correlation (P < 0.001) between P uptake and root mass, length and root hair cylinder volume (RHCV; estimated using root diameter, root hair length and root length) at all levels of P supply. An exception was the RHCV of B. hordeaceus, where expected P uptake was not achieved due to a markedly reduced root length at low-P supply. For the Rytidosperma species, morphological plasticity for specific root length, root mass fraction and root hair length ranged from 1.5-fold to 2.7-fold between high- and low-P supply. However, across all species and P levels no single root morphological trait was identified for universally increasing the size of the root system and P uptake.

Conclusions

Fast-growing species took up more P as a result of an overall larger root mass, greater root length and larger RHCV.

     

Adjunctive corticosteroid therapy for inpatients with <Emphasis Type="Italic">Mycoplasma pneumoniae</Emphasis> pneumonia

Background

There is conflicting evidence regarding the benefit of adjunctive corticosteroid therapy in patients with Mycoplasma pneumoniae pneumonia. We hypothesised that corticosteroid therapy could reduce mortality and length of stay (LOS) in such patients.

Methods

Adult patients with M. pneumoniae pneumonia from January 2010 to December 2013 were identified from the Japanese Diagnosis Procedure Combination inpatient database. The effects of low-dose and high-dose corticosteroid therapies on mortality, LOS, drug costs and hyperglycaemia requiring insulin treatment were evaluated using propensity score analyses.

Results

Eligible patients (n?=?2228) from 630 hospitals were divided into no-corticosteroid (n?=?1829), low-dose corticosteroid (n?=?267) and high-dose corticosteroid (n?=?132) groups. The propensity score-matched pairs were generated from no-corticoid and low-dose corticoid groups (251 pairs), or no-corticoid and high-dose corticosteroid groups (120 pairs). Adjunctive corticosteroid therapy did not decrease 30-day mortality. In addition, both low-dose and high-dose corticosteroid therapies were associated with increases in LOS. Furthermore, hyperglycaemia requiring insulin treatment and drug cost increased with corticosteroid use.

Conclusions

Adjunctive treatment with low-dose or high-dose corticosteroids may not be beneficial in M. pneumoniae pneumonia.

     

SNP markers associated with body size and pelt length in American mink (<Emphasis Type="Italic">Neovison vison</Emphasis>)

Background

Identification of genes underlying production traits is a key aim of the mink research community. Recent availability of genomic tools have opened the possibility for faster genetic progress in mink breeding. Availability of mink genome assembly allows genome-wide association studies in mink.

Results

In this study, we used genotyping-by-sequencing to obtain single nucleotide polymorphism (SNP) genotypes of 2496 mink. After multiple rounds of filtering, we retained 28,336 high quality SNPs and 2352 individuals for a genome-wide association study (GWAS). We performed the first GWAS for body weight, behavior, along with 10 traits related to fur quality in mink.

Conclusions

Combining association results with existing functional information of genes and mammalian phenotype databases, we proposed WWC3, MAP2K4, SLC7A1 and USP22 as candidate genes for body weight and pelt length in mink.

     

A constitutive expression system for <Emphasis Type="Italic">Pichia pastoris</Emphasis> based on the <Emphasis Type="Italic">PGK1</Emphasis> promoter

Objectives

To develop a new vector for constitutive expression in Pichia pastoris based on the endogenous glycolytic PGK1 promoter.

Results

P. pastoris plasmids bearing at least 415 bp of PGK1 promoter sequences can be used to drive plasmid integration by addition at this locus without affecting cell growth. Based on this result, a new P. pastoris integrative vector, pPICK2, was constructed bearing some features that facilitate protein production in this yeast: a ~620 bp PGK1 promoter fragment with three options of restriction sites for plasmid linearization prior to yeast transformation: a codon-optimized α-factor secretion signal, a new polylinker, and the kan marker for vector propagation in bacteria and selection of yeast transformants.

Conclusions

A new constitutive vector for P. pastoris represents an alternative platform for recombinant protein production and metabolic engineering purposes.

     

Extraction of natural products from bark of <Emphasis Type="Italic">Betula pendula</Emphasis> using ionic liquids

Introduction

Aqueous–methanol mixtures have successfully been applied to extract a broad range of metabolites from plant tissue. However, a certain amount of material remains insoluble.

Objectives

To enlarge the metabolic compendium, two ionic liquids were selected to extract the methanol insoluble part of trunk from Betula pendula.

Methods

The extracted compounds were analyzed by LC/MS and GC/MS.

Results

The results show that 1-butyl-3-methylimidazolium acetate (IL-Ac) predominantly resulted in fatty acids, whereas 1-ethyl-3-methylimidazolium tosylate (IL-Tos) mostly yielded phenolic structures. Interestingly, bark yielded more ionic liquid soluble metabolites compared to interior wood.

Conclusion

From this one can conclude that the application of ionic liquids may expand the metabolic snapshot.

     

Characterization of an aryl-alcohol oxidase from the plant saprophytic basidiomycete <Emphasis Type="Italic">Coprinopsis cinerea</Emphasis> with broad substrate specificity against aromatic alcohols

Objectives

The aim of the study was to obtain information about the enzymatic properties of aryl-alcohol oxidase from the plant saprophytic basidiomycete Coprinopsis cinerea (rCcAAO), which is classified into the auxiliary activities family 3 subfamily 2 (AA3_2).

Results

The gene encoding AAO from the plant saprophytic basidiomycete Coprinopsis cinerea (CcAAO) was cloned, and the recombinant CcAAO (rCcAAO) was heterologously expressed in the methylotrophic yeast Pichia pastoris. The purified rCcAAO showed significant activity not only against trans,trans-2,4-hexadien-1-ol but also against a broad range of aromatic alcohols including aromatic compounds that were reported to be poor substrates for known AAOs. Moreover, site-directed mutagenesis analysis demonstrated that mutants with substitutions from leucine to phenylalanine and tryptophan at position 416 exhibited decreases of activity for aromatic alcohols but still maintained the activity for trans,trans-2,4-hexadien-1-ol.

Conclusions

Leucine 416 in CcAAO contributes to the broad substrate specificity against various aromatic alcohols, which is useful for the production of hydrogen peroxide using this enzyme.

     

Study on the Clinical Features and Prognosis of <Emphasis Type="Italic">Penicilliosis marneffei</Emphasis> Without Human Immunodeficiency Virus Infection

Objective

To improve the diagnosis and treatment of Penicilliosis marneffei without human immunodeficiency virus infection.

Methods

Analyze and review the clinical features, diagnosis and treatment of six cases of P. marneffei without human immunodeficiency virus infection at The First Affiliated Hospital of Fujian Medical University.

Results

Two cases were diagnosed in the ENT Department, three cases in the respiratory department and one case in the dermatological department. Penicillium marneffei infection was confirmed by sputum culture, blood culture and tissue biopsy. After definite diagnosis, one refused further treatment, and others showed significant improvement.

Conclusion

Penicilliosis marneffei is insidious onset and easy to be escaped and misdiagnosed. To achieve early diagnosis and appropriate treatment, doubtful cases should be alerted for the diagnoses as P. marneffei.

     

<Emphasis Type="Italic">Sedum</Emphasis> root foraging in layered green roof substrates

Background and aims

Layered profiles of designed soils may provide long-term benefits for green roofs, provided the vegetation can exploit resources in the different layers. We aimed to quantify Sedum root foraging for water and nutrients in designed soils of different texture and layering.

Methods

In a controlled pot experiment we quantified the root foraging ability of the species Sedum album (L.) and S. rupestre (L.) in response to substrate structure (fine, coarse, layered or mixed), vertical fertiliser placement (top or bottom half of pot) and watering (5, 10 or 20 mm week^?1).

Results

Water availability was the main driver of plant growth, followed by substrate structure, while fertiliser placement only had marginal effects on plant growth. Root foraging ability was low to moderate, as also reflected in the low proportion of biomass allocated to roots (5–13%). Increased watering reduced the proportion of root length and root biomass in deeper layers.

Conclusions

Both S. album and S. rupestre had a low ability to exploit water and nutrients by precise root foraging in substrates of different texture and layering. Allocation of biomass to roots was low and showed limited flexibility even under water-deficient conditions.