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.

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.

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.

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.

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.

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.

Genomic dark matter: the reliability of short read mapping illustrated by the genome mappability score.

MOTIVATION: Genome resequencing and short read mapping are two of the primary tools of genomics and are used for many important applications. The current state-of-the-art in mapping uses the quality values and mapping quality scores to evaluate the reliability of the mapping. These attributes, however, are assigned to individual reads and do not directly measure the problematic repeats across the genome. Here, we present the Genome Mappability Score (GMS) as a novel measure of the complexity of resequencing a genome. The GMS is a weighted probability that any read could be unambiguously mapped to a given position and thus measures the overall composition of the genome itself. RESULTS: We have developed the Genome Mappability Analyzer to compute the GMS of every position in a genome. It leverages the parallelism of cloud computing to analyze large genomes, and enabled us to identify the 5-14% of the human, mouse, fly and yeast genomes that are difficult to analyze with short reads. We examined the accuracy of the widely used BWA/SAMtools polymorphism discovery pipeline in the context of the GMS, and found discovery errors are dominated by false negatives, especially in regions with poor GMS. These errors are fundamental to the mapping process and cannot be overcome by increasing coverage. As such, the GMS should be considered in every resequencing project to pinpoint the 'dark matter' of the genome, including of known clinically relevant variations in these regions. AVAILABILITY: The source code and profiles of several model organisms are available at http://gma-bio.sourceforge.net

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.

Sulfur-containing mixed-donor tribenzo-macrocycles and their endo- and exocyclic supramolecular silver(I) and copper(I) complexes.

New sulfur-containing mixed-donor macrocycles L(1)-L(4) with different donor sets and/or ring sizes (L(1), 20-membered O(3)S(2); L(2), 20-membered NO(2)S(2); L(3), 20-membered O(2)S(3); and L(4), 23-membered O(4)S(2)) were synthesized and structurally characterized by X-ray analysis. A comparative investigation of the coordination behavior of these macrocyclic ligands with silver(I) salts and copper(I) iodide is reported. The X-ray structures of seven complexes (1-7) have been determined, and a range of structural types and coordination modes, including mono- to multinuclear and endo- to exocyclic ones via M-S bonds, is shown to occur. In all seven complexes, the ether oxygens of each ring are unbound. Reactions of L(1)-L(4) with the silver(I) salts (ClO(4)(-), BF(4)(-), and NO(3)(-)) afforded four complexes with different topologies: the exocyclic 1-D coordination polymer [Ag(L(1))(ClO(4))](n) (1), the endocyclic 1:1 monomer complex [Ag(L(2))]ClO(4).CH(3)OH (2), the endo-/exocyclic 1-D coordination polymer {[Ag(2)(L(2))(CH(3)OH)](BF(4))(2)}(n) (3), and the cyclic dimer complex [Ag(2)(L(3))(2)(NO(3))(2)].2CH(3)CN (4). NMR titrations of L(1) and L(2) with silver(I) perchlorate were carried out to explore their complexation behavior in solution and also for comparison with the solid state structures. As well, reaction of the above ligands with copper(I) iodide afforded the exocyclic dimer, [Cu(2)I(2)(L(1))(2)].CH(2)Cl(2) (5), linked with diamond-type Cu-I(2)-Cu unit, the exocyclic monomer complex [CuI(L(3))].CH(2)Cl(2) (6), and the emissive 1-D poly(cyclic dimer) complex {[Cu(4)I(4)(L(4))(2)].CH(2)Cl(2)}(n) (7) linked by a cubane-type Cu(4)I(4) unit. The formation of both the discrete and continuous supramolecular complexes with copper(I) iodide is discussed in terms of a ligand-directed "sulfur-to-sulfur" distance effect.

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.

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.

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.

Reevaluating the capability of Taq DNA polymerase: long PCR amplification.

We tested the ability of Taq DNA polymerase (Taq) to amplify long DNA fragments and showed that, if the conditions were set properly, Taq could successfully perform the "long PCR" up to 24 kb. The conditions include: (1) longer primers, (2) a 2-step cycling, and (3) a "long buffer." We propose that the most important requirements are the survival rate of Taq at high temperatures and that of the primers against the 5' to 3' exonuclease activity of Taq.

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.

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.

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.

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.