Differentially Evolved Genes of Salmonella Pathogenicity Islands: Insights into the Mechanism of Host Specificity in Salmonella

Background

The species Salmonella enterica (S. enterica) includes many serovars that cause disease in avian and mammalian hosts. These serovars differ greatly in their host range and their degree of host adaptation. The host specificity of S. enterica serovars appears to be a complex phenomenon governed by multiple factors acting at different stages of the infection process, which makes identification of the cause/s of host specificity solely by experimental methods difficult.

Methodology/Principal Findings

In this study, we have employed a molecular evolution and phylogenetics based approach to identify genes that might play important roles in conferring host specificity to different serovars of S. enterica. These genes are ‘differentially evolved’ in different S. enterica serovars. This list of ‘differentially evolved’ genes includes genes that encode translocon proteins (SipD, SseC and SseD) of both Salmonella pathogenicity islands 1 and 2 encoded type three secretion systems, sptP, which encodes an effector protein that inhibits the mitogen-activated protein kinase pathway of the host cell, and genes which encode effector proteins (SseF and SifA) that are important in placing the Salmonella-containing vacuole in a juxtanuclear position.

Conclusions/Significance

Analysis of known functions of these ‘differentially evolved genes’ indicates that the products of these genes directly interact with the host cell and manipulate its functions and thereby confer host specificity, at least in part, to different serovars of S. enterica that are considered in this study.     

Age-Dependent Enterocyte Invasion and Microcolony Formation by Salmonella

The coordinated action of a variety of virulence factors allows Salmonella enterica to invade epithelial cells and penetrate the mucosal barrier. The influence of the age-dependent maturation of the mucosal barrier for microbial pathogenesis has not been investigated. Here, we analyzed Salmonella infection of neonate mice after oral administration. In contrast to the situation in adult animals, we observed spontaneous colonization, massive invasion of enteroabsorptive cells, intraepithelial proliferation and the formation of large intraepithelial microcolonies. Mucosal translocation was dependent on enterocyte invasion in neonates in the absence of microfold (M) cells. It further resulted in potent innate immune stimulation in the absence of pronounced neutrophil-dominated pathology. Our results identify factors of age-dependent host susceptibility and provide important insight in the early steps of Salmonella infection in vivo. We also present a new small animal model amenable to genetic manipulation of the host for the analysis of the Salmonella enterocyte interaction in vivo.     

The Stability of Complement-Mediated Bactericidal Activity in Human Serum against Salmonella

The complement cascade includes heat-labile proteins and care is required when handling serum in order to preserve its functional integrity. We have previously used a whole human serum bactericidal assay to show that antibody and an intact complement system are required in blood for killing of invasive isolates of Salmonella. The aim of the present study was to evaluate the conditions under which human serum can be stored and manipulated while maintaining complement integrity. Serum bactericidal activity against Salmonella was maintained for a minimum of 35 days when stored at 4°C, eight days at 22°C and 54 hours at 37°C. Up to three freeze-thaw cycles had no effect on the persistence of bactericidal activity and hemolytic complement assays confirmed no effect on complement function. Delay in the separation of serum for up to four days from clotted blood stored at 22°C did not affect bactericidal activity. Dilution of serum resulted in an increased rate of loss of bactericidal activity and so serum should be stored undiluted. These findings indicate that the current guidelines concerning manipulation and storage of human serum to preserve complement integrity and function leave a large margin for safety with regards to bactericidal activity against Salmonella. The study provides a scheme for determining the requirements for serum handling in relation to functional activity of complement in other systems.     

Value Addition in the Efficacy of Conventional Antibiotics by Nisin against Salmonella

Frequent and indiscriminate use of existing battery of antibiotics has led to the development of multi drug resistant (MDR) strains of pathogens. As decreasing the concentration of the antibiotic required to treat Salmonellosis might help in combating the development of resistant strains, the present study was designed to assess the synergistic effects, if any, of nisin, in combination with conventional anti-Salmonella antibiotics against Salmonella enterica serovar Typhimurium. Minimum inhibitory concentrations (MICs) of the selected antimicrobial agents were determined by micro and macro broth dilution assays. In-vitro synergy between the agents was evaluated by radial diffusion assay, fractional inhibitory concentration (FIC) index (checkerboard test) and time-kill assay. Scanning electron microscopy (SEM) was also performed to substantiate the effect of the combinations. In-vivo synergistic efficacy of the combinations selected on the basis of in-vitro results was also evaluated in the murine model, in terms of reduction in the number of Salmonellae in liver, spleen and intestine. Nisin-ampicillin and nisin-EDTA combinations were observed to have additive effects, whereas the combinations of nisin-ceftriaxone and nisin-cefotaxime were found to be highly synergistic against serovar Typhimurium as evident by checkerboard test and time-kill assay. SEM results revealed marked changes on the outer membrane of the bacterial cells treated with various combinations. In-vivo synergy was evident from the larger log unit decreases in all the target organs of mice treated with the combinations than in those treated with drugs alone. This study thus highlights that nisin has the potential to act in conjunction with conventional antibiotics at much lower MICs. These observations seem to be significant, as reducing the therapeutic concentrations of antibiotics may be a valuable strategy for avoiding/reducing the development of emerging antibiotic resistance. Value added potential of nisin in the efficacy of conventional antibiotics may thus be exploited not only against Salmonella but against other Gram-negative infections as well.     

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.     

Proto-Organism Kinetics: Evolutionary Dynamics of Lipid Aggregates with Genes and Metabolism

A synthetic proto-organism could be self-assembled by integrating a lipid proto-container with a proto-metabolic subsystem and a proto-genetic subsystem. This three-component system can use energy and nutrients by means of either redox or photo-chemical reactions, evolve its protogenome by means of template directed replication, and ultimately die. The evolutionary dynamics of the proto-organism depends crucially on the chemical kinetics of its sub-systems and on their interplay. In this work the template replication kinetics is investigated and it is found that the product inhibition inherent in the ligation-like replication process allows for coexistence of unrelated self-replicating proto-genes in the lipid surface layer. The combined catalytic effects from the proto-genes on the metabolic production rates determine the fate of the strain protocell.     

lac Repressor Is an Antivirulence Factor of Salmonella enterica: Its Role in the Evolution of Virulence in Salmonella

The genus Salmonella includes many pathogens of great medical and veterinary importance. Bacteria belonging to this genus are very closely related to those belonging to the genus Escherichia. lacZYA operon and lacI are present in Escherichia coli, but not in Salmonella enterica. It has been proposed that Salmonella has lost lacZYA operon and lacI during evolution. In this study, we have investigated the physiological and evolutionary significance of the absence of lacI in Salmonella enterica. Using murine model of typhoid fever, we show that the expression of LacI causes a remarkable reduction in the virulence of Salmonella enterica. LacI also suppresses the ability of Salmonella enterica to proliferate inside murine macrophages. Microarray analysis revealed that LacI interferes with the expression of virulence genes of Salmonella pathogenicity island 2. This effect was confirmed by RT-PCR and Western blot analysis. Interestingly, we found that SBG0326 of Salmonella bongori is homologous to lacI of Escherichia coli. Salmonella bongori is the only other species of the genus Salmonella and it lacks the virulence genes of Salmonella pathogenicity island 2. Overall, our results demonstrate that LacI is an antivirulence factor of Salmonella enterica and suggest that absence of lacI has facilitated the acquisition of virulence genes of Salmonella pathogenicity island 2 in Salmonella enterica making it a successful systemic pathogen.     

Independent Bottlenecks Characterize Colonization of Systemic Compartments and Gut Lymphoid Tissue by Salmonella

Vaccination represents an important instrument to control typhoid fever in humans and protects mice from lethal infection with mouse pathogenic serovars of Salmonella species. Mixed infections with tagged Salmonella can be used in combination with probabilistic models to describe the dynamics of the infection process. Here we used mixed oral infections with tagged Salmonella strains to identify bottlenecks in the infection process in naïve and vaccinated mice. We established a next generation sequencing based method to characterize the composition of tagged Salmonella strains which offers a fast and reliable method to characterise the composition of genome-tagged Salmonella strains. We show that initial colonization of Salmonella was distinguished by a non-Darwinian selection of few bacteria setting up the infection independently in gut associated lymphoid tissue and systemic compartments. Colonization of Peyer''s patches fuels the sustained spread of bacteria into mesenteric lymph nodes via dendritic cells. In contrast, infection of liver and spleen originated from an independent pool of bacteria. Vaccination only moderately reduced invasion of Peyer''s patches but potently uncoupled bacterial populations present in different systemic compartments. Our data indicate that vaccination differentially skews the capacity of Salmonella to colonize systemic and gut immune compartments and provide a framework for the further dissection of infection dynamics.     

The Small,Slow and Specialized CRISPR and Anti-CRISPR of Escherichia and Salmonella

Prokaryotes thrive in spite of the vast number and diversity of their viruses. This partly results from the evolution of mechanisms to inactivate or silence the action of exogenous DNA. Among these, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) are unique in providing adaptive immunity against elements with high local resemblance to genomes of previously infecting agents. Here, we analyze the CRISPR loci of 51 complete genomes of Escherichia and Salmonella. CRISPR are in two pairs of loci in Escherichia, one single pair in Salmonella, each pair showing a similar turnover rate, repeat sequence and putative linkage to a common set of cas genes. Yet, phylogeny shows that CRISPR and associated cas genes have different evolutionary histories, the latter being frequently exchanged or lost. In our set, one CRISPR pair seems specialized in plasmids often matching genes coding for the replication, conjugation and antirestriction machinery. Strikingly, this pair also matches the cognate cas genes in which case these genes are absent. The unexpectedly high conservation of this anti-CRISPR suggests selection to counteract the invasion of mobile elements containing functional CRISPR/cas systems. There are few spacers in most CRISPR, which rarely match genomes of known phages. Furthermore, we found that strains divergent less than 250 thousand years ago show virtually identical CRISPR. The lack of congruence between cas, CRISPR and the species phylogeny and the slow pace of CRISPR change make CRISPR poor epidemiological markers in enterobacteria. All these observations are at odds with the expectedly abundant and dynamic repertoire of spacers in an immune system aiming at protecting bacteria from phages. Since we observe purifying selection for the maintenance of CRISPR these results suggest that alternative evolutionary roles for CRISPR remain to be uncovered.     

F-Actin Dynamics in Neurospora crassa

This study demonstrates the utility of Lifeact for the investigation of actin dynamics in Neurospora crassa and also represents the first report of simultaneous live-cell imaging of the actin and microtubule cytoskeletons in filamentous fungi. Lifeact is a 17-amino-acid peptide derived from the nonessential Saccharomyces cerevisiae actin-binding protein Abp140p. Fused to green fluorescent protein (GFP) or red fluorescent protein (TagRFP), Lifeact allowed live-cell imaging of actin patches, cables, and rings in N. crassa without interfering with cellular functions. Actin cables and patches localized to sites of active growth during the establishment and maintenance of cell polarity in germ tubes and conidial anastomosis tubes (CATs). Recurrent phases of formation and retrograde movement of complex arrays of actin cables were observed at growing tips of germ tubes and CATs. Two populations of actin patches exhibiting slow and fast movement were distinguished, and rapid (1.2 μm/s) saltatory transport of patches along cables was observed. Actin cables accumulated and subsequently condensed into actin rings associated with septum formation. F-actin organization was markedly different in the tip regions of mature hyphae and in germ tubes. Only mature hyphae displayed a subapical collar of actin patches and a concentration of F-actin within the core of the Spitzenkörper. Coexpression of Lifeact-TagRFP and β-tubulin–GFP revealed distinct but interrelated localization patterns of F-actin and microtubules during the initiation and maintenance of tip growth.Actins are highly conserved proteins found in all eukaryotes and have an enormous variety of cellular roles. The monomeric form (globular actin, or G-actin) can self-assemble, with the aid of numerous actin-binding proteins (ABPs), into microfilaments (filamentous actin, or F-actin), which, together with microtubules, form the two major components of the fungal cytoskeleton. Numerous pharmacological and genetic studies of fungi have demonstrated crucial roles for F-actin in cell polarity, exocytosis, endocytosis, cytokinesis, and organelle movement (6, 7, 20, 34, 35, 51, 52, 59). Phalloidin staining, immunofluorescent labeling, and fluorescent-protein (FP)-based live-cell imaging have revealed three distinct subpopulations of F-actin-containing structures in fungi: patches, cables, and rings (1, 14, 28, 34, 60, 63, 64). Actin patches are associated with the plasma membrane and represent an accumulation of F-actin around endocytic vesicles (3, 26, 57). Actin cables are bundles of actin filaments stabilized with cross-linking proteins, such as tropomyosins and fimbrin, and are assembled by formins at sites of active growth, where they form tracks for myosin V-dependent polarized secretion and organelle transport (10, 16, 17, 27, 38, 47, 48). Cables, unlike patches, are absolutely required for polarized growth in the budding yeast Saccharomyces cerevisiae (34, 38). Contractile actomyosin rings are essential for cytokinesis in budding yeast, whereas in filamentous fungi, actin rings are less well studied but are known to be involved in septum formation (20, 28, 34, 39, 40).Actin cables and patches have been particularly well studied in budding yeast. However, there are likely to be important differences between F-actin architecture and dynamics in budding yeast and those in filamentous fungi, as budding yeasts display only a short period of polarized growth during bud formation, which is followed by isotropic growth over the bud surface (10). Sustained polarized growth during hyphal morphogenesis is a defining feature of filamentous fungi (21), making them attractive models for studying the roles of the actin cytoskeleton in cell polarization, tip growth, and organelle transport.In Neurospora crassa and other filamentous fungi, disruption of the actin cytoskeleton leads to rapid tip swelling, which indicates perturbation of polarized tip growth, demonstrating a critical role for F-actin in targeted secretion to particular sites on the plasma membrane (7, 22, 29, 56). Immunofluorescence studies of N. crassa have shown that F-actin localizes to hyphal tips as “clouds” and “plaques” (7, 54, 59). However, immunolabeling has failed to reveal actin cables in N. crassa and offers limited insights into F-actin dynamics. Live-cell imaging of F-actin architecture and dynamics has not been accomplished in N. crassa, yet it is expected to yield key insights into cell polarization, tip growth, and intracellular transport.We took advantage of a recently developed live-cell imaging probe for F-actin called Lifeact (43). Lifeact is a 17-amino-acid peptide derived from the N terminus of the budding yeast actin-binding protein Abp140 (5, 63) and has recently been demonstrated to be a universal live-cell imaging marker for F-actin in eukaryotes (43). Here, we report the successful application of fluorescent Lifeact fusion constructs for live-cell imaging of F-actin in N. crassa. We constructed two synthetic genes consisting of Lifeact fused to “synthetic” green fluorescent protein (sGFP) (S65T) (henceforth termed GFP) (12) or red fluorescent protein (TagRFP) (33) and expressed these constructs in various N. crassa strains. In all strain backgrounds, fluorescent Lifeact constructs clearly labeled actin patches, cables, and rings and revealed a direct association of F-actin structures with sites of cell polarization and active tip growth. Our results demonstrate the efficacy of Lifeact as a nontoxic live-cell imaging probe in N. crassa.     

Arabidopsis thaliana Genes Encoding Defense Signaling and Recognition Proteins Exhibit Contrasting Evolutionary Dynamics

The interplay between pathogen effectors, their host targets, and cognate recognition proteins provides various opportunities for antagonistic cycles of selection acting on plant and pathogen to achieve or abrogate resistance, respectively. Selection has previously been shown to maintain diversity in plant proteins involved in pathogen recognition and some of their cognate pathogen effectors. We analyzed the signatures of selection on 10 Arabidopsis thaliana genes encoding defense signal transduction proteins in plants, which are potential targets of pathogen effectors. There was insufficient evidence to reject neutral evolution for 6 genes encoding signaling components consistent with these proteins not being targets of effectors and/or indicative of constraints on their ability to coevolve with pathogen effectors. Functional constraints on effector targets may have provided the driving selective force for the evolution of guard proteins. PBS1, a known target of an effector, showed little variation but is known to be monitored by a variable guard protein. Evidence of selection maintaining diversity was present at NPR1, PAD4, and EDS1. Differences in the signatures of selection observed may reflect the numbers of effectors that target a particular protein, the presence or absence of a cognate guard protein, as well as functional constraints imposed by biochemical activities or interactions with plant proteins.     

The Bacterial Cytoskeleton Modulates Motility,Type 3 Secretion,and Colonization in Salmonella

Although there have been great advances in our understanding of the bacterial cytoskeleton, major gaps remain in our knowledge of its importance to virulence. In this study we have explored the contribution of the bacterial cytoskeleton to the ability of Salmonella to express and assemble virulence factors and cause disease. The bacterial actin-like protein MreB polymerises into helical filaments and interacts with other cytoskeletal elements including MreC to control cell-shape. As mreB appears to be an essential gene, we have constructed a viable ΔmreC depletion mutant in Salmonella. Using a broad range of independent biochemical, fluorescence and phenotypic screens we provide evidence that the Salmonella pathogenicity island-1 type three secretion system (SPI1-T3SS) and flagella systems are down-regulated in the absence of MreC. In contrast the SPI-2 T3SS appears to remain functional. The phenotypes have been further validated using a chemical genetic approach to disrupt the functionality of MreB. Although the fitness of ΔmreC is reduced in vivo, we observed that this defect does not completely abrogate the ability of Salmonella to cause disease systemically. By forcing on expression of flagella and SPI-1 T3SS in trans with the master regulators FlhDC and HilA, it is clear that the cytoskeleton is dispensable for the assembly of these structures but essential for their expression. As two-component systems are involved in sensing and adapting to environmental and cell surface signals, we have constructed and screened a panel of such mutants and identified the sensor kinase RcsC as a key phenotypic regulator in ΔmreC. Further genetic analysis revealed the importance of the Rcs two-component system in modulating the expression of these virulence factors. Collectively, these results suggest that expression of virulence genes might be directly coordinated with cytoskeletal integrity, and this regulation is mediated by the two-component system sensor kinase RcsC.     

Evolutionary Dynamics of Insertion Sequences in Relation to the Evolutionary Histories of the Chromosome and Symbiotic Plasmid Genes of Rhizobium etli Populations

Insertion sequences (IS) are mobile genetic elements that are distributed in many prokaryotes. In particular, in the genomes of the symbiotic nitrogen-fixing bacteria collectively known as rhizobia, IS are fairly abundant in plasmids or chromosomal islands that carry the genes needed for symbiosis. Here, we report an analysis of the distribution and genetic conservation of the IS found in the genome of Rhizobium etli CFN42 in a collection of 87 Rhizobium strains belonging to populations with different geographical origins. We used PCR to generate presence/absence profiles of the 39 IS found in R. etli CFN42 and evaluated whether the IS were located in consistent genomic contexts. We found that the IS from the symbiotic plasmid were frequently present in the analyzed strains, whereas the chromosomal IS were observed less frequently. We then examined the evolutionary dynamics of these strains based on a population genetic analysis of two chromosomal housekeeping genes (glyA and dnaB) and three symbiotic sequences (nodC and the two IS elements). Our results indicate that the IS contained within the symbiotic plasmid have a higher degree of genomic context conservation, lower nucleotide diversity and genetic differentiation, and fewer recombination events than the chromosomal housekeeping genes. These results suggest that the R. etli populations diverged recently in Mexico, that the symbiotic plasmid also had a recent origin, and that the IS elements have undergone a process of cyclic infection and expansion.Insertion sequences (IS) are the smallest transposable elements found in prokaryotes (usually less than 3 kb in size). They encode a transposase and may also encode small hypothetical proteins (4, 9). IS are distinguished by their ability to move within prokaryotic replicons, including both the chromosome and plasmids, and to copy themselves into various genomic sites. In this manner, IS elements can inactivate or alter the expression of adjacent genes (4). When IS occur in two or more identical copies within a genome, they can participate in various types of genetic rearrangements (e.g., duplications, inversions, and deletions), suggesting that IS may play an important role in the evolution of their hosts by promoting genomic plasticity (34). Due to these evolutionary dynamics, the diversity and distribution of IS elements differ greatly between taxa and even within strains of the same species (27).Various theories have been proposed to explain the evolution of IS elements in laboratory model strains and environmental bacterial populations (8, 18, 25, 29). Two main hypotheses seek to explain how these elements are maintained over the long term in their host genomes. The first proposes that they occasionally generate beneficial mutations and therefore may represent a selective advantage to their hosts (34). The second suggests that IS elements are genomic parasites that are maintained by their high rate of transposition and might be disseminated among different bacterial lineages by horizontal gene transfer (HGT). Data supporting the second hypothesis have shown that some IS elements may transpose at high rates upon entering a new host (42). After the initial infection, however, purifying selection may continuously remove these elements from the genome. Thus, IS may undergo an infection-expansion-extinction cycle that allows them to remain in different bacterial populations within the gene pool (42). These two hypotheses are not contradictory, and the evolutionary dynamics and distribution of IS may differ greatly depending on several factors, including (most notably) the rate of transposition and HGT, as well as selective pressures, population size, and the host''s habitat (6, 18, 21, 25, 27, 29).In the nitrogen-fixing symbiotic bacteria of the genera Rhizobium, Sinorhizobium, Mesorhizobium, Bradyrhizobium (of the alphaproteobacteria), Cupriavidus, and Burkholderia (of the betaproteobacteria), IS are particularly abundant in symbiotic plasmids (pSym) and symbiotic chromosomal islands (SI) (2, 12, 14, 19, 20, 43). SI and pSym include most of the genes needed to establish symbiosis in the roots of leguminous plants through nodule formation and nitrogen fixation (11). It is generally believed that these elements entered the rhizobial genomes through HGT (39, 40). Comparative genomic analyses have shown that both pSym and SI are highly variable, with the exception of a common set of genes encoding factors critical to nitrogen fixation (nif) and nodulation (nod) (5, 14). SI and pSym have been found to have lower GC contents and different codon usages than the corresponding chromosomal and nonsymbiotic plasmid sequences, suggesting that they were recently acquired by HGT.Some of these symbiotic elements, such as in pSym of Rhizobium etli CFN42 and the SI of Mesorhizobium loti, are conjugative and mobile (30, 32). Genomic analysis of R. etli CFN42 revealed the presence of 39 IS belonging to different families (14); they were found in the chromosome (11 IS); the 371-kb symbiotic plasmid (13 IS); the smaller 192-kb conjugative plasmid, p42a (13 IS); and two other plasmids, p42b and p42c (2 IS). Interestingly, this particular strain shows no evidence of IS disrupting open reading frames (ORFs) or having transpositional activity. However, another 42 incomplete IS may be found in the chromosome, pSym, and the conjugative plasmid; these incomplete sequences are truncated or contain stop codons in their coding sequences.Here, we focused on the dynamics and distribution of IS in different populations of the nitrogen-fixing symbiont R. etli. Since the maintenance of IS in bacterial species might depend on their transpositional activities and horizontal transfer rates, the identification of IS in the same genomic contexts across different strains of the same species could provide new insights into their persistence and divergence over short evolutionary periods. To examine the evolutionary dynamics of IS in natural populations of R. etli, we characterized the distributions, genomic contexts, and sequence variations of IS in isolates of R. etli from three populations with different origins, as well as in some other Rhizobium species. More specifically, we used PCR to generate presence/absence profiles of the 39 IS found in R. etli CFN42 in a collection of 87 strains representing different geographical sites and a gradient of domestication of the bacterial host, the common bean (Phaseolus vulgaris). We also evaluated whether the IS were conserved in the same genomic context relative to their position in R. etli CFN42 and determined the nucleotide sequences of two IS found in most of the isolates. Several population genetic tests applied to these IS, another pSym gene (nodC), and two chromosomal housekeeping genes (dnaB and glyA) suggest that these two IS elements have been inherited vertically and represent recent components of the R. etli gene pool. Finally, the present study strongly suggests that symbiotic plasmids have a recent origin within the R. etli populations.     

Evolutionary Dynamics of HTLV-I

Using mathematical models to describe the in vivo dynamics of HTLV-I infection, an explanation is offered for the slow rate of evolution of HTLV-I relative to HIV-1. In agreement with experimental findings, it is assumed that cell activation is required for successful replication in T helper cells and that HTLV-I induces a significant degree of bystander activation. It is found that the rate of evolution of HTLV-I is limited by the restricted availability of activated uninfected T cells, both at high and low proviral loads. This limits the within-host sequence diversity of HTLV-I and may therefore account for the slow rate of evolution of the virus in the population. Specific differences in the in vivo dynamics of HTLV-I and HIV-1 are identified which may account for the discrepancy in the rate of evolution of these two retroviruses. Received: 7 September 1999 / Accepted: 6 December 1999     

Evolutionary Diversification of New Caledonian Araucaria

New Caledonia is a global biodiversity hotspot. Hypotheses for its biotic richness suggest either that the island is a ‘museum’ for an old Gondwana biota or alternatively it has developed following relatively recent long distance dispersal and in situ radiation. The conifer genus Araucaria (Araucariaceae) comprises 19 species globally with 13 endemic to this island. With a typically Gondwanan distribution, Araucaria is particularly well suited to testing alternative biogeographic hypotheses concerning the origins of New Caledonian biota. We derived phylogenetic estimates using 11 plastid and rDNA ITS2 sequence data for a complete sampling of Araucaria (including multiple accessions of each of the 13 New Caledonian Araucaria species). In addition, we developed a dataset comprising 4 plastid regions for a wider taxon sample to facilitate fossil based molecular dating. Following statistical analyses to identify a credible and internally consistent set of fossil constraints, divergence times estimated using a Bayesian relaxed clock approach were contrasted with geological scenarios to explore the biogeographic history of Araucaria. The phylogenetic data resolve relationships within Araucariaceae and among the main lineages in Araucaria, but provide limited resolution within the monophyletic New Caledonian species group. Divergence time estimates suggest a Late Cretaceous-Cenozoic radiation of extant Araucaria and a Neogene radiation of the New Caledonian lineage. A molecular timescale for the evolution of Araucariaceae supports a relatively recent radiation, and suggests that earlier (pre-Cenozoic) fossil types assigned to Araucaria may have affinities elsewhere in Araucariaceae. While additional data will be required to adequately resolve relationships among the New Caledonian species, their recent origin is consistent with overwater dispersal following Eocene emersion of New Caledonia but is too old to support a single dispersal from Australia to Norfolk Island for the radiation of the Pacific Araucaria sect. Eutacta clade.     

Divergent Evolutionary and Expression Patterns between Lineage Specific New Duplicate Genes and Their Parental Paralogs in Arabidopsis thaliana

Gene duplication is an important mechanism for the origination of functional novelties in organisms. We performed a comparative genome analysis to systematically estimate recent lineage specific gene duplication events in Arabidopsis thaliana and further investigate whether and how these new duplicate genes (NDGs) play a functional role in the evolution and adaption of A. thaliana. We accomplished this using syntenic relationship among four closely related species, A. thaliana, A. lyrata, Capsella rubella and Brassica rapa. We identified 100 NDGs, showing clear origination patterns, whose parental genes are located in syntenic regions and/or have clear orthologs in at least one of three outgroup species. All 100 NDGs were transcribed and under functional constraints, while 24% of the NDGs have differential expression patterns compared to their parental genes. We explored the underlying evolutionary forces of these paralogous pairs through conducting neutrality tests with sequence divergence and polymorphism data. Evolution of about 15% of NDGs appeared to be driven by natural selection. Moreover, we found that 3 NDGs not only altered their expression patterns when compared with parental genes, but also evolved under positive selection. We investigated the underlying mechanisms driving the differential expression of NDGs and their parents, and found a number of NDGs had different cis-elements and methylation patterns from their parental genes. Overall, we demonstrated that NDGs acquired divergent cis-elements and methylation patterns and may experience sub-functionalization or neo-functionalization influencing the evolution and adaption of A. thaliana.