Conserved stromal-immune cell circuits secure B cell homeostasis and function.

B cell zone reticular cells (BRCs) form stable microenvironments that direct efficient humoral immunity with B cell priming and memory maintenance being orchestrated across lymphoid organs. However, a comprehensive understanding of systemic humoral immunity is hampered by the lack of knowledge of global BRC sustenance, function and major pathways controlling BRC-immune cell interactions. Here we dissected the BRC landscape and immune cell interactome in human and murine lymphoid organs. In addition to the major BRC subsets underpinning the follicle, including follicular dendritic cells, PI16+ RCs were present across organs and species. As well as BRC-produced niche factors, immune cell-driven BRC differentiation and activation programs governed the convergence of shared BRC subsets, overwriting tissue-specific gene signatures. Our data reveal that a canonical set of immune cell-provided cues enforce bidirectional signaling programs that sustain functional BRC niches across lymphoid organs and species, thereby securing efficient humoral immunity.

A model of human neural networks reveals NPTX2 pathology in ALS and FTLD.

Human cellular models of neurodegeneration require reproducibility and longevity, which is necessary for simulating age-dependent diseases. Such systems are particularly needed for TDP-43 proteinopathies1, which involve human-specific mechanisms2-5 that cannot be directly studied in animal models. Here, to explore the emergence and consequences of TDP-43 pathologies, we generated induced pluripotent stem cell-derived, colony morphology neural stem cells (iCoMoNSCs) via manual selection of neural precursors6. Single-cell transcriptomics and comparison to independent neural stem cells7 showed that iCoMoNSCs are uniquely homogenous and self-renewing. Differentiated iCoMoNSCs formed a self-organized multicellular system consisting of synaptically connected and electrophysiologically active neurons, which matured into long-lived functional networks (which we designate iNets). Neuronal and glial maturation in iNets was similar to that of cortical organoids8. Overexpression of wild-type TDP-43 in a minority of neurons within iNets led to progressive fragmentation and aggregation of the protein, resulting in a partial loss of function and neurotoxicity. Single-cell transcriptomics revealed a novel set of misregulated RNA targets in TDP-43-overexpressing neurons and in patients with TDP-43 proteinopathies exhibiting a loss of nuclear TDP-43. The strongest misregulated target encoded the synaptic protein NPTX2, the levels of which are controlled by TDP-43 binding on its 3' untranslated region. When NPTX2 was overexpressed in iNets, it exhibited neurotoxicity, whereas correcting NPTX2 misregulation partially rescued neurons from TDP-43-induced neurodegeneration. Notably, NPTX2 was consistently misaccumulated in neurons from patients with amyotrophic lateral sclerosis and frontotemporal lobar degeneration with TDP-43 pathology. Our work directly links TDP-43 misregulation and NPTX2 accumulation, thereby revealing a TDP-43-dependent pathway of neurotoxicity.

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes.

Although transcriptomics data is typically used to analyze mature spliced mRNA, recent attention has focused on jointly investigating spliced and unspliced (or precursor-) mRNA, which can be used to study gene regulation and changes in gene expression production. Nonetheless, most methods for spliced/unspliced inference (such as RNA velocity tools) focus on individual samples, and rarely allow comparisons between groups of samples (e.g. healthy vs. diseased). Furthermore, this kind of inference is challenging, because spliced and unspliced mRNA abundance is characterized by a high degree of quantification uncertainty, due to the prevalence of multi-mapping reads, ie reads compatible with multiple transcripts (or genes), and/or with both their spliced and unspliced versions. Here, we present DifferentialRegulation, a Bayesian hierarchical method to discover changes between experimental conditions with respect to the relative abundance of unspliced mRNA (over the total mRNA). We model the quantification uncertainty via a latent variable approach, where reads are allocated to their gene/transcript of origin, and to the respective splice version. We designed several benchmarks where our approach shows good performance, in terms of sensitivity and error control, vs. state-of-the-art competitors. Importantly, our tool is flexible, and works with both bulk and single-cell RNA-sequencing data. DifferentialRegulation is distributed as a Bioconductor R package.

DESpace: spatially variable gene detection via differential expression testing of spatial clusters.

MOTIVATION: Spatially resolved transcriptomics (SRT) enables scientists to investigate spatial context of mRNA abundance, including identifying spatially variable genes (SVGs), i.e. genes whose expression varies across the tissue. Although several methods have been proposed for this task, native SVG tools cannot jointly model biological replicates, or identify the key areas of the tissue affected by spatial variability. RESULTS: Here, we introduce DESpace, a framework, based on an original application of existing methods, to discover SVGs. In particular, our approach inputs all types of SRT data, summarizes spatial information via spatial clusters, and identifies spatially variable genes by performing differential gene expression testing between clusters. Furthermore, our framework can identify (and test) the main cluster of the tissue affected by spatial variability; this allows scientists to investigate spatial expression changes in specific areas of interest. Additionally, DESpace enables joint modeling of multiple samples (i.e. biological replicates); compared to inference based on individual samples, this approach increases statistical power, and targets SVGs with consistent spatial patterns across replicates. Overall, in our benchmarks, DESpace displays good true positive rates, controls for false positive and false discovery rates, and is computationally efficient. AVAILABILITY AND IMPLEMENTATION: DESpace is freely distributed as a Bioconductor R package at https://bioconductor.org/packages/DESpace.

Ten simple rules for computational biologists collaborating with wet lab researchers.

Computational biologists are frequently engaged in collaborative data analysis with wet lab researchers. These interdisciplinary projects, as necessary as they are to the scientific endeavor, can be surprisingly challenging due to cultural differences in operations and values. In this Ten Simple Rules guide, we aim to help dry lab researchers identify sources of friction and provide actionable tools to facilitate respectful, open, transparent, and rewarding collaborations.

Phase I Trial Characterizing the Pharmacokinetic Profile of N-803, a Chimeric IL-15 Superagonist, in Healthy Volunteers.

The oncotherapeutic promise of IL-15, a potent immunostimulant, is limited by a short serum t 1/2 The fusion protein N-803 is a chimeric IL-15 superagonist that has a >20-fold longer in vivo t 1/2 versus IL-15. This phase 1 study characterized the pharmacokinetic (PK) profile and safety of N-803 after s.c. administration to healthy human volunteers. Volunteers received two doses of N-803, and after each dose, PK and safety were assessed for 9 d. The primary endpoint was the N-803 PK profile, the secondary endpoint was safety, and immune cell levels and immunogenicity were measures of interest. Serum N-803 concentrations peaked 4 h after administration and declined with a t 1/2 of ∼20 h. N-803 did not cause treatment-emergent serious adverse events (AEs) or grade ≥3 AEs. Injection site reactions, chills, and pyrexia were the most common AEs. Administration of N-803 was well tolerated and accompanied by proliferation of NK cells and CD8+ T cells and sustained increases in the number of NK cells. Our results suggest that N-803 administration can potentiate antitumor immunity.

TCF/LEF dependent and independent transcriptional regulation of Wnt/ß-catenin target genes.

During canonical Wnt signalling, the activity of nuclear ß-catenin is largely mediated by the TCF/LEF family of transcription factors. To challenge this view, we used the CRISPR/Cas9 genome editing approach to generate HEK 293T cell clones lacking all four TCF/LEF genes. By performing unbiased whole transcriptome sequencing analysis, we found that a subset of ß-catenin transcriptional targets did not require TCF/LEF factors for their regulation. Consistent with this finding, we observed in a genome-wide analysis that ß-catenin occupied specific genomic regions in the absence of TCF/LEF Finally, we revealed the existence of a transcriptional activity of ß-catenin that specifically appears when TCF/LEF factors are absent, and refer to this as ß-catenin-GHOST response. Collectively, this study uncovers a previously neglected modus operandi of ß-catenin that bypasses the TCF/LEF transcription factors.

Validation of hypermethylated DNA regions found in colorectal cancers as potential aging-independent biomarkers of precancerous colorectal lesions.

BACKGROUND: We previously identified 16,772 colorectal cancer-associated hypermethylated DNA regions that were also detectable in precancerous colorectal lesions (preCRCs) and unrelated to normal mucosal aging. We have now conducted a study to validate 990 of these differentially methylated DNA regions (DMRs) in a new series of preCRCs. METHODS: We used targeted bisulfite sequencing to validate these 990 potential biomarkers in 59 preCRC tissue samples (41 conventional adenomas, 18 sessile serrated lesions), each with a patient-matched normal mucosal sample. Based on differential DNA methylation tests, a panel of candidate DMRs was chosen on a subset of our cohort and then validated on the remaining part of our cohort and two publicly available datasets with respect to their stratifying potential between preCRCs and normal mucosa. RESULTS: Strong statistical significance for the difference in methylation levels was observed across the full set of 990 investigated DMRs. From these, a selected candidate panel of 30 DMRs correctly identified 58/59 tumors (area under the receiver operating curve: 0.998). CONCLUSIONS: These validated DNA hypermethylation markers can be exploited to develop more accurate noninvasive colorectal tumor screening assays.

MiR-CLIP reveals iso-miR selective regulation in the miR-124 targetome.

Many microRNAs regulate gene expression via atypical mechanisms, which are difficult to discern using native cross-linking methods. To ascertain the scope of non-canonical miRNA targeting, methods are needed that identify all targets of a given miRNA. We designed a new class of miR-CLIP probe, whereby psoralen is conjugated to the 3p arm of a pre-microRNA to capture targetomes of miR-124 and miR-132 in HEK293T cells. Processing of pre-miR-124 yields miR-124 and a 5'-extended isoform, iso-miR-124. Using miR-CLIP, we identified overlapping targetomes from both isoforms. From a set of 16 targets, 13 were differently inhibited at mRNA/protein levels by the isoforms. Moreover, delivery of pre-miR-124 into cells repressed these targets more strongly than individual treatments with miR-124 and iso-miR-124, suggesting that isomirs from one pre-miRNA may function synergistically. By mining the miR-CLIP targetome, we identified nine G-bulged target-sites that are regulated at the protein level by miR-124 but not isomiR-124. Using structural data, we propose a model involving AGO2 helix-7 that suggests why only miR-124 can engage these sites. In summary, access to the miR-124 targetome via miR-CLIP revealed for the first time how heterogeneous processing of miRNAs combined with non-canonical targeting mechanisms expand the regulatory range of a miRNA.

T-cell acute leukaemia exhibits dynamic interactions with bone marrow microenvironments.

It is widely accepted that complex interactions between cancer cells and their surrounding microenvironment contribute to disease development, chemo-resistance and disease relapse. In light of this observed interdependency, novel therapeutic interventions that target specific cancer stroma cell lineages and their interactions are being sought. Here we studied a mouse model of human T-cell acute lymphoblastic leukaemia (T-ALL) and used intravital microscopy to monitor the progression of disease within the bone marrow at both the tissue-wide and single-cell level over time, from bone marrow seeding to development/selection of chemo-resistance. We observed highly dynamic cellular interactions and promiscuous distribution of leukaemia cells that migrated across the bone marrow, without showing any preferential association with bone marrow sub-compartments. Unexpectedly, this behaviour was maintained throughout disease development, from the earliest bone marrow seeding to response and resistance to chemotherapy. Our results reveal that T-ALL cells do not depend on specific bone marrow microenvironments for propagation of disease, nor for the selection of chemo-resistant clones, suggesting that a stochastic mechanism underlies these processes. Yet, although T-ALL infiltration and progression are independent of the stroma, accumulated disease burden leads to rapid, selective remodelling of the endosteal space, resulting in a complete loss of mature osteoblastic cells while perivascular cells are maintained. This outcome leads to a shift in the balance of endogenous bone marrow stroma, towards a composition associated with less efficient haematopoietic stem cell function. This novel, dynamic analysis of T-ALL interactions with the bone marrow microenvironment in vivo, supported by evidence from human T-ALL samples, highlights that future therapeutic interventions should target the migration and promiscuous interactions of cancer cells with the surrounding microenvironment, rather than specific bone marrow stroma, to combat the invasion by and survival of chemo-resistant T-ALL cells.

Censcyt: censored covariates in differential abundance analysis in cytometry.

BACKGROUND: Innovations in single cell technologies have lead to a flurry of datasets and computational tools to process and interpret them, including analyses of cell composition changes and transition in cell states. The diffcyt workflow for differential discovery in cytometry data consist of several steps, including preprocessing, cell population identification and differential testing for an association with a binary or continuous covariate. However, the commonly measured quantity of survival time in clinical studies often results in a censored covariate where classical differential testing is inapplicable. RESULTS: To overcome this limitation, multiple methods to directly include censored covariates in differential abundance analysis were examined with the use of simulation studies and a case study. Results show that multiple imputation based methods offer on-par performance with the Cox proportional hazards model in terms of sensitivity and error control, while offering flexibility to account for covariates. The tested methods are implemented in the R package censcyt as an extension of diffcyt and are available at https://bioconductor.org/packages/censcyt . CONCLUSION: Methods for the direct inclusion of a censored variable as a predictor in GLMMs are a valid alternative to classical survival analysis methods, such as the Cox proportional hazard model, while allowing for more flexibility in the differential analysis.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells.

BACKGROUND: Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown. RESULTS: The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14-29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D. CONCLUSION: Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

ARPEGGIO: Automated Reproducible Polyploid EpiGenetic GuIdance workflOw.

BACKGROUND: Whole genome duplication (WGD) events are common in the evolutionary history of many living organisms. For decades, researchers have been trying to understand the genetic and epigenetic impact of WGD and its underlying molecular mechanisms. Particular attention was given to allopolyploid study systems, species resulting from an hybridization event accompanied by WGD. Investigating the mechanisms behind the survival of a newly formed allopolyploid highlighted the key role of DNA methylation. With the improvement of high-throughput methods, such as whole genome bisulfite sequencing (WGBS), an opportunity opened to further understand the role of DNA methylation at a larger scale and higher resolution. However, only a few studies have applied WGBS to allopolyploids, which might be due to lack of genomic resources combined with a burdensome data analysis process. To overcome these problems, we developed the Automated Reproducible Polyploid EpiGenetic GuIdance workflOw (ARPEGGIO): the first workflow for the analysis of epigenetic data in polyploids. This workflow analyzes WGBS data from allopolyploid species via the genome assemblies of the allopolyploid's parent species. ARPEGGIO utilizes an updated read classification algorithm (EAGLE-RC), to tackle the challenge of sequence similarity amongst parental genomes. ARPEGGIO offers automation, but more importantly, a complete set of analyses including spot checks starting from raw WGBS data: quality checks, trimming, alignment, methylation extraction, statistical analyses and downstream analyses. A full run of ARPEGGIO outputs a list of genes showing differential methylation. ARPEGGIO was made simple to set up, run and interpret, and its implementation ensures reproducibility by including both package management and containerization. RESULTS: We evaluated ARPEGGIO in two ways. First, we tested EAGLE-RC's performance with publicly available datasets given a ground truth, and we show that EAGLE-RC decreases the error rate by 3 to 4 times compared to standard approaches. Second, using the same initial dataset, we show agreement between ARPEGGIO's output and published results. Compared to other similar workflows, ARPEGGIO is the only one supporting polyploid data. CONCLUSIONS: The goal of ARPEGGIO is to promote, support and improve polyploid research with a reproducible and automated set of analyses in a convenient implementation. ARPEGGIO is available at https://github.com/supermaxiste/ARPEGGIO .

Bias, robustness and scalability in single-cell differential expression analysis.

Many methods have been used to determine differential gene expression from single-cell RNA (scRNA)-seq data. We evaluated 36 approaches using experimental and synthetic data and found considerable differences in the number and characteristics of the genes that are called differentially expressed. Prefiltering of lowly expressed genes has important effects, particularly for some of the methods developed for bulk RNA-seq data analysis. However, we found that bulk RNA-seq analysis methods do not generally perform worse than those developed specifically for scRNA-seq. We also present conquer, a repository of consistently processed, analysis-ready public scRNA-seq data sets that is aimed at simplifying method evaluation and reanalysis of published results. Each data set provides abundance estimates for both genes and transcripts, as well as quality control and exploratory analysis reports.

CD8+ T cells retain protective functions despite sustained inhibitory receptor expression during Epstein-Barr virus infection in vivo.

Epstein Barr virus (EBV) is one of the most ubiquitous human pathogens in the world, persistently infecting more than 90% of the adult human population. It drives some of the strongest human CD8+ T cell responses, which can be observed during symptomatic primary infection known as infectious mononucleosis (IM). Despite high viral loads and prolonged CD8+ T cell stimulation during IM, EBV enters latency and is under lifelong immune control in most individuals that experience this disease. We investigated whether changes in T cell function, as frequently characterized by PD-1 up-regulation, occur during IM due to the prolonged exposure to high antigen levels. We readily detected the expansion of PD-1 positive CD8+ T cells together with high frequencies of Tim-3, 2B4, and KLRG1 expression during IM and in mice with reconstituted human immune system components (huNSG mice) that had been infected with a high dose of EBV. These PD-1 positive CD8+ T cells, however, retained proliferation, cytokine production, and cytotoxic abilities. Multiple subsets of CD8+ T cells expanded during EBV infection, including PD-1+Tim-3+KLRG1+ cells that express CXCR5 and TCF-1 germinal center homing and memory markers, and may also contain BATF3. Moreover, blocking the PD-1 axis compromised EBV specific immune control and resulted in virus-associated lymphomagenesis. Finally, PD-1+, Tim-3+, and KLRG1+ CD8+ T cell expansion coincided with declining viral loads during low dose EBV infection. These findings suggest that EBV infection primes PD-1 positive CD8+ T cell populations that rely on this receptor axis for the efficient immune control of this ubiquitous human tumor virus.

Circulating neutrophil subsets in advanced lung cancer patients exhibit unique immune signature and relate to prognosis.

The accumulation of circulating low-density neutrophils (LDN) has been described in cancer patients and associated with tumor-supportive properties, as opposed to the high-density neutrophils (HDN). Here we aimed to evaluate the clinical significance of circulating LDN in lung cancer patients, and further assessed its diagnostic vs prognostic value. Using mass cytometry (CyTOF), we identified major subpopulations within the circulating LDN/HDN subsets and determined phenotypic modulations of these subsets along tumor progression. LDN were highly enriched in the low-density (LD) fraction of advanced lung cancer patients (median 7.0%; range 0.2%-80%, n = 64), but not in early stage patients (0.7%; 0.05%-6%; n = 35), healthy individuals (0.8%; 0%-3.5%; n = 15), or stable chronic obstructive pulmonary disease (COPD) patients (1.2%; 0.3%-7.4%, n = 13). Elevated LDN (>10%) remarkably related with poorer prognosis in late stage patients. We identified three main neutrophil subsets which proportions are markedly modified in cancer patients, with CD66b+ /CD10low /CXCR4+ /PDL1inter subset almost exclusively found in advanced lung cancer patients. We found substantial variability in subsets between patients, and demonstrated that HDN and LDN retain a degree of inherent spontaneous plasticity. Deep phenotypic characterization of cancer-related circulating neutrophils and their modulation along tumor progression is an important advancement in understanding the role of myeloid cells in lung cancer.

The Spinal Transcriptome after Cortical Stroke: In Search of Molecular Factors Regulating Spontaneous Recovery in the Spinal Cord.

In response to cortical stroke and unilateral corticospinal tract degeneration, compensatory sprouting of spared corticospinal fibers is associated with recovery of skilled movement in rodents. To date, little is known about the molecular mechanisms orchestrating this spontaneous rewiring. In this study, we provide insights into the molecular changes in the spinal cord tissue after large ischemic cortical injury in adult female mice, with a focus on factors that might influence the reinnervation process by contralesional corticospinal neurons. We mapped the area of cervical gray matter reinnervation by sprouting contralesional corticospinal axons after unilateral photothrombotic stroke of the motor cortex in mice using anterograde tracing. The mRNA profile of this reinnervation area was analyzed using whole-genome sequencing to identify differentially expressed genes at selected time points during the recovery process. Bioinformatic analysis revealed two phases of processes: early after stroke (4-7 d post-injury), the spinal transcriptome is characterized by inflammatory processes, including phagocytic processes as well as complement cascade activation. Microglia are specifically activated in the denervated corticospinal projection fields in this early phase. In a later phase (28-42 d post-injury), biological processes include tissue repair pathways with upregulated genes related to neurite outgrowth. Thus, the stroke-denervated spinal gray matter, in particular its intermediate laminae, represents a growth-promoting environment for sprouting corticospinal fibers originating from the contralesional motor cortex. This dataset provides a solid starting point for future studies addressing key elements of the post-stroke recovery process, with the goal to improve neuroregenerative treatment options for stroke patients.SIGNIFICANCE STATEMENT We show that the molecular changes in the spinal cord target tissue of the stroke-affected corticospinal tract are mainly defined by two phases: an early inflammatory phase during which microglia are specifically activated in the target area of reinnervating corticospinal motor neurons; and a late phase during which growth-promoting factors are upregulated which can influence the sprouting response, arborization, and synapse formation. By defining for the first time the endogenous molecular machinery in the stroke-denervated cervical spinal gray matter with a focus on promotors of axon growth through the growth-inhibitory adult CNS, this study will serve as a basis to address novel neuroregenerative treatment options for chronic stroke patients.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics.

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

Faithful mRNA splicing depends on the Prp19 complex subunit faint sausage and is required for tracheal branching morphogenesis in Drosophila.

Morphogenesis requires the dynamic regulation of gene expression, including transcription, mRNA maturation and translation. Dysfunction of the general mRNA splicing machinery can cause surprisingly specific cellular phenotypes, but the basis for these effects is not clear. Here, we show that the Drosophila faint sausage (fas) locus, which is implicated in epithelial morphogenesis and has previously been reported to encode a secreted immunoglobulin domain protein, in fact encodes a subunit of the spliceosome-activating Prp19 complex, which is essential for efficient pre-mRNA splicing. Loss of zygotic fas function globally impairs the efficiency of splicing, and is associated with widespread retention of introns in mRNAs and dramatic changes in gene expression. Surprisingly, despite these general effects, zygotic fas mutants show specific defects in tracheal cell migration during mid-embryogenesis when maternally supplied splicing factors have declined. We propose that tracheal branching, which relies on dynamic changes in gene expression, is particularly sensitive for efficient spliceosome function. Our results reveal an entry point to study requirements of the splicing machinery during organogenesis and provide a better understanding of disease phenotypes associated with mutations in general splicing factors.

The DNA hypermethylation phenotype of colorectal cancer liver metastases resembles that of the primary colorectal cancers.

BACKGROUND: Identifying molecular differences between primary and metastatic colorectal cancers-now possible with the aid of omics technologies-can improve our understanding of the biological mechanisms of cancer progression and facilitate the discovery of novel treatments for late-stage cancer. We compared the DNA methylomes of primary colorectal cancers (CRCs) and CRC metastases to the liver. Laser microdissection was used to obtain epithelial tissue (10 to 25 × 106 µm2) from sections of fresh-frozen samples of primary CRCs (n = 6), CRC liver metastases (n = 12), and normal colon mucosa (n = 3). DNA extracted from tissues was enriched for methylated sequences with a methylCpG binding domain (MBD) polypeptide-based protocol and subjected to deep sequencing. The performance of this protocol was compared with that of targeted enrichment for bisulfite sequencing used in a previous study of ours. RESULTS: MBD enrichment captured a total of 322,551 genomic regions (249.5 Mb or ~ 7.8% of the human genome), which included over seven million CpG sites. A few of these regions were differentially methylated at an expected false discovery rate (FDR) of 5% in neoplastic tissues (primaries: 0.67%, i.e., 2155 regions containing 279,441 CpG sites; liver metastases: 1%, i.e., 3223 regions containing 312,723 CpG sites) as compared with normal mucosa samples. Most of the differentially methylated regions (DMRs; 94% in primaries; 70% in metastases) were hypermethylated, and almost 80% of these (1882 of 2396) were present in both lesion types. At 5% FDR, no DMRs were detected in liver metastases vs. primary CRC. However, short regions of low-magnitude hypomethylation were frequent in metastases but rare in primaries. Hypermethylated DMRs were far more abundant in sequences classified as intragenic, gene-regulatory, or CpG shelves-shores-island segments, whereas hypomethylated DMRs were equally represented in extragenic (mainly, open-sea) and intragenic (mainly, gene bodies) sequences of the genome. Compared with targeted enrichment, MBD capture provided a better picture of the extension of CRC-associated DNA hypermethylation but was less powerful for identifying hypomethylation. CONCLUSIONS: Our findings demonstrate that the hypermethylation phenotype in CRC liver metastases remains similar to that of the primary tumor, whereas CRC-associated DNA hypomethylation probably undergoes further progression after the cancer cells have migrated to the liver.