Supplementary Methods and Table from Protein Expression Profiling Identifies Subclasses of Breast Cancer and Predicts Prognosis

Supplementary Table 3 from Gene Expression Profiling for Molecular Characterization of Inflammatory Breast Cancer and Prediction of Response to Chemotherapy

Supplementary Table 2 from Gene Expression Profiling for Molecular Characterization of Inflammatory Breast Cancer and Prediction of Response to Chemotherapy

Supplementary Table 1 from Gene Expression Profiling for Molecular Characterization of Inflammatory Breast Cancer and Prediction of Response to Chemotherapy

The impact of tuberculosis and of anti-tuberculosis therapy on composition and modification of human lung microbiota has been the object of several investigations. However, no clear outcome has been presented so far and the relationship between M. tuberculosis pulmonary infection and the resident lung microbiota remains vague. In this work we describe the results obtained from a multicenter study of the microbiota of sputum samples from patients with tuberculosis or unrelated lung diseases and healthy donors recruited in Switzerland, Italy and Bangladesh, with the ultimate goal of discovering a microbiota-based biomarker associated with tuberculosis. Bacterial 16S rDNA amplification, high-throughput sequencing and extensive bioinformatic analyses revealed patient-specific flora and high variability in taxon abundance. No common signature could be identified among the individuals enrolled except for minor differences which were not consistent among the different geographical settings. Moreover, anti-tuberculosis therapy did not cause any important variation in microbiota diversity, thus precluding its exploitation as a biomarker for the follow up of tuberculosis patients undergoing treatment.

There is a growing need for unbiased clustering methods, ideally automated. We have developed a topology-based analysis tool called Two-Tier Mapper (TTMap) to detect subgroups in global gene expression datasets and identify their distinguishing features. First, TTMap discerns and adjusts for highly variable features in the control group and identifies outliers. Second, the deviation of each test sample from the control group in a high-dimensional space is computed and the test samples are clustered in a global and local network using a new topological algorithm based on Mapper. Validation of TTMap on both synthetic and biological datasets shows that it outperforms current clustering methods in sensitivity and stability; clustering is not affected by removal of samples from the control group, choice of normalization nor subselection of data. There is no user induced bias because all parameters are data-driven. Datasets can readily be combined into one analysis. TTMap reveals hitherto undetected gene expression changes in mouse mammary glands related to hormonal changes during the estrous cycle. This illustrates the ability to extract information from highly variable biological samples and its potential for personalized medicine.

Mouse liver regeneration after partial hepatectomy involves cells in the remaining tissue synchronously entering the cell division cycle. We have used this system and H3K4me3, Pol II and Pol III profiling to characterize adaptations in Pol III transcription. Our results broadly define a class of genes close to H3K4me3 and Pol II peaks, whose Pol III occupancy is high and stable, and another class, distant from Pol II peaks, whose Pol III occupancy strongly increases after partial hepatectomy. Pol III regulation in the liver thus entails both highly expressed housekeeping genes and genes whose expression can adapt to increased demand.

Proteins are a matter of dual nature. As a physical object, a protein molecule is a folded chain of amino acids with multifarious biochemistry. But it is also an instantiation along an evolutionary trajectory determined by the function performed by the protein within a hierarchy of interwoven interaction networks of the cell, the organism and the population. A physical theory of proteins therefore needs to unify both aspects, the biophysical and the evolutionary. Specifically, it should provide a model of how the DNA gene is mapped into the functional phenotype of the protein. We review several physical approaches to the protein problem, focusing on a mechanical framework which treats proteins as evolvable condensed matter: Mutations introduce localized perturbations in the gene, which are translated to localized perturbations in the protein matter. A natural tool to examine how mutations shape the phenotype are Green's functions. They map the evolutionary linkage among mutations in the gene (termed epistasis) to cooperative physical interactions among the amino acids in the protein. We discuss how the mechanistic view can be applied to examine basic questions of protein evolution and design.

How trafficking pathways and organelle abundance adapt in response to metabolic and physiological changes is still mysterious, although a few transcriptional regulators of organellar biogenesis have been identified in recent years. We previously found that the Wnt signaling directly controls lipid droplet formation, linking the cell storage capacity to the established functions of Wnt in development and differentiation. In the present paper, we report that Wnt-induced lipid droplet biogenesis does not depend on the canonical TCF/LEF transcription factors. Instead, we find that TFAP2 family members mediate the pro-lipid droplet signal induced by Wnt3a, leading to the notion that the TFAP2 transcription factor may function as a “master” regulator of lipid droplet biogenesis.

GETPrime (http://bbcftools.epfl.ch/getprime) is a database with a web frontend providing geneand transcript-specific, pre-computed qPCR primer pairs. The primers have been optimized for genome-wide specificity and for allowing the selective amplification of one or several splice variants of most known genes. To ease selection, primers have also been ranked according to defined criteria such as genome-wide specificity (with BLAST), amplicon size, and isoform coverage. Here, we report a major upgrade (2.0) of the database: eight new species (yeast, chicken, macaque, chimpanzee, rat, platypus, pufferfish, and Anolis carolinensis) now complement the five already included in the previous version (human, mouse, zebrafish, fly, and worm). Furthermore, the genomic reference has been updated to Ensembl v81 (while keeping earlier versions for backward compatibility) as a result of re-designing the back-end database and automating the import of relevant sections of the Ensembl database in species-independent fashion. This also allowed us to map known polymorphisms to the primers (on average three per primer for human), with the aim of reducing experimental error when targeting specific strains or individuals. Another consequence is that the inclusion of future Ensembl releases and other species has now become a relatively straightforward task.

Changes in cellular sterol species and concentrations can have profound effects on the transcriptional profile. In yeast, mutants defective in sterol biosynthesis show a wide range of changes in transcription, including a coinduction of anaerobic genes and ergosterol biosynthesis genes, biosynthesis of basic amino acids, and several stress genes. However the mechanisms underlying these changes are unknown. We identified mutations in the SAGA complex, a coactivator of transcription, which abrogate the ability to carry out most of these sterol-dependent transcriptional changes. In the erg3 mutant, the SAGA complex increases its occupancy time on many of the induced ergosterol and anaerobic gene promoters, increases its association with several relevant transcription factors and the SWI/SNF chromatin remodeling complex, and surprisingly, associates with an endocytic protein, Rvs167p, suggesting a moonlighting function for this protein in the sterol-regulated induction of the heat shock protein, HSP42 and HSP102, mRNAs.

Copy number variants (CNVs) are major contributors to genomic imbalance disorders. Phenotyping of 137 unrelated deletion and reciprocal duplication carriers of the distal 16p11.2 220 kb BP2-BP3 interval showed that these rearrangements are associated with autism spectrum disorders and mirror phenotypes of obesity/underweight and macrocephaly/microcephaly. Such phenotypes were previously associated with rearrangements of the non-overlapping proximal 16p11.2 600 kb BP4-BP5 interval. These two CNV-prone regions at 16p11.2 are reciprocally engaged in complex chromatin looping, as successfully confirmed by 4C-seq, fluorescence in situ hybridization and Hi-C, as well as coordinated expression and regulation of encompassed genes. We observed that genes differentially expressed in 16p11.2 BP4-BP5 CNV carriers are concomitantly modified in their chromatin interactions, suggesting that disruption of chromatin interplays could participate in the observed phenotypes. We also identified cis - and trans -acting chromatin contacts to other genomic regions previously associated with analogous phenotypes. For example, we uncovered that individuals with reciprocal rearrangements of the trans -contacted 2p15 locus similarly display mirror phenotypes on head circumference and weight. Our results indicate that chromosomal contacts’ maps could uncover functionally and clinically related genes.

BACKGROUND:Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors.METHODS:We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes.RESULTS:Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1 -/- mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1 -/- mouse model.CONCLUSIONS:These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

The SIB Swiss Institute of Bioinformatics (www.isb-sib.ch) provides world-class bioinformatics databases, software tools, services and training to the international life science community in academia and industry. These solutions allow life scientists to turn the exponentially growing amount of data into knowledge. Here, we provide an overview of SIB's resources and competence areas, with a strong focus on curated databases and SIB's most popular and widely used resources. In particular, SIB's Bioinformatics resource portal ExPASy features over 150 resources, including UniProtKB/Swiss-Prot, ENZYME, PROSITE, neXtProt, STRING, UniCarbKB, SugarBindDB, SwissRegulon, EPD, arrayMap, Bgee, SWISS-MODEL Repository, OMA, OrthoDB and other databases, which are briefly described in this article.

In human, mutations in bicaudal C1 (BICC1), an RNA binding protein, have been identified in patients with kidney dysplasia. Deletion of Bicc1 in mouse leads to left-right asymmetry randomization and renal cysts. Here, we show that BICC1 is also expressed in both the pancreatic progenitor cells that line the ducts during development, and in the ducts after birth, but not in differentiated endocrine or acinar cells. Genetic inactivation of Bicc1 leads to ductal cell over-proliferation and cyst formation. Transcriptome comparison between WT and Bicc1 KO pancreata, before the phenotype onset, reveals that PKD2 functions downstream of BICC1 in preventing cyst formation in the pancreas. Moreover, the analysis highlights immune cell infiltration and stromal reaction developing early in the pancreas of Bicc1 knockout mice. In addition to these functions in duct morphogenesis, BICC1 regulates NEUROG3+ endocrine progenitor production. Its deletion leads to a late but sustained endocrine progenitor decrease, resulting in a 50% reduction of endocrine cells. We show that BICC1 functions downstream of ONECUT1 in the pathway controlling both NEUROG3+ endocrine cell production and ductal morphogenesis, and suggest a new candidate gene for syndromes associating kidney dysplasia with pancreatic disorders, including diabetes.

A significant challenge of functional genomics is to develop methods for genome-scale acquisition and analysis of cell biological data. Here, we present an integrated method that combines genome-wide genetic perturbation of Saccharomyces cerevisiae with high-content screening to facilitate the genetic description of sub-cellular structures and compartment morphology. As proof of principle, we used a Rad52-GFP marker to examine DNA damage foci in similar to 20 million single cells from similar to 5,000 different mutant backgrounds in the context of selected genetic or chemical perturbations. Phenotypes were classified using a machine learning-based automated image analysis pipeline. 345 mutants were identified that had elevated numbers of DNA damage foci, almost half of which were identified only in sensitized backgrounds. Subsequent analysis of Vid22, a protein implicated in the DNA damage response, revealed that it acts together with the Sgs1 helicase at sites of DNA damage and preferentially binds G-quadruplex regions of the genome. This approach is extensible to numerous other cell biological markers and experimental systems.

Previous studies indicate that eukaryotic promoters display a stereotypical chromatin landscape characterized by a well-positioned +1 nucleosome near the transcription start site and an upstream -1 nucleosome that together demarcate a nucleosome-free (or -depleted) region. Here we present evidence that there are two distinct types of promoters distinguished by the resistance of the -1 nucleosome to micrococcal nuclease digestion. These different architectures are characterized by two sequence motifs that are broadly deployed at one set of promoters where a nuclease-sensitive ("fragile") nucleosome forms, but concentrated in a narrower, nucleosome-free region at all other promoters. The RSC nucleosome remodeler acts through the motifs to establish stable +1 and -1 nucleosome positions, while binding of a small set of general regulatory (pioneer) factors at fragile nucleosome promoters plays a key role in their destabilization. We propose that the fragile nucleosome promoter architecture is adapted for regulation of highly expressed, growth-related genes.

During the past 20 years, aerobic granular sludge (AGS) has been extensively studied with the aim to develop an attractive alternative to conventional activated sludge for wastewater treatment. The phosphate-accumulating organism (PAO) Candidatus ‘Accumulibacter phosphatis’ is often found with significant abundance in AGS as well as in other enhanced biological phosporus removal (EPBR) systems. Although they have never been isolated in pure culture, members of this bacterial genus appear to be genetically and physiologically more diversified than initially expected. Impaired EBPR performances observed in lab- and full-scale reactors have often been correlated to a decrease in Accumulibacter populations. This phenomenon has mainly been linked to the presence of bacterial competitors such as glycogen-accumulating organisms (GAO), and bacteriophages have only rarely been suspected to be responsible for this depletion (1). In the present study, the metagenome of 46 individual granules from a lab-scale AGS sequencing batch reactor was sequenced. The results showed a surprisingly variable relative abundance of Accumulibacter populations amongst the different granules that could only be partially explained by the “phenotype” of these granules. A co-occurrence analysis revealed a strong negative correlation between the number of Accumulibacter sequencing reads with the relative abundance of two bacteriophages, namely EBPR podovirus 1 (EPV1) and EBPR podovirus 3 (EPV3), that have been previously detected in a lab-scale EBPR reactor (2). These results suggest that these phages are the major reason for the variability of Accumulibacter relative abundance in the sampled granules which raises the question whether the Accumulibacter populations in the different granules have different sensitivities towards these phages. (1) Barr et al., 2010, Fems Microbiology Ecology, 631-642 (2) Skennerton et al., 2011, Plos One

Intestinal helminths are potent regulators of their host's immune system and can ameliorate inflammatory diseases such as allergic asthma. In the present study we have assessed whether this anti-inflammatory activity was purely intrinsic to helminths, or whether it also involved crosstalk with the local microbiota. We report that chronic infection with the murine helminth Heligmosomoides polygyrus bakeri (Hpb) altered the intestinal habitat, allowing increased short chain fatty acid (SCFA) production. Transfer of the Hpb-modified microbiota alone was sufficient to mediate protection against allergic asthma. The helminth-induced anti-inflammatory cytokine secretion and regulatory T cell suppressor activity that mediated the protection required the G protein-coupled receptor (GPR)-41. A similar alteration in the metabolic potential of intestinal bacterial communities was observed with diverse parasitic and host species, suggesting that this represents an evolutionary conserved mechanism of host-microbe-helminth interactions.