Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

A genome-wide identification of genes potentially associated with host specificity of <Emphasis Type="Italic">Brucella</Emphasis> species

Brucella species are facultative intracellular pathogenic α-Proteobacteria that can cause brucellosis in humans and domestic animals. The clinical and veterinary importance of the bacteria has led to well established studies on the molecular mechanisms of Brucella infection of host organisms. However, to date, no genome-wide study has scanned for genes related to the host specificity of Brucella spp. The majority of bacterial genes related to specific environmental adaptations such as host specificity are well-known to have evolved under positive selection pressure. We thus detected signals of positive selection for individual orthologous genes among Brucella genomes and identified genes related to host specificity. We first determined orthologous sets from seven completely sequenced Brucella genomes using the Reciprocal Best Hits (RBH). A maximum likelihood analysis based on the branch-site test was accomplished to examine the presence of positive selection signals, which was subsequently confirmed by phylogenetic analysis. Consequently, 12 out of 2,033 orthologous genes were positively selected by specific Brucella lineages, each of which belongs to a particular animal host. Extensive literature reviews revealed that half of these computationally identified genes are indeed involved in Brucella host specificity. We expect that this genome-wide approach based on positive selection may be reliably used to screen for genes related to environmental adaptation of a particular species and that it will provide a set of appropriate candidate genes.     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

Perspectives of genomic diversification and molecular recombination towards R-gene evolution in plants

Plants are under strong evolutionary pressure in developing new and noble R genes to recognize pathogen avirulence (avr) determinants and bring about stable defense for generation after generations. Duplication, sequence variation by mutation, disparity in the length and structure of leucine rich repeats etc., causes tremendous variations within and among R genes in a plant thereby developing diverse recognitional specificity suitable enough for defense against new pathogens. Recent studies on genome sequencing, diversity and population genetics in different plants have thrown new insights on the molecular evolution of these genes. Tandem and segmental duplication are important factors in R gene abundance as inferred from the distribution of major nucleotide binding site-leucine rich repeats (NBS-LRRs) type R-genes in plant genomes. Likewise, R-gene evolution is also thought to be facilitated by cluster formation thereby causing recombination and sequence exchange and resulting in haplotypic diversity. Population studies have further proven that balancing selection is responsible for the maintenance of allelic diversity in R genes. In this review, we emphasize and discuss on improved perspectives towards the molecular mechanisms and selection pressure responsible for the evolution of NBS-LRR class resistance genes in plants.     

Molecular evolution of the meiotic recombination pathway in mammals

Meiotic recombination shapes evolution and helps to ensure proper chromosome segregation in most species that reproduce sexually. Recombination itself evolves, with species showing considerable divergence in the rate of crossing‐over. However, the genetic basis of this divergence is poorly understood. Recombination events are produced via a complicated, but increasingly well‐described, cellular pathway. We apply a phylogenetic comparative approach to a carefully selected panel of genes involved in the processes leading to crossovers—spanning double‐strand break formation, strand invasion, the crossover/non‐crossover decision, and resolution—to reconstruct the evolution of the recombination pathway in eutherian mammals and identify components of the pathway likely to contribute to divergence between species. Eleven recombination genes, predominantly involved in the stabilization of homologous pairing and the crossover/non‐crossover decision, show evidence of rapid evolution and positive selection across mammals. We highlight TEX11 and associated genes involved in the synaptonemal complex and the early stages of the crossover/non‐crossover decision as candidates for the evolution of recombination rate. Evolutionary comparisons to MLH1 count, a surrogate for the number of crossovers, reveal a positive correlation between genome‐wide recombination rate and the rate of evolution at TEX11 across the mammalian phylogeny. Our results illustrate the power of viewing the evolution of recombination from a pathway perspective.     

Positive selection in AvrP4 avirulence gene homologues across the genus Melampsora

Pathogen genes involved in interactions with their plant hosts are expected to evolve under positive Darwinian selection or balancing selection. In this study a single copy avirulence gene, AvrP4, in the plant pathogen Melampsora lini, was used to investigate the evolution of such a gene across species. Partial translation elongation factor 1-alpha sequences were obtained to establish phylogenetic relationships among the Melampsora species. We amplified AvrP4 homologues from species pathogenic on hosts from different plant families and orders, across the inferred phylogeny. Translations of the AvrP4 sequences revealed a predicted signal peptide and towards the C-terminus of the protein, six identically spaced cysteines were identified in all sequences. Maximum likelihood analysis of synonymous versus non-synonymous substitution rates indicated that positive selection played a role in the evolution of the gene during the diversification of the genus. Fourteen codons under significant positive selection reside in the C-terminal 28 amino acid region, suggesting that this region interacts with host molecules in most sequenced accessions. Selection pressures on the gene may be either due to the pathogenicity or avirulence function of the gene or both.     

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.     

Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent,Treponema pallidum

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.     

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Conjugation genes are common throughout the genus Rickettsia and are transmitted horizontally

Rickettsia are endosymbionts of arthropods, some of which are vectored to vertebrates where they cause disease. Recently, it has been found that some Rickettsia strains harbour conjugative plasmids and others encode some conjugative machinery within the bacterial genome. We investigated the distribution of these conjugation genes in a phylogenetically diverse collection of Rickettsia isolated from arthropods. We found that these genes are common throughout the genus and, in stark contrast to other genes in the genome, conjugation genes are frequently horizontally transmitted between strains. There is no evidence to suggest that these genes are preferentially transferred between phylogenetically related strains, which is surprising given that closely related strains infect similar host species. In addition to detecting patterns of horizontal transmission between diverse Rickettsia species, these findings have implications for the evolution of pathogenicity, the evolution of Rickettsia genomes and the genetic manipulation of intracellular bacteria.     

Characterization and evaluation of an arbitrary primed Polymerase Chain Reaction (PCR) product for the specific detection of Brucella species

Laboratory detection of Brucella is based largely on bacterial isolation and phenotypic characterization. These methods are lengthy and labor-intensive and have been associated with a heightened risk of laboratory-acquired infection. Antibody based indirect detection methods also suffer from limitations in proper diagnosis of the organism. To overcome these problems, nucleic acid amplification has been explored for rapid detection and confirmation of the presence of Brucella spp. PCR-based diagnostics is useful for screening large populations of livestock to identify infected individuals and confirms the presence of the pathogen. Random Amplification of Polymorphic DNA (RAPD) was performed and identified a 1.3 kb PCR fragment specifically amplifiable from DNA isolated from Brucella. A BLAST search revealed no significant homology with the reported sequences from species other than the members of Brucella. The isolated fragment seems to be a part of d-alanine–d-alanine ligase gene in Brucella sp. Translational BLAST revealed certain degree of homology of this sequence with orthologs of this gene reported from other microbial species at the deduced amino acid level. The sequence information was used to develop PCR based assays to detect Brucella sp. from various samples. The minimum detection limit of Brucella from blood and milk samples spiked with Brucella DNA was found to be 1 ng/ml and 10 ng/ml, respectively. In conclusion, we demonstrated that the PCR based detection protocol was successfully used for the detection of Brucella from various organs and spiked samples of diseased sheep. Diagnosis of Brucellosis by PCR based method reported in this study is relatively rapid, specific and simple.     

Genome sequence and comparative analysis of a putative entomopathogenic Serratia isolated from Caenorhabditis briggsae

Background

Entomopathogenic associations between nematodes in the genera Steinernema and Heterorhabdus with their cognate bacteria from the bacterial genera Xenorhabdus and Photorhabdus, respectively, are extensively studied for their potential as biological control agents against invasive insect species. These two highly coevolved associations were results of convergent evolution. Given the natural abundance of bacteria, nematodes and insects, it is surprising that only these two associations with no intermediate forms are widely studied in the entomopathogenic context. Discovering analogous systems involving novel bacterial and nematode species would shed light on the evolutionary processes involved in the transition from free living organisms to obligatory partners in entomopathogenicity.

Results

We report the complete genome sequence of a new member of the enterobacterial genus Serratia that forms a putative entomopathogenic complex with Caenorhabditis briggsae. Analysis of the 5.04 MB chromosomal genome predicts 4599 protein coding genes, seven sets of ribosomal RNA genes, 84 tRNA genes and a 64.8 KB plasmid encoding 74 genes. Comparative genomic analysis with three of the previously sequenced Serratia species, S. marcescens DB11 and S. proteamaculans 568, and Serratia sp. AS12, revealed that these four representatives of the genus share a core set of ~3100 genes and extensive structural conservation. The newly identified species shares a more recent common ancestor with S. marcescens with 99 % sequence identity in rDNA sequence and orthology across 85.6 % of predicted genes. Of the 39 genes/operons implicated in the virulence, symbiosis, recolonization, immune evasion and bioconversion, 21 (53.8 %) were present in Serratia while 33 (84.6 %) and 35 (89 %) were present in Xenorhabdus and Photorhabdus EPN bacteria respectively.

Conclusion

The majority of unique sequences in Serratia sp. SCBI (South African Caenorhabditis briggsae Isolate) are found in ~29 genomic islands of 5 to 65 genes and are enriched in putative functions that are biologically relevant to an entomopathogenic lifestyle, including non-ribosomal peptide synthetases, bacteriocins, fimbrial biogenesis, ushering proteins, toxins, secondary metabolite secretion and multiple drug resistance/efflux systems. By revealing the early stages of adaptation to this lifestyle, the Serratia sp. SCBI genome underscores the fact that in EPN formation the composite end result – killing, bioconversion, cadaver protection and recolonization- can be achieved by dissimilar mechanisms. This genome sequence will enable further study of the evolution of entomopathogenic nematode-bacteria complexes.

Electronic supplementary material

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

Genome size correlates with reproductive fitness in seed beetles

The ultimate cause of genome size (GS) evolution in eukaryotes remains a major and unresolved puzzle in evolutionary biology. Large-scale comparative studies have failed to find consistent correlations between GS and organismal properties, resulting in the ‘C-value paradox’. Current hypotheses for the evolution of GS are based either on the balance between mutational events and drift or on natural selection acting upon standing genetic variation in GS. It is, however, currently very difficult to evaluate the role of selection because within-species studies that relate variation in life-history traits to variation in GS are very rare. Here, we report phylogenetic comparative analyses of GS evolution in seed beetles at two distinct taxonomic scales, which combines replicated estimation of GS with experimental assays of life-history traits and reproductive fitness. GS showed rapid and bidirectional evolution across species, but did not show correlated evolution with any of several indices of the relative importance of genetic drift. Within a single species, GS varied by 4–5% across populations and showed positive correlated evolution with independent estimates of male and female reproductive fitness. Collectively, the phylogenetic pattern of GS diversification across and within species in conjunction with the pattern of correlated evolution between GS and fitness provide novel support for the tenet that natural selection plays a key role in shaping GS evolution.     

Genetic structure and contrasting selection pattern at two major histocompatibility complex genes in wild house mouse populations

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.