Calypso is an easy-to-use online software suite that allows non-expert users to mine, interpret and compare taxonomic information from metagenomic or 16S rDNA datasets. Calypso has a focus on multivariate statistical approaches that can identify complex environment-microbiome associations. The software enables quantitative visualizations, statistical testing, multivariate analysis, supervised learning, factor analysis, multivariable regression, network analysis and diversity estimates. Comprehensive help pages, tutorials and videos are provided via a wiki page.Availability and ImplementationThe web-interface is accessible via http://cgenome.net/calypso/. The software is programmed in Java, PERL and R and the source code is available from Zenodo (https://zenodo.org/record/50931). The software is freely available for non-commercial users.Supplementary information Supplementary data are available at Bioinformatics online.

Chikungunya virus (CHIKV) is an arthritogenic alphavirus causing epidemics of acute and chronic arthritic disease. Herein we describe a comprehensive RNA-Seq analysis of feet and lymph nodes at peak viraemia (day 2 post infection), acute arthritis (day 7) and chronic disease (day 30) in the CHIKV adult wild-type mouse model. Genes previously shown to be up-regulated in CHIKV patients were also up-regulated in the mouse model. CHIKV sequence information was also obtained with up to ≈8% of the reads mapping to the viral genome; however, no adaptive viral genome changes were apparent. Although day 2, 7 and 30 represent distinct stages of infection and disease, there was a pronounced overlap in up-regulated host genes and pathways. Type I interferon response genes (IRGs) represented up to ≈50% of up-regulated genes, even after loss of type I interferon induction on days 7 and 30. Bioinformatic analyses suggested a number of interferon response factors were primarily responsible for maintaining type I IRG induction. A group of genes prominent in the RNA-Seq analysis and hitherto unexplored in viral arthropathies were granzymes A, B and K. Granzyme A-/- and to a lesser extent granzyme K-/-, but not granzyme B-/-, mice showed a pronounced reduction in foot swelling and arthritis, with analysis of granzyme A-/- mice showing no reductions in viral loads but reduced NK and T cell infiltrates post CHIKV infection. Treatment with Serpinb6b, a granzyme A inhibitor, also reduced arthritic inflammation in wild-type mice. In non-human primates circulating granzyme A levels were elevated after CHIKV infection, with the increase correlating with viral load. Elevated granzyme A levels were also seen in a small cohort of human CHIKV patients. Taken together these results suggest granzyme A is an important driver of arthritic inflammation and a potential target for therapy.Trial Registration: ClinicalTrials.gov NCT00281294

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.

We sequenced 8 melanoma exomes to identify novel somatic mutations in metastatic melanoma. Focusing on the MAP3K family, we found that 24% of melanoma cell lines have mutations in the protein-coding regions of either MAP3K5 or MAP3K9. Structural modelling predicts that mutations in the kinase domain may affect the activity and regulation of MAP3K5/9 protein kinases. The position of the mutations and loss of heterozygosity of MAP3K5 and MAP3K9 in 85% and 67% of melanoma samples, respectively, together suggest that the mutations are likely inactivating. In vitro kinase assay shows reduction in kinase activity in MAP3K5 I780F and MAP3K9 W333X mutants. Overexpression of MAP3K5 or MAP3K9 mutant in HEK293T cells reduces phosphorylation of downstream MAP kinases. Attenuation of MAP3K9 function in melanoma cells using siRNA leads to increased cell viability after temozolomide treatment, suggesting that decreased MAP3K pathway activity can lead to chemoresistance in melanoma.

BackgroundSmchd1 is an epigenetic modifier essential for X chromosome inactivation: female embryos lacking Smchd1 fail during midgestational development. Male mice are less affected by Smchd1-loss, with some (but not all) surviving to become fertile adults on the FVB/n genetic background. On other genetic backgrounds, all males lacking Smchd1 die perinatally. This suggests that, in addition to being critical for X inactivation, Smchd1 functions to control the expression of essential autosomal genes.ResultsUsing genome-wide microarray expression profiling and RNA-seq, we have identified additional genes that fail X inactivation in female Smchd1 mutants and have identified autosomal genes in male mice where the normal expression pattern depends upon Smchd1. A subset of genes in the Snrpn imprinted gene cluster show an epigenetic signature and biallelic expression consistent with loss of imprinting in the absence of Smchd1. In addition, single nucleotide polymorphism analysis of expressed genes in the placenta shows that the Igf2r imprinted gene cluster is also disrupted, with Slc22a3 showing biallelic expression in the absence of Smchd1. In both cases, the disruption was not due to loss of the differential methylation that marks the imprint control region, but affected genes remote from this primary imprint controlling element. The clustered protocadherins (Pcdhα, Pcdhβ, and Pcdhγ) also show altered expression levels, suggesting that their unique pattern of random combinatorial monoallelic expression might also be disrupted.ConclusionsSmchd1 has a role in the expression of several autosomal gene clusters that are subject to monoallelic expression, rather than being restricted to functioning uniquely in X inactivation. Our findings, combined with the recent report implicating heterozygous mutations of SMCHD1 as a causal factor in the digenically inherited muscular weakness syndrome facioscapulohumeral muscular dystrophy-2, highlight the potential importance of Smchd1 in the etiology of diverse human diseases.