Tiling microarray analysis of rice chromosome 10 to identify the transcriptome and relate its expression to chromosomal architecture

Background  

Sequencing and annotation of the genome of rice (Oryza sativa) have generated gene models in numbers that top all other fully sequenced species, with many lacking recognizable sequence homology to known genes. Experimental evaluation of these gene models and identification of new models will facilitate rice genome annotation and the application of this knowledge to other more complex cereal genomes.     

Origin and Evolution of Rickettsial Plasmids

Background

Rickettsia species are strictly intracellular bacteria that have undergone a reductive genomic evolution. Despite their allopatric lifestyle, almost half of the 26 currently validated Rickettsia species have plasmids. In order to study the origin, evolutionary history and putative roles of rickettsial plasmids, we investigated the evolutionary processes that have shaped 20 plasmids belonging to 11 species, using comparative genomics and phylogenetic analysis between rickettsial, microbial and non-microbial genomes.

Results

Plasmids were differentially present among Rickettsia species. The 11 species had 1 to 4 plasmid (s) with a size ranging from 12 kb to 83 kb. We reconstructed pRICO, the last common ancestor of the current rickettsial plasmids. pRICO was vertically inherited mainly from Rickettsia/Orientia chromosomes and diverged vertically into a single or multiple plasmid(s) in each species. These plasmids also underwent a reductive evolution by progressive gene loss, similar to that observed in rickettsial chromosomes, possibly leading to cryptic plasmids or complete plasmid loss. Moreover, rickettsial plasmids exhibited ORFans, recent gene duplications and evidence of horizontal gene transfer events with rickettsial and non-rickettsial genomes mainly from the α/γ-proteobacteria lineages. Genes related to maintenance and plasticity of plasmids, and to adaptation and resistance to stress mostly evolved under vertical and/or horizontal processes. Those involved in nucleotide/carbohydrate transport and metabolism were under the influence of vertical evolution only, whereas genes involved in cell wall/membrane/envelope biogenesis, cycle control, amino acid/lipid/coenzyme and secondary metabolites biosynthesis, transport and metabolism underwent mainly horizontal transfer events.

Conclusion

Rickettsial plasmids had a complex evolution, starting with a vertical inheritance followed by a reductive evolution associated with increased complexity via horizontal gene transfer as well as gene duplication and genesis. The plasmids are plastic and mosaic structures that may play biological roles similar to or distinct from their co-residing chromosomes in an obligate intracellular lifestyle.     

Rickettsial evolution in the light of comparative genomics

Rickettsia are best known as strictly intracellular vector‐borne bacteria that cause mild to severe diseases in humans and other animals. Recent advances in molecular tools and biological experiments have unveiled a wide diversity of Rickettsia spp. that include species with a broad host range and some species that act as endosymbiotic associates. Molecular phylogenies of Rickettsia spp. contain some ambiguities, such as the position of R. canadensis and relationships within the spotted fever group. In the modern era of genomics, with an ever‐increasing number of sequenced genomes, there is enhanced interest in the use of whole‐genome sequences to understand pathogenesis and assess evolutionary relationships among rickettsial species. Rickettsia have small genomes (1.1–1.5 Mb) as a result of reductive evolution. These genomes contain split genes, gene remnants and pseudogenes that, owing to the colinearity of some rickettsial genomes, may represent different steps of the genome degradation process. Genomics reveal extreme genome reduction and massive gene loss in highly vertebrate‐pathogenic Rickettsia compared to less virulent or endosymbiotic species. Information gleaned from rickettsial genomics challenges traditional concepts of pathogenesis that focused primarily on the acquisition of virulence factors. Another intriguing phenomenon about the reduced rickettsial genomes concerns the large fraction of non‐coding DNA and possible functionality of these “non‐coding” sequences, because of the high conservation of these regions. Despite genome streamlining, Rickettsia spp. contain gene families, selfish DNA, repeat palindromic elements and genes encoding eukaryotic‐like motifs. These features participate in sequence and functional diversity and may play a crucial role in adaptation to the host cell and pathogenesis. Genome analyses have identified a large fraction of mobile genetic elements, including plasmids, suggesting the possibility of lateral gene transfer in these intracellular bacteria. Phylogenetic analyses have identified several candidates for horizontal gene acquisition among Rickettsia spp. including tra, pat2, and genes encoding for the type IV secretion system and ATP/ADP translocase that may have been acquired from bacteria living in amoebae. Gene loss, gene duplication, DNA repeats and lateral gene transfer all have shaped rickettsial genome evolution. A comprehensive analysis of the entire genome, including genes and non‐coding DNA, will help to unlock the mysteries of rickettsial evolution and pathogenesis.     

Evolution of four gene families with patchy phylogenetic distributions: influx of genes into protist genomes

Background  

Lateral gene transfer (LGT) in eukaryotes from non-organellar sources is a controversial subject in need of further study. Here we present gene distribution and phylogenetic analyses of the genes encoding the hybrid-cluster protein, A-type flavoprotein, glucosamine-6-phosphate isomerase, and alcohol dehydrogenase E. These four genes have a limited distribution among sequenced prokaryotic and eukaryotic genomes and were previously implicated in gene transfer events affecting eukaryotes. If our previous contention that these genes were introduced by LGT independently into the diplomonad and Entamoeba lineages were true, we expect that the number of putative transfers and the phylogenetic signal supporting LGT should be stable or increase, rather than decrease, when novel eukaryotic and prokaryotic homologs are added to the analyses.     

The complete mitochondrial genome of the oriental fruit moth <Emphasis Type="Italic">Grapholita molesta</Emphasis> (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae) currently is one of the economically most destructive pest species of stone and pome fruits worldwide. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,776 bp long, with an A + T content of 81.24%, containing 37 typical animal mitochondrial genes and an A + T-rich region. All gene are arranged as hypothesized ancestral gene order of insects except for trnM, which was shuffled from 3′ downstream of trnQ to 5′ upstream of trnI. cox1 gene uses unusual CGA start codon, as that in all other sequenced lepidopteran mitochondrial genome. The secondary structures for the two rRNA genes were predicted. All helices typically present in insect mitochondrial rRNA genes are generated. A microsatellite sequence was inserted into the region of H2347 in rrnL in G. molesta and two other sequenced tortricid mitochondrial genomes, indicating that the insertion event in this helix might occurred anciently in family Tortricidae. All of the 22 typical animal tRNA genes have a typical cloverleaf structure except for trnS2, in which the D-stem pairings in the DHU arm are absent. An intergenic sequence is present between trnQ and nad2 as well as in other sequenced lepidopteran mitochondrial genomes, which was presumed to be a remnant of trnM gene and its boundary sequences after the duplication of trnM to the upstream of trnI in Lepidoptera. The A + T-rich region is 836 bp, containing six repeat sequences of “TTATTATTATTATTAAATA(G)TTT.”     

Sequencing and comparison of the Rickettsia genomes from the whitefly Bemisia tabaci Middle East Asia Minor I

The whitefly, Bemisia tabaci, harbors the primary symbiont ‘Candidatus Portiera aleyrodidarum’ and a variety of secondary symbionts. Among these secondary symbionts, Rickettsia is the only one that can be detected both inside and outside the bacteriomes. Infection with Rickettsia has been reported to influence several aspects of the whitefly biology, such as fitness, sex ratio, virus transmission and resistance to pesticides. However, mechanisms underlying these differences remain unclear, largely due to the lack of genomic information of Rickettsia. In this study, we sequenced the genome of two Rickettsia strains isolated from the Middle East Asia Minor 1 (MEAM1) species of the B. tabaci complex in China and Israel. Both Rickettsia genomes were of high coding density and AT‐rich, containing more than 1000 coding sequences, much larger than that of the coexisted primary symbiont, Portiera. Moreover, the two Rickettsia strains isolated from China and Israel shared most of the genes with 100% identity and only nine genes showed sequence differences. The phylogenetic analysis using orthologs shared in the genus, inferred the proximity of Rickettsia in MEAM1 and Rickettsia bellii. Functional analysis revealed that Rickettsia was unable to synthesize amino acids required for complementing the whitefly nutrition. Besides, a type IV secretion system and a number of virulence‐related genes were detected in the Rickettsia genome. The presence of virulence‐related genes might benefit the symbiotic life of the bacteria, and hint on potential effects of Rickettsia on whiteflies. The genome sequences of Rickettsia provided a basis for further understanding the function of Rickettsia in whiteflies.     

Comparative BAC end sequence analysis of tomato and potato reveals overrepresentation of specific gene families in potato

Background  

Tomato (Solanum lycopersicon) and potato (S. tuberosum) are two economically important crop species, the genomes of which are currently being sequenced. This study presents a first genome-wide analysis of these two species, based on two large collections of BAC end sequences representing approximately 19% of the tomato genome and 10% of the potato genome.     

OrthoParaMap: Distinguishing orthologs from paralogs by integrating comparative genome data and gene phylogenies

Background  

In eukaryotic genomes, most genes are members of gene families. When comparing genes from two species, therefore, most genes in one species will be homologous to multiple genes in the second. This often makes it difficult to distinguish orthologs (separated through speciation) from paralogs (separated by other types of gene duplication). Combining phylogenetic relationships and genomic position in both genomes helps to distinguish between these scenarios. This kind of comparison can also help to describe how gene families have evolved within a single genome that has undergone polyploidy or other large-scale duplications, as in the case of Arabidopsis thaliana – and probably most plant genomes.     

Complete plastid genome sequences of Drimys, Liriodendron, and Piper : implications for the phylogenetic relationships of magnoliids

Background  

The magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content.     

Using complementary approaches to identify trans‐domain nuclear gene transfers in the extremophile Galdieria sulphuraria (Rhodophyta)

Identification of horizontal gene transfers (HGTs) has primarily relied on phylogenetic tree based methods, which require a rich sampling of sequenced genomes to ensure a reliable inference. Because the success of phylogenetic approaches depends on the breadth and depth of the database, researchers usually apply stringent filters to detect only the most likely gene transfers in the genomes of interest. One such study focused on a highly conservative estimate of trans‐domain gene transfers in the extremophile eukaryote, Galdieria sulphuraria (Galdieri) Merola (Rhodophyta), by applying multiple filters in their phylogenetic pipeline. This led to the identification of 75 inter‐domain acquisitions from Bacteria or Archaea. Because of the evolutionary, ecological, and potential biotechnological significance of foreign genes in algae, alternative approaches and pipelines complementing phylogenetics are needed for a more comprehensive assessment of HGT. We present here a novel pipeline that uncovered 17 novel foreign genes of prokaryotic origin in G. sulphuraria, results that are supported by multiple lines of evidence including composition‐based, comparative data, and phylogenetics. These genes encode a variety of potentially adaptive functions, from metabolite transport to DNA repair.     

Hon-yaku: a biology-driven Bayesian methodology for identifying translation initiation sites in prokaryotes

Background  

Computational prediction methods are currently used to identify genes in prokaryote genomes. However, identification of the correct translation initiation sites remains a difficult task. Accurate translation initiation sites (TISs) are important not only for the annotation of unknown proteins but also for the prediction of operons, promoters, and small non-coding RNA genes, as this typically makes use of the intergenic distance. A further problem is that most existing methods are optimized for Escherichia coli data sets; applying these methods to newly sequenced bacterial genomes may not result in an equivalent level of accuracy.     

Application of quantitative measures of gene order similarity to phylogenetic reconstructions (exemplified by bacteria of the genus <Emphasis Type="Italic">Rickettsia</Emphasis>)

Data reflecting evolutionary changes in chromosomal gene order can be used for phylogenetic reconstructions along with the results of nucleotide sequence comparison. By the example of bacteria of the genus Rickettsia, we have shown that phylogenetic reconstructions based on quantitative estimates of the similarity and cladistic analysis of gene order data, may, in some cases, amend and fill up classical phylogenetic trees. When applied, these approaches enabled us to substantiate the hypothesis that Rickettsia felis species had split before the typhus (R. typhi, R. prowazekii) and spotted fever (R. connorii) group divergence and thus R. felis does not belong to the latter group. In general, rickettsias evolved towards increasing intracellular parasitic specialization. Five Rickettsia species whose genomes have been sequenced and annotated completely actually form an evolutionary series R. bellii—R. felis—R. conorii—R. prowazekii—R. typhi. Within this series, a reduction in genome size and rapid decrease of genome rearrangement rates (genome plasticity loss) gradually occur.     

Rickettsia Phylogenomics: Unwinding the Intricacies of Obligate Intracellular Life

Background

Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular α-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs).

Methodology/Principal Findings

We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (∼1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions.

Conclusion/Significance

Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets.     

Lineage-specific variations of congruent evolution among DNA sequences from three genomes,and relaxed selective constraints on <Emphasis Type="Italic">rbc</Emphasis>L in <Emphasis Type="Italic">Cryptomonas</Emphasis> (Cryptophyceae)

Background  

Plastid-bearing cryptophytes like Cryptomonas contain four genomes in a cell, the nucleus, the nucleomorph, the plastid genome and the mitochondrial genome. Comparative phylogenetic analyses encompassing DNA sequences from three different genomes were performed on nineteen photosynthetic and four colorless Cryptomonas strains. Twenty-three rbc L genes and fourteen nuclear SSU rDNA sequences were newly sequenced to examine the impact of photosynthesis loss on codon usage in the rbc L genes, and to compare the rbc L gene phylogeny in terms of tree topology and evolutionary rates with phylogenies inferred from nuclear ribosomal DNA (concatenated SSU rDNA, ITS2 and partial LSU rDNA), and nucleomorph SSU rDNA.     

The genus Cyclops (Copepoda,Cyclopoida) in Europe

Although copepods of the genus Cyclops are among the most common and dominant plankton taxa of lakes in the northern temperate zone, their taxonomy is still unclear. We analysed an extensive array of Cyclops populations from Europe by means of molecular methods and evaluated morphological characters. Altogether, 68 populations of Cyclops species were sampled, assigned to morphospecies and sequenced for the 12S rRNA gene. Selected populations of each morphospecies were additionally sequenced for three mitochondrial (16S rRNA, cytochrome b, COI) and two nuclear genes (18S rRNA, ITS1) and analysed for micromorphological traits. Our analysis revealed fifteen lineages that can be regarded as separate species. Thirteen of these match currently accepted species, while the remaining two lineages were distinct from the other described species. Thus, their taxonomic status is open to further studies. Besides taxonomy, our study gives new insights into the ecology, distribution and phylogenetic relationships of these species. Finally, a set of morphological traits was selected to facilitate identification.     

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

Principal Findings

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

Conclusions/Significance

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.     

Uncovering Cryptic Diversity in Two Amoebozoan Species Using Complete Mitochondrial Genome Sequences

The Amoebozoa are a major eukaryotic lineage that encompasses a wide range of amoeboid organisms. The group is understudied from a systematic perspective: molecular tools have only been applied in the last 15 yr. Hence, there is an undersampling of both genes and taxa in the group especially compared to plants, animals, and fungi. Here, we present the complete mitochondrial genomes of two ubiquitous and abundant morpho‐species (Acanthamoeba castellanii and Vermamoeba vermiformis). Both have mitochondrial genomes of close relatives previously available, enabling insights into recent divergences at a genomic scale, while simultaneously offering comparisons with divergence estimates obtained from traditionally used single genes, SSU rDNA and cox1. The newly sequenced mt genomes are significantly divergent from their previously sequenced conspecifics (A. castellannii 16.4% divergence at nucleotide level and 10.4% amino acid; V. vermiformis 21.6% and 13.1%, respectively), while divergence at the small subunit ribosomal DNA is below 1% within both species. Morphological analyses determined that these lineages are indistinguishable from their previously sequenced counterparts. Phylogenetic reconstructions using 26 mt genes also indicate a level of divergence that is comparable to divergence among species, while reconstructions using the small subunit ribosomal DNA (SSU rDNA) do not. In addition, we demonstrate that between closely related taxa, there are high levels of synteny, which can be explored for primer design to obtain larger fragments than the traditional barcoding genes. We conclude that, although most systematic work has relied on SSU, this gene alone can severely underestimate diversity. Thus, we suggest that the mt genome emerges as an alternative for unraveling the lower level phylogenetic relationships of Amoebozoa.     

Draft genome of the American Eel (Anguilla rostrata)

Freshwater eels (Anguilla sp.) have large economic, cultural, ecological and aesthetic importance worldwide, but they suffered more than 90% decline in global stocks over the past few decades. Proper genetic resources, such as sequenced, assembled and annotated genomes, are essential to help plan sustainable recoveries by identifying physiological, biochemical and genetic mechanisms that caused the declines or that may lead to recoveries. Here, we present the first sequenced genome of the American eel. This genome contained 305 043 contigs (N50 = 7397) and 79 209 scaffolds (N50 = 86 641) for a total size of 1.41 Gb, which is in the middle of the range of previous estimations for this species. In addition, protein‐coding regions, including introns and flanking regions, are very well represented in the genome, as 95.2% of the 458 core eukaryotic genes and 98.8% of the 248 ultra‐conserved subset were represented in the assembly and a total of 26 564 genes were annotated for future functional genomics studies. We performed a candidate gene analysis to compare three genes among all three freshwater eel species and, congruent with the phylogenetic relationships, Japanese eel (A. japanica) exhibited the most divergence. Overall, the sequenced genome presented in this study is a crucial addition to the presently available genetic tools to help guide future conservation efforts of freshwater eels.