Emergence of saliva protein genes in the secretory calcium-binding phosphoprotein (SCPP) locus and accelerated evolution in primates.

Genes within the secretory calcium-binding phosphoprotein (SCPP) family evolved in conjunction with major evolutionary milestones: the formation of a calcified skeleton in vertebrates, the emergence of tooth enamel in fish, and the introduction of lactation in mammals. The SCPP gene family also contains genes expressed primarily and abundantly in human saliva. Here, we explored the evolution of the saliva-related SCPP genes by harnessing currently available genomic and transcriptomic resources. Our findings provide insights into the expansion and diversification of SCPP genes, notably identifying previously undocumented convergent gene duplications. In primate genomes, we found additional duplication and diversification events that affected genes coding for proteins secreted in saliva. These saliva-related SCPP genes exhibit signatures of positive selection in the primate lineage while the other genes in the same locus remain conserved. We found that regulatory shifts and gene turnover events facilitated the accelerated gain of salivary expression. Collectively, our results position the SCPP gene family as a hotbed of evolutionary innovation, suggesting the potential role of dietary and pathogenic pressures in the adaptive diversification of the saliva composition in primates, including humans.

Ancient AMY1 gene duplications primed the amylase locus for adaptive evolution upon the onset of agriculture.

Starch digestion is a cornerstone of human nutrition. The amylase enzyme, which digests starch, plays a key role in starch metabolism. Indeed, the copy number of the human amylase gene has been associated with metabolic diseases and adaptation to agricultural diets. Previous studies suggested that duplications of the salivary amylase gene are of recent origin. In the course of characterizing 51 distinct amylase haplotypes across 98 individuals employing long-read DNA sequencing and optical mapping methods, we detected four 31mers linked to duplication of the amylase locus. Analyses with these 31mers suggest that the first duplication of the amylase locus occurred more than 700,000 years ago before the split between modern humans and Neanderthals. After the original duplication events, amplification of the AMY1 genes likely occurred via nonallelic homologous recombination in a manner that consistently results in an odd number of copies per chromosome. These findings suggest that amylase haplotypes may have been primed for bursts of natural-selection associated duplications that coincided with the incorporation of starch into human diets.

Gene communities in co-expression networks across different tissues.

With the recent availability of tissue-specific gene expression data, e.g., provided by the GTEx Consortium, there is interest in comparing gene co-expression patterns across tissues. One promising approach to this problem is to use a multilayer network analysis framework and perform multilayer community detection. Communities in gene co-expression networks reveal groups of genes similarly expressed across individuals, potentially involved in related biological processes responding to specific environmental stimuli or sharing common regulatory variations. We construct a multilayer network in which each of the four layers is an exocrine gland tissue-specific gene co-expression network. We develop methods for multilayer community detection with correlation matrix input and an appropriate null model. Our correlation matrix input method identifies five groups of genes that are similarly co-expressed in multiple tissues (a community that spans multiple layers, which we call a generalist community) and two groups of genes that are co-expressed in just one tissue (a community that lies primarily within just one layer, which we call a specialist community). We further found gene co-expression communities where the genes physically cluster across the genome significantly more than expected by chance (on chromosomes 1 and 11). This clustering hints at underlying regulatory elements determining similar expression patterns across individuals and cell types. We suggest that KRTAP3-1, KRTAP3-3, and KRTAP3-5 share regulatory elements in skin and pancreas. Furthermore, we find that CELA3A and CELA3B share associated expression quantitative trait loci in the pancreas. The results indicate that our multilayer community detection method for correlation matrix input extracts biologically interesting communities of genes.

Human subsistence and signatures of selection on chemosensory genes.

Chemosensation (olfaction, taste) is essential for detecting and assessing foods, such that dietary shifts elicit evolutionary changes in vertebrate chemosensory genes. The transition from hunting and gathering to agriculture dramatically altered how humans acquire food. Recent genetic and linguistic studies suggest agriculture may have precipitated olfactory degeneration. Here, we explore the effects of subsistence behaviors on olfactory (OR) and taste (TASR) receptor genes among rainforest foragers and neighboring agriculturalists in Africa and Southeast Asia. We analyze 378 functional OR and 26 functional TASR genes in 133 individuals across populations in Uganda (Twa, Sua, BaKiga) and the Philippines (Agta, Mamanwa, Manobo) with differing subsistence histories. We find no evidence of relaxed selection on chemosensory genes in agricultural populations. However, we identify subsistence-related signatures of local adaptation on chemosensory genes within each geographic region. Our results highlight the importance of culture, subsistence economy, and drift in human chemosensory perception.

Gene communities in co-expression networks across different tissues.

With the recent availability of tissue-specific gene expression data, e.g., provided by the GTEx Consortium, there is interest in comparing gene co-expression patterns across tissues. One promising approach to this problem is to use a multilayer network analysis framework and perform multilayer community detection. Communities in gene co-expression networks reveal groups of genes similarly expressed across individuals, potentially involved in related biological processes responding to specific environmental stimuli or sharing common regulatory variations. We construct a multilayer network in which each of the four layers is an exocrine gland tissue-specific gene co-expression network. We develop methods for multilayer community detection with correlation matrix input and an appropriate null model. Our correlation matrix input method identifies five groups of genes that are similarly co-expressed in multiple tissues (a community that spans multiple layers, which we call a generalist community) and two groups of genes that are co-expressed in just one tissue (a community that lies primarily within just one layer, which we call a specialist community). We further found gene co-expression communities where the genes physically cluster across the genome significantly more than expected by chance (on chromosomes 1 and 11). This clustering hints at underlying regulatory elements determining similar expression patterns across individuals and cell types. We suggest that KRTAP3-1, KRTAP3-3, and KRTAP3-5 share regulatory elements in skin and pancreas. Furthermore, we find that CELA3A and CELA3B share associated expression quantitative trait loci in the pancreas. The results indicate that our multilayer community detection method for correlation matrix input extracts biologically interesting communities of genes.

A paleogenome from a Holocene individual supports genetic continuity in Southeast Alaska.

Many specifics of the population histories of the Indigenous peoples of North America remain contentious owing to a dearth of physical evidence. Only few ancient human genomes have been recovered from the Pacific Northwest Coast, a region increasingly supported as a coastal migration route for the initial peopling of the Americas. Here, we report paleogenomic data from the remains of a ∼3,000-year-old female individual from Southeast Alaska, named Tatóok yík yées sháawat (TYYS). Our results demonstrate at least 3,000 years of matrilineal genetic continuity in Southeast Alaska, and that TYYS is most closely related to ancient and present-day northern Pacific Northwest Coast Indigenous Americans. We find no evidence of Paleo-Inuit (represented by Saqqaq) ancestry in present-day or ancient Pacific Northwest peoples. Instead, our analyses suggest the Saqqaq genome harbors Northern Native American ancestry. This study sheds further light on the human population history of the northern Pacific Northwest Coast.

SREBP activation contributes to fatty acid accumulations in necroptosis.

Necroptosis is a type of programmed cell death. It is characterized by membrane permeabilization and is associated with the release of intracellular components due to compromised membrane integrity which induces a strong inflammatory response. We recently showed that the accumulation of very long chain fatty acids (VLCFAs) contributes to membrane permeabilization during necroptosis. However, the mechanisms that result in the accumulation of these cytotoxic lipids remain unknown. Using comparative transcriptomics and digital PCR validations, we found that several target genes of sterol regulatory element-binding proteins (SREBPs) were upregulated during necroptosis, suggesting that they might be responsible for the accumulation of VLCFA in this process. We demonstrated that activation of SREBPs during necroptosis exacerbates the permeability of the plasma membrane and cell death. Consistent with these observations, targeting sterol regulatory element-binding protein cleavage-activating protein (SCAP), a protein involved in SREBP activation, reversed the accumulation of VLCFAs, and restored cell death and membrane permeabilization during necroptosis. Collectively, our results highlight a role for SREBP in regulating lipid changes during necroptosis and suggest SREBP-mediated lipid remodeling as a potential target for therapeutics to reduce membrane permeabilization during necroptosis.

Amniotes co-opt intrinsic genetic instability to protect germ-line genome integrity.

Unlike PIWI-interacting RNA (piRNA) in other species that mostly target transposable elements (TEs), >80% of piRNAs in adult mammalian testes lack obvious targets. However, mammalian piRNA sequences and piRNA-producing loci evolve more rapidly than the rest of the genome for unknown reasons. Here, through comparative studies of chickens, ducks, mice, and humans, as well as long-read nanopore sequencing on diverse chicken breeds, we find that piRNA loci across amniotes experience: (1) a high local mutation rate of structural variations (SVs, mutations ≥ 50 bp in size); (2) positive selection to suppress young and actively mobilizing TEs commencing at the pachytene stage of meiosis during germ cell development; and (3) negative selection to purge deleterious SV hotspots. Our results indicate that genetic instability at pachytene piRNA loci, while producing certain pathogenic SVs, also protects genome integrity against TE mobilization by driving the formation of rapid-evolving piRNA sequences.

Balancing selection on genomic deletion polymorphisms in humans.

A key question in biology is why genomic variation persists in a population for extended periods. Recent studies have identified examples of genomic deletions that have remained polymorphic in the human lineage for hundreds of millennia, ostensibly owing to balancing selection. Nevertheless, genome-wide investigation of ancient and possibly adaptive deletions remains an imperative exercise. Here, we demonstrate an excess of polymorphisms in present-day humans that predate the modern human-Neanderthal split (ancient polymorphisms), which cannot be explained solely by selectively neutral scenarios. We analyze the adaptive mechanisms that underlie this excess in deletion polymorphisms. Using a previously published measure of balancing selection, we show that this excess of ancient deletions is largely owing to balancing selection. Based on the absence of signatures of overdominance, we conclude that it is a rare mode of balancing selection among ancient deletions. Instead, more complex scenarios involving spatially and temporally variable selective pressures are likely more common mechanisms. Our results suggest that balancing selection resulted in ancient deletions harboring disproportionately more exonic variants with GWAS (genome-wide association studies) associations. We further found that ancient deletions are significantly enriched for traits related to metabolism and immunity. As a by-product of our analysis, we show that deletions are, on average, more deleterious than single nucleotide variants. We can now argue that not only is a vast majority of common variants shared among human populations, but a considerable portion of biologically relevant variants has been segregating among our ancestors for hundreds of thousands, if not millions, of years.

The persistence of versions of genes that cause severe disease in human populations has long perplexed scientists. It is common for many versions of a gene to exist. But scientists expect that over time natural selection will eliminate versions of genes harmful to human health. Sometimes, there are good reasons that a disease-causing gene may persist. For example, having two copies of a particular gene variant causes a condition, called sickle cell disease. But having one sickle cell-causing copy of the gene and one non-disease-causing copy protects against malaria. As a result, the version of the gene that causes sickle cell is more common in people from areas where malaria is prevalent despite the risks to people who end up with two copies. Scientists call this phenomenon balancing selection because trade-offs in the gene's benefits and risks cause it to persist in the population. Aqil et al. show that balancing selection has likely caused many ancient gene variants to persist in human populations. In the experiments, Aqil et al. scoured the genomes of hundreds of modern humans from around the world and four groups of ancient human ancestors, including Neanderthals and Denisovans. The experiments looked for structural changes in genes, like deletions, that date back to more than 700,000 years ago ­ before modern humans split from their ancestors. They found large numbers of such ancient genes in modern humans. Using computer modeling, Aqil et al. showed that these ancient genes likely persist because of balancing selection. Many of these ancient genes regulate the immune response and metabolism. These genes may protect against infectious diseases outbreaks and starvation, which have occurred periodically throughout human history. But these same genes may cause immune or metabolic diseases in modern humans not currently facing these threats. The experiments show how such biological trade-offs have shaped human evolution and reveal that modern human populations, regardless of race or region of origin, share the same genetic variation that already our ancestors carried within them.

Genomic evidence for ancient human migration routes along South America's Atlantic coast.

An increasing body of archaeological and genomic evidence has hinted at a complex settlement process of the Americas by humans. This is especially true for South America, where unexpected ancestral signals have raised perplexing scenarios for the early migrations into different regions of the continent. Here, we present ancient human genomes from the archaeologically rich Northeast Brazil and compare them to ancient and present-day genomic data. We find a distinct relationship between ancient genomes from Northeast Brazil, Lagoa Santa, Uruguay and Panama, representing evidence for ancient migration routes along South America's Atlantic coast. To further add to the existing complexity, we also detect greater Denisovan than Neanderthal ancestry in ancient Uruguay and Panama individuals. Moreover, we find a strong Australasian signal in an ancient genome from Panama. This work sheds light on the deep demographic history of eastern South America and presents a starting point for future fine-scale investigations on the regional level.

An evolutionary trade-off between host immunity and metabolism drives fatty liver in male mice.

Adaptations to infectious and dietary pressures shape mammalian physiology and disease risk. How such adaptations affect sex-biased diseases remains insufficiently studied. In this study, we show that sex-dependent hepatic gene programs confer a robust (~300%) survival advantage for male mice during lethal bacterial infection. The transcription factor B cell lymphoma 6 (BCL6), which masculinizes hepatic gene expression at puberty, is essential for this advantage. However, protection by BCL6 protein comes at a cost during conditions of dietary excess, which result in overt fatty liver and glucose intolerance in males. Deleting hepatic BCL6 reverses these phenotypes but markedly lowers male survival during infection, thus establishing a sex-dependent trade-off between host defense and metabolic systems. Our findings offer strong evidence that some current sex-biased diseases are rooted in ancient evolutionary trade-offs between immunity and metabolism.

A mechanism of gene evolution generating mucin function.

How novel gene functions evolve is a fundamental question in biology. Mucin proteins, a functionally but not evolutionarily defined group of proteins, allow the study of convergent evolution of gene function. By analyzing the genomic variation of mucins across a wide range of mammalian genomes, we propose that exonic repeats and their copy number variation contribute substantially to the de novo evolution of new gene functions. By integrating bioinformatic, phylogenetic, proteomic, and immunohistochemical approaches, we identified 15 undescribed instances of evolutionary convergence, where novel mucins originated by gaining densely O-glycosylated exonic repeat domains. Our results suggest that secreted proteins rich in proline are natural precursors for acquiring mucin function. Our findings have broad implications for understanding the role of exonic repeats in the parallel evolution of new gene functions, especially those involving protein glycosylation.

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants.

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.

Sex-specific phenotypic effects and evolutionary history of an ancient polymorphic deletion of the human growth hormone receptor.

The common deletion of the third exon of the growth hormone receptor gene (GHRd3) in humans is associated with birth weight, growth after birth, and time of puberty. However, its evolutionary history and the molecular mechanisms through which it affects phenotypes remain unresolved. We present evidence that this deletion was nearly fixed in the ancestral population of anatomically modern humans and Neanderthals but underwent a recent adaptive reduction in frequency in East Asia. We documented that GHRd3 is associated with protection from severe malnutrition. Using a novel mouse model, we found that, under calorie restriction, Ghrd3 leads to the female-like gene expression in male livers and the disappearance of sexual dimorphism in weight. The sex- and diet-dependent effects of GHRd3 in our mouse model are consistent with a model in which the allele frequency of GHRd3 varies throughout human evolution as a response to fluctuations in resource availability.

A Single-Organelle Optical Omics Platform for Cell Science and Biomarker Discovery.

Research in fundamental cell biology and pathology could be revolutionized by developing the capacity for quantitative molecular analysis of subcellular structures. To that end, we introduce the Ramanomics platform, based on confocal Raman microspectrometry coupled to a biomolecular component analysis algorithm, which together enable us to molecularly profile single organelles in a live-cell environment. This emerging omics approach categorizes the entire molecular makeup of a sample into about a dozen of general classes and subclasses of biomolecules and quantifies their amounts in submicrometer volumes. A major contribution of our study is an attempt to bridge Raman spectrometry with big-data analysis in order to identify complex patterns of biomolecules in a single cellular organelle and leverage discovery of disease biomarkers. Our data reveal significant variations in organellar composition between different cell lines. We also demonstrate the merits of Ramanomics for identifying diseased cells by using prostate cancer as an example. We report large-scale molecular transformations in the mitochondria, Golgi apparatus, and endoplasmic reticulum that accompany the development of prostate cancer. Based on these findings, we propose that Ramanomics datasets in distinct organelles constitute signatures of cellular metabolism in healthy and diseased states.

Protein acylation by saturated very long chain fatty acids and endocytosis are involved in necroptosis.

Necroptosis is a form of cell death characterized by receptor-interacting protein kinase activity and plasma membrane permeabilization via mixed-lineage kinase-like protein (MLKL). This permeabilization is responsible for the inflammatory properties of necroptosis. We previously showed that very long chain fatty acids (VLCFAs) are functionally involved in necroptosis, potentially through protein fatty acylation. Here, we define the scope of protein acylation by saturated VLCFAs during necroptosis. We show that MLKL and phosphoMLKL, key for membrane permeabilization, are exclusively acylated during necroptosis. Reducing the levels of VLCFAs decreases their membrane recruitment, suggesting that acylation by VLCFAs contributes to their membrane localization. Acylation of phosphoMLKL occurs downstream of phosphorylation and oligomerization and appears to be, in part, mediated by ZDHHC5 (a palmitoyl transferase). We also show that disruption of endosomal trafficking increases cell viability during necroptosis, possibly by preventing recruitment, or removal, of phosphoMLKL from the plasma membrane.

Evolutionary context of psoriatic immune skin response.

The skin is vital for protecting the body and perceiving external stimuli in the environment. Ability to adapt between environments is in part based on skin phenotypic plasticity, indicating evolved homeostasis between skin and environment. This homeostasis reflects the greater relationship between the body and the environment, and disruptions in this balance may lead to accumulation of susceptibility factors for autoimmune conditions like psoriasis. In this study, we examined the relationship between rapid, lineage-specific evolution of human skin and formation of psoriatic skin responses at the transcriptome level. We collected skin tissue biopsies from individuals with psoriasis and compared gene expression in psoriatic plaques to non-plaque psoriatic skin. We then compared these data with non-psoriatic skin transcriptome data from multiple primate species. We found 67 genes showing human-specific skin expression that are also differentially regulated in psoriatic skin; these genes are significantly enriched for skin barrier function, immunity and neuronal development. We identified six gene clusters with differential expression in the context of human evolution and psoriasis, suggesting underlying regulatory mechanisms in these loci. Human and psoriasis-specific enrichment of neuroimmune genes shows the importance of the ongoing evolved homeostatic relationship between skin and external environment. These results have implications for both evolutionary medicine and public health, using transcriptomic data to acknowledge the importance of an individual's surroundings on their overall health.

Functional Specialization of Human Salivary Glands and Origins of Proteins Intrinsic to Human Saliva.

Salivary proteins are essential for maintaining health in the oral cavity and proximal digestive tract, and they serve as potential diagnostic markers for monitoring human health and disease. However, their precise organ origins remain unclear. Through transcriptomic analysis of major adult and fetal salivary glands and integration with the saliva proteome, the blood plasma proteome, and transcriptomes of 28+ organs, we link human saliva proteins to their source, identify salivary-gland-specific genes, and uncover fetal- and adult-specific gene repertoires. Our results also provide insights into the degree of gene retention during gland maturation and suggest that functional diversity among adult gland types is driven by specific dosage combinations of hundreds of transcriptional regulators rather than by a few gland-specific factors. Finally, we demonstrate the heterogeneity of the human acinar cell lineage. Our results pave the way for future investigations into glandular biology and pathology, as well as saliva's use as a diagnostic fluid.

Loss-of-function tolerance of enhancers in the human genome.

Previous studies have surveyed the potential impact of loss-of-function (LoF) variants and identified LoF-tolerant protein-coding genes. However, the tolerance of human genomes to losing enhancers has not yet been evaluated. Here we present the catalog of LoF-tolerant enhancers using structural variants from whole-genome sequences. Using a conservative approach, we estimate that individual human genomes possess at least 28 LoF-tolerant enhancers on average. We assessed the properties of LoF-tolerant enhancers in a unified regulatory network constructed by integrating tissue-specific enhancers and gene-gene interactions. We find that LoF-tolerant enhancers tend to be more tissue-specific and regulate fewer and more dispensable genes relative to other enhancers. They are enriched in immune-related cells while enhancers with low LoF-tolerance are enriched in kidney and brain/neuronal stem cells. We developed a supervised learning approach to predict the LoF-tolerance of all enhancers, which achieved an area under the receiver operating characteristics curve (AUROC) of 98%. We predict 3,519 more enhancers would be likely tolerant to LoF and 129 enhancers that would have low LoF-tolerance. Our predictions are supported by a known set of disease enhancers and novel deletions from PacBio sequencing. The LoF-tolerance scores provided here will serve as an important reference for disease studies.

Analysis of Haplotypic Variation and Deletion Polymorphisms Point to Multiple Archaic Introgression Events, Including from Altai Neanderthal Lineage.

The time, extent, and genomic effect of the introgressions from archaic humans into ancestors of extant human populations remain some of the most exciting venues of population genetics research in the past decade. Several studies have shown population-specific signatures of introgression events from Neanderthals, Denisovans, and potentially other unknown hominin populations in different human groups. Moreover, it was shown that these introgression events may have contributed to phenotypic variation in extant humans, with biomedical and evolutionary consequences. In this study, we present a comprehensive analysis of the unusually divergent haplotypes in the Eurasian genomes and show that they can be traced back to multiple introgression events. In parallel, we document hundreds of deletion polymorphisms shared with Neanderthals. A locus-specific analysis of one such shared deletion suggests the existence of a direct introgression event from the Altai Neanderthal lineage into the ancestors of extant East Asian populations. Overall, our study is in agreement with the emergent notion that various Neanderthal populations contributed to extant human genetic variation in a population-specific manner.