title: The venom-gland transcriptome of the eastern coral snake (Micrurus fulvius) reveals high venom complexity in the intragenomic evolution of venoms.
abstract: BACKGROUND
  Snake venom is shaped by the ecology and evolution of venomous species, and signals of positive selection in toxins have been consistently documented, reflecting the role of venoms as an ecologically critical phenotype. New World coral snakes (Elapidae) are represented by three genera and over 120 species and subspecies that are capable of causing significant human morbidity and mortality, yet coral-snake venom composition is poorly understood in comparison to that of Old World elapids. High-throughput sequencing is capable of identifying thousands of loci, while providing characterizations of expression patterns and the molecular evolutionary forces acting within the venom gland.
  
  
  RESULTS
  We describe the de novo assembly and analysis of the venom-gland transcriptome of the eastern coral snake (Micrurus fulvius). We identified 1,950 nontoxin transcripts and 116 toxin transcripts. These transcripts accounted for 57.1% of the total reads, with toxins accounting for 45.8% of the total reads. Phospholipases A(2) and three-finger toxins dominated expression, accounting for 86.0% of the toxin reads. A total of 15 toxin families were identified, revealing venom complexity previously unknown from New World coral snakes. Toxins exhibited high levels of heterozygosity relative to nontoxins, and overdominance may favor gene duplication leading to the fixation of advantageous alleles. Phospholipase A(2) expression was uniformly distributed throughout the class while three-finger toxin expression was dominated by a handful of transcripts, and phylogenetic analyses indicate that toxin divergence may have occurred following speciation. Positive selection was detected in three of the four most diverse toxin classes, suggesting that venom diversification is driven by recurrent directional selection.
  
  
  CONCLUSIONS
  We describe the most complete characterization of an elapid venom gland to date. Toxin gene duplication may be driven by heterozygote advantage, as the frequency of polymorphic toxin loci was significantly higher than that of nontoxins. Diversification among toxins appeared to follow speciation reflecting species-specific adaptation, and this divergence may be directly related to dietary shifts and is suggestive of a coevolutionary arms race.
authors.affiliation: Department of Biological Science, Florida State University, Tallahassee, FL 32306-4295, USA.
authors.forename: Mark J
authors.initials: MJ
authors.lastname: Margres
authors.forename: Karalyn
authors.initials: K
authors.lastname: Aronow
authors.forename: Jacob
authors.initials: J
authors.lastname: Loyacano
authors.forename: Darin R
authors.initials: DR
authors.lastname: Rokyta
mesh_terms: D000818:Animals; D004546:Elapid Venoms / Q000235:genetics*; D017815:Elapidae; D019143:Evolution, Molecular; D020440:Gene Duplication; D020869:Gene Expression Profiling; D049810:Genetic Speciation; D023281:Genomics; D006579:Heterozygote; D045726:Metalloproteases / Q000235:genetics; D009498:Neurotoxins / Q000235:genetics; D054467:Phospholipases A2 / Q000235:genetics; D017421:Sequence Analysis
chemical_list: D004546:Elapid Venoms; D009498:Neurotoxins; C053972:micrurus venom; D054467:Phospholipases A2; D045726:Metalloproteases
references.citation: BMC Genomics. 2012 Apr 24;13:149
references.pmid: 22530590
references.citation: BMC Mol Biol. 2007 Dec 20;8:115
references.pmid: 18096037
references.citation: Toxicon. 2009 May;53(6):693-7
references.pmid: 19673084
references.citation: J Bacteriol. 2006 Feb;188(3):1134-42
references.pmid: 16428417
references.citation: Genetics. 2000 May;155(1):431-49
references.pmid: 10790415
references.citation: Comput Appl Biosci. 1997 Oct;13(5):555-6
references.pmid: 9367129
references.citation: Toxicon. 2012 Jan;59(1):132-42
references.pmid: 22108621
references.citation: Comp Biochem Physiol A Mol Integr Physiol. 2005 Jan;140(1):125-34
references.pmid: 15664321
references.citation: PLoS One. 2011;6(6):e20997
references.pmid: 21731638
references.citation: Mol Biol Evol. 2002 Jan;19(1):49-57
references.pmid: 11752189
references.citation: BMC Genomics. 2012 Jul 16;13:312
references.pmid: 23025625
references.citation: Toxicon. 2000 May;38(5):713-21
references.pmid: 10673162
references.citation: BMC Evol Biol. 2008 Feb 28;8:70
references.pmid: 18307759
references.citation: Toxicon. 2010 Feb-Mar;55(2-3):250-5
references.pmid: 19647760
references.citation: Gene. 2000 Dec 30;261(1):43-52
references.pmid: 11164036
references.citation: Toxicon. 2011 Apr;57(5):657-71
references.pmid: 21255598
references.citation: J Struct Funct Genomics. 2003;3(1-4):201-12
references.pmid: 12836699
references.citation: Toxicon. 2011 Feb;57(2):323-31
references.pmid: 21184770
references.citation: Comp Biochem Physiol C Toxicol Pharmacol. 2001 Mar;128(3):425-56
references.pmid: 11255115
references.citation: Proc Biol Sci. 2009 Jul 7;276(1666):2443-9
references.pmid: 19364745
references.citation: Toxicon. 2010 Jul;55(8):1463-74
references.pmid: 20227433
references.citation: BMC Genomics. 2009 Nov 30;10:564
references.pmid: 19948012
references.citation: Cell Mol Life Sci. 2008 Dec;65(24):4039-54
references.pmid: 18979207
references.citation: Bioinformatics. 2001 Aug;17(8):754-5
references.pmid: 11524383
references.citation: J Cell Biol. 2002 Dec 9;159(5):893-902
references.pmid: 12460986
references.citation: Toxicon. 2009 Mar 1;53(3):375-82
references.pmid: 19135468
references.citation: Toxicon. 1987;25(5):547-54
references.pmid: 3617089
references.citation: BMC Bioinformatics. 2009 Dec 03;10 Suppl 15:S2
references.pmid: 19958512
references.citation: BMC Genomics. 2009 Mar 16;10:112
references.pmid: 19291316
references.citation: PLoS Negl Trop Dis. 2010 Mar 09;4(3):e622
references.pmid: 20231886
references.citation: Mol Biol Evol. 1999 Mar;16(3):372-82
references.pmid: 10331263
references.citation: BMC Evol Biol. 2007 Jan 18;7:2
references.pmid: 17233905
references.citation: Bioinformatics. 2003 Aug 12;19(12):1572-4
references.pmid: 12912839
references.citation: J Mol Evol. 2003 Jul;57(1):110-29
references.pmid: 12962311
references.citation: J Proteomics. 2011 Aug 24;74(9):1810-25
references.pmid: 21803179
references.citation: Nature. 2011 Nov 16;479(7373):410-4
references.pmid: 22094702
references.citation: J Proteomics. 2011 Aug 24;74(9):1795-809
references.pmid: 21515432
references.citation: J Proteome Res. 2011 Apr 1;10(4):1816-27
references.pmid: 21280576
references.citation: Genetics. 2008 Aug;179(4):2319-24
references.pmid: 18689880
references.citation: J Proteomics. 2011 Dec 21;75(2):655-67
references.pmid: 21963438
references.citation: J Biol Chem. 2009 Apr 10;284(15):9885-91
references.pmid: 19208624
references.citation: Biochem J. 1971 Feb;121(4):609-12
references.pmid: 4940047
references.citation: J Mol Evol. 2008 Feb;66(2):151-66
references.pmid: 18253686
references.citation: Nucleic Acids Res. 1994 Nov 11;22(22):4673-80
references.pmid: 7984417
references.citation: Genetica. 1991;84(1):5-11
references.pmid: 1874440
references.citation: PLoS One. 2010 Jul 28;5(7):e11840
references.pmid: 20676378
references.citation: Adv Genet. 2002;46:451-83
references.pmid: 11931235