Low affinity binding site clusters confer hox specificity and regulatory robustness.

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.

Dense and pleiotropic regulatory information in a developmental enhancer.

Changes in gene regulation underlie much of phenotypic evolution1. However, our understanding of the potential for regulatory evolution is biased, because most evidence comes from either natural variation or limited experimental perturbations2. Using an automated robotics pipeline, we surveyed an unbiased mutation library for a developmental enhancer in Drosophila melanogaster. We found that almost all mutations altered gene expression and that parameters of gene expression-levels, location, and state-were convolved. The widespread pleiotropic effects of most mutations may constrain the evolvability of developmental enhancers. Consistent with these observations, comparisons of diverse Drosophila larvae revealed apparent biases in the phenotypes influenced by the enhancer. Developmental enhancers may encode a higher density of regulatory information than has been appreciated previously, imposing constraints on regulatory evolution.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila.

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Evolution of the tan locus contributed to pigment loss in Drosophila santomea: a response to Matute et al.

We have shown previously that the loss of abdominal pigmentation in D. santomea relative to its sister species D. yakuba resulted, in part, from cis-regulatory mutations at the tan locus. Matute et al. claim, based solely upon extrapolation from genetic crosses of D. santomea and D. melanogaster, a much more divergent species, that at least four X chromosome regions but not tan are responsible for pigmentation differences. Here, we provide additional evidence from introgressions of D. yakuba genes into D. santomea that support a causative role for tan in the loss of pigmentation and present analyses that contradict Matute et al.'s claims. We discuss how the choice of parental species and other factors affect the ability to identify loci responsible for species divergence, and we affirm that all of our previously reported results and conclusions stand.

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.

Natural courtship song variation caused by an intronic retroelement in an ion channel gene.

Animal species display enormous variation for innate behaviours, but little is known about how this diversity arose. Here, using an unbiased genetic approach, we map a courtship song difference between wild isolates of Drosophila simulans and Drosophila mauritiana to a 966 base pair region within the slowpoke (slo) locus, which encodes a calcium-activated potassium channel. Using the reciprocal hemizygosity test, we confirm that slo is the causal locus and resolve the causal mutation to the evolutionarily recent insertion of a retroelement in a slo intron within D. simulans. Targeted deletion of this retroelement reverts the song phenotype and alters slo splicing. Like many ion channel genes, slo is expressed widely in the nervous system and influences a variety of behaviours; slo-null males sing little song with severely disrupted features. By contrast, the natural variant of slo alters a specific component of courtship song, illustrating that regulatory evolution of a highly pleiotropic ion channel gene can cause modular changes in behaviour.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola.

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Accurate and sensitive quantification of protein-DNA binding affinity.

Transcription factors (TFs) control gene expression by binding to genomic DNA in a sequence-specific manner. Mutations in TF binding sites are increasingly found to be associated with human disease, yet we currently lack robust methods to predict these sites. Here, we developed a versatile maximum likelihood framework named No Read Left Behind (NRLB) that infers a biophysical model of protein-DNA recognition across the full affinity range from a library of in vitro selected DNA binding sites. NRLB predicts human Max homodimer binding in near-perfect agreement with existing low-throughput measurements. It can capture the specificity of the p53 tetramer and distinguish multiple binding modes within a single sample. Additionally, we confirm that newly identified low-affinity enhancer binding sites are functional in vivo, and that their contribution to gene expression matches their predicted affinity. Our results establish a powerful paradigm for identifying protein binding sites and interpreting gene regulatory sequences in eukaryotic genomes.

Functional regulatory evolution outside of the minimal even-skipped stripe 2 enhancer.

Transcriptional enhancers are regions of DNA that drive precise patterns of gene expression. Although many studies have elucidated how individual enhancers can evolve, most of this work has focused on what are called 'minimal' enhancers, the smallest DNA regions that drive expression that approximates an aspect of native gene expression. Here, we explore how the Drosophila erecta even-skipped (eve) locus has evolved by testing its activity in the divergent D. melanogaster genome. We found, as has been reported previously, that the D. erecta eve stripe 2 enhancer (eveS2) fails to drive appreciable expression in D. melanogaster However, we found that a large transgene carrying the entire D. erecta eve locus drives normal eve expression, including in stripe 2. We performed a functional dissection of the region upstream of the D. erecta eveS2 region and found multiple Zelda motifs that are required for normal expression. Our results illustrate how sequences outside of minimal enhancer regions can evolve functionally through mechanisms other than changes in transcription factor-binding sites that drive patterning.

Experimental and statistical reevaluation provides no evidence for Drosophila courtship song rhythms.

From 1980 to 1992, a series of influential papers reported on the discovery, genetics, and evolution of a periodic cycling of the interval between Drosophila male courtship song pulses. The molecular mechanisms underlying this periodicity were never described. To reinitiate investigation of this phenomenon, we previously performed automated segmentation of songs but failed to detect the proposed rhythm [Arthur BJ, et al. (2013) BMC Biol 11:11; Stern DL (2014) BMC Biol 12:38]. Kyriacou et al. [Kyriacou CP, et al. (2017) Proc Natl Acad Sci USA 114:1970-1975] report that we failed to detect song rhythms because (i) our flies did not sing enough and (ii) our segmenter did not identify many of the song pulses. Kyriacou et al. manually annotated a subset of our recordings and reported that two strains displayed rhythms with genotype-specific periodicity, in agreement with their original reports. We cannot replicate this finding and show that the manually annotated data, the original automatically segmented data, and a new dataset provide no evidence for either the existence of song rhythms or song periodicity differences between genotypes. Furthermore, we have reexamined our methods and analysis and find that our automated segmentation method was not biased to prevent detection of putative song periodicity. We conclude that there is no evidence for the existence of Drosophila courtship song rhythms.

The genetic causes of convergent evolution.

The evolution of phenotypic similarities between species, known as convergence, illustrates that populations can respond predictably to ecological challenges. Convergence often results from similar genetic changes, which can emerge in two ways: the evolution of similar or identical mutations in independent lineages, which is termed parallel evolution; and the evolution in independent lineages of alleles that are shared among populations, which I call collateral genetic evolution. Evidence for parallel and collateral evolution has been found in many taxa, and an emerging hypothesis is that they result from the fact that mutations in some genetic targets minimize pleiotropic effects while simultaneously maximizing adaptation. If this proves correct, then the molecular changes underlying adaptation might be more predictable than has been appreciated previously.

Identification of loci that cause phenotypic variation in diverse species with the reciprocal hemizygosity test.

The reciprocal hemizygosity test is a straightforward genetic test that can positively identify genes that have evolved to contribute to a phenotypic difference between strains or between species. The test involves a comparison between hybrids that are genetically identical throughout the genome except at the test locus, which is rendered hemizygous for alternative alleles from the two parental strains. If the two reciprocal hemizygotes display different phenotypes, then the two parental alleles must have evolved. New methods for targeted mutagenesis will allow application of the reciprocal hemizygosity test in many organisms. This review discusses the principles, advantages, and limitations of the test.

An unsupervised method for quantifying the behavior of paired animals.

Behaviors involving the interaction of multiple individuals are complex and frequently crucial for an animal's survival. These interactions, ranging across sensory modalities, length scales, and time scales, are often subtle and difficult to characterize. Contextual effects on the frequency of behaviors become even more difficult to quantify when physical interaction between animals interferes with conventional data analysis, e.g. due to visual occlusion. We introduce a method for quantifying behavior in fruit fly interaction that combines high-throughput video acquisition and tracking of individuals with recent unsupervised methods for capturing an animal's entire behavioral repertoire. We find behavioral differences between solitary flies and those paired with an individual of the opposite sex, identifying specific behaviors that are affected by social and spatial context. Our pipeline allows for a comprehensive description of the interaction between two individuals using unsupervised machine learning methods, and will be used to answer questions about the depth of complexity and variance in fruit fly courtship.

Morphological evolution caused by many subtle-effect substitutions in regulatory DNA.

Morphology evolves often through changes in developmental genes, but the causal mutations, and their effects, remain largely unknown. The evolution of naked cuticle on larvae of Drosophila sechellia resulted from changes in five transcriptional enhancers of shavenbaby (svb), a transcript of the ovo locus that encodes a transcription factor that governs morphogenesis of microtrichiae, hereafter called 'trichomes'. Here we show that the function of one of these enhancers evolved through multiple single-nucleotide substitutions that altered both the timing and level of svb expression. The consequences of these nucleotide substitutions on larval morphology were quantified with a novel functional assay. We found that each substitution had a relatively small phenotypic effect, and that many nucleotide changes account for this large morphological difference. In addition, we observed that the substitutions had non-additive effects. These data provide unprecedented resolution of the phenotypic effects of substitutions and show how individual nucleotide changes in a transcriptional enhancer have caused morphological evolution.

TALE-mediated modulation of transcriptional enhancers in vivo.

We tested whether transcription activator-like effectors (TALEs) could mediate repression and activation of endogenous enhancers in the Drosophila genome. TALE repressors (TALERs) targeting each of the five even-skipped (eve) stripe enhancers generated repression specifically of the focal stripes. TALE activators (TALEAs) targeting the eve promoter or enhancers caused increased expression primarily in cells normally activated by the promoter or targeted enhancer, respectively. This effect supports the view that repression acts in a dominant fashion on transcriptional activators and that the activity state of an enhancer influences TALE binding or the ability of the VP16 domain to enhance transcription. In these assays, the Hairy repression domain did not exhibit previously described long-range transcriptional repression activity. The phenotypic effects of TALER and TALEA expression in larvae and adults are consistent with the observed modulations of eve expression. TALEs thus provide a novel tool for detection and functional modulation of transcriptional enhancers in their native genomic context.

Phenotypic robustness conferred by apparently redundant transcriptional enhancers.

Genes include cis-regulatory regions that contain transcriptional enhancers. Recent reports have shown that developmental genes often possess multiple discrete enhancer modules that drive transcription in similar spatio-temporal patterns: primary enhancers located near the basal promoter and secondary, or 'shadow', enhancers located at more remote positions. It has been proposed that the seemingly redundant activity of primary and secondary enhancers contributes to phenotypic robustness. We tested this hypothesis by generating a deficiency that removes two newly discovered enhancers of shavenbaby (svb, a transcript of the ovo locus), a gene encoding a transcription factor that directs development of Drosophila larval trichomes. At optimal temperatures for embryonic development, this deficiency causes minor defects in trichome patterning. In embryos that develop at both low and high extreme temperatures, however, absence of these secondary enhancers leads to extensive loss of trichomes. These temperature-dependent defects can be rescued by a transgene carrying a secondary enhancer driving transcription of the svb cDNA. Finally, removal of one copy of wingless, a gene required for normal trichome patterning, causes a similar loss of trichomes only in flies lacking the secondary enhancers. These results support the hypothesis that secondary enhancers contribute to phenotypic robustness in the face of environmental and genetic variability.

Conserved regulatory architecture underlies parallel genetic changes and convergent phenotypic evolution.

Similar morphological, physiological, and behavioral features have evolved independently in different species, a pattern known as convergence. It is known that morphological convergence can occur through changes in orthologous genes. In some cases of convergence, cis-regulatory changes generate parallel modifications in the expression patterns of orthologous genes. Our understanding of how changes in cis-regulatory regions contribute to convergence is hampered, usually, by a limited understanding of the global cis-regulatory structure of the evolving genes. Here we examine the genetic causes of a case of precise phenotypic convergence between Drosophila sechellia and Drosophila ezoana, species that diverged ~40 Mya. Previous studies revealed that changes in multiple transcriptional enhancers of shavenbaby (svb, a transcript of the ovo locus) caused phenotypic evolution in the D. sechellia lineage. It has also been shown that the convergent phenotype of D. ezoana was likely caused by cis-regulatory evolution of svb. Here we show that the large-scale cis-regulatory architecture of svb is conserved between these Drosophila species. Furthermore, we show that the D. ezoana orthologs of the evolved D. sechellia enhancers have also evolved expression patterns that correlate precisely with the changes in the phenotype. Our results suggest that phenotypic convergence resulted from multiple noncoding changes that occurred in parallel in the D. sechellia and D. ezoana lineages.

Reported Drosophila courtship song rhythms are artifacts of data analysis.

BACKGROUND: In a series of landmark papers, Kyriacou, Hall, and colleagues reported that the average inter-pulse interval of Drosophila melanogaster male courtship song varies rhythmically (KH cycles), that the period gene controls this rhythm, and that evolution of the period gene determines species differences in the rhythm's frequency. Several groups failed to recover KH cycles, but this may have resulted from differences in recording chamber size. RESULTS: Here, using recording chambers of the same dimensions as used by Kyriacou and Hall, I found no compelling evidence for KH cycles at any frequency. By replicating the data analysis procedures employed by Kyriacou and Hall, I found that two factors--data binned into 10-second intervals and short recordings--imposed non-significant periodicity in the frequency range reported for KH cycles. Randomized data showed similar patterns. CONCLUSIONS: All of the results related to KH cycles are likely to be artifacts of binning data from short songs. Reported genotypic differences in KH cycles cannot be explained by this artifact and may have resulted from the use of small sample sizes and/or from the exclusion of samples that did not exhibit song rhythms.

Multiplexed shotgun genotyping for rapid and efficient genetic mapping.

We present a new approach to genotyping based on multiplexed shotgun sequencing that can identify recombination breakpoints in a large number of individuals simultaneously at a resolution sufficient for most mapping purposes, such as quantitative trait locus (QTL) mapping and mapping of induced mutations. We first describe a simple library construction protocol that uses just 10 ng of genomic DNA per individual and makes the approach accessible to any laboratory with standard molecular biology equipment. Sequencing this library results in a large number of sequence reads widely distributed across the genomes of multiplexed bar-coded individuals. We develop a Hidden Markov Model to estimate ancestry at all genomic locations in all individuals using these data. We demonstrate the utility of the approach by mapping a dominant marker allele in D. simulans to within 105 kb of its true position using 96 F1-backcross individuals genotyped in a single lane on an Illumina Genome Analyzer. We further demonstrate the utility of our method by genetically mapping more than 400 previously unassembled D. simulans contigs to linkage groups and by evaluating the quality of targeted introgression lines. At this level of multiplexing and divergence between strains, our method allows estimation of recombination breakpoints to a median of 38-kb intervals. Our analysis suggests that higher levels of multiplexing and/or use of strains with lower levels of divergence are practicable.

Multi-channel acoustic recording and automated analysis of Drosophila courtship songs.

BACKGROUND: Drosophila melanogaster has served as a powerful model system for genetic studies of courtship songs. To accelerate research on the genetic and neural mechanisms underlying courtship song, we have developed a sensitive recording system to simultaneously capture the acoustic signals from 32 separate pairs of courting flies as well as software for automated segmentation of songs. RESULTS: Our novel hardware design enables recording of low amplitude sounds in most laboratory environments. We demonstrate the power of this system by collecting, segmenting and analyzing over 18 hours of courtship song from 75 males from five wild-type strains of Drosophila melanogaster. Our analysis reveals previously undetected modulation of courtship song features and extensive natural genetic variation for most components of courtship song. Despite having a large dataset with sufficient power to detect subtle modulations of song, we were unable to identify previously reported periodic rhythms in the inter-pulse interval of song. We provide detailed instructions for assembling the hardware and for using our open-source segmentation software. CONCLUSIONS: Analysis of a large dataset of acoustic signals from Drosophila melanogaster provides novel insight into the structure and dynamics of species-specific courtship songs. Our new system for recording and analyzing fly acoustic signals should therefore greatly accelerate future studies of the genetics, neurobiology and evolution of courtship song.