The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.

Introduction {#s1}
============

*Neurospora crassa* is a filamentous fungus that has served for nearly 50 years as the model organism for studying circadian pathways and the molecular clock in eukaryotes [@pone.0112706-Baker1]. It is classified under the phylum Ascomycota and is widely distributed in nature [@pone.0112706-Baker1], typically found growing on vegetation killed by fire [@pone.0112706-Jacobson1]. As such, it is able to synthesize and secrete high levels of all main enzyme types involved in lignocellulosic biomass degradation and is also able to convert pentose sugars, cellulose polymers, and agro-industrial residues to ethanol [@pone.0112706-Dogaris1]. Its genome is predicted to contain twice as many cellulase genes in comparison the current state-of-the-art industrial strain used in biofuel production, *Trichoderma reesei* [@pone.0112706-Dogaris1]. Additionally, recent work has shown that the variety of molecular, genetic, and biochemical techniques available for *N. crassa* can advance the analyses of fungal deconstruction of plant biomass [@pone.0112706-Glass1]. Therefore, great interest in the use of *N. crassa* for biotechnological applications has sparked in recent years.

Gene expression profiling is an effective means of studying an organism's response to its environment, a necessary component in identifying the genes and regulatory mechanisms associated with biological processes in *N. crassa* such as circadian pathways and important biotechnological applications such as cellulase degradation. Recent work in our lab has utilized *N. crassa* maintained in continuous, long-term bioreactor cultures in an effort to understand the triggers of cellulase production. Most *N. crassa* studies are based on either racetube or batch culture conditions. Understanding the complexity of *N. crassa* growth in bioreactors is lacking and more importantly, the identification of optimal reference genes for molecular analysis needs to be addressed.

Reverse-transcription quantitative PCR (RT-qPCR) is a high-throughput means of studying the expression of a large suite of genes over time and under different conditions. With this technique, two quantification strategies are available to measure the level of expressed genes: absolute and relative quantification. Absolute quantification relates the fluorescent PCR signal of the sample to input copy number using an external calibration curve, while relative quantification measures the relative change in mRNA expression levels [@pone.0112706-Pfaffl1]. Relative expression is based on normalization of target gene expression to an internal control gene, often referred to as a "reference" or "housekeeping" gene, that should not vary in the cells under investigation or in response to environmental or experimental conditions [@pone.0112706-Vandesompele1]. The selection of a suitable reference gene is critical for gene expression studies using relative quantification [@pone.0112706-Pfaffl1]. However, these reference genes, thought to be constitutively expressed and minimally regulated, are commonly used for expression profiling in RT-qPCR assays without *a priori* conformation that mRNA levels do not change under the experimental conditions being investigated. The use of reference genes necessitates that they be validated for the specific experimental set-up and it is typically necessary to choose more than one [@pone.0112706-Bustin1]. Work by Vandesompele *et al* [@pone.0112706-Vandesompele1] demonstrated that the typical procedure in which a single reference gene is used for normalization can lead to relatively large errors [@pone.0112706-Vandesompele1], and it is now common to employ a combination of genes. Thus, in order to understand the metabolic processes when *N. crassa* is grown under different environmental conditions, the selection and evaluation of reference genes must be addressed.

*N. crassa* is a morphologically complex multicellular organism, possessing at least 28 distinct cell types [@pone.0112706-Bistis1] between the vegetative (asexual) and sexual phases. In bioreactors, it is maintained in the vegetative phase. The primary cells occurring in the vegetative phase include hyphae, macro- and microconidia, all of which are comprised of several different types [@pone.0112706-Bistis1]. One of the biological differences among these cell types is the number of nuclei per cell: the macroconidia are multinucleated, possessing on average between 3--6 nuclei per cell, while microconidia possess only a single nuclei ([Figure 1](#pone-0112706-g001){ref-type="fig"}). This life cycle trait must be taken into account when performing gene expression studies, particularly under continuous, long-term culture conditions such as those conducted in bioreactors. In general, our bioreactors are dominated by clusters of aggregated colonies dispersed throughout the vessel, similar in cell type composition to the mats observed in batch cultures yet comprised of more discrete units. Therefore, it is necessary to determine whether changes in gene expression are the result of true fluctuations in transcripts, or whether an increase or decrease in transcripts (i.e. gene expression) is actually a reflection in gene copy number fluctuation, which would occur when cells are a mix of macro- and micro-conidia. One means of assessing gene expression changes is to determine the number of transcripts normalized to a certain value such as cell number, volume of sample, or mass. This number can then be related to the gene copy number normalized to the same units to provide a ratio of the transcript to gene ratio. These data standardize the transcripts to the gene copy number and so determine whether the gene is changing expression, or whether transcript variation is a result of gene variation. This type of situation applies to cultures that occur as multiple cell types during which some morphologies are multinucleated.

![Generalized lifestyle of *N. crassa* during which different cell types occur.\
Center image illustrates bioreactor set-up. Within the bioreactors utilized in this study, *N. crassa* was present as both multi- and uninucleated cells.](pone.0112706.g001){#pone-0112706-g001}

While numerous gene expression studies have been performed with *N. crassa*, little to no information is provided as to reference gene validation, and the extent of gene expression fluctuation as a result of changes in gene copy number (i.e. as cells progress from hyphae to macroconidia) has not been addressed. Additionally, genes that are commonly used as reference genes in other organisms have been shown to be under clock control in *N. crassa*. In general, the actin gene is frequently used as a normalizer in gene expression studies [@pone.0112706-Sun1]--[@pone.0112706-Tralau1], though in most cases there is no mention of stability verification. Therefore, the overall objective of this study was to examine a suite of genes for the potential to serve as appropriate reference genes under a variety of environmental conditions during growth in *N. crassa*. We evaluated a set of 12 genes under different environmental conditions using a combination of both relative and absolute quantification. Using this combination allowed for a more comprehensive assessment of fluctuations in expression, as the transcript copy number per mL was calculated and its ratio to gene copy number determined. Additionally, the well-established software programs BestKeeper [@pone.0112706-Pfaffl2] and NormFinder [@pone.0112706-Andersen1] were used to determine expression stability and rank of all genes to assess their suitability as reference genes.

Materials and Methods {#s2}
=====================

Organism Maintenance {#s2a}
--------------------

*N. crassa* cultures were maintained on agar slants (1.5% w/v sucrose with Vogel's medium [@pone.0112706-Vogel1]) at −20°C for up to 1 month. Slants were prepared following standard protocols. Briefly, the medium comprised of Vogel's salts, 1.5% sucrose, and 1.5% agar was autoclaved and cooled to ≤50°C whereupon the trace minerals and biotin solutions (1.0×10^−5^ M) were added. Agar was distributed among test tubes and cooled at an angle, resulting in a "slanted" agar medium. Strain X-661 (*ras-1bd* mutant with the *frq:luc* transcriptional reporter) was used for these studies.

Bioreactor culture conditions and cell harvesting {#s2b}
-------------------------------------------------

New slants were inoculated from the frozen stock culture at the start of each experiment. *N. crassa* X-661 grew on slants for 2--3 days at 30°C in constant darkness. Conidia were isolated from slants by washing the slant culture with ∼1 mL of Vogel's medium [@pone.0112706-Vogel1] lacking a carbon source. The collected material was centrifuged (10,000 rpm, 1 min) and the resulting conidial pellet was washed 3 times with Vogel's medium (no C-source). The final pellet was re-suspended and inoculated into 60 mL of Vogel's medium with 1% glucose, 0.2% Junlon, (100,000 MW, Sigma Aldrich) and 0.008% Antifoam 204 (Sigma Aldrich). The culture was grown at 30°C with shaking for 24 hrs with light. This culture was inoculated into a BioFlo/CelliGen 115, 1.3 L bioreactor ([Figure 2](#pone-0112706-g002){ref-type="fig"}) containing the same medium composition and grown in batch for 24 hrs with light. Continuous flow into and export from the culture was started as the cells reached the exponential growth phase. When the cells reached steady state (∼11 hrs later) the light remained on for 11 hrs and then was shut off. At this point, the cells were grown in one of three conditions: continuous all-dark conditions, under light/dark cycling, or under temperature flux (25°C/37°C) in constant darkness. The light/dark and temperature flux oscillated with 22-hr periods (11 hrs light, 11 hrs dark) in keeping with the period length of *N. crassa* circadian oscillations. During the experiments, *N. crassa* cultures were maintained in steady-state conditions, with pH 5.5 and 10% dissolved O~2~. The O~2~ was bubbled through the medium at 1.0 L/min. A dynamic agitation program was used to keep the cells suspended in bioreactor cultures. The program controlled the rotation rate of the propeller in the reactor and was set up as a 10-min cycle using the BioCommand Batch Control Plus BioProcessing software (New Brunswick Scientific, \#M1326-0051, revision B). The propeller rotated at 1,000 rpm for 1 min, 100 rpm for 1 min, and then 400 rpm for 8 min. This program was started ∼1--2 hrs following bioreactor inoculation. The temperature was either maintained at 25°C or cycled (with 11-hr periods) between 25°C and 37°C. The bioreactors were equipped with an external heating jacket and an internal cooling finger and temperature probe to accurately maintain the temperature. A temperature cycling program (BioCommand software) was used for varied temperature experiments. The growth state of the cells (i.e., log phase, steady-state) was monitored using the CO~2~ output. A CO~2~ gas sensor (Vernier) was used to detect CO~2~ produced by the culture. The CO~2~ data was recorded using the LoggerLite software (version 1.6.1). Samples were removed under darkroom light conditions (safelight) and the 1.5 mL aliquots were flash frozen in centrifuge tubes and stored at −80°C. Biological replicates, in the form of individual bioreactors, were performed for each type of environmental condition, with a total of 8--10 time points collected per condition.

![Schematic of bioreactor design and operation.](pone.0112706.g002){#pone-0112706-g002}

DNA extraction {#s2c}
--------------

Genomic DNA was extracted from biological replicates using the DNeasy Plant Mini kit (Qiagen, Valencia, CA) following the manufacturer's instructions, with slight modifications. Briefly, cells from 1 mL of a bioreactor sample were disrupted using a mortar and pestle under liquid nitrogen. Following addition of 400 µL Buffer AP1 and 4 µL RNase A, the slurry was ground an additional 1--2 min. The entire volume was then transferred to a 1.5 mL microfuge tube, vortexed vigorously, and 400 µL of this slurry removed for DNA extraction. The remaining steps followed the manufacturer's protocol. DNA was eluted in nuclease-free water rather than Buffer AE. DNA concentration and purity were assessed using the Nanodrop-2000 (Thermo Fisher Scientific, Wilmington, DE). DNA used in quantitative PCR (qPCR) assays was diluted 1∶5 in nuclease-free water prior to use.

RNA isolation and cDNA synthesis {#s2d}
--------------------------------

Total RNA was extracted frozen from a 1 mL bioreactor sample using the RNeasy Plant Mini kit (Qiagen, Valencia, CA) following the manufacturer's protocol, with slight modifications. Briefly, cells were disrupted using a mortar and pestle under liquid nitrogen. Four hundred fifty microliters of buffer RLC was added and the slurry mixed until pipetteable. The sample was vortexed vigorously and 450 µL added to the QIAshredder spin column. Subsequent steps followed the manufacturer's protocol. The DNase digestion was extended to 1 hr. An additional spin at top speed (16,000 rpm) for 1.5 min was included immediately prior to elution in RNase-free water. Total RNA concentration and purity was measured using the Nanodrop-2000.

Total RNA was converted to cDNA using the High-Capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA). The protocol was optimized to include less template RNA and a larger overall reaction volume, as preliminary screening with several of the potential reference genes revealed greater fluctuations in C~T~s when using the recommended 2 µg in a 20 µL volume (data not shown). The following protocol was used for all reverse transcriptase reactions: total RNA was first diluted 1∶5 in nuclease-free water, and 500 ng added to the reaction. Each reaction contained: 2 µL 20x enzyme mix, 2x RT buffer, a volume equivalent to 500 ng total RNA, and brought to a final volume of 40 µL with nuclease-free water. Controls consisted of all components except the enzyme, in which case an equal volume of nuclease-free water was added. Reactions were incubated in a 2720 thermocycler (Applied Biosystems, Foster City, CA) under the following conditions: 37°C for 60 min followed by heating to 95°C for 5 min. Prior to use in qPCR, samples were again diluted 1∶5 in nuclease-free water.

Candidate reference gene primer design and validation {#s2e}
-----------------------------------------------------

Quantitative PCR assays were developed for a suite of candidate reference genes ([Table 1](#pone-0112706-t001){ref-type="table"}) based on the SYBR Green chemistry [@pone.0112706-Morrison1]. For each assay, primers were designed that targeted an ca. 100--130 bp region of the coding sequence of each gene ([Table 1](#pone-0112706-t001){ref-type="table"}). When possible, primers were designed to span exon-exon boundaries. Both coding and genomic sequences of each gene were obtained from the *N. crassa* database located on the Broad Institute website (<http://www.broadinstitute.org/annotation/genome/neurospora/GeneIndex.html>). Primers were designed using the Primer3 software (<http://bioinfo.ut.ee/primer3/>), with an annealing temperature of 60°C. Primers were examined for potential secondary structures using Mfold [@pone.0112706-Zuker1]. The potential for self and cross dimer formation was examined for each of the primer sets using an online oligonucleotide analysis tool (<https://www.operon.com/oligos/toolkit.php>). Self-complementarity was determined using the oligonucleotide properties calculator (<http://www.basic.northwestern.edu/biotools/oligocalc.html>). Primers were obtained from Invitrogen (Carlsbad, CA). Primer sets were optimized over concentrations spanning 150--300 nM, with an annealing/extension temperature of 60°C. Product specificity was determined via melt curve analysis. The PCR efficiency of each assay was determined using LinRegPCR [@pone.0112706-Ramakers1].

10.1371/journal.pone.0112706.t001

###### List of candidate reference genes.

![](pone.0112706.t001){#pone-0112706-t001-1}

  GeneName    Accession \#             Description                              Primer Sequence(5′-3′)                    AmpliconLength   T~M~    PCR *E*   FinalConcentration
  ---------- -------------- --------------------------------- ---------------------------------------------------------- ---------------- ------- --------- --------------------
  *act*         NCU04173               beta-actin                               F:GATCGGTATGGGCCAGAAGG                         110         79.10    1.91           300 nM
                                                                                R:CGTCCCAGTTGGTAACGACA                                                      
  *acl*         NCU06783            ATP citrate lyase          F:AGATCCTGATCCCCGTCGACC[\*](#nt101){ref-type="table-fn"}        120         80.55    1.92           300 nM
                                                               R:AACCTCGGCAGAGGTGGCATC[\*](#nt101){ref-type="table-fn"}                                     
  *adk*         NCU00414            adenosine kinase                            F:TCCTCCACGATGCTGTCAAG                         148         83.82    1.86           300 nM
                                                                                R:TCAAGGTCGTAGTGGTTGGC                                                      
  *asl*        NCU09119.7    ATP synthaselyase gamma subunit                    F:TACCAACGCCAAGGACATCG                         114         81.35    1.89           300 nM
                                                                                R:GAGTAGTCGCCCTTGAGCTG                                                      
  *btl*         NCU04054              beta-tubulin                              F:CCACTTCTTCATGGTCGGCT                         103         83.21    1.83           150 nM
                                                                                R:CTTGGGGTCGAACATCTGCT                                                      
  *frh*         NCU03363       Frq-interactingRNA helicase                      F:TAGCAAGAATGCCAGGTGGG                         139         82.35    1.92           300 nM
                                                                                R:GATCTCGTCAACGGCTTTGC                                                      
  *l6*          NCU02707      60S subunitribosomal protein                      F:AGAAGGAGGTTTCCAGCAGC                         148         82.25    1.90           300 nM
                                                                                R:CACTTCATCTCGTGGGGCTT                                                      
  *tbp*        NCU044770         TATA binding boxprotein                        F:GGTGCCAAGTCCGAAGATGA                         136         79.38    1.94           300 nM
                                                                                R:TGGGGAACTTGATGTCGCAA                                                      
  *trk*        AJ009758.1             K transporter                           F:AGTACACTATACCCTCATCATTG                        110         78.40    1.89           300 nM
                                                                                R:CCCCAGCAGCGAAGAAGAG                                                       
  *vma1*        J03955.1        vacuolar ATPase subunit 1                       F:CCAACAAAATGGCGCCGAG                          111         80.45    1.89           300 nM
                                                                               R:CCAACTCATACATAGCAACACC                                                     
  *vma2*        J03956.1        vacuolar ATPasesubunit 2                        F:GTCGTCCAGGTCTTCGAGG                          127         78.90    1.89           150 nM
                                                                               R:TGCCGGAACCATCAAAGATACG                                                     
  *vma3*       NCU01332.7       vacuolar ATPasesubunit 3                        F:ACAGGACCACTACGCTCTCT                         114         80.08    1.90           300 nM
                                                                                R:ACACCAGCATCACCGACAAT                                                      

\*Primers from [@pone.0112706-Wang1].

Construction of external standards for absolute quantification {#s2f}
--------------------------------------------------------------

The expression of a subset of the candidate reference genes was also assessed via absolute quantification of the ratio of transcripts to gene copies per mL of sample. This was achieved though the creation of external calibration curves. Plasmid-based standards were constructed for the actin (*act*), ATP synthase lyase γ subumit (*asl*), adenosine kinase (*adk*), and vacuolar ATPase subunit 3 (*vma3*) genes and used for absolute quantification of both transcript and gene copy number. Standards were created via amplification from genomic DNA using gene-specific primers ([Table S1](#pone.0112706.s010){ref-type="supplementary-material"}). DNA was used rather than RNA as it has been shown that external calibration curves constructed with DNA rather than RNA serve as better models for mRNA quantification due to higher sensitivity, increased quantification range, higher reproducibility, and is more stable than the RNA calibration curve [@pone.0112706-Pfaffl3]. PCRs were conducted in 25 µL reaction volumes with the PCR Reagent System (Invitrogen, Carlsbad, CA) in a 2720 thermocycler (Applied Biosystems, Foster City, CA) and contained (as final concentration): 1x PCR buffer (minus Mg), 2 mM MgCl~2~, 0.25 mM each dNTP, 200 nM each forward and reverse primer, 0.5 U Taq polymerase, and brought to a final volume of 25 µL with nuclease-free water. Thermocycling conditions consisted of 94°C, 3 min; 35 cycles of 94°C, 45 s, 51°C (48°C for *act*) for 30 s, 72°C, 1 min, and a final extension at 72°C, 7 min. PCR products were visualized by gel electrophoresis to confirm product amplification of the proper size. The resulting amplicons were cloned into the pCR4 vector using the TOPO TA cloning kit (Invitrogen, Carlsbad, CA) and transformed into chemically competent TOP10 *E. coli* following the manufacturer's instructions. Cells were spread onto Miller's Luria Broth (LB) plates containing 50 µg mL^−1^ kanamycin (LBkan~50~) and incubated overnight at 37°C. Colonies were screened for inserts of the expected size via the gene-specific PCRs described above. Colonies producing a band of the expected size were grown overnight in 2 mL LBkan~50~ broth. Plasmids were isolated using the QIAprep Spin Miniprep kit (Qiagen, Valencia, CA) and sequenced using the M13 forward and reverse primers by GeneWiz, Inc. (Germantown, MD). Sequence identity was confirmed using the standard nucleotide BLAST function on the NCBI website, with organism specificity set to "Neurospora." The number of copies per µL was calculated using the equation: X g µL^−1^DNA/\[PCR amplicon + plasmid length×660\])×6.022×10^23^. The plasmid stock solution was diluted in nuclease-free water to a working concentration of 5×10^8^ copies per µL. Serial 10-fold dilutions were used to create external calibration curves that spanned eight orders of magnitude ranging from 5×10^7^ to 5×10^−1^ copies per µL. Two µL of each dilution were used per reaction, yielding a standard curve that ranged from 1×10^1^ to 1×10^8^ copies. The standard curve for each assay was obtained by plotting the log of the calculated copy number against the cycle at which fluorescence for that sample crossed the threshold value (cycle threshold, C~T~) ([Figures S1](#pone.0112706.s001){ref-type="supplementary-material"}, [S2](#pone.0112706.s002){ref-type="supplementary-material"}, [S3](#pone.0112706.s003){ref-type="supplementary-material"}, and [S4](#pone.0112706.s004){ref-type="supplementary-material"}). The slope, y-intercept, r^2^ value, and PCR efficiency of each assay can be found in the Supporting Information ([Figures S1](#pone.0112706.s001){ref-type="supplementary-material"}, [S2](#pone.0112706.s002){ref-type="supplementary-material"}, [S3](#pone.0112706.s003){ref-type="supplementary-material"}, and [S4](#pone.0112706.s004){ref-type="supplementary-material"}).

Quantitative PCR {#s2g}
----------------

The Fast SYBR Green Master Mix (Applied Biosystems, Foster City, CA) was used for all qPCR assays. Each reaction contained: 10 µL 2x Fast SYBR Green master mix, 150--300 nM each forward and reverse primer, 5 ng (2 µL) template cDNA, and brought to a final volume of 20 µL with nuclease-free water. Reactions were conducted in 96-well plates on a ViiA 7 with the Fast 96-well block system (Applied Biosystems, Foster City, CA). The following protocol was used for all assays: an initial 20 s incubation at 95°C, followed by 40 cycles of 95°C for 1 s and 60°C for 20 s, followed by a melt curve analysis of 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s to determine product specificity. All qPCR reactions were performed in triplicate using Applied Biosystems MicroAmp Fast 96-well reaction plates (0.1 mL with barcode) sealed with MicroAmp optical adhesive film. No-template controls were also included in each amplification run to monitor for contamination. Reactions were recorded and analyzed using the Applied Biosystems ViiA 7 System software. The raw C~T~ values and intra-assay variations for both relative and absolute quantification assays are provided in [Data S1](#pone.0112706.s013){ref-type="supplementary-material"} and [Data S2](#pone.0112706.s014){ref-type="supplementary-material"}.

Data analysis {#s2h}
-------------

Three types of environmental conditions were used to determine the appropriateness of a suite of genes for use as reference genes: total darkness, in which samples were collected over a time span of 70 hrs; light/dark cycling, in which samples were collected from alternating periods of light and dark over the span of ca. 70 hrs; and temperature fluctuations in darkness, in which samples were collected from alternating cycles of 25°C and 37°C over the span of ca. 70 hrs. Expression levels of all twelve candidate reference genes were initially determined from their C~T~ values. The expression stability of all candidate reference genes under each of the different conditions and among all conditions combined was assessed using the BestKeeper v. 1 and Normfinder software programs. For data input into NormFinder, C~T~s of each gene at each time point were converted to relative quantities using the equation Q = *E* ^ΔC^ ~T~, where ΔC~T~ = C~T~ highest abundant sample -- C~T~ sample, *E* = PCR efficiency, and Q = relative quantity.

The number of transcripts and gene copies per mL were calculated via external calibration curves for *act*, *adk*, *asl*, and *vma3*. Transcript or gene copy numbers from total RNA or DNA extracted from 1 mL bioreactor samples were calculated automatically from the regression lines by the ViiA 7 System software. Gene copies were calculated using an equation described previously [@pone.0112706-Ritalahti1], in which gene copies per mL sample = (gene copies per reaction mix x volume of DNA \[µL\])/(DNA per reaction mix \[µL\]×mL sample used) and adjusted to account for the 1∶5 initial dilution factor. The following equation was used to calculate transcript copies per mL: \[(copies per µL per reaction mix×RNA dilution factor)/(volume µL total RNA added to RT reaction/RT reaction volume µL×cDNA dilution factor\]×vol total RNA eluted (µL) ([Data S2](#pone.0112706.s014){ref-type="supplementary-material"} and [Data S3](#pone.0112706.s015){ref-type="supplementary-material"}). The transcript to gene ratio was determined by dividing transcript copies per mL by gene copies per mL. The Shapiro-Wilkes test was used to test for normality of both the ratio and C~T~ values ([Table S2](#pone.0112706.s011){ref-type="supplementary-material"}). If both data sets passed the test for normality, the Pearson product moment correlation test was applied to examine the strength of the correlation between the two. If one of the data sets failed the test for normality, the Spearman rank order correlation test was applied.

Results and Discussion {#s3}
======================

Quantitative PCR assay design, optimization, and general characteristics {#s3a}
------------------------------------------------------------------------

Typical reference genes include GAPDH, albumin, actins, tubulins, and small and large subunit rRNA [@pone.0112706-Pfaffl1], as they are present in all nucleated cell types and are necessary for basic cell survival. However, numerous studies, across different organisms, have shown that the typical reference genes are regulated and vary under different conditions [@pone.0112706-Pfaffl1]. Based on the experimental evidence thus far, it appears that no reference gene retains sufficient overall expression stability suitable for all contexts and that the best candidates differ between tissue type [@pone.0112706-Hruz1] and environmental condition. Therefore, it is recommended that the individual investigator choose a reference gene that is most stably expressed under their specific conditions, such as tissue type or stage of development, and should be impervious to perturbations such as external stimuli, disease, or genetic modification [@pone.0112706-Pfaffl1], [@pone.0112706-Hruz1]--[@pone.0112706-Gutierrez1].

Several approaches exist for qPCR data normalization. One approach recently implemented has been to select genes from a genome-wide background derived from large sets of microarray data. This approach has been used for *Arabadopsis*, multiple human cell lines and cancers, and mouse [@pone.0112706-Hruz1]. However, while *N. crassa* serves as a model organism for studying eukaryotic circadian pathways, it is not included in this database. Additionally, studies utilizing whole-transcriptome profiling such as microarrays and RNA sequencing are able to select appropriate reference genes based on the transcriptomic profiling. For example, genes encoding an adenosine kinase and the ATP synthase gamma chain have been used as reference genes in the filamentous fungus *Penicillium glabrum* under conditions of heat stress, as their expression was shown via initial microarray analysis to be unaffected by this condition [@pone.0112706-Nevarez1]. Additionally, a recent study evaluated 12 potential reference genes in the yeast *S. cerevisiae* based on pooled data from publicly-available transcriptome studies. The majority of selected genes were not previously reported as reference genes and the results revealed that these newly-identified genes outperformed the commonly-used reference genes [@pone.0112706-CankorurCetinkaya1]. There is extremely limited gene expression data with *N. crassa* growth in bioreactors, such that large data sets from which to select potential housekeeping genes do not exist. Also lacking are publicly-available whole-transcriptome data sets such as has been used for other organisms when selecting reference genes. Therefore, potential reference genes were selected based on the following features: previous transcriptomic-based studies from *N. crassa* (*acl*) [@pone.0112706-Wang1] or other fungi (*adk*, *asl*) [@pone.0112706-Nevarez1] in which the expression of these genes remained unchanged; via selection from the Broad Institute *N. crassa* database (<http://www.broadinstitute.org/annotation/genome/neurospora/MultiHome.html>) with genes for which no expression data existed (*tbp*), or from previous studies of genes whose results pertaining to their function and regulation indicated the potential for use as a reference gene (*vma1*, *vma2*, *vma3*, trk) [@pone.0112706-Wechser1], [@pone.0112706-Rivetta1] Additionally, *N. crassa* genes employed as reference genes in previous studies were also evaluated here (*act*, *btl*, *l6*, *frh*) [@pone.0112706-IshShalom1], albeit with newly-designed primers for use in the *btl* and *frh* assays.

The compliance of the RT-qPCR assays with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments [@pone.0112706-Bustin2]) guidelines is shown in the MIQE checklist ([Table S3](#pone.0112706.s012){ref-type="supplementary-material"}). All assays developed in this study were specific for their intended gene target, as indicated by a single peak in the melt curve analysis ([Figure S5](#pone.0112706.s005){ref-type="supplementary-material"}). Additionally, all possessed similar PCR efficiencies ([Table 1](#pone-0112706-t001){ref-type="table"}), so that if all had displayed equal stability, they could have been used interchangeably in relative quantification methods. A preliminary assessment of the 12 reference genes showed that average C~T~ values, when assimilating the data from all genes, ranged from 21.47--30.34 under light/dark fluctuations ([Figure 3A](#pone-0112706-g003){ref-type="fig"}), 22.53--30.72 under all dark conditions ([Figure 3B](#pone-0112706-g003){ref-type="fig"}), and 19.93--28.73 under temperature flux ([Figure 3C](#pone-0112706-g003){ref-type="fig"}). Based on this initial analysis, the two genes displaying the highest average C~T~ values, *frh* and *trk*, were discarded from further analysis. This was due to the recommendation of reference genes being expressed at comparable levels to both other reference genes and target genes [@pone.0112706-Bustin3]. Additionally, *frh* has been shown to be under clock control [@pone.0112706-Dong1].

![Boxplot of C~T~ values of candidate reference genes during light/dark cycling, temperature fluctuation, and continuous darkness.\
A line across the box depicts the median. The box indicates the 25th and 75th percentiles. Whiskers represent the 90th and 10th percentiles.](pone.0112706.g003){#pone-0112706-g003}

Reference gene stability under light/dark conditions {#s3b}
----------------------------------------------------

The expression stability of candidate reference genes was initially estimated based on calculated variations (SD, standard deviation of the C~T~, and CV, the coefficient of variance expressed as a percentage on the C~T~ level) using BestKeeper. With this program, reference genes can be ordered from the most stably expressed, exhibiting the lowest variation, to the least stable. In addition to stability (SD and CV) values, additional key parameters include the Pearson correlation coefficient (*r*), in which pair-wise correlation analyses are performed to estimate inter-gene relationships among all reference genes, and their associated probability (*p*) values. Thus, *r* and *p* describe the relation between the BestKeeper index and each gene. The candidate reference genes were then ranked as to their suitability under one environmental condition and across all three using the program NormFinder. This program utilizes a model-based approach that, in addition to assessing overall expression variation, allows for the samples to be assigned to different groups, allowing for the intra- and intergroup expression variation of the candidate genes to be assessed; this grouping feature was used when assessing reference gene stability under light/dark and temperature cycling. It was not employed for all-dark growth at constant temperature. The comprehensive results of the BestKeeper and NormFinder programs are provided in [Data S4](#pone.0112706.s016){ref-type="supplementary-material"} and [Data S5](#pone.0112706.s017){ref-type="supplementary-material"}, respectively.

The greatest variability among the three environmental conditions was observed under light/dark fluctuation, as evidenced by both the range in C~T~s and the SD values obtained with BestKeeper ([Table 2](#pone-0112706-t002){ref-type="table"}). In general, most genes displayed overall stability with respect to light/dark cycling. Three genes (*btl*, *asl*, and *vma1*) displayed the greatest expression stability (SD = 1.10, 1.13, and 1.21, respectively), while *l6* and *acl* were the least stable (SD = 1.41). However, based on the Pearson correlation coefficients (*r*), they correlated well with the other genes and thus were retained for further analysis per previous methods [@pone.0112706-Pfaffl2]. Overall, genes displayed similar Pearson correlation coefficients, ranging between 0.936--0.997. NormFinder identified *vma1* as the most stable gene ([Table 2](#pone-0112706-t002){ref-type="table"}), with the best combination of two genes being *vma1* and *asl*. The similarity in stability values between *vma1* and *asl* was also reflected in the intragroup variation, which was 0.001 in the light for both. However, *vma1* expression was less influenced under dark conditions than *asl*, with an intragroup variation of 0.002 versus the 0.016 calculated for *asl* ([Table 2](#pone-0112706-t002){ref-type="table"}).

10.1371/journal.pone.0112706.t002

###### Descriptive statistics and stability values of candidate reference genes under conditions of light/dark cycling (n = 8).

![](pone.0112706.t002){#pone-0112706-t002-2}

                              *btl*       *asl*      *vma1*       *act*      *vma2*       *tbp*       *adk*      *vma3*       *acl*       *l6*
  ------------------------ ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- -----------
  **geo Mean \[C~T~\]**       24.60       22.72       26.24       21.15       24.82       25.21       24.77       23.33       25.85       22.00
  **ar Mean \[C~T~\]**        24.63       22.75       26.28       21.19       24.86       25.25       24.81       23.38       25.90       22.06
  **min \[C~T~\]**            22.16       21.44       24.85       19.60       23.21       23.88       22.84       21.62       23.73       20.51
  **max \[C~T~\]**            26.76       24.91       28.59       23.47       27.05       27.86       26.89       25.75       28.17       25.14
  **SD \[±C~T~\]**          **1.10**    **1.13**    **1.21**    **1.24**    **1.24**    **1.26**    **1.34**    **1.38**    **1.41**    **1.41**
  **CV \[% C~T~\]**           4.482       4.967       4.608       5.834       4.990       5.004       5.381       5.881       5.449       6.397
  ***r***                     0.936       0.997       0.997       0.988       0.995       0.987       0.970       0.985       0.954       0.985
  ***p***                     0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001
  **NF Stability value**    **0.194**   **0.070**   **0.067**   **0.154**   **0.114**   **0.169**   **0.175**   **0.252**   **0.278**   **0.274**
  **Int Var - Light**         0.062       0.001       0.001       0.033       0.001       0.157       0.140       0.001       0.134       0.342
  **Int Var - dark**          0.105       0.016       0.002       0.039       0.002       0.008       0.032       0.026       0.145       0.028

NF = NormFinder, Int Var = Intragroup Variation as calculated by Normfinder.

Reference gene stability under all-dark conditions {#s3c}
--------------------------------------------------

Based on the SD as calculated by BestKeeper, all genes were stable during complete darkness, as the SD for all genes was \<1 ([Table 3](#pone-0112706-t003){ref-type="table"}). The only exception to this was *vma3*, in which the SD = 1.07; however, the Pearson correlation coefficient was acceptable, and in fact higher than some of the other genes that displayed a lower SD than *vma3*. Genes displaying the lowest SD, and thus the greatest stability, were *vma1*, *btl*, and *tbp*. All genes were significantly correlated to the BestKeeper index (p\<0.001), and the Pearson correlation coefficient was very similar among genes. In general, the analysis showed a strong correlation for all 10 genes. NormFinder identified *tbp*, *vma2*, and *vma1* as being the most stable of the genes under conditions of complete darkness, with *tbp* and *vma2* recording almost identical stability values (0.141 and 0.143, respectively) ([Table 3](#pone-0112706-t003){ref-type="table"}).

10.1371/journal.pone.0112706.t003

###### Descriptive statistics and stability values of candidate reference genes under all-dark conditions (n = 10).

![](pone.0112706.t003){#pone-0112706-t003-3}

                             *vma1*       *btl*       *tbp*      *vma2*       *acl*       *adk*       *asl*       *act*       *l6*       *vma3*
  ------------------------ ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- -----------
  **geo Mean \[C~T~\]**       27.07       24.84       25.22       26.15       27.16       26.46       24.96       22.50       23.18       24.90
  **ar Mean \[C~T~\]**        27.08       24.85       25.24       26.17       27.17       26.48       24.98       22.53       23.21       24.93
  **min \[C~T~\]**            25.44       23.32       23.20       24.27       25.11       25.15       23.40       20.43       21.14       22.97
  **max \[C~T~\]**            28.22       26.19       26.58       27.63       28.79       27.78       26.33       24.29       25.09       26.55
  **SD \[±C~T~\]**          **0.489**   **0.584**   **0.699**   **0.729**   **0.739**   **0.777**   **0.899**   **0.932**   **0.989**   **1.074**
  **CV \[% CP\]**             1.807       2.349       2.768       2.787       2.718       2.936       3.600       4.134       4.260       4.308
  ***r***                     0.931       0.886       0.98        0.986       0.943       0.734       0.934       0.983       0.990       0.971
  ***p***                     0.001       0.001       0.001       0.001       0.001       0.016       0.001       0.001       0.001       0.001
  **NF Stability Value**    **0.169**   **0.181**   **0.141**   **0.144**   **0.211**   **0.499**   **0.398**   **0.268**   **0.271**   **0.261**

NF = NormFinder.

Reference gene stability during temperature fluctuation {#s3d}
-------------------------------------------------------

All genes were stably expressed under conditions of temperature fluctuations, as reflected by the SD ranging from 0.501--0.865. Overall, SDs ranged between 0.501--0.741, with the exceptions of *l6* (0.796) and *acl* (0.865). Based on the SD, *asl* (0.501), *btl* (0.616), and *vma1* (0.622) exhibited the greatest stability ([Table 4](#pone-0112706-t004){ref-type="table"}). As had been observed for all-dark conditions, all genes correlated very strongly with the index, with most displaying *r*\>0.9. Under conditions of fluctuating temperature, Normfinder identified *vma1* as the most stable, with the best combination being *vma1* and *act* ([Table 4](#pone-0112706-t004){ref-type="table"}). *vma1* displayed extremely low variation under 25°C (0.0001) with an increase in variance observed at 37°C for both *vma1* and *act* (0.026 and 0.018, respectively).

10.1371/journal.pone.0112706.t004

###### Descriptive statistics and stability values of candidate reference genes under temperature fluctuations (n = 12).

![](pone.0112706.t004){#pone-0112706-t004-4}

                              *asl*       *btl*      *vma1*       *adk*       *act*       *tbp*      *vma2*      *vma3*       *l6*        *acl*
  ------------------------ ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- -----------
  geo Mean \[C~T~\]           22.51       24.60       25.55       24.93       19.91       24.19       25.52       23.74       20.70       25.55
  ar Mean \[C~T~\]            22.52       24.62       25.56       24.95       19.93       24.20       25.54       23.76       20.72       25.57
  min \[C~T~\]                21.37       23.31       24.27       23.49       18.60       22.94       23.83       22.03       19.04       23.85
  max \[C~T~\]                23.87       26.01       26.69       26.66       21.38       25.68       27.10       25.37       22.64       27.36
  **SD \[±C~T~\]**          **0.501**   **0.616**   **0.622**   **0.673**   **0.675**   **0.693**   **0.718**   **0.741**   **0.796**   **0.865**
  CV \[% C~T~\]               2.225       2.504       2.435       2.698       3.387       2.863       2.811       3.117       3.841       3.381
  *r*                         0.953       0.865       0.970       0.792       0.978       0.906       0.837       0.998       0.944       0.950
  *p*                         0.001       0.001       0.001       0.002       0.001       0.001       0.001       0.001       0.001       0.001
  **NF Stability value**    **0.133**   **0.175**   **0.073**   **0.262**   **0.087**   **0.133**   **0.175**   **0.090**   **0.165**   **0.190**
  Int Var 25°C                0.022       0.074       0.000       0.012       0.015       0.006       0.115       0.011       0.153       0.066
  Int Var 37°C                0.006       0.041       0.026       0.266       0.018       0.093       0.017       0.003       0.050       0.271

NF = NormFinder, Int Var = Intragroup Variation as calculated by Normfinder.

Reference gene stability across all conditions {#s3e}
----------------------------------------------

Assimilating the results of both BestKeeper and NormFinder under each individual condition, *vma1* always ranked among the top three genes in terms of stability; *asl* and *btl* also consistently ranked in the top three in at least one of the analyses. Also, in keeping with their biological role as subunits of the same enzyme, *vma1*, *vma2*, and, under several conditions, *vma3* ranked closely together in terms of stability, suggesting they can be used interchangeably. However, in order to accurately identify the most stable genes under all three conditions, all data were combined and analyzed with both BestKeeper and Normfinder. The results of BestKeeper revealed that 4 of the 10 genes displayed acceptable stability across all conditions (reflected by a SD\<1): *btl*, *tbp*, *vma2*, and *vma1* ([Table 5](#pone-0112706-t005){ref-type="table"}). Conversely, the two genes frequently used for normalization, *l6* and *act*, displayed the greatest variance ([Table 5](#pone-0112706-t005){ref-type="table"}). The most appropriate gene for use across all conditions as identified via NormFinder was *vma1*, with the best combination consisting of *vma1* and *asl*, with a combined stability value of 0.089. *tbp* and *vma2* were the next most-stable genes among the 10, in agreement with the results obtained with BestKeeper ([Table 5](#pone-0112706-t005){ref-type="table"}).

10.1371/journal.pone.0112706.t005

###### Descriptive statistics and stability values of candidate reference genes under light/dark cycling, all-dark growth, and temperature fluctuations (n = 30).

![](pone.0112706.t005){#pone-0112706-t005-5}

                              *btl*       *tbp*      *vma2*      *vma1*       *adk*      *vma3*       *acl*       *asl*       *act*       *l6*
  ------------------------ ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- ----------- -----------
  geo Mean \[C~T~\]           24.68       24.80       25.54       26.23       25.39       24.01       26.16       23.35       21.08       21.85
  ar Mean \[C~T~\]            24.70       24.83       25.57       26.26       25.42       24.05       26.19       23.40       21.13       21.91
  min \[C~T~\]                22.16       22.94       23.21       24.27       22.84       21.62       23.73       21.37       18.60       19.04
  max \[C~T~\]                26.76       27.86       27.63       28.59       27.78       26.55       28.79       26.33       24.29       25.14
  **SD \[±C~T~\]**          **0.750**   **0.927**   **0.956**   **0.982**   **1.088**   **1.137**   **1.193**   **1.257**   **1.341**   **1.415**
  CV \[% C~T~\]               3.036       3.733       3.739       3.741       4.280       4.729       4.555       5.370       6.344       6.459
  *r*                         0.810       0.903       0.880       0.973       0.892       0.951       0.961       0.939       0.955       0.962
  *p*                         0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001       0.001
  **NF Stability value**    **0.191**   **0.137**   **0.178**   **0.129**   **0.230**   **0.243**   **0.279**   **0.148**   **0.185**   **0.218**
  Int Var: Dark Flux          0.129       0.003       0.007       0.004       0.022       0.019       0.129       0.023       0.040       0.025
  Int Var: Light Flux         0.089       0.136       0.001       0.001       0.150       0.001       0.166       0.002       0.025       0.300
  Int Var: All Dark           0.040       0.024       0.020       0.032       0.254       0.058       0.051       0.156       0.066       0.067
  Int Var: 25°C               0.074       0.006       0.115       0.000       0.012       0.011       0.066       0.022       0.015       0.153
  Int Var: 37°C               0.041       0.093       0.017       0.026       0.266       0.003       0.271       0.006       0.018       0.050

NF = NormFinder, Int Var = Intragroup Variation as calculated by Normfinder.

Absolute quantification of transcript and gene copy numbers {#s3f}
-----------------------------------------------------------

In bioreactors, *N. crassa* exists as different vegetative cell types, some of which are multi-nucleated. Based on the parameters [@pone.0112706-Vandesompele1], [@pone.0112706-Pfaffl2] defined for appropriate reference gene selection, all 10 candidate reference genes displayed fluctuations in expression. This prompted the question: Were fluctuating C~T~s due to changes in gene expression, or were gene (DNA) copy numbers changing, so that if transcript levels remained relatively stable, it was manifested as expression fluctuations? Therefore, expression changes in a subset of the candidate genes was assessed using absolute quantification. The transcript (RNA) and gene (DNA) copy numbers of *act*, *adk*, *asl*, and *vma3*, were calculated and used as an additional means of examining gene expression. These four were selected based on combination of factors since they are involved in different biological processes within the cell and ranged in their stability values under the conditions examined. For example, *asl* frequently ranked as one of the most stable, while *adk* was one of the least stable. External standards were created and used to calculate both transcript and gene copy numbers per mL of sample in order to determine the ratio of transcript copies to gene copies over time (representative image in [Figure 4](#pone-0112706-g004){ref-type="fig"}, plots of transcript and gene copy numbers for all genes under all conditions [Figures S6](#pone.0112706.s006){ref-type="supplementary-material"}, [S7](#pone.0112706.s007){ref-type="supplementary-material"}, [S8](#pone.0112706.s008){ref-type="supplementary-material"}, and [S9](#pone.0112706.s009){ref-type="supplementary-material"}). This ratio indicated whether C~T~ fluctuations were due to changes in transcript abundance, indicative of changes in gene expression, or whether it was a result of changes in gene copy number, such as could occur in multinucleated cells.

![Representative figure of data used to obtain transcript-to-gene ratios.\
*act* all-dark transcript and gene copies per mL of sample were calculated via absolute quantification using external calibration curves. Error bars represent the standard deviation derived from triplicate qPCRs from biological replicates.](pone.0112706.g004){#pone-0112706-g004}

Plotting both the RNA: DNA ratio and C~T~ values over time indicate that C~T~ changes are due to fluctuations in transcript copy number ([Figure 5](#pone-0112706-g005){ref-type="fig"}). The Pearson Product Correlation test was applied to the ratio and corresponding C~T~ of each time point for each gene to test the hypothesis that a negative correlation exists between ratio and C~T~ (i.e. as the ratio increases, the C~T~ decreases). With all four genes under all three conditions, a strong negative correlation was found to exist between ratio and C~T~ values, demonstrating that C~T~ changes serve as a reliable reflection of transcript, and not gene copy number fluctuations ([Table 6](#pone-0112706-t006){ref-type="table"}). For example, under conditions of light/dark cycling, the *act* C~T~ value decreased at the second time point ([Figure 5A](#pone-0112706-g005){ref-type="fig"}), indicative of an increase in transcripts. If this were due to changes in gene (DNA) copy number, the RNA: DNA ratio would not increase. However, an increase in ratio indicates an increase in transcript abundance. As statistically significant negative correlations were found to exist between the transcript: gene ratio and C~T~ values among the four genes under all conditions, these data demonstrate that even with the presence of different cell types, relative quantification is an acceptable method for measuring gene expression changes in *N. crassa* during growth in bioreactors.

![Correlation between gene expression and C~T~ values as derived from absolute quantification of transcript and gene copy numbers.\
Standard curves were created for the *act*, *adk*, *asl*, and *vma3* genes and used to calculate the transcript-to-gene ratios under conditions of light/dark cycling, temperature fluctuation, and continuous darkness. Top panel depicts light/dark cycling, with yellow highlight indicating light on. Bottom panel displays conditions of temperature flux between 25°C and 37°C, with 37°C illustrated by orange highlight.](pone.0112706.g005){#pone-0112706-g005}

10.1371/journal.pone.0112706.t006

###### Pearson Product Correlation values obtained for ratio and CT of *act*, *adk*, *asl*, and *vma3* under each condition.

![](pone.0112706.t006){#pone-0112706-t006-6}

  Gene/Condition                          *r*                     p-value
  -------------------- ----------------------------------------- ---------
  *act* Light/Dark                      −0.747                     0.033
  *act* All Dark                        −0.868                     0.001
  *act* Temperature     −0.685[\*](#nt106){ref-type="table-fn"}    0.012
  *asl* Light/Dark                      −0.880                     0.004
  *asl* All Dark                        −0.765                     0.010
  *asl* Temperature                     −0.721                     0.008
  *adk* Light/Dark                      −0.924                     0.001
  *adk* All Dark        −0.636[\*](#nt106){ref-type="table-fn"}    0.040
  *adk*Temperature                      −0.721                     0.008
  *vma3* Light/Dark                     −0.921                     0.001
  *vma3* All Dark       −0.939[\*](#nt106){ref-type="table-fn"}    0.000
  *vma3* Temperature                    −0.748                     0.004

\*indicates Spearman rank value.

Conclusions {#s4}
===========

The work presented here mirrors the results obtained in recent studies with multiple organisms in which the most appropriate reference genes for use in RT-qPCR were not those commonly employed when using the relative quantification method, underscoring the importance of proper reference gene evaluation and selection. The increased use of RT-qPCR coupled with the realization of the importance of selecting the appropriate reference gene has led to dedicated studies on reference gene selection in both model and non-model organisms, including rats [@pone.0112706-Langnaese1], [@pone.0112706-Lardizabal1], zooplankton [@pone.0112706-Spanier1], fruit crops [@pone.0112706-Zhong1] (*Litchi chinensis*), and insects [@pone.0112706-Bansal1]--[@pone.0112706-Shen3]. In some cases, no single gene was found to be stable under all conditions, with stability affected by both endogenous and exogenous factors such as cell type, developmental stage, diet, and environmental stressors [@pone.0112706-Shen2], [@pone.0112706-Shen3]. Additionally, reference gene studies have been undertaken in recent years with multiple genera of filamentous fungi. A study evaluating potential reference genes in *Trichoderma reesei*, similar in industrial importance to *N. crassa* due to its high secretory capacity for hydrolases, identified *sar1*, encoding a small GTPase, as most stable, whereas the gene coding for actin (*act*), did not rank among the best validated ones [@pone.0112706-Steiger1]. In the filamentous fungus *Fusarium graminearum*, an economically-important plant pathogen, evaluation of 15 genes previously identified as housekeeping genes or those selected from the whole transcriptome sequencing data under different culture conditions found the most appropriate reference genes to vary based on growth stage and/or toxin production [@pone.0112706-Kim1] and were comprised of genes coding for a cell surface flocculin, ubiquitin C-terminal hydrolase, the eukaryotic translation initiation factor 1 alpha, and the mitochondrial ribosomal protein S16.

In this study, we identified and validated new reference genes for use in *N. crassa* expression studies under multiple environmental conditions, including long-term growth in constant darkness, light/dark cycling, and temperature flux in constant darkness. In addition to different environmental conditions, we assessed the time-course variation in reference gene expression by sampling at frequent intervals over a period of \>60 hours. These diverse conditions and continuous culture conditions span primary developmental (mycelia, macro- and microconidia) and/or physiological stages, thus increasing the reliability of the reference genes selected in this study. Overall, the results demonstrate that, under the culture conditions described here, *vma1*, *vma2*, *tbp*, and *btl* serve as the most appropriate reference genes for use in qPCR. We also demonstrated that *l6*, and to a lesser extent *act*, often used in the relative quantification method, are less suitable as reference genes in comparison to multiple other genes.

Supporting Information {#s5}
======================

###### 

**Standard curve of** ***act*** **assay.** A plasmid-based standard was constructed for the actin gene via amplification from genomic DNA using actin-specific primers. The number of copies per µL was calculated using the equation: X g µL-1 DNA/\[PCR amplicon+plasmid length×660\])×6.022×1023. Serial 10-fold dilutions were used to create external calibration curves that spanned eight orders of magnitude ranging from 1×10^1^ to 1×10^8^ copies. The standard curve for each assay was obtained by plotting the log of the calculated copy number against the cycle at which fluorescence for that sample crossed the threshold value. The slope, y-intercept, r^2^ value, and PCR efficiency of the actin assay is provided in the figure.

(TIF)

###### 

Click here for additional data file.

###### 

**Standard curve of** ***adk*** **assay, constructed as described for the act assay.**

(TIF)

###### 

Click here for additional data file.

###### 

**Standard curve of** ***asl*** **assay, constructed as described for the act assay.**

(TIF)

###### 

Click here for additional data file.

###### 

**Standard curve of** ***vma3*** **assay, constructed as described for the act assay.**

(TIF)

###### 

Click here for additional data file.

###### 

**Melt curve profiles of the 12 reference gene assays evaluated in this study.**

(TIF)

###### 

Click here for additional data file.

###### 

***act*** **transcript and gene copies over time under conditions of light/dark cycling, temperature flux, and continuous darkness.** Transcript and gene copies per mL of sample were calculated via absolute quantification using external calibration curves as described in the text. Error bars represent the standard deviation derived from triplicate qPCRs from biological replicates.

(TIF)

###### 

Click here for additional data file.

###### 

***adk*** **transcript and gene copies over time under conditions of light/dark cycling, temperature flux, and continuous darkness.** Transcript and gene copies per mL of sample were calculated via absolute quantification using external calibration curves as described in the text. Error bars represent the standard deviation derived from triplicate qPCRs from biological replicates.

(TIF)

###### 

Click here for additional data file.

###### 

***asl*** **transcript and gene copies over time under conditions of light/dark cycling, temperature flux, and continuous darkness.** Transcript and gene copies per mL of sample were calculated via absolute quantification using external calibration curves as described in the text. Error bars represent the standard deviation derived from triplicate qPCRs from biological replicates.

(TIF)

###### 

Click here for additional data file.

###### 

***vma3*** **transcript and gene copies over time under conditions of light/dark cycling, temperature flux, and continuous darkness.** Transcript and gene copies per mL of sample were calculated via absolute quantification using external calibration curves as described in the text. Error bars represent the standard deviation derived from triplicate qPCRs from biological replicates.

(TIF)

###### 

Click here for additional data file.

###### 

**Primers used in construction of** ***act*** **,** ***adk*** **,** ***asl*** **, and** ***vma3*** **external standards for absolute quantification.**

(PDF)

###### 

Click here for additional data file.

###### 

**Normality assessment of transcript: gene ratio and C~T~ values of** ***act*** **,** ***adk*** **,** ***asl*** **, and** ***vma3*** **genes by Shapiro-Wilkes test.**

(PDF)

###### 

Click here for additional data file.

###### 

**MIQE checklist.**

(XLS)

###### 

Click here for additional data file.

###### 

**Raw C~T~ values and intra-assay variations for 12 reference genes assessed using relative quantification.**

(XLSX)

###### 

Click here for additional data file.

###### 

**Raw C~T~ values and intra-assay variations for reference genes assessed using absolute quantification.**

(XLSX)

###### 

Click here for additional data file.

###### 

**Sample description, total RNA concentration, dilution factor, and volume added to reverse transcription reaction.**

(XLSX)

###### 

Click here for additional data file.

###### 

**Comprehensive output of BestKeeper analysis.**

(XLSX)

###### 

Click here for additional data file.

###### 

**Comprehensive output of NormFinder analysis.**

(XLSX)

###### 

Click here for additional data file.

We thank the National Research Council for K.D.C.'s and A.L.C.'s postdoctoral research associateships. We also thank Kwangwon Lee (Rutgers University) and Christian Hong (University of Cincinnati) for providing the *N. crassa* strain, initial batch culture techniques, and meaningful discussions. The views expressed are those of the authors and do not reflect the official policy or position of the Department of Defense or the U.S. Government.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: KDC LAF JCB. Performed the experiments: KDC LAF RKP ALC ERP JCB. Analyzed the data: KDC. Contributed reagents/materials/analysis tools: KDC LAF JCB. Contributed to the writing of the manuscript: KDC JCB.
