The 19 genomes of Drosophila: a BAC library resource for genus-wide and genome-scale comparative evolutionary research

The role of fruitless in specifying courtship behaviors across divergent Drosophila species

Sex-specific behaviors are critical for reproduction and species survival. The sex-specifically spliced transcription factor fruitless (fru) helps establish male courtship behaviors in invertebrates. Forcing male-specific fru (fruM) splicing in Drosophila melanogaster females produces male-typical behaviors while disrupting female-specific behaviors. However, whether fru's joint role in specifying male and inhibiting female behaviors is conserved across species is unknown. We used CRISPR-Cas9 to force FruM expression in female Drosophila virilis, a species in which males and females produce sex-specific songs. In contrast to D. melanogaster, in which one fruM allele is sufficient to generate male behaviors in females, two alleles are needed in D. virilis females. D. virilis females expressing FruM maintain the ability to sing female-typical song as well as lay eggs, whereas D. melanogaster FruM females cannot lay eggs. These results reveal potential differences in fru function between divergent species and underscore the importance of studying diverse behaviors and species for understanding the genetic basis of sex differences.

Genomic resources of broomcorn millet: demonstration and application of a high-throughput BAC mapping pipeline

Background

With high-efficient water-use and drought tolerance, broomcorn millet has emerged as a candidate for food security. To promote its research process for molecular breeding and functional research, a comprehensive genome resource is of great importance.

Results

Herein, we constructed a BAC library for broomcorn millet, generated BAC end sequences based on the clone-array pooled shotgun sequencing strategy and Illumina sequencing technology, and integrated BAC clones into genome by a novel pipeline for BAC end profiling. The BAC library consisted of 76,023 clones with an average insert length of 123.48 Kb, covering about 9.9-fold of the 850 Mb genome. Of 9216 clones tested using our pipeline, 8262 clones were mapped on the broomcorn millet cultivar longmi4 genome. These mapped clones covered 308 of the 829 gaps left by the genome. To our knowledge, this is the only BAC resource for broomcorn millet.

Conclusions

We constructed a high-quality BAC libraray for broomcorn millet and designed a novel pipeline for BAC end profiling. BAC clones can be browsed and obtained from our website (http://eightstarsbio.com/gresource/JBrowse-1.16.5/index.html). The high-quality BAC clones mapped on genome in this study will provide a powerful genomic resource for genome gap filling, complex segment sequencing, FISH, functional research and genetic engineering of broomcorn millet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-01003-z.

Chemosensation and Evolution of Drosophila Host Plant Selection

The ability to respond to chemosensory cues is critical for survival of most organisms. Among insects, Drosophila melanogaster has the best characterized olfactory system, and the availability of genome sequences of 30 Drosophila species provides an ideal scenario for studies on evolution of chemosensation. Gene duplications of chemoreceptor genes allow for functional diversification of the rapidly evolving chemoreceptor repertoire. Although some species of the genus Drosophila are generalists for host plant selection, rapid evolution of olfactory receptors, gustatory receptors, odorant-binding proteins, and cytochrome P450s has enabled diverse host specializations of different members of the genus. Here, I review diversification of the chemoreceptor repertoire among members of the genus Drosophila along with co-evolution of detoxification mechanisms that may have enabled occupation of diverse host plant ecological niches.

Evolution of Reproductive Behavior

Behaviors associated with reproduction are major contributors to the evolutionary success of organisms and are subject to many evolutionary forces, including natural and sexual selection, and sexual conflict. Successful reproduction involves a range of behaviors, from finding an appropriate mate, courting, and copulation, to the successful production and (in oviparous animals) deposition of eggs following mating. As a consequence, behaviors and genes associated with reproduction are often under strong selection and evolve rapidly. Courtship rituals in flies follow a multimodal pattern, mediated through visual, chemical, tactile, and auditory signals. Premating behaviors allow males and females to assess the species identity, reproductive state, and condition of their partners. Conflicts between the "interests" of individual males, and/or between the reproductive strategies of males and females, often drive the evolution of reproductive behaviors. For example, seminal proteins transmitted by males often show evidence of rapid evolution, mediated by positive selection. Postmating behaviors, including the selection of oviposition sites, are highly variable and Drosophila species span the spectrum from generalists to obligate specialists. Chemical recognition features prominently in adaptation to host plants for feeding and oviposition. Selection acting on variation in pre-, peri-, and postmating behaviors can lead to reproductive isolation and incipient speciation. Response to selection at the genetic level can include the expansion of gene families, such as those for detecting pheromonal cues for mating, or changes in the expression of genes leading to visual cues such as wing spots that are assessed during mating. Here, we consider the evolution of reproductive behavior in Drosophila at two distinct, yet complementary, scales. Some studies take a microevolutionary approach, identifying genes and networks involved in reproduction, and then dissecting the genetics underlying complex behaviors in D. melanogaster Other studies take a macroevolutionary approach, comparing reproductive behaviors across the genus Drosophila and how these might correlate with environmental cues. A full synthesis of this field will require unification across these levels.

Optimal Scaling of Critical Size for Metamorphosis in the Genus Drosophila

Juveniles must reach a critical body size to become a mature adult. Molecular determinants of critical size have been studied, but the evolutionary importance of critical size is still unclear. Here, using nine fly species, we show that interspecific variation in organism size can be explained solely by species-specific critical size. The observed variation in critical size quantitatively agrees with the interspecific scaling relationship predicted by the life history model, which hypothesizes that critical size mediates an energy allocation switch between juvenile and adult tissues. The mechanism underlying critical size scaling is explained by an inverse relationship between growth duration and growth rate, which cancels out their contributions to the final size. Finally, we show that evolutionary changes in growth duration can be traced back to the scaling of ecdysteroid hormone dynamics. We conclude that critical size adaptively optimizes energy allocation, and has a central role in organism size determination.

The Drosophila Dot Chromosome: Where Genes Flourish Amidst Repeats

The F element of the Drosophila karyotype (the fourth chromosome in Drosophila melanogaster) is often referred to as the "dot chromosome" because of its appearance in a metaphase chromosome spread. This chromosome is distinct from other Drosophila autosomes in possessing both a high level of repetitious sequences (in particular, remnants of transposable elements) and a gene density similar to that found in the other chromosome arms, ∼80 genes distributed throughout its 1.3-Mb "long arm." The dot chromosome is notorious for its lack of recombination and is often neglected as a consequence. This and other features suggest that the F element is packaged as heterochromatin throughout. F element genes have distinct characteristics (e.g, low codon bias, and larger size due both to larger introns and an increased number of exons), but exhibit expression levels comparable to genes found in euchromatin. Mapping experiments show the presence of appropriate chromatin modifications for the formation of DNaseI hypersensitive sites and transcript initiation at the 5' ends of active genes, but, in most cases, high levels of heterochromatin proteins are observed over the body of these genes. These various features raise many interesting questions about the relationships of chromatin structures with gene and chromosome function. The apparent evolution of the F element as an autosome from an ancestral sex chromosome also raises intriguing questions. The findings argue that the F element is a unique chromosome that occupies its own space in the nucleus. Further study of the F element should provide new insights into chromosome structure and function.

Drosophila melanogaster as a Model for Diabetes Type 2 Progression

Drosophila melanogaster has been used as a very versatile and potent model in the past few years for studies in metabolism and metabolic disorders, including diabetes types 1 and 2. Drosophila insulin signaling, despite having seven insulin-like peptides with partially redundant functions, is very similar to the human insulin pathway and has served to study many different aspects of diabetes and the diabetic state. Yet, very few studies have addressed the chronic nature of diabetes, key for understanding the full-blown disease, which most studies normally explore. One of the advantages of having Drosophila mutant viable combinations at different levels of the insulin pathway, with significantly reduced insulin pathway signaling, is that the abnormal metabolic state can be studied from the onset of the life cycle and followed throughout. In this review, we look at the chronic nature of impaired insulin signaling. We also compare these results to the results gleaned from vertebrate model studies.

Genomic changes following the reversal of a Y chromosome to an autosome in Drosophila pseudoobscura

Robertsonian translocations resulting in fusions between sex chromosomes and autosomes shape karyotype evolution by creating new sex chromosomes from autosomes. These translocations can also reverse sex chromosomes back into autosomes, which is especially intriguing given the dramatic differences between autosomes and sex chromosomes. To study the genomic events following a Y chromosome reversal, we investigated an autosome‐Y translocation in Drosophila pseudoobscura. The ancestral Y chromosome fused to a small autosome (the dot chromosome) approximately 10–15 Mya. We used single molecule real‐time sequencing reads to assemble the D. pseudoobscura dot chromosome, including this Y‐to‐dot translocation. We find that the intervening sequence between the ancestral Y and the rest of the dot chromosome is only ∼78 Kb and is not repeat‐dense, suggesting that the centromere now falls outside, rather than between, the fused chromosomes. The Y‐to‐dot region is 100 times smaller than the D. melanogaster Y chromosome, owing to changes in repeat landscape. However, we do not find a consistent reduction in intron sizes across the Y‐to‐dot region. Instead, deletions in intergenic regions and possibly a small ancestral Y chromosome size may explain the compact size of the Y‐to‐dot translocation.

Genetic and Transgenic Reagents for Drosophila simulans, D. mauritiana, D. yakuba, D. santomea, and D. virilis

Species of the Drosophila melanogaster species subgroup, including the species D. simulans, D. mauritiana, D. yakuba, and D. santomea, have long served as model systems for studying evolution. However, studies in these species have been limited by a paucity of genetic and transgenic reagents. Here, we describe a collection of transgenic and genetic strains generated to facilitate genetic studies within and between these species. We have generated many strains of each species containing mapped piggyBac transposons including an enhanced yellow fluorescent protein (EYFP) gene expressed in the eyes and a ϕC31 attP site-specific integration site. We have tested a subset of these lines for integration efficiency and reporter gene expression levels. We have also generated a smaller collection of other lines expressing other genetically encoded fluorescent molecules in the eyes and a number of other transgenic reagents that will be useful for functional studies in these species. In addition, we have mapped the insertion locations of 58 transposable elements in D. virilis that will be useful for genetic mapping studies.

Genome engineering: Drosophila melanogaster and beyond

A central challenge in investigating biological phenomena is the development of techniques to modify genomic DNA with nucleotide precision that can be transmitted through the germ line. Recent years have brought a boon in these technologies, now collectively known as genome engineering. Defined genomic manipulations at the nucleotide level enable a variety of reverse engineering paradigms, providing new opportunities to interrogate diverse biological functions. These genetic modifications include controlled removal, insertion, and substitution of genetic fragments, both small and large. Small fragments up to a few kilobases (e.g., single nucleotide mutations, small deletions, or gene tagging at single or multiple gene loci) to large fragments up to megabase resolution can be manipulated at single loci to create deletions, duplications, inversions, or translocations of substantial sections of whole chromosome arms. A specialized substitution of chromosomal portions that presumably are functionally orthologous between different organisms through syntenic replacement, can provide proof of evolutionary conservation between regulatory sequences. Large transgenes containing endogenous or synthetic DNA can be integrated at defined genomic locations, permitting an alternative proof of evolutionary conservation, and sophisticated transgenes can be used to interrogate biological phenomena. Precision engineering can additionally be used to manipulate the genomes of organelles (e.g., mitochondria). Novel genome engineering paradigms are often accelerated in existing, easily genetically tractable model organisms, primarily because these paradigms can be integrated in a rigorous, existing technology foundation. The Drosophila melanogaster fly model is ideal for these types of studies. Due to its small genome size, having just four chromosomes, the vast amount of cutting-edge genetic technologies, and its short life-cycle and inexpensive maintenance requirements, the fly is exceptionally amenable to complex genetic analysis using advanced genome engineering. Thus, highly sophisticated methods developed in the fly model can be used in nearly any sequenced organism. Here, we summarize different ways to perform precise inheritable genome engineering using integrases, recombinases, and DNA nucleases in the D. melanogaster. For further resources related to this article, please visit the WIREs website.

Immune-directed support of rich microbial communities in the gut has ancient roots

The animal gut serves as a primary location for the complex host-microbe interplay that is essential for homeostasis and may also reflect the types of ancient selective pressures that spawned the emergence of immunity in metazoans. In this review, we present a phylogenetic survey of gut host-microbe interactions and suggest that host defense systems arose not only to protect tissue directly from pathogenic attack but also to actively support growth of specific communities of mutualists. This functional dichotomy resulted in the evolution of immune systems much more tuned for harmonious existence with microbes than previously thought, existing as dynamic but primarily cooperative entities in the present day. We further present the protochordate Ciona intestinalis as a promising model for studying gut host-bacterial dialogue. The taxonomic position, gut physiology and experimental tractability of Ciona offer unique advantages in dissecting host-microbe interplay and can complement studies in other model systems.

Ecological genomics of host shifts in Drosophila mojavensis

Advances in next-generation sequencing technologies have liberated our dependency on model laboratory species for answering genomic and transcriptomic level questions. These new techniques have dramatically expanded our breadth of study organisms and have allowed the analysis of species from diverse ecological environments. One such species is the cactophilic Drosophila mojavensis that inhabits the deserts of western North America. These insects feed and develop in the necrotic cacti, feeding largely on the microflora of the necrotic plant tissues. Drosophila mojavensis is composed of four geographically and ecologically separated populations. Each population (Baja California peninsula, mainland Sonoran Desert, Mojave Desert and Santa Catalina Island) utilizes the necrotic tissues of distinct cactus species. The differences in the nutritional and chemical composition of the necroses include a set of toxic compounds to which resident population must adapt. These ecological differences have facilitated many of the life history, behavior, physiological and genetic differences between the cactus host populations. Genomic resources have allowed investigators to examine the genomic and transcriptional level changes associated with the local adaptation of the four D. mojavensis populations, thereby providing further understanding of the genetic mechanism of adaptation and its role in the divergence of ecologically distinct populations.

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

A comparison of the signalling properties of two tyramine receptors from Drosophila

In invertebrates, the phenolamines, tyramine and octopamine, mediate many functional roles usually associated with the catecholamines, noradrenaline and adrenaline, in vertebrates. The α- and β-adrenergic classes of insect octopamine receptor are better activated by octopamine than tyramine. Similarly, the Tyramine 1 subgroup of receptors (or Octopamine/Tyramine receptors) are better activated by tyramine than octopamine. However, recently, a new Tyramine 2 subgroup of receptors was identified, which appears to be activated highly preferentially by tyramine. We examined immunocytochemically the ability of CG7431, the founding member of this subgroup from Drosophila melanogaster, to be internalized in transfected Chinese hamster ovary (CHO) cells by different agonists. It was only internalized after activation by tyramine. Conversely, the structurally related receptor, CG16766, was internalized by a number of biogenic amines, including octopamine, dopamine, noradrenaline, adrenaline, which also were able to elevate cyclic AMP levels. Studies with synthetic agonists and antagonists confirm that CG16766 has a different pharmacological profile to that of CG7431. Species orthologues of CG16766 were only found in Drosophila species, whereas orthologues of CG7431 could be identified in the genomes of a number of insect species. We propose that CG16766 represents a new group of tyramine receptors, which we have designated the Tyramine 3 receptors.

Dissecting sources of quantitative gene expression pattern divergence between Drosophila species

Gene expression patterns can diverge between species due to changes in a gene's regulatory DNA or changes in the proteins, e.g., transcription factors (TFs), that regulate the gene. We developed a modeling framework to uncover the sources of expression differences in blastoderm embryos of three Drosophila species, focusing on the regulatory circuit controlling expression of the hunchback (hb) posterior stripe. Using this framework and cellular-resolution expression measurements of hb and its regulating TFs, we found that changes in the expression patterns of hb's TFs account for much of the expression divergence. We confirmed our predictions using transgenic D. melanogaster lines, which demonstrate that this set of orthologous cis-regulatory elements (CREs) direct similar, but not identical, expression patterns. We related expression pattern differences to sequence changes in the CRE using a calculation of the CRE's TF binding site content. By applying this calculation in both the transgenic and endogenous contexts, we found that changes in binding site content affect sensitivity to regulating TFs and that compensatory evolution may occur in circuit components other than the CRE.

A pair of new BAC and BIBAC vectors that facilitate BAC/BIBAC library construction and intact large genomic DNA insert exchange

BACKGROUND

Large-insert BAC and BIBAC libraries are important tools for structural and functional genomics studies of eukaryotic genomes. To facilitate the construction of BAC and BIBAC libraries and the transfer of complete large BAC inserts into BIBAC vectors, which is desired in positional cloning, we developed a pair of new BAC and BIBAC vectors.

RESULTS

The new BAC vector pIndigoBAC536-S and the new BIBAC vector BIBAC-S have the following features: 1) both contain two 18-bp non-palindromic I-SceI sites in an inverted orientation at positions that flank an identical DNA fragment containing the lacZ selection marker and the cloning site. Large DNA inserts can be excised from the vectors as single fragments by cutting with I-SceI, allowing the inserts to be easily sized. More importantly, because the two vectors contain different antibiotic resistance genes for transformant selection and produce the same non-complementary 3' protruding ATAA ends by I-SceI that suppress self- and inter-ligations, the exchange of intact large genomic DNA inserts between the BAC and BIBAC vectors is straightforward; 2) both were constructed as high-copy composite vectors. Reliable linearized and dephosphorylated original low-copy pIndigoBAC536-S and BIBAC-S vectors that are ready for library construction can be prepared from the high-copy composite vectors pHZAUBAC1 and pHZAUBIBAC1, respectively, without the need for additional preparation steps or special reagents, thus simplifying the construction of BAC and BIBAC libraries. BIBAC clones constructed with the new BIBAC-S vector are stable in both E. coli and Agrobacterium. The vectors can be accessed through our website http://GResource.hzau.edu.cn.

CONCLUSIONS

The two new vectors and their respective high-copy composite vectors can largely facilitate the construction and characterization of BAC and BIBAC libraries. The transfer of complete large genomic DNA inserts from one vector to the other is made straightforward.

Genomic sequence around butterfly wing development genes: annotation and comparative analysis

BACKGROUND

Analysis of genomic sequence allows characterization of genome content and organization, and access beyond gene-coding regions for identification of functional elements. BAC libraries, where relatively large genomic regions are made readily available, are especially useful for species without a fully sequenced genome and can increase genomic coverage of phylogenetic and biological diversity. For example, no butterfly genome is yet available despite the unique genetic and biological properties of this group, such as diversified wing color patterns. The evolution and development of these patterns is being studied in a few target species, including Bicyclus anynana, where a whole-genome BAC library allows targeted access to large genomic regions.

METHODOLOGY/PRINCIPAL FINDINGS

We characterize ∼1.3 Mb of genomic sequence around 11 selected genes expressed in B. anynana developing wings. Extensive manual curation of in silico predictions, also making use of a large dataset of expressed genes for this species, identified repetitive elements and protein coding sequence, and highlighted an expansion of Alcohol dehydrogenase genes. Comparative analysis with orthologous regions of the lepidopteran reference genome allowed assessment of conservation of fine-scale synteny (with detection of new inversions and translocations) and of DNA sequence (with detection of high levels of conservation of non-coding regions around some, but not all, developmental genes).

CONCLUSIONS

The general properties and organization of the available B. anynana genomic sequence are similar to the lepidopteran reference, despite the more than 140 MY divergence. Our results lay the groundwork for further studies of new interesting findings in relation to both coding and non-coding sequence: 1) the Alcohol dehydrogenase expansion with higher similarity between the five tandemly-repeated B. anynana paralogs than with the corresponding B. mori orthologs, and 2) the high conservation of non-coding sequence around the genes wingless and Ecdysone receptor, both involved in multiple developmental processes including wing pattern formation.