Large Variable Region Research Articles

Novel Sulfolobus Fuselloviruses with Extensive Genomic Variations.

Fuselloviruses are among the most widespread and best-characterized archaeal viruses. They exhibit remarkable diversity, as the list of members of this family is rapidly growing. However, it has yet to be shown how a fuselloviral genome may undergo variation at the levels of both single nucleotides and sequence stretches. Here, we report the isolation and characterization of four novel spindle-shaped viruses, named Sulfolobus spindle-shaped viruses 19 to 22 (SSV19-22), from a hot spring in the Philippines. SSV19 is a member of the genus Alphafusellovirus, whereas SSV20-22 belong to the genus Betafusellovirus The genomes of SSV20-SSV22 are identical except for the presence of two large variable regions, as well as numerous sites of single-nucleotide polymorphisms (SNPs) unevenly distributed throughout the genomes and enriched in certain regions, including the gene encoding the putative end filament protein VP4. We show that coinfection of the host with SSV20 and SSV22 led to the formation of an SSV21-like virus, presumably through homologous recombination. In addition, large numbers of SNPs were identified in DNA sequences retrieved by PCR amplification targeting the SSV20-22 vp4 gene from the original enrichment culture, indicating the enormous diversity of SSV20-22-like viruses in the environment. The high variability of VP4 is consistent with its potential role in host recognition and binding by the virus.IMPORTANCE How a virus survives in the arms race with its host is an intriguing question. In this study, we isolated and characterized four novel fuselloviruses, named Sulfolobus spindle-shaped viruses 19 to 22 (SSV19-22). Interestingly, SSV20-22 differ primarily in two genomic regions and are apparently convertible through homologous recombination during coinfection. Moreover, sites of single-nucleotide polymorphism (SNP) were identified throughout the genomes of SSV20-22 and, notably, enriched in certain regions, including the gene encoding the putative end filament protein VP4, which is believed to be involved in host recognition and binding by the virus.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

Monthly maps of optimally interpolated in situ hydrography in the North Sea from 1948 to 2013

We used a statistical model mainly based on optimal or Gauss–Markov interpolation (OI) to produce monthly maps of in situ hydrography in the North Sea through 66years (from 1948 to 2013) with a resolution of 0.2°×0.2°. The in situ hydrography data are approximately 430,000 hydrographic profiles and aggregated thermosalinograph data from various international databases. Duplicates, outliers and vertical density instabilities were removed. Regions with poor OI estimates were replaced with a harmonic reconstruction arising from the most reliable OI estimates. Adjustments for vertical density stability were based on the standards of the World Ocean Atlas. We mapped at 54 depth levels through the water column, focusing on surface and bottom hydrography because this type of map is of particular interest for ecosystem and fisheries research. Average OI temperature and salinity expected errors at the surface are 0.3°C and 0.1, respectively. OI errors decrease with depth following decrement of signal and noise variances and apparently independent of the data amount (indicating a good data coverage). Alternative error estimates were obtained with the Median Absolute Deviation between our hydrography estimates and time series excluded from the analysis and are on average 0.3°C and 0.1 salinity units. While our product seems limited for analysis of variability on monthly and seasonal time scales, particularly in the regions of large variability, it is suitable for studies of inter-annual and decadal variability. A comparison with two alternative analyses (KLIWAS and SODA) is discussed. As direct application of our results, we present a new hydrographic climatology of the North Sea at various depths with an improved effective resolution.

Computational Refinement and Validation Protocol for Proteins with Large Variable Regions Applied to Model HIV Env Spike in CD4 and 17b Bound State

Envelope glycoprotein gp120 of HIV-1 possesses several variable regions; their precise structure has been difficult to establish. We report a new model of gp120, in complex with antibodies CD4 and 17b, complete with its variable regions. The model wasproduced by a computational protocol that uses cryo-electron microscopy (EM) maps, atomic-resolution structures of the core, and information on binding interactions. Our model has excellent fitwith EMD: 5020, is stereochemically and energetically favorable, and has the expected binding interfaces. Comparison of the ternary arrangement of the loops in this model with those bound to PGT122 and PGV04 suggested a possible motion ofthe V1V2 away from the CCR5 binding site andtoward CD4. Our study also revealed that theCD4-bound state of the V1V2 loop is not optimalfor gp120 bound with several neutralizing antibodies.

P60 Study of carotid doppler in patients with systemic lupus erythematosus and correlation with clinical variables

Different chemical reagents were used to study the tertiary structure of yeast tRNASer, a tRNA with a large variable region: ethylnitrosourea, which alkylates the phosphate groups; dimethylsulphate, which methylates N-7 of guanosine and N-3 of cytosine; and diethylpyrocarbonate, which modifies N-7 of adenine. The non-reactivity of N-3 of cytidine 47:1, 47:6, 47:7 and 47:8 and the reactivity of cytidine 47:3 confirms the existence of a variable stem of four base-pairs and a short variable loop of three residues. For the N-7 positions in purines, accessible residues are G1, G10, Gm18, G19, G30, I34, G35, A36, i6A37, G45, G47, G47:5, G47:9 and G73. The protection of N-7 atoms of residues G9, G15, A21, A22 and G47:9 reflects the tertiary folding. Strong phosphate protection was observed for P8 to P11, P20:1 to P22, P48 to P50 and for P59 and P60. A model was built on a PS300 graphic system on the basis of these data and its stereochemistry refined. While trying to keep most tertiary interactions, we adapted the tertiary folding of the known structures of tRNAAsp and tRNAPhe to the present sequence and solution data. The resulting model has the variable arm not far from the plane of the common L-shaped structure. A generalization of this model to other tRNAs with large variable regions is discussed.

Demonstration of the shorter flagellin (flaA) gene of urease-positive thermophilicCampylobacter isolated from the natural environment in Northern Ireland

The PCR amplicons (about 1450 bp in length) of flaA gene fragments of 11 isolates of urease-positive thermophilic Campylobacter (UPTC) isolated from the natural environment not including wild birds in Northern Ireland were demonstrated to be shorter than those of C. jejuni 81116 and six isolates of C. jejuni and C. coli (about 1700 bp) isolated in Northern Ireland and Japan. When the nucleotide lengths of the possible open reading frame (ORF) of the flaA genes were determined, those from the 11 UPTC isolates were estimated to be 1464-1503 bp, and those from the six C. jejuni and C. coli isolates and C. jejuni 81116 strain to be 1716-1728 bp. Nucleotide sequence and deduced amino acid sequence alignments of the possible ORFs demonstrated that the ORFs from the 11 UPTC isolates lack about 80 amino acid residues, mainly from the approximate residue numbers 390-470 of the large variable region in the flaA protein of the seven isolates of C. jejuni and C. coli, and do not have any internal termination codons. High amino acid sequence similarity of both amino- and carboxy-termini of the ORFs of the flaA gene was demonstrated between the 11 isolates of UPTC and the 7 isolates of C. jejuni and C. coli. The 11 UPTC isolates examined were strongly suggested to possess a shorter flaA gene without any internal termination codons.

A Stratified Diagnosis of the Indian Monsoon—Eurasian Snow Cover Relationship

Abstract This study is aimed at reconciling inconsistencies between modeling studies of the Indian monsoon–Eurasian snow cover relationship, which support the Blanford hypothesis, and recent diagnostic studies, which have largely failed to show its existence. Recently released version 2 NOAA/NESDIS satellite-based retrievals of snow cover are used. A focus is given to diagnosing (a) spatial and temporal complexity in the Eurasian snow cover distribution, (b) the role of ENSO in modulating the snow cover–all-India rainfall (AIR) relationship, and (c) the spatial distribution of the rainfall association within India. With these goals in mind, the climatological distribution of snow cover and its interannual variability are first assessed. While snow cover is largest in northern Eurasia, its variability is most pronounced in southwestern (SW) Asia between the Black and Caspian Seas, and across the northern Indian, Himalayan, and Tibetan Plateau (NIHT) regions. In regions of large variability, and for Eurasia...

Molecular cloning, nucleotide sequencing and characterization of the flagellin gene from isolates of urease-positive thermophilic Campylobacter

A primer pair which was expected to generate an amplicon of the estimated size (approximately 1700 base pair (bp)) of the flaA gene for Campylobacter jejuni amplified products of approximately 1450 bp for 33 of the 44 isolates of urease-positive thermophilic Campylobacter (UPTC). The primer pair, however, failed to amplify fragments for 11 isolates of UPTC, for all of the 12 isolates of urease-negative C. lari and for one isolate of C. coli. Nevertheless, it successfully amplified fragments of approximately 1700 bp for five isolates of C. jejuni and for nine isolates of C. coli. Thus, the fragments of the flaA gene of UPTC were shorter than those of C. jejuni and C. coli. After PCR amplification and nucleotide sequencing of the flaA genes from five UPTC NCTC isolates, the putative open reading frames (ORFs) were found to range from 1461 to 1479 bp. The amino acid and nucleotide sequence alignments demonstrated that the PCR clones contained the flaA gene; however, our data indicated that this locus was markedly shorter in the UPTC organisms examined, as they were approximately 85 amino acid residues shorter, mainly corresponding to approximate residue numbers 390–470 of the large variable region of C. jejuni 81116. Heterogeneity was indicated in the molecular mass of the flagellin purified from the isolates examined. Flagellin of UPTC was demonstrated to be genotypically and phenotypically smaller than those of C. jejuni.

The Dependence of Ensemble Dispersion on Analysis–Forecast Systems: Implications to Short-Range Ensemble Forecasting of Precipitation

Abstract The impact of differences in analysis–forecast systems on dispersion of an ensemble forecast is examined for a case of cyclogenesis. Changes in the dispersion properties between two 25-member ensemble forecasts with different cumulus parameterization schemes and different initial analyses are compared. The statistical significance of the changes is assessed. Error growth due to initial condition uncertainty depends significantly on the analysis–forecast system. Quantitative precipitation forecasts and probabilistic quantitative precipitation forecasts are extremely sensitive to the specification of physical parameterizations in the model. Regions of large variability tend to coincide with a high likelihood of parameterized convection. Analysis of other model fields suggests that those with relatively large energy in the mesoscale also exhibit highly significant differences in dispersion. The results presented here provide evidence that the combined effect of uncertainties in model physics and the...

The peculiar architectural framework of tRNASec is fully recognized by yeast AspRS.

The wild-type transcript of Escherichia coli tRNASec, characterized by a peculiar core architecture and a large variable region, was shown to be aspartylatable by yeast AspRS. Similar activities were found for tRNASec mutants with methionine, leucine, and tryptophan anticodons. The charging efficiency of these molecules was found comparable to that of a minihelix derived from tRNAAsp and is accounted for by the presence of the discriminator residue G73, which is a major aspartate identity determinant. Introducing the aspartate identity elements from the anticodon loop (G34, U35, C36, C38) into tRNASec transforms this molecule into an aspartate acceptor with kinetic properties identical to tRNAAsp. Expression of the aspartate identity set in tRNASec is independent of the size of its variable region. The functional study was completed by footprinting experiments with four different nucleases as structural probes. Protection patterns by AspRS of transplanted tRNASec and tRNAAsp were found similar. They are modified, particularly in the anticodon loop, upon changing the aspartate anticodon into that of methionine. Altogether, it appears that recognition of a tRNA by AspRS is more governed by the presence of the aspartate identity set than by the structural framework that carries this set.

A genetic and functional analysis of the unusually large variable region in the M.AluI DNA-(cytosine C5)-methyltransferase.

The M.AluI DNA-(cytosine C5)-methyltransferase (5mC methylase) acts on the sequence 5'-AGCT-3'. The amino acid sequences of known 5mC methylases contain ten conserved motifs, with a variable region between Motifs VIII and IX that contains one or more "target-recognizing domains" (TRDs) responsible for DNA sequence specificity. Monospecific 5mC methylases are believed to have only one TRD, while multispecific 5mC methylases have as many as five. M.AluI has the second-largest variable region of all known 5mC methylases, and sequence analysis reveals five candidate TRDs. In testing whether M.AluI is in fact monospecific it was found that AGCT methylation represents only 80-90% of the methylating activity of this enzyme, while control experiments with the enzyme M.HhaI gave no unexplained activity. Because individual TRDs can be deleted from multispecific methylases without general loss of activity, a series of insertion and deletion mutants of the M.AluI variable region were prepared. All deletions that removed more than single amino acids from the variable region caused significant loss of activity; a sensitive in vivo assay for methylase activity based on McrBC restriction suggested that the central portion of the variable region is particularly important. In some cases, multispecific methylases can accommodate a TRD from another multispecific methylase, thereby acquiring an additional specificity. When TRDs were moved from a multispecific methylase into two different locations in the variable region of M.AluI, all hybrid enzymes had greatly reduced activity and no new specificities. M.AluI thus behaves in most respects as a monospecific methylase despite the remarkable size of its variable region.

Selection of aminoacyl-tRNAs at sense codons: the size of the tRNA variable loop determines whether the immediate 3' nucleotide to the codon has a context effect.

Codon context can affect translational efficiency by several molecular mechanisms. The base stacking interactions between a codon-anticodon complex and the neighboring nucleotide immediately 3' can facilitate translation by amber suppressors and the tRNA structure is also known to modulate the sensitivity to context. In this study the relative rates of aminoacyl-tRNA selection were measured at four sense codons (UGG, CUC, UUC and UCA), in all four 3' nucleotide contexts, through direct competition with a programmed frameshift at a site derived from the release factor 2 gene. Two codons (UGG and UUC) are read by tRNAs with small variable regions and their rates of aminoacyl-tRNA selection correlated with the potential base stacking strength of the 3' neighboring nucleotide. The other two codons (CUC and UCA) are read by tRNAs with large variable regions and the rate of selection of the aminoacyl-tRNAs in these cases varied little among the four contexts. Re-examination of published data on amber suppression also revealed an inverse correlation between context sensitivity and the size of the variable region. Collectively the data suggest that a large variable loop in a tRNA decreases the influence of the 3' context on tRNA selection, probably by strengthening tRNA-ribosomal interactions.

Secondary structure at the 3′ terminal region of RNA coliphages: Comparison with tRNA

Secondary structure models for the 3′ non-coding region of the four groups of coliphage RNA are proposed based on comparative sequence analysis and on previously published data on the sensitivity of nucleotides in MS2 RNA to chemical modification and enzymes. We report the following observations. (1) In contrast to the coding regions, the structure at the 3′ terminus is characterized by stable regular helices. We note the occurrence of the loop sequences 5′-GUUCGC and 5′-CGAAAG, that are reported to confer exceptional stability to stem structures. These features are probably present to promote the segregation of mother and daughter strands during replication. (2) Comparison of homologous helices indicates that only those base pair substitutions are allowed that maintain the thermodynamic stability. (3) We have compared the structure of phage RNA with tRNA. Overall similarity is low, but one common element may exist. It is a quasi-continuous helix of 12 basepairs that could be the equivalent of the 12 basepair long coaxially stacked helix, formed by the TψC arm and the aminoacyl acceptor arm in tRNA. As in tRNA, this structure element starts after the fourth nucleotide from the 3′ end. (4) Phage RNA contains a large variable region of about 35 nucleotides bulging out from the quasi-continuous helix. We speculate that the variable loop in present-day tRNA could be the remnant of the variable region found in phage RNA. The variable region contains overlapping binding sites for the replicase enzyme and the maturation protein. This common binding site may serve as a switch from replication to packaging.

Solution structure of a tRNA with a large variable region: Yeast tRNA Ser

Different chemical reagents were used to study the tertiary structure of yeast tRNA Ser, a tRNA with a large variable region: ethylnitrosourea, which alkylates the phosphate groups; dimethylsulphate, which methylates N-7 of guanosine and N-3 of cytosine; and diethylpyrocarbonate, which modifies N-7 of adenine. The non-reactivity of N-3 of cytidine 47:1, 47:6, 47:7 and 47:8 and the reactivity of cytidine 47:3 confirms the existence of a variable stem of four base-pairs and a short variable loop of three residues. For the N-7 positions in purines, accessible residues are G1, G10, G m18, G19, G30, I34, G35, A36, i 6A37, G45, G47, G47:5, G47:9 and G73. The protection of N-7 atoms of residues G9, G15, A21, A22 and G47:9 reflects the tertiary folding. Strong phosphate protection was observed for P8 to P11, P20:1 to P22, P48 to P50 and for P59 and P60. A model was built on a PS300 graphic system on the basis of these data and its stereochemistry refined. While trying to keep most tertiary interactions, we adapted the tertiary folding of the known structures of tRNA Asp and tRNA Phe to the present sequence and solution data. The resulting model has the variable arm not far from the plane of the common L-shaped structure. A generalization of this model to other tRNAs with large variable regions is discussed.

The avian retrovirus env gene family: molecular analysis of host range and antigenic variants

The nucleotide sequence of the env gp85-coding domain from two avian sarcoma and leukosis retrovirus isolates was determined to identify host range and antigenic determinants. The predicted amino acid sequence of gp85 from a subgroup D virus isolate of the Schmidt-Ruppin strain of Rous sarcoma virus was compared with the previously reported sequences of subgroup A, B, C, and E avian sarcoma and leukosis retroviruses. Subgroup D viruses are closely related to the subgroup B viruses but have an extended host range that includes the ability to penetrate certain mammalian cells. There are 27 amino acid differences shared between the subgroup D sequence and three subgroup B sequences. At 16 of these sites, the subgroup D sequence is identical to the sequence of one or more of the other subgroup viruses (A, C, and E). The remaining 11 sites are specific to subgroup D and show some clustering in the two large variable regions that are thought to be major determinants of host range. Biological analysis of recombinant viruses containing a dominant selectable marker confirmed the role of the gp85-coding domain in determining the host range of the subgroup D virus in the infection of mammalian cells. We also compared the sequence of the gp85-coding domain from two subgroup A viruses, Rous-associated virus type 1 and a subgroup A virus of the Schmidt-Ruppin strain of Rous sarcoma virus. The comparison revealed 24 nonconservative amino acid changes, of which 6 result in changes in potential glycosylation sites. The positions of 10 amino acid differences are coincident with the positions of 10 differences found between two subgroup B virus env gene sequences. These 10 sites identify seven domains in the sequence which may constitute determinants of type-specific antigenicity. Using a molecular recombinant, we demonstrated that type-specific neutralization of two subgroup A viruses was associated with the gp85-coding domain of the virus.

Evolution of pandemic influenza virus strains.

It was found by the technique of molecular hybridization that the pandemic influenza virus strains of 1957 (H2N2) and 1968 (H3N2) evolved by reassortment of RNA segments from the foregoing pandemic strains, replacing four genes including those coding for haemagglutinin (HA) and neuraminidase (NA) (1957) or only that coding for the haemagglutinin (1968), respectively. The earlier pandemic strains from 1918-9 (as deduced from the swine influenza strain (Hsw1N1), which is assumed to be a survivor of the Spanish influenza), from 1933-4 (H0N1), and from 1947 (H1N1) were derived from each other through a number of point mutations only. The Russian strain from 1977 (H1N1) is genetically almost identical with the FW strain from 1950 (H1N1). In contrast to the conserved genes coding for the internal viral proteins, the genes coding for the viral surface glycoproteins consist of a relatively small highly conserved part which presumably is responsible for the functional integrity of the gene products, and the relatively large variable region which probably determines the antigenic property of their gene products. This structural feature of the surface glycoprotein genes explains the fast antigenic drift by mutations in the variable region and selection by the immune response of the host.

The genes coding for the surface glycoproteins of influenza A viruses contain a small conserved and a large variable region

The variable genes of influenza A viruses coding for hemagglutinin and neuraminidase consist of a relatively small highly conserved region, which is presumably involved in the functional integrity of the gene products, and a relatively large highly variable region, presumably involved in the immunological properties. This is demonstrated by melting profiles of homologous and heterologous RNA hybrid molecules. The results are discussed as a possible mechanism how serologically different influenza strains evolve by mutation in the variable part and selection by the immune system of the host.