Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis.

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

MHC class I expression dependent on bacterial infection and parental factors in whitefish embryos (Salmonidae).

Ecological conditions can influence not only the expression of a phenotype, but also the heritability of a trait. As such, heritable variation for a trait needs to be studied across environments. We have investigated how pathogen challenge affects the expression of MHC genes in embryos of the lake whitefish Coregonus palaea. In order to experimentally separate paternal (i.e. genetic) from maternal and environmental effects, and determine whether and how stress affects the heritable variation for MHC expression, embryos were produced in full-factorial in vitro fertilizations, reared singly, and exposed at 208 degree days (late-eyed stage) to either one of two strains of Pseudomonas fluorescens that differ in their virulence characteristics (one increased mortality, while both delayed hatching time). Gene expression was assessed 48 h postinoculation, and virulence effects of the bacterial infection were monitored until hatching. We found no evidence of MHC class II expression at this stage of development. MHC class I expression was markedly down-regulated in reaction to both pseudomonads. While MHC expression could not be linked to embryo survival, the less the gene was expressed, the earlier the embryos hatched within each treatment group, possibly due to trade-offs between immune function and developmental rate or further factors that affect both hatching timing and MHC expression. We found significant additive genetic variance for MHC class I expression in some treatments. That is, changes in pathogen pressures could induce rapid evolution in MHC class I expression. However, we found no additive genetic variance in reaction norms in our study population.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Pleiotropy in the melanocortin system: expression levels of this system are associated with melanogenesis and pigmentation in the tawny owl (Strix aluco).

The adaptive function of melanin-based coloration is a long-standing debate. A recent genetic model suggested that pleiotropy could account for covariations between pigmentation, behaviour, morphology, physiology and life history traits. We explored whether the expression levels of genes belonging to the melanocortin system (MC1R, POMC, PC1/3, PC2 and the antagonist ASIP), which have many pleiotropic effects, are associated with melanogenesis (through variation in the expression of the genes MITF, SLC7A11, TYR, TYRP1) and in turn melanin-based coloration. We considered the tawny owl (Strix aluco) because individuals vary continuously from light to dark reddish, and thus, colour variation is likely to stem from differences in the levels of gene expression. We measured gene expression in feather bases collected in nestlings at the time of melanin production. As expected, the melanocortin system was associated with the expression of melanogenic genes and pigmentation. Offspring of darker reddish fathers expressed PC1/3 to lower levels but tended to express PC2 to higher levels. The convertase enzyme PC1/3 cleaves the POMC prohormone to obtain ACTH, while the convertase enzyme PC2 cleaves ACTH to produce α-melanin-stimulating hormone (α-MSH). ACTH regulates glucocorticoids, hormones that modulate stress responses, while α-MSH induces eumelanogenesis. We therefore conclude that the melanocortin system, through the convertase enzymes PC1/3 and PC2, may account for part of the interindividual variation in melanin-based coloration in nestling tawny owls. Pleiotropy may thus account for the covariation between phenotypic traits involved in social interactions (here pigmentation) and life history, morphology, behaviour and physiology.

Phenotypic plasticity in salmonid embryos : characterizing pathogen-induced stress effects on heritable variation in traits and reaction norms

Les pressions écologiques peuvent varier tant en nature qu'en intensité dans le temps et l'espace. C'est pourquoi, un phénotype unique ne peut pas forcément conférer la meilleure valeur sélective. La plasticité phénotypique peut être un moyen de s'accommoder de cette situation, en augmentant globalement la tolérance aux changements environnementaux. Comme pour tout trait de caractère, une variation génétique doit persister pour qu'évoluent les traits plastiques dans une population donnée. Cependant, les pressions extérieures peuvent affecter l'héritabilité, et la direction de ces changements peut dépendre du caractère en question, de l'espèce mais aussi du type de stress. Dans la présente thèse, nous avons cherché à élucider les effets des pressions pathogéniques sur les phénotypes et la génétique quantitative de plusieurs traits plastiques chez les embryons de deux salmonidés, la palée (Coregonus palaea), et la truite de rivière (Salmo trutta). Les salmonidés se prêtent à de telles études du fait de leur extraordinaire variabilité morphologique, comportementale et des traits d'histoire de vie. Par ailleurs, avec le déclin des salmonidés dans le monde, il est important de savoir combien la variabilité génétique persiste dans les normes de réaction afin d'aider à prédire leur capacité à répondre aux changements de leur milieu. Nous avons observé qu'une augmentation de la croissance des communautés microbiennes symbiotiques entraînait une mortalité accrue et une éclosion précoce chez la palée, et dévoilait la variance génétique additive pour ces deux caractères (Chapitres 1-2). Bien qu'aucune variation génétique n'ait été trouvée pour les normes de réaction, nous avons observé une variabilité de la plasticité d'éclosion. Néanmoins, on a trouvé que les temps d'éclosion étaient corrélés entre les environnements, ce qui pourrait limiter l'évolution de la norme de réaction. Le temps d'éclosion des embryons est lié à la taille des géniteurs mâles, ce qui indique des effets pléiotropiques. Dans le Chapitre 3, nous avons montré qu'une interaction triple entre la souche bactérienne {Pseudomonas fluorescens}, l'état de dévelopement de l'hôte ainsi que ses gènes ont une influence sur la mortalité, le temps d'éclosion et la taille des alevins de la palée. Nous avons démontré qu'une variation génétique subsistait généralement dans les normes de réaction des temps d'éclosion, mais rarement pour la taille des alevins, et jamais pour la mortalité. Dans le même temps, nous avons exhibé que des corrélations entre environnements dépendaient des caractères phénotypiques, mais contrairement au Chapitre 2, nous n'avons pas trouvé de preuve de corrélations transgénérationnelles. Le Chapitre 4 complète le chapitre précédent, en se plaçant du point de vue moléculaire, et décrit comment le traitement d'embryons avec P. fluorescens s'est traduit par une régulation négative d'expression du CMH-I indépendemment de la souche bactérienne. Nous avons non seulement trouvé une variation génétique des caractères phénotypiques moyens, mais aussi de la plasticité. Les deux derniers chapitres traitent de l'investigation, chez la truite de rivière, des différences spécifiques entre populations pour des normes de réaction induites par les pathogènes. Dans le Chapitre 5, nous avons illustré que le métissage entre des populations génétiquement distinctes n'affectait en rien la hauteur ou la forme des normes de réaction d'un trait précoce d'histoire de vie suite au traitement pathogénique. De surcroît, en dépit de l'éclosion tardive et de la réduction de la taille des alevins, le traitement n'a pas modifié la variation héritable des traits de caractère. D'autre part, dans le Chapitre 6, nous avons démontré que le traitement d'embryons avec des stimuli contenus dans l'eau de conspécifiques infectés a entraîné des réponses propre à chaque population en terme de temps d'éclosion ; néanmoins, nous avons observé peu de variabilité génétique des normes de réaction pour ce temps d'éclosion au sein des populations. - Ecological stressors can vary in type and intensity over space and time, and as such, a single phenotype may not confer the highest fitness. Phenotypic plasticity can act as a means to accommodate this situation, increasing overall tolerance to environmental change. As with any trait, for plastic traits to evolve in a population, genetic variation must persist. However, environmental stress can alter trait heritability, and the direction of this shift can be trait, species, and stressor-dependent. In this thesis, we sought to understand the effects of pathogen stressors on the phenotypes and genetic architecture of several plastic traits in the embryos of two salmonids, the whitefish (Coregonus palaea), and the brown trout (Salmo trutta). Salmonids lend themselves to such studies because their extraordinary variability in morphological, behavioral, and life-history traits. Also, with declines in salmonids worldwide, knowing how much genetic variability persists in reaction norms may help predict their ability to respond to environmental change. We found that increasing growth of symbiotic microbial communities increased mortality and induced hatching in whitefish, and released additive genetic variance for both traits (Chapters 1-2). While no genetic variation was found for survival reaction norms, we did find variability in hatching plasticity. Nevertheless, hatching time was correlated across environments, which could constrain evolution of the reaction norm. Hatching time in the induced environment was also correlated to sire size, indicating pleiotropic effects. In Chapter 3 we report that a three-way interaction between bacterial strain (Pseudomonas fluorescens), host developmental stage, and host genetics impacted mortality, hatching time, and hatchling size in whitefish. We also showed that genetic variation generally persisted in hatching age reaction norms, but rarely for hatchling length, and never for mortality. At the same time, we demonstrated that cross-environmental correlations were trait-dependent, and unlike Chapter 2, we found no evidence of cross-generational correlations. Chapter 4 expands on the previous chapter, moving to the molecular level, and describes how treatment of embryos with P. fluorescens resulted in strain-independent downregulation of MHC class I. Genetic variation was evident not only in trait means, but also in plasticity. In the last two chapters, we investigated population level differences in pathogen- induced reaction norms in brown trout. In Chapter 5, we found that interbreeding between genetically distinct populations did not affect the elevation or shapes of the reaction norms of early life-history traits after pathogen challenge. Moreover, despite delaying hatching and reducing larval length, treatment produced no discernable shifts in heritable variation in traits. On the other hand, in Chapter 6, we found that treatment of embryos with water-borne cues from infected conspecifics elicited population-specific responses in terms of hatching time; however, we found little evidence of genetic variability in hatching reaction norms within populations. We have made considerable progress in understanding how pathogen stressors affect various early life-history traits in salmonid embryos. We have demonstrated that the effect of a particular stressor on heritable variation in these traits can vary according to the trait and species under consideration, in addition to the developmental stage of the host. Moreover, we found evidence of genetic variability in some, but not all reaction norms in whitefish and brown trout.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Inflammation-induced changes in expression and glycosylation of genetic variants of alpha 1-acid glycoprotein. Studies with human sera, primary cultures of human hepatocytes and transgenic mice.

The relative occurrence of genetic variants of human alpha 1-acid glycoprotein (AGP) in relation to changes in glycosylation was studied in sera of patients with burn injury, media of cytokine-treated primary cultures of human hepatocytes and Hep 3B cells, and sera of transgenic mice expressing the human AGP-A gene. It is concluded (i) that the glycosylation of AGP was not dependent on its genetic expression and (ii) that both the variants determined by the AGP-A gene as well as by the AGP-B/B' genes are increased after inflammation or treatment with interleukins 1 and 6.

Development of stable cell lines for production or regulated expression using matrix attachment regions.

One of the major hurdles of isolating stable, inducible or constitutive high-level producer cell lines is the time-consuming selection procedure. Given the variation in the expression levels of the same construct in individual clones, hundreds of clones must be isolated and tested to identify one or more with the desired characteristics. Various boundary elements (BEs), matrix attachment regions, and locus control regions (LCRs) were screened for their ability to augment the expression of heterologous genes in Chinese hamster ovary (CHO) cells. Of the chromatin elements assayed, the chicken lysozyme matrix-attachment region (MAR) was the only element to significantly increase stable reporter expression. We found that the use of the MAR increases the proportion of high-producing clones, thus reducing the number of clones that need to be screened. These benefits are observed both for constructs with MARs flanking the transgene expression cassette, as well as when constructs are co-transfected with the MAR on a separate plasmid. Moreover, the MAR was co-transfected with a multicomponent regulatable beta-galactosidase expression system in C2C12 cells and several clones exhibiting regulated expression were identified. Hence, MARs are useful in the development of stable cell lines for production or regulated expression.

Copy number variants, diseases and gene expression.

Copy number variation (CNV) has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals in humans and a number of model organisms, using bioinformatics or hybridization-based methods. These have allowed uncovering associations between copy number changes and complex diseases in whole-genome association studies, as well as identify new genomic disorders. At the genome-wide scale, however, the functional impact of CNV remains poorly studied. Here we review the current catalogs of CNVs, their association with diseases and how they link genotype and phenotype. We describe initial evidence which revealed that genes in CNV regions are expressed at lower and more variable levels than genes mapping elsewhere, and also that CNV not only affects the expression of genes varying in copy number, but also have a global influence on the transcriptome. Further studies are warranted for complete cataloguing and fine mapping of CNVs, as well as to elucidate the different mechanisms by which they influence gene expression.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

CleanEx: a database of heterogeneous gene expression data based on a consistent gene nomenclature and linked to an improved annotation system

Résumé: L'automatisation du séquençage et de l'annotation des génomes, ainsi que l'application à large échelle de méthodes de mesure de l'expression génique, génèrent une quantité phénoménale de données pour des organismes modèles tels que l'homme ou la souris. Dans ce déluge de données, il devient très difficile d'obtenir des informations spécifiques à un organisme ou à un gène, et une telle recherche aboutit fréquemment à des réponses fragmentées, voir incomplètes. La création d'une base de données capable de gérer et d'intégrer aussi bien les données génomiques que les données transcriptomiques peut grandement améliorer la vitesse de recherche ainsi que la qualité des résultats obtenus, en permettant une comparaison directe de mesures d'expression des gènes provenant d'expériences réalisées grâce à des techniques différentes. L'objectif principal de ce projet, appelé CleanEx, est de fournir un accès direct aux données d'expression publiques par le biais de noms de gènes officiels, et de représenter des données d'expression produites selon des protocoles différents de manière à faciliter une analyse générale et une comparaison entre plusieurs jeux de données. Une mise à jour cohérente et régulière de la nomenclature des gènes est assurée en associant chaque expérience d'expression de gène à un identificateur permanent de la séquence-cible, donnant une description physique de la population d'ARN visée par l'expérience. Ces identificateurs sont ensuite associés à intervalles réguliers aux catalogues, en constante évolution, des gènes d'organismes modèles. Cette procédure automatique de traçage se fonde en partie sur des ressources externes d'information génomique, telles que UniGene et RefSeq. La partie centrale de CleanEx consiste en un index de gènes établi de manière hebdomadaire et qui contient les liens à toutes les données publiques d'expression déjà incorporées au système. En outre, la base de données des séquences-cible fournit un lien sur le gène correspondant ainsi qu'un contrôle de qualité de ce lien pour différents types de ressources expérimentales, telles que des clones ou des sondes Affymetrix. Le système de recherche en ligne de CleanEx offre un accès aux entrées individuelles ainsi qu'à des outils d'analyse croisée de jeux de donnnées. Ces outils se sont avérés très efficaces dans le cadre de la comparaison de l'expression de gènes, ainsi que, dans une certaine mesure, dans la détection d'une variation de cette expression liée au phénomène d'épissage alternatif. Les fichiers et les outils de CleanEx sont accessibles en ligne (http://www.cleanex.isb-sib.ch/). Abstract: The automatic genome sequencing and annotation, as well as the large-scale gene expression measurements methods, generate a massive amount of data for model organisms. Searching for genespecific or organism-specific information througout all the different databases has become a very difficult task, and often results in fragmented and unrelated answers. The generation of a database which will federate and integrate genomic and transcriptomic data together will greatly improve the search speed as well as the quality of the results by allowing a direct comparison of expression results obtained by different techniques. The main goal of this project, called the CleanEx database, is thus to provide access to public gene expression data via unique gene names and to represent heterogeneous expression data produced by different technologies in a way that facilitates joint analysis and crossdataset comparisons. A consistent and uptodate gene nomenclature is achieved by associating each single gene expression experiment with a permanent target identifier consisting of a physical description of the targeted RNA population or the hybridization reagent used. These targets are then mapped at regular intervals to the growing and evolving catalogues of genes from model organisms, such as human and mouse. The completely automatic mapping procedure relies partly on external genome information resources such as UniGene and RefSeq. The central part of CleanEx is a weekly built gene index containing crossreferences to all public expression data already incorporated into the system. In addition, the expression target database of CleanEx provides gene mapping and quality control information for various types of experimental resources, such as cDNA clones or Affymetrix probe sets. The Affymetrix mapping files are accessible as text files, for further use in external applications, and as individual entries, via the webbased interfaces . The CleanEx webbased query interfaces offer access to individual entries via text string searches or quantitative expression criteria, as well as crossdataset analysis tools, and crosschip gene comparison. These tools have proven to be very efficient in expression data comparison and even, to a certain extent, in detection of differentially expressed splice variants. The CleanEx flat files and tools are available online at: http://www.cleanex.isbsib. ch/.

Trade-off between constitutive and inducible resistance against herbivores is only partially explained by gene expression and glucosinolate production.

The hypothesis that constitutive and inducible plant resistance against herbivores should trade-off because they use the same resources and impose costs to plant fitness has been postulated for a long time. Negative correlations between modes of deployment of resistance and defences have been observed across and within species in common garden experiments. It was therefore tested whether that pattern of resistance across genotypes follows a similar variation in patterns of gene expression and chemical defence production. Using the genetically tractable model Arabidopsis thaliana and different modes of induction, including the generalist herbivore Spodoptera littoralis, the specialist herbivore Pieris brassicae, and jasmonate application, constitutive and inducibility of resistance was measured across seven A. thaliana accessions that were previously selected based on constitutive levels of defence gene expression. According to theory, it was found that modes of resistance traded-off among accessions, particularly against S. littoralis, in which accessions investing in high constitutive resistance did not increase it substantially after attack and vice-versa. Accordingly, the average expression of eight genes involved in glucosinolate production negatively predicted larval growth across the seven accessions. Glucosinolate production and genes related to defence induction on healthy and herbivore-damaged plants were measured next. Surprisingly, only a partial correlation between glucosinolate production, gene expression, and the herbivore resistance results was found. These results suggest that the defence outcome of plants against herbivores goes beyond individual molecules or genes but stands on a complex network of interactions.

Learning-induced gene expression in the heads of two Nasonia species that differ in long-term memory formation.

BACKGROUND: Cellular processes underlying memory formation are evolutionary conserved, but natural variation in memory dynamics between animal species or populations is common. The genetic basis of this fascinating phenomenon is poorly understood. Closely related species of Nasonia parasitic wasps differ in long-term memory (LTM) formation: N. vitripennis will form transcription-dependent LTM after a single conditioning trial, whereas the closely-related species N. giraulti will not. Genes that were differentially expressed (DE) after conditioning in N. vitripennis, but not in N. giraulti, were identified as candidate genes that may regulate LTM formation. RESULTS: RNA was collected from heads of both species before and immediately, 4 or 24 hours after conditioning, with 3 replicates per time point. It was sequenced strand-specifically, which allows distinguishing sense from antisense transcripts and improves the quality of expression analyses. We determined conditioning-induced DE compared to naïve controls for both species. These expression patterns were then analysed with GO enrichment analyses for each species and time point, which demonstrated an enrichment of signalling-related genes immediately after conditioning in N. vitripennis only. Analyses of known LTM genes and genes with an opposing expression pattern between the two species revealed additional candidate genes for the difference in LTM formation. These include genes from various signalling cascades, including several members of the Ras and PI3 kinase signalling pathways, and glutamate receptors. Interestingly, several other known LTM genes were exclusively differentially expressed in N. giraulti, which may indicate an LTM-inhibitory mechanism. Among the DE transcripts were also antisense transcripts. Furthermore, antisense transcripts aligning to a number of known memory genes were detected, which may have a role in regulating these genes. CONCLUSION: This study is the first to describe and compare expression patterns of both protein-coding and antisense transcripts, at different time points after conditioning, of two closely related animal species that differ in LTM formation. Several candidate genes that may regulate differences in LTM have been identified. This transcriptome analysis is a valuable resource for future in-depth studies to elucidate the role of candidate genes and antisense transcription in natural variation in LTM formation.

Genetic Variation in the TP53 Pathway and Bladder Cancer Risk. A Comprehensive Analysis

Introduction: Germline variants in TP63 have been consistently associated with several tumors, including bladder cancer, indicating the importance of TP53 pathway in cancer genetic susceptibility. However, variants in other related genes, including TP53 rs1042522 (Arg72Pro), still present controversial results. We carried out an in depth assessment of associations between common germline variants in the TP53 pathway and bladder cancer risk. Material and Methods: We investigated 184 tagSNPs from 18 genes in 1,058 cases and 1,138 controls from the Spanish Bladder Cancer/EPICURO Study. Cases were newly-diagnosed bladder cancer patients during 1998–2001. Hospital controls were age-gender, and area matched to cases. SNPs were genotyped in blood DNA using Illumina Golden Gate and TaqMan assays. Cases were subphenotyped according to stage/grade and tumor p53 expression. We applied classical tests to assess individual SNP associations and the Least Absolute Shrinkage and Selection Operator (LASSO)-penalized logistic regression analysis to assess multiple SNPs simultaneously. Results: Based on classical analyses, SNPs in BAK1 (1), IGF1R (5), P53AIP1 (1), PMAIP1 (2), SERINPB5 (3), TP63 (3), and TP73 (1) showed significant associations at p-value#0.05. However, no evidence of association, either with overall risk or with specific disease subtypes, was observed after correction for multiple testing (p-value$0.8). LASSO selected the SNP rs6567355 in SERPINB5 with 83% of reproducibility. This SNP provided an OR = 1.21, 95%CI 1.05–1.38, p-value = 0.006, and a corrected p-value = 0.5 when controlling for over-estimation. Discussion: We found no strong evidence that common variants in the TP53 pathway are associated with bladder cancer susceptibility. Our study suggests that it is unlikely that TP53 Arg72Pro is implicated in the UCB in white Europeans. SERPINB5 and TP63 variation deserve further exploration in extended studies.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.