Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance

Background: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. Results: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. Conclusions: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

Detection of polymorphisms of the mtDNA control region of Caretta caretta (Testudines: Cheloniidae) by PCR-SSCP

Marine turtles are increasingly being threatened worldwide by anthropogenic activities. Better understanding of their life cycle, behavior and population structure is imperative for the design of adequate conservation strategies. The mtDNA control region is a fast-evolving matrilineal marker that has been employed in the study of marine turtle populations. We developed and tested a simple molecular tracing system for Caretta caretta mtDNA haplotypes by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Using this technique, we were able to distinguish the SSCP patterns of 18 individuals of the haplotypes CC-A4, CC-A24 and CCxLO, which are commonly found in turtles sampled on the Brazilian coast. When we analyzed 15 turtles with previously unknown sequences, we detected two other haplotypes, in addition to the other four. Based on DNA sequencing, they were identified as the CC-A17 and CC-A1 haplotypes. Further analyses were made with the sea turtles, Chelonia mydas (N = 8), Lepidochelys olivacea (N = 3) and Eretmochelys imbricata (N = 1), demonstrating that the PCR-SSCP technique is able to distinguish intra-and interspecific variation in the family Cheloniidae. We found that this technique can be useful for identifying sea turtle mtDNA haplotypes, reducing the need for sequencing.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

The Monocarboxylate Transporter 8 and L-Type Amino Acid Transporters 1 and 2 Are Expressed in Mouse Skeletons and in Osteoblastic MC3T3-E1 Cells

Background: Several plasma membrane transporters have been shown to mediate the cellular influx and/or efflux of iodothyronines, including the sodium-independent organic anion co-transporting polypeptide 1 (OATP1), the sodium taurocholate co-transporting polypeptide (NTCP), the L-type amino acid transporter 1 (LAT1) and 2 (LAT2), and the monocarboxylate transporter 8 (MCT8). The aim of this study was to investigate if the mRNAs of these transporters were expressed and regulated by thyroid hormone (TH) in mouse calvaria-derived osteoblastic MC3T3-E1 cells and in the fetal and postnatal bones of mice. Methods: The mRNA expression of the iodothyronine transporters was investigated with real-time polymerase chain reaction analysis in euthyroid and hypothyroid fetuses and litters of mice and in MC3T3-E1 cells treated with increasing doses of triiodothyronine (T(3); 10(-10) to 10(-6) M) or with 10(-8) M T(3) for 1-9 days. Results: MCT8, LAT1, and LAT2 mRNAs were detected in fetal and postnatal femurs and in MC3T3-E1 cells, while OATP1 and NTCP mRNAs were not. LAT1 and LAT2 mRNAs were not affected by TH status in vivo or in vitro or by the stage of bone development or osteoblast maturation (analyzed by the expression of osteocalcin and alkaline phosphatase, which are key markers of osteoblastic differentiation). In contrast, the femoral mRNA expression of MCT8 decreased significantly during post-natal development, whereas MCT8 mRNA expression increased as MC3T3-E1 cells differentiated. We also showed that MCT8 mRNA was up-regulated in the femur of hypothyroid animals, and that it was down-regulated by treatment with T(3) in MC3T3-E1 cells. Conclusions: This is the first study to demonstrate the mRNA expression of LAT1, LAT2, and MCT8 in the bone tissue of mice and in osteoblast-like cells. In addition, the pattern of MCT8 expression observed in vivo and in vitro suggests that MCT8 may be important to modulate TH effects on osteoblast differentiation and on bone development and metabolism.

High Iodine Concentration Attenuates RET/PTC3 Oncogene Activation in Thyroid Follicular Cells

Background: Papillary thyroid carcinoma (PTC) is frequently associated with a RET gene rearrangement that generates a RET/PTC oncogene. RET/PTC is a fusion of the tyrosine kinase domain of RET to the 50 portion of a different gene. This fusion results in a constitutively active MAPK pathway, which plays a key role in PTC development. The RET/PTC3 fusion is primarily associated with radiation-related PTC. Epidemiological studies show a lower incidence of PTC in radiation-exposed regions that are associated with an iodine-rich diet. Since the influence of excess iodine on the development of thyroid cancer is still unclear, the aim of this study is to evaluate the effect of high iodine concentrations on RET/PTC3-activated thyroid cells. Methods: PTC3-5 cells, a rat thyroid cell lineage harboring doxycycline-inducible RET/PTC3, were treated with 10(-3) M NaI. Cell growth was analyzed by cell counting and the MTT assay. The expression and phosphorylation state of MAPK pathway-related (Braf, Erk, pErk, and pRet) and thyroid-specific (natrium-iodide symporter [Nis] and thyroid-stimulating hormone receptor [Tshr]) proteins were analyzed by Western blotting. Thyroid-specific gene expression was further analyzed by quantitative reverse transcription (RT)-polymerase chain reaction. Results: A significant inhibition of proliferation was observed, along with no significant variation in cell death rate, in the iodine-treated cells. Further, iodine treatment attenuated the loss of Nis and Tshr gene and protein expression induced by RET/PTC3 oncogene induction. Finally, iodine treatment reduced Ret and Erk phosphorylation, without altering Braf and Erk expression. Conclusion: Our results indicate an antioncogenic role for excess iodine during thyroid oncogenic activation. These findings contribute to a better understanding of the effect of iodine on thyroid follicular cells, particularly how it may play a protective role during RET/PTC3 oncogene activation.

Complete Genome Viral Phylogenies Suggests the Concerted Evolution of Regulatory Cores and Accessory Satellites

We consider the concerted evolution of viral genomes in four families of DNA viruses. Given the high rate of horizontal gene transfer among viruses and their hosts, it is an open question as to how representative particular genes are of the evolutionary history of the complete genome. To address the concerted evolution of viral genes, we compared genomic evolution across four distinct, extant viral families. For all four viral families we constructed DNA-dependent DNA polymerase-based (DdDp) phylogenies and in addition, whole genome sequence, as quantitative descriptions of inter-genome relationships. We found that the history of the polymerase gene was highly predictive of the history of the genome as a whole, which we explain in terms of repeated, co-divergence events of the core DdDp gene accompanied by a number of satellite, accessory genetic loci. We also found that the rate of gene gain in baculovirus and poxviruses proceeds significantly more quickly than the rate of gene loss and that there is convergent acquisition of satellite functions promoting contextual adaptation when distinct viral families infect related hosts. The congruence of the genome and polymerase trees suggests that a large set of viral genes, including polymerase, derive from a phylogenetically conserved core of genes of host origin, secondarily reinforced by gene acquisition from common hosts or co-infecting viruses within the host. A single viral genome can be thought of as a mutualistic network, with the core genes acting as an effective host and the satellite genes as effective symbionts. Larger virus genomes show a greater departure from linkage equilibrium between core and satellites functions.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

High Prevalence of bla(CTX-M) Extended Spectrum Beta-Lactamase Genes in Klebsiella pneumoniae Isolates from a Tertiary Care Hospital: First report of bla(SHV-12), bla(SHV-31), bla(SHV-38), and bla(CTX-M-15) in Brazil

The aim of this study was to investigate the presence and prevalence of bla(TEM), bla(SHV), and bla(CTX-M) and bla(GES)-like genes, responsible for extended spectrum beta-lactamases (ESBLs) production in clinical isolates of Klebsiella pneumoniae collected from a Brazilian tertiary care hospital. Sixty-five ESBL producing K. pneumoniae isolates, collected between 2005 and 2007, were screened by polymerase chain reaction (PCR). Identification of bla genes was achieved by sequencing. Genotyping of ESBL producing K. pneumoniae was performed by the enterobacterial repetitive intergenic consensus-PCR with cluster analysis by the Dice coefficient. The presence of genes encoding ESBLs was confirmed in 59/65 (90.8%) isolates, comprising 20 bla(CTX-M-2), 14 bla(CTX-M-59), 12 bla(CTX-M-15), 9 bla(SHV-12), 1 bla(SHV-2), 1 bla(SHV-2a), 1 bla(SHV-5), and 1 bla(SHV-31) genes. The ESBL genes bla(SHV-12), bla(SHV-31), and bla(CTX-M-15), and the chromosome-encoded SHV-type beta-lactamase capable of hydrolyzing imipenem were detected in Brazil for the first time. The analysis of the enterobacterial repetitive intergenic consensus-PCR band patterns revealed a high rate of multiclonal bla(CTX-M) carrying K. pneumoniae isolates (70.8%), suggesting that dissemination of encoding plasmids is likely to be the major cause of the high prevalence of these genes among the K. pneumoniae isolates considered in this study.

Social Networks Shape the Transmission Dynamics of Hepatitis C Virus

Hepatitis C virus (HCV) infects 170 million people worldwide, and is a major public health problem in Brazil, where over 1% of the population may be infected and where multiple viral genotypes co-circulate. Chronically infected individuals are both the source of transmission to others and are at risk for HCV-related diseases, such as liver cancer and cirrhosis. Before the adoption of anti-HCV control measures in blood banks, this virus was mainly transmitted via blood transfusion. Today, needle sharing among injecting drug users is the most common form of HCV transmission. Of particular importance is that HCV prevalence is growing in non-risk groups. Since there is no vaccine against HCV, it is important to determine the factors that control viral transmission in order to develop more efficient control measures. However, despite the health costs associated with HCV, the factors that determine the spread of virus at the epidemiological scale are often poorly understood. Here, we sequenced partial NS5b gene sequences sampled from blood samples collected from 591 patients in Sao Paulo state, Brazil. We show that different viral genotypes entered Sao Paulo at different times, grew at different rates, and are associated with different age groups and risk behaviors. In particular, subtype 1b is older and grew more slowly than subtypes 1a and 3a, and is associated with multiple age classes. In contrast, subtypes 1a and 3b are associated with younger people infected more recently, possibly with higher rates of sexual transmission. The transmission dynamics of HCV in Sao Paulo therefore vary by subtype and are determined by a combination of age, risk exposure and underlying social network. We conclude that social factors may play a key role in determining the rate and pattern of HCV spread, and should influence future intervention policies.

Anti-Anopheles darlingi saliva antibodies as marker of Plasmodium vivax infection and clinical immunity in the Brazilian Amazon

Background: Despite governmental and private efforts on providing malaria control, this disease continues to be a major health threat. Thus, innovative strategies are needed to reduce disease burden. The malaria vectors, through the injection of saliva into the host skin, play important role on disease transmission and may influence malaria morbidity. This study describes the humoral immune response against Anopheles (An.) darlingi saliva in volunteers from the Brazilian Amazon and addresses the association between levels of specific antibodies and clinical presentation of Plasmodium (P.) vivax infection. Methods: Adult volunteers from communities in the Rondonia State, Brazil, were screened in order to assess the presence of P. vivax infection by light microscopy and nested PCR. Non-infected volunteers and individuals with symptomatic or symptomless infection were randomly selected and plasma collected. An. darlingi salivary gland sonicates (SGS) were prepared and used to measure anti-saliva antibody levels. Plasma interleukin (IL)-10 and interferon (IFN)-gamma levels were also estimated and correlated to anti-SGS levels. Results: Individuals infected with P. vivax presented higher levels of anti-SGS than non-infected individuals and antibody levels could discriminate infection. Furthermore, anti-saliva antibody measurement was also useful to distinguish asymptomatic infection from non-infection, with a high likelihood ratio. Interestingly, individuals with asymptomatic parasitaemia presented higher titers of anti-SGS and lower IFN-gamma/IL-10 ratio than symptomatic ones. In P. vivax-infected asymptomatic individuals, the IFN-gamma/IL-10 ratio was inversely correlated to anti-SGS titers, although not for while in symptomatic volunteers. Conclusion: The estimation of anti-An. darlingi antibody levels can indicate the probable P. vivax infection status and also could serve as a marker of disease severity in this region of Brazilian Amazon.

Hepatitis B Infection Is Associated with Asymptomatic Malaria in the Brazilian Amazon

Background: Areas that are endemic for malaria are also highly endemic for hepatitis B virus (HBV) infection. Nevertheless, it is unknown whether HBV infection modifies the clinical presentation of malaria. This study aimed to address this question. Methodology and Findings: An observational study of 636 individuals was performed in Rondonia, western Amazon, Brazil between 2006 and 2007. Active and passive case detections identified Plasmodium infection by field microscopy and nested Polymerase Chain Reaction (PCR). HBV infections were identified by serology and confirmed by real-time PCR. Epidemiological information and plasma cytokine profiles were studied. The data were analyzed using adjusted multinomial logistic regression. Plasmodium-infected individuals with active HBV infection were more likely to be asymptomatic (OR: 120.13, P < 0.0001), present with lower levels of parasitemia and demonstrate a decreased inflammatory cytokine profile. Nevertheless, co-infected individuals presented higher HBV viremia. Plasmodium parasitemia inversely correlated with plasma HBV DNA levels (r=-0.6; P=0.0003). Conclusion: HBV infection diminishes the intensity of malaria infection in individuals from this endemic area. This effect seems related to cytokine balance and control of inflammatory responses. These findings add important insights to the understanding of the factors affecting the clinical outcomes of malaria in endemic regions.

Identification of COL6A1 as a differentially expressed gene in human astrocytomas

Diffuse infiltrating gliomas are the most common tumors of the central nervous system. Gliomas are classified by the WHO according to their histopathological and clinical characteristics into four classes: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). Several genes have already been correlated with astrocytomas, but many others are yet to be uncovered. By analyzing the public SAGE data from 21 patients, comprising low malignant grade astrocytomas and glioblastomas, we found COL6A1 to be differentially expressed, confirming this finding by real time RT-PCR in 66 surgical samples. To the best of our knowledge, COL6A1 has never been described in gliomas. The expression of this gene has significantly different means when normal glia is compared with low-grade astrocytomas (grades I and II) and high-grade astrocytomas (grades III and IV), with a tendency to be greater in higher grade samples, thus rendering it a powerful tumor marker.

Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly in patients with heart failure

Background -: Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly and Thr164Ile were suggested to have an effect in heart failure. We evaluated these polymorphisms relative to clinical characteristics and prognosis of alarge cohort of patients with heart failure of different etiologies. Methods -: We studied 501 patients with heart failure of different etiologies. Mean age was 58 years (standard deviation 14.4 years), 298 (60%) were men. Polymorphisms were identified by polymerase chain reaction-restriction fragment length polymorphism. Results -: During the mean follow-up of 12.6 months (standard deviation 10.3 months), 188 (38%) patients died. Distribution of genotypes of polymorphism Arg16Gly was different relative to body mass index (chi(2) = 9.797; p = 0.04). Overall the probability of survival was not significantly predicted by genotypes of Gln27Glu, Arg16Gly, or Thr164Ile. Allele and haplotype analysis also did not disclose any significant difference regarding mortality. Exploratory analysis through classification trees pointed towards a potential association between the Gln27Glu polymorphism and mortality in older individuals. Conclusion -: In this study sample, we were not able to demonstrate an overall influence of polymorphisms Gln27Glu and Arg16Gly of beta-2 receptor gene on prognosis. Nevertheless, Gln27Glu polymorphism may have a potential predictive value in older individuals.

Identification of archaeal proteins that affect the exosome function in vitro

Background: The archaeal exosome is formed by a hexameric RNase PH ring and three RNA binding subunits and has been shown to bind and degrade RNA in vitro. Despite extensive studies on the eukaryotic exosome and on the proteins interacting with this complex, little information is yet available on the identification and function of archaeal exosome regulatory factors. Results: Here, we show that the proteins PaSBDS and PaNip7, which bind preferentially to poly-A and AU-rich RNAs, respectively, affect the Pyrococcus abyssi exosome activity in vitro. PaSBDS inhibits slightly degradation of a poly-rA substrate, while PaNip7 strongly inhibits the degradation of poly-A and poly-AU by the exosome. The exosome inhibition by PaNip7 appears to depend at least partially on its interaction with RNA, since mutants of PaNip7 that no longer bind RNA, inhibit the exosome less strongly. We also show that FITC-labeled PaNip7 associates with the exosome in the absence of substrate RNA. Conclusions: Given the high structural homology between the archaeal and eukaryotic proteins, the effect of archaeal Nip7 and SBDS on the exosome provides a model for an evolutionarily conserved exosome control mechanism.

Sequence-specific electrochemical detection of Alicyclobacillus acidoterrestris DNA using electroconductive polymer-modified fluorine tin oxide electrodes

This study outlines the quantification of low levels of Alicyclobacillus acidoterrestris in pure cultures, since this bacterium is not inactivated by pasteurization and may remain in industrialized foods and beverages. Electroconductive polymer-modified fluorine tin oxide (FTO) electrodes and multiple nanoparticle labels were used for biosensing. The detection of A. acidoterrestris in pure cultures was performed by reverse transcription polymerase chain reaction (RT-PCR) and the sensitivity was further increased by asymmetric nested RT-PCR using electrochemical detection for quantification of the amplicon. The quantification of nested RT-PCR products by Ag/Au-based electrochemical detection was able to detect 2 colony forming units per mL (CFU mL(-1)) of spores in pure culture and low detection and quantification limits (7.07 and 23.6 nM, respectively) were obtained for the target A. acidoterrestris on the electrochemical detection bioassay.