L1 sequence of a new human papillomavirus type-58 variant associated with cervical intraepithelial neoplasia

The present study on molecular characterization of a human papillomavirus (HPV) isolated in Central Brazil describes the L1 gene sequence from a new variant of HPV-58, the isolate Bsb-02. The sample was from a smear obtained from a woman with cervical intraepithelial neoplasia grade II. The whole L1 gene from isolate Bsb-02 was sequenced automatically, showing 99.1% nucleotide identity with the gene from the HPV-58 reference. The clustering between Bsb-02 and HPV-58 reference sequence was also supported by phylogenetic analysis. Fourteen nucleotide substitutions were observed: eight were synonymous and six were associated with amino acid substitutions. A10V and V144I have not been previously described. At GenBank, the only complete L1 sequence from HPV-58 in addition to the HPV-58 reference one is that of Bsb-02. These data provide information that may be relevant to HPV diagnosis and to rational vaccine strategies. HPV variants may also be associated with host immune responses and with the risk of cervical neoplasia.

Genotyping hepatitis C virus from hemodialysis patients in Central Brazil by line probe assay and sequence analysis

The present study examined the distribution of hepatitis C virus (HCV) genotypes and subtypes in a hemodialysis population in Goiás State, Central Brazil, and evaluated the efficiency of two genotyping methods: line probe assay (LiPA) based on the 5' noncoding region and nucleotide sequencing of the nonstructural 5B (NS5B) region of the genome. A total of 1095 sera were tested for HCV RNA by RT-nested PCR of the 5' noncoding region. The LiPA assay was able to genotype all 131 HCV RNA-positive samples. Genotypes 1 (92.4%) and 3 (7.6%) were found. Subtype 1a (65.7%) was the most prevalent, followed by subtypes 1b (26.7%) and 3a (7.6%). Direct nucleotide sequencing of 340 bp from the NS5B region was performed in 106 samples. The phylogenetic tree showed that 98 sequences (92.4%) were classified as genotype 1, subtypes 1a (72.6%) and 1b (19.8%), and 8 sequences (7.6%) as subtype 3a. The two genotyping methods gave concordant results within HCV genotypes and subtypes in 100 and 96.2% of cases, respectively. Only four samples presented discrepant results, with LiPA not distinguishing subtypes 1a and 1b. Therefore, HCV genotype 1 (subtype 1a) is predominant in hemodialysis patients in Central Brazil. By using sequence analysis of the NS5B region as a reference standard method for HCV genotyping, we found that LiPA was efficient at the genotype level, although some discrepant results were observed at the subtype level (sensitivity of 96.1% for subtype 1a and 95.2% for subtype 1b). Thus, analysis of the NS5B region permitted better discrimination between HCV subtypes, as required in epidemiological investigations.

Sequence change in the HS2-LCR and Gg-globin gene promoter region of sickle cell anemia patients

The fetal hemoglobin (HbF) levels and ßS-globin gene haplotypes of 125 sickle cell anemia patients from Brazil were investigated. We sequenced the Gg- and Ag-globin gene promoters and the DNase I-2 hypersensitive sites in the locus control regions (HS2-LCR) of patients with HbF level disparities as compared to their ßS haplotypes. Sixty-four (51.2%) patients had CAR/Ben genotype; 36 (28.8%) Ben/Ben; 18 (14.4%) CAR/CAR; 2 (1.6%) CAR/Atypical; 2 (1.6%) Ben/Cam; 1 (0.8%) CAR/Cam; 1 (0.8%) CAR/Arab-Indian, and 1 (0.8%) Sen/Atypical. The HS2-LCR sequence analyses demonstrated a c.-10.677G>A change in patients with the Ben haplotype and high HbF levels. The Gg gene promoter sequence analyses showed a c.-157T>C substitution shared by all patients, and a c.-222_-225del related to the Cam haplotype. These results identify new polymorphisms in the HS2-LCR and Gg-globin gene promoter. Further studies are required to determine the correlation between HbF synthesis and the clinical profile of sickle cell anemia patients.

Cloning, sequence analysis, and expression of cDNA coding for the major house dust mite allergen, Der f 1, in Escherichia coli

Our objective was to clone, express and characterize adult Dermatophagoides farinae group 1 (Der f 1) allergens to further produce recombinant allergens for future clinical applications in order to eliminate side reactions from crude extracts of mites. Based on GenBank data, we designed primers and amplified the cDNA fragment coding for Der f 1 by nested-PCR. After purification and recovery, the cDNA fragment was cloned into the pMD19-T vector. The fragment was then sequenced, subcloned into the plasmid pET28a(+), expressed in Escherichia coli BL21 and identified by Western blotting. The cDNA coding for Der f 1 was cloned, sequenced and expressed successfully. Sequence analysis showed the presence of an open reading frame containing 966 bp that encodes a protein of 321 amino acids. Interestingly, homology analysis showed that the Der p 1 shared more than 87% identity in amino acid sequence with Eur m 1 but only 80% with Der f 1. Furthermore, phylogenetic analyses suggested that D. pteronyssinus was evolutionarily closer to Euroglyphus maynei than to D. farinae, even though D. pteronyssinus and D. farinae belong to the same Dermatophagoides genus. A total of three cysteine peptidase active sites were found in the predicted amino acid sequence, including 127-138 (QGGCGSCWAFSG), 267-277 (NYHAVNIVGYG) and 284-303 (YWIVRNSWDTTWGDSGYGYF). Moreover, secondary structure analysis revealed that Der f 1 contained an a helix (33.96%), an extended strand (17.13%), a ß turn (5.61%), and a random coil (43.30%). A simple three-dimensional model of this protein was constructed using a Swiss-model server. The cDNA coding for Der f 1 was cloned, sequenced and expressed successfully. Alignment and phylogenetic analysis suggests that D. pteronyssinus is evolutionarily more similar to E. maynei than to D. farinae.

Promoter region sequence differences in the A and G gamma globin genes of Brazilian sickle cell anemia patients

Fetal hemoglobin (HbF), encoded by the HBG2 and HBG1 genes, is the best-known genetic modulator of sickle cell anemia, varying dramatically in concentration in the blood of these patients. This variation is partially associated with polymorphisms located in the promoter region of the HBG2 and HBG1 genes. In order to explore known and unknown polymorphisms in these genes, the sequences of their promoter regions were screened in sickle cell anemia patients and correlated with both their HbF levels and their βS-globin haplotypes. Additionally, the sequences were compared with genes from 2 healthy groups, a reference one (N = 104) and an Afro-descendant one (N = 98), to identify polymorphisms linked to the ethnic background.The reference group was composed by healthy individuals from the general population. Four polymorphisms were identified in the promoter region of HBG2 and 8 in the promoter region of HBG1 among the studied groups. Four novel single nucleotide polymorphisms (SNP) located at positions -324, -317, -309 and -307 were identified in the reference group. A deletion located between -396 and -391 in the HBG2 promoter region and the SNP -271 C→T in the HBG1 promoter region were associated with the Central African Republic βS-globin haplotype. In contrast, the -369 C→G and 309 A→G SNPs in the HBG2 promoter region were correlated to the Benin haplotype. The polymorphisms -396_-391 del HBG2, -369 SNP HBG2 and -271 SNP HBG1 correlated with HbF levels. Hence, we suggest an important role of HBG2 and HBG1 gene polymorphisms on the HbF synthesis.

Gene expression promoted by the SV40 DNA targeting sequence and the hypoxia-responsive element under normoxia and hypoxia

The main objective of the present study was to find suitable DNA-targeting sequences (DTS) for the construction of plasmid vectors to be used to treat ischemic diseases. The well-known Simian virus 40 nuclear DTS (SV40-DTS) and hypoxia-responsive element (HRE) sequences were used to construct plasmid vectors to express the human vascular endothelial growth factor gene (hVEGF). The rate of plasmid nuclear transport and consequent gene expression under normoxia (20% O2) and hypoxia (less than 5% O2) were determined. Plasmids containing the SV40-DTS or HRE sequences were constructed and used to transfect the A293T cell line (a human embryonic kidney cell line) in vitro and mouse skeletal muscle cells in vivo. Plasmid transport to the nucleus was monitored by real-time PCR, and the expression level of the hVEGF gene was measured by ELISA. The in vitro nuclear transport efficiency of the SV40-DTS plasmid was about 50% lower under hypoxia, while the HRE plasmid was about 50% higher under hypoxia. Quantitation of reporter gene expression in vitro and in vivo, under hypoxia and normoxia, confirmed that the SV40-DTS plasmid functioned better under normoxia, while the HRE plasmid was superior under hypoxia. These results indicate that the efficiency of gene expression by plasmids containing DNA binding sequences is affected by the concentration of oxygen in the medium.

New multilocus sequence typing of MRSA in São Paulo, Brazil

An increased incidence of nosocomial and community-acquired infections caused by methicillin-resistant Staphylococcus aureus (MRSA) has been observed worldwide. The molecular characterization of MRSA has played an important role in demonstrating the existence of internationally disseminated clones. The use of molecular biology methods in the surveillance programs has enabled the tracking of MRSA spread within and among hospitals. These data are useful to alert nosocomial infection control programs about the potential introduction of these epidemic clones in their areas. Four MRSA blood culture isolates from patients hospitalized at two hospitals in the city of São Paulo, Brazil, were analyzed; one of them was community acquired. The isolates were characterized as SCCmec, mecA and PVL by PCR, pulsed-field gel electrophoresis (PFGE) profile and molecular sequence typing (MLST) genotyping. The isolates presented type IV SCCmec, and none proved to be positive for PVL. The isolates showed a PFGE profile similar to the pediatric clone. MLST genotyping demonstrated that the isolates belonged to clonal complex 5 (CC5), showing a new yqiL allele gene, resulting in a new sequence typing (ST) (1176). Our results showed that strains of MRSA carrying a new ST are emerging in community and nosocomial infections, including bacteremia, in São Paulo, Brazil.

Sequence analysis of the 5′ third of glycoprotein C gene of South American bovine herpesviruses 1 and 5

Bovine herpesviruses 1 (BoHV-1) and 5 (BoHV-5) share high genetic and antigenic similarities, but exhibit marked differences in tissue tropism and neurovirulence. The amino-terminal region of glycoprotein C (gC), which is markedly different in each of the viruses, is involved in virus binding to cellular receptors and in interactions with the immune system. This study investigated the genetic and antigenic differences of the 5′ region of the gC (5′ gC) gene (amino-terminal) of South American BoHV-1 (n=19) and BoHV-5 (n=25) isolates. Sequence alignments of 374 nucleotides (104 amino acids) revealed mean similarity levels of 97.3 and 94.2% among BoHV-1 gC (gC1), respectively, 96.8 and 95.6% among BoHV-5 gC (gC5), and 62 and 53.3% between gC1 and gC5. Differences included the absence of 40 amino acid residues (27 encompassing predicted linear epitopes) scattered throughout 5′ gC1 compared to 5′ gC5. Virus neutralizing assays testing BoHV-1 and BoHV-5 antisera against each isolate revealed a high degree of cross-neutralization between the viruses, yet some isolates were neutralized at very low titers by heterologous sera, and a few BoHV-5 isolates reacted weakly with either sera. The virus neutralization differences observed within the same viral species, and more pronounced between BoHV-1 and BoHV-5, likely reflect sequence differences in neutralizing epitopes. These results demonstrate that the 5′ gC region is well conserved within each viral species but is divergent between BoHV-1 and BoHV-5, likely contributing to their biological and antigenic differences.

The consensus sequence of FAMLF alternative splice variants is overexpressed in undifferentiated hematopoietic cells

The familial acute myeloid leukemia related factor gene (FAMLF) was previously identified from a familial AML subtractive cDNA library and shown to undergo alternative splicing. This study used real-time quantitative PCR to investigate the expression of the FAMLF alternative-splicing transcript consensus sequence (FAMLF-CS) in peripheral blood mononuclear cells (PBMCs) from 119 patients with de novo acute leukemia (AL) and 104 healthy controls, as well as in CD34+cells from 12 AL patients and 10 healthy donors. A 429-bp fragment from a novel splicing variant of FAMLF was obtained, and a 363-bp consensus sequence was targeted to quantify total FAMLF expression. Kruskal-Wallis, Nemenyi, Spearman's correlation, and Mann-Whitney U-tests were used to analyze the data. FAMLF-CS expression in PBMCs from AL patients and CD34+ cells from AL patients and controls was significantly higher than in control PBMCs (P<0.0001). Moreover,FAMLF-CS expression in PBMCs from the AML group was positively correlated with red blood cell count (rs=0.317, P=0.006), hemoglobin levels (rs=0.210, P=0.049), and percentage of peripheral blood blasts (rs=0.256, P=0.027), but inversely correlated with hemoglobin levels in the control group (rs=–0.391, P<0.0001). AML patients with high CD34+ expression showed significantly higherFAMLF-CS expression than those with low CD34+ expression (P=0.041). Our results showed thatFAMLF is highly expressed in both normal and malignant immature hematopoietic cells, but that expression is lower in normal mature PBMCs.

Sequence-characterized amplified region (SCAR) technique used for variety discrimination in vinca (Catharanthus roseus (L.) G.Don)

Sequence-Characterized Amplified Region (SCAR) appears as a useful technique for genetic purity testing and variety discrimination, applicable to species in which some other techniques have failed. In particular, this technique is very attractive with species in which RAPD results were not consistent. The RAPD polymorphic bands were cloned, sequenced and from the sequence information, primers pairs for normal PCR were developed. Since the probability of obtaining successful SCAR primers from RAPD polymorphic bands was about 50%, a larger number of RAPD polymorphic bands are needed to develop sufficient SCAR primers for varietal discrimination in vinca. In addition, the efficiency of the SCAR technique is strongly affected by the quality of DNA extracted from seeds. The SCAR banding patterns obtained from vinca seed were consistent and repeatable making the results reliable for genetic purity testing and variety discrimination. The SCAR technique is simple, fast, relatively inexpensive and allows the use of DNA extracted from dry seeds, which is very important in a seed-quality evaluating program

Probing MST1 Kinase Sequence and Structure for Rational Drug Discovery (Targeting diabetes by blocking protein-protein interactions)

Apoptotic beta cell death is an underlying cause majorly for type I and to a lesser extent for type II diabetes. Recently, MST1 kinase was identified as a key apoptotic agent in diabetic condition. In this study, I have examined MST1 and closely related kinases namely, MST2, MST3 and MST4, aiming to tackle diabetes by exploring ways to selectively block MST1 kinase activity. The first investigation was directed towards evaluating possibilities of selectively blocking the ATP binding site of MST1 kinase that is essential for the activity of the enzymes. Structure and sequence analyses of this site however revealed a near absolute conservation between the MSTs and very few changes with other kinases. The observed residue variations also displayed similar physicochemical properties making it hard for selective inhibition of the enzyme. Second, possibilities for allosteric inhibition of the enzyme were evaluated. Analysis of the recognized allosteric site also posed the same problem as the MSTs shared almost all of the same residues. The third analysis was made on the SARAH domain, which is required for the dimerization and activation of MST1 and MST2 kinases. MST3 and MST4 lack this domain, hence selectivity against these two kinases can be achieved. Other proteins with SARAH domains such as the RASSF proteins were also examined. Their interaction with the MST1 SARAH domain were evaluated to mimic their binding pattern and design a peptide inhibitor that interferes with MST1 SARAH dimerization. In molecular simulations the RASSF5 SARAH domain was shown to strongly interact with the MST1 SARAH domain and possibly preventing MST1 SARAH dimerization. Based on this, the peptidic inhibitor was suggested to be based on the sequence of RASSF5 SARAH domain. Since the MST2 kinase also interacts with RASSF5 SARAH domain, absolute selectivity might not be achieved.

Continental margin architecture : the palynological signature of glacioeustasy /

Palynomorphs from two siliciclastic margins were examined to gain insights into continental margin architecture. Sea level change is thought to be one of the primary controls on continental margin architecture. Because Late Neogene glacioeustasy has been well studied marine sediments deposited during the Late Neogene were examined to test this concept. Cores from the outer shelf and upper slope were taken from the New Jersey margin in the western North Atlantic Ocean and from the Sunda Shelf margin in the South China Sea. Continental margin architecture is often described in a sequence stratigraphic context. One of the main goals of both coring projects was to test the theoretical sequence stratigraphic models developed by a research group at Exxon (e.g. Wilgus et al., 1988). Palynomorphs provide one of the few methods of inferring continental margin architecture in monotonous, siliciclastic marine sediments where calcareous sediments are rare (e.g. New Jersey margin). In this study theoretical models of the palynological signature expected in sediment packages deposited during the various increments of a glacioeustatic cycle were designed. These models were based on the modem palynomorph trends and taphonomic factors thought to control palynomorph distribution. Both terrestrial (pollen and spores) and marine (dinocysts) palynomorphs were examined. The palynological model was then compared with New Jersey margin and Sunda Shelf margin sediments. The predicted palynological trends provided a means of identifying a complete cycle of glacioeustatic change (Oxygen Isotope Stage 5e to present) in the uppermost 80 meters of sediment on the slope at the New Jersey margin. Sediment availability, not sea meters of sediment on the slope at the New Jersey margin. Sediment availability, not sea level change, is thought to be the major factor controlling margin architecture during the late Pleistocene here at the upper slope. This is likely a function of the glacial scouring of the continents which significantly increases sediment availability during glacial stages. The subaerially exposed continental shelf during the lowstand periods would have been subject to significant amounts of erosion fi:om the proglacial rivers flowing fi-om the southern regions of the ice-sheet. The slope site is non-depositional today and was also non-depositional during the last full interglacial period. The palynomorph data obtained fi-om the South China Sea indicate that the major difference between the New Jersey Margin sites and the Sunda Shelf margin sites is the variation in sediment supply and the rate of sediment accumulation. There was significantly less variation in sediment supply between glacial and interglacial periods and less overall sediment accumulation at the Sunda Shelf margin. The data presented here indicate that under certain conditions the theoretical palynological models allow the identification of individual sequence stratigraphic units and therefore, allow inferences regarding continental margin architecture. The major condition required in this approach is that a complete and reliable database of the contemporaneous palynomorphs be available.

Microfacies analysis of the lower paleocene beda formation in Western Concession 59 and Block 59F, Sirte Basin, Libya /

The rock sequence of the Tertiary Beda Formation of S. W. concession 59 and 59F block in Sirte Basin of Libya has been subdivided into twelve platformal carbonate microfacies. These microfacies are dominated by muddy carbonates, such as skeletal mudstones, wackestones, and packstones with dolomites and anhydrite. Rock textures, faunal assemblages and sedimentary structures suggest shallow, clear, warm waters and low to moderate energy conditions within the depositional shelf environment. The Beda Formation represents a shallowing-upward sequence typical of lagoonal and tidal flat environments marked at the top by sabkha and brackish-water sediments. Microfossils include benthonic foraminifera, such as miliolids, Nummulites, - oerculina and other smaller benthonics, in addition to dasycladacean algae, ostracods, molluscs, echinoderms, bryozoans and charophytes. Fecal pellets and pelloids, along with the biotic allochems, contributed greatly to the composition of the various microfacies. Dolomite, where present, is finely crystalline and an early replacement product. Anhydrite occurs as nodular, chickenwire and massive textures indicating supratidal sabkha deposition. Compaction, micr it i zat ion , dolomit izat ion , recrystallization, cementation, and dissolution resulted in alteration and obliteration of primary sedimentary structures of the Beda Formation microfacies. The study area is located in the Gerad Trough which developed as a NE-SW trending extensional graben. The Gerad trough was characterized by deep-shallow water conditions throughout the deposition of the Beda Formation sediments. The study area is marked by several horsts and grabens; as a result of extent ional tectonism. The area was tectonically active throughout the Tertiary period. Primary porosity is intergranular and intragranular, and secondary processes are characterized by dissolution, intercrystalline, fracture and fenestral features. Diagenesis, through solution leaching and dolomitization, contributed greatly to porosity development. Reservoir traps of the Beda Formation are characterized by normal fault blocks and the general reservoir characteristics/properties appear to be facies controlled.

Geochemistry and petrogenesis of the McElroy and Larder Lake assemblages, Abitibi Greenstone Belt, Northeastern Ontario

The McElroy and Larder Lake assemblages, located in the southern Abitibi Greenstone Belt are two late Archean metavolcanic sequences having markedly contrasting physical characteristics arid are separated from one another by a regional fault. An assemblage is an informal term which describes stratified volcanic and/or sedimentary rock units built during a specific time period in a similar depositional or volcanic setting and are commonly bounded by faults, unconformities or intrusions. The petrology and petrogenesis of these assemblages have been investigated to determine if a genetic link exists between the two adjacent assemblages. The McElroy assemblage is homoclinal sequence of evolved massive and pillowed fl.ows, which except for the basal unit represents a progressively fractionated volcanic pile. From the base to the top of the assemblage the lithologies include Fe-tholeiitic, dendritic flows; komatiite basaltic, ultramafic flows; Mg-tholeiitic, leucogabbro; Mg-tholeiitic, massive flows and Fe-tholeiitic, pillowed flows. Massive flows range from coarse grained to aphanitic and are commonly plagioclase glomerophyric. The Larder Lake assemblage consists of komatiitic, Mg-rich and Fe-rich tholeiitic basalts, structurally disrupted by folds and faults. Tholeiitic rocks in the Larder Lake assemblage range from aphanitic to coarse grained massive and pillowed flows. Komatiitic flows contain both spinifex and massive textures. Geochemical variability within both assemblages is attributed to different petrogenetic histories. The lithologies of the McElroy assemblage were derived by partial melting of a primitive mantle source followed by various degrees of crystal fractionation. Partial melting of a primitive mantle source generated the ultramafic flows and possibly other flows in the assemblage. Fractionation of ultramafic flows may have also produced the more evolved McElroy lithologies. The highly evolved, basal, dendritic flow may represent the upper unit 3 of a missing volcanic pile in which continued magmatism generated the remaining McElroy lithologies. Alternatively, the dendritic flows may represent a primary lava derived from a low degree (10-15%) partial melt of a primitive mantle source which was followed by continued partial melting to generate the ultramafic flows. The Larder Lake lithologies were derived by partial melting of a komatiitic source followed by gabbroic fractionation. The tectonic environment for both assemblages is interpreted to be an oceanic arc setting. The McElroy assemblage lavas were generated in a mature back arc setting whereas the Larder Lake lithologies were produced during the early stages of komatiitc crust subduction. This setting is consistent with previous models involving plate tectonic processes for the generation of other metavolcanic assemblages in the Abitibi Greenstone Belt.

Improving Short DNA Sequence Alignment with Parallel Computing

Variations in different types of genomes have been found to be responsible for a large degree of physical diversity such as appearance and susceptibility to disease. Identification of genomic variations is difficult and can be facilitated through computational analysis of DNA sequences. Newly available technologies are able to sequence billions of DNA base pairs relatively quickly. These sequences can be used to identify variations within their specific genome but must be mapped to a reference sequence first. In order to align these sequences to a reference sequence, we require mapping algorithms that make use of approximate string matching and string indexing methods. To date, few mapping algorithms have been tailored to handle the massive amounts of output generated by newly available sequencing technologies. In otrder to handle this large amount of data, we modified the popular mapping software BWA to run in parallel using OpenMPI. Parallel BWA matches the efficiency of multithreaded BWA functions while providing efficient parallelism for BWA functions that do not currently support multithreading. Parallel BWA shows significant wall time speedup in comparison to multithreaded BWA on high-performance computing clusters, and will thus facilitate the analysis of genome sequencing data.