Multiple Nucleotide Changes Research Articles

Conserved Pattern and Potential Role of Recurrent Deletions in SARS-CoV-2 Evolution.

ABSTRACTSARS-CoV-2 continues adapting to human hosts during the current worldwide pandemic since 2019. This virus evolves through multiple means, such as single nucleotide mutations and structural variations, which has brought great difficulty to disease prevention and control of COVID-19. Structural variation, including multiple nucleotide changes like insertions and deletions, has a greater impact relative to single nucleotide mutation on both genome structures and protein functions. In this study, we found that deletion occurred frequently in not only SARS-CoV-2 but also in other SARS-related coronaviruses. These deletions showed obvious location bias and formed 45 recurrent deletion regions in the viral genome. Some of these deletions showed proliferation advantages, including four high-frequency deletions (nsp6 Δ106-109, S Δ69-70, S Δ144, and Δ28271) that were detected in around 50% of SARS-CoV-2 genomes and other 19 median-frequency deletions. In addition, the association between deletions and the WHO reported variants of concern (VOC) and variants of interest (VOI) of SARS-CoV-2 indicated that these variants had a unique combination of deletion patterns. In the spike (S) protein, the deletions in SARS-CoV-2 were mainly in the N-terminal domain. Some deletions, such as S Δ144/145 and S Δ243-244, have been confirmed to block the binding sites of neutralizing antibodies. Overall, this study revealed a conservative regional pattern and the potential effect of some deletions in SARS-CoV-2 over the whole genome, providing important evidence for potential epidemic control and vaccine development.IMPORTANCE Mutations in SARS-CoV-2 were studied extensively, while only the structure variations on the spike protein were discussed well in previous studies. To study the role of structural variations in virus evolution, we described the distribution of structure variations on the whole genome. Conserved patterns were found of deletions among SARS-CoV-2, SARS-CoV-2-like, and SARS-CoV-like viruses. There were 45 recurrent deletion regions (RDRs) in SARS-CoV-2 generated through the integration of deleted positions. In these regions, four high-frequency deletions parallelly appeared in multiple strains. Furthermore, in the spike protein, the deletions in SARS-CoV-2 were mainly in the N-terminal domain, blocking the binding sites of some neutralizing antibodies, while the structural variations in SARS-related coronavirus were mainly in the N-terminal domain and receptor binding domain. The receptor binding domain is highly related to hosting recognition. The deletions in the receptor binding domain may play a role in host adaption.

Rapid screening for SARS-CoV-2 variants of concern in clinical and environmental samples using nested RT-PCR assays targeting key mutations of the spike protein

New SARS-CoV-2 mutations are constantly emerging, raising concerns of increased transmissibility, virulence or escape from host immune response.We describe a nested RT-PCR assay (~1500 bps) to detect multiple nucleotide changes resulting in key spike protein mutations distinctive of the major known circulating SARS-CoV-2 variants, including the three Variants of Concern (VOCs) 20I/501Y.V1 (United Kingdom), 20H/501Y.V2 (South Africa), and 20 J/501Y.V3 (Brazil), as well as the 20E.EU1 variant (Spain), the CAL.20C recently identified in California, and the mink-associated variant (GR, lineage B.1.1.298). Prior to application to field samples, the discriminatory potential of this PCR assay was explored using GISAID and Nextclade. To extend variant detection to challenging matrices such as sewage, where the amplification of long fragments is problematic, two short nested RT-PCR assays (~300 bps) were also designed, targeting portions of the region spanned by the long nested assay.The three newly-designed assays were then tested on field samples, including 31 clinical samples (7 fully-sequenced swab samples, and 24 uncharacterized ones) and 34 urban wastewater samples, some of which collected in areas where circulation of VOCs had been reported.The long assay successfully amplified 29 of the 31 swabs (93%), allowing the correct identification of variants 20I/501Y.V1 and 20E.EU1 present in the panel of previously characterized samples. The Spanish variant was detected in 14/24 of the uncharacterized samples as well. The sequences obtained using the short assays were consistent with those obtained with the long assay. Mutations characteristic of VOCs (UK and Brazilian variant) and of other variant (Spanish) were detected in sewage samples. To our knowledge, this is the first evidence of the presence of sequences harboring key mutations of 20I/501Y.V1 and 20 J/501Y.V3 in urban wastewaters, highlighting the potential contribution of wastewater surveillance to explore SARS-CoV-2 diversity.The developed nested RT-PCR assays can be used as an initial rapid screening test to select clinical samples containing mutations of interest. This can speed up diagnosis and optimize resources since it allows full genome sequencing to be done only on clinically relevant specimens. The assays can be also employed for a rapid and cost-effective detection of VOCs or other variants in sewage for the purposes of wastewater-based epidemiology. The approach proposed here can be used to better understand SARS-CoV-2 variant diversity, geographic distribution and impact worldwide.

Extra base hits: Widespread empirical support for instantaneous multiple-nucleotide changes

Despite many attempts to introduce evolutionary models that permit substitutions to instantly alter more than one nucleotide in a codon, the prevailing wisdom remains that such changes are rare and generally negligible or are reflective of non-biological artifacts, such as alignment errors. Codon models continue to posit that only single nucleotide change have non-zero rates. Here, we develop and test a simple hierarchy of codon-substitution models with non-zero evolutionary rates for only one-nucleotide (1H), one- and two-nucleotide (2H), or any (3H) codon substitutions. Using over 42, 000 empirical alignments, we find widespread statistical support for multiple hits: 61% of alignments prefer models with 2H allowed, and 23%-with 3H allowed. Analyses of simulated data suggest that these results are not likely to be due to simple artifacts such as model misspecification or alignment errors. Further modeling reveals that synonymous codon island jumping among codons encoding serine, especially along short branches, contributes significantly to this 3H signal. While serine codons were prominently involved in multiple-hit substitutions, there were other common exchanges contributing to better model fit. It appears that a small subset of sites in most alignments have unusual evolutionary dynamics not well explained by existing model formalisms, and that commonly estimated quantities, such as dN/dS ratios may be biased by model misspecification. Our findings highlight the need for continued evaluation of assumptions underlying workhorse evolutionary models and subsequent evolutionary inference techniques. We provide a software implementation for evolutionary biologists to assess the potential impact of extra base hits in their data in the HyPhy package and in the Datamonkey.org server.

In Silico Analysis of the Biodiversity and Conservation Status of Mitochondrial Cytochrome C Oxidase Subunit 1 (CO1) Gene of Taenia multiceps.

Taenia multiceps resides in the small intestine of carnivores such as dogs, foxes, woles, jackals, while Coenurus cerebralis which is the larval form usually settle in the central nervous system and spinal cord of intermediate hosts like sheep, cattle and goats. The aim of this study was to analyze the haplotype diversity, genetic variation and population structure of the mt-CO1 gene sequences of sheep isolates of T. multiceps had been submitted to GenBank from different countries. A total of 102 sequences from the mt-CO1 gene fragment belong to T. multiceps sheep isolates were used for bioinformatic analyses. Haplotype analysis, phylogenetic analysis and diversity, neutrality, fixation and gene flow analyses were applied to the sequences. As a result, 20 haplotypes together with different multiple nucleotide changes were determined after the sequence analysis. Trimmed fragment length was 337bp hereby 19 polymorphic areas, 12 of which were parsimony informative, were identified, and any insertion-deletion was found. The number of mutations between major haplotypes and the others range from one to nine. The highest (0.72) genetic differentiation (Fst value) was observed between Turkey and Egypt populations while the lowest (- 0.22) was reported from Greece. These findings are important in terms of showing the diversity of nucleotide variation in T. multiceps sheep isolates. This study serves as the basis for future large-scale studies on T. multiceps worldwide epidemiology, bioecology, geographic distribution and population structure.

Multiple Changes Underlie Allelic Divergence of CUP2 Between Saccharomyces Species.

Under the model of micromutationism, phenotypic divergence between species is caused by accumulation of many small-effect changes. While mapping the causal changes to single nucleotide resolution could be difficult for diverged species, genetic dissection via chimeric constructs allows us to evaluate whether a large-effect gene is composed of many small-effect nucleotide changes. In a previously described non-complementation screen, we found an allele difference of CUP2, a copper-binding transcription factor, underlies divergence in copper resistance between Saccharomyces cerevisiae and S. uvarum. Here, we tested whether the allele effect of CUP2 was caused by multiple nucleotide changes. By analyzing chimeric constructs containing four separate regions in the CUP2 gene, including its distal promoter, proximal promoter, DNA binding domain and transcriptional activation domain, we found that all four regions of the S. cerevisiae allele conferred copper resistance, with the proximal promoter showing the largest effect, and that both additive and epistatic effects are likely involved. These findings support a model of multiple changes underlying evolution and suggest an important role of both protein coding and cis-regulatory changes in evolution.

Molecular characterization and haplotypes of sheep and goat isolates of Cysticercus tenuicollis in Turkey.

Taenia hydatigena, is a common parasite, lived mostly in dogs and wild carnivores in its mature stage, and the larvae, Cysticercus tenuicollis, is found on ruminants and pigs. The aim of the current study was to determine the genetic diversity in 20 isolates of the sheep and goats. After the isolation of total genomic DNA from C. tenuicollis isolates, genetic characterization of the mitochondrial NADH dehydrogenase subunit 1 gene region was amplified using specific JB11-JB12 primers in PCR and the PCR products were sequenced and haplotype and genetic diversity analyses were utilized. As a result, multiple nucleotide changes were determined in the sequence analyses of the isolates leading to detection of 16 and 15 different haplotypes in sheep and goat samples, respectively. These findings are important in terms of showing the diversity of nucleotide variation in C. tenuicollis in Turkey.

Improved inference of site-specific positive selection under a generalized parametric codon model when there are multinucleotide mutations and multiple nonsynonymous rates

BackgroundAn excess of nonsynonymous substitutions, over neutrality, is considered evidence of positive Darwinian selection. Inference for proteins often relies on estimation of the nonsynonymous to synonymous ratio (ω = dN/dS) within a codon model. However, to ease computational difficulties, ω is typically estimated assuming an idealized substitution process where (i) all nonsynonymous substitutions have the same rate (regardless of impact on organism fitness) and (ii) instantaneous double and triple (DT) nucleotide mutations have zero probability (despite evidence that they can occur). It follows that estimates of ω represent an imperfect summary of the intensity of selection, and that tests based on the ω > 1 threshold could be negatively impacted.ResultsWe developed a general-purpose parametric (GPP) modelling framework for codons. This novel approach allows specification of all possible instantaneous codon substitutions, including multiple nonsynonymous rates (MNRs) and instantaneous DT nucleotide changes. Existing codon models are specified as special cases of the GPP model. We use GPP models to implement likelihood ratio tests for ω > 1 that accommodate MNRs and DT mutations. Through both simulation and real data analysis, we find that failure to model MNRs and DT mutations reduces power in some cases and inflates false positives in others. False positives under traditional M2a and M8 models were very sensitive to DT changes. This was exacerbated by the choice of frequency parameterization (GY vs. MG), with rates sometimes > 90% under MG. By including MNRs and DT mutations, accuracy and power was greatly improved under the GPP framework. However, we also find that over-parameterized models can perform less well, and this can contribute to degraded performance of LRTs.ConclusionsWe suggest GPP models should be used alongside traditional codon models. Further, all codon models should be deployed within an experimental design that includes (i) assessing robustness to model assumptions, and (ii) investigation of non-standard behaviour of MLEs. As the goal of every analysis is to avoid false conclusions, more work is needed on model selection methods that consider both the increase in fit engendered by a model parameter and the degree to which that parameter is affected by un-modelled evolutionary processes.

Bacterial DNA polymerases participate in oligonucleotide recombination.

Synthetic single-strand oligonucleotides (oligos) with homology to genomic DNA have proved to be highly effective for constructing designed mutations in targeted genomes, a process referred to as recombineering. The cellular functions important for this type of homologous recombination have yet to be determined. Towards this end, we have identified Escherichia coli functions that process the recombining oligo and affect bacteriophage λ Red-mediated oligo recombination. To determine the nature of oligo processing during recombination, each oligo contained multiple nucleotide changes: a single base change allowing recombinant selection, and silent changes serving as genetic markers to determine the extent of oligo processing during the recombination. Such oligos were often not incorporated into the host chromosome intact; many were partially degraded in the process of recombination. The position and number of these silent nucleotide changes within the oligo strongly affect both oligo processing and recombination frequency. Exonucleases, especially those associated with DNA Polymerases I and III, affect inheritance of the silent nucleotide changes in the oligos. We demonstrate for the first time that the major DNA polymerases (Pol I and Pol III) and DNA ligase are directly involved with oligo recombination.

Superiority of a mechanistic codon substitution model even for protein sequences in Phylogenetic analysis

BackgroundNucleotide and amino acid substitution tendencies are characteristic of each species, organelle, and protein family. Hence, various empirical amino acid substitution rate matrices have needed to be estimated for phylogenetic analysis: JTT, WAG, and LG for nuclear proteins, mtREV for mitochondrial proteins, cpREV10 and cpREV64 for chloroplast-encoded proteins, and FLU for influenza proteins. On the other hand, in a mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the ratio of fixation depending on the type of amino acid replacement, mutation rates and the strength of selective constraint on amino acids can be tailored to each protein family with additional 11 parameters. As a result, in the evolutionary analysis of codon sequences it outperforms codon substitution models equivalent to empirical amino acid substitution matrices. Is it superior even for amino acid sequences, among which synonymous substitutions cannot be identified?ResultsNucleotide mutations are assumed to occur independently of codon positions but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene with a linear function of a given estimate of selective constraints, which were estimated by maximizing the likelihood of an empirical amino acid or codon substitution frequency matrix, each of JTT, WAG, LG, and KHG. It is shown that the mechanistic codon substitution model with the assumption of equal codon usage yields better values of Akaike and Bayesian information criteria for all three phylogenetic trees of mitochondrial, chloroplast, and influenza-A hemagglutinin proteins than the empirical amino acid substitution models with mtREV, cpREV64, and FLU, which were designed specifically for those protein families, respectively. The variation of selective constraint across sites fits the datasets significantly better than variable codon mutation rates, confirming that substitution rate variations across sites detected by amino acid substitution models are caused primarily by the variation of selective constraint against amino acid substitutions rather than the variation of codon mutation rate.ConclusionsThe mechanistic codon substitution model is superior to amino acid substitution models even in the evolutionary analysis of protein sequences.

Advantages of a Mechanistic Codon Substitution Model for Evolutionary Analysis of Protein-Coding Sequences

BackgroundA mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated.ResultsThe present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.

Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes

BackgroundEmpirical substitution matrices represent the average tendencies of substitutions over various protein families by sacrificing gene-level resolution. We develop a codon-based model, in which mutational tendencies of codon, a genetic code, and the strength of selective constraints against amino acid replacements can be tailored to a given gene. First, selective constraints averaged over proteins are estimated by maximizing the likelihood of each 1-PAM matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution matrices. Then, selective constraints specific to given proteins are approximated as a linear function of those estimated from the empirical substitution matrices.ResultsAkaike information criterion (AIC) values indicate that a model allowing multiple nucleotide changes fits the empirical substitution matrices significantly better. Also, the ML estimates of transition-transversion bias obtained from these empirical matrices are not so large as previously estimated. The selective constraints are characteristic of proteins rather than species. However, their relative strengths among amino acid pairs can be approximated not to depend very much on protein families but amino acid pairs, because the present model, in which selective constraints are approximated to be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can provide a good fit to other empirical substitution matrices including cpREV for chloroplast proteins and mtREV for vertebrate mitochondrial proteins.Conclusions/SignificanceThe present codon-based model with the ML estimates of selective constraints and with adjustable mutation rates of nucleotide would be useful as a simple substitution model in ML and Bayesian inferences of molecular phylogenetic trees, and enables us to obtain biologically meaningful information at both nucleotide and amino acid levels from codon and protein sequences.

3P306 多塩基置換を考慮したコドン置換モデルにより評価したアミノ酸置換における選択的制約(生命情報科学-分子進化,第48回日本生物物理学会年会)

3P306 多塩基置換を考慮したコドン置換モデルにより評価したアミノ酸置換における選択的制約(生命情報科学-分子進化,第48回日本生物物理学会年会)

RelB and RelE of Escherichia coli Form a Tight Complex That Represses Transcription via the Ribbon–Helix–Helix Motif in RelB

RelB, the ribbon–helix–helix (RHH) repressor encoded by the relBE toxin–antitoxin locus of Escherichia coli, interacts with RelE and thereby counteracts the mRNA cleavage activity of RelE. In addition, RelB dimers repress the strong relBE promoter and this repression by RelB is enhanced by RelE; that is, RelE functions as a transcriptional co-repressor. RelB is a Lon protease substrate, and Lon is required both for activation of relBE transcription and for activation of the mRNA cleavage activity of RelE. Here we characterize the molecular interactions important for transcriptional control of the relBE model operon. Using an in vivo screen for relB mutants, we identified multiple nucleotide changes that map to important amino acid positions within the DNA-binding domain formed by the N-terminal RHH motif of RelB. Analysis of DNA binding of a subset of these mutant RHH proteins by gel-shift assays, transcriptional fusion assays and a structure model of RelB–DNA revealed amino acid residues making crucial DNA–backbone contacts within the operator (relO) DNA. Mutational and footprinting analyses of relO showed that RelB dimers bind on the same face of the DNA helix and that the RHH motif recognizes four 6-bp repeats within the bipartite binding site. The spacing between each half-site was found to be essential for cooperative interactions between adjacently bound RelB dimers stabilized by the co-repressor RelE. Kinetic and stoichiometric measurements of the interaction between RelB and RelE confirmed that the proteins form a high-affinity complex with a 2:1 stoichiometry. Lon degraded RelB in vitro and degradation was inhibited by RelE, consistent with the proposal that RelE protects RelB from proteolysis by Lon in vivo.

The coat protein of Tomato mosaic virus L11Y is associated with virus-induced chlorosis on infected tobacco plants

The L11Y strain of Tomato mosaic virus (ToMV) causes severe chlorosis on infected tobacco leaves. Sequencing analysis for the genome showed that L11Y contained multiple nucleotide changes and that some led to amino acid substitutions, when compared with that of the common L strain of ToMV. The chimeric virus, which has the CP of L11Y in the context of the L strain RNA genome, caused severe chlorosis on infected tobacco plants, suggesting that the CP of L11Y containing three amino acid changes (E33S, A86T and E97K) was the determinant of the chlorosis. Two of these amino acid changes (A86T and E97K) were associated with the induction of chlorosis when present together in the CP. Severe destruction and deformation of chloroplasts and the formation of discrete dark-staining materials adjacent to chloroplasts were observed with electron microscopy in L11Y-infected plants. Fewer virus particles accumulated in the cytoplasm of L11Y-infected plant cells. The level of accumulation of CP subgenomic RNA and CP in the infected protoplasts was similar between L and L11Y. Fewer virus particles accumulated in L11Y-infected protoplasts, and many of them were shorter-than-full-length.

An Empirical Codon Model for Protein Sequence Evolution

In the past, 2 kinds of Markov models have been considered to describe protein sequence evolution. Codon-level models have been mechanistic with a small number of parameters designed to take into account features, such as transition-transversion bias, codon frequency bias, and synonymous-nonsynonymous amino acid substitution bias. Amino acid models have been empirical, attempting to summarize the replacement patterns observed in large quantities of data and not explicitly considering the distinct factors that shape protein evolution. We have estimated the first empirical codon model (ECM). Previous codon models assume that protein evolution proceeds only by successive single nucleotide substitutions, but our results indicate that model accuracy is significantly improved by incorporating instantaneous doublet and triplet changes. We also find that the affiliations between codons, the amino acid each encodes and the physicochemical properties of the amino acids are main factors driving the process of codon evolution. Neither multiple nucleotide changes nor the strong influence of the genetic code nor amino acids' physicochemical properties form a part of standard mechanistic models and their views of how codon evolution proceeds. We have implemented the ECM for likelihood-based phylogenetic analysis, and an assessment of its ability to describe protein evolution shows that it consistently outperforms comparable mechanistic codon models. We point out the biological interpretation of our ECM and possible consequences for studies of selection.

Estimating the frequency of events that cause multiple-nucleotide changes.

Existing mathematical models of DNA sequence evolution assume that all substitutions derive from point mutations. There is, however, increasing evidence that larger-scale events, involving two or more consecutive sites, may also be important. We describe a model, denoted SDT, that allows for single-nucleotide, doublet, and triplet mutations. Applied to protein-coding DNA, the SDT model allows doublet and triplet mutations to overlap codon boundaries but still permits data to be analyzed using the simplifying assumption of independence of sites. We have implemented the SDT model for maximum-likelihood phylogenetic inference and have applied it to an alignment of mammalian globin sequences and to 258 other protein-coding sequence alignments from the Pandit database. We find the SDT model's inclusion of doublet and triplet mutations to be overwhelmingly successful in giving statistically significant improvements in fit of model to data, indicating that larger-scale mutation events do occur. Distributions of inferred parameter values over all alignments analyzed suggest that these events are far more prevalent than previously thought. Detailed consideration of our results and the absence of any known mechanism causing three adjacent nucleotides to be substituted simultaneously, however, leads us to suggest that the actual evolutionary events occurring may include still-larger-scale events, such as gene conversion, inversion, or recombination, or a series of rapid compensatory changes.

Hepatitis delta virus RNA encoding the large delta antigen cannot sustain replication due to rapid accumulation of mutations associated with RNA editing.

Hepatitis delta virus (HDV) contains two RNA species (HDV-S and HDV-L), which encode the small and large forms of hepatitis delta antigens (S- and L-HDAg), respectively. HDV-L RNA is a result of an RNA editing event occurring at an amber/W site of HDV-S RNA. RNA editing must be regulated to prevent premature and excessive accumulation of HDV-L RNA in the viral life cycle. In this study, we used an RNA transfection procedure to study the replication abilities of HDV-L and HDV-S RNA. While HDV-S led to robust RNA replication, HDV-L could not replicate even after 6 days following transfection. The failure of HDV-L to replicate was not due to insufficient amounts of S-HDAg, as identical results were obtained in a cell line that stably overexpresses S-HDAg. Also, it was not due to possible inhibition by L-HDAg, as HDV-S RNA replication was not affected when both HDV-L and HDV-S RNA were cotransfected. Further, when L-HDAg expression from HDV-L RNA was abolished by site-directed mutagenesis, the mutant HDV-L RNA also failed to replicate. Unexpectedly, when the kinetics of RNA replication was examined daily, HDV-L was found to replicate at a low level at the early time points (1 to 2 days posttransfection) but then lose this capability at later time points. Sequence analysis of the replicated HDV-L RNA at day 1 posttransfection showed that it had undergone multiple nucleotide changes, particularly in the region near the putative promoter region of HDV RNA replication. In contrast, very few mutations were found in HDV-S RNA. These results suggest that the editing at the amber/W site triggers a series of additional mutations which rapidly reduce the replication efficiency of the resultant HDV genome and thus help regulate the amount of HDV-L RNA in infected cells. They also explain why L-HDAg is not produced early in HDV infection, despite the fact that HDV-L RNA is present in the virion.

Functional conservation of the hydrophobic domain of polypeptide 3AB between human rhinovirus and poliovirus

In this study we exchanged portions of the poliovirus type 1 (PV1) hydrophobic domain within the membrane-associated polypeptide 3AB for the analogous sequences from human rhinovirus 14 (HRV14). The sequence exchanges were based upon a previous report in which the 22 amino acid hydrophobic region was subdivided into two domains, I and II, the latter of which was shown to be required for membrane association (J. Biol. Chem. 271 (1996), 26810). Using these divisions, the HRV14 sequences were cloned into the complete poliovirus type 1 cDNA sequence. RNAs transcribed from these cDNAs were transfected into HeLa cell monolayers and used in HeLa cell-free translation/replication assays. The data indicated that 3AB sequences from PV1 and HRV14 are interchangeable; however, the substitutions cause a range of significant RNA replication defects, and in some cases, protein processing defects. Following transfection of RNAs encoding the domain substitutions into HeLa cell monolayers, virus isolates were harvested, and the corresponding viral RNAs were sequenced. The sequence data revealed that for the carboxy-terminal domain substitutions (domain II), multiple nucleotide changes were identified in the first, second, and third positions of different codons. In addition, the data indicated that for one of the PV1/HRV14 chimeras to replicate, compensatory mutations within poliovirus protein 2B may be required.

Conserved subfamilies of the Drosophila HeT-A telomere-specific retrotransposon.

HeT-A, a major component of Drosophila telomeres, is the first retrotransposon proposed to have a vital cellular function. Unlike most retrotransposons, more than half of its genome is noncoding. The 3' end contains > 2.5 kb of noncoding sequence. Copies of HeT-A differ by insertions or deletions and multiple nucleotide changes, which initially led us to conclude that HeT-A noncoding sequences are very fluid. However, we can now report, on the basis of new sequences and further analyses, that most of these differences are due to the existence of a small number of conserved sequence subfamilies, not to extensive sequence change during each transposition event. The high level of sequence conservation within subfamilies suggests that they arise from a small number of replicatively active elements. All HeT-A subfamilies show preservation of two intriguing features. First, segments of extremely A-rich sequence form a distinctive pattern within the 3' noncoding region. Second, there is a strong strand bias of nucleotide composition: The DNA strand running 5' to 3' toward the middle of the chromosome is unusually rich in adenine and unusually poor in guanine. Although not faced with the constraints of coding sequences, the HeT-A 3' noncoding sequence appears to be under other evolutionary constraints, possibly reflecting its roles in the telomeres.

Influence of the 5' noncoding region of hepatitis A virus strain GBM on its growth in different cell lines.

Previous sequence analysis of consecutive passages of the hepatitis A virus (HAV) strain GBM/WT in human embryonic kidney cells (HEK cells), human embryonic lung fibroblasts (HFS cells) and in FRhK-4 cells (foetal rhesus monkey kidney cells) pointed to a host cell dependent cell culture adaptation of GBM/WT in HFS cells involving mutations in the 5' noncoding region (5'NCR). Multiple nucleotide changes occurred in the 5'NCR of the GBM genome after the cell line used for virus passage was changed from HEK cells to HFS cells. In contrast, no mutations in the 5'NCR occurred during the first 20 passages of GBM/WT in FRhK-4 cells. In order to analyse the influence of the 5'NCR on host cell specific adaptation of HAV strain GBM in different cell cultures, GBM/HM175 chimeras were constructed which contained 5'NCRs from different GBM variants by replacing the 5'NCR of the infectious clone pHAV/7. Parallel transfection assays in FRhK-4 and HFS cells, performed with transcripts from the chimeric GBM/HM175 constructs, showed that the 5'NCR of the GBM variant GBM/HFS is essential for virus growth in HFS cells. The GBM/HM175 chimeric RNA, which contained the 5'NCR of GBM/HFS, exclusively, was able to produce infectious virus after transfection of HFS cells. The growth of the different GBM/HM175 chimeras in FRhK-4 cells, in contrast, did not seem to be strongly influenced by a specific sequence of the 5'NCR.