The fine-scale architecture of structural variants in 17 mouse genomes.

BACKGROUND: Accurate catalogs of structural variants (SVs) in mammalian genomes are necessary to elucidate the potential mechanisms that drive SV formation and to assess their functional impact. Next generation sequencing methods for SV detection are an advance on array-based methods, but are almost exclusively limited to four basic types: deletions, insertions, inversions and copy number gains. RESULTS: By visual inspection of 100 Mbp of genome to which next generation sequence data from 17 inbred mouse strains had been aligned, we identify and interpret 21 paired-end mapping patterns, which we validate by PCR. These paired-end mapping patterns reveal a greater diversity and complexity in SVs than previously recognized. In addition, Sanger-based sequence analysis of 4,176 breakpoints at 261 SV sites reveal additional complexity at approximately a quarter of structural variants analyzed. We find micro-deletions and micro-insertions at SV breakpoints, ranging from 1 to 107 bp, and SNPs that extend breakpoint micro-homology and may catalyze SV formation. CONCLUSIONS: An integrative approach using experimental analyses to train computational SV calling is essential for the accurate resolution of the architecture of SVs. We find considerable complexity in SV formation; about a quarter of SVs in the mouse are composed of a complex mixture of deletion, insertion, inversion and copy number gain. Computational methods can be adapted to identify most paired-end mapping patterns.

Impact of normalization on miRNA microarray expression profiling.

Profiling miRNA levels in cells with miRNA microarrays is becoming a widely used technique. Although normalization methods for mRNA gene expression arrays are well established, miRNA array normalization has so far not been investigated in detail. In this study we investigate the impact of normalization on data generated with the Agilent miRNA array platform. We have developed a method to select nonchanging miRNAs (invariants) and use them to compute linear regression normalization coefficients or variance stabilizing normalization (VSN) parameters. We compared the invariants normalization to normalization by scaling, quantile, and VSN with default parameters as well as to no normalization using samples with strong differential expression of miRNAs (heart-brain comparison) and samples where only a few miRNAs are affected (by p53 overexpression in squamous carcinoma cells versus control). All normalization methods performed better than no normalization. Normalization procedures based on the set of invariants and quantile were the most robust over all experimental conditions tested. Our method of invariant selection and normalization is not limited to Agilent miRNA arrays and can be applied to other data sets including those from one color miRNA microarray platforms, focused gene expression arrays, and gene expression analysis using quantitative PCR.

Mapping the molecular characteristics of Brazilian human T-cell lymphotropic virus type 1 Env (gp46) and Pol amino acid sequences for vaccine design

This study was carried out to evaluate the molecular pattern of all available Brazilian human T-cell lymphotropic virus type 1 Env (n = 15) and Pol (n = 43) nucleotide sequences via epitope prediction, physico-chemical analysis, and protein potential sites identification, giving support to the Brazilian AIDS vaccine program. In 12 previously described peptides of the Env sequences we found 12 epitopes, while in 4 peptides of the Pol sequences we found 4 epitopes. The total variation on the amino acid composition was 9 and 17% for human leukocyte antigen (HLA) class I and class II Env epitopes, respectively. After analyzing the Pol sequences, results revealed a total amino acid variation of 0.75% for HLA-I and HLA-II epitopes. In 5 of the 12 Env epitopes the physico-chemical analysis demonstrated that the mutations magnified the antigenicity profile. The potential protein domain analysis of Env sequences showed the loss of a CK-2 phosphorylation site caused by D197N mutation in one epitope, and a N-glycosylation site caused by S246Y and V247I mutations in another epitope. Besides, the analysis of selection pressure have found 8 positive selected sites (w = 9.59) using the codon-based substitution models and maximum-likelihood methods. These studies underscore the importance of this Env region for the virus fitness, for the host immune response and, therefore, for the development of vaccine candidates.

A novel study of copy number variations in Hirschsprung disease using the multiple ligation-dependent probe amplification (MLPA) technique.

BACKGROUND Hirschsprung disease (HSCR) is a congenital malformation of the hindgut produced by a disruption in neural crest cell migration during embryonic development. HSCR has a complex genetic etiology and mutations in several genes, mainly the RET proto-oncogene, have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. METHODS In this study we have aimed to analyze the presence of CNVs in some HSCR genes (RET, EDN3, GDNF and ZFHX1B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. RESULTS Two alterations in the MLPA profiles of RET and EDN3 were detected, but a detailed inspection showed that the decrease in the corresponding dosages were due to point mutations affecting the hybridization probes regions. CONCLUSION Our results indicate that CNVs of the gene coding regions analyzed here are not a common molecular cause of Hirschsprung disease. However, further studies are required to determine the presence of CNVs affecting non-coding regulatory regions, as well as other candidate genes.

Altered expression of beta-galactosidase-1-like protein 3 (Glb1l3) in the retinal pigment epithelium (RPE)-specific 65-kDa protein knockout mouse model of Leber's congenital amaurosis

Purpose: In this study, we investigated the expression of the gene encoding beta-galactosidase (Glb)-1-like protein 3 (Glb1l3), a member of the glycosyl hydrolase 35 family, during retinal degeneration in the retinal pigment epithelium (RPE)-specific 65-kDa protein knockout (Rpe65(-/-)) mouse model of Leber congenital amaurosis (LCA). Additionally, we assessed the expression of the other members of this protein family, including beta-galactosidase-1 (Glb1), beta-galactosidase-1-like (Glb1l), and beta-galactosidase-1-like protein 2 (Glb1l2).Methods: The structural features of Glb1l3 were assessed using bioinformatic tools. mRNA expression of Glb-related genes was investigated by oligonucleotide microarray, real-time PCR, and reverse transcription (RT) -PCR. The localized expression of Glb1l3 was assessed by combined in situ hybridization and immunohistochemistry.Results: Glb1l3 was the only Glb-related member strongly downregulated in Rpe65(-/-) retinas before the onset and during progression of the disease. Glb1l3 mRNA was only expressed in the retinal layers and the RPE/choroid. The other Glb-related genes were ubiquitously expressed in different ocular tissues, including the cornea and lens. In the healthy retina, expression of Glb1l3 was strongly induced during postnatal retinal development; age-related increased expression persisted during adulthood and aging.Conclusions: These data highlight early-onset downregulation of Glb1l3 in Rpe65-related disease. They further indicate that impaired expression of Glb1l3 is mostly due to the absence of the chromophore 11-cis retinal, suggesting that Rpe65 deficiency may have many metabolic consequences in the underlying neuroretina.

Single cell analysis reveals similar functional competence of dominant and nondominant CD8 T-cell clonotypes.

Immune protection from infectious diseases and cancer is mediated by individual T cells of different clonal origin. Their functions are tightly regulated but not yet fully characterized. Understanding the contribution of each T cell will improve the prediction of immune protection based on laboratory assessment of T-cell responses. Here we developed techniques for simultaneous molecular and functional assessment of single CD8 T cells directly ex vivo. We studied two groups of patients with melanoma after vaccination with two closely related tumor antigenic peptides. Vaccination induced T cells with strong memory and effector functions, as found in virtually all T cells of the first patient group, and fractions of T cells in the second group. Interestingly, high functionality was not restricted to dominant clonotypes. Rather, dominant and nondominant clonotypes acquired equal functional competence. In parallel, this was also found for EBV- and CMV-specific T cells. Thus, the nondominant clonotypes may contribute similarly to immunity as their dominant counterparts.

Effectiveness of holistic neuropsychological rehabilitation for Spanish population with acquired brain injury measured using Rasch analysis.

INTRODUCTION The Rasch model is increasingly used in the field of rehabilitation because it improves the accuracy of measurements of patient status and their changes after therapy. OBJECTIVE To determine the long-term effectiveness of a holistic neuropsychological rehabilitation program for Spanish outpatients with acquired brain injury (ABI) using Rasch analysis. METHODS Eighteen patients (ten with long evolution - patients who started the program > 6 months after ABI- and eight with short evolution) and their relatives attended the program for 6 months. Patients' and relatives' answers to the European Brain Injury Questionnaire and the Frontal Systems Behavior Scale at 3 time points (pre-intervention. post-intervention and 12 month follow-up) were transformed into linear measures called logits. RESULTS The linear measures revealed significant improvements with large effects at the follow-up assessment on cognitive and executive functioning, social and emotional self-regulation, apathy and mood. At follow-up, the short evolution group achieved greater improvements in mood and cognitive functioning than the long evolution patients. CONCLUSIONS The program showed long-term effectiveness for most of the variables, and it was more effective for mood and cognitive functioning when patients were treated early. Relatives played a key role in the effectiveness of the rehabilitation program.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

Risk assessment of recurrence in sporadic retinoblastoma using a molecular-based algorithm.

PURPOSE: Most RB1 mutations are unique and distributed throughout the RB1 gene. Their detection can be time-consuming and the yield especially low in cases of conservatively-treated sporadic unilateral retinoblastoma (Rb) patients. In order to identify patients with true risk of developing Rb, and to reduce the number of unnecessary examinations under anesthesia in all other cases, we developed a universal sensitive, efficient and cost-effective strategy based on intragenic haplotype analysis. METHODS: This algorithm allows the calculation of the a posteriori risk of developing Rb and takes into account (a) RB1 loss of heterozygosity in tumors, (b) preferential paternal origin of new germline mutations, (c) a priori risk derived from empirical data by Vogel, and (d) disease penetrance of 90% in most cases. We report the occurrence of Rb in first degree relatives of patients with sporadic Rb who visited the Jules Gonin Eye Hospital, Lausanne, Switzerland, from January 1994 to December 2006 compared to expected new cases of Rb using our algorithm. RESULTS: A total of 134 families with sporadic Rb were enrolled; testing was performed in 570 individuals and 99 patients younger than 4 years old were identified. We observed one new case of Rb. Using our algorithm, the cumulated total a posteriori risk of recurrence was 1.77. CONCLUSIONS: This is the first time that linkage analysis has been validated to monitor the risk of recurrence in sporadic Rb. This should be a useful tool in genetic counseling, especially when direct RB1 screening for mutations leaves a negative result or is unavailable.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

Essential role of the N-terminal region of TFII-I in viability and behavior

Background: GTF2I codes for a general intrinsic transcription factor and calcium channel regulator TFII-I, with high and ubiquitous expression, and a strong candidate for involvement in the morphological and neuro-developmental anomalies of the Williams-Beuren syndrome (WBS). WBS is a genetic disorder due to a recurring deletion of about 1,55-1,83 Mb containing 25-28 genes in chromosome band 7q11.23 including GTF2I. Completed homozygous loss of either the Gtf2i or Gtf2ird1 function in mice provided additional evidence for the involvement of both genes in the craniofacial and cognitive phenotype. Unfortunately nothing is now about the behavioral characterization of heterozygous mice. Methods: By gene targeting we have generated a mutant mice with a deletion of the first 140 amino-acids of TFII-I. mRNA and protein expression analysis were used to document the effect of the study deletion. We performed behavioral characterization of heterozygous mutant mice to document in vivo implications of TFII-I in the cognitive profile of WBS patients. Results: Homozygous and heterozygous mutant mice exhibit craniofacial alterations, most clearly represented in homozygous condition. Behavioral test demonstrate that heterozygous mutant mice exhibit some neurobehavioral alterations and hyperacusis or odynacusis that could be associated with specific features of WBS phenotype. Homozygous mutant mice present highly compromised embryonic viability and fertility. Regarding cellular model, we documented a retarded growth in heterozygous MEFs respect to homozygous or wild-type MEFs. Conclusion: Our data confirm that, although additive effects of haploinsufficiency at several genes may contribute to the full craniofacial or neurocognitive features of WBS, correct expression of GTF2I is one of the main players. In addition, these findings show that the deletion of the fist 140 amino-acids of TFII-I altered it correct function leading to a clear phenotype, at both levels, at the cellular model and at the in vivo model.

The peroxisome proliferator-activated receptor alpha regulates amino acid metabolism.

The peroxisome proliferator-activated receptor alpha is a ligand-activated transcription factor that plays an important role in the regulation of lipid homeostasis. PPARalpha mediates the effects of fibrates, which are potent hypolipidemic drugs, on gene expression. To better understand the biological effects of fibrates and PPARalpha, we searched for genes regulated by PPARalpha using oligonucleotide microarray and subtractive hybridization. By comparing liver RNA from wild-type and PPARalpha null mice, it was found that PPARalpha decreases the mRNA expression of enzymes involved in the metabolism of amino acids. Further analysis by Northern blot revealed that PPARalpha influences the expression of several genes involved in trans- and deamination of amino acids, and urea synthesis. Direct activation of PPARalpha using the synthetic PPARalpha ligand WY14643 decreased mRNA levels of these genes, suggesting that PPARalpha is directly implicated in the regulation of their expression. Consistent with these data, plasma urea concentrations are modulated by PPARalpha in vivo. It is concluded that in addition to oxidation of fatty acids, PPARalpha also regulates metabolism of amino acids in liver, indicating that PPARalpha is a key controller of intermediary metabolism during fasting.

Hypertension prevalence among indigenous populations in Brazil: a systematic review with meta-analysis

Abstract OBJECTIVE Evaluating the evidence of hypertension prevalence among indigenous populations in Brazil through a systematic review and meta-analysis. METHODS A search was performed by two reviewers, with no restriction of date or language in the databases of PubMed, LILACS, SciELO, Virtual Health Library and Capes Journal Portal. Also, a meta-regression model was designed in which the last collection year of each study was used as a moderating variable. RESULTS 23 articles were included in the review. No hypertension was found in indigenous populations in 10 studies, and its prevalence was increasing and varied, reaching levels of up to 29.7%. Combined hypertension prevalence in Indigenous from the period of 1970 to 2014 was 6.2% (95% CI, 3.1% - 10.3%). In the regression, the value of the odds ratio was 1.12 (95% CI, 1.07 - 1.18; p <0.0001), indicating a 12% increase every year in the probability of an indigenous person presenting hypertension. CONCLUSION There has been a constant increase in prevalence despite the absence of hypertension in about half of the studies, probably due to changes in cultural, economic and lifestyle habits, resulting from indigenous interaction with non-indigenous society.