Ribosomal RNA Operon Research Articles

Ribosomal RNA phylogeny of bacterial and fungal pathogens of Agriotes wireworms

Wireworms, the polyphagous larvae of click beetles belonging to the genus Agriotes (Coleoptera: Elateridae), are severe and widespread agricultural pests affecting numerous crops. Biological control agents and methods for this general pest are highly solicited. In a screening for microbial Agriotes pathogens, an intracellular bacterium and a mitosporic fungus were isolated. Phylogenetic analysis based on ribosomal RNA operon sequences of both micro-organisms corroborated their previous morphology-based taxonomic classification. The bacterial pathogen has been assigned to the taxonomic genus Rickettsiella (Gammaproteobacteria) wherein it represents a new pathotype, ‘Rickettsiella agriotidis’, that appears most closely related to subjective synonyms of the nomenclatural type species, Rickettsiella popillae. The fungal pathogen has been shown to belong to the form-species Beauveria bassiana, i.e., an obligate anamorph related to the genus Cordyceps (Ascomycota: Hypocreales). Furthermore, the B. bassiana strain from Agriotes has been shown to be potentially susceptible to identification by gIi-diagnosis, i.e., a diagnostic method making use of the strain-specific presence of self-splicing group-I introns within the ribosomal RNA operons of certain hyphomycetous fungi.

Comparison of the ribosomal RNA operon from Texas Phoenix decline and lethal yellowing phytoplasmas

A comparison was made of the two palm yellows phytoplasmas affecting palms to determine if the entire ribosomal RNA operon portion of the phytoplasma genome, or portions thereof, could account for the observed palm host differences. Polymerase chain reaction (PCR) was used to amplify a 5.0 kb DNA fragment consisting of the entire ribosomal RNA operon from a subgroup 16SrIV-D phytoplasma that causes Texas Phoenix palm decline (TPD) in cabbage (Sabal palmetto) palm in west central Florida and from a subgroup 16SrIV-A phytoplasma that causes lethal yellowing (LY) in coconut (Cocos nucifera) palm in Jamaica. Before the PCR reaction, we sequenced by 454 sequencing a draft genome of the coconut LY phytoplasma, strain LYFL, that infects C. nucifera in Florida, and obtained from this draft sequence both copies of the entire ribosomal operon. Sequence analysis of the ribosomal RNA operons from both the LY and TPD phytoplasmas revealed the gene composition and orientation for the operons to be 5′16S rRNA-tRNAIle-23S rRNA-5S rRNA3′ and a tRNAVal3′ downstream of the 5S rRNA gene. Based on molecular comparisons using the sequences of the ribosomal RNA operon, the TPD (16SrIV-D) strain was 98 % similar to the LY (16SrIV-A) strains.

Bacterial communities in different sections of a municipal wastewater treatment plant revealed by 16S rDNA 454 pyrosequencing

In this study, we successfully demonstrated that 454 pyrosequencing was a powerful approach for investigating the bacterial communities in the activated sludge, digestion sludge, influent, and effluent samples of a full scale wastewater treatment plant treating saline sewage. For each sample, 18,808 effective sequences were selected and utilized to do the bacterial diversity and abundance analysis. In total, 2,455, 794, 1,667, and 1,932 operational taxonomic units were obtained at 3 % distance cutoff in the activated sludge, digestion sludge, influent, and effluent samples, respectively. The corresponding most dominant classes in the four samples are Alphaproteobacteria, Thermotogae, Deltaproteobacteria, and Gammaproteobacteria. About 67 % sequences in the digestion sludge sample were found to be affiliated with the Thermotogales order. Also, these sequences were assigned into a recently proposed genus Kosmotoga by the Ribosomal Database Project classifier. In the effluent sample, we found high abundance of Mycobacterium and Vibrio, which are genera containing pathogenic bacteria. Moreover, in this study, we proposed a method to differentiate the “gene percentage” and “cell percentage” by using Ribosomal RNA Operon Copy Number Database.Electronic supplementary materialThe online version of this article (doi:10.1007/s00253-012-4082-4) contains supplementary material, which is available to authorized users.

Tidying up international nucleotide sequence databases: ecological, geographical and sequence quality annotation of its sequences of mycorrhizal fungi.

Sequence analysis of the ribosomal RNA operon, particularly the internal transcribed spacer (ITS) region, provides a powerful tool for identification of mycorrhizal fungi. The sequence data deposited in the International Nucleotide Sequence Databases (INSD) are, however, unfiltered for quality and are often poorly annotated with metadata. To detect chimeric and low-quality sequences and assign the ectomycorrhizal fungi to phylogenetic lineages, fungal ITS sequences were downloaded from INSD, aligned within family-level groups, and examined through phylogenetic analyses and BLAST searches. By combining the fungal sequence database UNITE and the annotation and search tool PlutoF, we also added metadata from the literature to these accessions. Altogether 35,632 sequences belonged to mycorrhizal fungi or originated from ericoid and orchid mycorrhizal roots. Of these sequences, 677 were considered chimeric and 2,174 of low read quality. Information detailing country of collection, geographical coordinates, interacting taxon and isolation source were supplemented to cover 78.0%, 33.0%, 41.7% and 96.4% of the sequences, respectively. These annotated sequences are publicly available via UNITE (http://unite.ut.ee/) for downstream biogeographic, ecological and taxonomic analyses. In European Nucleotide Archive (ENA; http://www.ebi.ac.uk/ena/), the annotated sequences have a special link-out to UNITE. We intend to expand the data annotation to additional genes and all taxonomic groups and functional guilds of fungi.

Genomic analysis reveals Lactobacillus sanfranciscensis as stable element in traditional sourdoughs

Sourdough has played a significant role in human nutrition and culture for thousands of years and is still of eminent importance for human diet and the bakery industry. Lactobacillus sanfranciscensis is the predominant key bacterium in traditionally fermented sourdoughs.The genome of L. sanfranciscensis TMW 1.1304 isolated from an industrial sourdough fermentation was sequenced with a combined Sanger/454-pyrosequencing approach followed by gap closing by walking on fosmids. The sequencing data revealed a circular chromosomal sequence of 1,298,316 bp and two additional plasmids, pLS1 and pLS2, with sizes of 58,739 bp and 18,715 bp, which are predicted to encode 1,437, 63 and 19 orfs, respectively. The overall GC content of the chromosome is 34.71%. Several specific features appear to contribute to the ability of L. sanfranciscensis to outcompete other bacteria in the fermentation. L. sanfranciscensis contains the smallest genome within the lactobacilli and the highest density of ribosomal RNA operons per Mbp genome among all known genomes of free-living bacteria, which is important for the rapid growth characteristics of the organism. A high frequency of gene inactivation and elimination indicates a process of reductive evolution. The biosynthetic capacity for amino acids scarcely availably in cereals and exopolysaccharides reveal the molecular basis for an autochtonous sourdough organism with potential for further exploitation in functional foods. The presence of two CRISPR/cas loci versus a high number of transposable elements suggests recalcitrance to gene intrusion and high intrinsic genome plasticity.

A multilocus sequence analysis approach to the phylogeny and taxonomy of the Halobacteriales

Members of the order Halobacteriales are obligate extreme halophiles that belong to the domain Archaea. The classification of the Halobacteriales currently relies on a polyphasic approach, which integrates phenotypic, genotypic and chemotaxonomic characterization. However, the most utilized genetic marker for phylogeny, the 16S rRNA gene, has multiple drawbacks for use with the Halobacteriales: the species of many genera exhibit large intragenic differences between multiple ribosomal RNA operons, the gene is too conserved to discriminate reliably at the species level and it appears to be the most frequently recombined gene between closely related species. Moreover, the Halobacteriales is a rapidly expanding group due to recent successes at cultivating novel strains from a diverse set of hypersaline environments; a fast, reliable, inexpensive, portable molecular method for discriminating species is required for their investigation. Recently, multilocus sequence analysis (MLSA) has been shown to be an effective tool for strain identification and taxonomic designation, even for those taxa that experience frequent lateral gene transfer and homologous recombination. In this study, MLSA was utilized for evolutionary and taxonomic investigation of the Halobacteriales. Efficacy of the MLSA approach was tested across a hierarchical gradient using 52 halobacterial strains, representing 33 species (including names without standing in nomenclature) and 14 genera. A subset of 21 strains from the genus Haloarcula was analysed separately to test the sensitivity and relevance of the MLSA approach among closely related strains and species. The results demonstrated that MLSA differentiated individual strains, reliably grouped strains into species and species into genera and identified potential novel species and also family-like relationships. This study demonstrates that MLSA is a rapid and informative molecular method that will probably accommodate strain analysis at any taxonomic level within the Halobacteriales.

DNA topology of highly transcribed operons in Salmonella enterica serovar Typhimurium

Bacteria differ from eukaryotes by having the enzyme DNA gyrase, which catalyses the ATP-dependent negative supercoiling of DNA. Negative supercoils are essential for condensing chromosomes into an interwound (plectonemic) and branched structure known as the nucleoid. Topo-1 removes excess supercoiling in an ATP-independent reaction and works with gyrase to establish a topological equilibrium where supercoils move within 10 kb domains bounded by stochastic barriers along the sequence. However, transcription changes the stochastic pattern by generating supercoil diffusion barriers near the sites of gene expression. Using supercoil-dependent Tn3 and γδ resolution assays, we studied DNA topology upstream, downstream and across highly transcribed operons. Whenever two Res sites flanked efficiently transcribed genes, resolution was inhibited and the loss in recombination efficiency was proportional to transcription level. Ribosomal RNA operons have the highest transcription rates, and resolution assays at the rrnG and rrnH operons showed inhibitory levels 40-100 times those measured in low-transcription zones. Yet, immediately upstream and downstream of RNA polymerase (RNAP) initiation and termination sites, supercoiling characteristics were similar to poorly transcribed zones. We present a model that explains why RNAP blocks plectonemic supercoil movement in the transcribed track and suggests how gyrase and TopA control upstream and downstream transcription-driven supercoiling.

Extreme Reconfiguration of Plastid Genomes in the Angiosperm Family Geraniaceae: Rearrangements, Repeats, and Codon Usage

Geraniaceae plastid genomes (plastomes) have experienced a remarkable number of genomic changes. The plastomes of Erodium texanum, Geranium palmatum, and Monsonia speciosa were sequenced and compared with other rosids and the previously published Pelargonium hortorum plastome. Geraniaceae plastomes were found to be highly variable in size, gene content and order, repetitive DNA, and codon usage. Several unique plastome rearrangements include the disruption of two highly conserved operons (S10 and rps2-atpA), and the inverted repeat (IR) region in M. speciosa does not contain all genes in the ribosomal RNA operon. The sequence of M. speciosa is unusually small (128,787 bp); among angiosperm plastomes sequenced to date, only those of nonphotosynthetic species and those that have lost one IR copy are smaller. In contrast, the plastome of P. hortorum is the largest, at 217,942 bp. These genomes have experienced numerous gene and intron losses and partial and complete gene duplications. Some of the losses are shared throughout the family (e.g., trnT-GGU and the introns of rps16 and rpl16); however, other losses are homoplasious (e.g., trnG-UCC intron in G. palmatum and M. speciosa). IR length is also highly variable. The IR in P. hortorum was previously shown to be greatly expanded to 76 kb, and the IR is lost in E. texanum and reduced in G. palmatum (11 kb) and M. speciosa (7 kb). Geraniaceae plastomes contain a high frequency of large repeats (>100 bp) relative to other rosids. Within each plastome, repeats are often located at rearrangement end points and many repeats shared among the four Geraniaceae flank rearrangement end points. GC content is elevated in the genomes and also in coding regions relative to other rosids. Codon usage per amino acid and GC content at third position sites are significantly different for Geraniaceae protein-coding sequences relative to other rosids. Our findings suggest that relaxed selection and/or mutational biases lead to increased GC content, and this in turn altered codon usage. We propose that increases in genomic rearrangements, repetitive DNA, nucleotide substitutions, and GC content may be caused by relaxed selection resulting from improper DNA repair.

Global transport of thermophilic bacteria in atmospheric dust

Aerosols from dust storms generated in the Sahara-Sahel desert area of Africa are transported north over Europe and periodically result in dry dust precipitation in the Mediterranean region. Samples of dust collected in Turkey and Greece following two distinct desert storm events contained viable thermophilic organisms of the genus Geobacillus, namely G. thermoglucosidasius and G. thermodenitrificans, and the recently reclassified Aeribacillus pallidus (formerly Geobacillus pallidus). We present here evidence that African dust storms create an atmospheric bridge between distant geographical regions and that they are also probably the source of thermophilic geobacilli later deposited over northern Europe by rainfall or dust plumes themselves. The same organisms (99% similarity in the 16S rDNA sequence) were found in dust collected in the Mediterranean region and inhabiting cool soils in Northern Ireland. This study also contributes new insights to the taxonomic identification of Geobacillus sp. Attempts to identify these organisms using 16S rRNA gene sequences have revealed that they contain multiple and diverse copies of the ribosomal RNA operon (up to 10 copies with nine different sequences), which dictates care in interpreting data about the systematics of this genus.

Metabolic streamlining in an open-ocean nitrogen-fixing cyanobacterium

Nitrogen (N(2))-fixing marine cyanobacteria are an important source of fixed inorganic nitrogen that supports oceanic primary productivity and carbon dioxide removal from the atmosphere. A globally distributed, periodically abundant N(2)-fixing marine cyanobacterium, UCYN-A, was recently found to lack the oxygen-producing photosystem II complex of the photosynthetic apparatus, indicating a novel metabolism, but remains uncultivated. Here we show, from metabolic reconstructions inferred from the assembly of the complete UCYN-A genome using massively parallel pyrosequencing of paired-end reads, that UCYN-A has a photofermentative metabolism and is dependent on other organisms for essential compounds. We found that UCYN-A lacks a number of major metabolic pathways including the tricarboxylic acid cycle, but retains sufficient electron transport capacity to generate energy and reducing power from light. Unexpectedly, UCYN-A has a reduced genome (1.44 megabases) that is structurally similar to many chloroplasts and some bacteria, in that it contains inverted repeats of ribosomal RNA operons. The lack of biosynthetic pathways for several amino acids and purines suggests that this organism depends on other organisms, either in close association or in symbiosis, for critical nutrients. However, size fractionation experiments using natural populations have so far not provided evidence of a symbiotic association with another microorganism. The UCYN-A cyanobacterium is a paradox in evolution and adaptation to the marine environment, and is an example of the tight metabolic coupling between microorganisms in oligotrophic oceanic microbial communities.

Co-Orientation of Replication and Transcription Preserves Genome Integrity

In many bacteria, there is a genome-wide bias towards co-orientation of replication and transcription, with essential and/or highly-expressed genes further enriched co-directionally. We previously found that reversing this bias in the bacterium Bacillus subtilis slows replication elongation, and we proposed that this effect contributes to the evolutionary pressure selecting the transcription-replication co-orientation bias. This selection might have been based purely on selection for speedy replication; alternatively, the slowed replication might actually represent an average of individual replication-disruption events, each of which is counter-selected independently because genome integrity is selected. To differentiate these possibilities and define the precise forces driving this aspect of genome organization, we generated new strains with inversions either over ∼1/4 of the chromosome or at ribosomal RNA (rRNA) operons. Applying mathematical analysis to genomic microarray snapshots, we found that replication rates vary dramatically within the inverted genome. Replication is moderately impeded throughout the inverted region, which results in a small but significant competitive disadvantage in minimal medium. Importantly, replication is strongly obstructed at inverted rRNA loci in rich medium. This obstruction results in disruption of DNA replication, activation of DNA damage responses, loss of genome integrity, and cell death. Our results strongly suggest that preservation of genome integrity drives the evolution of co-orientation of replication and transcription, a conserved feature of genome organization.

Intracellular Yeasts in Endothelial Cells of a Great Blue Heron (Ardea Herodias)

Intracellular organisms in the endothelial cells of several organs of an adult great blue heron (Ardea herodias) were identified as a yeast in the family Saccharomycetales based on ultrastructural morphology and sequence data from the ribosomal RNA operon. Morphologically similar organisms of unknown identity have been described previously in Muscovy (Cairina moschata) and domestic (Anas platyrhynchos domestica) ducks.

Coordination of genomic structure and transcription by the main bacterial nucleoid‐associated protein HU

The histone-like protein HU is a highly abundant DNA architectural protein that is involved in compacting the DNA of the bacterial nucleoid and in regulating the main DNA transactions, including gene transcription. However, the coordination of the genomic structure and function by HU is poorly understood. Here, we address this question by comparing transcript patterns and spatial distributions of RNA polymerase in Escherichia coli wild-type and hupA/B mutant cells. We demonstrate that, in mutant cells, upregulated genes are preferentially clustered in a large chromosomal domain comprising the ribosomal RNA operons organized on both sides of OriC. Furthermore, we show that, in parallel to this transcription asymmetry, mutant cells are also impaired in forming the transcription foci-spatially confined aggregations of RNA polymerase molecules transcribing strong ribosomal RNA operons. Our data thus implicate HU in coordinating the global genomic structure and function by regulating the spatial distribution of RNA polymerase in the nucleoid.

Nonoptimal Microbial Response to Antibiotics Underlies Suppressive Drug Interactions

Suppressive drug interactions, in which one antibiotic can actually help bacterial cells to grow faster in the presence of another, occur between protein and DNA synthesis inhibitors. Here, we show that this suppression results from nonoptimal regulation of ribosomal genes in the presence of DNA stress. Using GFP-tagged transcription reporters in Escherichia coli, we find that ribosomal genes are not directly regulated by DNA stress, leading to an imbalance between cellular DNA and protein content. To test whether ribosomal gene expression under DNA stress is nonoptimal for growth rate, we sequentially deleted up to six of the seven ribosomal RNA operons. These synthetic manipulations of ribosomal gene expression correct the protein-DNA imbalance, lead to improved survival and growth, and completely remove the suppressive drug interaction. A simple mathematical model explains the nonoptimal regulation in different nutrient environments. These results reveal the genetic mechanism underlying an important class of suppressive drug interactions.

Characterization of a Phytoplasma Associated with Frogskin Disease in Cassava.

Cassava frogskin disease (CFSD) is an economically important root disease of cassava (Manihot esculenta) in Colombia and other South American countries, including Brazil, Venezuela, Peru, Costa Rica, and Panama. The roots of severely affected plants are thin, making them unsuitable for consumption. In Colombia, phytoplasma infections were confirmed in 35 of 39 genotypes exhibiting mild or severe CFSD symptoms either by direct or nested polymerase chain reaction (PCR) assays employing ribosomal (r)RNA operon primer pairs. The CFSD-associated phytoplasmas were identified as group 16SrIII strains by restriction fragment length polymorphism (RFLP) and sequence analyses of amplified rDNA products, and results were corroborated by PCRs employing group 16SrIII-specific rRNA gene or ribosomal protein (rp) gene primers. Collectively, RFLP analyses indicated that CFSD strains differed from all phytoplasmas described previously in group 16SrIII and, on this basis, the strains were tentatively assigned to new ribosomal and ribosomal protein subgroups 16SrIII-L and rpIII-H, respectively. This is the first molecular identification of a phytoplasma associated with CFSD in cassava in Colombia.

Fine tuning of the E. coli NusB:NusE complex affinity to BoxA RNA is required for processive antitermination

Phage λ propagation in Escherichia coli host cells requires transcription antitermination on the λ chromosome mediated by λN protein and four host Nus factors, NusA, B, E (ribosomal S10) and G. Interaction of E. coli NusB:NusE heterodimer with the single stranded BoxA motif of λnutL or λnutR RNA is crucial for this reaction. Similarly, binding of NusB:NusE to a BoxA motif is essential to suppress transcription termination in the ribosomal RNA (rrn) operons. We used fluorescence anisotropy to measure the binding properties of NusB and of NusB:NusE heterodimer to BoxA-containing RNAs differing in length and sequence. Our results demonstrate that BoxA is necessary and sufficient for binding. We also studied the gain-of-function D118N NusB mutant that allows λ growth in nusA1 or nusE71 mutants. In vivo λ burst-size determinations, CD thermal unfolding measurements and X-ray crystallography of this as well as various other NusB D118 mutants showed the importance of size and polarity of amino acid 118 for RNA binding and other interactions. Our work suggests that the affinity of the NusB:NusE complex to BoxA RNA is precisely tuned to maximize control of transcription termination.

Understanding the Differences between Genome Sequences of Escherichia coli B Strains REL606 and BL21(DE3) and Comparison of the E. coli B and K-12 Genomes

Each difference between the genome sequences of Escherichia coli B strains REL606 and BL21(DE3) can be interpreted in light of known laboratory manipulations plus a gene conversion between ribosomal RNA operons. Two treatments with 1-methyl-3-nitro-1-nitrosoguanidine in the REL606 lineage produced at least 93 single-base-pair mutations (∼90% GC-to-AT transitions) and 3 single-base-pair GC deletions. Two UV treatments in the BL21(DE3) lineage produced only 4 single-base-pair mutations but 16 large deletions. P1 transductions from K-12 into the two B lineages produced 317 single-base-pair differences and 9 insertions or deletions, reflecting differences between B DNA in BL21(DE3) and integrated restriction fragments of K-12 DNA inherited by REL606. Two sites showed selective enrichment of spontaneous mutations. No unselected spontaneous single-base-pair mutations were evident. The genome sequences revealed that a progenitor of REL606 had been misidentified, explaining initially perplexing differences. Limited sequencing of other B strains defined characteristic properties of B and allowed assembly of the inferred genome of the ancestral B of Delbrück and Luria. Comparison of the B and K-12 genomes shows that more than half of the 3793 proteins of their basic genomes are predicted to be identical, although ∼310 appear to be functional in either B or K-12 but not in both. The ancestral basic genome appears to have had ∼4039 coding sequences occupying ∼4.0 Mbp. Repeated horizontal transfer from diverged Escherichia coli genomes and homologous recombination may explain the observed variable distribution of single-base-pair differences. Fifteen sites are occupied by phage-related elements, but only six by comparable elements at the same site. More than 50 sites are occupied by IS elements in both B and K, 16 in common, and likely founding IS elements are identified. A signature of widespread cryptic phage P4-type mobile elements was identified. Complex deletions (dense clusters of small deletions and substitutions) apparently removed nonessential genes from ∼30 sites in the basic genomes.

Two Structurally Independent Domains of E. coli NusG Create Regulatory Plasticity via Distinct Interactions with RNA Polymerase and Regulators

NusG is a conserved regulatory protein that interacts with elongation complexes (ECs) of RNA polymerase, DNA, and RNA to modulate transcription in multiple and sometimes opposite ways. In Escherichia coli, NusG suppresses pausing and increases elongation rate, enhances termination by E. coli ρ and phage HK022 Nun protein, and promotes antitermination by λN and in ribosomal RNA operons. We report NMR studies that suggest that E. coli NusG consists of two largely independent N- and C-terminal structural domains, NTD and CTD, respectively. Based on tests of the functions of the NTD and CTD and variants of NusG in vivo and in vitro, we find that NTD alone is sufficient to suppress pausing and enhance transcript elongation in vitro. However, neither domain alone can enhance ρ-dependent termination or support antitermination, indicating that interactions of both domains with ECs are required for these processes. We propose that the two domains of NusG mediate distinct interactions with ECs: the NTD interacts with RNA polymerase and the CTD interacts with ρ and other regulators, providing NusG with different combinations of interactions to effect different regulatory outcomes.

Molecular characterization of the sequences of the 16S-23S rDNA internal spacer region (ISR) from isolates of Taylorella asinigenitalis

BackgroundSequence information on the 16S-23S rDNA internal spacer region (ISR) exhibits a large degree of sequence and length variation at both the genus and species levels. A primer pair for the amplification of 16S-23S rDNA ISR generated three amplicons for each of isolates of Taylorella asinigenitalis (UCD-1T, UK-1 and UK-2).FindingsFollowing TA cloning and sequencing, the three isolates of T. asinigenitalis were demonstrated to possess three ISR units of different lengths. Although the three corresponding ISRs (A, B and C) were identified to be identical to each other (UK-1 and UK-2 isolates), the ISRs shared approximately 95.3–98.9% nucleotide sequence similarities between the UCD-1T and UK-1/-2 isolates. A typical order of two intercistronic tRNA genes (5'-tRNAIle-tRNAAla-3') with the different nucleotide spacers [44 through 51 base pairs (bp)] in length was identified among the isolates. The consensus sequences of the antiterminators of boxB and boxA were also identified in all ISRs. Thus, three ISRs were identified for each isolate, and therefore, at least three distinctly different ribosomal RNA operons were suggested to occur in the genome of T. asinigenitalis. This was also confirmed by Southern hybridization procedure.ConclusionThe present study represents a dendrogram constructed based on the nucleotide sequence data of 16S-23S rDNA ISR for T. asinigenitalis, which may aid in the phylogenetic positioning of T. asinigenitalis within the genus Taylorella, and in the molecular discrimination of T. asinigenitalis.

The Complete Plastid Genome Sequence of the Secondarily Nonphotosynthetic Alga Cryptomonas paramecium: Reduction, Compaction, and Accelerated Evolutionary Rate

The cryptomonads are a group of unicellular algae that acquired photosynthesis through the engulfment of a red algal cell, a process called secondary endosymbiosis. Here, we present the complete plastid genome sequence of the secondarily nonphotosynthetic species Cryptomonas paramecium CCAP977/2a. The ∼78 kilobase pair (Kbp) C. paramecium genome contains 82 predicted protein genes, 29 transfer RNA genes, and a single pseudogene (atpF). The C. paramecium plastid genome is approximately 50 Kbp smaller than those of the photosynthetic cryptomonads Guillardia theta and Rhodomonas salina; 71 genes present in the G. theta and/or R. salina plastid genomes are missing in C. paramecium. The pet, psa, and psb photosynthetic gene families are almost entirely absent. Interestingly, the ribosomal RNA operon, present as inverted repeats in most plastid genomes (including G. theta and R. salina), exists as a single copy in C. paramecium. The G + C content (38%) is higher in C. paramecium than in other cryptomonad plastid genomes, and C. paramecium plastid genes are characterized by significantly different codon usage patterns and increased evolutionary rates. The content and structure of the C. paramecium plastid genome provides insight into the changes associated with recent loss of photosynthesis in a predominantly photosynthetic group of algae and reveals features shared with the plastid genomes of other secondarily nonphotosynthetic eukaryotes.