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.

It is unclear how regulatory genes establish neural circuits that compose sex-specific behaviors. The Drosophila melanogaster male courtship song provides a powerful model to study this problem. Courting males vibrate a wing to sing bouts of pulses and hums, called pulse and sine song, respectively. We report the discovery of male-specific thoracic interneurons—the TN1A neurons—that are required specifically for sine song. The TN1A neurons can drive the activity of a sex-non-specific wing motoneuron, hg1, which is also required for sine song. The male-specific connection between the TN1A neurons and the hg1 motoneuron is regulated by the sexual differentiation gene doublesex. We find that doublesex is required in the TN1A neurons during development to increase the density of the TN1A arbors that interact with dendrites of the hg1motoneuron. Our findings demonstrate how a sexual differentiation gene can build a sex-specific circuit motif by modulating neuronal arborization.

•Doublesex-expressing TN1 neurons are necessary and sufficient for the male sine song•A subclass of TN1 neurons, TN1A, contributes to the sine song•TN1A neurons are functionally coupled to a sine song motoneuron, hg1•Doublesex regulates the connectivity between the TN1A and hg1 neurons

It is unclear how developmental regulatory genes specify sex-specific behaviors. Shirangi et al. demonstrate that the Drosophila sexual differentiation gene doublesex encodes a sex-specific behavior—male song—by promoting the connectivity between the male-specific TN1A neurons and the sex-non-specific hg1 neurons, which are required for production of the song.

Interindividual differences in neuronal wiring may contribute to behavioral individuality and affect susceptibility to neurological disorders. To investigate the causes and potential consequences of wiring variation in Drosophila melanogaster, we focused on a hemilineage of ventral nerve cord interneurons that exhibits morphological variability. We find that late-born subclasses of the 12A hemilineage are highly sensitive to genetic and environmental variation. Neurons in the second thoracic segment are particularly variable with regard to two developmental decisions, whereas its segmental homologs are more robust. This variability "hotspot" depends on Ultrabithorax expression in the 12A neurons, indicating variability is cell-intrinsic and under genetic control. 12A development is more variable and sensitive to temperature in long-established laboratory strains than in strains recently derived from the wild. Strains with a high frequency of one of the 12A variants also showed a high frequency of animals with delayed spontaneous flight initiation, whereas other wing-related behaviors did not show such a correlation and were thus not overtly affected by 12A variation. These results show that neurodevelopmental robustness is variable and under genetic control in Drosophila and suggest that the fly may serve as a model for identifying conserved gene pathways that stabilize wiring in stressful developmental environments. Moreover, some neuronal lineages are variation hotspots and thus may be more amenable to evolutionary change.

Genes are regulated by transcription factors that bind to regions of genomic DNA called enhancers. Considerable effort is focused on identifying transcription factor binding sites, with the goal of predicting gene expression from DNA sequence. Despite this effort, general, predictive models of enhancer function are currently lacking. Here we combine quantitative models of enhancer function with manipulations using engineered transcription factors to examine the extent to which enhancer function can be controlled in a quantitatively predictable manner. Our models, which incorporate few free parameters, can accurately predict the contributions of ectopic transcription factor inputs. These models allow the predictable 'tuning' of enhancers, providing a framework for the quantitative control of enhancers with engineered transcription factors.

Multiple methods have been introduced over the past 30 years to identify the genomic insertion sites of transposable elements and other DNA elements that integrate into genomes. However, each of these methods suffer from limitations that can frustrate attempts to map multiple insertions in a single genome and to map insertions in genomes of high complexity that contain extensive repetitive DNA. I introduce a new method for transposon mapping that is simple to perform, can accurately map multiple insertions per genome, and generates long sequence reads that facilitate mapping to complex genomes. The method, called TagMap, for Tagmentation-based Mapping, relies on a modified Tn5 tagmentation protocol with a single tagmentation adaptor followed by PCR using primers specific to the tranposable element and the adaptor sequence. Several minor modifications to normal tagmentation reagents and protocols allow easy and rapid preparation of TagMap libraries. Short read sequencing starting from the adaptor sequence generates oriented reads that flank and are oriented toward the transposable element insertion site. The convergent orientation of adjacent reads at the insertion site allows straightforward prediction of the precise insertion site(s). A Linux shell script is provided to identify insertion sites from fastq files.

Transcription factor proteins regulate gene expression by binding to specific DNA regions. Most studies of transcription factor binding sites have focused on the highest affinity sites for each factor. There is abundant evidence, however, that binding sites with a range of affinities, including very low affinities, are critical to gene regulation. Here, we present the theoretical and experimental evidence for the importance of low-affinity sites in gene regulation and development. We also discuss the implications of the widespread use of low-affinity sites in eukaryotic genomes for robustness, precision, specificity, and evolution of gene regulation.

Rapid evolution of genitalia shape, a widespread phenomenon in animals with internal fertilization, offers the opportunity to dissect the genetic architecture of morphological evolution linked to sexual selection and speciation. Most quantitative trait loci (QTL) mapping studies of genitalia divergence have focused on Drosophila melanogaster and its three most closely related species, D. simulans, D. mauritiana, and D. sechellia, and have suggested that the genetic basis of genitalia evolution involves many loci. We report the first genetic study of male genitalia evolution between D. yakuba and D. santomea, two species of the D. melanogaster species subgroup. We focus on male ventral branches, which harm females during interspecific copulation. Using landmark-based geometric morphometrics, we characterized shape variation in parental species, F1 hybrids, and backcross progeny and show that the main axis of shape variation within the backcross population matches the interspecific variation between parental species. For genotyping, we developed a new molecular method to perform multiplexed shotgun genotyping (MSG), which allowed us to prepare genomic DNA libraries from 365 backcross individuals in a few days using little DNA. We detected only three QTL, one of which spans 2.7 Mb and exhibits a highly significant effect on shape variation that can be linked to the harmfulness of the ventral branches. We conclude that the genetic architecture of genitalia morphology divergence may not always be as complex as suggested by previous studies.

Male sexual characters are often among the first traits to diverge between closely related species and identifying the genetic basis of such changes can contribute to our understanding of their evolutionary history. However, little is known about the genetic architecture or the specific genes underlying the evolution of male genitalia. The morphology of the claspers, posterior lobes and anal plates exhibit striking differences between Drosophila mauritiana and Drosophila simulans. Using QTL and introgression-based high-resolution mapping, we identified several small regions on chromosome arms 3L and 3R that contribute to differences in these traits. However, we found that the loci underlying the evolution of clasper differences between these two species are independent from those that contribute to posterior lobe and anal plate divergence. Furthermore, while most of the loci affect each trait in the same direction and act additively, we also found evidence for epistasis between loci for clasper bristle number. In addition, we conducted an RNAi screen in D. melanogaster to investigate if positional and expression candidate genes located on chromosome 3L, are also involved in genital development. We found that six of these genes, including components of Wnt signaling and male-specific lethal 3 (msl3), regulate the development of genital traits consistent with the effects of the introgressed regions where they are located and that thus represent promising candidate genes for the evolution these traits.

Multilocus sequence typing (MLST) has become the preferred method for genotyping many biological species, and it is especially useful for analyzing haploid eukaryotes. MLST is rigorous, reproducible, and informative, and MLST genotyping has been shown to identify major phylogenetic clades, molecular groups, or subpopulations of a species, as well as individual strains or clones. MLST molecular types often correlate with important phenotypes. Conventional MLST involves the extraction of genomic DNA and the amplification by PCR of several conserved, unlinked gene sequences from a sample of isolates of the taxon under investigation. In some cases, as few as three loci are sufficient to yield definitive results. The amplicons are sequenced, aligned, and compared by phylogenetic methods to distinguish statistically significant differences among individuals and clades. Although MLST is simpler, faster, and less expensive than whole genome sequencing, it is more costly and time-consuming than less reliable genotyping methods (e.g. amplified fragment length polymorphisms). Here, we describe a new MLST method that uses next-generation sequencing, a multiplexing protocol, and appropriate analytical software to provide accurate, rapid, and economical MLST genotyping of 96 or more isolates in single assay. We demonstrate this methodology by genotyping isolates of the well-characterized, human pathogenic yeast Cryptococcus neoformans.

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.