Multiomic analysis of familial adenomatous polyposis reveals molecular pathways associated with early tumorigenesis.

Familial adenomatous polyposis (FAP) is a genetic disease causing hundreds of premalignant polyps in affected persons and is an ideal model to study transitions of early precancer states to colorectal cancer (CRC). We performed deep multiomic profiling of 93 samples, including normal mucosa, benign polyps and dysplastic polyps, from six persons with FAP. Transcriptomic, proteomic, metabolomic and lipidomic analyses revealed a dynamic choreography of thousands of molecular and cellular events that occur during precancerous transitions toward cancer formation. These involve processes such as cell proliferation, immune response, metabolic alterations (including amino acids and lipids), hormones and extracellular matrix proteins. Interestingly, activation of the arachidonic acid pathway was found to occur early in hyperplasia; this pathway is targeted by aspirin and other nonsteroidal anti-inflammatory drugs, a preventative treatment under investigation in persons with FAP. Overall, our results reveal key genomic, cellular and molecular events during the earliest steps in CRC formation and potential mechanisms of pharmaceutical prophylaxis.

Global loss of promoter-enhancer connectivity and rebalancing of gene expression during early colorectal cancer carcinogenesis.

Although three-dimensional (3D) genome architecture is crucial for gene regulation, its role in disease remains elusive. We traced the evolution and malignant transformation of colorectal cancer (CRC) by generating high-resolution chromatin conformation maps of 33 colon samples spanning different stages of early neoplastic growth in persons with familial adenomatous polyposis (FAP). Our analysis revealed a substantial progressive loss of genome-wide cis-regulatory connectivity at early malignancy stages, correlating with nonlinear gene regulation effects. Genes with high promoter-enhancer (P-E) connectivity in unaffected mucosa were not linked to elevated baseline expression but tended to be upregulated in advanced stages. Inhibiting highly connected promoters preferentially represses gene expression in CRC cells compared to normal colonic epithelial cells. Our results suggest a two-phase model whereby neoplastic transformation reduces P-E connectivity from a redundant state to a rate-limiting one for transcriptional levels, highlighting the intricate interplay between 3D genome architecture and gene regulation during early CRC progression.

Netrin G1 Ligand is a new stromal immunomodulator that promotes pancreatic cancer.

Understanding pancreatic cancer biology is fundamental for identifying new targets and for developing more effective therapies. In particular, the contribution of the stromal microenvironment to pancreatic cancer tumorigenesis requires further exploration. Here, we report the stromal roles of the synaptic protein Netrin G1 Ligand (NGL-1) in pancreatic cancer, uncovering its pro-tumor functions in cancer-associated fibroblasts and in immune cells. We observed that the stromal expression of NGL-1 inversely correlated with patients' overall survival. Moreover, germline knockout (KO) mice for NGL-1 presented decreased tumor burden, with a microenvironment that is less supportive of tumor growth. Of note, tumors from NGL-1 KO mice produced less immunosuppressive cytokines and displayed an increased percentage of CD8 + T cells than those from control mice, while preserving the physical structure of the tumor microenvironment. These effects were shown to be mediated by NGL-1 in both immune cells and in the local stroma, in a TGF-ß-dependent manner. While myeloid cells lacking NGL-1 decreased the production of immunosuppressive cytokines, NGL-1 KO T cells showed increased proliferation rates and overall polyfunctionality compared to control T cells. CAFs lacking NGL-1 were less immunosuppressive than controls, with overall decreased production of pro-tumor cytokines and compromised ability to inhibit CD8 + T cells activation. Mechanistically, these CAFs downregulated components of the TGF-ß pathway, AP-1 and NFAT transcription factor families, resulting in a less tumor-supportive phenotype. Finally, targeting NGL-1 genetically or using a functionally antagonistic small peptide phenocopied the effects of chemotherapy, while modulating the immunosuppressive tumor microenvironment (TME), rather than eliminating it. We propose NGL-1 as a new local stroma and immunomodulatory molecule, with pro-tumor roles in pancreatic cancer. Statement of Significance: Here we uncovered the pro-tumor roles of the synaptic protein NGL-1 in the tumor microenvironment of pancreatic cancer, defining a new target that simultaneously modulates tumor cell, fibroblast, and immune cell functions. This study reports a new pathway where NGL-1 controls TGF-ß, AP-1 transcription factor members and NFAT1, modulating the immunosuppressive microenvironment in pancreatic cancer. Our findings highlight NGL-1 as a new stromal immunomodulator in pancreatic cancer.

Rate of Pathogenic Germline Variants in Patients With Lung Cancer.

PURPOSE: Germline genetic testing (GGT) is now recommended for all patients diagnosed with ovarian or pancreatic cancer and for a large proportion of patients based solely on a diagnosis of colorectal or breast cancer. However, GGT is not yet recommended for all patients diagnosed with lung cancer (LC), primarily because of a lack of evidence that supports a significant frequency of identifying pathogenic germline variants (PGVs) in these patients. This study characterizes GGT results in a cohort of patients with LC. METHODS: We reviewed deidentified data for 7,788 patients with GGT (2015-2022). PGV frequencies were compared to a control cohort of unaffected individuals. GGT results were stratified by genomic ancestry, history of cancer, and PGV clinical actionability per current guidelines. RESULTS: Of all patients with LC, 14.9% (1,161/7,788) had PGVs. The rate was similar when restricted to patients with no cancer family history (FH) or personal history (PH) of other cancers (14.3%). PGVs were significantly enriched in BRCA2, ATM, CHEK2, BRCA1, and mismatch repair genes compared with controls. Patients of European (EUR) genomic ancestry had the highest PGV rate (18%) and variants of uncertain significance were significantly higher in patients of non-EUR genomic ancestry. Of the PGVs identified, 61.3% were in DNA damage repair (DDR) genes and 95% were clinically actionable. CONCLUSION: This retrospective study shows a LC diagnosis identifies patients with a significant likelihood of having a cancer-predisposing PGV across genomic ancestries. Enrichment of PGVs in DDR genes suggests that these PGVs may contribute to LC cancer predisposition. The frequency of PGVs among patients with LC did not differ significantly according to FH or PH of other cancers.

Epigenetic balance ensures mechanistic control of MLL amplification and rearrangement.

MLL/KMT2A amplifications and translocations are prevalent in infant, adult, and therapy-induced leukemia. However, the molecular contributor(s) to these alterations are unclear. Here, we demonstrate that histone H3 lysine 9 mono- and di-methylation (H3K9me1/2) balance at the MLL/KMT2A locus regulates these amplifications and rearrangements. This balance is controlled by the crosstalk between lysine demethylase KDM3B and methyltransferase G9a/EHMT2. KDM3B depletion increases H3K9me1/2 levels and reduces CTCF occupancy at the MLL/KMT2A locus, in turn promoting amplification and rearrangements. Depleting CTCF is also sufficient to generate these focal alterations. Furthermore, the chemotherapy doxorubicin (Dox), which associates with therapy-induced leukemia and promotes MLL/KMT2A amplifications and rearrangements, suppresses KDM3B and CTCF protein levels. KDM3B and CTCF overexpression rescues Dox-induced MLL/KMT2A alterations. G9a inhibition in human cells or mice also suppresses MLL/KMT2A events accompanying Dox treatment. Therefore, MLL/KMT2A amplifications and rearrangements are controlled by epigenetic regulators that are tractable drug targets, which has clinical implications.

TET1 and TDG Suppress Inflammatory Response in Intestinal Tumorigenesis: Implications for Colorectal Tumors With the CpG Island Methylator Phenotype.

BACKGROUND & AIMS: Aberrant DNA methylation is frequent in colorectal cancer (CRC), but underlying mechanisms and pathologic consequences are poorly understood. METHODS: We disrupted active DNA demethylation genes Tet1 and/or Tdg from ApcMin mice and characterized the methylome and transcriptome of colonic adenomas. Data were compared to human colonic adenocarcinomas (COAD) in The Cancer Genome Atlas. RESULTS: There were increased numbers of small intestinal adenomas in ApcMin mice expressing the TdgN151A allele, whereas Tet1-deficient and Tet1/TdgN151A-double heterozygous ApcMin colonic adenomas were larger with features of erosion and invasion. We detected reduction in global DNA hypomethylation in colonic adenomas from Tet1- and Tdg-mutant ApcMin mice and hypermethylation of CpG islands in Tet1-mutant ApcMin adenomas. Up-regulation of inflammatory, immune, and interferon response genes was present in Tet1- and Tdg-mutant colonic adenomas compared to control ApcMin adenomas. This up-regulation was also seen in murine colonic organoids and human CRC lines infected with lentiviruses expressing TET1 or TDG short hairpin RNA. A 127-gene inflammatory signature separated colonic adenocarcinomas into 4 groups, closely aligned with their microsatellite or chromosomal instability and characterized by different levels of DNA methylation and DNMT1 expression that anticorrelated with TET1 expression. Tumors with the CpG island methylator phenotype (CIMP) had concerted high DNMT1/low TET1 expression. TET1 or TDG knockdown in CRC lines enhanced killing by natural killer cells. CONCLUSIONS: Our findings reveal a novel epigenetic regulation, linked to the type of genomic instability, by which TET1/TDG-mediated DNA demethylation decreases methylation levels and inflammatory/interferon/immune responses. CIMP in CRC is triggered by an imbalance of methylating activities over demethylating activities. These mice represent a model of CIMP CRC.

Longitudinal linked-read sequencing reveals ecological and evolutionary responses of a human gut microbiome during antibiotic treatment.

Gut microbial communities can respond to antibiotic perturbations by rapidly altering their taxonomic and functional composition. However, little is known about the strain-level processes that drive this collective response. Here, we characterize the gut microbiome of a single individual at high temporal and genetic resolution through a period of health, disease, antibiotic treatment, and recovery. We used deep, linked-read metagenomic sequencing to track the longitudinal trajectories of thousands of single nucleotide variants within 36 species, which allowed us to contrast these genetic dynamics with the ecological fluctuations at the species level. We found that antibiotics can drive rapid shifts in the genetic composition of individual species, often involving incomplete genome-wide sweeps of pre-existing variants. These genetic changes were frequently observed in species without obvious changes in species abundance, emphasizing the importance of monitoring diversity below the species level. We also found that many sweeping variants quickly reverted to their baseline levels once antibiotic treatment had concluded, demonstrating that the ecological resilience of the microbiota can sometimes extend all the way down to the genetic level. Our results provide new insights into the population genetic forces that shape individual microbiomes on therapeutically relevant timescales, with potential implications for personalized health and disease.

Candidate variants in TUB are associated with familial tremor.

Essential tremor (ET) is the most common adult-onset movement disorder. In the present study, we performed whole exome sequencing of a large ET-affected family (10 affected and 6 un-affected family members) and identified a TUB p.V431I variant (rs75594955) segregating in a manner consistent with autosomal-dominant inheritance. Subsequent targeted re-sequencing of TUB in 820 unrelated individuals with sporadic ET and 630 controls revealed significant enrichment of rare nonsynonymous TUB variants (e.g. rs75594955: p.V431I, rs1241709665: p.Ile20Phe, rs55648406: p.Arg49Gln) in the ET cohort (SKAT-O test p-value = 6.20e-08). TUB encodes a transcription factor predominantly expressed in neuronal cells and has been previously implicated in obesity. ChIP-seq analyses of the TUB transcription factor across different regions of the mouse brain revealed that TUB regulates the pathways responsible for neurotransmitter production as well thyroid hormone signaling. Together, these results support the association of rare variants in TUB with ET.

Chromosome-level de novo assembly of the pig-tailed macaque genome using linked-read sequencing and HiC proximity scaffolding.

BACKGROUND: Macaque species share >93% genome homology with humans and develop many disease phenotypes similar to those of humans, making them valuable animal models for the study of human diseases (e.g., HIV and neurodegenerative diseases). However, the quality of genome assembly and annotation for several macaque species lags behind the human genome effort. RESULTS: To close this gap and enhance functional genomics approaches, we used a combination of de novo linked-read assembly and scaffolding using proximity ligation assay (HiC) to assemble the pig-tailed macaque (Macaca nemestrina) genome. This combinatorial method yielded large scaffolds at chromosome level with a scaffold N50 of 127.5 Mb; the 23 largest scaffolds covered 90% of the entire genome. This assembly revealed large-scale rearrangements between pig-tailed macaque chromosomes 7, 12, and 13 and human chromosomes 2, 14, and 15. We subsequently annotated the genome using transcriptome and proteomics data from personalized induced pluripotent stem cells derived from the same animal. Reconstruction of the evolutionary tree using whole-genome annotation and orthologous comparisons among 3 macaque species, human, and mouse genomes revealed extensive homology between human and pig-tailed macaques with regards to both pluripotent stem cell genes and innate immune gene pathways. Our results confirm that rhesus and cynomolgus macaques exhibit a closer evolutionary distance to each other than either species exhibits to humans or pig-tailed macaques. CONCLUSIONS: These findings demonstrate that pig-tailed macaques can serve as an excellent animal model for the study of many human diseases particularly with regards to pluripotency and innate immune pathways.

Multi-faceted epigenetic dysregulation of gene expression promotes esophageal squamous cell carcinoma.

Epigenetic landscapes can shape physiologic and disease phenotypes. We used integrative, high resolution multi-omics methods to delineate the methylome landscape and characterize the oncogenic drivers of esophageal squamous cell carcinoma (ESCC). We found 98% of CpGs are hypomethylated across the ESCC genome. Hypo-methylated regions are enriched in areas with heterochromatin binding markers (H3K9me3, H3K27me3), while hyper-methylated regions are enriched in polycomb repressive complex (EZH2/SUZ12) recognizing regions. Altered methylation in promoters, enhancers, and gene bodies, as well as in polycomb repressive complex occupancy and CTCF binding sites are associated with cancer-specific gene dysregulation. Epigenetic-mediated activation of non-canonical WNT/ß-catenin/MMP signaling and a YY1/lncRNA ESCCAL-1/ribosomal protein network are uncovered and validated as potential novel ESCC driver alterations. This study advances our understanding of how epigenetic landscapes shape cancer pathogenesis and provides a resource for biomarker and target discovery.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop.

BACKGROUND: Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing. RESULTS: Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences. CONCLUSIONS: This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

Lifelong physical activity is associated with promoter hypomethylation of genes involved in metabolism, myogenesis, contractile properties and oxidative stress resistance in aged human skeletal muscle.

Lifelong regular physical activity is associated with reduced risk of type 2 diabetes (T2D), maintenance of muscle mass and increased metabolic capacity. However, little is known about epigenetic mechanisms that might contribute to these beneficial effects in aged individuals. We investigated the effect of lifelong physical activity on global DNA methylation patterns in skeletal muscle of healthy aged men, who had either performed regular exercise or remained sedentary their entire lives (average age 62 years). DNA methylation was significantly lower in 714 promoters of the physically active than inactive men while methylation of introns, exons and CpG islands was similar in the two groups. Promoters for genes encoding critical insulin-responsive enzymes in glycogen metabolism, glycolysis and TCA cycle were hypomethylated in active relative to inactive men. Hypomethylation was also found in promoters of myosin light chain, dystrophin, actin polymerization, PAK regulatory genes and oxidative stress response genes. A cluster of genes regulated by GSK3ß-TCF7L2 also displayed promoter hypomethylation. Together, our results suggest that lifelong physical activity is associated with DNA methylation patterns that potentially allow for increased insulin sensitivity and a higher expression of genes in energy metabolism, myogenesis, contractile properties and oxidative stress resistance in skeletal muscle of aged individuals.

Piercing the dark matter: bioinformatics of long-range sequencing and mapping.

Several new genomics technologies have become available that offer long-read sequencing or long-range mapping with higher throughput and higher resolution analysis than ever before. These long-range technologies are rapidly advancing the field with improved reference genomes, more comprehensive variant identification and more complete views of transcriptomes and epigenomes. However, they also require new bioinformatics approaches to take full advantage of their unique characteristics while overcoming their complex errors and modalities. Here, we discuss several of the most important applications of the new technologies, focusing on both the currently available bioinformatics tools and opportunities for future research.

Hybrid assembly with long and short reads improves discovery of gene family expansions.

BACKGROUND: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. METHODS: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. RESULTS: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. CONCLUSION: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

The pineapple genome and the evolution of CAM photosynthesis.

Pineapple (Ananas comosus (L.) Merr.) is the most economically valuable crop possessing crassulacean acid metabolism (CAM), a photosynthetic carbon assimilation pathway with high water-use efficiency, and the second most important tropical fruit. We sequenced the genomes of pineapple varieties F153 and MD2 and a wild pineapple relative, Ananas bracteatus accession CB5. The pineapple genome has one fewer ancient whole-genome duplication event than sequenced grass genomes and a conserved karyotype with seven chromosomes from before the ρ duplication event. The pineapple lineage has transitioned from C3 photosynthesis to CAM, with CAM-related genes exhibiting a diel expression pattern in photosynthetic tissues. CAM pathway genes were enriched with cis-regulatory elements associated with the regulation of circadian clock genes, providing the first cis-regulatory link between CAM and circadian clock regulation. Pineapple CAM photosynthesis evolved by the reconfiguration of pathways in C3 plants, through the regulatory neofunctionalization of preexisting genes and not through the acquisition of neofunctionalized genes via whole-genome or tandem gene duplication.

Nasal immunization with M cell-targeting ligand-conjugated ApxIIA toxin fragment induces protective immunity against Actinobacillus pleuropneumoniae infection in a murine model.

Actinobacillus pleuropneumoniae is the causative agent of porcine pleuropneumonia and severe economic loss in the swine industry has been caused by the infection. Therefore, the development of an effective vaccine against the bacteria is necessary. ApxII toxin, among several virulence factors expressed by the bacteria, is considered to be a promising vaccine candidate because ApxII toxin not only accompanies cytotoxic and hemolytic activities, but is also expressed in all 15 serotypes of bacteria except serotypes 10 and 14. In this study, we identified the peptide ligand capable of targeting the ligand-conjugated ApxIIA #5 fragment antigen to nasopharynx-associated lymphoid tissue. It was found that nasal immunization with ligand-conjugated ApxIIA #5 induced efficient mucosal and systemic immune responses measured at the levels of antigen-specific antibodies, cytokine-secreting cells after antigen exposure, and antigen-specific lymphocyte proliferation. More importantly, the nasal immunization induced protective immunity against nasal challenge infection of the bacteria, which was confirmed by histopathological studies and bacterial clearance after challenge infection. Collectively, we confirmed that the ligand capable of targeting the ligand-conjugated antigen to nasopharynx-associated lymphoid tissue can be used as an effective nasal vaccine adjuvant to induce protective immunity against A. pleuropneumoniae infection.

Expression of the ATP-gated P2X7 Receptor on M Cells and Its Modulating Role in the Mucosal Immune Environment.

Interactions between microbes and epithelial cells in the gastrointestinal tract are closely associated with regulation of intestinal mucosal immune responses. Recent studies have highlighted the modulation of mucosal immunity by microbe-derived molecules such as ATP and short-chain fatty acids. In this study, we undertook to characterize the expression of the ATP-gated P2X7 receptor (P2X7R) on M cells and its role in gastrointestinal mucosal immune regulation because it was poorly characterized in Peyer's patches, although purinergic signaling via P2X7R and luminal ATP have been considered to play an important role in the gastrointestinal tract. Here, we present the first report on the expression of P2X7R on M cells and characterize the role of P2X7R in immune enhancement by ATP or LL-37.

Immune-Modulating Activity of Extract Prepared from Mycelial Culture of Chinese Caterpillar Mushroom, Ophiocordyceps sinensis (Ascomycetes).

Ophiocordyceps sinensis is a natural fungus that has been valued as a health food and traditional Chinese medicine for centuries. The fungus is parasitic and colonizes insect larva. Naturally occurring O. sinensis thrives at high altitude in cold and grassy alpine meadows on the Himalayan mountain ranges. Wild O. sinensis is becoming increasingly rare in its natural habitats, and its price is out of reach for clinical practice. For these reasons, development of a standardized alternative is a great focus of research to allow the use of O. sinensis as a medicine. To develop an alternative for wild O. sinensis, a refined standardized extract, CBG-CS-2, was produced by artificial fermentation and extraction of the mycelial strain Paecilomyces hepiali CBG-CS-1, which originated from wild O. sinensis. In this study, we analyzed the in vivo immune-modulating effect of CBG-CS-2 in mice. Oral administration of CBG-CS-2 supported splenocyte stimulation and enhanced Th1-type cytokine expression from the splenocytes. Importantly, the same treatment significantly enhanced the natural killer cell activity of the splenocytes. Finally, oral administration of CBG-CS-2 enhanced the potential for inflammatory responses. Together, these findings indicate that the mycelial culture extract prepared from O. sinensis exhibited immune-modulating activity and suggest its possible use in the treatment of diseases caused by abnormal immune function.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica.

BACKGROUND: The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate. RESULTS: Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the "pan-genome" of three divergent rice varieties and document several megabases of each genome absent in the other two. CONCLUSIONS: Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

SplitMEM: a graphical algorithm for pan-genome analysis with suffix skips.

MOTIVATION: Genomics is expanding from a single reference per species paradigm into a more comprehensive pan-genome approach that analyzes multiple individuals together. A compressed de Bruijn graph is a sophisticated data structure for representing the genomes of entire populations. It robustly encodes shared segments, simple single-nucleotide polymorphisms and complex structural variations far beyond what can be represented in a collection of linear sequences alone. RESULTS: We explore deep topological relationships between suffix trees and compressed de Bruijn graphs and introduce an algorithm, splitMEM, that directly constructs the compressed de Bruijn graph in time and space linear to the total number of genomes for a given maximum genome size. We introduce suffix skips to traverse several suffix links simultaneously and use them to efficiently decompose maximal exact matches into graph nodes. We demonstrate the utility of splitMEM by analyzing the nine-strain pan-genome of Bacillus anthracis and up to 62 strains of Escherichia coli, revealing their core-genome properties.