Evolution of a complex locus: exon gain, loss and divergence at the Gr39a locus in Drosophila

Taste evolution in an herbivorous drosophilid

Plant secondary metabolites pose a challenge for generalist herbivorous insects because they are not only potentially toxic, they also may trigger aversion. On the contrary, some highly specialized herbivorous insects evolved to use these same compounds as 'token stimuli' for unambiguous determination of their host plants. Two questions that emerge from these observations are how recently derived herbivores evolve to overcome this aversion to plant secondary metabolites and the extent to which they evolve increased attraction to these same compounds. In this study, we addressed these questions by focusing on the evolution of bitter taste preferences in the herbivorous drosophilid Scaptomyza flava, which is phylogenetically nested deep in the paraphyletic Drosophila. We measured behavioral and neural responses of S. flava and a set of non-herbivorous species representing a phylogenetic gradient (S. pallida, S. hsui, and D. melanogaster) towards host- and non-host derived bitter plant compounds. We observed that S. flava evolved a shift in bitter detection, rather than a narrow shift towards glucosinolates, the precursors of mustard-specific defense compounds. In a dye-based consumption assay, S. flava exhibited shifts in aversion toward the non-mustard bitter, plant-produced alkaloids caffeine and lobeline, and reduced aversion towards glucosinolates, whereas the non-herbivorous species each showed strong aversion to all bitter compounds tested. We then examined whether these changes in bitter preferences of S. flava could be explained by changes in sensitivity in the peripheral nervous system and compared electrophysiological responses from the labellar sensilla of S. flava, S. pallida, and D. melanogaster. Using scanning electron microscopy, we also created a map of labellar sensilla in S. flava and S. pallida. We assigned each sensillum to a functional sensilla class based on their morphology and initial response profiles to bitter and sweet compounds. Despite a high degree of conservation in the morphology and spatial placement of sensilla between S. flava and S. pallida, electrophysiological studies revealed that S. flava had reduced sensitivity to glucosinolates to varying degrees. We found this reduction only in I type sensilla. Finally, we speculate on the potential role that evolutionary genetic changes in gustatory receptors between S. pallida and S. flava may play in driving these patterns. Specifically, we hypothesize that the evolution of bitter receptors expressed in I type sensilla may have driven the reduced sensitivity observed in S. flava, and ultimately, its reduced bitter aversion. The S. flava system showcases the importance of reduced aversion to bitter defense compounds in relatively young herbivorous lineages, and how this may be achieved at the molecular and physiological level.

Genome-wide comparative identification and analysis of membrane-FADS-like superfamily genes in freshwater economic fishes

Membrane fatty acid desaturase (FADS)-like superfamily proteins (FADSs) are essential for the synthesis of unsaturated fatty acids (UFAs). Recently, studies on FADS in fishes have mostly focused on marine species, and a comprehensive analysis of the FADS superfamily, including the FADS, stearoyl-CoA desaturase (SCD), and sphingolipid delta 4-desaturase (DEGS) families, in freshwater economic fishes is urgently required. To this end, we conducted a thorough analysis of the number, gene/protein structure, chromosomal location, gene linkage map, phylogeny, and expression of the FADS superfamily. We identified 156 FADSs genes in the genome of 27 representative species. Notably, FADS1 and SCD5 were lost in most freshwater fish and other teleosts. All FADSs proteins contain 4 transmembrane helices and 2-3 amphipathic α-helices. FADSs in the same family are often linked on the same chromosome; moreover, FADS and SCD or DEGS are frequently collocated on the same chromosome. In addition, FADS, SCD, and DEGS family proteins share similar evolutionary patterns. Interestingly, FADS6, as a member of the FADS family, exhibits a similar gene structure and chromosome location to that of SCD family members, which may be the transitional form of FADS and SCD. This study shed light on the type, structure, and phylogenetic relationship of FADSs in freshwater fishes, offering a new perspective into the functional mechanism analysis of FADSs.

Utilizing paralogs for phylogenetic reconstruction has the potential to increase species tree support and reduce gene tree discordance in target enrichment data

The analysis of target enrichment data in phylogenetics lacks optimization toward using paralogs for phylogenetic reconstruction. We developed a novel approach of detecting paralogs and utilizing them for phylogenetic tree inference, by retrieving both ortho- and paralogous copies and creating orthologous alignments, from which the gene trees are built. We implemented this approach in ParalogWizard and demonstrate its performance in plant groups that underwent a whole genome duplication relatively recently: the subtribe Malinae (family Rosaceae), using Angiosperms353 as well as Malinae481 probes, the genus Oritrophium (family Asteraceae), using Compositae1061 probes, and the genus Amomum (family Zingiberaceae), using Zingiberaceae1180 probes. Discriminating between orthologs and paralogs reduced gene tree discordance and increased the species tree support in the case of the Malinae, but not for Oritrophium and Amomum. This may relate to the difference in the proportion of paralogous loci between the datasets, which was highest for the Malinae. Overall, retrieving paralogs for phylogenetic reconstruction following ParalogWizard has the potential to increase the species tree support and reduce gene tree discordance in target enrichment data, particularly if the proportion of paralogous loci is high.

Characterization of gustatory receptor 7 in the brown planthopper reveals functional versatility

Insect pests consume tastants as their necessary energy and nutrient sources. Gustatory receptors play important roles in insect life and can form within an extremely complicated regulatory network. However, there are still many gustatory genes that have a significant impact on insect physiology, but their functional mechanism is still unknown. Here, we purified and characterized a gustatory receptor (protein) coding gene, NlGr7, from the brown planthopper (BPH) Nilaparvata lugens, which is an important insect pest of rice. Our results revealed that NlGr7 has an active association with various ligands, such as lectins, lipids (phospho- and sphingolipid) and copper. The mass-spectrometry result showed that NlGr7 is a sugar receptor, and NlGr7 is validated by different types of insoluble polysaccharides and a varied range of tastants. Further, we observed that NlGr7-bound ATP hydrolysed on the ATPase activity assay, which indicated that NlGr7 may be associated with important biological functions in the BPH. Furthermore, an injection of NlGr7 (protein), into newly emerged female adults of BPH, showed the reduced vitellogenin in ovary. The important NlGr7 for chemoreception has now been characterized in the BPH. We showed that NlGr7 in the BPH is required for various protein-ligands, as well as protein-sugars interactions, and for regulation of fecundity marker to play crucial roles in this pest. This study will provide valuable information for further functional studies of chemoreception mechanisms in this important agricultural pest.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Background

The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.

Results

This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways.

Conclusions

The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00975-9.

Molecular Logic and Evolution of Bitter Taste in Drosophila

Taste systems detect a vast diversity of toxins, which are perceived as bitter. When a species adapts to a new environment, its taste system must adapt to detect new death threats. We deleted each of six commonly expressed bitter gustatory receptors (Grs) from Drosophila melanogaster. Systematic analysis revealed that requirements for these Grs differed for the same tastant in different neurons and for different tastants in the same neuron. Responses to some tastants in some neurons required four Grs, including Gr39a. Deletions also produced increased or novel responses, supporting a model of Gr-Gr inhibitory interactions. Coexpression of four Grs conferred several bitter responses to a sugar neuron. We then examined bitter coding in three other Drosophila species. We found major evolutionary shifts. One shift depended on the concerted activity of seven Grs. This work shows how the complex logic of bitter coding provides the capacity to detect innumerable hazards and the flexibility to adapt to new ones.

Expansion of a bitter taste receptor family in a polyphagous insect herbivore

The Insect taste system plays a central role in feeding behaviours and co-evolution of insect-host interactions. Gustatory receptors form the interface between the insect taste system and the environment. From genome and transcriptome sequencing we identified 197 novel gustatory receptor (GR) genes from the polyphagous pest Helicoverpa armigera. These GRs include a significantly expanded bitter receptor family (180 GRs) that could be further divided into three categories based on polypeptide lengths, gene structure and amino acid sequence. Type 1 includes 29 bitter Gr genes that possess introns. Type 2 includes 13 long intronless bitter Gr genes, while Type 3 comprises 131 short intronless bitter Gr genes. Calcium imaging analysis demonstrated that three Type 3 GRs (HarmGR35, HarmGR50 and HarmGR195) can be activated by a crude extract of cotton leaves. HarmGR195, a GR specifically and selectively expressed in adult tarsi, showed a specific response to proline, an amino acid widely present in plant tissues. We hypothesise that the expansion in the H. armigera GR family may be functionally tied to its polyphagous behavior. Understanding the molecular basis of polyphagy may provide opportunities for the development of new environmentally friendly pest control strategies.

Evaluation of Ancestral Sequence Reconstruction Methods to Infer Nonstationary Patterns of Nucleotide Substitution

Inference of gene sequences in ancestral species has been widely used to test hypotheses concerning the process of molecular sequence evolution. However, the approach may produce spurious results, mainly because using the single best reconstruction while ignoring the suboptimal ones creates systematic biases. Here we implement methods to correct for such biases and use computer simulation to evaluate their performance when the substitution process is nonstationary. The methods we evaluated include parsimony and likelihood using the single best reconstruction (SBR), averaging over reconstructions weighted by the posterior probabilities (AWP), and a new method called expected Markov counting (EMC) that produces maximum-likelihood estimates of substitution counts for any branch under a nonstationary Markov model. We simulated base composition evolution on a phylogeny for six species, with different selective pressures on G+C content among lineages, and compared the counts of nucleotide substitutions recorded during simulation with the inference by different methods. We found that large systematic biases resulted from (i) the use of parsimony or likelihood with SBR, (ii) the use of a stationary model when the substitution process is nonstationary, and (iii) the use of the Hasegawa-Kishino-Yano (HKY) model, which is too simple to adequately describe the substitution process. The nonstationary general time reversible (GTR) model, used with AWP or EMC, accurately recovered the substitution counts, even in cases of complex parameter fluctuations. We discuss model complexity and the compromise between bias and variance and suggest that the new methods may be useful for studying complex patterns of nucleotide substitution in large genomic data sets.

Dynamics of copy number variation in host races of the pea aphid

Copy number variation (CNV) makes a major contribution to overall genetic variation and is suspected to play an important role in adaptation. However, aside from a few model species, the extent of CNV in natural populations has seldom been investigated. Here, we report on CNV in the pea aphid Acyrthosiphon pisum, a powerful system for studying the genetic architecture of host-plant adaptation and speciation thanks to multiple host races forming a continuum of genetic divergence. Recent studies have highlighted the potential importance of chemosensory genes, including the gustatory and olfactory receptor gene families (Gr and Or, respectively), in the process of host race formation. We used targeted resequencing to achieve a very high depth of coverage, and thereby revealed the extent of CNV of 434 genes, including 150 chemosensory genes, in 104 individuals distributed across eight host races of the pea aphid. We found that CNV was widespread in our global sample, with a significantly higher occurrence in multigene families, especially in Ors. We also observed a decrease in the gene probability of being completely duplicated or deleted (CDD) with increase in coding sequence length. Genes with CDD variants were usually more polymorphic for copy number, especially in the P450 gene family where toxin resistance may be related to gene dosage. We found that Gr were overrepresented among genes discriminating host races, as were CDD genes and pseudogenes. Our observations shed new light on CNV dynamics and are consistent with CNV playing a role in both local adaptation and speciation.

New gene evolution: little did we know

Genes are perpetually added to and deleted from genomes during evolution. Thus, it is important to understand how new genes are formed and how they evolve to be critical components of the genetic systems that determine the biological diversity of life. Two decades of effort have shed light on the process of new gene origination and have contributed to an emerging comprehensive picture of how new genes are added to genomes, ranging from the mechanisms that generate new gene structures to the presence of new genes in different organisms to the rates and patterns of new gene origination and the roles of new genes in phenotypic evolution. We review each of these aspects of new gene evolution, summarizing the main evidence for the origination and importance of new genes in evolution. We highlight findings showing that new genes rapidly change existing genetic systems that govern various molecular, cellular, and phenotypic functions.

A sugar gustatory receptor identified from the foregut of cotton bollworm Helicoverpa armigera

Helicoverpa armigera (Hübner) is one of the most polyphagous and cosmopolitan pest species, the larvae of which feed on numerous important crops. The gustatory system is critical in guiding insect feeding behavior. Here, we identified a gustatory receptor from H. armigera, HaGR9, which shows high levels of identity to DmGR43a from Drosophila melanogaster and BmGR9 from Bombyx mori. Reverse transcriptase PCR (RT-PCR) revealed HaGR9 is highly expressed in larval foregut, with little or no expression in other chemosensory tissues. Membrane topology studies indicated that, like two previously studied B. mori GRs, BmGR8 and BmGR53, HaGR9 has an inverted topology relative to G protein-coupled receptors (GPCRs), an intracellular N-terminus and an extracellular C-terminus. Calcium imaging studies confirmed HaGR9 is a sugar receptor showing dose-dependent responses to D-galactose, D-maltose, and D-fructose. This highly-expressed foregut-specific gustatory receptor may contribute to the regulation of larval feeding behavior.

Brain sex differences and function of the fruitless gene in Drosophila

The fruitless (fru) gene in Drosophila plays a pivotal role in the formation of neural circuits underlying gender-specific behaviors. Specific labeling of fru expressing neurons has revealed a core circuit responsible for male courtship behavior.Females with a small number of masculinized neuronal clusters in their brain can initiate male-type courtship behavior. By examining the correlations between the masculinized neurons and behavioral gender type, a male-specific neuronal cluster,named P1, which coexpresses fru and double sex, was identified as a putative trigger center for male-type courtship behavior. P1 neurons extend dendrite to the lateral horn,where multimodal sensory inputs converge. Molecular studies suggest that fru determines the level of masculinization of neurons by orchestrating the transcription of a set of downstream genes, which remain to be identified.

Gr39a, a highly diversified gustatory receptor in Drosophila, has a role in sexual behavior

Sexual recognition among individuals is crucial for the reproduction of animals. In Drosophila, like in many other animals, pheromones are suggested to play an important role in conveying information about an individual, such as sex, maturity and mating status. Sex-specific cuticular hydrocarbon components are thought to be major sex pheromones in Drosophila, and are postulated to act through the gustatory system, since they are mostly non-volatile chemicals. However, very little is known about the molecular and neural bases of gustatory pheromone reception. So far, a few putative gustatory receptors, including Gr32a and Gr68a, have been implicated in courtship behavior. Here, we examine another putative gustatory receptor, Gr39a, which shares a cluster with both Gr32a and Gr68a in a molecular phylogeny of the gustatory receptor family, for its potential role in courtship behavior. The Gr39a gene produces four isoforms through alternative splicing of different 5'-most exons. A quantitative real-time PCR analysis showed that the expression levels of all four splice variants of Gr39a were reduced in a fly line in which a P element was inserted into the Gr39a locus. Homozygous and hemizygous males for the P-element insertion, as well as males in which Gr39a was knocked down by RNAi, all showed reduced courtship levels toward females. The courtship levels returned to normal when the P element was excised out. A close analysis of courtship behavior of the mutant males revealed that the average duration of a continuous courtship bout was significantly shorter in the mutants than in the wild type. The results suggest that Gr39a has a role in sustaining courtship behavior in males, possibly through the reception of a stimulating arrestant pheromone.

Shaping of Drosophila male courtship posture by a gustatory pheromone

Drosophila melanogaster males display stereotypical courtship ritual in courting potential mates. The fruitless (fru) gene plays a pivotal role in the generation of male-typical courtship behavior by instructing the formation of the central circuitry underlying this behavior during development. The fru gene expression can be monitored by a reporter, fruNP21, that labels approximately 800 neurons in the adult male brain. Among these fru-expressing neurons, a male-specific neural cluster called P1 initiated male-typical courtship behavior; when the P1 cluster was ectopically produced in the female brain, many such females displayed male-typical courtship behavior toward a target female. To further elucidate the central circuitry for courtship behavior, we examined central projections of sensory neurons that appear to perceive sex pheromones sensed in the foreleg gustatory receptors. The central terminals of these sensory neurons are located near dendrites of the fru-expressing interneurons named mAL. These results suggest that different subsets of fru-expressing neurons are involved in both the sensory integration of sexual cues and the activation of motor centers that generate motor outputs for courtship behavior in the Drosophila brain.