We address the problem of inferring the number of independently blinking fluorescent light emitters, when only their combined intensity contributions can be observed at each timepoint. This problem occurs regularly in light microscopy of objects that are smaller than the diffraction limit, where one wishes to count the number of fluorescently labelled subunits. Our proposed solution directly models the photo-physics of the system, as well as the blinking kinetics of the fluorescent emitters as a fully differentiable hidden Markov model. Given a trace of intensity over time, our model jointly estimates the parameters of the intensity distribution per emitter, their blinking rates, as well as a posterior distribution of the total number of fluorescent emitters. We show that our model is consistently more accurate and increases the range of countable subunits by a factor of two compared to current state-of-the-art methods, which count based on autocorrelation and blinking frequency, Further-more, we demonstrate that our model can be used to investigate the effect of blinking kinetics on counting ability, and therefore can inform experimental conditions that will maximize counting accuracy.

Although different animal species often exhibit extensive variation in many behaviors, typically scientists examine one or a small number of behaviors in any single study. Here, we propose a new framework to simultaneously study the evolution of many behaviors. We measured the behavioral repertoire of individuals from six species of fruit flies using unsupervised techniques and identified all stereotyped movements exhibited by each species. We then fit a Generalized Linear Mixed Model to estimate the intra- and inter-species behavioral covariances, and, by using the known phylogenetic relationships among species, we estimated the (unobserved) behaviors exhibited by ancestral species. We found that much of intra-specific behavioral variation has a similar covariance structure to previously described long-time scale variation in an individual's behavior, suggesting that much of the measured variation between individuals of a single species in our assay reflects differences in the status of neural networks, rather than genetic or developmental differences between individuals. We then propose a method to identify groups of behaviors that appear to have evolved in a correlated manner, illustrating how sets of behaviors, rather than individual behaviors, likely evolved. Our approach provides a new framework for identifying co-evolving behaviors and may provide new opportunities to study the mechanistic basis of behavioral evolution.

Regions of genomic DNA called enhancers encode binding sites for transcription factor proteins. Binding of activators and repressors increase and reduce transcription, respectively, but it is not understood how combinations of activators and repressors generate precise patterns of transcription during development. Here, we explore this problem using a fully synthetic transcriptional platform in Drosophila consisting of engineered transcription factor gradients and artificial enhancers. We found that binding sites for a transcription factor that makes DNA accessible are required together with binding sites for transcriptional activators to produce a functional enhancer. Only in this context can changes in the number of activator binding sites mediate quantitative control of transcription. Using an engineered transcriptional repressor gradient, we demonstrate that overlapping repressor and activator binding sites provide more robust repression and sharper expression boundaries than non-overlapping sites. This may explain why this common motif is observed in many developmental enhancers.

Wing dimorphisms have long served as models for examining the ecological and evolutionary tradeoffs associated with alternative phenotypes. Here, we investigated the genetic cause of the pea aphid () male wing dimorphism, wherein males exhibit one of two morphologies that differ in correlated traits that include the presence or absence of wings. We mapped this trait difference to a single genomic region and, using third generation, long-read sequencing, we identified a 120 kb insertion in the wingless allele. This insertion includes a duplicated gene, which is a strong candidate gene in the minimal mapped interval to cause the dimorphism. We found that both alleles were present prior to pea aphid biotype lineage diversification, we estimated that the insertion occurred millions of years ago, and we propose that both alleles have been maintained in the species, likely due to balancing selection.

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.

In an elaborate form of inter-species exploitation, many insects hijack plant development to induce novel plant organs called galls that provide the insect with a source of nutrition and a temporary home. Galls result from dramatic reprogramming of plant cell biology driven by insect molecules, but the roles of specific insect molecules in gall development have not yet been determined. Here, we study the aphid Hormaphis cornu, which makes distinctive "cone" galls on leaves of witch hazel Hamamelis virginiana. We found that derived genetic variants in the aphid gene determinant of gall color (dgc) are associated with strong downregulation of dgc transcription in aphid salivary glands, upregulation in galls of seven genes involved in anthocyanin synthesis, and deposition of two red anthocyanins in galls. We hypothesize that aphids inject DGC protein into galls and that this results in differential expression of a small number of plant genes. dgc is a member of a large, diverse family of novel predicted secreted proteins characterized by a pair of widely spaced cysteine-tyrosine-cysteine (CYC) residues, which we named BICYCLE proteins. bicycle genes are most strongly expressed in the salivary glands specifically of galling aphid generations, suggesting that they may regulate many aspects of gall development. bicycle genes have experienced unusually frequent diversifying selection, consistent with their potential role controlling gall development in a molecular arms race between aphids and their host plants.

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.

To establish functional connectivity between two candidate neurons that might form a circuit element, a common approach is to activate an optogenetic tool such as Chrimson in the candidate pre-synaptic neuron and monitor fluorescence of the calcium-sensitive indicator GCaMP in a candidate post-synaptic neuron. While performing such experiments, we found that low levels of leaky Chrimson expression can lead to strong artifactual GCaMP signals in presumptive postsynaptic neurons even when Chrimson is not intentionally expressed in any particular neurons. Withholding all-trans retinal, the chromophore required as a co-factor for Chrimson response to light, eliminates GCaMP signal but does not provide an experimental control for leaky Chrimson expression. Leaky Chrimson expression appears to be an inherent feature of current Chrimson transgenes, since artifactual connectivity was detected with Chrimson transgenes integrated into three different genomic locations (two insertions tested in larvae; a third insertion tested in the adult fly). These false-positive signals may complicate the interpretation of functional connectivity experiments. We illustrate how a no-Gal4 negative control improves interpretability of functional connectivity assays. We also propose a simple but effective procedure to identify experimental conditions that minimize potentially incorrect interpretations caused by leaky Chrimson expression.

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.

Chromosome inversions are of unique importance in the evolution of genomes and species because when heterozygous with a standard arrangement chromosome, they suppress meiotic crossovers within the inversion. In Drosophila species, heterozygous inversions also cause the interchromosomal effect, whereby the presence of a heterozygous inversion induces a dramatic increase in crossover frequencies in the remainder of the genome within a single meiosis. To date, the interchromosomal effect has been studied exclusively in species that also have high frequencies of inversions in wild populations. We took advantage of a recently developed approach for generating inversions in Drosophila simulans, a species that does not have inversions in wild populations, to ask if there is an interchromosomal effect. We used the existing chromosome 3R balancer and generated a new chromosome 2L balancer to assay for the interchromosomal effect genetically and cytologically. We found no evidence of an interchromosomal effect in D. simulans. To gain insight into the underlying mechanistic reasons, we qualitatively analyzed the relationship between meiotic double-strand break formation and synaptonemal complex assembly. We find that the synaptonemal complex is assembled prior to double-strand break formation as in D. melanogaster; however, we show that the synaptonemal complex is assembled prior to localization of the oocyte determination factor Orb, whereas in D. melanogaster, synaptonemal complex formation does not begin until Orb is localized. Together, our data show heterozygous inversions in D. simulans do not induce an interchromosomal effect and that there are differences in the developmental programming of the early stages of meiosis.