Eukaryotic 18S rDNA PCR amplification and high-throughput sequencing of antFOCE Biofilms

This metadata record contains an Excel spreadsheet with Operational Taxonomic Units (OTUs) gained from Eukaryotic 18S rDNA PCR amplification and high-throughput sequencing of samples from Biofilm slides deployed as part of the antFOCE experiment in the austral summer of 2014/15 at Casey station, East Antarctica. 
Refer to antFOCE report section 4.5.3 for deployment, sampling and analysis details.
https://data.aad.gov.au/metadata/records/AAS_4127_antFOCE_Project4127
Sampling design
2 trays of 8 horizontal standard glass microscope slides (72 x 25 mm) per chamber. Four of the glass slides were scored with a diamond pencil approximately 18 mm from the right hand end of the slide and deployed scored side up. The remaining four slides were unmodified. Slides were sampled at:
- Tmid - one tray per chamber / open plot. The sampled try was repopulated with fresh slides and redeployed
- Tend – 2 slides trays per chamber / open plot.
Sampling procedure
After 31 days deployment, 1 slide tray per chamber / open plot was sampled. At Tend both trays in each chamber / open plot were sampled. To minimize disturbance while being raised to the surface, each tray was removed from the tray holder by divers and placed in a seawater filled container with a lid. On the surface, slides were removed from the tray using ethanol sterilized forceps. The four unscoured slides per chamber / open plot were placed in a plastic microscope slide holder with a sealable lid. The scoured slides were placed individually in 70 ml plastic sample jars.  
Lab procedure - Casey
The slide holder (4 unscoured slides) from each chamber / open plot was frozen at -20C immediately upon return to the lab. The scoured slides were preserved in sea water containing 1% final concentration glutaraldehyde in separate jars.
Preservation Issue: Scoured slides were not refrigerated, either at Casey, during RTA or in Kingston before the 26th Nov 2015, when they were transferred to the 4C Cold Store. 
antFOCE Background 
The antFOCE experimental system was deployed in O’Brien Bay, approximately 5 kilometres south of Casey station, East Antarctica, in the austral summer of 2014/15. Surface and sub-surface (in water below the sea ice) infrastructure allowed controlled manipulation of seawater pH levels (reduced by 0.4 pH units below ambient) in 2 chambers placed on the sea floor over natural benthic communities. Two control chambers (no pH manipulation) and two open plots (no chambers, no pH manipulation) were also sampled to compare to the pH manipulated (acidified) treatment chambers.  
Details of the antFOCE experiment can be found in the report – "antFOCE 2014/15 – Experimental System, Deployment, Sampling and Analysis". This report and a diagram indicating how the various antFOCE data sets relate to each other are available at:
https://data.aad.gov.au/metadata/records/AAS_4127_antFOCE_Project4127
AntFOCE biofilm DNA methods
Laurence Clarke, Shane Powell, Bruce Deagle
DNA extraction
The biofilm was removed from the top of each slide with a cotton swab and DNA extracted directly from the swab using the MoBio PowerBiofilm DNA isolation kit following the manufacturer’s protocol. Extraction blanks were extracted in parallel to detect contamination. 
Eukaryotic 18S rDNA PCR amplification and high-throughput sequencing
DNA extracts were PCR-amplified in triplicate with primers designed to amplify 140-170 bp of eukaryotic 18S ribosomal DNA (Jarman et al. 2013). The forward primer was modified to improve amplification of protists.
Table 1. First and second round primers, including MID tags (Xs).
ILF_ProSSU3'F_X     TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG XXXXXX CACCGCCCGTCGCWMCTACCG
ILR_SSU3'R_Y    GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG XXXXXX GGTTCACCTACGGAAACCTTGTTACG
msqFX   AATGATACGGCGACCACCGAGATCTACAC XXXXXXXXXX TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG
msqRY   CAAGCAGAAGACGGCATACGAGAT XXXXXXXXXX GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG
PCR amplifications were performed in two rounds, the first to amplify the 18S region and add sample-specific multiplex-identifier (MID) tags and Illumina sequencing primers, the second to add the P5 and P7 sequencing adapters and additional MIDs.
Each reaction mix for the first PCR contained 0.1 µM each of forward and reverse primer, 0.2 µg/µL BSA, 0.2 U Phusion DNA polymerase in 1 x Phusion Master Mix (New England Biolabs, Ipswich, MA, USA) and 1 micro L DNA extract in a total reaction volume of 10 micro L. PCR thermal cycling conditions were initial denaturation at 98 degrees C for 30 secs, followed by 25 cycles of 98 degrees C for 5 secs, 67 degrees C for 20 secs and 72 degrees C for 20 secs, with a final extension at 72 degrees C for 5 min. Replicate PCR products were pooled then diluted 1:10 and Illumina sequencing adapters added in a second round of PCR using the same reaction mix and thermal cycling conditions as the first round, except the concentration of BSA was halved (0.1 micro g/micro L), and the number of cycles was reduced to 10 with an annealing temperature of 55 degrees C.
Products from each round of PCR were visualized on 2% agarose gels. Second round PCR products were pooled in equimolar ratios based on band intensity. The pooled products were purified using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA) and the concentration of the library measured using the Qubit dsDNA HS assay on a QUBIT 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA). The pool was diluted to 2 nM and paired-end reads generated on a MiSeq (Illumina, San Diego, CA, USA) with MiSeq Reagent Nano kit vs (300-cycles). 
Bacterial 16S rDNA PCR amplification and high-throughput sequencing
Bioinformatics
Reads were sorted by sample-specific MIDs added in the second round PCR using the MiSeq Reporter software. Fastq reads were merged using the -fastq_mergepairs command in USEARCH v8.0.1623 (Edgar 2010). Merged reads were sorted by "internal" 6 bp MID tags, and locus-specific primers trimmed with custom R scripts using the ShortRead package (Morgan et al. 2009), with only reads containing perfect matches to the expected MIDs and primers retained. Reads for all samples were dereplicated and global singletons discarded (-derep_fulllength -minuniquesize 2), and clustered into OTUs with the UPARSE algorithm (Edgar 2013) using the '-cluster_otus' command. Potentially chimeric reads were also discarded during this step. Reads for each sample were then assigned to OTUs (-usearch_global -id .97), and an OTU table generated using a custom R script.
Taxonomy was assigned to each OTU using MEGAN version 5.10.5 (Huson et al. 2011) based on 50 hits per OTU generated by BLASTN searches against the NCBI 'nt' database (downloaded August 2015). Default LCA parameters were used, except Min support = 1, Min score = 100, Top percent = 10. Alpha and beta-diversity analyses were performed based on a rarefied OTU table with QIIME v1.8.0 (alpha_rarefaction.py, beta_diversity_through_plots.py, Caporaso et al. 2010). 
References
Caporaso JG, Kuczynski J, Stombaugh J, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nature Methods 7, 335-336.
Huson DH, Mitra S, Ruscheweyh HJ, Weber N, Schuster SC (2011) Integrative analysis of environmental sequences using MEGAN4. Genome Research 21, 1552-1560.
Jarman SN, McInnes JC, Faux C, et al. (2013) Adelie penguin population diet monitoring by analysis of food DNA in scats. PLoS One 8, e82227.

Tags: ADS > Automated DNA Sequencer, AMD, AMD/AU, BIOFILMS, CEOS, CONTINENT > ANTARCTICA, EARTH SCIENCE > BIOSPHERE > ECOSYSTEMS > MARINE ECOSYSTEMS > BENTHIC, EARTH SCIENCE > OCEANS > OCEAN CHEMISTRY > CARBON DIOXIDE, EUKARYOTES, FIELD SURVEYS, GEOGRAPHIC REGION > POLAR, LABORATORY, MICROSCOPES, OCEAN > SOUTHERN OCEAN, OCEAN ACIDIFICATION
Organization: Australian Ocean Data Network