Exploring Adaptive Phenotypes for the Human Calcium-Sensing Receptor Polymorphism R990G.

Rainforest hunter-gatherers from Southeast Asia are characterized by specific morphological features including a particularly dark skin color (D), short stature (S), woolly hair (W), and the presence of steatopygia (S)-fat accumulation localized in the hips (DSWS phenotype). Based on previous evidence in the Andamanese population, we first characterized signatures of adaptive natural selection around the calcium-sensing receptor gene in Southeast Asian rainforest groups presenting the DSWS phenotype and identified the R990G substitution (rs1042636) as a putative adaptive variant for experimental follow-up. Although the calcium-sensing receptor has a critical role in calcium homeostasis by directly regulating the parathyroid hormone secretion, it is expressed in different tissues and has been described to be involved in many biological functions. Previous works have also characterized the R990G substitution as an activating polymorphism of the calcium-sensing receptor associated with hypocalcemia. Therefore, we generated a knock-in mouse for this substitution and investigated organismal phenotypes that could have become adaptive in rainforest hunter-gatherers from Southeast Asia. Interestingly, we found that mouse homozygous for the derived allele show not only lower serum calcium concentration but also greater body weight and fat accumulation, probably because of enhanced preadipocyte differentiation and lipolysis impairment resulting from the calcium-sensing receptor activation mediated by R990G. We speculate that such differential features in humans could have facilitated the survival of hunter-gatherer groups during periods of nutritional stress in the challenging conditions of the Southeast Asian tropical rainforests.

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.

Sleep Disturbances in At-Risk Mental States and First Episode of Psychosis: A Narrative Review on Interventions.

Sleep disturbances are a common yet often overlooked symptom of psychosis that can drastically affect the quality of life and well-being of those living with the condition. Sleep disorders are common in people diagnosed with schizophrenia and have significant negative effects on the clinical course of the illness and the functional outcomes and quality of life of patients. There is a limited number of studies addressing this question in first-episode psychosis (FEP). In this narrative review, we aimed to provide an overview of sleep disorders in populations with FEP and at-risk mental states (ARMS). The review was focused on the various treatments currently used for sleep disorders, including both non-pharmacological and pharmacological treatments. A total of 48 studies were included. We found that sleep disturbances are associated with attenuated psychotic symptoms and other psychopathological symptoms in ARMSs. The association of sleep disturbances with the transition to psychosis has been poorly investigated. Sleep disturbances have an impact on the quality of life and the psychopathological symptoms of people suffering from FEP. The non-pharmacological treatments include cognitive behavioral therapy for insomnia, bright light therapy, cognitive restructuring techniques, sleep restriction therapy, basic sleep hygiene education, and the provision of portable sleep trackers. Other treatments include antipsychotics in acute phases and melatonin. The early intervention in sleep disturbances may improve overall prognosis in emerging psychosis populations.

Identifying signatures of positive selection in human populations from North Africa.

Because of its location, North Africa (NA) has witnessed continuous demographic movements with an impact on the genomes of present-day human populations. Genomic data describe a complex scenario with varying proportions of at least four main ancestry components: Maghrebi, Middle Eastern-, European-, and West-and-East-African-like. However, the footprint of positive selection in NA has not been studied. Here, we compile genome-wide genotyping data from 190 North Africans and individuals from surrounding populations, investigate for signatures of positive selection using allele frequencies and linkage disequilibrium-based methods and infer ancestry proportions to discern adaptive admixture from post-admixture selection events. Our results show private candidate genes for selection in NA involved in insulin processing (KIF5A), immune function (KIF5A, IL1RN, TLR3), and haemoglobin phenotypes (BCL11A). We also detect signatures of positive selection related to skin pigmentation (SLC24A5, KITLG), and immunity function (IL1R1, CD44, JAK1) shared with European populations and candidate genes associated with haemoglobin phenotypes (HPSE2, HBE1, HBG2), other immune-related (DOCK2) traits, and insulin processing (GLIS3) traits shared with West and East African populations. Finally, the SLC8A1 gene, which codifies for a sodium-calcium exchanger, was the only candidate identified under post-admixture selection in Western NA.

Polygenic risk scores enhance prediction of body mass index increase in individuals with a first episode of psychosis.

BACKGROUND: Individuals with a first episode of psychosis (FEP) show rapid weight gain during the first months of treatment, which is associated with a reduction in general physical health. Although genetics is assumed to be a significant contributor to weight gain, its exact role is unknown. METHODS: We assembled a population-based FEP cohort of 381 individuals that was split into a Training (n = 224) set and a Validation (n = 157) set to calculate the polygenic risk score (PRS) in a two-step process. In parallel, we obtained reference genome-wide association studies for body mass index (BMI) and schizophrenia (SCZ) to examine the pleiotropic landscape between the two traits. BMI PRSs were added to linear models that included sociodemographic and clinical variables to predict BMI increase (∆BMI) in the Validation set. RESULTS: The results confirmed considerable shared genetic susceptibility for the two traits involving 449 near-independent genomic loci. The inclusion of BMI PRSs significantly improved the prediction of ∆BMI at 12 months after the onset of antipsychotic treatment by 49.4% compared to a clinical model. In addition, we demonstrated that the PRS containing pleiotropic information between BMI and SCZ predicted ∆BMI better at 3 (12.2%) and 12 months (53.2%). CONCLUSIONS: We prove for the first time that genetic factors play a key role in determining ∆BMI during the FEP. This finding has important clinical implications for the early identification of individuals most vulnerable to weight gain and highlights the importance of examining genetic pleiotropy in the context of medically important comorbidities for predicting future outcomes.

Uncovering Signals of Positive Selection in Peruvian Populations from Three Ecological Regions.

Peru hosts extremely diverse ecosystems which can be broadly classified into the following three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide single-nucleotide polymorphism-array data from 286 Peruvians living across the three ecoregions and analyzed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44, and ICAM1 in the coast; CBLB and PRDM1 in the rainforest; and BRD2, HLA-DOA, HLA-DPA1 regions in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

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.

Contribution of Evolutionary Selected Immune Gene Polymorphism to Immune-Related Disorders: The Case of Lymphocyte Scavenger Receptors CD5 and CD6.

Pathogens are one of the main selective pressures that ancestral humans had to adapt to. Components of the immune response system have been preferential targets of natural selection in response to such pathogen-driven pressure. In turn, there is compelling evidence showing that positively selected immune gene variants conferring increased resistance to past or present infectious agents are today associated with increased risk for autoimmune or inflammatory disorders but decreased risk of cancer, the other side of the same coin. CD5 and CD6 are lymphocytic scavenger receptors at the interphase of the innate and adaptive immune responses since they are involved in both: (i) microbial-associated pattern recognition; and (ii) modulation of intracellular signals mediated by the clonotypic antigen-specific receptor present in T and B cells (TCR and BCR, respectively). Here, we review available information on CD5 and CD6 as targets of natural selection as well as on the role of CD5 and CD6 variation in autoimmunity and cancer.

DDR1 methylation is associated with bipolar disorder and the isoform expression and methylation of myelin genes.

Aim: To investigate DDR1 methylation in the brains of bipolar disorder (BD) patients and its association with DDR1 mRNA levels and comethylation with myelin genes. Materials & methods: Genome-wide profiling of DNA methylation (Infinium MethylationEPIC BeadChip) corrected for glial composition and DDR1 gene expression analysis in the occipital cortices of individuals with BD (n = 15) and healthy controls (n = 15) were conducted. Results:DDR1 5-methylcytosine levels were increased and directly associated with DDR1b mRNA expression in the brains of BD patients. We also observed that DDR1 was comethylated with a group of myelin genes. Conclusion:DDR1 is hypermethylated in BD brain tissue and is associated with isoform expression. Additionally, DDR1 comethylation with myelin genes supports the role of this receptor in myelination.

The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm.

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.

The shared genetic architecture of schizophrenia, bipolar disorder and lifespan.

Psychiatric disorders such as Schizophrenia (SCZ) and Bipolar Disorder (BD) represent an evolutionary paradox, as they exhibit strong negative effects on fitness, such as decreased fecundity and early mortality, yet they persist at a worldwide prevalence of approximately 1%. Molecular mechanisms affecting lifespan, which may be widely common among complex diseases with fitness effects, can be studied by the integrated analysis of data from genome-wide association studies (GWAS) of human longevity together with any disease of interest. Here, we report the first of such studies, focusing on the genetic overlap-pleiotropy-between two psychiatric disorders with shortened lifespan, SCZ and BD, and human parental lifespan (PLS) as a surrogate of life expectancy. Our results are twofold: first, we demonstrate extensive polygenic overlap between SCZ and PLS and to a lesser extent between BD and PLS. Second, we identified novel loci shared between PLS and SCZ (n = 39), and BD (n = 8). Whereas most of the identified SCZ (66%) and BD (62%) pleiotropic risk alleles were associated with reduced lifespan, we also detected some antagonistic protective alleles associated to shorter lifespans. In fact, top-associated SNPs with SCZ seems to explain longevity variance explained (LVE) better than many other life-threatening diseases, including Type 2 diabetes and most cancers, probably due to a high overlap with smoking-related pathways. Overall, our study provides evidence of a genetic burden driven through premature mortality among people with SCZ, which can have profound implications for understanding, and potentially treating, the mortality gap associated with this psychiatric disorder.

Adaptive selection drives TRPP3 loss-of-function in an Ethiopian population.

TRPP3 (also called PKD2L1) is a nonselective, cation-permeable channel activated by multiple stimuli, including extracellular pH changes. TRPP3 had been considered a candidate for sour sensor in humans, due to its high expression in a subset of tongue receptor cells detecting sour, along with its membership to the TRP channel family known to function as sensory receptors. Here, we describe the functional consequences of two non-synonymous genetic variants (R278Q and R378W) found to be under strong positive selection in an Ethiopian population, the Gumuz. Electrophysiological studies and 3D modelling reveal TRPP3 loss-of-functions produced by both substitutions. R278Q impairs TRPP3 activation after alkalinisation by mislocation of H+ binding residues at the extracellular polycystin mucolipin domain. R378W dramatically reduces channel activity by altering conformation of the voltage sensor domain and hampering channel transition from closed to open state. Sour sensitivity tests in R278Q/R378W carriers argue against both any involvement of TRPP3 in sour detection and the role of such physiological process in the reported evolutionary positive selection past event.

The shaping of immunological responses through natural selection after the Roma Diaspora.

The Roma people are the largest transnational ethnic minority in Europe and can be considered the last human migration of South Asian origin into the continent. They left Northwest India approximately 1,000 years ago, reaching the Balkan Peninsula around the twelfth century and Romania in the fourteenth century. Here, we analyze whole-genome sequencing data of 40 Roma and 40 non-Roma individuals from Romania. We performed a genome-wide scan of selection comparing Roma, their local host population, and a Northwestern Indian population, to identify the selective pressures faced by the Roma mainly after they settled in Europe. We identify under recent selection several pathways implicated in immune responses, among them cellular metabolism pathways known to be rewired after immune stimulation. We validated the interaction between PIK3-mTOR-HIF-1α and cytokine response influenced by bacterial and fungal infections. Our results point to a significant role of these pathways for host defense against the most prevalent pathogens in Europe during the last millennium.

A New Risk Variant for Multiple Sclerosis at 11q23.3 Locus Is Associated with Expansion of CXCR5+ Circulating Regulatory T Cells.

Genome-wide association studies and meta-analysis have contributed to the identification of more than 200 loci associated with multiple sclerosis (MS). However, a proportion of MS heritability remains unknown. We aimed to uncover new genetic variants associated with MS and determine their functional effects. For this, we resequenced the exons and regulatory sequences of 14 MS risk genes in a cohort of MS patients and healthy individuals (n = 1,070) and attempted to validate a selection of signals through genotyping in an independent cohort (n = 5,138). We identified three new MS-associated variants at C-X-C motif chemokine receptor 5 (CXCR5), Ts translation elongation factor, mitochondrial (TSFM) and cytochrome P450 family 24 subfamily A member 1 (CYP24A1). Rs10892307 resulted in a new signal at the CXCR5 region that explains one of the associations with MS within the locus. This polymorphism and three others in high linkage disequilibrium mapped within regulatory regions. Of them, rs11602393 showed allele-dependent enhancer activity in the forward orientation as determined by luciferase reporter assays. Immunophenotyping using peripheral blood mononuclear cells from MS patients associated the minor allele of rs10892307 with increased percentage of regulatory T cells expressing CXCR5. This work reports a new signal for the CXCR5 MS risk locus and points to rs11602393 as the causal variant. The expansion of CXCR5+ circulating regulatory T cells induced by this variant could cause its MS association.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes.

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-ß compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

Macrophage-specific MHCII expression is regulated by a remote Ciita enhancer controlled by NFAT5.

MHCII in antigen-presenting cells (APCs) is a key regulator of adaptive immune responses. Expression of MHCII genes is controlled by the transcription coactivator CIITA, itself regulated through cell type-specific promoters. Here we show that the transcription factor NFAT5 is needed for expression of Ciita and MHCII in macrophages, but not in dendritic cells and other APCs. NFAT5-deficient macrophages showed defective activation of MHCII-dependent responses in CD4+ T lymphocytes and attenuated capacity to elicit graft rejection in vivo. Ultrasequencing analysis of NFAT5-immunoprecipitated chromatin uncovered an NFAT5-regulated region distally upstream of Ciita This region was required for CIITA and hence MHCII expression, exhibited NFAT5-dependent characteristics of active enhancers such as H3K27 acetylation marks, and required NFAT5 to interact with Ciita myeloid promoter I. Our results uncover an NFAT5-regulated mechanism that maintains CIITA and MHCII expression in macrophages and thus modulates their T lymphocyte priming capacity.

Sequence diversity of the Rh blood group system in Basques.

Basques show specific cultural, demographic, and genetic characteristics that have placed them as an isolated and unique population within Europe, such as their non-Indo-European language, Euskara. They have historically lived along the Western Pyrenees, between Spain and France, in one of the most important European glacial refugia during the Last Glacial Maximum. The most striking genetic characteristic is their highest frequency of the RhD blood group negative allele, a variant related to the hemolytic disease of the newborn. Both demographic and adaptive processes have been suggested as possible causes of the high frequency of RhD negative in Basques, but neither hypothesis has been clearly demonstrated. While previous studies on the Rh system in Basques have been mostly focused on serological and genotyping diversity, in this work we analyze genotyping and next generation sequencing data in order to provide a general framework of the genetic scenario of the system in Basques. In particular, we genotyped the most relevant variants of the system (D/d, E/e, and C/c), and sequenced three ~6 kb flanking regions surrounding the Rh genes in Basques and also in other populations for comparison. Our results are in agreement with previous studies, with Basques presenting the highest frequency of the RHD deletion (47.2%). Haplotype analyses of D/d, E/e, and C/c variants confirmed an association between the RhC allele, previously suggested to be under positive selection, and the RhD positive variant in non-sub-Saharan populations, including Basques. We also found extreme differentiation for the C/c variant when comparing sub-Saharan to non-sub-Saharan populations.

Evaluating the Genetics of Common Variable Immunodeficiency: Monogenetic Model and Beyond.

Common variable immunodeficiency (CVID) is the most frequent symptomatic primary immunodeficiency characterized by recurrent infections, hypogammaglobulinemia and poor response to vaccines. Its diagnosis is made based on clinical and immunological criteria, after exclusion of other diseases that can cause similar phenotypes. Currently, less than 20% of cases of CVID have a known underlying genetic cause. We have analyzed whole-exome sequencing and copy number variants data of 36 children and adolescents diagnosed with CVID and healthy relatives to estimate the proportion of monogenic cases. We have replicated an association of CVID to p.C104R in TNFRSF13B and reported the second case of homozygous patient to date. Our results also identify five causative genetic variants in LRBA, CTLA4, NFKB1, and PIK3R1, as well as other very likely causative variants in PRKCD, MAPK8, or DOCK8 among others. We experimentally validate the effect of the LRBA stop-gain mutation which abolishes protein production and downregulates the expression of CTLA4, and of the frameshift indel in CTLA4 producing expression downregulation of the protein. Our results indicate a monogenic origin of at least 15-24% of the CVID cases included in the study. The proportion of monogenic patients seems to be lower in CVID than in other PID that have also been analyzed by whole exome or targeted gene panels sequencing. Regardless of the exact proportion of CVID monogenic cases, other genetic models have to be considered for CVID. We propose that because of its prevalence and other features as intermediate penetrancies and phenotypic variation within families, CVID could fit with other more complex genetic scenarios. In particular, in this work, we explore the possibility of CVID being originated by an oligogenic model with the presence of heterozygous mutations in interacting proteins or by the accumulation of detrimental variants in particular immunological pathways, as well as perform association tests to detect association with rare genetic functional variation in the CVID cohort compared to healthy controls.