USING THE R PACKAGE crlmm FOR GENOTYPING AND COPY NUMBER ESTIMATION

Genotyping platforms such as Affymetrix can be used to assess genotype-phenotype as well as copy number-phenotype associations at millions of markers. While genotyping algorithms are largely concordant when assessed on HapMap samples, tools to assess copy number changes are more variable and often discordant. One explanation for the discordance is that copy number estimates are susceptible to systematic differences between groups of samples that were processed at different times or by different labs. Analysis algorithms that do not adjust for batch effects are prone to spurious measures of association. The R package crlmm implements a multilevel model that adjusts for batch effects and provides allele-specific estimates of copy number. This paper illustrates a workflow for the estimation of allele-specific copy number, develops markerand study-level summaries of batch effects, and demonstrates how the marker-level estimates can be integrated with complimentary Bioconductor software for inferring regions of copy number gain or loss. All analyses are performed in the statistical environment R. A compendium for reproducing the analysis is available from the author’s website (http://www.biostat.jhsph.edu/~rscharpf/crlmmCompendium/index.html).

Large deletions in the CFTR gene: clinics and genetics in Swiss patients with CF

Cystic fibrosis (CF) is the most common life-shortening autosomal recessive disorder in Caucasians, and is associated with at least one mutation on each CF transmembrane conductance regulator (CFTR) allele. Some patients, however, with only one identifiable point mutation carry on the other allele, a large deletion that is not detected by conventional screening methods. The overall frequency of large deletions in patients with CF is estimated to be 1-3%. Using the CFTR Multiplex Ligation dependent Probe Amplification Kit (MRC-Holland, Amsterdam, Netherlands) that allows the exact detection of copy numbers from all 27 exons in the CFTR gene, we screened 50 patients with only one identified mutation for large deletions in the CFTR gene. Each detected deletion was confirmed using our real-time polymerase chain reaction (PCR) assay and deletion-specific PCR reactions using junction fragment primers. We detected large deletions in eight patients (16%). These eight CF alleles belong to four different deletion types (CFTRindel2, CFTRdele14b-17b, CFTRdele17a-17b and CFTRdele 2-9) whereof the last is novel. Comparing detailed clinical data of all these patients with CF and the molecular genetic findings, we were able to elaborate criteria for deletion screenings and possible genotype-phenotype associations. In conclusion, we agree with other authors that deletion screenings should be implemented in routine genetic diagnostics of CF.

Detection of cryptic pathogenic copy number variations and constitutional loss of heterozygosity using high resolution SNP microarray analysis in 117 patients referred for cytogenetic analysis and impact on clinical practice

BACKGROUND: Microarray genome analysis is realising its promise for improving detection of genetic abnormalities in individuals with mental retardation and congenital abnormality. Copy number variations (CNVs) are now readily detectable using a variety of platforms and a major challenge is the distinction of pathogenic from ubiquitous, benign polymorphic CNVs. The aim of this study was to investigate replacement of time consuming, locus specific testing for specific microdeletion and microduplication syndromes with microarray analysis, which theoretically should detect all known syndromes with CNV aetiologies as well as new ones. METHODS: Genome wide copy number analysis was performed on 117 patients using Affymetrix 250K microarrays. RESULTS: 434 CNVs (195 losses and 239 gains) were found, including 18 pathogenic CNVs and 9 identified as "potentially pathogenic". Almost all pathogenic CNVs were larger than 500 kb, significantly larger than the median size of all CNVs detected. Segmental regions of loss of heterozygosity larger than 5 Mb were found in 5 patients. CONCLUSIONS: Genome microarray analysis has improved diagnostic success in this group of patients. Several examples of recently discovered "new syndromes" were found suggesting they are more common than previously suspected and collectively are likely to be a major cause of mental retardation. The findings have several implications for clinical practice. The study revealed the potential to make genetic diagnoses that were not evident in the clinical presentation, with implications for pretest counselling and the consent process. The importance of contributing novel CNVs to high quality databases for genotype-phenotype analysis and review of guidelines for selection of individuals for microarray analysis is emphasised.

Cardiac channel molecular autopsy: insights from 173 consecutive cases of autopsy-negative sudden unexplained death referred for postmortem genetic testing

OBJECTIVE To perform long QT syndrome and catecholaminergic polymorphic ventricular tachycardia cardiac channel postmortem genetic testing (molecular autopsy) for a large cohort of cases of autopsy-negative sudden unexplained death (SUD). METHODS From September 1, 1998, through October 31, 2010, 173 cases of SUD (106 males; mean ± SD age, 18.4 ± 12.9 years; age range, 1-69 years; 89% white) were referred by medical examiners or coroners for a cardiac channel molecular autopsy. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing, a comprehensive mutational analysis of the long QT syndrome susceptibility genes (KCNQ1, KCNH2, SCN5A, KCNE1, and KCNE2) and a targeted analysis of the catecholaminergic polymorphic ventricular tachycardia type 1-associated gene (RYR2) were conducted. RESULTS Overall, 45 putative pathogenic mutations absent in 400 to 700 controls were identified in 45 autopsy-negative SUD cases (26.0%). Females had a higher yield (26/67 [38.8%]) than males (19/106 [17.9%]; P<.005). Among SUD cases with exercise-induced death, the yield trended higher among the 1- to 10-year-olds (8/12 [66.7%]) compared with the 11- to 20-year-olds (4/27 [14.8%]; P=.002). In contrast, for those who died during a period of sleep, the 11- to 20-year-olds had a higher yield (9/25 [36.0%]) than the 1- to 10-year-olds (1/24 [4.2%]; P=.01). CONCLUSION Cardiac channel molecular autopsy should be considered in the evaluation of autopsy-negative SUD. Several interesting genotype-phenotype observations may provide insight into the expected yields of postmortem genetic testing for SUD and assist in selecting cases with the greatest potential for mutation discovery and directing genetic testing efforts.

Clinical and genetic characteristics of long QT syndrome

Long QT syndrome (LQTS) is an arrhythmogenic ion channel disorder characterized by severely abnormal ventricular repolarization, which results in prolongation of the electrocardiographic QT interval. The condition is associated with sudden cardiac death due to malignant ventricular arrhythmias similar in form to the hallmark torsade de pointes. Eleven years after the identification of the principle cardiac channels involved in the condition, hundreds of mutations in, to date, 10 genes have been associated with the syndrome. Genetic investigations carried out up until the present have shown that, although the severe form of the disease is sporadic, there are a number of common polymorphisms in genes associated with the condition that may confer susceptibility to the development of torsade de pointes in some individuals, particularly when specific drugs are being administered. Moreover, some polymorphisms have been shown to have regulatory properties that either enhance or counteract a particular mutation's impact. Understanding of the molecular processes underlying the syndrome has enabled treatment to be optimized and has led to better survival among sufferers, thereby demonstrating a key correspondence between genotype, phenotype and therapy. Despite these developments, a quarter of patients do not have mutations in the genes identified to date. Consequently, LQTS continues to be an area of active research. This article contains a summary of the main clinical and genetic developments concerning the syndrome that have taken place during the last decade.

New perspectives in long QT syndrome

Long QT Syndrome (LQTS) is a cardiac channelopathy characterized by prolonged ventricular repolarization and increased risk to sudden death secondary to ventricular dysrrhythmias. Was the first cardiac channelopathy described and is probably the best understood. After a decade of the sentinel identification of ion channel mutation in LQTS, genotype-phenotype correlations have been developed along with important improvement in risk stratification and genetic guided-treatment. Genetic screening has shown that LQTS is more frequent than expected and interestingly, ethnic specific polymorphism conferring increased susceptibility to drug induced QT prolongation and torsades de pointes have been identified. A better understanding of ventricular arrhythmias as an adverse effect of ion channel binding drugs, allow the development of more safety formulas and better control of this public health problem. Progress in understanding the molecular basis of LQTS has been remarkable; eight different genes have been identified, however still 25% of patients remain genotype-negative. This article is an overview of the main LQTS knowledge developed during the last years.

Transforming growth factor-β and inflammation in vascular (type IV) Ehlers-Danlos syndrome.

BACKGROUND Vascular Ehlers-Danlos syndrome (VEDS) causes reduced life expectancy because of arterial dissections/rupture and hollow organ rupture. Although the causative gene, COL3A1, was identified >20 years ago, there has been limited progress in understanding the disease mechanisms or identifying treatments. METHODS AND RESULTS We studied inflammatory and transforming growth factor-β (TGF-β) signaling biomarkers in plasma and from dermal fibroblasts from patients with VEDS. Analyses were done in terms of clinical disease severity, genotype-phenotype correlations, and body composition and fat deposition alterations. VEDS subjects had increased circulating TGF-β1, TGF-β2, monocyte chemotactic protein-1, C-reactive protein, intercellular adhesion molecule-1, vascular cell adhesion molecule-1, and leptin and decreased interleukin-8 versus controls. VEDS dermal fibroblasts secreted more TGF-β2, whereas downstream canonical/noncanonical TGF-β signaling was not different. Patients with COL3A1 exon skipping mutations had higher plasma intercellular adhesion molecule-1 and vascular cell adhesion molecule-1, and VEDS probands had abnormally high plasma C-reactive protein versus affected patients identified through family members before any disease manifestations. Patients with VEDS had higher mean platelet volumes, suggesting increased platelet turnover because of ongoing vascular damage, as well as increased regional truncal adiposity. CONCLUSIONS These findings suggest that VEDS is a systemic disease with a major inflammatory component. C-reactive protein is linked to disease state and may be a disease activity marker. No changes in downstream TGF-β signaling and increased platelet turnover suggest that chronic vascular damage may partially explain increased plasma TGF-β1. Finally, we found a novel role for dysregulated TGF-β2, as well as adipocyte dysfunction, as demonstrated through reduced interleukin-8 and elevated leptin in VEDS.

The Vacuolar H+-ATPase B1 Subunit Polymorphism p.E161K Associates with Impaired Urinary Acidification in Recurrent Stone Formers

Mutations in the vacuolar–type H+-ATPase B1 subunit gene ATP6V1B1 cause autosomal–recessive distal renal tubular acidosis (dRTA). We previously identified a single-nucleotide polymorphism (SNP) in the human B1 subunit (c.481G.A; p.E161K) that causes greatly diminished pump function in vitro. To investigate the effect of this SNP on urinary acidification, we conducted a genotype-phenotype analysis of recurrent stone formers in theDallas and Bern kidney stone registries. Of 555 patients examined, 32 (5.8%) were heterozygous for the p.E161K SNP, and the remaining 523 (94.2%) carried two wild–type alleles. After adjustment for sex, age, body mass index, and dietary acid and alkali intake, p.E161K SNP carriers had a nonsignificant tendency to higher urinary pH on a random diet (6.31 versus 6.09; P=0.09). Under an instructed low–Ca and low–Na diet, urinary pH was higher in p.E161K SNP carriers (6.56 versus 6.01; P,0.01). Kidney stones of p.E161K carriers were more likely to contain calcium phosphate than stones of wild-type patients. In acute NH4Cl loading, p.E161K carriers displayed a higher trough urinary pH (5.34 versus 4.89; P=0.01) than wild-type patients. Overall, 14.6% of wild-type patients and 52.4% of p.E161K carriers were unable to acidify their urine below pH 5.3 and thus, can be considered to have incomplete dRTA. In summary, our data indicate that recurrent stone formers with the vacuolar H+-ATPase B1 subunit p.E161K SNP exhibit a urinary acidification deficit with an increased prevalence of calcium phosphate– containing kidney stones. The burden of E161K heterozygosity may be a forme fruste of dRTA.

Cellular hyper-excitability caused by mutations that alter the activation process of voltage-gated sodium channels.

Voltage-gated sodium channels (Nav) are widely expressed as macro-molecular complexes in both excitable and non-excitable tissues. In excitable tissues, the upstroke of the action potential is the result of the passage of a large and rapid influx of sodium ions through these channels. NaV dysfunction has been associated with an increasingly wide range of neurological, muscular and cardiac disorders. The purpose of this review is to summarize the recently identified sodium channel mutations that are linked to hyper-excitability phenotypes and associated with the alteration of the activation process of voltage gated sodium channels. Indeed, several clinical manifestations that demonstrate an alteration of tissue excitability were recently shown to be strongly associated with the presence of mutations that affect the activation process of the Nav. These emerging genotype-phenotype correlations have expanded the clinical spectrum of sodium channelopathies to include disorders which feature a hyper-excitability phenotype that may or may not be associated with a cardiomyopathy. The p.I141V mutation in SCN4A and SCN5A, as well as its homologous p.I136V mutation in SCN9A, are interesting examples of mutations that have been linked to inherited hyperexcitability myotonia, exercise-induced polymorphic ventricular arrhythmias and erythromelalgia, respectively. Regardless of which sodium channel isoform is investigated, the substitution of the isoleucine to valine in the locus 141 induces similar modifications in the biophysical properties of the Nav by shifting the voltage-dependence of steady state activation toward more negative potentials.

Identification of genetic variants that influence atherosclerosis in young persons

Atherosclerosis is widely accepted as a complex genetic phenotype and is the usual cause of cardiovascular disease, the world’s leading killer. Genetic factors have been proven to be important risk contributors for atherosclerosis and much work has been done to identify promising candidates that might play a role in the development of atherosclerosis. It is well known that many independent replications are needed to unequivocally establish a valid genotype-phenotype association across different populations before the findings are extended to clinical settings and to the expensive follow-up studies designed to identify causal genetic variants. Aiming to replicate the association with atherosclerosis in the Pathobiological Determinants of Atherosclerosis in Youth (PDAY) study, we assessed the relationship of 32 atherosclerosis candidate SNPs to atherosclerosis in the PDAY cohort, consisting of AA and EA young people aged 15-34 years who died of non-medical causes. Two association studies, a whole sample study and a 1:1 matched case control study were performed by use of multiple linear regression and logistic regression analyses, respectively. For the whole sample association study, 32 SNPs among 2,650 individuals (1,369 AA and 1,281 EA) were tested for the association with six early atherosclerosis phenotypes: abdominal aorta fatty streaks, abdominal aorta raised lesions, right coronary artery fatty streaks, right coronary artery raised lesions, thoracic aorta fatty streaks, and thoracic aorta raised lesions. For the matched case-control association study, 337 case-control paired samples were included; cases were chosen with the highest total raised lesion scores from the studied population, while controls were randomly selected from individuals that had no raised lesions and matched to cases by age, gender and race. Sixteen SNPs in 13 genes were found to be significantly associated with atherosclerosis in at least one of the PDAY association studies. Among these 16 findings: eight SNPs (rs9579646, rs6053733, rs3849150, rs10499903, rs2148079, rs5073691, rs10116277, and rs17228212) successfully replicated previous results, six SNPs (rs17222814, rs10811661, rs7028570, rs7291467, rs16996148 and rs10401969) were reported as new findings exclusive to our study, the last two of the 16 SNPs, rs501120 and rs6922269, showed either intriguing or conflicting result. SNP rs17222814 in ALOX5AP and SNP rs3849150 in LRRC18 were consistently associated with atherosclerosis in both prior and the two PDAY association studies. SNP rs3849150 was also identified to be highly correlated with a non-synonymous coding SNP, rs17772611, which may damage the protein (polyphen score = 0.996), suggesting that SNP rs17772611 may be the causal functional variant.^ In conclusion, our study added more support for the association of these candidate genes with atherosclerosis. SNPs rs3849150 and rs17772611 of LRRC18, as well as SNP rs17222814 of ALOX5AP, were the most significant findings from our study, and may be ranked among the best for further study.^

Ferroportin Disease (Haemochromatosis-type IV): a case report

Introdution: Haemochromatosis-type IV, the ferroportin disease, is characterized by an autosomal-dominant transmission and early iron accumulation in macrophages. It is caused by mutations in the transmembrane iron exporter protein ferroportin1 (SLC40A1 gene). In form A (classic), ferroportin loss of function mutants are unable to export iron from cells leading to cellular iron accumulation with decreased availability of iron for serum transferrin (TS). We present a Portuguese rare clinical case of HH-IV. Materials and Methods: A 41-year-old woman with hyperferritinemia and normal TS. Causes of hyperferritinemia (inflammation, chronic alcohol consumption, metabolic syndrome, cell necrosis, non-alcoholic fatty liver disease and aceruloplasminemia) were assessed. Liver iron, evaluated by magnetic resonance imaging (MRI) was carried out. Screening for mutation in HFE and SCL40A1 genes were performed by Sanger sequencing. Baseline: Ferritin:708ng/ml; TS: 27%; MRI:85µmol/g; Hb:13,6g/dl. Therapy: weekly 450ml Therapeutic Phlebotomies (TP) until ferritin≤50ng/ml. Results: Hyperferritinemia comorbidities and common genetic mutations for haemochromatosis were negative. However, sequencing of the patient SLC40A1 gene has revealed the presence in heterozygosity of the variant c.238G>A; p.Gly80Ser. Due to low tolerance to TP, we adopted smaller phlebotomies every three weeks. Conclusion: This patient has a rare autosomal-dominant Ferroportin disease due to a mutated ferroportin which is predicted to be defective in iron cellular export. In agreement, she presents hyperferritinemia, with normal TS and liver iron overload. The genotype/phenotype association allowed to diagnosis this rare FD case. Although a mild form A, we decided to start TP. Her father also has been treated for iron overload.

Certain observations concerning the effects of epistasis on complex traits and the evolution of genomes.

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

High throughput determination of sequence-function relationships in protein and RNA

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

CYP2D6: a global analysis of phenotypic and genotypic variation in search of radical cure of Plasmodium vivax malaria

Thesis (Master's)--University of Washington, 2016-06

Methodological Considerations in the Study of Genetic Discrimination

The potential significance and dimensions of genetic discrimination have been described extensively in published literature, but epidemiological and verified case data are limited. Obtaining unbiased data from individuals about discrimination which has been based on erroneous or unjustifiable assumptions about their genetic predispositions poses unique challenges. Through review and discussion of research literature, we identify methodological considerations for collecting valid epidemiological data on genetic discrimination from individuals in the community; in particular, we consider issues which relate to sampling, selection and response. We identify issues to promote sound study design, with particular attention to verification of genetic discrimination, and highlight the importance of clinical and genetic knowledge of complex genotype-phenotype relationships.