A protein with amino acid sequence homology to bovine insulin is present in the legume Vigna unguiculata (cowpea)

Since the discovery of bovine insulin in plants, much effort has been devoted to the characterization of these proteins and elucidation of their functions. We report here the isolation of a protein with similar molecular mass and same amino acid sequence to bovine insulin from developing fruits of cowpea (Vigna unguiculata) genotype Epace 10. Insulin was measured by ELISA using an anti-human insulin antibody and was detected both in empty pods and seed coats but not in the embryo. The highest concentrations (about 0.5 ng/µg of protein) of the protein were detected in seed coats at 16 and 18 days after pollination, and the values were 1.6 to 4.0 times higher than those found for isolated pods tested on any day. N-terminal amino acid sequencing of insulin was performed on the protein purified by C4-HPLC. The significance of the presence of insulin in these plant tissues is not fully understood but we speculate that it may be involved in the transport of carbohydrate to the fruit.

Comparative sequence analysis of the P-, M- and L-coding region of the measles virus CAM-70 live attenuated vaccine strain

Measles virus is a highly contagious agent which causes a major health problem in developing countries. The viral genomic RNA is single-stranded, nonsegmented and of negative polarity. Many live attenuated vaccines for measles virus have been developed using either the prototype Edmonston strain or other locally isolated measles strains. Despite the diverse geographic origins of the vaccine viruses and the different attenuation methods used, there was remarkable sequence similarity of H, F and N genes among all vaccine strains. CAM-70 is a Japanese measles attenuated vaccine strain widely used in Brazilian children and produced by Bio-Manguinhos since 1982. Previous studies have characterized this vaccine biologically and genomically. Nevertheless, only the F, H and N genes have been sequenced. In the present study we have sequenced the remaining P, M and L genes (approximately 1.6, 1.4 and 6.5 kb, respectively) to complete the genomic characterization of CAM-70 and to assess the extent of genetic relationship between CAM-70 and other current vaccines. These genes were amplified using long-range or standard RT-PCR techniques, and the cDNA was cloned and automatically sequenced using the dideoxy chain-termination method. The sequence analysis comparing previously sequenced genotype A strains with the CAM-70 Bio-Manguinhos strain showed a low divergence among them. However, the CAM-70 strains (CAM-70 Bio-Manguinhos and a recently sequenced CAM-70 submaster seed strain) were assigned to a specific group by phylogenetic analysis using the neighbor-joining method. Information about our product at the genomic level is important for monitoring vaccination campaigns and for future studies of measles virus attenuation.

Rarity of DNA sequence alterations in the promoter region of the human androgen receptor gene

The human androgen receptor (AR) gene promoter lies in a GC-rich region containing two principal sites of transcription initiation and a putative Sp1 protein-binding site, without typical "TATA" and "CAAT" boxes. It has been suggested that mutations within the 5'untranslated region (5'UTR) may contribute to the development of prostate cancer by changing the rates of gene transcription and/or translation. In order to investigate this question, the aim of the present study was to search for the presence of mutations or polymorphisms at the AR-5'UTR in 92 prostate cancer patients, where histological diagnosis of adenocarcinoma was established in specimens obtained from transurethral resection or after prostatectomy. The AR-5'UTR was amplified by PCR from genomic DNA samples of the patients and of 100 healthy male blood donors, included as controls. Conformation-sensitive gel electrophoresis was used for DNA sequence alteration screening. Only one band shift was detected in one individual from the blood donor group. Sequencing revealed a new single nucleotide deletion (T) in the most conserved portion of the promoter region at position +36 downstream from the transcription initiation site I. Although the effect of this specific mutation remains unknown, its rarity reveals the high degree of sequence conservation of the human androgen promoter region. Moreover, the absence of detectable variation within the critical 5'UTR in prostate cancer patients indicates a low probability of its involvement in prostate cancer etiology.

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.

APOA1/C3/A4 gene cluster variability and lipid levels in Brazilian children

Genetic studies have suggested that polymorphisms of genes coding for apolipoproteins are significant determinants of serum lipoprotein and lipid levels in adults. However, only a few studies have investigated the association of these polymorphisms in children. Therefore, in the present investigation we studied the distribution of APOA1 -75 G>A, +83 C>T, APOC3 -482 C>T, -455 T>C and 3238 C>G, and APOA4 Q360H and T347S polymorphisms and their influence on plasma lipoprotein levels in children from a Brazilian northeastern admixed population. The seven polymorphic sites were genotyped in 414 children aged 5 to 15 years (mean 8.9 ± 2.9). The genotypes of the seven polymorphic sites were assessed by PCR-RFLP methods. The frequencies of the less common alleles were, in general, intermediate among parental populations, as expected. Strong linkage disequilibrium was detected between polymorphisms at the APOA1, APOC3 and APOA4 loci in this admixed population sample. Overall the genotype effects seen in adults were weaker or absent in children. The APOC3/-455 and APOA4 T347S variants showed significant effects on HDL cholesterol in girls (P = 0.033 and P = 0.016, respectively). Significantly higher plasma total (P = 0.003) and LDL cholesterol (P = 0.004) levels were observed in boys who were carriers of the 3238G allele at the APOC3/3238 C>G site. These results disclosed an overall absence of associations between these polymorphisms and lipids in children. This finding is not unexpected because expression of the effect of these polymorphisms might depend on the interaction with environmental variables both internal and external to the individual.

Heart rate variability in athletes and nonathletes at rest and during head-up tilt

The purpose of the present study was to determine if autonomic heart rate modulation, indicated by heart rate variability (HRV), differs during supine rest and head-up tilt (HUT) when sedentary and endurance-trained cyclists are compared. Eleven sedentary young men (S) and 10 trained cyclists (C) were studied. The volunteers were submitted to a dynamic ECG Holter to calculate HRV at rest and during a 70º HUT. The major aerobic capacity of athletes was expressed by higher values of at anaerobic threshold and peak conditions (P < 0.05). At rest the athletes had lower heart rates (P < 0.05) and higher values in the time domain of HRV compared with controls (SD of normal RR interval, SDNN, medians): 59.1 ms (S) vs 89.9 ms (C), P < 0.05. During tilt athletes also had higher values in the time domain of HRV compared with controls (SDNN, medians): 55.7 ms (S) vs 69.7 ms (C), P < 0.05. No differences in power spectral components of HRV at rest or during HUT were detected between groups. Based on the analysis of data by the frequency domain method, we conclude that in athletes the resting bradycardia seems to be much more related to changes in intrinsic mechanisms than to modifications in autonomic control. Also, HUT caused comparable changes in sympathetic and parasympathetic modulation of the sinus node in both groups.

A semi-automatic computerized method to measure baroreflex-mediated heart rate responses that reduces interobserver variability

Arterial baroreflex sensitivity estimated by pharmacological impulse stimuli depends on intrinsic signal variability and usually a subjective choice of blood pressure (BP) and heart rate (HR) values. We propose a semi-automatic method to estimate cardiovascular reflex sensitivity to bolus infusions of phenylephrine and nitroprusside. Beat-to-beat BP and HR time series for male Wistar rats (N = 13) were obtained from the digitized signal (sample frequency = 2 kHz) and analyzed by the proposed method (PRM) developed in Matlab language. In the PRM, time series were low-pass filtered with zero-phase distortion (3rd order Butterworth used in the forward and reverse direction) and presented graphically, and parameters were selected interactively. Differences between basal mean values and peak BP (deltaBP) and HR (deltaHR) values after drug infusions were used to calculate baroreflex sensitivity indexes, defined as the deltaHR/deltaBP ratio. The PRM was compared to the method traditionally (TDM) employed by seven independent observers using files for reflex bradycardia (N = 43) and tachycardia (N = 61). Agreement was assessed by Bland and Altman plots. Dispersion among users, measured as the standard deviation, was higher for TDM for reflex bradycardia (0.60 ± 0.46 vs 0.21 ± 0.26 bpm/mmHg for PRM, P < 0.001) and tachycardia (0.83 ± 0.62 vs 0.28 ± 0.28 bpm/mmHg for PRM, P < 0.001). The advantage of the present method is related to its objectivity, since the routine automatically calculates the desired parameters according to previous software instructions. This is an objective, robust and easy-to-use tool for cardiovascular reflex studies.

Epidemiological and functional implications of molecular variants of human papillomavirus

Human papillomavirus genomes are classified into molecular variants when they present more than 98% of similarity to the prototype sequence within the L1 gene. Comparative nucleotide sequence analyses of these viruses have elucidated some features of their phylogenetic relationship. In addition, human papillomavirus intratype variability has also been used as an important tool in epidemiological studies of viral transmission, persistence and progression to clinically relevant cervical lesions. Until the present, little has been published concerning the functional significance of molecular variants. It has been shown that nucleotide variability within the long control region leads to differences in the binding affinity of some cellular transcriptional factors and to the enhancement of the expression of E6 and E7 oncogenes. Furthermore, in vivo and in vitro studies revealed differences in E6 and E7 biochemical and biological properties among molecular variants. Nevertheless, further correlation with additional functional information is needed to evaluate the significance of genome intratypic variability. These results are also important for the development of vaccines and to determine the extent to which immunization with L1 virus-like particles of one variant could induce antibodies that cross-neutralize other variants.

Pipeline for macro- and microarray analyses

The pipeline for macro- and microarray analyses (PMmA) is a set of scripts with a web interface developed to analyze DNA array data generated by array image quantification software. PMmA is designed for use with single- or double-color array data and to work as a pipeline in five classes (data format, normalization, data analysis, clustering, and array maps). It can also be used as a plugin in the BioArray Software Environment, an open-source database for array analysis, or used in a local version of the web service. All scripts in PMmA were developed in the PERL programming language and statistical analysis functions were implemented in the R statistical language. Consequently, our package is a platform-independent software. Our algorithms can correctly select almost 90% of the differentially expressed genes, showing a superior performance compared to other methods of analysis. The pipeline software has been applied to 1536 expressed sequence tags macroarray public data of sugarcane exposed to cold for 3 to 48 h. PMmA identified thirty cold-responsive genes previously unidentified in this public dataset. Fourteen genes were up-regulated, two had a variable expression and the other fourteen were down-regulated in the treatments. These new findings certainly were a consequence of using a superior statistical analysis approach, since the original study did not take into account the dependence of data variability on the average signal intensity of each gene. The web interface, supplementary information, and the package source code are available, free, to non-commercial users at http://ipe.cbmeg.unicamp.br/pub/PMmA.

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.

Moderate-term reproducibility of heart rate variability during rest and light to moderate exercise in children

Previous studies have demonstrated the high reproducibility of heart rate variability (HRV) measures in adults while little information exists concerning HRV reproducibility in children. Subsequently, the aim of the current study was to examine the moderate-term reproducibility of heart rate and frequency domain measures of HRV during rest and light to moderate exercise in children. Ten healthy children (6 males, 4 females) aged between 7 and 12 years of age volunteered for this study with HRV recordings obtained during supine rest and three treadmill walking exercise work rates (≤60% maximum heart rate), initially and then 8 weeks later. Differences (P < 0.05) between variables were examined using paired t-tests or Wilcoxon signed rank tests while reliability and reproducibility were examined by intraclass correlation coefficients (ICC), coefficients of variation (CV), and mean bias ratio and ratio limits of agreement (LOA). Heart rate and all measures of HRV at rest and exercise were unchanged after 8 weeks. Significant ICC were documented primarily during rest (0.72-0.85) while weaker relationships (-0.02-0.87) were evident during exercise. A large range of CV was identified during rest (6-33%) and exercise (3-128%) while the ratio LOA were variable and substantial (1.04-2.73). Despite similar HRV over an 8-week period, variable ICC and sizable CV and ratio LOA indicate moderate to poor reproducibility of HRV in children, particularly during light to moderate exercise. Studies examining HRV in children should include age- or maturation stage-matched control participants to address the age-related change in HRV and inadequate HRV reliability.

p.F508del in a heterogeneous cystic fibrosis population from Minas Gerais, Brazil

Cystic fibrosis (CF) is the most common autosomal recessive disease of the Caucasian population. Among the various CF mutations, p.F508del is the most frequent, accounting for two-thirds of the global CF chromosomes, although showing great variability among populations. We have studied 115 unrelated CF patients from a mixed population of Minas Gerais (Brazil). To evaluate part of the DNA sequence of the cystic fibrosis transmembrane conductance regulator (CFTR) gene, blood DNA was obtained and PCR was performed using two pairs of primers that anneal to exons 10 and 24 of the CFTR gene. The PCR product was then submitted to automatic sequencing using the ABI PRISM 310 Genetic Analyzer. The p.F508del mutation was found in 50 (21.7%) of 230 unrelated CF alleles. Fifteen (13.0%) patients were homozygous for this mutation, while 20 (17.4%) were heterozygous; the remaining 80 (69.6%) patients did not carry the p.F508del mutation. Exon 24 sequence had no change in 75 (65.2%) patients, 21 (18.3%) had the sequence variation 4521G/A, 11 (9.6%) had a not yet described sequence variation 4407T/A and 8 (7.0%) patients had both sequence variations (4521G/A and 4407T/A). The polymorphism 4407T/A results in an amino acid modification from aspartic acid to glutamic acid, which will probably have no function effect in CFTR. This low p.F508del prevalence can be due to the variable ethnic origin of this population from Minas Gerais, which may have a high diversity of CF rare mutations.

Comparison of the Polar S810i monitor and the ECG for the analysis of heart rate variability in the time and frequency domains

The aim of the present study was to compare heart rate variability (HRV) at rest and during exercise using a temporal series obtained with the Polar S810i monitor and a signal from a LYNX® signal conditioner (BIO EMG 1000 model) with a channel configured for the acquisition of ECG signals. Fifteen healthy subjects aged 20.9 ± 1.4 years were analyzed. The subjects remained at rest for 20 min and performed exercise for another 20 min with the workload selected to achieve 60% of submaximal heart rate. RR series were obtained for each individual with a Polar S810i instrument and with an ECG analyzed with a biological signal conditioner. The HRV indices (rMSSD, pNN50, LFnu, HFnu, and LF/HF) were calculated after signal processing and analysis. The unpaired Student t-test and intraclass correlation coefficient were used for data analysis. No statistically significant differences were observed when comparing the values analyzed by means of the two devices for HRV at rest and during exercise. The intraclass correlation coefficient demonstrated satisfactory correlation between the values obtained by the devices at rest (pNN50 = 0.994; rMSSD = 0.995; LFnu = 0.978; HFnu = 0.978; LF/HF = 0.982) and during exercise (pNN50 = 0.869; rMSSD = 0.929; LFnu = 0.973; HFnu = 0.973; LF/HF = 0.942). The calculation of HRV values by means of temporal series obtained from the Polar S810i instrument appears to be as reliable as those obtained by processing the ECG signal captured with a signal conditioner.