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.

Characterizing Genetic Drivers of Lymphoma through High-Throughput Sequencing

The advent of next-generation sequencing, now nearing a decade in age, has enabled, among other capabilities, measurement of genome-wide sequence features at unprecedented scale and resolution.

In this dissertation, I describe work to understand the genetic underpinnings of non-Hodgkin’s lymphoma through exploration of the epigenetics of its cell of origin, initial characterization and interpretation of driver mutations, and finally, a larger-scale, population-level study that incorporates mutation interpretation with clinical outcome.

In the first research chapter, I describe genomic characteristics of lymphomas through the lens of their cells of origin. Just as many other cancers, such as breast cancer or lung cancer, are categorized based on their cell of origin, lymphoma subtypes can be examined through the context of their normal B Cells of origin, Naïve, Germinal Center, and post-Germinal Center. By applying integrative analysis of the epigenetics of normal B Cells of origin through chromatin-immunoprecipitation sequencing, we find that differences in normal B Cell subtypes are reflected in the mutational landscapes of the cancers that arise from them, namely Mantle Cell, Burkitt, and Diffuse Large B-Cell Lymphoma.

In the next research chapter, I describe our first endeavor into understanding the genetic heterogeneity of Diffuse Large B Cell Lymphoma, the most common form of non-Hodgkin’s lymphoma, which affects 100,000 patients in the world. Through whole-genome sequencing of 1 case as well as whole-exome sequencing of 94 cases, we characterize the most recurrent genetic features of DLBCL and lay the groundwork for a larger study.

In the last research chapter, I describe work to characterize and interpret the whole exomes of 1001 cases of DLBCL in the largest single-cancer study to date. This highly-powered study enabled sub-gene, gene-level, and gene-network level understanding of driver mutations within DLBCL. Moreover, matched genomic and clinical data enabled the connection of these driver mutations to clinical features such as treatment response or overall survival. As sequencing costs continue to drop, whole-exome sequencing will become a routine clinical assay, and another diagnostic dimension in addition to existing methods such as histology. However, to unlock the full utility of sequencing data, we must be able to interpret it. This study undertakes a first step in developing the understanding necessary to uncover the genomic signals of DLBCL hidden within its exomes. However, beyond the scope of this one disease, the experimental and analytical methods can be readily applied to other cancer sequencing studies.

Thus, this dissertation leverages next-generation sequencing analysis to understand the genetic underpinnings of lymphoma, both by examining its normal cells of origin as well as through a large-scale study to sensitively identify recurrently mutated genes and their relationship to clinical outcome.

The inoculum effect and band-pass bacterial response to periodic antibiotic treatment.

The inoculum effect (IE) refers to the decreasing efficacy of an antibiotic with increasing bacterial density. It represents a unique strategy of antibiotic tolerance and it can complicate design of effective antibiotic treatment of bacterial infections. To gain insight into this phenomenon, we have analyzed responses of a lab strain of Escherichia coli to antibiotics that target the ribosome. We show that the IE can be explained by bistable inhibition of bacterial growth. A critical requirement for this bistability is sufficiently fast degradation of ribosomes, which can result from antibiotic-induced heat-shock response. Furthermore, antibiotics that elicit the IE can lead to 'band-pass' response of bacterial growth to periodic antibiotic treatment: the treatment efficacy drastically diminishes at intermediate frequencies of treatment. Our proposed mechanism for the IE may be generally applicable to other bacterial species treated with antibiotics targeting the ribosomes.

Neutralization activity in a geographically diverse East London cohort of human immunodeficiency virus type 1-infected patients: clade C infection results in a stronger and broader humoral immune response than clade B infection.

The array of human immunodeficiency virus (HIV) subtypes encountered in East London, an area long associated with migration, is unusually heterogeneous, reflecting the diverse geographical origins of the population. In this study it was shown that viral subtypes or clades infecting a sample of HIV type 1 (HIV-1)-positive individuals in East London reflect the global pandemic. The authors studied the humoral response in 210 treatment-naïve chronically HIV-1-infected (>1 year) adult subjects against a panel of 12 viruses from six different clades. Plasmas from individuals infected with clade C, but also plasmas from clade A, and to a lesser degree clade CRF02_AG and CRF01_AE, were significantly more potent at neutralizing the tested viruses compared with plasmas from individuals infected with clade B. The difference in humoral robustness between clade C- and B-infected patients was confirmed in titration studies with an extended panel of clade B and C viruses. These results support the approach to develop an HIV-1 vaccine that includes clade C or A envelope protein (Env) immunogens for the induction of a potent neutralizing humoral response.

An integrated approach to the prediction of chemotherapeutic response in patients with breast cancer.

BACKGROUND: A major challenge in oncology is the selection of the most effective chemotherapeutic agents for individual patients, while the administration of ineffective chemotherapy increases mortality and decreases quality of life in cancer patients. This emphasizes the need to evaluate every patient's probability of responding to each chemotherapeutic agent and limiting the agents used to those most likely to be effective. METHODS AND RESULTS: Using gene expression data on the NCI-60 and corresponding drug sensitivity, mRNA and microRNA profiles were developed representing sensitivity to individual chemotherapeutic agents. The mRNA signatures were tested in an independent cohort of 133 breast cancer patients treated with the TFAC (paclitaxel, 5-fluorouracil, adriamycin, and cyclophosphamide) chemotherapy regimen. To further dissect the biology of resistance, we applied signatures of oncogenic pathway activation and performed hierarchical clustering. We then used mRNA signatures of chemotherapy sensitivity to identify alternative therapeutics for patients resistant to TFAC. Profiles from mRNA and microRNA expression data represent distinct biologic mechanisms of resistance to common cytotoxic agents. The individual mRNA signatures were validated in an independent dataset of breast tumors (P = 0.002, NPV = 82%). When the accuracy of the signatures was analyzed based on molecular variables, the predictive ability was found to be greater in basal-like than non basal-like patients (P = 0.03 and P = 0.06). Samples from patients with co-activated Myc and E2F represented the cohort with the lowest percentage (8%) of responders. Using mRNA signatures of sensitivity to other cytotoxic agents, we predict that TFAC non-responders are more likely to be sensitive to docetaxel (P = 0.04), representing a viable alternative therapy. CONCLUSIONS: Our results suggest that the optimal strategy for chemotherapy sensitivity prediction integrates molecular variables such as ER and HER2 status with corresponding microRNA and mRNA expression profiles. Importantly, we also present evidence to support the concept that analysis of molecular variables can present a rational strategy to identifying alternative therapeutic opportunities.

Gene expression signatures of radiation response are specific, durable and accurate in mice and humans.

BACKGROUND: Previous work has demonstrated the potential for peripheral blood (PB) gene expression profiling for the detection of disease or environmental exposures. METHODS AND FINDINGS: We have sought to determine the impact of several variables on the PB gene expression profile of an environmental exposure, ionizing radiation, and to determine the specificity of the PB signature of radiation versus other genotoxic stresses. Neither genotype differences nor the time of PB sampling caused any lessening of the accuracy of PB signatures to predict radiation exposure, but sex difference did influence the accuracy of the prediction of radiation exposure at the lowest level (50 cGy). A PB signature of sepsis was also generated and both the PB signature of radiation and the PB signature of sepsis were found to be 100% specific at distinguishing irradiated from septic animals. We also identified human PB signatures of radiation exposure and chemotherapy treatment which distinguished irradiated patients and chemotherapy-treated individuals within a heterogeneous population with accuracies of 90% and 81%, respectively. CONCLUSIONS: We conclude that PB gene expression profiles can be identified in mice and humans that are accurate in predicting medical conditions, are specific to each condition and remain highly accurate over time.

An integrated alcohol abuse and medical treatment model for patients with hepatitis C.

BACKGROUND: Patients with chronic hepatitis C virus (HCV) infection have high rates of alcohol consumption, which is associated with progression of fibrosis and lower response rates to HCV treatment. AIMS: This prospective cohort study examined the feasibility of a 24-week integrated alcohol and medical treatment to HCV-infected patients. METHODS: Patients were recruited from a hepatology clinic if they had an Alcohol Use Disorders Identification Test score >4 for women and >8 for men, suggesting hazardous alcohol consumption. The integrated model included patients receiving medical care and alcohol treatment within the same clinic. Alcohol treatment consisted of 6 months of group and individual therapy from an addictions specialist and consultation from a study team psychiatrist as needed. RESULTS: Sixty patients were initially enrolled, and 53 patients participated in treatment. The primary endpoint was the Addiction Severity Index (ASI) alcohol composite scores, which significantly decreased by 0.105 (41.7% reduction) between 0 and 3 months (P < 0.01) and by 0.128 (50.6% reduction) between 0 and 6 months (P < 0.01) after adjusting for covariates. Alcohol abstinence was reported by 40% of patients at 3 months and 44% at 6 months. Patients who did not become alcohol abstinent had reductions in their ASI alcohol composite scores from 0.298 at baseline to 0.219 (26.8% reduction) at 6 months (P = 0.08). CONCLUSION: This study demonstrated that an integrated model of alcohol treatment and medical care could be successfully implemented in a hepatology clinic with significant favorable impact on alcohol use and abstinence among patients with chronic HCV.

cAMP stimulates transcription of the beta 2-adrenergic receptor gene in response to short-term agonist exposure.

In addition to conveying cellular responses to an effector molecule, receptors are often themselves regulated by their effectors. We have demonstrated that epinephrine modulates both the rate of transcription of the beta 2-adrenergic receptor (beta 2AR) gene and the steady-state level of beta 2AR mRNA in DDT1MF-2 cells. Short-term (30 min) exposure to epinephrine (100 nM) stimulates the rate of beta 2AR gene transcription, resulting in a 3- to 4-fold increase in steady-state beta 2AR mRNA levels. These effects are mimicked by 1 mM N6,O2'-dibutyryladenosine 3',5'-cyclic monophosphate (Bt2cAMP) or foskolin but not by phorbol esters. The half-life of the beta 2AR mRNA after addition of actinomycin D (46.7 +/- 10.2 min; mean +/- SEM; n = 5) remained unchanged after 30 min of epinephrine treatment (46.8 +/- 10.6 min; mean +/- SEM; n = 4), indicating that a change in transcription rate is the predominant factor responsible for the increase of beta 2AR mRNA. Whereas brief exposure to epinephrine or Bt2cAMP does not significantly affect the total number of cellular beta 2ARs (assessed by ligand binding), continued exposure results in a gradual decline in beta 2AR number to approximately 20% (epinephrine) or approximately 45% (Bt2cAMP) of the levels in control cells by 24 hr. Similar decreases in agonist-stimulated adenylyl cyclase activity are observed. This loss of receptors with prolonged agonist exposure is accompanied by a 50% reduction in beta 2AR mRNA. Transfection of the beta 2AR promoter region cloned onto a reporter gene (bacterial chloramphenicol acetyltransferase) allowed demonstration of a 2- to 4-fold induction of transcription by agents that elevate cAMP levels, such as forskolin or phosphodiesterase inhibitors. These results establish the presence of elements within the proximal promoter region of the beta 2AR gene responsible for the transcriptional enhancing activity of cAMP and demonstrate that beta 2AR gene expression is regulated by a type of feedback mechanism involving the second messenger cAMP.

Rapid complete response of metastatic melanoma in a patient undergoing ipilimumab immunotherapy in the setting of active ulcerative colitis.

While blockade of the cytotoxic T-lymphocyte antigen-4 (CTLA-4) T cell regulatory receptor has become a commonly utilized strategy in the management of advanced melanoma, many questions remain regarding the use of this agent in patient populations with autoimmune disease. We present a case involving the treatment of a patient with stage IV melanoma and ulcerative colitis (UC) with anti-CTLA-4 antibody immunotherapy. Upon initial treatment, the patient developed grade III colitis requiring tumor necrosis factor-alpha (TNF-α) blocking antibody therapy, however re-treatment with anti-CTLA-4 antibody following a total colectomy resulted in a rapid complete response accompanied by the development of a tracheobronchitis, a previously described extra-intestinal manifestation of UC. This case contributes to the evolving literature on the use of checkpoint inhibitors in patients also suffering from autoimmune disease, supports future clinical trials investigating the use of these agents in patients with autoimmune diseases, and suggests that an understanding of the specific molecular pathways involved in a patient's autoimmune pathology may provide insight into the development of more effective novel combinatorial immunotherapeutic strategies.

Anti-hypotensive treatment and endothelin blockade synergistically antagonize exercise fatigue in rats under simulated high altitude.

Rapid ascent to high altitude causes illness and fatigue, and there is a demand for effective acute treatments to alleviate such effects. We hypothesized that increased oxygen delivery to the tissue using a combination of a hypertensive agent and an endothelin receptor A antagonist drugs would limit exercise-induced fatigue at simulated high altitude. Our data showed that the combination of 0.1 mg/kg ambrisentan with either 20 mg/kg ephedrine or 10 mg/kg methylphenidate significantly improved exercise duration in rats at simulated altitude of 4,267 m, whereas the individual compounds did not. In normoxic, anesthetized rats, ephedrine alone and in combination with ambrisentan increased heart rate, peripheral blood flow, carotid and pulmonary arterial pressures, breathing rate, and vastus lateralis muscle oxygenation, but under inspired hypoxia, only the combination treatment significantly enhanced muscle oxygenation. Our results suggest that sympathomimetic agents combined with endothelin-A receptor blockers offset altitude-induced fatigue in rats by synergistically increasing the delivery rate of oxygen to hypoxic muscle by concomitantly augmenting perfusion pressure and improving capillary conductance in the skeletal muscle. Our findings might therefore serve as a basis to develop an effective treatment to prevent high-altitude illness and fatigue in humans.

Inter-annual variability of precipitation constrains the production response of boreal Pinus sylvestris to nitrogen fertilization

© 2015 Published by Elsevier B.V.Tree growth resources and the efficiency of resource-use for biomass production determine the productivity of forest ecosystems. In nutrient-limited forests, nitrogen (N)-fertilization increases foliage [N], which may increase photosynthetic rates, leaf area index (L), and thus light interception (IC). The product of such changes is a higher gross primary production and higher net primary production (NPP). However, fertilization may also alter carbohydrate partitioning from below- to aboveground, increasing aboveground NPP (ANPP). We analyzed effects of long-term N-fertilization on NPP, and that of long-term carbon storing organs (NPPS) in a Pinus sylvestris forest on sandy soil, a wide-ranging forest type in the boreal region. We based our analyses on a combination of destructive harvesting, consecutive mensuration, and optical measurements of canopy openness. After eight-year fertilization with a total of 70gNm^-2, ANPP was 27±7% higher in the fertilized (F) relative to the reference (R) stand, but although L increased relative to its pre-fertilization values, IC was not greater than in R. On the seventh year after the treatment initiation, the increase of ANPP was matched by the decrease of belowground NPP (78 vs. 92gCm^-2yr^-1; ~17% of NPP) and, given the similarity of IC, suggests that the main effect of N-fertilization was changed carbon partitioning rather than increased canopy photosynthesis. Annual NPPS increased linearly with growing season temperature (T) in both treatments, with an upward shift of 70.2gCm^-2yr^-1 by fertilization, which also caused greater amount of unexplained variation (r²=0.53 in R, 0.21 in F). Residuals of the NPPS-T relationship of F were related to growing season precipitation (P, r²=0.48), indicating that T constrains productivity at this site regardless of fertility, while P is important in determining productivity where N-limitation is alleviated. We estimated that, in a growing season average T (11.5±1.0°C; 33-year-mean), NPPS response to N-fertilization will be nullified with P 31mm less than the mean (325±85mm), and would double with P 109mm greater than the mean. These results suggest that inter-annual variation in climate, particularly in P, may help explaining the reported large variability in growth responses to fertilization of pine stands on sandy soils. Furthermore, forest management of long-rotation systems, such as those of boreal and northern temperate forests, must consider the efficiency of fertilization in terms of wood production in the context of changes in climate predicted for the region.

The role of GRK6 in animal models of Parkinson's disease and L-DOPA treatment.

G protein-coupled Receptor Kinase 6 (GRK6) belongs to a family of kinases that phosphorylate GPCRs. GRK6 levels were found to be altered in Parkinson's Disease (PD) and D(2) dopamine receptors are supersensitive in mice lacking GRK6 (GRK6-KO mice). To understand how GRK6 modulates the behavioral manifestations of dopamine deficiency and responses to L-DOPA, we used three approaches to model PD in GRK6-KO mice: 1) the cataleptic response to haloperidol; 2) introducing GRK6 mutation to an acute model of absolute dopamine deficiency, DDD mice; 3) hemiparkinsonian 6-OHDA model. Furthermore, dopamine-related striatal signaling was analyzed by assessing the phosphorylation of AKT/GSK3β and ERK1/2. GRK6 deficiency reduced cataleptic behavior, potentiated the acute effect of L-DOPA in DDD mice, reduced rotational behavior in hemi-parkinsonian mice, and reduced abnormal involuntary movements induced by chronic L-DOPA. These data indicate that approaches to regulate GRK6 activity could be useful in modulating both therapeutic and side-effects of L-DOPA.

Adjunctive albuterol enhances the response to enzyme replacement therapy in late-onset Pompe disease.

Effective dosages for enzyme replacement therapy (ERT) in Pompe disease are much higher than for other lysosomal storage disorders, which has been attributed to low cation-independent mannose-6-phosphate receptor (CI-MPR) in skeletal muscle. We have previously demonstrated the benefit of increased CI-MPR-mediated uptake of recombinant human acid-α-glucosidase during ERT in mice with Pompe disease following addition of albuterol therapy. Currently we have completed a pilot study of albuterol in patients with late-onset Pompe disease already on ERT for >2 yr, who were not improving further. The 6-min walk test (6MWT) distance increased in all 7 subjects at wk 6 (30±13 m; P=0.002), wk 12 (34±14 m; P=0.004), and wk 24 (42±37 m; P=0.02), in comparison with baseline. Grip strength was improved significantly for both hands at wk 12. Furthermore, individual subjects reported benefits; e.g., a female patient could stand up from sitting on the floor much more easily (time for supine to standing position decreased from 30 to 11 s), and a male patient could readily swing his legs out of his van seat (hip abduction increased from 1 to 2+ on manual muscle testing). Finally, analysis of the quadriceps biopsies suggested increased CI-MPR at wk 12 (P=0.08), compared with baseline. With the exception of 1 patient who succumbed to respiratory complications of Pompe disease in the first week, only mild adverse events have been reported, including tremor, transient difficulty falling asleep, and mild urinary retention (requiring early morning voiding). Therefore, this pilot study revealed initial safety and efficacy in an open label study of adjunctive albuterol therapy in patients with late-onset Pompe disease who had been stable on ERT with no improvements noted over the previous several years.