Causes of Clonal Hematopoiesis: a Review.

PURPOSE OF REVIEW: Clonal hematopoiesis (CH) is an age-dependent process detectable using advanced sequencing technologies and is associated with multiple adverse health outcomes including cardiovascular disease and cancer. The purpose of this review is to summarize known causes of CH mutations and to identify key areas and considerations for future research on CH. RECENT FINDINGS: Studies have identified multiple potential causes of CH mutations including smoking, cancer therapies, cardiometabolic disease, inflammation, and germline risk factors. Additionally, large-scale studies have facilitated the identification of gene-specific effects of CH mutation risk factors that may have unique downstream health implications. For example, cancer therapies and sources of environmental radiation appear to cause CH through their impact on DNA damage repair genes. There is a growing body of evidence defining risk factors for CH mutations. Standardization in the identification of CH mutations may have important implications for future research. Additional studies in underrepresented populations and their diverse environmental exposures are needed to facilitate broad public health impact of the study of CH mutations.

Augmented Lipocalin-2 Is Associated with Chronic Obstructive Pulmonary Disease and Counteracts Lung Adenocarcinoma Development.

Rationale: Early pathogenesis of lung adenocarcinoma (LUAD) remains largely unknown. We found that, relative to wild-type littermates, the innate immunomodulator Lcn2 (lipocalin-2) was increased in normal airways from mice with knockout of the airway lineage gene Gprc5a (Gprc5a-/-) and that are prone to developing inflammation and LUAD. Yet, the role of LCN2 in lung inflammation and LUAD is poorly understood.Objectives: Delineate the role of Lcn2 induction in LUAD pathogenesis.Methods: Normal airway brushings, uninvolved lung tissues, and tumors from Gprc5a-/- mice before and after tobacco carcinogen exposure were analyzed by RNA sequencing. LCN2 mRNA was analyzed in public and in-house data sets of LUAD, lung squamous cancer (LUSC), chronic obstructive pulmonary disease (COPD), and LUAD/LUSC with COPD. LCN2 protein was immunohistochemically analyzed in a tissue microarray of 510 tumors. Temporal lung tumor development, gene expression programs, and host immune responses were compared between Gprc5a-/- and Gprc5a-/-/Lcn2-/- littermates.Measurements and Main Results:Lcn2 was progressively elevated during LUAD development and positively correlated with proinflammatory cytokines and inflammation gene sets. LCN2 was distinctively elevated in human LUADs, but not in LUSCs, relative to normal lungs and was associated with COPD among smokers and patients with LUAD. Relative to Gprc5a-/- mice, Gprc5a-/-/Lcn2-/- littermates exhibited significantly increased lung tumor development concomitant with reduced T-cell abundance (CD4+) and richness, attenuated antitumor immune gene programs, and increased immune cell expression of protumor inflammatory cytokines.Conclusions: Augmented LCN2 expression is a molecular feature of COPD-associated LUAD and counteracts LUAD development in vivo by maintaining antitumor immunity.

A genome-wide approach to children's aggressive behavior: The EAGLE consortium.

Individual differences in aggressive behavior emerge in early childhood and predict persisting behavioral problems and disorders. Studies of antisocial and severe aggression in adulthood indicate substantial underlying biology. However, little attention has been given to genome-wide approaches of aggressive behavior in children. We analyzed data from nine population-based studies and assessed aggressive behavior using well-validated parent-reported questionnaires. This is the largest sample exploring children's aggressive behavior to date (N = 18,988), with measures in two developmental stages (N = 15,668 early childhood and N = 16,311 middle childhood/early adolescence). First, we estimated the additive genetic variance of children's aggressive behavior based on genome-wide SNP information, using genome-wide complex trait analysis (GCTA). Second, genetic associations within each study were assessed using a quasi-Poisson regression approach, capturing the highly right-skewed distribution of aggressive behavior. Third, we performed meta-analyses of genome-wide associations for both the total age-mixed sample and the two developmental stages. Finally, we performed a gene-based test using the summary statistics of the total sample. GCTA quantified variance tagged by common SNPs (10-54%). The meta-analysis of the total sample identified one region in chromosome 2 (2p12) at near genome-wide significance (top SNP rs11126630, P = 5.30 × 10(-8) ). The separate meta-analyses of the two developmental stages revealed suggestive evidence of association at the same locus. The gene-based analysis indicated association of variation within AVPR1A with aggressive behavior. We conclude that common variants at 2p12 show suggestive evidence for association with childhood aggression. Replication of these initial findings is needed, and further studies should clarify its biological meaning. © 2015 Wiley Periodicals, Inc.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

A prospective study of the effects of breastfeeding and FADS2 polymorphisms on cognition and hyperactivity/attention problems.

Breastfeeding has been associated with improved cognitive functioning. There is a beneficial effect on IQ, and possibly on associated phenotypes such as attention problems. It has been suggested that the effect on IQ is moderated by polymorphisms in the FADS2 gene, which is involved in fatty acid metabolism. In this study we tested the relation between breastfeeding and FADS2 polymorphisms on the one hand and IQ, educational attainment, overactivity, and attention problems on the other hand. IQ at age 5, 7, 10, 12, and/or 18 (n = 1,313), educational attainment at age 12 (n = 1,857), overactive behavior at age 3 (n = 2,560), and attention problems assessed at age 7, 10, and 12 years (n = 2,479, n = 2,423, n = 2,226) were predicted by breastfeeding and two SNPs in FADS2 (rs174575 and rs1535). Analyses were performed using structural equation modeling. After correction for maternal education, a main effect of breastfeeding was found for educational attainment at age 12 and overactive behavior at age 3. For IQ, the effect of breastfeeding across age was marginally significant (P = 0.05) and amounted to 1.6 points after correcting for maternal education. Neither a main effect of the FADS2 polymorphisms nor an interaction with breastfeeding was detected for any of the phenotypes. This developmentally informed study confirms that breastfeeding is associated with higher educational attainment at age 12, less overactive behavior at age 3 and a trend toward higher IQ after correction for maternal education. In general, the benefits of breastfeeding were small and did not interact with SNPs in FADS2.

Combination of Sulindac and Bexarotene for Prevention of Intestinal Carcinogenesis in Familial Adenomatous Polyposis.

Familial adenomatous polyposis (FAP) is a hereditary colorectal cancer syndrome, which results in the development of hundreds of adenomatous polyps carpeting the gastrointestinal tract. NSAIDs have reduced polyp burden in patients with FAP and synthetic rexinoids have demonstrated the ability to modulate cytokine-mediated inflammation and WNT signaling. This study examined the use of the combination of an NSAID (sulindac) and a rexinoid (bexarotene) as a durable approach for reducing FAP colonic polyposis to prevent colorectal cancer development. Whole transcriptomic analysis of colorectal polyps and matched normal mucosa in a cohort of patients with FAP to identify potential targets for prevention in FAP was performed. Drug-dose synergism of sulindac and bexarotene in cell lines and patient-derived organoids was assessed, and the drug combination was tested in two different mouse models. This work explored mRNA as a potential predictive serum biomarker for this combination in FAP. Overall, transcriptomic analysis revealed significant activation of inflammatory and cell proliferation pathways. A synergistic effect of sulindac (300 µmol/L) and bexarotene (40 µmol/L) was observed in FAP colonic organoids with primary targeting of polyp tissue compared with normal mucosa. This combination translated into a significant reduction in polyp development in ApcMin/+ and ApcLoxP/+-Cdx2 mice. Finally, the reported data suggest miRNA-21 could serve as a predictive serum biomarker for polyposis burden in patients with FAP. These findings support the clinical development of the combination of sulindac and bexarotene as a treatment modality for patients with FAP. PREVENTION RELEVANCE: This study identified a novel chemopreventive regimen combining sulindac and bexarotene to reduce polyposis in patients with FAP using in silico tools, ex vivo, and in vivo models. This investigation provides the essential groundwork for moving this drug combination forward into a clinical trial.

Defining the Comprehensive Genomic Landscapes of Pancreatic Ductal Adenocarcinoma Using Real-World Endoscopic Aspiration Samples.

PURPOSE: Most patients with pancreatic ductal adenocarcinoma (PDAC) present with surgically unresectable cancer. As a result, endoscopic ultrasound-guided fine-needle aspiration (EUS-FNA) is the most common biospecimen source available for diagnosis in treatment-naïve patients. Unfortunately, these limited samples are often not considered adequate for genomic analysis, precluding the opportunity for enrollment on precision medicine trials. EXPERIMENTAL DESIGN: Applying an epithelial cell adhesion molecule (EpCAM)-enrichment strategy, we show the feasibility of using real-world EUS-FNA for in-depth, molecular-barcoded, whole-exome sequencing (WES) and somatic copy-number alteration (SCNA) analysis in 23 patients with PDAC. RESULTS: Potentially actionable mutations were identified in >20% of patients. Further, an increased mutational burden and higher aneuploidy in WES data were associated with an adverse prognosis. To identify predictive biomarkers for first-line chemotherapy, we developed an SCNA-based complexity score that was associated with response to platinum-based regimens in this cohort. CONCLUSIONS: Collectively, these results emphasize the feasibility of real-world cytology samples for in-depth genomic characterization of PDAC and show the prognostic potential of SCNA for PDAC diagnosis.

The Transcriptomic Landscape of Mismatch Repair-Deficient Intestinal Stem Cells.

Lynch syndrome is the most common cause of hereditary colorectal cancer and is secondary to germline alterations in one of four DNA mismatch repair (MMR) genes. Here we aimed to provide novel insights into the initiation of MMR-deficient (MMRd) colorectal carcinogenesis by characterizing the expression profile of MMRd intestinal stem cells (ISC). A tissue-specific MMRd mouse model (Villin-Cre;Msh2 LoxP/LoxP ) was crossed with a reporter mouse (Lgr5-EGFP-IRES-creERT2) to trace and isolate ISCs (Lgr5+) using flow cytometry. Three different ISC genotypes (Msh2-KO, Msh2-HET, and Msh2-WT) were isolated and processed for mRNA-seq and mass spectrometry, followed by bioinformatic analyses to identify expression signatures of complete MMRd and haplo-insufficiency. These findings were validated using qRT-PCR, IHC, and whole transcriptomic sequencing in mouse tissues, organoids, and a cohort of human samples, including normal colorectal mucosa, premalignant lesions, and early-stage colorectal cancers from patients with Lynch syndrome and patients with familial adenomatous polyposis (FAP) as controls. Msh2-KO ISCs clustered together with differentiated intestinal epithelial cells from all genotypes. Gene-set enrichment analysis indicated inhibition of replication, cell-cycle progression, and the Wnt pathway and activation of epithelial signaling and immune reaction. An expression signature derived from MMRd ISCs successfully distinguished MMRd neoplastic lesions of patients with Lynch syndrome from FAP controls. SPP1 was specifically upregulated in MMRd ISCs and colocalized with LGR5 in Lynch syndrome colorectal premalignant lesions and tumors. These results show that expression signatures of MMRd ISC recapitulate the initial steps of Lynch syndrome carcinogenesis and have the potential to unveil novel biomarkers of early cancer initiation. SIGNIFICANCE: The transcriptomic and proteomic profile of MMR-deficient intestinal stem cells displays a unique set of genes with potential roles as biomarkers of cancer initiation and early progression.

Single-Cell Transcriptomics of Pancreatic Cancer Precursors Demonstrates Epithelial and Microenvironmental Heterogeneity as an Early Event in Neoplastic Progression.

PURPOSE: Early detection of pancreatic ductal adenocarcinoma (PDAC) remains elusive. Precursor lesions of PDAC, specifically intraductal papillary mucinous neoplasms (IPMNs), represent a bona fide pathway to invasive neoplasia, although the molecular correlates of progression remain to be fully elucidated. Single-cell transcriptomics provides a unique avenue for dissecting both the epithelial and microenvironmental heterogeneities that accompany multistep progression from noninvasive IPMNs to PDAC. EXPERIMENTAL DESIGN: Single-cell RNA sequencing was performed through droplet-based sequencing on 5,403 cells from 2 low-grade IPMNs (LGD-IPMNs), 2 high-grade IPMNs (HGD-IPMN), and 2 PDACs (all surgically resected). RESULTS: Analysis of single-cell transcriptomes revealed heterogeneous alterations within the epithelium and the tumor microenvironment during the progression of noninvasive dysplasia to invasive cancer. Although HGD-IPMNs expressed many core signaling pathways described in PDAC, LGD-IPMNs harbored subsets of single cells with a transcriptomic profile that overlapped with invasive cancer. Notably, a proinflammatory immune component was readily seen in low-grade IPMNs, composed of cytotoxic T cells, activated T-helper cells, and dendritic cells, which was progressively depleted during neoplastic progression, accompanied by infiltration of myeloid-derived suppressor cells. Finally, stromal myofibroblast populations were heterogeneous and acquired a previously described tumor-promoting and immune-evading phenotype during invasive carcinogenesis. CONCLUSIONS: This study demonstrates the ability to perform high-resolution profiling of the transcriptomic changes that occur during multistep progression of cystic PDAC precursors to cancer. Notably, single-cell analysis provides an unparalleled insight into both the epithelial and microenvironmental heterogeneities that accompany early cancer pathogenesis and might be a useful substrate to identify targets for cancer interception.See related commentary by Hernandez-Barco et al., p. 2027.

Mutational Heterogeneity in APC and KRAS Arises at the Crypt Level and Leads to Polyclonality in Early Colorectal Tumorigenesis.

Purpose: The majority of genomic alterations causing intratumoral heterogeneity (ITH) in colorectal cancer are thought to arise during early stages of carcinogenesis as a burst but only after truncal mutations in APC have expanded a single founder clone. We have investigated if the initial source of ITH is consequent to multiple independent lineages derived from different crypts harboring distinct truncal APC and driver KRAS mutations, thus challenging the prevailing monoclonal monocryptal model.Experimental Design: High-depth next-generation sequencing and SNP arrays were performed in whole-lesion extracts of 37 familial adenomatous polyposis colorectal adenomas. Also, ultrasensitive genotyping of hotspot mutations of APC and KRAS was performed using nanofluidic PCRs in matched bulk biopsies (n = 59) and crypts (n = 591) from 18 adenomas and seven carcinomas and adjacent normal tissues.Results: Multiple co-occurring truncal APC and driver KRAS alterations were uncovered in whole-lesion extracts from adenomas and subsequently confirmed to belong to multiple clones. Ultrasensitive genotyping of bulk biopsies and crypts revealed novel undetected APC mutations that were prominent among carcinomas, whereas abundant wild-type APC crypts were detected in adenomas. KRAS mutational heterogeneity within crypts was evident in both adenomas and carcinomas with a higher degree of concordance between biopsy and crypt genotyping in carcinomas. Nonrandom heterogeneity among crypts was also observed.Conclusions: The striking degree of nonrandom intercrypt heterogeneity in truncal and driver gene mutations observed in adenomas and carcinomas is consistent with a polycryptal model derived from multiple independent initiation linages as the source of early ITH in colorectal carcinogenesis. Clin Cancer Res; 23(19); 5936-47. ©2017 AACR.

Pathway analysis of bladder cancer genome-wide association study identifies novel pathways involved in bladder cancer development.

Genome-wide association studies (GWAS) are designed to identify individual regions associated with cancer risk, but only explain a small fraction of the inherited variability. Alternative approach analyzing genetic variants within biological pathways has been proposed to discover networks of susceptibility genes with additional effects. The gene set enrichment analysis (GSEA) may complement and expand traditional GWAS analysis to identify novel genes and pathways associated with bladder cancer risk. We selected three GSEA methods: Gen-Gen, Aligator, and the SNP Ratio Test to evaluate cellular signaling pathways involved in bladder cancer susceptibility in a Texas GWAS population. The candidate genetic polymorphisms from the significant pathway selected by GSEA were validated in an independent NCI GWAS. We identified 18 novel pathways (P < 0.05) significantly associated with bladder cancer risk. Five of the most promising pathways (P ≤ 0.001 in any of the three GSEA methods) among the 18 pathways included two cell cycle pathways and neural cell adhesion molecule (NCAM), platelet-derived growth factor (PDGF), and unfolded protein response pathways. We validated the candidate polymorphisms in the NCI GWAS and found variants of RAPGEF1, SKP1, HERPUD1, CACNB2, CACNA1C, CACNA1S, COL4A2, SRC, and CACNA1C were associated with bladder cancer risk. Two CCNE1 variants, rs8102137 and rs997669, from cell cycle pathways showed the strongest associations; the CCNE1 signal at 19q12 has already been reported in previous GWAS. These findings offer additional etiologic insights highlighting the specific genes and pathways associated with bladder cancer development. GSEA may be a complementary tool to GWAS to identify additional loci of cancer susceptibility.

Comparison of marker types and map assumptions using Markov chain Monte Carlo-based linkage analysis of COGA data.

We performed multipoint linkage analysis of the electrophysiological trait ECB21 on chromosome 4 in the full pedigrees provided by the Collaborative Study on the Genetics of Alcoholism (COGA). Three Markov chain Monte Carlo (MCMC)-based approaches were applied to the provided and re-estimated genetic maps and to five different marker panels consisting of microsatellite (STRP) and/or SNP markers at various densities. We found evidence of linkage near the GABRB1 STRP using all methods, maps, and marker panels. Difficulties encountered with SNP panels included convergence problems and demanding computations.

Intelligence: shared genetic basis between Mendelian disorders and a polygenic trait.

Multiple inquiries into the genetic etiology of human traits indicated an overlap between genes underlying monogenic disorders (eg, skeletal growth defects) and those affecting continuous variability of related quantitative traits (eg, height). Extending the idea of a shared genetic basis between a Mendelian disorder and a classic polygenic trait, we performed an association study to examine the effect of 43 genes implicated in autosomal recessive cognitive disorders on intelligence in an unselected Dutch population (N=1316). Using both single-nucleotide polymorphism (SNP)- and gene-based association testing, we detected an association between intelligence and the genes of interest, with genes ELP2, TMEM135, PRMT10, and RGS7 showing the strongest associations. This is a demonstration of the relevance of genes implicated in monogenic disorders of intelligence to normal-range intelligence, and a corroboration of the utility of employing knowledge on monogenic disorders in identifying the genetic variability underlying complex traits.

Attention-deficit/hyperactivity disorder polygenic risk scores predict attention problems in a population-based sample of children.

OBJECTIVE: Clinically, attention-deficit/hyperactivity disorder (ADHD) is characterized by hyperactivity, impulsivity, and inattention and is among the most common childhood disorders. These same traits that define ADHD are variable in the general population, and the clinical diagnosis may represent the extreme end of a continuous distribution of inattentive and hyperactive behaviors. This hypothesis can be tested by assessing the predictive value of polygenic risk scores derived from a discovery sample of ADHD patients in a target sample from the general population with continuous scores of inattention and hyperactivity. In addition, the genetic overlap between ADHD and continuous ADHD scores can be tested across rater and age. METHOD: The Psychiatric Genomics Consortium has performed the largest genome-wide analysis (GWA) study of ADHD so far, including 5,621 clinical patients and 13,589 controls. The effects sizes of single nucleotide polymorphisms (SNPs) estimated in this meta-analysis were used to obtain individual polygenic risk scores in an independent population-based cohort of 2,437 children from the Netherlands Twin Register. The variance explained in Attention Problems (AP) scale scores by the polygenic risk scores was estimated by linear mixed modeling. RESULTS: The ADHD polygenic risk scores significantly predicted both parent and teacher ratings of AP in preschool- and school-aged children. CONCLUSION: These results indicate genetic overlap between a diagnosis of ADHD and AP scale scores across raters and age groups and provides evidence for a dimensional model of ADHD. Future GWA studies on ADHD can likely benefit from the inclusion of population-based cohorts and the analysis of continuous scores.

A genome-wide association meta-analysis of preschool internalizing problems.

OBJECTIVE: Preschool internalizing problems (INT) are highly heritable and moderately genetically stable from childhood into adulthood. Gene-finding studies are scarce. In this study, the influence of genome-wide measured single nucleotide polymorphisms (SNPs) was investigated in 3 cohorts (total N = 4,596 children) in which INT was assessed with the same instrument, the Child Behavior Checklist (CBCL). METHOD: First, genome-wide association (GWA) results were used for density estimation and genome-wide complex trait analysis (GCTA) to calculate the variance explained by all SNPs. Next, a fixed-effect inverse variance meta-analysis of the 3 GWA analyses was carried out. Finally, the overlap in results with prior GWA studies of childhood and adulthood psychiatric disorders and treatment responses was tested by examining whether SNPs associated with these traits jointly showed a significant signal for INT. RESULTS: Genome-wide SNPs explained 13% to 43% of the total variance. This indicates that the genetic architecture of INT mirrors the polygenic model underlying adult psychiatric traits. The meta-analysis did not yield a genome-wide significant signal but was suggestive for the PCSK2 gene located on chromosome 20p12.1. SNPs associated with other psychiatric disorders appeared to be enriched for signals with INT (λ = 1.26, p < .03). CONCLUSION: Our study provides evidence that INT is influenced by many common genetic variants, each with a very small effect, and that, even as early as age 3, genetic variants influencing INT overlap with variants that play a role in childhood and adulthood psychiatric disorders.

Newly identified loci that influence lipid concentrations and risk of coronary artery disease.

To identify genetic variants influencing plasma lipid concentrations, we first used genotype imputation and meta-analysis to combine three genome-wide scans totaling 8,816 individuals and comprising 6,068 individuals specific to our study (1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables) and 2,758 individuals from the Diabetes Genetics Initiative, reported in a companion study in this issue. We subsequently examined promising signals in 11,569 additional individuals. Overall, we identify strongly associated variants in eleven loci previously implicated in lipid metabolism (ABCA1, the APOA5-APOA4-APOC3-APOA1 and APOE-APOC clusters, APOB, CETP, GCKR, LDLR, LPL, LIPC, LIPG and PCSK9) and also in several newly identified loci (near MVK-MMAB and GALNT2, with variants primarily associated with high-density lipoprotein (HDL) cholesterol; near SORT1, with variants primarily associated with low-density lipoprotein (LDL) cholesterol; near TRIB1, MLXIPL and ANGPTL3, with variants primarily associated with triglycerides; and a locus encompassing several genes near NCAN, with variants strongly associated with both triglycerides and LDL cholesterol). Notably, the 11 independent variants associated with increased LDL cholesterol concentrations in our study also showed increased frequency in a sample of coronary artery disease cases versus controls.