Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations.

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Brain volume, energy balance, and cardiovascular health in two nonindustrial South American populations.

Little is known about brain aging or dementia in nonindustrialized environments that are similar to how humans lived throughout evolutionary history. This paper examines brain volume (BV) in middle and old age among two indigenous South American populations, the Tsimane and Moseten, whose lifestyles and environments diverge from those in high-income nations. With a sample of 1,165 individuals aged 40 to 94, we analyze population differences in cross-sectional rates of decline in BV with age. We also assess the relationships of BV with energy biomarkers and arterial disease and compare them against findings in industrialized contexts. The analyses test three hypotheses derived from an evolutionary model of brain health, which we call the embarrassment of riches (EOR). The model hypothesizes that food energy was positively associated with late life BV in the physically active, food-limited past, but excess body mass and adiposity are now associated with reduced BV in industrialized societies in middle and older ages. We find that the relationship of BV with both non-HDL cholesterol and body mass index is curvilinear, positive from the lowest values to 1.4 to 1.6 SDs above the mean, and negative from that value to the highest values. The more acculturated Moseten exhibit a steeper decrease in BV with age than Tsimane, but still shallower than US and European populations. Lastly, aortic arteriosclerosis is associated with lower BV. Complemented by findings from the United States and Europe, our results are consistent with the EOR model, with implications for interventions to improve brain health.

Characterization of Morreton virus as an oncolytic virotherapy platform for liver cancers.

BACKGROUND: Morreton virus (MORV) is an oncolytic Vesiculovirus , genetically distinct from vesicular stomatitis virus (VSV). AIM: To report that MORV induced potent cytopathic effects (CPEs) in cholangiocarcinoma (CCA) and hepatocellular carcinoma (HCC) in vitro models. APPROACH AND RESULTS: In preliminary safety analyses, high intranasal doses (up to 10 10 50% tissue culture infectious dose [TCID 50 ]) of MORV were not associated with significant adverse effects in immune competent, non-tumor-bearing mice. MORV was shown to be efficacious in a Hep3B hepatocellular cancer xenograft model but not in a CCA xenograft HuCCT1 model. In an immune competent, syngeneic murine CCA model, single intratumoral treatments with MORV (1 × 10 7 TCID 50 ) triggered a robust antitumor immune response leading to substantial tumor regression and disease control at a dose 10-fold lower than VSV (1 × 10 8 TCID 50 ). MORV led to increased CD8 + cytotoxic T cells without compensatory increases in tumor-associated macrophages and granulocytic or monocytic myeloid-derived suppressor cells. CONCLUSIONS: Our findings indicate that wild-type MORV is safe and can induce potent tumor regression via immune-mediated and immune-independent mechanisms in HCC and CCA animal models without dose limiting adverse events. These data warrant further development and clinical translation of MORV as an oncolytic virotherapy platform.

NeoScore Integrates Characteristics of the Neoantigen:MHC Class I Interaction and Expression to Accurately Prioritize Immunogenic Neoantigens.

Accurate prioritization of immunogenic neoantigens is key to developing personalized cancer vaccines and distinguishing those patients likely to respond to immune checkpoint inhibition. However, there is no consensus regarding which characteristics best predict neoantigen immunogenicity, and no model to date has both high sensitivity and specificity and a significant association with survival in response to immunotherapy. We address these challenges in the prioritization of immunogenic neoantigens by (1) identifying which neoantigen characteristics best predict immunogenicity; (2) integrating these characteristics into an immunogenicity score, the NeoScore; and (3) demonstrating a significant association of the NeoScore with survival in response to immune checkpoint inhibition. One thousand random and evenly split combinations of immunogenic and nonimmunogenic neoantigens from a validated dataset were analyzed using a regularized regression model for characteristic selection. The selected characteristics, the dissociation constant and binding stability of the neoantigen:MHC class I complex and expression of the mutated gene in the tumor, were integrated into the NeoScore. A web application is provided for calculation of the NeoScore. The NeoScore results in improved, or equivalent, performance in four test datasets as measured by sensitivity, specificity, and area under the receiver operator characteristics curve compared with previous models. Among cutaneous melanoma patients treated with immune checkpoint inhibition, a high maximum NeoScore was associated with improved survival. Overall, the NeoScore has the potential to improve neoantigen prioritization for the development of personalized vaccines and contribute to the determination of which patients are likely to respond to immunotherapy.

Synergistic combination of cytotoxic chemotherapy and cyclin-dependent kinase 4/6 inhibitors in biliary tract cancers.

BACKGROUND AND AIMS: Biliary tract cancers (BTCs) are uncommon, but highly lethal, gastrointestinal malignancies. Gemcitabine/cisplatin is a standard-of-care systemic therapy, but has a modest impact on survival and harbors toxicities, including myelosuppression, nephropathy, neuropathy, and ototoxicity. Whereas BTCs are characterized by aberrations activating the cyclinD1/cyclin-dependent kinase (CDK)4/6/CDK inhibitor 2a/retinoblastoma pathway, clinical use of CDK4/6 inhibitors as monotherapy is limited by lack of validated biomarkers, diffident preclinical efficacy, and development of acquired drug resistance. Emerging studies have explored therapeutic strategies to enhance the antitumor efficacy of CDK4/6 inhibitors by the combination with chemotherapy regimens, but their mechanism of action remains elusive. APPROACH AND RESULTS: Here, we report in vitro and in vivo synergy in BTC models, showing enhanced efficacy, reduced toxicity, and better survival with a combination comprising gemcitabine/cisplatin and CDK4/6 inhibitors. Furthermore, we demonstrated that abemaciclib monotherapy had only modest efficacy attributable to autophagy-induced resistance. Notably, triplet therapy was able to potentiate efficacy through elimination of the autophagic flux. Correspondingly, abemaciclib potentiated ribonucleotide reductase catalytic subunit M1 reduction, resulting in sensitization to gemcitabine. CONCLUSIONS: As such, these data provide robust preclinical mechanistic evidence of synergy between gemcitabine/cisplatin and CDK4/6 inhibitors and delineate a path forward for translation of these findings to preliminary clinical studies in advanced BTC patients.

The Pregnancy Pickle: Evolved Immune Compensation Due to Pregnancy Underlies Sex Differences in Human Diseases.

We hypothesize that, ancestrally, sex-specific immune modulation evolved to facilitate survival of the pregnant person in the presence of an invasive placenta and an immunologically challenging pregnancy - an idea we term the 'pregnancy compensation hypothesis' (PCH). Further, we propose that sex differences in immune function are mediated, at least in part, by the evolution of gene content and dosage on the sex chromosomes, and are regulated by reproductive hormones. Finally, we propose that changes in reproductive ecology in industrialized environments exacerbate these evolved sex differences, resulting in the increasing risk of autoimmune disease observed in females, and a counteracting reduction in diseases such as cancer that can be combated by heightened immune surveillance. The PCH generates a series of expectations that can be tested empirically and that may help to identify the mechanisms underlying sex differences in modern human diseases.

FGFR2-IIIb Expression by Immunohistochemistry Has High Specificity in Cholangiocarcinoma with FGFR2 Genomic Alterations.

BACKGROUND: FGFR2 genomic alterations are observed in 10-20% of cholangiocarcinoma (CCA). Although FGFR2 fusions are an important actionable target, FGFR2 protein expression has not been thoroughly characterized. AIMS: To evaluate FGFR2 protein expression in cholangiocarcinoma harboring FGFR2 genomic alterations. METHODS: FGFR2 protein expression was evaluated in 99 CCA cases with two different antibodies. FGFR2 genomic alterations were confirmed via next-generating sequencing (NGS) or FISH. Primary objective was to determine the specificity and sensitivity of FGFR2 immunohistochemistry staining for detecting FGFR2 genomic alterations. Secondary objectives included overall FGFR2 immunohistochemistry staining in CCA patients, and evaluation of whether FGFR2 expression correlates with clinical outcomes including overall survival (OS), progression-free survival (PFS), and time-to-tumor recurrence (TTR). RESULTS: Immunohistochemistry staining with two antibodies against FGFR2, FPR2-D, and clone 98706 showed high accuracy (78.7% and 91.9%) and specificity (82.9% and 97.7%), and moderate sensitivity (53.9% and 57.1%), respectively, when compared with the standard methods for detecting FGFR2 genomic alterations. In a median follow-up of 72 months, there were no statistically significant differences in OS, PFS, and TTR, for patients with positive or negative FGFR2 staining. CONCLUSION: FGFR2 protein expression by immunohistochemistry has high specificity and therefore could be used to imply the presence of FGFR2 genomic alterations in the context of a positive test. In the case of a negative test, NGS or FISH would be necessary to ascertain cases with FGFR2 genomic alterations.

High GILT Expression and an Active and Intact MHC Class II Antigen Presentation Pathway Are Associated with Improved Survival in Melanoma.

The MHC class I Ag presentation pathway in melanoma cells has a well-established role in immune-mediated destruction of tumors. However, the clinical significance of the MHC class II Ag presentation pathway in melanoma cells is less clear. In Ag-presenting cells, IFN-γ-inducible lysosomal thiol reductase (GILT) is critical for MHC class II-restricted presentation of multiple melanoma Ags. Although not expressed in benign melanocytes of nevi, GILT and MHC class II expression is induced in malignant melanocytes in a portion of melanoma specimens. Analysis of The Cancer Genome Atlas cutaneous melanoma data set showed that high GILT mRNA expression was associated with improved overall survival. Expression of IFN-γ, TNF-α, and IL-1ß was positively associated with GILT expression in melanoma specimens. These cytokines were capable of inducing GILT expression in human melanoma cells in vitro. GILT protein expression in melanocytes was induced in halo nevi, which are nevi undergoing immune-mediated regression, and is consistent with the association of GILT expression with improved survival in melanoma. To explore potential mechanisms of GILT's association with patient outcome, we investigated pathways related to GILT function and expression. In contrast to healthy skin specimens, in which the MHC class II pathway was nearly uniformly expressed and intact, there was substantial variation in the MHC class II pathway in the The Cancer Genome Atlas melanoma specimens. Both an active and intact MHC class II pathway were associated with improved overall survival in melanoma. These studies support a role for GILT and the MHC class II Ag presentation pathway in melanoma outcome.

Distinct molecular etiologies of male and female hepatocellular carcinoma.

BACKGROUND: Sex-differences in cancer occurrence and mortality are evident across tumor types; men exhibit higher rates of incidence and often poorer responses to treatment. Targeted approaches to the treatment of tumors that account for these sex-differences require the characterization and understanding of the fundamental biological mechanisms that differentiate them. Hepatocellular Carcinoma (HCC) is the second leading cause of cancer death worldwide, with the incidence rapidly rising. HCC exhibits a male-bias in occurrence and mortality, but previous studies have failed to explore the sex-specific dysregulation of gene expression in HCC. METHODS: Here, we characterize the sex-shared and sex-specific regulatory changes in HCC tumors in the TCGA LIHC cohort using combined and sex-stratified differential expression and eQTL analyses. RESULTS: By using a sex-specific differential expression analysis of tumor and tumor-adjacent samples, we uncovered etiologically relevant genes and pathways differentiating male and female HCC. While both sexes exhibited activation of pathways related to apoptosis and cell cycle, males and females differed in the activation of several signaling pathways, with females showing PPAR pathway enrichment while males showed PI3K, PI3K/AKT, FGFR, EGFR, NGF, GF1R, Rap1, DAP12, and IL-2 signaling pathway enrichment. Using eQTL analyses, we discovered germline variants with differential effects on tumor gene expression between the sexes. 24.3% of the discovered eQTLs exhibit differential effects between the sexes, illustrating the substantial role of sex in modifying the effects of eQTLs in HCC. The genes that showed sex-specific dysregulation in tumors and those that harbored a sex-specific eQTL converge in clinically relevant pathways, suggesting that the molecular etiologies of male and female HCC are partially driven by differential genetic effects on gene expression. CONCLUSIONS: Sex-stratified analyses detect sex-specific molecular etiologies of HCC. Overall, our results provide new insight into the role of inherited genetic regulation of transcription in modulating sex-differences in HCC etiology and provide a framework for future studies on sex-biased cancers.

Quantification of read species behavior within whole genome sequencing of cancer genomes for the stratification and visualization of genomic variation.

The cancer genome is abnormal genome, and the ability to monitor its sequence had undergone a technological revolution. Yet prognosis and diagnosis remain an expert-based decision, with only limited abilities to provide machine-based decisions. We introduce a heterogeneity-based method for stratifying and visualizing whole-genome sequencing (WGS) reads. This method uses the heterogeneity within WGS reads to markedly reduce the dimensionality of next-generation sequencing data; it is available through the tool HiBS (Heterogeneity-Based Subclassification) that allows cancer sample classification. We validated HiBS using >200 WGS samples from nine different cancer types from The Cancer Genome Atlas (TCGA). With HiBS, we show progress with two WGS related issues: (i) differentiation between normal (NB) and tumor (TP) samples based solely on the information structure of their WGS data, and (ii) identification of specific regions of chromosomal amplification/deletion and their association with tumor stage. By comparing results to those obtained through available WGS analyses tools, we demonstrate some of the novelties obtained by the approach implemented in HiBS and also show nearly perfect normal/tumor classification, used to identify known and unknown chromosomal aberrations. Finally, the HiBS index has been associated with breast cancer tumor stage.

CREBBP mutations in relapsed acute lymphoblastic leukaemia.

Relapsed acute lymphoblastic leukaemia (ALL) is a leading cause of death due to disease in young people, but the biological determinants of treatment failure remain poorly understood. Recent genome-wide profiling of structural DNA alterations in ALL have identified multiple submicroscopic somatic mutations targeting key cellular pathways, and have demonstrated substantial evolution in genetic alterations from diagnosis to relapse. However, DNA sequence mutations in ALL have not been analysed in detail. To identify novel mutations in relapsed ALL, we resequenced 300 genes in matched diagnosis and relapse samples from 23 patients with ALL. This identified 52 somatic non-synonymous mutations in 32 genes, many of which were novel, including the transcriptional coactivators CREBBP and NCOR1, the transcription factors ERG, SPI1, TCF4 and TCF7L2, components of the Ras signalling pathway, histone genes, genes involved in histone modification (CREBBP and CTCF), and genes previously shown to be targets of recurring DNA copy number alteration in ALL. Analysis of an extended cohort of 71 diagnosis-relapse cases and 270 acute leukaemia cases that did not relapse found that 18.3% of relapse cases had sequence or deletion mutations of CREBBP, which encodes the transcriptional coactivator and histone acetyltransferase CREB-binding protein (CREBBP, also known as CBP). The mutations were either present at diagnosis or acquired at relapse, and resulted in truncated alleles or deleterious substitutions in conserved residues of the histone acetyltransferase domain. Functionally, the mutations impaired histone acetylation and transcriptional regulation of CREBBP targets, including glucocorticoid responsive genes. Several mutations acquired at relapse were detected in subclones at diagnosis, suggesting that the mutations may confer resistance to therapy. These results extend the landscape of genetic alterations in leukaemia, and identify mutations targeting transcriptional and epigenetic regulation as a mechanism of resistance in ALL.

Tyrosine kinome sequencing of pediatric acute lymphoblastic leukemia: a report from the Children's Oncology Group TARGET Project.

One recently identified subtype of pediatric B-precursor acute lymphoblastic leukemia (ALL) has been termed BCR-ABL1-like or Ph-like because of similarity of the gene expression profile to BCR-ABL1 positive ALL suggesting the presence of lesions activating tyrosine kinases, frequent alteration of IKZF1, and poor outcome. Prior studies demonstrated that approximately half of these patients had genomic lesions leading to CRLF2 overexpression, with half of such cases harboring somatic mutations in the Janus kinases JAK1 and JAK2. To determine whether mutations in other tyrosine kinases might also occur in ALL, we sequenced the tyrosine kinome and downstream signaling genes in 45 high-risk pediatric ALL cases with either a Ph-like gene expression profile or other alterations suggestive of activated kinase signaling. Aside from JAK mutations and 1 FLT3 mutation, no somatic mutations were found in any other tyrosine kinases, suggesting that alternative mechanisms are responsible for activated kinase signaling in high-risk ALL.

Integrated cross-species transcriptional network analysis of metastatic susceptibility.

Metastatic disease is the proximal cause of mortality for most cancers and remains a significant problem for the clinical management of neoplastic disease. Recent advances in global transcriptional analysis have enabled better prediction of individuals likely to progress to metastatic disease. However, minimal overlap between predictive signatures has precluded easy identification of key biological processes contributing to the prometastatic transcriptional state. To overcome this limitation, we have applied network analysis to two independent human breast cancer datasets and three different mouse populations developed for quantitative analysis of metastasis. Analysis of these datasets revealed that the gene membership of the networks is highly conserved within and between species, and that these networks predicted distant metastasis free survival. Furthermore these results suggest that susceptibility to metastatic disease is cell-autonomous in estrogen receptor-positive tumors and associated with the mitotic spindle checkpoint. In contrast, nontumor genetics and pathway activities-associated stromal biology are significant modifiers of the rate of metastatic spread of estrogen receptor-negative tumors. These results suggest that the application of network analysis across species may provide a robust method to identify key biological programs associated with human cancer progression.

Addition of MR imaging features and genetic biomarkers strengthens glioblastoma survival prediction in TCGA patients.

PURPOSE: The purpose of our study was to assess whether a model combining clinical factors, MR imaging features, and genomics would better predict overall survival of patients with glioblastoma (GBM) than either individual data type. METHODS: The study was conducted leveraging The Cancer Genome Atlas (TCGA) effort supported by the National Institutes of Health. Six neuroradiologists reviewed MRI images from The Cancer Imaging Archive (http://cancerimagingarchive.net) of 102 GBM patients using the VASARI scoring system. The patients' clinical and genetic data were obtained from the TCGA website (http://www.cancergenome.nih.gov/). Patient outcome was measured in terms of overall survival time. The association between different categories of biomarkers and survival was evaluated using Cox analysis. RESULTS: The features that were significantly associated with survival were: (1) clinical factors: chemotherapy; (2) imaging: proportion of tumor contrast enhancement on MRI; and (3) genomics: HRAS copy number variation. The combination of these three biomarkers resulted in an incremental increase in the strength of prediction of survival, with the model that included clinical, imaging, and genetic variables having the highest predictive accuracy (area under the curve 0.679±0.068, Akaike's information criterion 566.7, P<0.001). CONCLUSION: A combination of clinical factors, imaging features, and HRAS copy number variation best predicts survival of patients with GBM.

Pathways of distinction analysis: a new technique for multi-SNP analysis of GWAS data.

Genome-wide association studies (GWAS) have become increasingly common due to advances in technology and have permitted the identification of differences in single nucleotide polymorphism (SNP) alleles that are associated with diseases. However, while typical GWAS analysis techniques treat markers individually, complex diseases (cancers, diabetes, and Alzheimers, amongst others) are unlikely to have a single causative gene. Thus, there is a pressing need for multi-SNP analysis methods that can reveal system-level differences in cases and controls. Here, we present a novel multi-SNP GWAS analysis method called Pathways of Distinction Analysis (PoDA). The method uses GWAS data and known pathway-gene and gene-SNP associations to identify pathways that permit, ideally, the distinction of cases from controls. The technique is based upon the hypothesis that, if a pathway is related to disease risk, cases will appear more similar to other cases than to controls (or vice versa) for the SNPs associated with that pathway. By systematically applying the method to all pathways of potential interest, we can identify those for which the hypothesis holds true, i.e., pathways containing SNPs for which the samples exhibit greater within-class similarity than across classes. Importantly, PoDA improves on existing single-SNP and SNP-set enrichment analyses, in that it does not require the SNPs in a pathway to exhibit independent main effects. This permits PoDA to reveal pathways in which epistatic interactions drive risk. In this paper, we detail the PoDA method and apply it to two GWAS: one of breast cancer and the other of liver cancer. The results obtained strongly suggest that there exist pathway-wide genomic differences that contribute to disease susceptibility. PoDA thus provides an analytical tool that is complementary to existing techniques and has the power to enrich our understanding of disease genomics at the systems-level.

Distinct sets of molecular characteristics define tumor-rejecting neoantigens.

Challenges in identifying tumor-rejecting neoantigens limit the efficacy of neoantigen vaccines to treat cancers, including cutaneous squamous cell carcinoma (cSCC). A minority of human cSCC tumors shared neoantigens, supporting the need for personalized vaccines. Using a UV-induced mouse cSCC model which recapitulated the mutational signature and driver mutations found in human disease, we found that CD8 T cells constrain cSCC. Two MHC class I neoantigens were identified that constrained cSCC growth. Compared to the wild-type peptides, one tumor-rejecting neoantigen exhibited improved MHC binding and the other had increased solvent accessibility of the mutated residue. Across known neoantigens that do not impact MHC binding, structural modeling of the peptide/MHC complexes indicated that increased solvent accessibility, which will facilitate TCR recognition of the neoantigen, distinguished tumor-rejecting from non-immunogenic neoantigens. This work reveals characteristics of tumor-rejecting neoantigens that may be of considerable importance in identifying optimal vaccine candidates in cSCC and other cancers.

Poor Oral Health Is Associated With Inflammation, Aortic Valve Calcification, and Brain Volume Among Forager-Farmers.

Poor oral health is associated with cardiovascular disease and dementia. Potential pathways include sepsis from oral bacteria, systemic inflammation, and nutritional deficiencies. However, in post-industrialized populations, links between oral health and chronic disease may be confounded because the lower socioeconomic exposome (poor diet, pollution, and low physical activity) often entails insufficient dental care. We assessed tooth loss, caries, and damaged teeth, in relation to cardiovascular and brain aging among the Tsimane, a subsistence population living a relatively traditional forager-horticulturalist lifestyle with poor dental health, but minimal cardiovascular disease and dementia. Dental health was assessed by a physician in 739 participants aged 40-92 years with cardiac and brain health measured by chest computed tomography (CT; n = 728) and brain CT (n = 605). A subset of 356 individuals aged 60+ were also assessed for dementia and mild cognitive impairment (n = 33 impaired). Tooth loss was highly prevalent, with 2.2 teeth lost per decade and a 2-fold greater loss in women. The number of teeth with exposed pulp was associated with higher inflammation, as measured by cytokine levels and white blood cell counts, and lower body mass index. Coronary artery calcium and thoracic aortic calcium were not associated with tooth loss or damaged teeth. However, aortic valve calcification and brain tissue loss were higher in those who had more teeth with exposed pulp. Overall, these results suggest that dental health is associated with indicators of chronic diseases in the absence of typical confounds, even in a population with low cardiovascular and dementia risk factors.

Increases in regional brain volume across two native South American male populations.

Industrialized environments, despite benefits such as higher levels of formal education and lower rates of infections, can also have pernicious impacts upon brain atrophy. Partly for this reason, comparing age-related brain volume trajectories between industrialized and non-industrialized populations can help to suggest lifestyle correlates of brain health. The Tsimane, indigenous to the Bolivian Amazon, derive their subsistence from foraging and horticulture and are physically active. The Moseten, a mixed-ethnicity farming population, are physically active but less than the Tsimane. Within both populations (N = 1024; age range = 46-83), we calculated regional brain volumes from computed tomography and compared their cross-sectional trends with age to those of UK Biobank (UKBB) participants (N = 19,973; same age range). Surprisingly among Tsimane and Moseten (T/M) males, some parietal and occipital structures mediating visuospatial abilities exhibit small but significant increases in regional volume with age. UKBB males exhibit a steeper negative trend of regional volume with age in frontal and temporal structures compared to T/M males. However, T/M females exhibit significantly steeper rates of brain volume decrease with age compared to UKBB females, particularly for some cerebro-cortical structures (e.g., left subparietal cortex). Across the three populations, observed trends exhibit no interhemispheric asymmetry. In conclusion, the age-related rate of regional brain volume change may differ by lifestyle and sex. The lack of brain volume reduction with age is not known to exist in other human population, highlighting the putative role of lifestyle in constraining regional brain atrophy and promoting elements of non-industrialized lifestyle like higher physical activity.