Optimal reference genes for RNA tissue analysis in small animal models of hemorrhagic fever viruses

Reverse-transcription quantitative polymerase chain reaction assays are frequently used to evaluate gene expression in animal model studies. Data analyses depend on normalization using a suitable reference gene (RG) to minimize effects of variation due to sample collection, sample processing, or experimental set-up. Here, we investigated the suitability of nine potential RGs in laboratory animals commonly used to study viral hemorrhagic fever infection. Using tissues (liver, spleen, gonad [ovary or testis], kidney, heart, lung, eye, brain, and blood) collected from naïve animals and those infected with Crimean–Congo hemorrhagic fever (mice), Nipah (hamsters), or Lassa (guinea pigs) viruses, optimal species-specific RGs were identified based on five web-based algorithms to assess RG stability. Notably, the Ppia RG demonstrated stability across all rodent tissues tested. Optimal RG pairs that include Ppia were determined for each rodent species (Ppia and Gusb for mice; Ppia and Hrpt for hamsters; and Ppia and Gapdh for guinea pigs). These RG pair assays were multiplexed with viral targets to improve assay turnaround time and economize sample usage. Finally, a pan-rodent Ppia assay capable of detecting Ppia across multiple rodent species was developed and successfully used in ecological investigations of field-caught rodents, further supporting its pan-species utility.

of viral hemorrhagic fever (VHF) to validate RGs for three highly pathogenic zoonoses: Crimean-Congo hemorrhagic fever virus (CCHFV) in mice 11 , Nipah virus (NiV) in hamsters 12 , and Lassa virus (LASV) in guinea pigs 13 .Multiple tissues (liver, spleen, gonad [ovary or testis], kidney, heart, lung, eye, brain, and blood) were evaluated both from naïve and infected animals to ensure the validity of the RG.To determine the optimal RG, five web-based algorithms-RefFinder, Bestkeeper, NormFinder, geNorm, and delta C t (ΔC t ) methods-were used to assess RG stability (Fig. 1).

Biosafety
All work with infectious virus or infected animals was conducted in a biosafety level 4 (BSL-4) laboratory at the Centers for Disease Control and Prevention (CDC) following established BSL-4 standard operating procedures approved by the Institutional Biosafety Committee.The CDC Institutional Biosafety Committee approved all  1) Selection and design of reference gene (RG) assays against potential optimal RGs, followed by assay validation.(2) Selection of sample RNA for determining RG stability.Samples were selected from prospective studies with study IDs, number of male (♂) and female (♀) animals required, and target organs indicated.(3) RG analysis methods.First, RG stability is assessed in each individual tissue, then across all tissue sets; finally, paired gene analysis was carried out to find optimal RG pairs.

Origin of tissue RNA
The study aimed to use RNA from tissues harvested and archived during previous experimental studies and from tissues harvested from in-house colony animals humanely euthanized for study-independent health indications (e.g., advanced chronic disease).In this study, naïve animal tissue RNA was selected from animals that did not receive virus challenge.Infected animal tissue RNA samples were selected from animals that succumbed to lethal disease (NiV, CCHFV, or LASV) or displayed high viral RNA loads in tissues at endpoint (NiV) after infection with well characterized lethal challenge doses.
Infected guinea pig tissue RNA was obtained from experiments (Study# 3129-1) in which strain 13/N guinea pigs (in-house CDC colony; male and female; 203-1195 days old) were inoculated SC in the interscapular region with 10 4 TCID 50 of recombinant LASV strain Josiah 17 (CDC Virharv# 813752; GenBank HQ688673; HQ688675), 500 µL total volume per animal.Naïve animals (Study# 3193-3) were taken from the in-house colony (male and female; 879-1345 days old) and euthanized for humane reasons (e.g., age, injury, or underlying conditions).At terminal or pre-determined endpoint, guinea pigs were euthanized by isoflurane inhalation and intracardiac administration of sodium pentoparbital followed by confirmation of cardiac arrest.

Animal husbandry
All animals were housed under climate-controlled conditions (range: 68-79 °C; 30-70% humidity) with a 12 h day/night cycle.All infected animals were housed in isolator caging systems containing a HEPA-filtered inlet and exhaust air supply.Both naïve and infected mice were group-housed (5-6 per cage) on corncob bedding (Bed-o'Cobs ® ¼", Anderson Lab Bedding; Care Fresh, Healthy Pet; Enviro-Dry, Shepherd Specialty Papers) with cotton nestlets in an isolator caging system (Tecniplast GM500, West Chester, PA, USA), and provided sterilized, commercially available mouse chow and water ad libitum.Naïve and infected Syrian hamsters were group-housed (4-6 per cage) on corncob bedding, with wood blocks, seeds (given weekly), and plastic hide-aways, in an isolator caging system (Tecniplast GR900), and provided rodent chow (Teklad Rodent Diet 18% protein, Envigo) and water ad libitum.Naïve guinea pigs were housed in floor pens in BSL-2 segregated by sex, in a climate-controlled laboratory with a 12:12 h light cycle.Animals were provided with ad libitum feed (Guinea Pig Diet 5025) and water, timothy hay daily, and fresh vegetables 3 times a week.Infected guinea pigs were singly housed on soft pellets (paper chip; Carefresh), with Enviro-Dri wood blocks, in an isolator caging system (Tecniplast Blue Line 1500U), provided chow and water ad libitum as above, and both timothy hay and fresh produce daily.This study is reported in accordance with ARRIVE guidelines.

PubMed query search for specifc reference genes
To determine a relative level of how commonly RGs are used in literature, a PubMed search was conducted for each chosen RG.Queries included the common names and abbreviations of each specifc RG (e.g., "18S" or "18SRNA" or "18S RNA") and "housekeeping gene" or "reference gene" (Fig. 2B).

RNA extraction
RNA was extracted from either homogenized tissue (~ 1 mm 3 sections) or 50 μL of whole blood in lithium heparin (mice and hamsters), or from 50 to 125 μL whole blood in EDTA (guinea pigs), using the MagMax-96 Total RNA isolation kit (Thermo Fisher Scientific) on a 96-well ABI MagMax extraction platform.Samples were treated with DNase (Lucigen) and eluted in 75 μL MagMax elution buffer.RNA was quantified using the Nanodrop 8000 spectrometer (Thermo Fisher Scientific).RNA purity was assessed and was found to have an average absorbance of 1.99 (A 260 /A 280 ) with a standard deviation of 0.17.Extracted RNA was kept at − 80 °C for long-term storage.

RT-qPCR
The equivalent of 100 ng RNA from each tissue type selected was analyzed by RT-qPCR to assess RG stability.These RT-qPCRs were performed using the OT-2 liquid handling platform (Opentrons) to reduce pipetting and volume errors.Assays were carried out using either in-house desgined primer sets (Tables 1 and 2) or a commercially available Eukaryotic 18S rRNA Endogenous Control (Thermo Fisher).RT-qPCRs were run on the CFX384 or CFX96 Real-Time PCR Detection systems (BioRad) with 10 μL reaction volume using the SuperScript III Platinum One-Step qRT-PCR kit (Thermo Fisher).Each PCR reaction was performed in duplicate.Thermal cycling conditions were as follows: 50 °C reverse transcription step for 15 min, 95 °C denaturation for 2 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s. Primer-probe set sequences and concentrations can be found in Tables 1 and 2.

Assay validation
Singleplex assay validation was performed on RNA harvested from the spleen of naïve hamsters, mice, or guinea pigs diluted serially tenfold from 100 ng/μL.In addition, singleplex and multiplex assay validations were performed on RNA harvested from the spleens of animals infected with CCHFV (mice), NiV (hamsters), or LASV (guinea pigs) diluted serially tenfold from 100 ng/μL.Efficiency calculations and inter-and intra-assay variability were determined from sample standard curves generated from at least five technical replicates, utilizing a minimum of three thermocyclers and three operators to account for both human and machine-derived variability.RT-qPCR efficiency was calculated from the mean of each dilution series, and the inter-and intra-assay variability were assessed using percentage correlation of variance (% CV).

Variability analysis
C t values were determined using the Bio-Rad CFX Manager 3.1 software by setting a standard baseline threshold level for each primer probe set.Analysis was carried out on mean C t values of duplicate reactions.The stability of each RG was assessed using the web-based RefFinder tool (http:// blooge.cn/ RefFi nder/?type= refer ence) 21,22 .This web-based platform also incorporates the following methods: comparative ΔC t method 23 , Bestkeeper 24 , NormFinder 25 , and geNorm 26 .Analyses were initially carried out using data sets from each individual tissue, followed by analysis across all tissue types.Analyses were carried out on naïve animals and infected animals separately, before combining for further evaluation.

RNA isolation and RT-qPCR from wild rodents
Field-caught rodents were identified by morphological characteristics or through molecular amplification of mitochondrial DNA cytochrome B gene.RNA was extracted from 100 µL rodent blood using TRI Reagent (Sigma-Aldrich) according to the manufacturer's instructions.RT-qPCRs were performed using the Luna Probe One-Step RT-qPCR kit (New England BioLabs) with 20 μL reaction volume containing 5 µL extracted RNA.Thermal cycling conditions were as follows: 55 °C reverse transcription step for 10 min, 95 °C denaturation for 1 min, followed by 45 cycles of 95 °C for 10 s and 60 °C for 30 s.Samples were considered positive if values obtained were less than 35 C t .

Statistical analyses
Data distribution was assessed using the Anderson-Darling normality test and confirmed to have normal distribution using a significance level of α = 0.05.Statistical analyses were carried out in GraphPad Prism v9.4.1.

Evaluation of eight species-specific reference gene assays
Commercial assays are available for select commonly used RGs, such as 18S.However, the availably of commercial assays for other candidate RGs is limited, especially for use across multiple rodent species.Here, we developed species-specific assays against 8 candidate RGs that are commonly used in the literature 1,23 (B2m, Tbp, Gapdh, Actb, Hrpt, Sdha, Ppia, and Gusb) (Fig. 2B).Efficacy and variance of these assays were validated by analyzing standard curves generated by tenfold serial dilution of total RNA harvested from the spleen of naïve animals (Figs.1-1 and 2C).Spleen RNA from naïve (uninfected) animals was chosen as the test material due to the consistently high concentrations of total RNA isolated from this tissue across all species.All species-specific singleplex assays were calculated to have efficiencies between 90 and 105% and the intra-and inter-assay variances were calculated to be < 3% CV and < 2% CV for all assays, respectively.Assay parameters are shown in Table 3.These in-house species-specific assays have now been evaluated to accurately amplify targets and therefore are appropriate to use to determine relative expression levels of RGs in tissues.

Determining relative expression levels of candidate reference genes for use on tissue RNA extracted from naïve or VHF-infected rodents
Once the candidate RGs were validated, we then determined expression of these 8 in-house species-specific RG assays and a pan-species commercial assay targeting 18S RNA in RNA isolated from liver, spleen, ovary, testis, kidney, heart, lung, eye, and brain (all species), plus mammary gland and bladder (guinea pigs).Furthermore, to ensure that disease did not alter RG expression, gene levels were evaluated in tissues from naïve (uninfected), infected, or both naïve and infected animals for each virus/rodent model pair: CCHFV-infected mice, NiVinfected hamsters, and LASV-infected guinea pigs.RG expression was analyzed by RT-qPCR using 100 ng RNA harvested from tissue.Sample size for each virus/rodent model pair and infection status was n = 4 (2 male and 2 female) for all tissues except gonad; for gonad tissue, n = 8 (4 males and 4 females).

Comparative stability of reference genes across tissue types in naïve and VHF-infected mice, hamsters, and guinea pigs
To determine the stability of RGs, we applied 5 commonly used analysis methods available through a web-based platform: RefFinder analysis 21,22 , the comparative ΔC t method 23 , Bestkeeper analysis 24 , Normfinder analysis 25 , and geNorm analysis 26 .These methods are readily utilized across studies of RG analysis and are freely available.These analyses provide a statistical basis to measure the variability of each gene's expression and output this as a value of stability.Here, the term stability will be used to describe this variability.A stable gene would be given a low score, indicating that it has minimal variability across the tissues and experimental conditions tested.RefFinder analysis uses a combination method in which the scores determined from each of the other analysis methods are weighted and the stability value is calculated from the geometric mean.Comparative ΔC t compares C t values and determines the relative expression of pairs of genes within the sample.Bestkeeper analyzes gene stability through the standard deviation of sample C t values.Normfinder provides a stability value (S), determined as the lowest variation between the inter-and intra-group variations.GeNorm analysis ranks gene stability based Table 3. Reference gene assay primer efficiency and variability.Assay parameters were determined from standard curves generated by tenfold dilutions of 100 ng/µL total RNA isolated from spleens of naïve mice, hamsters, or guinea pigs.Intra-and inter-assay variability was assessed using percentage correlation of variance (% CV).

Inter-assay variability (%)
Mouse  Due to the extensive range of tissues routinely examined in viral pathogenesis and medical countermeasure studies, we aimed to determine the most appropriate RG to implement across multiple tissue types in naïve and virus-infected animals for each of the previously investigated virus/rodent model pairs.Prior to assessing all tissues collectively, we first assessed the stability of each RG expression by tissue type to ensure that there were no tissue-specific outliers.We found that, occasionally, specific tissues in specific animal species (e.g., blood from mice and lungs from guinea pigs) had more variation in RG expression, but overall in individual tissues, irrespective of species or infection status, most RGs had acceptable stability scores (Fig. S2).Exceptions included B2m, which often ranked low amongst the tested RGs (particularly in mouse tissues), and 18S, which consistently demonstrated unstable expression scores across multiple tissue types and species using the 5 analysis methods.For 18S, this finding, and the previously noted wide range of C t values (Fig. S2, Table S1), indicated extensive variation and instability in gene expression; thus, 18S was removed from further analysis.
After ensuring consistency across tissue types, we carried out stability analyses utilizing C t values obtained collectively across all tissue RNAs from naïve and infected animals for each of the virus/rodent model pairs.Using this approach, we determined that the RGs demonstrating highest levels of stability across all tissues were: Ppia (determined by RefFinder, ΔC t , Normfinder, and geNorm) and Gapdh (Bestkeeper) in mice; Ppia (Ref-Finder, ΔC t , and Normfinder), Actb (Bestkeeper), and Hrpt (geNorm) in hamsters; and Ppia (RefFinder, ΔC t , and Normfinder), Tbp (Bestkeeper), and Gapdh (geNorm) in guinea pigs (Fig. 4).Ppia was the most stable RG when analyzed across multiple tissues and infection conditions among the rodent species examined here.We also note that other RGs demonstrating satisfactory stability scores could be used as acceptable alternates to Ppia (Fig. 4).

Optimal pairing of reference genes for tissue RNA analysis of naïve and VHF-infected rodents
We identified several RGs that may be appropriate for normalizing RT-qPCR data in these infection models when used alone; however, the use of multiple RGs is recommended to provide greater accuracy for RT-qPCR normalization than using a single RG 2,26 .Therefore, we identified RGs pairs using the same web-based analysis software as described above utilizing C t values obtained from the analyses of tissue RNA from naïve and infected animals for each virus/rodent model pair.As Ppia was determined to be the most suitable RG candidate across multiple species, we focused on RG pairs that included Ppia to provide a level of continuity across the different rodent models used in our studies (Fig. 5).

Validation of triplex assays for the detection of viral target and reference genes in VHF-infected rodent tissue samples
As animal studies can generate vast numbers of samples requiring analysis, we aimed to multiplex the optimal paired RG assays identified here with specific VHF viral targets to allow rapid sample processing and conservation of material.A pan-rodent Ppia assay was designed to detect Ppia across mouse, hamster, and guinea pig tissues, targeting a conserved region of this gene (Fig. 6A).The pan-Ppia assay demonstrated 90-105% efficiency across these tissues, with intra-and inter-assay variations of < 5% CV (Table 4).Ppia assays (species-specifc or pan-rodent) were combined with our previously determined optimal species RG assays and with assays targeting CCHFV, NiV, and LASV.
To detect CCHFV in mouse samples, the CCHFV nucleocapsid protein (NP) gene assay was combined with mouse Gusb and mouse Ppia or pan-rodent Ppia.To detect LASV in guinea pig samples, the LASV NP gene assay was combined with guinea pig Gapdh and guinea pig Ppia or pan-rodent Ppia.To detect NiV in hamster samples, the NiV nucleoprotein (N) gene assay was only combined with with hamster Ppia or pan-rodent Ppia, as attempts to multiplex with hamster Hrpt were unsuccessful due to primer-primer interactions.These multiplexed assays were validated by RT-qPCR on tenfold serially diluted RNA isolated from tissues of infected animals to generate a standard curve.Concurrently, singleplex assays were carried out for each target of interest to ensure that C t values differed < 1 C t between singleplex and multiplex assays (Fig. 6B).All single-and multiplexed assays were calculated to have efficiencies of 90-105%, and the intra-and inter-assay variations were calculated to be < 5% CV. Assay parameters can be seen in Table 4.

Pan-rodent Ppia assay can be used to detect rodent RNA in ecological field studies
To further the utility of the pan-rodent Ppia assay, we generated a sequence logo 28 of 26 rodents Ppia gene sequences from data available from Genbank (Table 2, Fig. S4A).The pan-rodent Ppia assay primers demonstrated a high level of complementarity to these rodent sequences.Therefore, we theorized they could be used as positive control RGs during rodent sampling in the field, for example, during epidemiological studies of rodentborne viruses.In a retrospective study of rodent viral reservoirs in Argentina, 31 blood specimens were analyzed representing the different rodent genotypes and their distribution in the country.Ppia-positive RT-qPCR results were obtained for Calomys boliviae, Holochilus chacarius, Necromys lasiurus, Oligoryzomys sp., (including O. nigripes, O. chacoensis, O. flavescens, and O. longicaudatus), and Oxymycterus rufus.The only rodent species that was negative on the pan-rodent Ppia assay was Akodon azarae, although other Akodon sp.(only identified to genus level) demonstrated positive results (Fig. S4B).These results suggest that the pan-rodent Ppia assay is able to detect Ppia across a wide range of rodents and may have further utility outside of the research laboratory.

Discussion
In vivo data regarding RG stability have only been described for a limited number of animal models of infectious disease thus far, including murine cytomegalovirus in the mouse model 9 ; Marek's disease virus infection in chickens 29 ; influenza in mice (with kidney yang deficiency syndrome) 30 ; and vesicular stomatitis virus in Syrian golden hamsters 31 .The impact of VHF infection on RG expression in animal models has not yet been reported, likely due to limited facilities and the time-consuming nature of conducting in vivo experiments in the high-containment environment.We aimed to utilize historical archived tissue RNA from CCHFV, NiV and LASV studies to identify RGs suitable for normalization when performing RT-qPCR analyses of RNA samples harvested from small animal models of VHF.Five web-based platforms (RefFinder, Bestkeeper, ∆C t , Normfinder,  Viral infections are known to modify the expression of host genes, including those often used as RGs.In lungs of influenza A-infected mice, expression of Gapdh and Actb becomes less stable than in lungs from naïve animals 30 ; we did not observe this destabilization during CCHFV infection.Additionally, assessment of RGs in naïve and VSV-infected Syrian hamster characterized B2m as more stable and Hrpt as less stable 31 , whereas we found Hrpt to be more stable and B2m less stable in naïve and NiV-infected Syrian hamsters.These differences underscore the importance of assessing RG stability under the appropriate experimental conditions due to differential modifications of host gene expression induced by different viruses.
Our work highlights how choice of RG alone may introduce discrepancy in data.Advances in animal model development are often hampered by experimental design-associated variability in outcome and virological indices.We found that commonly used RGs, including 18S and B2m, may be unsuitable for use in rodent models and potentially other species as well, as they exhibited large variations within and across tissues.This variability could result in further challenges when comparing data within or between experimental groups or between studies in which these RGs were used.In contrast, we established that the RG Ppia was stably expressed across all tested samples from naïve and infected mice, hamsters, and guinea pigs and serves as a suitable choice for use in these species.Ppia, though not as commonly used an RG as 18S, B2m, Gapdh, Actb, or Tbp 1 , has also been previously identified as a stably expressed gene in many small animal models, including mice 7,32-34 and rats 35 , and has been posited as a preferred RG in the analysis of human tissue 36,37 .Often a single RG is used in tissue analyses.To further increase confidence in samples analyzed with the use of Ppia, we determined optimal gene pairs to use as RGs for each species: Ppia and Gusb for mice, Ppia and Hrpt for hamsters, and Ppia and Gapdh for guinea pigs.
To our knowledge, this is the first report to investigate stable expression of RGs to identify appropriate controls for normalizing RT-qPCR data in studies of VHF infection in animal models.Furthermore, we demonstrate the validation of RG and viral target multiplex assays in each model species, providing a strong framework for further research into VHF pathogenesis and therapeutics evaluation.Finally, we describe development and use of Ppia as a pan-rodent assay for wide application in both laboratory animals and wild rodent species.Table 4. Multiplex reference gene assay primer efficiency and variability.Assay parameters were determined from standard curves generated by tenfold dilutions of 100 ng/µL total RNA isolated from spleens of CCHFVinfected mice, NiV-infected hamsters, or LASV-infected guinea pigs.Each assay was run as singleplex or was multiplexed (CCHFV and LASV: viral targets with species-specific RG and pan-rodent Ppia or species-specific Ppia; NiV: viral target with pan-rodent Ppia or species-specific Ppia).Intra-and inter-assay variability was assessed for both singleplex and multiplex assays using % CV.

Figure 1 .
Figure 1.Sample selection and analysis flow chart.(1) Selection and design of reference gene (RG) assays against potential optimal RGs, followed by assay validation.(2) Selection of sample RNA for determining RG stability.Samples were selected from prospective studies with study IDs, number of male (♂) and female (♀) animals required, and target organs indicated.(3) RG analysis methods.First, RG stability is assessed in each individual tissue, then across all tissue sets; finally, paired gene analysis was carried out to find optimal RG pairs.

Figure 2 .
Figure 2. Reference gene targets chosen for gene stability investigation to use in RT-qPCR normalization.(A) Reference gene function and expression.Data were collated from the GeneCards Suite 18 and the Mouse Genome Informatics Gene Expression Database 19 .(B) Number of results obtained from PubMed queries including the common names and abbreviations for each specifc RG (e.g., "18S" or "18SRNA" or "18S RNA") and "housekeeping gene" or "reference gene." (C) Standard curves were generated by making tenfold dilutions of 100 ng/µL total RNA isolated from the spleen of naïve mice, hamsters, or guinea pigs.Each point represents the mean of 8 independent RT-qPCR reactions, and error bars indicate the standard deviation from the mean.

Figure 3 .
Figure 3. Variation of reference gene expression across all tissues isolated from hamsters, mice, and guinea pigs.Expression profiles of 8 RGs in 100 ng of total RNA across all tissues collected from either naïve animals; animals infected with NiV-M (hamsters), CCHFV (mice), or LASV (guinea pigs); or naïve + infected combined for each of the virus/rodent models.Cycle threshold (C t ) values (n = 36-132) are presented as extended violin plots with median values indicated by a line.

Figure 4 .
Figure 4. Stability of reference genes in RNA harvested from hamster, mouse, and guinea pig tissues.Stability of RGs in 100 ng total RNA across multiple tissue types from mice (total replicates per gene: n = 76), hamsters (total replicates per gene: n = 80), and guinea pigs (total replicates per gene: n = 94) was assessed using the webbased algorithms RefFinder, Bestkeeper, ΔC t , Normfinder, and geNorm.Lower gene stability values indicate higher gene stability.Dashed line indicates the greatest acceptable stability value, where appropriate.

Figure 5 .
Figure 5. Stability of reference gene pairs in RNA harvested from hamster, mouse, and guinea pig tissues.Stability of RG pairs in 100 ng total RNA across multiple tissue types from mice (total replicates per gene: n = 76), hamsters (total replicates per gene: n = 80), and guinea pigs (total replicates per gene: n = 94) was assessed using the web-based algorithms RefFinder, Bestkeeper, ΔC t , Normfinder and geNorm.Lower gene stability values indicate higher gene stability.Dashed line indicates the greatest acceptable stability value, where appropriate.

Figure 6 .
Figure 6.Viral assays can be multiplexed with RG assays.(A) Multiple sequence alignment of Ppia coding sequences from mice (Genbank #NM_008907), hamsters (Genbank #XM_005086775), and guinea pigs (Genbank #XM_003465805) using ClustalW 27 .Asterisk (*) denotes conserved sequences.Locations of panspecies Ppia primers and probe are indicated.(B) Standard curves generated by tenfold dilutions of 100 ng/µL total RNA isolated from the spleen of CCHFV-infected mice, NiV-infected hamsters, or LASV-infected guinea pigs.Each assay was run as singleplex or was multiplexed (CCHFV and LASV: viral targets with species-specific RG and pan-rodent Ppia or species specific Ppia; NiV: viral target with pan-rodent Ppia or species-specific Ppia).Standard curves from both singleplex and multiplexed assays are shown.Each point represents the mean of 5-12 independent RT-qPCR reactions.

Table 1 .
Reference gene assay primer sequences for each rodent species.

Table 2 .
Pan -rodent Ppia and viral target multiplex assay primer sequences.GenBank accession numbers used to produce a consensus sequence for designing primers and probes against Ppia from multiple rodent species were NM_008907 (mouse), XM_005086775 (hamster), and XM_003465805 (guinea pig).Primer sequences are indicated in the 5′-3′ direction.The final concentration (conc.) of each individual primer or probe per reaction is indicated in micromolar (μM).