Variability in Estrogen-Metabolizing Genes and Their Association with Genomic Instability in Untreated Breast Cancer Patients and Healthy Women

In the present study, we investigated the relationship between polymorphisms in the estrogen-metabolizing genes CYP17, CYP1B1, CYP1A1, and COMT and genomic instability in the peripheral blood lymphocytes of 62 BC patients and 62 controls considering that increased or prolonged exposure to estrogen can damage the DNA molecule and increase the genomic instability process in breast tissue. Our data demonstrated increased genomic instability in BC patients and that individuals with higher frequencies of MN exhibited higher risk to BC when belonging Val/Met genotype of the COMT gene. We also observed that CYP17 and CYP1A1 polymorphisms can modify the risk to BC depending on the menopause status. We can conclude that the genetic background in estrogen metabolism pathway can modulate chromosome damage in healthy controls and patients and thereby influence the risk to BC. These findings suggest the importance to ally biomarkers of susceptibility and effects to estimate risk groups.

Tumor slices as a model to evaluate doxorubicin in vitro treatment and expression of trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 in canine mammary gland cancer

Background: In women with breast cancer submitted to neoadjuvant chemotherapy based in doxorubicin, tumor expression of groups of three genes (PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2) have classified them as responsive or resistant. We have investigated whether expression of these trios of genes could predict mammary carcinoma response in dogs and whether tumor slices, which maintain epithelial-mesenchymal interactions, could be used to evaluate drug response in vitro. Methods: Tumors from 38 dogs were sliced and cultured with or without doxorubicin 1 mu M for 24 h. Tumor cells were counted by two observers to establish a percentage variation in cell number, between slices. Based on these results, a reduction in cell number between treated and control samples >= 21.7%, arbitrarily classified samples, as drug responsive. Tumor expression of PRSS11, MTSS1, CLPTM1 and SMYD2, was evaluated by real time PCR. Relative expression results were then transformed to their natural logarithm values, which were spatially disposed according to the expression of trios of genes, comprising PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. Fisher linear discrimination test was used to generate a separation plane between responsive and non-responsive tumors. Results: Culture of tumor slices for 24 h was feasible. Nine samples were considered responsive and 29 non-responsive to doxorubicin, considering the pre-established cut-off value of cell number reduction = 21.7%, between doxorubicin treated and control samples. Relative gene expression was evaluated and tumor samples were then spatially distributed according to the expression of the trios of genes: PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. A separation plane was generated. However, no clear separation between responsive and non-responsive samples could be observed. Conclusion: Three-dimensional distribution of samples according to the expression of the trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 could not predict doxorubicin in vitro responsiveness. Short term culture of mammary gland cancer slices may be an interesting model to evaluate chemotherapy activity.

Sequential analysis of global gene expression profiles in immature and in vitro matured bovine oocytes: potential molecular markers of oocyte maturation

Background: Without intensive selection, the majority of bovine oocytes submitted to in vitro embryo production (IVP) fail to develop to the blastocyst stage. This is attributed partly to their maturation status and competences. Using the Affymetrix GeneChip Bovine Genome Array, global mRNA expression analysis of immature (GV) and in vitro matured (IVM) bovine oocytes was carried out to characterize the transcriptome of bovine oocytes and then use a variety of approaches to determine whether the observed transcriptional changes during IVM was real or an artifact of the techniques used during analysis. Results: 8489 transcripts were detected across the two oocyte groups, of which similar to 25.0% (2117 transcripts) were differentially expressed (p < 0.001); corresponding to 589 over-expressed and 1528 under-expressed transcripts in the IVM oocytes compared to their immature counterparts. Over expression of transcripts by IVM oocytes is particularly interesting, therefore, a variety of approaches were employed to determine whether the observed transcriptional changes during IVM were real or an artifact of the techniques used during analysis, including the analysis of transcript abundance in oocytes in vitro matured in the presence of a-amanitin. Subsets of the differentially expressed genes were also validated by quantitative real-time PCR (qPCR) and the gene expression data was classified according to gene ontology and pathway enrichment. Numerous cell cycle linked (CDC2, CDK5, CDK8, HSPA2, MAPK14, TXNL4B), molecular transport (STX5, STX17, SEC22A, SEC22B), and differentiation (NACA) related genes were found to be among the several over-expressed transcripts in GV oocytes compared to the matured counterparts, while ANXA1, PLAU, STC1and LUM were among the over-expressed genes after oocyte maturation. Conclusion: Using sequential experiments, we have shown and confirmed transcriptional changes during oocyte maturation. This dataset provides a unique reference resource for studies concerned with the molecular mechanisms controlling oocyte meiotic maturation in cattle, addresses the existing conflicting issue of transcription during meiotic maturation and contributes to the global goal of improving assisted reproductive technology.

Association of SNPs on CAPN1 and CAST genes with tenderness in Nellore cattle

We examined whether single-nucleotide polymorphisms (SNPs) in the calpain (CAPN) and calpastatin (CAST) genes, described from Bos primigenius taurus, are polymorphic in Nellore cattle. We also looked for a possible association of linkage disequilibrium of this polymorphism with tenderness of the longissimus dorsi muscle after 7, 14 and 21 days of postmortem aging in 638 purebred Nellore bulls. Meat tenderness was measured as Warner-Bratzler shear force. Additive and dominance effects were tested for SNPs of the three genotypic classes; the substitution effect was tested for SNPs with missing genotypic classes. Genotypic and gene frequencies were also calculated for the different SNPs. An increase in tenderness was observed from 7 to 21 days; the average values for shear force at 7, 14 and 21 days of aging were 5.92 +/- 0.06, 4.92 +/- 0.05, and 4.38 +/- 0.04 kg, respectively. All markers showed polymorphism, but there was no CC genotype for CAPN316, and few animals showed the AA genotype for CAPN530. The alleles CAPN4751, UOGCAST1, and WSUCAST were found to have additive and dominance effects for shear force at 7, 14 and 21 days, while CAPN316 showed a substitution effect for shear force at 7 and 21 days. An additive-by-additive epistatic interaction was observed between CAPN4751 and markers on the CAST gene. In conclusion, these markers should be considered for use in breeding programs.

Effects of polymorphisms of LHR and FSHR genes on sexual precocity in a Bos taurus x Bos indicus beef composite population

The purpose of the present research was to investigate the effects of polymorphisms of luteinizing hormone receptor (LHR) and follicle-stimulating hormone receptor (FSHR) genes, evaluated by polymerase chain reaction-restriction fragment length polymorphism in European-Zebu composite beef heifers from six different breed compositions. The polymorphism site analysis from digestion with HhaI and AluI restriction endonucleases allowed the genotype identification for LHR (TT, CT and CC) and FSHR (GG, CG and CC) genes. A high frequency of heterozygous animals was recorded in all breed compositions for both genes, except in two compositions for LHR. The probability of pregnancy (PP) at first breeding was used to evaluate the polymorphism effect on sexual precocity. The PP was analyzed as a binary trait, with a value of 1 (success) assigned to heifers that were diagnosed pregnant by rectal palpation and a value of 0 (failure) assigned to those that were not pregnant at that time. Heterozygous heifers showed a higher pregnancy rate (67 and 66% for LHR and FSHR genes, respectively), but no significant effects were observed for the genes studied (P=0.9188 and 0.8831 for LHR and FSHR, respectively) on the PP. These results do not justify the inclusion of LHR and FSHR restriction fragment length polymorphism markers in selection programs for sexual precocity in beef heifers. Nevertheless, these markers make possible the genotype characterization and may be used in additional studies to evaluate the genetic structure in other bovine populations.

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.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

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.

Identification and Phylogenetic Analysis of Heme Synthesis Genes in Trypanosomatids and Their Bacterial Endosymbionts

It has been known for decades that some insect-infecting trypanosomatids can survive in culture without heme supplementation while others cannot, and that this capability is associated with the presence of a betaproteobacterial endosymbiont in the flagellate's cytoplasm. However, the specific mechanisms involved in this process remained obscure. In this work, we sequence and phylogenetically analyze the heme pathway genes from the symbionts and from their hosts, as well as from a number of heme synthesis-deficient Kinetoplastida. Our results show that the enzymes responsible for synthesis of heme are encoded on the symbiont genomes and produced in close cooperation with the flagellate host. Our evidence suggests that this synergistic relationship is the end result of a history of extensive gene loss and multiple lateral gene transfer events in different branches of the phylogeny of the Trypanosomatidae.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

New genes of Xanthomonas citri subsp citri involved in pathogenesis and adaptation revealed by a transposon-based mutant library

Background: Citrus canker is a disease caused by the phytopathogens Xanthomonas citri subsp. citri, Xanthomonas fuscans subsp. aurantifolli and Xanthomonas alfalfae subsp. citrumelonis. The first of the three species, which causes citrus bacterial canker type A, is the most widely spread and severe, attacking all citrus species. In Brazil, this species is the most important, being found in practically all areas where citrus canker has been detected. Like most phytobacterioses, there is no efficient way to control citrus canker. Considering the importance of the disease worldwide, investigation is needed to accurately detect which genes are related to the pathogen-host adaptation process and which are associated with pathogenesis. Results: Through transposon insertion mutagenesis, 10,000 mutants of Xanthomonas citri subsp. citri strain 306 (Xcc) were obtained, and 3,300 were inoculated in Rangpur lime (Citrus limonia) leaves. Their ability to cause citrus canker was analyzed every 3 days until 21 days after inoculation; a set of 44 mutants showed altered virulence, with 8 presenting a complete loss of causing citrus canker symptoms. Sequencing of the insertion site in all 44 mutants revealed that 35 different ORFs were hit, since some ORFs were hit in more than one mutant, with mutants for the same ORF presenting the same phenotype. An analysis of these ORFs showed that some encoded genes were previously known as related to pathogenicity in phytobacteria and, more interestingly, revealed new genes never implicated with Xanthomonas pathogenicity before, including hypothetical ORFs. Among the 8 mutants with no canker symptoms are the hrpB4 and hrpX genes, two genes that belong to type III secretion system (TTSS), two hypothetical ORFS and, surprisingly, the htrA gene, a gene reported as involved with the virulence process in animal-pathogenic bacteria but not described as involved in phytobacteria virulence. Nucleic acid hybridization using labeled cDNA probes showed that some of the mutated genes are differentially expressed when the bacterium is grown in citrus leaves. Finally, comparative genomic analysis revealed that 5 mutated ORFs are in new putative pathogenicity islands. Conclusion: The identification of these new genes related with Xcc infection and virulence is a great step towards the understanding of plant-pathogen interactions and could allow the development of strategies to control citrus canker.

The Influence of Different Land Uses on the Structure of Archaeal Communities in Amazonian Anthrosols Based on 16S rRNA and amoA Genes

Soil from the Amazonian region is usually regarded as unsuitable for agriculture because of its low organic matter content and low pH; however, this region also contains extremely rich soil, the Terra Preta Anthrosol. A diverse archaeal community usually inhabits acidic soils, such as those found in the Amazon. Therefore, we hypothesized that this community should be sensitive to changes in the environment. Here, the archaeal community composition of Terra Preta and adjacent soil was examined in four different sites in the Brazilian Amazon under different anthropic activities. The canonical correspondence analysis of terminal restriction fragment length polymorphisms has shown that the archaeal community structure was mostly influenced by soil attributes that differentiate the Terra Preta from the adjacent soil (i.e., pH, sulfur, and organic matter). Archaeal 16S rRNA gene clone libraries indicated that the two most abundant genera in both soils were Candidatus nitrosphaera and Canditatus nitrosocaldus. An ammonia monoxygenase gene (amoA) clone library analysis indicated that, within each site, there was no significant difference between the clone libraries of Terra Preta and adjacent soils. However, these clone libraries indicated there were significant differences between sites. Quantitative PCR has shown that Terra Preta soils subjected to agriculture displayed a higher number of amoA gene copy numbers than in adjacent soils. On the other hand, soils that were not subjected to agriculture did not display significant differences on amoA gene copy numbers between Terra Preta and adjacent soils. Taken together, our findings indicate that the overall archaeal community structure in these Amazonian soils is determined by the soil type and the current land use.