Phylogenetic relationships of Drosophila melanogaster species group deduced from spacer regions of histone gene H2A-H2B

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

Genome-Wide Characterization and Analysis of Expression of the Histone Gene Family in Razor Clam, Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta), a bivalve species widely distributed in estuaries and mudflats, is often exposed to extreme environmental and microbial stresses. Histones are fundamental components of chromatin and play an important role in innate immunity, as demonstrated by its antimicrobial activities in clams. However, little attention has been paid to histones in bivalves. To fill this gap, we investigated the genomic distribution, structural characteristics, conserved motifs, and phylogenetic relationships of histones in S. constricta. A total of 114 histone genes were detected in the S. constricta genome, which were divided into 25 types in phylogenetic analysis. Among them, partial histones exhibited a tissue-dependent expression pattern, indicating that they may be involved in sustaining the homeostasis of organs/tissues in adult S. constricta. Furthermore, mRNA expression of certain histones changed significantly in S. constricta when infected with Vibrio parahaemolyticus, suggesting that histones play a role in the immune defense of S. constricta. All together, this study on histone genes in S. constricta not only greatly expands our knowledge of histone function in the clam, but also histone evolution in molluscs.

The chromosomes of Drosophila suzukii (Diptera: Drosophilidae): detailed photographic polytene chromosomal maps and in situ hybridization data

The spotted wing drosophila, D. suzukii, is a serious agricultural pest attacking a variety of soft fruits and vegetables. Although originating from East Asia it has recently invaded America and Europe raising major concern about its expansion potential and the consequent economic losses. Since cytogenetic information on the species is scarce, we report here the mitotic karyotype and detailed photographic maps of the salivary gland polytene chromosomes of D. suzukii. The mitotic metaphase complement contains three pairs of autosomes, one of which is dot-like, and one pair of heteromorphic (XX/XY) sex chromosomes. The salivary gland polytene complement consists of five long polytene arms, representing the two metacentric autosomes and the acrocentric X chromosome, and one very short polytene element, which corresponds to the dot-like autosome. Banding pattern as well as the most characteristic features and prominent landmarks of each polytene chromosome arm are presented and discussed. Furthermore, twelve gene markers have been mapped on the polytene chromosomes of D. suzukii by in situ hybridization. Their distribution pattern was found quite similar to that of D. melanogaster revealing conservation of synteny although the relative position within each chromosome arm for most of the genes differed significantly between D. suzukii and D. melanogaster. The chromosome information presented here is suitable for comparative cytogenetic studies and phylogenetic exploration, while it could also facilitate the assembly of the genome sequence and support the development of genetic tools for species-specific and environment-friendly biological control applications such as the sterile insect technique.

Morphological Markers for Cryopreservation in the Embryonic Development of Drosophila suzukii (Diptera: Drosophilidae)

Drosophila suzukii (Matsumura) (Diptera: Drosophilidae) is an invasive pest recently reported in Europe whose spread into new areas has caused severe economic damage to many agricultural crops. There are serious concerns about the currently available chemical insecticides because of their low efficacy in controlling the species and their environmental impact; so, several studies have focused on environmentally safe strategies. The sterile insect technique (SIT), which requires colony maintenance in laboratory and production of large numbers of live animals, can be utilized in pest management programs and could be integrated with other control strategies if the potential risks associated with the rearing and maintenance of the insect line under laboratory conditions are given sufficient attention. In this regard, the ability to cryobiologically preserve such stocks would be of substantial value. Important prerequisites for long-term cryopreservation are determination of the embryonic stages, identification of specific embryonic stages, and knowledge of development time. This paper describes the main visible markers for the different stages of embryonic development and determines the timing of development at 25°C. D. suzukii embryogenesis lasts 23-25 h at 25°C and can be divided into 17 stages defined by specific morphological markers. The point at which 50% of embryos are at Stage 14 and 50% are at Stage 15, the most tolerant stages for cryopreservation treatment, as ascertained for Drosophila melanogaster Meigen in prior studies, is reached in 14-15 h. The efficiency of this procedure might be impaired by the retention of eggs in the oviducts, making it impossible to determine the stage of embryonic development for ∼25% of laid eggs.

Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments

Type IIB restriction endonucleases are site-specific endonucleases that cut both strands of double-stranded DNA upstream and downstream of their recognition sequences. These restriction enzymes have recognition sequences that are generally interrupted and range from 5 to 7 bases long. They produce DNA fragments which are uniformly small, ranging from 21 to 33 base pairs in length (without cohesive ends). The fragments are generated from throughout the entire length of a genomic DNA providing an excellent fractional representation of the genome. In this study we simulated restriction enzyme digestions on 21 sequenced genomes of various Drosophila species using the predicted targets of 16 Type IIB restriction enzymes to effectively produce a large and arbitrary selection of loci from these genomes. The fragments were then used to compare organisms and to calculate the distance between genomes in pair-wise combination by counting the number of shared fragments between the two genomes. Phylogenetic trees were then generated for each enzyme using this distance measure and the consensus was calculated. The consensus tree obtained agrees well with the currently accepted tree for the Drosophila species. We conclude that multi-locus sub-genomic representation combined with next generation sequencing, especially for individuals and species without previous genome characterization, can accelerate studies of comparative genomics and the building of accurate phylogenetic trees.

DNA Barcoding and Molecular Phylogeny of Drosophila lini and Its Sibling Species

Drosophila lini and its two sibling species, D. ohnishii and D. ogumai, are hardly distinguishable from one another in morphology. These species are more or less reproductively isolated. The mitochondrial ND2 and COI-COII and the nuclear ITS1-ITS2 regions were sequenced to seek for the possibility of DNA barcoding and to reconstruct the phylogeny of them. The character-based approach for DNA barcoding detected some diagnostic nucleotides only for monophyletic D. ogumai, but no informative sites for the other two very closely species, D. lini and D. ohnishii, of which strains intermingled in the molecular phylogenetic trees. Thus, this study provides another case of limited applicability of DNA barcoding in species delineation, as in other cases of related Drosophila species. The molecular phylogenetic tree inferred from the concatenated sequences strongly supported the monophyly of the cluster of the three species, that is, the lini clade. We propose some hypotheses of evolutionary events in this clade.

Increasing the data size to accurately reconstruct the phylogenetic relationships between nine subgroups of the Drosophila melanogaster species group (Drosophilidae, Diptera)

Previous phylogenetic analyses of the melanogaster species group have led to conflicting hypotheses concerning their relationship; therefore the addition of new sequence data is necessary to discover the phylogeny of this species group. Here we present new data derived from 17 genes and representing 48 species to reconstruct the phylogeny of the melanogaster group. A variety of statistical tests, as well as maximum likelihood mapping analysis, were performed to estimate data quality, suggesting that all genes had a high degree of contribution to resolve the phylogeny. Individual locus was analyzed using maximum likelihood (ML), and the concatenated dataset (12,988 bp) were analyzed using partitioned maximum likelihood (ML) and Bayesian analyses. Separated analysis produced various phylogenetic relationships, however, phylogenetic topologies from ML and Bayesian analysis based on concatenated dataset, at the subgroup level, were completely identical to each other with high levels of support. Our results recovered three major clades: the ananassae subgroup, followed by the montium subgroup, the melanogaster subgroup and the oriental subgroups form the third monophyletic clade, in which melanogaster (takahashii, suzukii) forms one subclade and ficusphila [eugracilis (elegans, rhopaloa)] forms another. However, more data are necessary to determine the phylogenetic position of Drosophila lucipennis which proved difficult to place.

Histone and ribosomal RNA repetitive gene clusters of the boll weevil are linked in a tandem array

Histones are the major protein component of chromatin structure. The histone family is made up of a quintet of proteins, four core histones (H2A, H2B, H3 & H4) and the linker histones (H1). Spacers are found between the coding regions. Among insects this quintet of genes is usually clustered and the clusters are tandemly repeated. Ribosomal DNA contains a cluster of the rRNA sequences 18S, 5.8S and 28S. The rRNA genes are separated by the spacers ITS1, ITS2 and IGS. This cluster is also tandemly repeated. We found that the ribosomal RNA repeat unit of at least two species of Anthonomine weevils, Anthonomus grandis and Anthonomus texanus (Coleoptera: Curculionidae), is interspersed with a block containing the histone gene quintet. The histone genes are situated between the rRNA 18S and 28S genes in what is known as the intergenic spacer region (IGS). The complete reiterated Anthonomus grandis histone-ribosomal sequence is 16,248 bp.

Research on the karyotype and evolution of Drosophila melanogaster species group

Mitotic metaphase chromosomes of 34 species of Drosophila melanogaster species group were examined. Certain new karyotypes were described for the first time, and their evolutionary and interspecific genetic relationships among 8 subgroups of D. melanogaster species group were analyzed systematically. The results were as follows. The basic karyotype of elegans subgroup was type A. The karyotypes of eugracilis subgroup, melanogaster subgroup, and ficusphila subgroup were all type C. The karyotypes of takahashii subgroup and suzukii subgroup were both type C and type D. The montium subgroup had six kinds of karyotypes: types B, C, C', D, D', and E. The ananassae subgroup had three kinds of karyotypes: types F, G, and H. Thus, the melanogaster species group was classified into five pedigrees based on the diversity of these karyotypes: 1) elegans; 2) eugracilis-melanogaster-ficusphila; 3) takkahashii-suzukii; 4) montium; 5) ananassae. The above-mentioned results in karyotypic evolution were consistent with those of DNA sequence analysis reported by Yang except for the elegans subgroup and this subgroup was considered as the ancestral subgroup. Karyotype analysis of the same drosophila from different isofemale lines indicated that the same Drosophila from different places showed karyotypic variation which might be due to different geographical environment and evolutionary degree or interaction between the two factors.

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.

Genomic organization, nucleotide sequence analysis of the core histone genes cluster in Chlamys farreri and molecular evolution assessment of the H2A and H2B

This work represents the nucleotide sequence of the core histone gene cluster from scallop Chlamys farreri. The tandemly repeated unit of 5671 bp containing a copy of the four core histone genes H4, H2B, H2A and H3 was amplified and identified by the techniques of homology cloning and genomic DNA walking. All the histone genes in the cluster had the structures in their 3' flanking region which related to the evolution of histone gene expression patterns throughout the cell cycle, including two different termination signals, the hairpin structure and at least one AATAAA polyadenylation signal. In their 5' region, the transcription initiation sites with a conserved sequence of 5'-PyATTCPu-3' known as the CAP site were present in all genes except to H2B, generally 37-45 bp upstream of the start code. Canonical TATA and CAAT boxes were identified only in certain histone genes. In the case of the promoters of H2B and H2A genes, there was a 5'-GATCC-3' element, which had been found to be essential to start transcription at the appropriate site. After this element, in the promoter of H2B, there was another sequence, 5'-GGATCGAAACGTTC-3', which was similar to the consensus sequence of 5'-GGAATAAACGTATTC-3' corresponding to the H2B-specific promoter element. The presence of enhancer sequences (5'-TGATATATG-3') was identified from the H4 and H3 genes, matching perfectly with the consensus sequence defined for histone genes. There were several slightly more complex repetitive DNA in the intergene regions. The presence of the series of conserved sequences and reiterated sequences was consistent with the view that mollusc histone gene cluster arose by duplicating of an ancestral precursor histone gene, the birth-and-death evolution model with strong purifying selection enabled the histone cluster less variation and more conserved function. Meanwhile, the H2A and the H2B were demonstrated to be potential good marks for phylogenetic analysis. All the results will be contributed to the characterization of repeating histone gene families in molluscs.

Phylogenetic incongruence in the Drosophila melanogaster species group

Drosophila melanogaster and its close relatives are used extensively in comparative biology. Despite the importance of phylogenetic information for such studies, relationships between some melanogaster species group members are unclear due to conflicting phylogenetic signals at different loci. In this study, we use twelve nuclear loci (eleven coding and one non-coding) to assess the degree of phylogenetic incongruence in this model system. We focus on two nodes: (1) the node joining the Drosophila erecta-Drosophila orena, Drosophila melanogaster-Drosophila simulans, and Drosophila yakuba-Drosophila teissieri lineages, and (2) the node joining the lineages leading to the melanogaster, takahashii, and eugracilis subgroups. We find limited evidence for incongruence at the first node; our data, as well as those of several previous studies, strongly support monophyly of a clade consisting of D. erecta-D. orena and D. yakuba-D. teissieri. By contrast, using likelihood based tests of congruence, we find robust evidence for topological incongruence at the second node. Different loci support different relationships among the melanogaster, takahashii, and eugracilis subgroups, and the observed incongruence is not easily attributable to homoplasy, non-equilibrium base composition, or positive selection on a subset of loci. We argue that lineage sorting in the common ancestor of these three subgroups is the most plausible explanation for our observations. Such lineage sorting may lead to biased estimation of tree topology and evolutionary rates, and may confound inferences of positive selection.

Positive selection of Iris, a retroviral envelope-derived host gene in Drosophila melanogaster

Eukaryotic genomes can usurp enzymatic functions encoded by mobile elements for their own use. A particularly interesting kind of acquisition involves the domestication of retroviral envelope genes, which confer infectious membrane-fusion ability to retroviruses. So far, these examples have been limited to vertebrate genomes, including primates where the domesticated envelope is under purifying selection to assist placental function. Here, we show that in Drosophila genomes, a previously unannotated gene (CG4715, renamed Iris) was domesticated from a novel, active Kanga lineage of insect retroviruses at least 25 million years ago, and has since been maintained as a host gene that is expressed in all adult tissues. Iris and the envelope genes from Kanga retroviruses are homologous to those found in insect baculoviruses and gypsy and roo insect retroviruses. Two separate envelope domestications from the Kanga and roo retroviruses have taken place, in fruit fly and mosquito genomes, respectively. Whereas retroviral envelopes are proteolytically cleaved into the ligand-interaction and membrane-fusion domains, Iris appears to lack this cleavage site. In the takahashii/suzukii species groups of Drosophila, we find that Iris has tandemly duplicated to give rise to two genes (Iris-A and Iris-B). Iris-B has significantly diverged from the Iris-A lineage, primarily because of the “invention” of an intron de novo in what was previously exonic sequence. Unlike domesticated retroviral envelope genes in mammals, we find that Iris has been subject to strong positive selection between Drosophila species. The rapid, adaptive evolution of Iris is sufficient to unambiguously distinguish the phylogenies of three closely related sibling species of Drosophila (D. simulans, D. sechellia, and D. mauritiana), a discriminative power previously described only for a putative “speciation gene.” Iris represents the first instance of a retroviral envelope–derived host gene outside vertebrates. It is also the first example of a retroviral envelope gene that has been found to be subject to positive selection following its domestication. The unusual selective pressures acting on Iris suggest that it is an active participant in an ongoing genetic conflict. We propose a model in which Iris has “switched sides,” having been recruited by host genomes to combat baculoviruses and retroviruses, which employ homologous envelope genes to mediate infection.

Mobile genetic elements have made homes within eukaryotic (host) genomes for hundreds of millions of years. These include retroviruses that integrate into host genomes as an essential step in their life cycle. While most such integration events are likely to be either deleterious or of little consequence to the host, on rare occasions host genomes can preserve and exploit capabilities of mobile elements for their own function. Especially intriguing are instances where host genomes have chosen to retain the envelope genes of retroviruses; the same envelope genes are responsible for conferring infectious ability to retroviruses. Primates and rodent genomes each have domesticated retroviral envelope genes (called “syncytin” genes) for important roles in placental function.

Now, Harmit Malik and colleagues show that a similar, ancient domestication event has taken place within the fruit fly Drosophila melanogaster. They identify a gene, Iris, which was acquired from an envelope gene of insect retroviruses, and has been maintained as a host gene for more than 25 million years. Unexpectedly, the authors find that Iris continues to evolve rapidly whereas previous studies have shown that mammalian syncytin genes do not. They suggest a model in which the Iris gene has “switched sides,” from its original role in causing infections to its current role in preventing them.

Basal relationships in the Drosophila melanogaster species group

The Drosophila melanogaster species group is a popular model for evolutionary studies due to its morphological and ecological diversity and its inclusion of the model species D. melanogaster. However, phylogenetic relationships among major lineages within this species group remain controversial. In this report, the phylogeny of 10 species representing each of the well-supported monophyletic clades in the melanogaster group was studied using the sequences of 14 loci that together comprise 9493 nucleotide positions. Combined Bayesian analysis using gene-specific substitution models produced a 100% credible set of two trees. In the strict consensus of these trees, the ananassae subgroup branches first in the melanogaster species group, followed by the montium subgroup. The remaining lineages form a monophyletic clade in which D. ficusphila and D. elegans branch first, followed by D. biarmipes, D. eugracilis, and the melanogaster subgroup. This strongly supported phylogeny resolves most basal relationships in the melanogaster species group, and provides a framework that can be extended in the future to encompass more species.

The expansion of the metazoan microRNA repertoire

BACKGROUND

MicroRNAs have been identified as crucial regulators in both animals and plants. Here we report on a comprehensive comparative study of all known miRNA families in animals. We expand the MicroRNA Registry 6.0 by more than 1000 new homologs of miRNA precursors whose expression has been verified in at least one species. Using this uniform data basis we analyze their evolutionary history in terms of individual gene phylogenies and in terms of preservation of genomic nearness across species. This allows us to reliably identify microRNA clusters that are derived from a common transcript.

RESULTS

We identify three episodes of microRNA innovation that correspond to major developmental innovations: A class of about 20 miRNAs is common to protostomes and deuterostomes and might be related to the advent of bilaterians. A second large wave of innovations maps to the branch leading to the vertebrates. The third significant outburst of miRNA innovation coincides with placental (eutherian) mammals. In addition, we observe the expected expansion of the microRNA inventory due to genome duplications in early vertebrates and in an ancestral teleost. The non-local duplications in the vertebrate ancestor are predated by local (tandem) duplications leading to the formation of about a dozen ancient microRNA clusters.

CONCLUSION

Our results suggest that microRNA innovation is an ongoing process. Major expansions of the metazoan miRNA repertoire coincide with the advent of bilaterians, vertebrates, and (placental) mammals.

Cold and heat tolerance of drosophilid flies with reference to their latitudinal distributions

The relation between thermal tolerance and latitudinal distribution was studied with 30 drosophilid species collected from the cool-temperate region (Sapporo), the warm-temperate region (Tokyo and Kyoto) and the subtropical region (Iriomote island) in Japan. In addition, intraspecific variation was examined for five species collected from two localities. The subtropical strains of Scaptodrosophila coracina, Drosophila bizonata and D. daruma were less tolerant to cold than their temperate strains. However, the difference of cold tolerance between these two geographic strains was much smaller than the difference between the species restricted to the subtropical region and those occurring in the temperate region. In D. auraria and D. suzukii, no difference was observed in thermal tolerance between their cool- and warm-temperate strains. Thus, geographic variation in thermal tolerance within species was low or negligible. Interspecific comparisons by phylogenetic independent contrasts revealed that species which had the northern boundaries of their distributions at higher latitudes were generally more tolerant to cold than those which had their boundaries at lower latitudes. However, the data for some species did not agree with this trend. The use of man-protected warm places for overwintering, competition or predation would also affect their distributions. It also appeared that species which had their southern boundaries at higher latitudes were generally more cold-tolerant. The acquisition of cold tolerance may lower a fly's capacity to compete, survive or reproduce in warmer climates. On the other hand, no relation was observed between heat tolerance and latitudinal distribution. Heat tolerance was higher in species inhabiting openlands or the forest canopy than in those inhabiting the forest understorey.