Unravelling epigenetic mechanisms in Cerastoderma edule genome: a comparison of healthy and neoplastic cockles.

Cancer is a multifaceted genetic disease characterized by the acquisition of several essential hallmarks. Notably, certain cancers exhibit horizontal transmissibility, observed across mammalian species and diverse bivalves, the latter referred to as hemic neoplasia. Within this complex landscape, epigenetic mechanisms such as histone modifications and cytosine methylation emerge as fundamental contributors to the pathogenesis of these transmissible cancers. Our study delves into the epigenetic landscape of Cerastoderma edule, focusing on whole-genome methylation and hydroxymethylation profiles in heathy specimens and transmissible neoplasias by means of Nanopore long-read sequencing. Our results unveiled a global hypomethylation in the neoplastic specimens compared to their healthy counterparts, emphasizing the role of DNA methylation in these tumorigenic processes. Furthermore, we verified that intragenic CpG methylation positively correlated with gene expression, emphasizing its role in modulating transcription in healthy and neoplastic cockles, as also highlighted by some up-methylated oncogenic genes. Hydroxymethylation levels were significantly more elevated in the neoplastic samples, particularly within satellites and complex repeats, likely related to structural functions. Additionally, our analysis also revealed distinct methylation and activity patterns in retrotransposons, providing additional insights into bivalve neoplastic processes. Altogether, these findings contribute to understanding the epigenetic dynamics of bivalve neoplasias and shed light on the roles of DNA methylation and hydroxymethylation in tumorigenesis. Understanding these epigenetic alterations holds promise for advancing our broader understanding of cancer epigenetics.

A chromosome-scale fishing cat reference genome for the evaluation of potential germline risk variants.

The fishing cat, Prionailurus viverrinus, faces a population decline, increasing the importance of maintaining healthy zoo populations. Unfortunately, zoo-managed individuals currently face a high prevalence of transitional cell carcinoma (TCC), a form of bladder cancer. To investigate the genetics of inherited diseases among captive fishing cats, we present a chromosome-scale assembly, generate the pedigree of the zoo-managed population, reaffirm the close genetic relationship with the Asian leopard cat (Prionailurus bengalensis), and identify 7.4 million single nucleotide variants (SNVs) and 23,432 structural variants (SVs) from whole genome sequencing (WGS) data of healthy and TCC cats. Only BRCA2 was found to have a high recurrent number of missense mutations in fishing cats diagnosed with TCC when compared to inherited human cancer risk variants. These new fishing cat genomic resources will aid conservation efforts to improve their genetic fitness and enhance the comparative study of feline genomes.

Deciphering cell types by integrating scATAC-seq data with genome sequences.

The single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) technology provides insight into gene regulation and epigenetic heterogeneity at single-cell resolution, but cell annotation from scATAC-seq remains challenging due to high dimensionality and extreme sparsity within the data. Existing cell annotation methods mostly focus on the cell peak matrix without fully utilizing the underlying genomic sequence. Here we propose a method, SANGO, for accurate single-cell annotation by integrating genome sequences around the accessibility peaks within scATAC data. The genome sequences of peaks are encoded into low-dimensional embeddings, and then iteratively used to reconstruct the peak statistics of cells through a fully connected network. The learned weights are considered as regulatory modes to represent cells, and utilized to align the query cells and the annotated cells in the reference data through a graph transformer network for cell annotations. SANGO was demonstrated to consistently outperform competing methods on 55 paired scATAC-seq datasets across samples, platforms and tissues. SANGO was also shown to be able to detect unknown tumor cells through attention edge weights learned by the graph transformer. Moreover, from the annotated cells, we found cell-type-specific peaks that provide functional insights/biological signals through expression enrichment analysis, cis-regulatory chromatin interaction analysis and motif enrichment analysis.

Comprehensive profiling of L1 retrotransposons in mouse.

L1 elements are retrotransposons currently active in mammals. Although L1s are typically silenced in most normal tissues, elevated L1 expression is associated with a variety of conditions, including cancer, aging, infertility and neurological disease. These associations have raised interest in the mapping of human endogenous de novo L1 insertions, and a variety of methods have been developed for this purpose. Adapting these methods to mouse genomes would allow us to monitor endogenous in vivo L1 activity in controlled, experimental conditions using mouse disease models. Here, we use a modified version of transposon insertion profiling, called nanoTIPseq, to selectively enrich young mouse L1s. By linking this amplification step with nanopore sequencing, we identified >95% annotated L1s from C57BL/6 genomic DNA using only 200 000 sequencing reads. In the process, we discovered 82 unannotated L1 insertions from a single C57BL/6 genome. Most of these unannotated L1s were near repetitive sequence and were not found with short-read TIPseq. We used nanoTIPseq on individual mouse breast cancer cells and were able to identify the annotated and unannotated L1s, as well as new insertions specific to individual cells, providing proof of principle for using nanoTIPseq to interrogate retrotransposition activity at the single-cell level in vivo.

Emergence of enhancers at late DNA replicating regions.

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication.

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

On the genetic basis of tail-loss evolution in humans and apes.

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.

Autonomous transposons tune their sequences to ensure somatic suppression.

Transposable elements (TEs) are a major constituent of human genes, occupying approximately half of the intronic space. During pre-messenger RNA synthesis, intronic TEs are transcribed along with their host genes but rarely contribute to the final mRNA product because they are spliced out together with the intron and rapidly degraded. Paradoxically, TEs are an abundant source of RNA-processing signals through which they can create new introns1, and also functional2 or non-functional chimeric transcripts3. The rarity of these events implies the existence of a resilient splicing code that is able to suppress TE exonization without compromising host pre-mRNA processing. Here we show that SAFB proteins protect genome integrity by preventing retrotransposition of L1 elements while maintaining splicing integrity, via prevention of the exonization of previously integrated TEs. This unique dual role is possible because of L1's conserved adenosine-rich coding sequences that are bound by SAFB proteins. The suppressive activity of SAFB extends to tissue-specific, giant protein-coding cassette exons, nested genes and Tigger DNA transposons. Moreover, SAFB also suppresses LTR/ERV elements in species in which they are still active, such as mice and flies. A significant subset of splicing events suppressed by SAFB in somatic cells are activated in the testis, coinciding with low SAFB expression in postmeiotic spermatids. Reminiscent of the division of labour between innate and adaptive immune systems that fight external pathogens, our results uncover SAFB proteins as an RNA-based, pattern-guided, non-adaptive defence system against TEs in the soma, complementing the RNA-based, adaptive Piwi-interacting RNA pathway of the germline.

Sequencing of N6-methyl-deoxyadenosine at single-base resolution across the mammalian genome.

Although DNA N6-methyl-deoxyadenosine (6mA) is abundant in bacteria and protists, its presence and function in mammalian genomes have been less clear. We present Direct-Read 6mA sequencing (DR-6mA-seq), an antibody-independent method, to measure 6mA at base resolution. DR-6mA-seq employs a unique mutation-based strategy to reveal 6mA sites as misincorporation signatures without any chemical or enzymatic modulation of 6mA. We validated DR-6mA-seq through the successful mapping of the well-characterized G(6mA)TC motif in the E. coli DNA. As expected, when applying DR-6mA-seq to mammalian systems, we found that genomic DNA (gDNA) 6mA abundance is generally low in most mammalian tissues and cells; however, we did observe distinct gDNA 6mA sites in mouse testis and glioblastoma cells. DR-6mA-seq provides an enabling tool to detect 6mA at single-base resolution for a comprehensive understanding of DNA 6mA in eukaryotes.

Time-resolved, integrated analysis of clonally evolving genomes.

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes.

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Mouse genome rewriting and tailoring of three important disease loci.

Genetically engineered mouse models (GEMMs) help us to understand human pathologies and develop new therapies, yet faithfully recapitulating human diseases in mice is challenging. Advances in genomics have highlighted the importance of non-coding regulatory genome sequences, which control spatiotemporal gene expression patterns and splicing in many human diseases1,2. Including regulatory extensive genomic regions, which requires large-scale genome engineering, should enhance the quality of disease modelling. Existing methods set limits on the size and efficiency of DNA delivery, hampering the routine creation of highly informative models that we call genomically rewritten and tailored GEMMs (GREAT-GEMMs). Here we describe 'mammalian switching antibiotic resistance markers progressively for integration' (mSwAP-In), a method for efficient genome rewriting in mouse embryonic stem cells. We demonstrate the use of mSwAP-In for iterative genome rewriting of up to 115 kb of a tailored Trp53 locus, as well as for humanization of mice using 116 kb and 180 kb human ACE2 loci. The ACE2 model recapitulated human ACE2 expression patterns and splicing, and notably, presented milder symptoms when challenged with SARS-CoV-2 compared with the existing K18-hACE2 model, thus representing a more human-like model of infection. Finally, we demonstrated serial genome writing by humanizing mouse Tmprss2 biallelically in the ACE2 GREAT-GEMM, highlighting the versatility of mSwAP-In in genome writing.

Whole-genome doubling in tissues and tumors.

The overwhelming majority of proliferating somatic human cells are diploid, and this genomic state is typically maintained across successive cell divisions. However, failures in cell division can induce a whole-genome doubling (WGD) event, in which diploid cells transition to a tetraploid state. While some WGDs are developmentally programmed to produce nonproliferative tetraploid cells with specific cellular functions, unscheduled WGDs can be catastrophic: erroneously arising tetraploid cells are ill-equipped to cope with their doubled cellular and chromosomal content and quickly become genomically unstable and tumorigenic. Deciphering the genetics that underlie the genesis, physiology, and evolution of whole-genome doubled (WGD+) cells may therefore reveal therapeutic avenues to selectively eliminate pathological WGD+ cells.

Comparative y-Narratives Inscribed Throughout Genomes of Man and Dog.

We have until now focused solely on the non-coding, more precisely the non-protein-coding (npc), part of DNA of man and dog in the search for hidden y-texts written by means of y-words - spelled by nucleotides A, C, G, and T and delimited by stop-codons. In this paper we use the same methods to analyse the whole human and canine genome, but we divide the genome into the genetic part, the naturally occurring sequence of exons, and the non-protein-coding genome according to definitions. By use of the y-text-finder we determine the number of zipf-qualified and a-qualified texts hidden in each of these parts. We present the actual methods and procedures, and the results in twelve figures, six for Homo sapiens sapiens and six for Canis lupus familiaris. Results show that there are lots of y-texts in the genetic part of the genome just as there are in the npc-genome. There is even a non-negligible number of ?-texts hidden in the sequence of exons. In addition, we show how many genes we find included in or overlapping zipf-qualified and a-qualified y-texts in the one-stranded DNA of man and dog. We assume that all this information represents the cell's total ability to behave in all of life's situations and discuss briefly ?-text reading and disease aetiology; carcinogenesis are also discussed.

CRISPR/Cas-mediated genome editing in mice for the development of drug delivery mechanism.

BACKGROUND: To manipulate particular locations in the bacterial genome, researchers have recently resorted to a group of unique sequences in bacterial genomes that are responsible for safeguarding bacteria against bacteriophages. Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) are two such systems, each of which consists of an RNA component and an enzyme component. METHODS AND RESULTS: This review focuses primarily on how CRISPR/Cas9 technology can be used to make models to study human diseases in mice. Creating RNA molecules that direct endonucleases to a specific position in the genome are crucial for achieving a specific genetic modification. CRISPR/Cas9 technology has allowed scientists to edit the genome with greater precision than ever before. Researchers can use knock-in and knock-out methods to model human diseases such as Neurological, cardiovascular disease, and cancer. CONCLUSIONS: In terms of developing innovative methods to discover ailments for diseases/disorders, improved CRISPR/Cas9 technology will provide easier access to valuable novel animal models.

Genome-wide in silico analysis leads to identification of deleterious L290V mutation in RBBP5 gene in Bos indicus.

Genome-wide deleterious mutations were identified in zebu cattle (Bos indicus) using in silico approach. The ddRAD sequence data of Sahiwal cattle were annotated and aligned with the cattle reference genome (ARS-UCD1.2). A total of 279,383 SNPs were identified at Read Depth10, which were further filtered to 692 missense SNPs. These SNPs were further analyzed, for functional consequences, by using Variant Effect Predictor, PolyPhen, PROVEAN, and PANTHER tools. A total of 18 SNPs, were finally identified as deleterious, and among these, 12 SNPs were mapped on nine different genes. ERRAT, ProSA-web, Project HOPE, TM-Align, and YASSARA tools, further confirmed the protein malfunctioning of one missense (L290V) mutation of Retinoblastoma binding protein-5 (RBBP5) gene, transcribing a cell cycle regulatory protein and associated with Retinoblastoma in human. This derived bioinformatics pipeline may be useful for preliminarily identifying the deleterious DNA mutations in livestock, specifically in absence of any genetic disease records.

A new method to synthesize multiple gRNA libraries and functional mapping of mammalian H3K4me3 regions.

Genetic screening based on the clustered regularly interspaced palindromic repeat (CRISPR) system has been indicated to be a powerful tool for identifying regulatory genes or cis-elements. However, when applying CRISPR screens to pinpoint functional elements at particular loci, a large number of guide RNA (gRNA) spacers may be required to achieve saturated coverage. Here, we present a controlled template-dependent elongation (CTDE) method relying on reversible terminators to synthesize gRNA libraries with genomic regions of interest. By applying this approach to H3K4me3 chromatin immunoprecipitation (ChIP)-derived DNA of mammalian cells, mega-sized gRNA libraries were synthesized in a tissue-specific manner, with which we conducted screening experiments to annotate essential sites for cell proliferation. Additionally, we confirmed that an essential site within the intron of LINC00339 regulates its own mRNA and that LINC00339 is a novel regulator of the cell cycle that maintains HepG2 proliferation. The CTDE method has the potential to be automated with high efficiency at low cost, and will be widely used to identify functional elements in mammalian genomes.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs.

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed's strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

Genome-wide identification glutathione-S-transferase gene superfamily in Daphnia pulex and its transcriptional response to nanoplastics.

Glutathione S-transferases (GSTs) are key multifunctional phase II detoxification enzymes involved in the regulation of growth, development, and stress responses. However, the knowledge of GSTs in the model invertebrate organism Daphnia pulex at the genomic level remains limited. In the present study, 35 GST genes were identified in D. pulex (Dp-GST), belonging to eight subfamilies, with the sigma, mu, and delta/epsilon subfamilies constituting approximately 29 %, 20 %, and 20 % of the GST superfamily, respectively. Chromosome tandem duplication of genes within the same subfamily was observed, which may be the main force driving GST expansion in D. pulex. DpGST genes showed different expression patterns in response to nanoplastic exposure for 96 h and 21 days. Some homologous GST genes in D. pulex showed similar expression patterns in response to nanoplastic exposure, likely owing to their unique motifs. For example, motif 9 is found in all delta/epsilon GST genes, whereas motifs 1, 2, 3, 5, and 7 are highly conserved in sigma GST genes. The characterization of D. pulex GSTs extending the knowledge of GST-mediated environmental contaminants, especially nanoplastics.

Generation and analysis of context-specific genome-scale metabolic models derived from single-cell RNA-Seq data.

Single-cell RNA sequencing combined with genome-scale metabolic models (GEMs) has the potential to unravel the differences in metabolism across both cell types and cell states but requires new computational methods. Here, we present a method for generating cell-type-specific genome-scale models from clusters of single-cell RNA-Seq profiles. Specifically, we developed a method to estimate the minimum number of cells required to pool to obtain stable models, a bootstrapping strategy for estimating statistical inference, and a faster version of the task-driven integrative network inference for tissues algorithm for generating context-specific GEMs. In addition, we evaluated the effect of different RNA-Seq normalization methods on model topology and differences in models generated from single-cell and bulk RNA-Seq data. We applied our methods on data from mouse cortex neurons and cells from the tumor microenvironment of lung cancer and in both cases found that almost every cell subtype had a unique metabolic profile. In addition, our approach was able to detect cancer-associated metabolic differences between cancer cells and healthy cells, showcasing its utility. We also contextualized models from 202 single-cell clusters across 19 human organs using data from Human Protein Atlas and made these available in the web portal Metabolic Atlas, thereby providing a valuable resource to the scientific community. With the ever-increasing availability of single-cell RNA-Seq datasets and continuously improved GEMs, their combination holds promise to become an important approach in the study of human metabolism.