Signatures of Evolutionary Adaptation in Quantitative Trait Loci Influencing Trace Element Homeostasis in Liver

Association between Multiple Trace Elements, Executive Function, and Cognitive Impairment with No Dementia in Older Adults

Many studies suggest a significant association between individual essential trace elements (ETEs) and cognitive impairment in older adults, but evidence of the synchronized effect of multiple ETEs on cognitive function is lacking. We investigated the association between multiple ETEs, cognitive impairment with no dementia (CIND), and executive function in older Korean adults, using the Bayesian kernel machine regression (BKMR) model. Three hundred and thirty-six older adults were included as the study population and classified as the CIND and control groups. Blood manganese (Mn), copper (Cu), zinc (Zn), selenium (Se), and molybdenum (Mo) were measured as relevant ETEs. The frontal/executive tests included digit symbol coding (DSC), the Korean color word Stroop test (K-CWST), a controlled oral word association test (COWAT), and a trial-making test (TMT). Overall, the BKMR showed a negative association between multiple ETEs and the odds of CIND. Mn was designated as the most dominant element associated with the CIND (PIP = 0.6184), with a U-shaped relationship. Cu and Se levels were positively associated with the K-CWST percentiles (β = 31.78; 95% CI: 13.51, 50.06) and DSC percentiles (β = 25.10; 95% CI: 7.66, 42.53), respectively. Our results suggest that exposure to multiple ETEs may be linked to a protective mechanism against cognitive impairment in older adults.

Human genetic adaptation related to cellular zinc homeostasis

SLC30A9 encodes a ubiquitously zinc transporter (ZnT9) and has been consistently suggested as a candidate for positive selection in humans. However, no direct adaptive molecular phenotype has been demonstrated. Our results provide evidence for directional selection operating in two major complementary haplotypes in Africa and East Asia. These haplotypes are associated with differential gene expression but also differ in the Met50Val substitution (rs1047626) in ZnT9, which we show is found in homozygosis in the Denisovan genome and displays accompanying signatures suggestive of archaic introgression. Although we found no significant differences in systemic zinc content between individuals with different rs1047626 genotypes, we demonstrate that the expression of the derived isoform (ZnT9 50Val) in HEK293 cells shows a gain of function when compared with the ancestral (ZnT9 50Met) variant. Notably, the ZnT9 50Val variant was found associated with differences in zinc handling by the mitochondria and endoplasmic reticulum, with an impact on mitochondrial metabolism. Given the essential role of the mitochondria in skeletal muscle and since the derived allele at rs1047626 is known to be associated with greater susceptibility to several neuropsychiatric traits, we propose that adaptation to cold may have driven this selection event, while also impacting predisposition to neuropsychiatric disorders in modern humans.

Spatial and temporal diversity of positive selection on shared haplotypes at the PSCA locus among worldwide human populations

Selection on standing genetic variation is important for rapid local genetic adaptation when the environment changes. We report that, for the prostate stem cell antigen (PSCA) gene, different populations have different target haplotypes, even though haplotypes are shared among populations. The C-C-A haplotype, whereby the first C is located at rs2294008 of PSCA and is a low risk allele for gastric cancer, has become a target of positive selection in Asia. Conversely, the C-A-G haplotype carrying the same C allele has become a selection target mainly in Africa. However, Asian and African share both haplotypes, consistent with the haplotype divergence time (170 kya) prior to the out-of-Africa dispersal. The frequency of C-C-A/C-A-G is 0.344/0.278 in Asia and 0.209/0.416 in Africa. Two-dimensional site frequency spectrum analysis revealed that the extent of intra-allelic variability of the target haplotype is extremely small in each local population, suggesting that C-C-A or C-A-G is under ongoing hard sweeps in local populations. From the time to the most recent common ancestor (TMRCA) of selected haplotypes, the onset times of positive selection were recent (3-55 kya), concurrently with population subdivision from a common ancestor. Additionally, estimated selection coefficients from ABC analysis were up to ~3%, similar to those at other loci under recent positive selection. Phylogeny of local populations and TMRCA of selected haplotypes revealed that spatial and temporal switching of positive selection targets is a unique and novel feature of ongoing selection at PSCA. This switching may reflect the potential of rapid adaptability to distinct environments.

Association between whole blood essential trace elements and cognitive function in older adults

BACKGROUND

Essential trace elements (ETEs) are essential nutrients for keeping the nervous system functioning. Associations between ETEs and cognitive function are still inconclusive and limited.

OBJECTIVES

We aimed to investigate the individual and joint associations between ETEs and cognitive function among older adults.

METHODS

A population (N = 2181) at an average age≥ 65 from Yiwu cohort in China was available for this study. Whole blood chromium (Cr), selenium (Se), manganese (Mn), and copper (Cu) concentrations were measured by inductively coupled plasma mass spectrometry (ICP-MS). Cognitive function was assessed using the Mini-Mental State Examination (MMSE), consisting of five specific cognitive domains: orientation, registry, attention and calculation, recall, and language and praxis. Linear regression, restricted cubic spline (RCS) analysis, and Bayesian kernel machine regression (BKMR) were used to analyze the individual and joint associations between ETEs and cognitive function.

RESULTS

The association between Cr and MMSE score presented an inverted-U shape (Q3 versus Q1: β = 0.774, 95 % CI: 0.297, 1.250; Q4 versus Q1: β = 0.481, 95 % CI: 0.006, 0.956); and Cr was especially associated with the registry, recall, and language and praxis. Per IQR (36.32 μg/L) increase of Se was positively associated with the MMSE score (β = 0.497, 95 % CI: 0.277, 0.717) and all five cognitive domains. The BKMR showed that the dose-response association between Se and cognitive function increased initially and then decreased with increasing Se concentration when fixed the other ETEs in median. ETEs mixture was positively associated with cognitive function, and Se (posterior inclusion probabilities, PIPs = 0.915) was the most important contributor within the ETEs mixture.

CONCLUSIONS

The nonlinear association between Cr and cognitive function suggested further exploration of an appropriate concentration range for ETEs. A positive association between mixed ETEs and cognitive function is a reminder that their joint association should be considered. Further prospective studies or intervention studies are warranted to validate our findings in the future.

Metal Exposure Promotes Colorectal Tumorigenesis via the Aberrant N6-Methyladenosine Modification of ATP13A3

Element contamination, including that from heavy metals, is associated with gastrointestinal tumorigenesis, but the effects and mechanisms of crucial element exposure associated with colorectal cancer remain unclear. We profiled 56 elements by ICP-MS and used logistic regression, LASSO, BKMR, and GAM to identify colorectal cancer-relevant elements. A series of biochemical experiments were performed to demonstrate the cytotoxicity and the mechanisms of malignant transformation after metal exposure. Using an elementomics approach, we first found that the metal thallium (Tl) was positively correlated with many toxic metals and was associated with a significantly increased risk of colorectal cancer. Acute exposure to Tl induced cytotoxicity and cell death by accelerating the generation of reactive oxygen species and DNA damage. Chronic exposure to Tl led to the inhibition of cell death and thereby induced the malignant transformation of normal colon cells and xenograft tumor formation in nude mice. Furthermore, we describe the first identification of a significant metal quantitative trait locus for the novel colorectal cancer susceptibility locus rs1511625 near ATP13A3. Mechanistically, Tl increased the level of aberrant N⁶-methyladenosine (m⁶A) modification of ATP13A3 via the METLL3/METTL14/ALKBH5-ATP13A3 axis to promote colorectal tumorigenesis. This study provides a basis for the development of public health strategies for reducing metal exposure among populations at high risk for colorectal cancer.

Prospects of microgreens as budding living functional food: Breeding and biofortification through OMICS and other approaches for nutritional security

Nutrient deficiency has resulted in impaired growth and development of the population globally. Microgreens are considered immature greens (required light for photosynthesis and growing medium) and developed from the seeds of vegetables, legumes, herbs, and cereals. These are considered "living superfood/functional food" due to the presence of chlorophyll, beta carotene, lutein, and minerals like magnesium (Mg), Potassium (K), Phosphorus (P), and Calcium (Ca). Microgreens are rich at the nutritional level and contain several phytoactive compounds (carotenoids, phenols, glucosinolates, polysterols) that are helpful for human health on Earth and in space due to their anti-microbial, anti-inflammatory, antioxidant, and anti-carcinogenic properties. Microgreens can be used as plant-based nutritive vegetarian foods that will be fruitful as a nourishing constituent in the food industryfor garnish purposes, complement flavor, texture, and color to salads, soups, flat-breads, pizzas, and sandwiches (substitute to lettuce in tacos, sandwich, burger). Good handling practices may enhance microgreens'stability, storage, and shelf-life under appropriate conditions, including light, temperature, nutrients, humidity, and substrate. Moreover, the substrate may be a nutritive liquid solution (hydroponic system) or solid medium (coco peat, coconut fiber, coir dust and husks, sand, vermicompost, sugarcane filter cake, etc.) based on a variety of microgreens. However integrated multiomics approaches alongwith nutriomics and foodomics may be explored and utilized to identify and breed most potential microgreen genotypes, biofortify including increasing the nutritional content (macro-elements:K, Ca and Mg; oligo-elements: Fe and Zn and antioxidant activity) and microgreens related other traits viz., fast growth, good nutritional values, high germination percentage, and appropriate shelf-life through the implementation of integrated approaches includes genomics, transcriptomics, sequencing-based approaches, molecular breeding, machine learning, nanoparticles, and seed priming strategiesetc.

Associations Between Mild Cognitive Impairment and Whole Blood Zinc and Selenium in the Elderly Cohort

Some studies have shown that an imbalance in trace element homeostasis can lead to cognitive dysfunction, but data are lacking. The purpose of this study was to investigate the association between whole blood zinc (Zn), selenium (Se), copper-zinc ratio (Cu/Zn), copper-selenium ratio (Cu/Se), and zinc-selenium ratio (Zn/Se) and mild cognitive impairment (MCI) in elderly Chinese individuals. The study was based on the Elderly Health and Controlled Environmental Factors Cohort in Lu'an, Anhui Province, China, from June to September 2016. The cognitive function of the elderly was determined by the Mini-Mental State Examination (MMSE) and activities of daily living (ADL) scales. The concentrations of Zn, Cu, and Se in the whole blood were measured by inductively coupled plasma-mass spectrometry (ICP-MS). Binary logistic regression was used to analyze the associations between trace elements and MCI. A total of 1006 participants with an average age of 71.70 years old were included in this study. Compared with healthy people, MCI patients had higher whole blood Zn levels and lower Se levels, and Cu/Zn, Cu/Se, and Zn/Se were also significantly different. Binary logistic regression analysis showed that Zn, Cu/Se, and Zn/Se exposure in the third tertile was associated with an increased risk of MCI, while Se exposure in the third tertile was associated with a reduced risk of MCI. After adjustment for sex, age, marital status, BMI, and living status, whole blood Zn, Se, Cu/Zn, Cu/Se, and Zn/Se were significantly associated with MCI risk, especially in elderly women.

Understanding signatures of positive natural selection in human zinc transporter genes

Zinc is an essential micronutrient with a tightly regulated systemic and cellular homeostasis. In humans, some zinc transporter genes (ZTGs) have been previously reported as candidates for strong geographically restricted selective sweeps. However, since zinc homeostasis is maintained by the joint action of 24 ZTGs, other more subtle modes of selection could have also facilitated human adaptation to zinc availability. Here, we studied whether the complete set of ZTGs are enriched for signals of positive selection in worldwide populations and population groups from South Asia. ZTGs showed higher levels of genetic differentiation between African and non-African populations than would be randomly expected, as well as other signals of polygenic selection outside Africa. Moreover, in several South Asian population groups, ZTGs were significantly enriched for SNPs with unusually extended haplotypes and displayed SNP genotype-environmental correlations when considering zinc deficiency levels in soil in that geographical area. Our study replicated some well-characterized targets for positive selection in East Asia and sub-Saharan Africa, and proposes new candidates for follow-up in South Asia (SLC39A5) and Africa (SLC39A7). Finally, we identified candidate variants for adaptation in ZTGs that could contribute to different disease susceptibilities and zinc-related human health traits.

Tracing the Evolution of Human Gene Regulation and Its Association with Shifts in Environment

As humans populated the world, they adapted to many varying environmental factors, including climate, diet, and pathogens. Because many of these adaptations were mediated by multiple noncoding variants with small effects on gene regulation, it has been difficult to link genomic signals of selection to specific genes, and to describe the regulatory response to selection. To overcome this challenge, we adapted PrediXcan, a machine learning method for imputing gene regulation from genotype data, to analyze low-coverage ancient human DNA (aDNA). First, we used simulated genomes to benchmark strategies for adapting PrediXcan to increase robustness to incomplete data. Applying the resulting models to 490 ancient Eurasians, we found that genes with the strongest divergent regulation among ancient populations with hunter-gatherer, pastoralist, and agricultural lifestyles are enriched for metabolic and immune functions. Next, we explored the contribution of divergent gene regulation to two traits with strong evidence of recent adaptation: dietary metabolism and skin pigmentation. We found enrichment for divergent regulation among genes proposed to be involved in diet-related local adaptation, and the predicted effects on regulation often suggest explanations for known signals of selection, for example, at FADS1, GPX1, and LEPR. In contrast, skin pigmentation genes show little regulatory change over a 38,000-year time series of 2,999 ancient Europeans, suggesting that adaptation mainly involved large-effect coding variants. This work demonstrates that combining aDNA with present-day genomes is informative about the biological differences among ancient populations, the role of gene regulation in adaptation, and the relationship between genetic diversity and complex traits.

Pathogenic Variants in Selenoproteins and Selenocysteine Biosynthesis Machinery

Selenium is incorporated into selenoproteins as the 21st amino acid selenocysteine (Sec). There are 25 selenoproteins encoded in the human genome, and their synthesis requires a dedicated machinery. Most selenoproteins are oxidoreductases with important functions in human health. A number of disorders have been associated with deficiency of selenoproteins, caused by mutations in selenoprotein genes or Sec machinery genes. We discuss mutations that are known to cause disease in humans and report their allele frequencies in the general population. The occurrence of protein-truncating variants in the same genes is also presented. We provide an overview of pathogenic variants in selenoproteins genes from a population genomics perspective.

Ferroptosis-related gene signature associates with immunity and predicts prognosis accurately in patients with osteosarcoma

Osteosarcoma has been the most common malignant bone tumor in children and adolescents, while the 5‐y survival of osteosarcoma patients gained no significant improvement over the past decades. This study aimed to explore the role of ferroptosis‐related genes (FRGs) in the development and prognosis of osteosarcoma. The datasets of osteosarcoma patients including RNA sequencing data and clinical information were acquired from the TRGET and Gene Expression Omnibus (GEO) databases. The identification of molecular subgroups with different FRG expression patterns was achieved through nonnegative matrix factorization (NMF) clustering. The prognostic model was constructed using the least absolute shrinkage and selection operator (LASSO) algorithm and multivariate Cox regression analysis. The ESTIMATE algorithm was applied for determining the stromal score, immune score, ESTIMA score, and tumor purity of osteosarcoma patients. Functional analyses including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, Gene Set Enrichment Analysis (GSEA), and Gene Set Variation Analysis (GSVA) were conducted to explore the underlying mechanisms in the development and prognosis of osteosarcoma. Two molecular subgroups with different FRGs expression patterns were identified. The molecular subgroups with higher immune score and more active immune status showed better prognostic survival. On the basis of FRGs, a prognostic model and a nomogram integrating clinical characteristics were constructed and their prediction efficiency for osteosarcoma prognosis were well validated. Gene functional enrichment analysis showed that these differentially expressed FRGs were mainly enriched in immunity‐related signaling pathways, indicating that FRGs may affect the development and prognosis of osteosarcoma by regulating the immune microenvironment. The expression profiles of FRGs were closely related to the immunity status and prognostic survival of osteosarcoma patients. The interaction between ferroptosis and immunity in the development of osteosarcoma could provide a new insight into the exploration of molecular mechanisms and targeted therapies of osteosarcoma patients.

Knowledge-based analyses reveal new candidate genes associated with risk of hepatitis B virus related hepatocellular carcinoma

BACKGROUND

Recent genome-wide association studies (GWASs) have suggested several susceptibility loci of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) by statistical analysis at individual single-nucleotide polymorphisms (SNPs). However, these loci only explain a small fraction of HBV-related HCC heritability. In the present study, we aimed to identify additional susceptibility loci of HBV-related HCC using advanced knowledge-based analysis.

METHODS

We performed knowledge-based analysis (including gene- and gene-set-based association tests) on variant-level association p-values from two existing GWASs of HBV-related HCC. Five different types of gene-sets were collected for the association analysis. A number of SNPs within the gene prioritized by the knowledge-based association tests were selected to replicate genetic associations in an independent sample of 965 cases and 923 controls.

RESULTS

The gene-based association analysis detected four genes significantly or suggestively associated with HBV-related HCC risk: SLC39A8, GOLGA8M, SMIM31, and WHAMMP2. The gene-set-based association analysis prioritized two promising gene sets for HCC, cell cycle G1/S transition and NOTCH1 intracellular domain regulates transcription. Within the gene sets, three promising candidate genes (CDC45, NCOR1 and KAT2A) were further prioritized for HCC. Among genes of liver-specific expression, multiple genes previously implicated in HCC were also highlighted. However, probably due to small sample size, none of the genes prioritized by the knowledge-based association analyses were successfully replicated by variant-level association test in the independent sample.

CONCLUSIONS

This comprehensive knowledge-based association mining study suggested several promising genes and gene-sets associated with HBV-related HCC risks, which would facilitate follow-up functional studies on the pathogenic mechanism of HCC.

SLC39A8 is a risk factor for schizophrenia in Uygur Chinese: a case-control study

BACKGROUND

Schizophrenia is a severe mental disease with high morbidity and heritability. The SLC39A8 gene is located in 4q24 and encodes a protein that transports many metal ions. Multiple previous studies found that one of the most pleiotropic single nucleotide polymorphisms (SNPs) in SLC39A8, rs13107325, is associated with schizophrenia in the European population. However, the polymorphism of this locus is rare in other populations. In China, the Han Chinese and the Uygur Chinese are two ethnic populations that originate from different races.

METHODS

A case-control study was conducted with 983 schizophrenia cases and 1230 healthy controls of the Chinese Uygur population. To validate the most promising SNP, meta-analyses were conducted with the Han Chinese and the European PGC2 data sets reported previously.

RESULTS

A susceptible locus, rs10014145 (pallele = 0.014, pallele = 0.098 after correction; pgenotype = 0.004, pgenotype = 0.032 after correction) was identified in case-control study of the Chinese Uygur population. Further, the association between rs10014145 and schizophrenia was supported by a meta-analysis of Han and Uygur Chinese samples (pooled OR [95% CI] =1.10 [1.03-1.17], Z = 2.73, p = 0.006). The association between rs10014145 and schizophrenia was not significant in a meta-analysis of combined Chinese and European samples (pooled OR [95% CI] =1.07 [1.00-1.14], Z = 1.88, and p = 0.06). In addition, the "CCAC" haplotype of rs4698844-rs233814-rs13114343-rs151394 was significantly associated with schizophrenia in Uygur Chinese (P = 0.003, corrected p = 0.012).

CONCLUSIONS

The results of this study support that SLC39A8 is a susceptible gene for schizophrenia in the populations of Han Chinese and Uygur Chinese in China, further studies are suggested to validate the association.

Is there adaptation in the human genome for taste perception and phase I biotransformation?

Background

During the modern human expansion, new environmental pressures may have driven adaptation, especially in genes related to the perception of ingested substances and their detoxification. Consequently, positive (adaptive) selection may have occurred in genes related to taste, and in those related to the CYP450 system due to its role in biotransformation of potentially toxic compounds. A total of 91 genes (taste receptors and CYP450 superfamily) have been studied using Hierarchical Boosting, a powerful combination of different selection tests, to detect signatures of recent positive selection in three continental human populations: Northern Europeans (CEU), East Asians (CHB) and Africans (YRI). Analyses have been refined with selection analyses of the 26 populations of 1000 Genomes Project Phase 3.

Results

Genes related to taste perception have not been positively selected in the three continental human populations. This finding suggests that, contrary to results of previous studies, different allele frequencies among populations in genes such as TAS2R38 and TAS2R16 are not due to positive selection but to genetic drift. CYP1 and CYP2 genes, also previously considered to be under positive selection, did not show signatures of selective sweeps. However, three genes belonging to the CYP450 system have been identified by the Hierarchical Boosting as positively selected: CYP3A4 and CYP3A43 in CEU, and CYP27A1 in CHB.

Conclusions

No main adaptive differences are found in known taste receptor genes among the three continental human populations studied. However, there are important genetic adaptations in the cytochrome P450 system related to the Out of Africa expansion of modern humans. We confirmed that CYP3A4 and CYP3A43 are under selection in CEU, and we report for the first time CYP27A1 to be under positive selection in CHB.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1366-7) contains supplementary material, which is available to authorized users.

Using phenome-wide association to investigate the function of a schizophrenia risk locus at SLC39A8

While nearly all common genomic variants associated with schizophrenia have no known function, one corresponds to a missense variant associated with change in efficiency of a metal ion transporter, ZIP8, coded by SLC39A8. This variant has been linked to a range of phenotypes and is believed to be under recent selection pressure, but its impact on health is poorly understood. We sought to understand phenotypic implications of this variant in a large genomic biobank using an unbiased phenome-wide approach. Specifically, we generated 50 topics based on diagnostic codes using latent Dirichlet allocation, and examined them for association with the risk variant. Then, any significant topics were further characterized by examining association with individual diagnostic codes contributing to the topic. Among 50 topics, 1 was associated at an experiment-wide significance threshold (beta = 0.003, uncorrected p = 0.00049), comprising predominantly brain-related codes, including intracranial hemorrhage, cerebrovascular disease, and delirium/dementia. These results suggest that a functional variant previously associated with schizophrenia risk also increases liability to cerebrovascular disease. They further illustrate the utility of a topic-based approach to phenome-wide association.

The highly pleiotropic gene SLC39A8 as an opportunity to gain insight into the molecular pathogenesis of schizophrenia

There is a long way from the initial discovery of a genome-wide significant signal to mechanistic understanding of the association. Identification of the gene and causal polymorphism usually requires an extensive additional effort. The schizophrenia genome-wide significant locus at 4q24 may be a rare exception to this pattern. As discussed in this review, the association at this locus is most probably driven by a functional missense variant at the metal cations transporter SLC39A8. The variant, rs13107325, is almost exclusive of European populations and is one of the most pleiotropic variants of the genome, being associated at genome-wide significant level with several additional traits, such as body mass index, Crohn's disease, blood pressure related-traits, and serum levels of manganese, N-terminal pro-B-type natriuretic peptide and HDL-cholesterol. SLC39A8 seems to be subject to recent natural selection in Europeans. It is almost ubiquitously expressed and its physiological role is beginning to be elucidated, mainly in relation to immunity. This manuscript presents arguments in favor of the rs13107325 variant as the functional variant responsible for the association of this locus with schizophrenia, reviews the genetic associations with this gene, the evidences of natural selection on the gene, and the known aspects about its structure and physiological functions. Finally, some hypotheses about putative mechanisms for its association with schizophrenia are presented based on this knowledge, including impaired immunity/inflammation, interference with glutamatergic neurotransmission, homeostasis of essential metals in brain, such as iron, zinc or manganese, or neurotoxicity by heavy metals, such as cadmium or lead.