Functional impact of genetic variation on gene expression

Numerous links between genetic variants and phenotypes are known and genome-wide association studies dramatically increased the number of genetic variants associated with traits during the last decade. However, how changes in the DNA perturb the molecular mechanisms and impact on the phenotype of an organism remains elusive. Studies suggest that many traitassociated variants are in the non-coding region of the genome and probably act through regulation of gene expression. During my thesis I investigated how genetic variants affect gene expression through gene regulatory mechanisms. The first chapter was a collaborative project with a pharmaceutical company, where we investigated genome-wide copy number variation (CNVs) among Cynomolgus monkeys (Macaca fascicularis) used in pharmaceutical studies, and associated them to changes in gene expression. We found substantial copy number variation and identified CNVs linked to tissue-specific expression changes of proximal genes. The second and third chapters focus on genetic variation in humans and its effects on gene regulatory mechanisms and gene expression. The second chapter studies two human trios, where the allelic effects of genetic variation on genome-wide gene expression, protein-DNA binding and chromatin modifications were investigated. We found abundant allele specific activity across all measured molecular phenotypes and show extended coordinated behavior among them. In the third chapter, we investigated the impact of genetic variation on these phenotypes in 47 unrelated individuals. We found that chromatin phenotypes are organized into local variable modules, often linked to genetic variation and gene expression. Our results suggest that chromatin variation emerges as a result of perturbations of cis-regulatory elements by genetic variants, leading to gene expression changes. The work of this thesis provides novel insights into how genetic variation impacts gene expression by perturbing regulatory mechanisms. -- De nombreux liens entre variations génétiques et phénotypes sont connus. Les études d'association pangénomique ont considérablement permis d'augmenter le nombre de variations génétiques associées à des phénotypes au cours de la dernière décennie. Cependant, comprendre comment ces changements perturbent les mécanismes moléculaires et affectent le phénotype d'un organisme nous échappe encore. Des études suggèrent que de nombreuses variations, associées à des phénotypes, sont situées dans les régions non codantes du génome et sont susceptibles d'agir en modifiant la régulation d'expression des gènes. Au cours de ma thèse, j'ai étudié comment les variations génétiques affectent les niveaux d'expression des gènes en perturbant les mécanismes de régulation de leur expression. Le travail présenté dans le premier chapitre est un projet en collaboration avec une société pharmaceutique. Nous avons étudié les variations en nombre de copies (CNV) présentes chez le macaque crabier (Macaca fascicularis) qui est utilisé dans les études pharmaceutiques, et nous les avons associées avec des changements d'expression des gènes. Nous avons découvert qu'il existe une variabilité substantielle du nombre de copies et nous avons identifié des CNVs liées aux changements d'expression des gènes situés dans leur voisinage. Ces associations sont présentes ou absentes de manière spécifique dans certains tissus. Les deuxième et troisième chapitres se concentrent sur les variations génétiques dans les populations humaines et leurs effets sur les mécanismes de régulation des gènes et leur expression. Le premier se penche sur deux trios humains, père, mère, enfant, au sein duquel nous avons étudié les effets alléliques des variations génétiques sur l'expression des gènes, les liaisons protéine-ADN et les modifications de la chromatine. Nous avons découvert que l'activité spécifique des allèles est abondante abonde dans tous ces phénotypes moléculaires et nous avons démontré que ces derniers ont un comportement coordonné entre eux. Dans le second, nous avons examiné l'impact des variations génétiques de ces phénotypes moléculaires chez 47 individus, sans lien de parenté. Nous avons observé que les phénotypes de la chromatine sont organisés en modules locaux, qui sont liés aux variations génétiques et à l'expression des gènes. Nos résultats suggèrent que la variabilité de la chromatine est due à des variations génétiques qui perturbent des éléments cis-régulateurs, et peut conduire à des changements dans l'expression des gènes. Le travail présenté dans cette thèse fournit de nouvelles pistes pour comprendre l'impact des différentes variations génétiques sur l'expression des gènes à travers les mécanismes de régulation.

Expression of the Mycobacterium bovis P36 gene in Mycobacterium smegmatis and the baculovirus/insect cell system

In the present study we evaluated different systems for the expression of mycobacterial antigen P36 secreted by Mycobacterium bovis. P36 was detected by Western blot using a specific antiserum. The P36 gene was initially expressed in E. coli, under the control of the T7 promoter, but severe proteolysis prevented its purification. We then tried to express P36 in M. smegmatis and insect cells. For M. smegmatis, we used three different plasmid vectors differing in copy number and in the presence of a promoter for expression of heterologous proteins. P36 was detected in the cell extract and culture supernatant in both expression systems and was recognized by sera from M. bovis-infected cattle. To compare the expression level and compartmentalization, the MPB70 antigen was also expressed. The highest production was reached in insect cell supernatants. In conclusion, M. smegmatis and especially the baculovirus expression system are good choices for the production of proteins from pathogenic mycobacteria for the development of mycobacterial vaccines and diagnostic reagents.

Nuclear organisation and epigenetic regulation of gene expression in Dictyostelium discoideum

A series of vectors for the over-expression of tagged proteins in Dictyostelium were designed, constructed and tested. These vectors allow the addition of an N- or C-terminal tag (GFP, RFP, 3xFLAG, 3xHA, 6xMYC and TAP) with an optimized polylinker sequence and no additional amino acid residues at the N or C terminus. Different selectable markers (Blasticidin and gentamicin) are available as well as an extra chromosomal version; these allow copy number and thus expression level to be controlled, as well as allowing for more options with regard to complementation, co- and super-transformation. Finally, the vectors share standardized cloning sites, allowing a gene of interest to be easily transfered between the different versions of the vectors as experimental requirements evolve. The organisation and dynamics of the Dictyostelium nucleus during the cell cycle was investigated. The centromeric histone H3 (CenH3) variant serves to target the kinetochore to the centromeres and thus ensures correct chromosome segregation during mitosis and meiosis. A number of Dictyostelium histone H3-domain containing proteins as GFP-tagged fusions were expressed and it was found that one of them functions as CenH3 in this species. Like CenH3 from some other species, Dictyostelium CenH3 has an extended N-terminal domain with no similarity to any other known proteins. The targeting domain, comprising α-helix 2 and loop 1 of the histone fold is required for targeting CenH3 to centromeres. Compared to the targeting domain of other known and putative CenH3 species, Dictyostelium CenH3 has a shorter loop 1 region. The localisation of a variety of histone modifications and histone modifying enzymes was examined. Using fluorescence in situ hybridisation (FISH) and CenH3 chromatin-immunoprecipitation (ChIP) it was shown that the six telocentric centromeres contain all of the DIRS-1 and most of the DDT-A and skipper transposons. During interphase the centromeres remain attached to the centrosome resulting in a single CenH3 cluster which also contains the putative histone H3K9 methyltransferase SuvA, H3K9me3 and HP1 (heterochromatin protein 1). Except for the centromere cluster and a number of small foci at the nuclear periphery opposite the centromeres, the rest of the nucleus is largely devoid of transposons and heterochromatin associated histone modifications. At least some of the small foci correspond to the distal telomeres, suggesting that the chromosomes are organised in a Rabl-like manner. It was found that in contrast to metazoans, loading of CenH3 onto Dictyostelium centromeres occurs in late G2 phase. Transformation of Dictyostelium with vectors carrying the G418 resistance cassette typically results in the vector integrating into the genome in one or a few tandem arrays of approximately a hundred copies. In contrast, plasmids containing a Blasticidin resistance cassette integrate as single or a few copies. The behaviour of transgenes in the nucleus was examined by FISH, and it was found that low copy transgenes show apparently random distribution within the nucleus, while transgenes with more than approximately 10 copies cluster at or immediately adjacent to the centromeres in interphase cells regardless of the actual integration site along the chromosome. During mitosis the transgenes show centromere-like behaviour, and ChIP experiments show that transgenes contain the heterochromatin marker H3K9me2 and the centromeric histone variant H3v1. This clustering, and centromere-like behaviour was not observed on extrachromosomal transgenes, nor on a line where the transgene had integrated into the extrachromosomal rDNA palindrome. This suggests that it is the repetitive nature of the transgenes that causes the centromere-like behaviour. A Dictyostelium homolog of DET1, a protein largely restricted to multicellular eukaryotes where it has a role in developmental regulation was identified. As in other species Dictyostelium DET1 is nuclear localised. In ChIP experiments DET1 was found to bind the promoters of a number of developmentally regulated loci. In contrast to other species where it is an essential protein, loss of DET1 is not lethal in Dictyostelium, although viability is greatly reduced. Loss of DET1 results in delayed and abnormal development with enlarged aggregation territories. Mutant slugs displayed apparent cell type patterning with a bias towards pre-stalk cell types.

Oscheius tipulae as an example of eEF1A gene diversity in nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.

Screening de polimorfismos e análise comparativa de alterações no número de copias do gene tp53 em gliomas

As neoplasias astrocitárias correspondem a cerca de (60%) dos tumores do sistema nervoso central, sendo que atingem principalmente adultos numa fase altamente produtiva da vida e com evoluções pouco satisfatórias apesar dos tratamentos disponíveis. Um melhor entendimento de sua biologia molecular se faz necessário na tentativa de compreender sua evolução e melhor planejar e tratamento, assim como na busca de novas terapias. Este trabalho teve como objetivo analisar modificações no gene TP53 com relação ao número de cópias e polimorfismos nos éxons de 4 a 11, considerados hotspots para mutações. Um total de 14 amostras de diferentes graus de malignidade foram analisadas por experimentos de FISH interfásico com sondas loco-específicas do gene TP53 e centroméricas para o cromossomo 17, e também pela técnica de SSCP para o screening de polimorfismos dos éxons 4-10. Foram comparados os resultados obtidos entre tumores de graus I e II (benignos) com aqueles de graus III e IV (malignos). Os resultados referentes às sondas loco-específicas (gene TP53 e centrômero do cromossomo 17) mostraram que a ocorrência de deleções ou amplificações, apesar de importantes estatisticamente em relação aos núcleos com número de marcações normais, não apresentou correlação com idade, sexo ou grau de malignação. Entretanto, as alterações foram encontradas com maior freqüência nos paciente portadores de astrocitomas de grau intermediário (III). A técnica de SSCP revelou polimorfismos nos éxons 5, 7 e 10, e apesar de não estarem associados à malignidade tumoral, os casos polimórficos corresponderam aos pacientes com menor sobrevida após tratamentos, sugerindo, uma associação entre mutações nesses éxons e uma maior agressividade tumoral.

Evolutionary History of the PER3 Variable Number of Tandem Repeats (VNTR): Idiosyncratic Aspect of Primate Molecular Circadian Clock

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Alterações epigenéticas e do número de cópias do gene SFN em carcinomas mamários

Pós-graduação em Ciências Biológicas (Genética) - IBB

Organização genômima e análise da expressão do gene hnRNP Q-like (Heterogeneous nuclear ribonucleoprotein A-like) retroinserido no cromossomo B de Astatotilapia latifasciata

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Obtenção e caracterização molecular e fisiológica de plantas de soja contendo o gene AtGolS2 sob déficit hídrico

Pós-graduação em Agronomia (Genética e Melhoramento de Plantas) - FCAV

Detection of exon deletions within an entire gene (CFTR) by relative quantification on the LightCycler

BACKGROUND: Cystic fibrosis (CF) is associated with at least 1 pathogen point sequence variant on each CFTR allele. Some symptomatic patients, however, have only 1 detectable pathogen sequence variant and carry, on the other allele, a large deletion that is not detected by conventional screening methods. METHODS: For relative quantitative real-time PCR detection of large deletions in the CFTR gene, we designed DNA-specific primers for each exon of the gene and primers for a reference gene (beta2-microglobulin). For PCR we used a LightCycler system (Roche) and calculated the gene-dosage ratio of CFTR to beta2-microglobulin. We tested the method by screening all 27 exons in 3 healthy individuals and 2 patients with only 1 pathogen sequence variant. We then performed specific deletion screenings in 10 CF patients with known large deletions and a blinded analysis in which we screened 24 individuals for large deletions by testing 8 of 27 exons. RESULTS: None of the ratios for control samples were false positive (for deletions or duplications); moreover, for all samples from patients with known large deletions, the calculated ratios for deleted exons were close to 0.5. In addition, the results from the blinded analysis demonstrated that our method can also be used for the screening of single individuals. CONCLUSIONS: The LightCycler assay allows reliable and rapid screening for large deletions in the CFTR gene and detects the copy number of all 27 exons.

Gene duplication: a drive for phenotypic diversity and cause of human disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

Large genomic fibrillin-1 (FBN1) gene deletions provide evidence for true haploinsufficiency in Marfan syndrome

Mutations in the FBN1 gene are the major cause of Marfan syndrome (MFS), an autosomal dominant connective tissue disorder, which displays variable manifestations in the cardiovascular, ocular, and skeletal systems. Current molecular genetic testing of FBN1 may miss mutations in the promoter region or in other noncoding sequences as well as partial or complete gene deletions and duplications. In this study, we tested for copy number variations by successively applying multiplex ligation-dependent probe amplification (MLPA) and the Affymetrix Human Mapping 500 K Array Set, which contains probes for approximately 500,000 single-nucleotide polymorphisms (SNPs) across the genome. By analyzing genomic DNA of 101 unrelated individuals with MFS or related phenotypes in whom standard genetic testing detected no mutation, we identified FBN1 deletions in two patients with MFS. Our high-resolution approach narrowed down the deletion breakpoints. Subsequent sequencing of the junctional fragments revealed the deletion sizes of 26,887 and 302,580 bp, respectively. Surprisingly, both deletions affect the putative regulatory and promoter region of the FBN1 gene, strongly indicating that they abolish transcription of the deleted allele. This expectation of complete loss of function of one allele, i.e. true haploinsufficiency, was confirmed by transcript analyses. Our findings not only emphasize the importance of screening for large genomic rearrangements in comprehensive genetic testing of FBN1 but, importantly, also extend the molecular etiology of MFS by providing hitherto unreported evidence that true haploinsufficiency is sufficient to cause MFS.

The human granulocyte-macrophage colony- stimulating factor gene is autonomously regulated in vivo by an inducible tissue-specific enhancer

The granulocyte-macrophage colony-stimulating factor (GM-CSF) gene is part of a cytokine gene cluster and is directly linked to a conserved upstream inducible enhancer. Here we examined the in vitro and in vivo functions of the human GM-CSF enhancer and found that it was required for the correctly regulated expression of the GM-CSF gene. An inducible DNase I-hypersensitive site appeared within the enhancer in cell types such as T cells, myeloid cells, and endothelial cells that express GM-CSF, but not in nonexpressing cells. In a panel of transfected cells the human GM-CSF enhancer was activated in a tissue-specific manner in parallel with the endogenous gene. The in vivo function of the enhancer was examined in a transgenic mouse model that also addressed the issue of whether the GM-CSF locus was correctly regulated in isolation from other segments of the cytokine gene cluster. After correction for copy number the mean level of human GM-CSF expression in splenocytes from 11 lines of transgenic mice containing a 10.5-kb human GM-CSF transgene was indistinguishable from mouse GM-CSF expression (99% ± 56% SD). In contrast, a 9.8-kb transgene lacking just the enhancer had a significantly reduced (P = 0.004) and more variable level of activity (29% ± 89% SD). From these studies we conclude that the GM-CSF enhancer is required for the correct copy number-dependent expression of the human GM-CSF gene and that the GM-CSF gene is regulated independently from DNA elements associated with the closely linked IL-3 gene or other members of the cytokine gene cluster.

Homology-dependent Gene Silencing in Paramecium

Microinjection at high copy number of plasmids containing only the coding region of a gene into the Paramecium somatic macronucleus led to a marked reduction in the expression of the corresponding endogenous gene(s). The silencing effect, which is stably maintained throughout vegetative growth, has been observed for all Paramecium genes examined so far: a single-copy gene (ND7), as well as members of multigene families (centrin genes and trichocyst matrix protein genes) in which all closely related paralogous genes appeared to be affected. This phenomenon may be related to posttranscriptional gene silencing in transgenic plants and quelling in Neurospora and allows the efficient creation of specific mutant phenotypes thus providing a potentially powerful tool to study gene function in Paramecium. For the two multigene families that encode proteins that coassemble to build up complex subcellular structures the analysis presented herein provides the first experimental evidence that the members of these gene families are not functionally redundant.

Long-term expression of γ-globin mRNA in mouse erythrocytes from retrovirus vectors containing the human γ-globin gene fused to the ankyrin-1 promoter

Gene therapy for patients with hemoglobin disorders has been hampered by the inability of retrovirus vectors to transfer globin genes and their cis-acting regulatory sequences into hematopoietic stem cells without rearrangement. In addition, the expression from intact globin gene vectors has been variable in red blood cells due to position effects and retrovirus silencing. We hypothesized that by substituting the globin gene promoter for the promoter of another gene expressed in red blood cells, we could generate stable retrovirus vectors that would express globin at sufficient levels to treat hemoglobinopathies. Recently, we have shown that the human ankyrin (Ank) gene promoter directs position-independent, copy number-dependent expression of a linked γ-globin gene in transgenic mice. We inserted the Ank/Aγ-globin gene into retrovirus vectors that could transfer one or two copies of the Ank/Aγ-globin gene to target cells. Both vectors were stable, transferring only intact proviral sequences into primary mouse hematopoietic stem cells. Expression of Ank/Aγ-globin mRNA in mature red blood cells was 3% (single copy) and 8% (double copy) of the level of mouse α-globin mRNA. We conclude that these novel retrovirus vectors may be valuable for treating a variety of red cell disorders by gene replacement therapy including severe β-thalassemia if the level of expression can be further increased.