Reference gene selection in human periodontal ligament cells

Background: Real-time quantitative PCR (qPCR) is a highly sensitive and specific method which is used extensively for determining gene expression profiles in a variety of cell and tissue types. In order to obtain accurate and reliable gene expression quantification, qPCR data are generally normalised against so-called reference or housekeeping genes. Ideally, reference genes should have abundant and stable RNA transcriptomes under the experimental conditions employed. However, reference genes are often selected rather arbitrarily and indeed some have been shown to have variable expression in a variety of in vitro experimental conditions.
Objective: The objective of the current study was to investigate reference gene expression in human periodontal ligament (PDL) cells in response to treatment with lipopolysaccharide (LPS).
Method: Primary human PDL cells were grown in Dulbecco’s Modified Eagle Medium with L-glutamine supplemented with 10% fetal bovine serum, 100UI/ml penicillin and 100µg/ml streptomycin. RNA was isolated using the RNeasy Mini Kit (Qiagen) and reverse transcribed using the QuantiTect Reverse Transcription Kit (Qiagen). The expression of a total of 19 reference genes was studied in the presence and absence of LPS treatment using the Roche Reference Gene Panel. Data were analysed using NormFinder and Bestkeeper validation programs.
Results: Treatment of human PDL cells with LPS resulted in changes in expression of several commonly used reference genes, including GAPDH. On the other hand the reference genes β-actin, G6PDH and 18S were identified as stable genes following LPS treatment.
Conclusion: Many of the reference genes studied were robust to LPS treatment (up to 100 ng/ml). However several commonly employed reference genes, including GAPDH varied with LPS treatment, suggesting they would not be ideal candidates for normalisation in qPCR gene expression studies.

Lancefield group C and G streptococci in human infection : molecular typing, virulence and antimicrobial susceptibility

Tese de doutoramento, Ciências e Tecnologias da Saúde (Microbiologia), Universidade de Lisboa, Faculdade de Medicina, 2014

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Characterizing the Expression of Cytochrome P450s in Breast Cancer Cells

Une résistance aux agents anticancéreux utilisés dans le traitement du cancer du sein est souvent associée à un échec de traitement. Des variations dans le devenir des agents anticancéreux dans l’organisme, sont des facteurs pouvant expliquer des phénomènes de résistance. Notre but était d’évaluer l’impact des isoenzymes du CYP450s, dans le métabolisme local des agents anticancéreux. Notre premier objectif était de valider un gène rapporteur pour nos analyses de PCR en temps réel. Pour ce faire, nous avons criblé l’expression de 6 gènes rapporteurs dans 23 lignées cellulaires. NUP-214 a été démontré comme étant le gène rapporteur le plus stable avec un écart-type de seulement 0.55 Ct. Notre deuxième objectif était de déterminer le niveau d’expression des ARNm de 19 isoformes du CYP450 dans plusieurs lignées cellulaires du cancer du sein. Les ARNm des CYP450s ont démontré une très grande variabilité entre les lignées cellulaires. Les isoformes CYP1B1 et CYP2J2 démontrent l’expression la plus importante pour la majorité des lignées. Notre troisième objectif était d’évaluer la corrélation entre l’expression des isoformes des CYP450s et leur activité métabolique en utilisant les substrats spécifiques du CYP1B1 et 2J2, 7-éthoxyrésorufine et ébastine, respectivement. Une forte corrélation (r2=0.99) fut observée entre l’activité métabolique vis-à-vis l’ébastine et l’expression du CYP2J2. De même, le métabolisme du 7-éthoxyrésorufine était fortement corrélé (r2=0.98) avec l’expression du CYP1B1. En résumé, ces résultats suggèrent que le métabolisme local des agents anticancéreux pourrait significativement moduler le devenir des agents anticancéreux dans l’organisme, et pourrait être ainsi, une source de résistance.

Soft topographic map for clustering and classification of bacteria

In this work a new method for clustering and building a topographic representation of a bacteria taxonomy is presented. The method is based on the analysis of stable parts of the genome, the so-called “housekeeping genes”. The proposed method generates topographic maps of the bacteria taxonomy, where relations among different type strains can be visually inspected and verified. Two well known DNA alignement algorithms are applied to the genomic sequences. Topographic maps are optimized to represent the similarity among the sequences according to their evolutionary distances. The experimental analysis is carried out on 147 type strains of the Gammaprotebacteria class by means of the 16S rRNA housekeeping gene. Complete sequences of the gene have been retrieved from the NCBI public database. In the experimental tests the maps show clusters of homologous type strains and present some singular cases potentially due to incorrect classification or erroneous annotations in the database.

Metallothionein I and II gene expression as a response to metal contamination in bank vole Clethrionomys glareolus

The expression of two metallothionein genes (Mt-I and Mt-II) in the liver, kidney, and gonad of bank voles collected at four metal-contaminated sites (Cd, Zn, Pb, and Fe) were measured using the quantitative real-time PCR method (QPCR). Relative Mt gene expression was calculated by applying a normalization factor (NF) using the expression of two housekeeping genes, ribosomal 18S and beta-actin. Relative Mt expression in tissues of animals from contaminated sites was up to 54.8-fold higher than those from the reference site for Mt-I and up to 91.6-fold higher for Mt-II. Mt-II gene expression in the livers of bank voles from contaminated sites was higher than Mt-I gene expression. Inversely, Mt-II expression in the kidneys of voles was lower than Mt-I expression. Positive correlations between cadmium levels in the tissues and Mt-I were obtained in all studied tissues. Zinc, which undergoes homeostatic regulation, correlated positively with both Mt-I and Mt-II gene expression only in the kidney. Results showed that animals living in chronically contaminated environments intensively activate detoxifying mechanisms such as metallothionein expression. This is the first time that QPCR techniques to measure MT gene expression have been applied to assess the impact of environmental metal pollution on field collected bank voles.

Plasmodium falciparum Accompanied the Human Expansion out of Africa

Plasmodium falciparum is distributed throughout the tropics and is responsible for an estimated 230 million cases of malaria every year, with a further 1.4 billion people at risk of infection [1-3]. Little is known about the genetic makeup of P. falciparum populations, despite variation in genetic diversity being a key factor in morbidity, mortality, and the success of malaria control initiatives. Here we analyze a worldwide sample of 519 P. falciparum isolates sequenced for two housekeeping genes (63 single nucleotide polymorphisms from around 5000 nucleotides per isolate). We observe a strong negative correlation between within-population genetic diversity and geographic distance from sub-Saharan Africa (R(2) = 0.95) over Africa, Asia, and Oceania. In contrast, regional variation in transmission intensity seems to have had a negligible impact on the distribution of genetic diversity. The striking geographic patterns of isolation by distance observed in P. falciparum mirror the ones previously documented in humans [4-7] and point to a joint sub-Saharan African origin between the parasite and its host. Age estimates for the expansion of P. falciparum further support that anatomically modern humans were infected prior to their exit out of Africa and carried the parasite along during their colonization of the world.

Frequency distribution of TATA box and extension sequences on human promoters

TATA box is one of the most important transcription factor binding sites. But the exact sequences of TATA box are still not very clear yet. In this study, we conducted a dedicated analysis on the frequency distribution of TATA Box and its extension sequences on human promoters. Sixteen TATA elements derived from TATA Box motif, TATAWAWN, were classified into three distribution patterns: peak, bottom-peak and bottom. Fourteen TATA extension sequences (up to two base extensions) were predicted to be the new TATA Box elements because of their high motif factors, which indicate their statistical significance. Statistical analysis on the promoters of mouse, zebrafish and drosophila melanogaster verified seven of these elements. It was also observed that the distribution of TATA elements on the promoters of housekeeping genes are very similar with their distribution on the promoters of tissue specific genes in human.

Análise por sequências multilocus de Xanthomonas fuscans subsp. aurantifolii

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular Evaluation of Developmental Competence of Oocytes Collected In Vivo from Buffalo and Bovine Heifers during Winter and Summer

Buffaloes and bovines are polyestrous and seasonal or annual livestock, respectively, that show reduced fertility during heat stress. To investigate whether reduced fertility is related to oocyte competence in both species, immature oocytes from buffalo and bovine heifers were collected during winter and summer and subjected to molecular analyses. In each season, heifers of both species had their follicular wave emergence synchronized with a standard protocol (Ferreira et al., 2011). Before being subjected to ovum pick up (OPU), cutaneous (CT; degrees C) and rectal (RT; degrees C) temperatures and respiratory rate (RR; breaths/min) were measured. Oocytes' RNA was extracted to evaluate the expression of target genes related to mtDNA replication/transcription (PPARGC1A, TFAM and MT-CO1), apoptosis (BAX and BCL2) and HS (HSP90AA1 and HSPA1AB). ACTB, HIST1H2AG and GAPDH were initially chosen as housekeeping genes. In buffaloes, CT (35.0 +/- 0.4 vs 23.8 +/- 0.5), RT (38.7 +/- 0.1 vs 38.0 +/- 0) and RR (21.3 +/- 1.2 vs 15.4 +/- 1.1) were higher during summer than winter. However, in bovine heifers, RT (38.7 +/- 0.1 vs 38.6 +/- 0.1) and RR (44.8 +/- 1.5 vs 40.6 +/- 1.5) were similar in both seasons, while CT (31.6 +/- 0.3 vs 30.2 +/- 0.3) was increased during summer. Reduced expression of ACTB, HIST1H2AG and GAPDH was evidenced during summer, disqualifying them as housekeeping genes. Similarly, the expression of all target genes was reduced during summer in oocytes of both species. In summary, physiological responses to heat stress seem to be more intense in buffalo than bovine heifers. However, in both species, negative effects of heat stress upon oocyte quality occur at the molecular level and affects genes related to several biological functions.

Curtobacterium flaccumfaciens PV. flaccumfaciens: aspectos epidemiológicos e variabilidade genética

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Rapporti fra fattori ambientali e proteine di parete in Bifidobacterium

The normal gut microbiota has several important functions in host physiology and metabolism, and plays a key role in health and disease. Bifidobacteria, which are indigenous components of gastrointestinal microbiota, may play an important role in maintaining the well-being of the host although its precise function is very difficult to study. Its physiological and biochemical activities are controlled by many factors, particularly diet and environment. Adherence and colonization capacity are considered as contributing factors for immune modulation, pathogen exclusion, and enhanced contact with the mucosa. In this way, bifidobacteria would fortify the microbiota that forms an integral part of the mucosal barrier and colonization resistance against pathogens. Bifidobacteria are not only subjected to stressful conditions in industrial processes, but also in nature, where the ability to respond quickly to stress is essential for survival. Bifidobacteria, like other microorganisms, have evolved sensing systems for/and defences against stress that allow them to withstand harsh conditions and sudden environmental changes. Bacterial stress responses rely on the coordinated expression of genes that alter various cellular processes and structures (e.g. DNA metabolism, housekeeping genes, cell-wall proteins, membrane composition) and act in concert to improve bacterial stress tolerance. The integration of these stress responses is accomplished by regulatory networks that allow the cell to react rapidly to various and sometimes complex environmental changes. This work examined the effect of important stressful conditions, such as changing pH and osmolarity, on the biosynthesis of cell wall proteins in B. pseudolongum subsp. globosum. These environmental factors all influence heavily the expression of BIFOP (BIFidobacterial Outer Proteins) in the cell-wall and can have an impact in the interaction with host. Also evidence has been collected linking the low concentration of sugar in the culture medium with the presence or absence of extracromosomal DNA.

Gene expression analysis in ‘Candidatus Phytoplasma mali’-resistant and -susceptible Malus genotypes

Apple proliferation (AP) disease is the most important graft-transmissible and vector-borne disease of apple in Europe. ‘Candidatus Phytoplasma mali’ (Ca. P. mali) is the causal agent of AP. Apple (Malus x domestica) and other Malus species are the only known woody hosts. In European apple orchards, the cultivars are mainly grafted on one rootstock, M. x domestica cv. M9. M9 like all other M. x domestica cultivars is susceptible to ‘Ca. P. mali’. Resistance to AP was found in the wild genotype Malus sieboldii (MS) and in MS-derived hybrids but they were characterised by poor agronomic value. The breeding of a new rootstock carrying the resistant and the agronomic traits was the major aim of a project of which this work is a part. The objective was to shed light into the unknown resistance mechanism. The plant-phytoplasma interaction was studied by analysing differences between the ‘Ca. P. mali’-resistant and -susceptible genotypes related to constitutively expressed genes or to induced genes during infection. The cDNA-Amplified Fragment Length Polymorphism (cDNA-AFLP) technique was employed in both approaches. Differences related to constitutively expressed genes were identified between two ‘Ca. P. mali’-resistant hybrid genotypes (4551 and H0909) and the ‘Ca. P. mali’-susceptible M9. 232 cDNA-AFLP bands present in the two resistant genotypes but absent in the susceptible one were isolated but several different products associated to each band were found. Therefore, two different macroarray hybridisation experiments were performed with the cDNA-AFLP fragments yielding 40 sequences encoding for genes of unknown function or a wide array of functions including plant defence. In the second approach, individuation and analysis of the induced genes was carried out exploiting an in vitro system in which healthy and ‘Ca. P. mali’-infected micropropagated plants were maintained under controlled conditions. Infection trials using in vitro grafting of ‘Ca. P. mali’ showed that the resistance phenotype could be reproduced in this system. In addition, ex vitro plants were generated as an independent control of the genes differentially expressed in the in vitro plants. The cDNA-AFLP analysis in in vitro plants yielded 63 bands characterised by over-expression in the infected state of both the H0909 and MS genotypes. The major part (37 %) of the associated sequences showed homology with products of unknown function. The other genes were involved in plant defence, energy transport/oxidative stress response, protein metabolism and cellular growth. Real-time qPCR analysis was employed to validate the differential expression of the genes individuated in the cDNA-AFLP analysis. Since no internal controls were available for the study of the gene expression in Malus, an analysis on housekeeping genes was performed. The most stably expressed genes were the elongation factor-1 α (EF1) and the eukaryotic translation initiation factor 4-A (eIF4A). Twelve out of 20 genes investigated through qPCR were significantly differentially expressed in at least one genotype either in in vitro plants or in ex vitro plants. Overall, about 20% of the genes confirmed their cDNA-AFLP expression pattern in M. sieboldii or H0909. On the contrary, 30 % of the genes showed down-regulation or were not differentially expressed. For the remaining 50 % of the genes a contrasting behaviour was observed. The qPCR data could be interpreted as follows: the phytoplasma infection unbalance photosynthetic activity and photorespiration down-regulating genes involved in photosynthesis and in the electron transfer chain. As result, and in contrast to M. x domestica genotypes, an up-regulation of genes of the general response against pathogens was found in MS. These genes involved the pathway of H2O2 and the production of secondary metabolites leading to the hypothesis that a response based on the accumulation of H2O2 in MS would be at the base of its resistance. This resembles a phenomenon known as “recovery” where the spontaneous remission of the symptoms is observed in old susceptible plants but occurring in a stochastic way while the resistance in MS is an inducible but stable feature. As additional product of this work three cDNA-AFLP-derived markers were developed which showed independent distribution among the seedlings of two breeding progenies and were associated to a genomic region characteristic of MS. These markers will contribute to the development of molecular markers for the resistance as well as to map the resistance on the Malus genome.

Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains

Mycobacterium abscessus, Mycobacterium bolletii, and Mycobacterium massiliense (Mycobacterium abscessus sensu lato) are closely related species that currently are identified by the sequencing of the rpoB gene. However, recent studies show that rpoB sequencing alone is insufficient to discriminate between these species, and some authors have questioned their current taxonomic classification. We studied here a large collection of M. abscessus (sensu lato) strains by partial rpoB sequencing (752 bp) and multilocus sequence analysis (MLSA). The final MLSA scheme developed was based on the partial sequences of eight housekeeping genes: argH, cya, glpK, gnd, murC, pgm, pta, and purH. The strains studied included the three type strains (M. abscessus CIP 104536(T), M. massiliense CIP 108297(T), and M. bolletii CIP 108541(T)) and 120 isolates recovered between 1997 and 2007 in France, Germany, Switzerland, and Brazil. The rpoB phylogenetic tree confirmed the existence of three main clusters, each comprising the type strain of one species. However, divergence values between the M. massiliense and M. bolletii clusters all were below 3% and between the M. abscessus and M. massiliense clusters were from 2.66 to 3.59%. The tree produced using the concatenated MLSA gene sequences (4,071 bp) also showed three main clusters, each comprising the type strain of one species. The M. abscessus cluster had a bootstrap value of 100% and was mostly compact. Bootstrap values for the M. massiliense and M. bolletii branches were much lower (71 and 61%, respectively), with the M. massiliense cluster having a fuzzy aspect. Mean (range) divergence values were 2.17% (1.13 to 2.58%) between the M. abscessus and M. massiliense clusters, 2.37% (1.5 to 2.85%) between the M. abscessus and M. bolletii clusters, and 2.28% (0.86 to 2.68%) between the M. massiliense and M. bolletii clusters. Adding the rpoB sequence to the MLSA-concatenated sequence (total sequence, 4,823 bp) had little effect on the clustering of strains. We found 10/120 (8.3%) isolates for which the concatenated MLSA gene sequence and rpoB sequence were discordant (e.g., M. massiliense MLSA sequence and M. abscessus rpoB sequence), suggesting the intergroup lateral transfers of rpoB. In conclusion, our study strongly supports the recent proposal that M. abscessus, M. massiliense, and M. bolletii should constitute a single species. Our findings also indicate that there has been a horizontal transfer of rpoB sequences between these subgroups, precluding the use of rpoB sequencing alone for the accurate identification of the two proposed M. abscessus subspecies.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.