Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect
- NSF-PAR ID:
- 10015875
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Molecular Ecology
- Volume:
- 25
- Issue:
- 8
- ISSN:
- 0962-1083
- Page Range / eLocation ID:
- p. 1769-1784
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Summary Spirodela polyrhiza is a fast‐growing aquatic monocot with highly reduced morphology, genome size and number of protein‐coding genes. Considering these biological features of Spirodela and its basal position in the monocot lineage, understanding its genome architecture could shed light on plant adaptation and genome evolution. Like many draft genomes, however, the 158‐Mb Spirodela genome sequence has not been resolved to chromosomes, and important genome characteristics have not been defined. Here we deployed rapid genome‐wide physical maps combined with high‐coverage short‐read sequencing to resolve the 20 chromosomes of Spirodela and to empirically delineate its genome features. Our data revealed a dramatic reduction in the number of therDNA repeat units in Spirodela to fewer than 100, which is even fewer than that reported for yeast. Consistent with its unique phylogenetic position, smallRNA sequencing revealed 29 Spirodela‐specific microRNA , with only two being shared withElaeis guineensis (oil palm) andMusa balbisiana (banana). CombiningDNA methylation data and smallRNA sequencing enabled the accurate prediction of 20.5% long terminal repeats (LTR s) that doubled the previous estimate, and revealed a high Solo:IntactLTR ratio of 8.2. Interestingly, we found that Spirodela has the lowest globalDNA methylation levels (9%) of any plant species tested. Taken together our results reveal a genome that has undergone reduction, likely through eliminating non‐essential protein coding genes,rDNA andLTR s. In addition to delineating the genome features of this unique plant, the methodologies described and large‐scale genome resources from this work will enable future evolutionary and functional studies of this basal monocot family. -
Abstract There is mounting evidence that, across taxa, females breeding in competitive environments tend to allocate more testosterone to their offspring prenatally and these offspring typically have more aggressive and faster‐growing phenotypes. To date, no study has determined the mechanisms mediating this maternal effect's influence on offspring phenotype. However, levels of estrogen receptor alpha (
ER α ) gene expression are linked to differences in early growth and aggression; thus, maternal hormones may alter gene regulation, perhaps viaDNA methylation, ofER α in offspring during prenatal development. We performed a pilot study to examine natural variation in testosterone allocation to offspring through egg yolks in wild Eastern Bluebirds (Sialia sialis ) in varying breeding densities and percentDNA methylation ofCG dinucleotides in theER α promoter in offspring brain regions associated with growth and behavior. We hypothesized that breeding density would be positively correlated with yolk testosterone, and prenatal exposure to maternal‐derived yolk testosterone would be associated with greater offspring growth and decreasedER α promoter methylation. Yolk testosterone concentration was positively correlated with breeding density, nestling growth rate, and percentDNA methylation of one out of five investigated CpG sites (site 3) in the diencephalonER α promoter, but none in the telencephalon (n = 10). PercentDNA methylation of diencephalon CpG site 3 was positively correlated with growth rate. These data suggest a possible role for epigenetics in mediating the effects of the maternal environment on offspring phenotype. Experimentally examining this mechanism with a larger sample size in future studies may help elucidate a prominent way in which animals respond to their environment. Further, by determining the mechanisms that mediate maternal effects, we can begin to understand the potential for the heritability of these mechanisms and the impact that maternal effects are capable of producing at an evolutionary scale. -
Summary In plants, 24 nucleotide long heterochromatic si
RNA s (het‐siRNA s) transcriptionally regulate gene expression byRNA ‐directedDNA methylation (RdDM ). The biogenesis of most het‐siRNA s depends on the plant‐specificRNA polymeraseIV (PolIV ), andARGONAUTE 4 (AGO 4) is a major het‐siRNA effector protein. Through genome‐wide analysis ofsRNA ‐seq data sets, we found that is required for the accumulation of a small subset of het‐siAGO 4RNA s. The accumulation of ‐dependent het‐siAGO 4RNA s also requires several factors known to participate in the effector portion of the RdDM pathway, includingRNA POLYMERASE V (POL V),DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM 2) andSAWADEE HOMEODOMAIN HOMOLOGUE 1 (SHH 1). Like manyAGO proteins,AGO 4 is an endonuclease that can ‘slice’RNA s. We found that a slicing‐defectiveAGO 4 was unable to fully recover dependent het‐siAGO 4‐RNA accumulation fromago4 mutant plants. Collectively, our data suggest that ‐dependent siAGO 4RNA s are secondary siRNA s dependent on the prior activity of the RdDM pathway at certain loci. -
Population genomics of divergence among extreme and intermediate color forms in a polymorphic insect
Abstract Geographic variation in insect coloration is among the most intriguing examples of rapid phenotypic evolution and provides opportunities to study mechanisms of phenotypic change and diversification in closely related lineages. The bumble bee
Bombus bifarius comprises two geographically disparate color groups characterized by red‐banded and black‐banded abdominal pigmentation, but with a range of spatially and phenotypically intermediate populations across western North America. Microsatellite analyses have revealed thatB. bifarius in the USA are structured into two major groups concordant with geography and color pattern, but also suggest ongoing gene flow among regional populations. In this study, we better resolve the relationships among major color groups to better understand evolutionary mechanisms promoting and maintaining such polymorphism. We analyze >90,000 and >25,000 single‐nucleotide polymorphisms derived from transcriptome (RNA seq) and double digest restriction site associatedDNA sequencing (ddRAD ), respectively, in representative samples from spatial and color pattern extremes inB. bifarius as well as phenotypic and geographic intermediates. Both ddRAD andRNA seq data illustrate substantial genome‐wide differentiation of the red‐banded (eastern) color form from both black‐banded (western) and intermediate (central) phenotypes and negligible differentiation among the latter populations, with no obvious admixture among bees from the two major lineages. Results thus indicate much stronger background differentiation amongB. bifarius lineages than expected, highlighting potential challenges for revealing loci underlying color polymorphism from population genetic data alone. These findings will have significance for resolving taxonomic confusion in this species and in future efforts to investigate color‐pattern evolution inB. bifarius and other polymorphic bumble bee species. -
Abstract Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (
RE ) selection in double‐digest restriction‐site‐associatedDNA sequencing (ddRAD seq) by generating reduced representation genome‐wide data using four differentRE combinations. Our expectation was thatRE selections targeting longer, more complex restriction sites would recover fewer loci thanRE with shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in thepy program. The 6‐base pair cutterRAD Eco combined with methylated site‐specific 4‐base pair cutterRI MspI produced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRAD seq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable.