Dataset Information

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus.

ABSTRACT:

Background

Bungarus multicinctus, from which a classical Chinese medicine is produced, is known as the most venomous land snake in the world, but the chromatin organization and transcription factor activity during venom replenishment progress have not been explored yet. This study aimed to determine the roles of chromatin structure in toxin activity via bioinformatics and experimental validation.

Methods

Chromosome conformation capture (Hi-C) analysis was used to examine interactions among chromosomes and identify different scales of chromatin during envenomation in B. multicinctus. Correlations between epigenetic modifications and chromatin structure were verified through ChIP-seq analysis. RNA-seq was used to validate the influence of variation in chromatin structure and gene expression levels on venom production and regulation.

Results

Our results suggested that intra-chromosomal interactions are more intense than inter-chromosomal interactions among the control group, 3-day group of venom glands and muscles. Through this, we found that compartmental transition was correlated with chromatin interactions. Interestingly, the up-regulated genes in more compartmental switch regions reflect the function of toxin activity. Topologically associated domain (TAD) boundaries enriched with histone modifications are associated with different distributions of genes and the expression levels. Toxin-coding genes in the same loop are highly expressed, implying that the importance of epigenetic regulation during envenomination. On a smaller scale, the epigenetic markers affect transcriptional regulation by controlling the recruitment/inhibition of transcription initiation complexes.

Conclusions

Chromatin structure and epigenetic modifications could play a vital status role in the mechanisms of venom regulation in B. multicinctus.

SUBMITTER: Liao X

PROVIDER: S-EPMC8447776 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundKraits (genus Bungarus) and cobras (genus Naja) are two representative toxic genera of elapids in the old world. Although they are closely related genera and both of their venoms are very toxic, the compositions of their venoms are very different. To unveil their detailed venoms and their evolutionary patterns, we constructed venom gland cDNA libraries and genomic bacterial artificial chromosome (BAC) libraries for Bungarus multicinctus and Naja atra, respectively. We sequenced about 1500 cDNA clones for each of the venom cDNA libraries and screened BAC libraries of the two snakes by blot analysis using four kinds of toxin probes; i.e., three-finger toxin (3FTx), phospholipase A2 (PLA2), kunitz-type protease inhibitor (Kunitz), and natriuretic peptide (NP).ResultsIn total, 1092 valid expressed sequences tags (ESTs) for B. multicinctus and 1166 ESTs for N. atra were generated. About 70% of these ESTs can be annotated as snake toxin transcripts. 3FTx (64.5%) and β bungarotoxin (25.1%) comprise the main toxin classes in B. multicinctus, while 3FTx (95.8%) is the dominant toxin in N. atra. We also observed several less abundant venom families in B. multicinctus and N. atra, such as PLA2, C-type lectins, and Kunitz. Peculiarly a cluster of NP precursors with tandem NPs was detected in B. multicinctus. A total of 71 positive toxin BAC clones in B. multicinctus and N. atra were identified using four kinds of toxin probes (3FTx, PLA2, Kunitz, and NP), among which 39 3FTx-positive BACs were sequenced to reveal gene structures of 3FTx toxin genes.ConclusionsBased on the toxin ESTs and 3FTx gene sequences, the major components of B. multicinctus venom transcriptome are neurotoxins, including long chain alpha neurotoxins (α-ntx) and the recently originated β bungarotoxin, whereas the N. atra venom transcriptome mainly contains 3FTxs with cytotoxicity and neurotoxicity (short chain α-ntx). The data also revealed that tandem duplications contributed the most to the expansion of toxin multigene families. Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (dN/dS) indicates that not only multigene toxin families but also other less abundant toxins might have been under rapid diversifying evolution.