Antigenic and genetic typing of Whataroa viruses in Australia

We recently characterized three novel alphaviruses isolated from mosquitoes captured in New South Wales, Australia. Initial cross-neutralization studies revealed antigenic similarity to the Sindbis virus (SINV)-like Whataroa virus (WHAV), heretofore found only in New Zealand. Nucleotide sequence analysis showed that the WHAV-Iike viruses shared >99% nucleotide sequence similarity with each other, and 96-97% similarity with prototype WHAV. Enzyme-linked immunosorbent assay reactions of a panel of monoclonal antibodies to SINV showed that the novel WHAV-Iike viruses displayed identical binding patterns and were antigenically distinct from all SINV isolates examined. Although these viruses displayed a similar binding pattern to prototype WHAV, three monoclonal antibodies discriminated them from the New Zealand virus. Our results suggest that these novel alphaviruses are antigenic variants of WHAV and represent the first reported isolations of this virus from outside New Zealand. The monoclonal antibodies used in this study will be useful for typing new SINV and SINV-like isolates.

Molecular characterization and cancer risk associated with BRCA1 and BRCA2 splice site variants identified in multiple-case breast cancer families

Genetic screening of women from multiple-case breast cancer families and other research-based endeavors have identified an extensive collection of germline variations of BRCA1 and BRCA2 that can be classified as deleterious and have clinical relevance. For some variants, such as those in the conserved intronic splice site regions which are highly likely to alter splicing, it is not possible to classify them based on the identified DNA sequence variation alone. We studied 11 multiple-case breast cancer families carrying seven distinct splice site region genetic alterations in BRCA1 or BRCA2 (BRCA1, c.IVS6-2delA, c.IVS9-2A>C, c.IVS4-1G>T, c.IVS20+1G>A and BRCA2, c.IVS17-1G>C, c.IVS20+1G>A, c.IVS7-1G>A) and applied SpliceSiteFinder to predict possible changes in efficiency of splice donor and acceptor sites, characterized the transcripts, and estimated the average age-specific cumulative risk (penetrance) using a modified segregation analysis. SpliceSiteFinder predicted and we identified transcipts that illustrated that all variants caused exon skipping, and all but two led to frameshifts. The risks of breast cancer to age 70 yrs, averaged over all variants, over BRCA1 variants alone, and over BRCA2 variants alone, were 73% (95% confidence interval 47-93), 64% (95%CI 28-96) and 79% (95%CI 48-98) respectively (all P

Investigation of dyslexia and SLI risk variants in reading- and language-impaired subjects

Dyslexia (or reading disability) and specific language impairment (or SLI) are common childhood disorders that show considerable co-morbidity and diagnostic overlaps and have been suggested to share some genetic aetiology. Recently, genetic risk variants have been identified for SLI and dyslexia enabling the direct evaluation of possible shared genetic influences between these disorders. In this study we investigate the role of variants in these genes (namely MRPL19/C20RF3, ROBO1, DCDC2, KIAA0319, DYX1C1, CNTNAP2, ATP2C2 and CMIP) in the aetiology of SLI and dyslexia. We perform case–control and quantitative association analyses using measures of oral and written language skills in samples of SLI and dyslexic families and cases. We replicate association between KIAA0319 and DCDC2 and dyslexia and provide evidence to support a role for KIAA0319 in oral language ability. In addition, we find association between reading-related measures and variants in CNTNAP2 and CMIP in the SLI families.

Circadian gene variants and susceptibility to type 2 diabetes:a pilot study

Disruption of endogenous circadian rhythms has been shown to increase the risk of developing type 2 diabetes, suggesting that circadian genes might play a role in determining disease susceptibility. We present the results of a pilot study investigating the association between type 2 diabetes and selected single nucleotide polymorphisms (SNPs) in/near nine circadian genes. The variants were chosen based on their previously reported association with prostate cancer, a disease that has been suggested to have a genetic link with type 2 diabetes through a number of shared inherited risk determinants.

Common variants of the TCF7L2 gene are associated with increased risk of type 2 diabetes mellitus in a UK-resident South Asian population

Background - Recent studies have implicated variants of the transcription factor 7-like 2 (TCF7L2) gene in genetic susceptibility to type 2 diabetes mellitus in several different populations. The aim of this study was to determine whether variants of this gene are also risk factors for type 2 diabetes development in a UK-resident South Asian cohort of Punjabi ancestry. Methods - We genotyped four single nucleotide polymorphisms (SNPs) of TCF7L2 (rs7901695, rs7903146, rs11196205 and rs12255372) in 831 subjects with diabetes and 437 control subjects. Results - The minor allele of each variant was significantly associated with type 2 diabetes; the greatest risk of developing the disease was conferred by rs7903146, with an allelic odds ratio (OR) of 1.31 (95% CI: 1.11 – 1.56, p = 1.96 × 10-3). For each variant, disease risk associated with homozygosity for the minor allele was greater than that for heterozygotes, with the exception of rs12255372. To determine the effect on the observed associations of including young control subjects in our data set, we reanalysed the data using subsets of the control group defined by different minimum age thresholds. Increasing the minimum age of our control subjects resulted in a corresponding increase in OR for all variants of the gene (p ≤ 1.04 × 10-7). Conclusion - Our results support recent findings that TCF7L2 is an important genetic risk factor for the development of type 2 diabetes in multiple ethnic groups.

Assessing Methods for Assigning SNPs to Genes in Gene-Based Tests of Association Using Common Variants

Gene-based tests of association are frequently applied to common SNPs (MAF>5%) as an alternative to single-marker tests. In this analysis we conduct a variety of simulation studies applied to five popular gene-based tests investigating general trends related to their performance in realistic situations. In particular, we focus on the impact of non-causal SNPs and a variety of LD structures on the behavior of these tests. Ultimately, we find that non-causal SNPs can significantly impact the power of all gene-based tests. On average, we find that the “noise” from 6–12 non-causal SNPs will cancel out the “signal” of one causal SNP across five popular gene-based tests. Furthermore, we find complex and differing behavior of the methods in the presence of LD within and between non-causal and causal SNPs. Ultimately, better approaches for a priori prioritization of potentially causal SNPs (e.g., predicting functionality of non-synonymous SNPs), application of these methods to sequenced or fully imputed datasets, and limited use of window-based methods for assigning inter-genic SNPs to genes will improve power. However, significant power loss from non-causal SNPs may remain unless alternative statistical approaches robust to the inclusion of non-causal SNPs are developed.

Genetic markers, birth characteristics, and childhood leukemia risk

The cause for childhood acute lymphoblastic leukemia (ALL) remains unknown, but male gender is a risk factor, and among ethnicities, Hispanics have the highest risk. In this dissertation, we explored correlations among genetic polymorphisms, birth characteristics, and the risk of childhood ALL in a multi-ethnic sample in 161 cases and 231 controls recruited contemporaneously (2007-2012) in Houston, TX. We first examined three lymphoma risk markers, since lymphoma and ALL both stem from lymphoid cells. Of these, rs2395185 showed a risk association in non-Hispanic White males (OR=2.8, P=0.02; P interaction=0.03 for gender), but not in Hispanics. We verified previously known risk associations to validate the case-control sample. Mutations of HFE (C282Y, H63D) were genotyped to test whether iron-regulatory gene (IRG) variants known to elevate iron levels increase childhood ALL risk. Being positive for either polymorphism yielded only a modestly elevated OR in males, which increased to 2.96 (P=0.01) in the presence of a particular transferrin receptor (TFRC) genotype for rs3817672 (Pinteraction=0.04). SNP rs3817672 itself showed an ethnicity-specific association (P interaction=0.02 for ethnicity). We then examined additional IRG SNPs (rs422982, rs855791, rs733655), which showed risk associations in males (ORs=1.52 to 2.60). A polygenic model based on the number of polymorphic alleles in five IRG SNPs revealed a linear increase in risk (OR=2.00 per incremental change; P=0.002). Having three or more alleles compared with none was associated with increased risk in males (OR=4.12; P=0.004). Significant risk associations with childhood ALL was found with birth length (OR=1.18 per inch, P=0.04), high birth weight (>4,000g) (OR=1.93, P=0.01), and with gestational age (OR=1.10 per week, P=0.04). We observed a negative correlation between HFE SNP rs9366637 and gestational age (P=0.005), again, stronger in males ( P=0.001) and interacting with TFRC (P interaction=0.05). Our results showed that (i) ALL risk markers do not show universal associations across ethnicities or between genders, (ii) IRG SNPs modify ALL risk presumably by their effects on iron levels, (iii) a negative correlation between an HFE SNP and gestational age exists, which implicates an iron-related mechanism. The results suggest that currently unregulated supplemental iron intake may have implications on childhood ALL development.

Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Functional Evaluation of Causal Mutations Identified in Human Genetic Studies

Human genetics has been experiencing a wave of genetic discoveries thanks to the development of several technologies, such as genome-wide association studies (GWAS), whole-exome sequencing, and whole genome sequencing. Despite the massive genetic discoveries of new variants associated with human diseases, several key challenges emerge following the genetic discovery. GWAS is known to be good at identifying the locus associated with the patient phenotype. However, the actually causal variants responsible for the phenotype are often elusive. Another challenge in human genetics is that even the causal mutations are already known, the underlying biological effect might remain largely ambiguous. Functional evaluation plays a key role to solve these key challenges in human genetics both to identify causal variants responsible for the phenotype, and to further develop the biological insights from the disease-causing mutations.

We adopted various methods to characterize the effects of variants identified in human genetic studies, including patient genetic and phenotypic data, RNA chemistry, molecular biology, virology, and multi-electrode array and primary neuronal culture systems. Chapter 1 is a broader introduction for the motivation and challenges for functional evaluation in human genetic studies, and the background of several genetics discoveries, such as hepatitis C treatment response, in which we performed functional characterization.

Chapter 2 focuses on the characterization of causal variants following the GWAS study for hepatitis C treatment response. We characterized a non-coding SNP (rs4803217) of IL28B (IFNL3) in high linkage disequilibrium (LD) with the discovery SNP identified in the GWAS. In this chapter, we used inter-disciplinary approaches to characterize rs4803217 on RNA structure, disease association, and protein translation.

Chapter 3 describes another avenue of functional characterization following GWAS focusing on the novel transcripts and proteins identified near the IL28B (IFNL3) locus. It has been recently speculated that this novel protein, which was named IFNL4, may affect the HCV treatment response and clearance. In this chapter, we used molecular biology, virology, and patient genetic and phenotypic data to further characterize and understand the biology of IFNL4. The efforts in chapter 2 and 3 provided new insights to the candidate causal variant(s) responsible for the GWAS for HCV treatment response, however, more evidence is still required to make claims for the exact causal roles of these variants for the GWAS association.

Chapter 4 aims to characterize a mutation already known to cause a disease (seizure) in a mouse model. We demonstrate the potential use of multi-electrode array (MEA) system for the functional characterization and drug testing on mutations found in neurological diseases, such as seizure. Functional characterization in neurological diseases is relatively challenging and available systematic tools are relatively limited. This chapter shows an exploratory research and example to establish a system for the broader use for functional characterization and translational opportunities for mutations found in neurological diseases.

Overall, this dissertation spans a range of challenges of functional evaluations in human genetics. It is expected that the functional characterization to understand human mutations will become more central in human genetics, because there are still many biological questions remaining to be answered after the explosion of human genetic discoveries. The recent advance in several technologies, including genome editing and pluripotent stem cells, is also expected to make new tools available for functional studies in human diseases.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.

Identification of potentially pathogenic variants in candidate breast and ovarian cancer susceptibility genes in BRCAX patients

Mutations within the BRCA1 and BRCA2 genes account for approximately 20% of hereditary breast cancers, with a further 10%–15% being attributable to rare mutations in moderate-risk genes and common variants in low-risk genes. The genes harbouring mutations in the remaining ∼65% of hereditary breast cancers are unknown. The identification of mutation carriers in hereditary breast and ovarian cancer (hboc) families is critical for determining who is most at risk of developing the disease and therefore who should be offered risk-reducing procedures or more intensive screening, or both.

Many of the high- and moderate-risk genes for hereditary breast cancers encode proteins that work in concert to maintain genomic stability and in dna damage signalling and repair. A novel BRCA1 protein complex identified within the research group whose target genes are involved in dna repair provided novel candidates for hboc susceptibility genes. These 12 candidate genes were sequenced in a cohort of 675 affected individuals from the Kathleen Cunningham Foundation Consortium for Research into Familial Breast Cancer (kConFab) with hereditary breast or ovarian cancer, but with no mutations in known susceptibility genes (BRCAx patients). This analysis identified 20 individuals (each from a different BRCAx family) with different potentially pathogenic variants across 6 of the candidate hboc susceptibility genes. The family members of each BRCAx index case were tested for the presence of the specific mutation identified in the proband to examine segregation with disease. To further expand on the potential role of the novel candidate hboc susceptibility genes identified in this study, the genetic variation of a second cohort of 520 Northern Irish BRCAx patients is being characterized using a 61-gene panel.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10(-5), Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10(-06)/Pfemales: 3.45 × 10(-07)/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease.

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimer's disease (LOAD)1, 2. These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

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

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