EpiDiP/NanoDiP: a versatile unsupervised machine learning edge computing platform for epigenomic tumour diagnostics.

DNA methylation analysis based on supervised machine learning algorithms with static reference data, allowing diagnostic tumour typing with unprecedented precision, has quickly become a new standard of care. Whereas genome-wide diagnostic methylation profiling is mostly performed on microarrays, an increasing number of institutions additionally employ nanopore sequencing as a faster alternative. In addition, methylation-specific parallel sequencing can generate methylation and genomic copy number data. Given these diverse approaches to methylation profiling, to date, there is no single tool that allows (1) classification and interpretation of microarray, nanopore and parallel sequencing data, (2) direct control of nanopore sequencers, and (3) the integration of microarray-based methylation reference data. Furthermore, no software capable of entirely running in routine diagnostic laboratory environments lacking high-performance computing and network infrastructure exists. To overcome these shortcomings, we present EpiDiP/NanoDiP as an open-source DNA methylation and copy number profiling suite, which has been benchmarked against an established supervised machine learning approach using in-house routine diagnostics data obtained between 2019 and 2021. Running locally on portable, cost- and energy-saving system-on-chip as well as gpGPU-augmented edge computing devices, NanoDiP works in offline mode, ensuring data privacy. It does not require the rigid training data annotation of supervised approaches. Furthermore, NanoDiP is the core of our public, free-of-charge EpiDiP web service which enables comparative methylation data analysis against an extensive reference data collection. We envision this versatile platform as a useful resource not only for neuropathologists and surgical pathologists but also for the tumour epigenetics research community. In daily diagnostic routine, analysis of native, unfixed biopsies by NanoDiP delivers molecular tumour classification in an intraoperative time frame.

Clinical, histopathological and molecular features of dedifferentiated melanomas: An EORTC Melanoma Group Retrospective Analysis.

PURPOSE: Dedifferentiated melanoma (DedM) poses significant diagnostic challenges. We aimed to investigate the clinical, histopathological and molecular features of DedM. Methylation signature (MS) and copy number profiling (CNP) were carried out in a subgroup of cases. PATIENTS AND METHODS: A retrospective series of 78 DedM tissue samples from 61 patients retrieved from EORTC (European Organisation for Research and Treatment of Cancer) Melanoma Group centres were centrally reviewed. Clinical and histopathological features were retrieved. In a subgroup of patients, genotyping through Infinium Methylation microarray and CNP analysis was carried out. RESULTS: Most patients (60/61) had a metastatic DedM showing most frequently an unclassified pleomorphic, spindle cell, or small round cell morphology akin to undifferentiated soft tissue sarcoma, rarely associated with heterologous elements. Overall, among 20 successfully analysed tissue samples from 16 patients, we found retained melanoma-like MS in only 7 tissue samples while a non-melanoma-like MS was observed in 13 tissue samples. In two patients from whom multiple specimens were analysed, some of the samples had a preserved cutaneous melanoma MS while other specimens exhibited an epigenetic shift towards a mesenchymal/sarcoma-like profile, matching the histological features. In these two patients, CNP was largely identical across all analysed specimens, in line with their common clonal origin, despite significant modification of their epigenome. CONCLUSIONS: Our study further highlights that DedM represents a real diagnostic challenge. While MS and genomic CNP may help pathologists to diagnose DedM, we provide proof-of-concept that dedifferentiation in melanoma is frequently associated with epigenetic modifications.

Global network reorganization during dynamic adaptations of Bacillus subtilis metabolism.

Adaptation of cells to environmental changes requires dynamic interactions between metabolic and regulatory networks, but studies typically address only one or a few layers of regulation. For nutritional shifts between two preferred carbon sources of Bacillus subtilis, we combined statistical and model-based data analyses of dynamic transcript, protein, and metabolite abundances and promoter activities. Adaptation to malate was rapid and primarily controlled posttranscriptionally compared with the slow, mainly transcriptionally controlled adaptation to glucose that entailed nearly half of the known transcription regulation network. Interactions across multiple levels of regulation were involved in adaptive changes that could also be achieved by controlling single genes. Our analysis suggests that global trade-offs and evolutionary constraints provide incentives to favor complex control programs.

Alterations of pre-mRNA splicing in human inflammatory bowel disease.

Alternative pre-mRNA splicing is regarded as a pivotal mechanism for generating proteome diversity and complexity from a limited inventory of mammalian genes. Aberrant splicing has been described as a predisposing factor for a number of diseases, but very little is known about its role in chronic inflammation. In this study, we systematically screened 149 splicing factors and 145 potential intron retention events for occurrence and differential expression in inflammatory bowel diseases (IBD). As a result, we identified 47 splicing factors and 33 intron retention events that were differentially regulated in mucosal tissue of IBD patients at transcript level. Despite the fact that Crohn's disease and ulcerative colitis, two subtypes of IBD, share the expression patterns of splicing factors and intron retention events in the majority of cases, we observed significant differences. To investigate these subtype-specific changes in detail we determined the expression levels of seven splicing factors (DUSP11, HNRPAB, HNRPH3, SLU7, SFR2IP, SFPQ, SF3B14) and three intron retention events (PARC, IER3, FGD2) in a cohort of 165 patients with inflammatory diseases of the colon (120 with IBD) and 30 healthy controls by real time PCR (TaqMan). This study demonstrates the potential impact of regulated splicing factors on subsequent regulated intron retention in the pathogenesis of chronic inflammation, exemplified by IBD.

Positioning of necrotic lobular intraepithelial neoplasias (LIN, grade 3) within the sequence of breast carcinoma progression.

Lobular intraepithelial neoplasia Grade 3 (LIN3) is a recently recognized variant of intraepithelial lobular neoplasia (LIN) of the breast composed of either uniform, generally small cells with massive lobular distension, pleomorphic cells, signet-ring cells, or any cell type with necrosis. In contrast to classic forms of LIN, there is no consensus on therapeutic strategies for LIN3. In part this is due to the paucity of molecular data that could assist in defining the relationship of LIN3 to classic LIN and carcinomas. In this study we have employed array comparative genomic hybridization to determine the patterns of chromosomal aberrations in nine LIN3 lesions. By comparison to array CGH data of 13 classic LIN lesions, we demonstrate that classic LIN and LIN3 share several recurrent changes, in particular gains of 1q and losses of 16q. Both aberrations are known to appear early in tumorigenesis and to be associated with good prognosis. However, apart from this overlap, there were a number of karyotypic features that were observed exclusively in LIN3. Clearly, this lesion was characterized by a significantly higher number of DNA copy number changes (9 vs. 31 on average), a considerable complexity of chromosomal rearrangements with more than 16 breakpoints in one chromosome and overlapping high copy amplifications encompassing a number of known oncogenes. Our data suggest that, at the genetic level, LIN3 represents a highly advanced lesion with considerable resemblance to carcinomas and, therefore, might represent the transition state from an intraepithelial neoplasm to breast carcinoma.

Tbc1d1 mutation in lean mouse strain confers leanness and protects from diet-induced obesity.

We previously identified Nob1 as a quantitative trait locus for high-fat diet-induced obesity and diabetes in genome-wide scans of outcross populations of obese and lean mouse strains. Additional crossbreeding experiments indicated that Nob1 represents an obesity suppressor from the lean Swiss Jim Lambert (SJL) strain. Here we identify a SJL-specific mutation in the Tbc1d1 gene that results in a truncated protein lacking the TBC Rab-GTPase-activating protein domain. TBC1D1, which has been recently linked to human obesity, is related to the insulin signaling protein AS160 and is predominantly expressed in skeletal muscle. Knockdown of TBC1D1 in skeletal muscle cells increased fatty acid uptake and oxidation, whereas overexpression of TBC1D1 had the opposite effect. Recombinant congenic mice lacking TBC1D1 showed reduced body weight, decreased respiratory quotient, increased fatty acid oxidation and reduced glucose uptake in isolated skeletal muscle. Our data strongly suggest that mutation of Tbc1d1 suppresses high-fat diet-induced obesity by increasing lipid use in skeletal muscle.

Small heat shock proteins prevent aggregation of citrate synthase and bind to the N-terminal region which is absent in thermostable forms of citrate synthase.

Citrate synthase (CS) is often used in chaperone assays since this thermosensitive enzyme aggregates at moderately increased temperatures. Small heat shock proteins (sHsps) are molecular chaperones specialized in preventing the aggregation of other proteins, termed substrate proteins, under conditions of transient heat stress. To investigate the mechanism whereby sHsps bind to and stabilize a substrate protein, we here used peptide array screening covering the sequence of porcine CS (P00889). Strong binding of sHsps was detected to a peptide corresponding to the most N-terminal alpha-helix in CS (amino acids Leu(13) to Gln(27)). The N-terminal alpha-helices in the CS dimer intertwine with the C-terminus in the other subunit and together form a stem-like structure which is protruding from the CS dimer. This stem-like structure is absent in thermostable forms of CS from thermophilic archaebacteria like Pyrococcus furiosus and Sulfolobus solfatacarium. These data therefore suggest that thermostabilization of thermosensitive CS by sHsps is achieved by stabilization of the C- and N-terminae in the protruding thermosensitive softspot, which is absent in thermostable forms of the CS dimer.

Array CGH identifies reciprocal 16p13.1 duplications and deletions that predispose to autism and/or mental retardation.

Autism and mental retardation (MR) are often associated, suggesting that these conditions are etiologically related. Recently, array-based comparative genomic hybridization (array CGH) has identified submicroscopic deletions and duplications as a common cause of MR, prompting us to search for such genomic imbalances in autism. Here we describe a 1.5-Mb duplication on chromosome 16p13.1 that was found by high-resolution array CGH in four severe autistic male patients from three unrelated families. The same duplication was identified in several variably affected and unaffected relatives. A deletion of the same interval was detected in three unrelated patients with MR and other clinical abnormalities. In one patient we revealed a further rearrangement of the 16p13 imbalance that was not present in his unaffected mother. Duplications and deletions of this 1.5-Mb interval have not been described as copy number variants in the Database of Genomic Variants and have not been identified in >600 individuals from other cohorts examined by high-resolution array CGH in our laboratory. Thus we conclude that these aberrations represent recurrent genomic imbalances which predispose to autism and/or MR.

Characterization of a 5.3 Mb deletion in 15q14 by comparative genomic hybridization using a whole genome "tiling path" BAC array in a girl with heart defect, cleft palate, and developmental delay.

High-resolution array CGH utilizing sets of overlapping BAC and PAC clones ("tiling path") covering the whole genome is a powerful novel tool for fast detection of submicroscopic chromosome deletions or duplications. We describe the successful application of a submegabase resolution whole genome "tiling path" BAC array to confirm and characterize a de novo interstitial deletion of chromosome 15. The deletion has a size of 5.3 Mb and is located within chromosome band 15q14, distal to the Prader-Willi/Angelman region. The affected girl had a heart defect, cleft palate, recurrent infections, and developmental delay. In contrast to GTG banding, array CGH determined the exact number of deleted genes and thus allowed the identification of candidate genes for cleft palate (GREM1, CX36, MEIS2), congenital heart defect (ACTC, GREM1, CX36, MEIS2), and mental retardation (ARHGAP11A, CHRNA7, CHRM5).

Differential binding studies applying functional protein microarrays and surface plasmon resonance.

A variety of different in vivo and in vitro technologies provide comprehensive insights in protein-protein interaction networks. Here we demonstrate a novel approach to analyze, verify and quantify putative interactions between two members of the S100 protein family and 80 recombinant proteins derived from a proteome-wide protein expression library. Surface plasmon resonance (SPR) using Biacore technology and functional protein microarrays were used as two independent methods to study protein-protein interactions. With this combined approach we were able to detect nine calcium-dependent interactions between Arg-Gly-Ser-(RGS)-His6 tagged proteins derived from the library and GST-tagged S100B and S100A6, respectively. For the protein microarray affinity-purified proteins from the expression library were spotted onto modified glass slides and probed with the S100 proteins. SPR experiments were performed in the same setup and in a vice-versa approach reversing analytes and ligands to determine distinct association and dissociation patterns of each positive interaction. Besides already known interaction partners, several novel binders were found independently with both detection methods, albeit analogous immobilization strategies had to be applied in both assays.

Age-specific hormonal decline is accompanied by transcriptional changes in human sebocytes in vitro.

The importance of hormones in endogenous aging has been displayed by recent studies performed on animal models and humans. To decipher the molecular mechanisms involved in aging we maintained human sebocytes at defined hormone-substituted conditions that corresponded to average serum levels of females from 20 (f20) to 60 (f60) years of age. The corresponding hormone receptor expression was demonstrated by reverse transcription-polymerase chain reaction (RT-PCR), Western blotting and immunocytochemistry. Cells at f60 produced significantly lower lipids than at f20. Increased mRNA and protein levels of c-Myc and increased protein levels of FN1, which have been associated with aging, were detected in SZ95 sebocytes at f60 compared to those detected at f20 after 5 days of treatment. Expression profiling employing a cDNA microarray composed of 15 529 cDNAs identified 899 genes with altered expression levels at f20 vs. f60. Confirmation of gene regulation was performed by real-time RT-PCR. The functional annotation of these genes according to the Gene Ontology identified pathways related to mitochondrial function, oxidative stress, ubiquitin-mediated proteolysis, cell cycle, immune responses, steroid biosynthesis and phospholipid degradation - all hallmarks of aging. Twenty-five genes in common with those identified in aging kidneys and several genes involved in neurodegenerative diseases were also detected. This is the first report describing the transcriptome of human sebocytes and its modification by a cocktail of hormones administered in age-specific levels and provides an in vitro model system, which approximates some of the hormone-dependent changes in gene transcription that occur during aging in humans.

Recent advances of protein microarrays.

Technological innovations and novel applications have greatly advanced the field of protein microarrays. Over the past two years, different types of protein microarrays have been used for serum profiling, protein abundance determinations, and identification of proteins that bind DNA or small compounds. However, considerable development is still required to ensure common quality standards and to establish large content repertoires. Here, we summarize applications available to date and discuss recent technological achievements and efforts on standardization.

Primary differentiation in the human blastocyst: comparative molecular portraits of inner cell mass and trophectoderm cells.

The primary differentiation event during mammalian development occurs at the blastocyst stage and leads to the delineation of the inner cell mass (ICM) and the trophectoderm (TE). We provide the first global mRNA expression data from immunosurgically dissected ICM cells, TE cells, and intact human blastocysts. Using a cDNA microarray composed of 15,529 cDNAs from known and novel genes, we identify marker transcripts specific to the ICM (e.g., OCT4/POU5F1, NANOG, HMGB1, and DPPA5) and TE (e.g., CDX2, ATP1B3, SFN, and IPL), in addition to novel ICM- and TE-specific expressed sequence tags. The expression patterns suggest that the emergence of pluripotent ICM and TE cell lineages from the morula is controlled by metabolic and signaling pathways, which include inter alia, WNT, mitogen-activated protein kinase, transforming growth factor-beta, NOTCH, integrin-mediated cell adhesion, phosphatidylinositol 3-kinase, and apoptosis. These data enhance our understanding of the first step in human cellular differentiation and, hence, the derivation of both embryonic stem cells and trophoblastic stem cells from these lineages.

High throughput identification of potential Arabidopsis mitogen-activated protein kinases substrates.

Mitogen-activated protein kinase (MAPK) cascades are universal and highly conserved signal transduction modules in eucaryotes, including plants. These protein phosphorylation cascades link extracellular stimuli to a wide range of cellular responses. However, the underlying mechanisms are so far unknown as information about phosphorylation substrates of plant MAPKs is lacking. In this study we addressed the challenging task of identifying potential substrates for Arabidopsis thaliana mitogen-activated protein kinases MPK3 and MPK6, which are activated by many environmental stress factors. For this purpose, we developed a novel protein microarray-based proteomic method allowing high throughput study of protein phosphorylation. We generated protein microarrays including 1,690 Arabidopsis proteins, which were obtained from the expression of an almost nonredundant uniclone set derived from an inflorescence meristem cDNA expression library. Microarrays were incubated with MAPKs in the presence of radioactive ATP. Using a threshold-based quantification method to evaluate the microarray results, we were able to identify 48 potential substrates of MPK3 and 39 of MPK6. 26 of them are common for both kinases. One of the identified MPK6 substrates, 1-aminocyclopropane-1-carboxylic acid synthase-6, was just recently shown as the first plant MAPK substrate in vivo, demonstrating the potential of our method to identify substrates with physiological relevance. Furthermore we revealed transcription factors, transcription regulators, splicing factors, receptors, histones, and others as candidate substrates indicating that regulation in response to MAPK signaling is very complex and not restricted to the transcriptional level. Nearly all of the 48 potential MPK3 substrates were confirmed by other in vitro methods. As a whole, our approach makes it possible to shortlist candidate substrates of mitogen-activated protein kinases as well as those of other protein kinases for further analysis. Follow-up in vivo experiments are essential to evaluate their physiological relevance.

Substrate-dependent regulation of anaerobic degradation pathways for toluene and ethylbenzene in a denitrifying bacterium, strain EbN1.

Anaerobic biodegradation of toluene and ethylbenzene is of environmental concern and biochemical interest due to toxicity and novel reactions, respectively. The denitrifying strain EbN1 is unique in anaerobically degrading both alkylbenzenes via different pathways which converge at benzoyl coenzyme A. The organization of genes involved in both pathways was only recently determined for strain EbN1. In the present study, global expression analysis (DNA microarray and proteomics) indicated involvement of several thus-far-unknown proteins in the degradation of both alkylbenzenes. For example, orf68 and orf57, framing the ebd operon, are implicated in ethylbenzene degradation, and the ebA1932 and ebA1936 genes, located 7.2 kb upstream of the bbs operon, are implicated in toluene degradation. In addition, expression studies were now possible on the level of the complete pathways. Growth experiments demonstrated that degradative capacities for toluene and ethylbenzene could be simultaneously induced, regardless of the substrate used for adaptation. Regulation was studied at the RNA (real-time reverse transcription-PCR and DNA microarray) and protein (two-dimensional-difference gel electrophoresis) level by using cells adapted to anaerobic growth with benzoate, toluene, ethylbenzene, or a mixture of toluene and ethylbenzene. Expression of the two toluene-related operons (bss and bbs) was specifically induced in toluene-adapted cells. In contrast, genes involved in anaerobic ethylbenzene degradation were induced in ethylbenzene- and toluene-adapted cells, suggesting that toluene may act as a gratuitous inducer. In agreement with the predicted sequential regulation of the ethylbenzene pathway, Ebd proteins (encoding subunits of ethylbenzene dehydrogenase) were formed in ethylbenzene- but not in acetophenone-adapted cells, while Apc proteins (subunits of predicted acetophenone carboxylase) were formed under both conditions.

Cross-species hybridisation of human and bovine orthologous genes on high density cDNA microarrays.

BACKGROUND: Cross-species gene-expression comparison is a powerful tool for the discovery of evolutionarily conserved mechanisms and pathways of expression control. The usefulness of cDNA microarrays in this context is that broad areas of homology are compared and hybridization probes are sufficiently large that small inter-species differences in nucleotide sequence would not affect the analytical results. This comparative genomics approach would allow a common set of genes within a specific developmental, metabolic, or disease-related gene pathway to be evaluated in experimental models of human diseases. The objective of this study was to investigate the feasibility and reproducibility of cross-species analysis employing a human cDNA microarray as probe. RESULTS: As a proof of principle, total RNA derived from human and bovine fetal brains was used as a source of labelled targets for hybridisation onto a human cDNA microarray composed of 349 characterised genes. Each gene was spotted 20 times representing 6,980 data points thus enabling highly reproducible spot quantification. Employing high stringency hybridisation and washing conditions, followed by data analysis, revealed slight differences in the expression levels and reproducibility of the signals between the two species. We also assigned each of the genes into three expression level categories- i.e. high, medium and low. The correlation co-efficient of cross hybridisation between the orthologous genes was 0.94. Verification of the array data by semi-quantitative RT-PCR using common primer sequences enabled co-amplification of both human and bovine transcripts. Finally, we were able to assign gene names to previously uncharacterised bovine ESTs. CONCLUSIONS: Results of our study demonstrate the harnessing and utilisation power of comparative genomics and prove the feasibility of using human microarrays to facilitate the identification of co-expressed orthologous genes in common tissues derived from different species.

Systematic delineation of a calmodulin peptide interaction.

We present a comprehensive profile of amino acid side-chain constraints in a calmodulin (CaM) peptide complex. These data were obtained from the analysis of calmodulin binding to an array of all single substitution analogues as well as N- and C-terminal truncations of the skMLCK derived M13 peptide ligand. The experimentally derived binding data were evaluated with respect to the known 3D-structure of the CaM/M13 complex. Besides an almost perfect agreement between the measured affinities and the structural data, the unexpected high-affine Asn5Ala variant of the M13(*) peptide described by Montigiani et al. could be verified. In contrast to other reports our data clearly support the postulate of the minor and major hydrophobic anchors of this calcium dependent interaction.

Multiplexed sorting of libraries on libraries: a novel method for empirical protein design by affinity-driven phage enrichment on synthetic peptide arrays.

Chemically synthesized peptide arrays on planar cellulose carriers are proposed as libraries of ligands suitable for the multiplexed simultaneous capture of peptide-specific acceptor proteins from a large randomly mutagenized library of acceptor proteins presented on bacteriophage M13 particles. This experimental set-up can be exploited to rapidly screen for individual new, distinct binding partners from two complementary libraries (two-dimensional screening). The technical feasibility of this empirical protein design approach was demonstrated with calmodulin as an aceptor protein using an array of mastoparan variants for multiplexed phage affinity enrichment.