Tissue-Specific cis-Regulatory Divergence Implicates eloF in Inhibiting Interspecies Mating in Drosophila

Regulatory changes in the fatty acid elongase eloF underlie the evolution of sex-specific pheromone profiles in Drosophila prolongata

BACKGROUND

Pheromones play a key role in regulating sexual behavior throughout the animal kingdom. In Drosophila and other insects, many cuticular hydrocarbons (CHCs) are sexually dimorphic, and some are known to perform pheromonal functions. However, the genetic control of sex-specific CHC production is poorly understood outside of the model species D. melanogaster. A recent evolutionary change is found in D. prolongata, which, compared to its closest relatives, shows greatly increased sexual dimorphism in both CHCs and the chemosensory system responsible for their perception. A key transition involves a male-specific increase in the proportion of long-chain CHCs.

RESULTS

Perfuming D. prolongata females with the male-biased long-chain CHCs reduces copulation success, suggesting that these compounds function as sex pheromones. The evolutionary change in CHC profiles correlates with a male-specific increase in the expression of multiple genes involved in CHC biosynthesis, including fatty acid elongases, reductases and other key enzymes. In particular, elongase F, which is responsible for producing female-specific pheromones in D. melanogaster, is strongly upregulated in D. prolongata males compared both to females and to males of the sibling species. Mutations in eloF reduce the amount of long-chain CHCs, resulting in a partial feminization of pheromone profiles in D. prolongata males. Transgenic experiments show that sex-biased expression of eloF is caused in part by a putative transposable element honghaier insertion in its regulatory region.

CONCLUSIONS

These results show that cis-regulatory changes in the eloF gene, along with other changes in the CHC synthesis pathway, contribute to the evolution of sexual communication.

Olfactory projection neuron rewiring in the brain of an ecological specialist

Animal behaviors can differ greatly between closely related species. These behavioral changes are frequently linked to sensory system modifications, but central brain cell-type alterations might also be involved. Here, we develop advanced genetic tools to compare homologous central neurons in Drosophila sechellia, an ecological specialist, with the generalist Drosophila melanogaster. Through systematic morphological analysis of olfactory projection neurons (PNs), we reveal that the global anatomy of these second-order neurons is conserved. However, high-resolution, quantitative comparisons identify a striking case of convergent rewiring of PNs in two olfactory pathways critical for D. sechellia's host location. Calcium imaging and labeling of pre-synaptic sites in these evolved D. sechellia PNs indicate that species-specific connections with third-order partners are formed. This work demonstrates that peripheral sensory evolution is accompanied by selective wiring changes in the central brain to facilitate ecological specialization and paves the way to compare other cell types throughout the nervous system.

Two major-effect loci influence interspecific mating in females of the sibling species, Drosophila simulans and D. sechellia

Secondary contact between incompletely isolated species can produce a wide variety of outcomes. The vinegar flies Drosophila simulans and D. sechellia diverged on islands in the Indian Ocean and are currently separated by partial pre- and postzygotic barriers. The recent discovery of hybridization between D. simulans and D. sechellia in the wild presents an opportunity to monitor the prevalence of alleles that influence hybridization between these sibling species. We therefore sought to identify those loci in females that affect interspecific mating, and we adapted a two-choice assay to capture female mate choice and female attractiveness simultaneously. We used shotgun sequencing to genotype female progeny of reciprocal F1 backcrosses at high resolution and performed QTL analysis. We found 2 major-effect QTL in both backcrosses, one on either arm of the third chromosome that each account for 32-37% of the difference in phenotype between species. The QTL of both backcrosses overlap and may each be alternate alleles of the same locus. Genotypes at these 2 loci followed an assortative mating pattern with D. simulans males but not D. sechellia males, which mated most frequently with females that were hybrid at both loci. These data reveal how different allele combinations at 2 major loci may promote isolation and hybridization in the same species pair. Identification of these QTLs is an important step toward understanding how the genetic architecture of mate selection may shape the outcome of secondary contact.

Regulatory Changes in the Fatty Acid Elongase eloF Underlie the Evolution of Sex-specific Pheromone Profiles in Drosophila prolongata

Pheromones play a key role in regulating sexual behavior throughout the animal kingdom. In Drosophila and other insects, many cuticular hydrocarbons (CHCs) are sexually dimorphic, and some are known to perform pheromonal functions. However, the genetic control of sex-specific CHC production is not understood outside of the model species D. melanogaster. A recent evolutionary change is found in D. prolongata, which, compared to its closest relatives, shows greatly increased sexual dimorphism in both CHCs and the chemosensory system responsible for their perception. A key transition involves a male-specific increase in the proportion of long-chain CHCs. Perfuming D. prolongata females with the male-biased CHCs reduces copulation success, suggesting that these compounds function as sex pheromones. The evolutionary change in CHC profiles correlates with a male-specific increase in the expression of multiple genes involved in CHC biosynthesis, including fatty acid elongases and reductases and other key enzymes. In particular, elongase F, which is responsible for producing female-specific pheromones in D. melanogaster, is strongly upregulated in D. prolongata males compared both to females and to males of the sibling species. Induced mutations in eloF reduce the amount of long-chain CHCs, resulting in a partial feminization of pheromone profiles in D. prolongata males while having minimal effect in females. Transgenic experiments show that sex-biased expression of eloF is caused in part by a putative transposable element insertion in its regulatory region. These results reveal one of the genetic mechanisms responsible for a recent evolutionary change in sexual communication.

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

Cell-type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell-type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell-type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell-type-specific cis-regulatory changes. We find that cell-type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell-type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Allele-specific expression reveals genetic drivers of tissue regeneration in mice

In adult mammals, skin wounds typically heal by scarring rather than through regeneration. In contrast, "super-healer" Murphy Roths Large (MRL) mice have the unusual ability to regenerate ear punch wounds; however, the molecular basis for this regeneration remains elusive. Here, in hybrid crosses between MRL and non-regenerating mice, we used allele-specific gene expression to identify cis-regulatory variation associated with ear regeneration. Analyzing three major cell populations (immune, fibroblast, and endothelial), we found that genes with cis-regulatory differences specifically in fibroblasts were associated with wound-healing pathways and also co-localized with quantitative trait loci for ear wound-healing. Ectopic treatment with one of these proteins, complement factor H (CFH), accelerated wound repair and induced regeneration in typically fibrotic wounds. Through single-cell RNA sequencing (RNA-seq), we observed that CFH treatment dramatically reduced immune cell recruitment to wounds, suggesting a potential mechanism for CFH's effect. Overall, our results provide insights into the molecular drivers of regeneration with potential clinical implications.

Cell type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell type-specific cis-regulatory changes. We find that cell type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Odor-regulated oviposition behavior in an ecological specialist

Colonization of a novel ecological niche can require, or be driven by, evolution of an animal's behaviors promoting their reproductive success. We investigated the evolution and sensory basis of oviposition in Drosophila sechellia, a close relative of Drosophila melanogaster that exhibits extreme specialism for Morinda citrifolia noni fruit. D. sechellia produces fewer eggs than other drosophilids and lays these almost exclusively on noni substrates. We show that visual, textural and social cues do not explain this species-specific preference. By contrast, we find that loss of olfactory input in D. sechellia, but not D. melanogaster, essentially abolishes egg-laying, suggesting that olfaction gates gustatory-driven noni preference. Noni odors are detected by redundant olfactory pathways, but we discover a role for hexanoic acid and the cognate Ionotropic receptor 75b (Ir75b) in odor-evoked oviposition. Through receptor exchange in D. melanogaster, we provide evidence for a causal contribution of odor-tuning changes in Ir75b to the evolution of D. sechellia's oviposition behavior.

Accounting for cis-regulatory constraint prioritizes genes likely to affect species-specific traits

Measuring allele-specific expression in interspecies hybrids is a powerful way to detect cis-regulatory changes underlying adaptation. However, it remains difficult to identify genes most likely to explain species-specific traits. Here, we outline a simple strategy that leverages population-scale allele-specific RNA-seq data to identify genes that show constrained cis-regulation within species yet show divergence between species. Applying this strategy to data from human-chimpanzee hybrid cortical organoids, we identify signatures of lineage-specific selection on genes related to saccharide metabolism, neurodegeneration, and primary cilia. We also highlight cis-regulatory divergence in CUX1 and EDNRB that may shape the trajectory of human brain development.

Speciation and adaptation research meets genome editing

Understanding the genetic basis of reproductive isolation and adaptive traits in natural populations is one of the fundamental goals in evolutionary biology. Genome editing technologies based on CRISPR-Cas systems and site-specific recombinases have enabled us to modify a targeted genomic region as desired and thus to conduct functional analyses of target loci, genes and mutations even in non-conventional model organisms. Here, we review the technical properties of genome editing techniques by classifying them into the following applications: targeted gene knock-out for investigating causative gene functions, targeted gene knock-in of marker genes for visualizing expression patterns and protein functions, precise gene replacement for identifying causative alleles and mutations, and targeted chromosomal rearrangement for investigating the functional roles of chromosomal structural variations. We describe examples of their application to demonstrate functional analysis of naturally occurring genetic variations and discuss how these technologies can be applied to speciation and adaptation research. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Genetic and genomic architecture of species-specific cuticular hydrocarbon variation in parasitoid wasps

Cuticular hydrocarbons (CHCs) serve two fundamental functions in insects: protection against desiccation and chemical signalling. How the interaction of genes shapes CHC profiles, which are essential for insect survival, adaptation and reproductive success, is still poorly understood. Here we investigate the genetic and genomic basis of CHC biosynthesis and variation in parasitoid wasps of the genus Nasonia. We mapped 91 quantitative trait loci (QTL) explaining the variation of a total of 43 CHCs in F₂ hybrid males from interspecific crosses between three Nasonia species. To identify candidate genes, we localized orthologues of CHC biosynthesis-related genes in the Nasonia genomes. We discovered multiple genomic regions where the location of QTL coincides with the location of CHC biosynthesis-related candidate genes. Most conspicuously, on a region close to the centromere of chromosome 1, multiple CHC biosynthesis-related candidate genes co-localize with several QTL explaining variation in methyl-branched alkanes. The genetic underpinnings behind this compound class are not well understood so far, despite their high potential for encoding chemical information as well as their prevalence in hymenopteran CHC profiles. Our study considerably extends our knowledge on the genetic architecture governing this important compound class, establishing a model for methyl-branched alkane genetics in the Hymenoptera in general.

Phenotypic disruption of cuticular hydrocarbon production in hybrids between sympatric species of Hawaiian picture-wing Drosophila

Interspecies hybrids can express phenotypic traits far outside the range of parental species. The atypical traits of hybrids provide insight into differences in the factors that regulate the expression of these traits in the parental species. In some cases, the unusual phenotypic traits of hybrids can lead to phenotypic dysfunction with hybrids experiencing reduced survival or reproduction. Cuticular hydrocarbons (CHCs) in insects are important phenotypic traits that serve several functions, including desiccation resistance and pheromones for mating. We used gas chromatography mass spectrometry to investigate the differences in CHC production between two closely related sympatric Hawaiian picture-wing Drosophila species, Drosophila heteroneura and D. silvestris, and their F1 and backcross hybrid offspring. CHC profiles differed between males of the two species, with substantial sexual dimorphism in D. silvestris but limited sexual dimorphism in D. heteroneura. Surprisingly, F1 hybrids did not produce three CHCs, and the abundances of several other CHCs occurred outside the ranges present in the two parental species. Backcross hybrids produced all CHCs with greater variation than observed in F1 or parental species. Overall, these results suggest that the production of CHCs was disrupted in F1 and backcross hybrids, which may have important consequences for their survival or reproduction.

Drosophila sechellia: A Genetic Model for Behavioral Evolution and Neuroecology

Defining the mechanisms by which animals adapt to their ecological niche is an important problem bridging evolution, genetics, and neurobiology. We review the establishment of a powerful genetic model for comparative behavioral analysis and neuroecology, Drosophila sechellia. This island-endemic fly species is closely related to several cosmopolitan generalists, including Drosophila melanogaster, but has evolved extreme specialism, feeding and reproducing exclusively on the noni fruit of the tropical shrub Morinda citrifolia. We first describe the development and use of genetic approaches to facilitate genotype/phenotype associations in these drosophilids. Next, we survey the behavioral, physiological, and morphological adaptations of D. sechellia throughout its life cycle and outline our current understanding of the genetic and cellular basis of these traits. Finally, we discuss the principles this knowledge begins to establish in the context of host specialization, speciation, and the neurobiology of behavioral evolution and consider open questions and challenges in the field.

A shift to shorter cuticular hydrocarbons accompanies sexual isolation among Drosophila americana group populations

Because sensory signals often evolve rapidly, they could be instrumental in the emergence of reproductive isolation between species. However, pinpointing their specific contribution to isolating barriers, and the mechanisms underlying their divergence, remains challenging. Here, we demonstrate sexual isolation due to divergence in chemical signals between two populations of Drosophila americana (SC and NE) and one population of D. novamexicana, and dissect its underlying phenotypic and genetic mechanisms. Mating trials revealed strong sexual isolation between Drosophila novamexicana males and SC Drosophila americana females, as well as more moderate bi-directional isolation between D. americana populations. Mating behavior data indicate SC D. americana males have the highest courtship efficiency and, unlike males of the other populations, are accepted by females of all species. Quantification of cuticular hydrocarbon (CHC) profiles-chemosensory signals that are used for species recognition and mate finding in Drosophila-shows that the SC D. americana population differs from the other populations primarily on the basis of compound carbon chain-length. Moreover, manipulation of male CHC composition via heterospecific perfuming-specifically perfuming D. novamexicana males with SC D. americana males-abolishes their sexual isolation from these D. americana females. Of a set of candidates, a single gene-elongase CG17821-had patterns of gene expression consistent with a role in CHC differences between species. Sequence comparisons indicate D. novamexicana and our Nebraska (NE) D. americana population share a derived CG17821 truncation mutation that could also contribute to their shared "short" CHC phenotype. Together, these data suggest an evolutionary model for the origin and spread of this allele and its consequences for CHC divergence and sexual isolation in this group.

Pure species discriminate against hybrids in the Drosophila melanogaster species subgroup

Introgression, the exchange of alleles between species, is a common event in nature. This transfer of alleles between species must happen through fertile hybrids. Characterizing the traits that cause defects in hybrids illuminates how and when gene flow is expected to occur. Inviability and sterility are extreme examples of fitness reductions but are not the only type of defects in hybrids. Some traits specific to hybrids are more subtle but are important to determine their fitness. In this report, we study whether F1 hybrids between two species pairs of Drosophila are as attractive as the parental species. We find that in both species pairs, the sexual attractiveness of the F1 hybrids is reduced and that pure species discriminate strongly against them. We also find that the cuticular hydrocarbon (CHC) profile of the female hybrids is intermediate between the parental species. Perfuming experiments show that modifying the CHC profile of the female hybrids to resemble pure species improves their chances of mating. Our results show that behavioral discrimination against hybrids might be an important component of the persistence of species that can hybridize.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

Contact-Chemosensory Evolution Underlying Reproductive Isolation in Drosophila Species

The main theme of the review is how changes in pheromone biochemistry and the sensory circuits underlying pheromone detection contribute to mate choice and reproductive isolation. The review focuses primarily on gustatory and non-volatile signals in Drosophila. Premating isolation is prevalent among closely related species. In Drosophila, preference for conspecifics against other species in mate choice underlies premating isolation, and such preference relies on contact chemosensory communications between a female and male along with other biological factors. For example, although D. simulans and D. melanogaster are sibling species that yield hybrids, their premating isolation is maintained primarily by the contrasting effects of 7,11-heptacosadiene (7,11-HD), a predominant female pheromone in D. melanogaster, on males of the two species: it attracts D. melanogaster males and repels D. simulans males. The contrasting preference for 7,11-HD in males of these two species is mainly ascribed to opposite effects of 7,11-HD on neural activities in the courtship decision-making neurons in the male brain: 7,11-HD provokes both excitatory and inhibitory inputs in these neurons and differences in the balance between the two counteracting inputs result in the contrasting preference for 7,11-HD, i.e., attraction in D. melanogaster and repulsion in D. simulans. Introduction of two double bonds is a key step in 7,11-HD biosynthesis and is mediated by the desaturase desatF, which is active in D. melanogaster females but transcriptionally inactivated in D. simulans females. Thus, 7,11-HD biosynthesis diversified in females and 7,11-HD perception diversified in males, yet it remains elusive how concordance of the changes in the two sexes was attained in evolution.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.

Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic species.

The complex genetic architecture of male mate choice evolution between Drosophila species

Mate choice behaviors are among the most important reproductive isolating barriers in many animals. Little is known about the genetic basis of reproductively isolating behaviors, but examples to date provide evidence that they can have a simple genetic basis. However, it is unclear if these results indicate that individual genes with large effects are common, or are instead due to ascertainment biases. Here, we present the results of a QTL mapping study for the most important behavioral isolating barrier between Drosophila simulans and D. sechellia: male mate choice. Our QTL results initially suggested that differences in male mate choice may be due to a couple loci with large effects. However, as we divided the largest-effect QTL using stable introgression strains, we found evidence of multiple interacting loci. We further find that separate regions of the genome control different aspects of male choice. Taken together, our results suggest that the genetic architecture of mate choice behavior, in this case, is more complex than QTL mapping suggested, highlighting potential challenges to future mapping studies. We discuss the implications of these results as they relate to signal-receiver coevolution, mate choice, and reproductive isolation.

Molecular evolution of gland cell types and chemical interactions in animals

Across the Metazoa, the emergence of new ecological interactions has been enabled by the repeated evolution of exocrine glands. Specialized glands have arisen recurrently and with great frequency, even in single genera or species, transforming how animals interact with their environment through trophic resource exploitation, pheromonal communication, chemical defense and parental care. The widespread convergent evolution of animal glands implies that exocrine secretory cells are a hotspot of metazoan cell type innovation. Each evolutionary origin of a novel gland involves a process of 'gland cell type assembly': the stitching together of unique biosynthesis pathways; coordinated changes in secretory systems to enable efficient chemical release; and transcriptional deployment of these machineries into cells constituting the gland. This molecular evolutionary process influences what types of compound a given species is capable of secreting, and, consequently, the kinds of ecological interactions that species can display. Here, we discuss what is known about the evolutionary assembly of gland cell types and propose a framework for how it may happen. We posit the existence of 'terminal selector' transcription factors that program gland function via regulatory recruitment of biosynthetic enzymes and secretory proteins. We suggest ancestral enzymes are initially co-opted into the novel gland, fostering pleiotropic conflict that drives enzyme duplication. This process has yielded the observed pattern of modular, gland-specific biosynthesis pathways optimized for manufacturing specific secretions. We anticipate that single-cell technologies and gene editing methods applicable in diverse species will transform the study of animal chemical interactions, revealing how gland cell types are assembled and functionally configured at a molecular level.

Evolution of sexually dimorphic pheromone profiles coincides with increased number of male-specific chemosensory organs in Drosophila prolongata

Binary communication systems that involve sex-specific signaling and sex-specific signal perception play a key role in sexual selection and in the evolution of sexually dimorphic traits. The driving forces and genetic changes underlying such traits can be investigated in systems where sex-specific signaling and perception have emerged recently and show evidence of potential coevolution. A promising model is found in Drosophila prolongata, which exhibits a species-specific increase in the number of male chemosensory bristles. We show that this transition coincides with recent evolutionary changes in cuticular hydrocarbon (CHC) profiles. Long-chain CHCs that are sexually monomorphic in the closest relatives of D. prolongata (D. rhopaloa, D. carrolli, D. kurseongensis, and D. fuyamai) are strongly male-biased in this species. We also identify an intraspecific female-limited polymorphism, where some females have male-like CHC profiles. Both the origin of sexually dimorphic CHC profiles and the female-limited polymorphism in D. prolongata involve changes in the relative amounts of three mono-alkene homologs, 9-tricosene, 9-pentacosene, and 9-heptacosene, all of which share a common biosynthetic origin and point to a potentially simple genetic change underlying these traits. Our results suggest that pheromone synthesis may have coevolved with chemosensory perception and open the way for reconstructing the origin of sexual dimorphism in this communication system.

Birth-and-Death Evolution of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Diversification of Cuticular Hydrocarbon Synthesis in Drosophila

The birth-and-death evolutionary model proposes that some members of a multigene family are phylogenetically stable and persist as a single copy over time, whereas other members are phylogenetically unstable and undergo frequent duplication and loss. Functional studies suggest that stable genes are likely to encode essential functions, whereas rapidly evolving genes reflect phenotypic differences in traits that diverge rapidly among species. One such class of rapidly diverging traits are insect cuticular hydrocarbons (CHCs), which play dual roles in chemical communications as short-range recognition pheromones as well as protecting the insect from desiccation. Insect CHCs diverge rapidly between related species leading to ecological adaptation and/or reproductive isolation. Because the CHC and essential fatty acid biosynthetic pathways share common genes, we hypothesized that genes involved in the synthesis of CHCs would be evolutionary unstable, whereas those involved in fatty acid-associated essential functions would be evolutionary stable. To test this hypothesis, we investigated the evolutionary history of the fatty acyl-CoA reductases (FARs) gene family that encodes enzymes in CHC synthesis. We compiled a unique data set of 200 FAR proteins across 12 Drosophila species. We uncovered a broad diversity in FAR content which is generated by gene duplications, subsequent gene losses, and alternative splicing. We also show that FARs expressed in oenocytes and presumably involved in CHC synthesis are more unstable than FARs from other tissues. Taken together, our study provides empirical evidence that a comparative approach investigating the birth-and-death evolution of gene families can identify candidate genes involved in rapidly diverging traits between species.

Transcriptional regulatory divergence underpinning species-specific learned vocalization in songbirds

Learning of most motor skills is constrained in a species-specific manner. However, the proximate mechanisms underlying species-specific learned behaviors remain poorly understood. Songbirds acquire species-specific songs through learning, which is hypothesized to depend on species-specific patterns of gene expression in functionally specialized brain regions for vocal learning and production, called song nuclei. Here, we leveraged two closely related songbird species, zebra finch, owl finch, and their interspecific first-generation (F1) hybrids, to relate transcriptional regulatory divergence between species with the production of species-specific songs. We quantified genome-wide gene expression in both species and compared this with allele-specific expression in F1 hybrids to identify genes whose expression in song nuclei is regulated by species divergence in either cis- or trans-regulation. We found that divergence in transcriptional regulation altered the expression of approximately 10% of total transcribed genes and was linked to differential gene expression between the two species. Furthermore, trans-regulatory changes were more prevalent than cis-regulatory and were associated with synaptic formation and transmission in song nucleus RA, the avian analog of the mammalian laryngeal motor cortex. We identified brain-derived neurotrophic factor (BDNF) as an upstream mediator of trans-regulated genes in RA, with a significant correlation between individual variation in BDNF expression level and species-specific song phenotypes in F1 hybrids. This was supported by the fact that the pharmacological overactivation of BDNF receptors altered the expression of its trans-regulated genes in the RA, thus disrupting the learned song structures of adult zebra finch songs at the acoustic and sequence levels. These results demonstrate functional neurogenetic associations between divergence in region-specific transcriptional regulation and species-specific learned behaviors.

Genomics Analysis of L-DOPA Exposure in Drosophila sechellia

Drosophila sechellia is a dietary specialist fruit fly that evolved from a generalist ancestor to specialize on the toxic fruit of Morinda citrifolia This species pair has been the subject of numerous studies where the goal has largely been to determine the genetic basis of adaptations associated with host specialization. Because one of the most striking features of M. citrifolia fruit is the production of toxic volatile compounds that kill insects, most genomic studies in D. sechellia to date have focused on gene expression responses to the toxic compounds in its food. In this study, we aim to identify new genes important for host specialization by profiling gene expression response to 3,4-dihydroxyphenylalanine (L-DOPA). Recent work found it to be highly abundant in M. citrifolia, critical for reproductive success of D. sechellia, and supplementation of diet with the downstream pathway product dopamine can influence toxin resistance phenotypes in related species. Here we used a combination of functional genetics and genomics techniques to identify new genes that are important for D. sechellia ecological adaptation to this new niche. We show that L-DOPA exposure can affect toxin resistance phenotypes, identify genes with plastic responses to L-DOPA exposure, and functionally test an identified candidate gene. We found that knock-down of Esterase 6 (Est6) in a heterologous species alters toxin resistance suggesting Est6 may play an important role in D. sechellia host specialization.

Speciation: On the Scent of Mate Discrimination Genes

Animals use smell to recognize individuals from their own species and find suitable mates. A study of female chemical cues in two species of fruit flies uses a creative genetic strategy to identify an allele that is involved in species recognition and may play an important role in keeping these species apart in nature.