Coevolution between Male and Female Genitalia in the Drosophila melanogaster Species Subgroup

In contrast to male genitalia that typically exhibit patterns of rapid and divergent evolution among internally fertilizing animals, female genitalia have been less well studied and are generally thought to evolve slowly among closely-related species. As a result, few cases of male-female genital coevolution have been documented. In Drosophila, female copulatory structures have been claimed to be mostly invariant compared to male structures. Here, we re-examined male and female genitalia in the nine species of the D. melanogaster subgroup. We describe several new species-specific female genital structures that appear to coevolve with male genital structures, and provide evidence that the coevolving structures contact each other during copulation. Several female structures might be defensive shields against apparently harmful male structures, such as cercal teeth, phallic hooks and spines. Evidence for male-female morphological coevolution in Drosophila has previously been shown at the post-copulatory level (e.g., sperm length and sperm storage organ size), and our results provide support for male-female coevolution at the copulatory level.     

Evolution of DNA in Heterochromatin: the Drosophila Melanogaster Sibling Species Subgroup as a Resource

The Drosophila melanogasterspecies subgroup is a closely-knit collection of eight sibling species whose relationships are well defined. These species are too close for most evolutionary studies of euchromatic genes but are ideal to investigate the major changes that occur to DNA in heterochromatin over short periods during evolution. For example, it is not known whether the locations of genes in heterochromatin are conserved over this time. The 18S and 28S ribosomal RNA genes can be considered as genuine heterochromatic genes. In D. melanogasterthe rRNA genes are located at two sites, one each on the X and Y chromosome. In the other seven sibling species, rRNA genes are also located on the sex chromosomes but the positions often vary significantly, particularly on the Y. Furthermore, rDNA has been lost from the Y chromosome of both D. simulansand D. sechellia, presumably after separation of the line leading to present-day D. mauritiana.We conclude that changes to chromosomal position and copy number of rDNA arrays occur over much shorter evolutionary timespans than previously thought. In these respects the rDNA behaves more like the tandemly repeated satellite DNAs than euchromatic genes. This revised version was published online in July 2006 with corrections to the Cover Date.     

Comparative Genomics of Mitochondrial DNA in Members of the Drosophila melanogaster Subgroup

In this study, a comparative genomics approach is employed to investigate the forces that shape evolutionary change in the mitochondrial DNA (mtDNA) of members of the Drosophila melanogaster subgroup. This approach facilitates differentiation of the patterns of variation resulting from processes acting at a higher level from those acting on a single gene. The mitochondrial genomes of three isofemale lines of D. simulans (siI, -II, and -III), two of D. melanogaster (Oregon R and a line from Zimbabwe), and D. mauritiana (maI and -II), and one of D. sechellia were sequenced and compared with that derived from D. yakuba. Data presented here indicate that at least three broad mechanisms shape the evolutionary dynamics of mtDNA in these taxa. The first set of mechanisms is intrinsic to the molecule. Dominant processes may be interpreted as selection for an increased rate of replication of the mtDNA molecule, biases in DNA repair, and differences in the pattern of nucleotide substitution among strands. In the genes encoded on the major strand (62% of the coding DNA) changes to or from C predominate, whereas on the minor changes to or from G predominate. The second set of mechanisms affects distinct lineages. There are evolutionary rate differences among lineages, possibly owing to population demographic changes or changes in mutational biases. This is supported by the heterogeneity found in synonymous, nonsynonymous, and silent substitutions. The third set of mechanisms differentially affects distinct genes. A maximum-likelihood sliding-window analysis detected four disjunct regions that have a significantly different nucleotide substitution process from that derived from the complete sequence. These data show the potential for comparative genomics to tease apart subtle forces that shape the evolution of DNA. Received: 30 July 1999 / Accepted: 16 March 2000     

中国特有果蝇 Drosophila curviceps 种亚组在Drosophila immigrans种组中的种系发生地位

果蝇immigrans种组中的curviceps种亚组是1992年新建立的中国特有果蝇类群。该种亚组中的物种主要分布在中国大陆和台湾。目前除了形态学水平的研究外,还没有其他证据支持建立该种亚组的合理性及其起源和种系发生地位。为了在DNA分子水平上探讨果蝇curviceps种亚组在果蝇immigrans种组中的种系发生地位,从而为今后更深入地研究中国特有果蝇,甚至为果蝇亚属的进化遗传学提供理论依据,测定了immigrans种组5个种亚组（nasuta、immigrans、hypocausta、quadrilineata、curviceps)中12个代表物种的rDNA的ITS1和部分Adh基因的序列。其中ITS1序列的长度为513－587bp,共有191个信息位点;Adh基因片段的长度在714－747bp之间,共99个信息位点。考虑到单个分子提供的信息较少,将两个分子的序列综合起来,组成一个较长的复合序列。分别根据ITS1,Adh和两个分子的复合序列排比（Alignment)结果,和最大简约法和邻接法构建分子系统树,其中根据复合序列构建的系统树与形态学研究结果最为一致。分子树显示curviceps种亚组的特种确定单独形成一个分枝,为种亚组级的分类阶元,支持了形态学将其建立为一个新种亚组。根据Kimura距离,估算了复合分子的替换速率约为每百万年1．48％,进而计算出5个种亚组的分 歧年代。结合各物种的地理分布,推测了immigrans种组的进化历史：curviceps种亚组与quadrilineata种亚组的亲缘关系最近,主要分布在中国南部的温带地区。它们之间的分歧时间大约为3．4百万年,是最年轻的两个种亚组。主要分布在苏门答腊及附近的热带地区的hypocausta种亚组的物种是最早分化出来的,与其他种亚组的分歧时间约为9．2百万年。该结果与形态学和生物地理学研究相吻合。值得一提是的,目前归属仍存在争议的物种D.neohypocausta,在分子系统树中与hypocausta种亚组的物种相距较远,而与immiagrasn种亚组的关系较近,但分枝置信度较低（＜50％）。由于还缺乏其他方面的证据,因此D.neohypocausta的归属有待今后的研究来作定论。     

Evolution of late-life fecundity in Drosophila melanogaster

Late-life fecundity has been shown to plateau at late ages in Drosophila analogously to late-life mortality rates. In this study, we test an evolutionary theory of late life based on the declining force of natural selection that can explain the occurrence of these late-life plateaus in Drosophila. We also examine the viability of eggs laid by late-age females and test a population genetic mechanism that may be involved in the evolution of late-life fecundity: antagonistic pleiotropy. Together these experiments demonstrate that (i) fecundity plateaus at late ages, (ii) plateaus evolve according to the age at which the force of natural selection acting on fecundity reaches zero, (iii) eggs laid by females in late life are viable and (iv) antagonistic pleiotropy is involved in the evolution of late-life fecundity. This study further supports the evolutionary theory of late life based on the age-specific force of natural selection.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Identification of One Intron Loss and Phylogenetic Evolution of Dfak Gene in the Drosophila melanogaster Species Group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5–10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.     

Evolution of the olfactory code in the Drosophila melanogaster subgroup

The Drosophila melanogaster subgroup has been the focus of numerous studies about evolution. We address the question of how the olfactory code has evolved among the nine sister species. By using in vivo electrophysiological measurements, so called single-cell recordings, we have established the ligand affinity of a defined subset of olfactory receptor neurons (ORNs) across all nine species. We show that the olfactory code as relayed by the investigated subset of ORNs is conserved to a striking degree. Distinct shifts in the code have occurred only within the simulans clade. However, these shifts are restricted to an altered tuning profile of the same single ORN type in all three of the simulans siblings and a more drastic change unique to D. sechellia, involving a complete loss of one sensillum type in favour of another. The alterations observed in D. sechellia may represent a novel host-specific adaptation to its sole host, morinda fruit (Morinda citrifolia). The overall high degree of similarity of the code within the subgroup is intriguing when considering the great variety in distributions as well as in habitat and host choice of the siblings, factors that could greatly affect the olfactory system.     

Vertical Transmission of the Retrotransposable Elements R1 and R2 during the Evolution of the Drosophila Melanogaster Species Subgroup

R1 and R2 are non-long-terminal repeat retrotransposable elements that insert into specific sequences of insect 28S ribosomal RNA genes. These elements have been extensively described in Drosophila melanogaster. To determine whether these elements have been horizontally or vertically transmitted, we characterized R1 and R2 elements from the seven other members of the melanogaster species subgroup by genomic blotting and nucleotide sequencing. Each species was found to have homogeneous families of R1 and R2 elements with the exception of erecta and orena, which have no R2 elements. The DNA sequences of multiple R1 and R2 copies from each species indicated nucleotide divergence within each species averaged only 0.48% for R1 and 0.35% for R2, well below the level of divergence among the species. Most copies of R1 and R2 (40 of 47) sequenced from the seven species were potentially functional, as indicated by the absence of premature termination codons or translational frameshifts that would destroy the open reading frame of the element. The sequence relationships of both the R1 and R2 elements from the various members of the melanogaster subgroup closely followed that of the species phylogeny, suggesting that R1 and R2 have been stably maintained by vertical transmission since the origin of this species subgroup 17-20 million years ago. The remarkable stability of R1 and R2, compared to what has been suggested for transposable elements that insert at multiple locations in these same species, may be due to their unique specificity for sites in the rRNA gene locus. Under low copy number conditions, when it is essential for any mobile element to transpose, the insertion specificities of R1 and R2 ensure uniform developmentally regulated target sites that can be occupied with little or no detrimental effect on the host.     

果蝇nasuta亚群求爱歌的种间识别与进化遗传学研究

果蝇ｎａｓｕｔａ亚群由１４个种、亚种和分类群组成,广泛分布于印度－太平洋区域。本文首次记录了ｎａｓｕｔａ亚群种的求爱歌,测量了脉冲歌时域模式的参数：脉冲串间隔（ＩＢＩ）、脉冲间隔（ＩＰＩ）、脉冲串时间长度（ＰＴＬ）、每个脉冲串的脉冲数（ＰＮ）、脉冲时间长度（ＰＬ）、波动周期时间长度（ＣＬ）。采用计算机声谱分析技术,作出求爱歌信号的三维数字功率谱图,进行频率分析。发现Ｄ．ｐｕｌａｕｎａ和Ｔａｘｏｎ－Ｆ不发出求爱歌声信号,视觉在交配中可能起重要作用。对其余种、亚种和分类群的求爱歌分析表明,ｎａｓｕｔａ亚群种的求爱歌分为脉冲歌和正弦歌。对部分种的正反交Ｆ１求爱歌分析表明,脉冲歌时域参数,如ＩＰＩ平均值为Ｘ染色体连锁或常染色体多基因控制,正弦歌频率偏向母方。根据不同种、亚种和分类群脉冲歌的时域模式构建ｎａｓｕｔａ亚群的系统树,对亚群中不同种、亚种和分类群的亲缘关系进行讨论。     

Lineage-Specific Evolution of the Complex Nup160 Hybrid Incompatibility Between Drosophila melanogaster and Its Sister Species

Two genes encoding protein components of the nuclear pore complex Nup160 and Nup96 cause lethality in F₂-like hybrid genotypes between Drosophila simulans and Drosophila melanogaster. In particular, D. simulans Nup160 and Nup96 each cause inviability when hemizygous or homozygous in species hybrids that are also hemizygous (or homozygous) for the D. melanogaster X chromosome. The hybrid lethality of Nup160, however, is genetically complex, depending on one or more unknown additional factors in the autosomal background. Here we study the genetics and evolution of Nup160-mediated hybrid lethality in three ways. First, we test for variability in Nup160-mediated hybrid lethality within and among the three species of the D. simulans clade— D. simulans, D. sechellia, and D. mauritiana. We show that the hybrid lethality of Nup160 is fixed in D. simulans and D. sechellia but absent in D. mauritiana. Second, we explore how the hybrid lethality of Nup160 depends on other loci in the autosomal background. We find that D. simulans Nup160-mediated hybrid lethality does not depend on the presence of D. melanogaster Nup96, and we find that D. simulans and D. mauritiana are functionally differentiated at Nup160 as well as at other autosomal factor(s). Finally, we use population genetics data to show that Nup160 has experienced histories of recurrent positive selection both before and after the split of the three D. simulans clade species ∼240,000 years ago. Our genetic results suggest that a hybrid lethal Nup160 allele evolved before the split of the three D. simulans clade species, whereas the other autosomal factor(s) evolved more recently.     

Molecular Characterization and Evolution of the Amylase Multigene Family of Drosophila ananassae

Drosophila ananassae is known to produce numerous alpha-amylase variants. We have cloned seven different Amy genes in an African strain homozygous for the AMY1,2,3,4 electrophoretic pattern. These genes are organized as two main clusters: the first one contains three intronless copies on the 2L chromosome arm, two of which are tandemly arranged. The other cluster, on the 3L arm, contains two intron-bearing copies. The amylase variants AMY1 and AMY2 have been assigned to the intronless cluster, and AMY3 and AMY4 to the second one. The divergence of coding sequences between clusters is moderate (6.1% in amino acids), but the flanking regions are very different, which could explain their differential regulation. Within each cluster, coding and noncoding regions are conserved. Two very divergent genes were also cloned, both on chromosome 3L, but very distant from each other and from the other genes. One is the Amyrel homologous (41% divergent), the second one, Amyc1 (21.6% divergent) is unknown outside the D. ananassae subgroup. These two genes have unknown functions. Received: 30 May 2000 / Accepted: 17 July 2000     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.     

Evolution of the Gypsy endogenous retrovirus in the Drosophila melanogaster subgroup

We conducted a phylogenetic survey of the endogenous retrovirus Gypsy in the eight species of the Drosophila melanogaster subgroup. A 362-bp fragment from the integrase gene (int) was amplified, cloned, and sequenced. Phylogenetic relationships of the elements isolated from independent clones were compared with the host phylogeny. Our results indicate that two main lineages of Gypsy exist in the melanogaster subgroup and that vertical and horizontal transmission have played a crucial role in the evolution of this insect endogenous retrovirus.     

Strong Purifying Selection on the Odysseus Gene in Two Clades of Sibling Species of the Drosophila montium Species Subgroup

The Odysseus (OdsH) gene was duplicated from its ancestral neuron-expressed gene, unc-4, and then evolved very rapidly under strong positive Darwinian selection as a speciation gene causing hybrid-male sterility between closely related species of the Drosophila simulans clade. Has OdsH also experienced similar positive selection between Drosophila sibling species other than those of the simulans clade? We cloned and sequenced OdsH and unc-4 from two clades of the Drosophila montium species subgroup, the Drosophila lini and the Drosophila kikkawai clades. The ratios of Ka/Ks for OdsH were remarkably low between sibling species of these two clades, suggesting that OdsH has been subjected to strong purifying selection in these two clades.     

Molecular Evolution of the Duplicated Amy Locus in the Drosophila Melanogaster Species Subgroup: Concerted Evolution Only in the Coding Region and an Excess of Nonsynonymous Substitutions in Speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a ``specialist' species for host plants and probably also occurs in the process of speciation in general.