Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Gene expression variation between distinct areas of breast cancer measured from paraffin-embedded tissue cores

BACKGROUND: Diagnosis and prognosis in breast cancer are mainly based on histology and immunohistochemistry of formalin-fixed, paraffin-embedded (FFPE) material. Recently, gene expression analysis was shown to elucidate the biological variance between tumors and molecular markers were identified that led to new classification systems that provided better prognostic and predictive parameters. Archived FFPE samples represent an ideal source of tissue for translational research, as millions of tissue blocks exist from routine diagnostics and from clinical studies. These should be exploited to provide clinicians with more accurate prognostic and predictive information. Unfortunately, RNA derived from FFPE material is partially degraded and chemically modified and reliable gene expression measurement has only become successful after implementing novel and optimized procedures for RNA isolation, demodification and detection. METHODS: In this study we used tissue cylinders as known from the construction of tissue microarrays. RNA was isolated with a robust protocol recently developed for RNA derived from FFPE material. Gene expression was measured by quantitative reverse transcription PCR. RESULTS: Sixteen tissue blocks from 7 patients diagnosed with multiple histological subtypes of breast cancer were available for this study. After verification of appropriate localization, sufficient RNA yield and quality, 30 tissue cores were available for gene expression measurement on TaqMan(R) Low Density Arrays (16 invasive ductal carcinoma (IDC), 8 ductal carcinoma in situ (DCIS) and 6 normal tissue), and 14 tissue cores were lost. Gene expression values were used to calculate scores representing the proliferation status (PRO), the estrogen receptor status and the HER2 status. The PRO scores measured from entire sections were similar to PRO scores determined from IDC tissue cores. Scores determined from normal tissue cores consistently revealed lower PRO scores than cores derived from IDC or DCIS of the same block or from different blocks of the same patient. CONCLUSION: We have developed optimized protocols for RNA isolation from histologically distinct areas. RNA prepared from FFPE tissue cores is suitable for gene expression measurement by quantitative PCR. Distinct molecular scores could be determined from different cores of the same tumor specimen.

High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity

Background: Ecosystems worldwide are suffering the consequences of anthropogenic impact. The diverse ecosystem of coral reefs, for example, are globally threatened by increases in sea surface temperatures due to global warming. Studies to date have focused on determining genetic diversity, the sequence variability of genes in a species, as a proxy to estimate and predict the potential adaptive response of coral populations to environmental changes linked to climate changes. However, the examination of natural gene expression variation has received less attention. This variation has been implicated as an important factor in evolutionary processes, upon which natural selection can act. Results: We acclimatized coral nubbins from six colonies of the reef-building coral Acropora millepora to a common garden in Heron Island (Great Barrier Reef, GBR) for a period of four weeks to remove any site-specific environmental effects on the physiology of the coral nubbins. By using a cDNA microarray platform, we detected a high level of gene expression variation, with 17% (488) of the unigenes differentially expressed across coral nubbins of the six colonies (jsFDR-corrected, p < 0.01). Among the main categories of biological processes found differentially expressed were transport, translation, response to stimulus, oxidation-reduction processes, and apoptosis. We found that the transcriptional profiles did not correspond to the genotype of the colony characterized using either an intron of the carbonic anhydrase gene or microsatellite loci markers. Conclusion: Our results provide evidence of the high inter-colony variation in A. millepora at the transcriptomic level grown under a common garden and without a correspondence with genotypic identity. This finding brings to our attention the importance of taking into account natural variation between reef corals when assessing experimental gene expression differences. The high transcriptional variation detected in this study is interpreted and discussed within the context of adaptive potential and phenotypic plasticity of reef corals. Whether this variation will allow coral reefs to survive to current challenges remains unknown.

Variation of genetic expression during development, revealed by esterase patterns in Aedes aegypti (Diptera, Culicidae)

Polyacrylamide gel electrophoresis was used to analyze esterase patterns during development of Aedes aegypti from the cities of Marília and São José do Rio Preto (SJRP), Brazil. The zymograms showed a total of 23 esterase bands, 22 of which were in the specimens from Marília and 19 in those from SJRP. These esterase bands were considered to be the product of 23 alleles distributed tentatively in eight genetic loci. Most of the alleles were developmentally regulated. The larval stage expressed the greatest number of them (19 alleles, from the eight loci, in Marília; and 17 alleles, from seven loci, in SJRP). The pupal stage expressed 10 alleles from seven loci, in both populations, and the adult stage expressed 8 alleles from five and six loci in SJRP and Marília, respectively. Some alleles that were active in every stage were developmentally controlled at the level of expression (amount of product). A single allele was constitutively and highly expressed, in larvae, pupae, and adults, in both populations. Differences in esterase synthesis among stages are probably due to regulatory mechanisms acting in agreement with the requirements of a variable number of processes in which esterases are involved. The larval stage is the most active in developmental processes and shows very intense intake of food and very high mobility. These features may demand increased esterase production at that stage. Comparison of the two populations examined showed (besides the existence of alleles that they do not share) that they exhibit differences in the control of expression of other alleles. Such findings may reflect genetic differences between founders in each population, but the possibility of involvement of the intensive use of insecticides in SJRP is also discussed.

Extensive natural variation for cellular hydrogen peroxide release is genetically controlled.

Natural variation in DNA sequence contributes to individual differences in quantitative traits. While multiple studies have shown genetic control over gene expression variation, few additional cellular traits have been investigated. Here, we investigated the natural variation of NADPH oxidase-dependent hydrogen peroxide (H(2)O(2) release), which is the joint effect of reactive oxygen species (ROS) production, superoxide metabolism and degradation, and is related to a number of human disorders. We assessed the normal variation of H(2)O(2) release in lymphoblastoid cell lines (LCL) in a family-based 3-generation cohort (CEPH-HapMap), and in 3 population-based cohorts (KORA, GenCord, HapMap). Substantial individual variation was observed, 45% of which were associated with heritability in the CEPH-HapMap cohort. We identified 2 genome-wide significant loci of Hsa12 and Hsa15 in genome-wide linkage analysis. Next, we performed genome-wide association study (GWAS) for the combined KORA-GenCord cohorts (n = 279) using enhanced marker resolution by imputation (>1.4 million SNPs). We found 5 significant associations (p<5.00×10-8) and 54 suggestive associations (p<1.00×10-5), one of which confirmed the linked region on Hsa15. To replicate our findings, we performed GWAS using 58 HapMap individuals and ∼2.1 million SNPs. We identified 40 genome-wide significant and 302 suggestive SNPs, and confirmed genome signals on Hsa1, Hsa12, and Hsa15. Genetic loci within 900 kb from the known candidate gene p67phox on Hsa1 were identified in GWAS in both cohorts. We did not find replication of SNPs across all cohorts, but replication within the same genomic region. Finally, a highly significant decrease in H(2)O(2) release was observed in Down Syndrome (DS) individuals (p<2.88×10-12). Taken together, our results show strong evidence of genetic control of H(2)O(2) in LCL of healthy and DS cohorts and suggest that cellular phenotypes, which themselves are also complex, may be used as proxies for dissection of complex disorders.

Genomic architecture of sickle cell disease clinical variation in children from West Africa : a case-control study design

Contexte : L’anémie falciforme ou drépanocytose est un problème de santé important, particulièrement pour les patients d’origine africaine. La variation phénotypique de l’anémie falciforme est problématique pour le suivi et le traitement des patients. L’architecture génomique responsable de cette variabilité est peu connue. Principe : Mieux saisir la contribution génétique de la variation clinique de cette maladie facilitera l’identification des patients à risque de développer des phénotypes sévères, ainsi que l’adaptation des soins. Objectifs : L’objectif général de cette thèse est de combler les lacunes relatives aux connaissances sur l’épidémiologie génomique de l’anémie falciforme à l’aide d’une cohorte issue au Bénin. Les objectifs spécifiques sont les suivants : 1) caractériser les profils d’expressions génomiques associés à la sévérité de l’anémie falciforme ; 2) identifier des biomarqueurs de la sévérité de l’anémie falciforme ; 3) identifier la régulation génétique des variations transcriptionelles ; 4) identifier des interactions statistiques entre le génotype et le niveau de sévérité associé à l’expression ; 5) identifier des cibles de médicaments pour améliorer l’état des patients atteints d’anémie falciforme. Méthode : Une étude cas-témoins de 250 patients et 61 frères et soeurs non-atteints a été menée au Centre de Prise en charge Médical Intégré du Nourrisson et de la Femme Enceinte atteints de Drépanocytose, au Bénin entre février et décembre 2010. Résultats : Notre analyse a montré que des profils d’expressions sont associés avec la sévérité de l’anémie falciforme. Ces profils sont enrichis de génes des voies biologiques qui contribuent à la progression de la maladie : l’activation plaquettaire, les lymphocytes B, le stress, l’inflammation et la prolifération cellulaire. Des biomarqueurs transcriptionnels ont permis de distinguer les patients ayant des niveaux de sévérité clinique différents. La régulation génétique de la variation de l’expression des gènes a été démontrée et des interactions ont été identifiées. Sur la base de ces résultats génétiques, des cibles de médicaments sont proposées. Conclusion: Ce travail de thèse permet de mieux comprendre l’impact de la génomique sur la sévérité de l’anémie falciforme et ouvre des perspectives de développement de traitements ciblés pour améliorer les soins offerts aux patients.

Analysis of expression quantitative trait loci in D. melanogaster

Expression Quantitative Trait Loci (eQTL) analysis allows for the identification of genetic variation associated with variation in gene expression. It is often unclear however, which of the associated variants are causal, and by what mechanism. Integrating functional genomic data with eQTL data can provide insight into the impact of natural variation in the population, and the nature of the transcriptional machinery itself. In this thesis, I integrate functional genomic data with eQTL data derived from both 5’ CAGE and 3’ TagXseq expression assays, in developing embryos. I first use both datasets to analyse the transcription landscape in embryonic D., melanogaster, and then carry out an analysis of sequence motifs associated with transcription factor binding sites, promoters, and 3’ polyadenylation sites. Finally, I integrate functional genomic data, including these novel sequence motifs, to shed light on the mechanisms of gene expression variation in D.,melanogaster. I am able to demonstrate that some variants effecting gene regulation in Drosophila are found within haplotypes which buffer their effects.

Variation in protein expression levels in cell populations

Dissertation presented to obtain the Ph.D degree in Systems Biology

Variation d'expression des récepteurs hormonaux et d'HER-2 dans l'évolution du cancer du sein : quelles implications en pratique clinique [Hormone receptors and HER-2 changes during breast cancer progression: clinical implications].

Breast cancer remains a major public health problem. Even if there is an increase in this cancer curability, metastatic breast cancer remains a lethal disease in the vast majority of cases. Therapeutic advances in the chemotherapeutic and targeted therapies fields induced an increase in survival, however the proportion of long survivors remains low. Phenotypic instability, an early process initiated during tumour progression, and continued on the metastatic stage of the disease, can be one of the putative hypotheses explaining these results. An increasing amount of scientific data are pledging for a reanalysis of the phenotypic profile regarding hormone receptors and HER-2 status of metastatic lesions in order to identify drugable targets and allow individualisation of the treatment of these metastatic breast cancer patients. Phenotypic changes between the primary tumour and the paired metastatic lymph nodes are a challenging pitfall, raising the question of which site has to be assessed in the adjuvant treatment decision process. This article presents a comprehensive analysis of the frequency of theses phenotypic changes altogether with new modalities to evaluate this phenotypic status.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications.

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this "neighboring effect", we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples. We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype. GEO SERIES ACCESSION NUMBER: GSE33784, GSE33867.

Copy number variation modifies expression profiles over developmental time

Abstract : A preliminary understanding of the phenotypic effect of copy number variation (CNV) of DNA segments is emerging. These rearrangements were shown to influence, in a somewhat dose-dependent manner, the expression of genes mapping within them. They were also shown to modify the expression of genes located on their Hanks, sometimes at great distance. Here, we demonstrate by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, 'we find that some brain- expressed genes mapping within CNVS appear to be under compensatory loops only at specific time-points, indicating that the effect of CNVS on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time-course expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but its timing of expression as well. Résume : Nous commençons à comprendre les effets phénotypiques liés aux séquences d'ADN qui changent de nombre de copies d'un individu a l'autre. Des travaux précédents ont montré que ces variante de nombre de copies (CNVS) avaient une influence sur l'expression non seulement des gènes se trouvant dans le réarrangement, mais aussi sur ceux se trouvant à une certaine distance. Le présent travail étudie ces effets à différents stades du développement de la souris allant d'un embryon de deux semaines à la souris adulte. Nous avons observé que certains gènes exprimés dans le cerveau semblent soumis à un contrôle plus strict a certaines étapes du développement suggèrent que l'effet des CNVs est modulé différemment au cours de la vie. Notre travail sur trois souches différentes de souris a permis de montrer que les gènes ayant un profil d'expression différent dans le temps entre souches sont enrichis en gènes se trouvant dans des CNVs. Ceci nous amène à penser que les CNVs ont, non seulement une influence sur le niveau d'expression des gènes, mais aussi sur les moments durant lesquels ils seront exprimés. Résumé pour un large public : De nombreuses maladies sont dues soit a un gain (on parle alors de duplication) soit à une perte de matériel génétique (il s'agit dune délétion). Bien que les recherches visant à identifier les mécanismes moléculaires liés à ces réarrangements de notre génome progressent continuellement, la plupart des causes des maladies génétiques restent à élucider. Certaines parties de notre génome sont présentes en un nombre de copies qui diffère d'un individu à l'autre sans pour autant provoquer une ou des maladies. Ces segments d'ADN qui varient en nombre sont appelés Copy Number Variant (CNVs). Ils couvrent environ 12% de notre matériel génétique. Des études menées sur différents modèles animaux ont montré que les CNVs avaient une influence aussi bien sur les gènes qui sont a l'intérieur des CNVs que sur ceux qui sont dans leur voisinage. Ce travail étudie l'effet des CNVs à travers différents stades du développement de la souris. Nous avons démontré que les segments d'ADN qui varient en nombre de copies ont des effets variables selon le stade auxquels ils sont mesurés. Ainsi, les CNVs ont non seulement un impact sur l'expression des gènes présents dans ces régions et dans leur voisinage, mais influencent également leurs profils d'expression au cours du temps.

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

Copy number variation modifies expression time courses.

A preliminary understanding into the phenotypic effect of DNA segment copy number variation (CNV) is emerging. These rearrangements were demonstrated to influence, in a somewhat dose-dependent manner, the expression of genes that map within them. They were also shown to modify the expression of genes located on their flanks and sometimes those at a great distance from their boundary. Here we demonstrate, by monitoring these effects at multiple life stages, that these controls over expression are effective throughout mouse development. Similarly, we observe that the more specific spatial expression patterns of CNV genes are maintained through life. However, we find that some brain-expressed genes mapping within CNVs appear to be under compensatory loops only at specific time points, indicating that the effect of CNVs on these genes is modulated during development. Notably, we also observe that CNV genes are significantly enriched within transcripts that show variable time courses of expression between strains. Thus, modifying the copy number of a gene may potentially alter not only its expression level, but also the timing of its expression.