Searching for a Needle in a Haystack: Cas9-Targeted Nanopore Sequencing and DNA Methylation Profiling of Full-Length Glutenin Genes in a Big Cereal Genome

Sequencing and epigenetic profiling of target genes in plants are important tasks with various applications ranging from marker design for plant breeding to the study of gene expression regulation. This is particularly interesting for plants with big genome size for which whole-genome sequencing can be time-consuming and costly. In this study, we asked whether recently proposed Cas9-targeted nanopore sequencing (nCATS) is efficient for target gene sequencing for plant species with big genome size. We applied nCATS to sequence the full-length glutenin genes (Glu-1Ax, Glu-1Bx and Glu-1By) and their promoters in hexaploid triticale (X Triticosecale, AABBRR, genome size is 24 Gb). We showed that while the target gene enrichment per se was quite high for the three glutenin genes (up to 645×), the sequencing depth that was achieved from two MinION flowcells was relatively low (5–17×). However, this sequencing depth was sufficient for various tasks including detection of InDels and single-nucleotide variations (SNPs), read phasing and methylation profiling. Using nCATS, we uncovered SNP and InDel variation of full-length glutenin genes providing useful information for marker design and deciphering of variation of individual Glu-1By alleles. Moreover, we demonstrated that glutenin genes possess a ‘gene-body’ methylation epigenetic profile with hypermethylated CDS part and hypomethylated promoter region. The obtained information raised an interesting question on the role of gene-body methylation in glutenin gene expression regulation. Taken together, our work disclosures the potential of the nCATS approach for sequencing of target genes in plants with big genome size.

A merger between compatible but divergent genomes supports allopolyploidization in the Brassicaceae family

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species rarely occur in nature mainly due to reproductive barriers but how such hurdles can be overcome is largely unknown. xBrassicoraphanus is a fertile intergeneric allopolyploid synthesized between Brassica rapa and Raphanus sativus in the Brassicaceae family. Genomes of B. rapa and R. sativus are diverged enough to suppress synapsis formation between non-homologous progenitor chromosomes during meiosis, and we found that both genomes reside in the single nucleus of xBrassicoraphanus without genome loss or rearrangement. Expressions of syntenic orthologs identified in B. rapa and R. sativus were adjusted to a hybrid nuclear environment of xBrassicoraphanus, which necessitates reconfiguration of transcription network by rewiring cis-trans interactions. B. rapa coding sequences have a higher level of gene-body methylation than R. sativus, and such methylation asymmetry is maintained in xBrassicoraphanus. B. rapa-originated transposable elements were transcriptionally silenced in xBrassicoraphanus, rendered by gain of CHG methylation in trans via small RNAs derived from the same sequences of R. sativus subgenome. Our work proposes that not only transcription compatibility but also a certain extent of genome divergence supports hybrid genome stabilization, which may explain great diversification and expansion of angiosperms during evolution.

Gene Body Methylation Confers Transcription Robustness in Mangroves During Long-Term Stress Adaptation

Whether induced epigenetic changes contribute to long-term adaptation remains controversial. Recent studies indicate that environmentally cued changes in gene body methylation (gbM) can facilitate acclimatization. However, such changes are often associated with genetic variation and their contribution to long-term stress adaptation remains unclear. Using whole-genome bisulfite sequencing, we examined evolutionary gains and losses of gbM in mangroves that adapted to extreme intertidal environments. We treated mangrove seedlings with salt stress, and investigated expression changes in relation with stress-induced or evolutionarily-acquired gbM changes. Evolution and function of gbM was compared with that of genetic variation. Mangroves gained much more gbM than their terrestrial relatives, mainly through convergent evolution. Genes that convergently gained gbM during evolution are more likely to become methylated in response to salt stress in species where they are normally not marked. Stress-induced and evolutionarily convergent gains of gbM both correlate with reduction in expression variation, conferring genome-wide expression robustness under salt stress. Moreover, convergent gbM evolution is uncoupled with convergent sequence evolution. Our findings suggest that transgenerational inheritance of acquired gbM helps environmental canalization of gene expression, facilitating long-term stress adaptation of mangroves in the face of a severe reduction in genetic diversity.

ASXL1 mutations are associated with distinct epigenomic alterations that lead to sensitivity to venetoclax and azacytidine

The BCL2-inhibitor, Venetoclax (VEN), has shown significant anti-leukemic efficacy in combination with the DNMT-inhibitor, Azacytidine (AZA). To explore the mechanisms underlying the selective sensitivity of mutant leukemia cells to VEN and AZA, we used cell-based isogenic models containing a common leukemia-associated mutation in the epigenetic regulator ASXL1. KBM5 cells with CRISPR/Cas9-mediated correction of the ASXL1G710X mutation showed reduced leukemic growth, increased myeloid differentiation, and decreased HOXA and BCL2 gene expression in vitro compared to uncorrected KBM5 cells. Increased expression of the anti-apoptotic gene, BCL2, was also observed in bone marrow CD34+ cells from ASXL1 mutant MDS patients compared to CD34+ cells from wild-type MDS cases. ATAC-sequencing demonstrated open chromatin at the BCL2 promoter in the ASXL1 mutant KBM5 cells. BH3 profiling demonstrated increased dependence of mutant cells on BCL2. Upon treatment with VEN, mutant cells demonstrated increased growth inhibition. In addition, genome-wide methylome analysis of primary MDS samples and isogenic cell lines demonstrated increased gene-body methylation in ASXL1 mutant cells, with consequently increased sensitivity to AZA. These data mechanistically link the common leukemia-associated mutation ASXL1 to enhanced sensitivity to VEN and AZA via epigenetic upregulation of BCL2 expression and widespread alterations in DNA methylation.

Repression of CHROMOMETHYLASE 3 prevents epigenetic collateral damage in Arabidopsis

DNA methylation has evolved to silence mutagenic transposable elements (TEs) while typically avoiding the targeting of endogenous genes. Mechanisms that prevent DNA methyltransferases from ectopically methylating genes are expected to be of prime importance during periods of dynamic cell cycle activities including plant embryogenesis. However, virtually nothing is known regarding how DNA methyltransferase activities are precisely regulated during embryogenesis to prevent the induction of potentially deleterious and mitotically stable genic epimutations. Here, we report that microRNA-mediated repression of CHROMOMETHYLASE 3 (CMT3) and the chromatin features that CMT3 prefers help prevent ectopic methylation of thousands of genes during embryogenesis that can persist for weeks afterwards. Our results are also consistent with CMT3-induced ectopic methylation of promoters or bodies of genes undergoing transcriptional activation reducing their expression. Therefore, the repression of CMT3 prevents epigenetic collateral damage on endogenous genes. We also provide a model that may help reconcile conflicting viewpoints regarding the functions of gene-body methylation that occurs in nearly all flowering plants.

Integrated transcriptomics and epigenomics reveal chamber-specific and species-specific characteristics of human and mouse hearts

DNA methylation, chromatin accessibility, and gene expression represent different levels information in biological process, but a comprehensive multiomics analysis of the mammalian heart is lacking. Here, we applied nucleosome occupancy and methylome sequencing, which detected DNA methylation and chromatin accessibility simultaneously, as well as RNA-seq, for multiomics analysis of the 4 chambers of adult and fetal human hearts, and adult mouse hearts. Our results showed conserved region-specific patterns in the mammalian heart at transcriptome and DNA methylation level. Adult and fetal human hearts showed distinct features in DNA methylome, chromatin accessibility, and transcriptome. Novel long noncoding RNAs were identified in the human heart, and the gene expression profiles of major cardiovascular diseases associated genes were displayed. Furthermore, cross-species comparisons revealed human-specific and mouse-specific differentially expressed genes between the atria and ventricles. We also reported the relationship among multiomics and found there was a bell-shaped relationship between gene-body methylation and expression in the human heart. In general, our study provided comprehensive spatiotemporal and evolutionary insights into the regulation of gene expression in the heart.

Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant–pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.