Despite the initial benefits of treating HER2-amplified breast cancer patients with the tyrosine kinase inhibitor lapatinib, resistance inevitably develops. Here we report that lapatinib induces the degradation of the nuclear receptor ERRα, a master regulator of cellular metabolism, and that the expression of ERRα is restored in lapatinib-resistant breast cancer cells through reactivation of mTOR signalling. Re-expression of ERRα in resistant cells triggers metabolic adaptations favouring mitochondrial energy metabolism through increased glutamine metabolism, as well as ROS detoxification required for cell survival under therapeutic stress conditions. An ERRα inverse agonist counteracts these metabolic adaptations and overcomes lapatinib resistance in a HER2-induced mammary tumour mouse model. This work reveals a molecular mechanism by which ERRα-induced metabolic reprogramming promotes survival of lapatinib-resistant cancer cells and demonstrates the potential of ERRα inhibition as an effective adjuvant therapy in poor outcome HER2-positive breast cancer.

Genome-wide demethylation and remethylation of DNA during early embryogenesis is essential for development. Imprinted germline differentially methylated domains (gDMDs) established by sex-specific methylation in either male or female germ cells, must escape these dynamic changes and sustain precise inheritance of both methylated and unmethylated parental alleles. To identify other, gDMD-like sequences with the same epigenetic inheritance properties, we used a modified embryonic stem (ES) cell line that emulates the early embryonic demethylation and remethylation waves. Transient DNMT1 suppression revealed gDMD-like sequences requiring continuous DNMT1 activity to sustain a highly methylated state. Remethylation of these sequences was also compromised in vivo in a mouse model of transient DNMT1 loss in the preimplantation embryo. These novel regions, possessing heritable epigenetic features similar to imprinted-gDMDs are required for normal physiological and developmental processes and when disrupted are associated with disorders such as cancer and autism spectrum disorders. This study presents new perspectives on DNA methylation heritability during early embryo development that extend beyond conventional imprinted-gDMDs.

We present a new update to MetaboAnalyst (version 4.0) for comprehensive metabolomic data analysis, interpretation, and integration with other omics data. Since the last major update in 2015, MetaboAnalyst has continued to evolve based on user feedback and technological advancements in the field. For this year's update, four new key features have been added to MetaboAnalyst 4.0, including: (1) real-time R command tracking and display coupled with the release of a companion MetaboAnalystR package; (2) a MS Peaks to Pathways module for prediction of pathway activity from untargeted mass spectral data using the mummichog algorithm; (3) a Biomarker Meta-analysis module for robust biomarker identification through the combination of multiple metabolomic datasets and (4) a Network Explorer module for integrative analysis of metabolomics, metagenomics, and/or transcriptomics data. The user interface of MetaboAnalyst 4.0 has been reengineered to provide a more modern look and feel, as well as to give more space and flexibility to introduce new functions. The underlying knowledgebases (compound libraries, metabolite sets, and metabolic pathways) have also been updated based on the latest data from the Human Metabolome Database (HMDB). A Docker image of MetaboAnalyst is also available to facilitate download and local installation of MetaboAnalyst. MetaboAnalyst 4.0 is freely available at http://metaboanalyst.ca.

Genome-wide comparisons of transcription factor binding sites in different species can be used to evaluate evolutionary constraints that shape gene regulatory circuits and to understand how the interaction between transcription factors shapes their binding landscapes over evolution.

Structural variations (SVs) contribute significantly to the variability of the human genome and extensive genomic rearrangements are a hallmark of cancer. While genomic DNA paired-end-tag (DNA-PET) sequencing is an attractive approach to identify genomic SVs, the current application of PET sequencing with short insert size DNA can be insufficient for the comprehensive mapping of SVs in low complexity and repeat-rich genomic regions. We employed a recently developed procedure to generate PET sequencing data using large DNA inserts of 10-20 kb and compared their characteristics with short insert (1 kb) libraries for their ability to identify SVs. Our results suggest that although short insert libraries bear an advantage in identifying small deletions, they do not provide significantly better breakpoint resolution. In contrast, large inserts are superior to short inserts in providing higher physical genome coverage for the same sequencing cost and achieve greater sensitivity, in practice, for the identification of several classes of SVs, such as copy number neutral and complex events. Furthermore, our results confirm that large insert libraries allow for the identification of SVs within repetitive sequences, which cannot be spanned by short inserts. This provides a key advantage in studying rearrangements in cancer, and we show how it can be used in a fusion-point-guided-concatenation algorithm to study focally amplified regions in cancer.

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor alpha (ER-alpha) in the human genome. We found that most high-confidence remote ER-alpha-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-alpha functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Nephron number varies widely between 0.3 and 1.3 million per kidney in humans. During fetal life, the rate of nephrogenesis is influenced by local retinoic acid (RA) level such that even moderate maternal vitamin A deficiency limits the final nephron number in rodents. Inactivation of genes in the RA pathway causes renal agenesis in mice; however, the impact of retinoids on human kidney development is unknown. To resolve this, we tested for associations between variants of genes involved in RA metabolism (ALDH1A2, CYP26A1, and CYP26B1) and kidney size among normal newborns. Homozygosity for a common (1 in 5) variant, rs7169289(G), within an Sp1 transcription factor motif of the ALDH1A2 gene, showed a significant 22% increase in newborn kidney volume when adjusted for body surface area. Infants bearing this allele had higher umbilical cord blood RA levels compared to those with homozygous wild-type ALDH1A2 rs7169289(A) alleles. Furthermore, the effect of the rs7169289(G) variant was evident in subgroups with or without a previously reported hypomorphic RET 1476(A) proto-oncogene allele that is critical in determining final nephron number. As maternal vitamin A deficiency is widespread in developing countries and may compromise availability of retinol for fetal RA synthesis, our study suggests that the ALDH1A2 rs7169289(G) variant might be protective for such individuals.

Familial clustering of colorectal cancer occurs in 15-20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer.

Epilepsy will affect nearly 3% of people at some point during their lifetime. Previous copy number variants (CNVs) studies of epilepsy have used array-based technology and were restricted to the detection of large or exonic events. In contrast, whole-genome sequencing (WGS) has the potential to more comprehensively profile CNVs but existing analytic methods suffer from limited accuracy. We show that this is in part due to the non-uniformity of read coverage, even after intra-sample normalization. To improve on this, we developed PopSV, an algorithm that uses multiple samples to control for technical variation and enables the robust detection of CNVs. Using WGS and PopSV, we performed a comprehensive characterization of CNVs in 198 individuals affected with epilepsy and 301 controls. For both large and small variants, we found an enrichment of rare exonic events in epilepsy patients, especially in genes with predicted loss-of-function intolerance. Notably, this genome-wide survey also revealed an enrichment of rare non-coding CNVs near previously known epilepsy genes. This enrichment was strongest for non-coding CNVs located within 100 Kbp of an epilepsy gene and in regions associated with changes in the gene expression, such as expression QTLs or DNase I hypersensitive sites. Finally, we report on 21 potentially damaging events that could be associated with known or new candidate epilepsy genes. Our results suggest that comprehensive sequence-based profiling of CNVs could help explain a larger fraction of epilepsy cases.

Childhood brain tumors have suspected prenatal origins. To identify vulnerable developmental states, we generated a single-cell transcriptome atlas of >65,000 cells from embryonal pons and forebrain, two major tumor locations. We derived signatures for 191 distinct cell populations and defined the regional cellular diversity and differentiation dynamics. Projection of bulk tumor transcriptomes onto this dataset shows that WNT medulloblastomas match the rhombic lip-derived mossy fiber neuronal lineage and embryonal tumors with multilayered rosettes fully recapitulate a neuronal lineage, while group 2a/b atypical teratoid/rhabdoid tumors may originate outside the neuroectoderm. Importantly, single-cell tumor profiles reveal highly defined cell hierarchies that mirror transcriptional programs of the corresponding normal lineages. Our findings identify impaired differentiation of specific neural progenitors as a common mechanism underlying these pediatric cancers and provide a rational framework for future modeling and therapeutic interventions.

The paternal environment has been linked to infertility and negative outcomes. Such effects may be transmitted via sperm through histone modifications. To date, in-depth profiling of the sperm chromatin in men has been limited. Here, we use deep sequencing to characterize the sperm profiles of histone H3 lysine 4 tri-methylation (H3K4me3) and DNA methylation in a representative reference population of 37 men. Our analysis reveals that H3K4me3 is localized throughout the genome and at genes for fertility and development. Remarkably, enrichment is also found at regions that escape epigenetic reprogramming in primordial germ cells, embryonic enhancers, and short-interspersed nuclear elements (SINEs). There is significant overlap in H3K4me3 and DNA methylation throughout the genome, suggesting a potential interplay between these marks previously reported to be mutually exclusive in sperm. Comparisons made between H3K4me3 marked regions in sperm and the embryonic transcriptome suggest an influence of paternal chromatin on embryonic gene expression.

The functional effect of a gene edit by designer nucleases depends on the DNA repair outcome at the targeted locus. While non-homologous end joining (NHEJ) repair results in various mutations, microhomology-mediated end joining (MMEJ) creates precise deletions based on the alignment of flanking microhomologies (µHs). Recently, the sequence context surrounding nuclease-induced double strand breaks (DSBs) has been shown to predict repair outcomes, for which µH plays an important role. Here, we survey naturally occurring human deletion variants and identify that 11 million or 57% are flanked by µHs, covering 88% of protein-coding genes. These biologically relevant mutations are candidates for precise creation in a template-free manner by MMEJ repair. Using CRISPR-Cas9 in human induced pluripotent stem cells (hiPSCs), we efficiently create pathogenic deletion mutations for demonstrable disease models with both gain- and loss-of-function phenotypes. We anticipate this dataset and gene editing strategy to enable functional genetic studies and drug screening.

The sperm DNA methylation landscape is unique and critical for offspring health. If gamete-derived DNA methylation escapes reprograming in early embryos, epigenetic defects in sperm may be transmitted to the next generation. Current techniques to assess sperm DNA methylation show bias toward CpG-dense regions and do not target areas of dynamic methylation, those predicted to be environmentally sensitive and tunable regulatory elements.

Cancer stem cells are critical for cancer initiation, development, and treatment resistance. Our understanding of these processes, and how they relate to glioblastoma heterogeneity, is limited. To overcome these limitations, we performed single-cell RNA sequencing on 53586 adult glioblastoma cells and 22637 normal human fetal brain cells, and compared the lineage hierarchy of the developing human brain to the transcriptome of cancer cells. We find a conserved neural tri-lineage cancer hierarchy centered around glial progenitor-like cells. We also find that this progenitor population contains the majority of the cancer's cycling cells, and, using RNA velocity, is often the originator of the other cell types. Finally, we show that this hierarchal map can be used to identify therapeutic targets specific to progenitor cancer stem cells. Our analyses show that normal brain development reconciles glioblastoma development, suggests a possible origin for glioblastoma hierarchy, and helps to identify cancer stem cell-specific targets.

La-related protein 1 (LARP1) has been identified as a key translational inhibitor of terminal oligopyrimidine (TOP) mRNAs downstream of the nutrient sensing protein kinase complex, mTORC1. LARP1 exerts this inhibitory effect on TOP mRNA translation by binding to the mRNA cap and the adjacent 5'TOP motif, resulting in the displacement of the cap-binding protein eIF4E from TOP mRNAs. However, the involvement of additional signaling pathway in regulating LARP1-mediated inhibition of TOP mRNA translation is largely unexplored. In the present study, we identify a second nutrient sensing kinase GCN2 that converges on LARP1 to control TOP mRNA translation. Using chromatin-immunoprecipitation followed by massive parallel sequencing (ChIP-seq) analysis of activating transcription factor 4 (ATF4), an effector of GCN2 in nutrient stress conditions, in WT and GCN2 KO mouse embryonic fibroblasts, we determined that LARP1 is a GCN2-dependent transcriptional target of ATF4. Moreover, we identified GCN1, a GCN2 activator, participates in a complex with LARP1 on stalled ribosomes, suggesting a role for GCN1 in LARP1-mediated translation inhibition in response to ribosome stalling. Therefore, our data suggest that the GCN2 pathway controls LARP1 activity via two mechanisms: ATF4-dependent transcriptional induction of LARP1 mRNA and GCN1-mediated recruitment of LARP1 to stalled ribosomes.

Despite the availability of highly efficacious vaccines, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lacks effective drug treatment, which results in a high rate of mortality. To address this therapeutic shortcoming, we applied a systems biology approach to the study of patients hospitalized with severe COVID. We show that, at the time of hospital admission, patients who were equivalent on the clinical ordinal scale displayed significant differential monocyte epigenetic and transcriptomic attributes between those who would survive and those who would succumb to COVID-19. We identified messenger RNA metabolism, RNA splicing, and interferon signaling pathways as key host responses overactivated by patients who would not survive. Those pathways are prime drug targets to reduce mortality of critically ill patients with COVID-19, leading us to identify tacrolimus, zotatifin, and nintedanib as three strong candidates for treatment of severely ill patients at the time of hospital admission.

Québec was the Canadian province most impacted by COVID-19, with 401,462 cases as of September 24th, 2021, and 11,347 deaths due mostly to a very severe first pandemic wave. In April 2020, we assembled the Coronavirus Sequencing in Québec (CoVSeQ) consortium to sequence SARS-CoV-2 genomes in Québec to track viral introduction events and transmission within the province.

Here, we exploit a deep serological profiling strategy coupled with an integrated, computational framework for the analysis of SARS-CoV-2 humoral immune responses. Applying a high-density peptide array (HDPA) spanning the entire proteomes of SARS-CoV-2 and endemic human coronaviruses allowed identification of B cell epitopes and relate them to their evolutionary and structural properties. We identify hotspots of pre-existing immunity and identify cross-reactive epitopes that contribute to increasing the overall humoral immune response to SARS-CoV-2. Using a public dataset of over 38,000 viral genomes from the early phase of the pandemic, capturing both inter- and within-host genetic viral diversity, we determined the evolutionary profile of epitopes and the differences across proteins, waves, and SARS-CoV-2 variants. Lastly, we show that mutations in spike and nucleocapsid epitopes are under stronger selection between than within patients, suggesting that most of the selective pressure for immune evasion occurs upon transmission between hosts.

ZMYM2 is a transcriptional repressor whose role in development is largely unexplored. We found that Zmym2-/- mice show embryonic lethality by E10.5. Molecular characterization of Zmym2-/- embryos revealed two distinct defects. First, they fail to undergo DNA methylation and silencing of germline gene promoters, resulting in widespread upregulation of germline genes. Second, they fail to methylate and silence the evolutionarily youngest and most active LINE element subclasses in mice. Zmym2-/- embryos show ubiquitous overexpression of LINE-1 protein as well as aberrant expression of transposon-gene fusion transcripts. ZMYM2 homes to sites of PRC1.6 and TRIM28 complex binding, mediating repression of germline genes and transposons respectively. In the absence of ZMYM2, hypermethylation of histone 3 lysine 4 occurs at target sites, creating a chromatin landscape unfavourable for establishment of DNA methylation. ZMYM2-/- human embryonic stem cells also show aberrant upregulation and demethylation of young LINE elements, indicating a conserved role in repression of active transposons. ZMYM2 is thus an important new factor in DNA methylation patterning in early embryonic development.