A long-standing question concerns the role of Z-DNA in transcription. Here we use a deep learning approach DeepZ that predicts Z-flipons based on DNA sequence, structural properties of nucleotides and omics data. We examined Z-flipons that are conserved between human and mouse genomes after generating whole-genome Z-flipon maps and then validated them by orthogonal approaches based on high resolution chemical mapping of Z-DNA and the transformer algorithm Z-DNABERT. For human and mouse, we revealed similar pattern of transcription factors, chromatin remodelers, and histone marks associated with conserved Z-flipons. We found significant enrichment of Z-flipons in alternative and bidirectional promoters associated with neurogenesis genes. We show that conserved Z-flipons are associated with increased experimentally determined transcription reinitiation rates compared to promoters without Z-flipons, but without affecting elongation or pausing. Our findings support a model where Z-flipons engage Transcription Factor E and impact phenotype by enabling the reset of preinitiation complexes when active, and the suppression of gene expression when engaged by repressive chromatin complexes.

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

Kaposi\'s sarcoma herpesvirus (KSHV) ORF34 plays a significant role as a component of the viral pre-initiation complex (vPIC), which is indispensable for late gene expression across beta- and gammaherpesviruses. Although the key role of ORF34 within the vPIC and its function as a hub protein have been recognized, further clarification regarding its specific contribution to vPIC functionality and interactions with other components is required. This study employed a deep learning algorithm-assisted structural model of ORF34, revealing highly conserved amino acid residues across human beta- and gammaherpesviruses localized in structured domains. Thus, we engineered ORF34 alanine-scanning mutants by substituting conserved residues with alanine. These mutants were evaluated for their ability to interact with other vPIC factors and restore viral production in cells harboring the ORF34-deficient KSHV-BAC. Our experimental results highlight the crucial role of the four cysteine residues conserved in ORF34: a tetrahedral arrangement consisting of a pair of C-Xn-C consensus motifs. This suggests the potential incorporation of metal cations in interacting with ORF24 and ORF66 vPIC components, facilitating late gene transcription, and promoting overall virus production by capturing metal cations. In summary, our findings underline the essential role of conserved cysteines in KSHV ORF34 for effective vPIC assembly and viral replication, thereby enhancing our understanding of the complex interplay between the vPIC components.
OBJECTIVE: The initiation of late gene transcription is universally conserved across the beta- and gammaherpesvirus families. This process employs a viral pre-initiation complex (vPIC), which is analogous to a cellular PIC. Although KSHV ORF34 is a critical factor for viral replication and is a component of the vPIC, the specifics of vPIC formation and the essential domains crucial for its function remain unclear. Structural predictions suggest that the four conserved cysteines (C170, C175, C256, and C259) form a tetrahedron that coordinates the metal cation. We investigated the role of these conserved amino acids in interactions with other vPIC components, late gene expression, and virus production to demonstrate for the first time that these cysteines are pivotal for such functions. This discovery not only deepens our comprehensive understanding of ORF34 and vPIC dynamics but also lays the groundwork for more detailed studies on herpesvirus replication mechanisms in future research.

Coronavirus infectious disease 2019 (COVID-19), caused by severe acute respiratory virus type 2 (SARS-CoV-2), has caused a global public health crisis. As an RNA virus, the high gene mutability of SARS-CoV-2 poses significant challenges to the development of broad-spectrum vaccines and antiviral therapeutics. There remains a lack of specific therapeutics directly targeting SARS-CoV-2. With the ability to efficiently inhibit the expression of target genes in a sequence-specific way, small interfering RNA (siRNA) therapy has exhibited significant potential in antiviral and other disease treatments. In this work, we presented a highly effective self-assembled siRNA nanoparticle targeting multiple highly conserved regions of SARS-CoV-2. The siRNA sequences targeting viral conserved regions were first screened and evaluated by their thermodynamic features, off-target effects, and secondary structure toxicities. RNA motifs including siRNA sequences were then designed and self-assembled into siRNA nanoparticles. These siRNA nanoparticles demonstrated remarkable uniformity and stability and efficiently entered cells directly through cellular endocytic pathways. Moreover, these nanoparticles effectively inhibited the replication of SARS-CoV-2, exhibiting a superior inhibitory effect compared to free siRNA. These results demonstrated that these self-assembled siRNA nanoparticles targeting highly conserved regions of SARS-CoV-2 represent highly effective antiviral candidates for the treatment of infections, and are promisingly effective against current and future viral variants.

Puumala orthohantavirus (PUUV) is an emerging zoonotic virus endemic to Europe and Russia that causes nephropathia epidemica, a mild form of hemorrhagic fever with renal syndrome (HFRS). There are limited options for treatment and diagnosis of orthohantavirus infection, making the search for potential immunogenic candidates crucial. In the present work, various bioinformatics tools were employed to design conserved immunogenic peptides containing multiple epitopes of PUUV nucleocapsid protein. Eleven conserved peptides (90% conservancy) of the PUUV nucleocapsid protein were identified. Three conserved peptides containing multiple T and B cell epitopes were selected using a consensus epitope prediction algorithm. Molecular docking using the HPEP dock server demonstrated strong binding interactions between the epitopes and HLA molecules (ten alleles for each class I and II HLA). Moreover, an analysis of population coverage using the IEDB database revealed that the identified peptides have over 90% average population coverage across six continents. Molecular docking and simulation analysis reveal a stable interaction with peptide constructs of chosen immunogenic peptides and Toll-like receptor-4. These computational analyses demonstrate selected peptides\' immunogenic potential, which needs to be validated in different experimental systems.

Enterovirus genomic replication initiates at a predicted RNA cloverleaf (5\'CL) at the 5\' end of the RNA genome. The 5\'CL contains one stem (SA) and three stem-loops (SLB, SLC, SLD). Here, we present an analysis of 5\'CL conservation and divergence for 209 human health-related serotypes from the enterovirus genus, including enterovirus and rhinovirus species. Phylogenetic analysis indicates six distinct 5\'CL serotypes that only partially correlate with the species definition. Additional findings include that 5\'CL sequence conservation is higher between the EV species than between the RV species, the 5\'CL of EVA and EVB are nearly identical, and RVC has the lowest 5\'CL conservation. Regions of high conservation throughout all species include SA and the loop and nearby bases of SLB, which is consistent with known protein interactions at these sites. In addition to the known protein binding site for the Poly-C binding protein in the loop of SLB, other conserved consecutive cytosines in the stems of SLB and SLC provide additional potential interaction sites that have not yet been explored. Other sites of conservation, including the predicted bulge of SLD and other conserved stem, loop, and junction regions, are more difficult to explain and suggest additional interactions or structural requirements that are not yet fully understood. This more intricate understanding of sequence and structure conservation and variability in the 5\'CL may assist in the development of broad-spectrum antivirals against a wide range of enteroviruses, while better defining the range of virus isotypes expected to be affected by a particular antiviral.

Cytosolic sulfotransferases (SULTs) are Phase 2 drug-metabolizing enzymes that catalyze the conjugation of sulfonate to endogenous and xenobiotic compounds, increasing their hydrophilicity and excretion from cells. To date, 13 human SULTs have been identified and classified into five families. SULT4A1 mRNA encodes two variants: (1) the wild type, encoding a 284 amino acid, ~33 kDa protein, and (2) an alternative spliced variant resulting from a 126 bp insert between exon 6 and 7, which introduces a premature stop codon that enhances nonsense-mediated decay. SULT4A1 is classified as an SULT based on sequence and structural similarities, including PAPS-domains, active-site His, and the dimerization domain; however, the catalytic pocket lid \'Loop 3\' size is not conserved. SULT4A1 is uniquely expressed in the brain and localized in the cytosol and mitochondria. SULT4A1 is highly conserved, with rare intronic polymorphisms that have no outward manifestations. However, the SULT4A1 haplotype is correlated with Phelan-McDermid syndrome and schizophrenia. SULT4A1 knockdown revealed potential SULT4A1 functions in photoreceptor signaling and knockout mice display hampered neuronal development and behavior. Mouse and yeast models revealed that SULT4A1 protects the mitochondria from endogenously and exogenously induced oxidative stress and stimulates cell division, promoting dendritic spines\' formation and synaptic transmission. To date, no physiological enzymatic activity has been associated with SULT4A1.

Ultraconserved elements were discovered two decades ago, arbitrarily defined as sequences that are identical over a length ≥ 200 bp in the human, mouse, and rat genomes. The definition was subsequently extended to sequences ≥ 100 bp identical in at least three of five mammalian genomes (including dog and cow), and shown to have undergone rapid expansion from ancestors in fish and strong negative selection in birds and mammals. Since then, many more genomes have become available, allowing better definition and more thorough examination of ultraconserved element distribution and evolutionary history. We developed a fast and flexible analytical pipeline for identifying ultraconserved elements in multiple genomes, dedUCE, which allows manipulation of minimum length, sequence identity, and number of species with a detectable ultraconserved element according to specified parameters. We suggest an updated definition of ultraconserved elements as sequences ≥ 100 bp and ≥97% sequence identity in ≥50% of placental mammal orders (12,813 ultraconserved elements). By mapping ultraconserved elements to ∼200 species, we find that placental ultraconserved elements appeared early in vertebrate evolution, well before land colonization, suggesting that the evolutionary pressures driving ultraconserved element selection were present in aquatic environments in the Cambrian-Devonian periods. Most (>90%) ultraconserved elements likely appeared after the divergence of gnathostomes from jawless predecessors, were largely established in sequence identity by early Sarcopterygii evolution-before the divergence of lobe-finned fishes from tetrapods-and became near fixed in the amniotes. Ultraconserved elements are mainly located in the introns of protein-coding and noncoding genes involved in neurological and skeletomuscular development, enriched in regulatory elements, and dynamically expressed throughout embryonic development.

Magnetoreceptive biology as a field remains relatively obscure; compared with the breadth of species believed to sense magnetic fields, it remains under-studied. Here, we present grounds for the expansion of magnetoreception studies among teleosts. We begin with the electromagnetic perceptive gene (EPG) from Kryptopterus vitreolus and expand to identify 72 teleosts with homologous proteins containing a conserved three-phenylalanine (3F) motif. Phylogenetic analysis provides insight as to how EPG may have evolved over time and indicates that certain clades may have experienced a loss of function driven by different fitness pressures. One potential factor is water type with freshwater fish significantly more likely to possess the functional motif version (FFF), and saltwater fish to have the non-functional variant (FXF). It was also revealed that when the 3F motif from the homologue of Brachyhypopomus gauderio (B.g.) is inserted into EPG-EPG(B.g.)-the response (as indicated by increased intracellular calcium) is faster. This indicates that EPG has the potential to be engineered to improve upon its response and increase its utility to be used as a controller for specific outcomes.

Gene expression in Arabidopsis is regulated by more than 1,900 transcription factors (TFs), which have been identified genome-wide by the presence of well-conserved DNA-binding domains. Activator TFs contain activation domains (ADs) that recruit coactivator complexes; however, for nearly all Arabidopsis TFs, we lack knowledge about the presence, location and transcriptional strength of their ADs1. To address this gap, here we use a yeast library approach to experimentally identify Arabidopsis ADs on a proteome-wide scale, and find that more than half of the Arabidopsis TFs contain an AD. We annotate 1,553 ADs, the vast majority of which are, to our knowledge, previously unknown. Using the dataset generated, we develop a neural network to accurately predict ADs and to identify sequence features that are necessary to recruit coactivator complexes. We uncover six distinct combinations of sequence features that result in activation activity, providing a framework to interrogate the subfunctionalization of ADs. Furthermore, we identify ADs in the ancient AUXIN RESPONSE FACTOR family of TFs, revealing that AD positioning is conserved in distinct clades. Our findings provide a deep resource for understanding transcriptional activation, a framework for examining function in intrinsically disordered regions and a predictive model of ADs.