Sequence analysis, chromosomal distribution and long-range organization show that rapid turnover of new and old pBuM satellite DNA repeats leads to different patterns of variation in seven species of the Drosophila buzzatii cluster

We aimed to study patterns of variation and factors influencing the evolutionary dynamics of a satellite DNA, pBuM, in all seven Drosophila species from the buzzatii cluster (repleta group). We analyzed 117 alpha pBuM-1 (monomer length 190 bp) and 119 composite alpha/beta (370 bp) pBuM-2 repeats and determined the chromosome location and long-range organization on DNA fibers of major sequence variants. Such combined methodologies in the study of satDNAs have been used in very few organisms. In most species, concerted evolution is linked to high copy number of pBuM repeats. Species presenting low-abundance and scattered distributed pBuM repeats did not undergo concerted evolution and maintained part of the ancestral inter-repeat variability. The alpha and alpha/beta repeats colocalized in heterochromatic regions and were distributed on multiple chromosomes, with notable differences between species. High-resolution FISH revealed array sizes of a few kilobases to over 0.7 Mb and mutual arrangements of alpha and alpha/beta repeats along the same DNA fibers, but with considerable changes in the amount of each variant across species. From sequence, chromosomal and phylogenetic data, we could infer that homogenization and amplification events involved both new and ancestral pBuM variants. Altogether, the data on the structure and organization of the pBuM satDNA give insights into genome evolution including mechanisms that contribute to concerted evolution and diversification.

Long Terminal Repeat Sequence Analysis of HTLV-2 Molecular Variants Identified in Southern Brazil

In Brazil, human T-lymphotropic virus type 2 (HTLV-2) is endemic in Amerindians and epidemic in intravenous drug users (IDUs). The long terminal repeat (LTR) is the most divergent genomic region of HTLV-2, therefore useful to characterize subtypes. Nucleotide sequence and restriction fragment length polymorphism (RFLP) analysis of LTR genomic segments of fourteen HTLV-2 strains isolated from HIV-infected patients of Londrina, Southern Brazil, were carried out. Molecular analysis disclosed that all HTLV-2 strains belonged to 2a subtype, and RFLP detected the presence of the a4, a5, and a6 subgroups according to Switzer's nomenclature. RFLP correlated with nucleotide sequence, and phylogenetic analysis clustered HTLV-2 sequences of IDUs into subgroups a5 and a6. HTLV-2 sequences from individuals of sexual risk factor clustered into the a4 subgroup. These results extend the knowledge of the genetic diversity of HTLV-2 circulating in Brazil and provide insights into HTLV-2 transmission and virus movement in this geographic area.

Genome Survey Sequence Analysis and Identification of Homologs of Major Surface Protease (gp63) Genes in Trypanosoma rangeli

In this study, 222 genome survey sequences were generated for Trypanosoma rangeli strain P07 isolated from an opossum (Didelphis albiventris) in Minas Gerais State, Brazil. T. rangeli sequences were compared by BLASTX (Basic Local Alignment Search Tool X) analysis with the assembled contigs of Leishmania braziliensis, Leishmania infantum, Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Results revealed that 82% (182/222) of the sequences were associated with predicted proteins described, whereas 18% (40/222) of the sequences did not show significant identity with sequences deposited in databases, suggesting that they may represent T. rangeli-specific sequences. Among the 182 predicted sequences, 179 (80.6%) had the highest similarity with T. cruzi, 2 (0.9%) with T. brucei, and 1 (0.5%) with L. braziliensis. Computer analysis permitted the identification of members of various gene families described for trypanosomatids in the genome of T. rangeli, such as trans-sialidases, mucin-associated surface proteins, and major surface proteases (MSP or gp63). This is the first report identifying sequences of the MSP family in T. rangeli. Multiple sequence alignments showed that the predicted MSP of T. rangeli presented the typical characteristics of metalloproteases, such as the presence of the HEXXH motif, which corresponds to a region previously associated with the catalytic site of the enzyme, and various cysteine and proline residues, which are conserved among MSPs of different trypanosomatid species. Reverse transcriptase-polymerase chain reaction analysis revealed the presence of MSP transcripts in epimastigote forms of T. rangeli.

Characterization by restriction fragment length polymorphism and sequence analysis of field and vaccine strains of infectious laryngotracheitis virus involved in severe outbreaks

At the end of 2002 and throughout 2003, there was a severe outbreak of infectious laryngotracheitis (ILT) in an intensive production area of commercial hens in the Sao Paulo State of Brazil. ILT virus was isolated from 28 flocks, and 21 isolates were genotyped by polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) using four genes and eight restriction enzymes, and by partial sequencing of the infected cell protein 4 (ICP4) and thymidine kinase (TK) genes. Three groups resulted from the combinations of PCR-RFLP patterns: 19 field isolates formed Group I, and the remaining two isolates together with the chicken embryo origin (CEO) vaccine strains formed Group II. Group III comprised the tissue-culture origin (TCO) vaccine strain by itself. The PCR-RFLP results agreed with the sequencing results of two ICP4 gene fragments. The ICP4 gene sequence analysis showed that the 19 field isolates classified into Group I by RFLP-PCR were identical among themselves, but were different to the TCO and CEO vaccines. The two Group II isolates could not be distinguished from one of the CEO vaccines. The nucleotide and amino acid sequence analyses discriminated between the Brazilian and non-Brazilian isolates, as well as between the TCO and CEO vaccines. Sequence analysis of the TK gene enabled classification of the field isolates (Group I) as virulent and non-vaccine. This work shows that the severe ILT outbreak was caused by a highly virulent, non-vaccine strain.

Analysis of the genome of Spodoptera frugiperda nucleopolyhedrovirus (SfMNPV-19) and of the high genomic heterogeneity in group II nucleopolyhedroviruses

The genome of the most virulent among 22 Brazilian geographical isolates of Spodoptera frugiperda nucleopolyhedrovirus, isolate 19 (SfMNPV-1 9), was completely sequenced and shown to comprise 132 565 bp and 141 open reading frames (ORFs). A total of 11 ORFs with no homology to genes in the GenBank database were found. Of those, four had typical baculovirus; promoter motifs and polyadenylation sites. Computer-simulated restriction enzyme cleavage patterns of SfMNPV-1 9 were compared with published physical maps of other SfMNPV isolates. Differences were observed in terms of the restriction profiles and genome size. Comparison of SfMNPV-1 9 with the sequence of the SfMNPV isolate 3AP2 indicated that they differed due to a 1427 bp deletion, as well as by a series of smaller deletions and point mutations. The majority of genes of SfMNPV-1 9 were conserved in the closely related Spodoptera exigua NPV (SeMNPV) and Agrotis segetum NPV (AgseMNPV-A), but a few regions experienced major changes and rearrangements. Synthenic maps for the genomes of group 11 NPVs revealed that gene collinearity was observed only within certain clusters. Analysis of the dynamics of gene gain and loss along the phylogenetic tree of the NPVs showed that group 11 had only five defining genes and supported the hypothesis that these viruses form ten highly divergent ancient lineages. Crucially, more than 60% of the gene gain events followed a power-law relation to genetic distance among baculoviruses, indicative of temporal organization in the gene accretion process.

High genetic diversity in field isolates of Trypanosoma theileri assessed by analysis of cathepsin L-like sequences disclosed multiple and new genotypes infecting cattle in Thailand

In this study, we describe the first survey in Thailand of Trypanosoma theileri, a widespread and prevalent parasite of cattle that is transmitted by tabanid flies. Investigation of 210 bovine blood samples of Thai cattle from six farms by hematocrit centrifuge technique (HCT) revealed 14 samples with trypanosomes morphologically compatible to T. theileri. Additional animals were positive for T. theileri by PCR based on the Cathepsin L-like sequence (TthCATL-PCR) despite negative by HCT, indicating cryptic infections. Results revealed a prevalence of 26 +/- 15% (95% CI) of T. theileri infection. Additionally, 12 samples positive for T. theileri were detected in cattle from other 11 farms. From a total of 30 blood samples positive by HCT and/or PCR from 17 farms, seven were characterized to evaluate the genetic polymorphism of T. theileri through sequence analysis of PCR-amplified CATL DNA sequences. All CATL sequences of T. theileri from Thai cattle clustered with sequences of the previously described phylogenetic lineages TthI and TthII, supporting only two major lineages of T. theileri in cattle around the world. However, 11 of the 29 CATL sequences analyzed showed to be different, disclosing an unexpectedly large polymorphic genetic repertoire, with multiple genotypes of T. theileri not previously described in other countries circulating in Thai cattle. (C) 2011 Elsevier B.V. All rights reserved.

FAM5C Contributes to Aggressive Periodontitis

Aggressive periodontitis is characterized by a rapid and severe periodontal destruction in young systemically healthy subjects. A greater prevalence is reported in Africans and African descendent groups than in Caucasians and Hispanics. We first fine mapped the interval 1q24.2 to 1q31.3 suggested as containing an aggressive periodontitis locus. Three hundred and eighty-nine subjects from 55 pedigrees were studied. Saliva samples were collected from all subjects, and DNA was extracted. Twenty-one single nucleotide polymorphisms were selected and analyzed by standard polymerase chain reaction using TaqMan chemistry. Non-parametric linkage and transmission distortion analyses were performed. Although linkage results were negative, statistically significant association between two markers, rs1935881 and rs1342913, in the FAM5C gene and aggressive periodontitis (p = 0.03) was found. Haplotype analysis showed an association between aggressive periodontitis and the haplotype A-G (rs1935881-rs1342913; p = 0.009). Sequence analysis of FAM5C coding regions did not disclose any mutations, but two variants in conserved intronic regions of FAM5C, rs57694932 and rs10494634, were found. However, these two variants are not associated with aggressive periodontitis. Secondly, we investigated the pattern of FAM5C expression in aggressive periodontitis lesions and its possible correlations with inflammatory/immunological factors and pathogens commonly associated with periodontal diseases. FAM5C mRNA expression was significantly higher in diseased versus healthy sites, and was found to be correlated to the IL-1 beta, IL-17A, IL-4 and RANKL mRNA levels. No correlations were found between FAM5C levels and the presence and load of red complex periodontopathogens or Aggregatibacter actinomycetemcomitans. This study provides evidence that FAM5C contributes to aggressive periodontitis.

Multidimensional cluster stability analysis from a Brazilian Bradyrhizobium sp RFLP/PCR data set

The taxonomy of the N(2)-fixing bacteria belonging to the genus Bradyrhizobium is still poorly refined, mainly due to conflicting results obtained by the analysis of the phenotypic and genotypic properties. This paper presents an application of a method aiming at the identification of possible new clusters within a Brazilian collection of 119 Bradryrhizobium strains showing phenotypic characteristics of B. japonicum and B. elkanii. The stability was studied as a function of the number of restriction enzymes used in the RFLP-PCR analysis of three ribosomal regions with three restriction enzymes per region. The method proposed here uses Clustering algorithms with distances calculated by average-linkage clustering. Introducing perturbations using sub-sampling techniques makes the stability analysis. The method showed efficacy in the grouping of the species B. japonicum and B. elkanii. Furthermore, two new clusters were clearly defined, indicating possible new species, and sub-clusters within each detected cluster. (C) 2008 Elsevier B.V. All rights reserved.

Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production

Bioethanol is a biofuel produced mainly from the fermentation of carbohydrates derived from agricultural feedstocks by the yeast Saccharomyces cerevisiae. One of the most widely adopted strains is PE-2, a heterothallic diploid naturally adapted to the sugar cane fermentation process used in Brazil. Here we report the molecular genetic analysis of a PE-2 derived diploid (JAY270), and the complete genome sequence of a haploid derivative (JAY291). The JAY270 genome is highly heterozygous (similar to 2 SNPs/kb) and has several structural polymorphisms between homologous chromosomes. These chromosomal rearrangements are confined to the peripheral regions of the chromosomes, with breakpoints within repetitive DNA sequences. Despite its complex karyotype, this diploid, when sporulated, had a high frequency of viable spores. Hybrid diploids formed by outcrossing with the laboratory strain S288c also displayed good spore viability. Thus, the rearrangements that exist near the ends of chromosomes do not impair meiosis, as they do not span regions that contain essential genes. This observation is consistent with a model in which the peripheral regions of chromosomes represent plastic domains of the genome that are free to recombine ectopically and experiment with alternative structures. We also explored features of the JAY270 and JAY291 genomes that help explain their high adaptation to industrial environments, exhibiting desirable phenotypes such as high ethanol and cell mass production and high temperature and oxidative stress tolerance. The genomic manipulation of such strains could enable the creation of a new generation of industrial organisms, ideally suited for use as delivery vehicles for future bioenergy technologies.

Identification of protein coding regions using the modified Gabor-wavelet transform

An important topic in genomic sequence analysis is the identification of protein coding regions. In this context, several coding DNA model-independent methods based on the occurrence of specific patterns of nucleotides at coding regions have been proposed. Nonetheless, these methods have not been completely suitable due to their dependence on an empirically predefined window length required for a local analysis of a DNA region. We introduce a method based on a modified Gabor-wavelet transform (MGWT) for the identification of protein coding regions. This novel transform is tuned to analyze periodic signal components and presents the advantage of being independent of the window length. We compared the performance of the MGWT with other methods by using eukaryote data sets. The results show that MGWT outperforms all assessed model-independent methods with respect to identification accuracy. These results indicate that the source of at least part of the identification errors produced by the previous methods is the fixed working scale. The new method not only avoids this source of errors but also makes a tool available for detailed exploration of the nucleotide occurrence.

Expression analysis of putative vitellogenin and lipophorin receptors in honey bee (Apis mellifera L.) queens and workers

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth. (C) 2008 Elsevier Ltd. All rights reserved.

The p.A2215D Thyroglobulin Gene Mutation Leads to Deficient Synthesis and Secretion of the Mutated Protein and Congenital Hypothyroidism with Wide Phenotype Variation

Context: Thyroglobulin (TG) is a large glycoprotein and functions as a matrix for thyroid hormone synthesis. TG gene mutations give rise to goitrous congenital hypothyroidism (CH) with considerable phenotype variation. Objectives: The aim of the study was to report the genetic screening of 15 patients with CH due to TG gene mutations and to perform functional analysis of the p. A2215D mutation. Design: Clinical evaluation and DNA sequencing of the TG gene were performed in all patients. TG expression was analyzed in the goitrous tissue of one patient. Human cells were transfected with expression vectors containing mutated and wild-type human TG cDNA. Results: All patients had an absent rise of serum TG after stimulation with recombinant human TSH. Sequence analysis revealed three previously described mutations (p. A2215D, p. R277X, and g. IVS30 + 1G > T), and two novel mutations (p. Q2142X and g. IVS46-1G > A). Two known (g. IVS30 + 1G/p. A2215D and p. A2215D/p. R277X) and one novel (p. R277X/g. IVS46-1G > A) compound heterozygous constellations were also identified. Functional analysis indicated deficiency in TG synthesis, reduction of TG secretion, and retention of the mutant TG within the cell, leading to an endoplasmic reticulum storage disease, whereas small amounts of mutant TG were still secreted within the cell system. Conclusion: All studied patients were either homozygous or heterozygous for TG gene mutations. Two novel mutations have been detected, and we show that TG mutation p. A2215D promotes the retention of TG within the endoplasmic reticulum and reduces TG synthesis and secretion, causing mild hypothyroidism. In the presence of sufficient iodine supply, some patients with TG mutations are able to compensate the impaired hormonogenesis and generate thyroid hormone. (J Clin Endocrinol Metab 94: 2938-2944, 2009)

Molecular characterization of the Hepatitis B virus genotypes in Colombia: A Bayesian inference on the genotype F

Hepatitis B is a worldwide health problem affecting about 2 billion people and more than 350 million are chronic carriers of the virus. Nine HBV genotypes (A to I) have been described. The geographical distribution of HBV genotypes is not completely understood due to the limited number of samples from some parts of the world. One such example is Colombia, in which few studies have described the HBV genotypes. In this study, we characterized HBV genotypes in 143 HBsAg-positive volunteer blood donors from Colombia. A fragment of 1306 bp partially comprising HBsAg and the DNA polymerase coding regions (S/POL) was amplified and sequenced. Bayesian phylogenetic analyses were conducted using the Markov Chain Monte Carlo (MCMC) approach to obtain the maximum clade credibility (MCC) tree using BEAST v.1.5.3. Of all samples, 68 were positive and 52 were successfully sequenced. Genotype F was the most prevalent in this population (77%) - subgenotypes F3 (75%) and Fib (2%). Genotype G (7.7%) and subgenotype A2 (15.3%) were also found. Genotype G sequence analysis suggests distinct introductions of this genotype in the country. Furthermore, we estimated the time of the most recent common ancestor (TMRCA) for each HBV/F subgenotype and also for Colombian F3 sequences using two different datasets: (i) 77 sequences comprising 1306 bp of S/POL region and (ii) 283 sequences comprising 681 bp of S/POL region. We also used two other previously estimated evolutionary rates: (i) 2.60 x 10(-4) s/s/y and (ii) 1.5 x 10(-5) s/s/y. Here we report the HBV genotypes circulating in Colombia and estimated the TMRCA for the four different subgenotypes of genotype F. (C) 2010 Elsevier B.V. All rights reserved.

Conservation of hepatitis C virus nonstructural protein 3 amino acid sequence in viral isolates during liver transplantation

As a consequence of selective pressure exerted by the immune response during hepatitis C virus (HCV) infection, a high rate of nucleotide mutations in the viral genome is observed which leads to the emergence of viral escape mutants. The aim of this study was to evaluate the evolution of the amino acid (aa) sequence of the HCV nonstructural protein 3 (NS3) in viral isolates after liver transplantation. Six patients with HCV-induced liver disease undergoing liver transplantation (LT) were followed up for sequence analysis. Hepatitis C recurrence was observed in all patients after LT. The rate of synonymous (dS) nucleotide substitutions was much higher than that of nonsynonymous (dN) ones in the NS3 encoding region. The high values of the dS/dN ratios suggest no sustained adaptive evolution selection pressure and, therefore, absence of specific NS3 viral populations. Clinical genotype assignments were supported by phylogenetic analysis. Serial samples from each patient showed lower mean nucleotide genetic distance when compared with samples of the same HCV genotype and subtype. The NS3 samples studied had an N-terminal aa sequence with several differences as compared with reference ones, mainly in genotype 1b-infected patients. After LT, as compared with the sequences before, a few reverted aa substitutions and several established aa substitutions were observed at the N-terminal of NS3. Sites described to be involved in important functions of NS3, notably those of the catalytic triad and zinc binding, remained unaltered in terms of aa sequence. Rare or frequent aa substitutions occurred indiscriminately in different positions. Several cytotoxic T lymphocyte epitopes described for HCV were present in our 1b samples. Nevertheless, the deduced secondary structure of the NS3 protease showed a few alterations in samples from genotype 3a patients, but none were seen in 1b cases. Our data, obtained from patients under important selective pressure during LT, show that the NS3 protease remains well conserved, mainly in HCV 3a patients. It reinforces its potential use as an antigenic candidate for further studies aiming at the development of a protective immune response.

GP5+/6+ SYBR Green methodology for simultaneous screening and quantification of human papillomavirus

Background: Detection and quantification of human papillomavirus (HPV) may help in predicting the evolution of HPV infection and progression of associated lesions. Objectives: We propose a novel protocol using consensus primers GP5+/6+ in a SYBR Green quantitative real-time (Q-RT) polymerase chain reaction (PCR). The strategy permits screening for HPV infection and viral load quantification simultaneously. Study design: DNA from 153 archived cervical samples, previously tested for HPV detection by GP5+/6+ PCR and typed by EIA-RLB (enzyme immunoassay-reverse line blot) or sequence analysis, was analysed using SYBR Green Q-RT PCR. Melting temperature assay (T(m)) and cycle threshold (C(t)) were used to evaluate HPV positivity and viral load. The T(m) in the range of 77-82 degrees C was considered to be positive for HPV-DNA. HPV results generated through GP5+/6+ conventional PCR were considered the gold standard against which sensitivity and specificity of our assay were measured. Results: Out of 104 HPV positive samples, 100 (96.2%) were also determined as positive by SYBR Green Q-RT PCR; of the 49 HPV-negative samples, all were determined as negative. There was an excellent positivity agreement (K = 0.94) between the SYBR Green Q-RT and the previous methods employed. The specificity and sensitivity were 100% and 96.2%, respectively. Comparison of SYBR Green Q-RT and TaqMan oligo-probe technologies gave an excellent concordance (pc = 0.95) which validated the proposed strategy. Conclusions: We propose a sensitive and easy-to-perform technique for HPV screening and viral load quantification simultaneously. (C) 2009 Elsevier B.V. All rights reserved.