High Yap and Mll1 promote a persistent regenerative cell state induced by Notch signaling and loss of p53.

Specified intestinal epithelial cells reprogram and contribute to the regeneration and renewal of the epithelium upon injury. Mutations that deregulate such renewal processes may contribute to tumorigenesis. Using intestinal organoids, we show that concomitant activation of Notch signaling and ablation of p53 induce a highly proliferative and regenerative cell state, which is associated with increased levels of Yap and the histone methyltransferase Mll1. The induced signaling system orchestrates high proliferation, self-renewal, and niche-factor-independent growth, and elevates the trimethylation of histone 3 at lysine 4 (H3K4me3). We demonstrate that Yap and Mll1 are also elevated in patient-derived colorectal cancer (CRC) organoids and control growth and viability. Our data suggest that Notch activation and p53 ablation induce a signaling circuitry involving Yap and the epigenetic regulator Mll1, which locks cells in a proliferative and regenerative state that renders them susceptible for tumorigenesis.

TAp73 is a central transcriptional regulator of airway multiciliogenesis.

Motile multiciliated cells (MCCs) have critical roles in respiratory health and disease and are essential for cleaning inhaled pollutants and pathogens from airways. Despite their significance for human disease, the transcriptional control that governs multiciliogenesis remains poorly understood. Here we identify TP73, a p53 homolog, as governing the program for airway multiciliogenesis. Mice with TP73 deficiency suffer from chronic respiratory tract infections due to profound defects in ciliogenesis and complete loss of mucociliary clearance. Organotypic airway cultures pinpoint TAp73 as necessary and sufficient for basal body docking, axonemal extension, and motility during the differentiation of MCC progenitors. Mechanistically, cross-species genomic analyses and complete ciliary rescue of knockout MCCs identify TAp73 as the conserved central transcriptional integrator of multiciliogenesis. TAp73 directly activates the key regulators FoxJ1, Rfx2, Rfx3, and miR34bc plus nearly 50 structural and functional ciliary genes, some of which are associated with human ciliopathies. Our results position TAp73 as a novel central regulator of MCC differentiation.

Single-nucleus RNA sequencing of plant tissues using a nanowell-based system.

Single-cell genomics provides unprecedented potential for research on plant development and environmental responses. Here, we introduce a generic procedure for plant nucleus isolation combined with nanowell-based library preparation. Our method enables the transcriptome analysis of thousands of individual plant nuclei. It serves as an alternative to the use of protoplast isolation, which is currently a standard methodology for plant single-cell genomics, although it can be challenging for some plant tissues. We show the applicability of our nucleus isolation method by using different plant materials from different species. The potential of our single-nucleus RNA sequencing method is shown through the characterization of transcriptomes of seedlings and developing flowers from Arabidopsis thaliana. We evaluated the transcriptome dynamics during the early stages of anther development, identified stage-specific activities of transcription factors regulating this process, and predicted potential target genes of these transcription factors. Our nucleus isolation procedure can be applied in different plant species and tissues, thus expanding the toolkit for plant single-cell genomics experiments.

Regional gene expression patterns are associated with task-specific brain activation during reward and emotion processing measured with functional MRI.

The exploration of the spatial relationship between gene expression profiles and task-evoked response patterns known to be altered in neuropsychiatric disorders, for example depression, can guide the development of more targeted therapies. Here, we estimated the correlation between human transcriptome data and two different brain activation maps measured with functional magnetic resonance imaging (fMRI) in healthy subjects. Whole-brain activation patterns evoked during an emotional face recognition task were associated with topological mRNA expression of genes involved in cellular transport. In contrast, fMRI activation patterns related to the acceptance of monetary rewards were associated with genes implicated in cellular localization processes, metabolism, translation, and synapse regulation. An overlap of these genes with risk genes from major depressive disorder genome-wide association studies revealed the involvement of the master regulators TCF4 and PAX6 in emotion and reward processing. Overall, the identification of stable relationships between spatial gene expression profiles and fMRI data may reshape the prospects for imaging transcriptomics studies.

Dissecting Herpes Simplex Virus 1-Induced Host Shutoff at the RNA Level.

Herpes simplex virus 1 (HSV-1) induces a profound host shutoff during lytic infection. The virion host shutoff (vhs) protein plays a key role in this process by efficiently cleaving host and viral mRNAs. Furthermore, the onset of viral DNA replication is accompanied by a rapid decline in host transcriptional activity. To dissect relative contributions of both mechanisms and elucidate gene-specific host transcriptional responses throughout the first 8 h of lytic HSV-1 infection, we used transcriptome sequencing of total, newly transcribed (4sU-labeled) and chromatin-associated RNA in wild-type (WT) and Δvhs mutant infection of primary human fibroblasts. Following virus entry, vhs activity rapidly plateaued at an elimination rate of around 30% of cellular mRNAs per hour until 8 h postinfection (p.i.). In parallel, host transcriptional activity dropped to 10 to 20%. While the combined effects of both phenomena dominated infection-induced changes in total RNA, extensive gene-specific transcriptional regulation was observable in chromatin-associated RNA and was surprisingly concordant between WT and Δvhs infections. Both induced strong transcriptional upregulation of a small subset of genes that were poorly expressed prior to infection but already primed by H3K4me3 histone marks at their promoters. Most interestingly, analysis of chromatin-associated RNA revealed vhs-nuclease-activity-dependent transcriptional downregulation of at least 150 cellular genes, in particular of many integrin adhesome and extracellular matrix components. This was accompanied by a vhs-dependent reduction in protein levels by 8 h p.i. for many of these genes. In summary, our study provides a comprehensive picture of the molecular mechanisms that govern cellular RNA metabolism during the first 8 h of lytic HSV-1 infection.IMPORTANCE The HSV-1 virion host shutoff (vhs) protein efficiently cleaves both host and viral mRNAs in a translation-dependent manner. In this study, we model and quantify changes in vhs activity, as well as virus-induced global loss of host transcriptional activity, during productive HSV-1 infection. In general, HSV-1-induced alterations in total RNA levels were dominated by these two global effects. In contrast, chromatin-associated RNA depicted gene-specific transcriptional changes. This revealed highly concordant transcriptional changes in WT and Δvhs infections, confirmed DUX4 as a key transcriptional regulator in HSV-1 infection, and identified vhs-dependent transcriptional downregulation of the integrin adhesome and extracellular matrix components. The latter explained seemingly gene-specific effects previously attributed to vhs-mediated mRNA degradation and resulted in a concordant loss in protein levels by 8 h p.i. for many of the respective genes.

The Transcription Factor EB (TFEB) Sensitizes the Heart to Chronic Pressure Overload.

The transcription factor EB (TFEB) promotes protein degradation by the autophagy and lysosomal pathway (ALP) and overexpression of TFEB was suggested for the treatment of ALP-related diseases that often affect the heart. However, TFEB-mediated ALP induction may perturb cardiac stress response. We used adeno-associated viral vectors type 9 (AAV9) to overexpress TFEB (AAV9-Tfeb) or Luciferase-control (AAV9-Luc) in cardiomyocytes of 12-week-old male mice. Mice were subjected to transverse aortic constriction (TAC, 27G; AAV9-Luc: n = 9; AAV9-Tfeb: n = 14) or sham (AAV9-Luc: n = 9; AAV9-Tfeb: n = 9) surgery for 28 days. Heart morphology, echocardiography, gene expression, and protein levels were monitored. AAV9-Tfeb had no effect on cardiac structure and function in sham animals. TAC resulted in compensated left ventricular hypertrophy in AAV9-Luc mice. AAV9-Tfeb TAC mice showed a reduced LV ejection fraction and increased left ventricular diameters. Morphological, histological, and real-time PCR analyses showed increased heart weights, exaggerated fibrosis, and higher expression of stress markers and remodeling genes in AAV9-Tfeb TAC compared to AAV9-Luc TAC. RNA-sequencing, real-time PCR and Western Blot revealed a stronger ALP activation in the hearts of AAV9-Tfeb TAC mice. Cardiomyocyte-specific TFEB-overexpression promoted ALP gene expression during TAC, which was associated with heart failure. Treatment of ALP-related diseases by overexpression of TFEB warrants careful consideration.

Epigenetic alterations in longevity regulators, reduced life span, and exacerbated aging-related pathology in old father offspring mice.

Advanced age is not only a major risk factor for a range of disorders within an aging individual but may also enhance susceptibility for disease in the next generation. In humans, advanced paternal age has been associated with increased risk for a number of diseases. Experiments in rodent models have provided initial evidence that paternal age can influence behavioral traits in offspring animals, but the overall scope and extent of paternal age effects on health and disease across the life span remain underexplored. Here, we report that old father offspring mice showed a reduced life span and an exacerbated development of aging traits compared with young father offspring mice. Genome-wide epigenetic analyses of sperm from aging males and old father offspring tissue identified differentially methylated promoters, enriched for genes involved in the regulation of evolutionarily conserved longevity pathways. Gene expression analyses, biochemical experiments, and functional studies revealed evidence for an overactive mTORC1 signaling pathway in old father offspring mice. Pharmacological mTOR inhibition during the course of normal aging ameliorated many of the aging traits that were exacerbated in old father offspring mice. These findings raise the possibility that inherited alterations in longevity pathways contribute to intergenerational effects of aging in old father offspring mice.

CRISPLD1: a novel conserved target in the transition to human heart failure.

Heart failure is a major health problem worldwide with a significant morbidity and mortality rate. Although studied extensively in animal models, data from patients at the compensated disease stage are lacking. We sampled myocardium biopsies from aortic stenosis patients with compensated hypertrophy and moderate heart failure and used transcriptomics to study the transition to failure. Sequencing and comparative analysis of analogous samples of mice with transverse aortic constriction identified 25 candidate genes with similar regulation in response to pressure overload, reflecting highly conserved molecular processes. The gene cysteine-rich secretory protein LCCL domain containing 1 (CRISPLD1) is upregulated in the transition to failure in human and mouse and its function is unknown. Homology to ion channel regulatory toxins suggests a role in Ca2+ cycling. CRISPR/Cas9-mediated loss-of-function leads to dysregulated Ca2+ handling in human-induced pluripotent stem cell-derived cardiomyocytes. The downregulation of prohypertrophic, proapoptotic and Ca2+-signaling pathways upon CRISPLD1-KO and its upregulation in the transition to failure implicates a contribution to adverse remodeling. These findings provide new pathophysiological data on Ca2+ regulation in the transition to failure and novel candidate genes with promising potential for therapeutic interventions.

Organ-specific small non-coding RNA responses in domestic (Sudani) ducks experimentally infected with highly pathogenic avian influenza virus (H5N1).

The duck represents an important reservoir of influenza viruses for transmission to other avian and mammalian hosts, including humans. The increased pathogenicity of the recently emerging clades of highly pathogenic avian influenza (HPAI) viruses of the H5N1 subtype in ducks features systemic viral spread and organ-to-organ variation in viral transcription and tissue damage. We previously reported that experimental infection of Sudani ducks (Cairina moschata) with an Egyptian HPAI (H5N1) virus (clade 2.2.1.2) features high viral replication and severe tissue damage in lung, but lower viral replication and only mild histological changes in brain. Little is known about the involvement of miRNA in organ-specific responses to H5N1 viruses in ducks, and involvement of the other classes of small noncoding RNA (sncRNA) has not been investigated so far. Following RNA sequencing, we have annotated the duck sncRNome and compared global expression changes of the four major sncRNA classes (miRNAs, piRNAs, snoRNAs, snRNAs) between duck lung and brain during a 120 h time course of infection with this HPAI strain. We find major organ-specific differences in miRNA, piRNA and snoRNA populations even before infection and substantial reprogramming of all sncRNA classes throughout infection, which was less pronounced in brain. Pathway prediction analysis of miRNA targets revealed enrichment of inflammation-, infection- and apoptosis-related pathways in lung, but enrichment of metabolism-related pathways (including tryptophan metabolism) in brain. Thus, organ-specific differences in sncRNA responses may contribute to differences in viral replication and organ damage in ducks infected with isolates from this emerging HPAI clade, and likely other strains.

The landscape of human mutually exclusive splicing.

Mutually exclusive splicing of exons is a mechanism of functional gene and protein diversification with pivotal roles in organismal development and diseases such as Timothy syndrome, cardiomyopathy and cancer in humans. In order to obtain a first genomewide estimate of the extent and biological role of mutually exclusive splicing in humans, we predicted and subsequently validated mutually exclusive exons (MXEs) using 515 publically available RNA-Seq datasets. Here, we provide evidence for the expression of over 855 MXEs, 42% of which represent novel exons, increasing the annotated human mutually exclusive exome more than fivefold. The data provide strong evidence for the existence of large and multi-cluster MXEs in higher vertebrates and offer new insights into MXE evolution. More than 82% of the MXE clusters are conserved in mammals, and five clusters have homologous clusters in Drosophila Finally, MXEs are significantly enriched in pathogenic mutations and their spatio-temporal expression might predict human disease pathology.

Fructooligosaccharide intake promotes epigenetic changes in the intestinal mucosa in growing and ageing rats.

PURPOSE: The aim of this study was to investigate the relationship between fructooligosaccharide (FOS) intake at different life stages of Wistar rats and its stimulatory effects on intestinal parameters. METHODS: Recently weaned and ageing female rats were divided into growing and ageing treatments, which were fed diets that partially replaced sucrose with FOS for 12 weeks. RESULTS: Dietary FOS intake induced a significant increase in the numbers of Bifidobacterium and Lactobacillus in growing rats. FOS intake was associated with increased butyric acid levels and a reduced pH of the caecal contents at both ages. Differential gene expression patterns were observed by microarray analysis of growing and ageing animals fed the FOS diet. A total of 133 genes showed detectable changes in expression in the growing rats, while there were only 19 gene expression changes in ageing rats fed with FOS. CONCLUSION: These results suggest that dietary FOS intake may be beneficial for some parameters of intestinal health in growing rats.

Severe DCM phenotype of patient harboring RBM20 mutation S635A can be modeled by patient-specific induced pluripotent stem cell-derived cardiomyocytes.

The ability to generate patient-specific induced pluripotent stem cells (iPSCs) provides a unique opportunity for modeling heart disease in vitro. In this study, we generated iPSCs from a patient with dilated cardiomyopathy (DCM) caused by a missense mutation S635A in RNA-binding motif protein 20 (RBM20) and investigated the functionality and cell biology of cardiomyocytes (CMs) derived from patient-specific iPSCs (RBM20-iPSCs). The RBM20-iPSC-CMs showed abnormal distribution of sarcomeric α-actinin and defective calcium handling compared to control-iPSC-CMs, suggesting disorganized myofilament structure and altered calcium machinery in CMs of the RBM20 patient. Engineered heart muscles (EHMs) from RBM20-iPSC-CMs showed that not only active force generation was impaired in RBM20-EHMs but also passive stress of the tissue was decreased, suggesting a higher visco-elasticity of RBM20-EHMs. Furthermore, we observed a reduced titin (TTN) N2B-isoform expression in RBM20-iPSC-CMs by demonstrating a reduction of exon skipping in the PEVK region of TTN and an inhibition of TTN isoform switch. In contrast, in control-iPSC-CMs both TTN isoforms N2B and N2BA were expressed, indicating that the TTN isoform switch occurs already during early cardiogenesis. Using next generation RNA sequencing, we mapped transcriptome and splicing target profiles of RBM20-iPSC-CMs and identified different cardiac gene networks in response to the analyzed RBM20 mutation in cardiac-specific processes. These findings shed the first light on molecular mechanisms of RBM20-dependent pathological cardiac remodeling leading to DCM. Our data demonstrate that iPSC-CMs coupled with EHMs provide a powerful tool for evaluating disease-relevant functional defects and for a deeper mechanistic understanding of alternative splicing-related cardiac diseases.

SNooPer: a machine learning-based method for somatic variant identification from low-pass next-generation sequencing.

BACKGROUND: Next-generation sequencing (NGS) allows unbiased, in-depth interrogation of cancer genomes. Many somatic variant callers have been developed yet accurate ascertainment of somatic variants remains a considerable challenge as evidenced by the varying mutation call rates and low concordance among callers. Statistical model-based algorithms that are currently available perform well under ideal scenarios, such as high sequencing depth, homogeneous tumor samples, high somatic variant allele frequency (VAF), but show limited performance with sub-optimal data such as low-pass whole-exome/genome sequencing data. While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. RESULTS: For these reasons, we developed SNooPer, a versatile machine learning approach that uses Random Forest classification models to accurately call somatic variants in low-depth sequencing data. SNooPer uses a subset of variant positions from the sequencing output for which the class, true variation or sequencing error, is known to train the data-specific model. Here, using a real dataset of 40 childhood acute lymphoblastic leukemia patients, we show how the SNooPer algorithm is not affected by low coverage or low VAFs, and can be used to reduce overall sequencing costs while maintaining high specificity and sensitivity to somatic variant calling. When compared to three benchmarked somatic callers, SNooPer demonstrated the best overall performance. CONCLUSIONS: While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. The flexibility of SNooPer's random forest protects against technical bias and systematic errors, and is appealing in that it does not rely on user-defined parameters. The code and user guide can be downloaded at https://sourceforge.net/projects/snooper/ .

Identification of candidate genes for drought tolerance in coffee by high-throughput sequencing in the shoot apex of different Coffea arabica cultivars.

BACKGROUND: Drought is a widespread limiting factor in coffee plants. It affects plant development, fruit production, bean development and consequently beverage quality. Genetic diversity for drought tolerance exists within the coffee genus. However, the molecular mechanisms underlying the adaptation of coffee plants to drought are largely unknown. In this study, we compared the molecular responses to drought in two commercial cultivars (IAPAR59, drought-tolerant and Rubi, drought-susceptible) of Coffea arabica grown in the field under control (irrigation) and drought conditions using the pyrosequencing of RNA extracted from shoot apices and analysing the expression of 38 candidate genes. RESULTS: Pyrosequencing from shoot apices generated a total of 34.7 Mbp and 535,544 reads enabling the identification of 43,087 clusters (41,512 contigs and 1,575 singletons). These data included 17,719 clusters (16,238 contigs and 1,575 singletons) exclusively from 454 sequencing reads, along with 25,368 hybrid clusters assembled with 454 sequences. The comparison of DNA libraries identified new candidate genes (n = 20) presenting differential expression between IAPAR59 and Rubi and/or drought conditions. Their expression was monitored in plagiotropic buds, together with those of other (n = 18) candidates genes. Under drought conditions, up-regulated expression was observed in IAPAR59 but not in Rubi for CaSTK1 (protein kinase), CaSAMT1 (SAM-dependent methyltransferase), CaSLP1 (plant development) and CaMAS1 (ABA biosynthesis). Interestingly, the expression of lipid-transfer protein (nsLTP) genes was also highly up-regulated under drought conditions in IAPAR59. This may have been related to the thicker cuticle observed on the abaxial leaf surface in IAPAR59 compared to Rubi. CONCLUSIONS: The full transcriptome assembly of C. arabica, followed by functional annotation, enabled us to identify differentially expressed genes related to drought conditions. Using these data, candidate genes were selected and their differential expression profiles were confirmed by qPCR experiments in plagiotropic buds of IAPAR59 and Rubi under drought conditions. As regards the genes up-regulated under drought conditions, specifically in the drought-tolerant IAPAR59, several corresponded to orphan genes but also to genes coding proteins involved in signal transduction pathways, as well as ABA and lipid metabolism, for example. The identification of these genes should help advance our understanding of the genetic determinism of drought tolerance in coffee.

Oasis: online analysis of small RNA deep sequencing data.

UNLABELLED: Oasis is a web application that allows for the fast and flexible online analysis of small-RNA-seq (sRNA-seq) data. It was designed for the end user in the lab, providing an easy-to-use web frontend including video tutorials, demo data and best practice step-by-step guidelines on how to analyze sRNA-seq data. Oasis' exclusive selling points are a differential expression module that allows for the multivariate analysis of samples, a classification module for robust biomarker detection and an advanced programming interface that supports the batch submission of jobs. Both modules include the analysis of novel miRNAs, miRNA targets and functional analyses including GO and pathway enrichment. Oasis generates downloadable interactive web reports for easy visualization, exploration and analysis of data on a local system. Finally, Oasis' modular workflow enables for the rapid (re-) analysis of data. AVAILABILITY AND IMPLEMENTATION: Oasis is implemented in Python, R, Java, PHP, C++ and JavaScript. It is freely available at http://oasis.dzne.de. CONTACT: stefan.bonn@dzne.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hoxa9 collaborates with E2A-PBX1 in mouse B cell leukemia in association with Flt3 activation and decrease of B cell gene expression.

BACKGROUND: The fusion protein E2A-PBX1 induces pediatric B cell leukemia in human. Previously, we reported oncogenic interactions between homeobox (Hox) genes and E2A-PBX1 in murine T cell leukemia. A proviral insertional mutagenesis screen with our E2A-PBX1 B cell leukemia mouse model identified Hoxa genes as potential collaborators to E2A-PBX1. Here we studied whether Hoxa9 could enhance E2A-PBX1 leukemogenesis. RESULTS: We show that Hoxa9 confers a proliferative advantage to E2A-PBX1 B cells. Transplantation experiments with E2A-PBX1 transgenic B cells overexpressing Hoxa9 isolated from bone marrow chimeras showed that Hoxa9 accelerates the generation of E2A-PBX1 B cell leukemia, but Hoxa9 is unable to transform B cells alone. Quantitative-reverse transcriptase polymerase chain reaction analysis demonstrated a strong repression of B cell specific genes in these E2A-PBX1/Hoxa9 leukemias in addition to Flt3 activation, indicating inhibition of B cell differentiation in combination with enhanced proliferation. Overexpression of Hoxa9 in established E2A-PBX1 mouse leukemic B cells resulted in a growth advantage in vitro, which was also characterized by an enhanced expression of Flt3. CONCLUSIONS: we show for the first time that Hoxa9 collaborates with E2A-PBX1 in the oncogenic transformation of B cells in a mouse model that involves Flt3 signaling, which is potentially relevant to human disease.

Large-scale analysis of differential gene expression in coffee genotypes resistant and susceptible to leaf miner-toward the identification of candidate genes for marker assisted-selection.

BACKGROUND: A successful development of herbivorous insects into plant tissues depends on coordination of metabolic processes. Plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense response. To understand the transcriptional basis of this response, we compare gene expression profiles of two coffee genotypes, susceptible and resistant to leaf miner (Leucoptera coffella). A total of 22000 EST sequences from the Coffee Genome Database were selected for a microarray analysis. Fluorescence probes were synthesized using mRNA from the infested and non-infested coffee plants. Array hybridization, scanning and data normalization were performed using Nimble Scan® e ArrayStar® platforms. Genes with foldchange values +/-2 were considered differentially expressed. A validation of 18 differentially expressed genes was performed in infected plants using qRT-PCR approach. RESULTS: The microarray analysis indicated that resistant plants differ in gene expression profile. We identified relevant transcriptional changes in defense strategies before insect attack. Expression changes (>2.00-fold) were found in resistant plants for 2137 genes (1266 up-regulated and 873 down-regulated). Up-regulated genes include those responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and maintenance. Also, our analyses indicated that differential expression profiles between resistant and susceptible genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. Validation of selected genes pointed to four selected genes as suitable candidates for markers in assisted-selection of novel cultivars. CONCLUSIONS: Our results show evidences that coffee defense responses against leaf-miner attack are balanced with other cellular functions. Also analyses suggest a major metabolic reconfiguration that highlights the complexity of this response.

In situ Patch-seq analysis of microglia reveals a lack of stress genes as found in FACS-isolated microglia.

We applied the patch-seq technique to harvest transcripts from individual microglial cells from cortex, hippocampus and corpus callosum of acute brain slices from adult mice. After recording membrane currents with the patch-clamp technique, the cytoplasm was collected via the pipette and underwent adapted SMART-seq2 preparation with subsequent sequencing. On average, 4138 genes were detected in 113 cells from hippocampus, corpus callosum and cortex, including microglia markers such as Tmem119, P2ry12 and Siglec-H. Comparing our dataset to previously published single cell mRNA sequencing data from FACS-isolated microglia indicated that two clusters of cells were absent in our patch-seq dataset. Pathway analysis of marker genes in FACS-specific clusters revealed association with microglial activation and stress response. This indicates that under normal conditions microglia in situ lack transcripts associated with a stress-response, and that the microglia-isolation procedure by mechanical dissociation and FACS triggers the expression of genes related to activation and stress.

Joint genotype inference with germline and somatic mutations.

The joint sequencing of related genomes has become an important means to discover rare variants. Normal-tumor genome pairs are routinely sequenced together to find somatic mutations and their associations with different cancers. Parental and sibling genomes reveal de novo germline mutations and inheritance patterns related to Mendelian diseases.Acute lymphoblastic leukemia (ALL) is the most common paediatric cancer and the leading cause of cancer-related death among children. With the aim of uncovering the full spectrum of germline and somatic genetic alterations in childhood ALL genomes, we conducted whole-exome re-sequencing on a unique cohort of over 120 exomes of childhood ALL quartets, each comprising a patient's tumor and matched-normal material, and DNA from both parents. We developed a general probabilistic model for such quartet sequencing reads mapped to the reference human genome. The model is used to infer joint genotypes at homologous loci across a normal-tumor genome pair and two parental genomes.We describe the algorithms and data structures for genotype inference, model parameter training. We implemented the methods in an open-source software package (QUADGT) that uses the standard file formats of the 1000 Genomes Project. Our method's utility is illustrated on quartets from the ALL cohort.

Global analyses of Ceratocystis cacaofunesta mitochondria: from genome to proteome.

BACKGROUND: The ascomycete fungus Ceratocystis cacaofunesta is the causal agent of wilt disease in cacao, which results in significant economic losses in the affected producing areas. Despite the economic importance of the Ceratocystis complex of species, no genomic data are available for any of its members. Given that mitochondria play important roles in fungal virulence and the susceptibility/resistance of fungi to fungicides, we performed the first functional analysis of this organelle in Ceratocystis using integrated "omics" approaches. RESULTS: The C. cacaofunesta mitochondrial genome (mtDNA) consists of a single, 103,147-bp circular molecule, making this the second largest mtDNA among the Sordariomycetes. Bioinformatics analysis revealed the presence of 15 conserved genes and 37 intronic open reading frames in C. cacaofunesta mtDNA. Here, we predicted the mitochondrial proteome (mtProt) of C. cacaofunesta, which is comprised of 1,124 polypeptides - 52 proteins that are mitochondrially encoded and 1,072 that are nuclearly encoded. Transcriptome analysis revealed 33 probable novel genes. Comparisons among the Gene Ontology results of the predicted mtProt of C. cacaofunesta, Neurospora crassa and Saccharomyces cerevisiae revealed no significant differences. Moreover, C. cacaofunesta mitochondria were isolated, and the mtProt was subjected to mass spectrometric analysis. The experimental proteome validated 27% of the predicted mtProt. Our results confirmed the existence of 110 hypothetical proteins and 7 novel proteins of which 83 and 1, respectively, had putative mitochondrial localization. CONCLUSIONS: The present study provides the first partial genomic analysis of a species of the Ceratocystis genus and the first predicted mitochondrial protein inventory of a phytopathogenic fungus. In addition to the known mitochondrial role in pathogenicity, our results demonstrated that the global function analysis of this organelle is similar in pathogenic and non-pathogenic fungi, suggesting that its relevance in the lifestyle of these organisms should be based on a small number of specific proteins and/or with respect to differential gene regulation. In this regard, particular interest should be directed towards mitochondrial proteins with unknown function and the novel protein that might be specific to this species. Further functional characterization of these proteins could enhance our understanding of the role of mitochondria in phytopathogenicity.