Chan, Ho Ching. / Thesis M.Phil. Chinese University of Hong Kong 2016. / Includes bibliographical references (leaves 159-166). / Abstracts also in Chinese. / Title from PDF title page (viewed on 15, September, 2016).

Stochastic fluctuations in mRNA and protein copy number between cells are inevitable during the process gene expression, even when cells carry identical chromosomes. Such fluctuations are able to impact the phenotypic fate of the cell, and are known to have greater impact when the copy number of the molecule involved is low. Additionally, up to 50% of proteins in Escherichia coli are present in the cell at a level of 10 molecules per cell or fewer (Taniguchi et al. 2010). As such, quantification of low copy number gene expression products and their distribution in cellular populations is key in understanding the process of gene expression. Currently, there are few techniques that allow investigation with the single cell and single molecule resolution required to study low copy number gene expression products. This work presents a novel method for protein quantification at the single molecule level, Quantitative HaloTag-TMR labelling, and uses the technique to quantify the absolute numbers of the low copy number RecB, RecC and RecD subunits of the bacterial DNA repair enzyme RecBCD, finding each subunit is present at between two and eight molecules per cell with mean numbers per cell of 4.9, 4.7 and 4.5 respectively. Additionally single molecule mRNA FISH was used to quantify the mRNA levels of recB and recD within cells, with means of 0.21 and 0.31 mRNA per cell being observed respectively. Finally this work presents a new method for use detecting both mRNA and protein simultaneously in individual cells by combining the HaloTag and FISH protocols to give HaloFISH. This work introduces two novel techniques that allow for single cell examination of gene expression, and investigates RecBCD expression at the single molecule level.

Clefts of the lip and/or palate (CL/P) are one of the most common birth defects in the world occurring about every 1 in 700 live births. Individuals with non-syndromic clefting (NSCL/P) account for about 70% of all cleft cases and exhibit a cleft only whereas syndromic occurrences (SCL/P) include additional cognitive or structural abnormalities. Linkage, genome-wide association, candidate gene, animal model, sequencing and copy number variant (CNV) analyses have been used to study CL/P and have established that it is a heterogeneous, complex disorder. However, the impact of identified sequence variants on protein structure and the contribution of structural genetic variation to CL/P remains poorly understood. In our first analysis we reassessed the phenotype of a 30-year-old individual of SCL/P and noticed phenotypic overlap with Hartsfield syndrome, a rare syndrome resulting from sequence variants in Fibroblast growth factor 1 (FGFR1). We sequenced the coding region of FGFR1 and identified a novel, de novo variant. Due to the fact sequence variants in FGFR1 contribute to multiple syndromes encompassing a wide phenotypic spectrum, we performed an extensive literature search to record every published sequence variant of FGFR1 and mapped it to the protein structure by disease and phenotype. Although no statistically significant protein domain-phenotype correlations were identified, many regions neared significance. This work stresses the need for systematic, comprehensive phenotyping of patients and provides a method for assessing the impact of the location of sequence variants within the 3D structure of the protein. Although rare and common CNVs have been identified in individuals with CL/P, prior to our work no large-scale studies of rare CNVs for the identification of novel clefting genes had been performed. For our second set of analyses, we conducted two such studies, first focusing on a smaller cohort of 140 individuals with NSCL/P from the Philippines to establish our informatic and functional validation pipeline. We used whole-genome tiling arrays to assess rare deletions overlapping genes not previously implicated in clefting, and identified one deletion overlapping Isthmin1 (ISM1) and a deletion just 3’ of the gene in a second affected individual. Functional validation of Ism1 in Xenopus laevis showed strong expression in structures necessary for craniofacial development, and morpholino and CRISPR/Cas9 knockdown of Ism1 resulted in a median cleft lip in some embryos, establishing ISM1 as a novel craniofacial patterning gene. We then expanded our study and assessed genomic CNVs in 1021 individuals with NSCL/P and 81 individuals with SCL/P, finding no differences in CNV number, load or burden between these groups. We also identified 8 putative clefting genes overlapped by deletions in two or more individuals but at a rare (< 1% frequency) in the cohort. Functional validation of these genes using CRISPR/Cas9 in zebrafish and Xenopus tropicalis is currently underway. This work has identified a novel sequence variant leading to the diagnosis of Hartsfield syndrome in an individual with SCL/P, developed an innovative method for assessing the impact of sequence variation on protein structure, improved our understanding of the contribution of CNVs to SCL/P and NSCL/P and identified several putative novel clefting loci which may help explain a portion of the missing heritability of CL/P.

In recent years, microarray technology has revealed the widespread presence of submicroscopic deletions and duplications throughout the human genome termed copy number variants (CNVs). CNVs have a profound effect on gene expression and are an important source of normal genetic variation. In addition, a small proportion of CNVs contribute to genetically simple and complex disease. This thesis focuses on the identification of pathogenic CNVs contributing to the etiology of diseases with "missing heritability" using a well-planned study design individually tailored to each disease cohort to optimize CNV detection and interpretation. We performed a genome-wide analysis for CNVs in five disease cohorts with genetic etiology: autism, age-related macular degeneration (AMD), glaucoma, clubfoot, and Bardet-Biedl syndrome (BBS). Our results indicate that CNVs likely account for a proportion of cases for each disease cohort reported in this thesis. Approximately 20% of our cohort of individuals with autism from trio pedigrees harbors a CNV known to confer risk to develop autism and we identified other novel and rare variants that may play a role in autism pathogenesis. We also characterized a duplication of 2p25.3 identified in two male half-siblings with autism and determined that their mother was somatic mosaic for the duplication. Our work provides evidence that this novel CNV disrupting the genes PXDN and MYT1L are the autism-causing mutation in this pedigree. A comparative cases experimental design was used in the study of AMD and glaucoma. While no common "risk CNVs" were identified for either eye disorder, we did identify several rare overlapping CNVs disrupting genes known to play a role in the eye that may confer risk to disease in a small proportion of individuals. In a fourth genetically complex disease, clubfoot, we identified a duplication of 17q23.2 disrupting the genes TBX4, NACA2, and BRIP1 that segregates with the autosomal dominant clubfoot phenotype in a large pedigree with 16 affected individuals. In addition, the duplication is within the linkage interval identified for this family. We also applied microarray technology to analyze the genomes of individuals with BBS, an autosomal recessive disorder, for the presence of CNVs in known BBS genes as well as CNVs that elucidate novel candidate genes for BBS. From 34 BBS patients with an unidentified mutation, we observed one CNV, a heterozygous deletion of BBS10, unmasking a BBS10 frameshift mutation. A promising BBS candidate gene also emerged from our studies, implicated by an intragenic deletion of the gene MARK3 predicted to result in a frameshift and premature truncation of the protein. Functional studies utilizing antisense morpholino gene knockdown in the zebrafish provide additional evidence that MARK3 is a BBS gene as knockdown of zebrafish mark3 results in a Kupffer's Vesicle defect and a melanosome transport delay, two cardinal BBS phenotypes in the zebrafish. In addition to identifying CNVs involved in disease, the work outlined in this thesis provides valuable insight into the study design and interpretation of a genome-wide analysis of CNV. This includes the appropriate use of controls and publicly available control databases, methods for enriching for CNVs in a patient cohort to maximize efficiency and discovery, and the importance of analyzing all patient cohorts with heritable disease for the presence of CNVs disrupting known disease genes and CNVs that implicate novel genetic candidates. As the reliability and resolution of CNV detection continue to improve, allowing detection of > 1,000 CNVs in each individual genome, it becomes more important than ever to have a well-defined study design for both the detection and interpretation of CNVs.

The 16p11.2 copy-number variant (CNV) represents a well-characterized, high-risk factor for autism spectrum disorder that additionally predisposes deletion carriers (16pdel) to increased head circumference, known as macrocephaly. The 16p11.2 CNV consists of 29 known genes, many of which are associated with neurobiological processes relevant for macrocephaly such as cell proliferation and apoptosis, differentiation and cell growth. Our lab’s previous work has demonstrated that induced pluripotent stem cell (iPSC)-derived neurons from 16pdel carriers show altered cellular morphology related to growth, which include increased soma size, total dendritic length and dendritic complexity. However, specific CNV genes responsible for these phenotypes have not been established. Here, we investigate the relationship between three 16p11.2 genes and the observed cellular phenotypes. We differentiated neurons from control iPSC-derived neural progenitor cells (NPCs) and used short hairpin RNA (shRNA) to reduce the expression of these CNV genes: KCTD13, MAPK3 and C16ORF53. We then assessed neuronal morphology by evaluating soma size, total dendritic length and dendritic complexity. We demonstrate that knocking down KCTD13 and C16ORF53 increases soma size and total dendrite length, respectively, similar to that observed in 16pdel iPSC-derived neurons. For this reason, we speculate that these genes may have a role in cell growth and might underlie macrocephaly. Thus, our study investigates genes in the 16p11.2 CNV that contribute to neuronal morphology, which may have a role in influencing brain size.

Glaucoma is the second leading cause of irreversible blindness in the United States and is the leading cause of blindness in African Americans. Cupping or excavation of the optic nerve, which sends the visual signal from the photoreceptors in the eye to the brain, is a chief feature of glaucoma. A similar excavated appearance of the optic nerve is also the primary clinical sign of other congenital malformations of the eye including optic nerve head coloboma, optic pit, and morning glory disc anomaly collectively termed cavitary optic disc anomaly (CODA). Clinical similarities between CODA and glaucoma have suggested that these conditions may have overlapping pathophysiology. Although risk factors are known, such as the elevated intraocular pressure (IOP) observed in some glaucoma subjects, the biological pathways and molecular events that lead to excavation of the optic disc in glaucoma and in CODA are incompletely understood, which has hindered efforts to improve diagnosis and treatment of these diseases. Consequently, there is a critical need to clarify the biological mechanisms that lead to excavation of the optic nerve, which will lead to improvements in our understanding of these important disease processes. Because of their similar clinical phenotypes and the limited therapy geared at lowering IOP in glaucoma patients, our central hypothesis is that genes involved in Mendelian forms of CODA would also be involved in a subset of glaucoma cases and may provide insight into glaucomatous optic neuropathy. The purpose of my research project has been to identify and functionally characterize the gene that causes congenital autosomal dominant CODA in a multiplex family with 17 affected members. The gene that causes CODA was previously mapped to chromosome 12q14 and following screening of candidate genes within the region that did not yield any plausible coding sequence mutations, a triplication of a 6KB segment of DNA upstream of the matrix metalloproteinase 19 (MMP19) gene was subsequently identified using comparative genomic hybridization arrays and qPCR. This copy number variation (CNV) was present in all affected family members but absent in unaffected family members, a panel of 78 normal control subjects, and the Database of Genomic Variants. In a case-control study of singleton CODA subjects, CNVs were also detected; we detected the same 6KB triplication in 1 of 24 subjects screened. This subject was part of another 3-generation autosomal dominant CODA pedigree where affected members each have the same CNV identified in the larger CODA pedigree. A separate case-control study with 172 glaucoma cases (primary open angle glaucoma = 84, normal tension glaucoma = 88) was evaluated for MMP19 CNVs, however none were detected. Although our cohort of CODA patients is small limiting our ability to accurately determine the proportion of CODA caused by MMP19 mutations, our data indicates that the MMP19 CNV is not an isolated case and additional CODA subjects may have MMP19 defects. Because of the location of the CNV, we evaluated its effect on downstream gene expression with luciferase reporter gene assays. These assays revealed that the 6KB sequence spanned by the CNV in CODA subjects functioned as a transcriptional enhancer; in particular, a 773bp segment had a strong positive influence (8-fold higher) on downstream gene expression. MMP19, a largely understudied gene, was further characterized by expression studies in the optic nerve and retina. Using frozen sections from normal donor eyes, we demonstrated that MMP19 is predominantly localized to the optic nerve head in the lamina cribrosa region with moderate labeling in the postlaminar region, and weak labeling in the prelaminar region and retina. We also evaluated MMP19 expression in relation to the cell types that populate the optic nerve such as astrocytes and retinal ganglion cells. The pattern of expression is consistent with MMP19 being a secreted protein accumulating in the extracellular spaces and basement membranes of the optic nerve. Our studies have identified the first gene associated with CODA and future research is focused on recapitulating CODA phenotypes in animal models and assessing the mechanism of MMP19 involvement during development.

Malignant glioma have poor prognosis resulting mainly from high level of cell proliferation and invasion and resistance to conventional therapy. Identification of novel targets that are critical elements in gliomagenesis may help improve therapeutic outcome. Using genome-wide explorations of a comprehensive glioma specimen population, we identified whole gain of chromosome 19 as one of the major chromosomal aberrations in high grade glioma that correlates to patients' outcomes. Our analysis revealed for the first time NOTCH3 as one of the most significant gene amplifications mapped to chromosome 19. This amplification is associated with worse outcome compared to tumors with non-amplified locus. NOTCH signaling pathway is essential for cell proliferation, stem cell maintenance and differentiation and its deregulation has been reported in several human cancers. NOTCHs are key positive regulators of cell-cell interactions, angiogenesis, cell adhesion and stem cell niche development which have been shown to play critical roles in gliomagenesis and glioma drug resistance. Our objective is to determine NOTCH3 molecular roles in glioma pathogenesis and aggressiveness. Here we show for the first time that NOTCH3 plays a role in glioma cell proliferation, cell migration, invasion and apoptosis. We also found a NOTCH3 glioma addiction phenomenon. Therefore, our study uncovers, for the first time, the prognostic value and the oncogenic function of NOTCH3 in gliomagenesis and supports NOTCH3 as a promising target of therapy in high grade glioma. Our studies allow the identification of a subset of population that may benefit from GSI-based therapies. This may lead to the design of novel strategies to improve therapeutic outcome of patients with glioma by establishing medical and scientific basis for personalized medicine.

Background: The loss of photosynthesis has occurred often in eukaryotic evolution, even more than its acquisition, which occurred at least nine times independently and which generated the evolution of the supergroups Archaeplastida, Rhizaria, Chromalveolata and Excavata. This secondary loss of autotrophic capability is essential to explain the evolution of eukaryotes and the high diversity of protists, which has been severely underestimated until recently. However, the ecological and evolutionary scenarios behind this evolutionary ‘‘step back’’ are still largely unknown. Methodology/Principal Findings: Using a dynamic model of heterotrophic and mixotrophic flagellates and two types of prey, large bacteria and ultramicrobacteria, we examine the influence of DOC concentration, mixotroph’s photosynthetic growth rate, and external limitations of photosynthesis on the coexistence of both types of flagellates. Our key premises are: large bacteria grow faster than small ones at high DOC concentrations, and vice versa; and heterotrophic flagellates are more efficient than the mixotrophs grazing small bacteria (both empirically supported). We show that differential efficiency in bacteria grazing, which strongly depends on cell size, is a key factor to explain the loss of photosynthesis in mixotrophs (which combine photosynthesis and bacterivory) leading to purely heterotrophic lineages. Further, we show in what conditions an heterotroph mutant can coexist, or even out-compete, its mixotrophic ancestor, suggesting that bacterivory and cell size reduction may have been major triggers for the diversification of eukaryotes. Conclusions/Significance: Our results suggest that, provided the mixotroph’s photosynthetic advantage is not too large, the (small) heterotroph will also dominate in nutrient-poor environments and will readily invade a community of mixotrophs and bacteria, due to its higher efficiency exploiting the ultramicrobacteria. As carbon-limited conditions were presumably widespread throughout Earth history, such a scenario may explain the numerous transitions from phototrophy to mixotrophy and further to heterotrophy within virtually all major algal lineages. We challenge prevailing concepts that affiliated the evolution of phagotrophy with eutrophic or strongly light-limited environments only.

Familial Pancreatic Cancer (FPC) is a cancer syndrome characterized by clustering of pancreatic cancer in families, but most FPC cases do not have a known genetic etiology. Understanding genetic predisposition to pancreatic cancer is important for improving screening as well as treatment. The central aim of this thesis is to identify candidate susceptibility genes for FPC, and I used three approaches of increasing resolution. First, based on a candidate-gene approach, I hypothesized that BRCA1 is inactivated by loss-of-heterozygosity in pancreatic adenocarcinoma of germline mutation carriers. I demonstrated that 5/7 pancreatic tumors from BRCA1-mutation carriers show LOH, compared to only 1/9 sporadic tumors, suggesting that BRCA1 inactivation is involved in tumorigenesis in germline mutation carriers. Second, I hypothesized that the germline genomes of FPC subjects differ in copy-number profile from healthy genomes, and that regions affected by rare deletions or duplications in FPC subjects overlap candidate tumor-suppressors or oncogenes. I found no significant difference in the global copy-number profile of FPC and control genomes, but I identified 93 copy-number variable genomic regions unique to FPC subjects, overlapping 88 genes of which several have functional roles in cancer development. I investigated one duplication to sequence the breakpoints, but I found that this duplication did not segregate with disease in the affected family. Third, I hypothesized that in a family with multiple pancreatic cancer patients, genes containing rare variants shared by the affected members constitute susceptibility genes. Using next-generation sequencing to capture most bases in coding regions of the genome, I interrogated the germline exome of three relatives who died of pancreatic cancer and a relative who is healthy at advanced age. I identified a short-list of nine candidate genes with unreported mutations shared by the three affected relatives and absent in the unaffected relative, of which a few had functional relevance to tumorigenesis. I performed Sanger sequencing to screen an unrelated cohort of approximately 70 FPC patients for mutations in the top two candidate genes, but I found no additional rare variants in those genes. In conclusion, I present a list of candidate FPC susceptibility genes for further validation and investigation in future studies.

Familial Pancreatic Cancer (FPC) is a cancer syndrome characterized by clustering of pancreatic cancer in families, but most FPC cases do not have a known genetic etiology. Understanding genetic predisposition to pancreatic cancer is important for improving screening as well as treatment. The central aim of this thesis is to identify candidate susceptibility genes for FPC, and I used three approaches of increasing resolution. First, based on a candidate-gene approach, I hypothesized that BRCA1 is inactivated by loss-of-heterozygosity in pancreatic adenocarcinoma of germline mutation carriers. I demonstrated that 5/7 pancreatic tumors from BRCA1-mutation carriers show LOH, compared to only 1/9 sporadic tumors, suggesting that BRCA1 inactivation is involved in tumorigenesis in germline mutation carriers. Second, I hypothesized that the germline genomes of FPC subjects differ in copy-number profile from healthy genomes, and that regions affected by rare deletions or duplications in FPC subjects overlap candidate tumor-suppressors or oncogenes. I found no significant difference in the global copy-number profile of FPC and control genomes, but I identified 93 copy-number variable genomic regions unique to FPC subjects, overlapping 88 genes of which several have functional roles in cancer development. I investigated one duplication to sequence the breakpoints, but I found that this duplication did not segregate with disease in the affected family. Third, I hypothesized that in a family with multiple pancreatic cancer patients, genes containing rare variants shared by the affected members constitute susceptibility genes. Using next-generation sequencing to capture most bases in coding regions of the genome, I interrogated the germline exome of three relatives who died of pancreatic cancer and a relative who is healthy at advanced age. I identified a short-list of nine candidate genes with unreported mutations shared by the three affected relatives and absent in the unaffected relative, of which a few had functional relevance to tumorigenesis. I performed Sanger sequencing to screen an unrelated cohort of approximately 70 FPC patients for mutations in the top two candidate genes, but I found no additional rare variants in those genes. In conclusion, I present a list of candidate FPC susceptibility genes for further validation and investigation in future studies.