Mathematical models for the tumour control probability (TCP) are used to estimate the expected success of radiation treatment protocols of cancer. There are several TCP models in the literature, from the simplest (Poissonian TCP) to the well-advanced stochastic birth–death processes. Simple and complex models often make the same predictions. Hence, here, we present a systematic study where we compare six of these TCP models: the Poisson TCP, the Zaider–Minerbo TCP, a Monte Carlo TCP and their corresponding cell cycle (two-compartment) models. Several clinical non-uniform treatment protocols for prostate cancer are employed to evaluate these models. These include fractionated external beam radiotherapies, and high and low dose rate brachytherapies. We find that in realistic treatment scenarios, all one-compartment models and all two-compartment models give basically the same results. A difference occurs between one-compartment and two-compartment models due to reduced radiosensitivity of quiescent cells.We find that care must be taken for the right choice of parameters, such as the radiosensitivities α and β and the hazard function h. Typically, different hazard functions are used for fractionated treatment (fractionated survival fraction) and for brachytherapies (Lea–Catcheside protraction factor). We were able to combine these two approaches into one 'effective' hazard function. Based on our results, we can recommend the use of the Poissonian TCP for everyday treatment planning. More complicated models should only be used when absolutely necessary.
Over 400 variants in the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) are CF-causing. CFTR modulators target variants to improve lung function, but marked variability in response exists and current therapies do not address all CF-causing variants highlighting unmet needs. Alternative epithelial ion channel/transporters such as SLC26A9 could compensate for CFTR dysfunction, providing therapeutic targets that may benefit all individuals with CF. We investigate the relationship between rs7512462, a marker of SLC26A9 activity, and lung function pre- and post-treatment with CFTR modulators in Canadian and US CF cohorts, in the general population, and in those with chronic obstructive pulmonary disease (COPD). Rs7512462 CC genotype is associated with greater lung function in CF individuals with minimal function variants (for which there are currently no approved therapies; p = 0.008); and for gating (p = 0.033) and p.Phe508del/ p.Phe508del (p = 0.006) genotypes upon treatment with CFTR modulators. In parallel, human nasal epithelia with CC and p.Phe508del/p.Phe508del after Ussing chamber analysis of a combination of approved and experimental modulator treatments show greater CFTR function (p = 0.0022). Beyond CF, rs7512462 is associated with peak expiratory flow in a meta-analysis of the UK Biobank and Spirometa Consortium (p = 2.74 × 10
Abstract Variable number of tandem repeats (VNTRs) are major source of genetic variation in human. However due to their repetitive nature and large size, it is challenging to genotype them by short-read sequencing. Therefore, there is limited understanding of how they contribute to complex traits such as cystic fibrosis (CF) lung function. Genome-wide association study (GWAS) of CF lung disease identified two independent signals near SLC9A3 displaying a high density of VNTRs and CpG islands. Here, we used long-read (PacBio) phased sequence (N=58) to identify the boundaries and lengths of 49 common (frequency >2%) VNTRs in the region. Subsequently, associations of the VNTRs with gene expression were investigated in CF nasal epithelia using RNA sequencing (N=46). Two VNTRs tagged by the two GWAS signals and overlapping CpG islands were independently associated with SLC9A3 expression in CF nasal epithelia. The two VNTRs together explained 24% of SLC9A3 gene expression variation. One of them was also associated with TPPP expression. We then showed that the VNTR lengths can be estimated with good accuracy in short-read sequence in a subset of individuals with data on both long (PacBio) and short-read (10X Genomics) technologies (N=52). VNTR lengths were then estimated in the Genotype-Tissue Expression project (GTEx) and their association with gene expression was investigated. Both VNTRs were associated with SLC9A3 expression in multiple non-CF GTEx tissues including lung. The results confirm that VNTRs can explain substantial variation in gene expression and be responsible for GWAS signals, and highlight the critical role of long-read sequencing.
Cystic Fibrosis (CF) exhibits morbidity in several organs, including progressive lung disease in all patients and intestinal obstruction at birth (meconium ileus) in ~15%. Individuals with the same causal CFTR mutations show variable disease presentation which is partly attributed to modifier genes. With >6,500 participants from the International CF Gene Modifier Consortium, genome-wide association investigation identified a new modifier locus for meconium ileus encompassing ATP12A on chromosome 13 (min p = 3.83x10-10); replicated loci encompassing SLC6A14 on chromosome X and SLC26A9 on chromosome 1, (min p<2.2x10-16, 2.81x10−11, respectively); and replicated a suggestive locus on chromosome 7 near PRSS1 (min p = 2.55x10-7). PRSS1 is exclusively expressed in the exocrine pancreas and was previously associated with non-CF pancreatitis with functional characterization demonstrating impact on PRSS1 gene expression. We thus asked whether the other meconium ileus modifier loci impact gene expression and in which organ. We developed and applied a colocalization framework called the Simple Sum (SS) that integrates regulatory and genetic association information, and also contrasts colocalization evidence across tissues or genes. The associated modifier loci colocalized with expression quantitative trait loci (eQTLs) for ATP12A (p = 3.35x10-8), SLC6A14 (p = 1.12x10-10) and SLC26A9 (p = 4.48x10-5) in the pancreas, even though meconium ileus manifests in the intestine. The meconium ileus susceptibility locus on chromosome X appeared shifted in location from a previously identified locus for CF lung disease severity. Using the SS we integrated the lung disease association locus with eQTLs from nasal epithelia of 63 CF participants and demonstrated evidence of colocalization with airway-specific regulation of SLC6A14 (p = 2.3x10-4). Cystic Fibrosis is realizing the promise of personalized medicine, and identification of the contributing organ and understanding of tissue specificity for a gene modifier is essential for the next phase of personalizing therapeutic strategies.
Background Pseudomonas aeruginosa is a common pathogen that contributes to progressive lung disease in cystic fibrosis (CF). Genetic factors other than CF-causing CFTR (CF transmembrane conductance regulator) variations contribute ∼85% of the variation in chronic P. aeruginosa infection age in CF according to twin studies, but the susceptibility loci remain unknown. Our objective is to advance understanding of the genetic basis of host susceptibility to P. aeruginosa infection. Materials and methods We conducted a genome-wide association study of chronic P. aeruginosa infection age in 1037 Canadians with CF. We subsequently assessed the genetic correlation between chronic P. aeruginosa infection age and lung function through polygenic risk score (PRS) analysis and inferred their causal relationship through bidirectional Mendelian randomisation analysis. Results Two novel genome-wide significant loci with lead single nucleotide polymorphisms (SNPs) rs62369766 (chr5p12; p=1.98×10 −8 ) and rs927553 (chr13q12.12; p=1.91×10 −8 ) were associated with chronic P. aeruginosa infection age. The rs62369766 locus was validated using an independent French cohort (n=501). Furthermore, the PRS constructed from CF lung function-associated SNPs was significantly associated with chronic P. aeruginosa infection age (p=0.002). Finally, our analysis presented evidence for a causal effect of lung function on chronic P. aeruginosa infection age (β=0.782 years, p=4.24×10 −4 ). In the reverse direction, we observed a moderate effect (β=0.002, p=0.012). Conclusions We identified two novel loci that are associated with chronic P. aeruginosa infection age in individuals with CF. Additionally, we provided evidence of common genetic contributors and a potential causal relationship between P. aeruginosa infection susceptibility and lung function in CF. Therapeutics targeting these genetic factors may delay the onset of chronic infections, which account for significant remaining morbidity in CF.
Abstract Motivation Integration of next generation sequencing data (NGS) across different research studies can improve the power of genetic association testing by increasing sample size and can obviate the need for sequencing controls. Unfortunately, if differential genotype uncertainty across studies is not accounted for, combining data sets can also produce spurious association results. The robust variance score statistic (RVS) for genetic association of rare and common variants has been shown to effectively adjust for bias caused by the differences in read depth in case-control genetic association studies when the two groups were sequenced using different experimental designs. To enable consortium research, the aggregation of several data sets for genetic association analysis of quantitative and binary traits with covariate adjustment is required, and we developed the Variant Integration Kit for NGS (VikNGS) that expands the functionality of RVS (vRVS) for this purpose. Results VikNGS is a fast and computationally efficient cross-platform software package that provides an implementation for vRVS, as well as conventional rare and common variant genotype-based association analysis approaches. The package includes a graphical user interface that contains power simulation functionality and data visualization tools. Availability and Implementation The VikNGS package can be downloaded at http://www.tcag.ca/tools/index.html Documentation can be found at https://VikNGSdocs.readthedocs.io/en/latest/ Contact lisa.strug@sickkids.ca Supplementary information Supplementary data are available at Bioinformatics online.
The normal tissue complication probability (NTCP) is a measure for the estimated side effects of a given radiation treatment schedule. Here we use a stochastic logistic birth–death process to define an organ-specific and patient-specific NTCP. We emphasize an asymptotic simplification which relates the NTCP to the solution of a logistic differential equation. This framework is based on simple modelling assumptions and it prepares a framework for the use of the NTCP model in clinical practice. As example, we consider side effects of prostate cancer brachytherapy such as increase in urinal frequency, urinal retention and acute rectal dysfunction.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
