Continent-island boundary and environment-shaped evolution in the marine amphipod Ampithoe marcuzzii complex (Crustacea: Eumalacostraca: Ampithoidae).

Marine amphipods are crustaceans that lack a larval phase and consequently have low dispersion rates. Despite that, these crustaceans present a remarkable ability to be transported by rafting on natural floating substrata, especially macroalgae, where they find shelter, food and a mating ground. The species Ampithoe marcuzzii is widely distributed throughout the western Atlantic Ocean. Here, it was used as a model to study seascape genomics and phylogeography in invertebrates with low dispersion capacities. We anticipated that the lineages would present isolation-by-distance patterns. However, surface currents and other abiotic variables could facilitate connectivity among distant sites. Based on mitochondrial and nuclear genes, SNPs, and environmental associations, we observed the presence of a species complex within A. marcuzzii, separating mainland and insular populations. Each species showed an independent evolutionary history, with a strong latitudinal population structure and evidence of isolation-by-distance and isolation-by-environment, characterizing the 'continent' species. Historical expansion and environmental variables were observed associated with the southeastern population, and ecological niche modeling corroborated the region as a paleorefuge. Conversely, populations from 'islands' presented complicated evolutionary histories, with closer localities genetically isolated and distant localities connected. These findings indicate that insular populations with low dispersion capacity might be more susceptible to spatial connectivity by floating substrata and to changes in surface currents. In contrast, mainland populations might be more vulnerable to local climate changes due to lack of gene flow.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication.

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Comparing iDNA from mosquitoes and flies to survey mammals in a semi-controlled Neotropical area.

Ingested-derived DNA (iDNA) from insects represents a powerful tool for assessing vertebrate diversity because insects are easy to sample, have a diverse diet and are widely distributed. Because of these advantages, the use of iDNA for detecting mammals has gained increasing attention. Here we aimed to compare the effectiveness of mosquitoes and flies to detect mammals with a small sampling effort in a semi-controlled area, a zoo that houses native and non-native species. We compared mosquitoes and flies regarding the number of mammal species detected, the amount of mammal sequence reads recovered, and the flight distance range for detecting mammals. We also verified if the combination of two mini-barcodes (12SrRNA and 16SrRNA) would perform better than either mini-barcode alone to inform local mammal biodiversity from iDNA. To capture mosquitoes and flies, we distributed insect traps in eight sampling points during 5 days. We identified 43 Operational Taxonomic Units from 10 orders, from the iDNA of 17 mosquitoes and 46 flies. There was no difference in the number of species recovered per individual insect between mosquitoes and flies, but the number of flies captured was higher, resulting in more mammal species recovered by flies. Eight species were recorded exclusively by mosquitoes and 20 by flies, suggesting that using both samplers would allow a more comprehensive screening of the biodiversity. The maximum distance recorded was 337 m for flies and 289 m for mosquitoes, but the average range distance did not differ between insect groups. Our assay proved to be efficient for mammal detection, considering the high number of species detected with a reduced sampling effort.

Loss of S-nitrosoglutathione reductase disturbs phytohormone homeostasis and regulates shoot side branching and fruit growth in tomato.

S-Nitrosoglutathione plays a central role in nitric oxide (NO) homeostasis, and S-nitrosoglutathione reductase (GSNOR) regulates the cellular levels of S-nitrosoglutathione across kingdoms. Here, we investigated the role of endogenous NO in shaping shoot architecture and controlling fruit set and growth in tomato (Solanum lycopersicum). SlGSNOR silencing promoted shoot side branching and led to reduced fruit size, negatively impacting fruit yield. Greatly intensified in slgsnor knockout plants, these phenotypical changes were virtually unaffected by SlGSNOR overexpression. Silencing or knocking out of SlGSNOR intensified protein tyrosine nitration and S-nitrosation and led to aberrant auxin production and signaling in leaf primordia and fruit-setting ovaries, besides restricting the shoot basipetal polar auxin transport stream. SlGSNOR deficiency triggered extensive transcriptional reprogramming at early fruit development, reducing pericarp cell proliferation due to restrictions on auxin, gibberellin, and cytokinin production and signaling. Abnormal chloroplast development and carbon metabolism were also detected in early-developing NO-overaccumulating fruits, possibly limiting energy supply and building blocks for fruit growth. These findings provide new insights into the mechanisms by which endogenous NO fine-tunes the delicate hormonal network controlling shoot architecture, fruit set, and post-anthesis fruit development, emphasizing the relevance of NO-auxin interaction for plant development and productivity.

Coevolutionary Implications of Obligate Commensalism in Sea Turtles: the Case of the Genus Hyachelia Barnard, 1967 (Crustacea, Amphipoda).

Obligate commensalism in the marine environment and its evolutionary role are still poorly understood. Although sea turtles may serve as ideal substrates for epibionts, within amphipods, only the genus Hyachelia evolved in obligate commensalism with turtles. Here, we report a new host record for Hyachelia lowryi on the hawksbill turtle and describe a larger distribution of the genus in the Atlantic Ocean on green and loggerhead turtles. Hyachelia spp. were sampled from nesting sites of Caretta caretta and feeding grounds of Eretmochelys imbricata and Chelonia mydas along the Brazilian coast. Insights regarding the coevolution of this remarkable genus with its hosts based on molecular analyses are inferred based on mitochondrial (COI) and nuclear (18SrRNA) genes using new and previously available sequences from the infraorder Talitrida. Divergence times for Hyachelia are around the Cretaceous (~127.66 Mya), corresponding to an ancient origin and in agreement with modern green turtle (Chelonioidea) radiation. Later, diversification of Hyachelia species is dated at about 26 Mya, suggesting a coevolutionary association between amphipods and Carettini/Chelonini sea turtles.

Molecular evolution and signatures of selective pressures on Bos, focusing on the Nelore breed (Bos indicus).

Evolutionary history leads to genome changes over time, especially for species that have experienced intense selective pressures over a short period. Here, we investigated the genomic evolution of Bos species by searching for potential selection signatures, focusing on Nelore, an economically relevant cattle breed in Brazil. We assessed the genomic processes determining the molecular evolution across Nelore and thirteen other related taxa by evaluating (i) amino acid sequence conservation, (ii) the dN/dS ratio, and (iii) gene families' turnover rate (λ). Low conserved regions potentially associated with fatty acid metabolism seem to reflect differences in meat fat content in taxa with different evolutionary histories. All Bos species presented genes under positive selection, especially B. indicus and Nelore, which include transport protein cobalamin, glycolipid metabolism, and hormone signaling. These findings could be explained by constant selective pressures to obtain higher immune resistance and efficient metabolism. The gene contraction rate across the Nelore + B. indicus branch was almost nine times higher than that in other lineages (λ = 0.01043 vs. 0.00121), indicating gene losses during the domestication process. Amino acid biosynthesis, reproductive and innate immune system-related pathways were associated with genes recognized within the most frequent rapidly evolving gene families and in genes under positive selection, supporting the substantial relevance of such traits from a domestication perspective. Our data provide new insights into how the genome may respond to intense artificial selection in distinct taxa, and reinforces the presence of selective pressures on traits potentially relevant for future animal breeding investments.

Evolution, Expression Patterns, and Distribution of Novel Ribbon Worm Predatory and Defensive Toxins.

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

The state of Medusozoa genomics: current evidence and future challenges.

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Systems Biology Applied to the Study of Papaya Fruit Ripening: The Influence of Ethylene on Pulp Softening.

Papaya is a fleshy fruit that undergoes fast ethylene-induced modifications. The fruit becomes edible, but the fast pulp softening is the main factor that limits the post-harvest period. Papaya fast pulp softening occurs due to cell wall disassembling coordinated by ethylene triggering that massively expresses pectinases. In this work, RNA-seq analysis of ethylene-treated and non-treated papayas enabled a wide transcriptome overview that indicated the role of ethylene during ripening at the gene expression level. Several families of transcription factors (AP2/ERF, NAC, and MADS-box) were differentially expressed. ACO, ACS, and SAM-Mtase genes were upregulated, indicating a high rate of ethylene biosynthesis after ethylene treatment. The correlation among gene expression and physiological data demonstrated ethylene treatment can indeed simulate ripening, and regulation of changes in fruit color, aroma, and flavor could be attributed to the coordinated expression of several related genes. Especially about pulp firmness, the identification of 157 expressed genes related to cell wall metabolism demonstrated that pulp softening is accomplished by a coordinated action of several different cell wall-related enzymes. The mechanism is different from other commercially important fruits, such as strawberry, tomato, kiwifruit, and apple. The observed behavior of this new transcriptomic data confirms ethylene triggering is the main event that elicits fast pulp softening in papayas.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation.

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Multifaceted roles of nitric oxide in tomato fruit ripening: NO-induced metabolic rewiring and consequences for fruit quality traits.

Nitric oxide (NO) has been implicated as part of the ripening regulatory network in fleshy fruits. However, very little is known about the simultaneous action of NO on the network of regulatory events and metabolic reactions behind ripening-related changes in fruit color, taste, aroma and nutritional value. Here, we performed an in-depth characterization of the concomitant changes in tomato (Solanum lycopersicum) fruit transcriptome and metabolome associated with the delayed-ripening phenotype caused by NO supplementation at the pre-climacteric stage. Approximately one-third of the fruit transcriptome was altered in response to NO, including a multilevel down-regulation of ripening regulatory genes, which in turn restricted the production and tissue sensitivity to ethylene. NO also repressed hydrogen peroxide-scavenging enzymes, intensifying nitro-oxidative stress and S-nitrosation and nitration events throughout ripening. Carotenoid, tocopherol, flavonoid and ascorbate biosynthesis were differentially affected by NO, resulting in overaccumulation of ascorbate (25%) and flavonoids (60%), and impaired lycopene production. In contrast, the biosynthesis of compounds related to tomato taste (sugars, organic acids, amino acids) and aroma (volatiles) was slightly affected by NO. Our findings indicate that NO triggers extensive transcriptional and metabolic rewiring at the early ripening stage, modifying tomato antioxidant composition with minimal impact on fruit taste and aroma.

Shedding the Light on Litopenaeus vannamei Differential Muscle and Hepatopancreas Immune Responses in White Spot Syndrome Virus (WSSV) Exposure.

White Spot Syndrome Virus (WSSV) is one of the main threats to farming Litopenaeus vannamei, the most important crustacean commercialized in aquaculture worldwide. Here, we performed RNA-seq analyses in hepatopancreas and muscle from WSSV-negative (healthy) and WSSV-positive (unhealthy) L. vannamei, previously exposed to the virus, to obtain new insights about the molecular basis of resistance to WSSV. We detected 71% of our reads mapped against the recently described L. vannamei genome. This is the first report mapping RNA-seq transcripts from shrimps exposed to WSSV against the species reference genome. Differentially expressed gene (DEG) analyses were performed for four independent comparisons, and 13,338 DEGs were identified. When the redundancies and isoforms were disregarded, we observed 8351 and 6514 DEGs, respectively. Interestingly, after crossing the data, we detected a common set of DEGs for hepatopancreas and healthy shrimps, as well as another one for muscle and unhealthy shrimps. Our findings indicate that genes related to apoptosis, melanization, and the Imd pathway are likely to be involved in response to WSSV, offering knowledge about WSSV defense in shrimps exposed to the virus but not infected. These data present potential to be applied in further genetic studies in penaeids and other farmed shrimp species.

Role of neuroendocrine modulation and biochemistry in the sepsis in Piaractus mesopotamicus.

Sepsis is a systemic process with multifactorial pathophysiology that affects most animal species. It is responsible for high rates of morbidity and mortality. This work aimed to study the biochemical and neuroendocrine changes of the sepsis process in Piaractus mesopotamicus after Aeromonas hydrophila inoculation analyzing changes in blood leukocyte and differences in neuroendocrine-biochemical modulation using RNA-seq. Fish showed hypercortisolemia, inhibition of glucose absorption, followed by hypocortisolemia and then hyperglycemia. Thyroid hormones (T3 and T4) showed immediate decrease in serum and T4 increased 6 h post-inoculation (HPI). Sepsis-induced hormonal alterations triggered changes in the metabolic pathways increasing protein and lipid catabolism, use of transient anaerobic glycolysis and liver injury. A reference transcriptome was constructed based on blood leukocytes from P. mesopotamicus. The assembly resulted in total 266,272 contigs with a N50 of 2786 bp. There was a reorganization of plasma membrane of leukocytes at the beginning of the septic process with increased expression of neuroendocrine receptors and with continuous flow of neurotransmitters, hormones and solutes with compensatory regulation at 6 HPI. Three and nine HPI seemed to be critical, the expression of a number of transcription factors was increased, including the modulatory DEGs related to glucocorticoid and thyroid hormones induced and suppressed (FDR < 0.05). Neuroendocrine modulation can regulate leukocytes and biochemical parameters of peripheral blood, being important sources for the study of the pathophysiology of sepsis. These finding highlights the importance of further studies focusing on biochemical-neuroendocrine changes in blood leukocytes and systemic sepsis.

C4 and crassulacean acid metabolism within a single leaf: deciphering key components behind a rare photosynthetic adaptation.

Although biochemically related, C4 and crassulacean acid metabolism (CAM) systems are expected to be incompatible. However, Portulaca species, including P. oleracea, operate C4 and CAM within a single leaf, and the mechanisms behind this unique photosynthetic arrangement remain largely unknown. Here, we employed RNA-seq to identify candidate genes involved exclusively or shared by C4 or CAM, and provided an in-depth characterization of their transcript abundance patterns during the drought-induced photosynthetic transitions in P. oleracea. Data revealed fewer candidate CAM-specific genes than those recruited to function in C4 . The putative CAM-specific genes were predominantly involved in night-time primary carboxylation reactions and malate movement across the tonoplast. Analysis of gene transcript-abundance regulation and photosynthetic physiology indicated that C4 and CAM coexist within a single P. oleracea leaf under mild drought conditions. Developmental and environmental cues were shown to regulate CAM expression in stems, whereas the shift from C4 to C4 -CAM hybrid photosynthesis in leaves was strictly under environmental control. Moreover, efficient starch turnover was identified as part of the metabolic adjustments required for CAM operation in both organs. These findings provide insights into C4 /CAM connectivity and compatibility, contributing to a deeper understanding of alternative ways to engineer CAM into C4 crop species.

Landscape genomics to the rescue of a tropical bee threatened by habitat loss and climate change.

Habitat degradation and climate change are currently threatening wild pollinators, compromising their ability to provide pollination services to wild and cultivated plants. Landscape genomics offers powerful tools to assess the influence of landscape modifications on genetic diversity and functional connectivity, and to identify adaptations to local environmental conditions that could facilitate future bee survival. Here, we assessed range-wide patterns of genetic structure, genetic diversity, gene flow, and local adaptation in the stingless bee Melipona subnitida, a tropical pollinator of key biological and economic importance inhabiting one of the driest and hottest regions of South America. Our results reveal four genetic clusters across the species' full distribution range. All populations were found to be under a mutation-drift equilibrium, and genetic diversity was not influenced by the amount of reminiscent natural habitats. However, genetic relatedness was spatially autocorrelated and isolation by landscape resistance explained range-wide relatedness patterns better than isolation by geographic distance, contradicting earlier findings for stingless bees. Specifically, gene flow was enhanced by increased thermal stability, higher forest cover, lower elevations, and less corrugated terrains. Finally, we detected genomic signatures of adaptation to temperature, precipitation, and forest cover, spatially distributed in latitudinal and altitudinal patterns. Taken together, our findings shed important light on the life history of M. subnitida and highlight the role of regions with large thermal fluctuations, deforested areas, and mountain ranges as dispersal barriers. Conservation actions such as restricting long-distance colony transportation, preserving local adaptations, and improving the connectivity between highlands and lowlands are likely to assure future pollination services.

Revisiting metazoan phylogeny with genomic sampling of all phyla.

Proper biological interpretation of a phylogeny can sometimes hinge on the placement of key taxa-or fail when such key taxa are not sampled. In this light, we here present the first attempt to investigate (though not conclusively resolve) animal relationships using genome-scale data from all phyla. Results from the site-heterogeneous CAT + GTR model recapitulate many established major clades, and strongly confirm some recent discoveries, such as a monophyletic Lophophorata, and a sister group relationship between Gnathifera and Chaetognatha, raising continued questions on the nature of the spiralian ancestor. We also explore matrix construction with an eye towards testing specific relationships; this approach uniquely recovers support for Panarthropoda, and shows that Lophotrochozoa (a subclade of Spiralia) can be constructed in strongly conflicting ways using different taxon- and/or orthologue sets. Dayhoff-6 recoding sacrifices information, but can also reveal surprising outcomes, e.g. full support for a clade of Lophophorata and Entoprocta + Cycliophora, a clade of Placozoa + Cnidaria, and raising support for Ctenophora as sister group to the remaining Metazoa, in a manner dependent on the gene and/or taxon sampling of the matrix in question. Future work should test the hypothesis that the few remaining uncertainties in animal phylogeny might reflect violations of the various stationarity assumptions used in contemporary inference methods.

Low-intensity resistance training with partial blood flow restriction and high-intensity resistance training induce similar changes in skeletal muscle transcriptome in elderly humans.

We aimed to investigate the mechanisms underlying muscle growth after 12 weeks of resistance training performed with blood flow restriction (RT-BFR) and high-intensity resistance training (HRT) in older individuals. Participants were allocated into the following groups: HRT, RT-BFR, or a control group. High-throughput transcriptome sequencing was performed by the Illumina HiSeq 2500 platform. HRT and RT-BFR presented similar increases in the quadriceps femoris cross-sectional area, and few genes were differently expressed between interventions. The small differences in gene expression between interventions suggest that similar mechanisms may underpin training-induced muscle growth.

Data from proteomic analysis of bovine Longissimus dorsi muscle associated with intramuscular fat content.

The proteomic data presented in this article are associated with the research article entitled "Longissimus dorsi muscle label-free quantitative proteomic reveals biological mechanisms associated with intramuscular fat deposition" published in Journal of Proteomics [1]. In this article, we characterized the proteomic profile of bovine Longissimus dorsi muscle from Nelore steers and identified differentially abundant proteins associated with the intramuscular fat (IMF) content. An integrated transcriptome-assisted label-free quantitative proteomic approach by High Definition Mass Spectrometry (HDMSE) was employed to identify and quantify the proteins. A functional enrichment analysis using the differentially abundant proteins list was performed to understand the biological processes involved in IMF deposition. Moreover, to explore and clarify the biological mechanisms that influence IMF content, the mRNA data for the same trait from Cesar and collaborators [2] obtained by RNA-sequencing technology was compared with proteomic data. The mRNA data is deposited in the European Nucleotide Archive (ENA) repository (EMBL-EBI), under accession PRJEB13188.

Identification of SNPs potentially related to immune responses and growth performance in Litopenaeus vannamei by RNA-seq analyses.

Litopenaeus vannamei is one of the most important shrimp species for worldwide aquaculture. Despite this, little genomic information is available for this penaeid and other closely related taxonomic crustaceans. Consequently, genes, proteins and their respective polymorphisms are poorly known for these species. In this work, we used the RNA sequencing technology (RNA-seq) in L. vannamei shrimp evaluated for growth performance, and exposed to the White Spot Syndrome Virus (WSSV), in order to investigate the presence of Single Nucleotide Polymorphisms (SNPs) within genes related to innate immunity and growth, both features of great interest for aquaculture activity. We analyzed individuals with higher and lower growth rates; and infected (unhealthy) and non-infected (healthy), after exposure to WSSV. Approximately 7,000 SNPs were detected in the samples evaluated for growth, being 3,186 and 3,978 exclusive for individuals with higher and lower growth rates, respectively. In the animals exposed to WSSV we found about 16,300 unique SNPs, in which 9,338 were specific to non-infected shrimp, and 7,008 were exclusive to individuals infected with WSSV and symptomatic. In total, we describe 4,312 unigenes containing SNPs. About 60% of these unigenes returned GO blastX hits for Biological Process, Molecular Function and Cellular Component ontologies. We identified 512 KEGG unique KOs distributed among 275 pathways, elucidating the majority of metabolism roles related to high protein metabolism, growth and immunity. These polymorphisms are all located in coding regions, and certainly can be applied in further studies involving phenotype expression of complex traits, such as growth and immunity. Overall, the set of variants raised herein enriches the genomic databases available for shrimp, given that SNPs originated from nextgen are still rare for this relevant crustacean group, despite their huge potential of use in genomic selection approaches.