Validated predictive modelling of the environmental resistome

Multi-drug-resistant bacteria pose a significant threat to public health. The role of the environment in the overall rise in antibiotic-resistant infections and risk to humans is largely unknown. This study aimed to evaluate drivers of antibiotic-resistance levels across the River Thames catchment, model key biotic, spatial and chemical variables and produce predictive models for future risk assessment. Sediment samples from 13 sites across the River Thames basin were taken at four time points across 2011 and 2012. Samples were analysed for class 1 integron prevalence and enumeration of third-generation cephalosporin-resistant bacteria. Class 1 integron prevalence was validated as a molecular marker of antibiotic resistance; levels of resistance showed significant geospatial and temporal variation. The main explanatory variables of resistance levels at each sample site were the number, proximity, size and type of surrounding wastewater-treatment plants. Model 1 revealed treatment plants accounted for 49.5% of the variance in resistance levels. Other contributing factors were extent of different surrounding land cover types (for example, Neutral Grassland), temporal patterns and prior rainfall; when modelling all variables the resulting model (Model 2) could explain 82.9% of variations in resistance levels in the whole catchment. Chemical analyses correlated with key indicators of treatment plant effluent and a model (Model 3) was generated based on water quality parameters (contaminant and macro- and micro-nutrient levels). Model 2 was beta tested on independent sites and explained over 78% of the variation in integron prevalence showing a significant predictive ability. We believe all models in this study are highly useful tools for informing and prioritising mitigation strategies to reduce the environmental resistome.     

Hierarchical classification of functionally equivalent genes in prokaryotes

Functional classification of genes represents a fundamental problem to many biological studies. Most of the existing classification schemes are based on the concepts of homology and orthology, which were originally introduced to study gene evolution but might not be the most appropriate for gene function prediction, particularly at high resolution level. We have recently developed a scheme for hierarchical classification of genes (HCGs) in prokaryotes. In the HCG scheme, the functional equivalence relationships among genes are first assessed through a careful application of both sequence similarity and genomic neighborhood information; and genes are then classified into a hierarchical structure of clusters, where genes in each cluster are functionally equivalent at some resolution level, and the level of resolution goes higher as the clusters become increasingly smaller traveling down the hierarchy. The HCG scheme is validated through comparisons with the taxonomy of the prokaryotic genomes, Clusters of Orthologous Groups (COGs) of genes and the Pfam system. We have applied the HCG scheme to 224 complete prokaryotic genomes, and constructed a HCG database consisting of a forest of 5339 multi-level and 15 770 single-level trees of gene clusters covering approximately 93% of the genes of these 224 genomes. The validation results indicate that the HCG scheme not only captures the key features of the existing classification schemes but also provides a much richer organization of genes which can be used for functional prediction of genes at higher resolution and to help reveal evolutionary trace of the genes.     

Relationship of aerobic and anaerobic parameters with 400 m front crawl swimming performance

The aims of the present study were to investigate the relationship of aerobic and anaerobic parameters with 400 m performance, and establish which variable better explains long distance performance in swimming. Twenty-two swimmers (19.1±1.5 years, height 173.9±10.0 cm, body mass 71.2±10.2 kg; 76.6±5.3% of 400 m world record) underwent a lactate minimum test to determine lactate minimum speed (LMS) (i.e., aerobic capacity index). Moreover, the swimmers performed a 400 m maximal effort to determine mean speed (S400m), peak oxygen uptake (V.O2PEAK) and total anaerobic contribution (C_ANA). The C_ANA was assumed as the sum of alactic and lactic contributions. Physiological parameters of 400 m were determined using the backward extrapolation technique (V.O2PEAK and alactic contributions of C_ANA) and blood lactate concentration analysis (lactic anaerobic contributions of C_ANA). The Pearson correlation test and backward multiple regression analysis were used to verify the possible correlations between the physiological indices (predictor factors) and S400m (independent variable) (p < 0.05). Values are presented as mean ± standard deviation. Significant correlations were observed between S400m (1.4±0.1 m·s^-1) and LMS (1.3±0.1 m·s^-1; r = 0.80), V.O2PEAK (4.5±3.9 L·min^-1; r = 0.72) and C_ANA (4.7±1.5 L·O₂; r= 0.44). The best model constructed using multiple regression analysis demonstrated that LMS and V.O2PEAK explained 85% of the 400 m performance variance. When backward multiple regression analysis was performed, C_ANA lost significance. Thus, the results demonstrated that both aerobic parameters (capacity and power) can be used to predict 400 m swimming performance.     

Improving the performance of protein threading using insertion/deletion frequency arrays

As a protein evolves, not every part of the amino acid sequence has an equal probability of being deleted or for allowing insertions, because not every amino acid plays an equally important role in maintaining the protein structure. However, the most prevalent models in fold recognition methods treat every amino acid deletion and insertion as equally probable events. We have analyzed the alignment patterns for homologous and analogous sequences to determine patterns of insertion and deletion, and used that information to determine the statistics of insertions and deletions for different amino acids of a target sequence. We define these patterns as insertion/deletion (indel) frequency arrays (IFAs). By applying IFAs to the protein threading problem, we have been able to improve the alignment accuracy, especially for proteins with low sequence identity. We have also demonstrated that the application of this information can lead to an improvement in fold recognition.