Sensory nerves unlock the TOLL-7 gate for cancer spread.

Cancers hijack the nervous system for growth and spread. Thus, disrupting neuron-cancer crosstalk holds promise for blocking metastasis. Recently, Padmanaban et al. reported new therapeutic targets and showed that breast cancer cells activate sensory neurons to secrete the neuropeptide substance P (SP), leading to single-strand (ss)RNA release and noncanonical Toll-like receptor (TLR)7 signaling that drives metastasis.

Diversity and impact of single-stranded RNA viruses in Czech Heterobasidion populations.

Heterobasidion annosum sensu lato comprises some of the most devastating pathogens of conifers. Exploring virocontrol as a potential strategy to mitigate economic losses caused by these fungi holds promise for the future. In this study, we conducted a comprehensive screening for viruses in 98 H. annosum s.l. specimens from different regions of Czechia aiming to identify viruses inducing hypovirulence. Initial examination for dsRNA presence was followed by RNA-seq analyses using pooled RNA libraries constructed from H. annosum and Heterobasidion parviporum, with diverse bioinformatic pipelines employed for virus discovery. Our study uncovered 25 distinct ssRNA viruses, including two ourmia-like viruses, one mitovirus, one fusarivirus, one tobamo-like virus, one cogu-like virus, one bisegmented narna-like virus and one segment of another narna-like virus, and 17 ambi-like viruses, for which hairpin and hammerhead ribozymes were detected. Coinfections of up to 10 viruses were observed in six Heterobasidion isolates, whereas another six harbored a single virus. Seventy-three percent of the isolates analyzed by RNA-seq were virus-free. These findings show that the virome of Heterobasidion populations in Czechia is highly diverse and differs from that in the boreal region. We further investigated the host effects of certain identified viruses through comparisons of the mycelial growth rate and proteomic analyses and found that certain tested viruses caused growth reductions of up to 22% and significant alterations in the host proteome profile. Their intraspecific transmission rates ranged from 0% to 33%. Further studies are needed to fully understand the biocontrol potential of these viruses in planta.IMPORTANCEHeterobasidion annosum sensu lato is a major pathogen causing significant damage to conifer forests, resulting in substantial economic losses. This study is significant as it explores the potential of using viruses (virocontrol) to combat these fungal pathogens. By identifying and characterizing a diverse array of viruses in H. annosum populations from Czechia, the research opens new avenues for biocontrol strategies. The discovery of 25 distinct ssRNA viruses, some of which reduce fungal growth and alter proteome profiles, suggests that these viruses could be harnessed to mitigate the impact of Heterobasidion. Understanding the interactions between these viruses and their fungal hosts is crucial for developing effective, environmentally friendly methods to protect conifer forests and maintain ecosystem health. This study lays the groundwork for future research on the application of mycoviruses in forest disease management.

Influenza A Vaccine Candidates Based on Virus-like Particles Formed by Coat Proteins of Single-Stranded RNA Phages Beihai32 and PQ465.

Efficient control of influenza A infection can potentially be achieved through the development of broad-spectrum recombinant vaccines based on conserved antigens. The extracellular domain of the transmembrane protein M2 of influenza A virus (M2e) is highly conserved but poorly immunogenic and needs to be fused to an adjuvant protein or carrier virus-like particles (VLPs) to increase immunogenicity and provide protection against infection. In this study, we obtained VLPs based on capsid proteins (CPs) of single-stranded RNA phages Beihai32 and PQ465 bearing the M2e peptides. Four copies of the M2e peptide were linked to the C-terminus of the CP of phage Beihai32 and to the N and C termini of the CP of phage PQ465. The hybrid proteins, being expressed in Escherichia coli, formed spherical VLPs of about 30 nm in size. Immunogold transmission electron microscopy showed that VLPs formed by the phage PQ465 CP with a C-terminal M2e fusion present the M2e peptide on the surface. Subcutaneous immunization of mice with VLPs formed by both CPs containing four copies of the M2e peptide at the C termini induced high levels of M2e-specific IgG antibodies in serum and provided mice with protection against lethal influenza A virus challenge. In the case of an N-terminal fusion of M2e with the phage PQ465 CP, the immune response against M2e was significantly lower. CPs of phages Beihai32 and PQ465, containing four copies of the M2e peptide at their C termini, can be used to develop recombinant influenza A vaccine.

CD4+ T cells reactive to a hybrid peptide from insulin-chromogranin A adopt a distinct effector fate and are pathogenic in autoimmune diabetes.

T cell-mediated islet destruction is a hallmark of autoimmune diabetes. Here, we examined the dynamics and pathogenicity of CD4+ T cell responses to four different insulin-derived epitopes during diabetes initiation in non-obese diabetic (NOD) mice. Single-cell RNA sequencing of tetramer-sorted CD4+ T cells from the pancreas revealed that islet-antigen-specific T cells adopted a wide variety of fates and required XCR1+ dendritic cells for their activation. Hybrid-insulin C-chromogranin A (InsC-ChgA)-specific CD4+ T cells skewed toward a distinct T helper type 1 (Th1) effector phenotype, whereas the majority of insulin B chain and hybrid-insulin C-islet amyloid polypeptide-specific CD4+ T cells exhibited a regulatory phenotype and early or weak Th1 phenotype, respectively. InsC-ChgA-specific CD4+ T cells were uniquely pathogenic upon transfer, and an anti-InsC-ChgA:IAg7 antibody prevented spontaneous diabetes. Our findings highlight the heterogeneity of T cell responses to insulin-derived epitopes in diabetes and argue for the feasibility of antigen-specific therapies that blunts the response of pathogenic CD4+ T cells causing autoimmunity.

Virus discovery in cultured portunid crabs: Genomic, phylogenetic, histopathological and microscopic characterization of a reovirus and a new bunyavirus.

Portunid crabs are distributed worldwide and highly valued in aquaculture. Viral infections are the main limiting factor for the survival of these animals and, consequently, for the success of commercial-scale cultivation. However, there is still a lack of knowledge about the viruses that infect cultured portunid crabs worldwide. Herein, the genome sequence and phylogeny of Callinectes sapidus reovirus 2 (CsRV2) are described, and the discovery of a new bunyavirus in Callinectes danae cultured in southern Brazil is reported. The CsRV2 genome sequence consists of 12 dsRNA segments (20,909 nt) encode 13 proteins. The predicted RNA-dependent RNA polymerase (RdRp) shows a high level of similarity with that of Eriocheir sinensis reovirus 905, suggesting that CsRV2 belongs to the genus Cardoreovirus. The CsRV2 particles are icosahedral, measuring approximately 65 nm in diameter, and exhibit typical non-turreted reovirus morphology. High throughput sequencing data revealed the presence of an additional putative virus genome similar to bunyavirus, called Callinectes danae Portunibunyavirus 1 (CdPBV1). The CdPBV1 genome is tripartite, consisting of 6,654 nt, 3,120 nt and 1,656 nt single-stranded RNA segments that each encode a single protein. Each segment has a high identity with European shore crab virus 1, suggesting that CdPBV1 is a new representative of the family Cruliviridae. The putative spherical particles of CdPBV1 measure ∼120 nm in diameter and present a typical bunyavirus morphology. The results of the histopathological analysis suggest that these new viruses can affect the health and, consequently, the survival of C. danae in captivity. Therefore, the findings reported here should be used to improve prophylactic and pathogen control practices and contribute to the development and optimization of the production of soft-shell crabs on a commercial scale in Brazil.

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction.

Booster doses are crucial against severe COVID-19, as rapid virus mutations and variant emergence prolong the pandemic crisis. The virus's quick evolution, short generation-time, and adaptive changes impact virulence and evolvability, helping predictions about variant of concerns' (VOCs') landscapes. Here, in this study, we used a new computational algorithm, to predict the mutational pattern in SARS-CoV-2 ssRNA, proteomics, structural identification, mutation stability, and functional correlation, as well as immune escape mechanisms. Interestingly, the sequence diversity of SARS Coronavirus-2 has demonstrated a predominance of G- > A and C- > U substitutions. The best validation statistics are explored here in seven homologous models of the expected mutant SARS-CoV-2 spike ssRNA and employed for hACE2 and IgG interactions. The interactome profile of SARS-CoV-2 spike with hACE2 and IgG revealed a strong correlation between phylogeny and divergence time. The systematic adaptation of SARS-CoV-2 spike ssRNA influences infectivity and immune escape. Data suggest higher propensity of Adenine rich sequence promotes MHC system avoidance, preferred by A-rich codons. Phylogenetic data revealed the evolution of SARS-CoV-2 lineages' epidemiology. Our findings may unveil processes governing the genesis of immune-resistant variants, prompting a critical reassessment of the coronavirus mutation rate and exploration of hypotheses beyond mechanical aspects.

Purification of Single-Stranded RNA Bacteriophages and Host Receptors for Structural Determination Using Cryo-Electron Microscopy.

Single-stranded RNA bacteriophages (ssRNA phages) are small viruses with a compact genome (~3-4 kb) that infect gram-negative bacteria via retractile pili. These phages have been applied in various fields since their discovery approximately 60 years ago. To understand their biology, it is crucial to analyze the structure of mature virions. Cryo-electron microscopy (cryo-EM) has been employed to determine the structures of two ssRNA phages, MS2 and Qß. This chapter presents a method for purifying these two phages and their receptor, the F-pilus, to allow examination using cryo-EM.

HIPPO: HIstogram-based Pseudo-POtential for scoring protein-ssRNA fragment-based docking poses.

BACKGROUND: The RNA-Recognition motif (RRM) is a protein domain that binds single-stranded RNA (ssRNA) and is present in as much as 2% of the human genome. Despite this important role in biology, RRM-ssRNA interactions are very challenging to study on the structural level because of the remarkable flexibility of ssRNA. In the absence of atomic-level experimental data, the only method able to predict the 3D structure of protein-ssRNA complexes with any degree of accuracy is ssRNA'TTRACT, an ssRNA fragment-based docking approach using ATTRACT. However, since ATTRACT parameters are not ssRNA-specific and were determined in 2010, there is substantial opportunity for enhancement. RESULTS: Here we present HIPPO, a composite RRM-ssRNA scoring potential derived analytically from contact frequencies in near-native versus non-native docking models. HIPPO consists of a consensus of four distinct potentials, each extracted from a distinct reference pool of protein-trinucleotide docking decoys. To score a docking pose with one potential, for each pair of RNA-protein coarse-grained bead types, each contact is awarded or penalised according to the relative frequencies of this contact distance range among the correct and incorrect poses of the reference pool. Validated on a fragment-based docking benchmark of 57 experimentally solved RRM-ssRNA complexes, HIPPO achieved a threefold or higher enrichment for half of the fragments, versus only a quarter with the ATTRACT scoring function. In particular, HIPPO drastically improved the chance of very high enrichment (12-fold or higher), a scenario where the incremental modelling of entire ssRNA chains from fragments becomes viable. However, for the latter result, more research is needed to make it directly practically applicable. Regardless, our approach already improves upon the state of the art in RRM-ssRNA modelling and is in principle extendable to other types of protein-nucleic acid interactions.

Antiviral discovery in toxic cyanobacteria: Low hanging fruit in the age of pandemics.

The power of novel vaccination technologies and their rapid development were elucidated clearly during the COVID-19 pandemic. At the same time, it also became clear that there is an urgent need to discover and manufacture new antivirals that target emerging viral threats. Toxic species of cyanobacteria produce a range of bioactive compounds that makes them good candidates for drug discovery. Nevertheless, few studies demonstrate the antiviral potential of cyanobacteria. This is partly due to the lack of specific and simple protocols designed for the rapid detection of antiviral activity in cyanobacteria and partly because specialized facilities for work with pathogenic viruses are few and far between. We therefore developed an easy method for the screening of cyanobacterial cultures for antiviral activity and used our private culture collection of non-pathogenic virus isolates to show that antiviral activity is a prominent feature in the cyanobacterium Microcystis aeruginosa. In this proof-of-concept study, we show that M. aeruginosa extracts from three different cyanobacterial strains delay infection of diatom-infecting single-stranded DNA and single-stranded RNA viruses by up to 2 days. Our work shows the ease with which cyanobacteria from culture collections can be screened for antiviral activity and highlights the potential of cyanobacteria as an excellent source for the discovery of novel antiviral compounds, warranting further investigation.

Ultraviolet germicidal irradiation: A prediction model to estimate UV-C-induced infectivity loss in single-strand RNA viruses.

Ultraviolet germicidal irradiation (UVGI) disinfection technology is effective in inactivating microorganisms. However, its performance can vary against different microorganisms due to their diverse structural and genomic features. Thus, rapid predictions of UV (254 nm) inactivation kinetics are essential, particularly for highly infectious emerging pathogens, such as SARS-CoV-2, during the extemporary COVID-19 pandemic. In this study, aiming at single-strand RNA (ssRNA) viruses, an improved genomic model was introduced to predict the UV inactivation kinetics of viral genomes using genome sequence data. First, the overall virus infectivity loss in an aqueous matrix was estimated as the sum of damage to both the entire genome and the protein capsid. Then, the "UV rate constant ratio of aerosol and liquid" was used to convert the UV rate constant for viruses in a liquid-based matrix to an airborne state. The prediction model underwent both quantitative and qualitative validation using experimental data from this study and the literature. Finally, with the goal of mitigating potential airborne transmission of ssRNA viruses in indoor environments, this paper summarizes existing in-duct UVGI system designs and evaluates their germicidal performance. The prediction model may serve as a preliminary tool to assess the effectiveness of a UVGI system for emerging or unculturable viruses or to estimate the required UV dose when designing such a system.

The minus strand of positive-sense RNA viruses encodes small proteins.

It is widely accepted that the minus strands of positive single-strand RNA (+ssRNA) viruses function as replication templates only. Gong et al. revealed that the minus strand of two unrelated +ssRNA viruses encodes proteins. This textbook-changing discovery calls for the reconsideration of the molecular mechanisms underlying the infection cycle of +ssRNA viruses.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Roots of Rice Seedlings under Nitrogen Deficiency.

Long non-coding RNAs (lncRNAs) regulate gene expression in eukaryotic organisms. Research suggests that lncRNAs may be involved in the regulation of nitrogen use efficiency in plants. In this study, we identified 1628 lncRNAs based on the transcriptomic sequencing of rice roots under low-nitrogen (LN) treatment through the implementation of an integrated bioinformatics pipeline. After 4 h of LN treatment, 50 lncRNAs and 373 mRNAs were significantly upregulated, and 17 lncRNAs and 578 mRNAs were significantly downregulated. After 48 h LN treatment, 43 lncRNAs and 536 mRNAs were significantly upregulated, and 42 lncRNAs and 947 mRNAs were significantly downregulated. Moreover, the interaction network among the identified lncRNAs and mRNAs was investigated and one of the LN-induced lncRNAs (lncRNA24320.6) was further characterized. lncRNA24320.6 was demonstrated to positively regulate the expression of a flavonoid 3'-hydroxylase 5 gene (OsF3'H5). The overexpression of lncRNA24320.6 was shown to improve nitrogen absorption and promote growth in rice seedlings under LN conditions. Our results provide valuable insights into the roles of lncRNAs in the rice response to nitrogen starvation.

Structural Studies of Bacteriophage Φ6 and Its Transformations during Its Life Cycle.

From the first isolation of the cystovirus bacteriophage Φ6 from Pseudomonas syringae 50 years ago, we have progressed to a better understanding of the structure and transformations of many parts of the virion. The three-layered virion, encapsulating the tripartite double-stranded RNA (dsRNA) genome, breaches the cell envelope upon infection, generates its own transcripts, and coopts the bacterial machinery to produce its proteins. The generation of a new virion starts with a procapsid with a contracted shape, followed by the packaging of single-stranded RNA segments with concurrent expansion of the capsid, and finally replication to reconstitute the dsRNA genome. The outer two layers are then added, and the fully formed virion released by cell lysis. Most of the procapsid structure, composed of the proteins P1, P2, P4, and P7 is now known, as well as its transformations to the mature, packaged nucleocapsid. The outer two layers are less well-studied. One additional study investigated the binding of the host protein YajQ to the infecting nucleocapsid, where it enhances the transcription of the large RNA segment that codes for the capsid proteins. Finally, I relate the structural aspects of bacteriophage Φ6 to those of other dsRNA viruses, noting the similarities and differences.

Immune responses elicited by ssRNA(-) oncolytic viruses in the host and in the tumor microenvironment.

Oncolytic viruses (OVs) are at the forefront of biologicals for cancer treatment. They represent a diverse landscape of naturally occurring viral strains and genetically modified viruses that, either as single agents or as part of combination therapies, are being evaluated in preclinical and clinical settings. As the field gains momentum, the research on OVs has been shifting efforts to expand our understanding of the complex interplay between the virus, the tumor and the immune system, with the aim of rationally designing more efficient therapeutic interventions. Nowadays, the potential of an OV platform is no longer defined exclusively by the targeted replication and cancer cell killing capacities of the virus, but by its contribution as an immunostimulator, triggering the transformation of the immunosuppressive tumor microenvironment (TME) into a place where innate and adaptive immunity players can efficiently engage and lead the development of tumor-specific long-term memory responses. Here we review the immune mechanisms and host responses induced by ssRNA(-) (negative-sense single-stranded RNA) viruses as OV platforms. We focus on two ssRNA(-) OV candidates: Newcastle disease virus (NDV), an avian paramyxovirus with one of the longest histories of utilization as an OV, and influenza A (IAV) virus, a well-characterized human pathogen with extraordinary immunostimulatory capacities that is steadily advancing as an OV candidate through the development of recombinant IAV attenuated platforms.

Recent Advances in Structural Studies of Single-Stranded RNA Bacteriophages.

Positive-sense single-stranded RNA (ssRNA) bacteriophages (phages) were first isolated six decades ago. Since then, extensive research has been conducted on these ssRNA phages, particularly those infecting E. coli. With small genomes of typically 3-4 kb that usually encode four essential proteins, ssRNA phages employ a straightforward infectious cycle involving host adsorption, genome entry, genome replication, phage assembly, and host lysis. Recent advancements in metagenomics and transcriptomics have led to the identification of ~65,000 sequences from ssRNA phages, expanding our understanding of their prevalence and potential hosts. This review article illuminates significant investigations into ssRNA phages, with a focal point on their structural aspects, providing insights into the various stages of their infectious cycle.

Sequencing and Partial Molecular Characterization of BAB-TMP, the Babeș Strain of the Fixed Rabies Virus Adapted for Multiplication in Cell Lines.

The rabies virus is a major zoonosis that causes severe nervous disease in humans, leading to paralysis and death. The world's second anti-rabies center was established in 1888 by Victor Babeș, in Bucharest, where an eponymous strain of rabies was isolated and used to develop a method for immunization. The Babeș strain of the rabies virus was used for over 100 years in Romania to produce a rabies vaccine for human use, based on animal nerve tissue, thus having a proven history of prophylactic use. The present study aimed to sequence the whole genome of the Babeș strain and to explore its genetic relationships with other vaccine strains as well as to characterize its relevant molecular traits. After being adapted for multiplication in cell lines and designated BAB-TMP, 99% of the viral genome was sequenced. The overall organization of the genome is similar to that of other rabies vaccine strains. Phylogenetic analysis indicated that the BAB-TMP strain is closely related to the Russian RV-97 vaccine strain, and both seem to have a common ancestor. The nucleoprotein gene of the investigated genome was the most conserved, and the glycoprotein showed several unique amino acid substitutions within the major antigenic sites and linear epitopes.

Molecular Population Genetics of Aspen Mosaic-Associated Virus in Finland and Sweden.

Aspen mosaic-associated virus (AsMaV) is a newly identified Emaravirus, in the family Fimoviridae, Bunyavirales, associated with mosaic symptoms in aspen trees (Populus tremula). Aspen trees are widely distributed in Europe and understanding the population structure of AsMaV may aid in the development of better management strategies. The virus genome consists of five negative-sense single-stranded RNA (-ssRNA) molecules. To investigate the genetic diversity and population parameters of AsMaV, different regions of the genome were amplified and analyzed and full-length sequence of the divergent isolates were cloned and sequenced. The results show that RNA3 or nucleoprotein is a good representative for studying genetic diversity in AsMaV. Developed RT-PCR-RFLP was able to identify areas with a higher number of haplotypes and could be applied for screening the large number of samples. In general, AsMaV has a conserved genome and based on the phylogenetic studies, geographical structuring was observed in AsMaV isolates from Sweden and Finland, which could be attributed to founder effects. The genome of AsMaV is under purifying selection but not distributed uniformly on genomic RNAs. Distant AsMaV isolates displayed amino acid sequence variations compared to other isolates, and bioinformatic analysis predicted potential post-translational modification sites in some viral proteins.

Wildlife nidoviruses: biology, epidemiology, and disease associations of selected nidoviruses of mammals and reptiles.

Wildlife is the source of many emerging infectious diseases. Several viruses from the order Nidovirales have recently emerged in wildlife, sometimes with severe consequences for endangered species. The order Nidovirales is currently classified into eight suborders, three of which contain viruses of vertebrates. Vertebrate coronaviruses (suborder Cornidovirineae) have been extensively studied, yet the other major suborders have received less attention. The aim of this minireview was to summarize the key findings from the published literature on nidoviruses of vertebrate wildlife from two suborders: Arnidovirineae and Tornidovirineae. These viruses were identified either during investigations of disease outbreaks or through molecular surveys of wildlife viromes, and include pathogens of reptiles and mammals. The available data on key biological features, disease associations, and pathology are presented, in addition to data on the frequency of infections among various host populations, and putative routes of transmission. While nidoviruses discussed here appear to have a restricted in vivo host range, little is known about their natural life cycle. Observational field-based studies outside of the mortality events are needed to facilitate an understanding of the virus-host-environment interactions that lead to the outbreaks. Laboratory-based studies are needed to understand the pathogenesis of diseases caused by novel nidoviruses and their evolutionary histories. Barriers preventing research progress include limited funding and the unavailability of virus- and host-specific reagents. To reduce mortalities in wildlife and further population declines, proactive development of expertise, technologies, and networks should be developed. These steps would enable effective management of future outbreaks and support wildlife conservation.

Inhibiting immunoregulatory amidase NAAA blocks ZIKV maturation in Human Neural Stem Cells.

Recent evidence suggests that lipids play a crucial role in viral infections beyond their traditional functions of supplying envelope and energy, and creating protected niches for viral replication. In the case of Zika virus (ZIKV), it alters host lipids by enhancing lipogenesis and suppressing ß-oxidation to generate viral factories at the endoplasmic reticulum (ER) interface. This discovery prompted us to hypothesize that interference with lipogenesis could serve as a dual antiviral and anti-inflammatory strategy to combat the replication of positive sense single-stranded RNA (ssRNA+) viruses. To test this hypothesis, we examined the impact of inhibiting N-Acylethanolamine acid amidase (NAAA) on ZIKV-infected human Neural Stem Cells. NAAA is responsible for the hydrolysis of palmitoylethanolamide (PEA) in lysosomes and endolysosomes. Inhibition of NAAA results in PEA accumulation, which activates peroxisome proliferator-activated receptor-α (PPAR-α), directing ß-oxidation and preventing inflammation. Our findings indicate that inhibiting NAAA through gene-editing or drugs moderately reduces ZIKV replication by approximately one log10 in Human Neural Stem Cells, while also releasing immature virions that have lost their infectivity. This inhibition impairs furin-mediated prM cleavage, ultimately blocking ZIKV maturation. In summary, our study highlights NAAA as a host target for ZIKV infection.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex.

BACKGROUND: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. RESULTS: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. CONCLUSIONS: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.