A large genomic insertion containing a duplicated follistatin gene is linked to the pea aphid male wing dimorphism

Problems with Paralogs: The Promise and Challenges of Gene Duplicates in Evo-Devo Research

Gene duplicates, or paralogs, serve as a major source of new genetic material and comprise seeds for evolutionary innovation. While originally thought to be quickly lost or nonfunctionalized following duplication, now a vast number of paralogs are known to be retained in a functional state. Daughter paralogs can provide robustness through redundancy, specialize via sub-functionalization, or neo-functionalize to play new roles. Indeed, the duplication and divergence of developmental genes have played a monumental role in the evolution of animal forms (e.g., Hox genes). Still, despite their prevalence and evolutionary importance, the precise detection of gene duplicates in newly sequenced genomes remains technically challenging and often overlooked. This presents an especially pertinent problem for evolutionary developmental biology, where hypothesis testing requires accurate detection of changes in gene expression and function, often in nontraditional model species. Frequently, these analyses rely on molecular reagents designed within coding sequences that may be highly similar in recently duplicated paralogs, leading to cross-reactivity and spurious results. Thus, care is needed to avoid erroneously assigning diverged functions of paralogs to a single gene, and potentially misinterpreting evolutionary history. This perspective aims to overview the prevalence and importance of paralogs and to shed light on the difficulty of their detection and analysis while offering potential solutions.

Insect Flight: State of the Field and Future Directions

The evolution of flight in an early winged insect ancestral lineage is recognized as a key adaptation explaining the unparalleled success and diversification of insects. Subsequent transitions and modifications to flight machinery, including secondary reductions and losses, also play a central role in shaping the impacts of insects on broadscale geographic and ecological processes and patterns in the present and future. Given the importance of insect flight, there has been a centuries-long history of research and debate on the evolutionary origins and biological mechanisms of flight. Here, we revisit this history from an interdisciplinary perspective, discussing recent discoveries regarding the developmental origins, physiology, biomechanics, and neurobiology and sensory control of flight in a diverse set of insect models. We also identify major outstanding questions yet to be addressed and provide recommendations for overcoming current methodological challenges faced when studying insect flight, which will allow the field to continue to move forward in new and exciting directions. By integrating mechanistic work into ecological and evolutionary contexts, we hope that this synthesis promotes and stimulates new interdisciplinary research efforts necessary to close the many existing gaps about the causes and consequences of insect flight evolution.

Genome assembly and transcriptomic analysis to elucidate the ability of Nasonovia ribisnigri to break host plant resistance

Aphid genomic resources enable the study of complex life history traits and provide information on vector biology, host adaption and speciation. The currant-lettuce aphid (Nasonovia ribisnigri (Hemiptera: Aphididae) (Mosley)) is a cosmopolitan pest of outdoor lettuce (Lactuca sativa (Asterales: Asteraceae) (Linnaeus)). Until recently, the use of resistant cultivars was an effective method for managing N. ribisnigri. A resistant cultivar containing a single gene (Nr-locus), introduced in the 1980s, conferred complete resistance to feeding. Overreliance of this Nr-locus in lettuce resulted in N. ribisnigri's ability to break resistance mechanism, with first reports during 2003. Our work attempts to understand which candidate gene(s) are associated with this resistance-breaking mechanism. We present two de novo draft assembles for N. ribisnigri genomes, corresponding to both avirulent (Nr-locus susceptible) and virulent (Nr-locus resistant) biotypes. Changes in gene expression of the two N. ribisnigri biotypes were investigated using transcriptomic analyses of RNA-sequencing (RNA-seq) data to understand the potential mechanisms of resistance to the Nr-locus in lettuce. The draft genome assemblies were 94.2% and 91.4% complete for the avirulent and virulent biotypes, respectively. Out of the 18,872 differentially expressed genes, a single gene/locus was identified in N. ribisnigri that was shared between two resistant-breaking biotypes. This locus was further explored and validated in Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) experiments and has predicted localisations in both the cytoplasm and nucleus. This is the first study to provide evidence that a single gene/locus is likely responsible for the ability of N. ribisnigri to overcome the Nr-locus resistance in the lettuce host.

A plant virus manipulates the long-winged morph of insect vectors

Wing dimorphism of insect vectors is a determining factor for viral long-distance dispersal and large-area epidemics. Although plant viruses affect the wing plasticity of insect vectors, the potential underlying molecular mechanisms have seldom been investigated. Here, we found that a planthopper-vectored rice virus, rice stripe virus (RSV), specifically induces a long-winged morph in male insects. The analysis of field populations demonstrated that the long-winged ratios of male insects are closely associated with RSV infection regardless of viral titers. A planthopper-specific and testis-highly expressed gene, Encounter, was fortuitously found to play a key role in the RSV-induced long-winged morph. Encounter resembles malate dehydrogenase in the sequence, but it does not have corresponding enzymatic activity. Encounter is upregulated to affect male wing dimorphism at early larval stages. Encounter is closely connected with the insulin/insulin-like growth factor signaling pathway as a downstream factor of Akt, of which the transcriptional level is activated in response to RSV infection, resulting in the elevated expression of Encounter. In addition, an RSV-derived small interfering RNA directly targets Encounter to enhance its expression. Our study reveals an unreported mechanism underlying the direct regulation by a plant virus of wing dimorphism in its insect vectors, providing the potential way for interrupting viral dispersal.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

MTG-Link: leveraging barcode information from linked-reads to assemble specific loci

BACKGROUND

Local assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data.

RESULTS

We present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata.

CONCLUSIONS

MTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package.

Aphid male wing polymorphisms are transient and have evolved repeatedly

Polymorphic phenotypes have long been used to examine the maintenance of genetic variation within and between species. Most studies have focused on persistent polymorphisms, which are retained across species boundaries, and their positive effects on speciation rates. Far less is known about the macroevolutionary impacts of more transient polymorphisms, which are also common. Here we investigated male wing polymorphisms in aphids. We estimated the phylogenetic history of wing states across species, along with several other traits that could affect wing evolution. We found that male wing polymorphisms are transient: they are found in only ~4% of extant species, but have likely evolved repeatedly across the phylogeny. We reason that the repeated evolution of transient polymorphisms might be facilitated by the existence of the asexual female wing plasticity, which is common across aphids, and would maintain the wing development program even in species with wingless males. We also discovered that male wingedness correlates positively with host plant alternation and host plant breadth, and that winged morphs and wing polymorphisms may be associated with higher speciation rates. Our results provide new evolutionary insights into this well-studied group and suggest that even transient polymorphisms may impact species diversification rates.

A novel gene REPTOR2 activates the autophagic degradation of wing disc in pea aphid

Wing dimorphism in insects is an evolutionarily adaptive trait to maximize insect fitness under various environments, by which the population could be balanced between dispersing and reproduction. Most studies concern the regulatory mechanisms underlying the stimulation of wing morph in aphids, but relatively little research addresses the molecular basis of wing loss. Here, we found that, while developing normally in winged-destined pea aphids, the wing disc in wingless-destined aphids degenerated 30-hr postbirth and that this degeneration was due to autophagy rather than apoptosis. Activation of autophagy in first instar nymphs reduced the proportion of winged aphids, and suppression of autophagy increased the proportion. REPTOR2, associated with TOR signaling pathway, was identified by RNA-seq as a differentially expressed gene between the two morphs with higher expression in the thorax of wingless-destined aphids. Further genetic analysis indicated that REPTOR2 could be a novel gene derived from a gene duplication event that occurred exclusively in pea aphids on autosome A1 but translocated to the sex chromosome. Knockdown of REPTOR2 reduced autophagy in the wing disc and increased the proportion of winged aphids. In agreement with REPTOR's canonical negative regulatory role of TOR on autophagy, winged-destined aphids had higher TOR expression in the wing disc. Suppression of TOR activated autophagy of the wing disc and decreased the proportion of winged aphids, and vice versa. Co-suppression of TOR and REPTOR2 showed that dsREPTOR2 could mask the positive effect of dsTOR on autophagy, suggesting that REPTOR2 acted as a key regulator downstream of TOR in the signaling pathway. These results revealed that the TOR signaling pathway suppressed autophagic degradation of the wing disc in pea aphids by negatively regulating the expression of REPTOR2.

Investigating the genetic basis of vertebrate dispersal combining RNA-seq, RAD-seq and quantitative genetics

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

The transcription factor Zfh1 acts as a wing-morph switch in planthoppers

Insect wing polyphenism is characterized by its ability to produce two or more distinct wing morphs from a single genotype in response to changing environments. However, the molecular basis of this phenomenon remains poorly understood. Here, we identified a zinc finger homeodomain transcription factor Zfh1 that acts as an upstream regulator for the development of long-winged (LW) or shorted-winged (SW) morphs in planthoppers. Knockdown of Zfh1 directs SW-destined nymphs to develop into LW morphs by down-regulating the transcriptional level of FoxO, a prominent downstream effector of the insulin/IGF signaling (IIS) pathway. The balance between transcriptional regulation via the Zfh1-FoxO cascade and post-translational regulation via the IIS-FoxO cascade provides a flexible regulatory mechanism for the development of alternative wing morphs. These findings help us understand how phenotypic diversity is generated by altering the activity of conserved proteins, and provide an extended framework for the evolution of wing morphological diversity in insects.

Nutrition-dependent juvenile hormone sensitivity promotes flight-muscle degeneration during the aphid dispersal-reproduction transition

In insects, the loss of flight typically involves a dispersal-reproduction transition, but the underlying molecular mechanisms remain poorly understood. In the parthenogenetic pea aphid Acyrthosiphon pisum, winged females undergo flight-muscle degeneration after flight and feeding on new host plants. Similarly, topical application of a juvenile hormone (JH) mimic to starved aphids also induces flight-muscle degeneration. We found that feeding preferentially upregulated the expression of the JH receptor gene Met and a JH-inducible gene, Kr-h1, in the flight muscles, and, thus, enhanced tissue-specific JH sensitivity and signaling. RNAi-mediated knockdown of Kr-h1 prevented flight-muscle degeneration. Likewise, blocking nutritional signals by pharmacological inhibition of the target of rapamycin complex 1 (TORC1) impaired JH sensitivity of the flight muscles in feeding aphids and subsequently delayed muscle degeneration. RNA-sequencing analysis revealed that enhanced JH signaling inhibited the transcription of genes involved in the tricarboxylic acid cycle, likely resulting in reduction of the energy supply, mitochondrial dysfunction and muscle-fiber breakdown. This study shows that nutrient-dependent hormone sensitivity regulates developmental plasticity in a tissue-specific manner, emphasizing a relatively underappreciated mechanism of hormone sensitivity in modulating hormone signaling.

Pea aphid winged and wingless males exhibit reproductive, gene expression, and lipid metabolism differences

Alternative, intraspecific phenotypes offer an opportunity to identify the mechanistic basis of differences associated with distinctive life history strategies. Wing dimorphic insects, in which both flight-capable and flight-incapable individuals occur in the same population, are particularly well-studied in terms of why and how the morphs trade off flight for reproduction. Yet despite a wealth of studies examining the differences between female morphs, little is known about male differences, which could arise from different causes than those acting on females. Here we examined reproductive, gene expression, and biochemical differences between pea aphid (Acyrthosiphon pisum) winged and wingless males. We find that winged males are competitively superior in one-on-one mating circumstances, but wingless males reach reproductive maturity faster and have larger testes. We suggest that males tradeoff increased local matings with concurrent possible inbreeding for outbreeding and increased ability to find mates. At the mechanistic level, differential gene expression between the morphs revealed a possible role for activin and insulin signaling in morph differences; it also highlighted genes not previously identified as being functionally important in wing polymorphism, such as genes likely involved in sperm production. Further, we find that winged males have higher lipid levels, consistent with their use as flight fuel, but we find no consistent patterns of different levels of activity among five enzymes associated with lipid biosynthesis. Overall, our analyses provide evidence that winged versus wingless males exhibit differences at the reproductive, gene expression, and biochemical levels, expanding the field's understanding of the functional aspects of morph differences.

Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish

Background

Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens.

Results

We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication.

Conclusions

Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01205-y.

Frequent Drivers, Occasional Passengers: Signals of Symbiont-Driven Seasonal Adaptation and Hitchhiking in the Pea Aphid, Acyrthosiphon pisum

Insects harbor a variety of maternally inherited bacterial symbionts. As such, variation in symbiont presence/absence, in the combinations of harbored symbionts, and in the genotypes of harbored symbiont species provide heritable genetic variation of potential use in the insects' adaptive repertoires. Understanding the natural importance of symbionts is challenging but studying their dynamics over time can help to elucidate the potential for such symbiont-driven insect adaptation. Toward this end, we studied the seasonal dynamics of six maternally transferred bacterial symbiont species in the multivoltine pea aphid (Acyrthosiphon pisum). Our sampling focused on six alfalfa fields in southeastern Pennsylvania, and spanned 14 timepoints within the 2012 growing season, in addition to two overwintering periods. To test and generate hypotheses on the natural relevance of these non-essential symbionts, we examined whether symbiont dynamics correlated with any of ten measured environmental variables from the 2012 growing season, including some of known importance in the lab. We found that five symbionts changed prevalence across one or both overwintering periods, and that the same five species underwent such frequency shifts across the 2012 growing season. Intriguingly, the frequencies of these dynamic symbionts showed robust correlations with a subset of our measured environmental variables. Several of these trends supported the natural relevance of lab-discovered symbiont roles, including anti-pathogen defense. For a seventh symbiont-Hamiltonella defensa-studied previously across the same study periods, we tested whether a reported correlation between prevalence and temperature stemmed not from thermally varying host-level fitness effects, but from selection on co-infecting symbionts or on aphid-encoded alleles associated with this bacterium. In general, such "hitchhiking" effects were not evident during times with strongly correlated Hamiltonella and temperature shifts. However, we did identify at least one time period in which Hamiltonella spread was likely driven by selection on a co-infecting symbiont-Rickettsiella viridis. Recognizing the broader potential for such hitchhiking, we explored selection on co-infecting symbionts as a possible driver behind the dynamics of the remaining six species. Out of twelve examined instances of symbiont dynamics unfolding across 2-week periods or overwintering spans, we found eight in which the focal symbiont underwent parallel frequency shifts under single infection and one or more co-infection contexts. This supported the idea that phenotypic variation created by the presence/absence of individual symbionts is a direct target for selection, and that symbiont effects can be robust under co-habitation with other symbionts. Contrastingly, in two cases, we found that selection may target phenotypes emerging from symbiont co-infections, with specific species combinations driving overall trends for the focal dynamic symbionts, without correlated change under single infection. Finally, in three cases-including the one described above for Hamiltonella-our data suggested that incidental co-infection with a (dis)favored symbiont could lead to large frequency shifts for "passenger" symbionts, conferring no apparent cost or benefit. Such hitchhiking has rarely been studied in heritable symbiont systems. We propose that it is more common than appreciated, given the widespread nature of maternally inherited bacteria, and the frequency of multi-species symbiotic communities across insects.

Genomics of altitude-associated wing shape in two tropical butterflies

Understanding how organisms adapt to their local environment is central to evolution. With new whole-genome sequencing technologies and the explosion of data, deciphering the genomic basis of complex traits that are ecologically relevant is becoming increasingly feasible. Here, we studied the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges. Heliconius butterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1,100 butterflies from 71 broods of H. erato and H. melpomene in common-garden conditions and showed that wing aspect ratio, that is, elongatedness, is highly heritable in both species and that elevation-associated wing aspect ratio differences are maintained. Genome-wide associations with a published data set of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing aspect ratio variation in the wild. We identified several genes that have roles in wing morphogenesis or wing aspect ratio variation in Drosophila flies, making them promising candidates for future studies. There was little evidence for molecular parallelism in the two species, with only one shared candidate gene, nor for a role of the four known colour pattern loci, except for optix in H. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing aspect ratio in two Heliconius species, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

Male-Biased microRNA Discovery in the Pea Aphid

Epigenetic mechanisms modulate gene expression levels during development, shaping how a single genome produces a diversity of phenotypes. Here, we begin to explore the epigenetic regulation of sexual dimorphism in pea aphids (Acyrthosiphon pisum) by focusing on microRNAs. Previous analyses of microRNAs in aphids have focused solely on females, so we performed deep sequencing of a sample containing early-stage males. We used this sample, plus samples from Genbank, to find 207 novel pea aphid microRNA coding loci. We localized microRNA loci to a chromosome-level assembly of the pea aphid genome and found that those on the X chromosome have lower overall expression compared to those on autosomes. We then identified a set of 19 putative male-biased microRNAs and found them enriched on the X chromosome. Finally, we performed protein-coding RNA-Seq of first instar female and male pea aphids to identify genes with lower expression in males. 10 of these genes were predicted targets of the 19 male-biased microRNAs. Our study provides the most complete set of microRNAs in the pea aphid to date and serves as foundational work for future studies on the epigenetic control of sexual dimorphism.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

[Phenotypic plasticity in insects]

Insects represent 85% of the animals. They have adapted to many environments and play a major role in ecosystems. Many insect species exhibit phenotypic plasticity. We here report on the mechanisms involved in phenotypic plasticity of different insects (aphids, migratory locust, map butterfly, honeybee) and also on the nutritional size plasticity in Drosophila and the plasticity of the wing eye-spots of the butterfly Bicyclus anynana. We also describe in more detail our work concerning the thermal plasticity of pigmentation in Drosophila. We have shown that the expression of the tan, yellow and Ddc genes, encoding enzymes of the melanin synthesis pathway, is modulated by temperature and that it is a consequence, at least in part, of the temperature-sensitive expression of the bab locus genes that repress them.