A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins.

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

Generation and Evolution of Neural Cell Types and Circuits: Insights from the Drosophila Visual System.

The Drosophila visual system has become a premier model for probing how neural diversity is generated during development. Recent work has provided deeper insight into the elaborate mechanisms that control the range of types and numbers of neurons produced, which neurons survive, and how they interact. These processes drive visual function and influence behavioral preferences. Other studies are beginning to provide insight into how neuronal diversity evolved in insects by adding new cell types and modifying neural circuits. Some of the most powerful comparisons have been those made to the Drosophila visual system, where a deeper understanding of molecular mechanisms allows for the generation of hypotheses about the evolution of neural anatomy and function. The evolution of new neural types contributes additional complexity to the brain and poses intriguing questions about how new neurons interact with existing circuitry. We explore how such individual changes in a variety of species might play a role over evolutionary timescales. Lessons learned from the fly visual system apply to other neural systems, including the fly central brain, where decisions are made and memories are stored.

Evolution of a central neural circuit underlies Drosophila mate preferences.

Courtship rituals serve to reinforce reproductive barriers between closely related species. Drosophila melanogaster and Drosophila simulans exhibit reproductive isolation, owing in part to the fact that D. melanogaster females produce 7,11-heptacosadiene, a pheromone that promotes courtship in D. melanogaster males but suppresses courtship in D. simulans males. Here we compare pheromone-processing pathways in D. melanogaster and D. simulans males to define how these sister species endow 7,11-heptacosadiene with the opposite behavioural valence to underlie species discrimination. We show that males of both species detect 7,11-heptacosadiene using homologous peripheral sensory neurons, but this signal is differentially propagated to P1 neurons, which control courtship behaviour. A change in the balance of excitation and inhibition onto courtship-promoting neurons transforms an excitatory pheromonal cue in D. melanogaster into an inhibitory cue in D. simulans. Our results reveal how species-specific pheromone responses can emerge from conservation of peripheral detection mechanisms and diversification of central circuitry, and demonstrate how flexible nodes in neural circuits can contribute to behavioural evolution.

FlyAtlas 2 in 2022: enhancements to the Drosophila melanogaster expression atlas.

FlyAtlas 2 (flyatlas2.org) is a database and web application for studying the expression of the genes of Drosophila melanogaster in different tissues of adults and larvae. It is based on RNA-Seq data, and incorporates both genes encoding proteins and microRNAs. We have now completed the population of the database with 13 tissues from both male and female adults, five sex-specific tissues, and eight larval tissues. Larval garland cell nephrocytes have also been included. Major enhancements have been made to the application. First, a facility has been added for a 'Profile' search for genes with a similar pattern of tissue expression as a query gene. This may help establish the function of genes for which this is currently unknown. Second, a facility has been added dedicated to the larval midgut, where the difference in gene expression in the five regions of different pH can be explored. A variety of further improvements to the interface are described.

How sexual selection can drive the evolution of costly sperm ornamentation.

Post-copulatory sexual selection (PSS), fuelled by female promiscuity, is credited with the rapid evolution of sperm quality traits across diverse taxa. Yet, our understanding of the adaptive significance of sperm ornaments and the cryptic female preferences driving their evolution is extremely limited. Here we review the evolutionary allometry of exaggerated sexual traits (for example, antlers, horns, tail feathers, mandibles and dewlaps), show that the giant sperm of some Drosophila species are possibly the most extreme ornaments in all of nature and demonstrate how their existence challenges theories explaining the intensity of sexual selection, mating-system evolution and the fundamental nature of sex differences. We also combine quantitative genetic analyses of interacting sex-specific traits in D. melanogaster with comparative analyses of the condition dependence of male and female reproductive potential across species with varying ornament size to reveal complex dynamics that may underlie sperm-length evolution. Our results suggest that producing few gigantic sperm evolved by (1) Fisherian runaway selection mediated by genetic correlations between sperm length, the female preference for long sperm and female mating frequency, and (2) longer sperm increasing the indirect benefits to females. Our results also suggest that the developmental integration of sperm quality and quantity renders post-copulatory sexual selection on ejaculates unlikely to treat male-male competition and female choice as discrete processes.

The NHL domain of BRAT is an RNA-binding domain that directly contacts the hunchback mRNA for regulation.

The Drosophila protein brain tumor (Brat) forms a complex with Pumilio (Pum) and Nanos (Nos) to repress hunchback (hb) mRNA translation at the posterior pole during early embryonic development. It is currently thought that complex formation is initiated by Pum, which directly binds the hb mRNA and subsequently recruits Nos and Brat. Here we report that, in addition to Pum, Brat also directly interacts with the hb mRNA. We identify Brat-binding sites distinct from the Pum consensus motif and show that RNA binding and translational repression by Brat do not require Pum, suggesting so far unrecognized Pum-independent Brat functions. Using various biochemical and biophysical methods, we also demonstrate that the NHL (NCL-1, HT2A, and LIN-41) domain of Brat, a domain previously believed to mediate protein-protein interactions, is a novel, sequence-specific ssRNA-binding domain. The Brat-NHL domain folds into a six-bladed ß propeller, and we identify its positively charged top surface as the RNA-binding site. Brat belongs to the functional diverse TRIM (tripartite motif)-NHL protein family. Using structural homology modeling, we predict that the NHL domains of all TRIM-NHL proteins have the potential to bind RNA, indicating that Brat is part of a conserved family of RNA-binding proteins.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses.

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

Whole genome experimental maps of DNA G-quadruplexes in multiple species.

Genomic maps of DNA G-quadruplexes (G4s) can help elucidate the roles that these secondary structures play in various organisms. Herein, we employ an improved version of a G-quadruplex sequencing method (G4-seq) to generate whole genome G4 maps for 12 species that include widely studied model organisms and also pathogens of clinical relevance. We identify G4 structures that form under physiological K+ conditions and also G4s that are stabilized by the G4-targeting small molecule pyridostatin (PDS). We discuss the various structural features of the experimentally observed G-quadruplexes (OQs), highlighting differences in their prevalence and enrichment across species. Our study describes diversity in sequence composition and genomic location for the OQs in the different species and reveals that the enrichment of OQs in gene promoters is particular to mammals such as mouse and human, among the species studied. The multi-species maps have been made publicly available as a resource to the research community. The maps can serve as blueprints for biological experiments in those model organisms, where G4 structures may play a role.

Gene regulatory network architecture in different developmental contexts influences the genetic basis of morphological evolution.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

Phylogenetic position of the Drosophila fima and dentissima lineages, and the status of the D. melanogaster species group.

The subgenus Sophophora of Drosophila, which includes D. melanogaster, is an important model for the study of molecular evolution, comparative genomics, and evolutionary developmental biology. Numerous phylogenetic studies have examined species relationships in the well-known melanogaster, obscura, willistoni, and saltans species groups, as well as the relationships among these clades. In contrast, other species groups of Sophophora have been relatively neglected and have not been subjected to molecular phylogenetic analysis. Here, we focus on the endemic African Drosophila fima and dentissima lineages. We find that both these clades fall within the broadly defined melanogaster species group, but are otherwise distantly related to each other. The new phylogeny supports pervasive divergent and convergent evolution of male-specific grasping structures (sex combs). We discuss the implications of these results for defining the boundaries of the melanogaster species group, and weigh the relative merits of "splitting" and "lumping" approaches to the taxonomy of this key model system.

Evaluating the genetic architecture of quantitative traits via selection followed by inbreeding.

The deleterious mutation model proposes that quantitative trait variation should be dominated by rare, partially recessive, deleterious mutations. Following artificial selection on a focal trait, the ratio of the difference in inbreeding effects between control and selected populations (ΔB), to the difference in trait means caused by directional selection (ΔM), can inform the extent to which deleterious mutations cause quantitative trait variation. Here, we apply the ΔB/ΔM ratio test to two quantitative traits (male mating success and body size) in Drosophila melanogaster. For both traits, ΔB/ΔM ratios suggested that intermediate-frequency alleles, rather than rare, partially recessive alleles (i.e. deleterious mutations), caused quantitative trait variation. We discuss these results in relation to viability data, exploring how differences between regimens in segregating (measured through inbreeding) and fixed (measured through population crosses) mutational load could affect the ratio test. Finally, we present simulations that test the statistical power of the ratio test, providing guidelines for future research.

Identification and characterization of stem-bulge RNAs in Drosophila melanogaster.

Non-coding Y RNAs and stem-bulge RNAs are homologous small RNAs in vertebrates and nematodes, respectively. They share a conserved function in the replication of chromosomal DNA in these two groups of organisms. However, functional homologues have not been found in insects, despite their common early evolutionary history. Here, we describe the identification and functional characterization of two sbRNAs in Drosophila melanogaster, termed Dm1 and Dm2. The genes coding for these two RNAs were identified by a computational search in the genome of D. melanogaster for conserved sequence motifs present in nematode sbRNAs. The predicted secondary structures of Dm1 and Dm2 partially resemble nematode sbRNAs and show stability in molecular dynamics simulations. Both RNAs are phylogenetically closer related to nematode sbRNAs than to vertebrate Y RNAs. Dm1, but not Dm2 sbRNA is abundantly expressed in D. melanogaster S2 cells and adult flies. Only Dm1, but not Dm2 sbRNA can functionally replace Y RNAs in a human cell-free DNA replication initiation system. Therefore, Dm1 is the first functional sbRNA described in insects, allowing future investigations into the physiological roles of sbRNAs in the genetically tractable model organism D. melanogaster.

The evolution of dorsal-ventral patterning mechanisms in insects.

The gene regulatory network (GRN) underpinning dorsal-ventral (DV) patterning of the Drosophila embryo is among the most thoroughly understood GRNs, making it an ideal system for comparative studies seeking to understand the evolution of development. With the emergence of widely applicable techniques for testing gene function, species with sequenced genomes, and multiple tractable species with diverse developmental modes, a phylogenetically broad and molecularly deep understanding of the evolution of DV axis formation in insects is feasible. Here, we review recent progress made in this field, compare our emerging molecular understanding to classical embryological experiments, and suggest future directions of inquiry.

Evolution of a tissue-specific splicing network.

Alternative splicing of precursor mRNA (pre-mRNA) is a strategy employed by most eukaryotes to increase transcript and proteomic diversity. Many metazoan splicing factors are members of multigene families, with each member having different functions. How these highly related proteins evolve unique properties has been unclear. Here we characterize the evolution and function of a new Drosophila splicing factor, termed LS2 (Large Subunit 2), that arose from a gene duplication event of dU2AF(50), the large subunit of the highly conserved heterodimeric general splicing factor U2AF (U2-associated factor). The quickly evolving LS2 gene has diverged from the splicing-promoting, ubiquitously expressed dU2AF(50) such that it binds a markedly different RNA sequence, acts as a splicing repressor, and is preferentially expressed in testes. Target transcripts of LS2 are also enriched for performing testes-related functions. We therefore propose a path for the evolution of a new splicing factor in Drosophila that regulates specific pre-mRNAs and contributes to transcript diversity in a tissue-specific manner.

Compaction of chromatin by diverse Polycomb group proteins requires localized regions of high charge.

Polycomb group (PcG) proteins are required for the epigenetic maintenance of developmental genes in a silent state. Proteins in the Polycomb-repressive complex 1 (PRC1) class of the PcG are conserved from flies to humans and inhibit transcription. One hypothesis for PRC1 mechanism is that it compacts chromatin, based in part on electron microscopy experiments demonstrating that Drosophila PRC1 compacts nucleosomal arrays. We show that this function is conserved between Drosophila and mouse PRC1 complexes and requires a region with an overrepresentation of basic amino acids. While the active region is found in the Posterior Sex Combs (PSC) subunit in Drosophila, it is unexpectedly found in a different PRC1 subunit, a Polycomb homolog called M33, in mice. We provide experimental support for the general importance of a charged region by predicting the compacting capability of PcG proteins from species other than Drosophila and mice and by testing several of these proteins using solution assays and microscopy. We infer that the ability of PcG proteins to compact chromatin in vitro can be predicted by the presence of domains of high positive charge and that PRC1 components from a variety of species conserve this highly charged region. This supports the hypothesis that compaction is a key aspect of PcG function.

An investigation of the molecular and biochemical basis underlying chlorantraniliprole-resistant Drosophila strains and their cross-resistance to other insecticides.

Chlorantraniliprole is an anthranilic diamide insecticide that binds to the insect ryanodine receptor (RyR) and induces an uncontrolled release of Ca2+ , resulting in paralysis and ultimately death of the target insects. Recently, it was reported that chlorantraniliprole-resistant diamondback moths, Plutella xylostella Linnaeus, have mutations in their RyR. In this study, we developed two different chlorantraniliprole-resistant Drosophila melanogaster strain. The resistance ratio (RR) of the low-concentration chlorantraniliprole-treated resistant (Low-Res) strain was 2.3, while that of the high-concentration chlorantraniliprole-treated resistant (High-Res) strain was 21.3. The LC 50 of the untreated control (Con) strain was 23.8~25.9 ppm, which was significantly higher than that reported for the susceptible diamondback moth (0.03~0.51 ppm). The high LC 50 of the Con may be because the helix S2 amino acid sequence of D. melanogaster RyR ( DmRyR) is identical to the I4790M mutation of the chlorantraniliprole-resistant diamondback moths, resulting in a lower binding affinity of DmRyR for chlorantraniliprole. Among the tested detoxification enzymes, the activity of esterase was significantly increased in the two Res strains, but glutathione S-transferases and acetylcholinesterase were significantly decreased in the two Res strains. The cross-resistance of the High-Res strain to other insecticides with different modes of actions (MoAs) revealed that the RRs of the neuronal acetylcholine receptor allosteric and competitive modulators were significantly increased, while those of the Na 2+ channel modulators were significantly reduced. Our studies showed that RRs against the same insecticide vary with the treatment concentration, and that RRs against other insecticides with different MoAs can be altered.

Parallel Gene Expression Differences between Low and High Latitude Populations of Drosophila melanogaster and D. simulans.

Gene expression variation within species is relatively common, however, the role of natural selection in the maintenance of this variation is poorly understood. Here we investigate low and high latitude populations of Drosophila melanogaster and its sister species, D. simulans, to determine whether the two species show similar patterns of population differentiation, consistent with a role for spatially varying selection in maintaining gene expression variation. We compared at two temperatures the whole male transcriptome of D. melanogaster and D. simulans sampled from Panama City (Panama) and Maine (USA). We observed a significant excess of genes exhibiting differential expression in both species, consistent with parallel adaptation to heterogeneous environments. Moreover, the majority of genes showing parallel expression differentiation showed the same direction of differential expression in the two species and the magnitudes of expression differences between high and low latitude populations were correlated across species, further bolstering the conclusion that parallelism for expression phenotypes results from spatially varying selection. However, the species also exhibited important differences in expression phenotypes. For example, the genomic extent of genotype × environment interaction was much more common in D. melanogaster. Highly differentiated SNPs between low and high latitudes were enriched in the 3' UTRs and CDS of the geographically differently expressed genes in both species, consistent with an important role for cis-acting variants in driving local adaptation for expression-related phenotypes.

Genetic architecture of natural variation in Drosophila melanogaster aggressive behavior.

Aggression is an evolutionarily conserved complex behavior essential for survival and the organization of social hierarchies. With the exception of genetic variants associated with bioamine signaling, which have been implicated in aggression in many species, the genetic basis of natural variation in aggression is largely unknown. Drosophila melanogaster is a favorable model system for exploring the genetic basis of natural variation in aggression. Here, we performed genome-wide association analyses using the inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and replicate advanced intercross populations derived from the most and least aggressive DGRP lines. We identified genes that have been previously implicated in aggressive behavior as well as many novel loci, including gustatory receptor 63a (Gr63a), which encodes a subunit of the receptor for CO2, and genes associated with development and function of the nervous system. Although genes from the two association analyses were largely nonoverlapping, they mapped onto a genetic interaction network inferred from an analysis of pairwise epistasis in the DGRP. We used mutations and RNAi knock-down alleles to functionally validate 79% of the candidate genes and 75% of the candidate epistatic interactions tested. Epistasis for aggressive behavior causes cryptic genetic variation in the DGRP that is revealed by changing allele frequencies in the outbred populations derived from extreme DGRP lines. This phenomenon may pertain to other fitness traits and species, with implications for evolution, applied breeding, and human genetics.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

A natural polymorphism alters odour and DEET sensitivity in an insect odorant receptor.

Blood-feeding insects such as mosquitoes are efficient vectors of human infectious diseases because they are strongly attracted by body heat, carbon dioxide and odours produced by their vertebrate hosts. Insect repellents containing DEET (N,N-diethyl-meta-toluamide) are highly effective, but the mechanism by which this chemical wards off biting insects remains controversial despite decades of investigation. DEET seems to act both at close range as a contact chemorepellent, by affecting insect gustatory receptors, and at long range, by affecting the olfactory system. Two opposing mechanisms for the observed behavioural effects of DEET in the gas phase have been proposed: that DEET interferes with the olfactory system to block host odour recognition and that DEET actively repels insects by activating olfactory neurons that elicit avoidance behaviour. Here we show that DEET functions as a modulator of the odour-gated ion channel formed by the insect odorant receptor complex. The functional insect odorant receptor complex consists of a common co-receptor, ORCO (ref. 15) (formerly called OR83B; ref. 16), and one or more variable odorant receptor subunits that confer odour selectivity. DEET acts on this complex to potentiate or inhibit odour-evoked activity or to inhibit odour-evoked suppression of spontaneous activity. This modulation depends on the specific odorant receptor and the concentration and identity of the odour ligand. We identify a single amino-acid polymorphism in the second transmembrane domain of receptor OR59B in a Drosophila melanogaster strain from Brazil that renders OR59B insensitive to inhibition by the odour ligand and modulation by DEET. Our data indicate that natural variation can modify the sensitivity of an odour-specific insect odorant receptor to odour ligands and DEET. Furthermore, they support the hypothesis that DEET acts as a molecular 'confusant' that scrambles the insect odour code, and provide a compelling explanation for the broad-spectrum efficacy of DEET against multiple insect species.