Early human fetal lung atlas reveals the temporal dynamics of epithelial cell plasticity.

Studying human fetal lungs can inform how developmental defects and disease states alter the function of the lungs. Here, we sequenced >150,000 single cells from 19 healthy human pseudoglandular fetal lung tissues ranging between gestational weeks 10-19. We capture dynamic developmental trajectories from progenitor cells that express abundant levels of the cystic fibrosis conductance transmembrane regulator (CFTR). These cells give rise to multiple specialized epithelial cell types. Combined with spatial transcriptomics, we show temporal regulation of key signalling pathways that may drive the temporal and spatial emergence of specialized epithelial cells including ciliated and pulmonary neuroendocrine cells. Finally, we show that human pluripotent stem cell-derived fetal lung models contain CFTR-expressing progenitor cells that capture similar lineage developmental trajectories as identified in the native tissue. Overall, this study provides a comprehensive single-cell atlas of the developing human lung, outlining the temporal and spatial complexities of cell lineage development and benchmarks fetal lung cultures from human pluripotent stem cell differentiations to similar developmental window.

Evaluation of single-cell RNA-seq clustering algorithms on cancer tumor datasets.

Tumors are complex biological entities that comprise cell types of different origins, with different mutational profiles and different patterns of transcriptional dysregulation. The exploration of data related to cancer biology requires careful analytical methods to reflect the heterogeneity of cell populations in cancer samples. Single-cell techniques are now able to capture the transcriptional profiles of individual cells. However, the complexity of RNA-seq data, especially in cancer samples, makes it challenging to cluster single-cell profiles into groups that reflect the underlying cell types. We have developed a framework for a systematic examination of single-cell RNA-seq clustering algorithms for cancer data, which uses a range of well-established metrics to generate a unified quality score and algorithm ranking. To demonstrate this framework, we examined clustering performance of 15 different single-cell RNA-seq clustering algorithms on eight different cancer datasets. Our results suggest that the single-cell RNA-seq clustering algorithms fall into distinct groups by performance, with the highest clustering quality on non-malignant cells achieved by three algorithms: Seurat, bigSCale and Cell Ranger. However, for malignant cells, two additional algorithms often reach a better performance, namely Monocle and SC3. Their ability to detect known rare cell types was also among the best, along with Seurat. Our approach and results can be used by a broad audience of practitioners who analyze single-cell transcriptomic data in cancer research.

MetaFusion: a high-confidence metacaller for filtering and prioritizing RNA-seq gene fusion candidates.

MOTIVATION: Current fusion detection tools use diverse calling approaches and provide varying results, making selection of the appropriate tool challenging. Ensemble fusion calling techniques appear promising; however, current options have limited accessibility and function. RESULTS: MetaFusion is a flexible metacalling tool that amalgamates outputs from any number of fusion callers. Individual caller results are standardized by conversion into the new file type Common Fusion Format. Calls are annotated, merged using graph clustering, filtered and ranked to provide a final output of high-confidence candidates. MetaFusion consistently achieves higher precision and recall than individual callers on real and simulated datasets, and reaches up to 100% precision, indicating that ensemble calling is imperative for high-confidence results. MetaFusion uses FusionAnnotator to annotate calls with information from cancer fusion databases and is provided with a Benchmarking Toolkit to calibrate new callers. AVAILABILITY AND IMPLEMENTATION: MetaFusion is freely available at https://github.com/ccmbioinfo/MetaFusion. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

EpigenCentral: Portal for DNA methylation data analysis and classification in rare diseases.

Epigenetic processes play a key role in regulating gene expression. Genetic variants that disrupt chromatin-modifying proteins are associated with a broad range of diseases, some of which have specific epigenetic patterns, such as aberrant DNA methylation (DNAm), which may be used as disease biomarkers. While much of the epigenetic research has focused on cancer, there is a paucity of resources devoted to neurodevelopmental disorders (NDDs), which include autism spectrum disorder and many rare, clinically overlapping syndromes. To address this challenge, we created EpigenCentral, a free web resource for biomedical researchers, molecular diagnostic laboratories, and clinical practitioners to perform the interactive classification and analysis of DNAm data related to NDDs. It allows users to search for known disease-associated patterns in their DNAm data, classify genetic variants as pathogenic or benign to assist in molecular diagnostics, or analyze patterns of differential methylation in their data through a simple web form. EpigenCentral is freely available at http://epigen.ccm.sickkids.ca/.

DNA Methylation Signature for EZH2 Functionally Classifies Sequence Variants in Three PRC2 Complex Genes.

Weaver syndrome (WS), an overgrowth/intellectual disability syndrome (OGID), is caused by pathogenic variants in the histone methyltransferase EZH2, which encodes a core component of the Polycomb repressive complex-2 (PRC2). Using genome-wide DNA methylation (DNAm) data for 187 individuals with OGID and 969 control subjects, we show that pathogenic variants in EZH2 generate a highly specific and sensitive DNAm signature reflecting the phenotype of WS. This signature can be used to distinguish loss-of-function from gain-of-function missense variants and to detect somatic mosaicism. We also show that the signature can accurately classify sequence variants in EED and SUZ12, which encode two other core components of PRC2, and predict the presence of pathogenic variants in undiagnosed individuals with OGID. The discovery of a functionally relevant signature with utility for diagnostic classification of sequence variants in EZH2, EED, and SUZ12 supports the emerging paradigm shift for implementation of DNAm signatures into diagnostics and translational research.

Prevalence and Clinical Features of Inflammatory Bowel Diseases Associated With Monogenic Variants, Identified by Whole-Exome Sequencing in 1000 Children at a Single Center.

BACKGROUND & AIMS: A proportion of infants and young children with inflammatory bowel diseases (IBDs) have subtypes associated with a single gene variant (monogenic IBD). We aimed to determine the prevalence of monogenic disease in a cohort of pediatric patients with IBD. METHODS: We performed whole-exome sequencing analyses of blood samples from an unselected cohort of 1005 children with IBD, aged 0-18 years (median age at diagnosis, 11.96 years) at a single center in Canada and their family members (2305 samples total). Variants believed to cause IBD were validated using Sanger sequencing. Biopsies from patients were analyzed by immunofluorescence and histochemical analyses. RESULTS: We identified 40 rare variants associated with 21 monogenic genes among 31 of the 1005 children with IBD (including 5 variants in XIAP, 3 in DOCK8, and 2 each in FOXP3, GUCY2C, and LRBA). These variants occurred in 7.8% of children younger than 6 years and 2.3% of children aged 6-18 years. Of the 17 patients with monogenic Crohn's disease, 35% had abdominal pain, 24% had nonbloody loose stool, 18% had vomiting, 18% had weight loss, and 5% had intermittent bloody loose stool. The 14 patients with monogenic ulcerative colitis or IBD-unclassified received their diagnosis at a younger age, and their most predominant feature was bloody loose stool (78%). Features associated with monogenic IBD, compared to cases of IBD not associated with a single variant, were age of onset younger than 2 years (odds ratio [OR], 6.30; P = .020), family history of autoimmune disease (OR, 5.12; P = .002), extra-intestinal manifestations (OR, 15.36; P < .0001), and surgery (OR, 3.42; P = .042). Seventeen patients had variants in genes that could be corrected with allogeneic hematopoietic stem cell transplantation. CONCLUSIONS: In whole-exome sequencing analyses of more than 1000 children with IBD at a single center, we found that 3% had rare variants in genes previously associated with pediatric IBD. These were associated with different IBD phenotypes, and 1% of the patients had variants that could be potentially corrected with allogeneic hematopoietic stem cell transplantation. Monogenic IBD is rare, but should be considered in analysis of all patients with pediatric onset of IBD.

Correction to: DNA methylation signature is prognostic of choroid plexus tumor aggressiveness.

After publication of the original article [1], authors have requested to add a 'J' as middle name for Richard Gilbertson. Hence, full name should be Richard J Gilbertson.

DNA methylation signature is prognostic of choroid plexus tumor aggressiveness.

BACKGROUND: Histological grading of choroid plexus tumors (CPTs) remains the best prognostic tool to distinguish between aggressive choroid plexus carcinoma (CPC) and the more benign choroid plexus papilloma (CPP) or atypical choroid plexus papilloma (aCPP); however, these distinctions can be challenging. Standard treatment of CPC is very aggressive and often leads to severe damage to the young child's brain. Therefore, it is crucial to distinguish between CPC and less aggressive entities (CPP or aCPP) to avoid unnecessary exposure of the young patient to neurotoxic therapy. To better stratify CPTs, we utilized DNA methylation (DNAm) to identify prognostic epigenetic biomarkers for CPCs. METHODS: We obtained DNA methylation profiles of 34 CPTs using the HumanMethylation450 BeadChip from Illumina, and the data was analyzed using the Illumina Genome Studio analysis software. Validation of differentially methylated CpG sites chosen as biomarkers was performed using pyrosequencing analysis on additional 22 CPTs. Sensitivity testing of the CPC DNAm signature was performed on a replication cohort of 61 CPT tumors obtained from Neuropathology, University Hospital Münster, Germany. RESULTS: Generated genome-wide DNAm profiles of CPTs showed significant differences in DNAm between CPCs and the CPPs or aCPPs. The prediction of clinical outcome could be improved by combining the DNAm profile with the mutational status of TP53. CPCs with homozygous TP53 mutations clustered as a group separate from those carrying a heterozygous TP53 mutation or CPCs with wild type TP53 (TP53-wt) and showed the worst survival outcome. Specific DNAm signatures for CPCs revealed AK1, PER2, and PLSCR4 as potential biomarkers for CPC that can be used to improve molecular stratification for diagnosis and treatment. CONCLUSIONS: We demonstrate that combining specific DNAm signature for CPCs with histological approaches better differentiate aggressive tumors from those that are not life threatening. These findings have important implications for future prognostic risk prediction in clinical disease management.

Don't brush off buccal data heterogeneity.

Buccal epithelial cells are among the most clinically accessible tissues and are increasingly being used to identify epigenetic disease patterns. However, substantial variation in buccal DNA methylation patterns indicates heterogeneity of cell types within and between samples, raising questions of data quality. We systematically estimated cell-type composition for a large collection of buccal and saliva samples from 11 published studies of DNA methylation. In these we identified numerous cases of buccal samples with questionable purity, which may be affected by sampling from individuals with neurodevelopmental disorders, and by the brushes used for sample collection. Further challenges are involved in comparisons with tissues such as saliva, in which buccal component varies widely. We propose a reference-based method of correcting for buccal purity that reduces unwanted variation while preserving cross-tissue differences. Our work demonstrates the wide variation of buccal quality in epigenetic studies and suggests a possible approach to overcome this issue.

CHARGE and Kabuki Syndromes: Gene-Specific DNA Methylation Signatures Identify Epigenetic Mechanisms Linking These Clinically Overlapping Conditions.

Epigenetic dysregulation has emerged as a recurring mechanism in the etiology of neurodevelopmental disorders. Two such disorders, CHARGE and Kabuki syndromes, result from loss of function mutations in chromodomain helicase DNA-binding protein 7 (CHD7LOF) and lysine (K) methyltransferase 2D (KMT2DLOF), respectively. Although these two syndromes are clinically distinct, there is significant phenotypic overlap. We therefore expected that epigenetically driven developmental pathways regulated by CHD7 and KMT2D would overlap and that DNA methylation (DNAm) alterations downstream of the mutations in these genes would identify common target genes, elucidating a mechanistic link between these two conditions, as well as specific target genes for each disorder. Genome-wide DNAm profiles in individuals with CHARGE and Kabuki syndromes with CHD7LOF or KMT2DLOF identified distinct sets of DNAm differences in each of the disorders, which were used to generate two unique, highly specific and sensitive DNAm signatures. These DNAm signatures were able to differentiate pathogenic mutations in these two genes from controls and from each other. Analysis of the DNAm targets in each gene-specific signature identified both common gene targets, including homeobox A5 (HOXA5), which could account for some of the clinical overlap in CHARGE and Kabuki syndromes, as well as distinct gene targets. Our findings demonstrate how characterization of the epigenome can contribute to our understanding of disease pathophysiology for epigenetic disorders, paving the way for explorations of novel therapeutics.

Role of STAT5 and epigenetics in lactation-associated upregulation of multidrug transporter ABCG2 in the mammary gland.

The multidrug resistance efflux transporter ATP-binding cassette subfamily G member 2 (ABCG2) is not only overexpressed in certain drug-resistant cancers but is also highly expressed in the mammary gland during lactation, carrying xenobiotics and nutrients into milk. We sought to investigate the molecular mechanisms involved in the upregulation of ABCG2 during lactation. Expression profiling of different mouse Abcg2 mRNA isoforms (E1a, E1b, and E1c) revealed that E1b is predominantly expressed and induced in the lactating mouse mammary gland. Despite this induction, analyses of CpG methylation status and published ChIP-seq datasets reveal that E1b promoter sequences in the virgin gland are already hypomethylated and marked with the open chromatin histone mark H3K4me2. Using a forced-weaning model to shut down lactation, we found that within 24 h there was a significant reduction in Abcg2 mRNA expression and a loss of signal transducer and activator of transcription-5 (STAT5) occupancy at the mouse Abcg2 gene. Luciferase reporter assays further showed that some of these STAT5-binding regions that contained interferon-γ-activated sequence (GAS) motifs function as an enhancer after prolactin treatment. We conclude that Abcg2 is already poised for expression in the virgin mammary gland and that STAT5 plays an important role in Abcg2 expression during lactation.

Intercellular network structure and regulatory motifs in the human hematopoietic system.

The hematopoietic system is a distributed tissue that consists of functionally distinct cell types continuously produced through hematopoietic stem cell (HSC) differentiation. Combining genomic and phenotypic data with high-content experiments, we have built a directional cell-cell communication network between 12 cell types isolated from human umbilical cord blood. Network structure analysis revealed that ligand production is cell type dependent, whereas ligand binding is promiscuous. Consequently, additional control strategies such as cell frequency modulation and compartmentalization were needed to achieve specificity in HSC fate regulation. Incorporating the in vitro effects (quiescence, self-renewal, proliferation, or differentiation) of 27 HSC binding ligands into the topology of the cell-cell communication network allowed coding of cell type-dependent feedback regulation of HSC fate. Pathway enrichment analysis identified intracellular regulatory motifs enriched in these cell type- and ligand-coupled responses. This study uncovers cellular mechanisms of hematopoietic cell feedback in HSC fate regulation, provides insight into the design principles of the human hematopoietic system, and serves as a foundation for the analysis of intercellular regulation in multicellular systems.

Expanding the landscape of chromatin modification (CM)-related functional domains and genes in human.

Chromatin modification (CM) plays a key role in regulating transcription, DNA replication, repair and recombination. However, our knowledge of these processes in humans remains very limited. Here we use computational approaches to study proteins and functional domains involved in CM in humans. We analyze the abundance and the pair-wise domain-domain co-occurrences of 25 well-documented CM domains in 5 model organisms: yeast, worm, fly, mouse and human. Results show that domains involved in histone methylation, DNA methylation, and histone variants are remarkably expanded in metazoan, reflecting the increased demand for cell type-specific gene regulation. We find that CM domains tend to co-occur with a limited number of partner domains and are hence not promiscuous. This property is exploited to identify 47 potentially novel CM domains, including 24 DNA-binding domains, whose role in CM has received little attention so far. Lastly, we use a consensus Machine Learning approach to predict 379 novel CM genes (coding for 329 proteins) in humans based on domain compositions. Several of these predictions are supported by very recent experimental studies and others are slated for experimental verification. Identification of novel CM genes and domains in humans will aid our understanding of fundamental epigenetic processes that are important for stem cell differentiation and cancer biology. Information on all the candidate CM domains and genes reported here is publicly available.