Highly accurate RNA and DNA sequencing: Application to longstanding questions in aging and cancer

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

Multivariate QTL linkage analysis suggests a QTL for platelet count on chromosome 19q

Platelet count is a highly heritable trait with genetic factors responsible for around 80% of the phenotypic variance. We measured platelet count longitudinally in 327 monozygotic and 418 dizygotic twin pairs at 12, 14 and 16 years of age. We also performed a genome-wide linkage scan of these twins and their families in an attempt to localize QTLs that influenced variation in platelet concentrations. Suggestive linkage was observed on chromosome 19q13.13-19q13.31 at 12 (LOD=2.12, P=0.0009), 14 (LOD=2.23, P=0.0007) and 16 (LOD=1.01, P=0.016) years of age and multivariate analysis of counts at all three ages increased the LOD to 2.59 (P=0.0003). A possible candidate in this region is the gene for glycoprotein VI, a receptor involved in platelet aggregation. Smaller linkage peaks were also seen at 2p, 5p, 5q, 10p and 15q. There was little evidence for linkage to the chromosomal regions containing the genes for thrombopoietin (3q27) and the thrombopoietin receptor (1q34), suggesting that polymorphisms in these genes do not contribute substantially to variation in platelet count between healthy individuals.

Interactive effects of MC1R and OCA2 on melanoma risk phenotypes

The relationships between MC1R gene variants and red hair, skin reflectance, degree of freckling and nevus count were investigated in 2331 adolescent twins, their sibs and parents in 645 twin families. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. Of nine MC1R variant alleles assayed, four common alleles were strongly associated with red hair and fair skin (Asp84Glu, Arg151Cys, Arg160Trp and Asp294His), with a further three alleles having low penetrance (Val60Leu, Val92Met and Arg163Gln). These variants were separately combined for the purposes of this analysis and designated as strong 'R' (OR=63.3; 95% CI 31.9-139.6) and weak 'r ' (OR=5.1; 95% CI 2.5-11.3) red hair alleles. Red-haired individuals are predominantly seen in the R/R and R/r groups with 67.1 and 10.8%, respectively. To assess the interaction of the brown eye color gene OCA2 on the phenotypic effects of variant MC1R alleles we included eye color as a covariate, and also genotyped two OCA2 SNPs (Arg305Trp and Arg419Gln), which were confirmed as modifying eye color. MC1R genotype effects on constitutive skin color, freckling and mole count were modified by eye color, but not genotype for these two OCA2 SNPs. This is probably due to the association of these OCA2 SNPs with brown/green not blue eye color. Amongst individuals with a R/R genotype (but not R/r), those who also had brown eyes had a mole count twice that of those with blue eyes. This suggests that other OCA2 polymorphisms influence mole count and remain to be described.

Aprataxin, a novel protein that protects against genotoxic stress

Ataxia-oculomotor apraxia (AOA1) is a neurological disorder with symptoms that overlap those of ataxia-telangiectasia, a syndrome characterized by abnormal responses to double-strand DNA breaks and genome instability. The gene mutated in AOA1, APTX, is predicted to code for a protein called aprataxin that contains domains of homology with proteins involved in DNA damage signalling and repair. We demonstrate that aprataxin is a nuclear protein, present in both the nucleoplasm and the nucleolus. Mutations in the APTX gene destabilize the aprataxin protein, and fusion constructs of enhanced green fluorescent protein and aprataxin, representing deletions of putative functional domains, generate highly unstable products. Cells from AOA1 patients are characterized by enhanced sensitivity to agents that cause single-strand breaks in DNA but there is no evidence for a gross defect in single-strand break repair. Sensitivity to hydrogen peroxide and the resulting genome instability are corrected by transfection with full-length aprataxin cDNA. We also demonstrate that aprataxin interacts with the repair proteins XRCC1, PARP-1 and p53 and that it co-localizes with XRCC1 along charged particle tracks on chromatin. These results demonstrate that aprataxin influences the cellular response to genotoxic stress very likely by its capacity to interact with a number of proteins involved in DNA repair.

Sequence variation in the proximity of IDE may impact age at onset of both Parkinson disease and Alzheimer disease

We recently reported that a linkage disequilibrium (LD) block on chromosome 10q encompassing the gene encoding insulin-degrading enzyme (IDE) harbors sequence variants that associate with Alzheimer disease (AD). Evidence also indicated effects upon a number of quantitative indices of AD severity, including age-at-onset (AAO). Since linkage of this immediate region to AAO has been shown in both AD and Parkinson disease (PD), we have explored the possibility that polymorphism within this LD block might also influence PD. Utilizing single nucleotide polymorphisms that delineate common haplotypes from this region, we observed significant evidence of association with AAO in an Australian PD case-control sample. Analyses were complemented with AAO data from two independent Swedish AD case samples, for which previously reported findings were replicated. Results were consistent between AD and PD, suggesting the presence of equivalent detrimental and protective alleles. These data highlight a genomic region in the proximity of IDE that may contribute to AD and PD in a similar manner.

The estrogen receptor 1 G594A polymorphism is associated with migraine susceptibility in two independent case/control groups

Migraine is a painful and debilitating disorder with a significant genetic component. Steroid hormones, in particular estrogen, have long been considered to play a role in migraine, as variations in hormone levels are associated with migraine onset in many sufferers of the disorder. Steroid hormones mediate their activity via hormone receptors, which have a wide tissue distribution. Estrogen receptors have been localized to the brain in regions considered to be involved in migraine pathogenesis. Hence it is possible that genetic variation in the estrogen receptor gene may play a role in migraine susceptibility. This study thus examined the estrogen receptor 1 (ESRalpha) gene for a potential role in migraine pathogenesis and susceptibility. A population-based cohort of 224 migraine sufferers and 224 matched controls were genotyped for the G594A polymorphism located in exon 8 of the ESR1 gene. Statistical analysis indicated a significant difference between migraineurs and non-migraineurs in both the allele frequencies (P=0.003) and genotype distributions (P=0.008) in this sample. An independent follow-up study was then undertaken using this marker in an additional population-based cohort of 260 migraine sufferers and 260 matched controls. This resulted in a significant association between the two groups with regard to allele frequencies (P=8x10(-6)) and genotype distributions (P=4x10(-5)). Our findings support the hypothesis that genetic variation in hormone receptors, in particular the ESR1 gene, may play a role in migraine.

Normal HPRT coding region in a male with gout due to HPRT deficiency (Brief Communication)

A deficiency of the enzyme hypoxanthine-guanine phosphoribosyltransferase (HPRT; EC 2.4.2.8) is associated with a spectrum of disease that ranges from gouty arthritis (OMIM 300323) to the more severe Lesch-Nyhan syndrome (OMIM 300322). To date, all cases of HPRT deficiency have shown a mutation within the HPRT cDNA. In the present study of an individual with gout due to HPRT deficiency, we found a normal HPRT cDNA sequence. This is the first study to provide an example of HPRT deficiency which appears to be due to a defect in the regulation of the gene. © 2005 Elsevier Inc. All rights reserved.

Altered cell surface expression of human MC1R variant receptor alleles associated with red hair and skin cancer risk

The human melanocortin-1 receptor gene (MC1R) encodes a G-protein coupled receptor that is primarily expressed on melanocytes, where it plays a key role in pigmentation regulation. Variant alleles are associated with red hair colour and fair skin, known as the RHC phenotype, as well as skin cancer risk. The R151C, R160W and D294H alleles, designated 'R', are strongly associated with the RHC phenotype and have been proposed to result in loss of function receptors due to impaired G-protein coupling. We recently provided evidence that the R151C and R160W variants can efficiently couple to G-proteins in response to alpha-melanocyte stimulating hormone. The possibility that altered cellular localization of the R151C and R160W variant receptors could underlie their association with RHC was therefore considered. Using immunofluorescence and ligand binding studies, we found that melanocytic cells exogenously or endogenously expressing MC1R show strong surface localization of the wild-type and D294H alleles but markedly reduced cell surface expression of the R151C and R160W receptors. In additional exogenous expression studies, the R variant D84E and the rare I155T variant, also demonstrated a significant reduction in plasma membrane receptor numbers. The V60L, V92M and R163Q weakly associated RHC alleles, designated 'r', were expressed with normal or intermediate cell surface receptor levels. These results indicate that reduced receptor coupling activity may not be the only contributing factor to the genetic association between the MC1R variants and the RHC phenotype, with MC1R polymorphisms now linked to a change in receptor localization.

Small regulatory RNAs in mammals

Mammalian cells harbor numerous small non-protein-coding RNAs, including small nucleolar RNAs (snoRNAs), microRNAs (miRNAs), short interfering RNAs (siRNAs) and small double-stranded RNAs, which regulate gene expression at many levels including chromatin architecture, RNA editing, RNA stability, translation, and quite possibly transcription and splicing. These RNAs are processed by multistep pathways from the introns and exons of longer primary transcripts, including protein-coding transcripts. Most show distinctive temporal- and tissue-specific expression patterns in different tissues, including embryonal stem cells and the brain, and some are imprinted. Small RNAs control a wide range of developmental and physiological pathways in animals, including hematopoietic differentiation, adipocyte differentiation and insulin secretion in mammals, and have been shown to be perturbed in cancer and other diseases. The extent of transcription of non-coding sequences and the abundance of small RNAs suggests the existence of an extensive regulatory network on the basis of RNA signaling which may underpin the development and much of the phenotypic variation in mammals and other complex organisms and which may have different genetic signatures from sequences encoding proteins.

Zebrafish as a model for caveolin-associated muscle disease; caveolin-3 is required for myofibril organization and muscle cell patterning

Caveolae are an abundant feature of many animal cells. However, the exact function of caveolae remains unclear. We have used the zebrafish, Danio rerio, as a system to understand caveolae function focusing on the muscle-specific caveolar protein, caveolin-3 (Cav3). We have identified caveolin-1 (alpha and beta), caveolin-2 and Cav3 in the zebrafish. Zebrafish Cav3 has 72% identity to human CAV3, and the amino acids altered in human muscle diseases are conserved in the zebrafish protein. During embryonic development, cav3 expression is apparent by early segmentation stages in the first differentiating muscle precursors, the adaxial cells and slightly later in the notochord. cav3 expression appears in the somites during mid-segmentation stages and then later in the pectoral fins and facial muscles. Cav3 and caveolae are located along the entire sarcolemma of late stage embryonic muscle fibers, whereas beta-dystroglycan is restricted to the muscle fiber ends. Down-regulation of Cav3 expression causes gross muscle abnormalities and uncoordinated movement. Ultrastructural analysis of isolated muscle fibers reveals defects in myoblast fusion and disorganized myofibril and membrane systems. Expression of the zebrafish equivalent to a human muscular dystrophy mutant, CAV3P104L, causes severe disruption of muscle differentiation. In addition, knockdown of Cav3 resulted in a dramatic up-regulation of eng1a expression resulting in an increase in the number of muscle pioneer-like cells adjacent to the notochord. These studies provide new insights into the role of Cav3 in muscle development and demonstrate its requirement for correct intracellular organization and myoblast fusion.

Expression profiling of purified mouse gonadal somatic cells during the critical time window of sex determination reveals novel candidate genes for human sexual dysgenesis syndromes

Despite the identification of SRY as the testis-determining gene in mammals, the genetic interactions controlling the earliest steps of male sex determination remain poorly understood. In particular, the molecular lesions underlying a high proportion of human XY gonadal dysgenesis, XX maleness and XX true hermaphroditism remain undiscovered. A number of screens have identified candidate genes whose expression is modulated during testis or ovary differentiation in mice, but these screens have used whole gonads, consisting of multiple cell types, or stages of gonadal development well beyond the time of sex determination. We describe here a novel reporter mouse line that expresses enhanced green fluorescent protein under the control of an Sf1 promoter fragment, marking Sertoli and granulosa cell precursors during the critical period of sex determination. These cells were purified from gonads of male and female transgenic embryos at 10.5 dpc (shortly after Sry transcription is activated) and 11.5 dpc (when Sox9 transcription begins), and their transcriptomes analysed using Affymetrix genome arrays. We identified 266 genes, including Dhh, Fgf9 and Ptgds, that were upregulated and 50 genes that were downregulated in 11.5 dpc male somatic gonad cells only, and 242 genes, including Fst, that were upregulated in 11.5 dpc female somatic gonad cells only. The majority of these genes are novel genes that lack identifiable homology, and several human orthologues were found to map to chromosomal loci implicated in disorders of sexual development. These genes represent an important resource with which to piece together the earliest steps of sex determination and gonad development, and provide new candidates for mutation searching in human sexual dysgenesis syndromes.

How lifetimes shape epigenotype within and across generations

Despite our detailed characterization of the human genome at the level of the primary DNA sequence, we are still far from understanding the molecular events underlying phenotypic variation. Epigenetic modifications to the DNA sequence and associated chromatin are known to regulate gene expression and, as such, are a significant contributor to phenotype. Studies of inbred mice and monozygotic twins show that variation in the epigenotype can be seen even between genetically identical individuals and that this, in some cases at least, is associated with phenotypic differences. Moreover, recent evidence suggests that the epigenome can be influenced by the environment and these changes can last a lifetime. However, we also know that epigenetic states in real-time are in continual flux and, as a result, the epigenome exhibits instability both within and across generations. We still do not understand the rules governing the establishment and maintenance of the epigenotype at any particular locus. The underlying DNA sequence itself and the sequence at unlinked loci (modifier loci) are certainly involved. Recent support for the existence of transgenerational epigenetic inheritance in mammals suggests that the epigenetic state of the locus in the previous generation may also play a role. Over the next decade, many of these processes will be better understood, heralding a greater capacity for us to correlate measurable molecular marks with phenotype and providing the opportunity for improved diagnosis and presymptomatic healthcare.