The Genomics of Colorectal Cancer in Populations with African and European Ancestry

Blacks have a higher incidence of colorectal cancer (CRC) and worse survival rates when compared to Whites. Comprehensive genomic profiling was performed in 46,140 colorectal adenocarcinoma cases. Ancestry-informative markers identified 5,301 patients of African descent (AFR) and 33,770 patients of European descent (EUR). AFR were younger, had fewer MSI-H tumors, and had significantly more frequent alterations in KRAS, APC, and PIK3CA. AFR had increased frequency of KRAS mutations specifically KRAS G12D and KRAS G13. There were no differences in rates of actionable kinase driver alterations (HER2, MET, NTRK, ALK, ROS1, RET). In patients with young onset CRC (<50 years), AFR and EUR had similar frequency of MSI-H and TMB-H tumors, and strikingly different trends in APC mutations by age, as well as differences in MAPK pathway alterations. These findings inform treatment decisions, impact prognosis, and underscore the need for model systems representative of our diverse US population.

Vitamin D Deficiency During Development Permanently Alters Liver Cell Composition and Function

Vitamin D, a fat-soluble vitamin, plays a critical role in calcium homeostasis, the immune system, and normal development. Many epidemiological cohort studies globally have found high prevalence rates of vitamin D deficiency and insufficiency, recognized as an important health issue that needs to be solved. In particular, reproductive age and pregnant women low in vitamin D status may confer risks of diseases like obesity on their offspring. While observational studies have suggested associations between prenatal vitamin D deficiency and metabolic phenotypes in offspring, not yet determined is whether prenatal vitamin D deficiency permanently alters the development of the liver, a major metabolic organ. We tested the histopathology and the transcriptomic profiles of livers from male C57BL/6J mice exposed to prenatal vitamin D deficiency through a maternal dietary intervention model. We found that prenatal vitamin D deficiency increases the prevalence of histopathological changes in the liver, and alters its gene expression profile. Cell subtype proportion analysis showed that the liver of prenatal vitamin D deficiency alters non-parenchymal cells of the liver, specifically macrophages, a subset of endothelial cells, and dendritic cells. Our results indicate the long-term memory of prenatal vitamin D deficiency exposure in the adult liver, a potential contributor to offspring health risks.

Master Transcription Regulators and Transcription Factors Regulate Immune-Associated Differences Between Patients of African and European Ancestry With Colorectal Cancer

Background and Aims

Individuals of African (AFR) ancestry have a higher incidence of colorectal cancer (CRC) than those of European (EUR) ancestry and exhibit significant health disparities. Previous studies have noted differences in the tumor microenvironment between AFR and EUR patients with CRC. However, the molecular regulatory processes that underpin these immune differences remain largely unknown.

Methods

Multiomics analysis was carried out for 55 AFR and 456 EUR patients with microsatellite-stable CRC using The Cancer Genome Atlas. We evaluated the tumor microenvironment by using gene expression and methylation data, transcription factor, and master transcriptional regulator analysis to identify the cell signaling pathways mediating the observed phenotypic differences.

Results

We demonstrate that downregulated genes in AFR patients with CRC showed enrichment for canonical pathways, including chemokine signaling. Moreover, evaluation of the tumor microenvironment showed that cytotoxic lymphocytes and neutrophil cell populations are significantly decreased in AFR compared with EUR patients, suggesting AFR patients have an attenuated immune response. We further demonstrate that molecules called “master transcriptional regulators” (MTRs) play a critical role in regulating the expression of genes impacting key immune processes through an intricate signal transduction network mediated by disease-associated transcription factors (TFs). Furthermore, a core set of these MTRs and TFs showed a positive correlation with levels of cytotoxic lymphocytes and neutrophils across both AFR and EUR patients with CRC, thus suggesting their role in driving the immune infiltrate differences between the two ancestral groups.

Conclusion

Our study provides an insight into the intricate regulatory landscape of MTRs and TFs that orchestrate the differences in the tumor microenvironment between patients with CRC of AFR and EUR ancestry.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

GenomeDiver: a platform for phenotype-guided medical genomic diagnosis

Purpose:

Making a diagnosis from clinical genomic sequencing requires well-structured phenotypic data to guide genotype interpretation. A patient’s phenotypic features can be documented using the Human Phenotype Ontology (HPO), generating terms used to prioritize genes potentially causing the patient’s disease. We have developed GenomeDiver to provide a user interface for clinicians that allows more effective collaboration with the clinical diagnostic laboratory, with the goal of improving the success of the diagnostic process.

Methods:

GenomeDiver uses genomic data to prompt reverse phenotyping of patients undergoing genetic testing, enriching the amount and quality of structured phenotype data for the diagnostic laboratory, and helping clinicians to explore and flag diseases potentially causing their patient’s presentation.

Results:

We show how GenomeDiver communicates the clinician’s informed insights to the diagnostic lab in the form of HPO terms for interpretation of genomic sequencing data. We describe our user-driven design process, the engineering of the software for efficiency, security and portability, and examples of the performance of GenomeDiver using genomic testing data.

Conclusions:

GenomeDiver is a first step in a new approach to genomic diagnostics that enhances laboratory-clinician interactions, with the goal of directly engaging clinicians to improve the outcome of genomic diagnostic testing.

Disproportionate Vitamin A Deficiency in Women of Specific Ethnicities Linked to Differences in Allele Frequencies of Vitamin A-Related Polymorphisms

Background: While the current national prevalence rate of vitamin A deficiency (VAD) is estimated to be less than 1%, it is suggested that it varies between different ethnic groups and races within the U.S. We assessed the prevalence of VAD in pregnant women of different ethnic groups and tested these prevalence rates for associations with the vitamin A-related single nucleotide polymorphism (SNP) allele frequencies in each ethnic group. Methods: We analyzed two independent datasets of serum retinol levels with self-reported ethnicities and the differences of allele frequencies of the SNPs associated with vitamin A metabolism between groups in publicly available datasets. Results: Non-Hispanic Black and Hispanic pregnant women showed high VAD prevalence in both datasets. Interestingly, the VAD prevalence for Hispanic pregnant women significantly differed between datasets (p = 1.973 × 10⁻¹⁰, 95%CI 0.04–0.22). Alleles known to confer the risk of low serum retinol (rs10882272 C and rs738409 G) showed higher frequencies in the race/ethnicity groups with more VAD. Moreover, minor allele frequencies of a set of 39 previously reported SNPs associated with vitamin A metabolism were significantly different between the populations of different ancestries than those of randomly selected SNPs (p = 0.030). Conclusions: Our analysis confirmed that VAD prevalence varies between different ethnic groups/races and may be causally associated with genetic variants conferring risk for low retinol levels. Assessing genetic variant information prior to performing an effective nutrient supplementation program will help us plan more effective food-based interventions.

Cell type-specific chromatin accessibility analysis in the mouse and human brain

The Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) is becoming popular in the neuroscience field where chromatin regulation is thought to be involved in neurodevelopment, activity-dependent gene regulation, hormonal and environmental responses, and pathophysiology of neuropsychiatric disorders. The advantages of using ATAC-seq include a small amount of material needed, fast protocol, and the ability to capture a range of gene regulatory elements with a single assay. With increasing interest in chromatin research, it is an imperative to have feasible, reliable assays that are compatible with a range of neuroscience study designs. Here we tested three protocols for neuronal chromatin accessibility analysis, including a varying brain tissue freezing method followed by fluorescence-activated nuclei sorting (FANS) and ATAC-seq. Our study shows that the cryopreservation method impacts the number of open chromatin regions identified from frozen brain tissue using ATAC-seq. However, we show that all protocols generate consistent and robust data and enable the identification of functional regulatory elements in neuronal cells. Our study implies that the broad biological interpretation of chromatin accessibility data is not significantly affected by the freezing condition. We also reveal additional challenges of doing chromatin analysis on post-mortem human brain tissue. Overall, ATAC-seq coupled with FANS is a powerful method to capture cell-type-specific chromatin accessibility information in mouse and human brain. Our study provides alternative brain preservation methods that generate high-quality ATAC-seq data while fitting in different study designs, and further encourages the use of this method to uncover the role of epigenetic (dys)regulation in the brain.

GUÍA: a digital platform to facilitate result disclosure in genetic counseling

Purpose

Use of genomic sequencing is increasing at a pace that requires technological solutions to effectively meet the needs of a growing patient population. We developed GUÍA, a web-based application, to enhance the delivery of genomic results and related clinical information to patients and families.

Methods

GUÍA development occurred in five overlapping phases: formative research, content development, stakeholder/community member input, user interface design, and web application development. Development was informed by formative qualitative research involving parents (N = 22) whose children underwent genomic testing. Participants enrolled in the NYCKidSeq pilot study (N = 18) completed structured feedback interviews post–result disclosure using GUÍA. Genetic specialists, researchers, patients, and community stakeholders provided their perspectives on GUÍA’s design to ensure technical, cultural, and literacy appropriateness.

Results

NYCKidSeq participants responded positively to the use of GUÍA to deliver their children’s results. All participants (N = 10) with previous experience with genetic testing felt GUÍA improved result disclosure, and 17 (94%) participants said the content was clear.

Conclusion

GUÍA communicates complex genomic information in an understandable and personalized manner. Initial piloting demonstrated GUÍA’s utility for families enrolled in the NYCKidSeq pilot study. Findings from the NYCKidSeq clinical trial will provide insight into GUÍA’s effectiveness in communicating results among diverse, multilingual populations.

The NYCKidSeq project: study protocol for a randomized controlled trial incorporating genomics into the clinical care of diverse New York City children

BACKGROUND

Increasingly, genomics is informing clinical practice, but challenges remain for medical professionals lacking genetics expertise, and in access to and clinical utility of genomic testing for minority and underrepresented populations. The latter is a particularly pernicious problem due to the historical lack of inclusion of racially and ethnically diverse populations in genomic research and genomic medicine. A further challenge is the rapidly changing landscape of genetic tests and considerations of cost, interpretation, and diagnostic yield for emerging modalities like whole-genome sequencing.

METHODS

The NYCKidSeq project is a randomized controlled trial recruiting 1130 children and young adults predominantly from Harlem and the Bronx with suspected genetic disorders in three disease categories: neurologic, cardiovascular, and immunologic. Two clinical genetic tests will be performed for each participant, either proband, duo, or trio whole-genome sequencing (depending on sample availability) and proband targeted gene panels. Clinical utility, cost, and diagnostic yield of both testing modalities will be assessed. This study will evaluate the use of a novel, digital platform (GUÍA) to digitize the return of genomic results experience and improve participant understanding for English- and Spanish-speaking families. Surveys will collect data at three study visits: baseline (0 months), result disclosure visit (ROR1, + 3 months), and follow-up visit (ROR2, + 9 months). Outcomes will assess parental understanding of and attitudes toward receiving genomic results for their child and behavioral, psychological, and social impact of results. We will also conduct a pilot study to assess a digital tool called GenomeDiver designed to enhance communication between clinicians and genetic testing labs. We will evaluate GenomeDiver's ability to increase the diagnostic yield compared to standard practices, improve clinician's ability to perform targeted reverse phenotyping, and increase the efficiency of genetic testing lab personnel.

DISCUSSION

The NYCKidSeq project will contribute to the innovations and best practices in communicating genomic test results to diverse populations. This work will inform strategies for implementing genomic medicine in health systems serving diverse populations using methods that are clinically useful, technologically savvy, culturally sensitive, and ethically sound.

TRIAL REGISTRATION

ClinicalTrials.gov NCT03738098 . Registered on November 13, 2018 Trial Sponsor: Icahn School of Medicine at Mount Sinai Contact Name: Eimear Kenny, PhD (Principal Investigator) Address: Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Pl., Box 1003, New York, NY 10029 Email: eimear.kenny@mssm.edu.

The shape of gene expression distributions matter: how incorporating distribution shape improves the interpretation of cancer transcriptomic data

BACKGROUND

In genomics, we often assume that continuous data, such as gene expression, follow a specific kind of distribution. However we rarely stop to question the validity of this assumption, or consider how broadly applicable it may be to all genes that are in the transcriptome. Our study investigated the prevalence of a range of gene expression distributions in three different tumor types from the Cancer Genome Atlas (TCGA).

RESULTS

Surprisingly, the expression of less than 50% of all genes was Normally-distributed, with other distributions including Gamma, Bimodal, Cauchy, and Lognormal also represented. Most of the distribution categories contained genes that were significantly enriched for unique biological processes. Different assumptions based on the shape of the expression profile were used to identify genes that could discriminate between patients with good versus poor survival. The prognostic marker genes that were identified when the shape of the distribution was accounted for reflected functional insights into cancer biology that were not observed when standard assumptions were applied. We showed that when multiple types of distributions were permitted, i.e. the shape of the expression profile was used, the statistical classifiers had greater predictive accuracy for determining the prognosis of a patient versus those that assumed only one type of gene expression distribution.

CONCLUSIONS

Our results highlight the value of studying a gene's distribution shape to model heterogeneity of transcriptomic data and the impact on using analyses that permit more than one type of gene expression distribution. These insights would have been overlooked when using standard approaches that assume all genes follow the same type of distribution in a patient cohort.

Variant Classification Concordance using the ACMG-AMP Variant Interpretation Guidelines across Nine Genomic Implementation Research Studies

Harmonization of variant pathogenicity classification across laboratories is important for advancing clinical genomics. The two CLIA-accredited Electronic Medical Record and Genomics Network sequencing centers and the six CLIA-accredited laboratories and one research laboratory performing genome or exome sequencing in the Clinical Sequencing Evidence-Generating Research Consortium collaborated to explore current sources of discordance in classification. Eight laboratories each submitted 20 classified variants in the ACMG secondary finding v.2.0 genes. After removing duplicates, each of the 158 variants was annotated and independently classified by two additional laboratories using the ACMG-AMP guidelines. Overall concordance across three laboratories was assessed and discordant variants were reviewed via teleconference and email. The submitted variant set included 28 P/LP variants, 96 VUS, and 34 LB/B variants, mostly in cancer (40%) and cardiac (27%) risk genes. Eighty-six (54%) variants reached complete five-category (i.e., P, LP, VUS, LB, B) concordance, and 17 (11%) had a discordance that could affect clinical recommendations (P/LP versus VUS/LB/B). 21% and 63% of variants submitted as P and LP, respectively, were discordant with VUS. Of the 54 originally discordant variants that underwent further review, 32 reached agreement, for a post-review concordance rate of 84% (118/140 variants). This project provides an updated estimate of variant concordance, identifies considerations for LP classified variants, and highlights ongoing sources of discordance. Continued and increased sharing of variant classifications and evidence across laboratories, and the ongoing work of ClinGen to provide general as well as gene- and disease-specific guidance, will lead to continued increases in concordance.

AptCompare: optimized de novo motif discovery of RNA aptamers via HTS-SELEX

SUMMARY

High-throughput sequencing can enhance the analysis of aptamer libraries generated by the Systematic Evolution of Ligands by EXponential enrichment. Robust analysis of the resulting sequenced rounds is best implemented by determining a ranked consensus of reads following the processing by multiple aptamer detection algorithms. While several such approaches have been developed to this end, their installation and implementation is problematic. We developed AptCompare, a cross-platform program that combines six of the most widely used analytical approaches for the identification of RNA aptamer motifs and uses a simple weighted ranking to order the candidate aptamers, all driven within the same GUI-enabled environment. We demonstrate AptCompare's performance by identifying the top-ranked candidate aptamers from a previously published selection experiment in our laboratory, with follow-up bench assays demonstrating good correspondence between the sequences' rankings and their binding affinities.

AVAILABILITY AND IMPLEMENTATION

The source code and pre-built virtual machine images are freely available at https://bitbucket.org/shiehk/aptcompare.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of a novel subgroup of endometrial cancer patients with loss of thyroid hormone receptor beta expression and improved survival

BACKGROUND

Endometrial cancer (EC) is the most common gynecologic cancer in women, and the incidence of EC has increased by about 1% per year in the U. S over the last 10 years. Although 5-year survival rates for early-stage EC are around 80%, certain subtypes of EC that lose nuclear hormone receptor (NHR) expression are associated with poor survival rates. For example, estrogen receptor (ER)-negative EC typically harbors a worse prognosis compared to ER-positive EC. The molecular basis for the loss of NHR expression in endometrial tumors and its contribution to poor survival is largely unknown. Furthermore, there are no tools to systematically identify tumors that lose NHR mRNA expression relative to normal tissue. The development of such an approach could identify sets of NHR-based biomarkers for classifying patients into subgroups with poor survival outcomes.

METHODS

Here, a new computational method, termed receptLoss, was developed for identifying NHR expression loss in endometrial cancer relative to adjacent normal tissue. When applied to gene expression data from The Cancer Genome Atlas (TCGA), receptLoss identified 6 NHRs that were highly expressed in normal tissue and exhibited expression loss in a subset of endometrial tumors.

RESULTS

Three of the six identified NHRs - estrogen, progesterone, and androgen receptors - that are known to lose expression in ECs were correctly identified by receptLoss. Additionally, a novel association was found between thyroid hormone receptor beta (THRB) expression loss, increased expression of miRNA-146a, and increased rates of 5-year survival in the EC TCGA patient cohort. THRB expression loss occurs independently of estrogen and progesterone expression loss, suggesting the discovery of a distinct, clinically-relevant molecular subgroup.

CONCLUSION

ReceptLoss is a novel, open-source software tool to systematically identify NHR expression loss in cancer. The application of receptLoss to endometrial cancer gene expression data identified THRB, a previously undescribed biomarker of survival in endometrial cancer. Applying receptLoss to expression data from additional cancer types could lead to the development of biomarkers of disease progression for patients with any other tumor type. ReceptLoss can be applied to expression data from additional cancer types with the goal of identifying biomarkers of differential survival.

Quantitative Kinetic Analyses of Histone Turnover Using Imaging and Flow Cytometry

Dynamic histone changes occur as a central part of chromatin regulation. Deposition of histone variants and post-translational modifications of histones are strongly associated with properties of chromatin status. Characterizing the kinetics of histone variants allows important insights into transcription regulation, chromatin maintenance and other chromatin properties. Here we provide a protocol of quantitative and sensitive approaches to test the timing of incorporation and dissociation of histones using a two-color SNAP-labeling system, labelling pre-existing and newly-incorporated histones distinctly. Together with cell cycle synchronization methods and cell cycle markers, this approach enables a pulse-chase analysis to determine the turnover of histone variants during the cell cycle, detected using imaging or flow cytometry methods at single cell resolution. As well as testing global histone turnover, cell cycle-dependent cellular localization of histone variants can be also addressed using imaging approaches.

Functional Genomics of the Pediatric Obese Asthma Phenotype Reveal Enrichment of Rho-GTPase Pathways

Rationale: Obesity-related asthma disproportionately affects minority children and is associated with nonatopic T-helper type 1 (Th1) cell polarized inflammation that correlates with pulmonary function deficits. Its underlying mechanisms are poorly understood.Objectives: To use functional genomics to identify cellular mechanisms associated with nonatopic inflammation in obese minority children with asthma.Methods: CD4⁺ (cluster of differentiation 4-positive) Th cells from 59 obese Hispanic and African American children with asthma and 61 normal-weight Hispanic and African American children with asthma underwent quantification of the transcriptome and DNA methylome and genotyping. Expression and methylation quantitative trait loci revealed the contribution of genetic variation to transcription and DNA methylation. Adjusting for Th-cell subtype proportions discriminated loci where transcription or methylation differences were driven by differences in subtype proportions from loci that were independently associated with obesity-related asthma.Measurements and Main Results: Obese children with asthma had more memory and fewer naive Th cells than normal-weight children with asthma. Differentially expressed and methylated genes and methylation quantitative trait loci in obese children with asthma, independent of Th-cell subtype proportions, were enriched in Rho-GTPase pathways. Inhibition of CDC42 (cell division cycle 42), one of the Rho-GTPases associated with Th-cell differentiation, was associated with downregulation of the IFNγ, but not the IL-4, gene. Differential expression of the RPS27L (40S ribosomal protein S27-like) gene, part of the p53/mammalian target of rapamycin pathway, was due to nonrandom distribution of expression quantitative trait loci variants between groups. Differentially expressed and/or methylated genes, including RPS27L, were associated with pulmonary function deficits in obese children with asthma.Conclusions: We found enrichment of Rho-GTPase pathways in obese asthmatic Th cells, identifying them as a novel therapeutic target for obesity-related asthma, a disease that is suboptimally responsive to current therapies.

A Cellular Stress Response Induced by the CRISPR-dCas9 Activation System Is Not Heritable Through Cell Divisions

The CRISPR-Cas9 system can be modified to perform "epigenetic editing" by utilizing the catalytically inactive (dead) Cas9 (dCas9) to recruit regulatory proteins to specific genomic locations. In prior studies, epigenetic editing with multimers of the transactivator VP16 and guide RNAs (gRNAs) was found to cause adverse cellular responses. These side effects may confound studies inducing new cellular properties, especially if the cellular responses are maintained through cell divisions-an epigenetic regulatory property. Here, we show how distinct components of this CRISPR-dCas9 activation system, particularly dCas9 with untargeted gRNAs, upregulate genes associated with transcriptional stress, defense response, and regulation of cell death. Our results highlight a previously undetected acute stress response to CRISPR-dCas9 components in human cells, which is transient and not maintained through cell divisions.

High-efficiency genomic editing in Epstein-Barr virus-transformed lymphoblastoid B cells using a single-stranded donor oligonucleotide strategy

While human lymphoblastoid cell lines represent a valuable resource for population genetic studies, they have usually been regarded as difficult for CRISPR-mediated genomic editing because of very inefficient DNA transfection and retroviral or lentiviral transduction in these cells, which becomes a substantial problem when multiple constructs need to be co-expressed. Here we describe a protocol using a single-stranded donor oligonucleotide strategy for 'scarless' editing in lymphoblastoid cells, yielding 12/60 (20%) of clones with homology-directed recombination, when rates of <5-10% are frequently typical for many other cell types. The protocol does not require the use of lentiviruses or stable transfection, permitting lymphoblastoid cell lines to be used for CRISPR-mediated genomic targeting and screening in population genetic studies.

Functional genetic variants can mediate their regulatory effects through alteration of transcription factor binding

Functional variants in the genome are usually identified by their association with local gene expression, DNA methylation or chromatin states. DNA sequence motif analysis and chromatin immunoprecipitation studies have provided indirect support for the hypothesis that functional variants alter transcription factor binding to exert their effects. In this study, we provide direct evidence that functional variants can alter transcription factor binding. We identify a multifunctional variant within the TBC1D4 gene encoding a canonical NFκB binding site, and edited it using CRISPR-Cas9 to remove this site. We show that this editing reduces TBC1D4 expression, local chromatin accessibility and binding of the p65 component of NFκB. We then used CRISPR without genomic editing to guide p65 back to the edited locus, demonstrating that this re-targeting, occurring ~182 kb from the gene promoter, is enough to restore the function of the locus, supporting the central role of transcription factors mediating the effects of functional variants.

Misidentification of MLL3 and other mutations in cancer due to highly homologous genomic regions

The MLL3 gene has been shown to be recurrently mutated in many malignancies including in families with acute myeloid leukemia. We demonstrate that many MLL3 variant calls made by exome sequencing are false positives due to misalignment to homologous regions, including a region on chr21, and can only be validated by long-range PCR. Numerous other recurrently mutated genes reported in COSMIC and TCGA databases have pseudogenes and cannot also be validated by conventional short read-based sequencing approaches. Genome-wide identification of pseudogene regions demonstrates that frequency of these homologous regions is increased with sequencing read lengths below 200 bps. To enable identification of poor quality sequencing variants in prospective studies, we generated novel genome-wide maps of regions with poor mappability that can be used in variant calling algorithms. Taken together, our findings reveal that pseudogene regions are a source of false-positive mutations in cancers.

Chromatin organization in the female mouse brain fluctuates across the oestrous cycle

Male and female brains differ significantly in both health and disease, and yet the female brain has been understudied. Sex-hormone fluctuations make the female brain particularly dynamic and are likely to confer female-specific risks for neuropsychiatric disorders. The molecular mechanisms underlying the dynamic nature of the female brain structure and function are unknown. Here we show that neuronal chromatin organization in the female ventral hippocampus of mouse fluctuates with the oestrous cycle. We find chromatin organizational changes associated with the transcriptional activity of genes important for neuronal function and behaviour. We link these chromatin dynamics to variation in anxiety-related behaviour and brain structure. Our findings implicate an immediate-early gene product, Egr1, as part of the mechanism mediating oestrous cycle-dependent chromatin and transcriptional changes. This study reveals extreme, sex-specific dynamism of the neuronal epigenome, and establishes a foundation for the development of sex-specific treatments for disorders such as anxiety and depression.

The molecular mechanisms underlying the dynamic nature of the female brain structure and function remain poorly understood. Here the authors characterise chromatin organization in the mouse female ventral hippocampus, finding it fluctuates with the oestrous cycle, and identify changes in chromatin organization associated with the transcription of genes important for neuronal function and behaviour.

Aneuvis: web-based exploration of numerical chromosomal variation in single cells

Background

Numerical chromosomal variation is a hallmark of populations of malignant cells. Identifying the factors that promote numerical chromosomal variation is important for understanding mechanisms of carcinogenesis. However, the ability to quantify and visualize differences in chromosome number between experimentally-defined groups (e.g. control vs treated) obtained from single-cell experiments is currently limited by the lack of user-friendly software.

Results

Aneuvis is a web application that allows users to determine whether numerical chromosomal variation exists between experimental treatment groups. The web interface allows users to upload molecular cytogenetic or processed single cell whole-genome sequencing data in a cell-by-chromosome matrix format and automatically generates visualizations and summary statistics that reflect the degree of numeric chromosomal variability.

Conclusions

Aneuvis is the first user-friendly web application to help researchers identify the genetic and environmental perturbations that promote numerical chromosomal variation. Aneuvis is freely available as a web application at https://dpique.shinyapps.io/aneuvis/ and the source code for the application is available at https://github.com/dpique/aneuvis.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2842-1) contains supplementary material, which is available to authorized users.

Retargeting of macroH2A following mitosis to cytogenetic-scale heterochromatic domains

How macroH2A, a histone variant involved in silencing gene expression, is inherited from parent to daughter cells is unclear. Using a combination of imaging, biochemical, and genomic approaches, Sato et al. describe how newly synthesized macroH2A is incorporated predominantly in the G1 phase of human mitosis, targeting heterochromatic regions.

The heritability of chromatin states through cell division is a potential contributor to the epigenetic maintenance of cellular memory of prior states. The macroH2A histone variant has properties of a regulator of epigenetic cell memory, including roles controlling gene silencing and cell differentiation. Its mechanisms of regional genomic targeting and maintenance through cell division are unknown. Here, we combined in vivo imaging with biochemical and genomic approaches to show that human macroH2A is incorporated into chromatin in the G1 phase of the cell cycle following DNA replication. The newly incorporated macroH2A retargets the same large heterochromatic domains where macroH2A was already enriched in the previous cell cycle. It remains heterotypic, targeting individual nucleosomes that do not already contain a macroH2A molecule. The pattern observed resembles that of a new deposition of centromeric histone variants during the cell cycle, indicating mechanistic similarities for macrodomain-scale regulation of epigenetic properties of the cell.

Insights from deconvolution of cell subtype proportions enhance the interpretation of functional genomic data

Cell subtype proportion variability between samples contributes significantly to the variation of functional genomic properties such as gene expression or DNA methylation. Although the impact of the variation of cell subtype composition on measured genomic quantities is recognized, and some innovative tools have been developed for the analysis of heterogeneous samples, most functional genomics studies using samples with mixed cell types still ignore the influence of cell subtype proportion variation, or just deal with it as a nuisance variable to be eliminated. Here we demonstrate how harvesting information about cell subtype proportions from functional genomics data can provide insights into cellular changes associated with phenotypes. We focused on two types of mixed cell populations, human blood and mouse kidney. Cell type prediction is well developed in the former, but not currently in the latter. Estimating the cellular repertoire is easier when a reference dataset from purified samples of all cell types in the tissue is available, as is the case for blood. However, reference datasets are not available for most other tissues, such as the kidney. In this study, we showed that the proportion of alterations attributable to changes in the cellular composition varies strikingly in the two disorders (asthma and systemic lupus erythematosus), suggesting that the contribution of cell subtype proportion changes to functional genomic properties can be disease-specific. We also showed that a reference dataset from a single-cell RNA-seq study successfully estimated the cell subtype proportions in mouse kidney and allowed us to distinguish altered cell subtype differences between two different knock-out mouse models, both of which had reported a reduced number of glomeruli compared to their wild-type counterparts. These findings demonstrate that testing for changes in cell subtype proportions between conditions can yield important insights in functional genomics studies.

Engineering a haematopoietic stem cell niche by revitalizing mesenchymal stromal cells

Haematopoietic stem cells (HSCs) are maintained by bone marrow (BM) niches in vivo^1,2, but the ability of niche cells to maintain HSCs ex vivo is markedly diminished. Expression of niche factors by Nestin-GFP⁺ mesenchymal-derived stromal cells (MSCs) is downregulated upon culture, suggesting that transcriptional rewiring may contribute to this reduced HSC maintenance potential. Using an RNA sequencing screen, we identified 5 transcription factors (Klf7, Ostf1, Xbp1, Irf3, Irf7) that restored HSC niche function in cultured BM-derived MSCs. These revitalized MSCs (rMSCs) exhibited enhanced synthesis of HSC niche factors while retaining their mesenchymal differentiation capacity. In contrast to HSCs co-cultured with control MSCs, HSCs expanded with rMSCs showed higher repopulation capacity and protected lethally irradiated recipient mice. Competitive reconstitution assays revealed ~7-fold expansion of functional HSCs by rMSCs. rMSCs prevented the accumulation of DNA damage in cultured HSCs, a hallmark of ageing and replication stress. Analysis of the reprogramming mechanisms uncovered a role for myocyte enhancer factor 2c (Mef2c) in the revitalization of MSCs. These results provide insight in the transcriptional regulation of the niche with implications for stem cell-based therapies.