Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector.

BACKGROUND: The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes. RESULTS: A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups. CONCLUSIONS: The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Direct and indirect genetic effects of a social supergene.

Indirect genetic effects describe phenotypic variation that results from differences in the genotypic composition of social partners. Such effects represent heritable sources of environmental variation in eusocial organisms because individuals are typically reared by their siblings. In the fire ant Solenopsis invicta, a social supergene exhibits striking indirect genetic effects on worker regulation of colony queen number, such that the genotypic composition of workers at the supergene determines whether colonies contain a single or multiple queens. We assessed the direct and indirect genetic effects of this supergene on gene expression in brains and abdominal tissues from laboratory-reared workers and compared these with previously published data from field-collected prereproductive queens. We found that direct genetic effects caused larger gene expression changes and were more consistent across tissue types and castes than indirect genetic effects. Indirect genetic effects influenced the expression of many loci but were generally restricted to the abdominal tissues. Further, indirect genetic effects were only detected when the genotypic composition of social partners differed throughout the development and adult life of focal workers, and were often only significant with relatively lenient statistical cutoffs. Our study provides insight into direct and indirect genetic effects of a social supergene on gene regulatory dynamics across tissues and castes in a complex society.

Simple inheritance, complex regulation: Supergene-mediated fire ant queen polymorphism.

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion-based supergene that spans more than 13 Mb of a "social chromosome," contains over 400 protein-coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA-sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene-associated gene upregulation, (b) allele-specific expression and (c) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion-mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene-associated gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.

Epigenetics in Insects: Genome Regulation and the Generation of Phenotypic Diversity.

Epigenetic inheritance is fundamentally important to cellular differentiation and developmental plasticity. In this review, we provide an introduction to the field of molecular epigenetics in insects. Epigenetic information is passed across cell divisions through the methylation of DNA, the modification of histone proteins, and the activity of noncoding RNAs. Much of our knowledge of insect epigenetics has been gleaned from a few model species. However, more studies of epigenetic information in traditionally nonmodel taxa will help advance our understanding of the developmental and evolutionary significance of epigenetic inheritance in insects. To this end, we also provide a brief overview of techniques for profiling and perturbing individual facets of the epigenome. Doing so in diverse cellular, developmental, and taxonomic contexts will collectively help shed new light on how genome regulation results in the generation of diversity in insect form and function.

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Divergent whole-genome methylation maps of human and chimpanzee brains reveal epigenetic basis of human regulatory evolution.

DNA methylation is a pervasive epigenetic DNA modification that strongly affects chromatin regulation and gene expression. To date, it remains largely unknown how patterns of DNA methylation differ between closely related species and whether such differences contribute to species-specific phenotypes. To investigate these questions, we generated nucleotide-resolution whole-genome methylation maps of the prefrontal cortex of multiple humans and chimpanzees. Levels and patterns of DNA methylation vary across individuals within species according to the age and the sex of the individuals. We also found extensive species-level divergence in patterns of DNA methylation and that hundreds of genes exhibit significantly lower levels of promoter methylation in the human brain than in the chimpanzee brain. Furthermore, we investigated the functional consequences of methylation differences in humans and chimpanzees by integrating data on gene expression generated with next-generation sequencing methods, and we found a strong relationship between differential methylation and gene expression. Finally, we found that differentially methylated genes are strikingly enriched with loci associated with neurological disorders, psychological disorders, and cancers. Our results demonstrate that differential DNA methylation might be an important molecular mechanism driving gene-expression divergence between human and chimpanzee brains and might potentially contribute to the evolution of disease vulnerabilities. Thus, comparative studies of humans and chimpanzees stand to identify key epigenomic modifications underlying the evolution of human-specific traits.

Evolution at two levels in fire ants: the relationship between patterns of gene expression and protein sequence evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Epigenetic inheritance and genome regulation: is DNA methylation linked to ploidy in haplodiploid insects?

Organisms show great variation in ploidy level. For example, chromosome copy number varies among cells, individuals and species. One particularly widespread example of ploidy variation is found in haplodiploid taxa, wherein males are typically haploid and females are typically diploid. Despite the prevalence of haplodiploidy, the regulatory consequences of having separate haploid and diploid genomes are poorly understood. In particular, it remains unknown whether epigenetic mechanisms contribute to regulatory compensation for genome dosage. To gain greater insights into the importance of epigenetic information to ploidy compensation, we examined DNA methylation differences among diploid queen, diploid worker, haploid male and diploid male Solenopsis invicta fire ants. Surprisingly, we found that morphologically dissimilar diploid males, queens and workers were more similar to one another in terms of DNA methylation than were morphologically similar haploid and diploid males. Moreover, methylation level was positively associated with gene expression for genes that were differentially methylated in haploid and diploid castes. These data demonstrate that intragenic DNA methylation levels differ among individuals of distinct ploidy and are positively associated with levels of gene expression. Thus, these results suggest that epigenetic information may be linked to ploidy compensation in haplodiploid insects. Overall, this study suggests that epigenetic mechanisms may be important to maintaining appropriate patterns of gene regulation in biological systems that differ in genome copy number.

Relaxed selection is a precursor to the evolution of phenotypic plasticity.

Phenotypic plasticity allows organisms to produce alternative phenotypes under different conditions and represents one of the most important ways by which organisms adaptively respond to the environment. However, the relationship between phenotypic plasticity and molecular evolution remains poorly understood. We addressed this issue by investigating the evolution of genes associated with phenotypically plastic castes, sexes, and developmental stages of the fire ant Solenopsis invicta. We first determined if genes associated with phenotypic plasticity in S. invicta evolved at a rapid rate, as predicted under theoretical models. We found that genes differentially expressed between S. invicta castes, sexes, and developmental stages all exhibited elevated rates of evolution compared with ubiquitously expressed genes. We next investigated the evolutionary history of genes associated with the production of castes. Surprisingly, we found that orthologs of caste-biased genes in S. invicta and the social bee Apis mellifera evolved rapidly in lineages without castes. Thus, in contrast to some theoretical predictions, our results suggest that rapid rates of molecular evolution may not arise primarily as a consequence of phenotypic plasticity. Instead, genes evolving under relaxed purifying selection may more readily adopt new forms of biased expression during the evolution of alternate phenotypes. These results suggest that relaxed selective constraint on protein-coding genes is an important and underappreciated element in the evolutionary origin of phenotypic plasticity.

The genome of the fire ant Solenopsis invicta.

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta, are also major pests. Here, we present a draft genome of S. invicta, assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.

The evolution of invertebrate gene body methylation.

DNA methylation of transcription units (gene bodies) occurs in the genomes of many animal and plant species. Phylogenetic persistence of gene body methylation implies biological significance; yet, the functional roles of gene body methylation remain elusive. In this study, we analyzed methylation levels of orthologs from four distantly related invertebrate species, including the honeybee, silkworm, sea squirt, and sea anemone. We demonstrate that in all four species, gene bodies distinctively cluster to two groups, which correspond to high and low methylation levels. This pattern resembles that of sequence composition arising from the mutagenetic effect of DNA methylation. In spite of this effect, our results show that protein sequences of genes targeted by high levels of methylation are conserved relative to genes lacking methylation. Our investigation identified many genes that either gained or lost methylation during the course of invertebrate evolution. Most of these genes appear to have lost methylation in the insect lineages we investigated, particularly in the honeybee. We found that genes that are methylated in all four invertebrate taxa are enriched for housekeeping functions related to transcription and translation, whereas the loss of DNA methylation occurred in genes whose functions include cellular signaling and reproductive processes. Overall, our study helps to illuminate the functional significance of gene body methylation and its impacts on genome evolution in diverse invertebrate taxa.

Characterization of gene expression patterns in response to an orthotospovirus infection between two diploid peanut species and their hybrid.

Tomato spotted wilt orthotospovirus (TSWV) transmitted by thrips causes significant yield loss in peanut (Arachis hypogaea L.) production. Use of peanut cultivars with moderate field resistance has been critical for TSWV management. However, current TSWV resistance is often not adequate, and the availability of sources of tetraploid resistance to TSWV is very limited. Allotetraploids derived by crossing wild diploid species could help introgress alleles that confer TSWV resistance into cultivated peanut. Thrips-mediated TSWV screening identified two diploids and their allotetraploid possessing the AA, BB, and AABB genomes Arachis stenosperma V10309, Arachis valida GK30011, and [A. stenosperma × A. valida]4x (ValSten1), respectively. These genotypes had reduced TSWV infection and accumulation in comparison with peanut of pure cultivated pedigree. Transcriptomes from TSWV-infected and non-infected samples from A. stenosperma, A. valida, and ValSten1 were assembled, and differentially expressed genes (DEGs) following TSWV infection were assessed. There were 3,196, 8,380, and 1,312 significant DEGs in A. stenosperma, A. valida, and ValSten1, respectively. A higher proportion of genes decreased in expression following TSWV infection for A. stenosperma and ValSten1, whereas a higher proportion of genes increased in expression following infection in A. valida. The number of DEGs previously annotated as defense-related in relation to abiotic and biotic stress was highest in A. valida followed by ValSten1 and A. stenosperma. Plant phytohormone and photosynthesis genes also were differentially expressed in greater numbers in A. valida followed by ValSten1 and A. stenosperma, with over half of those exhibiting decreases in expression.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera.

The recent, unexpected discovery of a functional DNA methylation system in the genome of the social bee Apis mellifera underscores the potential importance of DNA methylation in invertebrates. The extent of genomic DNA methylation and its role in A. mellifera remain unknown, however. Here we show that genes in A. mellifera can be divided into 2 distinct classes, one with low-CpG dinucleotide content and the other with high-CpG dinucleotide content. This dichotomy is explained by the gradual depletion of CpG dinucleotides, a well-known consequence of DNA methylation. The loss of CpG dinucleotides associated with DNA methylation also may explain the unusual mutational patterns seen in A. mellifera that lead to AT-rich regions of the genome. A detailed investigation of this dichotomy implicates DNA methylation in A. mellifera development. High-CpG genes, which are predicted to be hypomethylated in germlines, are enriched with functions associated with developmental processes, whereas low-CpG genes, predicted to be hypermethylated in germlines, are enriched with functions associated with basic biological processes. Furthermore, genes more highly expressed in one caste than another are overrepresented among high-CpG genes. Our results highlight the potential significance of epigenetic modifications, such as DNA methylation, in developmental processes in social insects. In particular, the pervasiveness of DNA methylation in the genome of A. mellifera provides fertile ground for future studies of phenotypic plasticity and genomic imprinting.

Sociality is linked to rates of protein evolution in a highly social insect.

Eusocial insects exhibit unparalleled levels of cooperation and dominate terrestrial ecosystems. The success of eusocial insects stems from the presence of specialized castes that undertake distinct tasks. We investigated whether the evolutionary transition to societies with discrete castes was associated with changes in protein evolution. We predicted that proteins with caste-biased gene expression would evolve rapidly due to reduced antagonistic pleiotropy. We found that queen-biased proteins of the honeybee Apis mellifera did indeed evolve rapidly, as predicted. However, worker-biased proteins exhibited slower evolutionary rates than queen-biased or nonbiased proteins. We suggest that distinct selective pressures operating on caste-biased genes, rather than a general reduction in pleiotropy, explain the observed differences in evolutionary rates. Our study highlights, for the first time, the interaction between highly social behavior and dynamics of protein evolution.

DNA methylation in insects: on the brink of the epigenomic era.

DNA methylation plays an important role in gene regulation in animals. However, the evolution and function of DNA methylation has only recently emerged as the subject of widespread study in insects. In this review we profile the known distribution of DNA methylation systems across insect taxa and synthesize functional inferences from studies of DNA methylation in insects and vertebrates. Unlike vertebrate genomes, which tend to be globally methylated, DNA methylation is primarily targeted to genes in insects. Nevertheless, mounting evidence suggests that a specialized role exists for genic methylation in the regulation of transcription, and possibly mRNA splicing, in both insects and mammals. Investigations in several insect taxa further reveal that DNA methylation is preferentially targeted to ubiquitously expressed genes and may play a key role in the regulation of phenotypic plasticity. We suggest that insects are particularly amenable to advancing our understanding of the biological functions of DNA methylation, because insects are evolutionarily diverse, display several lineage-specific losses of DNA methylation and possess tractable patterns of DNA methylation in moderately sized genomes.

Defense-Related Gene Expression Following an Orthotospovirus Infection Is Influenced by Host Resistance in Arachis hypogaea.

Planting resistant cultivars is the most effective tactic to manage the thrips-transmitted tomato spotted wilt orthotospovirus (TSWV) in peanut plants. However, molecular mechanisms conferring resistance to TSWV in resistant cultivars are unknown. In this study, transcriptomes of TSWV-susceptible (SunOleic 97R) and field-resistant (Tifguard) peanut cultivars with and without TSWV infection were assembled and differentially expressed genes (DEGs) were compared. There were 4605 and 2579 significant DEGs in SunOleic 97R and Tifguard, respectively. Despite the lower number of DEGs in Tifguard, an increased proportion of defense-related genes were upregulated in Tifguard than in the susceptible cultivar. Examples included disease resistance (R) proteins, leucine-rich repeats, stilbene synthase, dicer, and calmodulin. Pathway analysis revealed the increased downregulation of genes associated with defense and photosynthesis in the susceptible cultivar rather than in the resistant cultivar. These results suggest that essential physiological functions were less perturbed in the resistant cultivar than in the susceptible cultivar and that the defense response following TSWV infection was more robust in the resistant cultivar than in the susceptible cultivar.

Supergene Evolution: Recombination Finds a Way.

Supergenes are multiple linked genes that regulate complex, polymorphic traits, but little is known about their evolution. A new study of an ancient supergene in several ant species suggests that rare recombination events shape supergene evolution in surprising ways.

Evolution of a supergene that regulates a trans-species social polymorphism.

Supergenes are clusters of linked genetic loci that jointly affect the expression of complex phenotypes, such as social organization. Little is known about the origin and evolution of these intriguing genomic elements. Here we analyse whole-genome sequences of males from native populations of six fire ant species and show that variation in social organization is under the control of a novel supergene haplotype (termed Sb), which evolved by sequential incorporation of three inversions spanning half of a 'social chromosome'. Two of the inversions interrupt protein-coding genes, resulting in the increased expression of one gene and modest truncation in the primary protein structure of another. All six socially polymorphic species studied harbour the same three inversions, with the single origin of the supergene in their common ancestor inferred by phylogenomic analyses to have occurred half a million years ago. The persistence of Sb along with the ancestral SB haplotype through multiple speciation events provides a striking example of a functionally important trans-species social polymorphism presumably maintained by balancing selection. We found that while recombination between the Sb and SB haplotypes is severely restricted in all species, a low level of gene flux between the haplotypes has occurred following the appearance of the inversions, potentially mitigating the evolutionary degeneration expected at genomic regions that cannot freely recombine. These results provide a detailed picture of the structural genomic innovations involved in the formation of a supergene controlling a complex social phenotype.