Variance and covariances of the numbers of synonymous and nonsynonymous substitutions per site

Sign of APOBEC editing, purifying selection, frameshift, and in-frame nonsense mutations in the microevolution of lumpy skin disease virus

The lumpy skin disease virus (LSDV), which mostly affects ruminants and causes huge-economic loss, was endemic in Africa, caused outbreaks in the Middle East, and was recently detected in Russia, Serbia, Greece, Bulgaria, Kazakhstan, China, Taiwan, Vietnam, Thailand, and India. However, the role of evolutionary drivers such as codon selection, negative/purifying selection, APOBEC editing, and genetic variations such as frameshift and in-frame nonsense mutations in the LSDVs, which cause outbreaks in cattle in various countries, are still largely unknown. In the present study, a frameshift mutation in LSDV035, LSDV019, LSDV134, and LSDV144 genes and in-frame non-sense mutations in LSDV026, LSDV086, LSDV087, LSDV114, LSDV130, LSDV131, LSDV145, LSDV154, LSDV155, LSDV057, and LSDV081 genes were revealed among different clusters. Based on the available complete genome sequences, the prototype wild-type cluster-1.2.1 virus has been found in other than Africa only in India, the wild-type cluster-1.2.2 virus found in Africa were spread outside Africa, and the recombinant viruses spreading only in Asia and Russia. Although LSD viruses circulating in different countries form a specific cluster, the viruses detected in each specific country are distinguished by frameshift and in-frame nonsense mutations. Furthermore, the present study has brought to light that the selection pressure for codons usage bias is mostly exerted by purifying selection, and this process is possibly caused by APOBEC editing. Overall, the present study sheds light on microevolutions in LSDV, expected to help in future studies towards disturbed ORFs, epidemiological diagnostics, attenuation/vaccine reverts, and predicting the evolutionary direction of LSDVs.

The phylodynamics of emerging porcine deltacoronavirus (PDCoV) in Southeast Asia

Porcine deltacoronavirus (PDCoV), a recently emerging pathogen, causes diarrhea in pigs. A previous phylogenetic analysis based on spike genes suggested that PDCoV was divided into 3 different groups, including China, US, and Southeast Asia (SEA). SEA PDCoV, however, is genetically separated from China and US but shares a common ancestor. Its origin and evolution have yet been identified. Herein, phylodynamic analyses based on the full-length genome were performed to investigate the origin and evolution of SEA PDCoV. In the study, 18 full-length genome sequences of SEA PDCoV identified in 2013-2016 together with PDCoV from other regions were used in analyses. The results demonstrated that PDCoV was classified into 2 genogroups including G1 and G2. G1 is further evolved into G1a (China), G1b (US). G2 (SEA) group is further evolved into 3 clades, including SEA-1 (Thailand), SEA-2 (Vietnam), and SEA-2r (Vietnam recombinant) clades. The time to the most recent common ancestor (MRCA) of global PDCoV was estimated to be approximately 1989-1990 and possibly have been circulated in SEA more than a decade. SEA PDCoV is genetically diverse compared to China and US PDCoV. The substitution rate of SEA PDCoV was lower than those of China and US, but the recombination rate of SEA was higher. Recombination analyses revealed 4 potential recombinant events in SEA PDCoV, suggesting that they were derived from the same ancestor of China PDCoV. The SEA-2r subgroup was potentially recombinant between SEA-2 and US strains. In conclusion, major mechanisms driving the complex evolution and genetic diversity of SEA PDCoV were multiple introductions of exotic PDCoV strains followed by recombination. This article is protected by copyright. All rights reserved.

Sequence Diversity, Locus Structure, and Evolutionary History of the SpTransformer Genes in the Sea Urchin Genome

The generation of large immune gene families is often driven by evolutionary pressure exerted on host genomes by their pathogens, which has been described as the immunological arms race. The SpTransformer (SpTrf) gene family from the California purple sea urchin, Strongylocentrotus purpuratus, is upregulated upon immune challenge and encodes the SpTrf proteins that interact with pathogens during an immune response. Native SpTrf proteins bind both bacteria and yeast, and augment phagocytosis of a marine Vibrio, while a recombinant SpTrf protein (rSpTrf-E1) binds a subset of pathogens and a range of pathogen associated molecular patterns. In the sequenced sea urchin genome, there are four SpTrf gene clusters for a total of 17 genes. Here, we report an in-depth analysis of these genes to understand the sequence complexities of this family, its genomic structure, and to derive a putative evolutionary history for the formation of the gene clusters. We report a detailed characterization of gene structure including the intron type and UTRs with conserved transcriptional start sites, the start codon and multiple stop codons, and locations of polyadenylation signals. Phylogenetic and percent mismatch analyses of the genes and the intergenic regions allowed us to predict the last common ancestral SpTrf gene and a theoretical evolutionary history of the gene family. The appearance of the gene clusters from the theoretical ancestral gene may have been driven by multiple duplication and deletion events of regions containing SpTrf genes. Duplications and ectopic insertion events, indels, and point mutations in the exons likely resulted in the extant genes and family structure. This theoretical evolutionary history is consistent with the involvement of these genes in the arms race in responses to pathogens and suggests that the diversification of these genes and their encoded proteins have been selected for based on the survival benefits of pathogen binding and host protection.

Coding-Gene Coevolution Analysis of Rotavirus Proteins: A Bioinformatics and Statistical Approach

Rotavirus remains a major cause of diarrhea in infants and young children worldwide. The permanent emergence of new genotypes puts the potential effectiveness of vaccines under serious question. The distribution of unusual genotypes subject to viral fitness is influenced by interactions among viral proteins. The present work aimed at analyzing the genetic constellation and the coevolution of rotavirus coding genes for the available rotavirus genotypes. Seventy-two full genome sequences of different genetic constellations were analyzed using a genetic algorithm. The results revealed an extensive genome-wide covariance network among the 12 viral proteins. Altogether, the emergence of new genotypes represents a challenge to the outcome and success of vaccination and the coevolutionary analysis of rotavirus proteins may boost efforts to better understand the interaction networks of proteins during viral replication/transcription.

Molecular evolution of genes encoding allergen proteins in the peanuts genus Arachis: Structural and functional implications

Food allergies are severe immune responses to plant and animal products mediated by immunoglobulin E (IgE). Peanuts (Arachis hypogaea L.) are among the top 15 crops that feed the world. However, peanuts is among the “big eight food allergens”, and allergies induced by peanuts are a significant public health problem and a life-threatening concern. Targeted mutation studies in peanuts demonstrate that single residue alterations in these allergen proteins could result in substantial reduction in allergenicity. Knowledge of peanut allergen proteins is confined to the allotetraploid crop and its two progenitors. We explored frequencies and positions of natural mutations in the hyperallergenic homologues Ara h 2 and Ara h 6 in newly generated sequences for 24 Arachis wild species and the crop species, assessed potential mutational impact on allergenicity using immunoblots and structural modeling, and evaluated whether these mutations follow evolutionary trends. We uncovered a wealth of natural mutations, both substitutions and gaps, including the elimination of immunodominant epitopes in some species. These molecular alterations appear to be associated with substantial reductions in allergenicity. The study demonstrated that Ara h 2 and Ara h 6 follow contrasting modes of natural selection and opposing mutational patterns, particularly in epitope regions. Phylogenetic analysis revealed a progressive trend towards immunodominant epitope evolution in Ara h 2. The findings provide valuable insight into the interactions among mutations, protein structure and immune system response, thus presenting a valuable platform for future manipulation of allergens to minimize, treat or eliminate allergenicity. The study strongly encourages exploration of genepools of economically important plants in allergenicity research.

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

Genomic Polymorphism Associated with the Emergence of Virulent Isolates of Mycobacterium bovis in the Nile Delta

Mycobacterium bovis is responsible for bovine tuberculosis in both animals and humans. Despite being one of the most important global zoonotic disease, data related to the ecology and pathogenicity of bovine tuberculosis is scarce, especially in developing countries. In this report, we examined the dynamics of M. bovis transmission among dairy cattle in the Nile Delta of Egypt. Animals belonging to 27 herds from 7 governorates were tested by the Single Intradermal Comparative Skin Tuberculin (SICST), as a preliminary screen for the presence of bovine tuberculosis. Positive SICST reactors were identified in 3% of the animals spread among 40% of the examined herds. Post-mortem examination of slaughtered reactors confirmed the presence of both pulmonary and/or digestive forms of tuberculosis in > 50% of the examined animals. Targeted and whole-genome analysis of M. bovis isolates indicated the emergences of a predominant spoligotype (SB0268) between 2013–2015, suggesting a recent clonal spread of this isolate within the Nile Delta. Surprisingly, 2 isolates belonged to M. bovis BCG group, which are not allowed for animal vaccination in Egypt, while the rest of isolates belonged to the virulent M. bovis clonal complex European 2 present in Latin America and several European countries. Analysis of strain virulence in the murine model of tuberculosis indicated the emergence of a more virulent strain (MBE4) with a specific genotype. More analysis is needed to understand the molecular basis for successful spread of virulent isolates of bovine tuberculosis among animals and to establish genotype/phenotype association.

Mycobacterium tuberculosis Acquires Limited Genetic Diversity in Prolonged Infections, Reactivations and Transmissions Involving Multiple Hosts

Background:Mycobacterium tuberculosis (MTB) has limited ability to acquire variability. Analysis of its microevolution might help us to evaluate the pathways followed to acquire greater infective success. Whole-genome sequencing (WGS) in the analysis of the transmission of MTB has elucidated the magnitude of variability in MTB. Analysis of transmission currently depends on the identification of clusters, according to the threshold of variability (<5 SNPs) between isolates. Objective: We evaluated whether the acquisition of variability in MTB, was more frequent in situations which could favor it, namely intrapatient, prolonged infections or reactivations and interpatient transmissions involving multiple sequential hosts. Methods: We used WGS to analyze the accumulation of variability in sequential isolates from prolonged infections or translations from latency to reactivation. We then measured microevolution in transmission clusters with prolonged transmission time, high number of involved cases, simultaneous involvement of latency and active transmission. Results: Intrapatient and interpatient acquisition of variability was limited, within the ranges expected according to the thresholds of variability proposed, even though bursts of variability were observed. Conclusions: The thresholds of variability proposed for MTB seem to be valid in most circumstances, including those theoretically favoring acquisition of variability. Our data point to multifactorial modulation of microevolution, although further studies are necessary to elucidate the factors underlying this modulation.

OmpA protein sequence-based typing and virulence-associated gene profiles of Pasteurella multocida isolates associated with bovine haemorrhagic septicaemia and porcine pneumonic pasteurellosis in Thailand

Background

Pasteurella multocida is a Gram-negative bacterium that causes economically significant infections of a broad range of animal species. Pneumonic and septicaemic pasteurellosis caused by this bacterium remain important problems in pigs, cattle, and water buffaloes in Thailand. The aim of this study was to characterise the virulence-associated gene profiles and to develop an OmpA molecular typing scheme for classifying 191 bovine and porcine isolates of P. multocida collected between 1989 and 2012 in Thailand using polymerase chain reactions (PCRs), nucleotide sequencing, and sequence and structural bioinformatics analyses.

Results

PCR screening successfully characterised the profiles of 25 virulence-associated genes in all isolates. The gene profiles separated these isolates into bovine and porcine clusters based on eight genes (hgbB, hsf1, tadD, nanH, pfhA, plpE, pmHAS, and tbpA). Phylogenetic analyses of the nucleotide and protein sequences corresponding to the ompA gene, which encodes a major outer membrane surface protein, showed two major bovine and porcine clusters. Structural prediction and analysis of the dN/dS ratio revealed four hypervariable extracellular loops of the OmpA transmembrane domains. These four loops were used to develop an OmpA typing scheme. This scheme classified 186 isolates into five major loop sequence types (LST8, LST12, LST15, LST18, and LST19), consistent with the phylogenetic results. The loop regions of the bovine isolates were predicted to be more antigenic than those of the porcine isolates. Thus, molecular evolution of the OmpA proteins could be used to classify P. multocida isolates into different capsular types, host types, and, possibly, pathogenicity levels.

Conclusions

Together with the virulence-associated gene profiles, the typing reported in this work provides a better understanding of P. multocida virulence. Effective monitoring and potential strain-specific subunit vaccines could be developed based on these loop oligopeptides.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-017-1157-6) contains supplementary material, which is available to authorized users.

Cassava brown streak virus has a rapidly evolving genome: implications for virus speciation, variability, diagnosis and host resistance

Cassava is a major staple food for about 800 million people in the tropics and sub-tropical regions of the world. Production of cassava is significantly hampered by cassava brown streak disease (CBSD), caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). The disease is suppressing cassava yields in eastern Africa at an alarming rate. Previous studies have documented that CBSV is more devastating than UCBSV because it more readily infects both susceptible and tolerant cassava cultivars, resulting in greater yield losses. Using whole genome sequences from NGS data, we produced the first coalescent-based species tree estimate for CBSV and UCBSV. This species framework led to the finding that CBSV has a faster rate of evolution when compared with UCBSV. Furthermore, we have discovered that in CBSV, nonsynonymous substitutions are more predominant than synonymous substitution and occur across the entire genome. All comparative analyses between CBSV and UCBSV presented here suggest that CBSV may be outsmarting the cassava immune system, thus making it more devastating and harder to control.

Integrated sequence and immunology filovirus database at Los Alamos

The Ebola outbreak of 2013–15 infected more than 28 000 people and claimed more lives than all previous filovirus outbreaks combined. Governmental agencies, clinical teams, and the world scientific community pulled together in a multifaceted response ranging from prevention and disease control, to evaluating vaccines and therapeutics in human trials. As this epidemic is finally coming to a close, refocusing on long-term prevention strategies becomes paramount. Given the very real threat of future filovirus outbreaks, and the inherent uncertainty of the next outbreak virus and geographic location, it is prudent to consider the extent and implications of known natural diversity in advancing vaccines and therapeutic approaches. To facilitate such consideration, we have updated and enhanced the content of the filovirus portion of Los Alamos Hemorrhagic Fever Viruses Database. We have integrated and performed baseline analysis of all family Filoviridae sequences deposited into GenBank, with associated immune response data, and metadata, and we have added new computational tools with web-interfaces to assist users with analysis. Here, we (i) describe the main features of updated database, (ii) provide integrated views and some basic analyses summarizing evolutionary patterns as they relate to geo-temporal data captured in the database and (iii) highlight the most conserved regions in the proteome that may be useful for a T cell vaccine strategy.

Database URL: www.hfv.lanl.gov

Molecular evolution of gas cavity in [NiFeSe] hydrogenases resurrected in silico

Oxygen tolerance of selenium-containing [NiFeSe] hydrogenases (Hases) is attributable to the high reducing power of the selenocysteine residue, which sustains the bimetallic Ni–Fe catalytic center in the large subunit. Genes encoding [NiFeSe] Hases are inherited by few sulphate-reducing δ-proteobacteria globally distributed under various anoxic conditions. Ancestral sequences of [NiFeSe] Hases were elucidated and their three-dimensional structures were recreated in silico using homology modelling and molecular dynamic simulation, which suggested that deep gas channels gradually developed in [NiFeSe] Hases under absolute anaerobic conditions, whereas the enzyme remained as a sealed edifice under environmental conditions of a higher oxygen exposure risk. The development of a gas cavity appears to be driven by non-synonymous mutations, which cause subtle conformational changes locally and distantly, even including highly conserved sequence regions.

Changes in pathogenicity and immunogenicity of Mycoplasma mycoides subsp. mycoides strains revealed by comparative genomics analysis

Mycoplasma mycoides subsp. mycoides is the causative agent of contagious bovine pleuropneumonia. A pathogenic strain BEN-1 was isolated from bovine lung and underwent continuous passages in rabbits for 468 generations. During this process, the strain's strong virulence became weak and, gradually, it lost the ability to confer protective immunity in cattle but developed virulence in rabbits. In order to gain insight into the mechanisms behind the reduction in virulence and the loss of immunogenicity, we sequenced five representative strains of the BEN series, including the original strain (BEN-1), the strain generation that first acquired virulence in rabbits (BEN-50), the two vaccine strain generations (BEN-181 and BEN-326), and the strain generation showing the greatest loss of immunogenicity (BEN-468). The gene mutation rate in the four different propagation stages varied greatly, and over half of variations observed in each generation were removed during the propagation process. However, the variation maintained in the BEN-468 generation might contribute to its changes in virulence and immunogenicity. We thus identified 18 genes associated with host adaptation, six genes contributing to virulence in cattle, and 35 genes participating in conferring immunity in cattle. These findings might help us optimize the vaccine to obtain more effective immunization results.

Longitudinal Analysis of Cerebrospinal Fluid and Plasma HIV-1 Envelope Sequences Isolated From a Single Donor with HIV Asymptomatic Neurocognitive Impairment

OBJECTIVE

Combined antiretroviral treatment (cART) has changed the clinical presentation of HIV-associated neurocognitive disorders (HAND) to that of the milder forms of the disease. Asymptomatic neurocognitive impairment (ANI) is now more prevalent and is associated with increased morbidity and mortality risk in HIV-1-infected people. HIV-1 envelope (env) genetic heterogeneity has been detected within the central nervous system (CNS) of individuals with ANI. Changes within env determine co-receptor use, cellular tropism, and neuropathogenesis. We hypothesize that compartmental changes are associated with HIV-1 env C2V4 during ANI and sought to analyze paired HIV-1 env sequences from plasma and cerebrospinal fluid (CSF) of a female subject undergoing long-term cART.

METHODS

Paired plasma and CSF samples were collected at 12-month intervals and HIV-1 env C2V4 was cloned and sequenced.

RESULTS

Phylogenetic analysis of paired samples consistently showed genetic variants unique to the CSF. Phenotypic prediction showed CCR5 (R5) variants for all CSF-derived sequences and showed minor X4 variants (or dual-tropic) in the plasma at later time points. Viral compartmentalization was evident throughout the study, suggesting that the occurrence of distinctive env strains may contribute to the neuropathogenesis of HAND.

CONCLUSIONS

Our study provides new insights about the genetic characteristics within the C2V4 of HIV-1 env that persist after long-term cART and during the course of persistent ANI.

A simple method for estimating the strength of natural selection on overlapping genes

Overlapping genes, where one DNA sequence codes for two proteins with different reading frames, are not uncommon in viruses and cellular organisms. Estimating the direction and strength of natural selection acting on overlapping genes is important for understanding their functionality, origin, evolution, maintenance, and potential interaction. However, the standard methods for estimating synonymous (d_S) and nonsynonymous (d_N) nucleotide substitution rates are inapplicable here because a nucleotide change can be simultaneously synonymous and nonsynonymous when both reading frames involved are considered. We have developed a simple method that can estimate d_N/d_S and test for the action of natural selection in each relevant reading frame of the overlapping genes. Our method is an extension of the modified Nei-Gojobori method previously developed for nonoverlapping genes. We confirmed the reliability of our method using extensive computer simulation. Applying this method, we studied the longest human sense–antisense overlapping gene pair, LRRC8E and ENSG00000214248. Although LRRC8E (leucine-rich repeat containing eight family, member E) is known to regulate cell size, the function of ENSG00000214248 is unknown. Our analysis revealed purifying selection on ENSG00000214248 and suggested that it originated in the common ancestor of bony vertebrates.

A genome-wide identification of genes undergoing recombination and positive selection in Neisseria

Currently, there is particular interest in the molecular mechanisms of adaptive evolution in bacteria. Neisseria is a genus of gram negative bacteria, and there has recently been considerable focus on its two human pathogenic species N. meningitidis and N. gonorrhoeae. Until now, no genome-wide studies have attempted to scan for the genes related to adaptive evolution. For this reason, we selected 18 Neisseria genomes (14 N. meningitidis, 3 N. gonorrhoeae and 1 commensal N. lactamics) to conduct a comparative genome analysis to obtain a comprehensive understanding of the roles of natural selection and homologous recombination throughout the history of adaptive evolution. Among the 1012 core orthologous genes, we identified 635 genes with recombination signals and 10 genes that showed significant evidence of positive selection. Further functional analyses revealed that no functional bias was found in the recombined genes. Positively selected genes are prone to DNA processing and iron uptake, which are essential for the fundamental life cycle. Overall, the results indicate that both recombination and positive selection play crucial roles in the adaptive evolution of Neisseria genomes. The positively selected genes and the corresponding amino acid sites provide us with valuable targets for further research into the detailed mechanisms of adaptive evolution in Neisseria.

Analysis of human rotaviruses from a single location over an 18-year time span suggests that protein coadaption influences gene constellations

Rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses that cause severe gastroenteritis in children. In addition to an error-prone genome replication mechanism, RVs can increase their genetic diversity by reassorting genes during host coinfection. Such exchanges allow RVs to acquire advantageous genes and adapt in the face of selective pressures. However, reassortment may also impose fitness costs if it unlinks genes/proteins that have accumulated compensatory, coadaptive mutations and that operate best when kept together. To better understand human RV evolutionary dynamics, we analyzed the genome sequences of 135 strains (genotype G1/G3/G4-P[8]-I1-C1-R1-A1-N1-T1-E1-H1) that were collected at a single location in Washington, DC, during the years 1974 to 1991. Intragenotypic phylogenetic trees were constructed for each viral gene using the nucleotide sequences, thereby defining novel allele level gene constellations (GCs) and illuminating putative reassortment events. The results showed that RVs with distinct GCs cocirculated during the vast majority of the collection years and that some of these GCs persisted in the community unchanged by reassortment. To investigate the influence of protein coadaptation on GC maintenance, we performed a mutual information-based analysis of the concatenated amino acid sequences and identified an extensive covariance network. Unexpectedly, amino acid covariation was highest between VP4 and VP2, which are structural components of the RV virion that are not thought to directly interact. These results suggest that GCs may be influenced by the selective constraints placed on functionally coadapted, albeit noninteracting, viral proteins. This work raises important questions about mutation-reassortment interplay and its impact on human RV evolution.

IMPORTANCE

Rotaviruses are devastating human pathogens that cause severe diarrhea and kill >450,000 children each year. The virus can evolve by accumulating mutations and by acquiring new genes from other strains via a process called reassortment. However, little is known about the relationship between mutation accumulation and gene reassortment for rotaviruses and how it impacts viral evolution. In this study, we analyzed the genome sequences of human strains found in clinical fecal specimens that were collected at a single hospital over an 18-year time span. We found that many rotaviruses did not reassort their genes but instead maintained them as specific sets (i.e., constellations). By analyzing the encoded proteins, we discovered concurrent amino acid changes among them, which suggests that they are functionally coadapted to operate best when kept together. This study increases our understanding of how rotaviruses evolve over time in the human population.

Diversity and evolution of bacterial twin arginine translocase protein, TatC, reveals a protein secretion system that is evolving to fit its environmental niche

Background

The twin-arginine translocation (Tat) protein export system enables the transport of fully folded proteins across a membrane. This system is composed of two integral membrane proteins belonging to TatA and TatC protein families and in some systems a third component, TatB, a homolog of TatA. TatC participates in substrate protein recognition through its interaction with a twin arginine leader peptide sequence.

Methodology/Principal Findings

The aim of this study was to explore TatC diversity, evolution and sequence conservation in bacteria to identify how TatC is evolving and diversifying in various bacterial phyla. Surveying bacterial genomes revealed that 77% of all species possess one or more tatC loci and half of these classes possessed only tatC and tatA genes. Phylogenetic analysis of diverse TatC homologues showed that they were primarily inherited but identified a small subset of taxonomically unrelated bacteria that exhibited evidence supporting lateral gene transfer within an ecological niche. Examination of bacilli tatCd/tatCy isoform operons identified a number of known and potentially new Tat substrate genes based on their frequent association to tatC loci. Evolutionary analysis of these Bacilli isoforms determined that TatCy was the progenitor of TatCd. A bacterial TatC consensus sequence was determined and highlighted conserved and variable regions within a three dimensional model of the Escherichia coli TatC protein. Comparative analysis between the TatC consensus sequence and Bacilli TatCd/y isoform consensus sequences revealed unique sites that may contribute to isoform substrate specificity or make TatA specific contacts. Synonymous to non-synonymous nucleotide substitution analyses of bacterial tatC homologues determined that tatC sequence variation differs dramatically between various classes and suggests TatC specialization in these species.

Conclusions/Significance

TatC proteins appear to be diversifying within particular bacterial classes and its specialization may be driven by the substrates it transports and the environment of its host.

An evolutionary perspective on Elovl5 fatty acid elongase: comparison of Northern pike and duplicated paralogs from Atlantic salmon

Background

The ability to produce physiologically critical LC-PUFA from dietary fatty acids differs greatly among teleost species, and is dependent on the possession and expression of fatty acyl desaturase and elongase genes. Atlantic salmon, as a result of a recently duplicated genome, have more of these enzymes than other fish. Recent phylogenetic studies show that Northern pike represents the closest extant relative of the preduplicated ancestral salmonid. Here we characterise a pike fatty acyl elongase, elovl5, and compare it to Atlantic salmon elovl5a and elovl5b duplicates.

Results

Phylogenetic analyses show that Atlantic salmon paralogs are evolving symmetrically, and they have been retained in the genome by purifying selection. Heterologous expression in yeast showed that Northern pike Elovl5 activity is indistinguishable from that of the salmon paralogs, efficiently elongating C18 and C20 substrates. However, in contrast to salmon, pike elovl5 was predominantly expressed in brain with negligible expression in liver and intestine.

Conclusions

We suggest that the predominant expression of Elovl5b in salmon liver and Elovl5a in salmon intestine is an adaptation, enabled by genome duplication, to a diet rich in terrestrial invertebrates which are relatively poor in LC-PUFA. Pike have retained an ancestral expression profile which supports the maintenance of PUFA in the brain but, due to a highly piscivorous LC-PUFA-rich diet, is not required in liver and intestine. Thus, the characterisation of elovl5 in Northern pike provides insights into the evolutionary divergence of duplicated genes, and the ecological adaptations of salmonids which have enabled colonisation of nutrient poor freshwaters.

Characterization of the bipartite degron that regulates ubiquitin-independent degradation of thymidylate synthase

TS (thymidylate synthase) is a key enzyme in the de novo biosynthesis of dTMP, and is indispensable for DNA replication. Previous studies have shown that intracellular degradation of the human enzyme [hTS (human thymidylate synthase)] is mediated by the 26S proteasome, and occurs in a ubiquitin-independent manner. Degradation of hTS is governed by a degron that is located at the polypeptide's N-terminus that is capable of promoting the destabilization of heterologous proteins to which it is attached. The hTS degron is bipartite, consisting of two subdomains: an IDR (intrinsically disordered region) that is highly divergent among mammalian species, followed by a conserved amphipathic α-helix (designated hA). In the present report, we have characterized the structure and function of the hTS degron in more detail. We have conducted a bioinformatic analysis of interspecies sequence variation exhibited by the IDR, and find that its hypervariability is not due to diversifying (or positive) selection; rather, it has been subjected to purifying (or negative) selection, although the intensity of such selection is relaxed or weakened compared with that exerted on the rest of the molecule. In addition, we have verified that both subdomains of the hTS degron are required for full activity. Furthermore, their co-operation does not necessitate that they are juxtaposed, but is maintained when they are physically separated. Finally, we have identified a 'cryptic' degron at the C-terminus of hTS, which is activated by the N-terminal degron and appears to function only under certain circumstances; its role in TS metabolism is not known.

Evolution of the HIV-1 env gene in the Rag2-/- gammaC-/- humanized mouse model

Rag2(-/-) gamma(C)(-/-) mice transplanted with human hematopoietic stem cells (DKO-hu-HSC mice) mimic aspects of human infection with human immunodeficiency virus type 1 (HIV-1), including sustained viral replication and CD4(+) T-cell decline. However, the extent of HIV-1 evolution during long-term infection in these humanized mice, a key feature of the natural infection, has not been assessed fully. In this study, we examined the types of genotypic and phenotypic changes in the viral env gene that occur in the viral populations of DKO-hu-HSC mice infected with the CCR5-tropic isolate HIV-1(JRCSF) for up to 44 weeks. The mean rate of divergence of viral populations in mice was similar to that observed in a cohort of humans during a similar period of infection. Many amino acid substitutions were common across mice, including losses of N-linked glycosylation sites and substitutions in the CD4 binding site and in CD4-induced epitopes, indicating common selective pressures between mice. In addition, env variants evolved sensitivity to antibodies directed at V3, suggesting a more open conformation for Env. This phenotypic change was associated with increased CD4 binding efficiency and was attributed to specific amino acid substitutions. In one mouse, env variants emerged that exhibited a CXCR4-tropic phenotype. These sequences were compartmentalized in the mesenteric lymph node. In summary, viral populations in these mice exhibited dynamic behavior that included sequence evolution, compartmentalization, and the appearance of distinct phenotypic changes. Thus, humanized mice offer a useful model for studying evolutionary processes of HIV-1 in a complex host environment.

Genetic diversity and stability of the porA allele as a genetic marker in human Campylobacter infection

The major outer-membrane protein (MOMP) of Campylobacter jejuni and Campylobacter coli, encoded by the porA gene, is extremely genetically diverse. Conformational MOMP epitopes are important in host immunity, and variation in surface-exposed regions probably occurs as a result of positive immune selection during infection. porA diversity has been exploited in genotyping studies using highly discriminatory nucleotide sequences to identify potentially epidemiologically linked cases of human campylobacteriosis. To understand the overall nature and extent of porA diversity and stability in C. jejuni and C. coli we investigated sequences in isolates (n=584) obtained from a defined human population (approx. 600 000) over a defined time period (1 year). A total of 196 distinct porA variants were identified. Regions encoding putative extracellular loops were the most variable in both nucleotide sequence and length. Phylogenetic analysis identified three porA allele clusters that originated in (i) predominantly C. jejuni and a few C. coli, (ii) solely C. jejuni or (iii) predominantly C. coli and a few C. jejuni. The stability of porA within an individual human host was investigated using isolates cultured longitudinally from 64 sporadic cases, 27 of which had prolonged infection lasting between 5 and 98 days (the remainder having illness of normal duration, 0–4 days), and 20 cases from family outbreaks. Evidence of mutation was detected in two patients with prolonged illness. Despite demonstrable positive immune selection in these two unusual cases, the persistence of numerous variants within the population indicated that the porA allele is a valuable tool for use in extended typing schemes.

Selection at a single locus leads to widespread expansion of Toxoplasma gondii lineages that are virulent in mice

Pathogenicity differences among laboratory isolates of the dominant clonal North American and European lineages of Toxoplasma gondii are largely controlled by polymorphisms and expression differences in rhoptry secretory proteins (ROPs). However, the extent to which such differences control virulence in natural isolates of T. gondii, including those from more diverse genetic backgrounds, is uncertain. We elucidated the evolutionary history and functional consequences of diversification in the serine/threonine kinase ROP18, a major virulence determinant in the mouse model. We characterized the extent of sequence polymorphism and the evolutionary forces acting on ROP18 and several antigen-encoding genes within a large collection of natural isolates, comparing them to housekeeping genes and introns. Surprisingly, despite substantial genetic diversity between lineages, we identified just three principal alleles of ROP18, which had very ancient ancestry compared to other sampled loci. Expression and allelic differences between these three alleles of ROP18 accounted for much of the variation in acute mouse virulence among natural isolates. While the avirulent type III allele was the most ancient, intermediate virulent (type II) and highly virulent (type I) lineages predominated and showed evidence of strong selective pressure. Out-group comparison indicated that historical loss of an upstream regulatory element increased ROP18 expression, exposing it to newfound diversifying selection, resulting in greatly enhanced virulence in the mouse model and expansion of new lineages. Population sweeps are evident in many genomes, yet their causes and evolutionary histories are rarely known. Our results establish that up-regulation of expression and selection at ROP18 in T. gondii has resulted in three distinct alleles with widely different levels of acute virulence in the mouse model. Preservation of all three alleles in the wild indicates they are likely adaptations for different niches. Our findings demonstrate that sweeping changes in population structure can result from alterations in a single gene.

Net positive charge of HIV-1 CRF01_AE V3 sequence regulates viral sensitivity to humoral immunity

The third variable region (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope gp120 subunit participates in determination of viral infection coreceptor tropism and host humoral immune responses. Positive charge of the V3 plays a key role in determining viral coreceptor tropism. Here, we examined by bioinformatics, experimental, and protein modelling approaches whether the net positive charge of V3 sequence regulates viral sensitivity to humoral immunity. We chose HIV-1 CRF01_AE strain as a model virus to address the question. Diversity analyses using CRF01_AE V3 sequences from 37 countries during 1984 and 2005 (n = 1361) revealed that reduction in the V3's net positive charge makes V3 less variable due to limited positive selection. Consistently, neutralization assay using CRF01_AE V3 recombinant viruses (n = 30) showed that the reduction in the V3's net positive charge rendered HIV-1 less sensitive to neutralization by the blood anti-V3 antibodies. The especially neutralization resistant V3 sequences were the particular subset of the CCR5-tropic V3 sequences with net positive charges of +2 to +4. Molecular dynamics simulation of the gp120 monomers showed that the V3's net positive charge regulates the V3 configuration. This and reported gp120 structural data predict a less-exposed V3 with a reduced net positive charge in the native gp120 trimer context. Taken together, these data suggest a key role of the V3's net positive charge in the immunological escape and coreceptor tropism evolution of HIV-1 CRF01_AE in vivo. The findings have molecular implications for the adaptive evolution and vaccine design of HIV-1.

Clonal dissemination of extended-spectrum beta-lactamase (ESBL)-producing Klebsiella pneumoniae isolates in a Korean hospital

In this study, we investigated the molecular characteristics of extended-spectrum beta-lactamase (ESBL)-producing Klebsiella pneumoniae isolates that were recovered from an outbreak in a Korean hospital. A new multilocus sequence typing (MLST) scheme for K. pneumoniae based on five housekeeping genes was developed and was evaluated for 43 ESBL-producing isolates from an outbreak as well as 38 surveillance isolates from Korea and also a reference strain. Overall, a total of 37 sequence types (STs) and six clonal complexes (CCs) were identified among the 82 K. pneumoniae isolates. The result of MLST analysis was concordant with that of pulsedfield gel electrophoresis. Most of the outbreak isolates belonged to a certain clone (ST2), and they produced SHV-1 and CTX-M14 enzymes, which was a different feature from that of the K. pneumoniae isolates from other Korean hospitals (ST20 and SHV-12). We also found a different distribution of CCs between ESBL-producing and -nonproducing K. pneumoniae isolates. The MLST method we developed in this study could provide unambiguous and well-resolved data for the epidemiologic study of K. pneumoniae. The outbreak isolates showed different molecular characteristics from the other K. pneumoniae isolates from other Korean hospitals.

The venom gland transcriptome of the Desert Massasauga rattlesnake (Sistrurus catenatus edwardsii): towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea)

Background

Snake venoms are complex mixtures of pharmacologically active proteins and peptides which belong to a small number of superfamilies. Global cataloguing of the venom transcriptome facilitates the identification of new families of toxins as well as helps in understanding the evolution of venom proteomes.

Results

We have constructed a cDNA library of the venom gland of a threatened rattlesnake (a pitviper), Sistrurus catenatus edwardsii (Desert Massasauga), and sequenced 576 ESTs. Our results demonstrate a high abundance of serine proteinase and metalloproteinase transcripts, indicating that the disruption of hemostasis is a principle mechanism of action of the venom. In addition to the transcripts encoding common venom proteins, we detected two varieties of low abundance unique transcripts in the library; these encode for three-finger toxins and a novel toxin possibly generated from the fusion of two genes. We also observed polyadenylated ribosomal RNAs in the venom gland library, an interesting preliminary obsevation of this unusual phenomenon in a reptilian system.

Conclusion

The three-finger toxins are characteristic of most elapid venoms but are rare in viperid venoms. We detected several ESTs encoding this group of toxins in this study. We also observed the presence of a transcript encoding a fused protein of two well-characterized toxins (Kunitz/BPTI and Waprins), and this is the first report of this kind of fusion in a snake toxin transcriptome. We propose that these new venom proteins may have ancillary functions for envenomation. The presence of a fused toxin indicates that in addition to gene duplication and accelerated evolution, exon shuffling or transcriptional splicing may also contribute to generating the diversity of toxins and toxin isoforms observed among snake venoms. The detection of low abundance toxins, as observed in this and other studies, indicates a greater compositional similarity of venoms (though potency will differ) among advanced snakes than has been previously recognized.

Genome inventory and analysis of nuclear hormone receptors in Tetraodon nigroviridis

Nuclear hormone receptors (NRs) form a large superfamily of ligand-activated transcription factors, which regulate genes underlying a wide range of (patho) physiological phenomena. Availability of the full genome sequence of Tetraodon nigroviridis facilitated a genome wide analysis of the NRs in fish genome. Seventy one NRs were found in Tetraodon and were compared with mammalian and fish NR family members. In general, there is a higher representation of NRs in fish genomes compared to mammalian ones. They showed high diversity across classes as observed by phylogenetic analysis. Nucleotide substitution rates show strong negative selection among fish NRs except for pregnane x receptor (PxR), estrogen receptor (ER) and liver x receptor (LxR). This may be attributed to crucial role played by them in metabolism and detoxification of xenobiotic and endobiotic compounds and might have resulted in slight positive selection. Chromosomal mapping and pairwise comparisons of NR distribution in Tetraodon and humans led to the identification of nine syntenic NR regions, of which three are common among fully sequenced vertebrate genomes. Gene structure analysis shows strong conservation of exon structures among orthologoues. Whereas paralogous members show different splicing patterns with intron gain or loss and addition or substitution of exons played a major role in evolution of NR superfamily.

Comprehensive genetic variant discovery in the surfactant protein B gene

Completely penetrant mutations in the surfactant protein B gene (SFTPB) and >75% reduction of SFTPB expression disrupt pulmonary surfactant function and cause neonatal respiratory distress syndrome. To inform studies of genetic regulation of SFTPB expression, we created a catalogue of SFTPB variants by comprehensive resequencing from an unselected, population-based cohort (n = 1,116). We found an excess of low-frequency variation [81 SNPs and five small insertion/deletions (in/dels)]. Despite its small genomic size (9.7 kb), SFTPB was characterized by weak linkage disequilibrium (LD) and high haplotype diversity. Using the HapMap Yoruban and European populations, we identified a recombination hot spot that spans SFTPB, was not detectable in our focused resequencing data, and accounts for weak LD. Using homology-based software tools, we discovered no definitively damaging exonic variants. We conclude that excess low-frequency variation, intragenic recombination and lack of common disruptive exonic variants favor complete resequencing as the optimal approach for genetic association studies to identify regulatory SFTPB variants that cause neonatal respiratory distress syndrome in genetically diverse populations.

A rapid progressor-specific variant clone of simian immunodeficiency virus replicates efficiently in vivo only in the absence of immune responses

A subset of simian immunodeficiency virus (SIV)-infected macaques progresses rapidly to disease with transient SIV-specific immune responses and high viral loads. Unique SIV variants with convergent Env mutations evolve in these rapid progressor (RP) macaques. To address the pathogenic significance of RP-specific variants, we generated infectious molecular clones from the terminal-phase plasma of an RP macaque. Inoculation of macaques with a representative clone, SIVsmH635FC, resulted in a persistent viremia, comparable to that produced by pathogenic SIVsmE543-3, and a chronic disease with progressive loss of CD4(+) T cells. However, SIVsmH635FC did not reproduce the rapid-disease phenomenon. Molecular analyses of viruses from these macaques revealed rapid reversion to the wild-type SIVsmE543-3 sequence at two RP-specific sites and slower reversion at another three sites. SIVsmH635FC infection was not sufficient to cause rapid progression even following coinoculation with SIVsmE543-3, despite acute depletion of memory CD4(+) T cells. SIVsmH635FC competed efficiently during primary infection in the coinoculated macaques, but SIVsmE543-3 predominated after the development of SIV-specific immune responses. These data suggest that the replication fitness of the RP variant was similar to that of SIVsmE543-3 in a naïve host; however, SIVsmH635FC was at a disadvantage following the development of SIV-specific immune responses. Consistent with these findings, neutralization assays revealed that SIVsmH635FC was highly sensitive to neutralization but that the parental SIVsmE543-3 strain was highly resistant. This study suggests that the evolution of RP-specific variants is the result of replication in a severely immunocompromised host, rather than the direct cause of rapid progression.

Selecton 2007: advanced models for detecting positive and purifying selection using a Bayesian inference approach

Biologically significant sites in a protein may be identified by contrasting the rates of synonymous (K_s) and non-synonymous (K_a) substitutions. This enables the inference of site-specific positive Darwinian selection and purifying selection. We present here Selecton version 2.2 (http://selecton.bioinfo.tau.ac.il), a web server which automatically calculates the ratio between K_a and K_s (ω) at each site of the protein. This ratio is graphically displayed on each site using a color-coding scheme, indicating either positive selection, purifying selection or lack of selection. Selecton implements an assembly of different evolutionary models, which allow for statistical testing of the hypothesis that a protein has undergone positive selection. Specifically, the recently developed mechanistic-empirical model is introduced, which takes into account the physicochemical properties of amino acids. Advanced options were introduced to allow maximal fine tuning of the server to the user's specific needs, including calculation of statistical support of the ω values, an advanced graphic display of the protein's 3-dimensional structure, use of different genetic codes and inputting of a pre-built phylogenetic tree. Selecton version 2.2 is an effective, user-friendly and freely available web server which implements up-to-date methods for computing site-specific selection forces, and the visualization of these forces on the protein's sequence and structure.

High levels of subgenomic HCV plasma RNA in immunosilent infections

A genetic analysis of hepatitis C virus (HCV) in rare blood donors who remained HCV seronegative despite long-term high-level viremia revealed the chronic presence of HCV genomes with large in frame deletions in their structural genes. Full-length HCV genomes were only detected as minority variants. In one immunodeficiency virus (HIV) co-infected donor the truncated HCV genome transiently decreased in frequency concomitant with delayed seroconversion and re-emerged following partial seroreversion. The long-term production of heavily truncated HCV genomes in vivo suggests that these viruses retained the necessary elements for RNA replication while the deleted structural functions necessary for their spread in vivo was provided in trans by wild-type helper virus in co-infected cells. The absence of immunological pressure and a high viral load may therefore promote the emergence of truncated HCV subgenomic replicons in vivo.

Fitness landscape of human immunodeficiency virus type 1 protease quasispecies

Here we show, at a high resolution (1%), the human immunodeficiency virus type 1 (HIV-1) protease gene quasispecies landscape from three infected naïve individuals. A huge range of genetic configurations was found (67%, 71%, and 80% of the nucleotide clones from the three individuals, respectively, were different), and these configurations created a dense net that linked different parts of the viral population. Similarly, a vast diversity of different protease activities was also found. Importantly, 65% of the analyzed enzymes had detectable protease activity, and 11% of the minority individual variants showed similar or better fitness than the master (most abundant) enzyme, suggesting that the viral complexity in this genomic region does not exclusively depend on the enzyme's catalytic efficiency. Several high-fitness minority variants had only one substitution compared to the master sequence, supporting the possibility that the rugged HIV-1 protease quasispecies fitness landscape may be formed by a continuous network that can be traversed by single mutational steps without passing through defective or less-adapted proteins.

Nature of allelic sequence polymorphism at the KIR3DL3 locus

KIR3DL3 is a framework gene of the Leukocyte Receptor Complex, present in all individuals and haplotypes analysed to date. We describe 17 novel KIR3DL3 alleles, including seven single nucleotide polymorphic (SNP) positions within the coding region. Sequence variation within introns included a VNTR within intron 1. As KIR3DL3 mRNA is known to be expressed in decidual NK cells, we investigated the impact of KIR3DL3 allelic variation on pre-eclampsia. No statistical difference in allele frequency or polymorphism was observed between pre-eclampsia patient and control cohorts. Linkage disequilibrium (LD) analysis of exonic SNPs suggested that recombination may be a mechanism of generating sequence diversity within KIR3DL3. A potential recombination hotspot was located within intron 5. A strong LD was detected between polymorphism in exon 6 of KIR3DL3 and the KIR gene -2DL3 or -2DS2 loci, which define the centromeric end of two main haplotypes (A and B) of the KIR cluster. Comparison of primate KIR sequences indicated that the Ig domains of KIR3DL3 are highly conserved between chimpanzee, gorilla and humans. Investigation of KIR3DL3 dN/dS ratios indicated a greater level of synonymous mutations consistent with purifying selection, although positive selection was detected acting on two sites within the stem region.

Evolution of the uniquely adaptable lentiviral envelope in a natural reservoir host

Background

The ability of emerging pathogens to infect new species is likely related to the diversity of pathogen variants present in existing reservoirs and their degree of genomic plasticity, which determines their ability to adapt to new environments. Certain simian immunodeficiency viruses (SIVcpz, SIVsm) have demonstrated tremendous success in infecting new species, including humans, resulting in the HIV-1 and HIV-2 epidemics. Although SIV diversification has been studied on a population level, the essential substrates for cross-species transmission, namely SIV sequence diversity and the types and extent of viral diversification present in individual reservoir animals have not been elucidated. To characterize this intra-host SIV diversity, we performed sequence analyses of clonal viral envelope (env) V1V2 and gag p27 variants present in individual SIVsm-infected sooty mangabeys over time.

Results

SIVsm demonstrated extensive intra-animal V1V2 length variation and amino acid diversity (le38%), and continual variation in V1V2 N-linked glycosylation consensus sequence frequency and location. Positive selection was the predominant evolutionary force. Temporal sequence shifts suggested continual selection, likely due to evolving antibody responses. In contrast, gag p27 was predominantly under purifying selection. SIVsm V1V2 sequence diversification is at least as great as that in HIV-1 infected humans, indicating that extensive viral diversification in and of itself does not inevitably lead to AIDS.

Conclusion

Positive diversifying selection in this natural reservoir host is the engine that has driven the evolution of the uniquely adaptable SIV/HIV envelope protein. These studies emphasize the importance of retroviral diversification within individual host reservoir animals as a critical substrate in facilitating cross-species transmission.

SIVsm quasispecies adaptation to a new simian host

Despite the potential for infectious agents harbored by other species to become emerging human pathogens, little is known about why some agents establish successful cross-species transmission, while others do not. The simian immunodeficiency viruses (SIVs), certain variants of which gave rise to the human HIV-1 and HIV-2 epidemics, have demonstrated tremendous success in infecting new host species, both simian and human. SIVsm from sooty mangabeys appears to have infected humans on several occasions, and was readily transmitted to nonnatural Asian macaque species, providing animal models of AIDS. Here we describe the first in-depth analysis of the tremendous SIVsm quasispecies sequence variation harbored by individual sooty mangabeys, and how this diverse quasispecies adapts to two different host species—new nonnatural rhesus macaque hosts and natural sooty mangabey hosts. Viral adaptation to rhesus macaques was associated with the immediate amplification of a phylogenetically related subset of envelope (env) variants. These variants contained a shorter variable region 1 loop and lacked two specific glycosylation sites, which may be selected for during acute infection. In contrast, transfer of SIVsm to its natural host did not subject the quasispecies to any significant selective pressures or bottleneck. After 100 d postinfection, variants more closely representative of the source inoculum reemerged in the macaques. This study describes an approach for elucidating how pathogens adapt to new host species, and highlights the particular importance of SIVsm env diversity in enabling cross-species transmission. The replicative advantage of a subset of SIVsm variants in macaques may be related to features of target cells or receptors that are specific to the new host environment, and may involve CD4-independent engagement of a viral coreceptor conserved among primates.

Why do some infectious agents establish successful cross-species transmission while others do not? Despite the clear potential for diseases harbored by animals to become emerging human pathogens, this question remains unanswered. Certain simian immunodeficiency viruses (SIVs) responsible for the human HIV-1 and HIV-2 epidemics have succeeded in infecting new host species, including humans. This study provides clues to how an SIV adapts to a new host in an experimental cross-species transmission. Indeed, many emerging diseases are caused by highly mutation-prone RNA viruses like SIV, which exist not as a single species, but rather as a population of genetic variants within a single infection. The presence of numerous viral variants in an infected animal increases the chance that variants with the ability to enter into or multiply in a new host species are present. This study describes how an SIV population from a natural reservoir host, the sooty mangabey, adapts to a new monkey species, the rhesus macaque. A limited subset of SIV variants containing unique viral surface proteins appears well suited to multiply in the new host. This study documents how viral variation facilitates cross-species transmission, and highlights the particular importance of immunodeficiency virus envelope variants in infecting new hosts.

Functional constraints of influenza A virus epitopes limit escape from cytotoxic T lymphocytes

Viruses can exploit a variety of strategies to evade immune surveillance by cytotoxic T lymphocytes (CTL), including the acquisition of mutations in CTL epitopes. Also for influenza A viruses a number of amino acid substitutions in the nucleoprotein (NP) have been associated with escape from CTL. However, other previously identified influenza A virus CTL epitopes are highly conserved, including the immunodominant HLA-A*0201-restricted epitope from the matrix protein, M1(58-66). We hypothesized that functional constraints were responsible for the conserved nature of influenza A virus CTL epitopes, limiting escape from CTL. To assess the impact of amino acid substitutions in conserved epitopes on viral fitness and recognition by specific CTL, we performed a mutational analysis of CTL epitopes. Both alanine replacements and more conservative substitutions were introduced at various positions of different influenza A virus CTL epitopes. Alanine replacements for each of the nine amino acids of the M1(58-66) epitope were tolerated to various extents, except for the anchor residue at the second position. Substitution of anchor residues in other influenza A virus CTL epitopes also affected viral fitness. Viable mutant viruses were used in CTL recognition experiments. The results are discussed in the light of the possibility of influenza viruses to escape from specific CTL. It was speculated that functional constraints limit variation in certain epitopes, especially at anchor residues, explaining the conserved nature of these epitopes.

Direct inoculation of simian immunodeficiency virus from sooty mangabeys in black mangabeys (Lophocebus aterrimus): first evidence of AIDS in a heterologous African species and different pathologic outcomes of experimental infection

A unique opportunity for the study of the role of serial passage and cross-species transmission was offered by a series of experiments carried out at the Tulane National Primate Research Center in 1990. To develop an animal model for leprosy, three black mangabeys (BkMs) (Lophocebus aterrimus) were inoculated with lepromatous tissue that had been serially passaged in four sooty mangabeys (SMs) (Cercocebus atys). All three BkMs became infected with simian immunodeficiency virus from SMs (SIVsm) by day 30 postinoculation (p.i.) with lepromatous tissue. One (BkMG140) died 2 years p.i. from causes unrelated to SIV, one (BkMG139) survived for 10 years, whereas the third (BkMG138) was euthanized with AIDS after 5 years. Histopathology revealed a high number of giant cells in tissues from BkMG138, but no SIV-related lesions were found in the remaining two BkMs. Four-color immunofluorescence revealed high levels of SIVsm associated with both giant cells and T lymphocytes in BkMG138 and no detectable SIV in the remaining two. Serum viral load (VL) showed a significant increase (>1 log) during the late stage of the disease in BkMG138, as opposed to a continuous decline in VL in the remaining two BkMs. With the progression to AIDS, neopterin levels increased in BkMG138. This study took on new significance when phylogenetic analysis unexpectedly showed that all four serially inoculated SMs were infected with different SIVsm lineages prior to the beginning of the experiment. Furthermore, the strain infecting the BkMs originated from the last SM in the series. Therefore, the virus infecting BkMs has not been serially passaged. In conclusion, we present the first compelling evidence that direct cross-species transmission of SIV may induce AIDS in heterologous African nonhuman primate (NHP) species. The results showed that cross-species-transmitted SIVsm was well controlled in two of three BkMs for 2 and 10 years, respectively. Finally, this case of AIDS in an African monkey suggests that the dogma of SIV nonpathogenicity in African NHP hosts should be reconsidered.

Duplicated genes evolve slower than singletons despite the initial rate increase

BACKGROUND

Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection.

RESULTS

Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins.

CONCLUSIONS

These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

Molecular evolution of the dotA gene in Legionella pneumophila

The molecular evolution of dotA, which is related to the virulence of Legionella pneumophila, was investigated by comparing the sequences of 15 reference strains (serogroups 1 to 15). It was found that dotA has a complex mosaic structure. The whole dotA gene of Legionella pneumophila subsp. pneumophila serogroups 2, 6, and 12 has been transferred from Legionella pneumophila subsp. fraseri. A discrepancy was found between the trees inferred from the nucleotide and deduced amino acid sequences of dotA, which suggests that multiple hits, resulting in synonymous substitutions, have occurred. Gene phylogenies inferred from three different segments (the 5'-end region, the central, large periplasmic domain, and the 3'-end region) showed impressively dissimilar topologies. This was concordant with the sequence polymorphisms, indicating that each region has experienced an independent evolutionary history, and was evident even within the same domain of each strain. For example, the PP2 domain was found to have a heterogeneous structure, which led us hypothesize that the dotA gene of L. pneumophila may have originated from two or more different sources. Comparisons of synonymous and nonsynonymous substitutions demonstrated that the PP2 domain has been under strong selective pressure with respect to amino acid change. Split decomposition analysis also supported the intragenic recombination of dotA. Multiple recombinational exchange within the dotA gene, encoding an integral cytoplasmic membrane protein that is secreted, probably provided increased fitness in certain environmental niches, such as within a particular biofilm community or species of amoebae.

Unique pattern of convergent envelope evolution in simian immunodeficiency virus-infected rapid progressor macaques: association with CD4-independent usage of CCR5

The rate of disease development in simian immunodeficiency virus (SIV) infection of macaques varies considerably among individual macaques. While the majority of macaques inoculated with pathogenic SIV develop AIDS within a period of 1 to 2 years, a minority exhibit a rapid disease course characterized by absence or transience of humoral and cellular immune responses and high levels of virus replication with widespread dissemination of SIV in macrophages and multinucleated giant cells. The goal of this study was to examine viral evolution in three SIVsmE543-3-inoculated rapid progressors to determine the contribution of viral evolution to the development of rapid disease and the effect of the absence of immune pressure upon viral evolution. PCR was used to amplify and clone the entire SIV genome from tissues collected at necropsy, and the course of viral evolution was assessed by env sequences cloned from sequential plasma samples of one rapid progressor (RP) macaque. The majority of sequence changes in RP macaques occurred in the envelope gene. Substitutions were observed in all three animals at specific conserved residues in envelope, including loss of a glycosylation site in V1/V2, a D-to-N/V substitution in a highly conserved GDPE motif, and a P-to-V/H/T substitution in the V3 loop analog. A cell-cell fusion assay revealed that representative env clones utilized CCR5 as a coreceptor, independent of CD4. The selection of specific substitutions in envelope in RP macaques suggests novel selection pressures on virus in such animals and suggests that viral variants that evolve in these animals may play a role in disease progression.

The reverse transcriptase sequence of human immunodeficiency virus type 1 is under positive evolutionary selection within the central nervous system

The human immunodeficiency virus type 1 (HIV-1) enters the central nervous system (CNS) during the acute phase of infection and causes AIDS-related encephalitis and dementia in 30% of individuals. Previous studies show that HIV-1 sequences derived from the CNS of infected patients, including the sequence encoding reverse transcriptase (RT), are genetically distinct from sequences in other tissues. The hypothesis of the current study is that the RT sequence of HIV-1 is under positive selection within the CNS. Multiple alignments of non-CNS-derived and CNS-derived HIV-1 RT sequences were constructed using the ClustalW 1.8 program. The multiple alignments were analyzed with the Synonymous/Nonsynonymous Analysis Program. Codon positions 122-125, 135-149, and 166-212 of the CNS-derived RT sequences underwent a greater accumulation of nonsynonymous than synonymous substitutions, which was markedly different from the analysis results of the non-CNS-derived RT sequences. These residues are located in the finger and palm subdomains of the RT protein structure, which encodes the polymerase active site. The analysis of CNS-derived partial-length RT sequences that encompass these regions yielded similar results. A comparison of CNS-derived RT sequences to a non-CNS-derived RT consensus sequence revealed that a majority of the nonsynonymous substitutions resulted in a specific amino acid replacement. These results indicate that reverse transcriptase is under positive selection within the CNS. The amino acid replacements were visualized on a three-dimensional structure of HIV-1 RT using the Sybyl software suite. The protein structure analysis revealed that the amino acid replacements observed among the CNS-derived sequences occurred in areas of known structural and functional significance.

A comparative genomic method for computational identification of prokaryotic translation initiation sites

The ever growing number of completely sequenced prokaryotic genomes facilitates cross-species comparisons by genomic annotation algorithms. This paper introduces a new probabilistic framework for comparative genomic analysis and demonstrates its utility in the context of improving the accuracy of prokaryotic gene start site detection. Our frame work employs a product hidden Markov model (PROD-HMM) with state architecture to model the species-specific trinucleotide frequency patterns in sequences immediately upstream and downstream of a translation start site and to detect the contrasting non-synonymous (amino acid changing) and synonymous (silent) substitution rates that differentiate prokaryotic coding from intergenic regions. Depending on the intricacy of the features modeled by the hidden state architecture, intergenic, regulatory, promoter and coding regions can be delimited by this method. The new system is evaluated using a preliminary set of orthologous Pyrococcus gene pairs, for which it demonstrates an improved accuracy of detection. Its robustness is confirmed by analysis with cross-validation of an experimentally verified set of Escherichia coli K-12 and Salmonella thyphimurium LT2 orthologs. The novel architecture has a number of attractive features that distinguish it from previous comparative models such as pair-HMMs.

Population genetic structure of Legionella pneumophila inferred from RNA polymerase gene (rpoB) and DotA gene (dotA) sequences

The population structure of Legionella pneumophila was studied by using partial RNA polymerase gene (rpoB) and DotA gene (dotA) sequences. Trees inferred from rpoB sequences showed that two subspecies of L. pneumophila, Legionella pneumophila subsp. pneumophila and Legionella pneumophila subsp. fraseri, were clearly separated genetically. In both rpoB and dotA trees, 79 Korean isolates used in this study constituted six clonal populations, four of which (designated subgroups P-I to P-IV) were identified in L. pneumophila subsp. pneumophila and two of which (designated subgroups F-I and F-II) were identified in L. pneumophila subsp. fraseri. Although the relationships among subgroups were not identical, such subgrouping was congruent between the rpoB and dotA trees. Type strains of several serogroups did not belong to any subgroup, presumably because isolates similar to these strains were not present among our local sample of the population. There was evidence that horizontal gene transfer or recombination had occurred within L. pneumophila. Contrary to the phylogeny from rpoB and the taxonomic context, subgroups P-III and P-IV of L. pneumophila subsp. pneumophila proved to be closely related to those of L. pneumophila subsp. fraseri or showed a distinct clustering in the dotA tree. It can be inferred that dotA of subgroups P-III and P-IV has been transferred horizontally from other subspecies. The diverse distribution of serogroup 1 strains through the gene trees suggests that surface antigen-coding genes that determine serogroup can be exchanged. Thus, it can be inferred that genetic recombination has been important in the evolution of L. pneumophila.

Identification of human immunodeficiency virus type 1 subtype C Gag-, Tat-, Rev-, and Nef-specific elispot-based cytotoxic T-lymphocyte responses for AIDS vaccine design

The most severe human immunodeficiency virus type 1 (HIV-1) epidemic is occurring in southern Africa. It is caused by HIV-1 subtype C (HIV-1C). In this study we present the identification and analysis of cumulative cytotoxic T-lymphocyte (CTL) responses in the southern African country of Botswana. CTLs were shown to be an important component of the immune response to control HIV-1 infection. The definition of optimal and dominant epitopes across the HIV-1C genome that are targeted by CTL is critical for vaccine design. The characteristics of the predominant virus that causes the HIV-1 epidemic in a certain geographic area and also the genetic background of the population, through the distribution of common HLA class I alleles, might impact dominant CTL responses in the vaccinee and in the general population. The enzyme-linked immunospot (Elispot) gamma interferon assay has recently been shown to be a reliable tool to map optimal CTL epitopes, correlating well with other methods, such as intracellular staining, tetramer staining, and the classical chromium release assay. Using Elispot with overlapping synthetic peptides across Gag, Tat, Rev, and Nef, we analyzed HIV-1C-specific CTL responses of HIV-1-infected blood donors. Profiles of cumulative Elispot-based CTL responses combined with diversity and sequence consensus data provide an additional characterization of immunodominant regions across the HIV-1C genome. Results of the study suggest that the construction of a poly-epitope subtype-specific HIV-1 vaccine that includes multiple copies of immunodominant CTL epitopes across the viral genome, derived from predominant HIV-1 viruses, might be a logical approach to the design of a vaccine against AIDS.

Positive Darwinian selection after gene duplication in primate ribonuclease genes

Evolutionary mechanisms of origins of new gene function have been a subject of long-standing debate. Here we report a convincing case in which positive Darwinian selection operated at the molecular level during the evolution of novel function by gene duplication. The genes for eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) in primates belong to the ribonuclease gene family, and the ECP gene, whose product has an anti-pathogen function not displayed by EDN, was generated by duplication of the EDN gene about 31 million years ago. Using inferred nucleotide sequences of ancestral organisms, we showed that the rate of nonsynonymous nucleotide substitution was significantly higher than that of synonymous substitution for the ECP gene. This strongly suggests that positive Darwinian selection operated in the early stage of evolution of the ECP gene. It was also found that the number of arginine residues increased substantially in a short period of evolutionary time after gene duplication, and these amino acid changes probably produced the novel anti-pathogen function of ECP.

Evolution of the Borrelia burgdorferi outer surface protein OspC

The genes coding for outer surface protein OspC from 22 Borrelia burgdorferi strains isolated from patients with Lyme borreliosis were cloned and sequenced. For reference purposes, the 16S rRNA genes from 17 of these strains were sequenced after being cloned. The deduced OspC amino acid sequences were aligned with 12 published OspC sequences and revealed the presence of 48 conserved amino acids. On the basis of the alignment, OspC could be divided into an amino-terminal relatively conserved region and a relatively variable region in the central portion. The distance tree obtained divided the ospC sequences into three groups. The first group contained ospC alleles from all (n = 13) sensu stricto strains, the second group contained ospC alleles from seven Borrelia afzelii strains, and the third group contained ospC alleles from five B. afzelii and all (n = 9) Borrelia garinii strains. The ratio of the mean number of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions per site calculated for B. burgdorferi sensu stricto, B. garinii, and B. afzelii ospC alleles suggested that the polymorphism of OspC is due to positive selection favoring diversity at the amino acid level in the relatively variable region. On the basis of the comparison of 16S rRNA gene sequences, Borrelia hermsii is more closely related to B. afzelii than to B. burgdorferi sensu stricto and B. garinii. In contrast, the phylogenetic tree obtained for the B. hermsii variable major protein, Vmp33, and 18 OspC amino acid sequences suggested that Vmp33 and OspC from B. burgdorferi sensu stricto strains share a common evolutionary origin.