Evolution of yellow gene regulation and pigmentation in Drosophila

BACKGROUND: Changes in developmental gene expression are central to phenotypic evolution, but the genetic mechanisms underlying these changes are not well understood. Interspecific differences in gene expression can arise from evolutionary changes in cis-regulatory DNA and/or in the expression of trans-acting regulatory proteins, but few case studies have distinguished between these mechanisms. Here, we compare the regulation of the yellow gene, which is required for melanization, among distantly related Drosophila species with different pigment patterns and determine the phenotypic effects of divergent Yellow expression. RESULTS: Yellow expression has diverged among D. melanogaster, D. subobscura, and D. virilis and, in all cases, correlates with the distribution of black melanin. Species-specific Yellow expression patterns were retained in D. melanogaster transformants carrying the D. subobscura and D. virilis yellow genes, indicating that sequence evolution within the yellow gene underlies the divergence of Yellow expression. Evolutionary changes in the activity of orthologous cis-regulatory elements are responsible for differences in abdominal Yellow expression; however, cis-regulatory element evolution is not the sole cause of divergent Yellow expression patterns. Transformation of the D. melanogaster yellow gene into D. virilis altered its expression pattern, indicating that trans-acting factors that regulate the D. melanogaster yellow gene have also diverged between these two species. Finally, we found that the phenotypic effects of evolutionary changes in Yellow expression depend on epistatic interactions with other genes. CONCLUSIONS: Evolutionary changes in Yellow expression correlate with divergent melanin patterns and are a result of evolution in both cis- and trans-regulation. These changes were likely necessary for the divergence of pigmentation, but evolutionary changes in other genes were also required.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Recent origins of sperm genes in Drosophila

Newly created genes often acquire testis-specific or enhanced expression but neither the mechanisms responsible for this specificity nor the functional consequences of these evolutionary processes are well understood. Genomic analyses of the Drosophila melanogaster sperm proteome has identified 2 recently evolved gene families on the melanogaster lineage and 4 genes created by retrotransposition during the evolution of the melanogaster group that encode novel sperm components. The expanded Mst35B (protamine) and tektin gene families are the result of tandem duplication events with all family members displaying testis-specific expression. The Mst35B family encodes rapidly evolving protamines that display a robust signature of positive selection within the DNA-binding high-mobility group box consistent with functional diversification in genome repackaging during sperm nuclear remodeling. The Mst35B paralogs also reside in a significant regional cluster of testis-overexpressed genes. Tektins, known components of the axoneme, are encoded by 3 nearly identical X-linked genes, a finding consistent with very recent gene family expansion. In addition to localized duplication events, the evolution of the sperm proteome has also been driven by recent retrotransposition events resulting in Cdlc2, CG13340, Vha36, and CG4706. Cdlc2, CG13340, and Vha36 all display high levels of overexpression in the testis, and Cdlc2 and CG13340 reside within testis-overexpressed gene clusters. Thus, gene creation is a dynamic force in the evolution of sperm composition and possibly function, which further suggests that acquisition of molecular functionality in sperm may be an influential pathway in the fixation of new genes.     

Protein subcellular relocalization: a new perspective on the origin of novel genes

Gene duplication is considered to be the most important evolutionary process for generating novel genes. However, the mechanisms involved in the evolution of such genetic innovations remain unclear. There is compelling evidence to suggest that changing the subcellular location of a protein can also alter its function, and that diversity in subcellular targeting within gene families is common. Here, we introduce the idea that protein subcellular relocalization might be an important evolutionary mechanism for the origins of new genes.     

Population genetics of duplicated disease-defense genes,hm1 and hm2, in maize (Zea mays ssp. mays L.) and its wild ancestor (Zea mays ssp. parviglumis)

Plant defense genes are subject to nonneutral evolutionary dynamics. Here we investigate the evolutionary dynamics of the duplicated defense genes hm1 and hm2 in maize and its wild ancestor Zea mays ssp. parviglumis. Both genes have been shown to confer resistance to the fungal pathogen Cochliobolus carbonum race 1, but the effectiveness of resistance differs between loci. The genes also display different population histories. The hm1 locus has the highest nucleotide diversity of any gene yet sampled in the wild ancestor of maize, and it contains a large number of indel polymorphisms. There is no evidence, however, that high diversity in hm1 is a product of nonneutral evolution. In contrast, hm2 has very low nucleotide diversity in the wild ancestor of maize. The distribution of hm2 polymorphic sites is consistent with nonneutral evolution, as indicated by Tajima's D and other neutrality tests. In addition, one hm2 haplotype is more frequent than expected under the equilibrium neutral model, suggesting hitchhiking selection. Both defense genes retain >80% of the level of genetic variation in maize relative to the wild ancestor, and this level is similar to other maize genes that were not subject to artificial selection during domestication.     

Molecular evolution of plant immune system genes

Molecular population genetic studies are providing new perspectives on the evolution of genes that confer resistance to pathogens and herbivores. Here, we compare the evolutionary history of different components of the defense response (detection, signaling and response) and of genes with parallel function in plants and Drosophila. A review of the literature indicates that the dominant form of selection acting on defense genes (balancing, positive and purifying) differs among components of defense. Sampling of particular classes of genes and genes from non-model organisms, however, remains limited. Future studies combining molecular evolutionary analyses with ecological genetic and functional analyses should better reveal how natural selection has shaped defense gene evolution.     

Ephemeral Association Between Gene CG5762 and Hybrid Male Sterility in Drosophila Sibling Species

Interspecies divergence in regulatory pathways may result in hybrid male sterility (HMS) when dominance and epistatic interactions between alleles that are functional within one genome are disrupted in hybrid genomes. The identification of genes contributing to HMS and other hybrid dysfunctions is essential for understanding the origin of new species (speciation). Previously, we identified a panel of male-specific loci misexpressed in sterile male hybrids of Drosophila simulans and D. mauritiana relative to parental species. In the current work, we attempt to dissect the genetic associations between HMS and one of the genes, CG5762, a Drosophila-unique locus characterized by rapid sequence divergence within the genus, presumably driven by positive natural selection. CG5762 is underexpressed in sterile backcross males compared with their fertile brothers. In CG5762 heterozygotes, the D. mauritiana allele is consistently overexpressed on both the D. simulans and D. mauritiana backcross genomic background, suggesting a cis-acting regulation factor. There is a significant association between heterozygosity and HMS in hybrid males from early but not later backcross generations. Microsatellite markers spanning CG5762 fail to associate with HMS. These observations lead to a conclusion that CG5762 is not a causative factor of HMS. Although genetic linkage between CG5762 and a neighboring causative introgression cannot be ruled out, it seems that the pattern is most consistent with CG5762 participating in epistatic interactions that are disrupted in flies with HMS.     

Interspecies comparison of a gene pair with partially redundant function: the rst and kirre genes in D. virilis and D. melanogaster

Abstract The D. melanogaster rst and kirre genes encode two highly related immunoglobulin-like cell adhesion molecules that function redundantly during embryonic muscle development. The two genes appear to be derived from a common ancestor by gene duplication. Gene duplications have been proposed to be of major evolutionary significance since duplicated redundant sequences can accumulate mutations without detrimental effects for the organism and leave the duplicated genes free to assume novel functions. To address the issue of conservation of the duplicated sequences and their putative redundancy, as well as to identify putative functional divergence of the paralogs during drosophilid evolution, we performed an interspecies comparison of the rst and kirre genes from D. virilis and D. melanogaster. The D. virilis genome contains orthologues of both rst and kirre and hence the duplication took place before the split of the two lineages and has subsequently been conserved. However, whilst the Rst orthologues show a high degree of sequence similarity, this similarity is lower in Kirre orthologues. Especially the intracellular domains of D. virilis and D. melanogaster Kirre sequences are highly divergent: the D. virilis kirre gene lacks the 3′-most exon present in D. melanogaster, which contains motifs conserved between kirre and rst in D. melanogaster. Hence, while each of the two genes is highly conserved at the level of its exon-intron organization, the selection forces acting on the rst and kirre coding sequences are different. These findings are discussed in the light of general evolutionary mechanisms.     

Tandem Repetitive D Domains of the Sperm Ligand Zonadhesin Evolve Faster in the Paralogue Than in the Orthologue Comparison

Gene duplication is regarded as an important evolutionary mechanism creating genetic and phenotypic novelty. At the same time, the evolutionary mechanisms following gene duplication have been a subject of much debate. Here we analyze the sequence evolution of zonadhesin, a mammalian sperm ligand that binds to the oocyte zona pellucida in a species-specific manner. In pig, rabbit, and primates, precursor zonadhesin comprises, among others, one partial and four complete tandem repetitive D domains. The mouse precursor is distinguished by 20 additional partial D3 domains consisting of 120 amino acids each. This gene structure allows sequence comparison in both paralogues and orthologues. Detailed sequence analysis reveals that D domains evolve faster across paralogues than orthologues. Moreover, at the codon level, partial D3 paralogues of mouse show evidence of positive selection, whereas the corresponding orthologues do not. Individual posttranslational motif patterns and positive selection point to neofunctionalization of partial D3 paralogues of mouse, rather than subfunctionalization. However, as we found additional evidence for homogenization by partial gene conversion, sequence evolution of partial D3 paralogues of mouse might be better described as a combination of divergent and convergent evolution. So far, the divergence at the codon level has outbalanced the convergence at the level of smaller fragments. The probable driving force behind the evolutionary patterns observed is sexual selection. We finally discuss whether the functional determination influences the evolutionary regime acting on sperm ligands and egg receptors, respectively. [Reviewing Editor: Dr. Yves Van de Peer]     

Evolution of gene sequence in response to chromosomal location

Evolutionary forces acting on the repetitive DNA of heterochromatin are not constrained by the same considerations that apply to protein-coding genes. Consequently, such sequences are subject to rapid evolutionary change. By examining the Troponin C gene family of Drosophila melanogaster, which has euchromatic and heterochromatic members, we find that protein-coding genes also evolve in response to their chromosomal location. The heterochromatic members of the family show a reduced CG content and increased variation in DNA sequence. We show that the CG reduction applies broadly to the protein-coding sequences of genes located at the heterochromatin:euchromatin interface, with a very strong correlation between CG content and the distance from centric heterochromatin. We also observe a similar trend in the transition from telomeric heterochromatin to euchromatin. We propose that the methylation of DNA is one of the forces driving this sequence evolution.     

Evolution of a highly polymorphic human cytochrome P450 gene cluster: CYP2D6.

The CYP2D gene cluster on human chromosome 22 containing the functional cytochrome P450 gene CYP2D6 and two or three highly homologous pseudogenes is involved in a clinically important variation in the inactivation of drugs and environmental chemicals. Several mutant haplotypes of CYP2D6 have been identified by restriction analysis and by PCR-based allele-specific amplification. To understand the evolutionary sequence of mutational events as well as recently discovered interracial differences, we analyzed the arrangement of the CYP2D haplotype containing a common mutant allele of CYP2D6 associated with a XbaI 44-kb fragment. This haplotype contains four CYP2D genes instead of three. Comparison of the sequences of these genes with those of previously characterized haplotypes suggests that an early point mutation was followed by a crossover and a gene conversion event, the latter found preferentially in Caucasians. These data are consistent with the rapid evolution of this locus during "plant-animal warfare" with practical consequences for present-day defense of the organism against environmental adversity.     

Identification of a novel target of D/V signaling in Drosophila wing disc: Wg-independent function of the organizer

Growth and patterning during Drosophila wing development are mediated by signaling from its dorso-ventral (D/V) organizer. Wingless is expressed in the D/V boundary and functions as a morphogen to activate target genes at a distance. Wingless pathway and thereby D/V signaling is negatively regulated by the homeotic gene Ultrabithorax (Ubx) to mediate haltere development. In an enhancer-trap screen to identify genes that show differential expression between wing and haltere discs, we identified CG32062, which codes for a RNA-binding protein. In wing discs, CG32062 is expressed only in non-D/V cells. CG32062 expression in non-D/V cells is dependent on Notch-mediated signaling from the D/V boundary. However, CG32062 expression is independent of Wingless function, thus providing evidence for a second long-range signaling mechanism of the D/V organizer. In haltere discs, CG32062 is negatively regulated by Ubx. The non-cell autonomous nature of Ubx-mediated repression of CG32062 expression suggests that the novel component of D/V signaling is also negatively regulated during haltere specification.     

Evidences of lateral gene transfer between archaea and pathogenic bacteria

Acquisition of new genetic material through horizontal gene transfer has been shown to be an important feature in the evolution of many pathogenic bacteria. Changes in the genetic repertoire, occurring through gene acquisition and deletion, are the major events underlying the emergence and evolution of bacterial pathogens. However, horizontal gene transfer across the domains i.e. archaea and bacteria is not so common. In this context, we explore events of horizontal gene transfer between archaea and bacteria. In order to determine whether the acquisition of archaeal genes by lateral gene transfer is an important feature in the evolutionary history of the pathogenic bacteria, we have developed a scheme of stepwise eliminations that identifies archaeal-like genes in various bacterial genomes. We report the presence of 9 genes of archaeal origin in the genomes of various bacteria, a subset of which is also unique to the pathogenic members and are not found in respective non-pathogenic counterparts. We believe that these genes, having been retained in the respective genomes through selective advantage, have key functions in the organism’s biology and may play a role in pathogenesis.