Gene-Environment Interaction Effects on Behavioral Variation and Risk of Complex Disorders: The Example of Alcoholism and Other Psychiatric Disorders

There have been few replicated examples of genotype x environment interaction effects on behavioral variation or risk of psychiatric disorder. We review some of the factors that have made detection of genotype x environment interaction effects difficult, and show how genotype x shared environment interaction (GxSE) effects are commonly confounded with genetic parameters in data from twin pairs reared together. Historic data on twin pairs reared apart can in principle be used to estimate such GxSE effects, but have rarely been used for this purpose. We illustrate this using previously published data from the Swedish Adoption Twin Study of Aging (SATSA), which suggest that GxSE effects could account for as much as 25% of the total variance in risk of becoming a regular smoker. Since few separated twin pairs will be available for study in the future, we also consider methods for modifying variance components linkage analysis to allow for environmental interactions with linked loci.

The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

Evaluating plant breeding strategies by simulating gene action and dryland environment effects

Functional genomics is the systematic study of genome-wide effects of gene expression on organism growth and development with the ultimate aim of understanding how networks of genes influence traits. Here, we use a dynamic biophysical cropping systems model (APSIM-Sorg) to generate a state space of genotype performance based on 15 genes controlling four adaptive traits and then search this spice using a quantitative genetics model of a plant breeding program (QU-GENE) to simulate recurrent selection. Complex epistatic and gene X environment effects were generated for yield even though gene action at the trait level had been defined as simple additive effects. Given alternative breeding strategies that restricted either the cultivar maturity type or the drought environment type, the positive (+) alleles for 15 genes associated with the four adaptive traits were accumulated at different rates over cycles of selection. While early maturing genotypes were favored in the Severe-Terminal drought environment type, late genotypes were favored in the Mild-Terminal and Midseason drought environment types. In the Severe-Terminal environment, there was an interaction of the stay-green (SG) trait with other traits: Selection for + alleles of the SG genes was delayed until + alleles for genes associated with the transpiration efficiency and osmotic adjustment traits had been fixed. Given limitations in our current understanding of trait interaction and genetic control, the results are not conclusive. However, they demonstrate how the per se complexity of gene X gene X environment interactions will challenge the application of genomics and marker-assisted selection in crop improvement for dryland adaptation.

Rapid Evolution of the Sequences and Gene Repertoires

Proteins secreted to the extracellular environment or to the periphery of the cell envelope, the secretome, play essential roles in foraging, antagonistic and mutualistic interactions. We hypothesize that arms races, genetic conflicts and varying selective pressures should lead to the rapid change of sequences and gene repertoires of the secretome. The analysis of 42 bacterial pan-genomes shows that secreted, and especially extracellular proteins, are predominantly encoded in the accessory genome, i.e. among genes not ubiquitous within the clade. Genes encoding outer membrane proteins might engage more frequently in intra-chromosomal gene conversion because they are more often in multi-genic families. The gene sequences encoding the secretome evolve faster than the rest of the genome and in particular at non-synonymous positions. Cell wall proteins in Firmicutes evolve particularly fast when compared with outer membrane proteins of Proteobacteria. Virulence factors are over-represented in the secretome, notably in outer membrane proteins, but cell localization explains more of the variance in substitution rates and gene repertoires than sequence homology to known virulence factors. Accordingly, the repertoires and sequences of the genes encoding the secretome change fast in the clades of obligatory and facultative pathogens and also in the clades of mutualists and free-living bacteria. Our study shows that cell localization shapes genome evolution. In agreement with our hypothesis, the repertoires and the sequences of genes encoding secreted proteins evolve fast. The particularly rapid change of extracellular proteins suggests that these public goods are key players in bacterial adaptation.

Gold nanoparticle-DNA conjugates for oligonucleotide vectorization towards gene silencing

Dissertação para obtenção do Grau de Mestre em Genética Molecular e Biomedicina

Gene Expression Profiling and Association Studies Implicate the Neuregulin Signaling Pathway in Behçet's Disease Susceptibility

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene–gene interactions. These association findings support a role for the EGF/ErbB signaling pathway inBD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.

Variant Rett Syndrome in a Girl with a Pericentric X-Chromosome Inversion Leading to Epigenetic Changes and Overexpression of the MECP2 Gene

Rett syndrome is a neurodevelopmental disorder caused by mutations in the MECP2 gene. We investigated the genetic basis of disease in a female patient with a Rett-like clinical. Karyotype analysis revealed a pericentric inversion in the X chromosome -46,X,inv(X)(p22.1q28), with breakpoints in the cytobands where the MECP2 and CDKL5 genes are located. FISH analysis revealed that the MECP2 gene is not dislocated by the inversion. However, and in spite of a balanced pattern of X inactivation, this patient displayed hypomethylation and an overexpression of the MECP2 gene at the mRNA level in the lymphocytes (mean fold change: 2.55±0.38) in comparison to a group of control individuals; the expression of the CDKL5 gene was similar to that of controls (mean fold change: 0.98±0.10). No gains or losses were detected in the breakpoint regions encompassing known or suspected transcription regulatory elements. We propose that the de-regulation of MECP2 expression in this patient may be due to alterations in long-range genomic interactions caused by the inversion and hypothesize that this type of epigenetic de-regulation of the MECP2 may be present in other RTT-like patients.

Transcriptional regulation of the mannosyltransferase-encoding gene pigm in inherited glycosylphosphatidylinositol (GPI) deficiency

RESUMO:O glicosilfosfatidilinositol (GPI) é um complexo glicolipídico utlizado por dezenas de proteínas, o qual medeia a sua ancoragem à superfície da célula. Proteínas de superfície celular ancoradas a GPI apresentam várias funções essenciais para a manutenção celular. A deficiência na síntese de GPI é o que caracteriza principalmente a deficiência hereditária em GPI, um grupo de doenças autossómicas raras que resultam de mutações nos genes PIGA, PIGL, PIGM, PIGV, PIGN, PIGO e PIGT, os quais sao indispensáveis para a biossíntese do GPI. Uma mutação pontual no motivo rico em GC -270 no promotor de PIGM impede a ligação do factor de transcrição (FT) Sp1 à sua sequência de reconhecimento, impondo a compactação da cromatina, associada à hipoacetilação de histonas, e consequentemente, impedindo a transcrição de PIGM. Desta forma, a adição da primeira manose ao GPI é comprometida, a síntese de GPI diminui assim como as proteínas ligadas a GPI à superficie das células. Pacientes com Deficiência Hereditária em GPI-associada a PIGM apresentam trombose e epilesia, e ausência de hemólise intravascular e anemia, sendo que estas duas últimas características definem a Hemoglobinúria Paroxística Nocturna (HPN), uma doença rara causada por mutações no gene PIGA. Embora a mutação que causa IGD seja constitutiva e esteja presente em todos os tecidos, o grau de deficiência em GPI varia entre células do mesmo tecido e entre células de tecidos diferentes. Por exemplo nos granulócitos e linfócitos B a deficiência em GPI é muito acentuada mas nos linfócitos T, fibroblastos, plaquetas e eritrócitos é aproximadamente normal, daí a ausência de hemólise intravascular. Os eventos transcricionais que estão na base da expressão diferencial da âncora GPI nas células hematopoiéticas são desconhecidos e constituem o objectivo geral desta tese. Em primeiro lugar, os resultados demonstraram que os níveis de PIGM mRNA variam entre células primárias hematopoiéticas normais. Adicionalmente, a configuração dos nucleossomas no promotor de PIGM é mais compacta em células B do que em células eritróides e tal está correlacionado com os níveis de expressão de PIGM, isto é, inferior nas células B. A presença de vários motivos de ligação para o FT específico da linhagem megacariocítica-eritróide GATA-1 no promotor de PIGM sugeriu que GATA-1 desempenha um papel regulador na sua transcrição. Os resultados mostraram que muito possivelmente GATA-1 desempenha um papel repressor em vez de activador da expressão de PIGM. Resultados preliminares sugerem que KLF1, um factor de transcrição restritamente eritróide, regula a transcrição de PIGM independentemente do motivo -270GC. Em segundo lugar, a investigação do papel dos FTs Sp demonstrou que Sp1 medeia directamente a transcrição de PIGM em ambas as células B e eritróide. Curiosamente, ao contrário do que acontece nas células B, em que a transcrição de PIGM requer a ligação do FT geral Sp1 ao motivo -270GC, nas células eritróides Sp1 regula a transcrição de PIGM ao ligar-se a montante e não ao motivo -270GC. Para além disso, demonstrou-se que Sp2 não é um regulador directo da transcrição de PIGM quer nas células B quer nas células eritróides. Estes resultados explicam a ausência de hemólise intravascular nos doentes com IGD associada a PIGM, uma das principais características que define a HPN. Por último, resultados preliminares mostraram que a repressão da transcrição de PIGM devida à mutação patogénica -270C>G está associada com a diminuição da frequência de interacções genómicas em cis entre PIGM e os seus genes “vizinhos”, sugerindo adicionalmente que a regulação de PIGM e desses genes é partilhada. No seu conjunto, os resultados apresentados nesta tese contribuem para o conhecimento do controlo transcricional de um gene housekeeping, específico-detecido, por meio de FTs genéricos e específicos de linhagem.-------------ABSTRACTC: Glycosylphosphatidylinositol (GPI) is a complex glycolipid used by dozens of proteins for cell surface anchoring. GPI-anchored proteins have various functions that are essential for the cellular maintenance. Defective GPI biosynthesis is the hallmark of inherited GPI deficiency (IGD), a group of rare autosomal diseases caused by mutations in PIGA, PIGL, PIGM, PIGV, PIGN, PIGO and PIGT, all genes indispensable for GPI biosynthesis. A point mutation in the -270GC-rich box in the core promoter of PIGM disrupts binding of the transcription factor (TF) Sp1 to it, imposing nucleosome compaction associated with histone hypoacetylation, thus abrogating transcription of PIGM. As a consequence of PIGM transcriptional repression, addition of the first mannose residue onto the GPI core and thus GPI production are impaired; and expression of GPI-anchored proteins on the surface of cells is severely impaired. Patients with PIGM-associated IGD suffer from life-threatening thrombosis and epilepsy but not intravascular haemolysis and anaemia, two defining features of paroxysmal nocturnal haemoglobinuria (PNH), a rare disease caused by somatic mutations in PIGA. Although the disease-causing mutation in IGD is constitutional and present in all tissues, the degree of GPI deficiency is variable and differs between cells of the same and of different tissues. Accordingly, GPI deficiency is severe in granulocytes and B cells but mild in T cells, fibroblasts, platelets and erythrocytes, hence the lack of intravascular haemolysis.The transcriptional events underlying differential expression of GPI in the haematopoietic cells of PIG-M-associated IGD are not known and constitute the general aim of this thesis. Firstly, I found that PIGM mRNA levels are variable amongst normal primary haematopoietic cells. In addition, the nucleosome configuration in the promoter of PIGM is more compacted in B cells than in erythroid cells and this correlated with the levels of PIGM mRNA expression, i.e., lower in B cells. The presence of several binding sites for GATA-1, a mega-erythroid lineage-specific transcription factor (TF), at the PIGM promoter suggested that GATA-1 has a role on PIGM transcription. My results showed that GATA-1 in erythroid cells is most likely a repressor rather than an activator of PIGM expression. Preliminary data suggested that KLF1, an erythroid-specific TF, regulates PIGM transcription but independently of the -270GC motif. Secondly, investigation of the role of the Sp TFs showed that Sp1 directly mediates PIGM transcriptional regulation in both B and erythroid cells. However, unlike in B cells in which active PIGM transcription requires binding of the generic TF Sp1 to the -270GC-rich box, in erythroid cells, Sp1 regulates PIGM transcription by binding upstream of but not to the -270GC-rich motif. Additionally, I showed that Sp2 is not a direct regulator of PIGM transcription in B and erythroid cells. These findings explain lack of intravascular haemolysis in PIGM-associated IGD, a defining feature of PNH. Lastly, preliminary work shows that transcriptional repression of PIG-M by the pathogenic -270C>G mutation is associated with reduced frequency of in cis genomic interactions between PIGM and its neighbouring genes, suggesting a shared regulatory link between these genes and PIGM. Altogether, the results presented in this thesis provide novel insights into tissuespecific transcriptional control of a housekeeping gene by lineage-specific and generic TFs.

Evolution of the gene regulatory network controlling flower dorsoventral asymmetry

Tese de Doutoramento em Biologia de Plantas.

Elf-1 contributes to the function of the complex interleukin (IL)-2-responsive enhancer in the mouse IL-2 receptor alpha gene.

Lymphocytes regulate their responsiveness to IL-2 through the transcriptional control of the IL-2R alpha gene, which encodes a component of the high affinity IL-2 receptor. In the mouse IL-2R alpha gene this control is exerted via two regulatable elements, a promoter proximal region, and an IL-2-responsive enhancer (IL-2rE) 1.3 kb upstream. In vitro and in vivo functional analysis of the IL-2rE in the rodent thymic lymphoma-derived, CD4- CD8- cell line PC60 demonstrated that three separate elements, sites I, II, and III, were necessary for IL-2 responsiveness; these three sites demonstrate functional cooperation. Site III contains a consensus binding motif for members of the Ets family of transcription factors. Here we demonstrate that Elf-1, an Ets-like protein, binds to site III and participates in IL-2 responsiveness. In vitro site III forms a complex with a protein constitutively present in nuclear extracts from PC60 cells as well as from normal CD4- CD8- thymocytes. We have identified this molecule as Elf-1 according to a number of criteria. The complex possesses an identical electrophoretic mobility to that formed by recombinant Elf-1 protein and is super-shifted by anti-Elf-1 antibodies. Biotinylated IL-2rE probes precipitate Elf-1 from PC60 extracts provided site III is intact and both recombinant and PC60-derived proteins bind with the same relative affinities to different mutants of site III. In addition, by introducing mutations into the core of the site III Ets-like motif and comparing the corresponding effects on the in vitro binding of Elf-1 and the in vivo IL-2rE activity, we provide strong evidence that Elf-1 is directly involved in IL-2 responsiveness. The nature of the functional cooperativity observed between Elf-1 and the factors binding sites I and II remains unresolved; experiments presented here however suggest that this effect may not require direct interactions between the proteins binding these three elements.

Estrogen receptor regulates MyoD gene expression by preventing AP-1-mediated repression.

Cell growth and differentiation are opposite events in the myogenic lineage. Growth factors block the muscle differentiation program by inducing the expression of transcription factors that negatively regulate the expression of muscle regulatory genes like MyoD. In contrast, extracellular clues that induce cell cycle arrest promote MyoD expression and muscle differentiation. Thus, the regulation of MyoD expression is critical for muscle differentiation. Here we show that estrogen induces MyoD expression in mouse skeletal muscle in vivo and in dividing myoblasts in vitro by relieving the MyoD promoter from AP-1 negative regulation through a mechanism involving estrogen receptor/AP-1 protein-protein interactions but independent of the estrogen receptor DNA binding activity.

Pro-opiomelanocortin gene and melanin-based colour polymorphism in a reptile

Colour polymorphism is widespread among vertebrates and plays important roles in prey-predator interactions, thermoregulation, social competition, and sexual selection. However, the genetic mechanisms involved in colour variation have been studied mainly in domestic mammals and birds, whereas information on wild animals remains scarce. Interestingly, the pro-opiomelanocortin gene (POMC) gives rise to melanocortin hormones that trigger melanogenesis (by binding the melanocortin-1-receptor; Mc1r) and other physiological and behavioural functions (by binding the melanocortin receptors Mc1-5rs). Owing to its pleiotropic effect, the POMC gene could therefore account for the numerous covariations between pigmentation and other phenotypic traits. We screened the POMC and Mc1r genes in 107 wild asp vipers (Vipera aspis) that can exhibit four discrete colour morphs (two unpatterned morphs: concolor or melanistic; two patterned morphs: blotched or lined) in a single population. Our study revealed a correlation between a single nucleotide polymorphism situated within the 3-untranslated region of the POMC gene and colour variation, whereas Mc1r was not found to be polymorphic. To the best of our knowledge, we disclose for the first time a relationship between a mutation at the POMC gene and coloration in a wild animal, as well as a correlation between a genetic marker and coloration in a snake species. Interestingly, similar mutations within the POMC 3-untranslated region are linked to human obesity and alcohol and drug dependence. Combined with our results, this suggests that the 3-untranslated region of the POMC gene may play a role in its regulation in distant vertebrates.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Protein-Binding Microarray Analysis of Tumor Suppressor AP2α Target Gene Specificity.

Cheap and massively parallel methods to assess the DNA-binding specificity of transcription factors are actively sought, given their prominent regulatory role in cellular processes and diseases. Here we evaluated the use of protein-binding microarrays (PBM) to probe the association of the tumor suppressor AP2α with 6000 human genomic DNA regulatory sequences. We show that the PBM provides accurate relative binding affinities when compared to quantitative surface plasmon resonance assays. A PBM-based study of human healthy and breast tumor tissue extracts allowed the identification of previously unknown AP2α target genes and it revealed genes whose direct or indirect interactions with AP2α are affected in the diseased tissues. AP2α binding and regulation was confirmed experimentally in human carcinoma cells for novel target genes involved in tumor progression and resistance to chemotherapeutics, providing a molecular interpretation of AP2α role in cancer chemoresistance. Overall, we conclude that this approach provides quantitative and accurate assays of the specificity and activity of tumor suppressor and oncogenic proteins in clinical samples, interfacing genomic and proteomic assays.

cAMP prevents the glucose-mediated stimulation of GLUT2 gene transcription in hepatocytes.

Glucose homoeostasis necessitates the presence in the liver of the high Km glucose transporter GLUT2. In hepatocytes, we and others have demonstrated that glucose stimulates GLUT2 gene expression in vivo and in vitro. This effect is transcriptionally regulated and requires glucose metabolism within the hepatocytes. In this report, we further characterized the cis-elements of the murine GLUT2 promoter, which confers glucose responsiveness on a reporter gene coding the chloramphenicol acetyl transferase (CAT) gene. 5'-Deletions of the murine GLUT2 promoter linked to the CAT reporter gene were transfected into a GLUT2 expressing hepatoma cell line (mhAT3F) and into primary cultured rat hepatocytes, and subsequently incubated at low and high glucose concentrations. Glucose stimulates gene transcription in a manner similar to that observed for the endogenous GLUT2 mRNA in both cell types; the -1308 to -212 bp region of the promoter contains the glucose-responsive elements. Furthermore, the -1308 to -338 bp region of the promoter contains repressor elements when tested in an heterologous thymidine kinase promoter. The glucose-induced GLUT2 mRNA accumulation was decreased by dibutyryl-cAMP both in mhAT3F cells and in primary hepatocytes. A putative cAMP-responsive element (CRE) is localized at the -1074/-1068 bp region of the promoter. The inhibitory effect of cAMP on GLUT2 gene expression was observed in hepatocytes transfected with constructs containing this CRE (-1308/+49 bp fragment), as well as with constructs not containing the consensus CRE (-312/+49 bp fragment). This suggests that the inhibitory effect of cAMP is not mediated by the putative binding site located in the repressor fragment of the GLUT2 promoter. Taken together, these data demonstrate that the elements conferring glucose and cAMP responsiveness on the GLUT2 gene are located within the -312/+49 region of the promoter.