Negligible nuclear introgression despite complete mitochondrial capture between two species of chipmunks

The idea that species boundaries can be semipermeable to gene flow is now widely accepted but the evolutionary importance of introgressive hybridization remains unclear. Here we examine the genomic contribution of gene flow between two hybridizing chipmunk species, Tamias ruficaudus and T. amoenus. Previous studies have shown that ancient hybridization has resulted in complete fixation of introgressed T. ruficaudus mitochondrial DNA (mtDNA) in some populations of T. amoenus, but the extent of nuclear introgression is not known. We used targeted capture to sequence over 10,500 gene regions from multiple individuals of both species. We found that most of the nuclear genome is sorted between these species and that overall genealogical patterns do not show evidence for introgression. Our analysis rules out all but very minor levels of interspecific gene flow, indicating that introgressive hybridization has had little impact on the overall genetic composition of these species outside of the mitochondrial genome. Given that much of the evidence for introgression in animals has come from mtDNA, our results underscore that unraveling the importance introgressive hybridization during animal speciation will require a genome‐wide perspective that is still absent for many species.     

Molecular variability in the Celleporella hyalina (Bryozoa; Cheilostomata) species complex: evidence for cryptic speciation from complete mitochondrial genomes

The bryozoan Celleporella has been shown to be composed of multiple, often cryptic, lineages. We sequenced two complete mitochondrial (mt) genomes of the Celleporella hyalina species complex from Wales, UK and Norway (i) to determine genetic divergence at the complete mt genome level, and (ii) to design new molecular markers for examining the interrelationships amongst the major lineages. In addressing (i), we estimated genetic divergence at three levels: (a) nucleotide diversity (π), (b) genome size, and (c) gene order. Genes nad4L, nad6, and atp8 showed the highest levels of divergence, and rrnL, rrnS, and cox1 showed the lowest levels. Inter-genome nucleotide divergence of protein-coding and ribosomal RNA genes, measured as π, was 0.21. The two genomes differed substantially in size, with the Norwegian genome being 2,573 base pairs (bp) longer than the Welsh genome, 17,265 and 14,692 bp, respectively. This difference in size is attributable to long non-coding regions present in the Norwegian genome. Both genomes exhibit similar gene orders, except for the translocation of one transfer RNA (trnA). Considering the high nucleotide diversity, genome size difference and change in gene order, these mt genomes are considered sufficiently divergent to have originated from two distinct species. In addressing (ii) we designed PCR primers that flank the most conserved regions of the genome: 1,300 bp of cox1 and a contiguous 2,000 bp fragment of rrnL + rrnS. The primers have yielded products for tissue from Wales, Norway, New Zealand, Alaska and Chile and should provide useful tools in establishing species- and population-level diversity within the Celleporella complex.     

Contrasting rates of evolution in Pm3 loci from three wheat species and rice

The Pm3 gene from wheat confers resistance against powdery mildew and recent studies have shown that it is a member of a multigene family in the wheat genome. We compared genomic sequences ranging from 178 to 332 kb containing six Pm3-like genes and five gene fragments from orthologous loci in the A genome of wheat at three different ploidy levels. We found that the wheat Pm3 loci display an extremely dynamic evolution where sequence conservation is minimal between species and basically limited to very short sequences containing the genetic markers that define the orthology. The Pm3-like genes and their up- and downstream regions were reshuffled by multiple rearrangements, resulting in a complex mosaic of conserved and unique sequences. Comparison with rice showed that the known wheat Pm3-like genes represent only one branch of a large superfamily with several clusters in rice and suggests the presence of additional similar genes in the wheat genome. Estimates of divergence times and transposable-element insertions indicate that the Pm3 locus in wheat has undergone more drastic changes in its recent evolution than its counterpart in rice. This indicates that loci containing homologous resistance gene analogs can evolve at highly variable speeds in different species.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Exploring genome architecture through GOV: a WWW-based gene order visualizer

MOTIVATION: The past decade has seen extension in the methods of sequence analysis from single gene based to analyzing multiple genes and proteins simultaneously. Consequently, there is a need for software tools that will allow mining of these enormous datasets at genome level effectively. A key challenge is to make them user-friendly, available to a larger community and integrate with public domain software without much hassle. RESULTS: A web-based interactive computational tool is described for visualization and comparison of gene order from prokaryotic and selected viral genome data. Many intriguing similarities and differences in gene order of multiple genomes can be compared and revealed. The interface facilitates easy extraction of the nucleotide sequence of the gene of interest and BLAST analysis against GenBank at NCBI to provide insights into gene functions and orthologs of the gene in other species.     

Rapid evolution of a sexual reproduction gene in centric diatoms of the genus Thalassiosira

Sexual reproduction is commonly assumed to occur in the vast majority of diatoms due to the intimate association of this process with cell size control. Surprisingly, however, little is known about the impact of sexual events on diatom population dynamics. The Sig1 gene is strongly upregulated during sexual reproduction in the centric diatom Thalassiosira weissflogii and has been hypothesized to encode a protein involved in gamete recognition. In the present study, degenerate PCR primers were designed and used to amplify a portion of Sig1 from three closely related species in the cosmopolitan genus Thalassiosira, Thalassiosira oceanica, Thalassiosira guillardii, and Thalassiosira pseudonana. Identification of Sig1 in these three additional species facilitated development of this gene as a molecular marker for diatom sexual events. Examination of the new sequences indicated that multiple copies of Sig1 are probably present in the genome. Moreover, compared to the housekeeping gene beta-tubulin, the Sig1 genes of isolates of T. weissflogii collected from different regions of the Atlantic and Pacific oceans displayed high levels of divergence. The Sig1 genes of the four closely related Thalassiosira species also displayed high levels of sequence divergence compared to the levels observed with a second gene, Fcp, probably explaining why Sig1 could not be amplified from more distantly related species. The high levels of sequence divergence both within and between species suggest that Sig1 is rapidly evolving in a manner reminiscent of the manner observed in other genes that encode gamete recognition proteins. A simple model is presented for Sig1 evolution and the implications of such a rapidly evolving sexual reproduction gene for diatom speciation and population dynamics.     

Prokaryotic genomes: the emerging paradigm of genome-based microbiology

Comparative analysis of the complete sequences of seven bacterial and three archaeal genomes leads to the first generalizations of emerging genome-based microbiology. Protein sequences are, generally, highly conserved, with ∼70% of the gene products in bacteria and archaea containing ancient conserved regions. In contrast, there is little conservation of genome organization, except for a few essential operons. The most striking conclusions derived by comparison of multiple genomes from phylogenetically distant species are that the number of universally conserved gene families is very small and that multiple events of horizontal gene transfer and genome fusion are major forces in evolution.     

Phylogeny and molecular evolution of the DMC1 gene within the StH genome species in Triticeae (Poaceae)

To estimate the phylogeny and molecular evolution of a single-copy nuclear disrupted meiotic cDNA (DMC1) gene within the StH genome species, two DMC1 homoeologous sequences were isolated from nearly all the sampled StH genome species and were analyzed with those from seven diploid taxa representing the St and H genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) there is a close relationship among North American StH genome species; (2) the DMC1 gene sequences of the StH genome species from North America and Eurasia are evolutionarily distinct; (3) the StH genome polyploids have higher levels of sequence diversity in the St genome homoeolog than the H genome homoeolog; (4) the DMC1 sequence may evolve faster in the polyploid species than in the diploids; (5) high dN and dN/dS values in the St genome within polyploid species could be caused by low selective constraints or AT-biased mutation pressure. Our result provides some insight on evolutionary dynamics of duplicate DMC1 gene, the polyploidization events and phylogeny of the StH genome species.     

Variation and divergence of multilocus genome markers in the species of the genus Chironomus (Diptera, Chironomidae)

Variation and divergence patterns of the multilocus genome markers in twelve Chironomus species belonging to the plumosus and piger sibling-species groups were examined by use of polymerase chain reaction with random primers (RAPD method). The chironomid species showed high levels of the RAPD markers polymorphism. The genetic distances (GD) were determined between the species within the group of closely related species, as well as between the species from different groups. The estimates obtained characterized the divergence levels between the sibling species (GD = 0.739) and morphologically distinct species (GD = 0.935). Comparison of the variation and divergence levels of the RAPD markers with those for the other genome markers, namely, the enzyme-coding genes and chromosomes (gene linkage groups) have demonstrated different divergence rates for different genome components during speciation of Chironimids.     

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks ( Tamias ruficaudus and T. amoenus ). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus , T. a. canicaudus , with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus . Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.     

Phylogeny and molecular evolution of the Acc1 gene within the StH genome species in Triticeae (Poaceae)

To estimate the phylogeny and molecular evolution of a single-copy gene encoding plastid acetyl-CoA carboxylase (Acc1) within the StH genome species, two Acc1 homoeologous sequences were isolated from nearly all the sampled StH genome species and were analyzed with those from 35 diploid taxa representing 19 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) the StH genome species from the same areas or neighboring geographic regions are closely related to each other; (2) the Acc1 gene sequences of the StH genome species from North America and Eurasia are evolutionarily distinct; (3) Dasypyrum has contributed to the nuclear genome of Elymus repens and Elymus mutabilis; (4) the StH genome polyploids have higher levels of sequence diversity in the H genome homoeolog than the St genome homoeolog; and (5) the Acc1 sequence may evolve faster in the polyploid species than in the diploids. Our result provides some insight on evolutionary dynamics of duplicate Acc1 gene, the polyploidy speciation and phylogeny of the StH genome species.     

Cryptic sex? Estimates of genome exchange in unisexual mole salamanders (Ambystoma sp.)

Cryptic sex has been argued to explain the exceptional longevity of certain parthenogenetic vertebrate lineages, yet direct measurements of genetic exchange between sexual and apparently parthenogenetic forms are rare. Female unisexual mole salamanders (Ambystoma sp.) are the oldest known unisexual vertebrate lineage (~5 million years), and one hypothesis for their persistence is that allopolyploid female unisexuals periodically exchange haploid genomes ‘genome exchange’ during gynogenetic reproduction with males from sympatric sexual species. We test this hypothesis by using genome‐specific microsatellite DNA markers to estimate the rates of genome exchange between sexual males and unisexual females in two ponds in NE Ohio. We also test the prediction that levels of gene flow should be higher for ‘sympatric’ (sexual males present) genomes in unisexuals compared to ‘allopatric’ (sexual males absent) unisexual genomes. We used a model testing framework in the coalescent‐based program MIGRATE‐N to compare models where unidirectional gene flow is present and absent between sexual species and unisexuals. As predicted, our results show higher levels of gene flow between sexuals and sympatric unisexual genomes compared to lower (likely artefactual) levels of gene flow between sexuals and allopatric unisexual genomes. Our results provide direct evidence that genome exchange between sexual and unisexual Ambystoma occurs and demonstrate that the magnitude depends on which sexual species are present. The relatively high levels of gene flow suggest that unisexuals must be at a selective advantage over sexual forms so as to avoid extinction due to genetic swamping through genome exchange.     

Comparative genomics of duplicate γ-glutamyl transferase genes in teleosts: medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), green spotted pufferfish (Tetraodon nigroviridis), fugu (Takifugu rubripes), and zebrafish (Danio rerio)

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we examine multiple paralogous genes of γ-glutamyl transferase (GGT) in several distantly related teleost species including medaka, stickleback, green spotted pufferfish, fugu, and zebrafish. Through mining genome databases, we have identified multiple GGT orthologs. Duplicate (paralogous) GGT sequences for GGT1 (GGT1 a and b), GGTL1 (GGTL1 a and b), and GGTL3 (GGTL3 a and b) were identified for each species. Phylogenetic analysis suggests that GGTs are ancient proteins conserved across most metazoan phyla and those paralogous GGTs in teleosts likely arose from the serial 3R genome duplication events. A third GGTL1 gene (GGTL1c) was found in green spotted pufferfish; however, this gene is not present in medaka, stickleback, or fugu. Similarly, one or both paralogs of GGTL3 appear to have been lost in green spotted pufferfish, fugu, and zebrafish. Syntenic relationships were highly maintained between duplicated teleost chromosomes, among teleosts and across ray-finned (Actinopterygii) and lobe-finned (Sarcopterygii) species. To assess subfunction partitioning, six medaka GGT genes were cloned and assessed for developmental and tissue-specific expression. On the basis of these data, we propose a modification of the "duplication-degeneration-complementation" model of subfunction partitioning where quantitative differences rather than absolute differences in gene expression are observed between gene paralogs. Our results demonstrate that multiple GGT genes have been retained within teleost genomes. Questions remain, however, regarding the functional roles of multiple GGTs in these species.     

PGMapper: a web-based tool linking phenotype to genes

SUMMARY: With the availability of whole genome sequence in many species, linkage analysis, positional cloning and microarray are gradually becoming powerful tools for investigating the links between phenotype and genotype or genes. However, in these methods, causative genes underlying a quantitative trait locus, or a disease, are usually located within a large genomic region or a large set of genes. Examining the function of every gene is very time consuming and needs to retrieve and integrate the information from multiple databases or genome resources. PGMapper is a software tool for automatically matching phenotype to genes from a defined genome region or a group of given genes by combining the mapping information from the Ensembl database and gene function information from the OMIM and PubMed databases. PGMapper is currently available for candidate gene search of human, mouse, rat, zebrafish and 12 other species. AVAILABILITY: Available online at http://www.genediscovery.org/pgmapper/index.jsp.     

The Plasmodium Apicoplast Genome: Conserved Structure and Close Relationship of P. ovale to Rodent Malaria Parasites

Apicoplast, a nonphotosynthetic plastid derived from secondary symbiotic origin, is essential for the survival of malaria parasites of the genus Plasmodium. Elucidation of the evolution of the apicoplast genome in Plasmodium species is important to better understand the functions of the organelle. However, the complete apicoplast genome is available for only the most virulent human malaria parasite, Plasmodium falciparum. Here, we obtained the near-complete apicoplast genome sequences from eight Plasmodium species that infect a wide variety of vertebrate hosts and performed structural and phylogenetic analyses. We found that gene repertoire, gene arrangement, and other structural attributes were highly conserved. Phylogenetic reconstruction using 30 protein-coding genes of the apicoplast genome inferred, for the first time, a close relationship between P. ovale and rodent parasites. This close relatedness was robustly supported using multiple evolutionary assumptions and models. The finding suggests that an ancestral host switch occurred between rodent and human Plasmodium parasites.     

Evolutionary genomics of the Fox genes: Origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.