Genetic determinants of clinical heterogeneity in sickle cell disease

L’anémie falciforme est une maladie monogénique causée par une mutation dans le locus de la β-globine. Malgré le fait que l’anémie falciforme soit une maladie monogénique, cette maladie présente une grande hétérogénéité clinique. On présume que des facteurs environnementaux et génétiques contribuent à cette hétérogénéité. Il a été observé qu’un haut taux d’hémoglobine fœtale (HbF) diminuait la sévérité et la mortalité des patients atteints de l’anémie falciforme. Le but de mon projet était d’identifier des variations génétiques modifiant la sévérité clinique de l’anémie falciforme. Dans un premier temps, nous avons effectué la cartographie-fine de trois régions précédemment associées avec le taux d’hémoglobine fœtale. Nous avons ensuite effectué des études d’association pan-génomiques avec deux complications cliniques de l’anémie falciforme ainsi qu’avec le taux d’hémoglobine fœtale. Hormis les régions déjà identifiées comme étant associées au taux d’hémoglobine fœtale, aucun locus n’a atteint le niveau significatif de la puce de génotypage. Pour identifier des groupes de gènes modérément associés au taux d’hémoglobine fœtale qui seraient impliqués dans de mêmes voies biologiques, nous avons effectué une étude des processus biologiques. Finalement, nous avons effectué l’analyse de 19 exomes de patients Jamaïcains ayant des complications cliniques mineures de l’anémie falciforme. Compte tenu de la taille des cohortes de réplication disponibles, nous n’avons pas les moyens de valider statistiquement les variations identifiées par notre étude. Cependant, nos résultats fournissent de bons gènes candidats pour des études fonctionnelles et pour les réplications futures. Nos résultats suggèrent aussi que le β-hydroxybutyrate en concentration endogène pourraient influencer le taux d’hémoglobine fœtale. De plus, nous montrons que la cartographie-fine des régions associées par des études pan-génomiques peut identifier des signaux d’association additionnels et augmenter la variation héritable expliquée par cette région.

Hemorheological alterations in sickle cell anemia and their clinical consequences: the role of genetic modulators

Sickle cell anemia (SCA) is an autosomal recessive chronic hemolytic anemia, caused by homozygosity for the HBB:c.20A>T mutation. The disease presents with high clinical heterogeneity, stroke being the most devastating manifestation. This study aimed to identify genetic modulators of severe hemolysis and stroke risk in children with SCA, as well as understand their consequences at the hemorheological level. Sixty-six children with SCA were categorised according to their degree of cerebral vasculopathy (Stroke/Risk/Control). Relevant data were collected from patients’ medical records. Several polymorphic regions in genes related to vascular cell adhesion and tonus were characterized by molecular methodologies. Data analyses were performed using R software. Several in silico tools (e.g. TFBind, MatInspector) were applied to investigate the main variant consequences. Some genetic variants in vascular adhesion molecule-1 gene promoter and endothelial nitric oxide synthase gene were associated with higher levels of hemolysis and stroke events. They modify important transcription factor binding sites or disturb the corresponding protein structure/function. Our findings emphasize the relevance of the genetic variants in modulating the degree of hemolysis and development of cerebral vasculopathy due to their effect on gene expression, modification of protein biological activities related with erythrocyte/endothelial interactions and consequent hemorheological abnormalities in SCA.

Understanding the genetic basis of phenotype variability in individuals with neurocognitive disorders

Thesis (Ph.D.)--University of Washington, 2016-06

Evidence for a modifier of onset age in Huntington disease linked to the HD gene in 4p16

Huntington disease (HD) is a neurodegenerative disorder caused by the abnormal expansion of CAG repeats in the HD gene on chromosome 4p16.3. A recent genome scan for genetic modifiers of age at onset of motor symptoms (AO) in HD suggests that one modifier may reside in the region close to the HD gene itself. We used data from 535 HD participants of the New England Huntington cohort and the HD MAPS cohort to assess whether AO was influenced by any of the three markers in the 4p16 region: MSX1 (Drosophila homeo box homologue 1, formerly known as homeo box 7, HOX7), Delta2642 (within the HD coding sequence), and BJ56 (D4S127). Suggestive evidence for an association was seen between MSX1 alleles and AO, after adjustment for normal CAG repeat, expanded repeat, and their product term (model P value 0.079). Of the variance of AO that was not accounted for by HD and normal CAG repeats, 0.8% could be attributed to the MSX1 genotype. Individuals with MSX1 genotype 3/3 tended to have younger AO. No association was found between Delta2642 (P=0.44) and BJ56 (P=0.73) and AO. This study supports previous studies suggesting that there may be a significant genetic modifier for AO in HD in the 4p16 region. Furthermore, the modifier may be present on both HD and normal chromosomes bearing the 3 allele of the MSX1 marker.

Ancestral haplotype 8.1 and lung disease severity in European cystic fibrosis patients

The clinical course of cystic fibrosis (CF) lung disease varies between patients bearing identical CFTR mutations. This suggests that additional genetic modifiers may contribute to the pulmonary phenotype. The highly conserved ancestral haplotype 8.1 (8.1AH), carried by up to one quarter of Caucasians, comprises linked gene polymorphisms on chromosome 6 that play a key role in the inflammatory response: LTA +252A/G; TNF -308G/A, HSP70-2 +1267A/G and RAGE -429T/C. As inflammation is a key component inducing CF lung damage, we investigated whether the 8.1AH represents a lung function modifier in CF.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

Introduction: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

Methods: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

Results: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

Conclusions: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Copy number variation and Huntington’s disease

Tese de doutoramento, Ciências Biomédicas (Neurociências), Universidade de Lisboa, Faculdade de Medicina, 2014

Systems Genetics of Metabolism: The Use of the BXD Murine Reference Panel for Multiscalar Integration of Traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.

Cytokine gene variants and venous thrombotic risk in the BRATROS (BRAZILIAN THROMBOSIS STUDY)

Introduction: Venous thrombosis (VT) and inflammation are two closely related entities. In the present investigation we assessed whether there is a relation between genetic modifiers of the inflammatory response and the risk of VT. Materials and methods: 420 consecutive and unrelated patients with an objective diagnosis of deep VT and 420 matched controls were investigated. The frequencies of the following gene polymorphisms were determined in all subjects: TNF-α- 308 G/A, LT-α+ 252 A/G, IL-6-174 G/C, IL1-ra 86 bp VNTR, IL-10-1082 A/G and CD-31 125 C/G. Results: Overall odds ratio (OR) for VT related to TNF-α- 308 G/A, LT-α+ 252 A/G, IL-6-174 G/C, A1 allele (4 bp repeat) of the IL1-ra 86 bp VNTR, IL-10-1082 A/G and CD-31 125 C/G were respectively: 1.0 (CI95: 0.8-1.5), 1.3 (CI95: 1.0-1.7), 1.1 (CI95: 0.9-1.5), 1.6 (CI95: 1-2.5), 1.2 (CI95: 0.8-1.7) and 0.8 (CI95: 0.6-1.1). A possible interaction between polymorphisms was observed only for the co-inheritance of the mutant alleles of the LT-α+ 252 A/G and IL-10-1082 G/A polymorphisms (OR = 2; CI95: 1.1-3.8). The risk of VT conferred by factor V Leiden and FII G20210A was not substantially altered by co-inheritance with any of the cytokine gene polymorphisms. Conclusions: Cytokine gene polymorphisms here investigated did not significantly influence venous thrombotic risk. © 2006 Elsevier Ltd. All rights reserved.

Bacterial colitis increases susceptibility to oral prion disease

Dietary exposure to prion-contaminated materials has caused kuru and variant Creutzfeldt-Jakob disease in humans and transmissible spongiform encephalopathies (TSEs) in cattle, mink, and felines. The epidemiology of dietary prion infections suggests that host genetic modifiers and possibly exogenous cofactors may play a decisive role in determining disease susceptibility. However, few cofactors influencing susceptibility to prion infection have been identified. In the present study, we investigated whether colitis might represent one such cofactor. We report that moderate colitis caused by an attenuated Salmonella strain more than doubles the susceptibility of mice to oral prion infection and modestly accelerates the development of disease after prion challenge. The prion protein was up-regulated in intestines and mesenteric lymph nodes of mice with colitis, providing a possible mechanism for the effect of colitis on the pathogenesis of prion disease. Therefore, moderate intestinal inflammation at the time of prion exposure may constitute one of the elusive risk factors underlying the development of TSE.

Use of a hypomorphic allele of myogenin to analyze Myogenin-dependent processes in mouse development

Myogenin is a muscle-specific transcription factor essential for skeletal muscle differentiation. A severe reduction in the number of fused myotubes is seen in myogenin-null mice, and the expression of genes characteristic of differentiated skeletal muscle is reduced. Additionally, sternebrae defects are seen in myogenin-null mice, a secondary defect in the sternal cartilage precursors. Very little is known about the quantitative requirement for myogenin in muscle differentiation and thoracic skeletal development in vivo. In this thesis I describe experiments utilizing a mouse line harboring a hypomorphic allele of myogenin, generated by gene targeting techniques in embryonic stem cells. The nature of the hypomorphism was due to lowered levels of myogenin from this allele. In embryos homozygous for the hypomorphic allele, normal sternum formation and extensive muscle differentiation was observed. However, muscle hypoplasia and reduced muscle-specific gene expression were apparent in these embryos, and the mice were not viable after birth. These results suggest skeletal muscle differentiation is highly sensitive to the absolute amounts of myogenin, and reveal distinct threshold requirements for myogenin in skeletal muscle differentiation, sternum formation, and viability in vivo. The hypomorphic allele was utilized as a genetically sensitized background to identify other components of myogenin-mediated processes. Using a candidate gene approach I crossed null mutations in MEF2C and MRF4 into the hypomorphic background and examined whether these mutations affected muscle differentiation and skeleton formation in the myogenin hypomorph. Although MEF2C mutation did not affect any phenotypes seen in the hypomorphic background, MRF4 was observed to be an essential component of myogenin-mediated processes of thoracic skeletal development. Additionally, the hypomorphic allele was very sensitive to genetic effects, suggesting the existence of mappable genetic modifiers of the hypomorphic allele of myogenin. ^

Tumor burden and clonality in multiple intestinal neoplasia mouse/normal mouse aggregation chimeras

Aggregation chimeras were formed between C57BL/6 mice heterozygous for the Apcmin (Min) mutation and wild-type SWR mice, that differ in their Pla2g2a status, a modifier of Apcmin, and also in their resistance to intestinal polyp formation. Variation in the dolichos biflorus agglutinin-staining patterns of the intestines of these mouse strains was used to determine the chimeric composition of the intestine in individual mice and to examine the clonal composition of adenomas. Macroscopic adenoma numbers in chimeric mice were compared with the expected adenoma numbers based on the percentage of C57BL/6J-Apcmin/+ epithelium in individual mice. These results unexpectedly show that there was no apparent inhibitory effect of the SWR-derived (Pla2g2a wild-type) tissue on adenoma formation in the C57BL/6J-Apcmin/+ epithelium. This suggests that the main genetic modifiers of the Min phenotype act at a cellular or crypt-restricted level with no discernable systemic effect. All adenomas were seen to contain C57BL/6J-Apcmin/+-derived epithelium, confirming that the germ-line mutation of the mApc gene is necessary to initiate tumorigenesis in this model system, and that the mApc gene acts in a cell autonomous fashion.

A subset of skin tumor modifier loci determines survival time of tumor-bearing mice

Studies of mouse models of human cancer have established the existence of multiple tumor modifiers that influence parameters of cancer susceptibility such as tumor multiplicity, tumor size, or the probability of malignant progression. We have carried out an analysis of skin tumor susceptibility in interspecific Mus musculus/Mus spretus hybrid mice and have identified another seven loci showing either significant (six loci) or suggestive (one locus) linkage to tumor susceptibility or resistance. A specific search was carried out for skin tumor modifier loci associated with time of survival after development of a malignant tumor. A combination of resistance alleles at three markers [D6Mit15 (Skts12), D7Mit12 (Skts2), and D17Mit7 (Skts10)], all of which are close to or the same as loci associated with carcinoma incidence and/or papilloma multiplicity, is significantly associated with increased survival of mice with carcinomas, whereas the reverse combination of susceptibility alleles is significantly linked to early mortality caused by rapid carcinoma growth (χ2 = 25.22; P = 5.1 × 10−8). These data indicate that host genetic factors may be used to predict carcinoma growth rate and/or survival of individual backcross mice exposed to the same carcinogenic stimulus and suggest that mouse models may provide an approach to the identification of genetic modifiers of cancer survival in humans.

The fourth chromosome of Drosophila melanogaster: Interspersed euchromatic and heterochromatic domains

The small fourth chromosome of Drosophila melanogaster (3.5% of the genome) presents a puzzle. Cytological analysis suggests that the bulk of the fourth, including the portion that appears banded in the polytene chromosomes, is heterochromatic; the banded region includes blocks of middle repetitious DNA associated with heterochromatin protein 1 (HP1). However, genetic screens indicate 50–75 genes in this region, a density similar to that in other euchromatic portions of the genome. Using a P element containing an hsp70-white gene and a copy of hsp26 (marked with a fragment of plant DNA designated pt), we have identified domains that allow for full expression of the white marker (R domains), and others that induce a variegating phenotype (V domains). In the former case, the hsp26-pt gene shows an accessibility and heat-shock-inducible activity similar to that seen in euchromatin, whereas in the latter case, accessibility and inducible expression are reduced to levels typical of heterochromatin. Mapping by in situ hybridization and by hybridization of flanking DNA sequences to a collection of cosmid and bacterial artificial chromosome clones shows that the R domains (euchromatin-like) and V domains (heterochromatin-like) are interspersed. Examination of the effect of genetic modifiers on the variegating transgenes shows some differences among these domains. The results suggest that heterochromatic and euchromatic domains are interspersed and closely associated within this 1.2-megabase region of the genome.