In silico prediction of splice-altering single nucleotide variants in the human genome

In silico tools have been developed to predict variants that may have an impact on pre-mRNA splicing. The major limitation of the application of these tools to basic research and clinical practice is the difficulty in interpreting the output. Most tools only predict potential splice sites given a DNA sequence without measuring splicing signal changes caused by a variant. Another limitation is the lack of large-scale evaluation studies of these tools. We compared eight in silico tools on 2959 single nucleotide variants within splicing consensus regions (scSNVs) using receiver operating characteristic analysis. The Position Weight Matrix model and MaxEntScan outperformed other methods. Two ensemble learning methods, adaptive boosting and random forests, were used to construct models that take advantage of individual methods. Both models further improved prediction, with outputs of directly interpretable prediction scores. We applied our ensemble scores to scSNVs from the Catalogue of Somatic Mutations in Cancer database. Analysis showed that predicted splice-altering scSNVs are enriched in recurrent scSNVs and known cancer genes. We pre-computed our ensemble scores for all potential scSNVs across the human genome, providing a whole genome level resource for identifying splice-altering scSNVs discovered from large-scale sequencing studies.     

In Vitro Cultivation of the Hymenoptera Genetic Model,Nasonia

The wasp genus Nasonia is a genetic model with unique advantages for the study of interspecific differences, including haplodiploidy and interfertile species. However, as a parasitoid, Nasonia is confined within a fly host, thus restricting direct observations and manipulation of development over time. Here, we present the first in vitro cultivation method for this system that decouples Nasonia from its host, allowing continuous observations from embryo to adulthood. Using transwell plates and a simple Nasonia rearing medium, we demonstrate a technique that will significantly expand the utility of the Nasonia model.     

In Vitro Model of Tumor Cell Extravasation

Tumor cells that disseminate from the primary tumor and survive the vascular system can eventually extravasate across the endothelium to metastasize at a secondary site. In this study, we developed a microfluidic system to mimic tumor cell extravasation where cancer cells can transmigrate across an endothelial monolayer into a hydrogel that models the extracellular space. The experimental protocol is optimized to ensure the formation of an intact endothelium prior to the introduction of tumor cells and also to observe tumor cell extravasation by having a suitable tumor seeding density. Extravasation is observed for 38.8% of the tumor cells in contact with the endothelium within 1 day after their introduction. Permeability of the EC monolayer as measured by the diffusion of fluorescently-labeled dextran across the monolayer increased 3.8 fold 24 hours after introducing tumor cells, suggesting that the presence of tumor cells increases endothelial permeability. The percent of tumor cells extravasated remained nearly constant from1 to 3 days after tumor seeding, indicating extravasation in our system generally occurs within the first 24 hours of tumor cell contact with the endothelium.     

In Vivo Assessment of Bone Regeneration in Alginate/Bone ECM Hydrogels with Incorporated Skeletal Stem Cells and Single Growth Factors

The current study has investigated the use of decellularised, demineralised bone extracellular matrix (ECM) hydrogel constructs for in vivo tissue mineralisation and bone formation. Stro-1-enriched human bone marrow stromal cells were incorporated together with select growth factors including VEGF, TGF-β3, BMP-2, PTHrP and VitD3, to augment bone formation, and mixed with alginate for structural support. Growth factors were delivered through fast (non-osteogenic factors) and slow (osteogenic factors) release PLGA microparticles. Constructs of 5 mm length were implanted in vivo for 28 days within mice. Dense tissue assessed by micro-CT correlated with histologically assessed mineralised bone formation in all constructs. Exogenous growth factor addition did not enhance bone formation further compared to alginate/bone ECM (ALG/ECM) hydrogels alone. UV irradiation reduced bone formation through degradation of intrinsic growth factors within the bone ECM component and possibly also ECM cross-linking. BMP-2 and VitD3 rescued osteogenic induction. ALG/ECM hydrogels appeared highly osteoinductive and delivery of angiogenic or chondrogenic growth factors led to altered bone formation. All constructs demonstrated extensive host tissue invasion and vascularisation aiding integration and implant longevity. The proposed hydrogel system functioned without the need for growth factor incorporation or an exogenous inducible cell source. Optimal growth factor concentrations and spatiotemporal release profiles require further assessment, as the bone ECM component may suffer batch variability between donor materials. In summary, ALG/ECM hydrogels provide a versatile biomaterial scaffold for utilisation within regenerative medicine which may be tailored, ultimately, to form the tissue of choice through incorporation of select growth factors.     

In silico Platform for Prediction of N-, O- and C-Glycosites in Eukaryotic Protein Sequences

Glycosylation is one of the most abundant and an important post-translational modification of proteins. Glycosylated proteins (glycoproteins) are involved in various cellular biological functions like protein folding, cell-cell interactions, cell recognition and host-pathogen interactions. A large number of eukaryotic glycoproteins also have therapeutic and potential technology applications. Therefore, characterization and analysis of glycosites (glycosylated residues) in these proteins is of great interest to biologists. In order to cater these needs a number of in silico tools have been developed over the years, however, a need to get even better prediction tools remains. Therefore, in this study we have developed a new webserver GlycoEP for more accurate prediction of N-linked, O-linked and C-linked glycosites in eukaryotic glycoproteins using two larger datasets, namely, standard and advanced datasets. In case of standard datasets no two glycosylated proteins are more similar than 40%; advanced datasets are highly non-redundant where no two glycosites’ patterns (as defined in methods) have more than 60% similarity. Further, based on our results with several algorihtms developed using different machine-learning techniques, we found Support Vector Machine (SVM) as optimum tool to develop glycosite prediction models. Accordingly, using our more stringent and non-redundant advanced datasets, the SVM based models developed in this study achieved a prediction accuracy of 84.26%, 86.87% and 91.43% with corresponding MCC of 0.54, 0.20 and 0.78, for N-, O- and C-linked glycosites, respectively. The best performing models trained on advanced datasets were then implemented as a user-friendly web server GlycoEP (http://www.imtech.res.in/raghava/glycoep/). Additionally, this server provides prediction models developed on standard datasets and allows users to scan sequons in input protein sequences.     

The Angiopoietin-1 Variant COMP-Ang1 Enhances BMP2-Induced Bone Regeneration with Recruiting Pericytes in Critical Sized Calvarial Defects

Craniofacial bone defects are observed in a variety of clinical situations, and their reconstructions require coordinated coupling between angiogenesis and osteogenesis. In this study, we explored the effects of cartilage oligomeric matrix protein-angiopoietin 1 (COMP-Ang1), a synthetic and soluble variant of angiopoietin 1, on bone morphogenetic protein 2 (BMP2)-induced cranial bone regeneration, and recruitment and osteogenic differentiation of perivascular pericytes. A critical-size calvarial defect was created in the C57BL/6 mouse and COMP-Ang1 and/or BMP2 proteins were delivered into the defects with absorbable collagen sponges. After 3 weeks, bone regeneration was evaluated using micro-computed tomography and histologic examination. Pericyte recruitment into the defects was examined using immunofluorescence staining with anti-NG2 and anti-CD31 antibodies. In vitro recruitment and osteoblastic differentiation of pericyte cells were assessed with Boyden chamber assay, staining of calcified nodules, RT-PCR and Western blot analyses. Combined administration of COMP-Ang1 and BMP2 synergistically enhanced bone repair along with the increased population of CD31 (an endothelial cell marker) and NG2 (a specific marker of pericyte) positive cells. In vitro cultures of pericytes consistently showed that pericyte infiltration into the membrane pore of Boyden chamber was more enhanced by the combination treatment. In addition, the combination further increased the osteoblast-specific gene expression, including bone sialoprotein (BSP), osteocalcin (OCN) and osterix (OSX), phosphorylation of Smad/1/5/8, and mineralized nodule formation. COMP-Ang1 can enhance BMP2-induced cranial bone regeneration with increased pericyte recruitment. Combined delivery of the proteins might be a therapeutic strategy to repair cranial bone damage.     

Effect of Starvation on the Endocytic Pathway in Dictyostelium Cells

Dictyostelium discoideum amoebae have been used extensively to study the structure and dynamics of the endocytic pathway. Here, we show that while the general structure of the endocytic pathway is maintained in starved cells, its dynamics rapidly slow down. In addition, analysis of apm3 and lvsB mutants reveals that the functional organization of the endocytic pathway is profoundly modified upon starvation. Indeed, in these mutant cells, some of the defects observed in rich medium persist in starved cells, notably an abnormally slow transfer of endocytosed material between endocytic compartments. Other parameters, such as endocytosis of the fluid phase or the rate of fusion of postlysosomes to the cell surface, vary dramatically upon starvation. Studying the endocytic pathway in starved cells can provide a different perspective, allowing the primary (invariant) defects resulting from specific mutations to be distinguished from their secondary (conditional) consequences.Dictyostelium discoideum is a widely used model organism for studying the organization and function of the endocytic pathway. In Dictyostelium, the organization of the endocytic pathway is similar to that in higher eukaryotes. The pathway in Dictyostelium can be divided into four steps (see Fig. S1 in the supplemental material): uptake at the plasma membrane of particles and medium, transfer through early acidic endocytic compartments (lysosomes), passage into less acidic postlysosomes (PLs), and finally, exocytosis of undigested materials (17, 20). Thus, Dictyostelium recapitulates many of the functions of the endocytic pathway in mammalian cells, including some features observed in most cell types (lysosome biogenesis) and some observed only in specialized cells (phagocytosis, macropinocytosis, and lysosome secretion).Dictyostelium amoebae live in the soil, where they feed by ingesting and digesting other microorganisms. In addition, axenic laboratory strains can macropinocytose medium to ensure their growth. Accordingly, both in natural situations and in laboratory settings, the endocytic pathway plays a key role in the acquisition of nutrients by Dictyostelium cells. In agreement with this notion, several observations suggest that the physiology of the endocytic pathway is sensitive to nutrient availability. In particular, starvation induces secretion of lysosomal enzymes by an unknown mechanism (11). The morphology of the endocytic pathway is also sensitive to nutritional cues, as shown for example by the observation that formation of multilamellar endosomes is enhanced in cells fed with bacteria (18).Here, we analyzed the effect of starvation on the organization as well as the dynamics of the endocytic pathway. We found that, while the overall organization was not extensively modified in starved cells, the dynamics of endocytic compartments were altered. Moreover, analysis of two specific knockout mutants, the apm3 (6) and lvsB (8) strains, revealed that their phenotype was profoundly altered upon starvation, providing further insight about the role of Apm3 and LvsB in the endocytic pathway.     

In silico comparisons of lipid-related genes between Mycobacterium tuberculosis and BCG vaccine strains

Despite highly variable efficacy, BCG (Bacillus Calmette-Guérin) is the only vaccine available to prevent the tuberculosis (TB). Genomic heterogeneity between attenuated BCG strains and virulent Mycobacterium tuberculosis might help to explain this vaccine’s impaired capacity to induce long-term protection. Here, we investigate the lipid-related genes absent in attenuated BCG strains in order to correlate changes in both lipid metabolism and cell-wall lipid content to vaccine impairment. Whole genome sequences of M. tuberculosis H37Rv and the six most used BCG strains worldwide were aligned and the absent regions functionally categorized. Genomes of the BCG strains showed a total of 14 non-homologous lipid-related genes, including those belonging to mce3 operon, as well as the gene echaA1, which encodes an enoyl-CoA hydratase, and the genes encoding phospholipases PlcA, PlcB and PlcC. Taken together, the depletion of these M. tuberculosis H37Rv genomic regions were associated with marked alterations in lipid-related genes of BCG strains. Such alterations may indicate a dormant-like state and can be determining factors to the vaccine’s inability to induce long-term protection. These lipids can be further evaluated as an adjuvant to boost the current BCG-based vaccine.     

Size Does Matter: An Integrative In Vivo-In Silico Approach for the Treatment of Critical Size Bone Defects

Although bone has a unique restorative capacity, i.e., it has the potential to heal scarlessly, the conditions for spontaneous bone healing are not always present, leading to a delayed union or a non-union. In this work, we use an integrative in vivo - in silico approach to investigate the occurrence of non-unions, as well as to design possible treatment strategies thereof. The gap size of the domain geometry of a previously published mathematical model was enlarged in order to study the complex interplay of blood vessel formation, oxygen supply, growth factors and cell proliferation on the final healing outcome in large bone defects. The multiscale oxygen model was not only able to capture the essential aspects of in vivo non-unions, it also assisted in understanding the underlying mechanisms of action, i.e., the delayed vascularization of the central callus region resulted in harsh hypoxic conditions, cell death and finally disrupted bone healing. Inspired by the importance of a timely vascularization, as well as by the limited biological potential of the fracture hematoma, the influence of the host environment on the bone healing process in critical size defects was explored further. Moreover, dependent on the host environment, several treatment strategies were designed and tested for effectiveness. A qualitative correspondence between the predicted outcomes of certain treatment strategies and experimental observations was obtained, clearly illustrating the model''s potential. In conclusion, the results of this study demonstrate that due to the complex non-linear dynamics of blood vessel formation, oxygen supply, growth factor production and cell proliferation and the interactions thereof with the host environment, an integrative in silico-in vivo approach is a crucial tool to further unravel the occurrence and treatments of challenging critical sized bone defects.     

The Masquelet Technique for Membrane Induction and the Healing of Ovine Critical Sized Segmental Defects

The healing of critical sized segmental defects is an ongoing clinical problem. No method has achieved pre-eminence. The Masquelet technique is a relatively new innovation involving the induction of a fibrous tissue membrane around the bone defect site taking advantage of the body’s foreign body reaction to the presence of a polymethylmethacrylate (PMMA) spacer. The aim of this study was to investigate the properties and characteristics of this induced membrane and its effectiveness when used in conjunction with allograft or an allograft/autograft mix as filler materials in an ovine critical sized defect model. The resultant induced membrane was found to be effective in containing the graft materials in situ. It was demonstrated to be an organised pseudosynovial membrane which expressed bone morphogenic protein 2 (BMP2), transforming growth factor- beta (TGFβ), vascular endothelial growth factor (VEGF), von Willerbrand factor (vWF), interleukin 6 (IL-6) and interleukin 8 (IL-8). While more new bone growth was evident in the test groups compared to the controls animals at 12 weeks, the volumes were not statistically different and no defects were fully bridged. Of the two graft material groups, the allograft/autograft mix was shown to have a more rapid graft resorption rate than the allograft only group. While the Masquelet technique proved effective in producing a membrane to enclose graft materials, its ability to assist in the healing of critical sized segmental defects when compared to empty controls remained inconclusive.     

Caenorhabditis elegans: A Simple Nematode Infection Model for Penicillium marneffei

Penicillium marneffei, one of the most important thermal dimorphic fungi, is a severe threat to the life of immunocompromised patients. However, the pathogenic mechanisms of P. marneffei remain largely unknown. In this work, we developed a model host by using nematode Caenorhabditis elegans to investigate the virulence of P. marneffei. Using two P. marneffei clinical isolate strains 570 and 486, we revealed that in both liquid and solid media, the ingestion of live P. marneffei was lethal to C. elegans (P<0.001). Meanwhile, our results showed that the strain 570, which can produce red pigment, had stronger pathogenicity in C. elegans than the strain 486, which can’t produce red pigment (P<0.001). Microscopy showed the formation of red pigment and hyphae within C. elegans after incubation with P. marneffei for 4 h, which are supposed to be two contributors in nematodes killing. In addition, we used C. elegans as an in vivo model to evaluate different antifungal agents against P. marneffei, and found that antifungal agents including amphotericin B, terbinafine, fluconazole, itraconazole and voriconazole successfully prolonged the survival of nematodesinfected by P. marneffei. Overall, this alternative model host can provide us an easy tool to study the virulence of P. marneffei and screen antifungal agents.     

De novo GTP Biosynthesis Is Critical for Virulence of the Fungal Pathogen Cryptococcus neoformans

We have investigated the potential of the GTP synthesis pathways as chemotherapeutic targets in the human pathogen Cryptococcus neoformans, a common cause of fatal fungal meningoencephalitis. We find that de novo GTP biosynthesis, but not the alternate salvage pathway, is critical to cryptococcal dissemination and survival in vivo. Loss of inosine monophosphate dehydrogenase (IMPDH) in the de novo pathway results in slow growth and virulence factor defects, while loss of the cognate phosphoribosyltransferase in the salvage pathway yielded no phenotypes. Further, the Cryptococcus species complex displays variable sensitivity to the IMPDH inhibitor mycophenolic acid, and we uncover a rare drug-resistant subtype of C. gattii that suggests an adaptive response to microbial IMPDH inhibitors in its environmental niche. We report the structural and functional characterization of IMPDH from Cryptococcus, revealing insights into the basis for drug resistance and suggesting strategies for the development of fungal-specific inhibitors. The crystal structure reveals the position of the IMPDH moveable flap and catalytic arginine in the open conformation for the first time, plus unique, exploitable differences in the highly conserved active site. Treatment with mycophenolic acid led to significantly increased survival times in a nematode model, validating de novo GTP biosynthesis as an antifungal target in Cryptococcus.     

PCR Amplification-Independent Methods for Detection of Microbial Communities by the High-Density Microarray PhyloChip

Environmental microbial community analysis typically involves amplification by PCR, despite well-documented biases. We have developed two methods of PCR-independent microbial community analysis using the high-density microarray PhyloChip: direct hybridization of 16S rRNA (dirRNA) or rRNA converted to double-stranded cDNA (dscDNA). We compared dirRNA and dscDNA communities to PCR-amplified DNA communities using a mock community of eight taxa, as well as experiments derived from three environmental sample types: chromium-contaminated aquifer groundwater, tropical forest soil, and secondary sewage in seawater. Community profiles by both direct hybridization methods showed differences that were expected based on accompanying data but that were missing in PCR-amplified communities. Taxon richness decreased in RNA compared to that in DNA communities, suggesting a subset of 20% in soil and 60% in groundwater that is active; secondary sewage showed no difference between active and inactive populations. Direct hybridization of dscDNA and RNA is thus a viable alternative to PCR-amplified microbial community analysis, providing identification of the active populations within microbial communities that attenuate pollutants, drive global biogeochemical cycles, or proliferate disease states.     

In silico Analysis and Experimental Validation of Lignan Extracts from Kadsura longipedunculata for Potential 5-HT1AR Agonists

Objectives

Kadsura longipedunculata (KL) has been widely used for the treatment of insomnia in traditional Chinese medicine. The aim of this study was to explore the mechanism of the sedative and hypnotic effects of KL.

Materials and Methods

The content of KL was evaluated by HPLC-TOF-MS, and a potential target was found and used to construct its 3D structure to screen for potential ligands among the compounds in KL by using bioinformatics analysis, including similarity ensemble approach (SEA) docking, homology modeling, molecular docking and ligand-based pharmacophore. The PCPA-induced insomnia rat model was then applied to confirm the potential targets related to the sedative effects of KL by performing the forced swimming test (FST), the tail suspension test (TST) and the measurement of target-related proteins using western blotting and immunofluorescence.

Results

Bioinformatics analysis showed that most of lignan compounds in KL were optimal ligands for the 5-HT1A receptor (5-HT1AR), and they were found to be potential targets related to sedative effects; the main lignan content of KL extracts was characterized by HPLC-TOF-MS, with 7 proposed lignans detected. Administration of KL could significantly reduce FST and TST immobility time in the PCPA-induced 5HT-depleted insomnia rat model. The expressions of proteins related to the 5-HT1AR pathway were regulated by extracts of KL in a concentration-dependent manner, indicating that extracts of KL had 5-HT1AR agonist-like effects.

Conclusion

In silico analysis and experimental validation together demonstrated that lignan extracts from KL can target 5-HT1AR in insomniac rats, which could shed light on its use as a potential 5-HT1AR agonist drug.     

Effects of Prestretch on Neonatal Peripheral Nerve: An In Vitro Study

Background  Characterizing the biomechanical failure responses of neonatal peripheral nerves is critical in understanding stretch-related peripheral nerve injury mechanisms in neonates. Objective  This in vitro study investigated the effects of prestretch magnitude and duration on the biomechanical failure behavior of neonatal piglet brachial plexus (BP) and tibial nerves. Methods  BP and tibial nerves from 32 neonatal piglets were harvested and prestretched to 0, 10, or 20% strain for 90 or 300 seconds. These prestretched samples were then subjected to tensile loading until failure. Failure stress and strain were calculated from the obtained load-displacement data. Results  Prestretch magnitude significantly affected failure stress but not the failure strain. BP nerves prestretched to 10 or 20% strain, exhibiting significantly lower failure stress than those prestretched to 0% strain for both prestretch durations (90 and 300 seconds). Likewise, tibial nerves prestretched to 10 or 20% strain for 300 seconds, exhibiting significantly lower failure stress than the 0% prestretch group. An effect of prestretch duration on failure stress was also observed in the BP nerves when subjected to 20% prestretch strain such that the failure stress was significantly lower for 300 seconds group than 90 seconds group. No significant differences in the failure strains were observed. When comparing BP and tibial nerve failure responses, significantly higher failure stress was reported in tibial nerve prestretched to 20% strain for 300 seconds than BP nerve. Conclusion  These data suggest that neonatal peripheral nerves exhibit lower injury thresholds with increasing prestretch magnitude and duration while exhibiting regional differences.     

In Silico Investigations of the Anti-Catabolic Effects of Pamidronate and Denosumab on Multiple Myeloma-Induced Bone Disease

It is unclear whether the new anti-catabolic agent denosumab represents a viable alternative to the widely used anti-catabolic agent pamidronate in the treatment of Multiple Myeloma (MM)-induced bone disease. This lack of clarity primarily stems from the lack of sufficient clinical investigations, which are costly and time consuming. However, in silico investigations require less time and expense, suggesting that they may be a useful complement to traditional clinical investigations. In this paper, we aim to (i) develop integrated computational models that are suitable for investigating the effects of pamidronate and denosumab on MM-induced bone disease and (ii) evaluate the responses to pamidronate and denosumab treatments using these integrated models. To achieve these goals, pharmacokinetic models of pamidronate and denosumab are first developed and then calibrated and validated using different clinical datasets. Next, the integrated computational models are developed by incorporating the simulated transient concentrations of pamidronate and denosumab and simulations of their actions on the MM-bone compartment into the previously proposed MM-bone model. These integrated models are further calibrated and validated by different clinical datasets so that they are suitable to be applied to investigate the responses to the pamidronate and denosumab treatments. Finally, these responses are evaluated by quantifying the bone volume, bone turnover, and MM-cell density. This evaluation identifies four denosumab regimes that potentially produce an overall improved bone-related response compared with the recommended pamidronate regime. This in silico investigation supports the idea that denosumab represents an appropriate alternative to pamidronate in the treatment of MM-induced bone disease.     

Tyrosine 129 of the Murine Gammaherpesvirus M2 Protein Is Critical for M2 Function In Vivo

A common strategy shared by all known gammaherpesviruses is their ability to establish a latent infection in lymphocytes – predominantly in B cells. In immunocompromised patients, such as transplant recipients or AIDS patients, gammaherpesvirus infections can lead to the development of lymphoproliferative disease and lymphoid malignancies. The human gamma-herpesviruses, EBV and KSHV, encode proteins that are capable of modulating the host immune signaling machinery, thereby subverting host immune responses. Murine gamma-herpesvirus 68 (MHV68) infection of laboratory strains of mice has proven to be useful small-animal model that shares important pathogenic strategies with the human gamma-herpesviruses. The MHV68 M2 protein is known to manipulate B cell signaling and, dependent on route and dose of virus inoculation, plays a role in both the establishment of latency and virus reactivation. M2 contains two tyrosines that are targets for phosphorylation, and have been shown to interact with the B cell signaling machinery. Here we describe in vitro and in vivo studies of M2 mutants which reveals that while both tyrosines Y120 and Y129 are required for M2 induction of IL-10 expression from primary murine B cells in vitro, only Y129 is critical for reactivation from latency and plasma cell differentiation in vivo.