RNA profile diversity across arthropoda: guidelines, methodological artifacts, and expected outcomes

Abstract High-quality RNA is an important precursor for high-throughput RNA sequencing (RNAseq) and subsequent analyses. However, the primary metric used to assess RNA quality, the RNA Integrity Number (RIN), was developed based on model bacterial and vertebrate organisms. Though the phenomenon is not widely recognized, invertebrate 28S ribosomal RNA (rRNA) is highly prone to a form of denaturation known as gap deletion, in which the subunit collapses into two smaller fragments. In many nonmodel invertebrates, this collapse of the 28S subunit appears as a single band similar in size to the 18S rRNA subunit. This phenomenon is hypothesized to be commonplace among arthropods and is often misinterpreted as a “degraded” rRNA profile. The limited characterization of gap deletion in arthropods, a highly diverse group, as well as other nonmodel invertebrates, often biases RNA quality assessments. To test whether the collapse of 28S is a general pattern or a methodological artifact, we sampled more than half of the major lineages within Arthropoda. We found that the 28S collapse is present in ∼90% of the species sampled. Nevertheless, RNA profiles exhibit considerable diversity with a range of banding patterns. High-throughput RNAseq and subsequent assembly of high-quality transcriptomes from select arthropod species exhibiting collapsed 28S subunits further illustrates the limitations of current RIN proxies in accurately characterizing RNA quality in nonmodel organisms. Furthermore, we show that this form of 28S denaturation, which is often mistaken for true “degradation,” can occur at relatively low temperatures.


Introduction
Next-generation RNA sequencing (RNAseq) has become a revolutionary tool for RNA-based characterization studies. Similar to traditional quantification methods (i.e., quantitative polymerase chain reaction (PCR)), it can now be performed in most modern laboratories with relative ease and accuracy [1] while yielding more precise and robust datasets [2]. For this reason, RNAseq has become a powerful technique for studying an array of biological questions in both model and nonmodel organisms, including de novo transcriptome reconstruction and profiling, gene discovery and functional analyses, differential gene expression studies [3], and evolutionary investigations (e.g., [4]). Moreover, continued advancements to sequencing technologies and analysis software have expanded the reach and utility of RNAseq to include microscopic inputs (singe-cell RNAseq), enabling spatial mapping and tissue reconstruction [5,6].
Compared with DNA, total RNA extraction and sequencing poses nontrivial methodological challenges to maintain RNA stability. Due to its ephemeral nature, RNA is prone to cellular degradation [7] from the initial cellular lysis step to final visualization and quantitation. In addition, ribonucleases (RNases), ubiquitous enzymes present in cells and the environment can further degrade RNA if proper care is not taken during sample preservation, extraction, and downstream protocols [7]. To prevent degradation and loss of transcriptomic information, most extraction guidelines suggest using dedicated equipment in isolated areas and frequent, thorough RNase decontamination throughout the protocol. Ideally, samples should be flash frozen or preserved in specialized salt-saturated buffers (e.g., RNAlater V R ; [8]), when the use of fresh tissue is not feasible, and kept frozen at À80 C prior to extraction and post-isolation to prevent native degradation. These rigorous methodological demands portray the labile nature of RNA.
Historically, RNA was assumed to be of high quality if the ratio of large to small ribosomal subunits (28S:18S) was $2 or higher [9], though this method of assessment is sample consuming and imprecise [10]. Presently, there are two common methods used to test for RNA integrity and degradation after extraction: (i) visualization via denaturing (agarose) gel electrophoresis and (ii) fragment size chip electropherograms on a Bioanalyzer (Agilent). In particular, the Bioanalyzer step calculates the RNA integrity number (RIN) [11], now considered as the gold standard for successful library preparation and sequencing [10,[12][13][14][15]. The RIN model is estimated from RNA fragment distribution peaks and areas under the curve, and ranges from 1 to 10 (low to pristine quality, respectively). Sample RIN values of !7 are typically deemed as "high-quality" and appropriate for RNAseq (e.g., [16]), and often used as a cutoff for which sequencing facilities will guarantee results. Another common standard used to assess RNA integrity is ribosomal subunit size. Using rodent and human models, size has been estimated to be between 4400 and 5000 base pairs (bps) in the large 28S subunit and $2000 bp in the small 18S subunit [17]. Since RNA degradation can occur at any step in most protocols, it is paramount to determine appropriate reference standards and metrics for accurate RNA quality control [12].
Although relevant ribosomal RNA (rRNA) profile standards exist for several model organisms [18], most branches across the "Tree of Life" lack genomic references and resources. This paucity of genomic data is particularly visible across Arthropoda, a globally distributed group that encompasses >80% of animal diversity [19] and occupies a variety of niches from the poles to the tropics, high mountains, and deep-sea trenches. They fulfill significant roles in ecosystems (e.g., pollination and pests) and are of medical (i.e., disease vectors) and economic importance (e.g., fisheries) [20]. Due to their global relevance, the need for increased genomic references and reliable metrics to evaluate nucleic acid quality is crucial for the understanding of multiple evolutionary and ecological phenomena.
There are few recent studies that have evaluated the universality and suitability of "standard" RNA quality metrics for nonmodel organisms (e.g., [21]). Early observations of RNA profiles in insects documented a "hidden break" of the 28S subunit, also referred to as gap deletion, in which the 28S denatures and collapses (decouples) into two smaller fragments (a and b) [22,23] and a short, variable stretch of rRNA is removed or deleted [24,25). This produces an apparent "degraded" RNA profile where both 18S and 28S overlap in a single band. As RNAseq has become a common method for RNA investigations in recent years, more and more invertebrate rRNA profiles suggest that this "hidden break" -a form of 28S denaturation -is the norm rather than the exception [21,26,27].
Hidden breaks in RNA are correlated with the presence of a cleavage site characterized as an Uracil (U) and Adenine (A) or UAAU-rich segment in the rRNA loop [28]. High temperature exposure is believed to result in UAAU decoupling and the subsequent dissociation (collapse) of the 28S [21][22][23]. This dissociation was hypothesized to occur as Hydrogen-bonds in this AU-rich region melt -the strength of the bonds are comparable to the strength of the chemical changes associated with denaturing conditions [22]; the double bonds of the AU pairs are weaker relative to the triple bonds of Guanine-Cytosine (GC) associations [29]. Differences in base-stacking stability [30], secondary structures and enzyme-dependent cleavage events are also thought to play a role [25]. Likewise, mutations in this UAAU-rich segment are thought to prevent this dissociation [28]. In a recent study, McCarthy et al. [21] hypothesized that exposing RNA to heat at 70 C, mainly for electropherogram preparation in a Bioanalyzer, was responsible for the 28S subunit collapse. However, there are additional steps often used in RNA extraction protocols, such as gel electrophoresis and deoxyribonuclease (DNase) applications for DNA removal, which may also expose samples to temperatures high enough to trigger 28S cleavage. As DNase treatments can be performed at relatively low temperatures ($37 C) for enzymes that do not require deactivation at 70 C, and temperature can be regulated during electrophoresis (e.g., ice bath and chilled buffers), 70 C may be an exaggerated upper limit for this denaturation to occur. Given the UAAU loop region is sensitive to heat, DNase enzyme activation temperatures of 37 C may result in 28S denaturation prior to 70 C exposure, and possibly result in different "degraded" profiles. Moreover, wider taxon sampling may reveal variable processing of this gap (loop) region, including insertions/deletions of additional UAAU-rich segments, yielding additional RNA profiles.
Here we document the range of variation among RNA extraction profiles for the major lineages within Arthropoda, with dense taxonomic sampling. Our primary motivation is to disentangle natural patterns of subunit denaturing from true degradation and methodological artifacts. Data available from prior studies are a composite of multiple extraction protocols from different laboratories that can further obscure these artifacts. This study was done using standardized reagents and conditions to improve comparisons among taxa. High-throughput RNAseq and de novo transcriptome assembly for two representative arthropod species exhibiting seemingly denatured RNA profiles provide further support for high-quality RNA despite evidence of 28S subunit collapse. Findings from this effort will aid future nucleic acid quality assessments for arthropods and other nonmodel invertebrates and will provide methodological insight for laboratories beginning to conduct RNA-based research.

Specimen preservation and tissue harvesting
We collected live specimens from the Biscayne Bay Campus of Florida International University (FIU) in North Miami, FL, USA, and brought them back to the lab for preservation. Some specimens were ordered from a commercial vendor (Carolina Biology, NC, USA) or donated by colleagues (see Acknowledgements section). Smaller individuals (bee size and below) were flash frozen and preserved at À80 C. Larger individuals were placed in RNAlater (Life Technologies, Foster City, CA, USA) with an abdominal incision and then placed at À80 C until RNA extraction. All specimens/or tissues are cataloged and deposited as vouchers (specimen or tissue) in the Florida International Crustacean Collection at FIU (Table 1).

RNA extraction
We extracted RNA from 40 species across the four Arthropod subphyla. RNA was extracted from replicates where possible, including: the isopods Asellus aquaticus (n ¼ 3), Cubaris cf. murinus (n ¼ 3), amphipod Niphargus hrabei (n ¼ 3), cave shrimp Barbouria cubensis (n ¼ 3), silverfish Lepisma saccharina (n ¼ 2), the insects Brachymesia gravida (n ¼ 2) and Panorpa debilis (n ¼ 2), the millipede Anadenobolus monilicornis (n ¼ 4), and centipede Hemiscolopendra marginata (n ¼ 4). Though RNA quality tissue replicates were difficult to obtain for some species, all species were included in the study regardless as an effort to more thoroughly survey RNA profile diversity across Arthropoda. All stages of our RNA extraction protocol were carried out under a fume hood in a dedicated RNase-free room. Bleach (2% solution in distilled water) was applied to all surfaces, equipment, and reagent containers to remove RNases and secondarily sterilized with RNaseZAP V R (Life Technologies, Carlsbad, CA, USA). We isolated total RNA from varied tissues or body parts depending on the specimen size (see Table 1). For very small specimens, we used the entire body or pooled multiple individuals. Tissues were homogenized in 1 mL TRIZOL V R Reagent (Life Technologies, Carlsbad, CA, USA). Approximately 10-12 sterile 2 mm zirconia/ ceramic beads were used in a MiniBead Beater (BioSpec, Bartlesville, OK, USA) for homogenization in a 2 mL screwcap tube, followed by the addition of chloroform (0.2 mL). Samples were mixed (shaken) vigorously by hand for 15 s and incubated at room temperature for 2-3 min. Phase separation was carried out by centrifuging samples (12 000g for 15 min) at 4 C, following manufacturer's instructions, and carefully transferring only the top (clear) aqueous phase containing the RNA into new vials. RNA was then precipitated from the clear aqueous phase by mixing in ($100%) isopropyl alcohol (0.5 mL) and glycogen (0.1 lL). Samples were then incubated at room temperature for 10 min and centrifuged (12 000g for 15 min) at 4 C. The supernatant was discarded, the RNA pellet washed twice with freshly made 75% ethanol (EtOH) and the supernatant was again discarded; 1 mL EtOH was added before centrifugation at 7500g for 5 min at 4 C (X2). Before elution, precipitated RNA was treated with DNases (Clonetech, Mountain View, CA, USA) to remove genomic DNA (15-min incubation at 37 C) following manufacturer's guidelines. After repeating the isopropyl precipitation and ethanol washes (2Â), total RNA was eluted in 50 lL of RNase-free water (Sigma, St Louis, MO, USA) and stored at À80 C.
RNA integrity was determined initially on a 1% agarose gel (made with Tris base, acetic acid and Edetic acid (EDTA) or TAE) treated with a 1% bleach solution and ran in cold buffer at 4 C to prevent RNA degradation during electrophoresis (e.g., 500 lL bleach diluted in 50 mL agarose TAE solution [36]). In order to determine the size of the 18S and 28S rRNA subunits and check for degraded samples, we added a RiboRuler RNA ladder (ThermoFisher, USA) to the agarose gel. The RNA concentration of each extraction was measured using a Qubit 2.0 fluorometer (Life Technologies, USA). Finally, we diluted a subsample from every extraction to 5 ng/lL and determined the RNA profile by measuring ribosomal number, size, and RIN number using a BioAnalyzer 2100 with RNA 6000 picochips (Agilent, USA) in the DNA Core Facility at FIU; all samples were processed without the heat exposure (70 C) typical of bioanalyzer protocols. Gloves were changed between each protocol step and RNaseZap was used liberally.

Distinguishing true RNA degradation from analytical artifacts
High temperature exposure is considered as one of the determining factors behind the 28S subunit collapse [21,22]. To test whether denaturation could occur at temperatures <70 C, we chose three representative lineages with differing rRNA profiles: millipede, centipede, and crab. We used a bacterial reference (total) RNA (Escherichia coli, Life Technologies, USA) as a noncleaved 28S positive control. RNA aliquots were taken from the three samples prior to the DNase treatment, precipitated, eluted, and kept at 4 C. We prepared a heat treatment gradient as recommended by the Bioanalyzer manufacturer, heating all samples in a MJ Research PTC-200 thermocycler, incubating for 2 min at a designated temperature using a nonheated lid, prior to cooling at 4 C. RNA from the four samples (including the positive control) were aliquoted into eight different temperature treatment sets (0.2 mL PCR strip tubes). Each set was placed across a temperature gradient ranging from 20 C to 90 C, at increments of 10 C. This range encompasses "room temperature" and standard DNase treatment temperatures, as well as temperatures 20 C higher than the recommended denaturing step in an attempt to force denaturation for comparison. After the eight tube sets were incubated at their corresponding temperature on the thermal gradient, we kept the samples at 4 C and stored them permanently at À80 C until they were run on both a 1% TAE bleach gel in cool buffer and Bioanalyzer picochips as described previously.

An array of RNA profile patterns within arthropoda
In total, we found six distinct patterns of 28S denaturation (Table 1), and banding patterns were consistent across replicates. In most arthropod lineages, we found a consistent "denatured" or "collapsed" RNA profile, where the 28S subunit is denatured into equally sized a-b fragments (Fig. 1). This dominant pattern was found in 22 species and appeared on the gel as a single thick band resulting from the collapse of 28S and subsequent overlap with 18S at around 1900 bp (Fig. 1a-f); a size estimate reported previously for the honey bee Apis mellifera (27). We detected two prominent bands in the fruit fly, cranefly (Diptera), pelagic shrimp, banded-coral shrimp, blue crab, and mudcrab (Decapoda), at $2000 and 1900 bp, respectively. Three well-defined bands were detected in river isopods (Isopoda), marine crabs (Decapoda), and the water flea (Cladocera) at $2200, 2000, and 1900 bp. The hermit crab (Decapoda) appeared to have four bands of similar size between 1800 and 2200 bp although it is possible that the presence of minor degradation may have confounded this profile. The only lineage with four distinct bands was the tiger centipede (Chilopoda) at $2000, 1800, 1000 and 800 bp. Land isopods (Isopoda) showed five distinct bands at $2200, 2000, 1900, 1700 and 600 bp though some degradation was present. Lastly, we found a nondenatured, two-band profile for cave shrimp (Decapoda) and bumblebee millipede (Diplopoda) -with subunit bands at $3900 bp and 1900 bp, similar to the model rRNA references [17] (Fig. 1). It is important to note that the rRNA fragments of some replicate lanes appear to have shifted on the picochip (e.g., Fig. 1f, lane 05) and would require manual adjustments for accurate size estimation (see Table 1), though the banding patterns remain consistent.
Some samples did show "true" degradation <1800 bp including the hermit crab Clibanarius tricolor and one replicate of the dragonfly Brachymesia gravida (Table 1 and Fig. 1), which appeared as multiple bands on the Bioanalyzer traces. Likewise, residual DNA contamination was evident in the mud crab Eurypanopeus depressus after the DNase treatment (Fig. 1e-l); it is possible that the DNase was not as effective because a higher volume of tissue was available and used for this particular extraction. Many samples also exhibited remnants of the original size of the 28S subunit below the 4000 bp mark, appearing as a faint band and a shorter peak on the electropherograms (Fig. 1). Mean sizes for the original 28S and the 18S peaks suggest that, in arthropods, RNA profiles with no denaturation can be measured at $3900 bp and 1900 bp, respectively. Several samples showed multiple faint bands beyond 1900 bp, which may imply some level of degradation and/or remnants of extracted 5/5.8S ribosomal subunits and tRNAs (Fig. 1).

Ribosomal 28S denaturation can occur at "low" temperatures
The temperature gradient revealed that samples with the 28S "hidden break" (blue centipede and blue crab) had at least a portion of their 28S subunit already collapsed at 20 C ( Fig. 2 and Supplementary Fig. S1). Blue crab RNA maintained some copies of its 28S subunit at $4000 bp until reaching 80 C. The positive bacterial control (E. coli) and the millipede (A. monilicornis) maintained 28S subunit integrity throughout the 20-90 C gradient, exhibiting a two-band profile reminiscent of the vertebrate rRNA model reference. The blue centipede (H. marginata) showed a single band throughout the entire gradient ( Fig. 2 and Supplementary Fig. S1). Some lanes on the electropherograms (Fig. 2) showed RNA degradation not found on the agarose gel ( Supplementary Fig. S1).

Discussion
One size does not fit all: diversity of RNA banding profiles across arthropoda There is a growing body of evidence suggesting that the current vertebrate RNA profile references are not reliable proxies to assess the quality of invertebrate RNA with seemingly "degraded" RNA profile patterns [21-23, 27, 28]. Our results suggest that 28S "collapse" or "denaturation" is a common occurrence across Arthropoda, but is likely a complex phenomenon resulting in diverse banding patterns, some of which are newly reported here. This renders model organism RNA profiles of limited use for various nonmodel taxa. Moreover, RINs could not be estimated for a majority of the samples exhibiting the denatured 28S rRNA subunits; paradoxically, RINs were only estimated for a few samples showing some level of "true" degradation [e.g., Fig. 1f  (L01-03)]. Therefore, RIN estimation may be impossible or biased if algorithms do not account for denaturation and collapsed ribosomal subunits in a variety of invertebrates though manual bioanalyzer corrections (i.e., higher anomaly thresholds upon detection of the additional peak(s) or signal) can partially reduce bias [11]. In addition to deviations from the common rRNA twobanded profile, our results further show that rRNA 28S in invertebrates is smaller than that in vertebrates, 3900 bp versus 4400-5000 bp, though the 18S sizes appear to be within a 100 bp differential margin. Bacterial RNA profile sizes are even smaller (E. coli 2900 bp for 23S and 1700 bp for 16S; Fig. 1e). These estimated size differences will allow for a more accurate and prompt assessment of rRNA in nonmodel invertebrates following RNA extraction, including the presence of bacterial and human contamination.
The diversity of RNA profiles, ranging from one to five bands ( Fig. 1 and Table 1), indicates that even though a oneband profile is the most frequently observed pattern among arthropods, individual lineages may vary. Some single bands within Insecta were also relatively thick, suggesting a possible artifact resulting from two high-quality overlapping bands of similar size. Denatured two-banded profiles, evident in several decapod crustaceans and insects [e.g., Fig.1a (L01), not to be mistaken with the banding profile of the vertebrate model reference], suggest that one of the two 28S fragments (a or b) is at least 100 bp longer than 18S and the other is very similar to 18S in length. The three-banded profile found among river isopods,  Table 1 These samples match those used in Supplementary Fig. S1. mud and blue crabs, and Daphnia water fleas, suggests a significant translocation of the UAAU loop with 400 bp separating the two 28S fragments. Likewise, the possibility exists that a second UAAU region, or underlying secondary structure, is creating an additional cleavage point. A four-banded profile in tiger centipedes also implies that complex secondary structures within the loop region are creating additional cleavage points, possibly making them more accessible to enzymes (e.g., RNases). Alternatively, a secondary UAAU-rich loop may be present, resulting in three 28S fragments (a, b, and a novel c) as tiger centipede RNA showed no signs of degradation (i.e., smears or a multitude of fragments). The third and fourth bands in tiger centipedes are also much smaller than bacterial rRNA subunits, rendering contamination from mitochondrial DNA or microorganisms unlikely. However, this four-banded pattern was not found in the blue centipede species (single band). This suggests that the four-banded pattern is not unanimous among Chilopoda. Similar to the RNA profiles observed in aphids [28,42], millipedes (Diplopoda), the sister clade to centipedes, maintained the integrity of the 28S subunit, suggesting this four-banded pattern could be a genus or speciesspecific phenomenon.
The observed variation in RNA profiles and 28S denaturation found in this study is comparable to other systems. For example, a majority of mammals have a typical nondenatured two-banded vertebrate profile, but rodents in the genus Ctenomys have a denatured 28S subunit [17]. Thus, an important distinction exists between subunit collapse/denaturing and degradation. The generation of high-quality de novo transcriptome assemblies from RNAseq data for two arthropod species exhibiting denatured RNA profiles provides additional support for high-quality RNA despite 28S subunit collapse [37]. Assembly statistics for both species compare favorably with previously published RNAseq studies on A. aquaticus [43] and other amphipods [44], as well as arthropod species with nondenatured RNA profiles [45]. Denaturation appears to be an arthropod-wide phenomenon with some lineage-specific exceptions, whereas degradation remains constant due to the ephemeral nature of RNA [46]. For researchers working on organisms with unknown rRNA profiles, we therefore suggest documenting (replicated) rRNA banding patterns to make accurate assessments about RNA integrity and avoid falsely labeling 28S rRNA denaturing as "true" degradation.
It's getting hot in here: the effect of temperature on 28S denaturation A prior study by McCarthy et al. [21] examining rRNA degradation in a few arthropod lineages theorized 28S cleavage would occur at temperatures around 70 C as part of the Bioanalyzer protocol. Our temperature-gradient experiments showed that 28S denaturing can occur across a broad array of temperatures in arthropods, including room temperature ($20 C). This can have important implications given human body temperatures are $37 C, and consequently excessive sample handling could contribute to unpredicted denaturing.
DNase treatments, which should be a standard practice in RNA extraction protocols, typically require heating samples to 37 C for enzyme activation. Therefore, 28S fragmentation will likely occur as part of the extraction protocol in all organisms with a UAAU-rich loop. Our findings also suggest that the large ribosomal subunit of some species may be resistant to denaturation as millipedes were able to maintain 28S integrity regardless of heat exposures to 90 C. This is possibly a result of a missing UAAU loop region similar to aphids (Hemiptera, Insecta; [28,42]). Our results further show that if samples are exposed to temperatures beyond 70 C for enzymes requiring heat deactivation, RNA integrity can be severely impacted. However, it appears that DNases can be largely removed with a second round of alcohol precipitation; this allowed us to avoid high temperature incubations by employing an enzyme brand that does not require temperature deactivation. Thus, we recommend avoiding enzymes or denaturing treatments that require high temperatures in order to ensure RNA integrity.

Denaturing gels versus electropherograms
Here, we show that both methods for evaluating RNA quality, denaturing gel electrophoresis and Bioanalyzer electropherograms, have their individual strengths while also being subjected to methodological artifacts. Though RIN values are unreliable in profiles where 28S is collapsed, electropherograms should still be used to assess RNA quality as they provide more accurate size estimates. Bleach agarose gels [36] can also be used for size estimation and to identify degradation or DNA contamination, in place of more toxic denaturing gels (e.g., Dimethyl sulfoxide (DMSO) and formaldehyde), to minimize time, effort, and costs. Furthermore, assessing RNA integrity via electrophoresis prior to diluting and running on a Bioanalyzer will allow researchers to more readily detect sources of degradation. For example, electropherograms for our gradient revealed degradation on a few lanes that were not initially present in the samples, indicating this occurred at some point during Bioanalyzer prep and/or transport. For these reasons, we recommend running samples on denaturing gels prior to or in conjunction with the bioanalyzer.

Arthropod extraction techniques
Although comparing multiple RNA extraction methods are outside the scope of this study, we found one in particular that worked well for arthropods and downstream RNAseq studies. We ultimately selected solvent-based (phenol-chloroform; e.g., TRIZOL) over column-based methods (e.g., QIAGEN and NUCLEOSPIN), as it was cost-effective while also yielding higher concentrations of high-molecular weight nucleic acids as well as small and micro-RNAs. Column membranes sequestered too much material and often resulted in the isolation of (1 lg of total RNA per sample. Obtaining high RNA yields becomes particularly important when working with small-bodied organisms, a challenge faced by many invertebrate researchers. Lastly, solvent-based methods prove to be reliable across taxa and laboratories as the RNA profiles reported here are consistent with past arthropod RNA studies (e.g., one band in spiders, [47]; three bands in crabs [48]).

Conclusions
Here we compared RNA profiles across Arthropoda from a diverse array of lineages including insects, crustaceans, arachnids, millipedes, and centipedes. We found that a majority of lineages displayed what would be incorrectly interpreted as a "degraded" RNA profile when compared with vertebrate references. Specifically, six distinct RNA banding patterns were revealed, providing insight into the high degree of variation that can be found across Arthropoda. However, 18S size (1900-2000 bp) remained constant among arthropods, making it a useful proxy for RNA integrity. Lastly, our study demonstrates that RNA is substantially labile and extraction protocol temperatures (>37 C) will denature 28S that contains the UAAU-rich loop region.

Data accessibility
Raw sequencing data are available on the National Center for Biotechnology Information (NCBI)'s Sequence Read Archive (SRA) database under Bioproject PRJNA476149.

Supplementary data
Supplementary data is available at Biology Methods and Protocols online.