Whole blood genome-wide gene expression profile in males after prolonged wakefulness and sleep recovery

Pellegrino R, Sunaga DY, Guindalini C, Martins RC, Mazzotti DR, Wei Z, Daye ZJ, Andersen ML, Tufik S. Whole blood genome-wide gene expression profile in males after prolonged wakefulness and sleep recovery. Physiol Genomics 44: 1003-1012, 2012. First published September 4, 2012; doi: 10.1152/physiolgenomics.00058.2012.-Although the specific functions of sleep have not been completely elucidated, the literature has suggested that sleep is essential for proper homeostasis. Sleep loss is associated with changes in behavioral, neurochemical, cellular, and metabolic function as well as impaired immune response. Using high-resolution microarrays we evaluated the gene expression profiles of healthy male volunteers who underwent 60 h of prolonged wakefulness (PW) followed by 12 h of sleep recovery (SR). Peripheral whole blood was collected at 8 am in the morning before the initiation of PW (Baseline), after the second night of PW, and one night after SR. We identified over 500 genes that were differentially expressed. Notably, these genes were related to DNA damage and repair and stress response, as well as diverse immune system responses, such as natural killer pathways including killer cell lectin-like receptors family, as well as granzymes and T-cell receptors, which play important roles in host defense. These results support the idea that sleep loss can lead to alterations in molecular processes that result in perturbation of cellular immunity, induction of inflammatory responses, and homeostatic imbalance. Moreover, expression of multiple genes was downregulated following PW and upregulated after SR compared with PW, suggesting an attempt of the body to re-establish internal homeostasis. In silico validation of alterations in the expression of CETN3, DNAJC, and CEACAM genes confirmed previous findings related to the molecular effects of sleep deprivation. Thus, the present findings confirm that the effects of sleep loss are not restricted to the brain and can occur intensely in peripheral tissues.

Genome-Wide Analysis in Brazilian Xavante Indians Reveals Low Degree of Admixture

Characterization of population genetic variation and structure can be used as tools for research in human genetics and population isolates are of great interest. The aim of the present study was to characterize the genetic structure of Xavante Indians and compare it with other populations. The Xavante, an indigenous population living in Brazilian Central Plateau, is one of the largest native groups in Brazil. A subset of 53 unrelated subjects was selected from the initial sample of 300 Xavante Indians. Using 86,197 markers, Xavante were compared with all populations of HapMap Phase III and HGDP-CEPH projects and with a Southeast Brazilian population sample to establish its population structure. Principal Components Analysis showed that the Xavante Indians are concentrated in the Amerindian axis near other populations of known Amerindian ancestry such as Karitiana, Pima, Surui and Maya and a low degree of genetic admixture was observed. This is consistent with the historical records of bottlenecks experience and cultural isolation. By calculating pair-wise F-st statistics we characterized the genetic differentiation between Xavante Indians and representative populations of the HapMap and from HGDP-CEPH project. We found that the genetic differentiation between Xavante Indians and populations of Ameridian, Asian, European, and African ancestry increased progressively. Our results indicate that the Xavante is a population that remained genetically isolated over the past decades and can offer advantages for genome-wide mapping studies of inherited disorders.

Genome wide scan for quantitative trait loci affecting tick resistance in cattle (Bos taurus × Bos indicus)

Abstract Background In tropical countries, losses caused by bovine tick Rhipicephalus (Boophilus) microplus infestation have a tremendous economic impact on cattle production systems. Genetic variation between Bos taurus and Bos indicus to tick resistance and molecular biology tools might allow for the identification of molecular markers linked to resistance traits that could be used as an auxiliary tool in selection programs. The objective of this work was to identify QTL associated with tick resistance/susceptibility in a bovine F2 population derived from the Gyr (Bos indicus) × Holstein (Bos taurus) cross. Results Through a whole genome scan with microsatellite markers, we were able to map six genomic regions associated with bovine tick resistance. For most QTL, we have found that depending on the tick evaluation season (dry and rainy) different sets of genes could be involved in the resistance mechanism. We identified dry season specific QTL on BTA 2 and 10, rainy season specific QTL on BTA 5, 11 and 27. We also found a highly significant genome wide QTL for both dry and rainy seasons in the central region of BTA 23. Conclusions The experimental F2 population derived from Gyr × Holstein cross successfully allowed the identification of six highly significant QTL associated with tick resistance in cattle. QTL located on BTA 23 might be related with the bovine histocompatibility complex. Further investigation of these QTL will help to isolate candidate genes involved with tick resistance in cattle.

Genome-wide metabolic (re-) annotation of Kluyveromyces lactis

Background: Even before having its genome sequence published in 2004, Kluyveromyces lactis had long been considered a model organism for studies in genetics and physiology. Research on Kluyveromyces lactis is quite advanced and this yeast species is one of the few with which it is possible to perform formal genetic analysis. Nevertheless, until now, no complete metabolic functional annotation has been performed to the proteins encoded in the Kluyveromyces lactis genome. Results: In this work, a new metabolic genome-wide functional re-annotation of the proteins encoded in the Kluyveromyces lactis genome was performed, resulting in the annotation of 1759 genes with metabolic functions, and the development of a methodology supported by merlin (software developed in-house). The new annotation includes novelties, such as the assignment of transporter superfamily numbers to genes identified as transporter proteins. Thus, the genes annotated with metabolic functions could be exclusively enzymatic (1410 genes), transporter proteins encoding genes (301 genes) or have both metabolic activities (48 genes). The new annotation produced by this work largely surpassed the Kluyveromyces lactis currently available annotations. A comparison with KEGG’s annotation revealed a match with 844 (~90%) of the genes annotated by KEGG, while adding 850 new gene annotations. Moreover, there are 32 genes with annotations different from KEGG. Conclusions: The methodology developed throughout this work can be used to re-annotate any yeast or, with a little tweak of the reference organism, the proteins encoded in any sequenced genome. The new annotation provided by this study offers basic knowledge which might be useful for the scientific community working on this model yeast, because new functions have been identified for the so-called metabolic genes. Furthermore, it served as the basis for the reconstruction of a compartmentalized, genome-scale metabolic model of Kluyveromyces lactis, which is currently being finished.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.

Allelic variations in the vitamin D receptor gene, insulin secretion and parents' heights are independently associated with height in obese children and adolescents

Polymorphisms in the VDR gene were reported to be associated with variations in intrauterine and postnatal growth and with adult height, but also with other traits that are strongly correlated such as the BMI, insulin sensitivity, insulin secretion and hyperglycemia. Here, we assessed the impact of VDR polymorphisms on body height and its interactions with obesity- and glucose tolerance-related traits in obese children and adolescents. We studied 173 prepubertal (Tanner's stage 1) and 146 pubertal (Tanner's stages 2-5) obese children who were referred for a weight-loss program. Three single nucleotide polymorphisms were genotyped: rs1544410 (BsmI), rs7975232 (ApaI) and rs731236 (TaqI). BsmI and TaqI genotypes were significantly associated with height in pubertal children, but the associations did not reach statistical significance in prepubertal children. In stepwise regression analyses, the lean body mass, insulin secretion, BsmI or TaqI genotypes and the father's and the mother's height were independently and positively associated with height in pubertal children. These covariables accounted for 46% of the trait variance. The height of homozygous carriers of the minor allele of BsmI was 0.65 z-scores (4 cm) higher than the height of homozygous carriers of the major allele (P=.0006). Haplotype analyses confirmed the associations of the minor alleles of BsmI and TaqI with increased height. In conclusion, VDR genotypes were significantly associated with height in pubertal obese children. The associations were independent from the effects of confounding traits, such as the body fat mass, insulin secretion, insulin sensitivity and glucose tolerance. (C) 2012 Elsevier Inc. All rights reserved.

Structural synaptic elements are differentially regulated in superior temporal cortex of schizophrenia patients

Inaccurate wiring and synaptic pathology appear to be major hallmarks of schizophrenia. A variety of gene products involved in synaptic neurotransmission and receptor signaling are differentially expressed in brains of schizophrenia patients. However, synaptic pathology may also develop by improper expression of intra- and extra-cellular structural elements weakening synaptic stability. Therefore, we have investigated transcription of these elements in the left superior temporal gyrus of 10 schizophrenia patients and 10 healthy controls by genome-wide microarrays (Illumina). Fourteen up-regulated and 22 downregulated genes encoding structural elements were chosen from the lists of differentially regulated genes for further qRT-PCR analysis. Almost all genes confirmed by this method were downregulated. Their gene products belonged to vesicle-associated proteins, that is, synaptotagmin 6 and syntaxin 12, to cytoskeletal proteins, like myosin 6, pleckstrin, or to proteins of the extracellular matrix, such as collagens, or laminin C3. Our results underline the pivotal roles of structural genes that control formation and stabilization of pre- and post-synaptic elements or influence axon guidance in schizophrenia. The glial origin of collagen or laminin highlights the close interrelationship between neurons and glial cells in establishment and maintenance of synaptic strength and plasticity. It is hypothesized that abnormal expression of these and related genes has a major impact on the pathophysiology of schizophrenia.

Loss of Interleukin-10 Signaling and Infantile Inflammatory Bowel Disease: Implications for Diagnosis and Therapy

BACKGROUND & AIMS: Homozygous loss of function mutations in interleukin-10 (IL10) and interleukin-10 receptors (IL10R) cause severe infantile (very early onset) inflammatory bowel disease (IBD). Allogeneic hematopoietic stem cell transplantation (HSCT) was reported to induce sustained remission in 1 patient with IL-10R deficiency. We investigated heterogeneity among patients with very early onset IBD, its mechanisms, and the use of allogeneic HSCT to treat this disorder. METHODS: We analyzed 66 patients with early onset IBD (younger than 5 years of age) for mutations in the genes encoding IL-10, IL-10R1, and IL-10R2. IL-10R deficiency was confirmed by functional assays on patients' peripheral blood mononuclear cells (immunoblot and enzyme-linked immunosorbent assay analyses). We assessed the therapeutic effects of standardized allogeneic HSCT. RESULTS: Using a candidate gene sequencing approach, we identified 16 patients with IL-10 or IL-10R deficiency: 3 patients had mutations in IL-10, 5 had mutations in IL-10R1, and 8 had mutations in IL-10R2. Refractory colitis became manifest in all patients within the first 3 months of life and was associated with perianal disease (16 of 16 patients). Extraintestinal symptoms included folliculitis (11 of 16) and arthritis (4 of 16). Allogeneic HSCT was performed in 5 patients and induced sustained clinical remission with a median follow-up time of 2 years. In vitro experiments confirmed reconstitution of IL-10R-mediated signaling in all patients who received the transplant. CONCLUSIONS: We identified loss of function mutations in IL-10 and IL-10R in patients with very early onset IBD. These findings indicate that infantile IBD patients with perianal disease should be screened for IL-10 and IL-10R deficiency and that allogeneic HSCT can induce remission in those with IL-10R deficiency.

Lack of evidence to support the association of a single IL28B genotype SNP rs12979860 with the HTLV-1 clinical outcomes and proviral load

Background: The Interleukin 28B (IL28B) rs12979860 polymorphisms was recently reported to be associated with the human T-cell leukemia virus type 1 (HTLV-1) proviral load (PvL) and the development of the HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Methods: In an attempt to examine this hypothesis, we assessed the association of the rs12979860 genotypes with HTLV-1 PvL levels and clinical status in 112 unrelated Brazilian subjects (81 HTLV-1 asymptomatic carriers, 24 individuals with HAM/TSP and 7 with Adult T cell Leukemia/Lymphoma (ATLL)). Results: All 112 samples were successfully genotyped and their PvLs compared. Neither the homozygote TT nor the heterozygote CT mutations nor the combination genotypes (TT/CT) were associated with a greater PvL. We also observed no significant difference in allele distribution between asymptomatic carriers and patients with HTLV-1 associated HAM/TSP. Conclusions: Our study failed to support the previously reported positive association between the IL28B rs12979860 polymorphisms and an increased risk of developing HAM/TSP in the Brazilian population.

Genetic analyses of smoking initiation, persistence, quantity, and age-at-onset of regular cigarette use in Brazilian families: the Baependi Heart Study

Background: The purpose of this study was to estimate the genetic influences on the initiation of cigarette smoking, the persistence, quantity and age-at-onset of regular cigarette use in Brazilian families. Methods: The data set consisted of 1,694 individuals enrolled in the Baependi Heart Study. The heritability and the heterogeneity in genetic and environmental variance components by gender were estimated from variance components approaches, using the SOLAR (Sequential Oligogenic Linkage Analysis Routines) computer package. The mixed-effects Cox model was used for the genetic analysis of the age-at onset of regular cigarette use. Results: The heritability estimates were high (> 50%) for smoking initiation and were intermediate, ranging from 23.4 to 31.9%, for smoking persistence and quantity. Significant evidence for heterogeneity in variance components by gender was observed for smoking initiation and age-at-onset of regular cigarette use. Genetic factors play an important role in the interindividual variation of these phenotypes in females, while in males there is a predominant environmental component, which could be explained by greater social influences in the initiation of tobacco use. Conclusions: Significant heritabilities were observed in smoking phenotypes for both males and females from the Brazilian population. These data add to the literature and are concordant with the notion of significant biological determination in smoking behavior. Samples from the Baependi Heart Study may be valuable for the mapping of genetic loci that modulate this complex biological trait.

MYLIP p.N342S polymorphism is not associated with lipid profile in the Brazilian population

Background: A recent study investigated the MYLIP region in the Mexican population in order to fine-map the actual susceptibility variants of this locus. The p.N342S polymorphism was identified as the underlying functional variant accounting for one of the previous signals of genome-wide association studies and the N342 allele was associated with higher cholesterol concentrations in Mexican dyslipidemic individuals. To date, there is no further evaluation on this genotype-phenotype association in the literature. In this scenario, and because of a possible pharmacotherapeutic target of dyslipidemia, the main aim of this study was to assess the influence of the MYLIP p.N342S polymorphism on lipid profile in Brazilian individuals. Methods: 1295 subjects of the general population and 1425 consecutive patients submitted to coronary angiography were selected. General characteristics, biochemical tests, blood pressures, pulse wave velocity, and coronary artery disease scores were analyzed. Genotypes for the MYLIP rs9370867 (p.N342S, c.G1025A) polymorphism were detected by high resolution melting analysis. Results: No association of the MYLIP rs9370867 genotypes with lipid profile, hemodynamic data, and coronary angiographic data was found. Analysis stratified by hyperlipidemia, gender, and ethnicity was also performed and the sub-groups presented similar results. In both general population and patient samples, the MYLIP rs9370867 polymorphism was differently distributed according to ethnicity. In the general population, subjects carrying GG genotypes had higher systolic blood pressure (BP), diastolic BP, and mean BP values (129.0 +/- 23.3; 84.9 +/- 14.6; 99.5 +/- 16.8 mmHg) compared with subjects carrying AA genotypes (123.7 +/- 19.5; 81.6 +/- 11.8; 95.6 +/- 13.6 mmHg) (p = 0.01; p = 0.02; p = 0.01, respectively), even after adjustment for covariates. However, in analysis stratified by ethnicity, this finding was not found and there is no evidence that the polymorphism influences BP. Conclusion: Our findings indicate that association studies involving this MYLIP variant can present distinct results according to the studied population. In this moment, further studies are needed to reaffirm if the MYLIP p.N342S polymorphism is functional or not, and to identify other functional markers within this gene.

Translocation capture sequencing: A method for high throughput mapping of chromosomal rearrangements

Chromosomal translocations require formation and joining of DNA double strand breaks (DSBs). These events disrupt the integrity of the genome and are involved in producing leukemias, lymphomas and sarcomas. Translocations are frequent, clonal and recurrent in mature B cell lymphomas, which bear a particularly high DNA damage burden by virtue of activation-induced cytidine deaminase (AID) expression. Despite the ubiquity of genomic rearrangements, the forces that underlie their genesis are not well understood. Here, we provide a detailed description of a new method for studying these events, translocation capture sequencing (TC-Seq). TC-Seq provides the means to document chromosomal rearrangements genome-wide in primary cells, and to discover recombination hotspots. Demonstrating its effectiveness, we successfully estimate the frequency of c-myc/IgH translocations in primary B cells, and identify hotspots of AID-mediated recombination. Furthermore. TC-Seq can be adapted to generate genome-wide rearrangement maps in any cell type and under any condition. (C) 2011 Elsevier B.V. All rights reserved.

Absence of Functional LIN28B Mutations in a Large Cohort of Patients with Idiopathic Central Precocious Puberty

Aim: To investigate LIN28B gene variants in children with idiopathic central precocious puberty (CPP). Patients and Methods: We studied 178 Brazilian children with CPP (171 girls, 16.8% familial cases). A large multiethnic group (1,599 subjects; Multiethnic Cohort, MEC) was used as control. DNA analysis and biochemical in vitro studies were performed. Results: A heterozygous LIN28B variant, p. H199R, was identified in a girl who developed CPP at 5.2 years. This variant was absent in 310 Brazilian control individuals, but it was found in the same allele frequency in women from the MEC cohort, independent of the age of menarche. Functional studies revealed that when ectopically expressed in cells, the mutant protein was capable of binding pre-let-7 microRNA and inhibiting let-7 expression to the same extent as wild-type Lin28B protein. Other rare LIN28B variants (p.P173P, c.198+32_33delCT, g.9575731A>C and c.-11C>T) were identified in CPP patients and controls. Therefore, no functional mutation was identified. Conclusion: In vitro studies revealed that the rare LIN28B p.H199R variant identified in a girl with CPP does not affect the Lin28B function in the regulation of let-7 expression. Although LIN28B SNPs were associated with normal pubertal timing, rare variations in this gene do not seem to be commonly involved in the molecular pathogenesis of CPP. Copyright (C) 2012 S. Karger AG, Basel

Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

Higher incidence of death in multi-vessel coronary artery disease patients associated with polymorphisms in chromosome 9p21

Background: We investigated whether 9p21 polymorphisms are associated with cardiovascular events in a group of 611 patients enrolled in the Medical, Angioplasty or Surgery Study II (MASS II), a randomized trial comparing treatments for patients with coronary artery disease (CAD) and preserved left ventricular function. Methods: The participants of the MASS II were genotyped for 9p21 polymorphisms (rs10757274, rs2383206, rs10757278 and rs1333049). Survival curves were calculated with the Kaplan-Meier method and compared with the log-rank statistic. We assessed the relationship between baseline variables and the composite end-point of death, death from cardiac causes and myocardial infarction using a Cox proportional hazards survival model. Results: We observed significant differences between patients within each polymorphism genotype group for baseline characteristics. The frequency of diabetes was lower in patients carrying GG genotype for rs10757274, rs2383206 and rs10757278 (29.4%, 32.8%, 32.0%) compared to patients carrying AA or AG genotypes (49.1% and 39.2%, p = 0.01; 52.4% and 40.1%, p = 0.01; 47.8% and 37.9%, p = 0.04; respectively). Significant differences in genotype frequencies between double and triple vessel disease patients were observed for the rs10757274, rs10757278 and rs1333049. Finally, there was a higher incidence of overall mortality in patients with the GG genotype for rs2383206 compared to patients with AA and AG genotypes (19.5%, 11.9%, 11.0%, respectively; p = 0.04). Moreover, the rs2383206 was still significantly associated with a 1.75-fold increased risk of overall mortality (p = 0.02) even after adjustment of a Cox multivariate model for age, previous myocardial infarction, diabetes, smoking and type of coronary anatomy. Conclusions: Our data are in accordance to previous evidence that chromosome 9p21 genetic variation may constitute a genetic modulator in the cardiovascular system in different scenarios. In patients with established CAD, we observed an association between the rs2383206 and higher incidence of overall mortality and death from cardiac causes in patients with multi-vessel CAD.