Assembly and comparison of two closely related Brassica napus genomes

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor‐bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor‐bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.     

Dictyochophyceae Plastid Genomes Reveal Unusual Variability in Their Organization

Dictyochophyceae (silicoflagellates) are unicellular freshwater and marine algae (Heterokontophyta, stramenopiles). Despite their abundance in global oceans and potential ecological significance, discovered in recent years, neither nuclear nor organellar genomes of representatives of this group were sequenced until now. Here, we present the first complete plastid genome sequences of Dictyochophyceae, obtained from four species: Dictyocha speculum, Rhizochromulina marina, Florenciella parvula and Pseudopedinella elastica. Despite their comparable size and genetic content, these four plastid genomes exhibit variability in their organization: plastid genomes of F. parvula and P. elastica possess conventional quadripartite structure with a pair of inverted repeats, R. marina instead possesses two direct repeats with the same orientation and D. speculum possesses no repeats at all. We also observed a number of unusual traits in the plastid genome of D. speculum, including expansion of the intergenic regions, presence of an intron in the otherwise non‐intron‐bearing psaA gene, and an additional copy of the large subunit of RuBisCO gene (rbcL), the last of which has never been observed in any plastid genome. We conclude that despite noticeable gene content similarities between the plastid genomes of Dictyochophyceae and their relatives (pelagophytes, diatoms), the number of distinctive features observed in this lineage strongly suggests that additional taxa require further investigation.     

Synechococcus: 3 billion years of global dominance

Cyanobacteria are among the most important primary producers on the Earth. However, the evolutionary forces driving cyanobacterial species diversity remain largely enigmatic due to both their distinction from macro‐organisms and an undersampling of sequenced genomes. Thus, we present a new genome of a Synechococcus‐like cyanobacterium from a novel evolutionary lineage. Further, we analyse all existing 16S rRNA sequences and genomes of Synechococcus‐like cyanobacteria. Chronograms showed extremely polyphyletic relationships in Synechococcus, which has not been observed in any other cyanobacteria. Moreover, most Synechococcus lineages bifurcated after the Great Oxidation Event, including the most abundant marine picoplankton lineage. Quantification of horizontal gene transfer among 70 cyanobacterial genomes revealed significant differences among studied genomes. Horizontal gene transfer levels were not correlated with ecology, genome size or phenotype, but were correlated with the age of divergence. All findings were synthetized into a novel model of cyanobacterial evolution, characterized by serial convergence of the features, that is multicellularity and ecology.     

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is ‘open’, with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.     

Comparative Whole-Genome Analysis of Clinical Isolates Reveals Characteristic Architecture of Mycobacterium tuberculosis Pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.     

Reconstruction of ancient genome and gene order from complete microbial genome sequences

Microbial genome sequences provide us with the fossil records for inferring their origination and evolution. Assuming that current microbial genomes are the evolutionary results of ancient genomes or fragments and the neighboring genes in ancient genomes are more likely neighbors in current genomes, in this paper we proposed a paleontological algorithm and assembled the orthologous gene groups from 66 complete and current microbial genome sequences into a pseudo-ancient genome, which consists of continuous fragments of various sizes. We performed bootstrap resampling and correlation analyses and the results showed that the assembled ancient genome and fragments are statistically significant and the genes of the same fragment are inherently related and likely derived from common ancestors. This method provides a new computational tool for studying microbial genome structure and evolution.     

Draft genomes and genomic divergence of two Lepidurus tadpole shrimp species (Crustacea,Branchiopoda, Notostraca)

Crustaceans of the order Notostraca (Branchiopoda) are distributed worldwide and are known for the remarkable morphological stasis between their extant and Permian fossil species. Moreover, these crustaceans show relevant ecological traits and a wide range of reproductive strategies. However, genomic studies on notostracans are fairly limited. Here, we present the genome sequences of two notostracan taxa, Lepidurus arcticus and Lepidurus apus lubbocki. Taking advantage of the small genome sizes (~0.11 pg) of these taxa, genomes were sequenced for one individual per species with one run on the Illumina HiSeq X platform. We finally assembled 73.2 Mbp (L. arcticus) and 90.3 Mbp (L. apus lubbocki) long genomes. Assemblies cover up to 84% of the estimated genome size, with a gene completeness >97% for both genomes. In total, 13%–16% of the assembled genomes consist of repeats, and based on read mapping, L. apus lubbocki shows a significantly lower transposable element content than L. arcticus. The analysis of 2,376 orthologous genes indicates an ~7% divergence between the two Lepidurus taxa, with a nucleotide substitution rate significantly lower than that of Daphnia taxa. K_a/K_s analysis suggests purifying selection in both branchiopod lineages, raising the question of whether the low substitution rate of Lepidurus is correlated with morphological conservation or is linked to specific biological traits. Our analysis demonstrates that, in these organisms, it is possible to obtain high‐quality draft genomes from single individuals with a relatively low sequencing effort. This result makes Lepidurus and Notostraca interesting models for genomic studies at taxonomic, ecological and evolutionary levels.     

Cryptic sex? Estimates of genome exchange in unisexual mole salamanders (Ambystoma sp.)

Cryptic sex has been argued to explain the exceptional longevity of certain parthenogenetic vertebrate lineages, yet direct measurements of genetic exchange between sexual and apparently parthenogenetic forms are rare. Female unisexual mole salamanders (Ambystoma sp.) are the oldest known unisexual vertebrate lineage (~5 million years), and one hypothesis for their persistence is that allopolyploid female unisexuals periodically exchange haploid genomes ‘genome exchange’ during gynogenetic reproduction with males from sympatric sexual species. We test this hypothesis by using genome‐specific microsatellite DNA markers to estimate the rates of genome exchange between sexual males and unisexual females in two ponds in NE Ohio. We also test the prediction that levels of gene flow should be higher for ‘sympatric’ (sexual males present) genomes in unisexuals compared to ‘allopatric’ (sexual males absent) unisexual genomes. We used a model testing framework in the coalescent‐based program MIGRATE‐N to compare models where unidirectional gene flow is present and absent between sexual species and unisexuals. As predicted, our results show higher levels of gene flow between sexuals and sympatric unisexual genomes compared to lower (likely artefactual) levels of gene flow between sexuals and allopatric unisexual genomes. Our results provide direct evidence that genome exchange between sexual and unisexual Ambystoma occurs and demonstrate that the magnitude depends on which sexual species are present. The relatively high levels of gene flow suggest that unisexuals must be at a selective advantage over sexual forms so as to avoid extinction due to genetic swamping through genome exchange.     

Genomic signatures and co‐occurrence patterns of the ultra‐small Saccharimonadia (phylum CPR/Patescibacteria) suggest a symbiotic lifestyle

The size of bacterial genomes is often associated with organismal metabolic capabilities determining ecological breadth and lifestyle. The recently proposed Candidate Phyla Radiation (CPR)/Patescibacteria encompasses mostly unculturable bacterial taxa with relatively small genome sizes with potential for co‐metabolism interdependencies. As yet, little is known about the ecology and evolution of CPR, particularly with respect to how they might interact with other taxa. Here, we reconstructed two novel genomes (namely, Candidatus Saccharibacter sossegus and Candidatus Chaer renensis) of taxa belonging to the class Saccharimonadia within the CPR/Patescibacteria using metagenomes obtained from acid mine drainage (AMD). By testing the hypothesis of genome streamlining or symbiotic lifestyle, our results revealed clear signatures of gene losses in these genomes, such as those associated with de novo biosynthesis of essential amino acids, nucleotides, fatty acids and cofactors. In addition, co‐occurrence analysis provided evidence supporting potential symbioses of these organisms with Hydrotalea sp. in the AMD system. Together, our findings provide a better understanding of the ecology and evolution of CPR/Patescibacteria and highlight the importance of genome reconstruction for studying metabolic interdependencies between unculturable Saccharimonadia representatives.     

Nucleomorph genome diversity and its phylogenetic implications in cryptomonad algae

The relationship between phylogeny and nucleomorph genome size was examined in 16 strains of cryptomonad algae using pulsed‐field gel electrophoresis, Southern hybridization and phylogenetic analyses. Our results suggest that all cryptomonads examined in this study contain three nucleomorph chromosomes and their total genome size ranges from 495 to 750 kb. In addition, we estimated the plastid genome size of the respective organisms. The plastid genomes of photosynthetic strains were approximately 120–160 kb in size, whereas the non‐photosynthetic Cryptomonas paramecium NIES715 possesses a genome of approximately 70 kb. Phylogenetic analysis of the nuclear small subunit ribosomal DNA (SSU rDNA) gene showed that nucleomorph genome size varies considerably within closely related strains. This result indicates that the reduction of nucleomorph genomes is a rapid phenomenon that occurred multiple times independently during cryptomonad evolution. The nucleomorph genome sizes of Cryptomonas rostratiformis NIES277 appeared to be approximately 495 kb. This is smaller than that of Guillardia theta CCMP327, which until now was thought to have the smallest known nucleomorph genome size among photosynthetic cryptomonads.     

Phylogenomics of Enterococcus faecalis from wild birds: new insights into host-associated differences in core and accessory genomes of the species

Wild birds have been suggested to be reservoirs of antimicrobial resistant and/or pathogenic Enterococcus faecalis (Efs) strains, but the scarcity of studies and available sequences limit our understanding of the population structure of the species in these hosts. Here, we analysed the clonal and plasmid diversity of 97 Efs isolates from wild migratory birds. We found a high diversity, with most sequence types (STs) being firstly described here, while others were found in other hosts including some predominant in poultry. We found that pheromone-responsive plasmids predominate in wild bird Efs while 35% of the isolates entirely lack plasmids. Then, to better understand the ecology of the species, the whole genome of fivestrains with known STs (ST82, ST170, ST16 and ST55) were sequenced and compared with all the Efs genomes available in public databases. Using several methods to analyse core and accessory genomes (AccNET, PLACNET, hierBAPS and PANINI), we detected differences in the accessory genome of some lineages (e.g. ST82) demonstrating specific associations with birds. Conversely, the genomes of other Efs lineages exhibited divergence in core and accessory genomes, reflecting different adaptive trajectories in various hosts. This pangenome divergence, horizontal gene transfer events and occasional epidemic peaks could explain the population structure of the species.     

Improving the genetic representation of rare taxa within complex microbial communities using DNA normalization methods

Complex microbial communities typically contain a large number of low abundance species, which collectively, comprise a considerable proportion of the community. This ‘rare biosphere’ has been speculated to contain keystone species and act as a repository of genomic diversity to facilitate community adaptation. Many environmental microbes are currently resistant to cultivation, and can only be accessed via culture‐independent approaches. To enhance our understanding of the role of the rare biosphere, we aimed to improve their metagenomic representation using DNA normalization methods, and assess normalization success via shotgun DNA sequencing. A synthetic metagenome was constructed from the genomic DNA of five bacterial species, pooled in a defined ratio spanning three orders of magnitude. The synthetic metagenome was fractionated and thermally renatured, allowing the most abundant sequences to hybridize. Double‐stranded DNA was removed either by hydroxyapatite chromatography, or by a duplex‐specific nuclease (DSN). The chromatographic method failed to enrich for the genomes present in low starting abundance, whereas the DSN method resulted in all genomes reaching near equimolar abundance. The representation of the rarest member was increased by approximately 450‐fold. De novo assembly of the normalized metagenome enabled up to 18.0% of genes from the rarest organism to be assembled, in contrast to the un‐normalized sample, where genes were not able to be assembled at the same sequencing depth. This study has demonstrated that the application of normalization methods to metagenomic samples is a powerful tool to enrich for sequences from rare taxa, which will shed further light on their ecological niches.     

Microbial ecology in the age of genomics and metagenomics: concepts, tools, and recent advances

Microbial ecology examines the diversity and activity of micro-organisms in Earth's biosphere. In the last 20 years, the application of genomics tools have revolutionized microbial ecological studies and drastically expanded our view on the previously underappreciated microbial world. This review first introduces the basic concepts in microbial ecology and the main genomics methods that have been used to examine natural microbial populations and communities. In the ensuing three specific sections, the applications of the genomics in microbial ecological research are highlighted. The first describes the widespread application of multilocus sequence typing and representational difference analysis in studying genetic variation within microbial species. Such investigations have identified that migration, horizontal gene transfer and recombination are common in natural microbial populations and that microbial strains can be highly variable in genome size and gene content. The second section highlights and summarizes the use of four specific genomics methods (phylogenetic analysis of ribosomal RNA, DNA-DNA re-association kinetics, metagenomics, and micro-arrays) in analysing the diversity and potential activity of microbial populations and communities from a variety of terrestrial and aquatic environments. Such analyses have identified many unexpected phylogenetic lineages in viruses, bacteria, archaea, and microbial eukaryotes. Functional analyses of environmental DNA also revealed highly prevalent, but previously unknown, metabolic processes in natural microbial communities. In the third section, the ecological implications of sequenced microbial genomes are briefly discussed. Comparative analyses of prokaryotic genomic sequences suggest the importance of ecology in determining microbial genome size and gene content. The significant variability in genome size and gene content among strains and species of prokaryotes indicate the highly fluid nature of prokaryotic genomes, a result consistent with those from multilocus sequence typing and representational difference analyses. The integration of various levels of ecological analyses coupled to the application and further development of high throughput technologies are accelerating the pace of discovery in microbial ecology.     

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.     

Tracing Patterns of Chloroplast Evolution in Euglenoids: Contributions from Colacium vesiculosum and Strombomonas acuminata (Euglenophyta)

The chloroplast genomes of two photosynthetic euglenoids, Colacium vesiculosum Ehrenberg (128,889 bp), and Strombomonas acuminata (Schmarda) Deflandre (144,167 bp) have been sequenced. These chloroplast genomes in combination with those of Euglena gracilis, Eutreptia viridis, and Eutreptiella gymnastica provide a snapshot of euglenoid chloroplast evolution allowing comparisons of gene content, arrangement, and expansion. The gene content of the five chloroplast genomes is very similar varying only in the presence or absence of, rrn5, roaA, psaI, psaM, rpoA, and two tRNAs. Large gene rearrangements have occurred within the C. vesiculosum and S. acuminata chloroplast genomes. Most of these rearrangements represent repositioning of entire operons rather than single genes. When compared with previously sequenced genomes, C. vesiculosum and S. acuminata chloroplast genomes more closely resemble the E. gracilis chloroplast genome in size of the genome, number of introns, and gene order than they do those of the Eutreptiales. Overall, the chloroplast genomes of these five species show an evolutionary trend toward increased intron number, a decrease in gene density, and substantial rearrangement of gene clusters.     

Selective Whole Genome Amplification for Resequencing Target Microbial Species from Complex Natural Samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host–microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10⁵-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2–9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.     

茎瘤芥根际一株产紫色杆菌素杜擀氏菌的分离鉴定及其基因组分析

【目的】对茎瘤芥根际微生物进行分离鉴定,分析微生物菌群构成,选择具有优良特性的菌株,评估其次级代谢产物合成能力,为茎瘤芥根际微生物多样性和菌种资源的挖掘利用奠定基础。【方法】采集重庆市涪陵区二渡村和邓家村的茎瘤芥根,分离培养根际微生物菌株,通过菌株形态观察和看家基因的序列分析,对菌株进行初步鉴定、归类和保存。选择具有优良性状的菌株,利用Pacbio RS II和Illumina HiSeq平台完成全基因组测序,通过antiSMASH分析评估其次级代谢产物合成潜力,克隆目的基因簇并进行异源表达和产物鉴定。【结果】分离得到256株微生物,初步鉴定120株;从中鉴定了一株产紫色杆菌素的杜擀氏菌BjR8,完成了基因组测序及分析,发现该菌基因组为一条环状染色体,全长7 205 593 bp,GC含量为64.67%,含有6 241个编码基因。生物信息学分析发现基因组含有9个次级代谢产物生物合成基因簇,其中7个基因簇与已知化合物编码基因簇同源性较低,说明该菌具有产生多种新型次级代谢产物的潜力;克隆得到紫色杆菌素生物合成基因簇,并在变铅青链霉菌TK23中完成了异源表达。【结论】从茎瘤芥根际分离得到25...     

A Tribute to Disorder in the Genome of the Bloom-Forming Freshwater Cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is one of the most common bloom-forming cyanobacteria in freshwater ecosystems worldwide. This species produces numerous secondary metabolites, including microcystins, which are harmful to human health. We sequenced the genomes of ten strains of M. aeruginosa in order to explore the genomic basis of their ability to occupy varied environments and proliferate. Our findings show that M. aeruginosa genomes are characterized by having a large open pangenome, and that each genome contains similar proportions of core and flexible genes. By comparing the GC content of each gene to the mean value of the whole genome, we estimated that in each genome, around 11% of the genes seem to result from recent horizontal gene transfer events. Moreover, several large gene clusters resulting from HGT (up to 19 kb) have been found, illustrating the ability of this species to integrate such large DNA molecules. It appeared also that all M. aeruginosa displays a large genomic plasticity, which is characterized by a high proportion of repeat sequences and by low synteny values between the strains. Finally, we identified 13 secondary metabolite gene clusters, including three new putative clusters. When comparing the genomes of Microcystis and Prochlorococcus, one of the dominant picocyanobacteria living in marine ecosystems, our findings show that they are characterized by having almost opposite evolutionary strategies, both of which have led to ecological success in their respective environments.     

Genomes of three methylotrophs from a single niche reveal the genetic and metabolic divergence of the methylophilaceae

The genomes of three representatives of the family Methylophilaceae, Methylotenera mobilis JLW8, Methylotenera versatilis 301, and Methylovorus glucosetrophus SIP3-4, all isolated from a single study site, Lake Washington in Seattle, WA, were completely sequenced. These were compared to each other and to the previously published genomes of Methylobacillus flagellatus KT and an unclassified Methylophilales strain, HTCC2181. Comparative analysis revealed that the core genome of Methylophilaceae may be as small as approximately 600 genes, while the pangenome may be as large as approximately 6,000 genes. Significant divergence between the genomes in terms of both gene content and gene and protein conservation was uncovered, including the varied presence of certain genes involved in methylotrophy. Overall, our data demonstrate that metabolic potentials can vary significantly between different species of Methylophilaceae, including organisms inhabiting the very same environment. These data suggest that genetic divergence among the members of this family may be responsible for their specialized and nonredundant functions in C? cycling, which in turn suggests means for their successful coexistence in their specific ecological niches.     

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.