A Novel Computational Method Identifies Intra- and Inter-Species Recombination Events in Staphylococcus aureus and Streptococcus pneumoniae

Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Mutation accumulation and fitness effects in hybridogenetic populations: a comparison to sexual and asexual systems

Background  

Female only unisexual vertebrates that reproduce by hybridogenesis show an unusual genetic composition. They are of hybrid origin but show no recombination between the genomes of their parental species. Instead, the paternal genome is discarded from the germline prior to meiosis, and gametes (eggs only) contain solely unrecombined maternal genomes. Hence hybridogens only transmit maternally inherited mutations. Hybridity is restored each generation by backcrossing with males of the sexual parental species whose genome was eliminated. In contrast, recombining sexual species propagate an intermixed pool of mutations derived from the maternal and paternal parts of the genome. If mutation rates are lower in female gametes than males, it raises the possibility for lower mutation accumulation in a hybridogenetic population, and consequently, higher population fitness than its sexual counterpart.     

How Hot Are Drosophila Hotspots? Examining Recombination Rate Variation and Associations with Nucleotide Diversity,Divergence, and Maternal Age in Drosophila pseudoobscura

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5kb regions of highly elevated recombination rates (10–100x the background rate) called “hotspots.” Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40Kb region of Drosophila pseudoobscura chromosome 2, with one 20kb interval assayed every 5Kb and the adjacent 20kb interval bisected into 10kb pieces. We found that recombination events across the 40kb stretch were relatively evenly distributed across each of the 5kb and 10kb intervals, rather than concentrated in a single 5kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.     

Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes

Background  

The bacterium Salmonella enterica includes a diversity of serotypes that cause disease in humans and different animal species. Some Salmonella serotypes show a broad host range, some are host restricted and exclusively associated with one particular host, and some are associated with one particular host species, but able to cause disease in other host species and are thus considered "host adapted". Five Salmonella genome sequences, representing a broad host range serotype (Typhimurium), two host restricted serotypes (Typhi [two genomes] and Paratyphi) and one host adapted serotype (Choleraesuis) were used to identify core genome genes that show evidence for recombination and positive selection.     

Genome sizes in afrotheria, xenarthra, euarchontoglires, and laurasiatheria

Topical literature and Web site databases provide genome sizes for approximately 4,000 animal species, invertebrates and vertebrates, 330 of which are mammals. We provide the genome size for 67 mammalian species, including 51 never reported before. Knowledge of genome size facilitates sequencing projects. The data presented here encompassed 5 Metatheria (order Didelphimorphia) and 62 Eutheria: 15 Xenarthra, 24 Euarchontoglires (Rodentia), as well as 23 Laurasiatheria (22 Chiroptera and 1 species from Perissodactyla). Already available karyotypes supplement the haploid nuclear DNA contents of the respective species. Thus, we established the first comprehensive set of genome size measurements for 15 Xenarthra species (armadillos) and for 12 house-mouse species; each group was previously represented by only one species. The Xenarthra exhibited much larger genomes than the modal 3 pg DNA known for mammals. Within the genus Mus, genome sizes varied between 2.98 pg and 3.68 pg. The 22 bat species we measured support the low 2.63 pg modal value for Chiroptera. In general, the genomes of Euarchontoglires and Laurasiatheria were found being smaller than those of (Afrotheria and) Xenarthra. Interspecific variation in genome sizes is discussed with particular attention to repetitive elements, which probably promoted the adaptation of extant mammals to their environment.     

Analysis of intragenic recombination at wx in rice: correlation between the molecular and genetic maps within the locus.

In plant genomes as well as other eukaryotic genomes, meiotic recombination does not occur uniformly. At the level of the gene, high recombination frequencies are often observed within genetic loci in maize, but this feature of intragenic recombination is not seen at the csr1 locus in Arabidopsis. These observations suggest that meiotic recombination in plant genomes varies considerably among species. In the present study we investigated meiotic recombination at the wx locus in rice. The mutation sites of wx mutants induced by ethyl methanesulfonate (EMS) treatment or gamma-ray irradiation and a spontaneous wx mutant were physically characterized, and the genetic distances between those wx mutation sites were estimated by pollen analysis. Based on these results, the recombination frequency at the wx locus in rice was estimated as 27.3 kb/cM, which was about 10 times higher than the average for the genome, suggesting that there was a radically different rate of meiotic recombination for intra- and intergenic regions in the rice genome.     

Recombination correlates with synaptonemal complex length and chromatin loop size in bovids—insights into mammalian meiotic chromosomal organization

Homologous chromosomes exchange genetic information through recombination during meiosis, a process that increases genetic diversity, and is fundamental to sexual reproduction. In an attempt to shed light on the dynamics of mammalian recombination and its implications for genome organization, we have studied the recombination characteristics of 112 individuals belonging to 28 different species in the family Bovidae. In particular, we analyzed the distribution of RAD51 and MLH1 foci during the meiotic prophase I that serve, respectively, as proxies for double-strand breaks (DSBs) which form in early stages of meiosis and for crossovers. In addition, synaptonemal complex length and meiotic DNA loop size were estimated to explore how genome organization determines DSBs and crossover patterns. We show that although the number of meiotic DSBs per cell and recombination rates observed vary between individuals of the same species, these are correlated with diploid number as well as with synaptonemal complex and DNA loop sizes. Our results illustrate that genome packaging, DSB frequencies, and crossover rates tend to be correlated, while meiotic chromosomal axis length and DNA loop size are inversely correlated in mammals. Moreover, axis length, DSB frequency, and crossover frequencies all covary, suggesting that these correlations are established in the early stages of meiosis.     

Genome size and recombination in angiosperms: a second look

Despite dramatic differences in genome size--and thus space for recombination to occur--previous workers found no correlation between recombination rate and genome size in flowering plants. Here I re-investigate these claims using phylogenetic comparative methods to test a large data set of recombination data in angiosperms. I show that genome size is significantly correlated with recombination rate across a wide sampling of species and that change in genome size explains a meaningful proportion ( approximately 20%) of variation in recombination rate. I show that the strength of this correlation is comparable with that of several characters previously linked to evolutionary change in recombination rate, but argue that consideration of processes of genome size change likely make the observed correlation a conservative estimate. And finally, although I find that recombination rate increases less than proportionally to change in genome size, several mechanistic and theoretical arguments suggest that this result is not unexpected.     

Relationships Among Powered Flight,Metabolic Rate,Body Mass,Genome Size,and the Retrotransposon Complement of Volant Birds

Avian genomes are of interest because the rapid metabolic rate associated with powered flight requires small cells which constrain genome size. Consequently, flying birds tend to have small genomes relative to other vertebrates such as mammals. It thus stands to reason that flying birds should have smaller genomes than ground-dwelling birds with lower metabolic rates. Small genomes could be condensed but uncompromised in a number of ways, including smaller intergenic intervals, shorter introns, and/or a reduced transposable element (TE) complement. We evaluated genome size in light of the orthologous TE complement among 41 flying (FY) and seven ground-dwelling (GD) bird species to determine if a preponderance of deletions in orthologous TEs might explain the compact genomes of flying birds with high metabolic rates. We measured, across multiple loci in all 48 species, the lengths of 50 contemporary orthologous chicken repeat 1 (CR1, a non-LTR retrotransposon) copies relative to inferred ancestral CR1 sequences. We found genome sizes in GD birds were not different than those in FY birds, but the mean lengths of orthologous CR1 loci were significantly shorter in FY birds than in GD birds. Moreover, we observed a negative correlation between basal metabolic rate and length of orthologous CR1 loci. Finally, we observed positive correlations between body mass and both genome sizes as well as length of orthologous CR1 loci, which we expected given that body mass correlates negatively with metabolic rates. Our results support the contention that metabolism helps shape the avian TE complement and thus indirectly contributes to the compact genomes of birds.     

An update on chloroplast genomes

Plant cells possess two more genomes besides the central nuclear genome: the mitochondrial genome and the chloroplast genome (or plastome). Compared to the gigantic nuclear genome, these organelle genomes are tiny and are present in high copy number. These genomes are less prone to recombination and, therefore, retain signatures of their age to a much better extent than their nuclear counterparts. Thus, they are valuable phylogenetic tools, giving useful information about the relative age and relatedness of the organisms possessing them. Unlike animal cells, mitochondrial genomes of plant cells are characterized by large size, extensive intramolecular recombination and low nucleotide substitution rates and are of limited phylogenetic utility. Chloroplast genomes, on the other hand, show resemblance to animal mitochondrial genomes in terms of phylogenetic utility and are more relevant and useful in case of plants. Conservation in gene order, content and lack of recombination make the plastome an attractive tool for plant phylogenetic studies. Their importance is reflected in the rapid increase in the availability of complete chloroplast genomes in the public databases. This review aims to summarize the progress in chloroplast genome research since its inception and tries to encompass all related aspects. Starting with a brief historical account, it gives a detailed account of the current status of chloroplast genome sequencing and touches upon RNA editing, ycfs, molecular phylogeny, DNA barcoding as well as gene transfer to the nucleus.     

Inheritance and recombination of mitochondrial genomes in plants, fungi and animals

It is generally assumed that mitochondrial genomes are uniparentally transmitted, homoplasmic and nonrecombining. However, these assumptions draw largely from early studies on animal mitochondrial DNA (mtDNA). In this review, we show that plants, animals and fungi are all characterized by episodes of biparental inheritance, recombination among genetically distinct partners, and selfish elements within the mitochondrial genome, but that the extent of these phenomena may vary substantially across taxa. We argue that occasional biparental mitochondrial transmission may allow organisms to achieve the best of both worlds by facilitating mutational clearance but continuing to restrict the spread of selfish genetic elements. We also show that methodological biases and disproportionately allocated study effort are likely to have influenced current estimates of the extent of biparental inheritance, heteroplasmy and recombination in mitochondrial genomes from different taxa. Despite these complications, there do seem to be discernible similarities and differences in transmission dynamics and likelihood of recombination of mtDNA in plant, animal and fungal taxa that should provide an excellent opportunity for comparative investigation of the evolution of mitochondrial genome dynamics.     

Connecting recombination,nucleotide diversity,and species divergence in Drosophila

The association between recombination rate and nucleotide diversity provides compelling evidence for the action of natural selection across much of the Drosophila melanogaster genome. This conclusion is further supported by the lack of association between recombination rate and nucleotide divergence between species. However, studies of other species, including other Drosophila, have not always yielded the same results. Our recent study measured these parameters within the D. pseudoobscura species group using next-generation sequencing and high-throughput genotyping technologies. We documented fine-scale variation in crossover rate within D. pseudoobscura, and we observed that crossover variation was strongly associated with nucleotide diversity only when measured at a fine-scale. We also observed associations between crossover rate and sequence differences between D. pseudoobscura and its close relatives. These latter associations could have been driven in part by mutagenic effects associated with double-strand break repair, but we cannot exclude the possibility that it results primarily from shared ancestral polymorphisms. Overall, this work strongly underscores the importance of scale in testing for associations of recombination rate with other parameters, and it brings us one small step closer to understanding the role of natural selection and other evolutionary forces in shaping divergence among genomes.     

New regulatory regions of Drosophila 412 retrotransposable element generated by recombination

There are no doubts that transposable elements (TEs) have greatly influenced genomes evolution. They have, however, evolved in different ways throughout mammals, plants, and invertebrates. In mammals they have been shown to be widely present but with low transposition activity; in plants they are responsible for large increases in genome size. In Drosophila, despite their low amount, transposition seems to be higher. Therefore, to understand how these elements have evolved in different genomes and how host genomes have proposed to go around them, are major questions on genome evolution. We analyzed sequences of the retrotransposable elements 412 in natural populations of the Drosophila simulans and D. melanogaster species that greatly differ in their amount of TEs. We identified new subfamilies of this element that were the result of mutation or insertion-deletion process, but also of interfamily recombinations. These new elements were well conserved in the D. simulans natural populations. The new regulatory regions produced by recombination could give rise to new elements able to overcome host control of transposition and, thus, become potential genome invaders.     

A reevaluation of rice mitochondrial evolution based on the complete sequence of male-fertile and male-sterile mitochondrial genomes

Plant mitochondrial genomes have features that distinguish them radically from their animal counterparts: a high rate of rearrangement, of uptake and loss of DNA sequences, and an extremely low point mutation rate. Perhaps the most unique structural feature of plant mitochondrial DNAs is the presence of large repeated sequences involved in intramolecular and intermolecular recombination. In addition, rare recombination events can occur across shorter repeats, creating rearrangements that result in aberrant phenotypes, including pollen abortion, which is known as cytoplasmic male sterility (CMS). Using next-generation sequencing, we pyrosequenced two rice (Oryza sativa) mitochondrial genomes that belong to the indica subspecies. One genome is normal, while the other carries the wild abortive-CMS. We find that numerous rearrangements in the rice mitochondrial genome occur even between close cytotypes during rice evolution. Unlike maize (Zea mays), a closely related species also belonging to the grass family, integration of plastid sequences did not play a role in the sequence divergence between rice cytotypes. This study also uncovered an excellent candidate for the wild abortive-CMS-encoding gene; like most of the CMS-associated open reading frames that are known in other species, this candidate was created via a rearrangement, is chimeric in structure, possesses predicted transmembrane domains, and coopted the promoter of a genuine mitochondrial gene. Our data give new insights into rice mitochondrial evolution, correcting previous reports.     

Genomic correlates of recombination rate and its variability across eight recombination maps in the western honey bee (Apis mellifera L.)

Background

Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results

Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions

The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1281-2) contains supplementary material, which is available to authorized users.     

Genome compaction and stability in microsporidian intracellular parasites

Microsporidian genomes are extraordinary among eukaryotes for their extreme reduction: although they are similar in form to other eukaryotic genomes, they are typically smaller than many prokaryotic genomes. At the same time, their rates of sequence evolution are among the highest for eukaryotic organisms. To explore the effects of compaction on nuclear genome evolution, we sequenced 685,000 bp of the Antonospora locustae genome (formerly Nosema locustae) and compared its organization with the recently completed genome of the human parasite Encephalitozoon cuniculi. Despite being very distantly related, the genomes of these two microsporidian species have retained an unexpected degree of synteny: 13% of genes are in the same context, and 30% of the genes were separated by a small number of short rearrangements. Microsporidian genomes are, therefore, paradoxically composed of rapidly evolving sequences harbored within a slowly evolving genome, although these two processes are sometimes considered to be coupled. Microsporidian genomes show that eukaryotic genomes (like genes) do not evolve in a clock-like fashion, and genome stability may result from compaction in addition to a lack of recombination, as has been traditionally thought to occur in bacterial and organelle genomes.     

Large-scale genomic 2D visualization reveals extensive CG-AT skew correlation in bird genomes

Background  

Bird genomes have very different compositional structure compared with other warm-blooded animals. The variation in the base skew rules in the vertebrate genomes remains puzzling, but it must relate somehow to large-scale genome evolution. Current research is inclined to relate base skew with mutations and their fixation. Here we wish to explore base skew correlations in bird genomes, to develop methods for displaying and quantifying such correlations at different scales, and to discuss possible explanations for the peculiarities of the bird genomes in skew correlation.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1672-4) contains supplementary material, which is available to authorized users.     

Low frequency of paleoviral infiltration across the avian phylogeny

Background

Mammalian genomes commonly harbor endogenous viral elements. Due to a lack of comparable genome-scale sequence data, far less is known about endogenous viral elements in avian species, even though their small genomes may enable important insights into the patterns and processes of endogenous viral element evolution.

Results

Through a systematic screening of the genomes of 48 species sampled across the avian phylogeny we reveal that birds harbor a limited number of endogenous viral elements compared to mammals, with only five viral families observed: Retroviridae, Hepadnaviridae, Bornaviridae, Circoviridae, and Parvoviridae. All nonretroviral endogenous viral elements are present at low copy numbers and in few species, with only endogenous hepadnaviruses widely distributed, although these have been purged in some cases. We also provide the first evidence for endogenous bornaviruses and circoviruses in avian genomes, although at very low copy numbers. A comparative analysis of vertebrate genomes revealed a simple linear relationship between endogenous viral element abundance and host genome size, such that the occurrence of endogenous viral elements in bird genomes is 6- to 13-fold less frequent than in mammals.

Conclusions

These results reveal that avian genomes harbor relatively small numbers of endogenous viruses, particularly those derived from RNA viruses, and hence are either less susceptible to viral invasions or purge them more effectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0539-3) contains supplementary material, which is available to authorized users.