Apparent versus true gene expression changes of three hypoxia-related genes in autopsy derived tissue and the importance of normalisation.

The aim of our work was to show how a chosen normal-isation strategy can affect the outcome of quantitative gene expression studies. As an example, we analysed the expression of three genes known to be upregulated under hypoxic conditions: HIF1A, VEGF and SLC2A1 (GLUT1). Raw RT-qPCR data were normalised using two different strategies: a straightforward normalisation against a single reference gene, GAPDH, using the 2(-ΔΔCt) algorithm and a more complex normalisation against a normalisation factor calculated from the quantitative raw data from four previously validated reference genes. We found that the two different normalisation strategies revealed contradicting results: normalising against a validated set of reference genes revealed an upregulation of the three genes of interest in three post-mortem tissue samples (cardiac muscle, skeletal muscle and brain) under hypoxic conditions. Interestingly, we found a statistically significant difference in the relative transcript abundance of VEGF in cardiac muscle between donors who died of asphyxia versus donors who died from cardiac death. Normalisation against GAPDH alone revealed no upregulation but, in some instances, a downregulation of the genes of interest. To further analyse this discrepancy, the stability of all reference genes used were reassessed and the very low expression stability of GAPDH was found to originate from the co-regulation of this gene under hypoxic conditions. We concluded that GAPDH is not a suitable reference gene for the quantitative analysis of gene expression in hypoxia and that validation of reference genes is a crucial step for generating biologically meaningful data.

A minimal i-motif stabilized by minor groove G:T:G:T tetrads

The repetitive DNA sequences found at telomeres and centromeres play a crucial role in the structure and function of eukaryotic chromosomes. This role may be related to the tendency observed in many repetitive DNAs to adopt non-canonical structures. Although there is an increasing recognition of the importance of DNA quadruplexes in chromosome biology, the co-existence of different quadruplex-forming elements in the same DNA structure is still a matter of debate. Here we report the structural study of the oligonucleotide d(TCGTTTCGT) and its cyclic analog d. Both sequences form dimeric quadruplex structures consisting of a minimal i-motif capped, at both ends, by a slipped minor groove-aligned G:T:G:T tetrad. These mini i-motifs, which do not exhibit the characteristic CD spectra of other i-motif structures, can be observed at neutral pH, although they are more stable under acidic conditions. This finding is particularly relevant since these oligonucleotide sequences do not contain contiguous cytosines. Importantly, these structures resemble the loop moiety adopted by an 11-nucleotide fragment of the conserved centromeric protein B (CENP-B) box motif, which is the binding site for the CENP-B.

A minimal i-motif stabilized by minor groove G:T:G:T tetrads

The repetitive DNA sequences found at telomeres and centromeres play a crucial role in the structure and function of eukaryotic chromosomes. This role may be related to the tendency observed in many repetitive DNAs to adopt non-canonical structures. Although there is an increasing recognition of the importance of DNA quadruplexes in chromosome biology, the co-existence of different quadruplex-forming elements in the same DNA structure is still a matter of debate. Here we report the structural study of the oligonucleotide d(TCGTTTCGT) and its cyclic analog d. Both sequences form dimeric quadruplex structures consisting of a minimal i-motif capped, at both ends, by a slipped minor groove-aligned G:T:G:T tetrad. These mini i-motifs, which do not exhibit the characteristic CD spectra of other i-motif structures, can be observed at neutral pH, although they are more stable under acidic conditions. This finding is particularly relevant since these oligonucleotide sequences do not contain contiguous cytosines. Importantly, these structures resemble the loop moiety adopted by an 11-nucleotide fragment of the conserved centromeric protein B (CENP-B) box motif, which is the binding site for the CENP-B.

Architectural features of bacterial genomes and their applications

Résumé -Caractéristiques architecturales des génomes bactériens et leurs applications Les bactéries possèdent généralement un seul chromosome circulaire. A chaque génération, ce chromosome est répliqué bidirectionnellement, par deux complexes enzymatiques de réplication se déplaçant en sens opposé depuis l'origine de réplication jusqu'au terminus, situé à l'opposé. Ce mode de réplication régit l'architecture du chromosome -l'orientation des gènes par rapport à la réplication, notamment - et est en grande partie à l'origine des pressions qui provoquent la variation de la composition en nucléotides du génome, hors des contraintes liées à la structure et à la fonction des protéines codées sur le chromosome. Le but de cette thèse est de contribuer à quantifier les effets de la réplication sur l'architecture chromosomique, en s'intéressant notamment aux gènes des ARN ribosomiques, cruciaux pour la bactérie. D'un autre côté, cette architecture est spécifique à l'espèce et donne ainsi une «identité génomique » aux gènes. Il est démontré ici qu'il est possible d'utiliser des marqueurs «naïfs » de cette identité pour détecter, notamment dans le génome du staphylocoque doré, des îlots de pathogénicité, qui concentrent un grand nombre de facteurs de virulence de la bactérie. Ces îlots de pathogénicité sont mobiles, et peuvent passer d'une bactérie à une autre, mais conservent durant un certain temps l'identité génomique de leur hôte précédent, ce qui permet de les reconnaître dans leur nouvel hôte. Ces méthodes simples, rapides et fiables seront de la plus haute importance lorsque le séquençage des génomes entiers sera rapide et disponible à très faible coût. Il sera alors possible d'analyser instantanément les déterminants pathogéniques et de résistance aux antibiotiques des agents pathogènes. Summary The bacterial genome is a highly organized structure, which may be referred to as the genome architecture, and is mainly directed by DNA replication. This thesis provides significant insights in the comprehension of the forces that shape bacterial chromosomes, different in each genome and contributing to confer them an identity. First, it shows the importance of the replication in directing the orientation of prokaryotic ribosomal RNAs, and how it shapes their nucleotide composition in a tax on-specific manner. Second, it highlights the pressure acting on the orientation of the genes in general, a majority of which are transcribed in the same direction as replication. Consequently, apparent infra-arm genome rearrangements, involving an exchange of the leading/lagging strands and shown to reduce growth rate, are very likely artifacts due to an incorrect contig assembly. Third, it shows that this genomic identity can be used to detect foreign parts in genomes, by establishing this identity for a given host and identifying the regions that deviate from it. This property is notably illustrated with Staphylococcus aureus: known pathogenicity islands and phages, and putative ancient pathogenicity islands concentrating many known pathogenicity-related genes are highlighted; the analysis also detects, incidentally, proteins responsible for the adhesion of S. aureus to the hosts' cells. In conclusion, the study of nucleotide composition of bacterial genomes provides the opportunity to better understand the genome-level pressures that shape DNA sequences, and to identify genes and regions potentially related to pathogenicity with fast, simple and reliable methods. This will be of crucial importance when whole-genome sequencing will be a rapid, inexpensive and routine tool.

IFN Stimulated Gene Expression in the Liver is a Better Predictor of Treatment Response in Chronic Hepatitis C than the IL28B (IFN lambda 3) Genotype

BACKGROUND: Therapy of chronic hepatitis C (CHC) with pegIFNα/ribavirin achieves a sustained virologic response (SVR) in ∼55%. Pre-activation of the endogenous interferon system in the liver is associated with non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. METHODS: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. RESULTS: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. CONCLUSIONS: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

Hyperammonemia: regulation of argininosuccinate synthetase and argininosuccinate lyase genes in aggregating cell cultures of fetal rat brain.

Hyperammonemia in the brain leads to poorly understood alterations of nitric oxide (NO) synthesis. Arginine, the substrate of nitric oxide synthases, might be recycled from the citrulline produced with NO by argininosuccinate synthetase (AS) and argininosuccinate lyase (AL). The regulation of AS and AL genes during hyperammonemia is unknown in the brain. We used brain cell aggregates cultured from dissociated telencephalic cortex of rat embryos to analyze the regulation of AS and AL genes in hyperammonemia. Using RNase protection assay and non-radioactive in situ hybridization on aggregate cryosections, we show that both AS and AL genes are induced in astrocytes but not in neurons of aggregates exposed to 5 mM NH4Cl. Our work suggests that the hyperammonemic brain might increase its recycling of citrulline to arginine.

Evolutionary footprint of coevolving positions in genes.

MOTIVATION: The analysis of molecular coevolution provides information on the potential functional and structural implication of positions along DNA sequences, and several methods are available to identify coevolving positions using probabilistic or combinatorial approaches. The specific nucleotide or amino acid profile associated with the coevolution process is, however, not estimated, but only known profiles, such as the Watson-Crick constraint, are usually considered a priori in current measures of coevolution. RESULTS: Here, we propose a new probabilistic model, Coev, to identify coevolving positions and their associated profile in DNA sequences while incorporating the underlying phylogenetic relationships. The process of coevolution is modeled by a 16 × 16 instantaneous rate matrix that includes rates of transition as well as a profile of coevolution. We used simulated, empirical and illustrative data to evaluate our model and to compare it with a model of 'independent' evolution using Akaike Information Criterion. We showed that the Coev model is able to discriminate between coevolving and non-coevolving positions and provides better specificity and specificity than other available approaches. We further demonstrate that the identification of the profile of coevolution can shed new light on the process of dependent substitution during lineage evolution.

Diversidad genética en cerdos criollos mexicanos con genes candidatos asociados a características productivas

El objetivo de este trabajo fue comparar la variabilidad, diversidad y distancias genéticas entre cerdos criollos, Pelón Mexicano (CPM) y Cuinos (CC), con Yorkshire, cuanto a los genes candidatos CAST, DECR1, HAL, HFABP4, LEP, LIPE, MCR4, MYOG, RN y CHX, a través de analysis por PCR-RFLP. Se evaluaron 180 cerdos: 59 CPM, 65 CC y 56 Yorkshire. Se analizaron las frecuencias génicas y genotípicas, heterocigosidad, distancias genéticas y árboles filogenéticos entre grupos raciales. Para CAST, DECR1, HFABP4, LEP, MCR4 y CHX las frecuencias génicas y genotípicas fueron diferentes al comparar las tres razas. En LIPE, los CC fueron iguales a los Yorkshire; en cuanto a MYOG, los CPM fueron iguales a los Yorkshire. No hubo diferencias entre poblaciones criollas y Yorkshire en las frecuencias génicas y genotípicas para HAL y RN. Los cerdos Yorkshire presentaron mayor frecuencia en alelos favorables para CAST, LIPE, MCR4 y MYOG, menor frecuencia de DECR1, HFABP4, CHX, y moderada en LEP. La heterocigosidad promedio para todos los genes fue mayor en CPM (0,42±0,05) y similar en CC (0,33±0,06) y Yorkshire (0,35±0,05). Al calcular distancias genéticas con todos los genes, los CC se encuentran más distantes de los Yorkshire.

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

The objective of this work was to validate, by quantitative PCR in real time (RT-qPCR), genes to be used as reference in studies of gene expression in soybean in drought-stressed trials. Four genes commonly used in soybean were evaluated: Gmβ-actin, GmGAPDH, GmLectin and GmRNAr18S. Total RNA was extracted from six samples: three from roots in a hydroponic system with different drought intensities (0, 25, 50, 75 and 100 minutes of water stress), and three from leaves of plants grown in sand with different soil moistures (15, 5 and 2.5% gravimetric humidity). The raw cycle threshold (Ct) data were analyzed, and the efficiency of each primer was calculated for an overall analysis of the Ct range among the different samples. The GeNorm application was used to evaluate the best reference gene, according to its stability. The GmGAPDH was the least stable gene, with the highest mean values of expression stability (M), and the most stable genes, with the lowest M values, were the Gmβ-actin and GmRNAr18S, when both root and leaves samples were tested. These genes can be used in RT-qPCR as reference gene for expression analysis.

Polymorphisms in DNA-repair genes in a cohort of prostate cancer patients from different areas in Spain: heterogeneity between populations as a confounding factor in association studies

Background: Differences in the distribution of genotypes between individuals of the same ethnicity are an important confounder factor commonly undervalued in typical association studies conducted in radiogenomics. Objective: To evaluate the genotypic distribution of SNPs in a wide set of Spanish prostate cancer patients for determine the homogeneity of the population and to disclose potential bias. Design, Setting, and Participants: A total of 601 prostate cancer patients from Andalusia, Basque Country, Canary and Catalonia were genotyped for 10 SNPs located in 6 different genes associated to DNA repair: XRCC1 (rs25487, rs25489, rs1799782), ERCC2 (rs13181), ERCC1 (rs11615), LIG4 (rs1805388, rs1805386), ATM (rs17503908, rs1800057) and P53 (rs1042522). The SNP genotyping was made in a Biotrove OpenArrayH NT Cycler. Outcome Measurements and Statistical Analysis: Comparisons of genotypic and allelic frequencies among populations, as well as haplotype analyses were determined using the web-based environment SNPator. Principal component analysis was made using the SnpMatrix and XSnpMatrix classes and methods implemented as an R package. Non-supervised hierarchical cluster of SNP was made using MultiExperiment Viewer. Results and Limitations: We observed that genotype distribution of 4 out 10 SNPs was statistically different among the studied populations, showing the greatest differences between Andalusia and Catalonia. These observations were confirmed in cluster analysis, principal component analysis and in the differential distribution of haplotypes among the populations. Because tumor characteristics have not been taken into account, it is possible that some polymorphisms may influence tumor characteristics in the same way that it may pose a risk factor for other disease characteristics. Conclusion: Differences in distribution of genotypes within different populations of the same ethnicity could be an important confounding factor responsible for the lack of validation of SNPs associated with radiation-induced toxicity, especially when extensive meta-analysis with subjects from different countries are carried out.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Identificação dos genes Ty-2 e Ty-3 de resistência a begomovírusem genótipos de tomateiro

O objetivo deste trabalho foi identificar a presença dos genes Ty-2 e Ty-3, de resistência a begomovírus, em acessos de tomateiro do Banco de Germoplasma de Hortaliças da Universidade Federal de Viçosa. Os oligonucleotídeos TO302 F/R e FLUW25 F/R foram utilizados em reações de PCR, para verificar a presença de marcadores relacionados aos genes Ty-2 e Ty-3, respectivamente. Observou-se a presença do gene Ty-2, em heterozigose na subamostra BGH-6881 (Solanum peruvianum), e do gene Ty-3, em homozigose nas subamostras BGH-6878, BGH-6897 (S. lycopersicum) e em heterozigose na subamostra BGH-6881. A identificação dos genes de resistência, com reações de PCR, representa um avanço para os programas de melhoramento de tomateiro no Brasil.

Expressão de genes da subfamília HD-Zip I em soja submetida à seca

O objetivo deste trabalho foi identificar genes candidatos da subfamília de fatores transcricionais HD-Zip I que contribuem para a tolerância à seca em soja. Foram avaliados trifólios de soja de cultivar tolerante (Embrapa 48) e suscetível à seca (BR 16), sob três níveis de deficit hídrico: ausência, moderado (-1,5 MPa) e severo (-3,0 MPa). Pela análise dos promotores, foi identificada a presença de possíveis elementos cis-regulatórios relacionados à resposta à seca, nos três genes avaliados (GmHB6, GmHB13 e GmHB21). No entanto, não houve padrão de distribuição específico associado à maior tolerância do genótipo à seca. Com a análise comparativa, foram identificados seis elementos cis-regulatórios potencialmente envolvidos na indução da expressão gênica sob seca. O gene GmHB13 foi exclusivamente induzido pela seca no genótipo tolerante, e o gene GmHB6 apresentou redução da expressão somente no genótipo suscetível. Já o gene GmHB21, apresentou aumento da expressão em ambos os genótipos. O gene GmHB13 é um importante elemento na regulação do mecanismo de tolerância à seca em soja, na cultivar tolerante Embrapa 48.

Identificação e caracterização de genes vip e cry coleóptero‑específicos em isolados de Bacillus thuringiensis

O objetivo deste trabalho foi identificar e caracterizar os genes cry3, vip1, vip2 e vip1/vip2 em uma coleção de 1.078 isolados de Bacillus thuringiensis potencialmente tóxicos para larvas de coleópteros. Foram utilizados pares de oligonucleotídeos iniciadores gerais obtidos a partir de regiões conservadas dos genes e do alinhamento de sequências consenso. Posteriormente, os isolados positivos foram caracterizados por meio da técnica de PCR‑RFLP, tendo-se utilizado enzimas de restrição específicas, para identificar novas subclasses de genes nos isolados. Cento e cinquenta e um isolados foram positivos para os genes avaliados, com maior frequência para o gene vip1/vip2 (139 isolados). Pela técnica de PCR‑RFLP, foram observados 14 perfis polimórficos, o que indica a presença de diferentes alelos e, consequentemente, de distintas subclasses desses genes.

Minor protease inhibitor mutations at baseline do not increase the risk for a virological failure in HIV-1 subtype B infected patients.

BACKGROUND: Minor protease inhibitor (PI) mutations often exist as polymorphisms in HIV-1 sequences from treatment-naïve patients. Previous studies showed that their presence impairs the antiretroviral treatment (ART) response. Evaluating these findings in a larger cohort is essential. METHODS: To study the impact of minor PI mutations on time to viral suppression and time to virological failure, we included patients from the Swiss HIV Cohort Study infected with HIV-1 subtype B who started first-line ART with a PI and two nucleoside reverse transcriptase inhibitors. Cox regression models were performed to compare the outcomes among patients with 0 and ≥ 1 minor PI mutation. Models were adjusted for baseline HIV-1 RNA, CD4 cell count, sex, transmission category, age, ethnicity, year of ART start, the presence of nucleoside reverse transcriptase inhibitor mutations, and stratified for the administered PIs. RESULTS: We included 1199 patients of whom 944 (78.7%) received a boosted PI. Minor PI mutations associated with the administered PI were common: 41.7%, 16.1%, 4.7% and 1.9% had 1, 2, 3 or ≥ 4 mutations, respectively. The time to viral suppression was similar between patients with 0 (reference) and ≥ 1 minor PI mutation (multivariable hazard ratio (HR): 1.1 [95% confidence interval (CI): 1.0-1.3], P = .196). The time to virological failure was also similar (multivariable HR:.9 [95% CI:.5-1.6], P = .765). In addition, the impact of each single minor PI mutation was analyzed separately: none was significantly associated with the treatment outcome. CONCLUSIONS: The presence of minor PI mutations at baseline has no effect on the therapy outcome in HIV infected individuals.