Selection pressure on Haemagglutinin genes of H9N2 influenza viruses from different hosts

Positive selection and differential selective pressure analyses were carried out to study Haemagglutinin (HA) genes of H9N2 influenza viruses from different hosts in this paper. Results showed that, although most positions in HAs were under neutral or purifying evolution, a few positions located in the antigenic regions and receptor binding sites were subject to positive selection and some of them were even positively selected at the population level. In addition, there were always some positions differentially selected for viruses from different hosts. Both selection pressure working on HA codons and positions differentially selected might account for the extension of the host range and adaptations to different hosts of H9N2 influenza viruses.     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsynonymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in viruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected positions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

Antiviral activity of pyridyl imidazolidinones against enterovirus 71 variants

Pyridyl imidazolidinone is a novel class of capsid binder which can inhibit enterovirus 71 (EV71). In this study, we tested the susceptibility of six recombinant viruses with different single-site mutations in VP1. Eleven modified pyridyl imidazolidinones were synthesized and used to probe the interaction between these compounds and the EV71 VP1 protein. We found that the D31N or E98K mutant viruses were susceptible to bulkier compounds, which suggested that mutations at these two sites in VP1 may widen the hydrophobic pocket of VP1, allowing bulkier compounds to enter and interfere VP1-receptor binding. Additionally, the Y116H mutant was more resistant to pyridyl imidazolidinone compounds containing a methyl group in the central position of the hydrophobic linker. When a trifluoromethyl group was substituted for the methyl group in the middle of the linker chain, the inhibitory effect was totally abolished in the Y116H mutant, suggesting that the interaction between Tyr (Y) 116 of VP1 and the central position of the linker chain of pyridyl imidazolodinone is very important for drug efficacy. A V192M mutant was resistant to most of the derivatives, indicating that residue 192 is a key mutation for resistance to pyridyl imidazolidinone.     

Recombinogenic activity of chimeric recA genes (Pseudomonas aeruginosa/Escherichia coli): a search for RecA protein regions responsible for this activity

The nature of selection on capsid genes of foot-and-mouth disease virus (FMDV) was characterized by examining the ratio of nonsynonymous to synonymous substitutions in 11 data sets of sequences obtained from six different serotypes of FMDV. Using a method of analysis that assigns each codon position to one of a number of estimated values of nonsynonymous to synonymous ratio, significant evidence of positive selection was identified in 5 data sets, operating at 1-7% of codon positions. Evidence of positive selection was identified in complete capsid sequences of serotypes A and C and in VP1 sequences of serotypes SAT 1 and 2. Sequences of serotype SAT-2 recovered from a persistently infected African buffalo also revealed evidence for positive selection. Locations of codons under positive selection coincide closely with those of antigenic sites previously identified with the use of monoclonal antibody escape mutants. The vast majority of codons are under mild to strong purifying selection. However, these results suggest that arising antigenic variants benefit from a selective advantage in their interaction with the immune system, either during the course of an infection or in transmission to individuals with previous exposure to antigen. Analysis of amino acid usage at sites under positive selection indicates that this selective advantage can be conferred by amino acid substitutions that share physicochemically similar properties.     

Molecular analysis of virulent determinants of enterovirus 71

Enterovirus 71 (EV71) is the most important causative agent of hand, foot and mouth disease (HFMD) in children. In most cases, it is a self-limiting illness. However some EV71 infectious cases can develop severe clinical outcomes, such as encephalitis, meningitis, poliomyelitis like paralysis, and even death. To identify the determinants of virulence, the deduced amino acid sequence of polyprotein and nucleotide sequence of 5'-NTR and 3'-NTR in 25 SC-EV71 strains (strains from severe cases) and 31 MC-EV71 strains (strains from mild cases) were analyzed. Results showed four amino acids on two positions (Gly(P710)/Gln(P710)/Arg(P710) and Glu(P729)) on the DE and EF loop of VP1, one (Lys(P930)) on the surface of protease 2A and four nucleotides on three positions (G(P272), U(P488) and A(P700)/U(P700)) in the 5'-NTR region are associated with EV71 virulent phenotype. Predicted secondary structure of RNA using the consensus sequence of 5'-NTR by RNAStructure showed the mutation of nucleotide at position 488 in strain BJ08-Z004-3 (position 491 in prototype strain BrCr) can result in the discrepancy of an additional pair of nucleotides and thus change the stability of the second structure of IRES. Fragment base content analysis showed that in the region 696 to 714 bp at the 5'-NTR, where the A(P700)/U(P700) was located, the nucleotide constitution ratios differed significantly between SC-EV71 and MC-EV71 strains. In conclusion, comparative genomic analysis showed that virulence of EV71 strains are mainly determined by the amino acids on two positions of VP1, one position of protease 2A and the nucleotides on three positions in 5'-NTR.     

Identification of cross-reactive and serotype 2-specific neutralization epitopes on VP3 of human rotavirus.

The group A rotaviruses are composed of at least seven serotypes. Serotype specificity is defined mainly by an outer capsid protein, VP7. In contrast, the other surface protein, VP3 (775 amino acids), appears to be associated with both serotype-specific and heterotypic immunity. To identify the cross-reactive and serotype-specific neutralization epitopes on VP3 of human rotavirus, we sequenced the VP3 gene of antigenic mutants resistant to each of seven anti-VP3 neutralizing monoclonal antibodies (N-MAbs) which exhibited heterotypic or serotype 2-specific reactivity, and we defined three distinct neutralization epitopes on VP3. The mutants sustained single amino acid substitutions at position 305, 392, 433, or 439. Amino acid position 305 was critical to epitope I, whereas amino acid position 433 was critical to epitope III. In contrast, epitope II appeared to be more dependent upon conformation and protein folding because both amino acid positions 392 and 439 appeared to be critical. These four positions clustered in a relatively limited area of VP5, the larger of the two cleavage products of VP3. At the positions where amino acid substitutions occurred, there was a correlation between amino acid sequence homology among different serotypes and the reactivity patterns of various viruses with the N-MAbs used for selection of mutants. A synthetic peptide (amino acids 296 to 313) which included the sequence of epitope I reacted with its corresponding N-MAb, suggesting that the region contains a sequential antigenic determinant. These data may prove useful in current efforts to develop vaccines against human rotavirus infection.     

Two proline residues are essential in the calcium-binding activity of rotavirus VP7 outer capsid protein.

Rotavirus maturation and stability of the outer capsid are calcium-dependent processes. It has been shown previously that the concentration of Ca2+-solubilizing outer capsid proteins from rotavirus particles is dependent on the virus strain. This property of viral particles has been associated with the gene coding for VP7 (gene 9). In this study the correlation between VP7 and resistance to low [Ca2+] was confirmed by analyzing the origin of gene 9 from reassortant viruses prepared under the selective pressure of low [Ca2+]. After chemical mutagenesis, we selected mutant viruses of the bovine strain RF that are more resistant to low [Ca2+]. The genes coding for the VP7 proteins of these independent mutants have been sequenced. Sequence analysis confirmed that these mutants are independent and revealed that all mutant VP7 proteins have proline 75 changed to leucine and have an outer capsid that solubilized at low [Ca2+]. The mutation of proline 279 to serine is found in all but two mutants. The phenotype of mutants having a single proline change can be distinguished from the phenotype of mutants having two proline changes. Sequence analysis showed that position 75 is in a region (amino acids 65 to 78) of great variability and that proline 75 is present in most of the bovine strains. In contrast, proline 279 is in a conserved region and is conserved in all the VP7 sequences in data banks. This region is rich in oxygenated residues that are correctly allocated in the metal-coordinating positions of the Ca2+-binding EF-hand structure pattern, suggesting that this region is important in the Ca2+ binding of VP7.     

The Role of VP1 Amino Acid Residue 145 of Enterovirus 71 in Viral Fitness and Pathogenesis in a Cynomolgus Monkey Model

Enterovirus 71 (EV71), a major causative agent of hand, foot, and mouth disease, occasionally causes severe neurological symptoms. We identified P-selectin glycoprotein ligand-1 (PSGL-1) as an EV71 receptor and found that an amino acid residue 145 in the capsid protein VP1 (VP1-145) defined PSGL-1-binding (PB) and PSGL-1-nonbinding (non-PB) phenotypes of EV71. However, the role of PSGL-1-dependent EV71 replication in neuropathogenesis remains poorly understood. In this study, we investigated viral replication, genetic stability, and the pathogenicity of PB and non-PB strains of EV71 in a cynomolgus monkey model. Monkeys were intravenously inoculated with cDNA-derived PB and non-PB strains of EV71, EV71-02363-EG and EV71-02363-KE strains, respectively, with two amino acid differences at VP1-98 and VP1-145. Mild neurological symptoms, transient lymphocytopenia, and inflammatory cytokine responses, were found predominantly in the 02363-KE-inoculated monkeys. During the early stage of infection, viruses were frequently detected in clinical samples from 02363-KE-inoculated monkeys but rarely in samples from 02363-EG-inoculated monkeys. Histopathological analysis of central nervous system (CNS) tissues at 10 days postinfection revealed that 02363-KE induced neuropathogenesis more efficiently than that induced by 02363-EG. After inoculation with 02363-EG, almost all EV71 variants detected in clinical samples, CNS, and non-CNS tissues, possessed a G to E amino acid substitution at VP1-145, suggesting a strong in vivo selection of VP1-145E variants and CNS spread presumably in a PSGL-1-independent manner. EV71 variants with VP1-145G were identified only in peripheral blood mononuclear cells in two out of four 02363-EG-inoculated monkeys. Thus, VP1-145E variants are mainly responsible for the development of viremia and neuropathogenesis in a non-human primate model, further suggesting the in vivo involvement of amino acid polymorphism at VP1-145 in cell-specific viral replication, in vivo fitness, and pathogenesis in EV71-infected individuals.     

Enterovirus 71 Binding to PSGL-1 on Leukocytes: VP1-145 Acts as a Molecular Switch to Control Receptor Interaction

Some strains of enterovirus 71 (EV71), but not others, infect leukocytes by binding to a specific receptor molecule: the P-selectin glycoprotein ligand-1 (PSGL-1). We find that a single amino acid residue within the capsid protein VP1 determines whether EV71 binds to PSGL-1. Examination of capsid sequences of representative EV71 strains revealed that the PSGL-1-binding viruses had either a G or a Q at residue 145 within the capsid protein VP1 (VP1-145G or Q), whereas PSGL-1-nonbinding viruses had VP1-145E. Using site-directed mutagenesis we found that PSGL-1-binding strains lost their capacity to bind when VP1-145G/Q was replaced by E; conversely, nonbinding strains gained the capacity to bind PSGL-1 when VP1-145E was replaced with either G or Q. Viruses with G/Q at VP1-145 productively infected a leukocyte cell line, Jurkat T-cells, whereas viruses with E at this position did not. We previously reported that EV71 binds to the N-terminal region of PSGL-1, and that binding depends on sulfated tyrosine residues within this region. We speculated that binding depends on interaction between negatively charged sulfate groups and positively charged basic residues in the virus capsid. VP1-145 on the virus surface is in close proximity to conserved lysine residues at VP1-242 and VP1-244. Comparison of recently published crystal structures of EV71 isolates with either Q or E at VP1-145 revealed that VP1-145 controls the orientation of the lysine side-chain of VP1-244: with VP1-145Q the lysine side chain faces outward, but with VP1-145E, the lysine side chain is turned toward the virus surface. Mutation of VP1-244 abolished virus binding to PSGL-1, and mutation of VP1-242 greatly reduced binding. We propose that conserved lysine residues on the virus surface are responsible for interaction with sulfated tyrosine residues at the PSGL-1 N-terminus, and that VP1-145 acts as a switch, controlling PSGL-1 binding by modulating the exposure of VP1-244K.     

A Novel Universal Neutralizing Monoclonal Antibody against Enterovirus 71 That Targets the Highly Conserved “Knob” Region of VP3 Protein

Hand, foot and mouth disease caused by enterovirus 71(EV71) leads to the majority of neurological complications and death in young children. While putative inactivated vaccines are only now undergoing clinical trials, no specific treatment options exist yet. Ideally, EV71 specific intravenous immunoglobulins could be developed for targeted treatment of severe cases. To date, only a single universally neutralizing monoclonal antibody against a conserved linear epitope of VP1 has been identified. Other enteroviruses have been shown to possess major conformational neutralizing epitopes on both the VP2 and VP3 capsid proteins. Hence, we attempted to isolate such neutralizing antibodies against conformational epitopes for their potential in the treatment of infection as well as differential diagnosis and vaccine optimization. Here we describe a universal neutralizing monoclonal antibody that recognizes a conserved conformational epitope of EV71 which was mapped using escape mutants. Eight escape mutants from different subgenogroups (A, B2, B4, C2, C4) were rescued; they harbored three essential mutations either at amino acid positions 59, 62 or 67 of the VP3 protein which are all situated in the “knob” region. The escape mutant phenotype could be mimicked by incorporating these mutations into reverse genetically engineered viruses showing that P59L, A62D, A62P and E67D abolish both monoclonal antibody binding and neutralization activity. This is the first conformational neutralization epitope mapped on VP3 for EV71.     

Recombinant VP1 protein expressed in Pichia pastoris induces protective immune responses against EV71 in mice

Human enterovirus 71 (EV71) is one of the major causative agents of hand, foot and mouth disease and is also associated with serious neurological diseases in children. Currently, there are no effective antiviral drugs or vaccines against EV71 infection. VP1, one of the major immunogenic capsid proteins of EV71, is widely considered to be the candidate antigen for an EV71 vaccine. In this study, VP1 of EV71 was expressed as a secretory protein with an N-terminal histidine tag in the methylotrophic yeast Pichia pastoris, and purified by Ni–NTA affinity chromatography. Immunogenicity and vaccine efficacy of the recombinant VP1 were assessed in mouse models. The results showed that the recombinant VP1 could efficiently induce anti-VP1 antibodies in BALB/c mice, which were able to neutralize EV71 viruses in an in vitro neutralization assay. Passive protection of neonatal mice further confirmed the prophylactic efficacy of the antisera from VP1 vaccinated mice. Furthermore, VP1 vaccination induced strong lymphoproliferative and Th1 cytokine responses. Taken together, our study demonstrated that the yeast-expressed VP1 protein retained good immunogenicity and was a potent EV71 vaccine candidate.     

Structural features of the scaffold interaction domain at the N terminus of the major capsid protein (VP5) of herpes simplex virus type 1

Protein-protein interactions drive the assembly of the herpes simplex virus type 1 capsid. A key interaction occurs between the C terminus of the scaffold protein and the N terminus of the major capsid protein (VP5). Results from alanine-scanning mutagenesis of hydrophobic residues in the N terminus of VP5 revealed seven residues (I27, L35, F39, L58, L65, L67, and L71) that reside in two predicted alpha helices (helix 1(22-42) and helix 2(58-72)) that are important for this bimolecular interaction. The goal of the present study was to further characterize the VP5 scaffold interaction domain (SID). Amino acids at the seven positions were replaced with L, M, V or P (I27); I, M, V, or P (L35, L58, L65, L67, and L71); and H, W, Y, or L (F39). Replacement with a hydrophobic side chain did not affect the interaction with scaffold protein in yeast cells or the ability of a virus specifying the mutation from replicating in cells. The mutation to the proline side chain abolished the interaction in all cases and was lethal for virus replication. Mutant viruses with proline substitutions in helix 1(22-42) at positions 27 and 35 assembled large open capsid shells that did not attain closure. Proline substitutions in helix 2(58-72) at either position 59, 65, or 67 abolished the accumulation of VP5 protein, and, at 58 and 71, although VP5 did accumulate, capsid shells were not assembled. Thus, the second SID, SID2, is highly structured, and this alpha helix (helix 2(58-72)) is likely involved in capsomere-capsomere interactions during shell accretion. Conserved glycine G59 in helix 2(58-72) was also mutated. G59 may act as a flexible "hinge" in helix 2(58-72) because decreasing the movement of this side chain by replacement with valine impaired capsid assembly. Thus, the N terminus of VP5 and the alpha helices embedded in this domain, as in the capsid shell proteins of some double-stranded DNA phages, are a key regulator of shell accretion and stabilization.     

BoLA class I allele diversity and polymorphism in a herd of cattle

Major histocompatibility complex class I genes are among the most polymorphic genes characterized. The high level of polymorphism is essential for generating host immune responses. In humans, three distinct genomic loci encode human leukocyte antigen (HLA) class I genes, allowing individuals to express up to six different HLA class I molecules. In cattle, the number of distinct genomic loci are currently at least six, and the number of different bovine leukocyte antigens (BoLA) class I molecules that are expressed in individual animals are variable. The extent of allele variation within the cattle population is unknown. In this study, the number and variety of BoLA class I sequences expressed by 36 individuals were determined from full-length BoLA class I cDNA clones. Twenty distinct BoLA class I alleles were identified, with only four being previously reported. The number of expressed BoLA class I alleles in individual animals ranged between one and four, with none of the animals having an identical complement of BoLA class I molecules. Variation existed in the number of BoLA class I alleles expressed as well as the composition of expressed alleles, however, several BoLA class I alleles were found in multiple individual animals. Polymorphic amino acid sites were analyzed for positive and negative selection using the ADAPTSITE program. In the antigen recognition sites (ARS), there were eight positions that were predicted to be under positive selection and three positions that were predicted to be under negative selection from 62 positions. In contrast, for non-antigen recognition sites (non-ARS), there were three positions that were predicted to be under positive selection and 20 that were predicted to be under negative selection from 278, indicating that positive selection of amino acids occurs at a greater frequency within the antigen recognition sites.     

Display of VP1 on the surface of baculovirus and its immunogenicity against heterologous human enterovirus 71 strains in mice

Background

Human Enterovirus 71 (EV71) is a common cause of hand, foot and mouth disease (HFMD) in young children. It is often associated with severe neurological diseases and has caused high mortalities in recent outbreaks across the Asia Pacific region. Currently, there is no effective vaccine and antiviral agents available against EV71 infections. VP1 is one of the major immunogenic capsid protein of EV71 and plays a crucial role in viral infection. Antibodies against VP1 are important for virus neutralization.

Methodology/Principal Finding

In the present study, infectious EV71 viruses were generated from their synthetic complementary DNA using the human RNA polymerase I reverse genetics system. Secondly, the major immunogenic capsid protein (VP1) of EV71-Fuyang (subgenogroup C4) was displayed on the surface of recombinant baculovirus Bac-Pie1-gp64-VP1 as gp64 fusion protein under a novel White Spot Syndrome Virus (WSSV) immediate early ie1 promoter. Baculovirus expressed VP1 was able to maintain its structural and antigenic conformity as indicated by immunofluorescence assay and western blot analysis. Interestingly, our results with confocal microscopy revealed that VP1 was able to localize on the plasma membrane of insect cells infected with recombinant baculovirus. In addition, we demonstrated with transmission electron microscopy that baculovirus successfully acquired VP1 from the insect cell membrane via the budding process. After two immunizations in mice, Bac-Pie1-gp64-VP1 elicited neutralization antibody titer of 1∶64 against EV71 (subgenogroup C4) in an in vitro neutralization assay. Furthermore, the antisera showed high cross-neutralization activities against all 11 subgenogroup EV71 strains.

Conclusion

Our results illustrated that Bac-Pie1-gp64-VP1 retained native epitopes of VP1 and acted as an effective EV71 vaccine candidate which would enable rapid production without any biosafety concerns.     

GeXP多重基因表达遗传分析系统在手足口病病原分型检测中的应用

利用GeXP多重基因表达遗传分析系统,建立一种多重逆转录-聚合酶链反应(RT-PCR)方法,同时检测引起手足口病的9种常见的人肠道病毒—人肠道病毒71型(HEV71)、柯萨奇病毒A组(CVA)16、4、5、9、10型和柯萨奇病毒B组(CVB)1、3、5型。优化多重反应体系中针对5’UTR区的肠道病毒通用引物和11对针对9种血清型人肠道病毒VP1区的特异性引物的浓度比例,分别以病毒细胞培养物和阳性粪便标本来验证多重反应体系的特异性,以TCID50定量的细胞培养物和克隆质粒体外转录的RNA梯度稀释液来检测多重检测体系的灵敏度。结果表明,优化后的多重检测体系,可扩增出人肠道病毒共有的保守片段的和型特异性片段,HEV71和CVA16细胞培养物的检测下限为100.5TCID50/μL,并可在103copies/μL水平同时、特异地检测出9种病毒RNA。该方法灵敏度高、特异性强,可快速对大量临床样本进行高通量检测,用于手足口病的分子流行病学调查。     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsyn-onymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in vi-ruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected posi-tions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

Translation of the minor capsid protein of a calicivirus is initiated by a novel termination-dependent reinitiation mechanism

Caliciviruses represent a family of positive strand RNA viruses responsible for a variety of syndromes in man and animals. VP10, a minor structural protein of the calicivirus rabbit hemorrhagic disease virus, is encoded in the small 3'-terminal open reading frame (ORF) 2 and is translated with an efficiency of approximately 20% of the preceding ORF1. The presence of the ORF1 termination codon is crucial for VP10 expression. Translation of VP10 starts at an AUG codon located at positions -5 to -3 of the ORF1 termination codon. However, VP10 was also expressed in the absence of an AUG initiation codon. The majority of ORF1 could be deleted or replaced by different sequences without significant influence on VP10 expression as long as translation terminated at the given position. The RNA sequence of the 3'-terminal 84 nucleotides of ORF1 but not the encoded peptide was found to be crucial for VP10 expression. In contrast, nearly the entire ORF2 could be replaced by a foreign sequence without abrogation of its translation. Accordingly, VP10 is expressed in a translation termination/reinitiation process that is particular because it is independent of an AUG translational start codon and requires the presence of a sequence element upstream of the initiation site.     

miR-142-3p is essential for hematopoiesis and affects cardiac cell fate in zebrafish

Hand, foot, and mouth disease (HFMD) is a common viral illness in young children. HFMD is caused by viruses belonging to the enterovirus genus of the picornavirus family. Recently, enterovirus 71 (EV71) has emerged as a virulent agent for HFMD with severe clinical outcomes. In the current report, we conducted a pilot antigen engineering study to optimize the expression and immunogenicity of subunit VP1 antigen for the design of EV71 vaccines. DNA immunization was adopted as a simple technical approach to test different designs of VP1 antigens without the need to express VP1 protein in vitro first. Our studies indicated that the expression and immunogenicity of VP1 protein can be improved with alternated VP1 antigen designs. Data presented in the current report revealed novel pathways to optimize the design of VP1 antigen-based EV71 vaccines.     

Nucleotide sequence variation of the envelope protein gene identifies two distinct genotypes of yellow fever virus.

The evolution of yellow fever virus over 67 years was investigated by comparing the nucleotide sequences of the envelope (E) protein genes of 20 viruses isolated in Africa, the Caribbean, and South America. Uniformly weighted parsimony algorithm analysis defined two major evolutionary yellow fever virus lineages designated E genotypes I and II. E genotype I contained viruses isolated from East and Central Africa. E genotype II viruses were divided into two sublineages: IIA viruses from West Africa and IIB viruses from America, except for a 1979 virus isolated from Trinidad (TRINID79A). Unique signature patterns were identified at 111 nucleotide and 12 amino acid positions within the yellow fever virus E gene by signature pattern analysis. Yellow fever viruses from East and Central Africa contained unique signatures at 60 nucleotide and five amino acid positions, those from West Africa contained unique signatures at 25 nucleotide and two amino acid positions, and viruses from America contained such signatures at 30 nucleotide and five amino acid positions in the E gene. The dissemination of yellow fever viruses from Africa to the Americas is supported by the close genetic relatedness of genotype IIA and IIB viruses and genetic evidence of a possible second introduction of yellow fever virus from West Africa, as illustrated by the TRINID79A virus isolate. The E protein genes of American IIB yellow fever viruses had higher frequencies of amino acid substitutions than did genes of yellow fever viruses of genotypes I and IIA on the basis of comparisons with a consensus amino acid sequence for the yellow fever E gene. The great variation in the E proteins of American yellow fever virus probably results from positive selection imposed by virus interaction with different species of mosquitoes or nonhuman primates in the Americas.     

Molecular evolution of the transmembrane domains of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are a superfamily of integral membrane proteins vital for signaling and are important targets for pharmaceutical intervention in humans. Previously, we identified a group of ten amino acid positions (called key positions), within the seven transmembrane domain (7TM) interhelical region, which had high mutual information with each other and many other positions in the 7TM. Here, we estimated the evolutionary selection pressure at those key positions. We found that the key positions of receptors for small molecule natural ligands were under strong negative selection. Receptors naturally activated by lipids had weaker negative selection in general when compared to small molecule-activated receptors. Selection pressure varied widely in peptide-activated receptors. We used this observation to predict that a subgroup of orphan GPCRs not under strong selection may not possess a natural small-molecule ligand. In the subgroup of MRGX1-type GPCRs, we identified a key position, along with two non-key positions, under statistically significant positive selection.