Gene duplication and speciation in Drosophila: evidence from the Odysseus locus

Selection with two alleles of X-linkage and its application to the fitness component analysis of OdsH in Drosophila

In organisms with the XY sex-determination system, there is an imbalance in the inheritance and transmission of the X chromosome between males and females. Unlike an autosomal allele, an X-linked recessive allele in a female will have phenotypic effects on its male counterpart. Thus, genes located on the X chromosome are of particular interest to researchers in molecular evolution and genetics. Here we present a model for selection with two alleles of X-linkage to understand fitness components associated with genes on the X chromosome. We apply this model to the fitness analysis of an X-linked gene, OdsH (16D), in the fruit fly Drosophila melanogaster. The function of OdsH is involved in sperm production and the gene is rapidly evolving under positive selection. Using site-directed gene targeting, we generated functional and defective OdsH variants tagged with the eye-color marker gene white. We compare the allele frequency changes of the two OdsH variants, each directly competing against a wild-type OdsH allele in concurrent but separate experimental populations. After 20 generations, the two genetically modified OdsH variants displayed a 40% difference in allele frequencies, with the functional OdsH variant demonstrating an advantage over the defective variant. Using maximum likelihood estimation, we determined the fitness components associated with the OdsH alleles in males and females. Our analysis revealed functional aspects of the fitness determinants associated with OdsH, and that sex-specific fertility and viability consequences both contribute to selection on an X-linked gene.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Structural Variants and Speciation: Multiple Processes at Play

Research on the genomic architecture of speciation has increasingly revealed the importance of structural variants (SVs) that affect the presence, abundance, position, and/or direction of a nucleotide sequence. SVs include large chromosomal rearrangements such as fusion/fissions and inversions and translocations, as well as smaller variants such as duplications, insertions, and deletions (CNVs). Although we have ample evidence that SVs play a key role in speciation, the underlying mechanisms differ depending on the type and length of the SV, as well as the ecological, demographic, and historical context. We review predictions and empirical evidence for classic processes such as underdominance due to meiotic aberrations and the coupling effect of recombination suppression before exploring how recent sequencing methodologies illuminate the prevalence and diversity of SVs. We discuss specific properties of SVs and their impact throughout the genome, highlighting that multiple processes are at play, and possibly interacting, in the relationship between SVs and speciation.

The role of gene presence-absence variations on genetic incompatibility in Asian rice

Genetic incompatibilities are widespread between species. However, it remains unclear whether they all originated after population divergence as suggested by the Bateson-Dobzhansky-Muller model, and if not, what is their prevalence and distribution within populations. The gene presence-absence variations (PAVs) provide an opportunity for investigating gene-gene incompatibility. Here, we searched for the repulsion of coexistence between gene PAVs to identify the negative interaction of gene functions separately in two Oryza sativa subspecies. Many PAVs are involved in subspecies-specific negative epistasis and segregate at low-to-intermediate frequencies in focal subspecies but at low or high frequencies in the other subspecies. Incompatible PAVs are enriched in two functional groups, defense response and protein phosphorylation, which are associated with plant immunity and consistent with autoimmunity being a known mechanism of hybrid incompatibility in plants. Genes in the two enriched functional groups are older and seldom directly interact with each other. Instead, they interact with other younger gene PAVs with diverse functions. Our results illustrate the landscape of genetic incompatibility at gene PAVs in rice, where many incompatible pairs have already segregated as polymorphisms within subspecies, and many are novel negative interactions between older defense-related genes and younger genes with diverse functions.

X-chromosome target specificity diverged between dosage compensation mechanisms of two closely related Caenorhabditis species

An evolutionary perspective enhances our understanding of biological mechanisms. Comparison of sex determination and X-chromosome dosage compensation mechanisms between the closely related nematode species Caenorhabditis briggsae (Cbr) and Caenorhabditis elegans (Cel) revealed that the genetic regulatory hierarchy controlling both processes is conserved, but the X-chromosome target specificity and mode of binding for the specialized condensin dosage compensation complex (DCC) controlling X expression have diverged. We identified two motifs within Cbr DCC recruitment sites that are highly enriched on X: 13 bp MEX and 30 bp MEX II. Mutating either MEX or MEX II in an endogenous recruitment site with multiple copies of one or both motifs reduced binding, but only removing all motifs eliminated binding in vivo. Hence, DCC binding to Cbr recruitment sites appears additive. In contrast, DCC binding to Cel recruitment sites is synergistic: mutating even one motif in vivo eliminated binding. Although all X-chromosome motifs share the sequence CAGGG, they have otherwise diverged so that a motif from one species cannot function in the other. Functional divergence was demonstrated in vivo and in vitro. A single nucleotide position in Cbr MEX can determine whether Cel DCC binds. This rapid divergence of DCC target specificity could have been an important factor in establishing reproductive isolation between nematode species and contrasts dramatically with the conservation of target specificity for X-chromosome dosage compensation across Drosophila species and for transcription factors controlling developmental processes such as body-plan specification from fruit flies to mice.

Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

The Importance of Gene Duplication and Domain Repeat Expansion for the Function and Evolution of Fertilization Proteins

The process of gene duplication followed by gene loss or evolution of new functions has been studied extensively, yet the role gene duplication plays in the function and evolution of fertilization proteins is underappreciated. Gene duplication is observed in many fertilization protein families including Izumo, DCST, ZP, and the TFP superfamily. Molecules mediating fertilization are part of larger gene families expressed in a variety of tissues, but gene duplication followed by structural modifications has often facilitated their cooption into a fertilization function. Repeat expansions of functional domains within a gene also provide opportunities for the evolution of novel fertilization protein. ZP proteins with domain repeat expansions are linked to species-specificity in fertilization and TFP proteins that experienced domain duplications were coopted into a novel sperm function. This review outlines the importance of gene duplications and repeat domain expansions in the evolution of fertilization proteins.

Chromosome-level genome assembly of Bactrocera dorsalis reveals its adaptation and invasion mechanisms

Bactrocera dorsalis is an invasive polyphagous pest causing considerable ecological and economic damage worldwide. We report a high-quality chromosome-level genome assembly and combine various transcriptome data to explore the molecular mechanisms of its rapid adaptation to new environments. The expansions of the DDE transposase superfamily and key gene families related to environmental adaptation and enrichment of the expanded and unique gene families in metabolism and defence response pathways explain its environmental adaptability. The relatively high but not significantly different expression of heat-shock proteins, regardless of the environmental conditions, suggests an intrinsic mechanism underlying its adaptation to high temperatures. The mitogen-activated protein kinase pathway plays a key role in adaptation to new environments. The prevalence of duplicated genes in its genome explains the diversity in the B. dorsalis complex. These findings provide insights into the genetic basis of the invasiveness and diversity of B. dorsalis, explaining its rapid adaptation and expansion.

Jiang et al. sequence the genome of Bactrocera dorsalis, a destructive and invasive agricultural pest. Insights from this chromosome-level assembly shed light on molecular adaptations that allow for the global invasion and expansion of this pest.

Functional Diversification of Chromatin on Rapid Evolutionary Timescales

Repeat-enriched genomic regions evolve rapidly and yet support strictly conserved functions like faithful chromosome transmission and the preservation of genome integrity. The leading resolution to this paradox is that DNA repeat-packaging proteins evolve adaptively to mitigate deleterious changes in DNA repeat copy number, sequence, and organization. Exciting new research has tested this model of coevolution by engineering evolutionary mismatches between adaptively evolving chromatin proteins of one species and the DNA repeats of a close relative. Here, we review these innovative evolution-guided functional analyses. The studies demonstrate that vital, chromatin-mediated cellular processes, including transposon suppression, faithful chromosome transmission, and chromosome retention depend on species-specific versions of chromatin proteins that package species-specific DNA repeats. In many cases, the ever-evolving repeats are selfish genetic elements, raising the possibility that chromatin is a battleground of intragenomic conflict.

How Important Are Structural Variants for Speciation?

Understanding the genetic basis of reproductive isolation is a central issue in the study of speciation. Structural variants (SVs); that is, structural changes in DNA, including inversions, translocations, insertions, deletions, and duplications, are common in a broad range of organisms and have been hypothesized to play a central role in speciation. Recent advances in molecular and statistical methods have identified structural variants, especially inversions, underlying ecologically important traits; thus, suggesting these mutations contribute to adaptation. However, the contribution of structural variants to reproductive isolation between species—and the underlying mechanism by which structural variants most often contribute to speciation—remain unclear. Here, we review (i) different mechanisms by which structural variants can generate or maintain reproductive isolation; (ii) patterns expected with these different mechanisms; and (iii) relevant empirical examples of each. We also summarize the available sequencing and bioinformatic methods to detect structural variants. Lastly, we suggest empirical approaches and new research directions to help obtain a more complete assessment of the role of structural variants in speciation.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

The Genetic Architecture of Post-Zygotic Reproductive Isolation Between Anopheles coluzzii and An. quadriannulatus

The Anopheles gambiae complex is comprised of eight morphologically indistinguishable species and has emerged as a model system for the study of speciation genetics due to the rapid radiation of its member species over the past two million years. Male hybrids between most An. gambiae complex species pairs are sterile, and some genotype combinations in hybrid males cause inviability. We investigated the genetic basis of hybrid male inviability and sterility between An. coluzzii and An. quadriannulatus by measuring segregation distortion and performing a QTL analysis of sterility in a backcross population. Hybrid males were inviable if they inherited the An. coluzzii X chromosome and were homozygous at one or more loci in 18.9 Mb region of chromosome 3. The An. coluzzii X chromosome has a disproportionately large effect on hybrid sterility when introgressed into an An. quadriannulatus genetic background. Additionally, an epistatic interaction between the An. coluzzii X and a 1.12 Mb, pericentric region of the An. quadriannulatus 3L chromosome arm has a statistically significant contribution to the hybrid sterility phenotype. This same epistatic interaction occurs when the An. coluzzii X is introgressed into the genetic background of An. arabiensis, the sister species of An. quadriannulatus, suggesting that this may represent one of the first Dobzhansky-Muller incompatibilities to evolve early in the radiation of the Anopheles gambiae species complex. We describe the additive effects of each sterility QTL, epistatic interactions between them, and genes within QTL with protein functions related to mating behavior, reproduction, spermatogenesis, and microtubule morphogenesis, whose divergence may contribute to post-zygotic reproductive isolation between An. coluzzii and An. quadriannulatus.

Speciation and changes in male gene expression in Drosophila

It has long been acknowledged that changes in the regulation of gene expression may account for major organismal differences. However, we still do not fully understand how changes in gene expression evolve and how do such changes influence organisms' differences. We are even less aware of the impact such changes might have in restricting gene flow between species. Here, we focus on studies of gene expression and speciation in the Drosophila model. We review studies that have identified gene interactions in post-mating reproductive isolation and speciation, particularly those that modulate male gene expression. We also address studies that have experimentally manipulated changes in gene expression to test their effect in post-mating reproductive isolation. We highlight the need for a more in-depth analysis of the role of selection causing disrupted gene expression of such candidate genes in sterile/inviable hybrids. Moreover, we discuss the relevance to incorporate more routinely assays that simultaneously evaluate the potential effects of environmental factors and genetic background in modulating plastic responses in male genes and their potential role in speciation.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Genomic Structural Variations Within Five Continental Populations of Drosophila melanogaster

Chromosomal structural variations (SV) including insertions, deletions, inversions, and translocations occur within the genome and can have a significant effect on organismal phenotype. Some of these effects are caused by structural variations containing genes. Large structural variations represent a significant amount of the genetic diversity within a population. We used a global sampling of Drosophila melanogaster (Ithaca, Zimbabwe, Beijing, Tasmania, and Netherlands) to represent diverse populations within the species. We used long-read sequencing and optical mapping technologies to identify SVs in these genomes. Among the five lines examined, we found an average of 2,928 structural variants within these genomes. These structural variations varied greatly in size and location, included many exonic regions, and could impact adaptation and genomic evolution.

Novel Genomic and Evolutionary Perspective of Cyanobacterial tRNAs

Transfer RNA (tRNA) plays a central role in protein synthesis and acts as an adaptor molecule between an mRNA and an amino acid. A tRNA has an L-shaped clover leaf-like structure and contains an acceptor arm, D-arm, D-loop, anti-codon arm, anti-codon loop, variable loop, Ψ-arm and Ψ-loop. All of these arms and loops are important in protein translation. Here, we aimed to delineate the genomic architecture of these arms and loops in cyanobacterial tRNA. Studies from tRNA sequences from 61 cyanobacterial species showed that, except for few tRNAs (tRNA^Asn, tRNA^Leu, tRNA^Gln, and tRNA^Met), all contained a G nucleotide at the 1st position in the acceptor arm. tRNA^Leu and tRNA^Met did not contain any conserved nucleotides at the 1st position whereas tRNA^Asn and tRNA^Gln contained a conserved U¹ nucleotide. In several tRNA families, the variable region also contained conserved nucleotides. Except for tRNA^Met and tRNA^Glu, all other tRNAs contained a conserved A nucleotide at the 1st position in the D-loop. The Ψ-loop contained a conserved U¹-U²-C³-x-A⁵-x-U⁷ sequence, except for tRNA^Gly, tRNA^Ala, tRNA^Val, tRNA^Phe, tRNA^Thr, and tRNA^Gln in which the U⁷ nucleotide was not conserved. However, in tRNA^Asp, the U⁷ nucleotide was substituted with a C⁷ nucleotide. Additionally, tRNA^Arg, tRNA^Gly, and tRNA^Lys of cyanobacteria contained a group I intron within the anti-codon loop region. Maximum composite likelihood study on the transition/transversion of cyanobacterial tRNA revealed that the rate of transition was higher than the rate of transversion. An evolutionary tree was constructed to understand the evolution of cyanobacterial tRNA and analyses revealed that cyanobacterial tRNA may have evolved polyphyletically with high rate of gene loss.

The comparative landscape of duplications in Heliconius melpomene and Heliconius cydno

Gene duplications can facilitate adaptation and may lead to interpopulation divergence, causing reproductive isolation. We used whole-genome resequencing data from 34 butterflies to detect duplications in two Heliconius species, Heliconius cydno and Heliconius melpomene. Taking advantage of three distinctive signals of duplication in short-read sequencing data, we identified 744 duplicated loci in H. cydno and H. melpomene and evaluated the accuracy of our approach using single-molecule sequencing. We have found that duplications overlap genes significantly less than expected at random in H. melpomene, consistent with the action of background selection against duplicates in functional regions of the genome. Duplicate loci that are highly differentiated between H. melpomene and H. cydno map to four different chromosomes. Four duplications were identified with a strong signal of divergent selection, including an odorant binding protein and another in close proximity with a known wing colour pattern locus that differs between the two species.

Comparative genomics reveals adaptive evolution of Asian tapeworm in switching to a new intermediate host

Taenia saginata, Taenia solium and Taenia asiatica (beef, pork and Asian tapeworms, respectively) are parasitic flatworms of major public health and food safety importance. Among them, T. asiatica is a newly recognized species that split from T. saginata via an intermediate host switch ∼1.14 Myr ago. Here we report the 169- and 168-Mb draft genomes of T. saginata and T. asiatica. Comparative analysis reveals that high rates of gene duplications and functional diversifications might have partially driven the divergence between T. asiatica and T. saginata. We observe accelerated evolutionary rates, adaptive evolutions in homeostasis regulation, tegument maintenance and lipid uptakes, and differential/specialized gene family expansions in T. asiatica that may favour its hepatotropism in the new intermediate host. We also identify potential targets for developing diagnostic or intervention tools against human tapeworms. These data provide new insights into the evolution of Taenia parasites, particularly the recent speciation of T. asiatica.

Only one of the three Taenia species causing taeniasis in humans was previously sequenced. Here the authors provide draft genomes of Taenia saginata and Taenia asiatica, analyse genome evolution of all three species, and identify potential targets for developing diagnostic markers or intervention tools.

Neighboring genes for DNA-binding proteins rescue male sterility in Drosophila hybrids

Crosses between closely related animal species often result in male hybrids that are sterile, and the molecular and functional basis of genetic factors for hybrid male sterility is of great interest. Here, we report a molecular and functional analysis of HMS1, a region of 9.2 kb in chromosome 3 of Drosophila mauritiana, which results in virtually complete hybrid male sterility when homozygous in the genetic background of sibling species Drosophila simulans. The HMS1 region contains two strong candidate genes for the genetic incompatibility, agt and Taf1 Both encode unrelated DNA-binding proteins, agt for an alkyl-cysteine-S-alkyltransferase and Taf1 for a subunit of transcription factor TFIID that serves as a multifunctional transcriptional regulator. The contribution of each gene to hybrid male sterility was assessed by means of germ-line transformation, with constructs containing complete agt and Taf1 genomic sequences as well as various chimeric constructs. Both agt and Taf1 contribute about equally to HMS1 hybrid male sterility. Transgenes containing either locus rescue sterility in about one-half of the males, and among fertile males the number of offspring is in the normal range. This finding suggests compensatory proliferation of the rescued, nondysfunctional germ cells. Results with chimeric transgenes imply that the hybrid incompatibilities result from interactions among nucleotide differences residing along both agt and Taf1 Our results challenge a number of preliminary generalizations about the molecular and functional basis of hybrid male sterility, and strongly reinforce the role of DNA-binding proteins as a class of genes contributing to the maintenance of postzygotic reproductive isolation.

A Proteome-wide Fission Yeast Interactome Reveals Network Evolution Principles from Yeasts to Human

Here, we present FissionNet, a proteome-wide binary protein interactome for S. pombe, comprising 2,278 high-quality interactions, of which ∼ 50% were previously not reported in any species. FissionNet unravels previously unreported interactions implicated in processes such as gene silencing and pre-mRNA splicing. We developed a rigorous network comparison framework that accounts for assay sensitivity and specificity, revealing extensive species-specific network rewiring between fission yeast, budding yeast, and human. Surprisingly, although genes are better conserved between the yeasts, S. pombe interactions are significantly better conserved in human than in S. cerevisiae. Our framework also reveals that different modes of gene duplication influence the extent to which paralogous proteins are functionally repurposed. Finally, cross-species interactome mapping demonstrates that coevolution of interacting proteins is remarkably prevalent, a result with important implications for studying human disease in model organisms. Overall, FissionNet is a valuable resource for understanding protein functions and their evolution.

Spermatogenesis Drives Rapid Gene Creation and Masculinization of the X Chromosome in Stalk-Eyed Flies (Diopsidae)

Throughout their evolutionary history, genomes acquire new genetic material that facilitates phenotypic innovation and diversification. Developmental processes associated with reproduction are particularly likely to involve novel genes. Abundant gene creation impacts the evolution of chromosomal gene content and general regulatory mechanisms such as dosage compensation. Numerous studies in model organisms have found complex and, at times contradictory, relationships among these genomic attributes highlighting the need to examine these patterns in other systems characterized by abundant sexual selection. Therefore, we examined the association among novel gene creation, tissue-specific gene expression, and chromosomal gene content within stalk-eyed flies. Flies in this family are characterized by strong sexual selection and the presence of a newly evolved X chromosome. We generated RNA-seq transcriptome data from the testes for three species within the family and from seven additional tissues in the highly dimorphic species, Teleopsis dalmanni. Analysis of dipteran gene orthology reveals dramatic testes-specific gene creation in stalk-eyed flies, involving numerous gene families that are highly conserved in other insect groups. Identification of X-linked genes for the three species indicates that the X chromosome arose prior to the diversification of the family. The most striking feature of this X chromosome is that it is highly masculinized, containing nearly twice as many testes-specific genes as expected based on its size. All the major processes that may drive differential sex chromosome gene content—creation of genes with male-specific expression, development of male-specific expression from pre-existing genes, and movement of genes with male-specific expression—are elevated on the X chromosome of T. dalmanni. This masculinization occurs despite evidence that testes expressed genes do not achieve the same levels of gene expression on the X chromosome as they do on the autosomes.

Genome-Wide Identification of Calcium Dependent Protein Kinase Gene Family in Plant Lineage Shows Presence of Novel D-x-D and D-E-L Motifs in EF-Hand Domain

Calcium ions are considered ubiquitous second messengers in eukaryotic signal transduction pathways. Intracellular Ca(2+) concentration are modulated by various signals such as hormones and biotic and abiotic stresses. Modulation of Ca(2+) ion leads to stimulation of calcium dependent protein kinase genes (CPKs), which results in regulation of gene expression and therefore mediates plant growth and development as well as biotic and abiotic stresses. Here, we reported the CPK gene family of 40 different plant species (950 CPK genes) and provided a unified nomenclature system for all of them. In addition, we analyzed their genomic, biochemical and structural conserved features. Multiple sequence alignment revealed that the kinase domain, auto-inhibitory domain and EF-hands regions of regulatory domains are highly conserved in nature. Additionally, the EF-hand domains of higher plants were found to contain four D-x-D and two D-E-L motifs, while lower eukaryotic plants had two D-x-D and one D-x-E motifs in their EF-hands. Phylogenetic analysis showed that CPK genes are clustered into four different groups. By studying the CPK gene family across the plant lineage, we provide the first evidence of the presence of D-x-D motif in the calcium binding EF-hand domain of CPK proteins.

The genomic architecture of population divergence between subspecies of the European rabbit

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

Primate genome architecture influences structural variation mechanisms and functional consequences

Although nucleotide resolution maps of genomic structural variants (SVs) have provided insights into the origin and impact of phenotypic diversity in humans, comparable maps in nonhuman primates have thus far been lacking. Using massively parallel DNA sequencing, we constructed fine-resolution genomic structural variation maps in five chimpanzees, five orang-utans, and five rhesus macaques. The SV maps, which are comprised of thousands of deletions, duplications, and mobile element insertions, revealed a high activity of retrotransposition in macaques compared with great apes. By comparison, nonallelic homologous recombination is specifically active in the great apes, which is correlated with architectural differences between the genomes of great apes and macaque. Transcriptome analyses across nonhuman primates and humans revealed effects of species-specific whole-gene duplication on gene expression. We identified 13 gene duplications coinciding with the species-specific gain of tissue-specific gene expression in keeping with a role of gene duplication in the promotion of diversification and the acquisition of unique functions. Differences in the present day activity of SV formation mechanisms that our study revealed may contribute to ongoing diversification and adaptation of great ape and Old World monkey lineages.

Chromosome inversions, genomic differentiation and speciation in the African malaria mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.

Novel genes from formation to function

The study of the evolution of novel genes generally focuses on the formation of new coding sequences. However, equally important in the evolution of novel functional genes are the formation of regulatory regions that allow the expression of the genes and the effects of the new genes in the organism as well. Herein, we discuss the current knowledge on the evolution of novel functional genes, and we examine in more detail the youngest genes discovered. We examine the existing data on a very recent and rapidly evolving cluster of duplicated genes, the Sdic gene cluster. This cluster of genes is an excellent model for the evolution of novel genes, as it is very recent and may still be in the process of evolving.

Simulated evolution of protein-protein interaction networks with realistic topology

We model the evolution of eukaryotic protein-protein interaction (PPI) networks. In our model, PPI networks evolve by two known biological mechanisms: (1) Gene duplication, which is followed by rapid diversification of duplicate interactions. (2) Neofunctionalization, in which a mutation leads to a new interaction with some other protein. Since many interactions are due to simple surface compatibility, we hypothesize there is an increased likelihood of interacting with other proteins in the target protein’s neighborhood. We find good agreement of the model on 10 different network properties compared to high-confidence experimental PPI networks in yeast, fruit flies, and humans. Key findings are: (1) PPI networks evolve modular structures, with no need to invoke particular selection pressures. (2) Proteins in cells have on average about 6 degrees of separation, similar to some social networks, such as human-communication and actor networks. (3) Unlike social networks, which have a shrinking diameter (degree of maximum separation) over time, PPI networks are predicted to grow in diameter. (4) The model indicates that evolutionarily old proteins should have higher connectivities and be more centrally embedded in their networks. This suggests a way in which present-day proteomics data could provide insights into biological evolution.

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.

How can satellite DNA divergence cause reproductive isolation? Let us count the chromosomal ways

Satellites are one of the most enigmatic parts of the eukaryotic genome. These highly repetitive, noncoding sequences make up as much as half or more of the genomic content and are known to play essential roles in chromosome segregation during meiosis and mitosis, yet they evolve rapidly between closely related species. Research over the last several decades has revealed that satellite divergence can serve as a formidable reproductive barrier between sibling species. Here we highlight several key studies on Drosophila and other model organisms demonstrating deleterious effects of satellites and their rapid evolution on the structure and function of chromosomes in interspecies hybrids. These studies demonstrate that satellites can impact chromosomes at a number of different developmental stages and through distinct cellular mechanisms, including heterochromatin formation. These findings have important implications for how loci that cause postzygotic reproductive isolation are viewed.

Functional evidence that a recently evolved Drosophila sperm-specific gene boosts sperm competition

In many species, both morphological and molecular traits related to sex and reproduction evolve faster in males than in females. Ultimately, rapid male evolution relies on the acquisition of genetic variation associated with differential reproductive success. Many newly evolved genes are associated with novel functions that might enhance male fitness. However, functional evidence of the adaptive role of recently originated genes in males is still lacking. The Sperm dynein intermediate chain multigene family, which encodes a Sperm dynein intermediate chain presumably involved in sperm motility, originated from complex genetic rearrangements in the lineage that leads to Drosophila melanogaster within the last 5.4 million years since its split from Drosophila simulans. We deleted all the members of this multigene family resident on the X chromosome of D. melanogaster by chromosome engineering and found that, although the deletion does not result in a reduction of progeny number, it impairs the competence of the sperm in the presence of sperm from wild-type males. Therefore, the Sperm dynein intermediate chain multigene family contributes to the differential reproductive success among males and illustrates precisely how quickly a new gene function can be incorporated into the genetic network of a species.

Chromatin evolution and molecular drive in speciation

Are there biological generalities that underlie hybrid sterility or inviability? Recently, around a dozen "speciation genes" have been identified mainly in Drosophila, and the biological functions of these genes are revealing molecular generalities. Major cases of hybrid sterility and inviability seem to result from chromatin evolution and molecular drive in speciation. Repetitive satellite DNAs within heterochromatin, especially at centromeres, evolve rapidly through molecular drive mechanisms (both meiotic and centromeric). Chromatin-binding proteins, therefore, must also evolve rapidly to maintain binding capability. As a result, chromatin binding proteins may not be able to interact with chromosomes from another species in a hybrid, causing hybrid sterility and inviability.

Gene duplication and ectopic gene conversion in Drosophila

The evolutionary impact of gene duplication events has been a theme of Drosophila genetics dating back to the Morgan School. While considerable attention has been placed on the genetic novelties that duplicates are capable of introducing, and the role that positive selection plays in their early stages of duplicate evolution, much less attention has been given to the potential consequences of ectopic (non-allelic) gene conversion on these evolutionary processes. In this paper we consider the historical origins of ectopic gene conversion models and present a synthesis of the current Drosophila data in light of several primary questions in the field.

Novel acid phosphatase in Candida glabrata suggests selective pressure and niche specialization in the phosphate signal transduction pathway

Evolution through natural selection suggests unnecessary genes are lost. We observed that the yeast Candida glabrata lost the gene encoding a phosphate-repressible acid phosphatase (PHO5) present in many yeasts including Saccharomyces cerevisiae. However, C. glabrata still had phosphate starvation-inducible phosphatase activity. Screening a C. glabrata genomic library, we identified CgPMU2, a member of a three-gene family that contains a phosphomutase-like domain. This small-scale gene duplication event could allow for sub- or neofunctionalization. On the basis of phylogenetic and biochemical characterizations, CgPMU2 has neofunctionalized to become a broad range, phosphate starvation-regulated acid phosphatase, which functionally replaces PHO5 in this pathogenic yeast. We determined that CgPmu2, unlike ScPho5, is not able to hydrolyze phytic acid (inositol hexakisphosphate). Phytic acid is present in fruits and seeds where S. cerevisiae grows, but is not abundant in mammalian tissues where C. glabrata grows. We demonstrated that C. glabrata is limited from an environment where phytic acid is the only source of phosphate. Our work suggests that during evolutionary time, the selection for the ancestral PHO5 was lost and that C. glabrata neofunctionalized a weak phosphatase to replace PHO5. Convergent evolution of a phosphate starvation-inducible acid phosphatase in C. glabrata relative to most yeast species provides an example of how small changes in signal transduction pathways can mediate genetic isolation and uncovers a potential speciation gene.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization

Background

Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates.

Results

Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species.

Conclusions

Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Altered heterochromatin binding by a hybrid sterility protein in Drosophila sibling species

Hybrid sterility of the heterogametic sex is one of the first postzygotic reproductive barriers to evolve during speciation, yet the molecular basis of hybrid sterility is poorly understood. We show that the hybrid male sterility gene Odysseus-site homeobox (OdsH) encodes a protein that localizes to evolutionarily dynamic loci within heterochromatin and leads to their decondensation. In Drosophila mauritiana x Drosophila simulans male hybrids, OdsH from D. mauritiana (OdsHmau) acts as a sterilizing factor by associating with the heterochromatic Y chromosome of D. simulans, whereas D. simulans OdsH (OdsHsim) does not. Characterization of sterile hybrid testes revealed that OdsH abundance and localization in the premeiotic phases of spermatogenesis differ between species. These results reveal that rapid heterochromatin evolution affects the onset of hybrid sterility.

Accelerated evolution of the Prdm9 speciation gene across diverse metazoan taxa

The onset of prezygotic and postzygotic barriers to gene flow between populations is a hallmark of speciation. One of the earliest postzygotic isolating barriers to arise between incipient species is the sterility of the heterogametic sex in interspecies' hybrids. Four genes that underlie hybrid sterility have been identified in animals: Odysseus, JYalpha, and Overdrive in Drosophila and Prdm9 (Meisetz) in mice. Mouse Prdm9 encodes a protein with a KRAB motif, a histone methyltransferase domain and several zinc fingers. The difference of a single zinc finger distinguishes Prdm9 alleles that cause hybrid sterility from those that do not. We find that concerted evolution and positive selection have rapidly altered the number and sequence of Prdm9 zinc fingers across 13 rodent genomes. The patterns of positive selection in Prdm9 zinc fingers imply that rapid evolution has acted on the interface between the Prdm9 protein and the DNA sequences to which it binds. Similar patterns are apparent for Prdm9 zinc fingers for diverse metazoans, including primates. Indeed, allelic variation at the DNA–binding positions of human PRDM9 zinc fingers show significant association with decreased risk of infertility. Prdm9 thus plays a role in determining male sterility both between species (mouse) and within species (human). The recurrent episodes of positive selection acting on Prdm9 suggest that the DNA sequences to which it binds must also be evolving rapidly. Our findings do not identify the nature of the underlying DNA sequences, but argue against the proposed role of Prdm9 as an essential transcription factor in mouse meiosis. We propose a hypothetical model in which incompatibilities between Prdm9-binding specificity and satellite DNAs provide the molecular basis for Prdm9-mediated hybrid sterility. We suggest that Prdm9 should be investigated as a candidate gene in other instances of hybrid sterility in metazoans.

Speciation, the process by which one species splits into two, involves reproductive barriers between previously interbreeding populations. The question of how speciation occurs has rightly occupied the attention of biologists since before Darwin's “On the Origin of Species.” Studies of recently diverged species have revealed the presence of hybrid sterility genes (colloquially referred to as “speciation genes”), alleles of which are associated with sterility of interspecies hybrids. Mouse Prdm9 is the only known such gene in vertebrate animals. Here we report that the Prdm9 protein has evolved extremely rapidly in its DNA-binding domain, comprising an array of “zinc fingers.” This suggests that hybrid sterility may arise from a mismatch between the DNA-binding specificity of Prdm9 and rapidly evolving DNA. We propose that Prdm9 binds to satellite-DNA repeats evolving rapidly within and between different species. Prdm9 evolution is unusual because other hybrid sterility genes appear only to evolve rapidly in isolated bursts, whereas Prdm9 has evolved rapidly over 700 million years, in many rodent species, diverse primates and other metazoans. This leads to the tantalizing possibility that Prdm9 may have served as a “speciation gene” on other occasions in metazoan evolution, a possibility that will now need to be investigated.

Gene duplication, transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.

Independent losses of function in a polyphenol oxidase in rice: differentiation in grain discoloration between subspecies and the role of positive selection under domestication

Asian rice (Oryza sativa) cultivars originated from wild rice and can be divided into two subspecies by several criteria, one of which is the phenol reaction (PHR) phenotype. Grains of indica cultivars turn brown in a phenol solution that accelerates a similar process that occurs during prolonged storage. By contrast, the grains of japonica do not discolor. This distinction may reflect the divergent domestication of these two subspecies. The PHR is controlled by a single gene, Phr1; here, we report the cloning of Phr1, which encodes a polyphenol oxidase. The Phr1 gene is indeed responsible for the PHR phenotype, as transformation with a functional Phr1 can complement a PHR negative cultivar. Phr1 is defective in all japonica lines but functional in nearly all indica and wild strains. Phylogenetic analysis showed that the defects in Phr1 arose independently three times. The multiple recent origins and rapid spread of phr1 in japonica suggest the action of positive selection, which is further supported by several population genetic tests. This case may hence represent an example of artificial selection driving the differentiation among domesticated varieties.

Gene family evolution across 12 Drosophila genomes

Comparison of whole genomes has revealed large and frequent changes in the size of gene families. These changes occur because of high rates of both gene gain (via duplication) and loss (via deletion or pseudogenization), as well as the evolution of entirely new genes. Here we use the genomes of 12 fully sequenced Drosophila species to study the gain and loss of genes at unprecedented resolution. We find large numbers of both gains and losses, with over 40% of all gene families differing in size among the Drosophila. Approximately 17 genes are estimated to be duplicated and fixed in a genome every million years, a rate on par with that previously found in both yeast and mammals. We find many instances of extreme expansions or contractions in the size of gene families, including the expansion of several sex- and spermatogenesis-related families in D. melanogaster that also evolve under positive selection at the nucleotide level. Newly evolved gene families in our dataset are associated with a class of testes-expressed genes known to have evolved de novo in a number of cases. Gene family comparisons also allow us to identify a number of annotated D. melanogaster genes that are unlikely to encode functional proteins, as well as to identify dozens of previously unannotated D. melanogaster genes with conserved homologs in the other Drosophila. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this genomic revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

Though comparative genome sequencing has revealed vast similarities in the total number of genes contained within closely related species, this similarity hides enormous complexities in the identity and number of constituent proteins. Species can differ in their complement of genes through both gene duplication and loss. Here we investigated the gain and loss of genes from the genomes of 12 fully sequenced Drosophila (fruit flies). We find high rates of gain and loss in all species and estimate that approximately one new gene is gained or lost every 60,000 years. We also find several hundred cases of extremely rapid gene turnover, with dozens of genes gained or lost in only a few million years. The highest turnover in gene number occurs in genes involved in sex and reproduction. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this evolutionary revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

Drosophila biology in the genomic age

Over the course of the past century, flies in the family Drosophilidae have been important models for understanding genetic, developmental, cellular, ecological, and evolutionary processes. Full genome sequences from a total of 12 species promise to extend this work by facilitating comparative studies of gene expression, of molecules such as proteins, of developmental mechanisms, and of ecological adaptation. Here we review basic biological and ecological information of the species whose genomes have recently been completely sequenced in the context of current research.