Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha‐proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein‐coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome‐wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad‐scale genome streamlining via loss of protein‐coding genes as well as noncoding, parasitic DNA elements.     

Retention and loss of RNA interference pathways in trypanosomatid protozoans

RNA interference (RNAi) pathways are widespread in metaozoans but the genes required show variable occurrence or activity in eukaryotic microbes, including many pathogens. While some Leishmania lack RNAi activity and Argonaute or Dicer genes, we show that Leishmania braziliensis and other species within the Leishmania subgenus Viannia elaborate active RNAi machinery. Strong attenuation of expression from a variety of reporter and endogenous genes was seen. As expected, RNAi knockdowns of the sole Argonaute gene implicated this protein in RNAi. The potential for functional genetics was established by testing RNAi knockdown lines lacking the paraflagellar rod, a key component of the parasite flagellum. This sets the stage for the systematic manipulation of gene expression through RNAi in these predominantly diploid asexual organisms, and may also allow selective RNAi-based chemotherapy. Functional evolutionary surveys of RNAi genes established that RNAi activity was lost after the separation of the Leishmania subgenus Viannia from the remaining Leishmania species, a divergence associated with profound changes in the parasite infectious cycle and virulence. The genus Leishmania therefore offers an accessible system for testing hypothesis about forces that may select for the loss of RNAi during evolution, such as invasion by viruses, changes in genome plasticity mediated by transposable elements and gene amplification (including those mediating drug resistance), and/or alterations in parasite virulence.     

Evolution of microbial genomes: sequence acquisition and loss

We present models describing the acquisition and deletion of novel sequences in populations of microorganisms. We infer that most novel sequences are neutral. Thus, sequence duplications and gene transfer between organisms sharing the same environment are rarely expected to generate adaptive functions. Two classes of models are considered: (1) a homogeneous population with constant size, and (2) an island model in which the population is subdivided into patches that are in contact through slow migration. Distributions of gene frequencies are derived in a Moran model with overlapping generations. We find that novel, neutral or near-neutral coding sequences in microorganisms will not be fixed globally because they offer large target sizes for mutations and because the populations are so large. At most, such genes may have a transient presence in only a small fraction of the population. Consequently, a microbial population is expected to have a very large diversity of transient neutral gene content. Only sequences that are under strong selection, globally or in individual patches, can be expected to persist. We suggest that genome size is maintained in microorganisms by a quasi-steady state mechanism in which random fluctuations in the effective acquisition and deletion rates result in genome sizes that vary from patch to patch. We assign the genomic identity of a global population to those genes that are required for the participation of patches in the genetic sweeps that maintain the genomic coherence of the population. In contrast, we stress the influence of sequence loss on the isolation and the divergence (speciation) of novel patches from a global population.     

Origin and evolution of the Amyrel gene in the α-amylase multigene family of Diptera

Alpha-amylase genes often form multigene families in living organisms. In Diptera, a remote paralog, Amyrel, had been discovered in Drosophila, where this gene is currently used as a population and phylogenetic marker. The putative encoded protein has about 40% divergence with the classical amylases. We have searched the presence of the paralog in other families of Diptera to track its origin and understand its evolution. Amyrel was detected in a number of families of Muscomorpha (Brachycera-Cyclorrapha), suggesting an origin much older than previously thought. It has not been found elsewhere to date, and it is absent from the Anopheles gambiae genome. The intron–exon structures of the genes found so far suggest that the ancestral gene (before the duplication which gave rise to Amyrel) had two introns, and that subsequent, repeated and independent loss of one or both introns occurred in some Muscomorpha families. It seems that the Amyrel protein has experienced specific amino acid substitutions in regions generally well conserved in amylases, raising the possibility of peculiar, functional adaptations of this protein.     

Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Evidence suggests that duplicated genes are retained at a much higher rate than originally thought and that functional divergence of gene copies is a major factor promoting their retention in the genome. We find that two Drosophila testes-specific alpha4 proteasome subunit genes (alpha4-t1 and alpha4-t2) have a higher polymorphism within species and are significantly more diverged between species than the somatic alpha4 gene. Our data suggest that following gene duplication, the alpha4-t1 gene experienced relaxed selective constraints, whereas the alpha4-t2 gene experienced positive selection acting on several codons. We report significant heterogeneity in evolutionary rates among all three paralogs at homologous codons, indicating that functional divergence has coincided with genic divergence. Reproductive subfunctionalization may allow for a more rapid evolution of reproductive traits and a greater specialization of testes function. Our data add to the increasing evidence that duplicated genes experience lower selective constraints and in some cases positive selection following duplication. Newly duplicated genes that are freer from selective constraints may provide a mechanism for developing new interactions and a pathway for the evolution of new genes.     

Intron loss mediated structural dynamics and functional differentiation of the polygalacturonase gene family in land plants

The polygalacturonase (PG) gene family is one of the largest gene families in plants. PGs are involved in various plant development steps. The evolutionary processes accounting for the functional divergence and the specialized functions of PGs in land plants are unclear. Whole sets of PG genes were retrieved from the genome web sites of model organisms in algae and land plants. The number of PG genes was expanded by lineage-specific manner with the biological complexity of the organism. Differentiation of PGs was related with phylogenetic hierarchy such as presence of rhamno-PGs from algae to plants, endo- and exo-PGs in land plants, exo-PGs in flowering plants. Gene structure analysis revealed that land plant PG genes resulted from differential intron gain and loss, with the latter event predominating. Differential intron losses partitioned the PGs into separate clades to be expressed differentially during plant development. Intron position and phase were not conserved between PGs of algae and land plants but conserved among PG genes of land plants from moss to vascular plants, indicating that the current introns in the PGs in land plants appeared after the split between unicellular algae and multicelluar land plants. The results demonstrate that the functional divergence and differentiation of PGs in land plants is attributable to intron losses.     

Comparable Rates of Gene Loss and Functional Divergence after Genome Duplications Early in Vertebrate Evolution

Duplicated genes are an important source of new protein functions and novel developmental and physiological pathways. Whereas most models for fate of duplicated genes show that they tend to be rapidly lost, models for pathway evolution suggest that many duplicated genes rapidly acquire novel functions. Little empirical evidence is available, however, for the relative rates of gene loss vs. divergence to help resolve these contradictory expectations. Gene families resulting from genome duplications provide an opportunity to address this apparent contradiction. With genome duplication, the number of duplicated genes in a gene family is at most 2(n), where n is the number of duplications. The size of each gene family, e.g., 1, 2, 3, . . . , 2(n), reflects the patterns of gene loss vs. functional divergence after duplication. We focused on gene families in humans and mice that arose from genome duplications in early vertebrate evolution and we analyzed the frequency distribution of gene family size, i.e., the number of families with two, three or four members. All the models that we evaluated showed that duplicated genes are almost as likely to acquire a new and essential function as to be lost through acquisition of mutations that compromise protein function. An explanation for the unexpectedly high rate of functional divergence is that duplication allows genes to accumulate more neutral than disadvantageous mutations, thereby providing more opportunities to acquire diversified functions and pathways.     

Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets

In many animals, male seminal proteins influence gamete interactions and fertilization ability and are probably involved in barriers to gene flow between diverging lineages. Here we use a proteomic approach to identify seminal proteins that are transferred to females during copulation and that may be involved in fertilization barriers between two hybridizing field crickets (Gryllus firmus and Gryllus pennsylvanicus). Analyses of patterns of divergence suggest that much of the field cricket genome has remained undifferentiated following the evolution of reproductive isolation. By contrast, seminal protein genes are highly differentiated. Tests of selection reveal that positive selection is likely to be responsible for patterns of differentiation. Together, our observations suggest that some of the loci encoding seminal proteins may indeed play a role in fertilization barriers in field crickets.     

After the duplication: gene loss and adaptation in Saccharomyces genomes

The ancient duplication of the Saccharomyces cerevisiae genome and subsequent massive loss of duplicated genes is apparent when it is compared to the genomes of related species that diverged before the duplication event. To learn more about the evolutionary effects of the duplication event, we compared the S. cerevisiae genome to other Saccharomyces genomes. We demonstrate that the whole genome duplication occurred before S. castellii diverged from S. cerevisiae. In addition to more accurately dating the duplication event, this finding allowed us to study the effects of the duplication on two separate lineages. Analyses of the duplication regions of the genomes indicate that most of the duplicated genes (approximately 85%) were lost before the speciation. Only a small amount of paralogous gene loss (4-6%) occurred after speciation. On the other hand, S. castellii appears to have lost several hundred genes that were not retained as duplicated paralogs. These losses could be related to genomic rearrangements that reduced the number of chromosomes from 16 to 9. In addition to S. castellii, other Saccharomyces sensu lato species likely diverged from S. cerevisiae after the duplication. A thorough analysis of these species will likely reveal other important outcomes of the whole genome duplication.     

Assessing the performance of single-copy genes for recovering robust phylogenies

Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.     

Complex prokaryotic genome structure: rapid evolution of chromosome II

Although many bacteria with two chromosomes have been sequenced, the roles of such complex genome structuring are still unclear. To uncover levels of chromosome I (CI) and chromosome II (CII) sequence divergence, Mauve 2.2.0 was used to align the CI- and CII-specific sequences of bacteria with complex genome structuring in two sets of comparisons: the first set was conducted among the CI and CII of bacterial strains of the same species, while the second set was conducted among the CI and CII of species in Alphaproteobacteria that possess two chromosomes. The analyses revealed a rapid evolution of CII-specific DNA sequences compared with CI-specific sequences in a majority of organisms. In addition, levels of protein divergence between CI-specific and CII-specific genes were determined using phylogenetic analyses and confirmed the DNA alignment findings. Analysis of synonymous and nonsynonymous substitutions revealed that the structural and functional constraints on CI and CII genes are not significantly different. Also, horizontal gene transfer estimates in selected organisms demonstrated that CII in many species has acquired higher levels of horizontally transferred segments than CI. In summary, rapid evolution of CII may perform particular roles for organisms such as aiding in adapting to specialized niches.     

The evolutionary history of human and chimpanzee Y-chromosome gene loss

Recent studies have suggested that gene gain and loss may contribute significantly to the divergence between humans and chimpanzees. Initial comparisons of the human and chimpanzee Y-chromosomes indicate that chimpanzees have a disproportionate loss of Y-chromosome genes, which may have implications for the adaptive evolution of sex-specific as well as reproductive traits, especially because one of the genes lost in chimpanzees is critically involved in spermatogenesis in humans. Here we have characterized Y-chromosome sequences in gorilla, bonobo, and several chimpanzee subspecies for 7 chimpanzee gene-disruptive mutations. Our analyses show that 6 of these gene-disruptive mutations predate chimpanzee-bonobo divergence at approximately 1.8 MYA, which indicates significant Y-chromosome change in the chimpanzee lineage relatively early in the evolutionary divergence of humans and chimpanzees.     

Mitochondrial genome organization and cytoplasmic male sterility in plants

Plant mitochondrial genomes are much larger and more complex than those of other eukaryotic organisms. They contain a very active recombination system and have a multipartite genome organization with a master circle resolving into two or more subgenomic circles by recombination through repeated sequences. Their protein coding capacity is very low and is comparable to that of animal and fungal systems. Several subunits of mitochondrial functional complexes, a complete set of tRNAs and 26S, 18S and 5S rRNAs are coded by the plant mitochondrial genome. The protein coding genes contain group II introns. The organelle genome contains stretches of DNA sequences homologous to chloroplast DNA. It also contains actively transcribed DNA sequences having open reading frames. Plasmid like DNA molecules are found in mitochondria of some plants Cytoplasmic male sterility in plants, characterized by failure to produce functional pollen grains, is a maternally inherited trait. This phenomenon has been found in many species of plants and is conveniently used for hybrid plant production. The genetic determinants for cytoplasmic male sterility reside in the mitochondrial genome. Some species of plants exhibit more than one type of cytoplasmic male sterility. Several nuclear genes are known to control expression of cytoplasmic male sterility. Different cytoplasmic male sterility types are distinguished by their specific nuclear genes(rfs) which restore pollen fertility. Cytoplasmic male sterility types are also characterized by mitochondrial DNA restriction fragment length polymorphism patterns, variations in mitochondrial RNAs, differences in protein synthetic profiles, differences in sensitivity to fungal toxins and insecticides, presence of plasmid DNAs or RNAs and also presence of certain unique sequences in the genome. Recently nuclear male sterility systems based on (i) over expression of agrobacterialrol C gene and (ii) anther specific expression of an RNase gene have been developed in tobacco andBrassica by genetic engineering methods.     

Structure of the mouse glial fibrillary acidic protein gene: implications for the evolution of the intermediate filament multigene family.

We report the complete sequence of the gene encoding mouse glial fibrillary acidic protein (GFAP), the intermediate filament (IF) protein specific to astrocytes. The 9.8 kb gene includes nine exons separated by introns ranging in size from 0.2 to 2.5 kb. A comparison of the organization of the GFAP gene with that of genes encoding other IF proteins reveals that the structure of IF genes is highly conserved in spite of considerable divergence at the amino acid level. Thus, most of the evolutionary events leading to the placement of introns in IF genes must have occurred prior to the duplication and subsequent divergence of IF genes from a presumptive common ancestral sequence. The conserved gene organization is unrelated to structural features of IF proteins. A curious feature of the GFAP gene is the large number of repeated sequences found in the introns. Six tracts of reiterated di- or trinucleotides are present, plus tandem repeats of two different novel sequences. One repeat is unique to the GFAP gene; the other occurs elsewhere in the mouse genome, although at relatively low frequency.     

Evolution by leaps: gene duplication in bacteria

Background  

Sequence related families of genes and proteins are common in bacterial genomes. In Escherichia coli they constitute over half of the genome. The presence of families and superfamilies of proteins suggest a history of gene duplication and divergence during evolution. Genome encoded protein families, their size and functional composition, reflect metabolic potentials of the organisms they are found in. Comparing protein families of different organisms give insight into functional differences and similarities.