Highly parallel genetic approaches in Mendelian and complex diseases

The recent advance in high-throughput sequencing and genotyping protocols allows rapid investigation of Mendelian and complex diseases on a scale not previously been possible. In my thesis research I took advantage of these modern techniques to study retinitis pigmentosa (RP), a rare inherited disease characterized by progressive loss of photoreceptors and leading to blindness; and hypertension, a common condition affecting 30% of the adult population. Firstly, I compared the performance of different next generation sequencing (NGS) platforms in the sequencing of the RP-linked gene PRPF31. The gene contained a mutation in an intronic repetitive element, which presented difficulties for both classic sequencing methods and NGS. We showed that all NGS platforms are powerful tools to identify rare and common DNA variants, also in case of more complex sequences. Moreover, we evaluated the features of different NGS platforms that are important in re-sequencing projects. The main focus of my thesis was then to investigate the involvement of pre-mRNA splicing factors in autosomal dominant RP (adRP). I screened 5 candidate genes in a large cohort of patients by using long-range PCR as enrichment step, followed by NGS. We tested two different approaches: in one, all target PCRs from all patients were pooled and sequenced as a single DNA library; in the other, PCRs from each patient were separated within the pool by DNA barcodes. The first solution was more cost-effective, while the second one allowed obtaining faster and more accurate results, but overall they both proved to be effective strategies for gene screenings in many samples. We could in fact identify novel missense mutations in the SNRNP200 gene, encoding an essential RNA helicase for splicing catalysis. Interestingly, one of these mutations showed incomplete penetrance in one family with adRP. Thus, we started to study the possible molecular causes underlying phenotypic differences between asymptomatic and affected members of this family. For the study of hypertension, I joined a European consortium to perform genome-wide association studies (GWAS). Thanks to the use of very informative genotyping arrays and of phenotipically well-characterized cohorts, we could identify a novel susceptibility locus for hypertension in the promoter region of the endothelial nitric oxide synthase gene (NOS3). Moreover, we have proven the direct causality of the associated SNP using three different methods: 1) targeted resequencing, 2) luciferase assay, and 3) population study. - Le récent progrès dans le Séquençage à haut Débit et les protocoles de génotypage a permis une plus vaste et rapide étude des maladies mendéliennes et multifactorielles à une échelle encore jamais atteinte. Durant ma thèse de recherche, j'ai utilisé ces nouvelles techniques de séquençage afin d'étudier la retinite pigmentale (RP), une maladie héréditaire rare caractérisée par une perte progressive des photorécepteurs de l'oeil qui entraine la cécité; et l'hypertension, une maladie commune touchant 30% de la population adulte. Tout d'abord, j'ai effectué une comparaison des performances de différentes plateformes de séquençage NGS (Next Generation Sequencing) lors du séquençage de PRPF31, un gène lié à RP. Ce gène contenait une mutation dans un élément répétable intronique, qui présentait des difficultés de séquençage avec la méthode classique et les NGS. Nous avons montré que les plateformes de NGS analysées sont des outils très puissants pour identifier des variations de l'ADN rares ou communes et aussi dans le cas de séquences complexes. De plus, nous avons exploré les caractéristiques des différentes plateformes NGS qui sont importantes dans les projets de re-séquençage. L'objectif principal de ma thèse a été ensuite d'examiner l'effet des facteurs d'épissage de pre-ARNm dans une forme autosomale dominante de RP (adRP). Un screening de 5 gènes candidats issus d'une large cohorte de patients a été effectué en utilisant la long-range PCR comme étape d'enrichissement, suivie par séquençage avec NGS. Nous avons testé deux approches différentes : dans la première, toutes les cibles PCRs de tous les patients ont été regroupées et séquencées comme une bibliothèque d'ADN unique; dans la seconde, les PCRs de chaque patient ont été séparées par code barres d'ADN. La première solution a été la plus économique, tandis que la seconde a permis d'obtenir des résultats plus rapides et précis. Dans l'ensemble, ces deux stratégies se sont démontrées efficaces pour le screening de gènes issus de divers échantillons. Nous avons pu identifier des nouvelles mutations faux-sens dans le gène SNRNP200, une hélicase ayant une fonction essentielle dans l'épissage. Il est intéressant de noter qu'une des ces mutations montre une pénétrance incomplète dans une famille atteinte d'adRP. Ainsi, nous avons commencé une étude sur les causes moléculaires entrainant des différences phénotypiques entre membres affectés et asymptomatiques de cette famille. Lors de l'étude de l'hypertension, j'ai rejoint un consortium européen pour réaliser une étude d'association Pangénomique ou genome-wide association study Grâce à l'utilisation de tableaux de génotypage très informatifs et de cohortes extrêmement bien caractérisées au niveau phénotypique, un nouveau locus lié à l'hypertension a été identifié dans la région promotrice du gène endothélial nitric oxide sinthase (NOS3). Par ailleurs, nous avons prouvé la cause directe du SNP associé au moyen de trois méthodes différentes: i) en reséquençant la cible avec NGS, ii) avec des essais à la luciférase et iii) une étude de population.

Genome-wide analysis of cancer/testis gene expression.

Cancer/Testis (CT) genes, normally expressed in germ line cells but also activated in a wide range of cancer types, often encode antigens that are immunogenic in cancer patients, and present potential for use as biomarkers and targets for immunotherapy. Using multiple in silico gene expression analysis technologies, including twice the number of expressed sequence tags used in previous studies, we have performed a comprehensive genome-wide survey of expression for a set of 153 previously described CT genes in normal and cancer expression libraries. We find that although they are generally highly expressed in testis, these genes exhibit heterogeneous gene expression profiles, allowing their classification into testis-restricted (39), testis/brain-restricted (14), and a testis-selective (85) group of genes that show additional expression in somatic tissues. The chromosomal distribution of these genes confirmed the previously observed dominance of X chromosome location, with CT-X genes being significantly more testis-restricted than non-X CT. Applying this core classification in a genome-wide survey we identified >30 CT candidate genes; 3 of them, PEPP-2, OTOA, and AKAP4, were confirmed as testis-restricted or testis-selective using RT-PCR, with variable expression frequencies observed in a panel of cancer cell lines. Our classification provides an objective ranking for potential CT genes, which is useful in guiding further identification and characterization of these potentially important diagnostic and therapeutic targets.

Genome-wide association study of co-occurring anxiety in major depression.

OBJECTIVES: Co-morbidity between depression and anxiety disorders is common. In this study we define a quantitative measure of anxiety by summating four anxiety items from the SCAN interview in a large collection of major depression (MDD) cases to identify genes contributing to this complex phenotype. METHODS: A total of 1522 MDD cases dichotomised according to those with at least one anxiety item scored (n = 1080) and those without anxiety (n = 442) were analysed, and also compared to 1588 healthy controls at a genome-wide level, to identify genes that may contribute to anxiety in MDD. RESULTS: For the quantitative trait, suggestive evidence of association was detected for two SNPs, and for the dichotomous anxiety present/absent ratings for three SNPs at genome-wide level. In the genome-wide analysis of MDD cases with co-morbid anxiety and healthy controls, two SNPs attained P values of < 5 × 10⁻⁶. Analysing candidate genes, P values ≤ 0.0005 were found with three SNPs for the quantitative trait and three SNPs for the dichotomous trait. CONCLUSIONS: This study provides an initial genome-wide assessment of possible genetic contribution to anxiety in MDD. Although suggestive evidence of association was found for several SNPs, our findings suggest that there are no common variants strongly associated with anxious depression.

Population pharmacokinetic analysis and pharmacogenetics of raltegravir in HIV-positive and healthy individuals.

The objectives of this study were to characterize raltegravir (RAL) population pharmacokinetics in HIV-positive (HIV(+)) and healthy individuals, identify influential factors, and search for new candidate genes involved in UDP glucuronosyltransferase (UGT)-mediated glucuronidation. The pharmacokinetic analysis was performed with NONMEM. Genetic association analysis was performed with PLINK using the relative bioavailability as the phenotype. Simulations were performed to compare once- and twice-daily regimens. A 2-compartment model with first-order absorption adequately described the data. Atazanavir, gender, and bilirubin levels influenced RAL relative bioavailability, which was 30% lower in HIV(+) than in healthy individuals. UGT1A9*3 was the only genetic variant possibly influencing RAL pharmacokinetics. The majority of RAL pharmacokinetic variability remains unexplained by genetic and nongenetic factors. Owing to the very large variability, trough drug levels might be very low under the standard dosing regimen, raising the question of a potential relevance of therapeutic drug monitoring of RAL in some situations.

A polymorphism in the leptin gene promoter is associated with anemia in patients with HIV disease.

To study factors associated with anemia and its effect on survival in HIV-infected persons treated with modern combined antiretroviral therapy (cART), we characterized the prevalence of anemia in the Veterans Aging Cohort Study (VACS) and used a candidate gene approach to identify proinflammatory gene single nucleotide polymorphisms (SNPs) associated with anemia in HIV disease. The study comprised 1597 HIV(+) and 865 HIV(-) VACS subjects with DNA, blood, and annotated clinical data available for analysis. Anemia was defined according to World Health Organization criteria (hemoglobin < 13 g/dL and < 12 g/dL in men and women, respectively). The prevalence of anemia in HIV(+) and HIV(-) subjects was 23.1% and 12.9%, respectively. Independent of HIV status, anemia was present in 23.4% and 8% in blacks and whites, respectively. Analysis of our candidate genes revealed that the leptin -2548 G/A SNP was associated with anemia in HIV(+), but not HIV(-), patients, with the AA and AG genotypes significantly predicting anemia (P < .003 and P < .039, respectively, logistic regression). This association was replicated in an independent cohort of HIV(+) women. Our study provides novel insight into the association between genetic variability in the leptin gene and anemia in HIV(+) individuals.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.

Genetics of essential tremor.

Essential tremor (ET) is a prevalent condition manifesting with progressive action tremor. Although ET was traditionally viewed as a sporadic disease, a significant proportion of cases report a positive family history of tremor. Autosomal dominant inheritance can be demonstrated in many families. Previously, genome-wide linkage studies in families mapped three loci for ET, hereditary essential tremor-1 (ETM1), ETM2 and ETM3. However, no causal mutation has been replicated in candidate genes within these loci, including dopamine D3 receptor (DRD3) and HS1-binding protein 3 (HS1BP3). Recently, the first genome-wide association study in ET followed by replication studies conducted in diverse populations identified a significant association between the leucine-rich repeat and Ig domain containing 1 gene (LINGO1) SNP rs9652490 and risk for ET Although further novel variants were indentified in LINGO1 and its paralog LINGO2 that may be associated with risk for ET, the pathogenic mechanisms involved remain elusive. Given the possibility that ET as a complex trait may be influenced by the combined effects of rare variants, novel high-throughput technologies sequencing all exons across the genome (exome sequencing) or the whole genome (genome sequencing) may become crucial in understanding/deciphering the genetic background of ET.

IROme, a New High-Throughput Molecular Tool for the Diagnosis of Inherited Retinal Dystrophies-A Price Comparison with Sanger Sequencing.

The molecular diagnosis of retinal dystrophies (RD) is difficult because of genetic and clinical heterogeneity. Previously, the molecular screening of genes was done one by one, sometimes in a scheme based on the frequency of sequence variants and the number of exons/length of the candidate genes. Payment for these procedures was complicated and the sequential billing of several genes created endless paperwork. We therefore evaluated the costs of generating and sequencing a hybridization-based DNA library enriched for the 64 most frequently mutated genes in RD, called IROme, and compared them to the costs of amplifying and sequencing these genes by the Sanger method. The production cost generated by the high-throughput (HT) sequencing of IROme was established at CHF 2,875.75 per case. Sanger sequencing of the same exons cost CHF 69,399.02. Turnaround time of the analysis was 3 days for IROme. For Sanger sequencing, it could only be estimated, as we never sequenced all 64 genes in one single patient. Sale cost for IROme calculated on the basis of the sale cost of one exon by Sanger sequencing is CHF 8,445.88, which corresponds to the sale price of 40 exons. In conclusion, IROme is cheaper and faster than Sanger sequencing and therefore represents a sound approach for the diagnosis of RD, both scientifically and economically. As a drop in the costs of HT sequencing is anticipated, target resequencing might become the new gold standard in the molecular diagnosis of RD.

Genome-wide association analyses identify 18 new loci associated with serum urate concentrations.

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

Genetic heterogeneity of the Coppock-like cataract: a mutation in CRYBB2 on chromosome 22q11.2.

PURPOSE: To identify the genetic defect for the Coppock-like cataract (CCL) affecting a Swiss family, which defect was unlinked to the chromosome 2q33-35 CCL locus. METHODS: A large family was characterized for linkage analysis by slit lamp examination or by the review of drawings made before cataract extraction. The affection status was attributed before genotyping, and the genotyping was masked to the affection status. Two-point and multipoint linkage analyses were performed using the MLINK and the LINKMAP components of the LINKAGE program package (ver. 5.1), respectively. Mutational analysis of candidate genes was performed by a combination of direct cycle sequencing and an amplification refractory mutation system assay. RESULTS: Ten individuals were affected with the CCL phenotype. The disease was autosomal dominant and appeared to be fully penetrant. A new CCL locus was identified on chromosome 22q11.2 within a 11.67-cM interval (maximum lod score [Zmax] = 4.14; theta = 0). Mutational analysis of the CRYBB2 candidate gene identified a disease-causing mutation in exon 6. This sequence change was identical with that previously described to be associated with the cerulean cataract, a clinically distinct entity. CONCLUSIONS: The CCL phenotype is genetically heterogeneous with a second gene on chromosome 22q11.2, CRYBB2. The CCL and the cerulean cataract are two distinct clinical entities associated with the same genetic defect. This work provides evidence for a modifier factor that influences cataract formation and that remains to be identified.

Proximal 11p deletion syndrome (P11pDS): additional evaluation of the clinical and molecular aspects.

The combination of multiple exostoses (EXT) and enlarged parietal foramina (foramina parietalia permagna, FPP) represent the main features of the proximal 11p deletion syndrome (P11pDS), a contiguous gene syndrome (MIM 601224) caused by an interstitial deletion on the short arm of chromosome 11. Here we present clinical aspects of two new P11pDS patients and the clinical follow-up of one patient reported in the original paper describing this syndrome. Recognised clinical signs include EXT, FPP, mental retardation, facial asymmetry, asymmetric calcification of coronary sutures, defective vision (severe myopia, nystagmus, strabismus), skeletal anomalies (small hands and feet, tapering fingers), heart defect, and anal stenosis. In addition fluorescence in situ hybridisation and molecular analysis were performed to gain further insight in potential candidate genes involved in P11pDS.

Fanconi-Bickel syndrome and autosomal recessive proximal tubulopathy with hypercalciuria (ARPTH) are allelic variants caused by GLUT2 mutations.

CONTEXT: Many inherited disorders of calcium and phosphate homeostasis are unexplained at the molecular level. OBJECTIVE: The objective of the study was to identify the molecular basis of phosphate and calcium abnormalities in two unrelated, consanguineous families. PATIENTS: The affected members in family 1 presented with rickets due to profound urinary phosphate-wasting and hypophosphatemic rickets. In the previously reported family 2, patients presented with proximal renal tubulopathy and hypercalciuria yet normal or only mildly increased urinary phosphate excretion. METHODS: Genome-wide linkage scans and direct nucleotide sequence analyses of candidate genes were performed. Transport of glucose and phosphate by glucose transporter 2 (GLUT2) was assessed using Xenopus oocytes. Renal sodium-phosphate cotransporter 2a and 2c (Npt2a and Npt2c) expressions were evaluated in transgenically rescued Glut2-null mice (tgGlut2-/-). RESULTS: In both families, genetic mapping and sequence analysis of candidate genes led to the identification of two novel homozygous mutations (IVS4-2A>G and R124S, respectively) in GLUT2, the gene mutated in Fanconi-Bickel syndrome, a rare disease usually characterized by renal tubulopathy, impaired glucose homeostasis, and hepatomegaly. Xenopus oocytes expressing the [R124S]GLUT2 mutant showed a significant reduction in glucose transport, but neither wild-type nor mutant GLUT2 facilitated phosphate import or export; tgGlut2-/- mice demonstrated a profound reduction of Npt2c expression in the proximal renal tubules. CONCLUSIONS: Homozygous mutations in the facilitative glucose transporter GLUT2, which cause Fanconi-Bickel syndrome, can lead to very different clinical and biochemical findings that are not limited to mild proximal renal tubulopathy but can include significant hypercalciuria and highly variable degrees of urinary phosphate-wasting and hypophosphatemia, possibly because of the impaired proximal tubular expression of Npt2c.

Autosomal-dominant distal myopathy associated with a recurrent missense mutation in the gene encoding the nuclear matrix protein, matrin 3.

Distal myopathies represent a heterogeneous group of inherited skeletal muscle disorders. One type of adult-onset, progressive autosomal-dominant distal myopathy, frequently associated with dysphagia and dysphonia (vocal cord and pharyngeal weakness with distal myopathy [VCPDM]), has been mapped to chromosome 5q31 in a North American pedigree. Here, we report the identification of a second large VCPDM family of Bulgarian descent and fine mapping of the critical interval. Sequencing of positional candidate genes revealed precisely the same nonconservative S85C missense mutation affecting an interspecies conserved residue in the MATR3 gene in both families. MATR3 is expressed in skeletal muscle and encodes matrin 3, a component of the nuclear matrix, which is a proteinaceous network that extends throughout the nucleus. Different disease related haplotype signatures in the two families provided evidence that two independent mutational events at the same position in MATR3 cause VCPDM. Our data establish proof of principle that the nuclear matrix is crucial for normal skeletal muscle structure and function and put VCPDM on the growing list of monogenic disorders associated with the nuclear proteome.

Mutations in FGF17, IL17RD, DUSP6, SPRY4, and FLRT3 Are Identified in Individuals with Congenital Hypogonadotropic Hypogonadism.

Congenital hypogonadotropic hypogonadism (CHH) and its anosmia-associated form (Kallmann syndrome [KS]) are genetically heterogeneous. Among the >15 genes implicated in these conditions, mutations in FGF8 and FGFR1 account for ∼12% of cases; notably, KAL1 and HS6ST1 are also involved in FGFR1 signaling and can be mutated in CHH. We therefore hypothesized that mutations in genes encoding a broader range of modulators of the FGFR1 pathway might contribute to the genetics of CHH as causal or modifier mutations. Thus, we aimed to (1) investigate whether CHH individuals harbor mutations in members of the so-called "FGF8 synexpression" group and (2) validate the ability of a bioinformatics algorithm on the basis of protein-protein interactome data (interactome-based affiliation scoring [IBAS]) to identify high-quality candidate genes. On the basis of sequence homology, expression, and structural and functional data, seven genes were selected and sequenced in 386 unrelated CHH individuals and 155 controls. Except for FGF18 and SPRY2, all other genes were found to be mutated in CHH individuals: FGF17 (n = 3 individuals), IL17RD (n = 8), DUSP6 (n = 5), SPRY4 (n = 14), and FLRT3 (n = 3). Independently, IBAS predicted FGF17 and IL17RD as the two top candidates in the entire proteome on the basis of a statistical test of their protein-protein interaction patterns to proteins known to be altered in CHH. Most of the FGF17 and IL17RD mutations altered protein function in vitro. IL17RD mutations were found only in KS individuals and were strongly linked to hearing loss (6/8 individuals). Mutations in genes encoding components of the FGF pathway are associated with complex modes of CHH inheritance and act primarily as contributors to an oligogenic genetic architecture underlying CHH.

Association of the hemochromatosis gene with pazopanib-induced transaminase elevation in renal cell carcinoma.

BACKGROUND & AIMS: Pazopanib has demonstrated clinical benefit in patients with advanced renal cell carcinoma (RCC) and is generally well tolerated. However, transaminase elevations have commonly been observed. This 2-stage study sought to identify genetic determinants of alanine transaminase (ALT) elevations in pazopanib-treated white patients with RCC.¦METHODS: Data from two separate clinical studies were used to examine the association of genetic polymorphisms with maximum on-treatment ALT levels.¦RESULTS: Of 6852 polymorphisms in 282 candidate genes examined in an exploratory dataset of 115 patients, 92 polymorphisms in 40 genes were significantly associated with ALT elevation (p<0.01). Two markers (rs2858996 and rs707889) in the HFE gene, which are not yet known to be associated with hemochromatosis, showed evidence for replication. Because of multiple comparisons, there was a 12% likelihood the replication occurred by chance. These two markers demonstrated strong linkage disequilibrium (r(2)=0.99). In the combined dataset, median (25-75th percentile) maximum ALT values were 1.2 (0.7-1.9), 1.1 (0.8-2.5), and 5.4 (1.9-7.6)×ULN for rs2858996 GG (n=148), GT (n=82), and TT (n=1 2) genotypes, respectively. All 12 TT patients had a maximum ALT>ULN, and 8 (67%) had ALT≥3×ULN. The odds ratio (95% CI) for ALT≥3×ULN for TT genotype was 39.7 (2.2-703.7) compared with other genotypes. As a predictor of ALT≥3×ULN, the TT genotype had a negative predictive value of 0.83 and positive predictive value of 0.67. No TT patients developed liver failure.¦CONCLUSIONS: The rs2858996/rs707889 polymorphisms in the HFE gene may be associated with reversible ALT elevation in pazo-panib-treated patients with RCC.