Codon Usage of Secondary Metabolite Genes in Desert Soils and Marine Sponges

This fileset contains the initial analysis of codon usage on secondary metabolite genes in desert soils (from Reddy et al 2012) and marine sponge microboim (from Trindade-Silva et al 2012). <strong>Why Does this Exist:</strong> I was wondering what were the similarites and differences in how microbial communties use secondary metabolite genes (genes that produce molecules that are sometimes useful to human as antibtiotics and as other things). My PhD work is focusing on looking at these genes across several different environments in New Mexico (caves, sides of cliffs, and springs). I thought that there might be a connection between environmental conditions and the types and use of these genes. In this data set I am looking at single domains of larger gene clusters. <strong>Methods:</strong> The Reddy et all dataset is available as an sra and fastq file under the identifier SRR342214 and the marine sponge KS domains are available under the accession numbers: JX012425:JX012657. The SRR342214 data set was mined with a custom python script to pull out all the barcoded KS domains and then remove domains with a size less then 150 nts. Codonw was used to calculate all condon bias indices and dinucleotide frequencies. The visualization of the data set was done in R studio using: ggplot2, gridExtrac, and FactoMineR. <strong>Questions Raised:</strong> <em>Does this single domain reflect the gene cluster?</em> Tenatively I say yes based on some whole gene cluster analysis I did. <em>How much is random GC mutation as opposed to other pressures acting on the NP genes in the community?</em> <em>Are there quantfiable difference between and within communites in ther codon usage of the KS domans?</em> With this small dataset we can see a difference between marine and desert soils. Will this hold up for my larger dataset? <em>Does this tell us something useful about how the communites share their genes within and what pressure select for certains types of codon bias in the NP genes?</em> <strong>References:</strong> <em>Codonw</em>, John Peden, Oxford University, available at http://bioweb.pasteur.fr/seqanal/interfaces/codonw.html <em>Natural product biosynthetic gene diversity in geographically distinct soil microbiomes.</em><br>Appl Environ Microbiol. 2012 May;78(10):3744-52. doi: 10.1128/AEM.00102-12. Epub 2012 Mar 16. Reddy BV, Kallifidas D, Kim JH, Charlop-Powers Z, Feng Z, Brady SF. <em>Polyketide synthase gene diversity within the microbiome of the sponge Arenosclera brasiliensis, endemic to the Southern Atlantic Ocean.</em> Appl Environ Microbiol. 2013 Mar;79(5):1598-605. doi: 10.1128/AEM.03354-12. Epub. 2012 Dec 28. Trindade-Silva AE, Rua CP, Andrade BG, Vicente AC, Silva GG, Berlinck RG, Thompson FL. R Core Team (2013). <em>R: A language and environment for statistical computing. R Foundation for Statistical Computing</em>, Vienna, Austria. URL http://www.R-project.org/. H. Wickham. <em>ggplot2: elegant graphics for data analysis</em>. Springer New York, 2009. Francois Husson, Julie Josse, Sebastien Le and Jeremy Mazet (2013). <em>FactoMineR: Multivariate Exploratory Data Analysis and Data Mining with R</em>. R package version 1.25. http://CRAN.R-project.org/package=FactoMineR