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.

Stable masking by H-2Dd cis ligand limits Ly49A relocalization to the site of NK cell/target cell contact.

Ly49A is an inhibitory receptor, which counteracts natural killer (NK) cell activation on the engagement with H-2D(d) (D(d)) MHC class I molecules (MHC-I) on target cells. In addition to binding D(d) on apposed membranes, Ly49A interacts with D(d) ligand expressed in the plane of the NK cells' membrane. Indeed, multivalent, soluble MHC-I ligand binds inefficiently to Ly49A unless the NK cells' D(d) complexes are destroyed. However, it is not known whether masked Ly49A remains constitutively associated with cis D(d) also during target cell interaction. Alternatively, it is possible that Ly49A has to be unmasked to significantly interact with its ligand on target cells. These two scenarios suggest distinct roles of Ly49A/D(d) cis interaction for NK cell function. Here, we show that Ly49A contributes to target cell adhesion and efficiently accumulates at synapses with D(d)-expressing target cells when NK cells themselves lack D(d). When NK cells express D(d), Ly49A no longer contributes to adhesion, and ligand-driven recruitment to the cellular contact site is strongly reduced. The destruction of D(d) complexes on NK cells, which unmasks Ly49A, is necessary and sufficient to restore Ly49A adhesive function and recruitment to the synapse. Thus, cis D(d) continuously sequesters a considerable fraction of Ly49A receptors, preventing efficient Ly49A recruitment to the synapse with D(d)+ target cells. The reduced number of Ly49A receptors that can functionally interact with D(d) on target cells explains the modest inhibitory capacity of Ly49A in D(d) NK cells. This property renders Ly49A NK cells more sensitive to react to diseased host cells.

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.

MHC Class II Transactivator Is an In Vivo Regulator of Osteoclast Differentiation and Bone Homeostasis Co-opted From Adaptive Immunity.

The molecular networks controlling bone homeostasis are not fully understood. The common evolution of bone and adaptive immunity encourages the investigation of shared regulatory circuits. MHC Class II Transactivator (CIITA) is a master transcriptional co-activator believed to be exclusively dedicated for antigen presentation. CIITA is expressed in osteoclast precursors, and its expression is accentuated in osteoporotic mice. We thus asked whether CIITA plays a role in bone biology. To this aim, we fully characterized the bone phenotype of two mouse models of CIITA overexpression, respectively systemic and restricted to the monocyte-osteoclast lineage. Both CIITA-overexpressing mouse models revealed severe spontaneous osteoporosis, as assessed by micro-computed tomography and histomorphometry, associated with increased osteoclast numbers and enhanced in vivo bone resorption, whereas osteoblast numbers and in vivo bone-forming activity were unaffected. To understand the underlying cellular and molecular bases, we investigated ex vivo the differentiation of mutant bone marrow monocytes into osteoclasts and immune effectors, as well as osteoclastogenic signaling pathways. CIITA-overexpressing monocytes differentiated normally into effector macrophages or dendritic cells but showed enhanced osteoclastogenesis, whereas CIITA ablation suppressed osteoclast differentiation. Increased c-fms and receptor activator of NF-κB (RANK) signaling underlay enhanced osteoclast differentiation from CIITA-overexpressing precursors. Moreover, by extending selected phenotypic and cellular analyses to additional genetic mouse models, namely MHC Class II deficient mice and a transgenic mouse line lacking a specific CIITA promoter and re-expressing CIITA in the thymus, we excluded MHC Class II expression and T cells from contributing to the observed skeletal phenotype. Altogether, our study provides compelling genetic evidence that CIITA, the molecular switch of antigen presentation, plays a novel, unexpected function in skeletal homeostasis, independent of MHC Class II expression and T cells, by exerting a selective and intrinsic control of osteoclast differentiation and bone resorption in vivo. © 2014 American Society for Bone and Mineral Research.

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.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

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.

Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

Functional diversity of bacterial genes associated with aromatic hydrocarbon degradation in anthropogenic dark earth of Amazonia

The objective of this work was to evaluate the catabolic gene diversity for the bacterial degradation of aromatic hydrocarbons in anthropogenic dark earth of Amazonia (ADE) and their biochar (BC). Functional diversity analyses in ADE soils can provide information on how adaptive microorganisms may influence the fertility of soils and what is their involvement in biogeochemical cycles. For this, clone libraries containing the gene encoding for the alpha subunit of aromatic ring-hydroxylating dioxygenases (α-ARHD bacterial gene) were constructed, totaling 800 clones. These libraries were prepared from samples of an ADE soil under two different land uses, located at the Caldeirão Experimental Station - secondary forest (SF) and agriculture (AG) -, and the biochar (SF_BC and AG_BC, respectively). Heterogeneity estimates indicated greater diversity in BC libraries; and Venn diagrams showed more unique operational protein clusters (OPC) in the SF_BC library than the ADE soil, which indicates that specific metabolic processes may occur in biochar. Phylogenetic analysis showed unidentified dioxygenases in ADE soils. Libraries containing functional gene encoding for the alpha subunit of the aromatic ring-hydroxylating dioxygenases (ARHD) gene from biochar show higher diversity indices than those of ADE under secondary forest and agriculture.