David L. Stern

Many insects manipulate plants by injecting effector proteins. In one extreme example of this molecular “hijacking,” Hormaphis cornu aphids inject bicycle proteins into Hamamelis virginiana , contributing to the development of novel organs called galls. Bicycle proteins share no amino acid sequence similarity with proteins of known function. Here, we report the crystal structures of two divergent bicycle proteins. Both proteins contain saposin-like folds: one with multiple disulfide bonds exhibits a swapped domain topology; the other has no disulfide bonds and possesses two distinct, tandem domains. To explore the structural evolution of bicycle proteins, we attempted to predict bicycle protein structures with Alphafold2 (AF2) and other deep learning programs. While AF2 did not recover the two experimental structures using existing databases, it succeeded when provided with multiple sequence alignments (MSAs) of protein sequences from newly sequenced closely related species. Using this approach, we generated 2,400 high-confidence bicycle protein predictions from seven aphid species. While all aphid bicycle proteins contain predicted saposin-like folds, they display a vast diversity of structural and physicochemical properties. While this diversity thwarts prediction of conserved functions encoded in structure, it suggests that bicycle proteins have evolved to target diverse plant processes and/or to evade plant immune surveillance. Our extension of AF2 with custom MSAs of proteins from closely related species provides a generalizable, powerful approach for predicting structures of rapidly evolving protein families.

Understanding how nervous systems generate coordinated movement requires precise measurement of body kinematics during natural behavior. The fruit fly, Drosophila , is a model organism with sophisticated behavior and well-studied neural circuits, but tracking fly movements in 3D remains challenging because of their teeny bodies, rapid movements, and frequent self-occlusions. Here we present a pipeline for markerless, full-body 3D pose estimation of fly terrestrial behavior, combining seven synchronized high-speed cameras to capture whole-body kinematics at 800 frames per second. We trained a hybrid 2D/3D deep learning model to track 50 keypoints, then refined them to produce anatomically feasible kinematic trajectories through a retargeting process that solved an inverse kinematics problem constrained by a biomechanical body model. Analysis of 3D kinematics revealed that flies perform grounded running across their full speed range, without transitioning between discrete gaits. Using multi-animal tracking, we found that courting males coordinate both wings during song and modulate body pitch to track the female's vertical position. Our open-source pipeline and 3D kinematic dataset of fly behavior provide a foundation for neuromechanical modeling and mechanistic studies of motor control in a genetically tractable model organism.

This protocol is designed for generating and preparing T7 Cas9 mRNA for co-injection into Drosophila. In vitro synthesis of capped mRNA is utilized in co-injections of many systems including Drosophila. This protocol utilizes plasmids with T7 expression system with a T7 promoter. It provides a double cleaning method that ensures a clean mRNA sample for embryo injections. While it specifically refers to Cas9 it may be utilized for other transcription of other plasmids uitlizing the T7 expression system such as T7 phiC31 for attP and attB integration.

Many insects manipulate plants by injecting effector proteins. In one extreme example of this molecular “hijacking”, Hormaphis cornu aphids inject bicycle proteins into Hamamelis virginiana , contributing to the development of novel organs called galls. Bicycle proteins share no amino acid sequence similarity with proteins of known function. Here, we report the crystal structures of two divergent bicycle proteins. Both proteins contain saposin-like folds: one with multiple disulfide bonds exhibits a swapped domain topology; the other has no disulfide bonds and possesses two distinct, tandem domains. To explore the structural evolution of bicycle proteins, we attempted to predict bicycle protein structures with Alphafold2 (AF2) and other deep learning programs. While AF2 did not recover the two experimental structures using existing databases, it succeeded when provided with multiple sequence alignments (MSAs) of protein sequences from newly sequenced closely related species. Using this approach, we generated 2400 high-confidence bicycle protein predictions from seven aphid species. While all aphid bicycle proteins contain predicted saposin-like folds, they display a vast diversity of structural and physicochemical properties. While this diversity thwarts prediction of conserved functions encoded in structure, it suggests that bicycle proteins have evolved to target diverse plant processes and/or to evade plant immune surveillance. Our extension of AF2 with custom MSAs of proteins from closely related species provides a generalizable, powerful approach for predicting structures of rapidly evolving protein families. Significance statement: Parasites introduce specialized “effector” proteins into hosts to suppress host immunity and to release nutrients. The molecular functions and structures of most effector proteins are unknown. Effector proteins often evolve rapidly and share no similarity with proteins of known function. Here, we demonstrate that machine learning algorithms can predict the structures of aphid “bicycle” effector proteins when supplemented with data from closely related species. We exploit this finding to generate predictions of 2400 bicycle protein structures. Aphid bicycle proteins exploit a common folding motif, yet exhibit topologically distinct structures that form separate structural clusters. Despite the clustering of these proteins in structure space, they occupy a nearly uniformly physicochemical space, suggesting that they encode a large diversity of molecular functions.

We have identified a Drosophila species which exhibits spontaneous and robust intermale sexual behavior. D. santomea males distinguish conspecific sexes but court both vigorously and seldom attack. Elevated intermale courtship stems from at least three evolutionarily derived pheromonal changes. In males, the sexually monomorphic cuticular pheromone 7-tricosene promotes rather than inhibits courtship and the courtship-inhibiting olfactory pheromone cVA is reduced 84-92% compared to close relatives. The third switch is in D. santomea females, where cVA suppresses rather than promotes sexual receptivity. Female cVA aversion and male cVA reduction may have co-evolved to maintain efficient intraspecific mating but prevent hybridization with the sympatric sibling species D . yakuba . High intermale courtship and low cVA also co-occur and appear selectively derived in a distant monomorphic species D. persimilis , implying pheromonal and social behavioral convergence. Changes in pheromone valence and levels may therefore explain the recent evolutionary emergence of intermale sexual behavior in Drosophila .

Two invasive adelgids are associated with widespread damage to several North American conifer species. Adelges tsugae, hemlock woolly adelgid, was introduced from Japan and reproduces parthenogenetically in North America, where it has rapidly decimated Tsuga canadensis and Tsuga caroliniana (eastern and Carolina hemlocks, respectively). Adelges abietis, eastern spruce gall adelgid, introduced from Europe, forms distinctive pineapple-shaped galls on several native spruce species. While not considered a major forest pest, it weakens trees and increases susceptibility to additional stressors. Broad-spectrum insecticides that are often used to control adelgid populations can have off-target impacts on beneficial insects. Whole-genome sequencing was performed on both species to aid in the development of targeted solutions that may minimize ecological impact. A. abietis was sequenced using barcoded linked-reads from 30 pooled individuals, with Hi-C scaffolding performed using data from a single individual collected from the same host plant. A. tsugae used long-read sequencing from pooled nymphs. The assembled A. tsugae and A. abietis genomes, pooled from several parthenogenetic females, are 220.75 and 253.16 Mb, respectively. Each consists of 8 autosomal chromosomes, as well as 2 sex chromosomes (X1/X2), supporting the XX-XO sex determination system. The genomes are over 96% complete based on BUSCO assessment. Genome annotation identified 11,424 and 12,060 protein-coding genes in A. tsugae and A. abietis, respectively. Comparative analysis of proteins across 31 hemipteran species and 14 arthropod outgroups identified 32,340 putative gene families. Gene family evolution analysis with CAFE revealed lineage-specific expansions in immune-related aminopeptidases (ERAP1) and juvenile hormone binding proteins (JHBP), contractions in juvenile hormone acid methyltransferases (JHAMT), and conservation of nicotinic acetylcholine receptors (nAChRs). These genes were explored as candidate families toward a long-term objective of developing adelgid-selective insecticides. Structural comparisons of proteins across 7 focal species (A. tsugae, A. abietis, Adelges cooleyi, Rhopalosiphum maidis, Apis mellifera, Danaus plexippus, and Drosophila melanogaster) revealed high conservation of nAChR and ERAP1, while JHAMT exhibited species-specific structural divergence. The potential of JHAMT as a lineage-specific target for pest control was explored through virtual screening of drugs and pesticides.

To perform most behaviors, animals must send commands from higher-order processing centers in the brain to premotor circuits that reside in ganglia distinct from the brain, such as the mammalian spinal cord or insect ventral nerve cord. How these circuits are functionally organized to generate the great diversity of animal behavior remains unclear. An important first step in unraveling the organization of premotor circuits is to identify their constituent cell types and create tools to monitor and manipulate these with high specificity to assess their functions. This is possible in the tractable ventral nerve cord of the fly. To generate such a toolkit, we used a combinatorial genetic technique (split-GAL4) to create 195 sparse transgenic driver lines targeting 196 individual cell types in the ventral nerve cord. These included wing and haltere motoneurons, modulatory neurons, and interneurons. Using a combination of behavioral, developmental, and anatomical analyses, we systematically characterized the cell types targeted in our collection. In addition, we identified correspondences between the cells in this collection and a recent connectomic data set of the ventral nerve cord. Taken together, the resources and results presented here form a powerful toolkit for future investigations of neuronal circuits and connectivity of premotor circuits while linking them to behavioral outputs.

Two invasive hemipteran adelgids are associated with widespread damage to several North American conifer species. Adelges tsugae, hemlock woolly adelgid, was introduced from Japan and reproduces parthenogenetically in North America, where it has rapidly decimated Tsuga canadensis and Tsuga caroliniana (the eastern and Carolina hemlocks, respectively). Adelges abietis , eastern spruce gall adelgid, introduced from Europe, forms distinctive pineapple-shaped galls on several native spruce species. While not considered a major forest pest, it weakens trees and increases susceptibility to additional stressors. Broad-spectrum insecticides that are often used to control adelgid populations can have off-target impacts on beneficial insects. Whole genome sequencing was performed on both species to aid in development of targeted solutions that may minimize ecological impact. Adelges abietis was sequenced using Illumina Linked-Read technology from 30 pooled individuals, with Hi-C scaffolding performed using data from a single individual collected from the same host plant. Adelges tsugae used Oxford Nanopore long-read sequencing from pooled nymphs. The assembled A. tsugae and A. abietis genomes, pooled from several parthenogenetic females, are 220.75 Mbp and 253.16 Mbp, respectively. Each consists of eight autosomal chromosomes, as well as two sex chromosomes (X1/X2), supporting the XX-XO sex determination system. The genomes are over 96% complete based on BUSCO assessment. Genome annotation identified 11,424 and 12,060 protein-coding genes in A. tsugae and A. abietis , respectively. Comparative analysis of proteins across 29 hemipteran species and 14 arthropod outgroups identified 31,666 putative gene families. Gene family evolution analysis with CAFE revealed lineage-specific expansions in immune-related aminopeptidases ( ERAP1 ) and juvenile hormone binding proteins ( JHBP ), contractions in juvenile hormone acid methyltransferases ( JHAMT ), and conservation of nicotinic acetylcholine receptors ( nAChR ). These genes were explored as candidate families towards a long-term objective of developing adelgid-selective insecticides. Structural comparisons of proteins across seven focal species ( Adelges tsugae , Adelges abietis , Adelges cooleyi , Rhopalosiphum maidis , Apis mellifera , Danaus plexippus , and Drosophila melanogaster ) revealed high conservation of nAChR and ERAP1, while JHAMT exhibited species-specific structural divergence. The potential of JHAMT as a lineage-specific target for pest control was explored through virtual drug and pesticide screening.

Motor systems implement diverse motor programs to pattern behavioral sequences, yet how different motor actions are controlled on a moment-by-moment basis remains unclear. Here, we investigated the neural circuit mechanisms underlying the control of distinct courtship songs in Drosophila. Courting males rapidly alternate between two types of song: pulse and sine. By recording calcium signals in the ventral nerve cord in singing flies, we found that one neural population is active during both songs, whereas an expanded neural population, which includes neurons from the first population, is active during pulse song. Brain recordings showed that this nested activation pattern is present in two descending pathways required for singing. Connectomic analysis reveals that these two descending pathways provide structured input to ventral nerve cord neurons in a manner consistent with their activation patterns. These results suggest that nested premotor circuit activity, directed by distinct descending signals, enables rapid switching between motor actions.

The evolutionary expansion of sensory neuron populations detecting important environmental cues is widespread, but functionally enigmatic. We investigated this phenomenon through comparison of homologous olfactory pathways of Drosophila melanogaster and its close relative Drosophila sechellia, an extreme specialist for Morinda citrifolia noni fruit. D. sechellia has evolved species-specific expansions in select, noni-detecting olfactory sensory neuron (OSN) populations, through multigenic changes. Activation and inhibition of defined proportions of neurons demonstrate that OSN number increases contribute to stronger, more persistent, noni-odour tracking behaviour. These expansions result in increased synaptic connections of sensory neurons with their projection neuron (PN) partners, which are conserved in number between species. Surprisingly, having more OSNs does not lead to greater odour-evoked PN sensitivity or reliability. Rather, pathways with increased sensory pooling exhibit reduced PN adaptation, likely through weakened lateral inhibition. Our work reveals an unexpected functional impact of sensory neuron population expansions to explain ecologically-relevant, species-specific behaviour.