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.     

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co‐occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next‐generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein‐encoding genes, with the Complex IV genes (coxI‐III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.     

High selectivity in symbiotic associations of lichenized ascomycetes and cyanobacteria

The selectivity of mycobionts and cyanobionts in lichen symbioses were examined. We analyzed symbiotic cyanobionts, collected from different sample sites, and compared them to free‐living cyanobacteria Nostoc. Cyanobionts were obtained from lichens assigned to the genera Pseudocyphellaria and Sticta, in particular. Multiple gene loci were screened and direct optimization was used in the phylogenetic analyses. We show that many lichen fungi are strongly selective towards their cyanobionts. Lichenized ascomycetes seem to be able to identify and choose a specific strain, species or a species group of Nostoc with which to associate. The present analyses also suggest that some of the Nostoc taxa may be specialized in symbiotic life with only lichenized ascomycetes. Despite the selectivity observed in fungi, there appears to be no coevolution between the partners. We have also discussed the problems of using the tRNA^Leu intron as a marker in phylogenetic analyses. © The Willi Hennig Society 2006.     

Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

Nuclear and mitochondrial subunits from the white shrimp Litopenaeus vannamei F0F1 ATP-synthase complex: cDNA sequence, molecular modeling, and mRNA quantification of atp9 and atp6

We studied for the first time the ATP-synthase complex from shrimp as a model to understand the basis of crustacean bioenergetics since they are exposed to endogenous processes as molting that demand high amount of energy. We analyzed the cDNA sequence of two subunits of the Fo sector from mitochondrial ATP-synthase in the white shrimp Litopenaeus vannamei. The nucleus encoded atp9 subunit presents a 773 bp sequence, containing a signal peptide sequence only observed in crustaceans, and the mitochondrial encoded atp6 subunit presents a sequence of 675 bp, and exhibits high identity with homologous sequences from invertebrate species. ATP9 and ATP6 protein structural models interaction suggest specific functional characteristics from both proteins in the mitochondrial enzyme. Differences in the steady-state mRNA levels of atp9 and atp6 from five different tissues correlate with tissue function. Moreover, significant changes in the mRNA levels of both subunits at different molt stages were detected. We discussed some insights about the enzyme structure and the regulation mechanisms from both ATP-synthase subunits related to the energy requirements of shrimp.     

The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)₁₂ element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.     

Mitochondrial sequences of <Emphasis Type="Italic">Seriatopora</Emphasis> corals show little agreement with morphology and reveal the duplication of a tRNA gene near the control region

The taxonomy of corals of the genus Seriatopora has not previously been studied using molecular sequence markers. As a first step toward a re-evaluation of species boundaries in this genus, mitochondrial sequence variability was analyzed in 51 samples collected from Okinawa, New Caledonia, and the Philippines. Four clusters of sequences were detected that showed little concordance with species currently recognized on a morphological basis. The most likely explanation is that the skeletal characters used for species identification are highly variable (polymorphic or phenotypically plastic); alternative explanations include introgression/hybridization, or deep coalescence and the retention of ancestral mitochondrial polymorphisms. In all individuals sequenced, two copies of trnW were found on either side of the atp8 gene near the putative D-loop, a novel mitochondrial gene arrangement that may have arisen from a duplication of the trnW-atp8 region followed by a deletion of one atp8.     

Organization of repetitive DNAs and the genomic regions carrying ribosomal RNA,<Emphasis Type="Italic"> cob</Emphasis>, and<Emphasis Type="Italic"> atp9</Emphasis> genes in the cucurbit mitochondrial genomes

Plants in the genus Cucumis (cucumber and melon) have the largest mitochondrial genomes known among all plants, due in part to the accumulation of repetitive DNAs of varying complexities. Recombination among these repetitive DNAs should produce highly rearranged mitochondrial genomes relative to the smaller mitochondrial genomes of related plants. We cloned and sequenced mitochondrial genomic regions near the rRNA, atp9 and cob genes from cucumber, melon, squash and watermelon (all members of the Cucurbitaceae family), and compared to the previously sequenced mitochondrial genomes of Arabidopsis thaliana and sugar beet to study the distribution and arrangement of coding and repetitive DNAs. Cucumber and melon had regions of concentrated repetitive DNAs spread throughout the sequenced regions; few repetitive DNAs were revealed in the mitochondrial genomes of A. thaliana, sugar beet, squash and watermelon. Recombination among these repetitive DNAs most likely produced unique arrangements of the rrn18 and rrn5 genes in the genus Cucumis. Cucumber mitochondrial DNA had more pockets of dispersed direct and inverted repeats than melon and the other plants, and we did not reveal repetitive sequences significantly contributing to mitochondrial genome expansion in both cucumber and melon.Disclaimer. Names are necessary to report factually on available data; however, the U.S. Department of Agriculture (USDA) neither guarantees nor warrants the standard of the product, and the use of the name by USDA implies no approval of the product to the exclusion of others that may also be suitable.Communicated by R. Hagemann     

Construction and characterization of a full-length cDNA library from mycobiont of <Emphasis Type="Italic">Endocarpon pusillum</Emphasis> (lichen-forming Ascomycota)

The first full-length cDNA library for lichenized fungi was constructed from cultured mycobiont of the arid desert lichen Endocarpon pusillum. Based on small-scale sequencing results, 111 genes of the lichenized fungi were identified for the first time, among which 11 genes shared no homology with any known fungal genes. Real-time PCR showed that the size of the mycobiont genome is 39.13 Mb and the copy number of ribosomal RNA gene repeat units is 43. The results of this study will be valuable for the ongoing lichen genome-sequencing project and the large-scale identification of functional genes from lichenized fungi.     

Comparative genomics of large mitochondria in placozoans

The first sequenced mitochondrial genome of a placozoan, Trichoplax adhaerens, challenged the conventional wisdom that a compact mitochondrial genome is a common feature among all animals. Three additional placozoan mitochondrial genomes representing highly divergent clades have been sequenced to determine whether the large Trichoplax mtDNA is a shared feature among members of the phylum Placozoa or a uniquely derived condition. All three mitochondrial genomes were found to be very large, 32- to 37-kb, circular molecules, having the typical 12 respiratory chain genes, 24 tRNAs, rnS, and rnL. They share with the Trichoplax mitochondrial genome the absence of atp8, atp9, and all ribosomal protein genes, the presence of several cox1 introns, and a large open reading frame containing an intron group I LAGLIDADG endonuclease domain. The differences in mtDNA size within Placozoa are due to variation in intergenic spacer regions and the presence or absence of long open reading frames of unknown function. Phylogenetic analyses of the 12 respiratory chain genes support the monophyly of Placozoa. The similarities in composition and structure between the three mitochondrial genomes reported here and that of Trichoplax's mtDNA suggest that their uncompacted state is a shared ancestral feature to other nonmetazoans while their gene content is a derived feature shared only among the Metazoa.     

宁夏贺兰山青海云杉根系共生真菌的分离与鉴定

[背景]青海云杉(Picea crassifolia)是中国特有植物,具有极高的生态价值,是西北地区重要的森林更新树种和荒山造林树种,其生长发育及其抗逆性与根系共生真菌多样性密切相关.[目的]从青海云杉根系分离并鉴定定殖的可培养共生真菌,阐明青海云杉根系可培养共生真菌的种类组成,为共生真菌在青海云杉育苗、造林及生态恢复...     

A genome scan of diversifying selection in Ophiocordyceps zombie‐ant fungi suggests a role for enterotoxins in co‐evolution and host specificity

Identification of the genes underlying adaptation sheds light on the biological functions targeted by natural selection. Searches for footprints of positive selection, in the form of rapid amino acid substitutions, and the identification of species‐specific genes have proved to be powerful approaches to identifying the genes involved in host specialization in plant‐pathogenic fungi. We used an evolutionary comparative genomic approach to identify genes underlying host adaptation in the ant‐infecting genus Ophiocordyceps, which manipulates ant behaviour. A comparison of the predicted genes in the genomes of species from three species complexes—O. unilateralis, O. australis and O. subramanianii—revealed an enrichment in pathogenesis‐associated functions, including heat‐labile enterotoxins, among species‐specific genes. Furthermore, these genes were overrepresented among those displaying significant footprints of positive selection. Other categories of genes suspected to be important for virulence and pathogenicity in entomopathogenic fungi (e.g., chitinases, lipases, proteases, core secondary metabolism genes) were much less represented, although a few candidate genes were found to evolve under positive selection. An analysis including orthologs from other entomopathogenic fungi in a broader context showed that positive selection on enterotoxins was specific to the ant‐infecting genus Ophiocordyceps. Together with previous studies reporting the overexpression of an enterotoxin during behavioural manipulation in diseased ants, our findings suggest that heat‐labile enterotoxins are important effectors in host adaptation and co‐evolution in the Ophiocordyceps entomopathogenic fungi.     

Mitochondrial gene diversity associated with the atp9 stop codon in natural populations of wild carrot (Daucus carota ssp. carota)

Mitochondrial genomes extracted from the wild populations of Daucus carota have been used as a genetic resource by breeders of cultivated carrot, yet little is known concerning the extent of their diversity in nature. Of special interest is an SNP in the putative stop codon of the mitochondrial gene atp9 that has been associated previously with male-sterile and male-fertile phenotypic variants. In this study, either the sequence or PCR/RFLP genotypes were obtained from the mitochondrial genes atp1, atp9, and cox1 found in D. carota individuals collected from 24 populations in the eastern United States. More than half of the 128 individuals surveyed had a CAA or AAA, rather than TAA, genotype at the position usually thought to function as an atp9 stop codon in this species. We also found no evidence for mitochondrial RNA editing (Cytosine to Uridine) of the CAA stop codon in either floral or leaf tissue. Evidence for intragenic recombination, as opposed to the more common intergenic recombination in plant mitochondrial genomes, in our data set is presented. Indel and SNP variants elsewhere in atp9, and in the other 2 genes surveyed, were nonrandomly associated with the 3 atp9 stop codon variants, though further analysis suggested that multilocus genotypic diversity had been enhanced by recombination. Overall the mitochondrial genetic diversity was only modestly structured among populations with an F(ST) of 0.34.     

Comparison of closely related,uncultivated Coxiella tick endosymbiont population genomes reveals clues about the mechanisms of symbiosis

Understanding the symbiotic interaction between Coxiella‐like endosymbionts (CLE) and their tick hosts is challenging due to lack of isolates and difficulties in tick functional assays. Here we sequenced the metagenome of a CLE population from wild Rhipicephalus sanguineus ticks (CRs) and compared it to the previously published genome of its close relative, CLE of R. turanicus (CRt). The tick hosts are closely related sympatric species, and their two endosymbiont genomes are highly similar with only minor differences in gene content. Both genomes encode numerous pseudogenes, consistent with an ongoing genome reduction process. In silico flux balance metabolic analysis (FBA) revealed the excess production of L‐proline for both genomes, indicating a possible proline transport from Coxiella to the tick. Additionally, both CR genomes encode multiple copies of the proline/betaine transporter, proP gene. Modelling additional Coxiellaceae members including other tick CLE, did not identify proline as an excreted metabolite. Although both CRs and CRt genomes encode intact B vitamin synthesis pathway genes, which are presumed to underlay the mechanism of CLE‐tick symbiosis, the FBA analysis indicated no changes for their products. Therefore, this study provides new testable hypotheses for the symbiosis mechanism and a better understanding of CLE genome evolution and diversity.     

Genome size‐dependent pcna gene copy number in dinoflagellates and molecular evidence of retroposition as a major evolutionary mechanism

Proliferating cell nuclear antigen (PCNA) plays critical roles in eukaryotic DNA replication and replication‐associated processes. It is typically encoded by one or two gene copies (pcna) in eukaryotic genomes. Recently reported higher copy numbers of pcna in some dinoflagellates raised a question of how this gene has uniquely evolved in this phylum. Through real‐time PCR quantification, we found a wide range of pcna copy number (2–287 copies) in 11 dinoflagellate species (n = 38), and a strong positive correlation between pcna copy number and genome size (log₁₀–log₁₀ transformed). Intraspecific pcna diverged up to 21% and are dominated by nonsynonymous substitutions, indicating strong purifying selection pressure on and hence functional necessity of this gene. By surveying pcna copy numbers in eukaryotes, we observed a genome size threshold at 4 pg DNA, above which more than two pcna copies are found. To examine whether retrotransposition is a mechanism of pcna duplication, we measured the copy number of retroposed pcna, taking advantage of the 22‐nt dinoflagellate‐specific spliced leader (DinoSL) capping the 5′ end of dinoflagellate nuclear‐encoded mRNAs, which would exist in the upstream region of a retroposed gene copy. We found that retroposed pcna copy number increased with total pcna copy number and genome size. These results indicate co‐evolution of dinoflagellate pcna copy number with genome size, and retroposition as a major mechanism of pcna duplication in dinoflagellates. Furthermore, we posit that the demand of faithful replication and maintenance of the large dinoflagellate genomes might have favored the preservation of the retroposed pcna as functional genes.     

Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids

Phenotypically diverse Lake Malawi cichlids exhibit similar genomes. The extensive sharing of genetic polymorphism among forms has both intrigued and frustrated biologists trying to understand the nature of diversity in this and other rapidly evolving systems. Shared polymorphism might result from hybridization and/or the retention of ancestrally polymorphic alleles. To examine these alternatives, we used new genomic tools to characterize genetic differentiation in widespread, geographically structured populations of Labeotropheus fuelleborni and Metriaclima zebra. These phenotypically distinct species share mitochondrial DNA (mtDNA) haplotypes and show greater mtDNA differentiation among localities than between species. However, Bayesian analysis of nuclear single nucleotide polymorphism (SNP) data revealed two distinct genetic clusters corresponding perfectly to morphologically diagnosed L. fuelleborni and M. zebra. This result is a function of the resolving power of the multi‐locus dataset, not a conflict between nuclear and mitochondrial partitions. Locus‐by‐locus analysis showed that mtDNA differentiation between species (F_CT) was nearly identical to the median single‐locus SNP F_CT. Finally, we asked whether there is evidence for gene flow at sites of co‐occurrence. We used simulations to generate a null distribution for the level of differentiation between co‐occurring populations of L. fuelleborni and M. zebra expected if there was no hybridization. The null hypothesis was rejected for the SNP data; populations that co‐occur at rock reef sites were slightly more similar than expected by chance, suggesting recent gene flow. The coupling of numerous independent markers with extensive geographic sampling and simulations utilized here provides a framework for assessing the prevalence of gene flow in recently diverged species.     

Nonneutral evolution of organelle genes in Silene vulgaris

Knowledge of mitochondrial gene evolution in angiosperms has taken a dramatic shift within the past decade, from universal slow rates of nucleotide change to a growing realization of high variation in rates among lineages. Additionally, evidence of paternal inheritance of plant mitochondria and recombination among mitochondrial genomes within heteroplasmic individuals has led to speculation about the potential for independent evolution of organellar genes. We report intraspecific mitochondrial and chloroplast sequence variation in a cosmopolitan sample of 42 Silene vulgaris individuals. There was remarkably high variation in two mitochondrial genes (atp1, atp9) and additional variation within a third gene (cob). Tests for patterns of nonneutral evolution were significant for atp1 and atp9, indicative of the maintenance of balanced polymorphisms. Two chloroplast genes (matK, ndhF) possessed less, but still high, variation and no divergence from neutral expectations. Phylogenetic patterns of organelle genes in both the chloroplast and mitochondria were incongruent, indicating the potential for independent evolutionary trajectories. Evidence indicated reassociation among cytoplasmic genomes and recombination between mitochondrial genes and within atp1, implying transient heteroplasmy in ancestral lineages. Although the mechanisms for long-term maintenance of mitochondrial polymorphism are currently unknown, frequency-dependent selection on linked cytoplasmic male sterility genes is a potential candidate.