TGF-ß2 Induces Ribosome Activity, Alters Ribosome Composition and Inhibits IRES-Mediated Translation in Chondrocytes.

Alterations in cell fate are often attributed to (epigenetic) regulation of gene expression. An emerging paradigm focuses on specialized ribosomes within a cell. However, little evidence exists for the dynamic regulation of ribosome composition and function. Here, we stimulated a chondrocytic cell line with transforming growth factor beta (TGF-ß2) and mapped changes in ribosome function, composition and ribosomal RNA (rRNA) epitranscriptomics. 35S Met/Cys incorporation was used to evaluate ribosome activity. Dual luciferase reporter assays were used to assess ribosomal modus. Ribosomal RNA expression and processing were determined by RT-qPCR, while RiboMethSeq and HydraPsiSeq were used to determine rRNA modification profiles. Label-free protein quantification of total cell lysates, isolated ribosomes and secreted proteins was done by LC-MS/MS. A three-day TGF-ß2 stimulation induced total protein synthesis in SW1353 chondrocytic cells and human articular chondrocytes. Specifically, TGF-ß2 induced cap-mediated protein synthesis, while IRES-mediated translation was not (P53 IRES) or little affected (CrPv IGR and HCV IRES). Three rRNA post-transcriptional modifications (PTMs) were affected by TGF-ß2 stimulation (18S-Gm1447 downregulated, 18S-ψ1177 and 28S-ψ4598 upregulated). Proteomic analysis of isolated ribosomes revealed increased interaction with eIF2 and tRNA ligases and decreased association of eIF4A3 and heterogeneous nuclear ribonucleoprotein (HNRNP)s. In addition, thirteen core ribosomal proteins were more present in ribosomes from TGF-ß2 stimulated cells, albeit with a modest fold change. A prolonged stimulation of chondrocytic cells with TGF-ß2 induced ribosome activity and changed the mode of translation. These functional changes could be coupled to alterations in accessory proteins in the ribosomal proteome.

Loss of Conserved rRNA Modifications in the Peptidyl Transferase Center Leads to Diminished Protein Synthesis and Cell Growth in Budding Yeast.

Ribosomal RNAs (rRNAs) are extensively modified during the transcription and subsequent maturation. Three types of modifications, 2'-O-methylation of ribose moiety, pseudouridylation, and base modifications, are introduced either by a snoRNA-driven mechanism or by stand-alone enzymes. Modified nucleotides are clustered at the functionally important sites, including peptidyl transferase center (PTC). Therefore, it has been hypothesised that the modified nucleotides play an important role in ensuring the functionality of the ribosome. In this study, we demonstrate that seven 25S rRNA modifications, including four evolutionarily conserved modifications, in the proximity of PTC can be simultaneously depleted without loss of cell viability. Yeast mutants lacking three snoRNA genes (snR34, snR52, and snR65) and/or expressing enzymatically inactive variants of spb1(D52A/E679K) and nop2(C424A/C478A) were constructed. The results show that rRNA modifications in PTC contribute collectively to efficient translation in eukaryotic cells. The deficiency of seven modified nucleotides in 25S rRNA resulted in reduced cell growth, cold sensitivity, decreased translation levels, and hyperaccurate translation, as indicated by the reduced missense and nonsense suppression. The modification m5C2870 is crucial in the absence of the other six modified nucleotides. Thus, the pattern of rRNA-modified nucleotides around the PTC is essential for optimal ribosomal translational activity and translational fidelity.

Assessment of different enrichment methods revealed the optimal approach to identify bovine circRnas.

Although circular RNAs (circRNAs) play important roles in regulating gene expression, the understanding of circRNAs in livestock animals is scarce due to the significant challenge to characterize them from a biological sample. In this study, we assessed the outcomes of bovine circRNA identification using six enrichment approaches with the combination of ribosomal RNAs removal (Ribo); linear RNAs degradation (R); linear RNAs and RNAs with structured 3' ends degradation (RTP); ribosomal RNAs coupled with linear RNAs elimination (Ribo-R); ribosomal RNA, linear RNAs and RNAs with poly (A) tailing elimination (Ribo-RP); and ribosomal RNA, linear RNAs and RNAs with structured 3' ends elimination (Ribo-RTP), respectively. RNA-sequencing analysis revealed that different approaches led to varied ratio of uniquely mapped reads, false-positive rate of identifying circRNAs, and the number of circRNAs per million clean reads (Padj <0.05). Out of 2,285 and 2,939 highly confident circRNAs identified in liver and rumen tissues, respectively, 308 and 260 were commonly identified from five methods, with Ribo-RTP method identified the highest number of circRNAs. Besides, 507 of 4,051 identified bovine highly confident circRNAs had shared splicing sites with human circRNAs. The findings from this work provide optimized methods to identify bovine circRNAs from cattle tissues for downstream research of their biological roles in cattle.

Decoding the ribosome's hidden language: rRNA modifications as key players in cancer dynamics and targeted therapies.

Ribosomal RNA (rRNA) modifications, essential components of ribosome structure and function, significantly impact cellular proteomics and cancer biology. These chemical modifications transcend structural roles, critically shaping ribosome functionality and influencing cellular protein profiles. In this review, the mechanisms by which rRNA modifications regulate both rRNA functions and broader cellular physiological processes are critically discussed. Importantly, by altering the translational output, rRNA modifications can shift the cellular equilibrium towards oncogenesis, thus playing a key role in cancer development and progression. Moreover, a special focus is placed on the functions of mitochondrial rRNA modifications and their aberrant expression in cancer, an area with profound implications yet largely uncharted. Dysregulation in these modifications can lead to metabolic dysfunction and apoptosis resistance, hallmark traits of cancer cells. Furthermore, the current challenges and future perspectives in targeting rRNA modifications are highlighted as a therapeutic approach for cancer treatment. In conclusion, rRNA modifications represent a frontier in cancer research, offering novel insights and therapeutic possibilities. Understanding and harnessing these modifications can pave the way for breakthroughs in cancer treatment, potentially transforming the approach to combating this complex disease.

The molecular machinery for maturation of primary mtDNA transcripts.

Human mitochondria harbour a circular, polyploid genome (mtDNA) encoding 11 messenger RNAs (mRNAs), two ribosomal RNAs (rRNAs) and 22 transfer RNAs (tRNAs). Mitochondrial transcription produces long, polycistronic transcripts that span almost the entire length of the genome, and hence contain all three types of RNAs. The primary transcripts then undergo a number of processing and maturation steps, which constitute key regulatory points of mitochondrial gene expression. The first step of mitochondrial RNA processing consists of the separation of primary transcripts into individual, functional RNA molecules and can occur by two distinct pathways. Both are carried out by dedicated molecular machineries that substantially differ from RNA processing enzymes found elsewhere. As a result, the underlying molecular mechanisms remain poorly understood. Over the last years, genetic, biochemical and structural studies have identified key players involved in both RNA processing pathways and provided the first insights into the underlying mechanisms. Here, we review our current understanding of RNA processing in mammalian mitochondria and provide an outlook on open questions in the field.

Directional selection, not the direction of selection, affects telomere length and copy number at ribosomal RNA loci.

Many fisheries exert directional selection on traits such as body size and growth rate. Whether directional selection impacts regions of the genome associated with traits related to growth is unknown. To address this issue, we characterised copy number variation in three regions of the genome associated with cell division, (1) telomeric DNA, (2) loci transcribed as ribosomal RNA (rDNA), and (3) mitochondrial DNA (mtDNA), in three selection lines of zebrafish reared at three temperatures (22 °C, 28 °C, and 34 °C). Selection lines differed in (1) the direction of selection (two lines experienced directional selection for large or small body size) and (2) whether they experienced any directional selection itself. Lines that had experienced directional selection were smaller, had lower growth rate, shorter telomeres, and lower rDNA copy number than the line that experiencing no directional selection. Neither telomere length nor rDNA copy number were affected by temperature. In contrast, mtDNA content increased at elevated temperature but did not differ among selection lines. Though directional selection impacts rDNA and telomere length, direction of such selection did not matter, whereas mtDNA acts as a stress marker for temperature. Future work should examine the consequences of these genomic changes in natural fish stocks.

Wickerhamiella lachancei f.a. sp. nov., a novel ascomycetous yeast species isolated from flowers of Lantana camara in India.

Two isolates representing a novel species of the genus Wickerhamiella were obtained in India from nectar of flowers of Lantana camara, an ornamental exotic species native to Central and South America. Phylogenetic analyses of the D1/D2 domain of the 26S large subunit (LSU) rRNA gene, internal transcribed spacer (ITS) region, and physiological characteristics, supported the recognition of the novel species, that we designate Wickerhamiella lachancei sp. nov (MycoBank no. MB851709), with MCC 9929T as the holotype and PYCC 10003T as the isotype. Considering pairwise sequence similarity, the type strain of the novel species differs from the type strain of the most closely related species, Wickerhamiella drosophilae CBS 8459T, by 16 nucleotide substitutions and two gaps (3.9 % sequence variation) in the D1/D2 region (560 bp compared) and 28 nucleotide substitutions and five gaps (7.22 % sequence variation) in the ITS region (444 bp compared).

Role of Rmd9p in 3'-end processing of mitochondrial 15S rRNA in Saccharomyces cerevisiae.

Ribosome biogenesis, involving processing/assembly of rRNAs and r-proteins is a vital process. In Saccharomyces cerevisiae mitochondria, ribosomal small subunit comprises 15S rRNA (15S). While the 15S 5'-end processing uses Ccm1p and Pet127p, the mechanisms of the 3'-end processing remain unclear. We reveal involvement of Rmd9p in safeguarding/processing 15S 3'-end. Rmd9p deficiency results in a cleavage at a position 183 nucleotides upstream of 15S 3'-end, and in the loss of the 3'-minor domain. Rmd9p binds to the sequences in the 3'-end region of 15S, and a genetic interaction between rmd9 and dss1 indicates that Rmd9p regulates/limits mtEXO activity during the 3'-end spacer processing.

The regulatory roles of small nucleolar RNAs within their host locus.

Small nucleolar RNAs (snoRNAs) are a class of conserved noncoding RNAs forming complexes with proteins to catalyse site-specific modifications on ribosomal RNA. Besides this canonical role, several snoRNAs are now known to regulate diverse levels of gene expression. While these functions are carried out in trans by mature snoRNAs, evidence has also been emerging of regulatory roles of snoRNAs in cis, either within their genomic locus or as longer transcription intermediates during their maturation. Herein, we review recent findings that snoRNAs can interact in cis with their intron to regulate the expression of their host gene. We also explore the ever-growing diversity of longer host-derived snoRNA extensions and their functional impact across the transcriptome. Finally, we discuss the role of snoRNA duplications into forging these new layers of snoRNA-mediated regulation, as well as their involvement in the genomic imprinting of their host locus.

A role for a Trypanosoma brucei cytosine RNA methyltransferase homolog in ribosomal RNA processing.

In Trypanosoma brucei, gene expression is primarily regulated posttranscriptionally making RNA metabolism critical. T. brucei has an epitranscriptome containing modified RNA bases. Yet, the identity of the enzymes catalyzing modified RNA base addition and the functions of the enzymes and modifications remain unclear. Homology searches indicate the presence of numerous T. brucei cytosine RNA methyltransferase homologs. One such homolog, TbNop2 was studied in detail. TbNop2 contains the six highly conserved motifs found in cytosine RNA methyltransferases and is evolutionarily related to the Nop2 protein family required for rRNA modification and processing. RNAi experiments targeting TbNop2 resulted in reduced levels of TbNop2 RNA and protein, and a cessation of parasite growth. Next generation sequencing of bisulfite-treated RNA (BS-seq) detected the presence of two methylation sites in the large rRNA; yet TbNop2 RNAi did not result in a significant reduction of methylation. However, TbNop2 RNAi resulted in the retention of 28S internal transcribed spacer RNAs, indicating a role for TbNop2 in rRNA processing.

The invasive land flatworm Arthurdendyus triangulatus has repeated sequences in the mitogenome, extra-long cox2 gene and paralogous nuclear rRNA clusters.

Using a combination of short- and long-reads sequencing, we were able to sequence the complete mitochondrial genome of the invasive 'New Zealand flatworm' Arthurdendyus triangulatus (Geoplanidae, Rhynchodeminae, Caenoplanini) and its two complete paralogous nuclear rRNA gene clusters. The mitogenome has a total length of 20,309 bp and contains repetitions that includes two types of tandem-repeats that could not be solved by short-reads sequencing. We also sequenced for the first time the mitogenomes of four species of Caenoplana (Caenoplanini). A maximum likelihood phylogeny associated A. triangulatus with the other Caenoplanini but Parakontikia ventrolineata and Australopacifica atrata were rejected from the Caenoplanini and associated instead with the Rhynchodemini, with Platydemus manokwari. It was found that the mitogenomes of all species of the subfamily Rhynchodeminae share several unusual structural features, including a very long cox2 gene. This is the first time that the complete paralogous rRNA clusters, which differ in length, sequence and seemingly number of copies, were obtained for a Geoplanidae.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis.

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

Using digital PCR to predict ciliate abundance from ribosomal RNA gene copy numbers.

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.

The DEAD-box ATPase Dbp10/DDX54 initiates peptidyl transferase center formation during 60S ribosome biogenesis.

DEAD-box ATPases play crucial roles in guiding rRNA restructuring events during the biogenesis of large (60S) ribosomal subunits, but their precise molecular functions are currently unknown. In this study, we present cryo-EM reconstructions of nucleolar pre-60S intermediates that reveal an unexpected, alternate secondary structure within the nascent peptidyl-transferase-center (PTC). Our analysis of three sequential nucleolar pre-60S intermediates reveals that the DEAD-box ATPase Dbp10/DDX54 remodels this alternate base pairing and enables the formation of the rRNA junction that anchors the mature form of the universally conserved PTC A-loop. Post-catalysis, Dbp10 captures rRNA helix H61, initiating the concerted exchange of biogenesis factors during late nucleolar 60S maturation. Our findings show that Dbp10 activity is essential for the formation of the ribosome active site and reveal how this function is integrated with subsequent assembly steps to drive the biogenesis of the large ribosomal subunit.

Initial study and phylogenetic analysis of hard ticks (Acari: Ixodidae) in Nantong, China along the route of avian migration.

The growing concern about migratory birds potentially spreading ticks due to global warming has become a significant issue. The city of Nantong in this study is situated along the East Asia-Australasian Flyway (EAAF), with numerous wetlands serving as roosting sites for migratory birds. We conducted an investigation of hard ticks and determined the phylogenetic characteristics of tick species in this city. We utilized three different genes for our study: the mitochondrial cytochrome oxidase subunit 1 (COX1) gene, the second internal transcribed spacer (ITS2), and the mitochondrial small subunit rRNA (12 S rRNA) gene. The predominant tick species were Haemaphysalis flava (H. flava) and Haemaphysalis longicornis (H. longicornis). Additionally, specimens of Haemaphysalis campanulata (H. campanulata) and Rhipicephalus sanguineus (R. sanguineus) were collected. The H. flava specimens in this study showed a close genetic relationship with those from inland provinces of China, as well as South Korea and Japan. Furthermore, samples of H. longicornis exhibited a close genetic relationship with those from South Korea, Japan, Australia, and the USA, as well as specific provinces in China. Furthermore, R. sanguineus specimens captured in Nantong showed genetic similarities with specimens from Egypt, Nigeria, and Argentina.

Analysis of the Overlength Main Noncoding Region in Metacarcinus magister (Decapoda: Brachyura) and a Phylogenetic Study of the Cancroidea Species.

Complete mitochondrial genomes (mitogenomes) can provide important information regarding the molecular evolution and phylogenetic relationships of marine invertebrates, especially in Brachyura. Only one Cancroidea species of mitogenomes has been sequenced before; in this research, the mitogenomic characteristics of Metacarcinus magister (Cancridae: Cancroidea) are newly studied. The length of the M. magister mitogenome was 48,820 bp, and it contained the typical 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. We performed a series of analyses on the characteristics of the mNCR of M. magister. The phylogenetics, life circumstances, and selective pressures were all analyzed to explain the formation of this length, which revealed the length of the M. magister mitogenome to be approximately three times greater than the normal length of Brachyuran mitogenomes. Phylogenetic analyses based on a dataset of 215 Decapodan mitogenomes indicated that all Eriphioidea crabs were clustered together as a group. Moreover, the rearrangement mechanism of the Cancroidea species was predicted to provide stronger evidence for the phylogenetic analysis. In general, the results obtained in this study will contribute to a better understanding of the cause of the unusual length of the M. magister mitogenome and provide new insights into the phylogeny of Brachyura.

Complete Mitochondrial Genome for Lucilia cuprina dorsalis (Diptera: Calliphoridae) from the Northern Territory, Australia.

The Australian sheep blowfly, Lucilia cuprina dorsalis, is a major sheep ectoparasite causing subcutaneous myiasis (flystrike), which can lead to reduced livestock productivity and, in severe instances, death of the affected animals. It is also a primary colonizer of carrion, an efficient pollinator, and used in maggot debridement therapy and forensic investigations. In this study, we report the complete mitochondrial (mt) genome of L. c. dorsalis from the Northern Territory (NT), Australia, where sheep are prohibited animals, unlike the rest of Australia. The mt genome is 15,943 bp in length, comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and a non-coding control region. The gene order of the current mt genome is consistent with the previously published L. cuprina mt genomes. Nucleotide composition revealed an AT bias, accounting for 77.5% of total mt genome nucleotides. Phylogenetic analyses of 56 species/taxa of dipterans indicated that L. c. dorsalis and L. sericata are the closest among all sibling species of the genus Lucilia, which helps to explain species evolution within the family Luciliinae. This study provides the first complete mt genome sequence for L. c. dorsalis derived from the NT, Australia to facilitate species identification and the examination of the evolutionary history of these blowflies.

A scalable and cost-efficient rRNA depletion approach to enrich RNAs for molecular biology investigations.

Transcriptomics analyses play pivotal roles in understanding the complex regulatory networks that govern cellular processes. The abundance of rRNAs, which account for 80%-90% of total RNA in eukaryotes, limits the detection and investigation of other transcripts. While mRNAs and long noncoding RNAs have poly(A) tails that are often used for positive selection, investigations of poly(A)- RNAs, such as circular RNAs, histone mRNAs, and small RNAs, typically require the removal of the abundant rRNAs for enrichment. Current approaches to deplete rRNAs for downstream molecular biology investigations are hampered by restrictive RNA input masses and high costs. To address these challenges, we developed rRNA Removal by RNaseH (rRRR), a method to efficiently deplete rRNAs from a wide range of human, mouse, and rat RNA inputs and of varying qualities at a cost 10- to 20-fold cheaper than other approaches. We used probe-based hybridization and enzymatic digestion to selectively target and remove rRNA molecules while preserving the integrity of non-rRNA transcripts. Comparison of rRRR to two commercially available approaches showed similar rRNA depletion efficiencies and comparable off-target effects. Our developed method provides researchers with a valuable tool for investigating gene expression and regulatory mechanisms across a wide range of biological systems at an affordable price that increases the accessibility for researchers to enter the field, ultimately advancing our understanding of cellular processes.

Oxytoxaceae are prorocentralean rather than peridinialean dinophytes and taxonomic clarification of heterotrophic Oxytoxum lohmannii (≡ "Amphidinium" crassum) by epitypification.

During evolution of Dinophyceae, size reduction of the episome has occurred in several lineages (including unarmoured Amphidiniales and armoured Prorocentrales). One such species is Amphidinium crassum, whose taxonomic identity is elusive though showing morphological similarities with Oxytoxaceae (currently placed in armoured Peridiniales). Plankton samples were taken at the type locality of A. crassum in Kiel Bight (Baltic Sea) in order to establish monoclonal strains. The protist material was examined in detail using light and electron microscopy, and a long (2984 bp) ribosomal RNA sequence gained was part of a taxon sample comprising 206 specimen vouchers and representing the known molecular diversity of Dinophyceae. Cells of A. crassum were ovoid and exhibited a plate pattern po, 4', 1a, 6'', 5c, 4s, 5''', 1''''. In the molecular phylogeny, the species seemed to belong neither to Amphidiniales nor to Peridiniales but to Prorocentrales and clustered with other representatives of Oxytoxaceae. The morphological diversity of Prorocentrales appears thus expanded, and the group may include a number of previously unrecognised representatives unusually having five postcingular and only a single antapical plate. The taxonomic identity of A. crassum is clarified by epitypification, and the species notably exhibits both an apical pore and an additional epithecal pore.

The complete mitochondrial genome of Chibiraga houshuaii (Lepidoptera, Limacodidae) and its phylogenetic implications.

Chibiraga is a mall East Asian genus in the family Limacodidae (slug-moths). The latter includes many agricultural pests. Mitochondrial genome analysis is an important tool for studying insect molecular identification and phylogenetics. However, there are very few mitogenome sequences available for Limacodidae species, and none for the genus Chibiraga at all. To explore the mitogenome features of Chibiraga and verify its phylogenetic position, the complete mitogenome of Chibiraga houshuaii was sequenced and annotated. The complete 15,487 bp genome encoded 37 mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region (CR). Most of the PCGs had typical ATN start codons and terminated with TAA or a single T residue. UUA (Leu2), AUU (Ile), UUU (Phe), AUA (Met) and AAU (Asn) were the five most frequently used codons. All tRNAs were folded into cloverleaf secondary structure, except for trnS1, which lacked the DHU arm. Phylogenetic analyses within the superfamily Zygaenoidea were performed based on multiple datasets from mitochondrial genes. The results showed that the families Phaudidae, Limacodidae and Zygaenidae were respectively recovered as monophyly; C. houshuaii was clustered in a clade with nettle type larvae in Limacodidae.