A signature of Neanderthal introgression on molecular mechanisms of environmental responses

Massively parallel approaches for characterizing noncoding functional variation in human evolution

The genetic differences underlying unique phenotypes in humans compared to our closest primate relatives have long remained a mystery. Similarly, the genetic basis of adaptations between human groups during our expansion across the globe is poorly characterized. Uncovering the downstream phenotypic consequences of these genetic variants has been difficult, as a substantial portion lies in noncoding regions, such as cis-regulatory elements (CREs). Here, we review recent high-throughput approaches to measure the functions of CREs and the impact of variation within them. CRISPR screens can directly perturb CREs in the genome to understand downstream impacts on gene expression and phenotypes, while massively parallel reporter assays can decipher the regulatory impact of sequence variants. Machine learning has begun to be able to predict regulatory function from sequence alone, further scaling our ability to characterize genome function. Applying these tools across diverse phenotypes, model systems, and ancestries is beginning to revolutionize our understanding of noncoding variation underlying human evolution.

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

Genotype × environment interactions in gene regulation and complex traits

Genotype × environment interactions (GxE) have long been recognized as a key mechanism underlying human phenotypic variation. Technological developments over the past 15 years have dramatically expanded our appreciation of the role of GxE in both gene regulation and complex traits. The richness and complexity of these datasets also required parallel efforts to develop robust and sensitive statistical and computational approaches. Although our understanding of the genetic architecture of molecular and complex traits has been maturing, a large proportion of complex trait heritability remains unexplained. Furthermore, there are increasing efforts to characterize the effect of environmental exposure on human health. We therefore review GxE in human gene regulation and complex traits, advocating for a comprehensive approach that jointly considers genetic and environmental factors in human health and disease.

Association of serum 25-hydroxy-vitamin D concentration and risk of mortality in cancer survivors in the United States

PURPOSE

Cancer survivors have a high risk of mortality, and vitamin D (VD) is associated with the risk of mortality. This study is aim to examine the impact of VD on mortality in cancer survivors.

METHODS

A prospective study was conducted using data from the National Health and Nutrition Examination Survey. Participants were obtained information on their baseline characteristics, dietary habits, comorbidities, lifestyle, and serum 25-hydroxy VD [25(OH)D] concentrations. The weighted Cox proportional hazard and competing risk regression models were used to estimate the hazard ratio and 95% confidence intervals (HR, 95% CI) of mortality for different serum 25(OH)D concentrations. Restricted cubic spline (RCS) curves were utilized to illustrate the dose-response relationship between serum 25(OH)D concentrations and mortality.

RESULTS

The study encompassed 2,495 participants with cancer diagnoses. Multivariate models indicated that, compared to serum 25(OH)D concentrations below 58.5 nmol/L, concentrations exceeding 81.6 nmol/L were associated with reduced HRs for all-cause mortality (HR = 0.70; 95% CI: 0.56-0.87), cardiovascular mortality (HR = 0.53; 95% CI: 0.32-0.86), and cancer-specific mortality (HR = 0.66; 95% CI: 0.45-0.99). RCS curves revealed "L-shaped" associations between serum 25(OH)D concentration and both all-cause and cancer-specific mortality, with threshold effects at 87.9 nmol/L and 84.6 nmol/L, respectively. Conversely, the relationship between serum 25(OH)D concentration and cardiovascular mortality exhibited a more linear pattern, with a threshold at 88.7 nmol/L. Subgroup analyses highlighted a gender-specific interaction that elevated serum 25(OH)D concentrations were significantly more protective against mortality in males than in females, especially regarding cancer-specific mortality (P-interaction = 0.009).

CONCLUSION

Elevated serum 25(OH)D concentrations were correlated with decreased risks of all-cause, cardiovascular, and cancer-specific mortality in cancer survivors, with benefit thresholds at 87.9, 88.7, and 84.6 nmol/L, respectively. These findings suggested that cancer survivors might benefit from higher vitamin D recommendations than the general population.

An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression

INTRODUCTION

Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear.

MATERIALS AND METHODS

Population genetics analysis and functional genomics work was performed for this locus.

RESULTS

Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human.

CONCLUSION

This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.

Evolutionary immuno-genetics of endoplasmic reticulum aminopeptidase II (ERAP2)

Endoplasmic reticulum aminopeptidase 2 (ERAP2) is a proteolytic enzyme involved in adaptive immunity. The ERAP2 gene is highly polymorphic and encodes haplotypes that confer resistance against lethal infectious diseases, but also increase the risk for autoimmune disorders. Identifying how ERAP2 influences susceptibility to these traits requires an understanding of the selective pressures that shaped and maintained allelic variation throughout human evolution. Our review discusses the genetic regulation of haplotypes and diversity in naturally occurring ERAP2 allotypes in the global population. We outline how these ERAP2 haplotypes evolved during human history and highlight the presence of Neanderthal DNA sequences in ERAP2 of modern humans. Recent evidence suggests that human adaptation during the last ~10,000 years and historic pandemics left a significant mark on the ERAP2 gene that determines susceptibility to infectious and inflammatory diseases today.

Leveraging massively parallel reporter assays for evolutionary questions

A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.

Denisovan introgression has shaped the immune system of present-day Papuans

Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.