Next Article in Journal
Directed Evolution of P450 BM3 towards Functionalization of Aromatic O-Heterocycles
Next Article in Special Issue
Transplacental Gene Delivery (TPGD) as a Noninvasive Tool for Fetal Gene Manipulation in Mice
Previous Article in Journal
Diagnostic Roles of Postmortem cTn I and cTn T in Cardiac Death with Special Regard to Myocardial Infarction: A Systematic Literature Review and Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantifying the Potential for Future Gene Therapy to Lower Lifetime Risk of Polygenic Late-Onset Diseases

by
Roman Teo Oliynyk
1,2
1
Centre for Computational Evolution, University of Auckland, Auckland 1010, New Zealand
2
Department of Computer Science, University of Auckland, Auckland 1010, New Zealand
Int. J. Mol. Sci. 2019, 20(13), 3352; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133352
Submission received: 21 June 2019 / Revised: 5 July 2019 / Accepted: 5 July 2019 / Published: 8 July 2019
(This article belongs to the Special Issue Genome Editing Therapies)

Abstract

:
Gene therapy techniques and genetic knowledge may sufficiently advance, within the next few decades, to support prophylactic gene therapy for the prevention of polygenic late-onset diseases. The risk of these diseases may, hypothetically, be lowered by correcting the effects of a subset of common low effect gene variants. In this paper, simulations show that if such gene therapy were to become technically possible; and if the incidences of the treated diseases follow the proportional hazards model with a multiplicative genetic architecture composed of a sufficient number of common small effect gene variants, then: (a) late-onset diseases with the highest familial heritability will have the largest number of variants available for editing; (b) diseases that currently have the highest lifetime risk, particularly those with the highest incidence rate continuing into older ages, will prove the most challenging cases in lowering lifetime risk and delaying the age of onset at a population-wide level; (c) diseases that are characterized by the lowest lifetime risk will show the strongest and longest-lasting response to such therapies; and (d) longer life expectancy is associated with a higher lifetime risk of these diseases, and this tendency, while delayed, will continue after therapy.

1. Introduction

In the past two decades, the human genome has been successfully sequenced. Whole genome sequencing (WGS) and genome-wide association studies (GWASs) of human genomes (as well as those of other organisms) have become an everyday occurrence [1]. Our knowledge of genetic variants, particularly the single nucleotide polymorphisms (SNPs) associated with susceptibility to diseases, has become deeper and more extensive.
Experimental gene therapy techniques, aimed at diseases caused by a single defective gene or a single SNP—the so-called Mendelian conditions—are being refined. Mendelian conditions cause high mortality and morbidity, but each of these conditions affects only a minute fraction of the population. As of June 2019, the OMIM Gene Map Statistics [2] compendium has listed 6436 phenotypic genetic conditions caused by 4102 gene mutations. This list includes a variety of conditions, with onsets ranging from very early to late. For example, type 1 diabetes mellitus is caused by single defects in the HLA-DQA1, HLA-DQB1, or HLA-DRB1 genes [3]. Early-onset Alzheimer’s disease is caused primarily by APP, PSEN1, or PSEN2 gene mutations and affects a relatively small proportion of the population, starting in their thirties, with the majority of mutation carriers being affected by the age of 65 [4]. In contrast, macular degeneration [5,6,7] is primarily caused by a small number of high-effect variants and manifests at a relatively old age. In some cases, individualized genetic diagnoses, where an SNP that needs to be edited can be specified precisely, are possible. Over the last two decades, 287 monogenic disease clinical trials have been conducted worldwide [8]. When the medical technology becomes available, individuals who receive treatment will be effectively cured and will have no need for concern about the single specific cause of their disease.
Polygenic or complex late-onset diseases (LODs) pose a more nuanced problem, and this study will focus on them. There are thousands of estimated gene variants or SNPs of typically small effect that, in combination, constitute the polygenic LOD risk of each individual [9,10]. These diseases include the old-age diseases that eventually affect most individuals and are exemplified by cardiovascular disease (particularly coronary artery disease (CAD)), cerebral stroke, type 2 diabetes (T2D), senile dementia, Alzheimer’s disease (AD), cancer, and osteoarthritis.
What distinguishes polygenic LODs from infectious diseases or from Mendelian genetic conditions is difficulty in terms of the concept of cure. The diseases of aging are primarily a consequence of an organism’s decline over time, leading to increased susceptibility to many LODs [11,12,13]. The combination of genetic liability, environmental factors, and the physiological decline of multiple organ systems leads to individual disease presentations [14]. Detrimental gene variants are exacerbating factors [15], compared to the average distribution of common gene variants that define human conditions, as they apply to polygenic LODs. The time of onset for each individual is modulated by genotype and environment [16]. While some individuals will be diagnosed at a relatively young age, others will not be diagnosed with a particular LOD during their lifetime [17]. According to the current consensus, a large number of common low-effect variants offer the likeliest explanation for the heritability of the majority of complex traits [18,19]. For example, in the cancers analyzed in this study, the fraction of all diagnoses that were attributed to highly detrimental inherited mutations was relatively low—it was estimated to explain heritability connected with 10%–14% of breast cancers [20,21], 10%–12% of prostate cancers [22,23,24,25], 5%–10% of colorectal cancers [26,27], and was assumed to be a relatively minor fraction for lung cancers [28,29,30]. For the majority of these cancers, liability is attributed to the common low-effect gene variants and environmental factors. The development of cancer is a multistage process, wherein individual variability in any tumorigenesis stage duration or liability may be influenced by hereditary predisposition, as well as environmental factors [31]. The level of susceptibility to the major polygenic LODs, and the difference between high-risk and low-risk individuals, may lie in a slightly higher- or lower-than-average fraction of detrimental gene variants. Certainly, the failure does not begin immediately prior to the age of diagnosis. For example, AD deterioration begins decades before symptoms first become noticeable [32]. A similar situation holds for cardiovascular disease [33,34] and cancer [35].
The best cure is prevention, and the time may be nearing when prophylactic gene therapy will be attempted for the prevention of complex polygenic diseases. Much scientific knowledge and technical expertise is required and many ethics questions will need to be settled before this level of prophylactic gene therapy can become possible. From an ethical perspective, as techniques have developed and the medical possibilities offered by gene therapy for improving health and preventing diseases have gradually materialized, its acceptance is becoming more widespread. This is exemplified by the findings of the U.S. Committee of the National Academies of Sciences, Engineering, and Medicine [36] in Human genome editing—Science, ethics, and governance, and the recommendations of the U.K. Nuffield Council on Bioethics [37] in Genome editing and human reproduction: Social and ethical issues, which considered germline editing as one possible application.
Computational techniques attempting to evaluate the effects of mutations or gene variants have been developed, although their accuracy needs to improve dramatically before they can become applicable to personalized human genetic evaluation or treatment [38]. Similarly, while extensive libraries of human SNPs have been compiled, including dbSNP, HapMap, SNPedia, and aggregating sites [39], the information is far from actionable as far as modifying multiple personalized SNPs is concerned. The ability to locate or be able to computationally estimate a complete set of the low-effect causal SNPs requires knowledge that may take decades to gain.
Gene editing technologies may also be a few decades away from the time when they can be used routinely, with the same low risk as applying an influenza vaccination, to modify a large number of gene variants distributed across the human genome. The latest gene editing technique, CRISPR-Cas9 [40], has supplemented and mostly replaced older technologies, such as zinc-finger nuclease (ZFN) [41] and transcription activator-like effector nuclease (TALEN) [42], although, for some applications, these older techniques continue to be more appropriate. While its selectivity and on-target precision have improved, CRISPR is still the most effective in gene knockdown operations. For modification and repair, only a small fraction of CRISPR operations—using homologous repair with a template or a sister chromatid sequence—succeed. A recent advance, reported by Smith et al. [43], proposed base editing with reduced DNA nicking, allowing for the simultaneous editing of >10,000 loci in human cells. CRISPR, which is only six years old, remains to be a rapidly developing technology that holds great promise. Synthetic genomics [44,45] could be another promising future technology. Synthetic genomics techniques could also help in developing the precise mapping of the effects of gene variants on disease phenotypes. If none of these approaches ultimately succeed in becoming reliable enough for the purposes of gene therapy, it is almost certain that a new, more suitable technique will be invented.
Changes in lifestyle and medical care, including the prevention and treatment of infectious diseases, have extended longevity over the last century, and this trend is projected to continue. This increased longevity is partly due to medical advances, helping people to live and function decades after first being diagnosed with historically deadly or debilitating illnesses. Preventive gene therapies may also become a future factor in prolonging health span. Actuarial science has tracked human mortality trends for centuries. The Gompertz–Makeham law of mortality, which was established more than 150 years ago, depicts an exponential increase in the rate of human mortality after the age of 30 [46,47]. While the parameters of the Gompertz–Makeham law continue to be adjusted, the principle remains valid. The apparent squaring of the mortality curve—the so-called compression of morbidity and mortality into older ages—implies that the maximum human lifespan is likely limited to about 120 years of age [48,49,50].
Within the next few decades, gene therapy techniques and genetic knowledge may sufficiently advance to support prophylactic gene therapy to prevent late-onset diseases. It may be timely to evaluate the extent of the effects that future gene therapies may have on delaying the onset of LODs or preventing them entirely.
The goal of this study is to establish how the proportional hazards model and multiplicative genetic architecture can be used to map the polygenic risk to hazard ratio of succumbing to common late-onset diseases with advancing age and apply this mapping to quantify the effects of hypothetical future prophylactic gene therapies. As its foundation, this study used earlier research [51], which reviewed epidemiology, heritability, and polygenic risk models, and developed a simulational basis for the analysis of eight of the most common diseases: AD; T2D; cerebral stroke; CA; and breast, prostate, colorectal, and lung cancers. Computer simulations in this study quantified the correlation between the aging process, the polygenic risk score (PRS), and the change in the hazard ratio with age—using as inputs the clinical incidence rate and familial heritability—and estimated the outcomes of hypothetical future prophylactic gene therapy on the lifetime risk and age of onset for these eight LODs, they also estimated the lifetime risk increase associated with longevity gains.

2. Results

2.1. Characteristics of the Aging Coefficient

The clinical incidence rate pattern was used to map the hazard ratio from the PRS of individuals diagnosed with an LOD at each year of age for eight LODs: AD; T2D; cerebral stroke; CAD; and breast, prostate, colorectal, and lung late-onset cancers. Knowledge of the hazard ratio for any given PRS value at any given age (the aging coefficient) is what is needed for further calculations and simulations modeling gene therapy and the effects of life expectancy increases. The aging coefficients were discovered for these LODs, as described in the Methods section “Aging coefficient discovery simulation steps”, and can be seen in Figure 1.
It is interesting to note how the magnitude of the aging coefficient changes with age for the analyzed diseases. The range of values spanned by the aging coefficient and its maximum value, seen in Figure 1A, is larger for AD than for all other analyzed LODs. This can be ascribed to two factors that distinguish AD: (1) The steepest rise of incidence rate and cumulative incidence; and (2) the highest heritability of the reviewed LODs, resulting in the highest variance according to our genetic architecture model. This larger variance results in more extreme values of high and low PRSs in the population and, therefore, the age coefficient multiplier is necessarily the smallest at young onset ages; with the opposite being true at older ages. T2D, CAD, and stroke show a comparatively moderate progression and maximum projected aging coefficient values, while cancers still show a smaller maximum magnitude. Lung cancer possesses both the lowest heritability and the lowest cumulative incidence and, consequently, the values discovered by simulation for the lung cancer aging coefficient (seen in Figure 1H) almost precisely match its incidence rate. In a limit case of an LOD characterized by a PRS that remains constant with age and a very low incidence rate, the aging coefficient should be identical to the incidence rate.

2.2. Longer Life Expectancy Corresponds to Increasing Lifetime Risk

The modeled increased life expectancy in the baseline scenario, without prophylactic gene therapy, is displayed in Figure 2. This analysis corresponds to the baseline incidence rate, represented by the blue line in Figure 3. All analyzed LODs show a significant lifetime risk increase with every five years of life expectancy extension. This is most prominent with AD, the incidence of which nearly triples with an extension of life expectancy of 15 years, while the incidence rate less than doubles for the rest of the LODs. The incidence rate density (as seen in Figure S1) shows the relative incidence increasing while the peak of incidence shifts toward older ages.

2.3. Lifetime Risk Estimates for Discrete Hazard Ratio Multiples

Having the aging coefficient, it is easy to recalculate lifetime incidence risk for a range of hazard ratios (HRs). Population mortality is one of the principal limiting factors on the lifetime risk of an LOD, and shifting the mortality curve (and, thus, emulating longer life expectancy) reveals how incidence rate and lifetime risk would change if life expectancy were to increase.
Figure S2 shows a grid display with HR ranging from 16.0 to 0.0625 and average life expectancy varying from baseline to extended up to 15 years. The results of these calculations show that the lifetime risk is proportionate to the HR, as long as lifetime risk is relatively low (below 50%) for all reviewed LODs, even though they display quite varied incidence patterns. Lowering the HR, through gene therapy or other means, implies proportionately lowering the lifetime risk, assuming life expectancy remains constant, as summarized in Table 1. For instance, lowering the HR four-fold is accompanied by a four-fold (or 400%) drop in lifetime risk. The only exception is AD, for which the lifetime risk decreases at a slower-than-proportionate rate—a possible explanation could be the near-exponential incidence rate increase to very advanced ages for AD, while all other LODs analyzed here can be approximated by a flattened or linear progression at a more moderate old age [51]. Figure S2 shows a dramatic increase in the lifetime risk of LODs with every five extra years of life expectancy. Evidently, with increased life expectancy, the projection approaches certainty for higher HRs and increasing life expectancies. This is most prominent for AD and T2D, which can be explained by the high heritability and high later-age prevalence of these LODs, leading to a significantly higher risk that individuals will become ill earlier. Furthermore, at an advanced age, the remaining lower-risk individuals are those who constitute the majority of incidence cases [51]. As a result, the modeled high-PRS individuals show improvement in delaying the statistical disease onset, while the lifetime risk may remain almost as high.
The application of the aging coefficient and life expectancy increases makes it simple to estimate the onset delay if the PRS were changed, as in the case of prophylactic gene therapy. Table 2 and Table 3 show the values of shifts in the onset delay on the cumulative incidence slope, at 30% of the lifetime risk and for the full lifetime risk, respectively. Emulating longer life expectancy raises lifetime risk, and we see a more complex picture than a mere proportionate decline in risk, depending on the incidence and heritability patterns for each LOD. The AD lifetime risk exceeded the baseline within approximately 1–4 years of longer life expectancy. It would take approximately 10–15 years of longer life for T2D, stroke, and CAD to approach or slightly exceed their respective baseline lifetime risks. In nearly all calculated scenarios, it took significantly longer than 15 years for cancers to exceed their baseline lifetime risks.
Reviewing the first four LODs in the first row and, to a lesser extent, AD and T2D in the second row of Table 3, it is noticeable that for high PRS values (i.e., 16.0 and 4.0) the response was lower than for the final row and also lower than for the slope onset delay in Table 2. Figure S2 shows that the likelihood of becoming ill is modeled as a near certainty for these values. The simulation results of aggregating outcomes for the population in the following section provide additional generalization and confirmation of these patterns.

2.4. Results of Simulated Gene Therapy Lowering Population PRS

The results of this simulation show an improved (lower) incidence rate after treatment, as is shown in Figure 3; and a corresponding decrease in lifetime risk, as presented in Figure 4. Figure 4 also serves as a good qualitative illustration of the results of the previous section. Comparing the baseline (blue line) to therapy, with life expectancy (green line) held constant, shows a significant improvement in both the lifetime risk and age of onset for all LODs. The lifetime risk is lower for all LODs and shows a delay of approximately a decade in the incidence rate curve for T2D, stroke, and CAD. AD benefited the least, and cancers showed the most significant improvement. Table 4 illuminates the results detailed earlier from a slightly different perspective. With life expectancy unchanged, the lifetime risk decreased by 30% or more for AD and T2D, by more than 66% for colorectal and lung cancer, and by 50% or more for the remaining four LODs.
Emulating longer life expectancy resulted in increased lifetime risk. The AD lifetime risk exceeded the baseline within approximately 3 years of longer life expectancy. It took 15 years of longer life for T2D, stroke, and CAD to approach or slightly exceed their respective baseline lifetime risks. Nevertheless, even with longer life expectancies, the onset remained delayed, compared to the baseline. All cancers remained far below their baseline lifetime risks, even with a 15-year longer life expectancy. The risk remained lower than the baseline for all analyzed cancers after an increase of 15 years in life expectancy, any LOD with a relatively high heritability and low prevalence should similarly benefit. The results are summarized in Table 4, along with relevant LOD statistics.

3. Discussion

For the purposes of this hypothetical treatise, it was assumed that it is possible to precisely identify individual gene variants and their detrimental or beneficial effects, then, use gene therapy to modify a large number of detrimental variants. Rather than analyzing arbitrary synthetic choices of heritability and disease incidence progressions, eight LODs were chosen as a case study. Using this approach allowed us to relate the findings to some of the highly prevalent LODs that cover the broad spectrum of heritability and disease incidence patterns and—while keeping in mind that the results are a model view, with each of the reviewed LODs certainly possessing deeper specific causal mechanisms—it allowed us to make generalizations about lifetime risk changes if the LOD risks were lowered by some intervention, in this case, by gene therapy. This hypothetical gene therapy model was applied to estimate what would happen to LOD progression as the population ages. Conceptually, gene therapy here does not consider additions of artificially designed genomic sequences, but rather, only corrections made to typically low-effect heterozygous in-population gene variants, that is, a correction of a detrimental variant to a naturally occurring neutral state. For the sake of simplicity, the model used SNP distributions, though the same would apply (albeit with a higher degree of complexity) to gene therapy using other gene variant types.
This study does not evaluate potential obstacles due to pleiotropy, defined for the purposes of gene therapy as the possible negative effect on other phenotypic features of any attempt to prevent an LOD by modifying a subset of SNPs [52,53]. The high-risk individual PRS is caused by numerous variants. In this model, these are normally distributed in the population. There is a relatively small difference in the absolute number of detrimental alleles between the population average and higher-risk individuals. Arguably, for the purpose of personalized prophylactic treatment, it will be possible to select a small fraction of variants from a large set of available choices (as seen in Table 4) that do not possess antagonistic pleiotropy, or perhaps even select SNPs that are agonistically pleiotropic with regard to some of the other LODs.
Applying the modeled aging coefficient to evaluate the impact of longer life expectancy on lifetime risk confirms the long-standing observation that aging itself is the predominant risk factor for many late-onset diseases and conditions. The calculations applying the discovered aging coefficient to the discrete hazard ratio values showed a delay in onset incidence for all analyzed LODs. The lifetime risk decreased in proportion with a decrease in hazard ratio, as long as the absolute value of lifetime risk remained low. With the introduction of an emulated life expectancy increase, the lifetime risk increased. The lifetime risk increase with age was most prominent for AD. In those countries with longer life expectancy, the lifetime risk of AD is usually higher, as was demonstrated by Wu et al. [54], using the example of Japan. These results confirm, once more, that if mortality from all causes is lower (resulting in a longer life expectancy), AD is an LOD that exhibits a rapid rise in advanced-age prevalence. It would be difficult to limit the prevalence of AD, which is delayed only by approximately 3 years with the modeled level of therapy. AD may require a higher number of gene edits, likely postponing the possibility of more effective treatment to a point even further in the future; yet, any improvement would be welcome. It is possible that a pharmaceutical intervention targeting a causal metabolic pathway or immune or inflammatory response may be more effective for AD, although past announcements that generated false hope regarding breakthroughs through these kinds of approaches are too numerous to cite.
The Framingham General Cardiovascular Risk Score included age as one of the major risk factors for stroke and CAD [55]. Boehme et al. [56] showed a similar pattern for T2D, which the results of the current study were in agreement with. For T2D, stroke, and CAD, lifetime risk will regain pretreatment baselines within 10 to 15 years of longer life, which is equivalent to delaying the average onset age of these LODs by as many years. Based on heritability and incidence rate combinations, prophylactic gene therapy holds the potential to bring significant and longer-lasting benefits for cancer prevention, even with a similar or smaller number of edited gene variants than for the more prevalent diseases. The potential limitation of this study is the possibility that GWAS (and other future techniques) will have difficulties in finding a sufficient number of common low-effect SNPs to decrease the disease liability to the level simulated in this research, or that gene-environment effects will not follow Cox’s proportional hazards model [57,58] for some of the late-onset polygenic diseases. The likeliest candidate is lung cancer, which has the lowest heritability and is the most environmentally affected of all cancers reviewed here. For lung cancer, addressing the polygenic risk of smoking [59], as well as genetically influenced carcinogenicity of smoking on an individual level [60] and environmental improvements may allow for similar amelioration of disease liability. Additionally, when such advanced gene therapy technologies become available, preventing monogenic, highly detrimental variants will be simple, and the combination of therapies can bring about even more substantial improvements in both individual and population-wide health outcomes.
Gene therapy simulation scenarios analyzing population statistics showed decreases in LOD incidence and delays in LOD onset. These simulations also showed the increase in lifetime risk with emulated longer life expectancy. Such estimates may be important for evaluating population health and well-being and the potential financial impact on healthcare systems. The estimates in this study, based on the proportional hazards model and multiplicative genetic architecture using the aging coefficient, allowed for an estimation of these effects accounting for a model genetic architecture of the LODs, rather than a more simplistic calculation based primarily on the statistical shape of the incidence rate progression. In a study, aptly titled “Projections of Alzheimer’s disease in the United States and the public health impact of delaying disease onset," Brookmeyer et al. [61], it was estimated that an intervention that achieved a two-fold AD hazard ratio decrease would shift the exponential rise curve of AD by five years, leading, in the long term, to a twofold decline in the cumulative incidence and prevalence of AD when accounting for mortality. The simulation reflecting age-related change in PRS distribution demonstrated that the positive effect on the lifetime risk of AD would be significantly lower than projected by the above study, in the case of preventative gene therapy. While AD has emerged as one of the most difficult diseases to prevent, LODs with low cumulative incidence, such as cancer, exhibit enduring improvement under this model.
Even though each LOD was analyzed independently in this study, prioritizing certain LODs for preventative therapy, in practice, could have a significant effect on other conditions not specifically targeted for treatment. For example, T2D is one of the diseases that causes the most comorbidities, accelerating the onset of cardiovascular and other diseases, sometimes by decades [56]. For this reason, preventative treatment of T2D could mean improvements in health or delays in the presentation of a range of LODs, either independently of or in addition to treating their specific gene variants.

4. Methods

4.1. Conceptual Summary

This section briefly summarizes the concepts of the earlier research [51] used as a foundation for this study. The subsequent sections will describe this study’s simulation flow and implementation.

4.1.1. Cox’s Proportional Hazards Model

According to Chatterjee et al. [58], the conditional age-specific incidence rate of the disease, I ( t | G ) , defined as the probability of developing the disease at a particular age t given that a subject has been disease-free until that age, can be modeled using Cox’s proportional hazards model [57]:
I ( t | G ) = I 0 ( t ) · e x p ( k b k G k ) ,
where G = ( G 1 , , G k ) is the multiplicative effect of a set of risk factors on the baseline hazard of the disease I 0 ( t ) . The set of age-independent variables in G could include genetic and environmental risk factors, as well as their interaction terms.
According to Chatterjee et al. [58], if it can be assumed that environmental risk factors are independent of the SNPs, the “post-GWAS epidemiological studies of gene-environment interactions have generally reported multiplicative joint associations between low-penetrant SNPs and environmental risk factors, with only a few exceptions”. This means that the polygenic score G = k b k G k , as the lifelong characteristic of each individual, is used multiplicatively with I 0 ( t ) , which encompasses environmental and aging effects.
The simulations in this study operate on model genetic architectures of the analyzed LODs, not a complete GWAS map of their experimentally discovered SNPs, because GWAS-discovered sets can explain only a fraction of the heritability of these LODs. These model genetic architecture SNPs are treated as “true” causal for disease liability and heritability variants, as discussed in Chatterjee et al. [58], contrary to GWAS-linked SNPs—and it is assumed that they can be accurately identified for the purposes of personalized gene therapy.

4.1.2. Allele Distribution Models

The allele distribution models were based on [62,63]. The allele scenarios were implemented in the simulations identically to those in the earlier research, see [51] for a comprehensive description.
The common allele low effect size genetic architecture model (Scenario A, Table 5) was expected to be most suitable for explaining the heritability of the analyzed LODs, and all results were reported using this scenario. For verification, all simulations were also performed using a more extreme rare allele medium effect size model (Scenario B, Table 5), and the results for lifetime risk were essentially identical to those in Scenario A. Any material difference in the results would have warranted additional investigation. As the results did not differ materially, separate figures are not presented for Scenario B—the corresponding simulation data are available in Supplementary Data. The number of variants needed for the Scenario A LODs is summarized in the Results section.

4.1.3. LOD Incidence Rate Functional Approximation

The incidence rate functional approximations and source data are detailed in Oliynyk [51]. The yearly incidence rate logistic and exponential regressions for the functional approximation of the LODs, based on the available clinical incidence statistics, are illustrated in Figure S3. The logistic approximations were used for all LODs except breast cancer, for which the exponential followed by linear regression more accurately approximated the incidence rate pattern.

4.2. The Aging Coefficient: Mapping PRS to Age-Dependent Probability of LODs

Earlier research [51] showed that the incidence rates for these LODs increase as individuals age and, statistically, the PRSs of individuals diagnosed at older ages decline. Following the multiplicative model of PRS and environmental effects, from a statistical perspective, two processes are balancing disease incidence—the aging or environmental effect increases the average susceptibility of the population, and this makes it more likely that older individuals will become ill. Individuals with higher PRSs are more likely to become ill and, with higher incidence rates, a larger fraction of these individuals have already become ill with every additional year of age and are therefore excluded from the high-risk pool for the following years of age. The aging process continues, as reflected in the incidence rate rising or remaining high in subsequent years.
Statistically, individual PRSs and environmental effects including aging affect the probability of a person becoming ill [58]. According to Cox’s proportional hazards model [57], in every year of age, each individual can be assigned a hazard ratio which describes that individual’s risk of being diagnosed with an LOD. The goal is to uncover the mapping between PRSs and individual hazard ratios on a yearly basis:
R u ( t ) A ( t ) · G u ,
where R u ( a g e ) is the hazard ratio of uth unaffected individual, and t is the individual’s age in years. The PRS G u remains constant for each individual for life, and the multiplier A ( t ) drives the age-related increase in LOD incidence— A ( t ) will also be called the “aging coefficient”.
It may be possible to map or “discover” A ( t ) by applying the yearly incidence rate in a population to the PRS distribution, based on an LOD genetic allele architecture, through the population simulation flow described in the next section. The discovered A ( t ) , then applied to the population simulation on a yearly basis, should precisely reproduce the initial LOD incidence rate pattern. Later, building a population with a modified PRS to emulate the effect of gene therapy and simulating aging of this population by applying A ( t ) , it will be possible to find the corresponding incidence rate and other resulting statistics.
In addition to the age-related change in the PRS distribution described above, the aging coefficient automatically incorporates other miscellaneous environmental effects that accumulate with age, as these are all reflected in the yearly incidence rate. The term “aging coefficient” is used throughout this publication, rather than “environmental effect”, to emphasize that it is an aggregate age-dependent parameter.
A useful parameter in the simulation and analysis is the incidence rate density D ( t ) , which depicts LOD incidence contributions relative to the initial population count at yearly increments, adjusted for mortality:
D ( t ) = I ( t ) · S ( t ) ,
where I ( t ) is the yearly incidence at age t, and S ( t ) is the survivor rate from the US Social Security Actuarial Life Table [64]. Integrating the area under the curve, or summing up discrete yearly values, corresponds to an LOD cumulative incidence C ( t ) —the limit or lifetime value C ( T ) equals the lifetime risk:
C ( t ) = t T D ( t ) .
The comprehensive description of the simulation steps, procedures, validation, and statistical analysis is available in Supplementary Materials Chapter S1.

4.3. Data Sources, Programming, and Equipment

The population mortality estimates from the [64] provide annual death probability and survivor numbers, up to 119 years of age, for both men and women.
Disease incidence data were obtained from the following sources: Alzheimer’s disease [61,65,66,67], type 2 diabetes [56], coronary artery disease and cerebral stroke [68], and cancers [69,70].
The simulations were performed on an Intel Xeon Gold 6154 CPU-based 36-core computer system with 288GB of RAM. The simulation was written in C++, and the source code can be found in Supplementary Data.
The final simulation data, additional plots, R scripts, and executables are also available in Supplementary Data. Intel Parallel Studio XE was used for multithreading support and the Boost C++ library for faster statistical functions—the executable can be built and function without these two libraries, with a corresponding slowdown in execution.

5. Conclusions

In this study, computer simulations mapped polygenic risk to the hazard ratio of being diagnosed with eight common LODs, based on their known heritability and incidence rates, under the proportional hazards model and multiplicative genetic architecture. The resulting mapping—the aging coefficient—enabled the researcher to quantify the population effects of the emulated prophylactic gene therapy, alongside longevity increases. Computer modeling and simulations deal with simplifications and generalizations of biological processes, and aim to make predictions about the behavior of the modeled systems when modifying parameters of a model, the conclusions of this study are made in such context. The conclusions of this study are contingent on progress in molecular genetics identifying a sufficient number of true causal SNPs for a particular LOD on an individual basis, gene editing technologies becoming capable to safely provide such a level of therapy, and prophylactic gene therapies successfully passing clinical trials and obtaining the approval of governmental agencies.
The intensive gene therapy simulated here could dramatically delay the average onset of the analyzed LODs and reduce the lifetime risk of the population. The simulations highlighted that the magnitude of familial heritability and cumulative incidence patterns distinguish the outcomes for the analyzed LODs when subjected to the same PRS decrease. This outcome can be characterized by the delay in LOD onset, that is, the estimate of the number of years it would take for each LOD to regain the pretreatment baseline level.
In summary, if gene therapy, as hypothesized here, were to become possible, and if the incidence of the treated diseases followed the proportional hazards model with multiplicative genetic architecture composed of a sufficient number of common low effect gene variants, then (a) late-onset diseases with the highest familial heritability will have the highest number of variants available for editing; (b) diseases with the highest current lifetime risk, particularly those with the highest incidence rate continuing into advanced age, will be the most resistant to attempts to lower the lifetime risk and delay the age of onset at a population level; (c) diseases that are characterized by the lowest lifetime risk will show the strongest and longest-lasting response to such therapies; and (d) longer life expectancy is associated with a higher lifetime risk of these diseases, and this tendency, while delayed, will continue after the therapy.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/13/3352/s1, Supplementary Document: A document Supplementary.PDF containing a supplementary methods chapter and figures. Supplementary Data: A zip file SupplementaryData.ZIP containing the simulation executable, the source code, R scripts, batch files, and simulation results.

Funding

This research received no external funding.

Acknowledgments

The author thanks Alexei J. Drummond at the University of Auckland for a number of helpful and challenging discussions.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ADAlzheimer’s disease
CADcoronary artery disease
GWASgenome-wide association study
HRhazard ratio
LODlate-onset disease
MAFminor allele frequency; customarily implying the ’effect allele frequency’
ORodds ratio
PRSpolygenic risk score
SNPsingle nucleotide polymorphism; in context of this study used synonymously with the term ’allele’
T2Dtype 2 diabetes

References

  1. Visscher, P.M.; Wray, N.R.; Zhang, Q.; Sklar, P.; McCarthy, M.I.; Brown, M.A.; Yang, J. 10 years of GWAS discovery: Biology, function, and translation. Am. J. Hum. Genet. 2017, 101, 5–22. [Google Scholar] [PubMed]
  2. OMIM. 2019. Available online: http://omim.org/statistics/geneMap (accessed on 2 June 2019).
  3. Steck, A.K.; Rewers, M.J. Genetics of type 1 diabetes. Clin. Chem. 2011, 57, 176–185. [Google Scholar] [CrossRef] [PubMed]
  4. Ghani, M.; Reitz, C.; St George-Hyslop, P.; Rogaeva, E. Genetic Complexity of Early-Onset Alzheimer’s Disease. In Neurodegenerative Diseases; Springer International Publishing: Cham, Switzerland, 2018; pp. 29–50. [Google Scholar]
  5. Jager, R.D.; Mieler, W.F.; Miller, J.W. Age-related macular degeneration. N. Engl. J. Med. 2008, 358, 2606–2617. [Google Scholar] [CrossRef] [PubMed]
  6. Manolio, T.A.; Collins, F.S.; Cox, N.J.; Goldstein, D.B.; Hindorff, L.A.; Hunter, D.J.; McCarthy, M.I.; Ramos, E.M.; Cardon, L.R.; Chakravarti, A. Finding the missing heritability of complex diseases. Nature 2009, 461, 747–753. [Google Scholar] [PubMed] [Green Version]
  7. Sobrin, L.; Ripke, S.; Yu, Y.; Fagerness, J.; Bhangale, T.R.; Tan, P.L.; Souied, E.H.; Buitendijk, G.H.; Merriam, J.E.; Richardson, A.J. Heritability and genome-wide association study to assess genetic differences between advanced age-related macular degeneration subtypes. Ophthalmology 2012, 119, 1874–1885. [Google Scholar] [PubMed]
  8. Ginn, S.L.; Amaya, A.K.; Alexander, I.E.; Edelstein, M.; Abedi, M.R. Gene therapy clinical trials worldwide to 2017: An update. J. Gene Med. 2018, 20, e3015. [Google Scholar] [PubMed]
  9. Eyre-Walker, A. Genetic architecture of a complex trait and its implications for fitness and genome-wide association studies. Proc. Natl. Acad. Sci. USA 2010, 107, 1752–1756. [Google Scholar]
  10. Yang, J.; Ferreira, T.; Morris, A.P.; Medland, S.E.; Madden, P.A.; Heath, A.C.; Martin, N.G.; Montgomery, G.W.; Weedon, M.N.; Loos, R.J. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat. Genet. 2012, 44, 369–375. [Google Scholar] [CrossRef]
  11. Ribezzo, F.; Shiloh, Y.; Schumacher, B. Systemic DNA damage responses in aging and diseases. In Seminars in Cancer Biology; Elsevier: Atlanta, GA, USA, 2016; Volume 37, pp. 26–35. [Google Scholar]
  12. Nelson, P.; Masel, J. Intercellular competition and the inevitability of multicellular aging. Proc. Natl. Acad. Sci. USA 2017. [Google Scholar] [CrossRef]
  13. Fedarko, N.S. Theories and Mechanisms of Aging. In Geriatric Anesthesiology; Springer International Publishing: Cham, Switzerland, 2018; pp. 19–25. [Google Scholar]
  14. Franceschi, C.; Garagnani, P.G.; Morsiani, C.; Conte, M.; Santoro, A.; Grignolio, A.; Monti, D.; Capri, M.; Salvioli, S. The continuum of aging and age-related diseases: Common mechanisms but different rates. Front. Med. 2018, 5, 61. [Google Scholar] [CrossRef]
  15. Gonzaga-Jauregui, C.; Lupski, J.R.; Gibbs, R.A. Human genome sequencing in health and disease. Annu. Rev. Med. 2012, 63, 35–61. [Google Scholar] [CrossRef] [PubMed]
  16. Yashin, A.I.; Arbeev, K.G.; Wu, D.; Arbeeva, L.; Kulminski, A.; Kulminskaya, I.; Akushevich, I.; Ukraintseva, S.V. How Genes Modulate Patterns of Aging-Related Changes on the Way to 100: Biodemographic Models and Methods in Genetic Analyses of Longitudinal Data. N. Am. Actuar J. 2016, 20, 201–232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Rzhetsky, A.; Wajngurt, D.; Park, N.; Zheng, T. Probing genetic overlap among complex human phenotypes. Proc. Natl. Acad. Sci. USA 2007, 104, 11694–11699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Anderson, C.A.; Soranzo, N.; Zeggini, E.; Barrett, J.C. Synthetic associations are unlikely to account for many common disease genome-wide association signals. PLoS Biol. 2011, 9, e1000580. [Google Scholar] [CrossRef] [PubMed]
  19. Yang, J.; Bakshi, A.; Zhu, Z.; Hemani, G.; Vinkhuyzen, A.A.; Lee, S.H.; Robinson, M.R.; Perry, J.R.; Nolte, I.M.; van Vliet-Ostaptchouk, J.V. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat. Genet. 2015, 47, 1114. [Google Scholar] [CrossRef]
  20. Risch, H.A.; McLaughlin, J.R.; Cole, D.E.; Rosen, B.; Bradley, L.; Fan, I.; Tang, J.; Li, S.; Zhang, S.; Shaw, P.A. Population BRCA1 and BRCA2 mutation frequencies and cancer penetrances: A kin–cohort study in Ontario, Canada. J. Natl. Cancer Inst. 2006, 98, 1694–1706. [Google Scholar] [CrossRef]
  21. Haley, B. Hereditary Breast Cancer: The Basics of BRCA and Beyond. Internal Medicine Grand Rounds; UT Southwestern Medical Center; 2016; Available online: https://utswmed-ir.tdl.org/handle/2152.5/3095 (accessed on 7 July 2019).
  22. Wu, X.; Gu, J. Heritability of prostate cancer: A tale of rare variants and common single nucleotide polymorphisms. Ann. Transl. Med. 2016, 4, 206. [Google Scholar] [CrossRef]
  23. Mancuso, N.; Rohland, N.; Rand, K.A.; Tandon, A.; Allen, A.; Quinque, D.; Mallick, S.; Li, H.; Stram, A.; Sheng, X. The contribution of rare variation to prostate cancer heritability. Nat. Genet. 2016, 48, 30. [Google Scholar] [CrossRef]
  24. Walsh, P.C. The search for the missing heritability of prostate cancer. Eur. Urol. 2017, 72, 657–659. [Google Scholar] [CrossRef]
  25. Lecarpentier, J.; Silvestri, V.; Kuchenbaecker, K.B.; Barrowdale, D.; Dennis, J.; McGuffog, L.; Soucy, P.; Leslie, G.; Rizzolo, P.; Navazio, A.S. Prediction of breast and prostate cancer risks in male BRCA1 and BRCA2 mutation carriers using polygenic risk scores. J. Clin. Oncol. 2017, 35, 2240–2250. [Google Scholar] [CrossRef]
  26. De Voer, R.M.; Hahn, M.M.; Weren, R.D.; Mensenkamp, A.R.; Gilissen, C.; van Zelst-Stams, W.A.; Spruijt, L.; Kets, C.M.; Zhang, J.; Venselaar, H. Identification of novel candidate genes for early-onset colorectal cancer susceptibility. PLoS Genet. 2016, 12, e1005880. [Google Scholar] [CrossRef] [PubMed]
  27. Graff, R.E.; Möller, S.; Passarelli, M.N.; Witte, J.S.; Skytthe, A.; Christensen, K.; Tan, Q.; Adami, H.O.; Czene, K.; Harris, J.R. Familial risk and heritability of colorectal cancer in the nordic twin study of cancer. Clin. Gastroenterol. Hepatol. 2017, 15, 1256–1264. [Google Scholar] [CrossRef] [PubMed]
  28. Weissfeld, J.L.; Lin, Y.; Lin, H.M.; Kurland, B.F.; Wilson, D.O.; Fuhrman, C.R.; Pennathur, A.; Romkes, M.; Nukui, T.; Yuan, J.M. Lung cancer risk prediction using common SNPs located in GWAS-identified susceptibility regions. J. Thorac. Oncol. 2015, 10, 1538–1545. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, D.C.; Wang, X. Tomorrow’s genome medicine in lung cancer. In Seminars in Cancer Biology; Elsevier: Shanghai, China, 2017; Volume 42, pp. 39–43. [Google Scholar]
  30. Kanwal, M.; Ding, X.J.; Cao, Y. Familial risk for lung cancer. Oncol. Lett. 2017, 13, 535–542. [Google Scholar] [CrossRef] [PubMed]
  31. Weinberg, R. The Biology of Cancer; Garland Science: New York, NY, USA, 2013. [Google Scholar]
  32. Alzheimer’s Association. 2017 Alzheimer’s disease facts and figures. Alzheimer’s Dement. 2017, 13, 325–373. [Google Scholar] [CrossRef]
  33. Lakatta, E.G.; Levy, D. Arterial and cardiac aging: Major shareholders in cardiovascular disease enterprises: Part I: Aging arteries: A “set up” for vascular disease. Circulation 2003, 107, 139–146. [Google Scholar] [CrossRef] [PubMed]
  34. Lakatta, E.G. So! What’s aging? Is cardiovascular aging a disease? J. Mol. Cell. Cardiol. 2015, 83, 1–13. [Google Scholar] [CrossRef]
  35. Mitchell, T.J.; Turajlic, S.; Rowan, A.; Nicol, D.; Farmery, J.H.; O’Brien, T.; Martincorena, I.; Tarpey, P.; Angelopoulos, N.; Yates, L.R. Timing the landmark events in the evolution of clear cell renal cell cancer: TRACERx renal. Cell 2018, 173, 611–623. [Google Scholar] [CrossRef]
  36. National Academies of Sciences, Engineering, and Medicine. Human Genome Editing: Science, Ethics, and Governance; National Academies Press: Washington, DC, USA, 2017. [Google Scholar]
  37. Nuffield Council on Bioethics. Genome Editing and Human Reproduction: Social and Ethical Issues; Nuffield Council on Bioethics: London, UK, 2018. [Google Scholar]
  38. Liu, X.; Wu, C.; Li, C.; Boerwinkle, E. dbNSFP v3. 0: A One-Stop Database of Functional Predictions and Annotations for Human Nonsynonymous and Splice-Site SNVs. Hum. Mutat. 2016, 37, 235–241. [Google Scholar] [CrossRef]
  39. Beck, T.; Hastings, R.K.; Gollapudi, S.; Free, R.C.; Brookes, A.J. GWAS Central: A comprehensive resource for the comparison and interrogation of genome-wide association studies. Eur. J. Hum. Genet. 2014, 22, 949. [Google Scholar] [CrossRef]
  40. Wright, A.V.; Nuñez, J.K.; Doudna, J.A. Biology and applications of CRISPR systems: Harnessing nature’s toolbox for genome engineering. Cell 2016, 164, 29–44. [Google Scholar] [CrossRef] [PubMed]
  41. Carroll, D. Genome engineering with zinc-finger nucleases. Genetics 2011, 188, 773–782. [Google Scholar] [CrossRef] [PubMed]
  42. Joung, J.K.; Sander, J.D. TALENs: A widely applicable technology for targeted genome editing. Nat. Rev. Mol. Cell Biol. 2013, 14, 49. [Google Scholar] [CrossRef] [PubMed]
  43. Smith, C.J.; Castanon, O.; Said, K.; Volf, V.; Khoshakhlagh, P.; Hornick, A.; Ferreira, R.; Wu, C.T.; Güell, M.; Garg, S.; et al. Enabling large-scale genome editing by reducing DNA nicking. bioRxiv 2019. [Google Scholar] [CrossRef]
  44. Thompson, D.; Aboulhouda, S.; Hysolli, E.; Smith, C.; Wang, S.; Castanon, O.; Church, G. The future of multiplexed eukaryotic genome engineering. ACS Chem. Biol. 2017, 13, 313–325. [Google Scholar] [CrossRef] [PubMed]
  45. Kohman, R.E.; Kunjapur, A.M.; Hysolli, E.; Wang, Y.; Church, G.M. From Designing the Molecules of Life to Designing Life: Future Applications Derived from Advances in DNA Technologies. Angew. Chem. 2018, 57, 4313–4328. [Google Scholar] [CrossRef] [PubMed]
  46. Makeham, W.M. On the law of mortality and construction of annuity tables. J. Inst. Actuar. 1860, 8, 301–310. [Google Scholar] [CrossRef]
  47. Missov, T.I.; Lenart, A. Gompertz–Makeham life expectancies: Expressions and applications. Theor. Popul. Biol. 2013, 90, 29–35. [Google Scholar] [CrossRef]
  48. Stallard, E. Compression of morbidity and mortality: New perspectives. N. Am. Actuar. J. 2016, 20, 341–354. [Google Scholar] [CrossRef]
  49. Geddes, L. Human age limit claim sparks debate. Nat. News 2016. [Google Scholar] [CrossRef]
  50. Barbi, E.; Lagona, F.; Marsili, M.; Vaupel, J.W.; Wachter, K.W. The plateau of human mortality: Demography of longevity pioneers. Science 2018, 360, 1459–1461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Oliynyk, R.T. Age-related late-onset disease heritability patterns and implications for genome-wide association studies. PeerJ 2019, 7, e7168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Stearns, F.W. One hundred years of pleiotropy: A retrospective. Genetics 2010, 186, 767–773. [Google Scholar] [CrossRef] [PubMed]
  53. Paaby, A.B.; Rockman, M.V. The many faces of pleiotropy. Trends Genet. 2013, 29, 66–73. [Google Scholar] [CrossRef] [PubMed]
  54. Wu, Y.T.; Beiser, A.S.; Breteler, M.M.; Fratiglioni, L.; Helmer, C.; Hendrie, H.C.; Honda, H.; Ikram, M.A.; Langa, K.M.; Lobo, A. The changing prevalence and incidence of dementia over time - current evidence. Nat. Rev. Neurol. 2017, 13, 327. [Google Scholar] [CrossRef] [PubMed]
  55. D’Agostino, R.B., Sr.; Vasan, R.S.; Pencina, M.J.; Wolf, P.A.; Massaro, J.M.; Kannel, W.B. General cardiovascular risk profile for use in primary care. Circulation 2008. [Google Scholar] [CrossRef] [PubMed]
  56. Boehme, M.W.; Buechele, G.; Frankenhauser-Mannuss, J.; Mueller, J.; Lump, D.; Boehm, B.O.; Rothenbacher, D. Prevalence, incidence and concomitant co-morbidities of type 2 diabetes mellitus in South Western Germany—A retrospective cohort and case control study in claims data of a large statutory health insurance. BMC Public Health 2015, 15, 855. [Google Scholar] [CrossRef]
  57. Cox, D. Regression Models and Life-Tables. J. R. Stat. Soc. Ser. B (Methodol.) 1972, 34, 187–220. [Google Scholar] [CrossRef]
  58. Chatterjee, N.; Shi, J.; García-Closas, M. Developing and evaluating polygenic risk prediction models for stratified disease prevention. Nat. Rev. Genet. 2016, 17, 392. [Google Scholar] [CrossRef]
  59. Belsky, D.W.; Moffitt, T.E.; Baker, T.B.; Biddle, A.K.; Evans, J.P.; Harrington, H.; Houts, R.; Meier, M.; Sugden, K.; Williams, B.; et al. Polygenic risk and the developmental progression to heavy, persistent smoking and nicotine dependence: Evidence from a 4-decade longitudinal study. JAMA Psychiatry 2013, 70, 534–542. [Google Scholar] [CrossRef]
  60. Hjelmborg, J.; Korhonen, T.; Holst, K.; Skytthe, A.; Pukkala, E.; Kutschke, J.; Harris, J.R.; Mucci, L.A.; Christensen, K.; Czene, K.; et al. Lung cancer, genetic predisposition and smoking: The Nordic Twin Study of Cancer. Thorax 2017, 72, 1021–2027. [Google Scholar] [CrossRef] [PubMed]
  61. Brookmeyer, R.; Gray, S.; Kawas, C. Projections of Alzheimer’s disease in the United States and the public health impact of delaying disease onset. Am. J. Public Health 1998, 88, 1337–1342. [Google Scholar] [PubMed]
  62. Noh, M.; Yip, B.; Lee, Y.; Pawitan, Y. Multicomponent variance estimation for binary traits in family-based studies. Genet. Epidemiol. 2006, 30, 37–47. [Google Scholar] [PubMed]
  63. Pawitan, Y.; Seng, K.C.; Magnusson, P.K. How many genetic variants remain to be discovered? PLoS ONE 2009, 4, e7969. [Google Scholar]
  64. Social Security Administration (US). 2014. Available online: https://www.ssa.gov/oact/STATS/table4c6.html (accessed on 2 June 2019).
  65. Edland, S.D.; Rocca, W.A.; Petersen, R.C.; Cha, R.H.; Kokmen, E. Dementia and Alzheimer disease incidence rates do not vary by sex in Rochester, Minn. Arch. Neurol. 2002, 59, 1589–1593. [Google Scholar] [CrossRef] [PubMed]
  66. Kokmen, E.; Chandra, V.; Schoenberg, B.S. Trends in incidence of dementing illness in Rochester, Minnesota, in three quinquennial periods, 1960–1974. Neurology 1988, 38, 975–980. [Google Scholar] [PubMed]
  67. Hebert, L.E.; Scherr, P.A.; Beckett, L.A.; Albert, M.S.; Pilgrim, D.M.; Chown, M.J.; Funkenstein, H.H.; Evans, D.A. Age-specific incidence of Alzheimer’s disease in a community population. JAMA 1995, 273, 1354–1359. [Google Scholar] [PubMed]
  68. Rothwell, P.; Coull, A.; Silver, L.; Fairhead, J.; Giles, M.; Lovelock, C.; Redgrave, J.; Bull, L.; Welch, S.; Cuthbertson, F. Population-based study of event-rate, incidence, case fatality, and mortality for all acute vascular events in all arterial territories (Oxford Vascular Study). Lancet 2005, 366, 1773–1783. [Google Scholar]
  69. Cancer Research UK. 2018. Available online: http://www.cancerresearchuk.org/health-professional/cancer-statistics-for-the-uk (accessed on 10 November 2018).
  70. Kuchenbaecker, K.B.; Hopper, J.L.; Barnes, D.R.; Phillips, K.A.; Mooij, T.M.; Roos-Blom, M.J.; Jervis, S.; Van Leeuwen, F.E.; Milne, R.L.; Andrieu, N. Risks of breast, ovarian, and contralateral breast cancer for BRCA1 and BRCA2 mutation carriers. JAMA 2017, 317, 2402–2416. [Google Scholar]
Figure 1. Aging coefficients reflect increases in late-onset disease (LOD) liability with age, based on the clinical incidence rate and model genetic architecture polygenic risk score (PRS). (A) Alzheimer’s disease (AD), (B) type 2 diabetes (T2D), (C) cerebral stroke, (D) coronary artery disease (CAD), (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer.
Figure 1. Aging coefficients reflect increases in late-onset disease (LOD) liability with age, based on the clinical incidence rate and model genetic architecture polygenic risk score (PRS). (A) Alzheimer’s disease (AD), (B) type 2 diabetes (T2D), (C) cerebral stroke, (D) coronary artery disease (CAD), (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer.
Ijms 20 03352 g001
Figure 2. Projected LOD cumulative incidence and lifetime risk increase for life expectancy increases of 5, 10, and 15 years (baseline without gene therapy). (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. This is the baseline scenario, without gene therapy or other health improvements for the plotted LOD. It represents the case where life expectancy increases due to causes other than the plotted LOD. Lifetime risk (lifetime cumulative incidence) corresponds to the lifetime (rightmost) values of the plots.
Figure 2. Projected LOD cumulative incidence and lifetime risk increase for life expectancy increases of 5, 10, and 15 years (baseline without gene therapy). (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. This is the baseline scenario, without gene therapy or other health improvements for the plotted LOD. It represents the case where life expectancy increases due to causes other than the plotted LOD. Lifetime risk (lifetime cumulative incidence) corresponds to the lifetime (rightmost) values of the plots.
Ijms 20 03352 g002
Figure 3. LOD incidence rate patterns, baseline and after emulated gene therapy. (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. All individuals in the population had emulated corrective gene therapy editing of, on average, 15 single nucleotide polymorphisms (SNPs) (corresponding to an odds ratio (OR) multiplier of 0.25).
Figure 3. LOD incidence rate patterns, baseline and after emulated gene therapy. (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. All individuals in the population had emulated corrective gene therapy editing of, on average, 15 single nucleotide polymorphisms (SNPs) (corresponding to an odds ratio (OR) multiplier of 0.25).
Ijms 20 03352 g003
Figure 4. Cumulative incidence and lifetime risk following therapy reducing the population risk fourfold. (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. All individuals in the population had emulated corrective gene therapy editing of, on average, 15 SNPs (corresponding to an OR multiplier of 0.25). Result of gene therapy and results in increases in life expectancy of 5, 10, and 15 years. Lifetime risk (lifetime cumulative incidence) corresponds to the lifetime (rightmost) values of the plots.
Figure 4. Cumulative incidence and lifetime risk following therapy reducing the population risk fourfold. (A) Alzheimer’s disease, (B) type 2 diabetes, (C) cerebral stroke, (D) coronary artery disease, (E) breast cancer, (F) prostate cancer, (G) colorectal cancer, and (H) lung cancer. All individuals in the population had emulated corrective gene therapy editing of, on average, 15 SNPs (corresponding to an OR multiplier of 0.25). Result of gene therapy and results in increases in life expectancy of 5, 10, and 15 years. Lifetime risk (lifetime cumulative incidence) corresponds to the lifetime (rightmost) values of the plots.
Ijms 20 03352 g004
Table 1. Lifetime risk corresponding to select discrete LOD PRSs (%).
Table 1. Lifetime risk corresponding to select discrete LOD PRSs (%).
Highly Prevalent LODsCancers
ADT2DStrokeCADBreastProstateColorectalLung
1.043.378.726.821.18.273.943.337.12
0.531.939.913.510.64.141.971.663.56
0.2522.020.06.755.302.070.9840.8321.78
Table 2. Estimate of LOD onset delay of incidence slope at 30% of lifetime risk (values in years).
Table 2. Estimate of LOD onset delay of incidence slope at 30% of lifetime risk (values in years).
Highly Prevalent LODsCancers
ADT2DStrokeCADBreastProstateColorectalLung
16.04.047101021171310
4.01.0312131431161620
1.00.25413111637211916
Table 3. Estimate of LOD onset delay for lifetime risk (values in years).
Table 3. Estimate of LOD onset delay for lifetime risk (values in years).
Highly Prevalent LODsCancers
ADT2DStrokeCADBreastProstateColorectalLung
16.04.01023>40192421
4.01.0291412>402834>40
1.00.253201714>402935>40
Table 4. LOD characteristics and lifetime risk in a range of life expectancy scenarios and with prophylactic gene therapy lowering population PRS fourfold.
Table 4. LOD characteristics and lifetime risk in a range of life expectancy scenarios and with prophylactic gene therapy lowering population PRS fourfold.
Highly Prevalent LODsCancers
ADT2DStrokeCADBreastProstateColorectalLung
Literature and clinical data:
Heritability0.7950.690.550.410.570.400.310.10
Max yearly incidence rate>20%2.5%4.4%3.6%<0.5%<0.8%<0.6%<0.6%
Genetic model SNP count3575212511756251250600400100
Lifetime risk, baseline + longer life:
+5 years life expectancy160%112%128%127%115%123%127%128%
+10 years life expectancy228%123%156%155%130%147%156%156%
+15 years life expectancy293%134%184%184%146%172%187%186%
Lifetime risk, odds ratio (OR) 0.25 therapy versus baseline:
Therapy, unchanged life expectancy70%67%44%50%30%41%31%27%
Therapy, +5 years life expectancy124%77%61%69%36%52%40%34%
Therapy, +10 years life expectancy191%88%81%89%41%65%51%43%
Therapy, +15 years life expectancy260%98%101%112%47%77%62%52%
LODs’ heritability and clinical incidence are discussed in the Methods section. The baseline is considered 100% for lifetime risk comparisons; 160% after a 5-year life expectancy extension, because AD indicates an increase in lifetime risk (LR) by 1.6 times; and 50% after gene therapy which means half the LR, compared to the baseline value. Genetic model single nucleotide polymorphism (SNP) count is the number needed for the common low-effect genetic architecture to achieve each LOD heritability.
Table 5. Genetic architecture scenarios with modeled allele frequencies and effect sizes.
Table 5. Genetic architecture scenarios with modeled allele frequencies and effect sizes.
ScenarioMAF RangeOR RangeMAF ValuesAllele OR Values
A. Common low0.073–0.4991.05–1.150.073, 0.18, 0.286, 0.393, 0.51.05, 1.075, 1.1, 1.125, 1.15
B. Rare medium0.0146–0.09981.28–2.010.0146, 0.036, 0.0572, 0.0785, 0.09981.28, 1.463, 1.645, 1.828, 2.01
To build the genetic architecture, minor allele frequencies (MAFs) and odds ratios (ORs) were chosen using 25 possible combinations of the values in the table. Following Pawitan et al. [63], the variants were assigned to individuals with frequencies proportionate to MAF p k for SNP k, producing, in accordance with the Hardy–Weinberg principle, three genotypes (AA, AB, or BB) for each SNP, with frequencies p k 2 , 2 p k ( 1 p k ) , and ( 1 p k ) 2 , respectively, resulting in the normal distribution of the individual risk ORs.

Share and Cite

MDPI and ACS Style

Oliynyk, R.T. Quantifying the Potential for Future Gene Therapy to Lower Lifetime Risk of Polygenic Late-Onset Diseases. Int. J. Mol. Sci. 2019, 20, 3352. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133352

AMA Style

Oliynyk RT. Quantifying the Potential for Future Gene Therapy to Lower Lifetime Risk of Polygenic Late-Onset Diseases. International Journal of Molecular Sciences. 2019; 20(13):3352. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133352

Chicago/Turabian Style

Oliynyk, Roman Teo. 2019. "Quantifying the Potential for Future Gene Therapy to Lower Lifetime Risk of Polygenic Late-Onset Diseases" International Journal of Molecular Sciences 20, no. 13: 3352. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133352

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop