Next Article in Journal
Rewriting Human History and Empowering Indigenous Communities with Genome Editing Tools
Next Article in Special Issue
Spatiotemporal Differentiation of Alpine Butterfly Parnassius glacialis (Papilionidae: Parnassiinae) in China: Evidence from Mitochondrial DNA and Nuclear Single Nucleotide Polymorphisms
Previous Article in Journal
Mitochondrial DNA: A Key Regulator of Anti-Microbial Innate Immunity
Previous Article in Special Issue
Developmental Transcriptomic Analysis of the Cave-Dwelling Crustacean, Asellus aquaticus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 1
10.3390/genes11010087
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider

by
Lindsay S. Miles
1,2,*,
Nadia A. Ayoub
3,
Jessica E. Garb
4,
Robert A. Haney
5 and
Brian C. Verrelli
1
1
Center for Life Sciences Education, Virginia Commonwealth University, Richmond, VA 23284, USA
2
Department of Biology, University of Toronto Mississauga, Mississauga, ON L5L 1C6, Canada
3
Department of Biology, Washington and Lee University, Lexington, VA 24450, USA
4
Department of Biological Sciences, University of Massachusetts Lowell, Lowell, MA 01854, USA
5
Department of Biology, Ball State University, Muncie, IN 47306, USA
*
Author to whom correspondence should be addressed.
Genes 2020, 11(1), 87; https://0-doi-org.brum.beds.ac.uk/10.3390/genes11010087
Submission received: 6 November 2019 / Revised: 6 December 2019 / Accepted: 7 January 2020 / Published: 12 January 2020
(This article belongs to the Special Issue Arthropod Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Due to their abundance and ability to invade diverse environments, many arthropods have become pests of economic and health concern, especially in urban areas. Transcriptomic analyses of arthropod ovaries have provided insight into life history variation and fecundity, yet there are few studies in spiders despite their diversity within arthropods. Here, we generated a de novo ovarian transcriptome from 10 individuals of the western black widow spider (Latrodectus hesperus), a human health pest of high abundance in urban areas, to conduct comparative ovarian transcriptomic analyses. Biological processes enriched for metabolism—specifically purine, and thiamine metabolic pathways linked to oocyte development—were significantly abundant in L. hesperus. Functional and pathway annotations revealed overlap among diverse arachnid ovarian transcriptomes for highly-conserved genes and those linked to fecundity, such as oocyte maturation in vitellogenin and vitelline membrane outer layer proteins, hormones, and hormone receptors required for ovary development, and regulation of fertility-related genes. Comparative studies across arachnids are greatly needed to understand the evolutionary similarities of the spider ovary, and here, the identification of ovarian proteins in L. hesperus provides potential for understanding how increased fecundity is linked to the success of this urban pest.

1. Introduction

Arthropods have long been used as models in studying life history variation and evolutionary adaptation [1,2,3]. Arguably one of the most important life history traits is fecundity, as it is closely tied to fitness, and studies of fecundity and population abundance highlight the arthropod literature [4,5,6]. Arthropods also dominate the literature as invasive and pest species as a result of their high fecundity and ability to rapidly adapt to different environments, which include the urban and agricultural areas that humans have created. A few human pest examples are the common flour beetle, due to their close affinity with grain storage [7,8,9]; aphids, of which over 80 global species destroy crops and ornamentals [10]; and termites that feed on wood in homes [11,12,13]. However, from an ecological and evolutionary perspective, some of the more interesting examples come from those arthropod pest species that are anthropogenic obligates to the extent that they have become threats to human health. For example, cockroaches, bedbugs, and even mosquito species [14,15,16] have adapted to these newly-formed environments, and their abundance implies they thrive within them. In this context, there is a great need to characterize the genotypic and phenotypic responses associated with pest species’ adaptation to human environments to prioritize management and health concerns [17].
With genomic and transcriptomic tools available for any species, comparative and functional analyses of gene expression can be used to investigate fitness-related traits in non-model arthropod pest species in human-altered environments. Ovaries fulfill several important functions tied to oocyte development, hormone secretion, and fertilization, and thus, transcriptomic analyses of these tissues can shed light onto factors driving fecundity. Adrian and Comeron [18] analyzed ovary transcript profiles of early- and late-meiosis development stages in Drosophila melanogaster and showed maternal effects differ across developmental stages. In doing so, they identified new transcripts that likely map to regions of heterochromatin, which is surprising given the very deep genomic sequencing of this model organism. Yang et al. [19] examined the ovary transcriptome in Portunus trituberculatus, the Japanese blue crab and the most widely fished species of crab in the world, and found significant variation in gene expression across sexual maturation stages that can help increase population size for this economically-relevant species. Finally, Uengwetwanita et al. [20] examined the ovary transcriptome in Penaeus monodon—the black tiger shrimp, and again a species of economic importance—to show that although new genes were identified, gene regulation followed similar pathways identified in vertebrates. These studies highlight how comparative transcriptome analyses reveal both unique and similar patterns of gene expression that can be used to understand the ecology and evolution of fecundity in arthropods of interest to humans.
The western black widow spider (Latrodectus hesperus) has become a pest species of human medical relevance due to its potent venom and unique characteristics as an urban adaptor [21,22,23,24]. Specifically, urban L. hesperus show higher fecundity and more dense aggregations than non-urban spiders [25,26,27]. Our previous sampling of thousands of genomewide single nucleotide polymorphisms (SNPs) from 21 urban and non-urban southwestern U.S. locales found higher genetic diversity within urban areas and higher gene flow among urban areas [28], a result that contrasts with the standard hypothesis associated with urbanization [29]. Additionally, our ‘popgraph’ network analysis identified specific urban ‘hubs’ that contribute more to L. hesperus population expansion and gene flow than others, and has implications for pest management [30]. These ecological and genetic patterns raise the question of what phenotypic variation exists for the genes and pathways associated with L. hesperus fecundity, not only between urban environments but also in comparison to other arachnid species. However, despite their unique phylogenetic and taxonomic diversity, there is a pronounced paucity of transcriptome information available overall for spiders, and these studies typically focus on venom and silk [31,32,33,34,35,36,37,38]. In fact, only two ovary studies have examined the biological and functional dynamics of transcript variation associated with fecundity [39,40].
Although tissue-specific transcriptome data exist for L. hesperus, there are no ovary-specific expression analyses and no annotated genome available, without which we struggle to understand factors underlying fecundity variation that drives population expansion in this urban pest. Here, we provide comparative ovarian transcriptomic analyses of the western black widow spider from tissue-specific RNA-seq data collected from multiple individuals from a population sample in the southwestern United States. Our overall objective is to construct a de novo ovarian transcriptome for L. hesperus to compare with other spiders and arthropods in identifying genes and pathways associated with ovarian expression and function. Addressing this objective will contribute to our understanding of the evolutionary diversity of arthropod ovaries in general, and will be specifically important to our understanding of the functional variation associated with fecundity in the western black widow spider as an urban pest model of human health concern.

2. Materials and Methods

2.1. Sample Collection and RNA Extraction

We collected 10 L. hesperus adult female spiders in June, 2017, during their breeding season, from an urban locale within the Phoenix (AZ, USA) city limits. Specifically, individuals were collected from a single area of infestation (GPS coordinates: 33.618972, −112.065231). As is typical of this species within urban areas (Miles et al., 2018a, b; [28,30]), individuals were found within their webs at high densities covering brick walls and rocky outcrops of artificial, man-made surfaces facing a major city street. Ovaries from each individual were dissected in 150 mM sodium chloride, 15 mM sodium citrate using forceps and a dissection microscope, and snap frozen in liquid nitrogen. Total RNA was initially isolated from the ovary tissue using Trizol (Invitrogen, Waltham, MA, USA) and chloroform extraction, then further purified using the RNeasy kit (Qiagen, Hilden, Germany), and any contaminating DNA was removed with on-column DNase treatment (Qiagen). The concentration and purity of the extracted RNA was quantified using an Agilent 2100 Bioanalyzer (Santa Clara, CA, USA).

2.2. Illumina Library Prep and Sequencing

The cDNA libraries were generated with the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA), followed by paired-end, 150 bp sequencing in one lane of HiSeq 4000 (Illumina). Reads were cleaned using Trim Galore! v.0.3.7 (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore) with FastQC v.0.11.2 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) that removed Illumina adaptors and low-quality reads.

2.3. De Novo Assembly of Ovary Transcriptome

There is currently no complete genome available for L hesperus or any closely related species. For this reason, and because our initial inspection of the data revealed very high transcript abundance and coverage compared to former transcriptome datasets of L hesperus, we took this opportunity to generate a de novo ovarian transcriptome. The clean reads were pooled for de novo assembly using default parameters in Trinity v.2.5.1 [41,42]. Benchmarking Universal Single-Copy Orthologs (BUSCO v3) was used with default parameters to assess the completeness of the transcriptome assembly using the ‘eukaryota odb9’ database, which contains 303 highly conserved genes [43,44]. Salmon v.1.55 [45] was used to quantify transcript abundance in the ovary transcriptome using the quasi-mapping algorithm.

2.4. Functional Annotation

Our methods are aligned with our objective to identify genes and pathways expressed in the ovary and compare these with other ovarian data in arachnids. At this time, we did not focus on analyses that attempt to identify those transcripts specific to or enriched in the ovary compared to other tissues. Instead, similar to other arthropod studies, we desired to first generate a complete dataset of genes and pathways expressed in the ovary. Although we expect this current approach will reveal a number of genes that are universally expressed across tissues (i.e., ‘housekeeping’ genes), this approach is vital for our current comparative analyses, as well as for our follow-up studies to investigate genes and pathways integral to western black widow spider ovarian evolution in urban areas. In fact, we expect genes involved in urban adaptation are linked to fecundity and ovarian pathways, but likely are not genes unique to the ovary. Finally, we also did not attempt to identify ‘unique transcripts’. In doing so, it would reduce the number of redundancies in the data, but it also will remove closely-related genes/orthologs (i.e., recent duplication) that are of interest down the road. Thus, here, we outline methods focused on analyses of transcript data initially interpreted as isoforms for an exhaustive dataset at this stage, and we evaluate our results in light of this approach in the Discussion.
The transcript isoforms were mapped to the NCBI nt database using blastn [46] in the Blast v2.7.1+ package (http://blast.ncbi.nlm.nih.gov) with e-value 1 (e-20), max_target_seq 1 and outfmt 6 parameters changed. Next, the transcript isoforms were mapped against the UniProt database to identify functional proteins using Blastx. Finally, excluding isoforms that had not been matched to UniProt, we mapped the remaining isoforms against the only two previously published arachnid ovary-specific transcriptomes, the Spanish funnel-web spider, Macrothele calpeiana [39], and the common house spider, Parasteatoda tepidariorum [40], to identify overlapping arachnid ovary genes.
To determine the transcripts with nearly full-length proteins and their corresponding UniProt BLAST hits, we first used Blast v2.7.1+ to determine the best matching proteins with the UniProt database, then we used a Trinity Perl script to extract the transcripts with >80% coverage of the UniProt proteins. Next, we extracted these transcripts and their corresponding UniProt BLAST hits to be used in the Gene Ontology (GO) analysis to identify genes in the transcriptome with high transcript abundance for the ‘biological process’, ‘cellular component’, and ‘molecular function’ categories in the ovary. The program PANTHER v.14.1 [47] was used to perform the GO analysis. We generated an annotation report using Trinotate v.3.2.0 [48], then uploaded the annotation results to Blast2GO [49] to perform the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, which identify potential functional pathways in the ovary.

3. Results

3.1. Sequencing and De Novo Assembly

For each of the 10 L. hesperus libraries, 35M-63M raw sequence reads were generated, and 98% of clean reads were retained after pre-processing (e.g., adaptor removal, quality trimming, ‘N’ removal). The assembly of transcripts resulted in a total of 690,256 unique transcript isoforms. The BUSCO analysis of the assembled ovary transcriptome resulted in 96% completeness, with 28% complete and single copy (N = 85 of 303 highly conserved genes), 68% complete and duplicate copy (N = 207 of 303), and 4% fragmented (N = 11 of 303). There are 5842 proteins that are represented by nearly full-length transcripts, having >80% alignment coverage and 4021 proteins that are covered by more than 90% of their protein lengths. Using a standard minimum threshold of five ‘transcripts per million transcripts’ (TPM), we found 64,580 transcript isoforms with this minimum threshold of expression in the ovary, and which were used to further characterize genes in our GO and KEGG analyses.

3.2. Functional Annotation and Classification

From the 64,580 transcript isoforms (noted above), there were 37,387 transcripts for the L. hesperus ovary that were assigned GO terms, with each of these transcripts assigned to one or more GO terms in the biological process, cellular component, and molecular function categories (Figure 1). This analysis finds that the GO terms in the L. hesperus ovary transcriptome are highly abundant for ‘biological processes’. Within this category, there were 4495 hits with cellular (33%) and metabolic (28%) processes as the most abundant (Supplementary Materials Figure S1, Table S1). For the cellular component category, there were 3392 hits with cell (49%) and organelle (29%) components as the most abundant. Finally, in the molecular function category, there were 3455 hits with catalytic activity (42%) and binding (33%) function as the most abundant.
The KEGG orthology (KO) assignments resulted in an annotation of 160 pathways that correspond to 14,878 transcripts in the ovary transcriptome (Table S2). The most abundant pathways are related to metabolism (Figure 2), with the two most representative pathways being purine metabolism and thiamine metabolism (Supplementary Materials Figure S2). Purine metabolism plays a role in cell growth or division, with expression in the ovary tied to the cell proliferation of oocytes [50,51]. Thiamine metabolism plays a role in the activity of key enzymes in cellular metabolism and has been linked to the meiotic maturation of oocytes [52,53].

3.3. Spider Ovary Genes of Interest

The UniProt BLAST identified 551 transcripts whose highest match was listed as an arthropod-annotated gene. Within these 551 top hits, 28 aligned to arachnid genes and 7 of these arachnid genes were associated with α-latrotoxins, which are toxin components found in Latrodectus spider venom [54]. Additionally, our BLAST search against M. calpeiana and P. tepidariorum ovarian transcriptomes resulted in an overlap of 905 transcripts among all three species, 440 shared between L. hesperus and M. calpeiana, and 82,310 shared between L. hesperus and P. tepidariorum (Figure 3). However, our L. hesperus transcriptome had an additional 80,184 transcripts that match to UniProt genes, but which did not match to transcripts in the other two arachnid ovary transcriptomes.
Although the vast majority of overlapping transcripts are of ‘housekeeping’ genes found in all cell types, Table 1 shows 15 genes identified in our transcriptome that overlapped with the 19 previously-identified fecundity-related genes in arthropods [55,56,57,58,59,60,61,62,63,64,65]. These genes include vitellogenin and vitelline membrane outer layer protein, which are involved in oocyte development [66]. Several hormones and hormone receptors are required for development in the ovary including 3-β-hydroxysteroid dehydrogenase, mandibular organ-inhibiting hormone, lutropin-choriogonadotropic hormone receptor, follicle-stimulating hormone receptor, and estrogen [67]. Lastly, genes that regulate vitellogenin and other fertility-related genes in the ovary include zinc-finger protein, phosphoglycerate kinase, carboxylesterase, protein geranylgeranyl transferase, SRY related HMG-Box-11, C-terminal-binding protein, fizzy, sex-lethal, serine-threonine kinase receptor, integrin protein, rab guanosine diphosphate/guanosine triphosphate exchange factor, and maternal protein pumilio [65]. Interestingly, only seven of these 15 fecundity-related genes (3 β-hydroxysteroid dehydrogenase, estrogen, zinc-finger protein, phosphoglycerate kinase, protein geranylgeranyl transferase, SRY related HMG-Box-11, C-terminal-binding protein) overlapped with the P. tepidariorum ovarian transcriptome, whereas none of them overlapped with M. calpeiana.

4. Discussion

The overall objective of this study was to construct a de novo ovarian transcriptome for L. hesperus to compare with other spiders and arthropods in identifying genes and pathways associated with ovarian expression and function. Because no genome is yet available for L. hesperus, a de novo assembly of the transcriptome was generated from a large dataset with transcript coverage. Specifically, our annotation rate was 54.5%, which is similar to other arthropod studies, even though these previous studies have an available genome. For example, the annotation rate is ~50% for several arthropods including Antheraea pernyi (50.8%), Helicoverpa armigera (50.8%), Helicoverpa assulta (54.0%), Spodoptera frugiperda (51.1%), and Athetis lepigone (41.5%) [68,69,70]. With only two spider ovarian transcriptomes assembled, L. hesperus has a significantly higher annotation rate than M. calpeiana (29.4%) [39] and P. tepidariorum (19.7%) [40]. These results reflect a successful de novo assembly and provide for comparisons with ovarian transcriptomes of other spiders.

4.1. Western Black Widow Spider Ovarian Gene and Pathway Categories

The main function of the ovary in spiders, and arthropods in general, is to develop mature oocytes and produce hormones needed for female reproduction. Thus, the identification of genes and pathways with transcripts of high abundance tied to biological processes involved in metabolism and cellular proliferation in the ovary is expected, and is a pattern that is consistent with other arthropod transcriptome studies [39,71,72,73]. Sub-groups of cellular processes and metabolic processes have the first and the second highest transcript count, respectively, for GO term annotations (Figure 1), which is consistent with other arthropod studies of ovarian transcriptomes.
While metabolic processes in general have been found as top KEGG pathways in other arthropods [39,72,74,75], this is the first study to identify purine and thiamine metabolism as the top pathways in an arthropod ovarian transcriptome. Purine metabolism has been well-documented as a key factor in cell growth and division during oocyte proliferation in mice [50,51]. Purine pathways are important in both follicular fluid and serum surrounding oocytes, which aid in the development of oocytes [76]. If purine metabolism is reduced, oocyte development is suppressed [76,77,78]. Purine metabolism is linked to oocyte quality, which can lead to lower fertility when there are changes to this metabolic pathway [51,76]. Similarly, thiamine metabolism has been identified as part of cellular metabolism and meiotic maturation of oocytes in mice [52,53]. Arthropods that are unable to sequester thiamine in their diets, such as tsetse flies and ticks, have reduced fertility when their symbiont—typically a bacteria that produces thiamine—is removed [79]. The reduced fertility in these arthropods is likely caused by vitamin B deficiencies that can inhibit meiotic maturation of oocytes [53]. While this study suggests these pathways may be major factors in ovarian development, investigating these patterns (such as vitamin deficiencies) in other closely-related spiders and within L. hesperus will determine whether they play a unique role in this species.

4.2. Species Comparisons of Spider Ovarian Transcriptomes

Although only two other spider ovary-specific transcriptome datasets are available, there are already a few interesting lessons to be learned here. We find that there is significantly more overlap among L. hesperus and P. tepidariorum than between L. hesperus and M. calpeiana. In fact, there is more overlap among all three transcriptomes than there is shared between L. hesperus and M. calpeiana. There are several potential explanations for these patterns, the first being the different datasets. Similar to our transcriptome, the P. tepidariorum data are generated from the ovaries of adult female spiders, and these data include all raw sequences, i.e., sequenced transcripts that blast to known and unknown genes [40]. The M. calpeiana data are also generated from the ovaries of adult female spiders, but these data are the result of a ‘subtractive’ approach, i.e., an attempt by the authors to enrich for tissue-specific transcripts at various steps [39]. As noted previously, this approach was specifically not taken by us to help guard against a bias towards ovary-specific transcripts, which may have less overlap with other spider and arthropod transcriptomes, as we speculate here.
Given the approach by Frias-Lopez [39], the overlap between L. hesperus and M. calpeiana may more reflect ‘ovary-specific’ transcripts in contrast to the L. hesperus and P. tepidariorum overlap. However, inspection of this overlap between L. hesperus and M. calpeiana finds the vast majority are standard housekeeping genes. In addition, we might expect the majority of genes in common with ovary development and oocyte production should be shared among the three species. However, over 98% of the three-species overlapping transcripts also map to housekeeping genes with >90% identity, and none of the arthropod fecundity-related genes exist in this three-species overlap. Instead, of the 19 fecundity-related genes, only 7 are shared between L. hesperus and P. tepidariorum and none are shared between L. hesperus and M. calpeiana. One consideration is that these genes are expressed differentially due to evolutionary changes in oocyte development. It is possible that if gene expression across oocyte developmental stages is examined (as in Drosophila) [18], that we find that these genes conserved in arthropods have diverse roles in egg maturation even among spiders.
Another explanation for these patterns in overlap among spider species is that L. hesperus and P. tepidariorum (from the Theridiidae family), are more phylogenetically related [80,81], whereas, the divergence time of M. calpeiana (Macrothelidae family) from the other two is estimated at over 300 million years [82]. Thus, given this significantly higher divergence, we expect that overlap among the three will largely reflect highly evolutionarily conserved genes, which is in fact what we see. This result alone does not imply that these housekeeping transcripts do not have ‘ovary-specific’ functions, especially given the curation of the M. calpeiana dataset; however, more analyses would be needed to test this hypothesis. For example, we may speculate that a proportion of the transcripts do not overlap because they have sufficiently diverged due to positive selection as a response to ovary-specific adaptation—a hypothesis we will address with molecular evolutionary analyses of these genes within and between species. These results make the obvious argument that more ovarian transcriptomes are needed from closely-related spiders to determine whether these annotated transcripts with no overlap are truly species-specific.

4.3. Western Black Widow Spider Genes Associated with Fecundity

Our analyses of the L. hesperus transcriptome identified 15 of 19 genes relevant for fecundity and the process of oocyte growth previously noted in arthropod ovaries [55,56,57,58,59,60,61,62,63,64,65]. These studies had identified serine-threonine kinase receptor, integrin protein, rab guanosine diphosphate/guanosine triphosphate exchange factor, and maternal protein pumillo as involved in arthropod fertility and vitellogenin regulation [65]. However, these four genes were not identified in our transcriptome. As noted previously, the lack of these transcripts in our dataset may less reflect the loss of the genes in L. hesperus, but instead that there is some temporal stage of expression or population variance. For example, there were several vitellogenin regulation genes that we did identify, and the lack of these four genes (listed above) in our transcriptome may indicate that not all of them are required to be expressed during oocyte development. Indeed, the hormonal control of vitellogenesis in spiders is not well understood, and there is a paucity of data on the hormones involved [83]. Finally, as previously noted, it is also possible that these genes were not identified because they have sufficiently diverged due to positive selection compared to other arthropod ovarian genes. There are a number of transcripts in L. hesperus that were annotated but do not overlap with the two species, and thus, our planned evolutionary analyses of transcripts initially annotated as ‘species-specific’ can address this hypothesis.
For each of the fecundity-related genes identified in L. hesperus, we identified evidence consistent with alternative splicing (Table 1). In fact, the most abundant represented gene is the zinc-finger protein, with 315,938 transcripts aligned. The zinc-finger protein is an abundant transcription factor in many transcriptomes, and its primary role is to regulate gene expression [84,85]. One of the most ubiquitous genes present in arthropod ovaries that is likely regulated by these zinc-finger transcription factors is vitellogenin, which is used in the process of oocyte development [66]. Our results indicate that these genes are also undergoing alternative splicing, as we find 184 transcripts of the vitellogenin gene, a pattern found in arthropods, and more specifically in other arachnids [19,86,87,88,89]. Vitellogenesis involves the synthesis and accumulation of proteins for eggs to develop [66]. Mated female spiders have ovaries in vitellogenesis and post-vitellogenesis phases, which would produce significant levels of vitellogenin proteins [90,91,92,93]. As our current study includes only gravid females, it would be of interest to examine different female stages (i.e., unmated, multiple matings, etc.) to detect potential adaptive expression variation in fecundity.

4.4. Fecundity and Human Pest Adaptation to Urban Environments

The western black widow spider is an example of a human pest that has thrived due to human-facilitated movement [28,29,30], and our identification of ovarian transcriptome patterns here has implications for its adaptation to urban environments. One interesting result was the identification of many isoforms associated with α-latrotoxin, which is a member of a family of highly-expressed neurotoxins found in Latrodectus spider venom that is specific to vertebrates [54]. Haney et al. [40] also identified these transcripts in their ovarian transcriptome for P. tepidariorum, and they were found here in the overlap with L. hesperus, but not between L. hesperus and M. calpeiana. Toxin-related transcripts have been identified in L. tredecimguttatus eggs, to which the authors speculated the eggs have their own toxic mechanism [94]. The ovaries collected in this study were from gravid females, with visible oocytes, and therefore these toxins may be produced in the yolk. Our sample was collected from an urban geographic area of Phoenix, AZ, and urban spiders have higher fecundity than non-urban spiders [25,26,27]. Thus, this additional toxin in the yolk may play a role in higher fitness of the urban egg sacs relative to non-urban ones. Future studies can test this hypothesis by measuring differential gene expression of these latrotoxins in ovaries between urban and non-urban populations to elucidate if this pattern reflects a unique response to the urban environment.
Although vitellogenin transcripts in the ovarian transcriptome are expected for an arthropod, it is unclear whether we may expect high expression of these genes to be associated with increased fecundity. Previous studies have suggested that urban L. hesperus produce significantly lower ‘quality’ eggs compared to non-urban spiders [25,26,95]; however, urban spiders have been shown to produce significantly more eggs per egg sac. Thus, it may be the case that selection in urban areas produces smaller, but more eggs compared to non-urban areas as a result of reduced interspecific competition and increased intraspecific competition [25,26,95]. In fact, the pattern seen here of genes linked to cell proliferation and vitellogenin and its regulation, could be consistent with higher fecundity in urban areas being adaptive. Future work can target these genes of interest between urban and non-urban populations to tease apart quality and quantity differences in reproductive success.
Much of what we have learned about how urban adapter pests perform and thrive in the urban environment is due to contrasts between populations, species, and phenotypes. This current study highlights the value and gaps in tissue-specific and comparative species analyses. For example, while we need more species transcriptomes to resolve whether annotated transcripts are species-specific, we also need more tissue-specific transcriptomes to determine whether genes and pathways reflect tissue-specific function. In fact, other than fecundity, there are other urban adapter phenotypes that co-evolve and contribute to urban adaptation. For example, previous studies have identified prey differences between urban and non-urban populations [27,96] that likely affect gene expression in food metabolism, venom production, and prey-capturing silk. Therefore, the analysis of transcriptome patterns across multiple tissues associated with these phenotypes across both urban and non-urban populations would reveal the genetic architecture involved in the adaptation of pest species, such as the western black widow spider, to human-altered environments.

Supplementary Materials

The following are available online: https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/1/87/s1, Table S1: GO terms and the number of L. hesperus ovarian transcripts per GO category; Table S2: KEGG pathway annotations for L. hesperus ovarian transcriptome; Figure S1: Gene Ontology (GO) level distribution of transcripts in the L. hesperus ovarian transcriptome, Figure S2: Purine and thiamine metabolism pathways in the L. hesperus ovarian transcriptome.

Author Contributions

Conceptualization, all authors; Formal analysis, L.S.M.; Funding acquisition, B.C.V.; Methodology, all authors; Project administration, B.C.V.; Supervision, B.C.V.; Writing—original draft, L.S.M. and B.C.V.; Writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

B.C.V. was partially supported by a National Science Foundation Research Coordination Network grant no. DEB-1840663.

Acknowledgments

We thank M.S. Rosenberg for assistance on data curation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lewontin, R.C. The Effects of Population Density and Composition on Viability in Drosophila melanogaster. Evolution 1955, 9, 27–41. [Google Scholar] [CrossRef]
  2. Dawson, P.S. Life History Strategy and Evolutionary History of Tribolium Flour Beetles. Evolution 1977, 31, 226–229. [Google Scholar] [CrossRef] [PubMed]
  3. Scriber, J.M.; Slansky, F. The nutritional ecology of immature insects. Annu. Rev. Entomol. 1981, 26, 183–211. [Google Scholar] [CrossRef]
  4. Leather, S.R. Size, Reproductive Potential and Fecundity in Insects: Things aren’t as Simple as They Seem. Oikos 1988, 51, 386–389. [Google Scholar] [CrossRef]
  5. Barnes, A.I.; Wigby, S.; Boone, J.M.; Partridge, L.; Chapman, T. Feeding, fecundity and lifespan in female Drosophila melanogaster. Proc. R. Soc. B Biol. Sci. 2008, 275, 1675–1683. [Google Scholar] [CrossRef] [Green Version]
  6. Lister, B.C.; Garcia, A. Climate-driven declines in arthropod abundance restructure a rainforest food web. Proc. Natl. Acad. Sci. USA 2018, 115, E10397–E10406. [Google Scholar] [CrossRef] [Green Version]
  7. Li, B.; Predel, R.; Neupert, S.; Hauser, F.; Tanaka, Y.; Cazzamali, G.; Williamson, M.; Arakane, Y.; Verleyen, P.; Schoofs, L.; et al. Genomics, transcriptomics, and peptidomics of neuropeptides and protein hormones in the red flour beetle Tribolium castaneum. Genome Res. 2008, 18, 113–122. [Google Scholar] [CrossRef] [Green Version]
  8. Tan, A.; Palli, S.R. Edysone receptor isoforms play distinct roles in controlling molting and metamorphosis in the red flour beetle, Tribolium castaneum. Mol. Cell. Endocrinol. 2008, 291, 42–49. [Google Scholar] [CrossRef] [Green Version]
  9. Shukla, J.N.; Palli, S.R. Doublesex target genes in the red flour beetle, Tribolium castaneum. Sci. Rep. 2012, 2, 948. [Google Scholar] [CrossRef] [Green Version]
  10. Van Emden, H.F.; Harrington, R. Aphids as Crop Pests; Cabi: Wallingford, UK, 2017. [Google Scholar]
  11. Constantino, R. The pest termites of South America: Taxonomy, distribution and status. J. Appl. Entomol. 2002, 126, 355–365. [Google Scholar] [CrossRef]
  12. Kumar, D.; Pardeshi, M. Biodiversity of termites in agro-ecosystem and relation between their niche breadth and pest status. J. Entomol. 2011, 8, 250–258. [Google Scholar] [CrossRef] [Green Version]
  13. Evans, T.A.; Iqbal, N. Termite (order Blattodea, infraorder Isoptera) baiting 20 years after commercial release. Pest Manag. Sci. 2015, 71, 897–906. [Google Scholar] [CrossRef] [PubMed]
  14. Byrne, K.; Nichols, R.A. Culex pipiens in London underground tunnels: Differentiation between surface and subterranean populations. Heredity 1999, 82, 7–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Vargo, E.L.; Crissman, J.R.; Booth, W.; Santangelo, R.G.; Mukha, D.V.; Schal, C. Hierarchical genetic analysis of German cockroach (Blattella germanica) populations from within buildings to across continents. PLoS ONE 2014, 9, e102321. [Google Scholar] [CrossRef] [Green Version]
  16. Booth, W.; Balvín, O.; Vargo, E.L.; Vilímová, J.; Schal, C. Host association drives genetic divergence in the bed bug, Cimex lectularius. Mol. Ecol. 2015, 24, 980–992. [Google Scholar] [CrossRef]
  17. Rivkin, L.R.; Santangelo, J.S.; Alberti, M.; Aronson, M.F.J.; de Keyzer, C.W.; Diamond, S.E.; Fortin, M.J.; Frazee, L.J.; Gorton, A.J.; Hendry, A.P.; et al. A roadmap for urban evolutionary ecology. Evol. Appl. 2019, 12, 384–398. [Google Scholar] [CrossRef] [Green Version]
  18. Adrian, A.B.; Comeron, J.M. The Drosophila early ovarian transcriptome provides insight to the molecular causes of recombination rate variation across genomes. BMC Genom. 2013, 14, 794. [Google Scholar] [CrossRef] [Green Version]
  19. Yang, Y.; Wang, J.; Han, T.; Liu, T.; Wang, C.; Xiao, J.; Mu, C.; Li, R.; Yu, F.; Shi, H. Ovarian transcriptome analysis of Portunus trituberculatus provides insights into genes expressed during phase III and IV development. PLoS ONE 2015, 10, e0138862. [Google Scholar] [CrossRef] [Green Version]
  20. Uengwetwanit, T.; Ponza, P.; Sangsrakru, D.; Wichadakul, D.; Ingsriswang, S.; Leelatanawit, R.; Klinbunga, S.; Tangphatsornruang, S.; Karoonuthaisiri, N. Transcriptome-based discovery of pathways and genes related to reproduction of the black tiger shrimp (Penaeus monodon). Mar. Genom. 2018, 37, 69–73. [Google Scholar] [CrossRef]
  21. Blair, R.B. Land use and avian species diversity along an urban gradient. Ecol. Appl. 1996, 6, 506–519. [Google Scholar] [CrossRef]
  22. McKinney, M.L.; Lockwood, J.L. Biotic homogenization: A few winners replacing many losers in the next mass extinction. Trends Ecol. Evol. 1999, 14, 450–453. [Google Scholar] [CrossRef]
  23. McIntyre, N.E.; Rango, J.; Fagan, W.F.; Faeth, S.H. Ground arthropod community structure in a heterogeneous urban environment. Landsc. Urban Plan. 2001, 52, 257–274. [Google Scholar] [CrossRef]
  24. Kark, S.; Iwaniuk, A.; Schalimtzek, A.; Banker, E. Living in the city: Can anyone become an ‘urban exploiter’? J. Biogeogr. 2007, 34, 638–651. [Google Scholar] [CrossRef]
  25. Johnson, J.C.; Trubl, P.J.; Miles, L.S. Black widows in an urban desert: City-living compromises spider fecundity and egg investment despite urban prey abundance. Am. Midl. Nat. 2012, 168, 333–341. [Google Scholar] [CrossRef]
  26. Johnson, J.C.; Miles, L.S.; Trubl, P.J.; Hagenmaier, A. Maternal effects on egg investment and offspring performance in black widow spiders. Anim. Behav. 2014, 91, 67–73. [Google Scholar] [CrossRef]
  27. Trubl, P.; Gburek, T.; Miles, L.; Johnson, J.C. Black widow spiders in an urban desert: Population variation in an arthropod pest across metropolitan Phoenix, AZ. Urban Ecosyst. 2012, 15, 599–609. [Google Scholar] [CrossRef]
  28. Miles, L.S.; Johnson, J.C.; Dyer, R.J.; Verrelli, B.C. Urbanization as a facilitator of gene flow in a human health pest. Mol. Ecol. 2018, 27, 3219–3230. [Google Scholar] [CrossRef]
  29. Miles, L.S.; Rivkin, L.R.; Johnson, M.T.J.; Munshi-South, J.; Verrelli, B.C. Gene flow and genetic drift in urban environments. Mol. Ecol. 2019, 28, 4138–4151. [Google Scholar] [CrossRef]
  30. Miles, L.S.; Dyer, R.J.; Verrelli, B.C. Urban hubs of connectivity: Contrasting patterns of gene flow within and among cities in the western black widow spider. Proc. R. Soc. B Biol. Sci. 2018, 285, 20181224. [Google Scholar] [CrossRef]
  31. Mattila, T.M.; Bechsgaard, J.S.; Hansen, T.T.; Schierup, M.H.; Bilde, T. Orthologous genes identified by transcriptome sequencing in the spider genus Stegodyphus. BMC Genom. 2012, 13, 70. [Google Scholar] [CrossRef] [Green Version]
  32. Posnien, N.; Zeng, V.; Schwager, E.E.; Pechmann, M.; Hilbrant, M.; Keefe, J.D.; Damen, W.G.; Prpic, N.M.; McGregor, A.P.; Extavour, C.G. A comprehensive reference transcriptome resource for the common house spider Parasteatoda tepidariorum. PLoS ONE 2014, 9, e104885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Sanggaard, K.W.; Bechsgaard, J.S.; Fang, X.; Duan, J.; Dyrlund, T.F.; Gupta, V.; Jiang, X.; Cheng, L.; Fan, D.; Feng, Y.; et al. Spider genomes provide insight into composition and evolution of venom and silk. Nat. Commun. 2014, 5, 3765. [Google Scholar] [CrossRef] [PubMed]
  34. Clarke, T.H.; Garb, J.E.; Hayashi, C.Y.; Haney, R.A.; Lancaster, A.K.; Corbett, S.; Ayoub, N.A. Multi-tissue transcriptomics of the black widow spider reveals expansions, co-options, and functional processes of the silk gland gene toolkit. BMC Genom. 2014, 15, 365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Clarke, T.H.; Garb, J.E.; Hayashi, C.Y.; Arensburger, P.; Ayoub, N.A. Spider transcriptomes identify ancient large-scale gene duplication event potentially important in silk gland evolution. Genome Biol. Evol. 2015, 7, 1856–1870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Clarke, T.H.; Garb, J.E.; Haney, R.A.; Chaw, R.C.; Hayashi, C.Y.; Ayoub, N.A. Evolutionary shifts in gene expression decoupled from gene duplication across functionally distinct spider silk glands. Sci. Rep. 2017, 7, 8393. [Google Scholar] [CrossRef] [PubMed]
  37. Haney, R.A.; Clarke, T.H.; Gadgil, R.; Fitzpatrick, R.; Hayashi, C.Y.; Ayoub, N.A.; Garb, J.E. Effects of gene duplication, positive selection, and shifts in gene expression on the evolution of the venom gland transcriptome in widow spiders. Genome Biol. Evol. 2016, 8, 228–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Garb, J.E.; Haney, R.A.; Schwager, E.E.; Gregorič, M.; Kuntner, M.; Agnarsson, I.; Blackledge, T.A. The transcriptome of Darwin’s bark spider silk glands predicts proteins contributing to dragline silk toughness. Commun. Biol. 2019, 2, 275. [Google Scholar] [CrossRef]
  39. Frías-López, C.; Almeida, F.C.; Guirao-Rico, S.; Vizueta, J.; Sánchez-Gracia, A.; Arnedo, M.A.; Rozas, J. Comparative analysis of tissue-specific transcriptomes in the funnel-web spider Macrothele calpeiana (Araneae, Hexathelidae). PeerJ 2015, 3, e1064. [Google Scholar] [CrossRef] [Green Version]
  40. Haney, R.A.; Garb, J.E.; Forsyth, F.S.; Matte, T. Alternative transcription at venom genes and its role as a complementary mechanism for the generation of venom complexity. Front. Ecol. Evol. 2019, 7, 85. [Google Scholar] [CrossRef] [Green Version]
  41. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644. [Google Scholar] [CrossRef] [Green Version]
  42. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494. [Google Scholar] [CrossRef] [PubMed]
  43. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [Green Version]
  44. Waterhouse, R.M.; Seppey, M.; Simão, F.A.; Manni, M.; Ioannidis, P.; Klioutchnikov, G.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 2017, 35, 543–548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Mi, H.; Huang, X.; Muruganujan, A.; Tang, H.; Mills, C.; Kang, D.; Thomas, P.D. PANTHER version 11: Expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2016, 45, D183–D189. [Google Scholar] [CrossRef] [Green Version]
  48. Bryant, D.M.; Johnson, K.; DiTommaso, T.; Tickle, T.; Couger, M.B.; Payzin-Dogru, D.; Lee, T.J.; Leigh, N.D.; Kuo, T.H.; Davis, F.G.; et al. A tissue-mapped axolotl de novo transcriptome enables identification of limb regeneration factors. Cell Rep. 2017, 18, 762–776. [Google Scholar] [CrossRef] [Green Version]
  49. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef]
  50. Downs, S.M.; Eppig, J.J. Induction of mouse oocyte maturation in vivo by perturbants of purine metabolism. Biol. Reprod. 1987, 36, 431–437. [Google Scholar] [CrossRef]
  51. Downs, S.M. Purine control of mouse oocyte maturation: Evidence that nonmetabolized hypoxanthine maintains meiotic arrest. Mol. Reprod. Dev. 1993, 35, 82–94. [Google Scholar] [CrossRef]
  52. Gillett, J.D. Variation in the time of release of the ovarian development hormone in Aedes aegypti. Nature 1957, 180, 656–657. [Google Scholar] [CrossRef] [PubMed]
  53. Tsuji, A.; Nakamura, T.; Shibata, K. Effects of mild and severe vitamin B1 deficiencies on the meiotic maturation of mice oocytes. Nutr. Metab. Insights 2017, 10. [Google Scholar] [CrossRef] [Green Version]
  54. Haney, R.A.; Ayoub, N.; Clarke, T.H.; Hayashi, C.Y.; Garb, J.E. Dramatic expansion of the black widow toxin arsenal uncovered by multi-tissue transcriptomics and venom proteomics. BMC Genom. 2014, 15, 366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Mikhailov, A.T.; Torrado, M. Carboxylesterases moonlight in the male reproductive tract: A functional shift pivotal for male fertility. Front. Biosci. 2000, 5, E53–E62. [Google Scholar] [PubMed] [Green Version]
  56. Wen, J.; Xie, J.; Liu, S.; Gui, J. Differential expression and characterization analysis of a new gene with WD domains in fish oogenesis. Sci. China C Life Sci. 2001, 44, 541–553. [Google Scholar] [CrossRef] [PubMed]
  57. Danshina, P.V.; Geyer, C.B.; Dai, Q.; Goulding, E.H.; Willis, W.D.; Kitto, G.B.; McCarrey, J.R.; Eddy, E.M.; O’Brien, D.A. Phosphoglycerate kinase 2 (PGK2) is essential for sperm function and male fertility in mice. Biol. Reprod. 2010, 82, 136–145. [Google Scholar] [CrossRef] [Green Version]
  58. Hashiyama, K.; Hayashi, Y.; Kobayashi, S. Drosophila sex lethal gene initiates female development in germline progenitors. Science 2011, 333, 885–888. [Google Scholar] [CrossRef]
  59. Van Doren, M. Determining sexual identity. Science 2011, 333, 829–830. [Google Scholar] [CrossRef]
  60. Wei, L.; Cheng, D.; Li, D.; Meng, M.; Peng, L.; Tang, L.; Pan, M.; Xiang, Z.; Xia, Q.; Lu, C. Identification and characterization of Sox genes in the silkworm, Bombyx mori. Mol. Biol. Rep. 2011, 38, 3573–3584. [Google Scholar] [CrossRef]
  61. Jiang, T.; Hou, C.C.; She, Z.Y.; Yang, W.X. The SOX gene family: Function and regulation in testis determination and male fertility maintenance. Mol. Biol. Rep. 2013, 40, 2187–2194. [Google Scholar] [CrossRef]
  62. Zhai, Y.; Zhang, J.; Sun, Z.; Dong, X.; He, Y.; Kang, K.; Liu, Z.; Zhang, W. Proteomic and transcriptomic analyses of fecundity in the brown planthopper Nilaparvata lugens (Stal). J. Proteome Res. 2013, 12, 5199–5212. [Google Scholar] [CrossRef] [PubMed]
  63. Zhai, Y.; Sun, Z.; Zhang, J.; Kang, K.; Chen, J.; Zhang, W. Activation of the TOR signalling pathway by glutamine regulates insect fecundity. Sci. Rep. 2015, 5, 10694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Sun, Z.X.; Zhai, Y.F.; Zhang, J.Q.; Kang, K.; Cai, J.H.; Fu, Y.; Qiu, J.Q.; Shen, J.W.; Zhang, W.Q. The genetic basis of population fecundity prediction across multiple field populations of Nilaparvata lugens. Mol. Ecol. 2015, 24, 771–784. [Google Scholar] [CrossRef] [PubMed]
  65. Qiu, J.; He, Y.; Zhang, J.; Kang, K.; Li, T.; Zhang, W. Discovery and functional identification of fecundity-related genes in the brown planthopper by large-scale RNA interference. Insect Mol. Biol. 2016, 25, 724–733. [Google Scholar] [CrossRef] [PubMed]
  66. Tufail, M.; Takeda, M. Molecular characteristics of insect vitellogenins. J. Insect Physiol. 2008, 54, 1447–1458. [Google Scholar] [CrossRef]
  67. Bateman, A.; Bennett, H.P. Granulins: The structure and function of an emerging family of growth factors. J. Endocrinol. 1998, 158, 145–151. [Google Scholar] [CrossRef] [Green Version]
  68. Li, L.T.; Zhu, Y.B.; Ma, J.F.; Li, Z.Y.; Dong, Z.P. An analysis of the Athetis lepigone transcriptome from four developmental stages. PLoS ONE 2013, 8, e73911. [Google Scholar] [CrossRef]
  69. Do Nascimento, A.R.; Fresia, P.; Cônsoli, F.L. Comparative transcriptome analysis of lufenuron-resistant and susceptible strains of Spodoptera frugiperda (Lepidoptera: Noctuidae). BMC Genom. 2015, 16, 985. [Google Scholar] [CrossRef] [Green Version]
  70. Zhang, J.; Wang, B.; Dong, S.; Cao, D.; Dong, J.; Walker, W.B.; Liu, Y.; Wang, G. Antennal transcriptome analysis and comparison of chemosensory gene families in two closely related noctuidae moths, Helicoverpa armigera and H. assulta. PLoS ONE 2015, 10, e0117054. [Google Scholar] [CrossRef]
  71. Rendón-Anaya, M.; Delaye, L.; Possani, L.D.; Herrera-Estrella, A. Global transcriptome analysis of the scorpion Centruroides noxius: New toxin families and evolutionary insights from an ancestral scorpion species. PLoS ONE 2012, 7, e43331. [Google Scholar] [CrossRef] [Green Version]
  72. Vizueta, J.; Frías-López, C.; Macías-Hernández, N.; Arnedo, M.A.; Sánchez-Gracia, A.; Rozas, J. Evolution of chemosensory gene families in arthropods: Insight from the first inclusive comparative transcriptome analysis across spider appendages. Genome Biol. Evol. 2017, 9, 178–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Bian, H.X.; Chen, D.B.; Zheng, X.X.; Ma, H.F.; Li, Y.P.; Li, Q.; Xia, R.X.; Wang, H.; Jiang, Y.R.; Liu, Y.Q.; et al. Transcriptomic analysis of the prothoracic gland from two lepidopteran insects, domesticated silkmoth Bombyx mori and wild silkmoth Antheraea pernyi. Sci. Rep. 2019, 9, 5313. [Google Scholar] [CrossRef] [PubMed]
  74. Li, C.; Weng, S.; Chen, Y.; Yu, X.; Lü, L.; Zhang, H.; He, J.; Xu, X. Analysis of Litopenaeus vannamei transcriptome using the next-generation DNA sequencing technique. PLoS ONE 2012, 7, e47442. [Google Scholar] [CrossRef] [PubMed]
  75. Kim, H.S.; Lee, B.Y.; Won, E.J.; Han, J.; Hwang, D.S.; Park, H.G.; Lee, J.S. Identification of xenobiotic biodegradation and metabolism-related genes in the copepod Tigriopus japonicus whole transcriptome analysis. Mar. Genom. 2015, 24, 207–208. [Google Scholar] [CrossRef] [PubMed]
  76. Chen, M.; Zhang, B.; Cai, S.; Zeng, X.; Ye, Q.; Mao, X.; Zhang, S.; Zeng, X.; Ye, C.; Qiao, S. Metabolic disorder of amino acids, fatty acids and purines reflects the decreases in oocyte quality and potential in sows. J. Proteom. 2019, 200, 134–143. [Google Scholar] [CrossRef]
  77. Sirard, M.A.; First, N.L. In vitro inhibition of oocyte nuclear maturation in the bovine. Biol. Reprod. 1988, 39, 229–234. [Google Scholar] [CrossRef]
  78. Carta, M.C.; Mattana, A.; Camici, M.; Allegrini, S.; Tozzi, M.G.; Sgarrella, F. Catabolism of exogenous deoxyinosine in cultured epithelial amniotic cells. Biochim. Biophys. Acta (BBA) Gen. Subj. 2001, 1528, 74–80. [Google Scholar] [CrossRef]
  79. Rio, R.V.; Attardo, G.M.; Weiss, B.L. Grandeur alliances: Symbiont metabolic integration and obligate arthropod hematophagy. Trends Parasitol. 2016, 32, 739–749. [Google Scholar] [CrossRef] [Green Version]
  80. Hormiga, G.; Griswold, C.E. Systematics, phylogeny, and evolution of orb-weaving spiders. Annu. Rev. Entomol. 2014, 59, 487–512. [Google Scholar] [CrossRef]
  81. Gendreau, K.L.; Haney, R.A.; Schwager, E.E.; Wierschin, T.; Stanke, M.; Richards, S.; Garb, J.E. House spider genome uncovers evolutionary shifts in the diversity and expression of black widow venom proteins associated with extreme toxicity. BMC Genom. 2017, 18, 178. [Google Scholar] [CrossRef] [Green Version]
  82. Garrison, N.L.; Rodriguez, J.; Agnarsson, I.; Coddington, J.A.; Griswold, C.E.; Hamilton, C.A.; Hedin, M.; Kocot, K.M.; Ledford, J.M.; Bond, J.E. Spider phylogenomics: Untangling the Spider Tree of Life. PeerJ 2016, 4, e1719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Sawadro, M.; Bednarek, A.; Babczyńska, A. The current state of knowledge on the neuroactive compounds that affect the development, mating and reproduction of spiders (Araneae) compared to insects. Invertebr. Neurosci. 2017, 2017, 17. [Google Scholar] [CrossRef]
  84. Nolte, R.T.; Conlin, R.M.; Harrison, S.C.; Brown, R.S. Differing roles for zinc fingers in DNA recognition: Structure of a six-finger transcription factor IIIA complex. Proc. Natl. Acad. Sci. USA 1998, 95, 2938–2943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Emerson, R.O.; Thomas, J.H. Adaptive evolution in zinc finger transcription factors. PLoS Genet. 2009, 5, e1000325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Thompson, D.M.; Khalil, S.M.S.; Jeffers, L.A.; Sonenshine, D.E.; Mitchel, R.D.; Osgood, C.J.; Michael, R.R. Sequence and the developmental and tissue-specific regulation of the first complete vitellogenin messenger RNA from ticks responsible for heme sequestration. Insect Biochem. Mol. Biol. 2007, 37, 363–374. [Google Scholar] [CrossRef] [PubMed]
  87. Wang, W.; Wu, X.; Liu, Z.; Zheng, H.; Cheng, Y. Insights into hepatopancreatic functions for nutrition metabolism and ovarian development in the crab Portunus trituberculatus: Gene discovery in the comparative transcriptome of different hepatopancreas stages. PLoS ONE 2014, 9, e84921. [Google Scholar] [CrossRef] [PubMed]
  88. Guo, J.; Wang, L.; Wu, H.; Cao, Y.; Xiao, R.; Lai, X.; Liu, J.; Yi, J.; Zhang, G. Molecular characterization and expression of vitellogenin genes from the wolf spider Pardosa pseudoannulata (Araneae: Lycosidae). Physiol. Entomol. 2018, 43, 295–305. [Google Scholar] [CrossRef]
  89. Bednarek, A.W.; Sawadro, M.K.; Nicewicz, Ł.; Babczyńska, A.I. Vitellogenins in the spider Parasteatoda tepidariorum–expression profile and putative hormonal regulation of vitellogenesis. BMC Dev. Biol. 2019, 19, 4. [Google Scholar] [CrossRef] [Green Version]
  90. Jean André, C.R. The ultrastructure of the vitelline body in the oocyte of the spider Tegenaria parietina. J. Biophys. Biochem. Cytol. 1957, 3, 977–1000. [Google Scholar] [CrossRef]
  91. Sotello, J.R.; Oscar, T.-C. Electron microscope study of the vitelline body of some spider oocytes. J. Biophys. Biochem. Cytol. 1957, 3, 301–317. [Google Scholar] [CrossRef] [Green Version]
  92. Osaki, H. Electron microscope studies on developing oocytes of the spider Plexippus paykulli. Annot. Zool. Jpn. 1972, 45, 187–200. [Google Scholar]
  93. Trabalon, M.; Bautz, A.M.; Moriniere, M.; Porcheron, P. Ovarian development and correlated changes in hemolymphatic ecdysteroid levels in two spiders, Coelotes terrestris and Tegenaria domestics (araneae, agelenidae). Gen. Comp. Endocrinol. 1992, 88, 128–136. [Google Scholar] [CrossRef]
  94. Yan, S.; Wang, X. Recent advances in research on widow spider venoms and toxins. Toxins 2015, 7, 5055–5067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Johnson, J.C.; Urcuyo, J.; Moen, C.; Stevens, D.R. Urban heat island conditions experienced by the Western black widow spider (Latrodectus hesperus): Extreme heat slows development but results in behavioral accommodations. PLoS ONE 2019, 14, e0220153. [Google Scholar] [CrossRef] [PubMed]
  96. Trubl, P.; Johnson, J.C. Ecological stoichiometry of the black widow spider and its prey from desert, urban and laboratory populations. J. Arid Environ. 2019, 163, 18–25. [Google Scholar] [CrossRef]
Genes 11 00087 g001 550
Figure 1. Gene Ontology (GO) analysis of the Latrodectus hesperus ovarian transcriptome according to the identified transcripts involvement in biological process, cellular component, and molecular function categories.
Figure 1. Gene Ontology (GO) analysis of the Latrodectus hesperus ovarian transcriptome according to the identified transcripts involvement in biological process, cellular component, and molecular function categories.
Genes 11 00087 g001
Genes 11 00087 g002 550
Figure 2. Kyoto Encyclopedia of Genes and Genomes (KEGG) Top 20 pathways for the L. hesperus ovarian transcriptome.
Figure 2. Kyoto Encyclopedia of Genes and Genomes (KEGG) Top 20 pathways for the L. hesperus ovarian transcriptome.
Genes 11 00087 g002
Genes 11 00087 g003 550
Figure 3. Transcripts overlapping in the ovarian transcriptomes among the three arachnids L. hesperus, Macrothele calpeia na, and Parasteatoda tepidariorum.
Figure 3. Transcripts overlapping in the ovarian transcriptomes among the three arachnids L. hesperus, Macrothele calpeia na, and Parasteatoda tepidariorum.
Genes 11 00087 g003
Table
Table 1. Arthropod fecundity-related genes present in the L. hesperus ovarian transcriptome.
Table 1. Arthropod fecundity-related genes present in the L. hesperus ovarian transcriptome.
Gene.FunctionCount *
Vitellogeninoocyte development181
Vitelline membrane outer layer proteinoocyte development174
Estrogenovary development3397
3 beta-hydroxysteroid dehydrogenase (3-β-HSD)ovary development405
Mandibular organ-inhibiting hormone (MOIH)ovary development100
Lutropin-choriogonadotropic hormone receptor (LSHR)ovary development30
Follicle-stimulating hormone receptor (FSHR)ovary development5
Zinc-Finger protein (ZFP)vitellogenin regulation315938
phosphoglycerate kinasevitellogenin regulation5424
carboxylesterasevitellogenin regulation2610
C-terminal-binding protein (CtBP)vitellogenin regulation584
protein geranylgeranyl transferasevitellogenin regulation297
fizzy (fzy)vitellogenin regulation244
SRY related HMG-Box-11 (SOX-11)vitellogenin regulation145
Sex-lethal (Sxl)vitellogenin regulation5
* The number of reads identified as the given name.

Share and Cite

MDPI and ACS Style

Miles, L.S.; Ayoub, N.A.; Garb, J.E.; Haney, R.A.; Verrelli, B.C. Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider. Genes 2020, 11, 87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11010087

AMA Style

Miles LS, Ayoub NA, Garb JE, Haney RA, Verrelli BC. Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider. Genes. 2020; 11(1):87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11010087

Chicago/Turabian Style

Miles, Lindsay S., Nadia A. Ayoub, Jessica E. Garb, Robert A. Haney, and Brian C. Verrelli. 2020. "Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider" Genes 11, no. 1: 87. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11010087

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop