Genomic Targets Research Articles

Evaluation of Interfering RNA Efficacy in Treating Hepatitis B: Is It Promising?

Background: Despite an existing safe and effective vaccine for hepatitis B virus (HBV), it is still a major public health concern. Nowadays, several drugs are used to treat chronic hepatitis B; however, full healing remains controversial. The viral covalently closed circular DNA (cccDNA) formed by HBV forms a major challenge in its treatment, as does the ability of HBV to integrate itself into the host genome, which enables infection reactivation. Interfering RNA (RNAi) is a gene-silencing post-transcriptional mechanism which forms as a promising alternative to treat chronic hepatitis B. The aim of the present review is to assess the evolution of hepatitis B treatment approaches based on using RNA interference. Methods: Data published between 2016 and 2023 in scientific databases (PubMed, PMC, LILACS, and Bireme) were assessed. Results: In total, 76,949 articles were initially identified and quality-checked, and 226 eligible reports were analyzed in depth. The main genomic targets, delivery systems, and major HBV therapy innovations are discussed in this review. This review reinforces the therapeutic potential of RNAi and identifies the need for conducting further studies to fill the remaining gaps between bench and clinical practice.

A Multiple Targeting Genome Editing System for Remodulation of Circulating Malignant Cells to Eliminate Cancer Immunosuppression and Restore Immune Responses.

Cancer immunotherapy, which aims to eliminate cancer immunosuppression and reactivate anticancer immunity, holds great promise in oncology treatments. However, it is challenging to accurately study the efficacy of immunotherapy based on human-derived cells through animal experiments due to xenogeneic immune rejection. Herein, a personalized and precise strategy to evaluate the effectiveness of immunotherapy using the blood samples of cancer patients is presented. Through the utilization of multiple cancer-targeting delivery system decorated with the epidermal growth factor receptor (EGFR)-specific aptamer CL4 and the AXL-specific aptamer GL21.T to achieve superior efficiency in delivering the genome editing plasmid for MUC1 knockout, effective modulation on the behavior of circulating malignant cells (CMCs) is realized. After genome editing, both mucin 1 (MUC1) and programmed death-ligand 1 (PD-L1) are significantly downregulated in CMCs. The elimination of immunosuppression results in markedly enhanced secretion of pro-inflammatory anticancer cytokines encompassing interleukins 2, 12, and 15 and interferon-γ by immune cells. The study not only provides a strategy to overcome immunosuppression but also yields critical insights for personalized immunotherapy approaches.

Insights into trait-association of selection signatures and adaptive eQTL in indigenous African cattle

BackgroundAfrican cattle represent a unique resource of genetic diversity in response to adaptation to numerous environmental challenges. Characterising the genetic landscape of indigenous African cattle and identifying genomic regions and genes of functional importance can contribute to targeted breeding and tackle the loss of genetic diversity. However, pinpointing the adaptive variant and determining underlying functional mechanisms of adaptation remains challenging.ResultsIn this study, we use selection signatures from whole-genome sequence data of eight indigenous African cattle breeds in combination with gene expression and quantitative trait loci (QTL) databases to characterise genomic targets of artificial selection and environmental adaptation and to identify the underlying functional candidate genes. In general, the trait-association analyses of selection signatures suggest the innate and adaptive immune system and production traits as important selection targets. For example, a large genomic region, with selection signatures identified for all breeds except N’Dama, was located on BTA27, including multiple defensin DEFB coding-genes. Out of 22 analysed tissues, genes under putative selection were significantly enriched for those overexpressed in adipose tissue, blood, lung, testis and uterus. Our results further suggest that cis-eQTL are themselves selection targets; for most tissues, we found a positive correlation between allele frequency differences and cis-eQTL effect size, suggesting that positive selection acts directly on regulatory variants.ConclusionsBy combining selection signatures with information on gene expression and QTL, we were able to reveal compelling candidate selection targets that did not stand out from selection signature results alone (e.g. GIMAP8 for tick resistance and NDUFS3 for heat adaptation). Insights from this study will help to inform breeding and maintain diversity of locally adapted, and hence important, breeds.

MHConstructor: a high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short-read, de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target capture short-read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short-read data. MHConstructor facilitates wide-spread access to high-quality, alignment-free MHC sequence analysis.

Real-time genomic characterization of pediatric acute leukemia using adaptive sampling.

Effective treatment of pediatric acute leukemia is dependent on accurate genomic classification, typically derived from a combination of multiple time-consuming and costly techniques such as flow cytometry, fluorescence in situ hybridization (FISH), karyotype analysis, targeted PCR, and microarrays. We investigated the feasibility of a comprehensive single-assay classification approach using long-read sequencing, with real-time genome target enrichment, to classify chromosomal abnormalities and structural variants characteristic of acute leukemia. We performed whole genome sequencing on DNA from diagnostic peripheral blood or bone marrow for 54 pediatric acute leukemia cases with diverse genomic subtypes. We demonstrated the characterization of known, clinically relevant karyotype abnormalities and structural variants concordant with standard-of-care clinical testing. Subtype-defining genomic alterations were identified in all cases following a maximum of forty-eight hours of sequencing. In 18 cases, we performed real-time analysis - concurrent with sequencing - and identified the driving alteration in as little as fifteen minutes (for karyotype) or up to six hours (for complex structural variants). Whole genome nanopore sequencing with adaptive sampling has the potential to provide detailed genomic classification of acute leukemia specimens with reduced cost and turnaround time compared to the current standard of care.

The multifaceted roles of phosphoethanolamine-modified lipopolysaccharides: from stress response and virulence to cationic antimicrobial resistance.

SUMMARYLipopolysaccharides (LPS) are an integral part of the outer membrane of Gram-negative bacteria and play essential structural and functional roles in maintaining membrane integrity as well as in stress response and virulence. LPS comprises a membrane-anchored lipid A group, a sugar-based core region, and an O-antigen formed by repeating oligosaccharide units. 3-Deoxy-D-manno-octulosonic acid-lipid A (Kdo2-lipid A) is the minimum LPS component required for bacterial survival. While LPS modifications are not essential, they play multifaceted roles in stress response and host-pathogen interactions. Gram-negative bacteria encode several distinct LPS-modifying phosphoethanolamine transferases (PET) that add phosphoethanolamine (pEtN) to lipid A or the core region of LPS. The pet genes differ in their genomic locations, regulation mechanisms, and modification targets of the encoded enzyme, consistent with their various roles in different growth niches and under varied stress conditions. The discovery of mobile colistin resistance genes, which represent lipid A-modifying pet genes that are encoded on mobile elements and associated with resistance to the last-resort antibiotic colistin, has led to substantial interest in PETs and pEtN-modified LPS over the last decade. Here, we will review the current knowledge of the functional diversity of pEtN-based LPS modifications, including possible roles in niche-specific fitness advantages and resistance to host-produced antimicrobial peptides, and discuss how the genetic and structural diversities of PETs may impact their function. An improved understanding of the PET group will further enhance our comprehension of the stress response and virulence of Gram-negative bacteria and help contextualize host-pathogen interactions.

Museomics reveal origins of East African Pleophylla forest chafers and Miocene forest connectivity

Here we present a nearly complete species-level phylogeny including 23 of the 25 known species of the forest-dwelling herbivorous scarab chafer beetle genus Pleophylla (Coleoptera: Scarabaeidae: Sericinae), based on the analysis of 950 nuclear genes (metazoan-level universal single-copy orthologs; mzl-USCOs). DNA sequences were obtained from freshly collected, ethanol-preserved samples and from dried museum specimens by target enrichment or genome shotgun sequencing. Alignment completeness of mzl-USCOs newly obtained here by target DNA enrichment of ethanol samples were very heterogenous and lower (29–62 %) than in Dietz et al. (2023a), while that of sequences recovered from dried samples was even lower (∼19 %). Alignment completeness of the sequences obtained from low coverage shotgun sequencing was highest (∼92 %), although the average coverage was much lower than for the target enrichment samples. We used the resulting phylogeny to reconstruct the historical biogeography of the group. To estimate a time-calibrated tree, we combined the mzl-USCO data of Pleophylla with a nucleotide alignment from an available transcriptomic dataset of Scarabaeoidea and used two different sets of secondary calibration points. Despite the problems associated with the capture rate of mzl-USCO sequences from museum specimens, we were able to infer a well-resolved phylogeny of the genus Pleophylla that also provided reliable estimates of the phylogenetic position of species for which we had little sequence data. Our study clearly identified South Africa as the geographic origin of Pleophylla. Timing and biogeographic history confirm a persistent fragmentation of forests since the Eocene. The occurrence of only one long-distance dispersal event from southern Africa to the Eastern African Arc even during the Miocene highlights the limited dispersal possibilities for these forest-adapted chafers, which do not seem to have had important northerly range expansions along hypothetical forest corridors during the Pleistocene.

ENL mutation and AML: a new model that reveals oncogenic condensate's function in leukemogenesis.

Precise regulation of gene expression is essential for proper development and the maintenance of homeostasis in organisms. Studies have shown that some transcriptional regulatory proteins influence gene expression through the formation of dynamic, locally concentrated assemblies known as condensates, while dysregulation of transcriptional condensates was associated with several cancers, such as Ewing sarcoma and AML [Wang Y etal. (2023) Nat Chem Biol 19, 1223-1234; Chandra B etal. (2022) Cancer Discov 12, 1152-1169]. Mutations in the histone acetylation "reader" eleven-nineteen-leukemia (ENL) have been shown to form discrete condensates at endogenous genomic targets, but it remains unclear how ENL mutations drive tumorigenesis and whether it is correlated with their condensate formation property. Liu etal. now show, using a conditional knock-in mouse model, that ENL YEATS domain mutation is a bona fide oncogenic driver for AML. This mutant ENL forms condensates in hematopoietic stem/progenitor cells at the genomic loci of key leukemogenic genes, including Meis1 and Hoxa cluster genes, and disrupting condensate formation via mutagenesis impairs its chromatin and oncogenic function. Furthermore, they show that small-molecule inhibition of the acetyl-binding activity displaces ENL mutant condensates from oncogenic target loci, and this inhibitor significantly impairs the onset and progression of AML driven by mutant ENL in vivo.

P53motifDB: integration of genomic information and tumor suppressor p53 binding motifs.

The tumor suppressor gene TP53 encodes the DNA binding transcription factor p53 and is one of the most commonly mutated genes in human cancer. Tumor suppressor activity requires binding of p53 to its DNA response elements and subsequent transcriptional activation of a diverse set of target genes. Despite decades of close study, the logic underlying p53 interactions with its numerous potential genomic binding sites and target genes is not yet fully understood. Here, we present a database of DNA and chromatin-based information focused on putative p53 binding sites in the human genome to allow users to generate and test new hypotheses related to p53 activity in the genome. Users can query genomic locations based on experimentally observed p53 binding, regulatory element activity, genetic variation, evolutionary conservation, chromatin modification state, and chromatin structure. We present multiple use cases demonstrating the utility of this database for generating novel biological hypotheses, such as chromatin-based determinants of p53 binding and potential cell type-specific p53 activity. All database information is also available as a precompiled sqlite database for use in local analysis or as a Shiny web application.

ASPYRE-Lung: validation of a simple, fast, robust and novel method for multi-variant genomic analysis of actionable NSCLC variants in FFPE tissue.

Genomic variant testing of tumors is a critical gateway for patients to access the full potential of personalized oncology therapeutics. Current methods such as next-generation sequencing are costly and challenging to interpret, while PCR assays are limited in the number of variants they can cover. We developed ASPYRE® (Allele-Specific PYrophosphorolysis REaction) technology to address the urgent need for rapid, accessible and affordable diagnostics informing actionable genomic target variants of a given cancer. The targeted ASPYRE-Lung panel for non-small cell carcinoma covers 114 variants in 11 genes (ALK, BRAF, EGFR, ERBB2, KRAS, RET, ROS1, MET & NTRK1/2/3) to robustly inform clinical management. The assay detects single nucleotide variants, insertions, deletions, and gene fusions from tissue-derived DNA and RNA simultaneously. We tested the limit of detection, specificity, analytical accuracy and analytical precision of ASPYRE-Lung using FFPE lung tissue samples from patients with non-small cell lung carcinoma, variant-negative FFPE tissue from healthy donors, and FFPE-based contrived samples with controllable variant allele fractions. The sensitivity of ASPYRE-Lung was determined to be ≤ 3% variant allele fraction for single nucleotide variants and insertions or deletions, 100 copies for fusions, and 200 copies for MET exon 14 skipping. The specificity was 100% with no false positive results. The analytical accuracy test yielded no discordant calls between ASPYRE-Lung and expected results for clinical samples (via orthogonal testing) or contrived samples, and results were replicable across operators, reagent lots, runs, and real-time PCR instruments with a high degree of precision. The technology is simple and fast, requiring only four reagent transfer steps using standard laboratory equipment (PCR and qPCR instruments) with analysis via a cloud-based algorithm. The ASPYRE-Lung assay has the potential to be transformative in facilitating access to rapid, actionable molecular profiling of tissue for patients with non-small cell carcinoma.

Regulatory insight for a Zn2Cys6 transcription factor controlling effector-mediated virulence in a fungal pathogen of wheat.

The regulation of virulence in plant-pathogenic fungi has emerged as a key area of importance underlying host infections. Recent work has highlighted individual transcription factors (TFs) that serve important roles. A prominent example is PnPf2, a member of the Zn2Cys6 family of fungal TFs, which controls the expression of effectors and other virulence-associated genes in Parastagonospora nodorum during infection of wheat. PnPf2 orthologues are similarly important for other major fungal pathogens during infection of their respective host plants, and have also been shown to control polysaccharide metabolism in model saprophytes. In each case, the direct genomic targets and associated regulatory mechanisms were unknown. Significant insight was made here by investigating PnPf2 through chromatin-immunoprecipitation (ChIP) and mutagenesis approaches in P. nodorum. Two distinct binding motifs were characterised as positive regulatory elements and direct PnPf2 targets identified. These encompass known effectors and other components associated with the P. nodorum pathogenic lifestyle, such as carbohydrate-active enzymes and nutrient assimilators. The results support a direct involvement of PnPf2 in coordinating virulence on wheat. Other prominent PnPf2 targets included TF-encoding genes. While novel functions were observed for the TFs PnPro1, PnAda1, PnEbr1 and the carbon-catabolite repressor PnCreA, our investigation upheld PnPf2 as the predominant transcriptional regulator characterised in terms of direct and specific coordination of virulence on wheat, and provides important mechanistic insights that may be conserved for homologous TFs in other fungi.

Zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters.

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of ortholog groups for individual protein-encoding genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of evolutionary statistics for each inferred ortholog group. Importantly, in comparison to tools for visual exploration of homologous relationships between gene clusters, zol can scale to thousands of gene cluster instances and produce detailed reports that are easy to digest. To showcase fai and zol, we apply them for: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-level genetic insights of two common BGCs in the fungal species Aspergillus flavus, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Making target sites in large structured RNAs accessible to RNA-cleaving DNAzymes through hybridization with synthetic DNA oligonucleotides.

The 10-23 DNAzyme is one of the most active DNA-based enzymes, and in theory, can be designed to target any purine-pyrimidine junction within an RNA sequence for cleavage. However, purine-pyrimidine junctions within a large, structured RNA (lsRNA) molecule of biological origin are not always accessible to 10-23, negating its general utility as an RNA-cutting molecular scissor. Herein, we report a generalizable strategy that allows 10-23 to access any purine-pyrimidine junction within an lsRNA. Using three large SARS-CoV-2 mRNA sequences of 566, 584 and 831 nucleotides in length as model systems, we show that the use of antisense DNA oligonucleotides (ASOs) that target the upstream and downstream regions flanking the cleavage site can restore the activity (kobs) of previously poorly active 10-23 DNAzyme systems by up to 2000-fold. We corroborated these findings mechanistically using in-line probing to demonstrate that ASOs reduced 10-23 DNAzyme target site structure within the lsRNA substrates. This approach represents a simple, efficient, cost-effective, and generalizable way to improve the accessibility of 10-23 to a chosen target site within an lsRNA molecule, especially where direct access to the genomic RNA target is necessary.

Characterisation of the erythrocyte invasion phenotype of FCB-2: A South American P. falciparum reference strain

The extent of parasite adaptive capability involved in erythrocyte invasion represents a significant challenge for the development of a Plasmodium falciparum vaccine. The parasite's geographical and populational origin may influence such adaptive behaviour; in vitro culture-adapted parasite strains are typically used for such studies. Previous studies have reported invasion phenotypes in strains from Africa and Asia and, to a lesser extent, from Latin America. This study was aimed at expanding the pool of characterised parasite strains from Latin America by describing the invasion phenotype of the P. falciparum Colombia Bogotá 2 (FCB2) strain. The FCB2 genome was sequenced and erythrocyte invasion ligand sequences were analysed and compared to other previously reported ones. RT-PCR was used for assessing Pfeba family erythrocyte invasion ligands and reticulocyte binding homologue (Pfrh) gene transcription. A flow cytometry-based erythrocyte invasion assay (using enzymatically-treated erythrocytes) was used for determining the FCB2 strain's invasion phenotype. The P. falciparum FCB2 genome sequence was analysed, bearing in mind that prolonged in vitro parasite culture may affect its genome sequence and, in some cases, lead to the deletion of certain genes; it was demonstrated that all erythrocyte invasion ligand gene sequences studied here were preserved. Comparative analysis showed that the target genome sequences were conserved whereas transcriptional analysis highlighted Pfebas and Pfrhs gene expression. Erythrocyte invasion analysis demonstrated that the FCB2 strain has a sialic acid-resistant invasion phenotype.

Robust genome and cell engineering via in vitro and in situ circularized RNAs.

Circularization can improve RNA persistence, yet simple and scalable approaches to achieve this are lacking. Here we report two methods that facilitate the pursuit of circular RNAs (cRNAs): cRNAs developed via in vitro circularization using group II introns, and cRNAs developed via in-cell circularization by the ubiquitously expressed RtcB protein. We also report simple purification protocols that enable high cRNA yields (40-75%) while maintaining low immune responses. These methods and protocols facilitate a broad range of applications in stem cell engineering as well as robust genome and epigenome targeting via zinc finger proteins and CRISPR-Cas9. Notably, cRNAs bearing the encephalomyocarditis internal ribosome entry enabled robust expression and persistence compared with linear capped RNAs in cardiomyocytes and neurons, which highlights the utility of cRNAs in these non-dividing cells. We also describe genome targeting via deimmunized Cas9 delivered as cRNA and a long-range multiplexed protein engineering methodology for the combinatorial screening of deimmunized protein variants that enables compatibility between persistence of expression and immunogenicity in cRNA-delivered proteins. The cRNA toolset will aid research and the development of therapeutics.

ARID1 is required to regulate and reinforce H3K9me2 in sperm cells in Arabidopsis

Heterochromatin de-condensation in companion gametic cells is conserved in both plants and animals. In plants, microspore undergoes asymmetric pollen mitosis (PMI) to produce a vegetative cell (VC) and a generative cell (GC). Subsequently, the GC undergoes pollen mitosis (PMII) to produce two sperm cells (SC). Consistent with heterochromatin de-condensation in the VC, H3K9me2, a heterochromatin mark, is barely detected in VC. However, how H3K9me2 is differentially regulated during pollen mitosis remains unclear. Here, we show that H3K9me2 is gradually evicted from the VC since PMI but remain unchanged in the GC and SC. ARID1, a pollen-specific transcription factor that facilitates PMII, promotes H3K9me2 maintenance in the GC/SC but slows down its eviction in the VC. The genomic targets of ARID1 mostly overlaps with H3K9me2 loci, and ARID1 recruits H3K9 methyltransferase SUVH6. Our results uncover that differential pattern of H3K9me2 between two cell types is regulated by ARID1 during pollen mitosis.

Maptcha: an efficient parallel workflow for hybrid genome scaffolding

BackgroundGenome assembly, which involves reconstructing a target genome, relies on scaffolding methods to organize and link partially assembled fragments. The rapid evolution of long read sequencing technologies toward more accurate long reads, coupled with the continued use of short read technologies, has created a unique need for hybrid assembly workflows. The construction of accurate genomic scaffolds in hybrid workflows is complicated due to scale, sequencing technology diversity (e.g., short vs. long reads, contigs or partial assemblies), and repetitive regions within a target genome.ResultsIn this paper, we present a new parallel workflow for hybrid genome scaffolding that would allow combining pre-constructed partial assemblies with newly sequenced long reads toward an improved assembly. More specifically, the workflow, called Maptcha, is aimed at generating long scaffolds of a target genome, from two sets of input sequences—an already constructed partial assembly of contigs, and a set of newly sequenced long reads. Our scaffolding approach internally uses an alignment-free mapping step to build a ⟨\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\langle $$\\end{document}contig,contig⟩\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rangle $$\\end{document} graph using long reads as linking information. Subsequently, this graph is used to generate scaffolds. We present and evaluate a graph-theoretic “wiring” heuristic to perform this scaffolding step. To enable efficient workload management in a parallel setting, we use a batching technique that partitions the scaffolding tasks so that the more expensive alignment-based assembly step at the end can be efficiently parallelized. This step also allows the use of any standalone assembler for generating the final scaffolds.ConclusionsOur experiments with Maptcha on a variety of input genomes, and comparison against two state-of-the-art hybrid scaffolders demonstrate that Maptcha is able to generate longer and more accurate scaffolds substantially faster. In almost all cases, the scaffolds produced by Maptcha are at least an order of magnitude longer (in some cases two orders) than the scaffolds produced by state-of-the-art tools. Maptcha runs significantly faster too, reducing time-to-solution from hours to minutes for most input cases. We also performed a coverage experiment by varying the sequencing coverage depth for long reads, which demonstrated the potential of Maptcha to generate significantly longer scaffolds in low coverage settings (1×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1\ imes $$\\end{document}–10×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$10\ imes $$\\end{document}).

SAMD1 suppresses epithelial-mesenchymal transition pathways in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) poses a significant threat due to its tendency to evade early detection, frequent metastasis, and the subsequent challenges in devising effective treatments. Processes that govern epithelial-mesenchymal transition (EMT) in PDAC hold promise for advancing novel therapeutic strategies. SAMD1 (SAM domain-containing protein 1) is a CpG island-binding protein that plays a pivotal role in the repression of its target genes. Here, we revealed that SAMD1 acts as a repressor of genes associated with EMT. Upon deletion of SAMD1 in PDAC cells, we observed significantly increased migration rates. SAMD1 exerts its effects by binding to specific genomic targets, including CDH2, encoding N-cadherin, which emerged as a driver of enhanced migration upon SAMD1 knockout. Furthermore, we discovered the FBXO11-containing E3 ubiquitin ligase complex as an interactor and negative regulator of SAMD1, which inhibits SAMD1 chromatin-binding genome-wide. High FBXO11 expression in PDAC is associated with poor prognosis and increased expression of EMT-related genes, underlining an antagonistic relationship between SAMD1 and FBXO11. In summary, our findings provide insights into the regulation of EMT-related genes in PDAC, shedding light on the intricate role of SAMD1 and its interplay with FBXO11 in this cancer type.

Capture of circulating metastatic cancer cell clusters from lung cancer patients can reveal unique genomic profiles and potential anti-metastatic molecular targets: A proof-of-concept study.

Metastasis remains the leading cause of cancer deaths worldwide and lung cancer, known for its highly metastatic progression, remains among the most lethal of malignancies. Lung cancer metastasis can selectively spread to multiple different organs, however the genetic and molecular drivers for this process are still poorly understood. Understanding the heterogeneous genomic profile of lung cancer metastases is considered key in identifying therapeutic targets that prevent its spread. Research has identified the key source for metastasis being clusters of cells rather than individual cancer cells. These clusters, known as metastatic cancer cell clusters (MCCCs) have been shown to be 100-fold more tumorigenic than individual cancer cells. Unfortunately, access to these primary drivers of metastases remains difficult and has limited our understanding of their molecular and genomic profiles. Strong evidence in the literature suggests that differentially regulated biological pathways in MCCCs can provide new therapeutic drug targets to help combat cancer metastases. In order to expand research into MCCCs and their role in metastasis, we demonstrate a novel, proof of principle technology, to capture MCCCs directly from patients' whole blood. Our platform can be readily tuned for different solid tumor types by combining a biomimicry-based margination effect coupled with immunoaffinity to isolate MCCCs. Adopting a selective capture approach based on overexpressed CD44 in MCCCs provides a methodology that preferentially isolates them from whole blood. Furthermore, we demonstrate a high capture efficiency of more than 90% when spiking MCCC-like model cell clusters into whole blood. Characterization of the captured MCCCs from lung cancer patients by immunofluorescence staining and genomic analyses, suggests highly differential morphologies and genomic profiles. This study lays the foundation to identify potential drug targets thus unlocking a new area of anti-metastatic therapeutics.

Olivar: towards automated variant aware primer design for multiplex tiled amplicon sequencing of pathogens

Tiled amplicon sequencing has served as an essential tool for tracking the spread and evolution of pathogens. Over 15 million complete SARS-CoV-2 genomes are now publicly available, most sequenced and assembled via tiled amplicon sequencing. While computational tools for tiled amplicon design exist, they require downstream manual optimization both computationally and experimentally, which is slow and costly. Here we present Olivar, a first step towards a fully automated, variant-aware design of tiled amplicons for pathogen genomes. Olivar converts each nucleotide of the target genome into a numeric risk score, capturing undesired sequence features that should be avoided. In a direct comparison with PrimalScheme, we show that Olivar has fewer mismatches overlapping with primers and predicted PCR byproducts. We also compare Olivar head-to-head with ARTIC v4.1, the most widely used primer set for SARS-CoV-2 sequencing, and show Olivar yields similar read mapping rates (~90%) and better coverage to the manually designed ARTIC v4.1 amplicons. We also evaluate Olivar on real wastewater samples and found that Olivar has up to 3-fold higher mapping rates while retaining similar coverage. In summary, Olivar automates and accelerates the generation of tiled amplicons, even in situations of high mutation frequency and/or density. Olivar is available online as a web application at https://olivar.rice.edu and can be installed locally as a command line tool with Bioconda. Source code, installation guide, and usage are available at https://github.com/treangenlab/Olivar.