Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Hypermutation of the inactive X chromosome is a frequent event in cancer.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

AI-driven discovery of blood xenobiotic biomarkers in neovascular age-related macular degeneration using iterative random forests.

PURPOSE: To investigate the xenobiotic profiles of patients with neovascular age-related macular degeneration (nAMD) undergoing anti-vascular endothelial growth factor (anti-VEGF) intravitreal therapy (IVT) to identify biomarkers indicative of clinical phenotypes through advanced AI methodologies. METHODS: In this cross-sectional observational study, we analyzed 156 peripheral blood xenobiotic features in a cohort of 46 nAMD patients stratified by choroidal neovascularization (CNV) control under anti-VEGF IVT. We employed Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) for measurement and leveraged an AI-driven iterative Random Forests (iRF) approach for robust pattern recognition and feature selection, aligning molecular profiles with clinical phenotypes. RESULTS: AI-augmented iRF models effectively refined the metabolite spectrum by discarding non-predictive elements. Perfluorooctanesulfonate (PFOS) and Ethyl ß-glucopyranoside were identified as significant biomarkers through this process, associated with various clinically relevant phenotypes. Unlike single metabolite classes, drug metabolites were distinctly correlated with subretinal fluid presence. CONCLUSIONS: This study underscores the enhanced capability of AI, particularly iRF, in dissecting complex metabolomic data to elucidate the xenobiotic landscape of nAMD and environmental impact on the disease. The preliminary biomarkers discovered offer promising directions for personalized treatment strategies, although further validation in broader cohorts is essential for clinical application.

A combined computational and functional approach identifies IGF2BP2 as a driver of chemoresistance in a wide array of pre-clinical models of colorectal cancer.

AIM: Chemoresistance is a major cause of treatment failure in colorectal cancer (CRC) therapy. In this study, the impact of the IGF2BP family of RNA-binding proteins on CRC chemoresistance was investigated using in silico, in vitro, and in vivo approaches. METHODS: Gene expression data from a well-characterized cohort and publicly available cross-linking immunoprecipitation sequencing (CLIP-Seq) data were collected. Resistance to chemotherapeutics was assessed in patient-derived xenografts (PDXs) and patient-derived organoids (PDOs). Functional studies were performed in 2D and 3D cell culture models, including proliferation, spheroid growth, and mitochondrial respiration analyses. RESULTS: We identified IGF2BP2 as the most abundant IGF2BP in primary and metastastatic CRC, correlating with tumor stage in patient samples and tumor growth in PDXs. IGF2BP2 expression in primary tumor tissue was significantly associated with resistance to selumetinib, gefitinib, and regorafenib in PDOs and to 5-fluorouracil and oxaliplatin in PDX in vivo. IGF2BP2 knockout (KO) HCT116 cells were more susceptible to regorafenib in 2D and to oxaliplatin, selumitinib, and nintedanib in 3D cell culture. Further, a bioinformatic analysis using CLIP data suggested stabilization of target transcripts in primary and metastatic tumors. Measurement of oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) revealed a decreased basal OCR and an increase in glycolytic ATP production rate in IGF2BP2 KO. In addition, real-time reverse transcriptase polymerase chain reaction (qPCR) analysis confirmed decreased expression of genes of the respiratory chain complex I, complex IV, and the outer mitochondrial membrane in IGF2BP2 KO cells. CONCLUSIONS: IGF2BP2 correlates with CRC tumor growth in vivo and promotes chemoresistance by altering mitochondrial respiratory chain metabolism. As a druggable target, IGF2BP2 could be used in future CRC therapy to overcome CRC chemoresistance.

Systemic Blood Proteome Patterns Reflect Disease Phenotypes in Neovascular Age-Related Macular Degeneration.

There is early evidence of extraocular systemic signals effecting function and morphology in neovascular age-related macular degeneration (nAMD). The prospective, cross-sectional BIOMAC study is an explorative investigation of peripheral blood proteome profiles and matched clinical features to uncover systemic determinacy in nAMD under anti-vascular endothelial growth factor intravitreal therapy (anti-VEGF IVT). It includes 46 nAMD patients stratified by the level of disease control under ongoing anti-VEGF treatment. Proteomic profiles in peripheral blood samples of every patient were detected with LC-MS/MS mass spectrometry. The patients underwent extensive clinical examination with a focus on macular function and morphology. In silico analysis includes unbiased dimensionality reduction and clustering, a subsequent annotation of clinical features, and non-linear models for recognition of underlying patterns. The model assessment was performed using leave-one-out cross validation. The findings provide an exploratory demonstration of the link between systemic proteomic signals and macular disease pattern using and validating non-linear classification models. Three main results were obtained: (1) Proteome-based clustering identifies two distinct patient subclusters with the smaller one (n = 10) exhibiting a strong signature for oxidative stress response. Matching the relevant meta-features on the individual patient's level identifies pulmonary dysfunction as an underlying health condition in these patients. (2) We identify biomarkers for nAMD disease features with Aldolase C as a putative factor associated with superior disease control under ongoing anti-VEGF treatment. (3) Apart from this, isolated protein markers are only weakly correlated with nAMD disease expression. In contrast, applying a non-linear classification model identifies complex molecular patterns hidden in a high number of proteomic dimensions determining macular disease expression. In conclusion, so far unconsidered systemic signals in the peripheral blood proteome contribute to the clinically observed phenotype of nAMD, which should be examined in future translational research on AMD.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

A framework for validating AI in precision medicine: considerations from the European ITFoC consortium.

BACKGROUND: Artificial intelligence (AI) has the potential to transform our healthcare systems significantly. New AI technologies based on machine learning approaches should play a key role in clinical decision-making in the future. However, their implementation in health care settings remains limited, mostly due to a lack of robust validation procedures. There is a need to develop reliable assessment frameworks for the clinical validation of AI. We present here an approach for assessing AI for predicting treatment response in triple-negative breast cancer (TNBC), using real-world data and molecular -omics data from clinical data warehouses and biobanks. METHODS: The European "ITFoC (Information Technology for the Future Of Cancer)" consortium designed a framework for the clinical validation of AI technologies for predicting treatment response in oncology. RESULTS: This framework is based on seven key steps specifying: (1) the intended use of AI, (2) the target population, (3) the timing of AI evaluation, (4) the datasets used for evaluation, (5) the procedures used for ensuring data safety (including data quality, privacy and security), (6) the metrics used for measuring performance, and (7) the procedures used to ensure that the AI is explainable. This framework forms the basis of a validation platform that we are building for the "ITFoC Challenge". This community-wide competition will make it possible to assess and compare AI algorithms for predicting the response to TNBC treatments with external real-world datasets. CONCLUSIONS: The predictive performance and safety of AI technologies must be assessed in a robust, unbiased and transparent manner before their implementation in healthcare settings. We believe that the consideration of the ITFoC consortium will contribute to the safe transfer and implementation of AI in clinical settings, in the context of precision oncology and personalized care.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris).

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

Decoding the regulatory landscape of medulloblastoma using DNA methylation sequencing.

Epigenetic alterations, that is, disruption of DNA methylation and chromatin architecture, are now acknowledged as a universal feature of tumorigenesis. Medulloblastoma, a clinically challenging, malignant childhood brain tumour, is no exception. Despite much progress from recent genomics studies, with recurrent changes identified in each of the four distinct tumour subgroups (WNT-pathway-activated, SHH-pathway-activated, and the less-well-characterized Group 3 and Group 4), many cases still lack an obvious genetic driver. Here we present whole-genome bisulphite-sequencing data from thirty-four human and five murine tumours plus eight human and three murine normal controls, augmented with matched whole-genome, RNA and chromatin immunoprecipitation sequencing data. This comprehensive data set allowed us to decipher several features underlying the interplay between the genome, epigenome and transcriptome, and its effects on medulloblastoma pathophysiology. Most notable were highly prevalent regions of hypomethylation correlating with increased gene expression, extending tens of kilobases downstream of transcription start sites. Focal regions of low methylation linked to transcription-factor-binding sites shed light on differential transcriptional networks between subgroups, whereas increased methylation due to re-normalization of repressed chromatin in DNA methylation valleys was positively correlated with gene expression. Large, partially methylated domains affecting up to one-third of the genome showed increased mutation rates and gene silencing in a subgroup-specific fashion. Epigenetic alterations also affected novel medulloblastoma candidate genes (for example, LIN28B), resulting in alternative promoter usage and/or differential messenger RNA/microRNA expression. Analysis of mouse medulloblastoma and precursor-cell methylation demonstrated a somatic origin for many alterations. Our data provide insights into the epigenetic regulation of transcription and genome organization in medulloblastoma pathogenesis, which are probably also of importance in a wider developmental and disease context.

Data and knowledge management in translational research: implementation of the eTRIKS platform for the IMI OncoTrack consortium.

BACKGROUND: For large international research consortia, such as those funded by the European Union's Horizon 2020 programme or the Innovative Medicines Initiative, good data coordination practices and tools are essential for the successful collection, organization and analysis of the resulting data. Research consortia are attempting ever more ambitious science to better understand disease, by leveraging technologies such as whole genome sequencing, proteomics, patient-derived biological models and computer-based systems biology simulations. RESULTS: The IMI eTRIKS consortium is charged with the task of developing an integrated knowledge management platform capable of supporting the complexity of the data generated by such research programmes. In this paper, using the example of the OncoTrack consortium, we describe a typical use case in translational medicine. The tranSMART knowledge management platform was implemented to support data from observational clinical cohorts, drug response data from cell culture models and drug response data from mouse xenograft tumour models. The high dimensional (omics) data from the molecular analyses of the corresponding biological materials were linked to these collections, so that users could browse and analyse these to derive candidate biomarkers. CONCLUSIONS: In all these steps, data mapping, linking and preparation are handled automatically by the tranSMART integration platform. Therefore, researchers without specialist data handling skills can focus directly on the scientific questions, without spending undue effort on processing the data and data integration, which are otherwise a burden and the most time-consuming part of translational research data analysis.

Low catalytic activity is insufficient to induce disease pathology in triosephosphate isomerase deficiency.

Triosephosphate isomerase (TPI) deficiency is a fatal genetic disorder characterized by hemolytic anemia and neurological dysfunction. Although the enzyme defect in TPI was discovered in the 1960s, the exact etiology of the disease is still debated. Some aspects indicate the disease could be caused by insufficient enzyme activity, whereas other observations indicate it could be a protein misfolding disease with tissue-specific differences in TPI activity. We generated a mouse model in which exchange of a conserved catalytic amino acid residue (isoleucine to valine, Ile170Val) reduces TPI specific activity without affecting the stability of the protein dimer. TPIIle170Val/Ile170Val mice exhibit an approximately 85% reduction in TPI activity consistently across all examined tissues, which is a stronger average, but more consistent, activity decline than observed in patients or symptomatic mouse models that carry structural defect mutant alleles. While monitoring protein expression levels revealed no evidence for protein instability, metabolite quantification indicated that glycolysis is affected by the active site mutation. TPIIle170Val/Ile170Val mice develop normally and show none of the disease symptoms associated with TPI deficiency. Therefore, without the stability defect that affects TPI activity in a tissue-specific manner, a strong decline in TPI catalytic activity is not sufficient to explain the pathological onset of TPI deficiency.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Synergistic cross-talk of hedgehog and interleukin-6 signaling drives growth of basal cell carcinoma.

Persistent activation of hedgehog (HH)/GLI signaling accounts for the development of basal cell carcinoma (BCC), a very frequent nonmelanoma skin cancer with rising incidence. Targeting HH/GLI signaling by approved pathway inhibitors can provide significant therapeutic benefit to BCC patients. However, limited response rates, development of drug resistance, and severe side effects of HH pathway inhibitors call for improved treatment strategies such as rational combination therapies simultaneously inhibiting HH/GLI and cooperative signals promoting the oncogenic activity of HH/GLI. In this study, we identified the interleukin-6 (IL6) pathway as a novel synergistic signal promoting oncogenic HH/GLI via STAT3 activation. Mechanistically, we provide evidence that signal integration of IL6 and HH/GLI occurs at the level of cis-regulatory sequences by co-binding of GLI and STAT3 to common HH-IL6 target gene promoters. Genetic inactivation of Il6 signaling in a mouse model of BCC significantly reduced in vivo tumor growth by interfering with HH/GLI-driven BCC proliferation. Our genetic and pharmacologic data suggest that combinatorial HH-IL6 pathway blockade is a promising approach to efficiently arrest cancer growth in BCC patients.

T cells become licensed in the lung to enter the central nervous system.

The blood­brain barrier (BBB) and the environment of the central nervous system (CNS) guard the nervous tissue from peripheral immune cells. In the autoimmune disease multiple sclerosis, myelin-reactive T-cell blasts are thought to transgress the BBB and create a pro-inflammatory environment in the CNS, thereby making possible a second autoimmune attack that starts from the leptomeningeal vessels and progresses into the parenchyma. Using a Lewis rat model of experimental autoimmune encephalomyelitis, we show here that contrary to the expectations of this concept, T-cell blasts do not efficiently enter the CNS and are not required to prepare the BBB for immune-cell recruitment. Instead, intravenously transferred T-cell blasts gain the capacity to enter the CNS after residing transiently within the lung tissues. Inside the lung tissues, they move along and within the airways to bronchus-associated lymphoid tissues and lung-draining mediastinal lymph nodes before they enter the blood circulation from where they reach the CNS. Effector T cells transferred directly into the airways showed a similar migratory pattern and retained their full pathogenicity. On their way the T cells fundamentally reprogrammed their gene-expression profile, characterized by downregulation of their activation program and upregulation of cellular locomotion molecules together with chemokine and adhesion receptors. The adhesion receptors include ninjurin 1, which participates in T-cell intravascular crawling on cerebral blood vessels. We detected that the lung constitutes a niche not only for activated T cells but also for resting myelin-reactive memory T cells. After local stimulation in the lung, these cells strongly proliferate and, after assuming migratory properties, enter the CNS and induce paralytic disease. The lung could therefore contribute to the activation of potentially autoaggressive T cells and their transition to a migratory mode as a prerequisite to entering their target tissues and inducing autoimmune disease.

Human Lineage-Specific Transcriptional Regulation through GA-Binding Protein Transcription Factor Alpha (GABPa).

A substantial fraction of phenotypic differences between closely related species are likely caused by differences in gene regulation. While this has already been postulated over 30 years ago, only few examples of evolutionary changes in gene regulation have been verified. Here, we identified and investigated binding sites of the transcription factor GA-binding protein alpha (GABPa) aiming to discover cis-regulatory adaptations on the human lineage. By performing chromatin immunoprecipitation-sequencing experiments in a human cell line, we found 11,619 putative GABPa binding sites. Through sequence comparisons of the human GABPa binding regions with orthologous sequences from 34 mammals, we identified substitutions that have resulted in 224 putative human-specific GABPa binding sites. To experimentally assess the transcriptional impact of those substitutions, we selected four promoters for promoter-reporter gene assays using human and African green monkey cells. We compared the activities of wild-type promoters to mutated forms, where we have introduced one or more substitutions to mimic the ancestral state devoid of the GABPa consensus binding sequence. Similarly, we introduced the human-specific substitutions into chimpanzee and macaque promoter backgrounds. Our results demonstrate that the identified substitutions are functional, both in human and nonhuman promoters. In addition, we performed GABPa knock-down experiments and found 1,215 genes as strong candidates for primary targets. Further analyses of our data sets link GABPa to cognitive disorders, diabetes, KRAB zinc finger (KRAB-ZNF), and human-specific genes. Thus, we propose that differences in GABPa binding sites played important roles in the evolution of human-specific phenotypes.

5-Fluorouracil affects assembly of stress granules based on RNA incorporation.

The antimetabolite 5-fluorouracil is a widely used chemotherapeutic for the treatment of several solid cancers. However, resistance to 5-fluorouracil remains a major drawback in its clinical use. In this study we report that treatment of HeLa cells with 5-fluorouracil resulted in de novo assembly of stress granules. Moreover, we revealed that stress granule assembly under stress conditions as well as disassembly is altered in cells treated with 5-fluorouracil. Notably, we discovered that RACK1, a protein mediating cell survival and apoptosis, is a component of 5-fluorouracil-induced stress granules. To explore the mode of action of 5-fluorouracil accountable for de novo stress granule assembly, we analyzed 5-fluorouracil metabolites and noticed that stress granule assembly is caused by RNA, not DNA incorporating 5-fluorouracil metabolites. Interestingly, we observed that other RNA incorporating drugs also cause assembly of stress granules. Thus, our results suggest that incorporation of chemotherapeutics into RNA may result in stress granule assembly with potential significance in chemoresistance.

ARH-seq: identification of differential splicing in RNA-seq data.

The computational prediction of alternative splicing from high-throughput sequencing data is inherently difficult and necessitates robust statistical measures because the differential splicing signal is overlaid by influencing factors such as gene expression differences and simultaneous expression of multiple isoforms amongst others. In this work we describe ARH-seq, a discovery tool for differential splicing in case-control studies that is based on the information-theoretic concept of entropy. ARH-seq works on high-throughput sequencing data and is an extension of the ARH method that was originally developed for exon microarrays. We show that the method has inherent features, such as independence of transcript exon number and independence of differential expression, what makes it particularly suited for detecting alternative splicing events from sequencing data. In order to test and validate our workflow we challenged it with publicly available sequencing data derived from human tissues and conducted a comparison with eight alternative computational methods. In order to judge the performance of the different methods we constructed a benchmark data set of true positive splicing events across different tissues agglomerated from public databases and show that ARH-seq is an accurate, computationally fast and high-performing method for detecting differential splicing events.

DNA-methylome analysis of mouse intestinal adenoma identifies a tumour-specific signature that is partly conserved in human colon cancer.

Aberrant CpG methylation is a universal epigenetic trait of cancer cell genomes. However, human cancer samples or cell lines preclude the investigation of epigenetic changes occurring early during tumour development. Here, we have used MeDIP-seq to analyse the DNA methylome of APC(Min) adenoma as a model for intestinal cancer initiation, and we present a list of more than 13,000 recurring differentially methylated regions (DMRs) characterizing intestinal adenoma of the mouse. We show that Polycomb Repressive Complex (PRC) targets are strongly enriched among hypermethylated DMRs, and several PRC2 components and DNA methyltransferases were up-regulated in adenoma. We further demonstrate by bisulfite pyrosequencing of purified cell populations that the DMR signature arises de novo in adenoma cells rather than by expansion of a pre-existing pattern in intestinal stem cells or undifferentiated crypt cells. We found that epigenetic silencing of tumour suppressors, which occurs frequently in colon cancer, was rare in adenoma. Quite strikingly, we identified a core set of DMRs, which is conserved between mouse adenoma and human colon cancer, thus possibly revealing a global panel of epigenetically modified genes for intestinal tumours. Our data allow a distinction between early conserved epigenetic alterations occurring in intestinal adenoma and late stochastic events promoting colon cancer progression, and may facilitate the selection of more specific clinical epigenetic biomarkers.

Combined ultra-low input mRNA and whole-genome sequencing of human embryonic stem cells.

BACKGROUND: Next Generation Sequencing has proven to be an exceptionally powerful tool in the field of genomics and transcriptomics. With recent development it is nowadays possible to analyze ultra-low input sample material down to single cells. Nevertheless, investigating such sample material often limits the analysis to either the genome or transcriptome. We describe here a combined analysis of both types of nucleic acids from the same sample material. METHODS: The method described enables the combined preparation of amplified cDNA as well as amplified whole-genome DNA from an ultra-low input sample material derived from a sub-colony of in-vitro cultivated human embryonic stem cells. cDNA is prepared by the application of oligo-dT coupled magnetic beads for mRNA capture, first strand synthesis and 3'-tailing followed by PCR. Whole-genome amplified DNA is prepared by Phi29 mediated amplification. Illumina sequencing is applied to short fragment libraries prepared from the amplified samples. RESULTS: We developed a protocol which enables the combined analysis of the genome as well as the transcriptome by Next Generation Sequencing from ultra-low input samples. The protocol was evaluated by sequencing sub-colony structures from human embryonic stem cells containing 150 to 200 cells. The method can be adapted to any available sequencing system. CONCLUSIONS: To our knowledge, this is the first report where sub-colonies of human embryonic stem cells have been analyzed both at the genomic as well as transcriptome level. The method of this proof of concept study may find useful practical applications for cases where only a limited number of cells are available, e.g. for tissues samples from biopsies, tumor spheres, circulating tumor cells and cells from early embryonic development. The results we present demonstrate that a combined analysis of genomic DNA and messenger RNA from ultra-low input samples is feasible and can readily be applied to other cellular systems with limited material available.

The direction of cross affects [corrected] obesity after puberty in male but not female offspring.

BACKGROUND: We investigated parent-of-origin and allele-specific expression effects on obesity and hepatic gene expression in reciprocal crosses between the Berlin Fat Mouse Inbred line (BFMI) and C57Bl/6NCrl (B6N). RESULTS: We found that F1-males with a BFMI mother developed 1.8 times more fat mass on a high fat diet at 10 weeks than F1-males of a BFMI father. The phenotype was detectable from six weeks on and was preserved after cross-fostering. RNA-seq data of liver provided evidence for higher biosynthesis and elongation of fatty acids (p = 0.00635) in obese male offspring of a BFMI mother versus lean offspring of a BFMI father. Furthermore, fatty acid degradation (p = 0.00198) and the peroxisome pathway were impaired (p = 0.00094). The circadian rhythm was affected as well (p = 0.00087). Among the highest up-regulated protein coding genes in obese males were Acot4 (1.82 fold, p = 0.022), Cyp4a10 (1.35 fold, p = 0.026) and Cyp4a14 (1.32 fold, p = 0.012), which hydroxylize fatty acids and which are known to be increased in liver steatosis. Obese males showed lower expression of the genetically imprinted and paternally expressed 3 (Peg3) gene (0.31 fold, p = 0.046) and higher expression of the androgen receptor (Ar) gene (2.38 fold, p = 0.068). Allelic imbalance was found for expression of ATP-binding cassette transporter gene Abca8b. Several of the differentially expressed genes contain estrogen response elements. CONCLUSIONS: Parent-of-origin effects during gametogenesis and/or fetal development in an obese mother epigenetically modify the transcription of genes that lead to enhanced fatty acid synthesis and impair ß-oxidation in the liver of male, but not female F1 offspring. Down-regulation of Peg3 could contribute to trigger this metabolic setting. At puberty, higher amounts of the androgen receptor and altered access to estrogen response elements in affected genes are likely responsible for male specific expression of genes that were epigenetically triggered. A suggestive lack of estrogen binding motifs was found for highly down-regulated genes in adult hepatocytes of obese F1 males (p = 0.074).