Sequence analysis of NS3 protease gene in clinical strains of hepatitis C virus

The amino terminal region of the non structural gene 3 (NS3) of hepatitis C virus (HCV) is a chymotripsinlike serine-protease responsible for cleavage of the non structural proteins of Hepatitis C virus (HCV). In order to investigate the genetic variation of this region, we developed a nested PCR to obtain NS3 protease sequences from 54 patients chronically infected with HCV genotypes 1a, 1b and 3, respectively. Comparison of nucleotide and amino acids sequences of NS3 protease domain with consensus sequence obtained within the same genotype, showed 3.73% nucleotide divergence and 1.64% amino acid divergence in isolates of genotype 3a, whereas isolates 1a exhibited 4.45% nucleotide and 4% amino acid change, respectively. Finally, NS3 sequence from 1b isolates revealed 6.47% nucleotide and 3.5 % aa changes. Comparison of consensus amino acid sequences derived from isolates 1a, 1b and 3, with the HCV prototypes showed a low amino acid sequence diversity. However, the consensus sequence of HCV genotype 3 isolates showed an amino acid changed from the prototype, that was located within a region important for enzyme structure and activity. These results indicated that the NS3 protease gene is highly conserved within the same HCV genotype. The domains involved in enzyme function were highly conserved in 1a and 1b strains, whereas consensus sequence of isolates 3a showed that the majority of these strains were not perfectly conserved in one of such regions. These findings altogether suggested that the NS3 protease enzyme of HCV may constitute an important target for antiviral therapy, but the NS3 protease variability of isolates 3 within a region that is a potential target for antiviral therapy could pose a problem for structure based drug development.     

Unusual distribution of mutations associated with serial bottleneck passages of human immunodeficiency virus type 1

Repeated bottleneck passages result in fitness losses of RNA viruses. In the case of human immunodeficiency virus type 1 (HIV-1), decreases in fitness after a limited number of plaque-to-plaque transfers in MT-4 cells were very drastic. Here we report an analysis of entire genomic nucleotide sequences of four HIV-1 clones derived from the same HIV-1 isolate and their low-fitness progeny following 7 to 15 plaque-to-plaque passages. Clones accumulated 4 to 28 mutations per genome, with dominance of A --> G and G --> A transitions (57% of all mutations) and 49% nonsynonymous replacements. One clone-but not three sibling clones-showed an overabundance of G --> A transitions, evidencing the highly stochastic nature of some types of mutational bias. The distribution of mutations along the genome was very unusual in that mutation frequencies in gag were threefold higher than in env. Particularly striking was the complete absence of replacements in the V3 loop of gp120, confirmed with partial nucleotide sequences of additional HIV-1 clones subjected to repeated bottleneck passages. The analyses revealed several amino acid replacements that have not been previously recorded among natural HIV-1 isolates and illustrate how evolution of an RNA virus genome, with regard to constant and variable regions, can be profoundly modified by alterations in population dynamics.     

Phylogenetic analysis of viral haemorrhagic septicaemia virus (VHSV) isolates from France (1971-1999)

The nucleotide sequences of a specific region of the glycoprotein gene were compared among 63 strains of viral haemorrhagic septicaemia virus (VHSV) isolated from fish in France between 1971 and 1999. The analysis was performed on a region corresponding to amino acids 238 to 331 of the glycoprotein gene, also designated the V2 region and previously shown to accumulate most of the mutations. The sequences of many VHSV isolates were found to be identical or very conserved. An isolate, designated L59X, obtained from elver in the Loire estuary, depicted a higher degree of divergence compared to the other French isolates. The deduced amino-acid sequences were analysed together with the results of neutralisation tests performed using monoclonal antibody 168m4 specific to serotype 1. Non-neutralised VHSV strains had mutations in the region corresponding to the previously described 168m4 epitope. Phylogenetic analysis showed that all the VHSV isolates studied, except L59X, belong to genotype I, previously described as containing VHSV strains isolated from continental Europe. Most of the VHSV isolates studied were found to be genetically related to one of the previously described VHSV strains representative of the major serotypes. Isolate L59X, which was the only French marine strain studied, was found to belong to genotype II, previously shown to encompass the VHSV strains isolated from the British Isles coastal waters. Overall there was a good correlation between the geographical origin of the studied isolates and their genetic characteristics.     

Rates and fitness consequences of new mutations in humans

The human mutation rate per nucleotide site per generation (μ) can be estimated from data on mutation rates at loci causing Mendelian genetic disease, by comparing putatively neutrally evolving nucleotide sequences between humans and chimpanzees and by comparing the genome sequences of relatives. Direct estimates from genome sequencing of relatives suggest that μ is about 1.1 × 10(-8), which is about twofold lower than estimates based on the human-chimp divergence. This implies that an average of ~70 new mutations arise in the human diploid genome per generation. Most of these mutations are paternal in origin, but the male:female mutation rate ratio is currently uncertain and might vary even among individuals within a population. On the basis of a method proposed by Kondrashov and Crow, the genome-wide deleterious mutation rate (U) can be estimated from the product of the number of nucleotide sites in the genome, μ, and the mean selective constraint per site. Although the presence of many weakly selected mutations in human noncoding DNA makes this approach somewhat problematic, estimates are U ≈ 2.2 for the whole diploid genome per generation and 0.35 for mutations that change an amino acid of a protein-coding gene. A genome-wide deleterious mutation rate of 2.2 seems higher than humans could tolerate if natural selection is "hard," but could be tolerated if selection acts on relative fitness differences between individuals or if there is synergistic epistasis. I argue that in the foreseeable future, an accumulation of new deleterious mutations is unlikely to lead to a detectable decline in fitness of human populations.     

Phylogenetic relationship of street rabies virus strains and their antigenic reactivity with antibodies induced by vaccine strains. I. Analysis of phylogenetic relationship of street rabies virus strains isolated in Poland

The aims of these studies were: genetic characteristic of street rabies virus strains isolated from different animal species in Poland and determination of phylogenetic relationships to reference laboratory strains of the street rabies viruses belonging to genotype 1 and 5. The variability of rabies isolates and their phylogenetic relationship were studied by comparing the nucleotide sequence of the virus genome fragment. The Polish strains of genotype 1 belong to four phylogenetic groups (NE, CE, NEE, EE) corresponding to four variants: fox-racoon dog (F-RD); European fox 1 (F1); European fox 2 (F2) and European fox 3 (F3). On the Polish territories there are no rabies strains representing the variant dog-wolf and typical for arctic fox variant. The similarity of nucleotide and amino acid sequences of street rabies strains belonging to genotype 1 and laboratory strain CVS is very high. It is about 91% similarity at nucleotide level and 95% at amino acid level. Rabies strain CVS is similar to genotype 5 bat strains (EBL 1) only in about 69% and 74% at nucleotide and amino acid level, respectively. The genetic divergence of rabies strains circulating in Poland raised the need of permanent epidemiological and virological surveillance. The genotype and variant of isolated strains should be determined (using PCR and RLFP methods).     

Interferon Resistance of Hepatitis C Virus Genotype 1b: Relationship to Nonstructural 5A Gene Quasispecies Mutations

A 40-amino-acid sequence located in the nonstructural 5A (NS5A) protein of hepatitis C virus genotype 1b (HCV-1b) was recently suggested to be the interferon sensitivity-determining region (ISDR), because HCV-1b strains with an ISDR amino acid sequence identical to that of the prototype strain HCV-J were found to be resistant to alpha interferon (IFN-α) whereas strains with amino acid substitutions were found to be sensitive (N. Enomoto, I. Sakuma, Y. Asahina, M. Kurosaki, T. Murakami, C. Yamamoto, N. Izumi, F. Marumo, and C. Sato, J. Clin. Invest. 96:224–230, 1995; N. Enomoto, I. Sakuma, Y. Asahina, M. Kurosaki, T. Murakami, C. Yamamoto, Y. Ogura, N. Izumi, F. Marumo, and C. Sato, N. Engl. J. Med. 334:77–81, 1996). We used single-strand conformation polymorphism (SSCP) analysis, combined with cloning and sequencing strategies, to characterize NS5A quasispecies in HCV-1b-infected patients and determine the relationships between pre- and posttreatment NS5A quasispecies mutations and the IFN-α sensitivity of HCV-1b. The serine residues involved in phosphorylation of NS5A protein were highly conserved both in the various patients and in quasispecies in a given patient, suggesting that phosphorylation is important in NS5A protein function. A hot spot for amino acid substitutions was found at positions 2217 to 2218; it could be the result of either strong selection pressure or tolerance to these amino acid replacements. The proportion of synonymous mutations was significantly higher than the proportion of nonsynonymous mutations, suggesting that genetic variability in the region studied was the result of high mutation rates and viral replication kinetics rather than of positive selection. Sustained HCV RNA clearance was associated with low viral load and low nucleotide sequence entropy, suggesting (i) that the replication kinetics when treatment is started plays a critical role in HCV-1b sensitivity to IFN-α and (ii) that HCV-1b resistance to IFN-α could be conferred by numerous and/or related mutations that could be patient specific and located at different positions throughout the viral genome and could allow escape variants to be selected by IFN-α-stimulated immune responses. No NS5A sequence appeared to be intrinsically resistant or sensitive to IFN-α, but the HCV-J sequence was significantly more frequent in nonresponder quasispecies than in sustained virological responder quasispecies, suggesting that the balance between NS5A quasispecies sequences in infected patients could have a subtle regulatory influence on HCV replication.     

In-silico genomic landscape characterization and evolution of SARS-CoV-2 variants isolated in India shows significant drift with high frequency of mutations

In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.     

Full-length sequence analysis of subacute sclerosing panencephalitis (SSPE) virus, a mutant of measles virus, isolated from brain tissues of a patient shortly after onset of SSPE

Subacute sclerosing panencephalitis (SSPE) virus, a measles virus (MeV) mutant, was isolated from brain tissues of a patient shortly after the clinical onset, and the entire viral genome was sequenced. The virus, named SSPE-Kobe-1, formed syncytia on B95a and Vero/SLAM cells without producing cell-free infectious virus particles, which is characteristic of SSPE virus. Phylogenetic analysis classified SSPE-Kobe-1 into genotype D3. When compared with an MeV field isolate of the same genotype (Ich-B strain), SSPE-Kobe-1 exhibited mutation rates of 0.8-1.6% at the nucleotide level in each of the proteincoding regions of the viral genome. It is noteworthy that the mutation rate of the M gene (1.2%) of SSPE-Kobe-1 was considerably lower than for other SSPE virus strains reported so far, but that the majority of the mutations (75%) were the uridine-to-cytidine biased hypermutation characteristic of the SSPE virus M gene. At the amino acid level, the viral proteins, such as N, P, C, V, M, F, H and L proteins, had point-mutations on 3, 7, 1, 4, 3, 9, 8 and 14 residues, respectively, compared with the Ich-B strain. In addition, the F and H proteins had mutated C-termini due to single-point mutations near or at the stop codons. Two of the three mutations in the M protein were Leu-to-Pro mutations, which are likely to affect the conformation and, therefore, the function of the protein. Because of the relatively small number of mutations, SSPE-Kobe-1 would be a useful tool to study genetic evolution of SSPE virus.     

Nucleotide substitutions in the acetolactate synthase genes of sulfonylurea-resistant biotypes of Monochoria vaginalis (Pontederiaceae)

Some point mutations in acetolactate synthase (ALS) confer resistance to ALS-inhibiting herbicides in weeds. To clarify the evolution of the herbicide resistance of Monochoria vaginalis, a weed in rice fields in Japan, the nucleotide sequences of four genes encoding ALS were surveyed in five sulfonylurea-resistant (SU-R) and five sulfonylurea-susceptible (SU-S) biotypes. In the ALS1 gene, two SU-R biotypes showed nucleotide substitutions changing Pro197 to Ser and Leu, respectively. In a different gene, ALS3, three other SU-R biotypes showed either of the two nonsynonymous nucleotide substitutions seen in ALS1. Only two biotypes geographically located distantly from each other shared the same mutation conferring SU resistance in the same gene. These patterns of nucleotide substitutions indicate that the SU-R phenotype was acquired independently by different biotypes. Nucleotide diversity values of the genes showing SU-R mutations were higher than those of ALS2 lacking any SU-R mutation and of a putative pseudogene, ALS4. This result suggests that the maintenance of nucleotide variability within target genes provides an opportunity for the evolution of SU-R phenotypes by herbicide-driven selection for mutations conferring resistance.     

Characterization of hypervariable regions in the putative envelope protein of hepatitis C virus.

We previously identified two hypervariable regions [HVR1 (27 amino acids) and HVR2 (7 amino acids)] in the putative envelope glycoprotein (gp70) by comparison of the amino acid sequences of many isolates of the HCV-II genotype. To understand the functional features of these HVRs, using the polymerase chain reaction we analyzed the rate of actual sequence variability in the region including HVR1 and HVR2 of HCV isolated successively at intervals of several months from two patients with chronic C-type hepatitis. In both patients, the amino acid sequence of HVR1, but not HVR2, was found to change dramatically during the observation period (about one amino acid per month). However, no alteration of the amino acid sequence of HVR1 of HCV was observed in a patient in the acute phase of chronic hepatitis. Restriction digestion analysis of sequence diversity showed that a HCV genome with a newly introduced mutation in HVR1 often became the predominant population at the next time of examination. Alterations of amino acids in HVR1 occurred sequentially in the two patients in the chronic phase. These findings suggest that mutations in HVR1 are involved in the mechanism of persistent chronic HCV infection.     

The Complete Genome Phylogeny of Geographically Distinct Dengue Virus Serotype 2 Isolates (1944-2013) Supports Further Groupings within the Cosmopolitan Genotype

Dengue virus serotype 2 (DENV-2) isolates have been implicated in deadly outbreaks of dengue fever (DF) and dengue hemorrhagic fever (DHF) in several regions of the world. Phylogenetic analysis of DENV-2 isolates collected from particular countries has been performed using partial or individual genes but only a few studies have examined complete whole-genome sequences collected worldwide. Herein, 50 complete genome sequences of DENV-2 isolates, reported over the past 70 years from 19 different countries, were downloaded from GenBank. Phylogenetic analysis was conducted and evolutionary distances of the 50 DENV-2 isolates were determined using maximum likelihood (ML) trees or Bayesian phylogenetic analysis created from complete genome nucleotide (nt) and amino acid (aa) sequences or individual gene sequences. The results showed that all DENV-2 isolates fell into seven main groups containing five previously defined genotypes. A Cosmopolitan genotype showed further division into three groups (C-I, C-II, and C-III) with the C-I group containing two subgroups (C-IA and C-IB). Comparison of the aa sequences showed specific mutations among the various groups of DENV-2 isolates. A maximum number of aa mutations was observed in the NS5 gene, followed by the NS2A, NS3 and NS1 genes, while the smallest number of aa substitutions was recorded in the capsid gene, followed by the PrM/M, NS4A, and NS4B genes. Maximum evolutionary distances were found in the NS2A gene, followed by the NS4A and NS4B genes. Based on these results, we propose that genotyping of DENV-2 isolates in future studies should be performed on entire genome sequences in order to gain a complete understanding of the evolution of various isolates reported from different geographical locations around the world.     

Hepatitis C virus protease gene diversity in patients coinfected with human immunodeficiency virus

The clonal variability of the hepatitis C virus (HCV) protease gene in 24 individuals with HCV genotypes 1a, 1b, 2b, and 3a who were coinfected with the human immunodeficiency virus was evaluated. Within-genotype variability at the nucleotide and amino acid levels ranged from 6.5 to 8.6% and 2.2 to 3.8%, respectively. After adjustments were made for correlation of intrapatient clonal variation, mixed-model analysis indicated that nucleotide and amino acid variability among patients with different genotypes did not differ significantly. However, within individual patients, clonal variability differed by up to 5.3% and 5.8% at the nucleotide and amino acid levels, respectively, and genotype 1a had significantly greater nucleotide variability than other genotypes (P = 0.01). Significant variability exists within HCV protease gene variants at the patient level and could affect the effectiveness of HCV protease inhibitors.     

Genome Sequencing and Phylogenetic Analysis of 39 Human Parainfluenza Virus Type 1 Strains Isolated from 1997–2010

Thirty-nine human parainfluenza type 1 (HPIV-1) genomes were sequenced from samples collected in Milwaukee, Wisconsin from 1997–2010. Following sequencing, phylogenetic analyses of these sequences plus any publicly available HPIV-1 sequences (from GenBank) were performed. Phylogenetic analysis of the whole genomes, as well as individual genes, revealed that the current HPIV-1 viruses group into three different clades. Previous evolutionary studies of HPIV-1 in Milwaukee revealed that there were two genotypes of HPIV-1 co-circulating in 1991 (previously described as HPIV-1 genotypes C and D). The current study reveals that there are still two different HPIV-1 viruses co-circulating in Milwaukee; however, both groups of HPIV-1 viruses are derived from genotype C indicating that genotype D may no longer be in circulation in Milwaukee. Analyses of genetic diversity indicate that while most of the genome is under purifying selection some regions of the genome are more tolerant of mutation. In the 40 HPIV-1 genomes sequenced in this study, the nucleotide sequence of the L gene is the most conserved while the sequence of the P gene is the most variable. Over the entire protein coding region of the genome, 81 variable amino acid residues were observed and as with nucleotide diversity, the P protein seemed to be the most tolerant of mutation (and contains the greatest proportion of non-synonymous to synonymous substitutions) while the M protein appears to be the least tolerant of amino acid substitution.     

Evolutionary coupling between the deleteriousness of gene mutations and the amount of non-coding sequences

The phenotypic effects of random mutations depend on both the architecture of the genome and the gene-trait relationships. Both levels thus play a key role in the mutational variability of the phenotype, and hence in the long-term evolutionary success of the lineage. Here, by simulating the evolution of organisms with flexible genomes, we show that the need for an appropriate phenotypic variability induces a relationship between the deleteriousness of gene mutations and the quantity of non-coding sequences maintained in the genome. The more deleterious the gene mutations, the shorter the intergenic sequences. Indeed, in a shorter genome, fewer genes are affected by rearrangements (duplications, deletions, inversions, translocations) at each replication, which compensates for the higher impact of each gene mutation. This spontaneous adjustment of genome structure allows the organisms to retain the same average fitness loss per replication, despite the higher impact of single gene mutations. These results show how evolution can generate unexpected couplings between distinct organization levels.     

拉米夫定联合阿德福韦酯治疗前后乙型肝炎病毒的全基因组序列分析

目的探究拉米夫定治疗反弹后联合阿德福韦酯治疗前后乙型肝炎全基因组序列变化。方法分别提取服用拉米夫定治疗24周反弹后和阿德福韦酯辅助治疗24周后的患者2份血清病毒核酸,用聚合酶链反应扩增核酸后进行全基因组测序分析。结果测序结果显示,共计有29个氨基酸发生了突变,其中,S区突变点有5个（17．2％）,C区突变点有12个（41．3％）,P区突变点有6个（20．6％）,X区突变点有6个（20．6％）,其中P区与拉米夫定的相关位点173和204位点发生了突变翻转,但服用阿德福韦后出现了与之相关的突变位点（181、214、236和237位点）。结论核苷酸药物的使用和HBV基因耐药突变密切相关,定期检测HBV基因突变对于合理使用核苷酸药物具有重要意义。     

小鼠Lrh-1基因与排卵数间相关性分析

目的 探讨肝受体类似物-1( liver receptor homolog-1,Lrh-1;NR5A2)基因序列差异与小鼠排卵数性状间关系.方法 根据GenBank中NM_001159769序列设计了5对引物扩增Lrh-1基因CDS区,用单链构象多态性(Single-strand conformation polymorphism,SSCP)分析小鼠Lrh-1基因序列差异,分析Lrh-1基因与小鼠排卵数间关系.结果 ①在5对引物中只有引物P2和P5扩增产物存在多态性,引物P2扩增产物存在两种基因型AA和AB型,其中AB型发生(C674A)突变,这种突变导致编码的140位氨基酸由谷氨酰胺变为赖氨酸.②引物P5扩增产物存在三种多态类型SSCP-1、SSCP-2和SSCP-3,SSCP-1和SSCP-3型核酸序列比SSCP-2型少25个碱基,SSCP-3型除缺失区段外,其他区域碱基序列与SSCP-2型有5处碱基突变,而与SSCP-1型核酸序列有34处碱基差异,经BLAST分析,SSCP-2型与NM_001159769序列相似,SSCP-1型与NM_001159769序列相差较多,而与NG_012313.1序列相似,通过氨基酸序列比对分析,SSCP-3型1652处的A→G的突变导致氨基酸由丙氨酸变为苏氨酸,1678位G→C的突变使原密码子TAC(酪氨酸Y)变为TAG(终止子).③引物P2 AA型(27.27±8.52)与AB型小鼠的平均排卵数(25.92±11.73)间差异无显著性(P＞0.05);引物P5 SSCP-2型平均排卵数(35.00±14.58)显著高于SSCP-3型平均排卵数(19.50±7.94) (P ＜0.05),SSCP-1型(26.2±8.18)与SSCP-2型、SSCP-3型平均排卵数间差异无显著性(P＞0.05).结论 明确了Lrh-1基因序列差异与小鼠排卵数性状间关系,为深入研究Lrh-1基因对动物生殖调控机理提供依据.     

The variability of the hepatitis B virus genome: statistical analysis and biological implications

A statistical analysis of the nucleotide sequence variability in 14 published hepatitis B virus (HBV) genomes was carried out using parametric and nonparametric methods. A parametric statistical model revealed that the different regions of the genome differed significantly in their variability. The conclusion was supported by a nonparametric kernel-density model of the HBV genome. Genes S, C, and P, region X, the precore region, and the pre-S2/pre-S1 regions were ranked in order of increasing variability. In many instances, conserved regions of the genome identified with sequences of known function in HBV biology. However, other characterized regions (such as pre-S) showed much variability despite the involvement of their encoded peptides in specific functions. Point mutations that may result in the formation of stop codons and amino acid changes may affect the clinical picture of HBV infection and may be reflected in atypical serological patterns.