GuiTope: an application for mapping random-sequence peptides to protein sequences

Analysis-ready VCF at Biobank scale using Zarr

Background

Variant Call Format (VCF) is the standard file format for interchanging genetic variation data and associated quality control metrics. The usual row-wise encoding of the VCF data model (either as text or packed binary) emphasises efficient retrieval of all data for a given variant, but accessing data on a field or sample basis is inefficient. Biobank scale datasets currently available consist of hundreds of thousands of whole genomes and hundreds of terabytes of compressed VCF. Row-wise data storage is fundamentally unsuitable and a more scalable approach is needed.

Results

We present the VCF Zarr specification, an encoding of the VCF data model using Zarr which makes retrieving subsets of the data much more efficient. Zarr is a cloud-native format for storing multi-dimensional data, widely used in scientific computing. We show how this format is far more efficient than standard VCF based approaches, and competitive with specialised methods for storing genotype data in terms of compression ratios and calculation performance. We demonstrate the VCF Zarr format (and the vcf2zarr conversion utility) on a subset of the Genomics England aggV2 dataset comprising 78,195 samples and 59,880,903 variants, with a 5X reduction in storage and greater than 300X reduction in CPU usage in some representative benchmarks.

Conclusions

Large row-encoded VCF files are a major bottleneck for current research, and storing and processing these files incurs a substantial cost. The VCF Zarr specification, building on widely-used, open-source technologies has the potential to greatly reduce these costs, and may enable a diverse ecosystem of next-generation tools for analysing genetic variation data directly from cloud-based object stores.

Data Independent Acquisition to Inform the Development of Targeted Proteomics Assays Using a Triple Quadrupole Mass Spectrometer

Mass spectrometry based targeted proteomics methods provide sensitive and high-throughput analysis of selected proteins. To develop a targeted bottom-up proteomics assay, peptides must be evaluated as proxies for the measurement of a protein or proteoform in a biological matrix. Candidate peptide selection typically relies on predetermined biochemical properties, data from semi-stochastic sampling, or by empirical measurements. These strategies require extensive testing and method refinement due to the difficulties associated with prediction of peptide response in the biological matrix of interest. Gas-phase fractionated (GPF) narrow window data-independent acquisition (DIA) aids in the development of reproducible selected reaction monitoring (SRM) assays by providing matrix-specific information on peptide detectability and quantification by mass spectrometry. To demonstrate the suitability of DIA data for selecting peptide targets, we reimplement a portion of an existing assay to measure 98 Alzheimer's disease proteins in cerebrospinal fluid (CSF). Peptides were selected from GPF-DIA based on signal intensity and reproducibility. The resulting SRM assay exhibits similar quantitative precision to published data, despite the inclusion of different peptides between the assays. This workflow enables development of new assays without additional up-front data acquisition, demonstrated here through generation of a separate assay for an unrelated set of proteins in CSF from the same dataset.

Genetic analysis and QTL mapping for pericarp thickness in maize (Zea mays L.)

Proper pericarp thickness protects the maize kernel against pests and diseases, moreover, thinner pericarp improves the eating quality in fresh corn. In this study, we aimed to investigate the dynamic changes in maize pericarp during kernel development and identified the major quantitative trait loci (QTLs) for maize pericarp thickness. It was observed that maize pericarp thickness first increased and then decreased. During the growth and formation stages, the pericarp thickness gradually increased and reached the maximum, after which it gradually decreased and reached the minimum during maturity. To identify the QTLs for pericarp thickness, a BC4F4 population was constructed using maize inbred lines B73 (recurrent parent with thick pericarp) and Baimaya (donor parent with thin pericarp). In addition, a high-density genetic map was constructed using maize 10 K SNP microarray. A total of 17 QTLs related to pericarp thickness were identified in combination with the phenotypic data. The results revealed that the heritability of the thickness of upper germinal side of pericarp (UG) was 0.63. The major QTL controlling UG was qPT1-1, which was located on chromosome 1 (212,215,145-212,948,882). The heritability of the thickness of upper abgerminal side of pericarp (UA) was 0.70. The major QTL controlling UA was qPT2-1, which was located on chromosome 2 (2,550,197-14,732,993). In addition, a combination of functional annotation, DNA sequencing analysis and quantitative real-time PCR (qPCR) screened two candidate genes, Zm00001d001964 and Zm00001d002283, that could potentially control maize pericarp thickness. This study provides valuable insights into the improvement of maize pericarp thickness during breeding.

What makes clocks tick? Characterizing developmental dynamics of adult epigenetic clock sites

DNA methylation (DNAm) at specific sites can be used to calculate 'epigenetic clocks', which in adulthood are used as indicators of age(ing). However, little is known about how these clock sites 'behave' during development and what factors influence their variability in early life. This knowledge could be used to optimize healthy aging well before the onset of age-related conditions. Here, we leveraged results from two longitudinal population-based cohorts (N=5,019 samples from 2,348 individuals) to characterize trajectories of adult clock sites from birth to early adulthood. We find that clock sites (i) diverge widely in their developmental trajectories, often showing non-linear change over time; (ii) are substantially more likely than non-clock sites to vary between individuals already from birth, differences that are predictive of DNAm variation at later ages; and (iii) show enrichment for genetic and prenatal environmental exposures, supporting an early-origins perspective to epigenetic aging.

Ensemble methods of rank-based trees for single sample classification with gene expression profiles

Building Single Sample Predictors (SSPs) from gene expression profiles presents challenges, notably due to the lack of calibration across diverse gene expression measurement technologies. However, recent research indicates the viability of classifying phenotypes based on the order of expression of multiple genes. Existing SSP methods often rely on Top Scoring Pairs (TSP), which are platform-independent and easy to interpret through the concept of "relative expression reversals". Nevertheless, TSP methods face limitations in classifying complex patterns involving comparisons of more than two gene expressions. To overcome these constraints, we introduce a novel approach that extends TSP rules by constructing rank-based trees capable of encompassing extensive gene-gene comparisons. This method is bolstered by incorporating two ensemble strategies, boosting and random forest, to mitigate the risk of overfitting. Our implementation of ensemble rank-based trees employs boosting with LogitBoost cost and random forests, addressing both binary and multi-class classification problems. In a comparative analysis across 12 cancer gene expression datasets, our proposed methods demonstrate superior performance over both the k-TSP classifier and nearest template prediction methods. We have further refined our approach to facilitate variable selection and the generation of clear, precise decision rules from rank-based trees, enhancing interpretability. The cumulative evidence from our research underscores the significant potential of ensemble rank-based trees in advancing disease classification via gene expression data, offering a robust, interpretable, and scalable solution. Our software is available at https://CRAN.R-project.org/package=ranktreeEnsemble .

Molecular cloning and functional analysis of 4-coumarate: CoA ligases from Marchantia paleacea and their roles in lignin and flavanone biosynthesis

Phenylpropanoids play important roles in plant physiology and the enzyme 4-coumarate: coenzyme A ligase (4CL) catalyzes the formation of thioesters. Despite extensive characterization in various plants, the functions of 4CLs in the liverwort Marchantia paleacea remain unknown. Here, four 4CLs from M. paleacea were isolated and functionally analyzed. Heterologous expression in Escherichia coli indicated the presence of different enzymatic activities in the four enzymes. Mp4CL1 and Mp4CL2 were able to convert caffeic, p-coumaric, cinnamic, ferulic, dihydro-p-coumaric, and 5-hydroxyferulic acids to their corresponding CoA esters, while Mp4CL3 and Mp4CL4 catalyzed none. Mp4CL1 transcription was induced when M. paleacea thalli were treated with methyl jasmonate (MeJA). The overexpression of Mp4CL1 increased the levels of lignin in transgenic Arabidopsis. In addition, we reconstructed the flavanone biosynthetic pathway in E. coli. The pathway comprised Mp4CL1, co-expressed with chalcone synthase (CHS) from different plant species, and the efficiency of biosynthesis was optimal when both the 4CL and CHS were obtained from the same species M. paleacea.

The human Toll-like receptor 2 (TLR2) response during pathogenic Leptospira infection

Background

Human innate immune responses are triggered through the interaction of human pattern recognition receptors and pathogen-associated molecular patterns. The role of toll-like receptor2 (TLR2) in mice innate immune response to leptospirosis is well established, while human studies are limited. The present study aimed to determine the TLR2 response among confirmed cases of leptospirosis.

Methodology/Principle findings

The study has two components. Clinically suspected patients of leptospirosis were confirmed using a previously validated qPCR assay. Total RNA was extracted from patients' RNA-stabilized whole blood samples. Human TLR2 gene expression (RT-qPCR) analysis was carried out using an exon-exon spanning primer pair, using CFX Maestro™ software. The first set of patient samples was used to calculate the Relative Normalized Expression (ΔΔCq value) of the TLR2 gene in comparison to a healthy control sample and normalized by the reference gene GAPDH (Glyceraldehyde-3-phosphate dehydrogenase). Secondly, recruited patient samples were subjected to TLR2 gene expression analysis and compared to healthy controls and normalized by the reference genes Beta-2-microglobulin(B2M), Hypoxanthine phosphoribosyltransferase 1 (HPRT 1).In the initial cohort of 64 confirmed leptospirosis cases, 18 were selected for human TLR2 gene expression analysis based on criteria of leptospiremia and RNA yield. Within this group, one individual exhibited a down-regulation of TLR2 gene (Expression/ΔΔCq=0.01352), whereas the remaining subjects presented no significant change in gene expression. In a subsequent cohort of 23 confirmed cases, 13 were chosen for similar analysis. Among these, three patients demonstrated down-regulation of TLR2 gene expression, with Expression/ΔΔCq values of 0.86574, 0.47200, and 0.28579, respectively. No TLR2 gene expression was noted in the other patients within this second group.

Conclusions

Our investigation into the acute phase of leptospirosis using human clinical samples has revealed a downregulation of TLR2 gene expression. This observation contrasts to the upregulation commonly reported in the majority of in-vitro and in-vivo studies of Leptospira infection. These preliminary findings prompt a need for further research to explore the mechanisms underlying TLR2's role in the pathogenesis of leptospirosis, which may differ in clinical settings compared to laboratory models.

Author Summary

The human immune system employs pattern recognition receptors like toll-like receptor 2 (TLR2) to detect and combat infections such as leptospirosis. While TLR2's role is well-documented in mice, its function in the human response to leptospirosis remains unclear. Our study evaluated TLR2 activity in patients with confirmed leptospirosis. We conducted a genetic analysis of blood samples from these patients, comparing TLR2 gene activity against healthy individuals, with standard reference genes for accuracy. Contrary to expectations and existing laboratory data, we observed a decrease in TLR2 activity in some patients. This suggests that human TLR2 responses in actual infections may diverge from established laboratory models. These findings indicate a need for further study to understand the human immune response to leptospirosis, which may significantly differ from that observed in controlled experimental settings.

Fine-tuning TrailMap: The utility of transfer learning to improve the performance of deep learning in axon segmentation of light-sheet microscopy images

Light-sheet microscopy has made possible the 3D imaging of both fixed and live biological tissue, with samples as large as the entire mouse brain. However, segmentation and quantification of that data remains a time-consuming manual undertaking. Machine learning methods promise the possibility of automating this process. This study seeks to advance the performance of prior models through optimizing transfer learning. We fine-tuned the existing TrailMap model using expert-labeled data from noradrenergic axonal structures in the mouse brain. By fine-tuning the final two layers of the neural network at a lower learning rate of the TrailMap model, we demonstrate an improved recall and an occasionally improved adjusted F1-score within our test dataset over using the originally trained TrailMap model.

Origins and Evolution of Novel Bacteroides in Captive Apes

Bacterial strains evolve in response to the gut environment of their hosts, with genomic changes that influence their interactions with hosts as well as with other members of the gut community. Great apes in captivity have acquired strains of Bacteroides xylanisolvens, which are common within gut microbiome of humans but not typically found other apes, thereby enabling characterization of strain evolution following colonization. Here, we isolate, sequence and reconstruct the history of gene gain and loss events in numerous captive-ape-associated strains since their divergence from their closest human-associated strains. We show that multiple captive-ape-associated B. xylanisolvens lineages have independently acquired gene complexes that encode functions related to host mucin metabolism. Our results support the finding of high genome fluidity in Bacteroides, in that several strains, in moving from humans to captive apes, have rapidly gained large genomic regions that augment metabolic properties not previously present in their relatives.

Epigenome-wide association study of incident type 2 diabetes in Black and White participants from the Atherosclerosis Risk in Communities Study

DNA methylation studies of incident type 2 diabetes in US populations are limited, and to our knowledge none included individuals of African descent living in the US. We performed an epigenome-wide association analysis of blood-based methylation levels at CpG sites with incident type 2 diabetes using Cox regression in 2,091 Black and 1,029 White individuals from the Atherosclerosis Risk in Communities study. At an epigenome-wide significance threshold of 10-7, we detected 7 novel diabetes-associated CpG sites in C1orf151 (cg05380846: HR= 0.89, p = 8.4 × 10-12), ZNF2 (cg01585592: HR= 0.88, p = 1.6 × 10-9), JPH3 (cg16696007: HR= 0.87, p = 7.8 × 10-9), GPX6 (cg02793507: HR= 0.85, p = 2.7 × 10-8 and cg00647063: HR= 1.20, p = 2.5 × 10-8), chr17q25 (cg16865890: HR= 0.8, p = 6.9 × 10-8), and chr11p15 (cg13738793: HR= 1.11, p = 7.7 × 10-8). The CpG sites at C1orf151, ZNF2, JPH3 and GPX6, were identified in Black adults, chr17q25 was identified in White adults, and chr11p15 was identified upon meta-analyzing the two groups. The CpG sites at JPH3 and GPX6 were likely associated with incident type 2 diabetes independent of BMI. All the CpG sites, except at JPH3, were likely consequences of elevated glucose at baseline. We additionally replicated known type 2 diabetes-associated CpG sites including cg19693031 at TXNIP, cg00574958 at CPT1A, cg16567056 at PLBC2, cg11024682 at SREBF1, cg08857797 at VPS25, and cg06500161 at ABCG1, 3 of which were replicated in Black adults at the epigenome-wide threshold. We observed modest increase in type 2 diabetes variance explained upon addition of the significantly associated CpG sites to a Cox model that included traditional type 2 diabetes risk factors and fasting glucose (increase from 26.2% to 30.5% in Black adults; increase from 36.9% to 39.4% in White adults). We examined if groups of proximal CpG sites were associated with incident type 2 diabetes using a gene-region specific and a gene-region agnostic differentially methylated region (DMR) analysis. Our DMR analyses revealed several clusters of significant CpG sites, including a DMR consisting of a previously discovered CpG site at ADCY7 and promoter regions of TP63 which were differentially methylated across all race groups. This study illustrates improved discovery of CpG sites/regions by leveraging both individual CpG site and DMR analyses in an unexplored population. Our findings include genes linked to diabetes in experimental studies (e.g., GPX6, JPH3, and TP63), and future gene-specific methylation studies could elucidate the link between genes, environment, and methylation in the pathogenesis of type 2 diabetes.

scMD: cell type deconvolution using single-cell DNA methylation references

The proliferation of single-cell RNA sequencing data has led to the widespread use of cellular deconvolution, aiding the extraction of cell type-specific information from extensive bulk data. However, those advances have been mostly limited to transcriptomic data. With recent development in single-cell DNA methylation (scDNAm), new avenues have been opened for deconvolving bulk DNAm data, particularly for solid tissues like the brain that lack cell-type references. Due to technical limitations, current scDNAm sequences represent a small proportion of the whole genome for each single cell, and those detected regions differ across cells. This makes scDNAm data ultra-high dimensional and ultra-sparse. To deal with these challenges, we introduce scMD (single cell Methylation Deconvolution), a cellular deconvolution framework to reliably estimate cell type fractions from tissue-level DNAm data. To analyze large-scale complex scDNAm data, scMD employs a statistical approach to aggregate scDNAm data at the cell cluster level, identify cell-type marker DNAm sites, and create a precise cell-type signature matrix that surpasses state-of-the-art sorted-cell or RNA-derived references. Through thorough benchmarking in several datasets, we demonstrate scMD's superior performance in estimating cellular fractions from bulk DNAm data. With scMD-estimated cellular fractions, we identify cell type fractions and cell type-specific differentially methylated cytosines associated with Alzheimer's disease.

rCRUX: A Rapid and Versatile Tool for Generating Metabarcoding Reference libraries in R

Key to making accurate taxonomic assignments are curated, comprehensive reference barcode databases. However, the generation and curation of such databases has remained challenging given the large and continuously growing volumes of DNA sequence data and novel reference barcode targets. Monitoring and research applications require a greater diversity of specialized gene regions and targeted taxa to meet taxonomic classification goals then are currently curated by professional staff. Thus, there is a growing need for an easy to implement tool that can generate comprehensive metabarcoding reference libraries for any bespoke locus. We address this need by reimagining CRUX from the Anacapa Toolkit and present the rCRUX package in R. The typical workflow involves searching for plausible seed amplicons (get_seeds_local() or get_seeds_remote()) by simulating in silico PCR to acquire seed sequences containing a user-defined primer set. Next these seeds are used to iteratively blast search seed sequences against a local NCBI formatted database using a taxonomic rank based stratified random sampling approach (blast_seeds()) that results in a comprehensive set of sequence matches. This database is dereplicated and cleaned (derep_and_clean_db()) by identifying identical reference sequences and collapsing the taxonomic path to the lowest taxonomic agreement across all matching reads. This results in a curated, comprehensive database of primer specific reference barcode sequences from NCBI. We demonstrate that rCRUX provides more comprehensive reference databases for the MiFish Universal Teleost 12S, Taberlet trnl, and fungal ITS locus than CRABS, METACURATOR, RESCRIPt, and ECOPCR reference databases. We then further demonstrate the utility of rCRUX by generating 16 reference databases for metabarcoding loci that lack dedicated reference database curation efforts. The rCRUX package provides a simple to use tool for the generation of curated, comprehensive reference databases for user-defined loci, facilitating accurate and effective taxonomic classification of metabarcoding and DNA sequence efforts broadly.

PASSerRank: Prediction of Allosteric Sites with Learning to Rank

Allostery plays a crucial role in regulating protein activity, making it a highly sought-after target in drug development. One of the major challenges in allosteric drug research is the identification of allosteric sites. In recent years, many computational models have been developed for accurate allosteric site prediction. Most of these models focus on designing a general rule that can be applied to pockets of proteins from various families. In this study, we present a new approach using the concept of Learning to Rank (LTR). The LTR model ranks pockets based on their relevance to allosteric sites, i.e., how well a pocket meets the characteristics of known allosteric sites. The model outperforms other common machine learning models with higher F1 score and Matthews correlation coefficient. After the training and validation on two datasets, the Allosteric Database (ASD) and CASBench, the LTR model was able to rank an allosteric pocket in the top 3 positions for 83.6% and 80.5% of test proteins, respectively. The trained model is available on the PASSer platform (https://passer.smu.edu) to aid in drug discovery research.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

Prediction and Design of Protease Enzyme Specificity Using a Structure-Aware Graph Convolutional Network

Site-specific proteolysis by the enzymatic cleavage of small linear sequence motifs is a key post-translational modification involved in physiology and disease. The ability to robustly and rapidly predict protease substrate specificity would also enable targeted proteolytic cleavage - editing - of a target protein by designed proteases. Current methods for predicting protease specificity are limited to sequence pattern recognition in experimentally-derived cleavage data obtained for libraries of potential substrates and generated separately for each protease variant. We reasoned that a more semantically rich and robust model of protease specificity could be developed by incorporating the three-dimensional structure and energetics of molecular interactions between protease and substrates into machine learning workflows. We present Protein Graph Convolutional Network (PGCN), which develops a physically-grounded, structure-based molecular interaction graph representation that describes molecular topology and interaction energetics to predict enzyme specificity. We show that PGCN accurately predicts the specificity landscapes of several variants of two model proteases: the NS3/4 protease from the Hepatitis C virus (HCV) and the Tobacco Etch Virus (TEV) proteases. Node and edge ablation tests identified key graph elements for specificity prediction, some of which are consistent with known biochemical constraints for protease:substrate recognition. We used a pre-trained PGCN model to guide the design of TEV protease libraries for cleaving two non-canonical substrates, and found good agreement with experimental cleavage results. Importantly, the model can accurately assess designs featuring diversity at positions not present in the training data. The described methodology should enable the structure-based prediction of specificity landscapes of a wide variety of proteases and the construction of tailor-made protease editors for site-selectively and irreversibly modifying chosen target proteins.

Investigating the Complexity of Gene Co-expression Estimation for Single-cell Data

With the rapid advance of single-cell RNA sequencing (scRNA-seq) technology, understanding biological processes at a more refined single-cell level is becoming possible. Gene co-expression estimation is an essential step in this direction. It can annotate functionalities of unknown genes or construct the basis of gene regulatory network inference. This study thoroughly tests the existing gene co-expression estimation methods on simulation datasets with known ground truth co-expression networks. We generate these novel datasets using two simulation processes that use the parameters learned from the experimental data. We demonstrate that these simulations better capture the underlying properties of the real-world single-cell datasets than previously tested simulations for the task. Our performance results on tens of simulated and eight experimental datasets show that all methods produce estimations with a high false discovery rate potentially caused by high-sparsity levels in the data. Finally, we find that commonly used pre-processing approaches, such as normalization and imputation, do not improve the co-expression estimation. Overall, our benchmark setup contributes to the co-expression estimator development, and our study provides valuable insights for the community of single-cell data analyses.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

Newcastle disease burden in Nepal and efficacy of Tablet I2 vaccine in commercial and backyard poultry production

Poultry (Gallus domesticus) farming plays an important role as an income generating enterprise in a developing country like Nepal, contributing more than 4% to the national Gross Domestic Product (GDP). Newcastle Disease (ND) is a major poultry disease affecting both commercial and backyard poultry production worldwide. There were more than 90 reported ND outbreaks in Nepal in 2018 with over 74,986 birds being affected. ND is responsible for over 7% of total poultry mortality in the country. Recent outbreaks of ND in 2021 affected many farms throughout Nepal and caused massive loss in poultry production. ND is caused by a single-stranded ribonucleic acid (RNA) virus that presents very similar clinical symptoms as Influenza A (commonly known as bird flu) adding much complexity to clinical disease identification and intervention. We conducted a nationwide ND and Influenza A (IA) prevalence study, collecting samples from representative commercial and backyard poultry farms from across the major poultry production hubs of Nepal. We used both serological and molecular assessments to determine disease exposure history and identification of strains of ND Virus (NDV). Of the 40 commercial farms tested, both NDV (n = 28, 70%) and IAV (n = 11, 27.5%) antibodies were detected in majority of the samples. In the backyard farms (n = 36), sero-prevalence of NDV and IAV were 17.5% (n = 7) and 7.5% (n = 3) respectively. Genotype II NDV was present in most of the commercial farms, which was likely due to live vaccine usage. We detected never reported Genotype I NDV in two backyard farm samples. Our investigation into 2021 ND outbreak implicated Genotype VII.2 NDV strain as the causative pathogen. Additionally, we developed a Tablet formulation of the thermostable I2-NDV vaccine (Ranigoldunga™) and assessed its efficacy on various (mixed) breeds of chicken (Gallus domesticus). Ranigoldunga™ demonstrated an overall efficacy >85% with a stability of 30 days at room temperature (25°C). The intraocularly administered vaccine was highly effective in preventing ND, including Genotype VII.2 NDV strain.

Identification of long non-coding RNA using single nucleotide epimutation analysis: a novel gene discovery approach

BACKGROUND

Long non-coding RNAs (lncRNAs) are involved in a variety of mechanisms related to tumorigenesis by functioning as oncogenes or tumor-suppressors or even harboring oncogenic and tumor-suppressing effects; representing a new class of cancer biomarkers and therapeutic targets. It is predicted that more than 35,000 ncRNA especially lncRNA are positioned at the intergenic regions of the human genome. Emerging research indicates that one of the key pathways controlling lncRNA expression and tissue specificity is epigenetic regulation.

METHODS

In the current article, a novel approach for lncRNA discovery based on the intergenic position of most lncRNAs and a single CpG site methylation level representing epigenetic characteristics has been suggested.

RESULTS

Using this method, a novel antisense lncRNA named LINC02892 presenting three transcripts without the capacity of coding a protein was found exhibiting nuclear, cytoplasmic, and exosome distributions.

CONCLUSION

The current discovery strategy could be applied to identify novel non-coding RNAs influenced by methylation aberrations.

NetRAX: accurate and fast maximum likelihood phylogenetic network inference

MOTIVATION

Phylogenetic networks can represent non-treelike evolutionary scenarios. Current, actively developed approaches for phylogenetic network inference jointly account for non-treelike evolution and incomplete lineage sorting (ILS). Unfortunately, this induces a very high computational complexity and current tools can only analyze small datasets.

RESULTS

We present NetRAX, a tool for maximum likelihood (ML) inference of phylogenetic networks in the absence of ILS. Our tool leverages state-of-the-art methods for efficiently computing the phylogenetic likelihood function on trees, and extends them to phylogenetic networks via the notion of 'displayed trees'. NetRAX can infer ML phylogenetic networks from partitioned multiple sequence alignments and returns the inferred networks in Extended Newick format. On simulated data, our results show a very low relative difference in Bayesian Information Criterion (BIC) score and a near-zero unrooted softwired cluster distance to the true, simulated networks. With NetRAX, a network inference on a partitioned alignment with 8000 sites, 30 taxa and 3 reticulations completes within a few minutes on a standard laptop.

AVAILABILITY AND IMPLEMENTATION

Our implementation is available under the GNU General Public License v3.0 at https://github.com/lutteropp/NetRAX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Understanding and Modulating Antibody Fine Specificity: Lessons from Combinatorial Biology

Combinatorial biology methods such as phage and yeast display, suitable for the generation and screening of huge numbers of protein fragments and mutated variants, have been useful when dissecting the molecular details of the interactions between antibodies and their target antigens (mainly those of protein nature). The relevance of these studies goes far beyond the mere description of binding interfaces, as the information obtained has implications for the understanding of the chemistry of antibody–antigen binding reactions and the biological effects of antibodies. Further modification of the interactions through combinatorial methods to manipulate the key properties of antibodies (affinity and fine specificity) can result in the emergence of novel research tools and optimized therapeutics.

Genomic, morphological, and biochemical analyses of a multi-metal resistant but multi-drug susceptible strain of Bordetella petrii from hospital soil

Contamination of soil by antibiotics and heavy metals originating from hospital facilities has emerged as a major cause for the development of resistant microbes. We collected soil samples surrounding a hospital effluent and measured the resistance of bacterial isolates against multiple antibiotics and heavy metals. One strain BMCSI 3 was found to be sensitive to all tested antibiotics. However, it was resistant to many heavy metals and metalloids like cadmium, chromium, copper, mercury, arsenic, and others. This strain was motile and potentially spore-forming. Whole-genome shotgun assembly of BMCSI 3 produced 4.95 Mb genome with 4,638 protein-coding genes. The taxonomic and phylogenetic analysis revealed it, to be a Bordetella petrii strain. Multiple genomic islands carrying mobile genetic elements; coding for heavy metal resistant genes, response regulators or transcription factors, transporters, and multi-drug efflux pumps were identified from the genome. A comparative genomic analysis of BMCSI 3 with annotated genomes of other free-living B. petrii revealed the presence of multiple transposable elements and several genes involved in stress response and metabolism. This study provides insights into how genomic reorganization and plasticity results in evolution of heavy metals resistance by acquiring genes from its natural environment.

Microevolution of Trypanosoma cruzi reveals hybridization and clonal mechanisms driving rapid genome diversification

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi, the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi, followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi.

Restriction of the Global IgM Repertoire in Antiphospholipid Syndrome

The typical anti-phospholipid antibodies (APLA) in the anti-phospholipid syndrome (APS) are reactive with the phospholipid-binding protein β2GPI as well as a growing list of other protein targets. The relation of APLA to natural antibodies and the fuzzy set of autoantigens involved provoked us to study the changes in the IgM repertoire in APS. To this end, peptides selected by serum IgM from a 7-residue linear peptide phage display library (PDL) were deep sequenced. The analysis was aided by a novel formal representation of the Igome (the mimotope set reflecting the IgM specificities) in the form of a sequence graph. The study involved women with APLA and habitual abortions (n=24) compared to age-matched clinically healthy pregnant women (n=20). Their pooled Igomes (297 028 mimotope sequences) were compared also to the global public repertoire Igome of pooled donor plasma IgM (n=2 796 484) and a set of 7-mer sequences found in the J regions of human immunoglobulins (n=4 433 252). The pooled Igome was represented as a graph connecting the sequences as similar as the mimotopes of the same monoclonal antibody. The criterion was based on previously published data. In the resulting graph, identifiable clusters of vertices were considered related to the footprints of overlapping antibody cross-reactivities. A subgraph based on the clusters with a significant differential expression of APS patients’ mimotopes contained predominantly specificities underrepresented in APS. The differentially expressed IgM footprints showed also an increased cross-reactivity with immunoglobulin J regions. The specificities underexpressed in APS had a higher correlation with public specificities than those overexpressed. The APS associated specificities were strongly related also to the human peptidome with 1 072 mimotope sequences found in 7 519 human proteins. These regions were characterized by low complexity. Thus, the IgM repertoire of the APS patients was found to be characterized by a significant reduction of certain public specificities found in the healthy controls with targets representing low complexity linear self-epitopes homologous to human antibody J regions.

Computational Methods for Automated Analysis of Malaria Parasite Using Blood Smear Images: Recent Advances

Malaria comes under one of the dangerous diseases in many countries. It is the primary reason for most of the causalities across the world. It is presently rated as a significant cause of the high mortality rate worldwide compared with other diseases that can be reduced significantly by its earlier detection. Therefore, to facilitate the early detection/diagnosis of malaria to reduce the mortality rate, an automated computational method is required with a high accuracy rate. This study is a solid starting point for researchers who want to look into automated blood smear analysis to detect malaria. In this paper, a comprehensive review of different computer-assisted techniques has been outlined as follows: (i) acquisition of image dataset, (ii) preprocessing, (iii) segmentation of RBC, and (iv) feature extraction and selection, and (v) classification for the detection of malaria parasites using blood smear images. This study will be helpful for: (i) researchers can inspect and improve the existing computational methods for early diagnosis of malaria with a high accuracy rate that may further reduce the interobserver and intra-observer variations; (ii) microbiologists to take the second opinion from the automated computational methods for effective diagnosis of malaria; and (iii) finally, several issues remain addressed, and future work has also been discussed in this work.

Parallel sequence tagging for concept recognition

BACKGROUND

Named Entity Recognition (NER) and Normalisation (NEN) are core components of any text-mining system for biomedical texts. In a traditional concept-recognition pipeline, these tasks are combined in a serial way, which is inherently prone to error propagation from NER to NEN. We propose a parallel architecture, where both NER and NEN are modeled as a sequence-labeling task, operating directly on the source text. We examine different harmonisation strategies for merging the predictions of the two classifiers into a single output sequence.

RESULTS

We test our approach on the recent Version 4 of the CRAFT corpus. In all 20 annotation sets of the concept-annotation task, our system outperforms the pipeline system reported as a baseline in the CRAFT shared task, a competition of the BioNLP Open Shared Tasks 2019. We further refine the systems from the shared task by optimising the harmonisation strategy separately for each annotation set.

CONCLUSIONS

Our analysis shows that the strengths of the two classifiers can be combined in a fruitful way. However, prediction harmonisation requires individual calibration on a development set for each annotation set. This allows achieving a good trade-off between established knowledge (training set) and novel information (unseen concepts).

Glutamic acid reshapes the plant microbiota to protect plants against pathogens

BACKGROUND

Plants in nature interact with other species, among which are mutualistic microorganisms that affect plant health. The co-existence of microbial symbionts with the host contributes to host fitness in a natural context. In turn, the composition of the plant microbiota responds to the environment and the state of the host, raising the possibility that it can be engineered to benefit the plant. However, technology for engineering the structure of the plant microbiome is not yet available.

RESULTS

The loss of diversity and reduction in population density of Streptomyces globisporus SP6C4, a core microbe, was observed coincident with the aging of strawberry plants. Here, we show that glutamic acid reshapes the plant microbial community and enriches populations of Streptomyces, a functional core microbe in the strawberry anthosphere. Similarly, in the tomato rhizosphere, treatment with glutamic acid increased the population sizes of Streptomyces as well as those of Bacillaceae and Burkholderiaceae. At the same time, diseases caused by species of Botrytis and Fusarium were significantly reduced in both habitats. We suggest that glutamic acid directly modulates the composition of the microbiome community.

CONCLUSIONS

Much is known about the structure of plant-associated microbial communities, but less is understood about how the community composition and complexity are controlled. Our results demonstrate that the intrinsic level of glutamic acid in planta is associated with the composition of the microbiota, which can be modulated by an external supply of a biostimulant. Video Abstract.

Identifying stress responsive genes using overlapping communities in co-expression networks

Background

This paper proposes a workflow to identify genes that respond to specific treatments in plants. The workflow takes as input the RNA sequencing read counts and phenotypical data of different genotypes, measured under control and treatment conditions. It outputs a reduced group of genes marked as relevant for treatment response. Technically, the proposed approach is both a generalization and an extension of WGCNA. It aims to identify specific modules of overlapping communities underlying the co-expression network of genes. Module detection is achieved by using Hierarchical Link Clustering. The overlapping nature of the systems’ regulatory domains that generate co-expression can be identified by such modules. LASSO regression is employed to analyze phenotypic responses of modules to treatment.

Results

The workflow is applied to rice (Oryza sativa), a major food source known to be highly sensitive to salt stress. The workflow identifies 19 rice genes that seem relevant in the response to salt stress. They are distributed across 6 modules: 3 modules, each grouping together 3 genes, are associated to shoot K content; 2 modules of 3 genes are associated to shoot biomass; and 1 module of 4 genes is associated to root biomass. These genes represent target genes for the improvement of salinity tolerance in rice.

Conclusions

A more effective framework to reduce the search-space for target genes that respond to a specific treatment is introduced. It facilitates experimental validation by restraining efforts to a smaller subset of genes of high potential relevance.

Trends in biological data integration for the selection of enzymes and transcription factors related to cellulose and hemicellulose degradation in fungi

Fungi are key players in biotechnological applications. Although several studies focusing on fungal diversity and genetics have been performed, many details of fungal biology remain unknown, including how cellulolytic enzymes are modulated within these organisms to allow changes in main plant cell wall compounds, cellulose and hemicellulose, and subsequent biomass conversion. With the advent and consolidation of DNA/RNA sequencing technology, different types of information can be generated at the genomic, structural and functional levels, including the gene expression profiles and regulatory mechanisms of these organisms, during degradation-induced conditions. This increase in data generation made rapid computational development necessary to deal with the large amounts of data generated. In this context, the origination of bioinformatics, a hybrid science integrating biological data with various techniques for information storage, distribution and analysis, was a fundamental step toward the current state-of-the-art in the postgenomic era. The possibility of integrating biological big data has facilitated exciting discoveries, including identifying novel mechanisms and more efficient enzymes, increasing yields, reducing costs and expanding opportunities in the bioprocess field. In this review, we summarize the current status and trends of the integration of different types of biological data through bioinformatics approaches for biological data analysis and enzyme selection.

Metapredict: a fast, accurate, and easy-to-use predictor of consensus disorder and structure

Intrinsically disordered proteins and protein regions make up a substantial fraction of many proteomes in which they play a wide variety of essential roles. A critical first step in understanding the role of disordered protein regions in biological function is to identify those disordered regions correctly. Computational methods for disorder prediction have emerged as a core set of tools to guide experiments, interpret results, and develop hypotheses. Given the multiple different predictors available, consensus scores have emerged as a popular approach to mitigate biases or limitations of any single method. Consensus scores integrate the outcome of multiple independent disorder predictors and provide a per-residue value that reflects the number of tools that predict a residue to be disordered. Although consensus scores help mitigate the inherent problems of using any single disorder predictor, they are computationally expensive to generate. They also necessitate the installation of multiple different software tools, which can be prohibitively difficult. To address this challenge, we developed a deep-learning-based predictor of consensus disorder scores. Our predictor, metapredict, utilizes a bidirectional recurrent neural network trained on the consensus disorder scores from 12 proteomes. By benchmarking metapredict using two orthogonal approaches, we found that metapredict is among the most accurate disorder predictors currently available. Metapredict is also remarkably fast, enabling proteome-scale disorder prediction in minutes. Importantly, metapredict is a fully open source and is distributed as a Python package, a collection of command-line tools, and a web server, maximizing the potential practical utility of the predictor. We believe metapredict offers a convenient, accessible, accurate, and high-performance predictor for single-proteins and proteomes alike.

Naive human B cells engage the receptor binding domain of SARS-CoV-2, variants of concern, and related sarbecoviruses

Exposure to a pathogen elicits an adaptive immune response aimed to control and eradicate. Interrogating the abundance and specificity of the naive B cell repertoire contributes to understanding how to potentially elicit protective responses. Here, we isolated naive B cells from 8 seronegative human donors targeting the SARS-CoV-2 receptor-binding domain (RBD). Single B cell analysis showed diverse gene usage with no restricted complementarity determining region lengths. We show that recombinant antibodies engage SARS-CoV-2 RBD, circulating variants, and pre-emergent coronaviruses. Representative antibodies signal in a B cell activation assay and can be affinity matured through directed evolution. Structural analysis of a naive antibody in complex with spike shows a conserved mode of recognition shared with infection-induced antibodies. Lastly, both naive and affinity-matured antibodies can neutralize SARS-CoV-2. Understanding the naive repertoire may inform potential responses recognizing variants or emerging coronaviruses enabling the development of pan-coronavirus vaccines aimed at engaging germline responses.

Isolation of antibody germline precursors targeting the receptor binding domain of coronaviruses.

Towards defining the Mechanisms of Alzheimer's disease based on a contextual analysis of molecular pathways

Alzheimer's disease (AD) is posing an increasingly profound problem to society. Our genuine understanding of the pathogenesis of AD is inadequate and as a consequence, diagnostic and therapeutic strategies are currently insufficient. The understandable focus of many studies is the identification of molecules with high diagnostic utility however the opportunity to obtain a further understanding of the mechanistic origins of the disease from such putative biomarkers is often overlooked. This study examines the involvement of biomarkers in AD to shed light on potential mechanisms and pathways through which they are implicated in the pathology of this devastating neurodegenerative disorder. The computational tools required to analyse ever-growing datasets in the context of AD are also discussed.

EARN: an ensemble machine learning algorithm to predict driver genes in metastatic breast cancer

BACKGROUND

Today, there are a lot of markers on the prognosis and diagnosis of complex diseases such as primary breast cancer. However, our understanding of the drivers that influence cancer aggression is limited.

METHODS

In this work, we study somatic mutation data consists of 450 metastatic breast tumor samples from cBio Cancer Genomics Portal. We use four software tools to extract features from this data. Then, an ensemble classifier (EC) learning algorithm called EARN (Ensemble of Artificial Neural Network, Random Forest, and non-linear Support Vector Machine) is proposed to evaluate plausible driver genes for metastatic breast cancer (MBCA). The decision-making strategy for the proposed ensemble machine is based on the aggregation of the predicted scores obtained from individual learning classifiers to be prioritized homo sapiens genes annotated as protein-coding from NCBI.

RESULTS

This study is an attempt to focus on the findings in several aspects of MBCA prognosis and diagnosis. First, drivers and passengers predicted by SVM, ANN, RF, and EARN are introduced. Second, biological inferences of predictions are discussed based on gene set enrichment analysis. Third, statistical validation and comparison of all learning methods are performed by some evaluation metrics. Finally, the pathway enrichment analysis (PEA) using ReactomeFIVIz tool (FDR < 0.03) for the top 100 genes predicted by EARN leads us to propose a new gene set panel for MBCA. It includes HDAC3, ABAT, GRIN1, PLCB1, and KPNA2 as well as NCOR1, TBL1XR1, SIRT4, KRAS, CACNA1E, PRKCG, GPS2, SIN3A, ACTB, KDM6B, and PRMT1. Furthermore, we compare results for MBCA to other outputs regarding 983 primary tumor samples of breast invasive carcinoma (BRCA) obtained from the Cancer Genome Atlas (TCGA). The comparison between outputs shows that ROC-AUC reaches 99.24% using EARN for MBCA and 99.79% for BRCA. This statistical result is better than three individual classifiers in each case.

CONCLUSIONS

This research using an integrative approach assists precision oncologists to design compact targeted panels that eliminate the need for whole-genome/exome sequencing. The schematic representation of the proposed model is presented as the Graphic abstract.

Enhancement of conformational B-cell epitope prediction using CluSMOTE

BACKGROUND

A conformational B-cell epitope is one of the main components of vaccine design. It contains separate segments in its sequence, which are spatially close in the antigen chain. The availability of Ag-Ab complex data on the Protein Data Bank allows for the development predictive methods. Several epitope prediction models also have been developed, including learning-based methods. However, the performance of the model is still not optimum. The main problem in learning-based prediction models is class imbalance.

METHODS

This study proposes CluSMOTE, which is a combination of a cluster-based undersampling method and Synthetic Minority Oversampling Technique. The approach is used to generate other sample data to ensure that the dataset of the conformational epitope is balanced. The Hierarchical DBSCAN algorithm is performed to identify the cluster in the majority class. Some of the randomly selected data is taken from each cluster, considering the oversampling degree, and combined with the minority class data. The balance data is utilized as the training dataset to develop a conformational epitope prediction. Furthermore, two binary classification methods, Support Vector Machine and Decision Tree, are separately used to develop model prediction and to evaluate the performance of CluSMOTE in predicting conformational B-cell epitope. The experiment is focused on determining the best parameter for optimal CluSMOTE. Two independent datasets are used to compare the proposed prediction model with state of the art methods. The first and the second datasets represent the general protein and the glycoprotein antigens respectively.

RESULT

The experimental result shows that CluSMOTE Decision Tree outperformed the Support Vector Machine in terms of AUC and Gmean as performance measurements. The mean AUC of CluSMOTE Decision Tree in the Kringelum and the SEPPA 3 test sets are 0.83 and 0.766, respectively. This shows that CluSMOTE Decision Tree is better than other methods in the general protein antigen, though comparable with SEPPA 3 in the glycoprotein antigen.

Antibody characterization using immunosignatures

Therapeutic monoclonal antibodies have the potential to work as biological therapeutics. OKT3, Herceptin, Keytruda and others have positively impacted healthcare. Antibodies evolved naturally to provide high specificity and high affinity once mature. These characteristics can make them useful as therapeutics. However, we may be missing characteristics that are not obvious. We present a means of measuring antibodies in an unbiased manner that may highlight therapeutic activity. We propose using a microarray of random peptides to assess antibody properties. We tested twenty-four different commercial antibodies to gain some perspective about how much information can be derived from binding antibodies to random peptide libraries. Some monoclonals preferred to bind shorter peptides, some longer, some preferred motifs closer to the C-term, some nearer the N-term. We tested some antibodies with clinical activity but whose function was blinded to us at the time. We were provided with twenty-one different monoclonal antibodies, thirteen mouse and eight human IgM. These antibodies produced a variety of binding patterns on the random peptide arrays. When unblinded, the antibodies with polyspecific binding were the ones with the greatest therapeutic activity. The protein target to these therapeutic monoclonals is still unknown but using common sequence motifs from the peptides we predicted several human and mouse proteins. The same five highest proteins appeared in both mouse and human lists.

Development and Application of Computational Methods in Phage Display Technology

BACKGROUND

Phage display is a powerful and versatile technology for the identification of peptide ligands binding to multiple targets, which has been successfully employed in various fields, such as diagnostics and therapeutics, drug-delivery and material science. The integration of next generation sequencing technology with phage display makes this methodology more productive. With the widespread use of this technique and the fast accumulation of phage display data, databases for these data and computational methods have become an indispensable part in this community. This review aims to summarize and discuss recent progress in the development and application of computational methods in the field of phage display.

METHODS

We undertook a comprehensive search of bioinformatics resources and computational methods for phage display data via Google Scholar and PubMed. The methods and tools were further divided into different categories according to their uses.

RESULTS

We described seven special or relevant databases for phage display data, which provided an evidence-based source for phage display researchers to clean their biopanning results. These databases can identify and report possible target-unrelated peptides (TUPs), thereby excluding false-positive data from peptides obtained from phage display screening experiments. More than 20 computational methods for analyzing biopanning data were also reviewed. These methods were classified into computational methods for reporting TUPs, for predicting epitopes and for analyzing next generation phage display data.

CONCLUSION

The current bioinformatics archives, methods and tools reviewed here have benefitted the biopanning community. To develop better or new computational tools, some promising directions are also discussed.

A general approach for predicting protein epitopes targeted by antibody repertoires using whole proteomes

Antibodies are essential to functional immunity, yet the epitopes targeted by antibody repertoires remain largely uncharacterized. To aid in characterization, we developed a generalizable strategy to predict antibody-binding epitopes within individual proteins and entire proteomes. Specifically, we selected antibody-binding peptides for 273 distinct sera out of a random library and identified the peptides using next-generation sequencing. To predict antibody-binding epitopes and the antigens from which these epitopes were derived, we tiled the sequences of candidate antigens into short overlapping subsequences of length k (k-mers). We used the enrichment over background of these k-mers in the antibody-binding peptide dataset to predict antibody-binding epitopes. As a positive control, we used this approach, termed K-mer Tiling of Protein Epitopes (K-TOPE), to predict epitopes targeted by monoclonal and polyclonal antibodies of well-characterized specificity, accurately recovering their known epitopes. K-TOPE characterized a commonly targeted antigen from Rhinovirus A, predicting four epitopes recognized by antibodies present in 87% of sera (n = 250). An analysis of 2,908 proteins from 400 viral taxa that infect humans predicted seven enterovirus epitopes and five Epstein-Barr virus epitopes recognized by >30% of specimens. Analysis of Staphylococcus and Streptococcus proteomes similarly predicted 22 epitopes recognized by >30% of specimens. Twelve of these common viral and bacterial epitopes agreed with previously mapped epitopes with p-values < 0.05. Additionally, we predicted 30 HSV2-specific epitopes that were 100% specific against HSV1 in novel and previously reported antigens. Experimentally validating these candidate epitopes could help identify diagnostic biomarkers, vaccine components, and therapeutic targets. The K-TOPE approach thus provides a powerful new tool to elucidate the organisms, antigens, and epitopes targeted by human antibody repertoires.

Multiscale modeling reveals angiogenesis-induced drug resistance in brain tumors and predicts a synergistic drug combination targeting EGFR and VEGFR pathways

BACKGROUND

Experimental studies have demonstrated that both the extracellular vasculature or microenvironment and intracellular molecular network (e.g., epidermal growth factor receptor (EGFR) signaling pathway) are important for brain tumor growth. Additionally, some drugs have been developed to inhibit EGFR signaling pathways. However, how angiogenesis affects the response of tumor cells to drug treatment has rarely been mechanistically studied. Therefore, a multiscale model is required to investigate such complex biological systems that contain interactions and feedback among multiple levels.

RESULTS

In this study, we developed a single cell-based multiscale spatiotemporal model to simulate vascular tumor growth and the drug response based on the vascular endothelial growth factor receptor (VEGFR) signaling pathway, the EGFR signaling pathway and the cell cycle as well as several microenvironmental factors that determine cell fate switches in a temporal and spatial context. By incorporating the EGFRI treatment effect, the model showed an interesting phenomenon in which the survival rate of tumor cells decreased in the early stage but rebounded in a later stage, revealing the emergence of drug resistance. Moreover, we revealed the critical role of angiogenesis in acquired drug resistance, since inhibiting blood vessel growth using a VEGFR inhibitor prevented the recovery of the survival rate of tumor cells in the later stage. We further investigated the optimal timing of combining VEGFR inhibition with EGFR inhibition and predicted that the drug combination targeting both the EGFR pathway and VEGFR pathway has a synergistic effect. The experimental data validated the prediction of drug synergy, confirming the effectiveness of our model. In addition, the combination of EGFR and VEGFR genes showed clinical relevance in glioma patients.

CONCLUSIONS

The developed multiscale model revealed angiogenesis-induced drug resistance mechanisms of brain tumors to EGFRI treatment and predicted a synergistic drug combination targeting both EGFR and VEGFR pathways with optimal combination timing. This study explored the mechanistic and functional mechanisms of the angiogenesis underlying tumor growth and drug resistance, which advances our understanding of novel mechanisms of drug resistance and provides implications for designing more effective cancer therapies.

A large-scale comparative assessment of methods for residue-residue contact prediction

Sequence-based prediction of residue-residue contact in proteins becomes increasingly more important for improving protein structure prediction in the big data era. In this study, we performed a large-scale comparative assessment of 15 locally installed contact predictors. To assess these methods, we collected a big data set consisting of 680 nonredundant proteins covering different structural classes and target difficulties. We investigated a wide range of factors that may influence the precision of contact prediction, including target difficulty, structural class, the alignment depth and distribution of contact pairs in a protein structure. We found that: (1) the machine learning-based methods outperform the direct-coupling-based methods for short-range contact prediction, while the latter are significantly better for long-range contact prediction. The consensus-based methods, which combine machine learning and direct-coupling methods, perform the best. (2) The target difficulty does not have clear influence on the machine learning-based methods, while it does affect the direct-coupling and consensus-based methods significantly. (3) The alignment depth has relatively weak effect on the machine learning-based methods. However, for the direct-coupling-based methods and consensus-based methods, the predicted contacts for targets with deeper alignment tend to be more accurate. (4) All methods perform relatively better on β and α + β proteins than on α proteins. (5) Residues buried in the core of protein structure are more prone to be in contact than residues on the surface (22 versus 6%). We believe these are useful results for guiding future development of new approach to contact prediction.

High-throughput sequencing of small RNAs revealed the diversified cold-responsive pathways during cold stress in the wild banana (Musa itinerans)

BACKGROUND

Cold stress is one of the most severe abiotic stresses affecting the banana production. Although some miRNAs have been identified, little is known about the role of miRNAs in response to cold stress in banana, and up to date, there is no report about the role of miRNAs in the response to cold stress in the plants of the cultivated or wild bananas.

RESULT

Here, a cold-resistant line wild banana (Musa itinerans) from China was used to profile the cold-responsive miRNAs by RNA-seq during cold stress. Totally, 265 known mature miRNAs and 41 novel miRNAs were obtained. Cluster analysis of differentially expressed (DE) miRNAs indicated that some miRNAs were specific for chilling or 0 °C treated responses, and most of them were reported to be cold-responsive; however, some were seldom reported to be cold-responsive in response to cold stress, e.g., miR395, miR408, miR172, suggesting that they maybe play key roles in response to cold stress. The GO and KEGG pathway enrichment analysis of DE miRNAs targets indicated that there existed diversified cold-responsive pathways, and miR172 was found likely to play a central coordinating role in response to cold stress, especially in the regulation of CK2 and the circadian rhythm. Finally, qPCR assays indicated the related targets were negatively regulated by the tested DE miRNAs during cold stress in the wild banana.

CONCLUSIONS

In this study, the profiling of miRNAs by RNA-seq in response to cold stress in the plants of the wild banana (Musa itinerans) was reported for the first time. The results showed that there existed diversified cold-responsive pathways, which provided insight into the roles of miRNAs during cold stress, and would be helpful for alleviating cold stress and cold-resistant breeding in bananas.

Adaptive optimal input design and parametric estimation of nonlinear dynamical systems: application to neuronal modeling

OBJECTIVE

Many physical models of biological processes including neural systems are characterized by parametric nonlinear dynamical relations between driving inputs, internal states, and measured outputs of the process. Fitting such models using experimental data (data assimilation) is a challenging task since the physical process often operates in a noisy, possibly non-stationary environment; moreover, conducting multiple experiments under controlled and repeatable conditions can be impractical, time consuming or costly. The accuracy of model identification, therefore, is dictated principally by the quality and dynamic richness of collected data over single or few experimental sessions. Accordingly, it is highly desirable to design efficient experiments that, by exciting the physical process with smart inputs, yields fast convergence and increased accuracy of the model.

APPROACH

We herein introduce an adaptive framework in which optimal input design is integrated with square root cubature Kalman filters (OID-SCKF) to develop an online estimation procedure that first, converges significantly quicker, thereby permitting model fitting over shorter time windows, and second, enhances model accuracy when only few process outputs are accessible. The methodology is demonstrated on common nonlinear models and on a four-area neural mass model with noisy and limited measurements. Estimation quality (speed and accuracy) is benchmarked against high-performance SCKF-based methods that commonly employ dynamically rich informed inputs for accurate model identification.

MAIN RESULTS

For all the tested models, simulated single-trial and ensemble averages showed that OID-SCKF exhibited (i) faster convergence of parameter estimates and (ii) lower dependence on inter-trial noise variability with gains up to around 1000 ms in speed and 81% increase in variability for the neural mass models. In terms of accuracy, OID-SCKF estimation was superior, and exhibited considerably less variability across experiments, in identifying model parameters of (a) systems with challenging model inversion dynamics and (b) systems with fewer measurable outputs that directly relate to the underlying processes.

SIGNIFICANCE

Fast and accurate identification therefore carries particular promise for modeling of transient (short-lived) neuronal network dynamics using a spatially under-sampled set of noisy measurements, as is commonly encountered in neural engineering applications.

Improving the interoperability of biomedical ontologies with compound alignments

Background

Ontologies are commonly used to annotate and help process life sciences data. Although their original goal is to facilitate integration and interoperability among heterogeneous data sources, when these sources are annotated with distinct ontologies, bridging this gap can be challenging. In the last decade, ontology matching systems have been evolving and are now capable of producing high-quality mappings for life sciences ontologies, usually limited to the equivalence between two ontologies. However, life sciences research is becoming increasingly transdisciplinary and integrative, fostering the need to develop matching strategies that are able to handle multiple ontologies and more complex relations between their concepts.

Results

We have developed ontology matching algorithms that are able to find compound mappings between multiple biomedical ontologies, in the form of ternary mappings, finding for instance that “aortic valve stenosis”(HP:0001650) is equivalent to the intersection between “aortic valve”(FMA:7236) and “constricted” (PATO:0001847). The algorithms take advantage of search space filtering based on partial mappings between ontology pairs, to be able to handle the increased computational demands. The evaluation of the algorithms has shown that they are able to produce meaningful results, with precision in the range of 60-92% for new mappings. The algorithms were also applied to the potential extension of logical definitions of the OBO and the matching of several plant-related ontologies.

Conclusions

This work is a first step towards finding more complex relations between multiple ontologies. The evaluation shows that the results produced are significant and that the algorithms could satisfy specific integration needs.

A novel non-parametric method for uncertainty evaluation of correlation-based molecular signatures: its application on PAM50 algorithm

Motivation

The PAM50 classifier is used to assign patients to the highest correlated breast cancer subtype irrespectively of the obtained value. Nonetheless, all subtype correlations are required to build the risk of recurrence (ROR) score, currently used in therapeutic decisions. Present subtype uncertainty estimations are not accurate, seldom considered or require a population-based approach for this context.

Results

Here we present a novel single-subject non-parametric uncertainty estimation based on PAM50's gene label permutations. Simulations results ( n  = 5228) showed that only 61% subjects can be reliably 'Assigned' to the PAM50 subtype, whereas 33% should be 'Not Assigned' (NA), leaving the rest to tight 'Ambiguous' correlations between subtypes. The NA subjects exclusion from the analysis improved survival subtype curves discrimination yielding a higher proportion of low and high ROR values. Conversely, all NA subjects showed similar survival behaviour regardless of the original PAM50 assignment. We propose to incorporate our PAM50 uncertainty estimation to support therapeutic decisions.

Availability and Implementation

Source code can be found in 'pbcmc' R package at Bioconductor.

Contacts

cristobalfresno@gmail.com or efernandez@bdmg.com.ar.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Emergence of Openness in Open-Source Projects: The Case of OpenEhR

The meaning of openness in open source is both intrinsically unstable and dynamic, and tends to fluctuate with time and context. We draw on a very particular open-source project primarily concerned with building rigorous clinical concepts to be used in electronic health records called openEHR. openEHR explains how openness is a concept that is purposely engaged with, and how, in this process of engagement, the very meaning of open matures and evolves within the project. Drawing on rich longitudinal data related to openEHR we theorise the evolving nature of openness and how this idea emerges through two intertwined processes of maturation and metamorphosis. While metamorphosis allows us to trace and interrogate the mutational evolution in openness, maturation analyses the small, careful changes crafted to build a very particular understanding of openness. Metamorphosis is less managed and controlled, whereas maturation is representative of highly precise work carried out in controlled form. Both processes work together in open-source projects and reinforce each other. Our study reveals that openness emerges and evolves in open-source projects where it can be understood to mean rigour; ability to participate; open implementation; and an open process. Our work contributes to a deepening in the theorisation of what it means to be an open-source project. The multiple and co-existing meanings of ‘open’ imply that open-source projects evolve in nonlinear ways where each critical meaning of openness causes a reflective questioning by the community of its continued status and existence.

LT adjuvant modulates epitope specificity and improves the efficacy of murine antibodies elicited by sublingual vaccination with the N-terminal domain of Streptococcus mutans P1

In this study, we evaluated the immunogenicity, protective efficacy and peptide-based immune signatures of antibodies raised in mice after sublingual immunization with a recombinant form of the P1 (aka AgI/II, PAc) adhesin (P1_39-512) of Streptococcus mutans, a major etiological agent of dental caries. Sublingual administration of P1_39-512 in combination with the mucosal adjuvant LTK4R (a derivative of heat-labile LT toxin) induced strong and long-lasting systemic and mucosal immune responses. Incorporation of the adjuvant resulted in an enhancement of the anti-adhesive and anti-colonization activity against S. mutans as evaluated both under in vitro and in vivo conditions. Incorporation of the adjuvant to the vaccine formulation also changed the epitope specificity of the induced antibodies as determined by immunological signatures of sera collected from vaccinated mice. Use of a peptide microarray library led to the identification of peptide targets recognized by antibodies in serum samples with enhanced anti-adhesive effects. Altogether, the results presented herein showed that the sublingual administration of a P1-based subunit vaccine represents a promising approach for the prevention of dental caries caused by S. mutans. In addition, the present study disclosed the role of adjuvants on the epitope specificity and functionality of antibodies raised by subunit vaccines.

Adjuvant-Mediated Epitope Specificity and Enhanced Neutralizing Activity of Antibodies Targeting Dengue Virus Envelope Protein

The heat-labile toxins (LT) produced by enterotoxigenic Escherichia coli display adjuvant effects to coadministered antigens, leading to enhanced production of serum antibodies. Despite extensive knowledge of the adjuvant properties of LT derivatives, including in vitro-generated non-toxic mutant forms, little is known about the capacity of these adjuvants to modulate the epitope specificity of antibodies directed against antigens. This study characterizes the role of LT and its non-toxic B subunit (LTB) in the modulation of antibody responses to a coadministered antigen, the dengue virus (DENV) envelope glycoprotein domain III (EDIII), which binds to surface receptors and mediates virus entry into host cells. In contrast to non-adjuvanted or alum-adjuvanted formulations, antibodies induced in mice immunized with LT or LTB showed enhanced virus-neutralization effects that were not ascribed to a subclass shift or antigen affinity. Nonetheless, immunosignature analyses revealed that purified LT-adjuvanted EDIII-specific antibodies display distinct epitope-binding patterns with regard to antibodies raised in mice immunized with EDIII or the alum-adjuvanted vaccine. Notably, the analyses led to the identification of a specific EDIII epitope located in the EF to FG loop, which is involved in the entry of DENV into eukaryotic cells. The present results demonstrate that LT and LTB modulate the epitope specificity of antibodies generated after immunization with coadministered antigens that, in the case of EDIII, was associated with the induction of neutralizing antibody responses. These results open perspectives for the more rational development of vaccines with enhanced protective effects against DENV infections.

DRNApred, fast sequence-based method that accurately predicts and discriminates DNA- and RNA-binding residues

Protein-DNA and protein-RNA interactions are part of many diverse and essential cellular functions and yet most of them remain to be discovered and characterized. Recent research shows that sequence-based predictors of DNA-binding residues accurately find these residues but also cross-predict many RNA-binding residues as DNA-binding, and vice versa. Most of these methods are also relatively slow, prohibiting applications on the whole-genome scale. We describe a novel sequence-based method, DRNApred, which accurately and in high-throughput predicts and discriminates between DNA- and RNA-binding residues. DRNApred was designed using a new dataset with both DNA- and RNA-binding proteins, regression that penalizes cross-predictions, and a novel two-layered architecture. DRNApred outperforms state-of-the-art predictors of DNA- or RNA-binding residues on a benchmark test dataset by substantially reducing the cross predictions and predicting arguably higher quality false positives that are located nearby the native binding residues. Moreover, it also more accurately predicts the DNA- and RNA-binding proteins. Application on the human proteome confirms that DRNApred reduces the cross predictions among the native nucleic acid binders. Also, novel putative DNA/RNA-binding proteins that it predicts share similar subcellular locations and residue charge profiles with the known native binding proteins. Webserver of DRNApred is freely available at http://biomine.cs.vcu.edu/servers/DRNApred/.

Equine skeletal muscle adaptations to exercise and training: evidence of differential regulation of autophagosomal and mitochondrial components

BACKGROUND

A single bout of exercise induces changes in gene expression in skeletal muscle. Regular exercise results in an adaptive response involving changes in muscle architecture and biochemistry, and is an effective way to manage and prevent common human diseases such as obesity, cardiovascular disorders and type II diabetes. However, the biomolecular mechanisms underlying such responses still need to be fully elucidated. Here we performed a transcriptome-wide analysis of skeletal muscle tissue in a large cohort of untrained Thoroughbred horses (n = 51) before and after a bout of high-intensity exercise and again after an extended period of training. We hypothesized that regular high-intensity exercise training primes the transcriptome for the demands of high-intensity exercise.

RESULTS

An extensive set of genes was observed to be significantly differentially regulated in response to a single bout of high-intensity exercise in the untrained cohort (3241 genes) and following multiple bouts of high-intensity exercise training over a six-month period (3405 genes). Approximately one-third of these genes (1025) and several biological processes related to energy metabolism were common to both the exercise and training responses. We then developed a novel network-based computational analysis pipeline to test the hypothesis that these transcriptional changes also influence the contextual molecular interactome and its dynamics in response to exercise and training. The contextual network analysis identified several important hub genes, including the autophagosomal-related gene GABARAPL1, and dynamic functional modules, including those enriched for mitochondrial respiratory chain complexes I and V, that were differentially regulated and had their putative interactions 're-wired' in the exercise and/or training responses.

CONCLUSION

Here we have generated for the first time, a comprehensive set of genes that are differentially expressed in Thoroughbred skeletal muscle in response to both exercise and training. These data indicate that consecutive bouts of high-intensity exercise result in a priming of the skeletal muscle transcriptome for the demands of the next exercise bout. Furthermore, this may also lead to an extensive 're-wiring' of the molecular interactome in both exercise and training and include key genes and functional modules related to autophagy and the mitochondrion.