Search PubMed⌕ Search

Biomedical subjects

Antonis Rokas

Publications and source records attributed to Antonis Rokas.

17 recordsLinked to original sources
PubMed2025

Discovery of additional ancient genome duplications in yeasts.

Whole genome duplication (WGD) has had profound macroevolutionary impacts on diverse lineages1,2, preceding adaptive radiations in vertebrates3-5, teleost fish6,7, and angiosperms8,9. In contrast to the many known ancient WGDs in animals10,11 and especially plants12-14, we are aware of evidence for only four in fungi15,16. The oldest of these occurred ~100 million years ago (mya) and is shared by ~60 extant Saccharomycetales species17,18, including the baker's yeast Saccharomyces cerevisiae (Fig. 1). Notably, this is the only known ancient WGD in the yeast subphylum Saccharomycotina. The dearth of ancient WGD events in fungi remains a mystery15. Some studies have suggested that fungal lineages that experience chromosome19 and genome15 duplication quickly go extinct, leaving no trace in the genomic record, while others contend that the lack of known WGD is due to an absence of data15,16. Under the second hypothesis, additional sampling and deeper sequencing of fungal genomes should lead to the discovery of more WGD events. Coupling hundreds of recently published genomes from nearly every described Saccharomycotina species with three additional long-read assemblies, we discovered three novel WGD events. While the functions of retained duplicate genes originating from these events are broad, they bear many similarities to the well-known WGD that occurred in the Saccharomycetales17. Our results suggest that WGD may be a more common evolutionary force in fungi than previously believed.

Convergent Evolution↗
PubMed2006

Bayesian estimation of concordance among gene trees.

Multigene sequence data have great potential for elucidating important and interesting evolutionary processes, but statistical methods for extracting information from such data remain limited. Although various biological processes may cause different genes to have different genealogical histories (and hence different tree topologies), we also may expect that the number of distinct topologies among a set of genes is relatively small compared with the number of possible topologies. Therefore evidence about the tree topology for one gene should influence our inferences of the tree topology on a different gene, but to what extent? In this paper, we present a new approach for modeling and estimating concordance among a set of gene trees given aligned molecular sequence data. Our approach introduces a one-parameter probability distribution to describe the prior distribution of concordance among gene trees. We describe a novel 2-stage Markov chain Monte Carlo (MCMC) method that first obtains independent Bayesian posterior probability distributions for individual genes using standard methods. These posterior distributions are then used as input for a second MCMC procedure that estimates a posterior distribution of gene-to-tree maps (GTMs). The posterior distribution of GTMs can then be summarized to provide revised posterior probability distributions for each gene (taking account of concordance) and to allow estimation of the proportion of the sampled genes for which any given clade is true (the sample-wide concordance factor). Further, under the assumption that the sampled genes are drawn randomly from a genome of known size, we show how one can obtain an estimate, with credibility intervals, on the proportion of the entire genome for which a clade is true (the genome-wide concordance factor). We demonstrate the method on a set of 106 genes from 8 yeast species.

Algorithms↗
PubMed2006

Regulation of body pigmentation by the Abdominal-B Hox protein and its gain and loss in Drosophila evolution.

Hox genes have been implicated in the evolution of many animal body patterns, but the molecular events underlying trait modification have not been elucidated. Pigmentation of the posterior male abdomen is a recently acquired trait in the Drosophila melanogaster lineage. Here, we show that the Abdominal-B (ABD-B) Hox protein directly activates expression of the yellow pigmentation gene in posterior segments. ABD-B regulation of pigmentation evolved through the gain of ABD-B binding sites in a specific cis-regulatory element of the yellow gene of a common ancestor of sexually dimorphic species. Within the melanogaster species group, male-specific pigmentation has subsequently been lost by at least three different mechanisms, including the mutational inactivation of a key ABD-B binding site in one lineage. These results demonstrate how Hox regulation of traits and target genes is gained and lost at the species level and have general implications for the evolution of body form at higher taxonomic levels.

Animals↗
PubMed2006

Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene.

The independent evolution of morphological similarities is widespread. For simple traits, such as overall body colour, repeated transitions by means of mutations in the same gene may be common. However, for more complex traits, the possible genetic paths may be more numerous; the molecular mechanisms underlying their independent origins and the extent to which they are constrained to follow certain genetic paths are largely unknown. Here we show that a male wing pigmentation pattern involved in courtship display has been gained and lost multiple times in a Drosophila clade. Each of the cases we have analysed (two gains and two losses) involved regulatory changes at the pleiotropic pigmentation gene yellow. Losses involved the parallel inactivation of the same cis-regulatory element (CRE), with changes at a few nucleotides sufficient to account for the functional divergence of one element between two sibling species. Surprisingly, two independent gains of wing spots resulted from the co-option of distinct ancestral CREs. These results demonstrate how the functional diversification of the modular CREs of pleiotropic genes contributes to evolutionary novelty and the independent evolution of morphological similarities.

Animals↗
PubMed2005

Animal evolution and the molecular signature of radiations compressed in time.

The phylogenetic relationships among most metazoan phyla remain uncertain. We obtained large numbers of gene sequences from metazoans, including key understudied taxa. Despite the amount of data and breadth of taxa analyzed, relationships among most metazoan phyla remained unresolved. In contrast, the same genes robustly resolved phylogenetic relationships within a major clade of Fungi of approximately the same age as the Metazoa. The differences in resolution within the two kingdoms suggest that the early history of metazoans was a radiation compressed in time, a finding that is in agreement with paleontological inferences. Furthermore, simulation analyses as well as studies of other radiations in deep time indicate that, given adequate sequence data, the lack of resolution in phylogenetic trees is a signature of closely spaced series of cladogenetic events.

Animals↗
PubMed2005

More genes or more taxa? The relative contribution of gene number and taxon number to phylogenetic accuracy.

The relative contribution of taxon number and gene number to accuracy in phylogenetic inference is a major issue in phylogenetics and of central importance to the choice of experimental strategies for the successful reconstruction of a broad sketch of the tree of life. Maximization of the number of taxa sampled is the strategy favored by most phylogeneticists, although its necessity remains the subject of debate. Vast increases in gene number are now possible due to advances in genomics, but large numbers of genes will be available for only modest numbers of taxa, raising the question of whether such genome-scale phylogenies will be robust to the addition of taxa. To examine the relative benefit of increasing taxon number or gene number to phylogenetic accuracy, we have developed an assay that utilizes the symmetric difference tree distance as a measure of phylogenetic accuracy. We have applied this assay to a genome-scale data matrix containing 106 genes from 14 yeast species. Our results show that increasing taxon number correlates with a slight decrease in phylogenetic accuracy. In contrast, increasing gene number has a significant positive effect on phylogenetic accuracy. Analyses of an additional taxon-rich data matrix from the same yeast clade show that taxon number does not have a significant effect on phylogenetic accuracy. The positive effect of gene number and the lack of effect of taxon number on phylogenetic accuracy are also corroborated by analyses of two data matrices from mammals and angiosperm plants, respectively. We conclude that, for typical data sets, the number of genes utilized may be a more important determinant of phylogenetic accuracy than taxon number.

Evolution, Molecular↗
PubMed2004

Parallel inactivation of multiple GAL pathway genes and ecological diversification in yeasts.

Understanding the evolutionary relationship between genome content and ecological niche is one of the fundamental challenges of biology. The distinct physiologies of yeast species provide a window into how genomes evolve in concert with niche. Although the enzymes of the well studied yeast galactose utilization pathway are present in all domains of life, we have found that multiple genes of the GAL pathway are absent from four yeast species that cannot use galactose. Whereas three species lack any trace of the pathway except a single gene, Saccharomyces kudriavzevii, a close relative of Saccharomyces cerevisiae, retains remnants of all seven dedicated GAL genes as syntenic pseudogenes, providing a rare glimpse of an entire pathway in the process of degeneration. An estimate of the timing of gene inactivation suggests that pathway degeneration began early in the lineage and proceeded rapidly. S. kudriavzevii exhibits several other divergent physiological properties that are associated with a shift in ecological niche. These results suggest that rapid and irreversible gene inactivation and pathway degeneration are associated with adaptation to new ecological niches in natural populations. Inactivated genes may generally serve as markers of specific functions made dispensable by recent adaptive shifts.

Adaptation, Physiological↗
PubMed2003

Genome-scale approaches to resolving incongruence in molecular phylogenies.

One of the most pervasive challenges in molecular phylogenetics is the incongruence between phylogenies obtained using different data sets, such as individual genes. To systematically investigate the degree of incongruence, and potential methods for resolving it, we screened the genome sequences of eight yeast species and selected 106 widely distributed orthologous genes for phylogenetic analyses, singly and by concatenation. Our results suggest that data sets consisting of single or a small number of concatenated genes have a significant probability of supporting conflicting topologies. By contrast, analyses of the entire data set of concatenated genes yielded a single, fully resolved species tree with maximum support. Comparable results were obtained with a concatenation of a minimum of 20 genes; substantially more genes than commonly used but a small fraction of any genome. These results have important implications for resolving branches of the tree of life.

Genes, Fungal↗
PubMed2003

Lifecycle closure, lineage sorting, and hybridization revealed in a phylogenetic analysis of European oak gallwasps (Hymenoptera: Cynipidae: Cynipini) using mitochondrial sequence data.

Oak gallwasps are cyclically parthenogenetic insects that induce a wide diversity of highly complex species- and generation-specific galls on oaks and other Fagaceae. Phylogenetic relationships within oak gallwasps remain to be established, while sexual and parthenogenetic generations of many species remain unpaired. Previous work on oak gallwasps has revealed substantial intra-specific variation, particularly between regions known to represent discrete Pleistocene glacial refuges. Here we use statistical phylogenetic inference methods on sequence data for a fragment of the mitochondrial cytochrome b gene to reconstruct the relationships among 62 oak gallwasp species. For 16 of these we also include 23 additional cytochrome b haplotype sequences from different Pleistocene refuge areas to test the effect of intra-specific variation on inter-specific phylogeny reconstruction. The reconstructed phylogenies show good intra-generic resolution and identify several conserved clades, but fail to reconstruct either very recent or very ancient divergences. Nine of the 16 species represented by multiple haplotypes are not monophyletic. The apparent discordance between the recovered gene tree and the current taxonomic classification can be explained through: (a) collapsing of some species currently known only from either a sexual or a parthenogenetic generation into a single cyclically parthenogenetic entity; (b) sorting of ancestral polymorphism in diverging lineages, and (c) horizontal transfer of haplotypes, perhaps due to hybridization within glacial refuges. Our conclusions emphasise the need for careful intra-specific sampling when reconstructing phylogenies for radiations of closely related species and imply that for certain taxonomic groups full phylogenetic resolution (using molecular markers) may not be attainable.

Animals↗
PubMed2003

Out of Anatolia: longitudinal gradients in genetic diversity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae.

Many studies have addressed the latitudinal gradients in intraspecific genetic diversity of European taxa generated during postglacial range expansion from southern refugia. Although Asia Minor is known to be a centre of diversity for many taxa, relatively few studies have considered its potential role as a Pleistocene refugium or a potential source for more ancient westward range expansion into Europe. Here we address these issues for an oak gallwasp, Andricus quercustozae (Hymenoptera: Cynipidae), whose distribution extends from Morocco along the northern coast of the Mediterranean through Turkey to Iran. We use sequence data for a fragment of the mitochondrial gene cytochrome b and allele frequency data for 12 polymorphic allozyme loci to answer the following questions: (1). which regions represent current centres of genetic diversity for A. quercustozae? Do eastern populations represent one refuge or several discrete glacial refugia? (2). Can we infer the timescale and sequence of the colonization processes linking current centres of diversity? Our results suggest that A. quercustozae was present in five distinct refugia (Iberia, Italy, the Balkans, southwestern Turkey and northeastern Turkey) with recent genetic exchange between Italy and Hungary. Genetic diversity is greatest in the Turkish refugia, suggesting that European populations are either (a). derived from Asia Minor, or (b). subject to more frequent population bottlenecks. Although Iberian populations show the lowest diversity for putatively selectively neutral markers, they have colonized a new oak host and represent a genetically and biologically discrete entity within the species.

Animals↗
PubMed2003

Comparing bootstrap and posterior probability values in the four-taxon case.

Assessment of the reliability of a given phylogenetic hypothesis is an important step in phylogenetic analysis. Historically, the nonparametric bootstrap procedure has been the most frequently used method for assessing the support for specific phylogenetic relationships. The recent employment of Bayesian methods for phylogenetic inference problems has resulted in clade support being expressed in terms of posterior probabilities. We used simulated data and the four-taxon case to explore the relationship between nonparametric bootstrap values (as inferred by maximum likelihood) and posterior probabilities (as inferred by Bayesian analysis). The results suggest a complex association between the two measures. Three general regions of tree space can be identified: (1) the neutral zone, where differences between mean bootstrap and mean posterior probability values are not significant, (2) near the two-branch corner, and (3) deep in the two-branch corner. In the last two regions, significant differences occur between mean bootstrap and mean posterior probability values. Whether bootstrap or posterior probability values are higher depends on the data in support of alternative topologies. Examination of star topologies revealed that both bootstrap and posterior probability values differ significantly from theoretical expectations; in particular, there are more posterior probability values in the range 0.85-1 than expected by theory. Therefore, our results corroborate the findings of others that posterior probability values are excessively high. Our results also suggest that extrapolations from single topology branch-length studies are unlikely to provide any general conclusions regarding the relationship between bootstrap and posterior probability values.

Computer Simulation↗
PubMed2002

A maximum-likelihood analysis of eight phylogenetic markers in gallwasps (Hymenoptera: Cynipidae): implications for insect phylogenetic studies.

We assessed the utility of eight DNA sequence markers (5.8S rDNA, 18S rDNA, 28S rDNA, ITS regions, long-wavelength opsin, elongation factor 1-alpha, cytochrome b, and cytochrome oxidase I) in reconstructing phylogenetic relationships at various levels of divergence in gallwasps (Hymenoptera: Cynipidae), using a set of eight exemplar taxa. We report sequence divergence values and saturation levels and compare phylogenetic results of these sequences analyzed both separately and combined to a well-corroborated morphological phylogeny. Likelihood ratio tests were used to find the best evolutionary model fitting each of the markers. The likelihood model best explaining the data is, for most loci, parameter rich, with strong A-T bias for mitochondrial loci and strong rate heterogeneity for the majority of loci. Our data suggest that 28S rDNA, elongation factor 1-alpha, and long-wavelength opsin may be potentially useful markers for the resolution of cynipid and other insect within-family-level divergences (circa 50-100 mya old), whereas mitochondrial loci and ITS regions are most useful for lower-level phylogenetics. In contrast, the 18S rDNA marker is likely to be useful for the resolution of above-family-level relationships.

Animals↗
PubMed2002

Evolutionary shifts between host oak sections and host-plant organs in Andricus gallwasps.

Gall-inducing insects have especially intimate interactions with their host plants and generally show great specificity with regard to both the host-plant species and the organ (e.g. flower, leaf) galled. However, the relative roles of shifts between host species and between host-plant organs in the diversification of gall-inducers are uncertain. We employ a novel and general maximum-likelihood approach to show that shifts between host-plant organs occur at a significantly greater rate than shifts between host oak sections in European Andricus gallwasps. This suggests that speciation has more often been associated with gall location shifts than with colonization of new host-plant species, and implies that it may be easier for gall-inducers to colonize new plant organs than new plant species. Andricus gallwasps have complex life cycles, with obligate alternation of sexual and parthenogenetic generations. Our phylogenetic analyses show that a life cycle with both generations galling white oaks (section Quercus) is ancestral, with a single shift of the sexual generation onto black oaks (section Cerris) to generate a clade with a novel host-alternating life cycle. This new life cycle provided the opportunity for further speciation, but may have also increased the risk of extinction of one or both generations by the demographic requirement for co-existence of both host-plant groups. In summary, it appears that Andricus gallwasp radiation may be a two-level process. Speciation events often involve shifts in gall location on the same host species. However, there are only so many ways to gall an oak, and rare shifts to new oak sections may contribute greatly to long-term diversification by opening up whole new adaptive zones.

Adaptation, Physiological↗
PubMedDate not supplied

Conflicting phylogenetic signals at the base of the metazoan tree.

A phylogenetic framework is essential for under-standing the origin and evolution of metazoan development. Despite a number of recent molecular studies and a rich fossil record of sponges and cnidarians, the evolutionary relationships of the early branching metazoan groups to each other and to a putative outgroup, the choanoflagellates, remain uncertain. This situation may be the result of the limited amount of phylogenetic information found in single genes and the small number of relevant taxa surveyed. To alleviate the effect of these analytical factors in the phylogenetic recons-truction of early branching metazoan lineages, we cloned multiple protein-coding genes from two choanoflagellates and diverse sponges, cnidarians, and a ctenophore. Comparisons of sequences for alpha-tubulin, beta-tubulin, elongation factor 2, HSP90, and HSP70 robustly support the hypothesis that choanoflagellates are closely affiliated with animals. However, analyses of single and concatenated amino acid sequences fail to resolve the relationships either between early branching metazoan groups or between Metazoa and choano-flagellates. We demonstrate that variable rates of evolution among lineages, sensitivity of the analyses to taxon selection, and conflicts in the phylogenetic signal contained in different amino acid sequences obscure the phylogenetic associations among the early branching Metazoa. These factors raise concerns about the ability to resolve the phylogenetic history of animals with molecular sequences. A consensus view of animal evolution may require investigations of genome-scale characters.

Animals↗