Sequence Elements Research Articles

CCAAT Promoter element regulates transgenerational expression of the MHC class I gene

Transgenerational gene expression depends on both underlying DNA sequences and epigenetic modifications. The latter, which can result in transmission of variegated gene expression patterns across multiple generations without DNA alterations, has been termed epigenetic inheritance and has been documented in plants, worms, flies and mammals. Whereas transcription factors binding to cognate DNA sequence elements regulate gene expression, the molecular basis for epigenetic inheritance has been linked to histone and DNA modifications and non-coding RNA. Here we report that mutation of the CCAAT box promoter element abrogates NF-Y binding and disrupts the stable transgenerational expression of an MHC class I transgene. Transgenic mice with a mutated CCAAT box in the MHC class I transgene display variegated expression of the transgene among littermates and progeny in multiple independently derived transgenic lines. After 4 generations, CCAAT mutant transgenic lines derived from a single founder stably displayed distinct patterns of expression. Histone modifications and RNA polymerase II binding correlate with expression of CCAAT mutant transgenic lines, whereas DNA methylation and nucleosome occupancy do not. Mutation of the CCAAT box also results in changes to CTCF binding and DNA looping patterns across the transgene that correlate with expression status. These studies identify the CCAAT promoter element as a regulator of stable transgenerational gene expression such that mutation of the CCAAT box results in variegated transgenerational inheritance. Considering that the CCAAT box is present in 30% of eukaryotic promoters, this study provides insights into how fidelity of gene expression patterns is maintained through multiple generations.

IProL: identifying DNA promoters from sequence information based on Longformer pre-trained model

Promoters are essential elements of DNA sequence, usually located in the immediate region of the gene transcription start sites, and play a critical role in the regulation of gene transcription. Its importance in molecular biology and genetics has attracted the research interest of researchers, and it has become a consensus to seek a computational method to efficiently identify promoters. Still, existing methods suffer from imbalanced recognition capabilities for positive and negative samples, and their recognition effect can still be further improved. We conducted research on E. coli promoters and proposed a more advanced prediction model, iProL, based on the Longformer pre-trained model in the field of natural language processing. iProL does not rely on prior biological knowledge but simply uses promoter DNA sequences as plain text to identify promoters. It also combines one-dimensional convolutional neural networks and bidirectional long short-term memory to extract both local and global features. Experimental results show that iProL has a more balanced and superior performance than currently published methods. Additionally, we constructed a novel independent test set following the previous specification and compared iProL with three existing methods on this independent test set.

Double-Helical Tiled Chain Structure of the Twist-Bend Liquid Crystal Phase in CB7CB

The twist-bend nematic liquid crystal phase is a three-dimensional fluid in which achiral bent molecules spontaneously form an orientationally ordered, macroscopically chiral, heliconical winding of a ten nanometer-scale pitch in the absence of positional ordering. Here, the structure of the twist-bend phase of the bent dimer CB7CB and its mixtures with 5CB is characterized, revealing a hidden invariance of the self-assembly of the twist-bend structure of CB7CB, such that over a wide range of concentrations and temperatures, the helix pitch and cone angle change as if the ground state for a pitch of the TB helix is an inextensible heliconical ribbon along the contour formed by following the local molecular long axis (the director). Remarkably, the distance along the length for a single turn of this helix is given by 2πRmol, where Rmol is the radius of bend curvature of a single all-trans CB7CB molecule. This relationship emerges from frustrated steric packing due to the bent molecular shape: space in the fluid that is hard to fill attracts the most flexible molecular subcomponents, a theme of nanosegregation that generates self-assembled, oligomer-like correlations of interlocking bent molecules in the form of a brickwork-like tiling of pairs of molecular strands into duplex double-helical chains. At higher temperatures in the twist-bend phase, the cone angle is small, the director contour is nearly along the helix axis z, and the duplex chains are sequences of biaxial elements formed by overlapping half-molecule pairs, with an approximately 45° rotation of the biaxis between each such element along the chain.

A comprehensive analysis of 3'UTRs in Caenorhabditis elegans.

3'Untranslated regions (3'UTRs) are essential portions of genes containing elements necessary for pre-mRNA 3'end processing and are involved in post-transcriptional gene regulation. Despite their importance, they remain poorly characterized in eukaryotes. Here, we have used a multi-pronged approach to extract and curate 3'UTR data from 11533 publicly available datasets, corresponding to the entire collection of Caenorhabditis elegans transcriptomes stored in the NCBI repository from 2009 to 2023. We have also performed high throughput cloning pipelines to identify and validate rare 3'UTR isoforms and incorporated and manually curated 3'UTR isoforms from previously published datasets. This updated C. elegans 3'UTRome (v3) is the most comprehensive resource in any metazoan to date, covering 97.4% of the 20362 experimentally validated protein-coding genes with refined and updated 3'UTR boundaries for 234893'UTR isoforms. We also used this novel dataset to identify and characterize sequence elements involved in pre-mRNA 3'end processing and update miRNA target predictions. This resource provides important insights into the 3'UTR formation, function, and regulation in eukaryotes.

Global evolutionary dynamics of virulence genes in ST11-KL47 carbapenem-resistant Klebsiella pneumoniae

ST11-KL47 is a hypervirulent carbapenem-resistant Klebsiella pneumoniae (CRKP) that is highly prevalent in China and poses a major public health risk. To investigate the evolutionary dynamics of virulence genes in this subclone, we analysed 78 sequenced isolates obtained from a long-term study across 29 centres from 17 cities in China. Virulence genes were located in large hybrid pNDM-Mar-like plasmids (length: ∼266 kilobases) rather than in classical pK2044-like plasmids. These hybrid plasmids, derived from the fusion of pK2044 and pNDM-Mar plasmids mediated by insertion sequence (IS) elements (such as ISKpn28 and IS26), integrated virulence gene fragments into the chromosome. Analysis of 217 sequences containing the special IncFIB (pNDM-Mar) replicon using public databases indicated that these plasmids typically contained T4SS-related and multiple antimicrobial resistance genes, were present in 24 countries, and were found in humans, animals, and the environment. Notably, the chromosomal integration of virulence genes was observed in strains across five countries across two continents. In vivo and in vitro models showed that the large hybrid plasmid increased the host fitness cost while increasing virulence. Conversely, virulence genes transferred to chromosomes resulted in increased fitness and lower virulence. In conclusion, virulence genes in the plasmids of ST11-KL47 CRKP are evolving, driven by adaptive negative selection, to enable vertical chromosomal inheritance along with conferring a survival advantage and low pathogenicity.

Comparative genomics of Deinococcus radiodurans: unveiling genetic discrepancies between ATCC 13939K and BAA-816 strains.

The Deinococcus genus is renowned for its remarkable resilience against environmental stresses, including ionizing radiation, desiccation, and oxidative damage. This resilience is attributed to its sophisticated DNA repair mechanisms and robust defense systems, enabling it to recover from extensive damage and thrive under extreme conditions. Central to Deinococcus research, the D. radiodurans strains ATCC BAA-816 and ATCC 13939 facilitate extensive studies into this remarkably resilient genus. This study focused on delineating genetic discrepancies between these strains by sequencing our laboratory's ATCC 13939 specimen (ATCC 13939K) and juxtaposing it with ATCC BAA-816. We uncovered 436 DNA sequence differences within ATCC 13939K, including 100 single nucleotide variations, 278 insertions, and 58 deletions, which could induce frameshifts altering protein-coding genes. Gene annotation revisions accounting for gene fusions and the reconciliation of gene lengths uncovered novel protein-coding genes and refined the functional categorizations of established ones. Additionally, the analysis pointed out genome structural variations due to insertion sequence (IS) elements, underscoring the D. radiodurans genome's plasticity. Notably, ATCC 13939K exhibited a loss of six ISDra2 elements relative to BAA-816, restoring genes fragmented by ISDra2, such as those encoding for α/β hydrolase and serine protease, and revealing new open reading frames, including genes imperative for acetoin decomposition. This comparative genomic study offers vital insights into the metabolic capabilities and resilience strategies of D. radiodurans.

Novel Loci (EIF4A2, ADIPOQ, TPRG1) for Triglyceride / High-density Lipoprotein Cholesterol Ratio Longitudinal Change (ΔTHR) among Subjects without Type 2 Diabetes: Evidence from the Long Life Family Study (LLFS) and the Framingham Heart Study (FHS) Offspring Cohort (OS).

Triglyceride (TG) /High density lipoprotein cholesterol (HDL-C) ratio (THR) represents a single surrogate predictor of hyperinsulinemia or insulin resistance that is associated with premature aging processes, risk of diabetes and increased mortality. To identify novel genetic loci for THR change over time (ΔTHR), we conducted genome-wide association study (GWAS) and genome-wide linkage scan (GWLS) among subjects of European ancestry who had complete data from two exams collected about seven years apart from the Long Life Family Study (LLFS, n=1384), a study with familial clustering of exceptional longevity in the US and Denmark. Subjects with diabetes or using medications for dyslipidemia were excluded from this analysis. ΔTHR was derived using growth curve modeling, and adjusted for age, sex, field centers, and principal components (PCs). GWAS was conducted using a linear mixed model accounted for familial relatedness. Our linkage scan was built on haplotype-based IBD estimation with 0.5 cM average spacing. Heritability of ΔTHR was moderate (46%). Our GWAS identified a significant locus at the LPL (p=1.58e-9) for ΔTHR; this gene locus has been reported before influencing baseline THR levels. Our GWLS found evidence for a significant linkage with a logarithm of the odds (LODs) exceeding 3 on 3q28 (LODs=4.1). Using a subset of 25 linkage enriched families (pedigree-specific LODs>0.1), we assessed sequence elements under 3q28 and identified two novel variants (EIF4A2/ADIPOQ-rs114108468, p=5e-6, MAF=1.8%; TPRG1-rs16864075, p=3e-6, MAF=8%; accounted for ~28% and ~29% of the linkage, respectively, and 57% jointly). While the former variant was associated with EIF4A2 (p=7e-5) / ADIPOQ (p=3.49e-2) RNA transcriptional levels, the latter variant was not associated with TPRG1 (p=0.23) RNA transcriptional levels. Replication in FHS OS observed modest effect of these loci on ΔTHR. Of 188 metabolites from 13 compound classes assayed in LLFS, we observed multiple metabolites (e.g., DG.38.5, PE.36.4, TG.58.3) that were significantly associated with the variants (p<3e-4). our linkage-guided sequence analysis approach permitted our discovery of two novel gene variants EIF4A2/ADIPOQ-rs114108468 and TPRG1-rs16864075 on 3q28 for ΔTHR among subjects without diabetes selected for exceptional survival and healthy aging.

Matrix characterisation of the hierarchy FiC n of bivaluated logics

In this paper, we define the family FiC n := { Fi C n } n ∈ ω of logics by means of semantics of bivaluations. This family is a variant of the family of paraconsistent bivaluated Ci u n -logics, defined and studied by J. Ciuciura. We will show that every logic in FiC n can be alternatively defined by means of finite matrices. This result arises from the characterisation of the truth-values of the involved matrices (relative to each Fi C n -logic) as being specific finite sequences of elements of the set 2 := { 0 , 1 } . Moreover, we will show in this paper that such characterisation is related to the well-known standard Fibonacci Sequence, which is presented here by means of its binary expansion.

Artificially designed synthetic promoter for a high level of salt induction using a cis-engineering approach

This work aimed to design a synthetic salt-inducible promoter using a cis-engineering approach. The designed promoter (PS) comprises a minimal promoter sequence for basal-level expression and upstream cis-regulatory elements (CREs) from promoters of salinity-stress-induced genes. The copy number, spacer lengths, and locations of CREs were manually determined based on their occurrence within native promoters. The initial activity profile of the synthesized PS promoter in transiently transformed N. tabacum leaves shows a seven-fold, five-fold, and four-fold increase in reporter GUS activity under salt, drought, and abscisic acid stress, respectively, at the 24-h interval, compared to the constitutive CaMV35S promoter. Analysis of gus expression in stable Arabidopsis transformants showed that the PS promoter induces over a two-fold increase in expression under drought or abscisic acid stress and a five-fold increase under salt stress at 24- and 48-h intervals, compared to the CaMV35S promoter. The promoter PS exhibits higher and more sustained activity under salt, drought, and abscisic acid stress compared to the constitutive CaMV35S.

Sensitive aptasensing of ATP based on a PAM site-regulated CRISPR/Cas12a activation.

The combination of CRISPR/Cas12a and functional DNA provides the possibility of constructing biosensors for detecting non-nucleic-acid targets. In the current study, the duplex protospacer adjacent motif (PAM) in the activator of CRISPR/Cas12a was used as a molecular switch, and a sensitive adenosine triphosphate (ATP) detection biosensor was constructed using an allosteric probe-conjugated PAM site formation in hybridization chain reaction (HCR) integrated with the CRISPR/Cas12a system (APF-CRISPR). In the absence of ATP, an aptamer-containing probe (AP) is in a stem-loop structure, which blocks the initiation of HCR. In the presence of ATP, the structure of AP is changed upon ATP binding, resulting in the release of the HCR trigger strand and the production of long duplex DNA with many PAM sites. Since the presence of a duplex PAM site is crucial for triggering the cleavage activity of CRISPR/Cas12a, the ATP-dependent formation of the PAM site in HCR products can initiate the FQ-reporter cleavage, allowing ATP quantification by measuring the fluorescent signals. By optimizing the sequence elements and detection conditions, the aptasensor demonstrated superior detection performance. The limit of detection (LOD) of the assay was estimated to be 1.16 nM, where the standard deviation of the blank was calculated based on six repeated measurements. The dynamic range of the detection was 25-750 nM, and the whole workflow of the assay was approximately 60min. In addition, the reliability and practicability of the aptasensor were validated by comparing it with a commercially available chemiluminescence kit for ATP detection in serum. Due to its high sensitivity, specificity, and reliable performance, the APF-CRISPR holds great potential in bioanalytical studies for ATP detection. In addition, we have provided a proof-of-principle for constructing a CRISPR/Cas12a-based aptasensor, in which the PAM is utilized to regulate Cas12a cleavage activity.

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants.

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

On a unified approach to homogeneous second-order linear difference equations with constant coefficients and some applications

AbstractIn this article, we establish a new closed formula for the solution of homogeneous second-order linear difference equations with constant coefficients by using matrix theory. This, in turn, gives new closed formulas concerning all sequences of this type such as the Fibonacci and Lucas sequences. Next, we show the main advantage of our formula which is based on the fact that future calculations rely on the previous ones and this is true from any desired starting point. As applications, we show that Binet’s formula, in this case, is valid for negative integers as well. Finally, new summation formulas relating elements of such sequences. As a conclusion, we present a formula for the sum of squares of Chebychev polynomials of the first and second kind.

The Analysis, Description, and Examination of the Maize LAC Gene Family's Reaction to Abiotic and Biotic Stress.

Laccase (LAC) is a diverse group of genes found throughout the plant genome essential for plant growth and the response to stress by converting monolignin into intricate lignin formations. However, a comprehensive investigation of maize laccase has not yet been documented. A bioinformatics approach was utilized in this research to conduct a thorough examination of maize (Zea mays L.), resulting in the identification and categorization of 22 laccase genes (ZmLAC) into six subfamilies. The gene structure and motifs of each subgroup were largely consistent. The distribution of the 22 LAC genes was uneven among the maize chromosomes, with the exception of chromosome 9. The differentiation of the genes was based on fragment replication, and the differentiation time was about 33.37 million years ago. ZmLAC proteins are primarily acidic proteins. There are 18 cis-acting elements in the promoter sequences of the maize LAC gene family associated with growth and development, stress, hormones, light response, and stress response. The analysis of tissue-specific expression revealed a high expression of the maize LAC gene family prior to the V9 stage, with minimal expression at post-V9. Upon reviewing the RNA-seq information from the publicly available transcriptome, it was discovered that ZmLAC5, ZmLAC10, and ZmLAC17 exhibited significant expression levels when exposed to various biotic and abiotic stress factors, suggesting their crucial involvement in stress responses and potential value for further research. This study offers an understanding of the functions of the LAC genes in maize's response to biotic and abiotic stress, along with a theoretical basis for comprehending the molecular processes at play.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130–160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation.

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Universally Quantitative Band-Selective Pure Shift NMR Spectroscopy.

NMR spectroscopy is often described as a quantitative analytical technique. Strictly, only the simple pulse-acquire experiment is universally quantitative, but the poor signal resolution of the 1H NMR pulse-acquie experiment frequently complicates quantitative analysis. Pure shift NMR techniques provide higher resolution, by reducing signal overlap, but they are susceptible to a variety of sources of site-dependent signal loss. Here, we introduce a new method that corrects for signal loss from such sources in band-selective pure shift NMR experiments, by performing different numbers of iterations of the same pulse sequence elements before acquisition to allow extrapolation back to the loss-free signal. We apply this method to both interferogram and semi-realtime acquisition modes, obtaining integrals within 1% of those acquired from a pulse-acquire experiment for a three-component mixture.

A chromosome-level genome assembly of the disco clam, Ctenoides ales.

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

A high-resolution map of functional miR-181 response elements in the thymus reveals the role of coding sequence targeting and an alternative seed match.

MicroRNAs (miRNAs) are critical post-transcriptional regulators in many biological processes. They act by guiding RNA-induced silencing complexes to miRNA response elements (MREs) in target mRNAs, inducing translational inhibition and/or mRNA degradation. Functional MREs are expected to predominantly occur in the 3' untranslated region and involve perfect base-pairing of the miRNA seed. Here, we generate a high-resolution map of miR-181a/b-1 (miR-181) MREs to define the targeting rules of miR-181 in developing murine Tcells. By combining a multi-omics approach with computational high-resolution analyses, we uncover novel miR-181 targets and demonstrate that miR-181 acts predominantly through RNA destabilization. Importantly, we discover an alternative seed match and identify a distinct set of targets with repeat elements in the coding sequence which are targeted by miR-181 and mediate translational inhibition. In conclusion, deep profiling of MREs in primary cells is critical to expand physiologically relevant targetomes and establish context-dependent miRNA targeting rules.

Sequence reliance of the Drosophila context-dependent transcription factor CLAMP.

Despite binding similar cis elements in multiple locations, a single transcription factor (TF) often performs context-dependent functions at different loci. How factors integrate cis sequence and genomic context is still poorly understood and has implications for off-target effects in genetic engineering. The Drosophila context-dependent TF chromatin-linked adaptor for male-specific lethal proteins (CLAMP) targets similar GA-rich cis elements on the X-chromosome and at the histone gene locus but recruits very different, locus-specific factors. We discover that CLAMP leverages information from both cis element and local sequence to perform context-specific functions. Our observations imply the importance of other cues, including protein-protein interactions and the presence of additional cofactors.

A study at the wildlife-livestock interface unveils the potential of feral swine as a reservoir for extended-spectrum β-lactamase-producing Escherichia coli

Wildlife is known to serve as carriers and sources of antimicrobial resistance (AMR). Due to their unrestricted movements and behaviors, they can spread antimicrobial resistant bacteria among livestock, humans, and the environment, thereby accelerating the dissemination of AMR. Extended-spectrum β-lactamase (ESBL)-producing Enterobacteriaceae is one of major concerns threatening human and animal health, yet transmission mechanisms at the wildlife-livestock interface are not well understood. Here, we investigated the mechanisms of ESBL-producing bacteria spreading across various hosts, including cattle, feral swine, and coyotes in the same habitat range, as well as from environmental samples over a two-year period. We report a notable prevalence and clonal dissemination of ESBL-producing E. coli in feral swine and coyotes, suggesting their persistence and adaptation within wildlife hosts. In addition, in silico studies showed that horizontal gene transfer, mediated by conjugative plasmids and insertion sequences elements, may play a key role in spreading the ESBL genes among these bacteria. Furthermore, the shared gut resistome of cattle and feral swine suggests the dissemination of antibiotic resistance genes at the wildlife-livestock interface. Taken together, our results suggest that feral swine may serve as a reservoir of ESBL-producing E. coli.