The genomic basis of adaptive evolution in threespine sticklebacks.

Marine stickleback fish have colonized and adapted to thousands of streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high-quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of twenty additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine-freshwater divergence. Our results indicate that reuse of globally shared standing genetic variation, including chromosomal inversions, has an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine-freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.

Myocardial NF-κB activation is essential for zebrafish heart regeneration.

Heart regeneration offers a novel therapeutic strategy for heart failure. Unlike mammals, lower vertebrates such as zebrafish mount a strong regenerative response following cardiac injury. Heart regeneration in zebrafish occurs by cardiomyocyte proliferation and reactivation of a cardiac developmental program, as evidenced by induction of gata4 regulatory sequences in regenerating cardiomyocytes. Although many of the cellular determinants of heart regeneration have been elucidated, how injury triggers a regenerative program through dedifferentiation and epicardial activation is a critical outstanding question. Here, we show that NF-κB signaling is induced in cardiomyocytes following injury. Myocardial inhibition of NF-κB activity blocks heart regeneration with pleiotropic effects, decreasing both cardiomyocyte proliferation and epicardial responses. Activation of gata4 regulatory sequences is also prevented by NF-κB signaling antagonism, suggesting an underlying defect in cardiomyocyte dedifferentiation. Our results implicate NF-κB signaling as a key node between cardiac injury and tissue regeneration.

A signature of pre-operative biomarkers of cellular senescence to predict risk of cardiac and kidney adverse events after cardiac surgery.

Objective: Understand the potential for pre-operative biomarkers of cellular senescence, a primary aging mechanism, to predict risk of cardiac surgery-associated adverse events. Methods: Biomarkers of senescence were assessed in blood samples collected prior to surgery in 331 patients undergoing CABG +/-valve repair or replacement. Patients were followed throughout the hospital stay and at a 30-day follow-up visit. Logistic regression models for pre-operative risk prediction were built for age-related clinical outcomes with high incidence including KDIGO-defined acute kidney injury (AKI), decline in eGFR ≥25% between pre-op and 30 days, and MACKE30, a composite endpoint of major adverse cardiac and kidney events at 30d. Results: AKI occurred in 19.9% of patients, persistent decline in kidney function at 30d occurred in 11.0%, and MACKE30 occurred in 13.4%. A network of six biomarkers of senescence (p16, p14, LAG3, CD244, CD28 and suPAR) were able to identify patients at risk for AKI (AUC 0.76), kidney decline at 30d (AUC 0.73), and MACKE30 (AUC 0.71). Comparing the top and bottom tertiles of senescence-based risk models, patients in the top tertile had 7.8 (3.3-8.4) higher odds of developing AKI, 4.5 (1.6-12.6) higher odds of developing renal decline at 30d, and 5.7 (2.1-15.6) higher odds of developing MACKE30. All models remained significant when adjusted for clinical variables. Patients with kidney function decline at 30d were largely non-overlapping and clinically distinct from those who experienced AKI, suggesting a different etiology. Typical clinical factors that predispose to AKI (e.g., age, CKD, surgery type) associated with AKI but not the 30d decline endpoint which was instead associated with new-onset atrial fibrillation. Conclusions: A six-member network of biomarkers of senescence, a fundamental mechanism of aging, can identify patients for risk of adverse kidney and cardiac events when measured pre-operatively.

Profiling an integrated network of cellular senescence and immune resilience measures in natural aging: a prospective multi-cohort study.

Background: Biological aging begins decades before the onset of age-related clinical conditions and is mediated by both cellular senescence and declining adaptive immune function. These processes are functionally related with the rate of senescent cell accumulation dependent upon a balance between induction and immune clearance. We previously showed that biomarkers in these domains can identify patients at-risk of surgery-related adverse events. Here, we describe evidence of clinical relevance in early aging and metabolic phenotypes in a general adult population. Methods: We enrolled a total of 482 participants (ages 25-90) into two prospective, cross-sectional healthy aging cohorts. Expression of biomarkers of adaptive immune function and cellular senescence (SapereX) was measured in CD3+ T cells isolated from peripheral blood. Findings: We established a network of biomarkers of adaptive immune function that correlate with cellular senescence and associate with early aging phenotypes. SapereX immune components associated with a decrease in CD4+ T cells, an increase in cytotoxic CD8+ T cells, and a loss of CD8+ naïve T cells (Pearson correlation 0.3-0.6). These components also associated with a metric of immune resilience, an ability to withstand antigen challenge and inflammation. In contrast, SapereX components were only weakly associated with GlycanAge (Pearson correlation 0.03-0.15) and commonly used DNA methylation clocks (Pearson correlation 0-0.25). Finally, SapereX biomarkers, in particular p16, were associated with chronic inflammation and metabolic dysregulation. Interpretation: Measurement of SapereX biomarkers may capture essential elements of the relationship between cellular senescence and dysregulated adaptive immune function and may provide a benchmark for clinically relevant health decisions.

A biomarker of aging, p16, predicts peripheral neuropathy in women receiving adjuvant taxanes for breast cancer.

Identifying patients at higher risk of chemotherapy-induced peripheral neuropathy (CIPN) is a major unmet need given its high incidence, persistence, and detrimental effect on quality of life. We determined if the expression of p16, a biomarker of aging and cellular senescence, predicts CIPN in a prospective, multi-center study of 152 participants enrolled between 2014 and 2018. Any women with newly diagnosed Stage I-III breast cancer scheduled to receive taxane-containing chemotherapy was eligible. The primary outcome was development of grade 2 or higher CIPN during chemotherapy graded by the clinician before each chemotherapy cycle (NCI-CTCAE v5 criteria). We measured p16 expression in peripheral blood T cells by qPCR before and at the end of chemotherapy. A multivariate model identified risk factors for CIPN and included taxane regimen type, p16Age Gap, a measure of discordance between chronological age and p16 expression, and p16 expression before chemotherapy. Participants with higher p16Age Gap-higher chronological age but lower p16 expression prior to chemotherapy - were at the highest risk. In addition, higher levels of p16 before treatment, regardless of patient age, conferred an increased risk of CIPN. Incidence of CIPN positively correlated with chemotherapy-induced increase in p16 expression, with the largest increase seen in participants with the lowest p16 expression before treatment. We have shown that p16 expression levels before treatment can identify patients at high risk for taxane-induced CIPN. If confirmed, p16 might help guide chemotherapy selection in early breast cancer.

The genetics of adaptive shape shift in stickleback: pleiotropy and effect size.

The distribution of effect sizes of genes underlying adaptation is unknown (Orr 2005). Are suites of traits that diverged under natural selection controlled by a few pleiotropic genes of large effect (major genes model), by many independently acting genes of small effect (infinitesimal model), or by a combination, with frequency inversely related to effect size (geometric model)? To address this we carried out a quantitative trait loci (QTL) study of a suite of 54 position traits describing body shapes of two threespine stickleback species: an ancestral Pacific marine form and a highly derived benthic species inhabiting a geologically young lake. About half of the 26 detected QTL affected just one coordinate and had small net effects, but several genomic regions affected multiple aspects of shape and had large net effects. The distribution of effect sizes followed the gamma distribution, as predicted by the geometric model of adaptation when detection limits are taken into account. The sex-determining chromosome region had the largest effect of any QTL. Ancestral sexual dimorphism was similar to the direction of divergence, and was largely eliminated during freshwater adaptation, suggesting that sex differences may provide variation upon which selection can act. Several shape QTL are linked to Eda, a major gene responsible for reduction of lateral body armor in freshwater. Our results are consistent with predictions of the geometric model of adaptation. Shape evolution in stickleback results from a few genes with large and possibly widespread effects and multiple genes of smaller effect.

Expression of p16INK4a is a biomarker of chondrocyte aging but does not cause osteoarthritis.

Cellular senescence drives a functional decline of numerous tissues with aging by limiting regenerative proliferation and/or by producing pro-inflammatory molecules known as the senescence-associated secretory phenotype (SASP). The senescence biomarker p16INK4a is a potent inhibitor of the cell cycle but is not essential for SASP production. Thus, it is unclear whether p16INK4a identifies senescence in hyporeplicative cells such as articular chondrocytes and whether p16INK4a contributes to pathologic characteristics of cartilage aging. To address these questions, we examined the role of p16INK4a in murine and human models of chondrocyte aging. We observed that p16INK4a mRNA expression was significantly upregulated with chronological aging in murine cartilage (~50-fold from 4 to 18 months of age) and in primary human chondrocytes from 57 cadaveric donors (r2  = .27, p < .0001). Human chondrocytes exhibited substantial replicative potential in vitro that depended on the activity of cyclin-dependent kinases 4 or 6 (CDK4/6), and proliferation was reduced in cells from older donors with increased p16INK4a expression. Moreover, increased chondrocyte p16INK4a expression correlated with several SASP transcripts. Despite the relationship between p16INK4a expression and these features of senescence, somatic inactivation of p16INK4a in chondrocytes of adult mice did not mitigate SASP expression and did not alter the rate of osteoarthritis (OA) with physiological aging or after destabilization of the medial meniscus. These results establish that p16INK4a expression is a biomarker of dysfunctional chondrocytes, but that the effects of chondrocyte senescence on OA are more likely driven by production of SASP molecules than by loss of chondrocyte replicative function.

Constraints on utilization of the EDA-signaling pathway in threespine stickleback evolution.

Many traits evolve in parallel in widely separated populations. The evolutionary radiation of threespine sticklebacks provides a powerful model for testing the molecular basis of parallel evolution in vertebrates. Although marine sticklebacks are completely covered with bony armor plates, most freshwater populations have dramatic reductions in plates. Recent genetic studies have shown that major changes in armor patterning are likely due to regulatory alterations in the gene encoding the secreted signaling molecule ectodysplasin (EDA). In mammals, mutations in many different components of the EDA-signaling pathway produce similar changes in hair, teeth, sweat glands, and dermal bones. To test whether other genes in the EDA pathway also control natural variation in armor plates, we identified and mapped stickleback EDA Receptor (EDAR), the EDAR-Associated Death Domain adaptor, Tumor Necrosis Factor Receptor (TNFR) SuperFamily member 19, its adaptor TNFR-Associated Factor 6, and the downstream regulator nuclear factor kappa B Essential Modulator (NEMO). In contrast to the diversity of genes underlying ectodermal dysplasia disease phenotypes in humans, none of these EDA pathway components map to chromosomes previously shown to modify armor plates in natural populations, though EDAR showed a small but significant effect on plate number. We further investigated whether these genes exhibit differences in copy number, target size, or genomic organization that might make them less suitable targets for evolutionary change. In comparison with EDA, all these genes have smaller surrounding noncoding (putative regulatory) regions, with fewer evolutionarily conserved regions. We suggest that the presence of highly modular cis-acting control sequences may be a key factor influencing the likelihood that particular genes will serve as the basis of major phenotypic changes in nature.

Induction of the neural crest: a multigene process.

In the embryo, the neural crest is an important population of cells that gives rise to diverse derivatives, including the peripheral nervous system and the craniofacial skeleton. Evolutionarily, the neural crest is of interest as an important innovation in vertebrates. Experimentally, it represents an excellent system for studying fundamental developmental processes, such as tissue induction. Classical embryologists have identified interactions between tissues that lead to neural crest formation. More recently, geneticists and molecular biologists have identified the genes that are involved in these interactions; this recent work has revealed that induction of the neural crest is a complex multistep process that involves many genes.

Identification and characterization of a calcium channel gamma subunit expressed in differentiating neurons and myoblasts.

Transient elevations of intracellular calcium (calcium transients) play critical roles in many developmental processes, including differentiation. Although the factors that regulate calcium transients are not clearly defined, calcium influx may be controlled by molecules interacting with calcium channels, including channel regulatory subunits. Here, we describe the chick gamma4 regulatory subunit (CACNG4), the first such subunit to be characterized in early development. CACNG4 is expressed early in the cranial neural plate, and later in the cranial and dorsal root ganglia; importantly, the timing of this later expression correlates precisely with the onset of neuronal differentiation. CACNG4 expression is also observed in nonneuronal tissues undergoing differentiation, specifically the myotome and a subpopulation of differentiating myoblasts in the limb bud. Finally, within the distal cranial ganglia, we show that CACNG4 is expressed in placode-derived cells (prospective neurons), but also, surprisingly, in neural crest-derived cells, previously shown to form only glia in this location; contrary to these previous results, we find that neural crest cells can form neurons in the distal ganglia. Given the proposed role of CACNG4 in modulating calcium channels and its expression in differentiating cells, we suggest that CACNG4 may promote differentiation via regulation of intracellular calcium levels.