Skip to main content
  • Facebook
  • Twitter
  • YouTube
  • LinkedIn
  • Google Plus
  • Other GSA Resources
    • Genetics Society of America
    • Genetics
    • Genes to Genomes: The GSA Blog
    • GSA Conferences
    • GeneticsCareers.org
  • Log in
G3: Genes | Genomes | Genetics

Main menu

  • HOME
  • ISSUES
    • Current Issue
    • Early Online
    • Archive
  • ABOUT
    • About the journal
    • Why publish with us?
    • Editorial board
    • Contact us
  • SERIES
    • All Series
    • Genomic Prediction
    • Multiparental Populations
  • ARTICLE TYPES
    • About Article Types
    • Genome Reports
    • Meeting Reports
    • Mutant Screen Reports
    • Software and Data Resources
  • PUBLISH & REVIEW
    • Scope & publication policies
    • Submission & review process
    • Article types
    • Prepare your manuscript
    • Submit your manuscript
    • After acceptance
    • Guidelines for reviewers
  • SUBSCRIBE
    • Email alerts
    • RSS feeds
  • Other GSA Resources
    • Genetics Society of America
    • Genetics
    • Genes to Genomes: The GSA Blog
    • GSA Conferences
    • GeneticsCareers.org

User menu

Search

  • Advanced search
G3: Genes | Genomes | Genetics

Advanced Search

  • HOME
  • ISSUES
    • Current Issue
    • Early Online
    • Archive
  • ABOUT
    • About the journal
    • Why publish with us?
    • Editorial board
    • Contact us
  • SERIES
    • All Series
    • Genomic Prediction
    • Multiparental Populations
  • ARTICLE TYPES
    • About Article Types
    • Genome Reports
    • Meeting Reports
    • Mutant Screen Reports
    • Software and Data Resources
  • PUBLISH & REVIEW
    • Scope & publication policies
    • Submission & review process
    • Article types
    • Prepare your manuscript
    • Submit your manuscript
    • After acceptance
    • Guidelines for reviewers
  • SUBSCRIBE
    • Email alerts
    • RSS feeds
Previous ArticleNext Article

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma

View ORCID ProfilePhilipp Brand, Nicholas Saleh, Hailin Pan, Cai Li, Karen M. Kapheim and Santiago R. Ramírez
G3: Genes, Genomes, Genetics September 1, 2017 vol. 7 no. 9 2891-2898; https://doi.org/10.1534/g3.117.043687
Philipp Brand
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Philipp Brand
  • For correspondence: pbrand@ucdavis.edu
Nicholas Saleh
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hailin Pan
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cai Li
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karen M. Kapheim
Department of Biology, Utah State University, Logan, Utah 84322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Santiago R. Ramírez
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
Loading

Abstract

Bees provide indispensable pollination services to both agricultural crops and wild plant populations, and several species of bees have become important models for the study of learning and memory, plant–insect interactions, and social behavior. Orchid bees (Apidae: Euglossini) are especially important to the fields of pollination ecology, evolution, and species conservation. Here we report the nuclear and mitochondrial genome sequences of the orchid bee Euglossa dilemma Bembé & Eltz. E. dilemma was selected because it is widely distributed, highly abundant, and it was recently naturalized in the southeastern United States. We provide a high-quality assembly of the 3.3 Gb genome, and an official gene set of 15,904 gene annotations. We find high conservation of gene synteny with the honey bee throughout 80 MY of divergence time. This genomic resource represents the first draft genome of the orchid bee genus Euglossa, and the first draft orchid bee mitochondrial genome, thus representing a valuable resource to the research community.

  • whole-genome assembly
  • corbiculate bee
  • orchid bee
  • eusocial
  • mitochondrial genome
  • Genome Report

Bees (Apoidea) are important models for the study of learning and memory (Menzel and Muller 1996), plant–insect interactions (Doetterl and Vereecken 2010), and the evolution of social behavior (Nowak et al. 2010; Woodard et al. 2011; Kapheim et al. 2015). Among the >20,000 bee species worldwide, lineages have evolved varied degrees of specialization on floral resources such as pollen, resins, and oils (Wcislo and Cane 1996; Michener 2007; Litman et al. 2011). These relationships are wide-ranging and have substantial impact on the health and function of natural and agricultural systems (Klein et al. 2007). Furthermore, several transitions from an ancestral solitary to a derived eusocial behavior have occurred within bees (Danforth 2002; Cardinal and Danforth 2011; Branstetter et al. 2017). Thus, bees provide unique opportunities to investigate the genetic underpinnings of multiple major ecological and evolutionary transitions. The repeated evolution of different behavioral phenotypes in bees, including foraging and social behavior, provides a natural experiment that allows for the determination of general as well as species-specific molecular genomic changes underlying phenotypic transitions. In order to capitalize on this potential, whole-genome sequences of a divergent array of bee species with different life histories are needed (Kapheim et al. 2015).

Orchid bees (Apidae; Euglossini) are among the most important pollinators of angiosperms in the neotropical region (Ramírez et al. 2002). While female orchid bees collect nectar, pollen, and resin for nest construction and brood-cell provisioning, male bees collect perfume compounds from floral and nonfloral sources (Vogel 1966; Whitten et al. 1993; Eltz et al. 1999; Roubik and Hanson 2004). These volatile compounds are used to concoct a species-specific perfume blend that is subsequently used during courtship display, presumably to attract conspecific females. This unique male scent-collecting behavior has recently been examined in a broad array of molecular ecological and evolutionary studies, focusing on phenotypic evolution, chemical communication, plant–insect mutualisms, and speciation (Eltz et al. 2008, 2011; Ramírez et al. 2011; Brand et al. 2015; Weber et al. 2016).

While most of the ∼220 species of orchid bees appear to be solitary, several species have transitioned to living in coordinated social groups (Garófalo 1985; Pech et al. 2008; Augusto and Garófalo 2009). Female Euglossa dilemma individuals, for example, can either form solitary nests and provision their own brood cells, or live in small groups where daughters remain in their natal nest and help their mother rear her offspring, instead of dispersing to found their own nest. The social E. dilemma nests (similar to the closely related E. viridissima; Pech et al. 2008) exhibit true division of labor, with subordinate daughters foraging for resources and the reproductively dominant mothers remaining in the nest and laying eggs (N. Saleh and S.R. Ramírez, personal observation). This facultative eusocial behavior represents an early stage in social evolution and makes E. dilemma well suited for studying the genetic mechanisms underlying the transition from solitary to eusocial behavior. While other facultative eusocial species evolved throughout the bee lineage, orchid bees have a unique taxonomic position (Cardinal and Danforth 2011). Orchid bees are part of the corbiculate bees, together with the honey bees, bumble bees, and stingless bees, three obligately eusocial bee lineages (Figure 1A). As the sister group to all other corbiculate bee lineages (Romiguier et al. 2016; Branstetter et al. 2017; Peters et al. 2017), orchid bees may provide key insights into the early stages of eusociality and the possible evolutionary trajectories that led to the obligate eusocial behavior exhibited by honey bees.

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Genomic features. (A) Phylogeny of the four corbiculate bee tribes with orchid bees as sistergroup to honey bees, stingless bees, and bumblebees (Romiguer et al. 2016). (B) K-mer distribution spectrum (k = 25) of genomic sequence reads. The positively skewed spectrum reveals a high abundance of a few k-mers, leading to an estimate of 87.7% repetitiveness of the E. dilemma genome. Red shows the k-mer spectrum before, and blue after error correction. (C) Genomic element density including genic and nongenic features as a fraction of the overall genome assembly length excluding stretches of N. 49.15% of the assembly could not be annotated with the selected methods. (D) Synteny between the E. dilemma and the honey bee (Apis mellifera) genome. In an analysis including E. dilemma scaffolds of ≥100 kb length, 83% showed ≥95% synteny to a single honey bee scaffold. Photographs in (A) are reproduced from Wikimedia under the CC BY-SA 3.0 license.

Here we present the draft genome of the orchid bee species E. dilemma. Using a combined paired-end and mate-pair library sequencing approach, we assembled 18% of the predicted 3.3 Gb genome, and annotated a high-quality gene set including 15,904 genes. In addition, we reconstructed three quarters of the mitochondrial genome with the help of transcriptome data, representing the first orchid bee mitogenome. These genomic resources will facilitate the genetic study of outstanding ecological and evolutionary questions, such as the evolution of resource preferences and the evolution of eusociality. Moreover, it provides an important genomic resource for a group of crucial neotropical pollinators, which are of specific concern for conservation biologists (Zimmermann et al. 2011; Suni and Brosi 2012; Soro et al. 2016; Suni 2016).

Materials and Methods

Genome sequencing and assembly

Nuclear genome:

Sequencing of the E. dilemma genome was based on six haploid male individuals collected at Fern Forest Nature Center in Broward County, FL (26°13’46.3’’N, 80°11’08.9’’W) in February 2011. This population was chosen based on its low nucleotide diversity resulting from a bottleneck during a single introduction to Southern Florida ∼15 yr ago (Skov and Wiley 2005; Pemberton and Wheeler 2006; Zimmermann et al. 2011). DNA was extracted from each bee independently and used for the construction of four paired-end (two 170 and 500 bp libraries, respectively) and four mate-pair (two 2 and 5 kb libraries, respectively) sequencing libraries. Next, the paired-end libraries were sequenced in 90 cycles and the mate-pair libraries in 49 cycles on an Illumina HiSeq2000. The resulting sequence data were run through fastuniq v1.1 (Xu et al. 2012) to remove polymerase chain reaction (PCR) duplicates and quality trimmed using trim_galore v0.3.7 (Babraham Bioinformatics). Subsequently, reads were used for de novo assembly with ALLPATHS-LG v51750 (Gnerre et al. 2011) and Soap-denovo2 (Luo et al. 2012) with varying settings. Gaps within scaffolds were closed using GapCloser v1.12 (Luo et al. 2012) for each assembly. ALLPATHS-LG with default settings resulted in the highest-quality assembly, based on assessments of annotation completeness (see below). This assembly (E. dilemma genome assembly v1.0) was used for all subsequent analyses. All other assemblies were excluded from analysis, but are available upon request.

The preprocessed reads were used for k-mer based genome size estimates. We used ALLPATHS-LG to produce and analyze the k-mer frequency spectrum (k = 25). Genome size was estimated on the basis of the consecutive length of all reads divided by the overall sequencing depth as (N × (L − K + 1) − B)/D = G, where N is the total number of reads, L is the single-read length, K is the k-mer length, B is the total count of low-frequency (frequency ≤3) k-mers that are most likely due to sequencing errors, D is the k-mer depth estimated from the k-mer frequency spectrum, and G is the genome size. In addition, we used the ALLPATHS-LG k-mer frequency spectrum to predict the repetitive fraction of the genome.

The quality of the genome assembly was assessed using standard N statistics and assembly completeness as measured by the CEGMA v2.5 (Parra et al. 2007) and BUSCO v1.1 (Simão et al. 2015) pipelines. CEGMA was run in default mode, whereas BUSCO was run with the arthropoda_odb9 OrthoDB database (Zdobnov et al. 2017) in genome mode.

We estimated the mean per-base genome coverage on the basis of the preprocessed reads and the estimated genome size as Embedded Image where R is the number of reads and L the mean read length of sequence library i, G is the estimated genome size, and C the resulting per-base coverage.

Mitogenome:

Initial attempts to reconstruct the mitochondrial genome from our whole-genome shotgun sequencing reads including read subsampling and exclusion of rare variants were only partially successful, due to high sequence variability of sequencing reads with similarity to mitochondrial loci (data not shown). In addition, we have observed that the amplification of mitochondrial DNA in standard PCR leads to a high level of polymorphic sites in E. dilemma and other orchid bees (P. Brand and S.R. Ramírez personal observation). Together, this suggests the presence of nuclear copies of the mitochondrial genome (NUMTs) that interfere with the assembly process and PCR amplification. Mitochondrial genes are expressed in almost every tissue in eukaryotes. We used this feature to reconstruct the mitochondrial genome as far as possible from RNA-Seq data. Therefore, we used available E. dilemma transcriptome assemblies in order to reconstruct the mitochondrial genome from cDNA (Brand et al. 2015). In order to find expressed mitochondrial genes represented in the transcriptome assembly of Brand et al. (2015), we used blastx with the honey bee mitochondrial genome as query (Crozier and Crozier 1993) and an E-value cutoff of 10E−12 (Altschul et al. 1990; Camacho et al. 2009). The contigs and scaffolds that were detected with this approach were annotated following Dietz et al. (2016). Briefly, we performed tblastn and tblastx searches with protein-coding genes and rRNA genes of the honey bee mitochondrial genome, respectively. All hits were used for manual gene annotation using Geneious v8.0.5 (Biomatters Ltd. 2012). Since the recovered mitochondrial mRNA scaffolds contained >1 gene, we searched and annotated intergenic tRNAs using ARWEN 1.2.3 (Laslett and Canbäck 2008) and tRNAscan-SE 1.21 (Lowe and Eddy 1997).

Genome annotation

Gene annotation:

Genes were annotated based on sequence homology and de novo gene predictions. The homology approach was based on the recently updated high-quality official gene set of the honey bee (OGS v3.2; Elsik et al. 2014). All honey bee original gene set (OGS) proteins were used in initial tblastn searches against all E. dilemma scaffolds with an E-value cutoff of 10E−4. All honey bee protein sequences with a blast hit to the E. dilemma genome assembly covering at least 50% of the protein were selected for homology-based annotation. The resulting set of honey bee proteins was used as input to exonerate v2.42.1 (Slater and Birney 2005) in order to annotate homologous open reading frames (ORFs) through accurate exon–intron boundary prediction for each scaffold. Exonerate was run with default settings and the minimum fraction of the possible optimal similarity per protein query set to 35%. In a second round, genes not annotated under the previous settings were rerun with minimum similarity set to 15%. In the case of overlapping annotations on the same strand (i.e., identical ORF orientation) resulting from honey bee queries with high similarity, we discarded all but one annotation with the best exonerate score (based on completeness and similarity). This approach proved feasible due to the relatedness of E. dilemma and the honey bee, as well as the high quality of the well curated honey bee OGS. For de novo gene prediction we used Augustus (Stanke et al. 2008) and SNAP (Korf 2004) trained on the honey bee, with the E. dilemma genome masked for repetitive regions (see below) as input. Only genes predicted by both programs were taken into account. Gene predictions with ≥85% sequence similarity to each other were discarded, to prevent the inclusion of putative unmasked transposable element derived genes in the official gene set. De novo predictions were added to the E. dilemma OGS if not annotated by the homology-based approach.

Repetitive element annotation:

Repetitive elements including tandem repeats, nuclear copies of the mitochondrial genome (NUMTs), and transposable elements (TEs) were annotated using multiple methods.

Tandem repeats:

We searched for micro- and mini-satellites (1–6 and 7–1000 bp motif length, respectively) in all scaffolds using Phobos 3.3.12 (Mayer 2010). We performed two independent runs for each class of tandem repeats with Phobos parameter settings following Leese et al. (2012) (gap score and mismatch score set to −4 and a minimum repeat score of 12; Leese et al. 2012).

NUMTs:

We annotated NUMTs using blastn runs with the partial mitochondrial genome (see above) as query and an E-value cutoff of 10E−4 as used in NUMT analyses of other insect genomes (Pamilo et al. 2007). This approach allowed us to find NUMTs with medium to high similarity to the actual transcriptome-based mitochondrial genome.

TEs:

In order to annotate TEs, we first ran RepeatModeler for de novo repeat element annotation and classification based on the genome annotation followed by RepeatMasker in order to detect the total fraction of repetitive elements present in the genome (Smit et al. 2016). In addition we used the k-mer based de novo repeat assembler REPdenovo (Chu et al. 2016) to identify highly repetitive genomic elements based on the generated short sequence reads. These two methods independently assess repetitive sequence content in genomic data using different approaches thus allowing for a more robust estimation of repetitive genome content. The incorporation of multiple de novo repeat detection pipelines was necessary due to the lack of a bee-specific repeat database. RepeatModeler v1.0.8 was run with default settings using the NCBI blast algorithm (Altschul et al. 1990) for repeat element detection. We used the resulting de novo repeat element annotations as a database for RepeatMasker v4.0.5 with Crossmatch v0.990329 as the search engine. We ran the analysis in sensitive mode in order to identify the total fraction of repeat elements in the assembly. We excluded low complexity regions and small RNA from the analysis (settings -nolow and -norna). REPdenovo v1.0 was run on all preprocessed genomic short reads in default mode with a minimum repeat frequency of 400× (i.e., the squared mean genome coverage of 20×; see Results and Discussion section). We then used Bowtie v2 (Langmead and Salzberg 2012) in the sensitive local alignment mode to map all reads to the resulting contigs. The mapping results were then used to calculate mean per-contig coverage with bamtools (Barnett et al. 2011). In order to estimate the fraction of the overall genome that corresponds to these highly abundant contigs, we divided the mean contig coverage by the respective contig length and the mean genome-wide coverage. The sum of the resulting normalized base pair counts divided by the estimated genome size was then used as an estimate of the fraction of the genome containing highly abundant sequences.

Genome structure

To analyze genome structure, we compared the genome-wide gene synteny of E. dilemma and the honey bee. We used the genomic locations of homologous genes (as determined above) of the honey bee and E. dilemma scaffolds of at least 100 kb length to build haplotype blocks with a minimum length of 1 kb. Haplotype blocks included the entire gene span as well as intergenic regions whenever two or more adjacent genes were homologous in both species. We discarded gene annotations from downstream analysis that were recovered as homologous to multiple genomic locations in either species. Furthermore, we excluded E. dilemma genes that were recovered as homologous to honey bee scaffolds belonging to unknown linkage groups.

Data availability

The E. dilemma genome assembly Edil_v1.0, the annotation, and the original gene set Edil_OGS_v1.0 are available for download via NCBI (Bioproject: PRJNA388474), Beebase (Elsik et al. 2016), the i5k NAL workspace (https://i5k.nal.usda.gov/euglossa-dilemma) (i5K Consortium 2013), and the Ramirez Lab website. The raw reads are available via NCBI (Bioproject: PRJNA388474). The published raw transcriptome sequence reads are available at the NCBI Sequence Read Archive (SRA: SRX765918) (Brand et al. 2015).

Results and Discussion

Whole-genome assembly

The E. dilemma genome assembly resulted in 22,698 scaffolds with an N50 scaffold length of 144 Kb and a total length of 588 Mb (Table 1). This represents 18% of the k-mer based estimated genome size of 3.2 Gb. Of all sequence reads, 68% aligned to the genome assembly, of which 56% aligned more than once. Further, the k-mer frequency spectrum based on all sequencing reads was strongly positively skewed indicating the presence of highly repetitive sequences in the read set (Figure 1B). Based on the k-mer frequency spectrum, 87.7% of the genome was estimated to be repetitive. This suggests that the genome of E. dilemma consists largely of highly repetitive sequences, explaining the low consecutive assembly length, the high assembly fragmentation, and the high fraction of sequence reads mapping multiple times to the assembly. The mean per-base coverage was estimated to be comparatively low in comparison to previous bee genome assemblies, with 19.7× based on the preprocessed reads and estimated genome size (Kocher et al. 2013, 95.65× coverage; Kapheim et al. 2015, 120–272.3× coverage). Total genomic GC content was 39.9%, and thus similar to previously sequenced bee genomes ranging between 32.7 and 41.5% (Table 1) (Kocher et al. 2013; Elsik et al. 2014; Kapheim et al. 2015).

View this table:
  • View inline
  • View popup
Table 1 E. dilemma genome assembly statistics in comparison to previously published bee genomes

Despite the fragmentation of the genome assembly representing <20% of the estimated genome size, CEGMA analysis revealed complete assemblies of 231 out of 248 core eukaryotic genes (93.2% completeness). Similarly, BUSCO analysis revealed that 1007 out of 1066 highly conserved arthropod genes were completely assembled (94.4% completeness). The BUSCO analysis detected the duplication of a fraction of 4.8% of the benchmark single-copy ortholog genes in the genome, which is similar to the fraction of arthropod single-copy orthologs found to be duplicated in the honey bee genome (6.9%, Weinstock et al. 2006).

Our gene prediction approach generated a comprehensive official gene set including 15,904 protein-coding genes (Table 1). Of these gene models, 11,139 were derived from homology-based predictions, representing 73% of the 15,314 honey bee genes used for annotation. These annotations are well within or exceeding previous bee genome assemblies and are similar to those reported for the other orchid bee genomes available (Table 1) (Kocher et al. 2013; Elsik et al. 2014; Kapheim et al. 2015; Park et al. 2015; Sadd et al. 2015).

The CEGMA and BUSCO analysis and the gene annotation results suggest that the gene-coding fraction of the E. dilemma genome was properly assembled, despite the large estimated genome size and comparatively low per-base sequencing coverage. However, genetic material obtained from natural populations as in our study can lead to the fragmentation of assemblies due to high nucleotide diversity (Kajitani et al. 2014). In addition, high genetic diversity in the underlying genetic material can lead to a high number of false duplicates due to multiple incorporation of divergent genomic regions in genome assemblies (Kelley and Salzberg 2010). Our BUSCO analysis suggests that the assembly did not produce an unusually high fraction of duplicated benchmark single-copy orthologs, indicating a relatively low abundance of false duplicates. The observed fragmentation in our assembly is thus likely to be primarily the result of repetitive genomic elements, and less likely the result of low coverage or high nucleotide diversity in the genetic material used for sequencing.

Overall, our results suggest that our approach was sufficient to produce a high-quality official gene set. The homology-based approach we used resulted in the majority of annotated genes in the official gene set with a known homology to honey bee genes (Supplemental Material, Table S1). This genomic resource will facilitate genome-wide expression studies including gene ontology analyses and comparative gene set analyses among insects.

Mitochondrial genome assembly

The recently published transcriptome assembly used for the reconstruction of the mitochondrial genome contained four scaffolds between 1222 and 4188 bp long with a total consecutive length of 11,128 bp (Figure 2). This corresponds to ∼75% of the estimated length of the mitochondrial genome, based on other corbiculate bee species (Crozier and Crozier 1993; Cha et al. 2007). The E. dilemma mitogenome fragments contained 5 out of 22 tRNAs, 11 out of 13 protein-coding genes of which two were only partially recovered, and the 16S rRNA gene. Within scaffolds all genes showed the known hymenopteran gene order and orientation, while the orientation of the five tRNAs detected was identical to those in the honey bee (Crozier and Crozier 1993; Cha et al. 2007). Attempts to complete the mitochondrial genome using the nuclear genome assembly yielded no improvement of the assembly (data not shown). Accordingly, the mitochondrial genome is entirely derived from transcriptomic sequences.

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

Mitochondrial genome reconstruction. The structure of the honey bee mitochondrial genome and information of the homologous reconstructed parts of the E. dilemma mitochondrial genome. Nonreconstructed parts of incompletely reconstructed genes are hatched.

The high success in mitochondrial gene reconstruction is likely due to the nature of the analyzed transcriptome data. Short intergenic regions as well as polycistronic mitochondrial mRNA likely lead to the assembly of multiple genes into single scaffolds. The A-T rich region is completely missing as well as the ND2 and 12S rRNA genes flanking the region in insect mitogenomes. This unrecovered region also contains a high number of tRNAs in the honey bee, which could explain the low number of recovered tRNAs in E. dilemma. While the partial mitochondrial genome assembly is only 75% complete, it represents the first mitogenome for the group of orchid bees and will thus be a valuable resource for future phylogenetic analyses within the lineage and between more distantly related bee taxa.

Repetitive elements

Tandem repeats:

We detected 76,001 microsatellite loci with a consecutive length of 2,291,067 bp. Mini-satellites with motif lengths from 7 to 1000 bp were less numerous in the genome (67,323 loci), totaling 13,343,515 bp. Accordingly, tandem repeats represent 3.86% of the genome assembly, suggesting that they contribute only a small proportion to the overall genome size (Figure 1C).

NUMTs:

We detected fragments with similarity to the draft mitochondrial genome on 129 scaffolds totaling a length of 150,670 bp. The fragments had a mean length of 764.8 bp and a mean similarity of 91.5% to the mitogenome. This suggests that these fragments are not derived from the mitochondrial genome and represent actual NUMTs. A total of 39 scaffolds carried multiple fragments with high similarity to the mitogenome with a concatenated length of up to 6566 bp, suggesting that respective NUMTs might have originated from larger fragments of the mitogenome. In total, only 0.04% of the whole-genome assembly had hits to the mitogenome (Figure 1C). This is likely an underestimate, due to the incompleteness of the reconstructed mitochondrial genome. Nevertheless, NUMTs likely represent only a small fraction of the whole nuclear E. dilemma genome. Previous analyses have shown a high density of NUMTs in the honey bee in comparison to other insect genomes totaling ∼0.1% of the overall genome size (Pamilo et al. 2007). Accordingly, given the NUMT content detected in E. dilemma, it is possible that a comparatively high NUMT density is a common feature of corbiculate bee genomes.

TEs:

The RepeatModeler analysis revealed a total of 566 repeat element families in the genome assembly of which 142 (25.1%) belonged to known TE families including 106 DNA transposons and 36 retroelements, while the remaining 424 (74.9%) repeat element families could not be classified into known TE families (Table 2 and Table S2). Using the 566 newly detected repeat element families as the input database for RepeatMasker, we annotated a total of 597,369 elements in the genome assembly of which 74,513 (12.5%) were derived from the 142 classified TE families. The remaining 522,856 (87.5%) elements belonged to the unclassified repeat element families. In total, all annotated repeat elements had a cumulative length of 163,384,833 bp corresponding to 38.48% of the total genome assembly. The majority was derived from unclassified (i.e., unknown) repeat element families corresponding to a total of 32.6% of the genome assembly (Figure 1C). Similarly, the read-based de novo repeat analysis in REPdenovo revealed 27,636 contigs derived from k-mers with a minimum coverage of 400× (i.e., the squared mean genome-wide coverage), which includes 831,433,228 bp (26%) of the estimated 3.3 Gb genome.

View this table:
  • View inline
  • View popup
Table 2 Transposable element repeat class analysis

The detected high fraction of the genome associated with repetitive element families in E. dilemma is not surprising given that large genome sizes are often associated with elevated TE activity and TE content. Similar patterns have been observed in the genomes of diverse lineages, ranging from unicellular eukaryotes to complex multicellular organisms like plants, invertebrates and vertebrates (Kidwell 2002). However, TEs are fast evolving and highly specific to their host lineages, which leads to large underestimates of genomic TE content in previously unstudied lineages (Chalopin et al. 2015; Platt et al. 2016). This likely explains the large fraction of unclassified repetitive element families that we detected in our genome assembly. Further, the remaining high fraction of unknown genome content (49.15%, Figure 1C) may have resulted from undetected repetitive elements with no similarity to elements from other genomes sequenced previously. The only publicly accessible bee repeat element families are derived from the honey bee, a species with a comparatively small genome (0.23 Gb) and low TE diversity and content (Weinstock et al. 2006; Kapheim et al. 2015). Our efforts to annotate TEs based on known honey bee elements did not improve the TE annotation (data not shown). Overall, our analysis suggests that a large proportion of the E. dilemma genome is repetitive. This result is similar to the results obtained in the orchid bee Eufriesea mexicana, which has an estimated repetitive genome content of 31% (Kapheim et al. 2015).

Genome structure

Of the 22,698 E. dilemma scaffolds, 580 were at least 100 kb in length and used for synteny analysis with the honey bee genome. A total of 356 of these scaffolds carried at least one gene annotation with known homology to the honey bee, and 329 of these E. dilemma scaffolds were homologous to honey bee scaffolds with known linkage group (LG) association (Table S1). Of these scaffolds, 272 (83%) showed ≥95% syntenic homology (Figure 1D). Overall, the detected syntenic linkage blocks cover 222 MB of scaffold length with homology to the honey bee, representing 85% of the 329 filtered scaffolds. This suggests that the genomic architecture is very similar between E. dilemma and the honey bee, representing a high level of conservation during the 80 MY since the two lineages diverged. Further, our results support a recent comparative analysis of the honey bee and the bumblebee genomes, which revealed high conservation of genomic synteny (Stolle et al. 2011). In comparison, in previous studies gene synteny was found to be less conserved in other insect groups. Extensive local shuffling of gene order beginning on the time scale of 20–40 MY evolutionary distance was described for dipterans, moths, and ants resulting in 60–70% genome-wide synteny after ∼60 MY of divergence time in flies and ants (Clark et al. 2007; d’Alençon et al. 2010; Obbard et al. 2012; Simola et al. 2013; Neafsey et al. 2015; Nygaard et al. 2016). Together, these results support a general pattern of slow evolution of gene synteny in corbiculate bees, independent of the fraction of repetitive genome content.

Conclusion

The genome assembly of the orchid bee E. dilemma that we present here is of high quality, despite its large genome size (estimated to be 3.3 Gb). The 15,904 gene annotations provide a comprehensive set of genes with known homology to the honey bee, facilitating future gene ontology and functional genomic analyses. While we were unable to annotate the mostly repetitive majority of the genome assembly with our approach, the provided sequence reads will be useful for future analyses of repetitive genetic elements in the genome. The nuclear and mitochondrial draft genomes represent a valuable genomic resource for the community of bee geneticists. This genomic resource will likely prove valuable in genetic and functional genomic analyses dealing with the ecology, evolution, and conservation of orchid bees. Furthermore, the genome of the facultatively eusocial E. dilemma will be helpful in the study of the evolution of eusociality, due to its taxonomic placement as the sister lineage to the three obligately eusocial corbiculate bee tribes including stingless bees, bumblebees, and honey bees.

Acknowledgments

We would like to thank Yannick Wurm and one anonymous reviewer for comments that improved the manuscript. We would like to thank Gene E. Robinson, Guojie Zhang, and the Beijing Genomics Institute for financial support for genome sequencing. P.B. was supported by a fellowship of the German academic exchange service (Deutscher Akademischer Austauschdienst, DAAD) for parts of the project. S.R.R. received support from the National Science Foundation (DEB-1457753). This work used the Vincent J. Coates Genomics Sequencing Laboratory at University of California, Berkeley, supported by NIH S10 Instrumentation Grants S10RR029668, S10RR027303, and OD018174.

Footnotes

  • Supplemental material is available online at www.g3journal.org/lookup/suppl/doi:10.1534/g3.117.043687/-/DC1.

  • Communicating editor: B. Oliver

  • Received May 24, 2017.
  • Accepted July 10, 2017.
  • Copyright © 2017 Brand et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Literature Cited

  1. ↵
    1. Altschul S. F.,
    2. Gish W.,
    3. Miller W.,
    4. Myers E. W.,
    5. Lipman D. J.
    , 1990 Basic local alignment search tool. J. Mol. Biol. 215: 403–410.
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Augusto S. C.,
    2. Garófalo C. A.
    , 2009 Bionomics and sociological aspects of Euglossa fimbriata (Apidae, Euglossini). Genet. Mol. Res. 8: 525–538.
    OpenUrlPubMed
  3. ↵
    1. Barnett D. W.,
    2. Garrison E. K.,
    3. Quinlan A. R.,
    4. Strömberg M. P.,
    5. Marth G. T.
    , 2011 BamTools: a C++ API and toolkit for analyzing and managing BAM files. Bioinformatics 27: 1691–1692.
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    1. Brand P.,
    2. Ramírez S. R.,
    3. Leese F.,
    4. Quezada-Euan J. J.,
    5. Tollrian R.,
    6. et al.
    , 2015 Rapid evolution of chemosensory receptor genes in a pair of sibling species of orchid bees (Apidae: Euglossini). BMC Evol. Biol. 15: 176.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Branstetter M. G.,
    2. Danforth B. N.,
    3. Pitts J. P.,
    4. Faircloth B. C.,
    5. Ward P. S.,
    6. et al.
    , 2017 Phylogenomic insights into the evolution of stinging wasps and the origins of ants and bees. Curr. Biol. 27: 1019–1025.
    OpenUrlCrossRef
  6. ↵
    1. Camacho C.,
    2. Coulouris G.,
    3. Avagyan V.,
    4. Ma N.,
    5. Papadopoulos J.,
    6. et al.
    , 2009 BLAST+: architecture and applications. BMC Bioinformatics 10: 421.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Cardinal S.,
    2. Danforth B. N.
    , 2011 The antiquity and evolutionary history of social behavior in bees. PLoS One 6: e21086.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cha S. Y.,
    2. Yoon H. J.,
    3. Lee E. M.,
    4. Yoon M. H.,
    5. Hwang J. S.,
    6. et al.
    , 2007 The complete nucleotide sequence and gene organization of the mitochondrial genome of the bumblebee, Bombus ignitus (Hymenoptera: Apidae). Gene 392: 206–220.
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    1. Chalopin D.,
    2. Naville M.,
    3. Plard F.,
    4. Galiana D.,
    5. Volff J.-N.
    , 2015 Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates. Genome Biol. Evol. 7: 567–580.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Chu C.,
    2. Nielsen R.,
    3. Wu Y.
    , 2016 REPdenovo: inferring de novo repeat motifs from short sequence reads. PLoS One 11: e0150719.
    OpenUrl
  11. ↵
    1. Clark A. G.,
    2. Eisen M. B.,
    3. Smith D. R.,
    4. Bergman C. M.,
    5. Oliver B.,
    6. et al.
    , 2007 Evolution of genes and genomes on the Drosophila phylogeny. Nature 450: 203–218.
    OpenUrlCrossRefPubMedWeb of Science
  12. ↵
    1. Crozier R. H.,
    2. Crozier Y. C.
    , 1993 The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics 133: 97–117.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. d’Alençon E.,
    2. Sezutsu H.,
    3. Legeai F.,
    4. Permal E.,
    5. Bernard-Samain S.,
    6. et al.
    , 2010 Extensive synteny conservation of holocentric chromosomes in Lepidoptera despite high rates of local genome rearrangements. Proc. Natl. Acad. Sci. USA 107: 7680–7685.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Danforth B. N.
    , 2002 Evolution of sociality in a primitively eusocial lineage of bees. Proc. Natl. Acad. Sci. USA 99: 286–290.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Dietz L.,
    2. Brand P.,
    3. Eschner L. M.,
    4. Leese F.
    , 2016 The mitochondrial genomes of the caddisflies Sericostoma personatum and Thremma gallicum (Insecta: Trichoptera). Mitochondrial DNA A DNA Mapp. Seq. Anal. 27: 3293–3294.
    OpenUrl
  16. ↵
    1. Doetterl S.,
    2. Vereecken N. J.
    , 2010 The chemical ecology and evolution of bee–flower interactions: a review and perspectives. Can. J. Zool. 88: 668–697.
    OpenUrlCrossRefWeb of Science
  17. ↵
    1. Elsik C. G.,
    2. Worley K. C.,
    3. Bennett A. K.,
    4. Beye M.,
    5. Camara F.,
    6. et al.
    , 2014 Finding the missing honey bee genes: lessons learned from a genome upgrade. BMC Genomics 15: 86.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Elsik C. G.,
    2. Tayal A.,
    3. Diesh C. M.,
    4. Unni D. R.,
    5. Emery M. L.,
    6. et al.
    , 2016 Hymenoptera genome database: integrating genome annotations in HymenopteraMine. Nucleic Acids Res. 44: D793–D800.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Eltz T.,
    2. Whitten W. M.,
    3. Roubik D. W.,
    4. Linsenmair K. E.
    , 1999 Fragrance collection, storage, and accumulation by individual male orchid bees. J. Chem. Ecol. 25: 157–176.
    OpenUrlCrossRefWeb of Science
  20. ↵
    1. Eltz T.,
    2. Zimmermann Y.,
    3. Pfeiffer C.,
    4. Pech J. R.,
    5. Twele R.,
    6. et al.
    , 2008 An olfactory shift is associated with male perfume differentiation and species divergence in orchid bees. Curr. Biol. 18: 1844–1848.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Eltz T.,
    2. Fritzsch F.,
    3. Pech J. R.,
    4. Zimmermann Y.,
    5. Ramirez S. R.,
    6. et al.
    , 2011 Characterization of the orchid bee Euglossa viridissima (Apidae: Euglossini) and a novel cryptic sibling species, by morphological, chemical, and genetic characters. Zool. J. Linnean Soc. 163: 1064–1076.
    OpenUrl
  22. ↵
    1. Garófalo C. A.
    , 1985 Social structure of Euglossa cordata nests (Hymenoptera: Apidae: Euglossini). Entomol. Gen. 11: 77–83.
    OpenUrl
  23. ↵
    1. Gnerre S.,
    2. MacCallum I.,
    3. Przybylski D.,
    4. Ribeiro F. J.,
    5. Burton J. N.,
    6. et al.
    , 2011 High-quality draft assemblies of mammalian genomes from massively parallel sequence data. Proc. Natl. Acad. Sci. USA 108: 1513–1518.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. i5K Consortium
    , 2013 The i5K initiative: advancing arthropod genomics for knowledge, human health, agriculture, and the environment. J. Hered. 104: 595–600.
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Kajitani R.,
    2. Toshimoto K.,
    3. Noguchi H.,
    4. Toyoda A.,
    5. Ogura Y.,
    6. et al.
    , 2014 Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 24: 1384–1395.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Kapheim K. M.,
    2. Pan H.,
    3. Li C.,
    4. Salzberg S. L.,
    5. Puiu D.,
    6. et al.
    , 2015 Social evolution. Genomic signatures of evolutionary transitions from solitary to group living. Science 348: 1139–1143.
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Kelley D. R.,
    2. Salzberg S. L.
    , 2010 Detection and correction of false segmental duplications caused by genome mis-assembly. Genome Biol. 11: R28.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Kidwell M. G.
    , 2002 Transposable elements and the evolution of genome size in eukaryotes. Genetica 115: 49–63.
    OpenUrlCrossRefPubMedWeb of Science
  29. ↵
    1. Klein A.-M.,
    2. Vaissière B. E.,
    3. Cane J. H.,
    4. Steffan-Dewenter I.,
    5. Cunningham S. A.,
    6. et al.
    , 2007 Importance of pollinators in changing landscapes for world crops. Proc. Biol. Sci. 274: 303–313.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Kocher S. D.,
    2. Li C.,
    3. Yang W.,
    4. Tan H.,
    5. Yi S. V.,
    6. et al.
    , 2013 The draft genome of a socially polymorphic halictid bee, Lasioglossum albipes. Genome Biol. 14: R142.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Korf I.
    , 2004 Gene finding in novel genomes. BMC Bioinformatics 5: 59.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Langmead B.,
    2. Salzberg S. L.
    , 2012 Fast gapped-read alignment with Bowtie 2. Nat. Methods 9: 357–359.
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    1. Laslett D.,
    2. Canbäck B.
    , 2008 ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics 24: 172–175.
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    1. Leese F.,
    2. Brand P.,
    3. Rozenberg A.,
    4. Mayer C.,
    5. Agrawal S.,
    6. et al.
    , 2012 Exploring Pandora’s box: potential and pitfalls of low coverage genome surveys for evolutionary biology. PLoS One 7: e49202.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Litman J. R.,
    2. Danforth B. N.,
    3. Eardley C. D.,
    4. Praz C. J.
    , 2011 Why do leafcutter bees cut leaves? New insights into the early evolution of bees. Proc. Biol. Sci. 278: 3593–3600.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Lowe T. M.,
    2. Eddy S. R.
    , 1997 tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25: 955–964.
    OpenUrlCrossRefPubMed
  37. ↵
    Luo, R., B. Liu, Y. Xie, Z. Li, W. Huang et al., 2012 SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience 1: 18 (erratum: GigaScience 4: 30).
  38. ↵
    Mayer, C., 2010 Phobos Version 3.3.12. A Tandem Repeat Search Program. 20 p. Available at: http://www.rub.de/spezzoo/cm/cm_phobos.htm. Accessed: September 30, 2016.
  39. ↵
    1. Menzel R.,
    2. Muller U.
    , 1996 Learning and memory in honeybees: from behavior to neural substrates. Annu. Rev. Neurosci. 19: 379–404.
    OpenUrlCrossRefPubMedWeb of Science
  40. ↵
    1. Michener C. D.
    , 2007 The Bees of the World, Ed. 2. The John Hopkins University Press, Baltimore, MD.
  41. ↵
    1. Neafsey D. E.,
    2. Waterhouse R. M.,
    3. Abai M. R.,
    4. Aganezov S. S.,
    5. Alekseyev M. A.,
    6. et al.
    , 2015 Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes. Science 347: 1258522.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Nowak M. A.,
    2. Tarnita C. E.,
    3. Wilson E. O.
    , 2010 The evolution of eusociality. Nature 466: 1057–1062.
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    1. Nygaard S.,
    2. Hu H.,
    3. Li C.,
    4. Schiøtt M.,
    5. Chen Z.,
    6. et al.
    , 2016 Reciprocal genomic evolution in the ant–fungus agricultural symbiosis. Nat. Commun. 7: 12233.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Obbard D. J.,
    2. Maclennan J.,
    3. Kim K.-W.,
    4. Rambaut A.,
    5. O’Grady P. M.,
    6. et al.
    , 2012 Estimating divergence dates and substitution rates in the Drosophila phylogeny. Mol. Biol. Evol. 29: 3459–3473.
    OpenUrlCrossRefPubMedWeb of Science
  45. ↵
    1. Pamilo P.,
    2. Viljakainen L.,
    3. Vihavainen A.
    , 2007 Exceptionally high density of NUMTs in the honeybee genome. Mol. Biol. Evol. 24: 1340–1346.
    OpenUrlCrossRefPubMedWeb of Science
  46. ↵
    1. Park D.,
    2. Jung J. W.,
    3. Choi B.-S.,
    4. Jayakodi M.,
    5. Lee J.,
    6. et al.
    , 2015 Uncovering the novel characteristics of Asian honey bee, Apis cerana, by whole genome sequencing. BMC Genomics 16: 1.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Parra G.,
    2. Bradnam K.,
    3. Korf I.
    , 2007 CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23: 1061–1067.
    OpenUrlCrossRefPubMedWeb of Science
  48. ↵
    Pech, M. E. C., W. de J. May-Itzá, L. A. M. Medina, and J. J. G. Quezada-Euan, 2008 Sociality in Euglossa (Euglossa) viridissima Friese (Hymenoptera, Apidae, Euglossini). Insectes Sociaux 55: 428–433.
  49. ↵
    1. Pemberton R. W.,
    2. Wheeler G. S.
    , 2006 Orchid bees don’t need orchids: evidence from the naturalization of an orchid bee in Florida. Ecology 87: 1995–2001.
    OpenUrlCrossRefPubMedWeb of Science
  50. ↵
    1. Peters R. S.,
    2. Krogmann L.,
    3. Mayer C.,
    4. Donath A.,
    5. Gunkel S.,
    6. et al.
    , 2017 Evolutionary history of the Hymenoptera. Curr. Biol. 27: 1013–1018.
    OpenUrlCrossRef
  51. ↵
    1. Platt R. N.,
    2. Blanco-Berdugo L.,
    3. Ray D. A.
    , 2016 Accurate transposable element annotation is vital when analyzing new genome assemblies. Genome Biol. Evol. 8: 403–410.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Ramírez S.,
    2. Dressler R. L.,
    3. Ospina M.
    , 2002 Abejas euglosinas (Hymenoptera: Apidae) de la región Neotropical: Listado de especies con notas sobre su biología. Biota Colombiana 3: 7–118.
    OpenUrl
  53. ↵
    1. Ramírez S. R.,
    2. Eltz T.,
    3. Fujiwara M. K.,
    4. Gerlach G.,
    5. Goldman-Huertas B.,
    6. et al.
    , 2011 Asynchronous diversification in a specialized plant-pollinator mutualism. Science 333: 1742–1746.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Romiguier J.,
    2. Cameron S. A.,
    3. Woodard S. H.,
    4. Fischman B. J.,
    5. Keller L.,
    6. et al.
    , 2016 Phylogenomics controlling for base compositional bias reveals a single origin of eusociality in corbiculate bees. Mol. Biol. Evol. 33: 670–678.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Roubik D. W.,
    2. Hanson P. E.
    , 2004 Orchid Bees of Tropical America: Biology and Field Guide. Instituto Nacional de Biodiversidad (INBio), Santo Domingo De Heredia, Costa Rica.
  56. ↵
    1. Sadd B. M.,
    2. Barribeau S. M.,
    3. Bloch G.,
    4. de Graaf D. C.,
    5. Dearden P.,
    6. et al.
    , 2015 The genomes of two key bumblebee species with primitive eusocial organization. Genome Biol. 16: 76.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Simão F. A.,
    2. Waterhouse R. M.,
    3. Ioannidis P.,
    4. Kriventseva E. V.,
    5. Zdobnov E. M.
    , 2015 BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31: 3210–3212.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Simola D. F.,
    2. Wissler L.,
    3. Donahue G.,
    4. Waterhouse R. M.,
    5. Helmkampf M.,
    6. et al.
    , 2013 Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality. Genome Res. 23: 1235–1247.
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Skov C.,
    2. Wiley J.
    , 2005 Establishment of the neotropical orchid bee Euglossa viridissima (Hymenoptera: Apidae) in Florida. Fla. Entomol. 88: 225–227.
    OpenUrlCrossRef
  60. ↵
    1. Slater G. S. C.,
    2. Birney E.
    , 2005 Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics 6: 31.
    OpenUrlCrossRefPubMed
    1. Smit A. F. A,
    2. Hubley R.,
    3. Green P.
    , 2013-2017 RepeatMasker Open-4.0.http://www.repeatmasker.org.
  61. ↵
    1. Soro A.,
    2. Quezada-Euan J. J. G.,
    3. Theodorou P.,
    4. Moritz R. F. A.,
    5. Paxton R. J.
    , 2016 The population genetics of two orchid bees suggests high dispersal, low diploid male production and only an effect of island isolation in lowering genetic diversity. Conserv. Genet. 18: 607–619.
    OpenUrl
  62. ↵
    1. Stanke M.,
    2. Diekhans M.,
    3. Baertsch R.,
    4. Haussler D.
    , 2008 Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24: 637–644.
    OpenUrlCrossRefPubMedWeb of Science
  63. ↵
    1. Stolle E.,
    2. Wilfert L.,
    3. Schmid-Hempel R.,
    4. Schmid-Hempel P.,
    5. Kube M.,
    6. et al.
    , 2011 A second generation genetic map of the bumblebee Bombus terrestris (Linnaeus, 1758) reveals slow genome and chromosome evolution in the Apidae. BMC Genomics 12: 48.
    OpenUrlCrossRefPubMed
  64. ↵
    1. Suni S. S.
    , 2016 Dispersal of the orchid bee Euglossa imperialis over degraded habitat and intact forest. Conserv. Genet. 18: 621–630.
    OpenUrl
  65. ↵
    1. Suni S. S.,
    2. Brosi B. J.
    , 2012 Population genetics of orchid bees in a fragmented tropical landscape. Conserv. Genet. 13: 323–332.
    OpenUrl
  66. ↵
    1. Vogel S.
    , 1966 Parfümsammelnde bienen als bestäuber von orchidaceen und Gloxinia. Österr. Bot. Z. 113: 302–361.
    OpenUrlCrossRef
  67. ↵
    Wcislo, W. T., and J. H. Cane, 1996 Floral resource utilization by solitary bees (Hymenoptera: Apoidea) and exploitation of their stored foods by natural enemies. Annu. Rev. Entomol. 41: 257–286.
  68. ↵
    1. Weber M. G.,
    2. Mitko L.,
    3. Eltz T.,
    4. Ramírez S. R.
    , 2016 Macroevolution of perfume signalling in orchid bees. Ecol. Lett. 19: 1314–1323.
    OpenUrl
  69. ↵
    1. Weinstock G. M.,
    2. Robinson G. E.,
    3. Gibbs R. A.,
    4. Worley K. C.,
    5. Evans J. D.,
    6. et al.
    , 2006 Insights into social insects from the genome of the honeybee Apis mellifera. Nature 443: 931–949.
    OpenUrlCrossRefPubMedWeb of Science
  70. ↵
    1. Whitten W. M.,
    2. Young A. M.,
    3. Stern D. L.
    , 1993 Nonfloral sources of chemicals that attract male euglossine bees (Apidae: Euglossini). J. Chem. Ecol. 19: 3017–3027.
    OpenUrlCrossRefWeb of Science
  71. ↵
    1. Woodard S. H.,
    2. Fischman B. J.,
    3. Venkat A.,
    4. Hudson M. E.,
    5. Varala K.,
    6. et al.
    , 2011 Genes involved in convergent evolution of eusociality in bees. Proc. Natl. Acad. Sci. USA 108: 7472–7477.
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Xu H.,
    2. Luo X.,
    3. Qian J.,
    4. Pang X.,
    5. Song J.,
    6. et al.
    , 2012 FastUniq: a fast de novo duplicates removal tool for paired short reads. PLoS One 7: e52249.
    OpenUrlCrossRefPubMed
  73. ↵
    1. Zdobnov E. M.,
    2. Tegenfeldt F.,
    3. Kuznetsov D.,
    4. Waterhouse R. M.,
    5. Simão F. A.,
    6. et al.
    , 2017 OrthoDB v9.1: cataloging evolutionary and functional annotations for animal, fungal, plant, archaeal, bacterial and viral orthologs. Nucleic Acids Res. 45: D744–D749.
    OpenUrlCrossRefPubMed
  74. ↵
    1. Zimmermann Y.,
    2. Schorkopf D. L. P.,
    3. Moritz R. F. A.,
    4. Pemberton R. W.,
    5. Quezada-Euan J. J. G.,
    6. et al.
    , 2011 Population genetic structure of orchid bees (Euglossini) in anthropogenically altered landscapes. Conserv. Genet. 12: 1183–1194.
    OpenUrl
View Abstract
Previous ArticleNext Article
Back to top

PUBLICATION INFORMATION

Volume 7 Issue 9, September 2017

G3: Genes|Genomes|Genetics: 7 (9)

ARTICLE CLASSIFICATION

Genome Report
View this article with LENS
Email

Thank you for sharing this G3: Genes | Genomes | Genetics article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma
(Your Name) has forwarded a page to you from G3: Genes | Genomes | Genetics
(Your Name) thought you would be interested in this article in G3: Genes | Genomes | Genetics.
Print
Alerts
Enter your email below to set up alert notifications for new article, or to manage your existing alerts.
SIGN UP OR SIGN IN WITH YOUR EMAIL
View PDF
Share

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma

View ORCID ProfilePhilipp Brand, Nicholas Saleh, Hailin Pan, Cai Li, Karen M. Kapheim and Santiago R. Ramírez
G3: Genes, Genomes, Genetics September 1, 2017 vol. 7 no. 9 2891-2898; https://doi.org/10.1534/g3.117.043687
Philipp Brand
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Philipp Brand
  • For correspondence: pbrand@ucdavis.edu
Nicholas Saleh
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hailin Pan
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cai Li
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karen M. Kapheim
Department of Biology, Utah State University, Logan, Utah 84322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Santiago R. Ramírez
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Citation

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma

View ORCID ProfilePhilipp Brand, Nicholas Saleh, Hailin Pan, Cai Li, Karen M. Kapheim and Santiago R. Ramírez
G3: Genes, Genomes, Genetics September 1, 2017 vol. 7 no. 9 2891-2898; https://doi.org/10.1534/g3.117.043687
Philipp Brand
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Philipp Brand
  • For correspondence: pbrand@ucdavis.edu
Nicholas Saleh
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616Graduate Group in Population Biology, University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hailin Pan
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cai Li
China National Genebank, Beijing Genomics Institute, Shenzhen 518083, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karen M. Kapheim
Department of Biology, Utah State University, Logan, Utah 84322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Santiago R. Ramírez
Department for Evolution and Ecology, Center for Population Biology, and University of California, Davis, California 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Related Articles

Cited By

More in this TOC Section

  • Long-Read Genome Sequence of the Sugar Beet Rhizosphere Mycoparasite Pythium oligandrum
  • Genome Assembly and Analysis of the North American Mountain Goat (Oreamnos americanus) Reveals Species-Level Responses to Extreme Environments
  • De Novo Genome Assembly of Populus simonii Further Supports That Populus simonii and Populus trichocarpa Belong to Different Sections
Show more Genome Report
  • Top
  • Article
    • Abstract
    • Materials and Methods
    • Results and Discussion
    • Acknowledgments
    • Footnotes
    • Literature Cited
  • Figures & Data
  • Supplemental
  • Info & Metrics

GSA

The Genetics Society of America (GSA), founded in 1931, is the professional membership organization for scientific researchers and educators in the field of genetics. Our members work to advance knowledge in the basic mechanisms of inheritance, from the molecular to the population level.

Online ISSN: 2160-1836

  • For Authors
  • For Reviewers
  • For Advertisers
  • Submit a Manuscript
  • Editorial Board
  • Press Releases

SPPA Logo

GET CONNECTED

RSS  Subscribe with RSS.

email  Subscribe via email. Sign up to receive alert notifications of new articles.

  • Facebook
  • Twitter
  • YouTube
  • LinkedIn
  • Google Plus

Copyright © 2019 by the Genetics Society of America

  • About G3
  • Terms of use
  • Permissions
  • Contact us
  • International access