This work was supported by the National Institutes of Health grant 5R21AI081041-02 and Georgetown University.

1. Larval Rearing of Two A. albopictus Groups to Adulthood
<ol>
	<li>Set two photoperiod cabinets with programmable lighting at 21 &#176;C for optimal diapause expression16 and approximately 80% relative humidity.
	<ol>
		<li>Program one cabinet for a 16L:8D light:dark cycle (a non-diapause inducing LD photoperiod). Set the second cabinet for an 8L:16D light:dark cycle (diapause inducing)13.</li>
		<li>Program &#8216;lights on&#8217; at the same time in both cabinets to synchronize circadian time between photoperiods.&#160;</li>
	</ol>
	</li>
	<li>Calculate the quantity of eggs needed to perform the experiment. Aim for 300-500 eggs per cage. At least three replicate diapause cages and three replicate non-diapause cages are needed for RNA generation. This totals to ca. 1,800-3,000 eggs per experiment.
	<ol>
		<li>Hatch the eggs by submerging egg papers into ca. 500 ml of deionized H2O.</li>
		<li>Add ca. 1 ml food slurry consisting of ground dog food and brine shrimp as previously described24. Cover container with mesh, and keep the mesh in place with a rubber band.</li>
		<li>Place in the LD photoperiod cabinet for ca. 24 hr.</li>
	</ol>
	</li>
	<li>Transfer hatched larvae to 10 x 10 x 2 cm Petri dishes filled with ca. 90 ml deionized H2O.
	<ol>
		<li>Maintain ca. 30 larvae per dish. Transfer larvae to clean dishes every 48&#8211;72 hr, for example every Monday, Wednesday, and Friday (M-W-F)24.</li>
		<li>Feed ca. 1 ml food slurry consisting of ground dog food and brine shrimp in deionized water every M-W-F as previously described24.&#160;</li>
	</ol>
	</li>
	<li>Set up three to four adult cages for each photoperiod treatment, where each cage comprises a biological replicate.
	<ol>
		<li>From opposite sides of 9.5 L buckets, cut out one 10-by-14 cm hole, and another hole with 15 cm diameter. Cover the first with mesh. Cut approximately one foot length of an orthopedic stocking, and glue one end around the inside of the other hole.</li>
		<li>Cut the foot end of the stocking off and knot shut &#8212; open only when access to the interior of the cage is needed. For the cage lid, cut out all of the interior, leaving only the rim, and replace the interior plastic with mesh24.</li>
		<li>Note the photoperiod, replicate number, cage start date, and other information relevant to the experiment with permanent marker on the side of the cage.&#160;</li>
	</ol>
	</li>
	<li>Line the bottom of the adult cages with wet filter paper. Dampen the filter paper with enough deionized H2O to increase local humidity in the cage, but avoid standing water, which can stimulate oviposition on the filter paper24. Check the filter paper daily for drying, re-wet when necessary.</li>
	<li>To produce sufficient eggs for an RNA library, include at least 100 females/9.5 L cage, with no more than 500 mosquitoes per cage.</li>
	<li>Collect pupae M-W-F and place in a small cup of clean H2O at a density of no more than 50 pupae per 25 ml H2O. Transfer the pupae cup to an adult cage. Place cages in the respective photoperiod cabinet &#8212; A. albopictus pupae are photosensitive15.</li>
	<li>Ensure daily that H2O in cups is clean and clear, and remove dead pupae, because build-up of dead pupae can cause mass mortality. Remove H2O cups after all pupae emerge.</li>
	<li>Place organic raisins on the top mesh of the cage to provide sugar for emerged adults. Monitor raisins, and change them every 3-5 days to prevent mold accumulation.</li>
</ol>
2. Maintenance of Adults to Allow Mating and Egg Production
<ol>
	<li>Maintain cages at high humidity (approx. 80%) by lining the cage bottom with a moist filter paper, and provide access to non-moldy raisins, as described above (Sections 1.5 and 1.9).</li>
</ol>
3. Blood-feeding
<ol>
	<li>Prepare to blood-feed females between two to six days after eclosion to ensure that females have been exposed to at least eight unambiguous short days before oviposition commences for nearly 100% diapause eggs14.</li>
	<li>Prepare the Hemotek Membrane Feeding system. Plug the feeding units into the power supply. Adjust the temperature of each unit to 37 &#176;C using the adjustment screw. Use an electronic thermometer and probe to measure the temperature of the feeding unit during calibration.</li>
	<li>Prepare the meal reservoir. Stretch a square of collagen feeding membrane over the aperture of the meal reservoir and secure it with an O-ring. Carefully pull the corners to remove wrinkles; trim the excess membrane with scissors. 
	NOTE: Collagen membrane may not work well for all mosquito species, and it may be necessary to try several types to find the optimal membrane if working with a species other than A. albopictus. Parafilm works well with Culex pipiens.
	<ol>
		<li>If blood is stored frozen, thaw at room temperature for at least 1 hr before using.</li>
		<li>Hold the reservoir so the membrane is facing down, unsupported, and the filling ports are facing up. Use a transfer pipette or syringe to fill the reservoir with approximately 3 ml of whole blood from chickens that has sodium citrate as an anti-coagulant. Seal filling ports with plastic plugs.</li>
	</ol>
	</li>
	<li>Attach the prepared reservoir to the feeder by screwing it onto the stud on the heat transfer plate on the bottom of the feeder. Invert the feeder and place it on top of the cage, membrane side down, so that mosquitoes can feed through the mesh of the cage. Keep the feeder on the cage for approximately 45 min to maximize feeding.</li>
</ol>
4. Stimulate Oviposition
<ol>
	<li>Four to five days post blood meal, equip each cage with a dark colored 50 ml cup lined with unbleached seed germination paper (egg paper) or textured non-bleached paper towel and fill halfway with deionized water9. If more than 250 mosquitoes are in a cage, use two cups. 
	NOTE: For small cages or single-female vials, &#8220;hay infusion&#8221; in oviposition containers may increase oviposition due to the odor of the microbial flora25.</li>
</ol>
5. Collect and Store Eggs
<ol>
	<li>Commence egg collection within 4-5 days of blood feeding because egg production typically peaks approximately five days after blood feeding, and then subsides over the next week.</li>
	<li>Vary egg collection frequency on necessities of the experiment. For general purposes, collect egg papers on an M-W-F schedule. Remove egg papers from each cage and replace with fresh paper. Place recently removed papers in Petri dishes and store in SD photoperiod cabinet to avoid confounding effects of egg storage.</li>
	<li>Allow egg papers to remain wet for 2 days post-oviposition to allow serosal cuticle formation, which increases egg desiccation resistance26.</li>
	<li>Approximately 48 hr post collection, dry eggs in open air. Dry the paper such that it is limp and slightly damp to the touch, but not so wet that the paper is dark from H2O or stimulates hatching of eggs. 
	NOTE: A 6.5&#8221; x 4&#8221; paper may take approximately 3.5 hr to dry. Be cautious not to over-dry egg papers, as this will result in egg desiccation27.</li>
	<li>Reserve additional eggs from both LD and SD photoperiods to assess diapause incidence and interpret the photoperiodic effect (see Measuring Diapause, Section 6).</li>
	<li>For long-term storage, keep egg papers at 21 &#176;C and approximately 80% humidity in Petri dishes. Keep Petri dishes in a Tupperware storage container with a flask of water to maintain local humidity as embryonic development takes four to five days at 21 &#176;C.</li>
</ol>
6. Measure Diapause Incidence
<ol>
	<li>Use additional reserved embryos (see Stimulate Oviposition, Section 4) that are 7&#8211;20 days old to quantify the diapause response.</li>
	<li>Record the number of eggs present on each egg paper.</li>
	<li>Stimulate eggs to hatch by completely submerging individual egg papers in a 90 ml Petri dish with approximately 80 ml deionized H2O. Add approximately 0.25 ml food slurry.</li>
	<li>After 24 hr tally the number of hatched first instar larvae. Place the Petri dish on a black surface to visualize larvae and place a light source on one side of dish. Larvae will move away from the light source, allowing for a clear tally of individual larvae. Remove individual larvae with a pipette while counting to prevent recounting individual larvae.</li>
	<li>Place egg papers in a new Petri dish and re-dry. Re-hatch eggs after ~1 week and again tally eggs hatched using the above method.</li>
	<li>Place egg papers with the remaining un-hatched eggs in new 90 ml Petri dishes with approximately 80 ml bleaching solution28. Ensure that the egg papers are completely submerged in the bleaching solution and leave under a fume hood overnight to avoid the odor of bleach. 
	NOTE: Bleaching solution can be stored for ~1 week at 4 &#176;C, but should otherwise be made fresh.</li>
	<li>Inspect eggs using a light microscope as the bleaching will clear the chorion and allow visualization of embryonated, un-hatched eggs. If the egg is embryonated, the egg will have an off-white color with eyes appearing as two small black dots opposite each other on the dorsal side. Tally the number of un-hatched, embryonated eggs13.</li>
	<li>Determine diapause incidence with the following formula: % diapause = no. embryonated un-hatched eggs / (no. hatched eggs + no. embryonated un-hatched eggs) x 100 13.</li>
</ol>
7. RNA Extraction from Eggs/pharate Larvae

NOTE: Use Trizol in a laminar flow hood.

<ol>
	<li>Brush mosquito eggs containing developing embryos or pharate larvae at distinct development time points from egg papers to glass grinders using a camel-hair brush. Grind the eggs in Trizol (1 ml&#160;per 50&#8211;100 mg of&#160;tissue) until completely pulverized. Use at least 400 eggs per library to yield sufficient RNA.
	<ol>
		<li>Alternatively, snap freeze eggs in liquid nitrogen and stored at -80 &#176;C in microcentrifuge tubes for up to a month before grinding in Trizol.</li>
	</ol>
	</li>
	<li>Perform RNA extraction in Trizol followed by isopropanol precipitation according to manufacturer&#8217;s instructions.</li>
	<li>Treat the bench with RNase decontamination solution or other agents to remove any residual nucleases to avoid RNA degradation.
	<ol>
		<li>Treat the extracted RNA with DNase. According to manufacturer&#8217;s instructions, incubate the RNA samples with DNase for 30 min at 37 &#176;C. Use 1 &#181;l DNase for up to 10 &#181;g of RNA in a 50 &#181;l reaction. Increase the amount of DNase if there are more than 10 &#181;g of RNA in one reaction.</li>
		<li>Inactivate DNase by adding 5 &#181;l suspended DNase inactivation reagent. Incubate 5 min at room temperature, mixing three times during incubation period (gentle vortexing).</li>
		<li>Centrifuge at 10,000 x g for 1.5 min. Transfer the supernatant containing the treated RNA samples to fresh tubes for subsequent steps.</li>
	</ol>
	</li>
	<li>Assess the quality of the total RNA samples by fluorometry. Send the samples to a specialty facility with proper instrument for this task. The facility will perform on-chip gel electrophoresis to determine the sizes of RNA species in the sample, visualized by fluorescent dye instilled in the chip. The results will be returned as an electropherogram.
	<ol>
		<li>Determine the integrity of total RNA samples by the presence or absence of degradation products, as evidenced by the presence of peaks between the 18S and 5S ribosomal RNA peaks on the resulting electropherogram (Figure 1B).</li>
	</ol>
	</li>
</ol>
8. RNA Sequencing
<ol>
	<li>Send total RNA samples with sufficiently high quality (Figure 1A) and quantity (usually &#62;3 &#956;g per library) to a commercial sequencing center for construction of enriched paired-end mRNA libraries and mRNA sequencing, following standard protocols.</li>
	<li>If more than one lane is being used for sequencing a single experiment, split individual libraries into two lanes for sequencing to account for technical variation among lanes during sequencing.</li>
</ol>
9. lllumina Read Cleaning
NOTE: Figure 2 summarizes the bioinformatics portion of this protocol. For a full list of all programs and resources used in the bioinformatics section of this protocol, refer to Table 1. In addition, Supplemental File 1 contains command line examples for each of the following bioinformatics protocol steps.
<ol>
	<li>Use ssaha229 (Table 1) to identify matches of 95% identity or higher to the NCBI UniVec Core database (Table 1), A. albopictus rRNA sequence (GenBank #L22060.1), and sequencing adapters (detailed command-line examples provided in Supplemental File 1). Remove read pairs with matches using Perl or a similar scripting tool, for example by adapting the provided Perl script (Supplemental File 2).</li>
	<li>Clean remaining reads with the SolexaQA package30 (Table 1; Supplemental File 1): trim regions with a phred score equivalent of less than 20 using the default settings of DynamicTrim.pl.
	<ol>
		<li>Remove reads shorter than 25 bp with LengthSort.pl on both forward and reverse reads simultaneously. Evaluate the quality of the cleaned fastq files with FastQC (Table 1) &#8212; in particular, verify that the per-base sequence quality and the per-sequence quality scores are above 20.</li>
	</ol>
	</li>
</ol>
10. Digital Normalization
<ol>
	<li>Perform one round of digital normalization on the cleaned reads using the khmer tool31 (Table 1; Supplemental File 1), specifically normalize-by-median.py (using k-mer size 20, a coverage cut-off of 20, and x = 1e10).</li>
	<li>Alternatively, if a machine with high RAM is available (hundreds of GBs), use Trinity&#39;s normalize_by_kmer_coverage.pl script (Table 1).</li>
</ol>
11. De Novo Transcriptome Assembly
<ol>
	<li>Obtain access to a computer or computer cluster with up to 256 Gb of RAM and 24 cpus, depending on the size of the assembly.</li>
	<li>Use Trinity32 (Table 1; Supplemental File 1) to assemble the digitally normalized read set into contigs. To reduce memory usage, use --min_kmer_cov 2.</li>
</ol>
12. Assembly Evaluation
<ol>
	<li>Run assemblathon_stats.pl from the Assemblathon2 project33 on the Trinity contig output. This script performs basic calculations relevant to evaluating assembly quality, such as number of scaffolds, N50, assembly composition, and more (Table 1; Supplemental File 1).</li>
</ol>
13. Annotation of the Assembled Transcriptome
<ol>
	<li>Perform Blastx (Table 1) of the assembly against a reference protein set; for mosquitoes, Drosophila melanogaster, Anopheles gambiae, Culex pipiens, and Aedes aegypti are suitable references. Specifically, format the reference protein fasta file for blast, followed by blastx (Supplemental File 1).</li>
</ol>
14. Map Reads to the Assembly Using RSEM34 (Table 1)
<ol>
	<li>Create a &#39;transcript-to-gene-map&#39; file, in which the first column contains the reference gene IDs, and the second column the contig IDs. In a spreadsheet editor, swap the first and second columns from the Blastx output, and write these columns to a .txt file. Use the LineBreak program to convert line breaks in the resulting .txt file to Unix format.</li>
	<li>Create a reference dataset from the transcriptome fasta file using the rsem-prepare-reference script, provided in the RSEM package (Supplemental File 1).</li>
	<li>Calculate the expression values separately for each library using the rsem-calculate-expression command, provided in the RSEM package (Supplemental File 1). As reads, use the paired fastq files resulting from the read-cleaning step (step 9.2).
	<ol>
		<li>If RNA from a biological replicate was split into two lanes for sequencing, include both fastq files in the expression calculation to generate a single file.</li>
	</ol>
	</li>
	<li>Convert the expression results from each library to a matrix easily processed by other programs using the provided script rsem-generate-data-matrix, provided in the RSEM package (Supplemental File 1).</li>
</ol>
15. Differential Expression Analysis
<ol>
	<li>Install R and EdgeR (Table 1).</li>
	<li>Use read.delim to load the RSEM results from step 14.3 (Supplemental File 1). If necessary, round the counts to the nearest integer. 
	NOTE: The EdgeR guide recommends limiting the dataset to genes with high enough expression to detect significance.</li>
	<li>To format the data for EdgeR, generate a DGEList object from the loaded data file (Supplemental File 1). Then, normalize the data using TMM normalization (Supplemental File 1). Estimate the common and tagwise dispersions of the data (Supplemental File 1).</li>
	<li>Identify differentially expressed genes with a Benjamini-Hochberg corrected p-value &#60;0.05 (Supplemental File 1). Plot the distribution of log-fold-change vs. abundance (Supplemental File 1).</li>
</ol>

<ol>
	<li>Bilyk, K. T., Cheng, C. H. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Model+of+gene+expression+in+extreme+cold+-+reference+transcriptome+for+the+high-Antarctic+cryopelagic+notothenioid+fish+Pagothenia+borchgrevinki.">Model of gene expression in extreme cold - reference transcriptome for the high-Antarctic cryopelagic notothenioid fish Pagothenia borchgrevinki.</a> BMC Genomics. 14, 634 (2013).</li><li>Chapman, M. A., Hiscock, S. J., Filatov, D. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genomic+divergence+during+speciation+driven+by+adaptation+to+altitude.">Genomic divergence during speciation driven by adaptation to altitude.</a> Mol. Biol. Evol. 30, 2553-2567 (2013).</li><li>Schwarz, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sympatric+ecological+speciation+meets+pyrosequencing:+sampling+the+transcriptome+of+the+apple+maggot+Rhagoletis+pomonella.">Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella.</a> BMC Genomics. 10, 633 (2009).</li><li>Barshis, D. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genomic+basis+for+coral+resilience+to+climate+change.">Genomic basis for coral resilience to climate change.</a> P. Natl Acad. Sci. USA. 110, 1387-1392 (2013).</li><li>Urbanski, J. M., Aruda, A., Armbruster, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+transcriptional+element+of+the+diapause+program+in+the+Asian+tiger+mosquito,+Aedes+albopictus,+identified+by+suppressive+subtractive+hybridization.">A transcriptional element of the diapause program in the Asian tiger mosquito, Aedes albopictus, identified by suppressive subtractive hybridization.</a> J. Insect Physiol. 56, 1147-1154 (2010).</li><li>Huang, Y. H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+duck+genome+and+transcriptome+provide+insight+into+an+avian+influenza+virus+reservoir+species.">The duck genome and transcriptome provide insight into an avian influenza virus reservoir species.</a> Nat. Genet. 45, 776-783 (2013).</li><li>Sessions, O. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Host+cell+transcriptome+profile+during+wild-type+and+attenuated+dengue+virus+infection.">Host cell transcriptome profile during wild-type and attenuated dengue virus infection.</a> PLoS Negl. Trop. Dis. 7 (3), 2107 (2013).</li><li>Poelchau, M. F., Reynolds, J. A., Elsik, C. G., Denlinger, D. L., Armbruster, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Deep+sequencing+reveals+complex+mechanisms+of+diapause+preparation+in+the+invasive+mosquito,+Aedes+albopictus.">Deep sequencing reveals complex mechanisms of diapause preparation in the invasive mosquito, Aedes albopictus.</a> P. R. Soc B. 280, (2013).</li><li>Poelchau, M. F., Reynolds, J. A., Elsik, C. G., Denlinger, D. L., Armbruster, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transcriptome+sequencing+as+a+platform+to+elucidate+molecular+components+of+the+diapause+response+in+Aedes+albopictus.">Transcriptome sequencing as a platform to elucidate molecular components of the diapause response in Aedes albopictus.</a> Physiol. Entomol. 38, 173-181 (2013).</li><li>Poelchau, M. F., Reynolds, J. A., Denlinger, D. L., Elsik, C. G., Armbruster, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+de+novo+transcriptome+of+the+Asian+tiger+mosquito,+Aedes+albopictus,+to+identify+candidate+transcripts+for+diapause+preparation.">A de novo transcriptome of the Asian tiger mosquito, Aedes albopictus, to identify candidate transcripts for diapause preparation.</a> BMC Genomics. 12, 619 (2011).</li><li>Benedict, M. Q., Levine, R. S., Hawley, W. A., Lounibos, L. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Spread+of+the+tiger:+Global+risk+of+invasion+by+the+mosquito+Aedes+albopictus.">Spread of the tiger: Global risk of invasion by the mosquito Aedes albopictus.</a> Vector-Borne Zoonot. 7, 76-85 (2007).</li><li>Lounibos, L. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Invasions+by+insect+vectors+of+human+disease.">Invasions by insect vectors of human disease.</a> Annu. Rev. Entomol. 47, 233-266 (2002).</li><li>Urbanski, J. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Rapid+adaptive+evolution+of+photoperiodic+response+during+invasion+and+range+expansion+across+a+climatic+gradient.">Rapid adaptive evolution of photoperiodic response during invasion and range expansion across a climatic gradient.</a> Am. Nat. 179, 490-500 (2012).</li><li>Lounibos, L. P., Escher, R. L., Lourenco-de-Oliveria, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Asymmetric+evolution+of+photoperiodic+diapause+in+temperate+and+tropical+invasive+populations+of+Aedes+albopictus+(Diptera+Culicidae).">Asymmetric evolution of photoperiodic diapause in temperate and tropical invasive populations of Aedes albopictus (Diptera Culicidae).</a> Ann. Entomol. Soc. Am. 96, 512-518 (2003).</li><li>Mori, A., Oda, T., Wada, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Studies on the egg diapause and overwintering of Aedes albopictus in Nagasaki. Trop. Med. 23, 79-90 (1981).</li><li>Pumpuni, C. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Factors influencing photoperiodic control of egg diapause in Aedes albopictus [dissertation]. , (1989).</li><li>Wang, R. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Observations+on+the+influence+of+photoperiod+on+egg+diapause+in+Aedes+albopictus+Skuse.">Observations on the influence of photoperiod on egg diapause in Aedes albopictus Skuse.</a> Acta Entomol. Sinica. 15, 75-77 (1966).</li><li>Sota, T., Mogi, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Survival-time+and+resistance+to+desiccation+of+diapause+and+non-diapause+eggs+of+temperate+Aedes+(Stegomyia)+mosquitoes.">Survival-time and resistance to desiccation of diapause and non-diapause eggs of temperate Aedes (Stegomyia) mosquitoes.</a> Entomol. Exp. Appl. 63, 155-161 (1992).</li><li>Reynolds, J. A., Poelchau, M. F., Rahman, Z., Armbruster, P. A., Denlinger, D. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transcript+profiling+reveals+mechanisms+for+lipid+conservation+during+diapause+in+the+mosquito,+Aedes+albopictus.">Transcript profiling reveals mechanisms for lipid conservation during diapause in the mosquito, Aedes albopictus.</a> J. Insect Physiol. 58, 966-973 (2012).</li><li>Andrewartha, H. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Diapause+in+relation+to+the+ecology+of+insects.">Diapause in relation to the ecology of insects.</a> Biol. Rev. 27, 50-107 (1952).</li><li>Danks, H. V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Insect+Dormancy:+An+Ecological+Perspective.">Insect Dormancy: An Ecological Perspective.</a> Biological Survey of Canada (Terrestrial Arthropods). , (1987).</li><li>Denlinger, D. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Regulation+of+diapause.">Regulation of diapause.</a> Ann. Rev. Entomol. 47, 93-122 (2002).</li><li>Rev Entomol, A. n. n. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> 59, 93-122 (2014).</li><li>Armbruster, P. A., Conn, J. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Geographic+variation+of+larval+growth+in+North+American+Aedes+albopictus+(Diptera).">Geographic variation of larval growth in North American Aedes albopictus (Diptera).</a> Culicidae). Ann. Entomol. Soc. Am. 99, 1234-1243 (2006).</li><li>Reiter, P., Amador, M. A., Colon, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Enhancement+of+the+CDC+ovitrap+with+hay+infusions+for+daily+monitoring+of+Aedes+aegypti+populations.">Enhancement of the CDC ovitrap with hay infusions for daily monitoring of Aedes aegypti populations.</a> J. Am. Mosquito Contr. Association. 7 (1), 52 (1991).</li><li>Rezende, G. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Embryonic+desiccation+resistance+in+Aedes+aegypti:+presumptive+role+of+the+chitinized+serosal+cuticle.">Embryonic desiccation resistance in Aedes aegypti: presumptive role of the chitinized serosal cuticle.</a> BMC Dev. Biol. 8, 182 (2008).</li><li>Munstermann, L., Crampton, J., Beard, C., C, L. o. u. i. s. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Care+and+maintenance+of+Aedes+mosquito+colonies.">Care and maintenance of Aedes mosquito colonies.</a> The Molecular Biology of Insect Disease Vectors. , 13-20 (1997).</li><li>Trpis, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+new+bleaching+and+decalcifying+method+for+general+use+in+zoology.">A new bleaching and decalcifying method for general use in zoology.</a> Can. J. Zoolog. 48, 892-893 (1970).</li><li>Ning, Z. M., Cox, A. J., Mullikin, J. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SSAHA:+A+fast+search+method+for+large+DNA+databases.">SSAHA: A fast search method for large DNA databases.</a> Genome Res. 11 (10), 1725-1729 (2001).</li><li>Cox, M. P., Peterson, D. A., Biggs, P. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SolexaQA:+At-a-glance+quality+assessment+of+Illumina+second-generation+sequencing+data.">SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data.</a> BMC Bioinformatics. 11, 485 (2010).</li><li>Grabherr, M. G., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Full-length+transcriptome+assembly+from+RNA-Seq+data+without+a+reference+genome.">Full-length transcriptome assembly from RNA-Seq data without a reference genome.</a> Nat. Biotechnol. 29 (7), 644-652 (2011).</li><li>Bradnam, K. R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Assemblathon+2:+evaluating+de+novo+methods+of+genome+assembly+in+three+vertebrate+species.">Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species.</a> GigaScience. 2, 10 (2013).</li><li>Li, B., Dewey, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RSEM:+accurate+transcript+quantification+from+RNA-Seq+data+with+or+without+a+reference+genome.">RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome.</a> BMC Bioinformatics. 12, 323 (2011).</li><li>White, B. N., De Lucca, F. L., Turner, R. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Preparation+and+analysis+of+RNA.">Preparation and analysis of RNA.</a> Analytical Biochemistry of Insects. , (1977).</li><li>Hawley, W. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+biology+of+Aedes+albopictus.">The biology of Aedes albopictus.</a> J. Am. Mosq. Contr. Assoc. 4, 1-39 (1988).</li><li>Dowling, Z., Ladeau, S. L., Armbruster, P., Biehler, D., Leisnham, P. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Socioeconomic+status+affect+mosquito+(Diptera:Culicidae)+larval+habitat+type+availability+and+infestation+level.">Socioeconomic status affect mosquito (Diptera:Culicidae) larval habitat type availability and infestation level.</a> J. Med Entomol. 50, 764-772 (2013).</li><li>Benedict, M. Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chapter+2,4,10.+Bloodfeeding:+Membrane+apparatuses+and+animals.">Chapter 2,4,10. Bloodfeeding: Membrane apparatuses and animals.</a> Methods in Anopheles Research. , (2010).</li><li>Das, S., Garver, S., Ramirez, J. R., Xi, Z., Dimopolous, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Protocol+for+dengue+infections+in+mosquitoes+(A.+aegypti)+and+infection+phenotype+determination.">Protocol for dengue infections in mosquitoes (A. aegypti) and infection phenotype determination.</a> J. Vis. Exp. 5 (220), (2007).</li><li>Goff, S. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+iPlant+collaborative:+cyberinfrastructure+for+plant+biology.">The iPlant collaborative: cyberinfrastructure for plant biology.</a> Plant Sci. 2 (34), (2011).</li><li>Robinson, M. D., McCarthy, D. J., Smyth, G. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=edgeR:+a+Bioconductor+package+for+differential+expression+analysis+of+digital+gene+expression+data.">edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.</a> Bioinformatics. 26, 139-140 (2010).</li><li>Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., Knight, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=UCHIME+improves+sensitivity+and+speed+of+chimera+detection.">UCHIME improves sensitivity and speed of chimera detection.</a> Bioinformatics. 27 (16), 2194-2220 (2011).</li><li>Goecks, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Galaxy:+a+comprehensive+approach+for+supporting+accessible,+reproducible,+and+transparent+computational+research+in+the+life+sciences.">Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences.</a> Genome Biol. 11 (8), 86 (2010).</li></ol>

The authors have nothing to disclose.

This protocol presents methods to discover differentially expressed genes due to photoperiodically induced diapause in A. albopictus. The protocol is significant in that it uniquely combines mosquito rearing and bioinformatics techniques to make all experimental aspects of a molecular physiology program accessible to novice users &#8212; in particular for those focusing on the photoperiodic diapause response. Existing methods, to our knowledge, do not provide as much detail in the rearing protocol &#8212; which ...

Rapid advances in next-generation sequencing (NGS) technologies are providing exciting opportunities to probe the molecular underpinnings of a wide range of genetically complex ecological adaptations in a broad diversity of non-model organisms1&#8211;3. This approach is extremely powerful because it establishes a basis for population and functional genomics studies of organisms with an especially interesting and/or well-described ecology or evolutionary history, as well as organisms of practical concern, such as agricultural pests and disease vectors. Thus, NGS technologies are leading to rapid advances in the fields of ecology and have the potential to address problems such as understanding the mechanistic bases of biological responses to rapid contemporary climate change4, the spread of invasive species5, and host-pathogen interactions6,7.
The extraordinary potential of NGS technologies for addressing basic and applied questions in ecology and evolutionary biology is in part due to the fact that these approaches can be applied to any organism at a moderate cost that is feasible for most research laboratories. Furthermore, these approaches provide genome-wide information without the requirement of a priori genetic resources such as a microarray chip or complete genome sequence. Nevertheless, to maximize the productivity of NGS experiments requires careful consideration of experimental design including issues such as the developmental timing and tissue-specificity of RNA sampling. Furthermore, the technical skills required to analyze the massive amounts of data produced by these experiments, often up to several hundred million DNA sequence reads, has been a particular challenge and has limited the widespread implementation of NGS approaches.
Recent RNA-Seq studies on the transcriptional bases of diapause in the invasive and medically important mosquito Aedes albopictus provide a useful example of some of the experimental protocols that can be employed to successfully apply NGS technology to studying the molecular basis of a complex ecological adaptation in a non-model organism8&#8211;10. A. albopictus is a highly invasive species that is native to Asia but has recently invaded North America, South America, Europe, and Africa11,12. Like many temperate insects, temperate populations of A. albopictus survive through winter by entering a type of dormancy referred to as photoperiodic diapause. In A. albopictus, exposure of pupal and adult females to short (autumnal) day lengths leads to the production of diapause eggs in which embryological development is completed, but the pharate larva inside the chorion of the egg enters a developmental arrest that renders the egg refractory to hatching stimulus15&#8211;17. Diapause eggs are more desiccation resistant5,18 and contain more total lipids19 than non-diapause eggs. Photoperiodic diapause in A. albopictus is thus a maternally controlled, adaptive phenotypic plasticity that is essential for surviving the harsh conditions of winter in temperate environments. Despite the well-understood ecological significance of photoperiodic diapause in a wide range of insects20,21, the molecular basis of this crucial adaptation is not well characterized in any insect22. In organisms such as A. albopictus that undergo an embryonic diapause at the pharate larval stage, it remains a particularly compelling challenge to understand how the photoperiodic signal received by the mother is passed to the offspring and persists through the course of embryonic development to cause arrest at the pharate larval stage.
This protocol describes mosquito rearing, experimental design and bioinformatics analyses for NGS experiments (transcriptome sequencing) performed to elucidate transcriptional components of photoperiodic diapause in A. albopictus. This protocol can be used for additional studies of diapause in A. albopictus, can be adapted to investigate diapause in other closely related species such as other aedine mosquitoes that undergo egg diapause23, and is also more generally relevant to employing NGS approaches to study the transcriptional bases of any complex adaptation in any insect.

<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Incubator - Model 818</td>
			<td style="width:116px;">Thermo-Scientific</td>
			<td align="right" style="width:109px;">3751</td>
			<td style="width:351px;">120V</td>
		</tr>
		<tr height="83">
			<td height="83" style="height:83px;width:209px;">Controlled environment room</td>
			<td>Thermax Scientific</td>
			<td style="width:109px;">N/A</td>
			<td style="width:351px;">Walk-in controlled environment room built to custom specifications by Thermax Scientific Products. A larger alternative to an incubator. http://thermmax.com/</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:209px;">Cool Fluorescent bulb</td>
			<td style="width:116px;">Philips</td>
			<td align="right">392183</td>
			<td style="width:351px;">4 Watt</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:209px;">Petri Dish 100mm x 20mm</td>
			<td style="width:116px;">Fisher</td>
			<td>08-772-E</td>
			<td style="width:351px;"></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:209px;">Filter Paper 20.5cm</td>
			<td style="width:116px;">Fisher</td>
			<td>09-803-6J</td>
			<td style="width:351px;"></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">9.5L Bucket</td>
			<td style="width:116px;">Plastican</td>
			<td style="width:109px;">Bway Products</td>
			<td style="width:351px;">http://www.bwayproducts.com/sites/portal/plastic-products/plastic-open-head-pails/117</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Utility Fabric-Mosquito Netting White</td>
			<td style="width:116px;">Joann</td>
			<td>10173292</td>
			<td style="width:351px;">http://www.joann.com/utility-fabric-mosquito-netting-white/10173292.html</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:209px;">Orthopedic stockings</td>
			<td>Albahealth</td>
			<td style="width:109px;">23650-040</td>
			<td style="width:351px;">product no. 081420</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Organic Raisins</td>
			<td style="width:116px;">Newman&#39;s Own</td>
			<td style="width:109px;">UPC: 884284040255</td>
			<td style="width:351px;"></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Oviposition cups (brown)</td>
			<td>Fisher Scientific</td>
			<td>03-007-52</td>
			<td style="width:351px;">The product is actually an amber 125 mL bottle that we saw the top off of.</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Recycled Paper Towels</td>
			<td style="width:116px;">Seventh Generation</td>
			<td style="width:109px;">30BPT120</td>
			<td style="width:351px;"></td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Modular Mates Square Tupperware Set</td>
			<td style="width:116px;">Tupperware</td>
			<td style="width:109px;"></td>
			<td style="width:351px;">http://order.tupperware.com/pls/htprod_www/coe$www.add_items</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Glass Grinder</td>
			<td style="width:116px;">Corning Incorporated</td>
			<td style="width:109px;">7727-2</td>
			<td style="width:351px;">These Tenbroeck tissue grinders break the eggs and release RNA into the TRI Reagent.</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">TRI Reagent</td>
			<td style="width:116px;">Sigma Aldrich</td>
			<td style="width:109px;">T9424</td>
			<td style="width:351px;">Apply 1ml TRI Reagent per 50-100mg of tissue. Caution - this reagent is toxic.</td>
		</tr>
		<tr height="62">
			<td height="62" style="height:62px;width:209px;">TURBO DNA-free</td>
			<td style="width:116px;">Ambion/Life Technologies</td>
			<td style="width:109px;">AM1907</td>
			<td style="width:351px;">This kit generates greater yield than traditional DNase treatment followed by phenol/chloroform cleanup, and it is simpler to use.</td>
		</tr>
		<tr height="104">
			<td height="104" style="height:104px;width:209px;">RNaseZap</td>
			<td style="width:116px;">Ambion/Life Technologies</td>
			<td style="width:109px;">AM9782</td>
			<td style="width:351px;">Apply liberally on the bench surfaces and any equipment that might be in contact with the RNA samples. The solution is slightly alkaline/corrosive, can cause irritation and is harmful when swallowed.</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">2100 Bioanalyzer</td>
			<td style="width:116px;">Agilent Technologies</td>
			<td style="width:109px;">G2939AA</td>
			<td style="width:351px;">Place up to 12 RNA samples on one chip.</td>
		</tr>
		<tr height="187">
			<td height="187" style="height:187px;width:209px;">Hemotek Membrane Feeder</td>
			<td style="width:116px;">Hemotek&#160;</td>
			<td style="width:109px;">5W1</td>
			<td style="width:351px;">This system&#160; provides 5 feeding stations that can be used simultaneously. Includes PS5 Power Unit and Power cord; 5 FUI Feeders + Meal Reservoirs and O-rings; Plastic Plugs, Hemotek collagen feeding membrane; Temperature setting tool; and Plug extracting tool. The company&#39;s mailing address is: Hemotek Ltd; Unit 5 Union Court; Alan Ramsbottom Way; Great Harwood; Lancashire, UK; BB6 7FD; tel: +44 1254 889 307.</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:209px;">Digital Thermometer and Probe</td>
			<td style="width:116px;">Hemotek&#160;</td>
			<td style="width:109px;">MT3KFU</td>
			<td style="width:351px;">MicroT3 thermometer and KFU probe. This is used to set the temperature of each FUI feeding unit.</td>
		</tr>
		<tr height="83">
			<td height="83" style="height:83px;width:209px;">Chicken Whole Blood, non-sterile with Sodium Citrate</td>
			<td style="width:116px;">Pel-Freez Biologicals</td>
			<td>33130-1</td>
			<td style="width:351px;">The 500 ml of blood were frozen and stored in 20 ml aliquots at -80 degrees C for up to 1 year.&#160; Thaw blood at room temperature for at least 1 h before using.</td>
		</tr>
	</tbody>
</table>

an experimental and bioinformatics protocol for rna-seq analyses of photoperiodic diapause in the asian tiger mosquito, aedes albopictus

Photoperiodic diapause is an important adaptation that allows individuals to escape harsh seasonal environments via a series of physiological changes, most notably developmental arrest and reduced metabolism. Global gene expression profiling via RNA-Seq can provide important insights into the transcriptional mechanisms of photoperiodic diapause. The Asian tiger mosquito, Aedes albopictus, is an outstanding organism for studying the transcriptional bases of diapause due to its ease of rearing, easily induced diapause, and the genomic resources available. This manuscript presents a general experimental workflow for identifying diapause-induced transcriptional differences in A. albopictus. Rearing techniques, conditions necessary to induce diapause and non-diapause development, methods to estimate percent diapause in a population, and RNA extraction and integrity assessment for mosquitoes are documented. A workflow to process RNA-Seq data from Illumina sequencers culminates in a list of differentially expressed genes. The representative results demonstrate that this protocol can be used to effectively identify genes differentially regulated at the transcriptional level in A. albopictus due to photoperiodic differences. With modest adjustments, this workflow can be readily adapted to study the transcriptional bases of diapause or other important life history traits in other mosquitoes.

Fluorometry of two representative RNA samples showed two bands at approximately 2,000 nt (Figure 1A, B). The insect 28S ribosomal RNA is comprised of two polynucleotide chains held together by hydrogen bonds, which are easily disrupted by brief heating or agents that disrupt hydrogen bonds35. The resulting two components are approximately the same size as the 18S ribosomal RNA. The second RNA sample showed high levels of degradation (Figure 1B).
Pho...

Watch this Scientific Journal Video about An Experimental and Bioinformatics Protocol for RNA-seq Analyses of Photoperiodic Diapause in the Asian Tiger Mosquito, Aedes albopictus at JoVE.com

An Experimental and Bioinformatics Protocol for RNA-seq Analyses of Photoperiodic Diapause in the Asian Tiger Mosquito, Aedes albopictus

RNA-Seq analyses are becoming increasingly important for identifying the molecular underpinnings of adaptive traits in non-model organisms. Here, a protocol to identify differentially expressed genes between diapause and non-diapause Aedes albopictus mosquitoes is described, from mosquito rearing, to RNA sequencing and bioinformatics analyses of RNA-Seq data.

An Experimental and Bioinformatics Protocol for RNA-seq Analyses of Photoperiodic Diapause in the Asian Tiger Mosquito, Aedes albopictus

Photoperiodic diapause is an important adaptation that allows individuals to escape harsh seasonal environments via a series of physiological changes, ...

an-experimental-bioinformatics-protocol-for-rna-seq-analyses

Microsoft Bing

Research

JoVE Journal

Biology

13.2K Views.  Georgetown University. RNA-Seq analyses are becoming increasingly important for identifying the molecular underpinnings of adaptive traits in non-model organisms. Here, a protocol to identify differentially expressed genes between diapause and non-diapause Aedes albopictus mosquitoes is described, from mosquito rearing, to RNA sequencing and bioinformatics analyses of RNA-Seq data. 

Video: An Experimental and Bioinformatics Protocol for RNA-seq Analyses of Photoperiodic Diapause in the Asian Tiger Mosquito, Aedes albopictus

Protocol for Mosquito Rearing (A. gambiae)

This protocol describes mosquito rearing in the insectary. The insectary rooms are maintained at 28&deg;C and ~80% humidity, with a 12 hr. day/night cycle. For this procedure, you'll need mosquito cages, 10% sterile sucrose solution, paper towels, beaker, whatman filter paper, glass feeders, human blood and serum, water bath, parafilm, distilled water, clean plastic trays, mosquito food (described below), mosquito net to cover the trays, vacuum, and a collection chamber to collect adults.

Protocol for Mosquito Rearing A. gambiae

This video illustrates the general techniques used to rear Anopheles gambiae in the laboratory. The methods for caring for laboratory mosquitoes are demonstrated through all stages of the organism's life cycle from larvae to pupae to blood-feeding adults.

This protocol describes mosquito rearing in the insectary. The insectary rooms are maintained at 28&deg;C and ~80% humidity, with a 12 hr. day/night ...

Video Bioinformatics Analysis of Human Embryonic Stem Cell Colony Growth

Because video data are complex and are comprised of many images, mining information from video material is difficult to do without the aid of computer software.  Video bioinformatics is a powerful quantitative approach for extracting spatio-temporal data from video images using computer software to perform dating mining and analysis. In this article, we introduce a video bioinformatics method for quantifying the growth of human embryonic stem cells (hESC) by analyzing time-lapse videos collected in a Nikon BioStation CT incubator equipped with a camera for video imaging. In our experiments, hESC colonies that were attached to Matrigel were filmed for 48 hours in the BioStation CT.  To determine the rate of growth of these colonies, recipes were developed using CL-Quant software which enables users to extract various types of data from video images.  To accurately evaluate colony growth, three recipes were created. The first segmented the image into the colony and background, the second enhanced the image to define colonies throughout the video sequence accurately, and the third measured the number of pixels in the colony over time.  The three recipes were run in sequence on video data collected in a BioStation CT to analyze the rate of growth of individual hESC colonies over 48 hours. To verify the truthfulness of the CL-Quant recipes, the same data were analyzed manually using Adobe Photoshop software. When the data obtained using the CL-Quant recipes and Photoshop were compared, results were virtually identical, indicating the CL-Quant recipes were truthful. The method described here could be applied to any video data to measure growth rates of hESC or other cells that grow in colonies. In addition, other video bioinformatics recipes can be developed in the future for other cell processes such as migration, apoptosis, and cell adhesion.  

Video bioinformatics is the automated processing, analysis, understanding, and data mining of biological spatio-temporal data extracted from microscopic videos. The purpose of this article is to demonstrate a method for measuring human embryonic stem cell colony growth using a video bioinformatics method.

Because video data are complex and are comprised of many images, mining information from video material is difficult to do without the aid of computer ...

A Protocol for Computer-Based Protein Structure and Function Prediction

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the understanding of their biological role. Although experimental methods can provide detailed information for a small fraction of these proteins, computational modeling is needed for the majority of protein molecules which are experimentally uncharacterized. The I-TASSER server is an on-line workbench for high-resolution modeling of protein structure and function. Given a protein sequence, a typical output from the I-TASSER server includes secondary structure prediction, predicted solvent accessibility of each residue, homologous template proteins detected by threading and structure alignments, up to five full-length tertiary structural models, and structure-based functional annotations for enzyme classification, Gene Ontology terms and protein-ligand binding sites. All the predictions are tagged with a confidence score which tells how accurate the predictions are without knowing the experimental data. To facilitate the special requests of end users, the server provides channels to accept user-specified inter-residue distance and contact maps to interactively change the I-TASSER modeling; it also allows users to specify any proteins as template, or to exclude any template proteins during the structure assembly simulations. The structural information could be collected by the users based on experimental evidences or biological insights with the purpose of improving the quality of I-TASSER predictions. The server was evaluated as the best programs for protein structure and function predictions in the recent community-wide CASP experiments. There are currently &gt;20,000 registered scientists from over 100 countries who are using the on-line I-TASSER server.

Guidelines for computer based structural and functional characterization of protein using the I-TASSER pipeline is described. Starting from query protein sequence, 3D models are generated using multiple threading alignments and iterative structural assembly simulations. Functional inferences are thereafter drawn based on matches to proteins with known structure and functions.

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the ...

RNA-seq Analysis of Transcriptomes in Thrombin-treated and Control Human Pulmonary Microvascular Endothelial Cells

The characterization of gene expression in cells via measurement of mRNA levels is a useful tool in determining how the transcriptional machinery of the cell is affected by external signals (e.g. drug treatment), or how cells differ between a healthy state and a diseased state. With the advent and continuous refinement of next-generation DNA sequencing technology, RNA-sequencing (RNA-seq) has become an increasingly popular method of transcriptome analysis to catalog all species of transcripts, to determine the transcriptional structure of all expressed genes and to quantify the changing expression levels of the total set of transcripts in a given cell, tissue or organism1,2 . RNA-seq is gradually replacing DNA microarrays as a preferred method for transcriptome analysis because it has the advantages of profiling a complete transcriptome, providing a digital type datum (copy number of any transcript) and not relying on any known genomic sequence3.
Here, we present a complete and detailed protocol to apply RNA-seq to profile transcriptomes in human pulmonary microvascular endothelial cells with or without thrombin treatment. This protocol is based on our recent published study entitled "RNA-seq Reveals Novel Transcriptome of Genes and Their Isoforms in Human Pulmonary Microvascular Endothelial Cells Treated with Thrombin,"4 in which we successfully performed the first complete transcriptome analysis of human pulmonary microvascular endothelial cells treated with thrombin using RNA-seq. It yielded unprecedented resources for further experimentation to gain insights into molecular mechanisms underlying thrombin-mediated endothelial dysfunction in the pathogenesis of inflammatory conditions, cancer, diabetes, and coronary heart disease, and provides potential new leads for therapeutic targets to those diseases.
The descriptive text of this protocol is divided into four parts. The first part describes the treatment of human pulmonary microvascular endothelial cells with thrombin and RNA isolation, quality analysis and quantification. The second part describes library construction and sequencing. The third part describes the data analysis. The fourth part describes an RT-PCR validation assay. Representative results of several key steps are displayed. Useful tips or precautions to boost success in key steps are provided in the Discussion section. Although this protocol uses human pulmonary microvascular endothelial cells treated with thrombin, it can be generalized to profile transcriptomes in both mammalian and non-mammalian cells and in tissues treated with different stimuli or inhibitors, or to compare transcriptomes in cells or tissues between a healthy state and a disease state.

This protocol presents a complete and detailed procedure to apply RNA-seq, a powerful next-generation DNA sequencing technology, to profile transcriptomes in human pulmonary microvascular endothelial cells with or without thrombin treatment. This protocol is generalizable to various cells or tissues affected by different reagents or disease states.

The characterization of gene expression in cells via measurement of mRNA levels is a useful tool in determining how the transcriptional machinery of the ...

Depletion of Ribosomal RNA for Mosquito Gut Metagenomic RNA-seq

The mosquito gut accommodates dynamic microbial communities across different stages of the insect's life cycle. Characterization of the genetic capacity and functionality of the gut community will provide insight into the effects of gut microbiota on mosquito life traits. Metagenomic RNA-Seq has become an important tool to analyze transcriptomes from various microbes present in a microbial community. Messenger RNA usually comprises only 1-3% of total RNA, while rRNA constitutes approximately 90%. It is challenging to enrich messenger RNA from a metagenomic microbial RNA sample because most prokaryotic mRNA species lack stable poly(A) tails. This prevents oligo d(T) mediated mRNA isolation. Here, we describe a protocol that employs sample derived rRNA capture probes to remove rRNA from a metagenomic total RNA sample. To begin, both mosquito and microbial small and large subunit rRNA fragments are amplified from a metagenomic community DNA sample. Then, the community specific biotinylated antisense ribosomal RNA probes are synthesized in vitro using T7 RNA polymerase. The biotinylated rRNA probes are hybridized to the total RNA. The hybrids are captured by streptavidin-coated beads and removed from the total RNA. This subtraction-based protocol efficiently removes both mosquito and microbial rRNA from the total RNA sample. The mRNA enriched sample is further processed for RNA amplification and RNA-Seq.

A ribosomal RNA (rRNA) depletion protocol was developed to enrich messenger RNA (mRNA) for RNA-seq of the mosquito gut metatranscriptome. Sample specific rRNA probes, which were used to remove rRNA via subtraction, were created from the mosquito and its gut microbes. Performance of the protocol can result in the removal of approximately 90-99% of rRNA.

The mosquito gut accommodates dynamic microbial communities across different stages of the insect's life cycle. Characterization of the genetic capacity ...

Physiological Recordings and RNA Sequencing of the Gustatory Appendages of the Yellow-fever Mosquito Aedes aegypti

Electrophysiological recording of action potentials from sensory neurons of mosquitoes provides investigators a glimpse into the chemical perception of these disease vectors. We have recently identified a bitter sensing neuron in the labellum of female Aedes aegypti that responds to DEET and other repellents, as well as bitter quinine, through direct electrophysiological investigation. These gustatory receptor neuron responses prompted our sequencing of total mRNA from both male and female labella and tarsi samples to elucidate the putative chemoreception genes expressed in these contact chemoreception tissues. Samples of tarsi were divided into pro-, meso- and metathoracic subtypes for both sexes. We then validated our dataset by conducting qRT-PCR on the same tissue samples and used statistical methods to compare results between the two methods. Studies addressing molecular function may now target specific genes to determine those involved in repellent perception by mosquitoes. These receptor pathways may be used to screen novel repellents towards disruption of host-seeking behavior to curb the spread of harmful viruses.

Physiological Recordings and RNA Sequencing of the Gustatory Appendages of the Yellow-fever Mosquito Aedes aegypti

Using two methods to estimate gene expression in the major gustatory appendages of Aedes aegypti, we have identified the set of genes putatively underlying the neuronal responses to bitter and repulsive compounds, as determined by electrophysiological examination.

Electrophysiological recording of action potentials from sensory neurons of mosquitoes provides investigators a glimpse into the chemical perception of ...

A Protocol for Laboratory Housing of Turquoise Killifish (Nothobranchius furzeri)

The development of husbandry practices in non-model laboratory fish used for experimental purposes has greatly benefited from the establishment of reference fish model systems, such as zebrafish and medaka. In recent years, an emerging fish &#8211; the turquoise killifish (Nothobranchius furzeri) &#8211; has been adopted by a growing number of research groups in the fields of biology of aging and ecology. With a captive life span of 4 - 8 months, this species is the shortest-lived vertebrate raised in captivity and allows the scientific community to test &#8211; in a short time &#8211; experimental interventions that can lead to alterations of the aging rate and life expectancy. Given the unique biology of this species, characterized by embryonic diapause, explosive sexual maturation, marked morphological and behavioral sexual dimorphism - and their relatively short adult life span -&#160;ad hoc husbandry practices are in urgent demand. This protocol reports a set of key husbandry measures that allow optimal turquoise killifish laboratory care, enabling the scientific community to adopt this species as a powerful laboratory animal model.

A Protocol for Laboratory Housing of Turquoise Killifish Nothobranchius furzeri

Laboratory housing of turquoise killifish can be scaled up&#160;to house and efficiently raise thousands of individual fish in a centralized water filtration system, employing the same infrastructure used for standard zebrafish facilities. Here we detail a list of standardized procedures that allow efficient killifish maintenance.

The development of husbandry practices in non-model laboratory fish used for experimental purposes has greatly benefited from the establishment of ...

Experimental Protocol for Using Drosophila As an Invertebrate Model System for Toxicity Testing in the Laboratory

Emergent properties and external factors (population-level and ecosystem-level interactions, in particular) play important roles in mediating ecologically-important endpoints, though they are rarely considered in toxicological studies. D. melanogaster is emerging as a toxicology model for the behavioral, neurological, and genetic impacts of toxicants, to name a few. More importantly, species in the genus Drosophila can be utilized as a model system for an integrative framework approach to incorporate emergent properties and answer ecologically-relevant questions in&#160;toxicology research. The aim of this paper is to provide a protocol for exposing species in the genus Drosophila to pollutants to be used as a model system for a range of phenotypic outputs and ecologically-relevant questions. More specifically, this protocol can be used to 1) link multiple biological levels of organization and understand the impact of toxicants on both individual- and population-level fitness; 2) test the impact of toxicants at different stages of developmental exposure; 3) test multigenerational and evolutionary implications of pollutants; and 4) test multiple contaminants and stressors simultaneously.

Experimental Protocol for Using Drosophila As an Invertebrate Model System for Toxicity Testing in the Laboratory

In this paper, we provide a detailed protocol for exposing species in the genus Drosophila to pollutants with the goal of studying the impact of exposure on a range of phenotypic outputs at different developmental stages and for more than one generation.

Emergent properties and external factors (population-level and ecosystem-level interactions, in particular) play important roles in mediating ...

A Bioinformatics Pipeline for Investigating Molecular Evolution and Gene Expression using RNA-seq

Distilling and reporting large datasets, such as whole genome or transcriptome data, is often a daunting task. One way to break down results is to focus on one or more gene families that are significant to the organism and study. In this protocol, we outline bioinformatic steps to generate a phylogeny and to quantify the expression of genes of interest. Phylogenetic trees can give insight into how genes are evolving within and between species as well as reveal orthology. These results can be enhanced using RNA-seq data to compare the expression of these genes in different individuals or tissues. Studies of molecular evolution and expression can reveal modes of evolution and conservation of gene function between species. The characterization of a gene family can serve as a springboard for future studies and can highlight an important gene family in a new genome or transcriptome paper.

The purpose of this protocol is to investigate the evolution and expression of candidate genes using RNA sequencing data.

Distilling and reporting large datasets, such as whole genome or transcriptome data, is often a daunting task. One way to break down results is to focus ...

A Step-by-Step Guide to Mosquito Electroantennography

Female mosquitoes are the deadliest animals on earth, claiming the lives of more than 1 million people every year due to pathogens they transmit when acquiring a blood-meal. To locate a host to feed on, mosquitoes rely on a wide range of sensory cues, including visual, mechanical, thermal, and olfactory. The study details a technique, electroantennography (EAG), that allows researchers to assess whether the mosquitoes can detect individual chemicals and blends of chemicals in a concentration-dependent manner. When coupled with gas-chromatography (GC-EAG), this technique allows to expose the antennae to a full headspace/complex mixture and determines which chemicals present in the sample of interest, the mosquito can detect. This is applicable to host body odors as well as plant floral bouquets or other ecologically relevant odors (e.g., oviposition sites odorants). Here, we described a protocol that permits long durations of preparation responsiveness time and is applicable to both female and male mosquitoes from multiple genera, including Aedes, Culex, Anopheles, and Toxorhynchites mosquitoes. As olfaction plays a major part in mosquito-host interactions and mosquito biology in general, EAGs and GC-EAG can reveal compounds of interest for the development of new disease vector control strategies (e.g., baits). Complemented with behavioral assays, the valence (e.g., attractant, repellent) of each chemical can be determined.

The present article details a step-by-step protocol for successful and low noise electroantennograms in several genera of mosquitoes, including both females and males.

Female mosquitoes are the deadliest animals on earth, claiming the lives of more than 1 million people every year due to pathogens they transmit when ...