Chemosensory function of Varroa gnathosoma: transcriptomic and proteomic analyses.

In this study, we evaluated the role of the gnathosoma (mouthparts) in chemosensing of the most devastating honey bee parasite, Varroa destructor mite. Through transcriptomic analysis, we compared the expression of putative chemosensory genes between the body parts containing the main chemosensory organs (the forelegs), gnathosoma and the rest of the body devoid of these two body parts. Furthermore, we checked the presence of chemosensory-related transcripts in the proteome of the gnathosoma. Our comparative transcriptomic analysis revealed the presence of 83 transcripts with known characteristic conserved domains belonging to eight chemosensory gene families in the three Varroa transcriptomes. Among these transcripts, 11 were significantly upregulated in the mite's forelegs, compared to 8 and 10 in the gnathosoma and body devoid of both organs, respectively. Whilst the gnathosoma and the forelegs share similar expression of some putative lipid carrier proteins, membrane-bound receptors, and associated proteins, they also differ in the expression profiles of some transcripts belonging to these protein families. This suggests two functional chemosensory organs that may differ in their chemosensory function according to specific characteristics of compounds they detect. Moreover, the higher expression of some chemosensory transcripts in the body devoid of forelegs and gnathosoma compared to the gnathosoma alone, may suggest the presence of additional function of these transcripts or alternatively presence of additional external or internal chemosensory organs. Insights into the functional annotation of a highly expressed gustatory receptor present in both organs using RNA interference (RNAi) are also revealed.

Automated phenol-chloroform extraction of high molecular weight genomic DNA for use in long-read single-molecule sequencing.

Background: Automation has increasingly become more commonplace in the research laboratory workspace. The introduction of articulated robotic arms allows the researcher more flexibility in the tasks a single piece of automated machinery can perform. We set out to incorporate automation in processing of genomic DNA organic extractions to increase throughput and limit researchers to the exposure of organic solvents. Methods: In order to automate the genome sequencing pipeline in our laboratory, we programmed a dual-arm anthropomorphic robot, the Robotic Biology Institute's Maholo LabDroid, to perform organic solvent-based genomic DNA extraction from cell lysates. To the best of our knowledge, this is the first time that automation of phenol-chloroform extraction has been reported. Results: We achieved routine extraction of high molecular weight genomic DNA (>100 kb) from diverse biological samples including algae cultured in sea water, bacteria, whole insects, and human cell lines. The results of pulse-field electrophoresis size analysis and the N50 sequencing metrics of reads obtained from Nanopore MinION runs verified the presence of intact DNA suitable for direct sequencing. Conclusions: We present the workflow that can be used to program similar robots and discuss the problems and solutions we encountered in developing the workflow. The protocol can be adapted to analogous methods such as RNA extraction, and there is ongoing work to incorporate further post-extraction steps such as library construction. This work shows the potential for automated robotic workflows to free molecular biological researchers from manual interventions in routine experimental work. A time-lapse movie of the entire automated run is included in this report.

Correction to "Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases".

[This corrects the article DOI: 10.1021/jacsau.1c00464.].

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases.

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

The CCR4-NOT deadenylase complex safeguards thymic positive selection by down-regulating aberrant pro-apoptotic gene expression.

A repertoire of T cells with diverse antigen receptors is selected in the thymus. However, detailed mechanisms underlying this thymic positive selection are not clear. Here we show that the CCR4-NOT complex limits expression of specific genes through deadenylation of mRNA poly(A) tails, enabling positive selection. Specifically, the CCR4-NOT complex is up-regulated in thymocytes before initiation of positive selection, where in turn, it inhibits up-regulation of pro-apoptotic Bbc3 and Dab2ip. Elimination of the CCR4-NOT complex permits up-regulation of Bbc3 during a later stage of positive selection, inducing thymocyte apoptosis. In addition, CCR4-NOT elimination up-regulates Dab2ip at an early stage of positive selection. Thus, CCR4-NOT might control thymocyte survival during two-distinct stages of positive selection by suppressing expression levels of pro-apoptotic molecules. Taken together, we propose a link between CCR4-NOT-mediated mRNA decay and T cell selection in the thymus.

Hidden proteome of synaptic vesicles in the mammalian brain.

Current proteomic studies clarified canonical synaptic proteins that are common to many types of synapses. However, proteins of diversified functions in a subset of synapses are largely hidden because of their low abundance or structural similarities to abundant proteins. To overcome this limitation, we have developed an "ultra-definition" (UD) subcellular proteomic workflow. Using purified synaptic vesicle (SV) fraction from rat brain, we identified 1,466 proteins, three times more than reported previously. This refined proteome includes all canonical SV proteins, as well as numerous proteins of low abundance, many of which were hitherto undetected. Comparison of UD quantifications between SV and synaptosomal fractions has enabled us to distinguish SV-resident proteins from potential SV-visitor proteins. We found 134 SV residents, of which 86 are present in an average copy number per SV of less than one, including vesicular transporters of nonubiquitous neurotransmitters in the brain. We provide a fully annotated resource of all categorized SV-resident and potential SV-visitor proteins, which can be utilized to drive novel functional studies, as we characterized here Aak1 as a regulator of synaptic transmission. Moreover, proteins in the SV fraction are associated with more than 200 distinct brain diseases. Remarkably, a majority of these proteins was found in the low-abundance proteome range, highlighting its pathological significance. Our deep SV proteome will provide a fundamental resource for a variety of future investigations on the function of synapses in health and disease.

Phylogenetic comparisons reveal mosaic histories of larval and adult shell matrix protein deployment in pteriomorph bivalves.

Molluscan shells are organo-mineral composites, in which the dominant calcium carbonate is intimately associated with an organic matrix comprised mainly of proteins and polysaccharides. However, whether the various shell matrix proteins (SMPs) date to the origin of hard skeletons in the Cambrian, or whether they represent later deployment through adaptive evolution, is still debated. In order to address this issue and to better understand the origins and evolution of biomineralization, phylogenetic analyses have been performed on the three SMP families, Von Willebrand factor type A (VWA) and chitin-binding domain-containing protein (VWA-CB dcp), chitobiase, and carbonic anhydrase (CA), which exist in both larval and adult shell proteomes in the bivalves, Crassostrea gigas and Pinctada fucata. In VWA-CB dcp and chitobiase, paralogs for larval and adult SMPs evolved before the divergence of these species. CA-SMPs have been taken as evidence for ancient origins of SMPs by their presumed indispensable function in biomineralization and ubiquitous distribution in molluscs. However, our results indicate gene duplications that gave rise to separate deployments as larval and adult CA-SMPs occurred independently in each lineage after their divergence, which is considerably more recent than hitherto assumed, supporting the "recent heritage and fast evolution" scenario for SMP evolution.

miR2118-dependent U-rich phasiRNA production in rice anther wall development.

Reproduction-specific small RNAs are vital regulators of germline development in animals and plants. MicroRNA2118 (miR2118) is conserved in plants and induces the production of phased small interfering RNAs (phasiRNAs). To reveal the biological functions of miR2118, we describe here rice mutants with large deletions of the miR2118 cluster. Our results demonstrate that the loss of miR2118 causes severe male and female sterility in rice, associated with marked morphological and developmental abnormalities in somatic anther wall cells. Small RNA profiling reveals that miR2118-dependent 21-nucleotide (nt) phasiRNAs in the anther wall are U-rich, distinct from the phasiRNAs in germ cells. Furthermore, the miR2118-dependent biogenesis of 21-nt phasiRNAs may involve the Argonaute proteins OsAGO1b/OsAGO1d, which are abundant in anther wall cell layers. Our study highlights the site-specific differences of phasiRNAs between somatic anther wall and germ cells, and demonstrates the significance of miR2118/U-phasiRNA functions in anther wall development and rice reproduction.

Rice Histone Propionylation and Generation of Chemically Derivatized Synthetic H3 and H4 Peptides for Identification of Acetylation Sites and Quantification.

Histone proteins are crucial in the study of chromatin dynamics owing to their wide-ranging implications in the regulation of gene expression. Modifications of histones are integral to these regulatory processes in concert with associated proteins, such as transcription factors and coactivators. One of the biochemical techniques available to enhance analysis of histone proteins is chemical derivatization using propionic anhydride. In this protocol, we describe the use of propionylation to efficiently derivatize acid-extracted histones from rice. We also synthesize H3 and H4 tryptic peptides, thus mimicking the nature of derivatized extracted peptides to aid in identification and quantification using targeted-mass spectrometry. Here we make available the masses of the precursor ions and the retention times (RT) of each synthesized peptide. These provide useful information to facilitate histone data analysis. Lastly, we note that we will distribute these synthetic peptides in nanomolar (nM) concentrations to those who wish to utilize them for assays and further experimental studies.

High-resolution archaellum structure reveals a conserved metal-binding site.

Many archaea swim by means of archaella. While the archaellum is similar in function to its bacterial counterpart, its structure, composition, and evolution are fundamentally different. Archaella are related to archaeal and bacterial type IV pili. Despite recent advances, our understanding of molecular processes governing archaellum assembly and stability is still incomplete. Here, we determine the structures of Methanococcus archaella by X-ray crystallography and cryo-EM The crystal structure of Methanocaldococcus jannaschii FlaB1 is the first and only crystal structure of any archaellin to date at a resolution of 1.5 Å, which is put into biological context by a cryo-EM reconstruction from Methanococcus maripaludis archaella at 4 Å resolution created with helical single-particle analysis. Our results indicate that the archaellum is predominantly composed of FlaB1. We identify N-linked glycosylation by cryo-EM and mass spectrometry. The crystal structure reveals a highly conserved metal-binding site, which is validated by mass spectrometry and electron energy-loss spectroscopy. We show in vitro that the metal-binding site, which appears to be a widespread property of archaellin, is required for filament integrity.

Dual Gene Repertoires for Larval and Adult Shells Reveal Molecules Essential for Molluscan Shell Formation.

Molluscan shells, mainly composed of calcium carbonate, also contain organic components such as proteins and polysaccharides. Shell organic matrices construct frameworks of shell structures and regulate crystallization processes during shell formation. To date, a number of shell matrix proteins (SMPs) have been identified, and their functions in shell formation have been studied. However, previous studies focused only on SMPs extracted from adult shells, secreted after metamorphosis. Using proteomic analyses combined with genomic and transcriptomic analyses, we have identified 31 SMPs from larval shells of the pearl oyster, Pinctada fucata, and 111 from the Pacific oyster, Crassostrea gigas. Larval SMPs are almost entirely different from those of adults in both species. RNA-seq data also confirm that gene expression profiles for larval and adult shell formation are nearly completely different. Therefore, bivalves have two repertoires of SMP genes to construct larval and adult shells. Despite considerable differences in larval and adult SMPs, some functional domains are shared by both SMP repertoires. Conserved domains include von Willebrand factor type A (VWA), chitin-binding (CB), carbonic anhydrase (CA), and acidic domains. These conserved domains are thought to play crucial roles in shell formation. Furthermore, a comprehensive survey of animal genomes revealed that the CA and VWA-CB domain-containing protein families expanded in molluscs after their separation from other Lophotrochozoan linages such as the Brachiopoda. After gene expansion, some family members were co-opted for molluscan SMPs that may have triggered to develop mineralized shells from ancestral, nonmineralized chitinous exoskeletons.

Organic and Peptidyl Constituents of Snake Venoms: The Picture Is Vastly More Complex Than We Imagined.

Small metabolites and peptides in 17 snake venoms (Elapidae, Viperinae, and Crotalinae), were quantified using liquid chromatography-mass spectrometry. Each venom contains >900 metabolites and peptides. Many small organic compounds are present at levels that are probably significant in prey envenomation, given that their known pharmacologies are consistent with snake envenomation strategies. Metabolites included purine nucleosides and their bases, neurotransmitters, neuromodulators, guanidino compounds, carboxylic acids, amines, mono- and disaccharides, and amino acids. Peptides of 2⁻15 amino acids are also present in significant quantities, particularly in crotaline and viperine venoms. Some constituents are specific to individual taxa, while others are broadly distributed. Some of the latter appear to support high anabolic activity in the gland, rather than having toxic functions. Overall, the most abundant organic metabolite was citric acid, owing to its predominance in viperine and crotaline venoms, where it chelates divalent cations to prevent venom degradation by venom metalloproteases and damage to glandular tissue by phospholipases. However, in terms of their concentrations in individual venoms, adenosine, adenine, were most abundant, owing to their high titers in Dendroaspis polylepis venom, although hypoxanthine, guanosine, inosine, and guanine all numbered among the 50 most abundant organic constituents. A purine not previously reported in venoms, ethyl adenosine carboxylate, was discovered in D. polylepis venom, where it probably contributes to the profound hypotension caused by this venom. Acetylcholine was present in significant quantities only in this highly excitotoxic venom, while 4-guanidinobutyric acid and 5-guanidino-2-oxopentanoic acid were present in all venoms.

Chitin-based barrier immunity and its loss predated mucus-colonization by indigenous gut microbiota.

Mammalian gut microbiota are integral to host health. However, how this association began remains unclear. We show that in basal chordates the gut space is radially compartmentalized into a luminal part where food microbes pass and an almost axenic peripheral part, defined by membranous delamination of the gut epithelium. While this membrane, framed with chitin nanofibers, structurally resembles invertebrate peritrophic membranes, proteome supports its affinity to mammalian mucus layers, where gut microbiota colonize. In ray-finned fish, intestines harbor indigenous microbes, but chitinous membranes segregate these luminal microbes from the surrounding mucus layer. These data suggest that chitin-based barrier immunity is an ancient system, the loss of which, at least in mammals, provided mucus layers as a novel niche for microbial colonization. These findings provide a missing link for intestinal immune systems in animals, revealing disparate mucosal environment in model organisms and highlighting the loss of a proven system as innovation.

The putative ceramide-conjugation protein Cwh43 regulates G0 quiescence, nutrient metabolism and lipid homeostasis in fission yeast.

Cellular nutrient states control whether cells proliferate, or whether they enter or exit quiescence. Here, we report characterizations of fission yeast temperature-sensitive (ts) mutants of the evolutionarily conserved transmembrane protein Cwh43, and explore its relevance to utilization of glucose, nitrogen source and lipids. GFP-tagged Cwh43 localizes at ER associated with the nuclear envelope and the plasma membrane, as in budding yeast. We found that cwh43 mutants failed to divide in low glucose and lost viability during quiescence under nitrogen starvation. In cwh43 mutants, comprehensive metabolome analysis demonstrated dramatic changes in marker metabolites that altered under low glucose and/or nitrogen starvation, although cwh43 cells apparently consumed glucose in the culture medium. Furthermore, we found that cwh43 mutant cells had elevated levels of triacylglycerols (TGs) and coenzyme A, and that they accumulated lipid droplets. Notably, TG biosynthesis was required to maintain cell division in the cwh43 mutant. Thus, Cwh43 affects utilization of glucose and nitrogen sources, as well as storage lipid metabolism. These results may fit a notion developed in budding yeast stating that Cwh43 conjugates ceramide to glycosylphosphatidylinositol (GPI)-anchored proteins and maintains integrity of membrane organization.

Genetic defects in SAPK signalling, chromatin regulation, vesicle transport and CoA-related lipid metabolism are rescued by rapamycin in fission yeast.

Rapamycin inhibits TOR (target of rapamycin) kinase, and is being used clinically to treat various diseases ranging from cancers to fibrodysplasia ossificans progressiva. To understand rapamycin mechanisms of action more comprehensively, 1014 temperature-sensitive (ts) fission yeast (Schizosaccharomyces pombe) mutants were screened in order to isolate strains in which the ts phenotype was rescued by rapamycin. Rapamycin-rescued 45 strains, among which 12 genes responsible for temperature sensitivity were identified. These genes are involved in stress-activated protein kinase (SAPK) signalling, chromatin regulation, vesicle transport, and CoA- and mevalonate-related lipid metabolism. Subsequent metabolome analyses revealed that rapamycin upregulated stress-responsive metabolites, while it downregulated purine biosynthesis intermediates and nucleotide derivatives. Rapamycin alleviated abnormalities in cell growth and cell division caused by sty1 mutants (Δsty1) of SAPK. Notably, in Δsty1, rapamycin reduced greater than 75% of overproduced metabolites (greater than 2× WT), like purine biosynthesis intermediates and nucleotide derivatives, to WT levels. This suggests that these compounds may be the points at which the SAPK/TOR balance regulates continuous cell proliferation. Rapamycin might be therapeutically useful for specific defects of these gene functions.

Coralsnake Venomics: Analyses of Venom Gland Transcriptomes and Proteomes of Six Brazilian Taxa.

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2-6 toxin classes that account for 91-99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A2 (PLA2s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA2s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1-2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6-9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the ß-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A2s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15-35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.

Diversity of Microbial Communities and Quantitative Chemodiversity in Layers of Marine Sediment Cores from a Causeway (Kaichu-Doro) in Okinawa Island, Japan.

Microbial community diversity and chemodiversity were investigated in marine sediments adjacent to the Okinawan "Kaichu-Doro" Causeway, which was constructed 46 years ago to connect a group of four islands (Henza-jima, Miyagi-jima, Ikei-jima, Hamahiga-jima) to the Okinawan main island. This causeway was not built on pilings, but by land reclamation; hence, it now acts as a long, thin peninsula. The construction of this causeway was previously shown to have influenced the surrounding marine ecosystem, causing ecosystem fragmentation and loss of water circulation. In this study, we collected sediment cores (n = 10) from five paired sites in 1 m water depths. Each pair of sites consisted of one site each on the immediate north and south sides of the causeway. Originally the members of each pair were much closer to each other (<150 m) than to other pairs, but now the members of each pair are isolated by the causeway. Each core was 60-80 cm long and was divided into 15-cm layers. We examined the vertical diversity of microbial communities and chemical compounds to determine the correlation between chemodiversity and microbial communities among marine sediment cores and layers. Principal coordinate analyses (PCoA) of detected compounds and of bacterial and archaeal operational taxonomic units (OTUs) revealed that the north and south sides of the causeway are relatively isolated, with each side having unique microbial OTUs. Additionally, some bacterial families (e.g., Acidaminobacteraceae, Rhizobiaceae, and Xanthomonadaceae) were found only on the south side of Kaichu-Doro. Interestingly, we found that the relative abundance of OTUs for some microbial families increased from top to bottom, but this was reversed in some other families. We conclude that the causeway has altered microbial community composition and metabolite profiles in marine sediments.

Polyamines as Snake Toxins and Their Probable Pharmacological Functions in Envenomation.

While decades of research have focused on snake venom proteins, far less attention has been paid to small organic venom constituents. Using mostly pooled samples, we surveyed 31 venoms (six elapid, six viperid, and 19 crotalid) for spermine, spermidine, putrescine, and cadaverine. Most venoms contained all four polyamines, although some in essentially trace quantities. Spermine is a potentially significant component of many viperid and crotalid venoms (≤0.16% by mass, or 7.9 µmol/g); however, it is almost completely absent from elapid venoms assayed. All elapid venoms contained larger molar quantities of putrescine and cadaverine than spermine, but still at levels that are likely to be biologically insignificant. As with venom purines, polyamines impact numerous physiological targets in ways that are consistent with the objectives of prey envenomation, prey immobilization via hypotension and paralysis. Most venoms probably do not contain sufficient quantities of polyamines to induce systemic effects in prey; however, local effects seem probable. A review of the pharmacological literature suggests that spermine could contribute to prey hypotension and paralysis by interacting with N-methyl-d-aspartate (NMDA) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors, nicotinic and muscarinic acetylcholine receptors, γ-Aminobutyric acid (GABA) receptors, blood platelets, ryanodine receptors, and Ca2+-ATPase. It also blocks many types of cation-permeable channels by interacting with negatively charged amino acid residues in the channel mouths. The site of envenomation probably determines which physiological targets assume the greatest importance; however, venom-induced liberation of endogenous, intracellular stores of polyamines could potentially have systemic implications and may contribute significantly to envenomation sequelae.

Developmental Progression in the Coral Acropora digitifera Is Controlled by Differential Expression of Distinct Regulatory Gene Networks.

Corals belong to the most basal class of the Phylum Cnidaria, which is considered the sister group of bilaterian animals, and thus have become an emerging model to study the evolution of developmental mechanisms. Although cell renewal, differentiation, and maintenance of pluripotency are cellular events shared by multicellular animals, the cellular basis of these fundamental biological processes are still poorly understood. To understand how changes in gene expression regulate morphogenetic transitions at the base of the eumetazoa, we performed quantitative RNA-seq analysis duringAcropora digitifera's development. We collected embryonic, larval, and adult samples to characterize stage-specific transcription profiles, as well as broad expression patterns. Transcription profiles reconstructed development revealing two main expression clusters. The first cluster grouped blastula and gastrula and the second grouped subsequent developmental time points. Consistently, we observed clear differences in gene expression between early and late developmental transitions, with higher numbers of differentially expressed genes and fold changes around gastrulation. Furthermore, we identified three coexpression clusters that represented discrete gene expression patterns. During early transitions, transcriptional networks seemed to regulate cellular fate and morphogenesis of the larval body. In late transitions, these networks seemed to play important roles preparing planulae for switch in lifestyle and regulation of adult processes. Although developmental progression inA. digitiferais regulated to some extent by differential coexpression of well-defined gene networks, stage-specific transcription profiles appear to be independent entities. While negative regulation of transcription is predominant in early development, cell differentiation was upregulated in larval and adult stages.

Snake venoms are integrated systems, but abundant venom proteins evolve more rapidly.

BACKGROUND: While many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown. RESULTS: Using next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A2 transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation. CONCLUSIONS: Given log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.