Genome Evolution Research Articles

Actinomycetota bioprospecting from ore-forming environments.

Natural products from Actinomycetota have served as inspiration for many clinically relevant therapeutics. Despite early triumphs in natural product discovery, the rate of unearthing new compounds has decreased, necessitating inventive approaches. One promising strategy is to explore environments where survival is challenging. These harsh environments are hypothesized to lead to bacteria developing chemical adaptations (e.g. natural products) to enable their survival. This investigation focuses on ore-forming environments, particularly fluoride mines, which typically have extreme pH, salinity and nutrient scarcity. Herein, we have utilized metagenomics, metabolomics and evolutionary genome mining to dissect the biodiversity and metabolism in these harsh environments. This work has unveiled the promising biosynthetic potential of these bacteria and has demonstrated their ability to produce bioactive secondary metabolites. This research constitutes a pioneering endeavour in bioprospection within fluoride mining regions, providing insights into uncharted microbial ecosystems and their previously unexplored natural products.

Single-cell RNA-seq of Drosophila miranda testis reveals the evolution and trajectory of germline sex chromosome regulation.

Although sex chromosomes have evolved from autosomes, they often have unusual regulatory regimes that are sex- and cell-type-specific such as dosage compensation (DC) and meiotic sex chromosome inactivation (MSCI). The molecular mechanisms and evolutionary forces driving these unique transcriptional programs are critical for genome evolution but have been, in the case of MSCI in Drosophila, subject to continuous debate. Here, we take advantage of the younger sex chromosomes in D. miranda (XR and the neo-X) to infer how former autosomes acquire sex-chromosome-specific regulatory programs using single-cell and bulk RNA sequencing and ribosome profiling, in a comparative evolutionary context. We show that contrary to mammals and worms, the X down-regulation through germline progression is most consistent with the shutdown of DC instead of MSCI, resulting in half gene dosage at the end of meiosis for all 3 X's. Moreover, lowly expressed germline and meiotic genes on the neo-X are ancestrally lowly expressed, instead of acquired suppression after sex linkage. For the young neo-X, DC is incomplete across all tissue and cell types and this dosage imbalance is rescued by contributions from Y-linked gametologs which produce transcripts that are translated to compensate both gene and protein dosage. We find an excess of previously autosomal testis genes becoming Y-specific, showing that the neo-Y and its masculinization likely resolve sexual antagonism. Multicopy neo-sex genes are predominantly expressed during meiotic stages of spermatogenesis, consistent with their amplification being driven to interfere with mendelian segregation. Altogether, this study reveals germline regulation of evolving sex chromosomes and elucidates the consequences these unique regulatory mechanisms have on the evolution of sex chromosome architecture.

Insights into chloroplast genome evolution in Rutaceae through population genomics

Chloroplast genomes, pivotal for understanding plant evolution, remain unexplored in Rutaceae, a family with key perennial crops like citrus. Leveraging next-generation sequencing data from 509 Rutaceae accessions across 15 species, we conducted a de novo assembly of 343 chloroplast genomes, unveiling a chloroplast variation map highlighting the heterogeneous evolution rates across genome regions. Notably, differences in chloroplast genome size primarily originate from large single-copy and small single-copy regions. Structural variants predominantly occurred in the single-copy region, with two insertions located at the single-copy and inverted repeat region boundary. Phylogenetic analysis, principal component analysis, and population genetic statistics confirmed the cohesive clustering of different Citrus species, reflecting evolutionary dynamics in Citrus diversification. Furthermore, a close chloroplast genetic affinity was revealed among Atalantia (previously regarded as primitive citrus), Clausena, and Murraya. Zanthoxylum formed a distinct group with heightened genetic diversity. Through expanding our analysis to include 34 published chloroplast genomes, we explored chloroplast gene selection, revealing divergent evolutionary trends in photosynthetic pathways. While Photosystem I and Photosystem II exhibited robust negative selection, indicating stability, the Nicotinamide adenine dinucleotide (NADH) dehydrogenase pathway demonstrated rapid evolution, which was indicative of environmental adaptation. Finally, we discussed the effects of gene length and GC content on chloroplast gene evolution. In conclusion, our study reveals the genetic characterization of chloroplast genomes during Rutaceae diversification, providing insights into the evolutionary history of this family.

Growth and Genome Features of Non-O1/O139 Vibrio cholerae Isolated from Three Species of Common Freshwater Fish

Vibrio cholerae is the etiological agent of cholera in humans. The bacterium is frequently detected in aquatic products worldwide. However, the current literature on the genome evolution of V. cholerae of aquatic animal origins is limited. Here, we firstly characterized the growth and genome features of V. cholerae isolates with different resistance phenotypes from three species of common freshwater fish. The results revealed that the non-O1/O139 V. cholerae isolates (n = 4) were halophilic and grew optimally at 2% NaCl and pH 8.0. Their draft genome sequences were 3.89 Mb–4.15 Mb with an average GC content of 47.35–47.63%. Approximately 3366–3561 genes were predicted to encode proteins, but 14.9–17.3% of them were of an unknown function. A number of strain-specific genes (n = 221–311) were found in the four V. cholerae isolates, 3 of which belonged to none of any of the known sequence types (STs). Several putative mobile genetic elements (MGEs) existed in the V. cholerae isolates, including genomic islands (n = 4–9), prophages (n = 0–3), integrons (n = 1–1), and insertion sequences (n = 0–3). Notably, CRISPR-Cas system arrays (n = 2–10) were found in the V. cholerae genomes, whereby the potential immunity defense system could be active. Comparative genomic analyses also revealed many putative virulence-associated genes (n = 106–122) and antibiotic resistance-related genes (n = 6–9). Overall, the results of this study demonstrate the bacterial broader-spectrum growth traits and fill prior gaps in the genomes of V. cholerae originating from freshwater fish.

Investigating the role of HSP90 in cancer cell phenotypic plasticity.

"What are the mechanisms driving tumor evolution under the selective pressure of chemotherapeutics?" The emerging importance of epigenetic gene regulation in cancer progression necessitates not only our understanding of which genes are potential targets but also what mechanisms are employed in targeting those genes. Understanding the mechanisms that promote the evolution of the normal genome and epigenome is central to understanding how cancer cells adapt to chemotherapy. Our previous investigations have shown that heat shock protein 90 (HSP90) has a critical role in epigenetic gene regulation through histone acetylation and phenotypic plasticity. We recently extended these results in an A549 lung cancer model to test the role of HSP90 in the plasticity of cells regarding multi-drug resistance and epithelial-to-mesenchymal transition phenotypes. HSP90 is over-expressed in multiple cancers with poor prognosis. We propose that inhibition of HSP90 results in lower phenotypic plasticity of cancer cells making them more susceptible to chemotherapeutic intervention. Here we review the context of our results in the broader field of evolution of these phenotypes.

Discovery and genome-wide characterization of a novel miniature inverted repeat transposable element reveal genome-specific distribution in Glycine.

Miniature inverted repeat transposable elements (MITEs) are a dynamic component responsible for genome evolution. Tourist MITEs are derived from and mobilized by elements from the harbinger superfamily. In this study, a novel family of Tourist-like MITE was characterized in wild soybean species Glycine falcata. The new GftoMITE1 was initially discovered as an insertional polymorphism of the centromere-specific histone H3 (CenH3) gene in G. falcata. Using polymerase chain reaction, cloning and sequencing approaches, we showed a high number of copies of the GftoMITE1 family. Extensive bioinformatic analyses revealed the genome-level distribution and locus-specific mapping of GftoMITE1 members in Glycine species. Our results provide the first extensive characterization of the GftoMITE1 family and contribute to the understanding of the evolution of MITEs in the Glycine genus. Genome-specific GftoMITE1 was prominent in perennial wild soybean species, but not in annual cultivated soybean (Glycine max) or its progenitor (Glycine soja). We discuss that the GftoMITE1 family reveals a single rapid amplification in G. falcata and could have potential implications for gene regulation and soybean breeding as an efficient genetic marker for germplasm utilization in the future.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation.

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Long range segmentation of prokaryotic genomes by gene age and functionality.

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kilobases that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa’s climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa’s climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

The genome of the Arctic snow alga Limnomonas spitsbergensis (Chlamydomonadales).

Snow algae are a diverse group of extremophilic microeukaryotes found on melting polar and alpine snowfields. They play an important role in the microbial ecology of the cryosphere, and their propagation on snow and ice surfaces may in part accelerate climate-induced melting of these systems. High-quality snow algae genomes are needed for studies on their unique physiology, adaptive mechanisms, and genome evolution under multiple forms of stress, including cold temperatures and intense sunlight. Here, we assembled and annotated the genome of Limnomonas spitsbergensis, a cryophilic biciliate green alga originally isolated from melting snow on Svalbard, in the Arctic. The L. spitsbergensis genome assembly is based primarily on the use of PacBio long reads and secondly Illumina short reads, with an assembly size of 260.248 Mb in 124 contigs. A combination of 3 alternative annotation strategies was used including protein homology, RNA-seq evidence, and PacBio full-length transcript isoforms. The best merged set of annotations identified 18,277 protein-coding genes, which were 95.2% complete based on Benchmarking Universal Single-Copy Orthologs analysis. We also provide the annotated mitogenome, which is a relatively large 77.942 kb circular mapping sequence containing extensive repeats. The L. spitsbergensis genome will provide a new resource for research on snow algae adaptation, behavior, and natural selection in unique, low-temperature terrestrial environments that are under threat from climate change.

Shared Features Underlying Compact Genomes and Extreme Habitat Use in Chironomid Midges.

Nonbiting midges (family Chironomidae) are found throughout the world in a diverse array of aquatic and terrestrial habitats, can often tolerate harsh conditions such as hypoxia or desiccation, and have consistently compact genomes. Yet we know little about the shared molecular basis for these attributes and how they have evolved across the family. Here, we address these questions by first creating high-quality, annotated reference assemblies for Tanytarsus gracilentus (subfamily Chironominae, tribe Tanytarsini) and Parochlus steinenii (subfamily Podonominae). Using these and other publicly available assemblies, we created a time-calibrated phylogenomic tree for family Chironomidae with outgroups from order Diptera. We used this phylogeny to test for features associated with compact genomes, as well as examining patterns of gene family evolution and positive selection that may underlie chironomid habitat tolerances. Our results suggest that compact genomes evolved in the common ancestor of Chironomidae and Ceratopogonidae and that this occurred mainly through reductions in noncoding regions (introns, intergenic sequences, and repeat elements). Significantly expanded gene families in Chironomidae included biological processes that may relate to tolerance of stressful environments, such as temperature homeostasis, carbohydrate transport, melanization defense response, and trehalose transport. We identified several positively selected genes in Chironomidae, notably sulfonylurea receptor, CREB-binding protein, and protein kinase D. Our results improve our understanding of the evolution of small genomes and extreme habitat use in this widely distributed group.

Surfing the Waves of SARS-CoV-2: Analysis of Viral Genome Variants Using an NGS Survey in Verona, Italy.

The availability of new technologies for deep sequencing, including next-generation sequencing (NGS), allows for the detection of viral genome variations. The epidemiological determination of SARS-CoV-2 viral genome changes during the pandemic waves displayed the genome evolution and subsequent onset of variants over time. These variants were often associated with a different impact on viral transmission and disease severity. We investigated, in a retrospective study, the trend of SARS-CoV-2-positive samples collected from the start of the Italian pandemic (January 2020) to June 2023. In addition, viral RNAs extracted from 938 nasopharyngeal swab samples were analyzed using NGS between February 2022 and June 2023. Sequences were analyzed with bioinformatic tools to identify lineages and mutations and for phylogenetic studies. Six pandemic waves were detected. In our samples, we predominantly detected BA.2, BQ.1, BA.5.1, BA.5.2, and, more recently, XBB.1 and its subvariants. The data describe the SARS-CoV-2 genome evolution involved in viral interactions with the host and the dynamics of specific genome mutations and deletions.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Emx2 underlies the development and evolution of marsupial gliding membranes

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium—an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Elimination of subtelomeric repeat sequences exerts little effect on telomere essential functions in Saccharomyces cerevisiae

Telomeres, which are chromosomal end structures, play a crucial role in maintaining genome stability and integrity in eukaryotes. In the baker’s yeast Saccharomyces cerevisiae, the X- and Y’-elements are subtelomeric repetitive sequences found in all 32 and 17 telomeres, respectively. While the Y’-elements serve as a backup for telomere functions in cells lacking telomerase, the function of the X-elements remains unclear. This study utilized the S. cerevisiae strain SY12, which has three chromosomes and six telomeres, to investigate the role of X-elements (as well as Y’-elements) in telomere maintenance. Deletion of Y’-elements (SY12YΔ), X-elements (SY12XYΔ+Y), or both X- and Y’-elements (SY12XYΔ) did not impact the length of the terminal TG1-3 tracks or telomere silencing. However, inactivation of telomerase in SY12YΔ, SY12XYΔ+Y, and SY12XYΔ cells resulted in cellular senescence and the generation of survivors. These survivors either maintained their telomeres through homologous recombination-dependent TG1-3 track elongation or underwent microhomology-mediated intra-chromosomal end-to-end joining. Our findings indicate the non-essential role of subtelomeric X- and Y’-elements in telomere regulation in both telomerase-proficient and telomerase-null cells and suggest that these elements may represent remnants of S. cerevisiae genome evolution. Furthermore, strains with fewer or no subtelomeric elements exhibit more concise telomere structures and offer potential models for future studies in telomere biology.

Elimination of subtelomeric repeat sequences exerts little effect on telomere essential functions in Saccharomyces cerevisiae.

Telomeres, which are chromosomal end structures, play a crucial role in maintaining genome stability and integrity in eukaryotes. In the baker's yeast Saccharomyces cerevisiae, the X- and Y'-elements are subtelomeric repetitive sequences found in all 32 and 17 telomeres, respectively. While the Y'-elements serve as a backup for telomere functions in cells lacking telomerase, the function of the X-elements remains unclear. This study utilized the S. cerevisiae strain SY12, which has three chromosomes and six telomeres, to investigate the role of X-elements (as well as Y'-elements) in telomere maintenance. Deletion of Y'-elements (SY12YΔ), X-elements (SY12XYΔ+Y), or both X- and Y'-elements (SY12XYΔ) did not impact the length of the terminal TG1-3 tracks or telomere silencing. However, inactivation of telomerase in SY12YΔ, SY12XYΔ+Y, and SY12XYΔ cells resulted in cellular senescence and the generation of survivors. These survivors either maintained their telomeres through homologous recombination-dependent TG1-3 track elongation or underwent microhomology-mediated intra-chromosomal end-to-end joining. Our findings indicate the non-essential role of subtelomeric X- and Y'-elements in telomere regulation in both telomerase-proficient and telomerase-null cells and suggest that these elements may represent remnants of S. cerevisiae genome evolution. Furthermore, strains with fewer or no subtelomeric elements exhibit more concise telomere structures and offer potential models for future studies in telomere biology.

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.

Pangenome analysis reveals transposon-driven genome evolution in cotton

BackgroundTransposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution.ResultsHere, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes.ConclusionsOur pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

The Draft Genome Sequences of 50 Salamander species (Caudata, Amphibia).

We present partial genome sequences of 50 salamander species (Urodela) from 10 genera and 4 families. These span nearly the entire range of genome sizes in salamanders, from ~14-130GB, the latter of which is among the largest of all animal genomes. Only three salamander genomes were available to this point, from Ambystomatidae (one species) and Salamandridae (two species from two genera), to which we have added Amphiumidae (one species), Plethodontidae (45 species from 6 genera), Proteidae (one species), and Sirenidae (three species from two genera). These span ~140 million years of evolutionary divergence, leaving only Cryptobranchidae, Hynobiidae, and Rhyacotritonidae as salamander families without genome assemblies. These data should facilitate additional future work on speciation and genome evolution, both within Urodela and across Animalia.

Deep learning model for predicting genetic diseases using DNA sequence data

The emergence of machine learning in the recent decade has excelled in determining new potential features and nonlinear relationships existing between the data derived from the DNA sequences of genetic diseases. Machine learning also enhances the process of handling data with maximum predicted variables compared to observations during the data mining process of prediction. In this context, our study presents a deep learning model for predicting Transcription Factor Binding Sites (TFBS) in DNA sequences, with a focus on features within genetic data associated with diseases. Transcription Factors (TFs) play a crucial role in modulating gene expression by binding to TFBS. The accurate prediction of TFBS is essential for understanding genome function and evolution. Thus, we develop an efficient deep learning model that considers TFBS prediction as a nucleotide-level binary classification task. In our proposed model, first we create an input matrix using the original DNA sequences. Next, we encode these DNA sequences using one-hot encoding, representing them as a sequence of numerical values. We then employ three convolutional layers, allowing our model to capture intricate patterns and motif features over a larger spatial range. To capture important features within the DNA sequence and to focus on them, we incorporate an attention layer. Finally, a dense layer, consisting of two fully connected layers and a dropout layer, calculates the probability of TF binding site occurrence based on the features learned by the proposed model. Our experimental results, using in-vivo datasets obtained from Chip-seq, demonstrate the superior performance of our proposed deep learning model in TFBS prediction compared to other existing state-of-the-art methods. The improvement in accuracy is due to additional layers of CNN and then an attention layer in the model. Thus, this result in a better performance of our approach in predicting the transcription factor binding sites and enhancing our understanding of gene regulation and genome function.