Regulation of body pigmentation by the Abdominal-B Hox protein and its gain and loss in Drosophila evolution

Molecular signature of primate astrocytes reveals pathways and regulatory changes contributing to human brain evolution

Astrocytes contribute to the development and regulation of the higher-level functions of the brain, the critical targets of evolution. However, how astrocytes evolve in primates is unsettled. Here, we obtain human, chimpanzee, and macaque induced pluripotent stem-cell-derived astrocytes (iAstrocytes). Human iAstrocytes are bigger and more complex than the non-human primate iAstrocytes. We identify new loci contributing to the increased human astrocyte. We show that genes and pathways implicated in long-range intercellular signaling are activated in the human iAstrocytes and partake in controlling iAstrocyte complexity. Genes downregulated in human iAstrocytes frequently relate to neurological disorders and were decreased in adult brain samples. Through regulome analysis and machine learning, we uncover that functional activation of enhancers coincides with a previously unappreciated, pervasive gain of "stripe" transcription factor binding sites. Altogether, we reveal the transcriptomic signature of primate astrocyte evolution and a mechanism driving the acquisition of the regulatory potential of enhancers.

Protocol for dissecting Drosophila pupae and visualizing RNA expression using hybridization chain reaction

Visualizing RNA expression in the Drosophila epidermis during pupal development is challenging because the tissue is fragile during early pupal stages and increasingly impermeable at later stages. Here, we present a protocol for tissue dissection and detection of RNA in situ. We describe steps for using the hybridization chain reaction (HCR) in early and late Drosophila pupal stages (and larval imaginal discs). This protocol facilitates the study of dynamically changing patterns in gene expression during pupal development.

Entangled and non-modular enhancer sequences producing independent spatial activities

The modularity of transcriptional enhancers is central to our understanding of morphological evolution, allowing specific changes to a gene expression pattern component, without affecting others. Enhancer modularity refers to physically separated stretches of regulatory sequence producing discrete spatiotemporal transcriptional activity. This concept stems from assays that test the sufficiency of a DNA segment to drive spatial reporter expression resembling that of the corresponding gene. Focusing on spatial patterns, it overlooks quantitative aspects of gene expression, underestimating the regulatory sequence actually required to reach full endogenous expression levels. Here, we show that five regulatory activities of the gene yellow in Drosophila, classically described as modular, result from extensively overlapping sequences, with broadly distributed regulatory information. Nevertheless, the independent regulatory activities of these entangled enhancers appear to be nucleated by specific segments that we called enhancer cores. Our work calls for a reappraisal of enhancer definition and properties, as well as of the consequences on regulatory evolution.

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

Characterization of Single-Cell Cis-regulatory Elements Informs Implications for Cell Differentiation

Cis-regulatory elements govern the specific patterns and dynamics of gene expression in cells during development, which are the fundamental mechanisms behind cell differentiation. However, the genomic characteristics of single-cell cis-regulatory elements closely linked to cell differentiation during development remain unclear. To explore this, we systematically analyzed ∼250,000 putative single-cell cis-regulatory elements obtained from snATAC-seq analysis of the developing mouse cerebellum. We found that over 80% of these single-cell cis-regulatory elements show pleiotropic effects, being active in 2 or more cell types. The pleiotropic degrees of proximal and distal single-cell cis-regulatory elements are positively correlated with the density and diversity of transcription factor binding motifs and GC content. There is a negative correlation between the pleiotropic degrees of single-cell cis-regulatory elements and their distances to the nearest transcription start sites, and proximal single-cell cis-regulatory elements display higher relevance strengths than distal ones. Furthermore, both proximal and distal single-cell cis-regulatory elements related to cell differentiation exhibit enhanced sequence-level evolutionary conservation, increased density and diversity of transcription factor binding motifs, elevated GC content, and greater distances from their nearest genes. Together, our findings reveal the general genomic characteristics of putative single-cell cis-regulatory elements and provide insights into the genomic and evolutionary mechanisms by which single-cell cis-regulatory elements regulate cell differentiation during development.

A genetic screen of transcription factors in the Drosophila melanogaster abdomen identifies novel pigmentation genes

Gene regulatory networks specify the gene expression patterns needed for traits to develop. Differences in these networks can result in phenotypic differences between organisms. Although loss-of-function genetic screens can identify genes necessary for trait formation, gain-of-function screens can overcome genetic redundancy and identify loci whose expression is sufficient to alter trait formation. Here, we leveraged transgenic lines from the Transgenic RNAi Project at Harvard Medical School to perform both gain- and loss-of-function CRISPR/Cas9 screens for abdominal pigmentation phenotypes. We identified measurable effects on pigmentation patterns in the Drosophila melanogaster abdomen for 21 of 55 transcription factors in gain-of-function experiments and 7 of 16 tested by loss-of-function experiments. These included well-characterized pigmentation genes, such as bab1 and dsx, and transcription factors that had no known role in pigmentation, such as slp2. Finally, this screen was partially conducted by undergraduate students in a Genetics Laboratory course during the spring semesters of 2021 and 2022. We found this screen to be a successful model for student engagement in research in an undergraduate laboratory course that can be readily adapted to evaluate the effect of hundreds of genes on many different Drosophila traits, with minimal resources.

The H3K79me3 methyl-transferase Grappa is involved in the establishment and thermal plasticity of abdominal pigmentation in Drosophila melanogaster females

Temperature sensitivity of abdominal pigmentation in Drosophila melanogaster females allows to investigate the mechanisms underlying phenotypic plasticity. Thermal plasticity of pigmentation is due to modulation of tan and yellow expression, encoding pigmentation enzymes. Furthermore, modulation of tan expression by temperature is correlated to the variation of the active histone mark H3K4me3 on its promoter. Here, we test the role of the DotCom complex, which methylates H3K79, another active mark, in establishment and plasticity of pigmentation. We show that several components of the DotCom complex are involved in the establishment of abdominal pigmentation. In particular, Grappa, the catalytic unit of this complex, plays opposite roles on pigmentation at distinct developmental stages. Indeed, its down-regulation from larval L2 to L3 stages increases female adult pigmentation, whereas its down-regulation during the second half of the pupal stage decreases adult pigmentation. These opposite effects are correlated to the regulation of distinct pigmentation genes by Grappa: yellow repression for the early role and tan activation for the late one. Lastly, reaction norms measuring pigmentation along temperature in mutants for subunits of the DotCom complex reveal that this complex is not only involved in the establishment of female abdominal pigmentation but also in its plasticity.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

The melanin pigment gene black mediates body pigmentation and courtship behaviour in the German cockroach Blattella germanica

Genes involved in melanin production directly impact insect pigmentation and can affect diverse physiology and behaviours. The role these genes have on sex behaviour, however, is unclear. In the present study, the crucial melanin pigment gene black was functionally characterised in an urban pest, the German cockroach, Blattella germanica. RNAi knockdown of B. germanica black (Bgblack) had no effect on survival, but did result in black pigmentation of the thoraxes, abdomens, heads, wings, legs, antennae, and cerci due to cuticular accumulation of melanin. Sex-specific variation in the pigmentation pattern was apparent, with females exhibiting darker coloration on the abdomen and thorax than males. Bgblack knockdown also resulted in wing deformation and negatively impacted the contact sex pheromone-based courtship behaviour of males. This study provides evidence for black function in multiple aspects of B. germanica biology and opens new avenues of exploration for novel pest control strategies.

Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

Mask exhibits trxG-like behavior and associates with H3K27ac marked chromatin

The Trithorax group (trxG) proteins counteract the repressive effect of Polycomb group (PcG) complexes and maintain transcriptional memory of active states of key developmental genes. Although chromatin structure and modifications appear to play a fundamental role in this process, it is not clear how trxG prevents PcG-silencing and heritably maintains an active gene expression state. Here, we report a hitherto unknown role of Drosophila Multiple ankyrin repeats single KH domain (Mask), which emerged as one of the candidate trxG genes in our reverse genetic screen. The genome-wide binding profile of Mask correlates with known trxG binding sites across the Drosophila genome. In particular, the association of Mask at chromatin overlaps with CBP and H3K27ac, which are known hallmarks of actively transcribed genes by trxG. Importantly, Mask predominantly associates with actively transcribed genes in Drosophila. Depletion of Mask not only results in the downregulation of trxG targets but also correlates with diminished levels of H3K27ac. The fact that Mask positively regulates H3K27ac levels in flies was also found to be conserved in human cells. Strong suppression of Pc mutant phenotype by mutation in mask provides physiological relevance that Mask contributes to the anti-silencing effect of trxG, maintaining expression of key developmental genes. Since Mask is a downstream effector of multiple cell signaling pathways, we propose that Mask may connect cell signaling with chromatin mediated epigenetic cell memory governed by trxG.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Changes in locus wide repression underlie the evolution of Drosophila abdominal pigmentation

Changes in gene regulation represent an important path to generate developmental differences affecting anatomical traits. Interspecific divergence in gene expression often results from changes in transcription-stimulating enhancer elements. While gene repression is crucial for precise spatiotemporal expression patterns, the relative contribution of repressive transcriptional silencers to regulatory evolution remains to be addressed. Here, we show that the Drosophila pigmentation gene ebony has mainly evolved through changes in the spatial domains of silencers patterning its abdominal expression. By precisely editing the endogenous ebony locus of D. melanogaster, we demonstrate the requirement of two redundant abdominal enhancers and three silencers that repress the redundant enhancers in a patterned manner. We observe a role for changes in these silencers in every case of ebony evolution observed to date. Our findings suggest that negative regulation by silencers likely has an under-appreciated role in gene regulatory evolution.

Evolution of modular and pleiotropic enhancers

Cis-regulatory elements (CREs), or enhancers, are segments of noncoding DNA that regulate the spatial and temporal expression of nearby genes. Sometimes, genes are expressed in more than one tissue, and this can be driven by two main types of CREs: tissue-specific "modular" CREs, where different CREs drive expression of the gene in the different tissues, or by "pleiotropic" CREs, where the same CRE drives expression in the different tissues. In this perspective, we will discuss some of the ways (i) modular and pleiotropic CREs might originate; (ii) propose that modular CREs might derive from pleiotropic CREs via a process of duplication, degeneration, and complementation (the CRE-DDC model); and (iii) propose that hotspot loci of evolution are associated with the origin of modular CREs belonging to any gene in a regulatory network.

Widespread cis- and trans-regulatory evolution underlies the origin, diversification, and loss of a sexually dimorphic fruit fly pigmentation trait

Changes in gene expression are a prominent feature of morphological evolution. These changes occur to hierarchical gene regulatory networks (GRNs) of transcription factor genes that regulate the expression of trait-building differentiation genes. While changes in the expression of differentiation genes are essential to phenotypic evolution, they can be caused by mutations within cis-regulatory elements (CREs) that drive their expression (cis-evolution) or within genes for CRE-interacting transcription factors (trans-evolution). Locating these mutations remains a challenge, especially when experiments are limited to one species that possesses the ancestral or derived phenotype. We investigated CREs that control the expression of the differentiation genes tan and yellow, the expression of which evolved during the gain, modification, and loss of dimorphic pigmentation among Sophophora fruit flies. We show these CREs to be necessary components of a pigmentation GRN, as deletion from Drosophila melanogaster (derived dimorphic phenotype) resulted in lost expression and lost male-specific pigmentation. We evaluated the ability of orthologous CRE sequences to drive reporter gene expression in species with modified (Drosophila auraria), secondarily lost (Drosophila ananassae), and ancestrally absent (Drosophila willistoni) pigmentation. We show that the transgene host frequently determines CRE activity, implicating trans-evolution as a significant factor for this trait's diversity. We validated the gain of dimorphic Bab transcription factor expression as a trans-change contributing to the dimorphic trait. Our findings suggest an amenability to change for the landscape of trans-regulators and begs for an explanation as to why this is so common compared to the evolution of differentiation gene CREs.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

A novel role for trithorax in the gene regulatory network for a rapidly evolving fruit fly pigmentation trait

Animal traits develop through the expression and action of numerous regulatory and realizator genes that comprise a gene regulatory network (GRN). For each GRN, its underlying patterns of gene expression are controlled by cis-regulatory elements (CREs) that bind activating and repressing transcription factors. These interactions drive cell-type and developmental stage-specific transcriptional activation or repression. Most GRNs remain incompletely mapped, and a major barrier to this daunting task is CRE identification. Here, we used an in silico method to identify predicted CREs (pCREs) that comprise the GRN which governs sex-specific pigmentation of Drosophila melanogaster. Through in vivo assays, we demonstrate that many pCREs activate expression in the correct cell-type and developmental stage. We employed genome editing to demonstrate that two CREs control the pupal abdomen expression of trithorax, whose function is required for the dimorphic phenotype. Surprisingly, trithorax had no detectable effect on this GRN's key trans-regulators, but shapes the sex-specific expression of two realizator genes. Comparison of sequences orthologous to these CREs supports an evolutionary scenario where these trithorax CREs predated the origin of the dimorphic trait. Collectively, this study demonstrates how in silico approaches can shed novel insights on the GRN basis for a trait's development and evolution.

The Genetic Mechanisms Underlying the Concerted Expression of the yellow and tan Genes in Complex Patterns on the Abdomen and Wings of Drosophila guttifera

How complex morphological patterns form is an intriguing question in developmental biology. However, the mechanisms that generate complex patterns remain largely unknown. Here, we sought to identify the genetic mechanisms that regulate the tan (t) gene in a multi-spotted pigmentation pattern on the abdomen and wings of Drosophila guttifera. Previously, we showed that yellow (y) gene expression completely prefigures the abdominal and wing pigment patterns of this species. In the current study, we demonstrate that the t gene is co-expressed with the y gene in nearly identical patterns, both transcripts foreshadowing the adult abdominal and wing melanin spot patterns. We identified cis-regulatory modules (CRMs) of t, one of which drives reporter expression in six longitudinal rows of spots on the developing pupal abdomen, while the second CRM activates the reporter gene in a spotted wing pattern. Comparing the abdominal spot CRMs of y and t, we found a similar composition of putative transcription factor binding sites that are thought to regulate the complex expression patterns of both terminal pigmentation genes y and t. In contrast, the y and t wing spots appear to be regulated by distinct upstream factors. Our results suggest that the D. guttifera abdominal and wing melanin spot patterns have been established through the co-regulation of y and t, shedding light on how complex morphological traits may be regulated through the parallel coordination of downstream target genes.

The regulation of a pigmentation gene in the formation of complex color patterns in Drosophila abdomens

Changes in the control of developmental gene expression patterns have been implicated in the evolution of animal morphology. However, the genetic mechanisms underlying complex morphological traits remain largely unknown. Here we investigated the molecular mechanisms that induce the pigmentation gene yellow in a complex color pattern on the abdomen of Drosophila guttifera. We show that at least five developmental genes may collectively activate one cis-regulatory module of yellow in distinct spot rows and a dark shade to assemble the complete abdominal pigment pattern of Drosophila guttifera. One of these genes, wingless, may play a conserved role in the early phase of spot pattern development in several species of the quinaria group. Our findings shed light on the evolution of complex animal color patterns through modular changes of gene expression patterns.

Genome-Wide Identification, Characterization, and Expression Analysis of Glutamate Receptor-like Gene (GLR) Family in Sugarcane

The plant glutamate receptor-like gene (GLR) plays a vital role in development, signaling pathways, and in its response to environmental stress. However, the GLR gene family has not been comprehensively and systematically studied in sugarcane. In this work, 43 GLR genes, including 34 in Saccharum spontaneum and 9 in the Saccharum hybrid cultivar R570, were identified and characterized, which could be divided into three clades (clade I, II, and III). They had different evolutionary mechanisms, the former was mainly on the WGD/segmental duplication, while the latter mainly on the proximal duplication. Those sugarcane GLR proteins in the same clade had a similar gene structure and motif distribution. For example, 79% of the sugarcane GLR proteins contained all the motifs, which proved the evolutionary stability of the sugarcane GLR gene family. The diverse cis-acting regulatory elements indicated that the sugarcane GLRs may play a role in the growth and development, or under the phytohormonal, biotic, and abiotic stresses. In addition, GO and KEGG analyses predicted their transmembrane transport function. Based on the transcriptome data, the expression of the clade III genes was significantly higher than that of the clade I and clade II. Furthermore, qRT-PCR analysis demonstrated that the expression of the SsGLRs was induced by salicylic acid (SA) treatment, methyl jasmonic acid (MeJA) treatment, and abscisic acid (ABA) treatment, suggesting their involvement in the hormone synthesis and signaling pathway. Taken together, the present study should provide useful information on comparative genomics to improve our understanding of the GLR genes and facilitate further research on their functions.

Transgenerational effects of gamma radiation dose and dose rate on Drosophila flies irradiated at an early embryonal stage

Ionizing radiation (IR) kills cells mainly through induction of DNA damages and the surviving cells may suffer from mutations. Transgenerational effects of IR are well documented, but the exact mechanisms underlying them are less well understood; they include induction of mutations in germ cells and epigenetic inheritance. Previously, effects in the offspring of mice and zebrafish exposed to IR have been reported. A few studies also showed indications of transgenerational effects of radiation in humans, particularly in nuclear power workers. In the present project, short- and long-term effects of low-dose-rate (LDR; 50 and 97 mGy/h) and high-dose-rate (HDR; 23.4, 47.1 and 495 Gy/h) IR in Drosophila embryos were investigated. The embryos were irradiated at different doses and dose rates and radiosensitivity at different developmental stages was investigated. Also, the survival of larvae, pupae and adults developed from embryos irradiated at an early stage (30 min after egg laying) were studied. The larval crawling and pupation height assays were applied to investigate radiation effects on larval locomotion and pupation behavior, respectively. In parallel, the offspring from 3 Gy irradiated early-stage embryos were followed up to 12 generations and abnormal phenotypes were studied. Acute exposure of embryos at different stages of development showed that the early stage embryo is the most sensitive. The effects on larval locomotion showed no significant differences between the dose rates but a significant decrease of locomotion activity above 7 Gy was observed. The results indicate that embryos exposed to the low dose rates have shorter eclosion times. At the same cumulative dose (1 up to 7 Gy), HDR is more embryotoxic than LDR. We also found a radiation-induced depigmentation on males (A5 segment of the dorsal abdomen, A5pig^-) that can be transmitted up to 12 generations. The phenomenon does not follow the classical Mendelian laws of segregation.

The chloroplast genome of Salix floderusii and characterization of chloroplast regulatory elements

Salix floderusii is a rare alpine tree species in the Salix genus. Unfortunately, no extensive germplasm identification, molecular phylogeny, and chloroplast genomics of this plant have been conducted. We sequenced the chloroplast (cp) genome of S. floderusii for the first time using second-generation sequencing technology. The cp genome was 155,540 bp long, including a large single-copy region (LSC, 84,401 bp), a small single-copy region (SSC, 16,221 bp), and inverted repeat regions (IR, 54,918 bp). A total of 131 genes were identified, including 86 protein genes, 37 tRNA genes, and 8 rRNA genes. The S. floderusii cp genome contains 1 complement repeat, 24 forward repeats, 17 palindromic repeats, and 7 reverse repeats. Analysis of the IR borders showed that the IRa and IRb regions of S. floderusii and Salix caprea were shorter than those of Salix cinerea, which may affect plastome evolution. Furthermore, four highly variable regions were found, including the rpl22 coding region, psbM/trnD-GUC non-coding region, petA/psbJ non-coding region, and ycf1 coding region. These high variable regions can be used as candidate molecular markers and as a reference for identifying future Salix species. In addition, phylogenetic analysis indicated that the cp genome of S. floderusii is sister to Salix cupularis and belongs to the Subgenus Vetrix. Genes (Sf-trnI, Sf-PpsbA, aadA, Sf-TpsbA, Sf-trnA) obtained via cloning were inserted into the pBluescript II SK (+) to yield the cp expression vectors, which harbored the selectable marker gene aadA. The results of a spectinomycin resistance test indicated that the cp expression vector had been successfully constructed. Moreover, the aadA gene was efficiently expressed under the regulation of predicted regulatory elements. The present study provides a solid foundation for establishing subsequent S. floderusii cp transformation systems and developing strategies for the genetic improvement of S. floderusii.

The role of the epidermis enhancer element in positive and negative transcriptional regulation of ebony in Drosophila melanogaster

The spatiotemporal regulation of gene expression is essential to ensure robust phenotypic outcomes. Pigmentation patterns in Drosophila are determined by pigments biosynthesized in the developing epidermis and the cis-regulatory elements of the genes involved in this process are well-characterized. Here, we report that the known primary epidermal enhancer is dispensable for the transcriptional activation of ebony (involved in light-colored pigment synthesis) in the developing epidermis of Drosophila melanogaster. The evidence was obtained by introducing an approximately 1 kbp deletion at the primary epidermal enhancer by genome editing. The effect of the primary epidermal enhancer deletion on pigmentation and on the endogenous expression pattern of a mCherry-fused ebony allele was examined in the abdomen. The expression levels of the mCherry-fused ebony in the primary epidermal enhancer-deleted strains were slightly higher than that of the control strain, indicating that the sequences outside the primary epidermal enhancer have an ability to drive an expression of this gene in the epidermis. Interestingly, the primary epidermal enhancer deletion resulted in a derepression of this gene in the dorsal midline of the abdominal tergites, where dark pigmentation is present in the wild-type individuals. This indicated that the primary epidermal enhancer fragment contains a silencer. Furthermore, the endogenous expression pattern of ebony in the 2 additional strains with partially deleted primary epidermal enhancer revealed that the silencer resides within a 351-bp fragment in the 5' portion of the primary epidermal enhancer. These results demonstrated that deletion assays combined with reporter assays are highly effective in detecting the presence of positively and negatively regulating sequences within and outside the focal cis-regulatory elements.

Micromanagement of Drosophila Post-Embryonic Development by Hox Genes

Hox genes function early in development to determine regional identity in animals. Consequently, the loss or gain of Hox gene expression can change this identity and cause homeotic transformations. Over 20 years ago, it was observed that the role of Hox genes in patterning animal body plans involves the fine-scale regulation of cell fate and identity during development, playing the role of 'micromanagers' as proposed by Michael Akam in key perspective papers. Therefore, as well as specifying where structures develop on animal bodies, Hox genes can help to precisely sculpt their morphology. Here, we review work that has provided important insights about the roles of Hox genes in influencing cell fate during post-embryonic development in Drosophila to regulate fine-scale patterning and morphology. We also explore how this is achieved through the regulation of Hox genes, specific co-factors and their complex regulation of hundreds of target genes. We argue that further investigating the regulation and roles of Hox genes in Drosophila post-embryonic development has great potential for understanding gene regulation, cell fate and phenotypic differentiation more generally.

Homeodomain proteins POU-M2, antennapedia and abdominal-B are involved in regulation of the segment-specific expression of the clip-domain serine protease gene CLIP13 in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs) play important roles in insect innate immunity and development. Our previous studies indicated that CLIP13, an epidermis-specific gene, was involved in cuticle remodeling during molting and metamorphosis in the silkworm, Bombyx mori. However, the transcriptional regulatory mechanism and regulatory pathways of CLIP13 remained unclear. In the present study, we investigated CLIP13 expression and the regulation pathway controlled by 20-hydroxyecdysone (20E) in the silkworm. At the transcriptional level, expression of CLIP13 exhibited pronounced spatial and temporal specificity in different regions of the epidermis; homeodomain transcription factors POU-M2, antennapedia (Antp), and abdominal-B (Abd-B) showed similar expression change trends as CLIP13 in the head capsule, thorax, and abdomen, respectively. Furthermore, results of cell transfection assays, electrophoretic mobility shift assays, and chromatin immunoprecipitation demonstrated that POU-M2, Antp, and Abd-B were involved in the transcriptional regulation of CLIP13 by directly binding to their cis-response elements in CLIP13 promoter. RNA interference-mediated silencing of POU-M2, Antp, and Abd-B led to a decrease of CLIP13 expression in the head capsule, the epidermis of the 1st to 3rd thoracic segments and the 7th to 10th abdominal segments, respectively. Consistent with CLIP13, 20E treatment significantly upregulated expression of POU-M2, Antp, and Abd-B in the silkworm epidermis. Taken together, these data suggest that 20E positively regulates transcription of CLIP13 via homeodomain proteins POU-M2, Antp, and Abd-B in different regions of the silkworm epidermis during metamorphosis, thus affecting the molting process. Our findings provide new insight into the functions of homeodomain transcription factors in insect molting.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

Rapid radiation in a highly diverse marine environment

Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.

Repression precedes independent evolutionary gains of a highly specific gene expression pattern

Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns.

Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.

Redundant enhancers in the iab-5 domain cooperatively activate Abd-B in the A5 and A6 abdominal segments of Drosophila

The Abdominal-B (Abd-B) gene belongs to the bithorax complex and its expression is controlled by four regulatory domains, iab-5, iab-6, iab-7 and iab-8, each of which is thought to be responsible for directing the expression of Abd-B in one of the abdominal segments from A5 to A8. A variety of experiments have supported the idea that BX-C regulatory domains are functionally autonomous and that each domain is both necessary and sufficient to orchestrate the development of the segment they specify. Unexpectedly, we discovered that this model does not always hold. Instead, we find that tissue-specific enhancers located in the iab-5 domain are required for the proper activation of Abd-B not only in A5 but also in A6. Our findings indicate that the functioning of the iab-5 and iab-6 domains in development of the adult cuticle A5 and A6 in males fit better with an additive model, much like that first envisioned by Ed Lewis.

Hox genes are essential for the development of eyespots in Bicyclus anynana butterflies

The eyespot patterns found on the wings of nymphalid butterflies are novel traits that originated first in hindwings and subsequently in forewings, suggesting that eyespot development might be dependent on Hox genes. Hindwings differ from forewings in the expression of Ultrabithorax (Ubx), but the function of this Hox gene in eyespot development as well as that of another Hox gene Antennapedia (Antp), expressed specifically in eyespots centers on both wings, are still unclear. We used CRISPR-Cas9 to target both genes in Bicyclus anynana butterflies. We show that Antp is essential for eyespot development on the forewings and for the differentiation of white centers and larger eyespots on hindwings, whereas Ubx is essential not only for the development of at least some hindwing eyespots but also for repressing the size of other eyespots. Additionally, Antp is essential for the development of silver scales in male wings. In summary, Antp and Ubx, in addition to their conserved roles in modifying serially homologous segments along the anterior-posterior axis of insects, have acquired a novel role in promoting the development of a new set of serial homologs, the eyespot patterns, in both forewings (Antp) and hindwings (Antp and Ubx) of B. anynana butterflies. We propose that the peculiar pattern of eyespot origins on hindwings first, followed by forewings, could be due to an initial co-option of Ubx into eyespot development followed by a later, partially redundant, co-option of Antp into the same network.

From Aedes to Zeugodacus: a review of dipteran body coloration studies regarding evolutionary developmental biology, pest control, and species discovery

Over the past two decades, evo-devo (evolution of development) studies have elucidated genetic mechanisms underlying novel dipteran body color patterns. Here we review the most recent developments, which show some departure from the model organism Drosophila melanogaster, leading the field into the investigation of more complex color patterns. We also discuss how the robust application of transgenic techniques has facilitated the study of many non-model pest species. Furthermore, we see that subtle pigmentation differences guide the discovery and description of new dipterans. Therefore, we argue that the existence of new field guides and the prevalence of pigmentation studies in non-model flies will enable scientists to adopt uninvestigated species into the lab, allowing them to study novel morphologies.

Variation in Pleiotropic Hub Gene Expression Is Associated with Interspecific Differences in Head Shape and Eye Size in Drosophila

Revealing the mechanisms underlying the breathtaking morphological diversity observed in nature is a major challenge in Biology. It has been established that recurrent mutations in hotspot genes cause the repeated evolution of morphological traits, such as body pigmentation or the gain and loss of structures. To date, however, it remains elusive whether hotspot genes contribute to natural variation in the size and shape of organs. As natural variation in head morphology is pervasive in Drosophila, we studied the molecular and developmental basis of differences in compound eye size and head shape in two closely related Drosophila species. We show differences in the progression of retinal differentiation between species and we applied comparative transcriptomics and chromatin accessibility data to identify the GATA transcription factor Pannier (Pnr) as central factor associated with these differences. Although the genetic manipulation of Pnr affected multiple aspects of dorsal head development, the effect of natural variation is restricted to a subset of the phenotypic space. We present data suggesting that this developmental constraint is caused by the coevolution of expression of pnr and its cofactor u-shaped (ush). We propose that natural variation in expression or function of highly connected developmental regulators with pleiotropic functions is a major driver for morphological evolution and we discuss implications on gene regulatory network evolution. In comparison to previous findings, our data strongly suggest that evolutionary hotspots are not the only contributors to the repeated evolution of eye size and head shape in Drosophila.

Many ways to make darker flies: Intra- and interspecific variation in Drosophila body pigmentation components

Body pigmentation is an evolutionarily diversified and ecologically relevant trait with substantial variation within and between species, and important roles in animal survival and reproduction. Insect pigmentation, in particular, provides some of the most compelling examples of adaptive evolution, including its ecological significance and genetic bases. Pigmentation includes multiple aspects of color and color pattern that may vary more or less independently, and can be under different selective pressures. We decompose Drosophila thorax and abdominal pigmentation, a valuable eco-evo-devo model, into distinct measurable traits related to color and color pattern. We investigate intra- and interspecific variation for those traits and assess its different sources. For each body part, we measured overall darkness, as well as four other pigmentation properties distinguishing between background color and color of the darker pattern elements that decorate each body part. By focusing on two standard D. melanogaster laboratory populations, we show that pigmentation components vary and covary in distinct manners depending on sex, genetic background, and temperature during development. Studying three natural populations of D. melanogaster along a latitudinal cline and five other Drosophila species, we then show that evolution of lighter or darker bodies can be achieved by changing distinct component traits. Our results paint a much more complex picture of body pigmentation variation than previous studies could uncover, including patterns of sexual dimorphism, thermal plasticity, and interspecific diversity. These findings underscore the value of detailed quantitative phenotyping and analysis of different sources of variation for a better understanding of phenotypic variation and diversification, and the ecological pressures and genetic mechanisms underlying them.

Redundant enhancers in the iab-5 domain cooperatively activate Abd-B in the A5 and A6 abdominal segments of Drosophila.. [DOI: 10.1101/2021.05.22.445252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The homeotic Abdominal-B (Abd-B) gene belongs to Bithorax complex and is regulated by four regulatory domains named iab-5, iab-6, iab-7 and iab-8, each of which is thought to be responsible for directing the expression of Abd-B in one of the abdominal segments from A5 to A8. It is assumed that male specific features of the adult cuticle in A5 is solely dependent on regulatory elements located in iab-5, while the regulatory elements in the iab-6 are both necessary and sufficient for the proper differentiation of the A6 cuticle. Unexpectedly, we found that this long held assumption is not correct. Instead, redundant tissue-specific enhancers located in the iab-5 domain are required for the proper activation of Abd-B not only in A5 but also in A6. Our study of deletions shows that the iab-5 initiator is essential for the functioning of the iab-5 enhancers in A5, as well as for the correct differentiation of A6. This requirement is circumvented by deletions that remove the initiator and most of the iab-5 regulatory domain sequences. While the remaining iab-5 enhancers are inactive in A5, they are activated in A6 and contribute to the differentiation of this segment. In this case, Abd-B stimulation by the iab-5 enhancers in A6 depends on the initiators in the iab-4 and iab-6 domains.

Summary Statement

In Drosophila, the segmental-specific expression of the homeotic gene Abdominal-B in the abdominal segments is regulated by autonomous regulatory domains. We demonstrated cooperation between these domains in activation of Abdominal-B.

Kinome-Wide RNAi Screen Uncovers Role of Ballchen in Maintenance of Gene Activation by Trithorax Group in Drosophila

Polycomb group (PcG) and trithorax group (trxG) proteins are evolutionary conserved factors that contribute to cell fate determination and maintenance of cellular identities during development of multicellular organisms. The PcG maintains heritable patterns of gene silencing while trxG acts as anti-silencing factors by conserving activation of cell type specific genes. Genetic and molecular analysis has revealed extensive details about how different PcG and trxG complexes antagonize each other to maintain cell fates, however, the cellular signaling components that contribute to the preservation of gene expression by PcG/trxG remain elusive. Here, we report an ex vivo kinome-wide RNAi screen in Drosophila aimed at identifying cell signaling genes that facilitate trxG in counteracting PcG mediated repression. From the list of trxG candidates, Ballchen (BALL), a histone kinase known to phosphorylate histone H2A at threonine 119 (H2AT119p), was characterized as a trxG regulator. The ball mutant exhibits strong genetic interactions with Polycomb (Pc) and trithorax (trx) mutants and loss of BALL affects expression of trxG target genes. BALL co-localizes with Trithorax on chromatin and depletion of BALL results in increased H2AK118 ubiquitination, a histone mark central to PcG mediated gene silencing. Moreover, BALL was found to substantially associate with known TRX binding sites across the genome. Genome wide distribution of BALL also overlaps with H3K4me3 and H3K27ac at actively transcribed genes. We propose that BALL mediated signaling positively contributes to the maintenance of gene activation by trxG in counteracting the repressive effect of PcG.

Transvection regulates the sex-biased expression of a fly X-linked gene

Sexual dimorphism in animals results from sex-biased gene expression patterns. These patterns are controlled by genetic sex determination hierarchies that establish the sex of an individual. Here we show that the male-biased wing expression pattern of the Drosophila biarmipes gene yellow, located on the X chromosome, is independent of the fly sex determination hierarchy. Instead, we find that a regulatory interaction between yellow alleles on homologous chromosomes (a process known as transvection) silences the activity of a yellow enhancer functioning in the wing. Therefore, this enhancer can be active in males (XY) but not in females (XX). This transvection-dependent enhancer silencing requires the yellow intron and the chromatin architecture protein Mod(mdg4). Our results suggest that transvection can contribute more generally to the sex-biased expression of X-linked genes.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.

[Phenotypic plasticity in insects]

Insects represent 85% of the animals. They have adapted to many environments and play a major role in ecosystems. Many insect species exhibit phenotypic plasticity. We here report on the mechanisms involved in phenotypic plasticity of different insects (aphids, migratory locust, map butterfly, honeybee) and also on the nutritional size plasticity in Drosophila and the plasticity of the wing eye-spots of the butterfly Bicyclus anynana. We also describe in more detail our work concerning the thermal plasticity of pigmentation in Drosophila. We have shown that the expression of the tan, yellow and Ddc genes, encoding enzymes of the melanin synthesis pathway, is modulated by temperature and that it is a consequence, at least in part, of the temperature-sensitive expression of the bab locus genes that repress them.

Bab2 Functions as an Ecdysone-Responsive Transcriptional Repressor during Drosophila Development

Drosophila development is governed by distinct ecdysone steroid pulses that initiate spatially and temporally defined gene expression programs. The translation of these signals into tissue-specific responses is crucial for metamorphosis, but the mechanisms that confer specificity to systemic ecdysone pulses are far from understood. Here, we identify Bric-à-brac 2 (Bab2) as an ecdysone-responsive transcriptional repressor that controls temporal gene expression during larval to pupal transition. Bab2 is necessary to terminate Salivary gland secretion (Sgs) gene expression, while premature Bab2 expression blocks Sgs genes and causes precocious salivary gland histolysis. The timely expression of bab2 is controlled by the ecdysone-responsive transcription factor Broad, and manipulation of EcR/USP/Broad signaling induces inappropriate Bab2 expression and termination of Sgs gene expression. Bab2 directly binds to Sgs loci in vitro and represses all Sgs genes in vivo. Our work characterizes Bab2 as a temporal regulator of somatic gene expression in response to systemic ecdysone signaling.

Transcriptome analysis reveals wingless regulates neural development and signaling genes in the region of wing pigmentation of a polka-dotted fruit fly

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Preexisting gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, to explore the genetic and molecular bases of the novel wing pigmentation pattern of a polka-dotted fruit fly (Drosophila guttifera), we performed de novo genome sequencing and transcriptome analyses. As a result, we comprehensively identified the genes associated with the pigmentation pattern. Furthermore, we revealed that 151 of these associated genes were positively or negatively regulated by wingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of a polka-dotted fruit fly might have emerged through multistep co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Copulatory function and development shape modular architecture of genitalia differently in males and females

Genitalia are multitasking structures whose development is mediated by numerous regulatory pathways. This multifactorial nature provides an avenue for multiple sources of selection. As a result, genitalia tend to evolve as modular systems comprising semi-independent subsets of structures, yet the processes that give rise to those patterns are still poorly understood. Here, we ask what are the relative roles of development and function in shaping modular patterns of genitalia within populations and across species of stink-bugs. We found that male genitalia are less integrated, more modular, and primarily shaped by functional demands. In contrast, females show higher integration, lower modularity, and a predominant role of developmental processes. Further, interactions among parts of each sex are more determinant to modularity than those between the sexes, and patterns of modularity are equivalent between and within species. Our results strongly indicate that genitalia have been subjected to sex-specific selection, although male and female genitalia are homologous and functionally associated. Moreover, modular patterns are seemingly constant in the evolutionary history of stink-bugs, suggesting a scenario of multivariate stabilizing selection within each sex. Our study demonstrates that interactions among genital parts of the same sex may be more fundamental to genital evolution than previously thought.

CRISPR/Cas9-mediated knockout of the eye pigmentation gene white leads to alterations in colour of head spots in the oriental fruit fly, Bactrocera dorsalis

The intensely studied white gene is widely used as a genetic marker in Drosophila melanogaster. Here, we cloned and characterized the white gene in an important pest of the fruit industry, Bactrocera dorsalis, to understand its functional role in pigmentation. We obtained BdWhite knockout strains, based on the wild-type strain, using the CRISPR/Cas9 genome editing system, and found that mutants lost pigmentation in the compound eye and their black head spots. We then examined differences in the expression levels of genes associated with melanin pigmentation between mutants and the wild-type strain using quantitative reverse transcription PCR. We found that transcription levels of the Bd-yellow1 were lower in the head of mutants than in the wild-type strain, and there were no significant differences in expression of the other six genes between mutants and the wild type. Since yellow is critical for melanin biosynthesis (Heinze et al., Scientific Reports. 2017;7:4582), the lower levels of expression of Bd-yellow1 in mutants led to reduced dark pigmentation in head spots. Our results provide the first evidence, to our knowledge, that white may play a functional role in cuticle pigmentation by affecting the expression of yellow.

The Hox gene Abdominal-B uses DoublesexF as a cofactor to promote neuroblast apoptosis in the Drosophila central nervous system

Highly conserved DM domain-containing transcription factors (Doublesex/MAB-3/DMRT1) are responsible for generating sexually dimorphic features. In the Drosophila central nervous system, a set of Doublesex (Dsx)-expressing neuroblasts undergo apoptosis in females whereas their male counterparts proliferate and give rise to serotonergic neurons crucial for adult mating behaviour. Our study demonstrates that the female-specific isoform of Dsx collaborates with Hox gene Abdominal-B (Abd-B) to bring about this apoptosis. Biochemical results suggest that proteins AbdB and Dsx interact through their highly conserved homeodomain and DM domain, respectively. This interaction is translated into a cooperative binding of the two proteins on the apoptotic enhancer in the case of females but not in the case of males, resulting in female-specific activation of apoptotic genes. The capacity of AbdB to use the sex-specific isoform of Dsx as a cofactor underlines the possibility that these two classes of protein are capable of cooperating in selection and regulation of target genes in a tissue- and sex-specific manner. We propose that this interaction could be a common theme in generating sexual dimorphism in different tissues across different species.

Modular cis-regulatory logic of yellow gene expression in silkmoth larvae

Colour patterns in butterflies and moths are crucial traits for adaptation. Previous investigations have highlighted genes responsible for pigmentation (ie yellow and ebony). However, the mechanisms by which these genes are regulated in lepidopteran insects remain poorly understood. To elucidate this, molecular studies involving dipterans have largely analysed the cis-regulatory regions of pigmentation genes and have revealed cis-regulatory modularity. Here, we used well-developed transgenic techniques in Bombyx mori and demonstrated that cis-regulatory modularity controls tissue-specific expression of the yellow gene. We first identified which body parts are regulated by the yellow gene via black pigmentation. We then isolated three discrete regulatory elements driving tissue-specific gene expression in three regions of B. mori larvae. Finally, we found that there is no apparent sequence conservation of cis-regulatory regions between B. mori and Drosophila melanogaster, and no expression driven by the regulatory regions of one species when introduced into the other species. Therefore, the trans-regulatory landscapes of the yellow gene differ significantly between the two taxa. The results of this study confirm that lepidopteran species use cis-regulatory modules to control gene expression related to pigmentation, and represent a powerful cadre of transgenic tools for studying evolutionary developmental mechanisms.

Changes throughout a Genetic Network Mask the Contribution of Hox Gene Evolution

Hox genes pattern the anterior-posterior axis of animals and are posited to drive animal body plan evolution, yet their precise role in evolution has been difficult to determine. Here, we identified evolutionary modifications in the Hox gene Abd-B that dramatically altered its expression along the body plan of Drosophila santomea. Abd-B is required for pigmentation in Drosophila yakuba, the sister species of D. santomea, and changes to Abd-B expression would be predicted to make large contributions to the loss of body pigmentation in D. santomea. However, manipulating Abd-B expression in current-day D. santomea does not affect pigmentation. We attribute this epistatic interaction to four other genes within the D. santomea pigmentation network, three of which have evolved expression patterns that do not respond to Abd-B. Our results demonstrate how body plans may evolve through small evolutionary steps distributed throughout Hox-regulated networks. Polygenicity and epistasis may hinder efforts to identify genes and mechanisms underlying macroevolutionary traits.

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.

A homeotic shift late in development drives mimetic color variation in a bumble bee

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.