11. |
- Mahmood, Syed Raza, et al.
(author)
-
β-actin mediated H3K27ac changes demonstrate the link between compartment switching and enhancer-dependent transcriptional regulation
- 2023
-
In: Genome Biology. - : Springer Science and Business Media LLC. - 1465-6906 .- 1474-760X. ; 24
-
Journal article (peer-reviewed)abstract
- Background: Recent work has demonstrated that three-dimensional genome organization is directly affected by changes in the levels of nuclear cytoskeletal proteins such as β-actin. The mechanisms which translate changes in 3D genome structure into changes in transcription, however, are not fully understood. Here, we use a comprehensive genomic analysis of cells lacking nuclear β-actin to investigate the mechanistic links between compartment organization, enhancer activity, and gene expression.Results: Using HiC-Seq, ATAC-Seq, and RNA-Seq, we first demonstrate that transcriptional and chromatin accessibility changes observed upon β-actin loss are highly enriched in compartment-switching regions. Accessibility changes within compartment switching genes, however, are mainly observed in non-promoter regions which potentially represent distal regulatory elements. Our results also show that β-actin loss induces widespread accumulation of the enhancer-specific epigenetic mark H3K27ac. Using the ABC model of enhancer annotation, we then establish that these epigenetic changes have a direct impact on enhancer activity and underlie transcriptional changes observed upon compartment switching. A complementary analysis of fibroblasts undergoing reprogramming into pluripotent stem cells further confirms that this relationship between compartment switching and enhancer-dependent transcriptional change is not specific to β-actin knockout cells but represents a general mechanism linking compartment-level genome organization to gene expression.Conclusions: We demonstrate that enhancer-dependent transcriptional regulation plays a crucial role in driving gene expression changes observed upon compartment-switching. Our results also reveal a novel function of nuclear β-actin in regulating enhancer function by influencing H3K27 acetylation levels.
|
|
12. |
- Malyukova, Alena, et al.
(author)
-
Sequential drug treatment targeting cell cycle and cell fate regulatory programs blocks non-genetic cancer evolution in acute lymphoblastic leukemia
- 2024
-
In: Genome Biology. - 1474-7596 .- 1474-760X. ; 25:1
-
Journal article (peer-reviewed)abstract
- Background: Targeted therapies exploiting vulnerabilities of cancer cells hold promise for improving patient outcome and reducing side-effects of chemotherapy. However, efficacy of precision therapies is limited in part because of tumor cell heterogeneity. A better mechanistic understanding of how drug effect is linked to cancer cell state diversity is crucial for identifying effective combination therapies that can prevent disease recurrence. Results: Here, we characterize the effect of G2/M checkpoint inhibition in acute lymphoblastic leukemia (ALL) and demonstrate that WEE1 targeted therapy impinges on cell fate decision regulatory circuits. We find the highest inhibition of recovery of proliferation in ALL cells with KMT2A-rearrangements. Single-cell RNA-seq and ATAC-seq of RS4;11 cells harboring KMT2A::AFF1, treated with the WEE1 inhibitor AZD1775, reveal diversification of cell states, with a fraction of cells exhibiting strong activation of p53-driven processes linked to apoptosis and senescence, and disruption of a core KMT2A-RUNX1-MYC regulatory network. In this cell state diversification induced by WEE1 inhibition, a subpopulation transitions to a drug tolerant cell state characterized by activation of transcription factors regulating pre-B cell fate, lipid metabolism, and pre-BCR signaling in a reversible manner. Sequential treatment with BCR-signaling inhibitors dasatinib, ibrutinib, or perturbing metabolism by fatostatin or AZD2014 effectively counteracts drug tolerance by inducing cell death and repressing stemness markers. Conclusions: Collectively, our findings provide new insights into the tight connectivity of gene regulatory programs associated with cell cycle and cell fate regulation, and a rationale for sequential administration of WEE1 inhibitors with low toxicity inhibitors of pre-BCR signaling or metabolism.
|
|
13. |
- Meadows, Jennifer, et al.
(author)
-
Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture
- 2023
-
In: Genome Biology. - : BioMed Central (BMC). - 1465-6906 .- 1474-760X. ; 24
-
Journal article (peer-reviewed)abstract
- Background: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 x data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.Results: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.Conclusions: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.
|
|
14. |
- Nordin, Anna, et al.
(author)
-
The CUT & RUN suspect list of problematic regions of the genome
- 2023
-
In: Genome Biology. - : BMC. - 1465-6906 .- 1474-760X. ; 24:1
-
Journal article (peer-reviewed)abstract
- BackgroundCleavage Under Targets and Release Using Nuclease (CUT & RUN) is an increasingly popular technique to map genome-wide binding profiles of histone modifications, transcription factors, and co-factors. The ENCODE project and others have compiled blacklists for ChIP-seq which have been widely adopted: these lists contain regions of high and unstructured signal, regardless of cell type or protein target, indicating that these are false positives. While CUT & RUN obtains similar results to ChIP-seq, its biochemistry and subsequent data analyses are different. We found that this results in a CUT & RUN-specific set of undesired high-signal regions.ResultsWe compile suspect lists based on CUT & RUN data for the human and mouse genomes, identifying regions consistently called as peaks in negative controls. Using published CUT & RUN data from our and other labs, we show that the CUT & RUN suspect regions can persist even when peak calling is performed with SEACR or MACS2 against a negative control and after ENCODE blacklist removal. Moreover, we experimentally validate the CUT & RUN suspect lists by performing reiterative negative control experiments in which no specific protein is targeted, showing that they capture more than 80% of the peaks identified.ConclusionsWe propose that removing these problematic regions can substantially improve peak calling in CUT & RUN experiments, resulting in more reliable datasets.
|
|
15. |
- Ren, Luyao, et al.
(author)
-
Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance
- 2023
-
In: Genome Biology. - : BioMed Central (BMC). - 1465-6906 .- 1474-760X. ; 24:1
-
Journal article (peer-reviewed)abstract
- BACKGROUND: Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome.RESULTS: We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data.CONCLUSIONS: The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.
|
|
16. |
- Sabin, Susanna, et al.
(author)
-
A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex
- 2020
-
In: Genome Biology. - : Springer Science and Business Media LLC. - 1474-7596 .- 1474-760X. ; 21:1
-
Journal article (peer-reviewed)abstract
- BACKGROUND: Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679). RESULTS: A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent. CONCLUSIONS: The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.
|
|
17. |
- Sahlin, Kristoffer, et al.
(author)
-
A survey of mapping algorithms in the long-reads era
- 2023
-
In: Genome Biology. - 1465-6906 .- 1474-760X. ; 24:1
-
Research review (peer-reviewed)abstract
- It has been over a decade since the first publication of a method dedicated entirely to mapping long-reads. The distinctive characteristics of long reads resulted in methods moving from the seed-and-extend framework used for short reads to a seed-and-chain framework due to the seed abundance in each read. The main novelties are based on alternative seed constructs or chaining formulations. Dozens of tools now exist, whose heuristics have evolved considerably. We provide an overview of the methods used in long-read mappers. Since they are driven by implementation-specific parameters, we develop an original visualization tool to understand the parameter settings (http:// bcazaux.polytech-lille.net/Minimap2/).
|
|
18. |
- Sahlin, Kristoffer
(author)
-
Strobealign : flexible seed size enables ultra-fast and accurate read alignment
- 2022
-
In: Genome Biology. - : Springer Science and Business Media LLC. - 1465-6906 .- 1474-760X. ; 23
-
Journal article (peer-reviewed)abstract
- Read alignment is often the computational bottleneck in analyses. Recently, several advances have been made on seeding methods for fast sequence comparison. We combine two such methods, syncmers and strobemers, in a novel seeding approach for constructing dynamic-sized fuzzy seeds and implement the method in a short-read aligner, strobealign. The seeding is fast to construct and effectively reduces repetitiveness in the seeding step, as shown using a novel metric E-hits. strobealign is several times faster than traditional aligners at similar and sometimes higher accuracy while being both faster and more accurate than more recently proposed aligners for short reads of lengths 150nt and longer. Availability: https://github.com/ksahlin/strobealign
|
|
19. |
- Sakthikumar, Sharadha, et al.
(author)
-
Whole-genome sequencing of glioblastoma reveals enrichment of non-coding constraint mutations in known and novel genes
- 2020
-
In: Genome Biology. - : Springer Science and Business Media LLC. - 1465-6906 .- 1474-760X. ; 21:1
-
Journal article (peer-reviewed)abstract
- BackgroundGlioblastoma (GBM) has one of the worst 5-year survival rates of all cancers. While genomic studies of the disease have been performed, alterations in the non-coding regulatory regions of GBM have largely remained unexplored. We apply whole-genome sequencing (WGS) to identify non-coding mutations, with regulatory potential in GBM, under the hypothesis that regions of evolutionary constraint are likely to be functional, and somatic mutations are likely more damaging than in unconstrained regions.ResultsWe validate our GBM cohort, finding similar copy number aberrations and mutated genes based on coding mutations as previous studies. Performing analysis on non-coding constraint mutations and their position relative to nearby genes, we find a significant enrichment of non-coding constraint mutations in the neighborhood of 78 genes that have previously been implicated in GBM. Among them, SEMA3C and DYNC1I1 show the highest frequencies of alterations, with multiple mutations overlapping transcription factor binding sites. We find that a non-coding constraint mutation in the SEMA3C promoter reduces the DNA binding capacity of the region. We also identify 1776 other genes enriched for non-coding constraint mutations with likely regulatory potential, providing additional candidate GBM genes. The mutations in the top four genes, DLX5, DLX6, FOXA1, and ISL1, are distributed over promoters, UTRs, and multiple transcription factor binding sites.ConclusionsThese results suggest that non-coding constraint mutations could play an essential role in GBM, underscoring the need to connect non-coding genomic variation to biological function and disease pathology.
|
|
20. |
- Sounart, Hailey, et al.
(author)
-
Dual spatially resolved transcriptomics for human host–pathogen colocalization studies in FFPE tissue sections
- 2023
-
In: Genome Biology. - : Springer Nature. - 1465-6906 .- 1474-760X. ; 24:1
-
Journal article (peer-reviewed)abstract
- Technologies to study localized host–pathogen interactions are urgently needed. Here, we present a spatial transcriptomics approach to simultaneously capture host and pathogen transcriptome-wide spatial gene expression information from human formalin-fixed paraffin-embedded (FFPE) tissue sections at a near single-cell resolution. We demonstrate this methodology in lung samples from COVID-19 patients and validate our spatial detection of SARS-CoV-2 against RNAScope and in situ sequencing. Host–pathogen colocalization analysis identified putative modulators of SARS-CoV-2 infection in human lung cells. Our approach provides new insights into host response to pathogen infection through the simultaneous, unbiased detection of two transcriptomes in FFPE samples.
|
|