Convergent adaptation of human lactase persistence in Africa and Europe

AncestryGrapher toolkit: Python command-line pipelines to visualize global- and local- ancestry inferences from the RFMIX version 2 software

SUMMARY

Admixture is a fundamental process that has shaped levels and patterns of genetic variation in human populations. RFMIX version 2 (RFMIX2) utilizes a robust modeling approach to identify the genetic ancestries in admixed populations. However, this software does not have a built-in method to visually summarize the results of analyses. Here, we introduce the AncestryGrapher toolkit, which converts the numerical output of RFMIX2 into graphical representations of global and local ancestry (i.e. the per-individual ancestry components and the genetic ancestry along chromosomes, respectively).

RESULTS

To demonstrate the utility of our methods, we applied the AncestryGrapher toolkit to visualize the global and local ancestry of individuals in the North African Mozabite Berber population from the Human Genome Diversity Panel. Our results showed that the Mozabite Berbers derived their ancestry from the Middle East, Europe, and sub-Saharan Africa (global ancestry). We also found that the population origin of ancestry varied considerably along chromosomes (local ancestry). For example, we observed variance in local ancestry in the genomic region on Chromosome 2 containing the regulatory sequence in the MCM6 gene associated with lactase persistence, a human trait tied to the cultural development of adult milk consumption. Overall, the AncestryGrapher toolkit facilitates the exploration, interpretation, and reporting of ancestry patterns in human populations.

AVAILABILITY AND IMPLEMENTATION

The AncestryGrapher toolkit is free and open source on https://github.com/alisi1989/RFmix2-Pipeline-to-plot.

Deep estimation of the intensity and timing of natural selection from ancient genomes

Leveraging past allele frequencies has proven to be key for identifying the impact of natural selection across time. However, this approach suffers from imprecise estimations of the intensity (s) and timing (T) of selection, particularly when ancient samples are scarce in specific epochs. Here, we aimed to bypass the computation of allele frequencies across arbitrarily defined past epochs and refine the estimations of selection parameters by implementing convolutional neural networks (CNNs) algorithms that directly use ancient genotypes sampled across time. Using computer simulations, we first show that genotype-based CNNs consistently outperform an approximate Bayesian computation (ABC) approach based on past allele frequency trajectories, regardless of the selection model assumed and the number of available ancient genotypes. When applying this method to empirical data from modern and ancient Europeans, we replicated the reported increased number of selection events in post-Neolithic Europe, independently of the continental subregion studied. Furthermore, we substantially refined the ABC-based estimations of s and T for a set of positively and negatively selected variants, including iconic cases of positive selection and experimentally validated disease-risk variants. Our CNN predictions support a history of recent positive and negative selection targeting variants associated with host defence against pathogens, aligning with previous work that highlights the significant impact of infectious diseases, such as tuberculosis, in Europe. These findings collectively demonstrate that detecting the footprints of natural selection on ancient genomes is crucial for unravelling the history of severe human diseases.

Genetic aspects of lactase deficiency in indigenous populations of Siberia

The ability to metabolize lactose in adulthood is associated with the persistence of lactase enzyme activity. In European populations, lactase persistence is determined mainly by the presence of the rs4988235-T variant in the MCM6 gene, which increases the expression of the LCT gene, encoding lactase. The highest rates of lactase persistence are characteristic of Europeans, and the lowest rates are found in East Asian populations. Analysis of published data on the distribution of the hypolactasia-associated variant rs4988235-C in the populations of Central Asia and Siberia showed that the frequency of this variant increases in the northeastern direction. The frequency of this allele is 87 % in Central Asia, 90.6 % in Southern Siberia, and 92.9 % in Northeastern Siberia. Consequently, the ability of the population to metabolize lactose decreases in the same geographical direction. The analysis of paleogenomic data has shown that the higher frequency of the rs4988235-T allele in populations of Central Asia and Southern Siberia is associated with the eastward spread of ancient populations of the Eastern European steppes, starting from the Bronze Age. The results of polymorphism analysis of exons and adjacent introns of the MCM6 and LCT genes in indigenous populations of Siberia indicate the possibility that polymorphic variants may potentially be related to lactose metabolism exist in East Asian populations. In East Asian populations, including Siberian ethnic groups, a ~26.5 thousand nucleotide pairs long region of the MCM6 gene, including a combination of the rs4988285-A, rs2070069-G, rs3087353-T, and rs2070068-A alleles, was found. The rs4988285 and rs2070069 loci are located in the enhancer region that regulates the activity of the LCT gene. Analysis of paleogenomic sequences showed that the genomes of Denisovans and Neanderthals are characterized by the above combination of alleles of the MCM6 gene. Thus, the haplotype discovered appears to be archaic. It could have been inherited from a common ancestor of modern humans, Neanderthals, and Denisovans, or it could have been acquired by hybridization with Denisovans or Neanderthals. The data obtained indicate a possible functional significance of archaic variants of the MCM6 gene.

Recurrent evolution and selection shape structural diversity at the amylase locus

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

Recent Research on the Human Biology of Pastoralists

Despite encroachment by agricultural systems and globalization, pastoral nomads maintain a robust presence in terms of numbers and subsistence activity. At the same time, increasing concern about climate change has promoted awareness that increased climatic fluctuation may push pastoral population past their capacity for resilience. The response of pastoralists to climate change has important implications for our evolutionary past and our increasingly problematic future. Yet, pastoralists have received less explicit attention than foragers as populations under consistent selective constraints including limited caloric intake, high levels of habitual activity, and high disease burdens. Additional factors include exposure to cold and high temperatures, as well as high altitude. Over the last 20 or so years, the use of new techniques for measuring energetics, including actigraphs and doubly labeled water have built on existing noninvasive sample collection for hormones, immune markers and genes to provide a more detailed picture of the human biology of pastoral populations. Here I consider recent work on pastoralists from Siberia and northern Europe, Africa, Asia, and South America. I survey what is known about maternal milk composition and infant health, childhood growth, lactase persistence, and adult energy expenditure and lactase persistence to build a picture of the pastoralist biological response to environmental conditions, including heat, cold, and high altitude. Where available I include information about population history because of its importance for selection. I end by outlining the impact of milk consumption and climate over the human life cycle and make suggestions for further research.

Bone Mineral Density and the Risk of Type-2 Diabetes in Postmenopausal Women: rs4988235 Polymorphism Associated with Lactose Intolerance Effects

Dairy products, a major source of calcium, demonstrate a number of beneficial effects, not only protecting against the development of osteoporosis (OP) but also suppressing the onset of type-2 diabetes (T2DM) and improving bone mineral density (BMD). Dairy consumption is closely linked to lactose tolerance. One of the genetic factors predisposing individuals to lactose intolerance is rs4988235 polymorphism of the MCM6 gene. The aim of this reported study was to analyse the relationship between the rs4988235 variant of the MCM6 gene and bone mineral density and the risk of type-2 diabetes in women after menopause.

METHODS

The study was conducted among 607 female patients in the postmenopausal period in whom bone densitometry and vitamin-D3 levels were assayed and genotyping of the rs4988235 polymorphism of MCM6 gene was performed. The obtained results were analysed for the presence of T2DM, obesity surrogates, medical data, and past medical history.

RESULTS

The distribution of genotype frequencies was consistent with the Hardy-Weinberg equilibrium (p > 0.050). Postmenopausal women with the GG homozygote of rs4988235 polymorphism consumed significantly less calcium (dairy), which was probably related to the observed lactose intolerance. The GG homozygote of women with rs4988235 polymorphism was significantly more likely to have T2DM relative to the A allele carriers (p = 0.023). GG homozygotes had significantly lower femoral-vertebral mineral density despite the significantly more frequent supplementation with calcium preparations (p = 0.010), vitamin D (p = 0.01), and anti-osteoporotic drugs (p = 0.040). The obtained results indicate a stronger loss of femoral-neck mineral density with age in the GG homozygotes relative to the A allele carriers (p = 0.038).

CONCLUSIONS

In the population of women after menopause, the carriage of the G allele of rs4988235 polymorphism of the MCM6 gene, i.e., among the patients with lactose intolerance, significantly increased the risk of developing T2DM and the loss of BMD.

The molecular basis of lactase persistence: Linking genetics and epigenetics

Lactase persistence (LP) - the genetic trait that determines the continued expression of the enzyme lactase into adulthood - has undergone recent, rapid positive selection since the advent of animal domestication and dairying in some human populations. While underlying evolutionary explanations have been widely posited and studied, the molecular basis of LP remains less so. This review considers the genetic and epigenetic bases of LP. Multiple single-nucleotide polymorphisms (SNPs) in an LCT enhancer in intron 13 of the neighbouring MCM6 gene are associated with LP. These SNPs alter binding of transcription factors (TFs) and likely prevent age-related increases in methylation in the enhancer, maintaining LCT expression into adulthood to cause LP. However, the complex relationship between the genetics and epigenetics of LP is not fully characterised, and the mode of action of methylation quantitative trait loci (meQTLs) (SNPs affecting methylation) generally remains poorly understood. Here, we examine published LP data to propose a model describing how methylation in the LCT enhancer is prevented in LP adults. We argue that this occurs through altered binding of the TF Oct-1 (encoded by the gene POU2F1) and neighbouring TFs GATA-6 (GATA6), HNF-3A (FOXA1) and c-Ets1 (ETS1) acting in concert. We therefore suggest a plausible new model for LCT downregulation in the context of LP, with wider relevance for future work on the mechanisms of other meQTLs.

Food-Intolerance Genetic Testing: A Useful Tool for the Dietary Management of Chronic Gastrointestinal Disorders

The rise in food intolerances and celiac disease, along with advanced diagnostic techniques, has prompted health professionals to seek effective and economical testing methods. This study evaluates combining genetic tests with routine carbohydrate-absorption breath tests to classify patients with chronic gastrointestinal disorders into therapeutic groups, enhancing dietary management and improving gut health and quality of life. Forty-nine patients with suspected carbohydrate intolerance underwent genetic testing for lactase non-persistence, hereditary fructose intolerance, and celiac disease risk. Simultaneously, breath tests assessed lactose and fructose absorption. The lactase non-persistence genotype appeared in 36.7% of cases, with one hereditary fructose-intolerance case in a heterozygous condition. Celiac disease risk markers (HLA-DQ2/8 haplotypes) were found in 49.0% of the population. Secondary lactose and/or fructose malabsorption was present in 67.3% of patients, with 66.1% of lactase non-persistence individuals showing secondary lactose malabsorption. Fructose malabsorption was prevalent in 45.8% of patients at risk for celiac disease. Two main treatment groups were defined based on genetic results, indicating primary and irreversible gastrointestinal disorder causes, followed by a sub-classification using breath test results. Genetic testing is a valuable tool for designing dietary management plans, avoiding unnecessary diet restrictions, and reducing recovery times.

Early survival in Atlantic salmon is associated with parental genotypes at loci linked to timing of maturation

Large effect loci often contain genes with critical developmental functions and potentially broad effects across life stages. However, their life stage-specific fitness consequences are rarely explored. In Atlantic salmon, variation in two large-effect loci, six6 and vgll3, is linked to age at maturity and several physiological and behavioral traits in early life. By genotyping the progeny of wild Atlantic salmon that were planted into natural streams with nutrient manipulations, we tested if genetic variation in these loci is associated with survival in early life. We found that higher early-life survival was linked to the genotype associated with late maturation in the vgll3, but with early maturation in the six6 locus. These effects were significant in high nutrients but not in low-nutrient streams. The differences in early survival were not explained by additive genetic effects in the offspring generation but by maternal genotypes in the six6 locus and by both parents' genotypes in the vgll3 locus. Our results suggest that indirect genetic effects of large-effect loci can be significant determinants of offspring fitness. This study demonstrates an intriguing case of how large-effect loci can exhibit complex fitness associations across life stages in the wild and indicates that predicting evolutionary dynamics is difficult.

Digesting the complex metabolic effects of diet on the host and microbiome

The past 50 years of interdisciplinary research in humans and model organisms has delivered unprecedented insights into the mechanisms through which diet affects energy balance. However, translating these results to prevent and treat obesity and its associated diseases remains challenging. Given the vast scope of this literature, we focus this Review on recent conceptual advances in molecular nutrition targeting the management of energy balance, including emerging dietary and pharmaceutical interventions and their interactions with the human gut microbiome. Notably, multiple current dietary patterns of interest embrace moderate-to-high fat intake or prioritize the timing of eating over macronutrient intake. Furthermore, the rapid expansion of microbiome research findings has complicated multiple longstanding tenets of nutrition while also providing new opportunities for intervention. Continued progress promises more precise and reliable dietary recommendations that leverage our growing knowledge of the microbiome, the changing landscape of clinical interventions, and our molecular understanding of human biology.

Lactase persistence in the Jordanian population: Potential effects of the Arabian Peninsula and Sahara's aridification

The single nucleotide polymorphism (SNP) -13910 C > T has proved a good predictor of the incidence of lactase persistence in Europe and South Asia. Yet, this is not the case in the Near East, although this region is a passageway between the two continents. Lactase persistence is associated with cattle breeding, which originated in the Fertile Crescent of the Near East and spread later during the Middle Neolithic throughout Europe. Here we analyzed five SNPs (-13915 T > G (rs41380347), -13910 C > T (rs4988235), -13907 C > G (rs41525747), -14009 T > G (rs869051967), and -14010 G > C (rs145946881)) in three Jordanian human groups, namely the Bedouins, Jordan valley farmers, and Jordanian urban people. The SNPs -14009 T > G and -14010 G > C were not detected in the sample, -13907 C > G was virtually non-existent, -13910 C > T showed low frequencies, and -13915 T > G exhibited salient frequencies. The estimated incidence of lactase persistence was lower in the urban population (16 %), intermediate in the Jordan Valley's farmer population (30 %), and higher among the Bedouins (62 %). In explaining our findings, we postulated climatic change brought about by the aridification episode of the Arabian Peninsula and the Sahara 4200 years ago. This climatic milestone caused the collapse of the Akkadian Empire and the Old Kingdom in Egypt. Also, it could have led to a drastic decline of cattle in the region, being replaced by the domestication of camels. Loss of traditional crops and increasing dependence on camel milk might have triggered local selective pressures, mainly associated with -13915 T > G and differentiated from the ones in Europe, associated with -13910 C > T.

Unveiling novel genetic variants in 370 challenging medically relevant genes using the long read sequencing data of 41 samples from 19 global populations

BACKGROUND

A large number of challenging medically relevant genes (CMRGs) are situated in complex or highly repetitive regions of the human genome, hindering comprehensive characterization of genetic variants using next-generation sequencing technologies. In this study, we employed long-read sequencing technology, extensively utilized in studying complex genomic regions, to characterize genetic alterations, including short variants (single nucleotide variants and short insertions and deletions) and copy number variations, in 370 CMRGs across 41 individuals from 19 global populations.

RESULTS

Our analysis revealed high levels of genetic variants in CMRGs, with 68.73% exhibiting copy number variations and 65.20% containing short variants that may disrupt protein function across individuals. Such variants can influence pharmacogenomics, genetic disease susceptibility, and other clinical outcomes. We observed significant differences in CMRG variation across populations, with individuals of African ancestry harboring the highest number of copy number variants and short variants compared to samples from other continents. Notably, 15.79% to 33.96% of short variants were exclusively detectable through long-read sequencing. While the T2T-CHM13 reference genome significantly improved the assembly of CMRG regions, thereby facilitating variant detection in these regions, some regions still lacked resolution.

CONCLUSION

Our results provide an important reference for future clinical and pharmacogenetic studies, highlighting the need for a comprehensive representation of global genetic diversity in the reference genome and improved variant calling techniques to fully resolve medically relevant genes.

Investigating food production-associated DNA methylation changes in paleogenomes: Lack of consistent signals beyond technical noise

The Neolithic transition introduced major diet and lifestyle changes to human populations across continents. Beyond well-documented bioarcheological and genetic effects, whether these changes also had molecular-level epigenetic repercussions in past human populations has been an open question. In fact, methylation signatures can be inferred from UDG-treated ancient DNA through postmortem damage patterns, but with low signal-to-noise ratios; it is thus unclear whether published paleogenomes would provide the necessary resolution to discover systematic effects of lifestyle and diet shifts. To address this we compiled UDG-treated shotgun genomes of 13 pre-Neolithic hunter-gatherers (HGs) and 21 Neolithic farmers (NFs) individuals from West and North Eurasia, published by six different laboratories and with coverage c.1×-58× (median = 9×). We used epiPALEOMIX and a Monte Carlo normalization scheme to estimate methylation levels per genome. Our paleomethylome dataset showed expected genome-wide methylation patterns such as CpG island hypomethylation. However, analyzing the data using various approaches did not yield any systematic signals for subsistence type, genetic sex, or tissue effects. Comparing the HG-NF methylation differences in our dataset with methylation differences between hunter-gatherers versus farmers in modern-day Central Africa also did not yield consistent results. Meanwhile, paleomethylome profiles did cluster strongly by their laboratories of origin. Using larger data volumes, minimizing technical noise and/or using alternative protocols may be necessary for capturing subtle environment-related biological signals from paleomethylomes.

Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans

The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each exhibit multiple duplications/deletions with mutation rates >10,000-fold the SNP mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern day populations with traditionally agricultural diets. Leveraging 533 ancient human genomes we find that duplication-containing haplotypes (i.e. haplotypes with more amylase gene copies than the ancestral haplotype) have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in West Eurasians. Together, our study highlights the potential impacts of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.

Comparative transcriptome and methylome of polar bears, giant and red pandas reveal diet-driven adaptive evolution

Epigenetic regulation plays an important role in the evolution of species adaptations, yet little information is available on the epigenetic mechanisms underlying the adaptive evolution of bamboo-eating in both giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens). To investigate the potential contribution of epigenetic to the adaptive evolution of bamboo-eating in giant and red pandas, we performed hepatic comparative transcriptome and methylome analyses between bamboo-eating pandas and carnivorous polar bears (Ursus maritimus). We found that genes involved in carbohydrate, lipid, amino acid, and protein metabolism showed significant differences in methylation and expression levels between the two panda species and polar bears. Clustering analysis of gene expression revealed that giant pandas did not form a sister group with the more closely related polar bears, suggesting that the expression pattern of genes in livers of giant pandas and red pandas have evolved convergently driven by their similar diets. Compared to polar bears, some key genes involved in carbohydrate metabolism and biological oxidation and cholesterol synthesis showed hypomethylation and higher expression in giant and red pandas, while genes involved in fat digestion and absorption, fatty acid metabolism, lysine degradation, resistance to lipid peroxidation and detoxification showed hypermethylation and low expression. Our study elucidates the special nutrient utilization mechanism of giant pandas and red pandas and provides some insights into the molecular mechanism of their adaptive evolution of bamboo feeding. This has important implications for the breeding and conservation of giant pandas and red pandas.

Emerging Parameters Justifying a Revised Quality Concept for Cow Milk

Milk has become a staple food product globally. Traditionally, milk quality assessment has been primarily focused on hygiene and composition to ensure its safety for consumption and processing. However, in recent years, the concept of milk quality has expanded to encompass a broader range of factors. Consumers now also consider animal welfare, environmental impact, and the presence of additional beneficial components in milk when assessing its quality. This shifting consumer demand has led to increased attention on the overall production and sourcing practices of milk. Reflecting on this trend, this review critically explores such novel quality parameters, offering insights into how such practices meet the modern consumer's holistic expectations. The multifaceted aspects of milk quality are examined, revealing the intertwined relationship between milk safety, compositional integrity, and the additional health benefits provided by milk's bioactive properties. By embracing sustainable farming practices, dairy farmers and processors are encouraged not only to fulfill but to anticipate consumer standards for premium milk quality. This comprehensive approach to milk quality underscores the necessity of adapting dairy production to address the evolving nutritional landscape and consumption patterns.

In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene

Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.

Sexually dimorphic traits are associated with subsistence strategy in African faces from the Sahel/Savannah belt

OBJECTIVES

Previous research revealed that in some African populations, food-production strategies are associated with facial shape. Nomadic pastoralists living in the African Sahel/Savannah belt have a different facial morphology than their sedentary neighbors. We investigated whether the lifestyle associated with a subsistence pattern has an impact on sexual dimorphism in the facial structure.

METHODS

We employed several methods from geometric morphometrics and demonstrated such effect in four ethnically distinct populations that share the same geographic space.

RESULTS

We show that the facial traits which correlate with a subsistence strategy are systematically associated with levels of facial sex-typicality. In particular, we found that faces with more pronounced pastoralist features have on average more masculine facial traits and that this effect is more pronounced in men than in women.

CONCLUSIONS

In general, though, the magnitude of overall facial dimorphism does not differ between pastoralists and farmers. Pastoralists (in contrast to farmers) tend to have a more masculine facial morphology but facial differences between the sexes are in both groups the same.

Prebiotic Strategies to Manage Lactose Intolerance Symptoms

Lactose intolerance, which affects about 65-75% of the world's population, is caused by a genetic post-weaning deficiency of lactase, the enzyme required to digest the milk sugar lactose, called lactase non-persistence. Symptoms of lactose intolerance include abdominal pain, bloating and diarrhea. Genetic variations, namely lactase persistence, allow some individuals to metabolize lactose effectively post-weaning, a trait thought to be an evolutionary adaptation to dairy consumption. Although lactase non-persistence cannot be altered by diet, prebiotic strategies, including the consumption of galactooligosaccharides (GOSs) and possibly low levels of lactose itself, may shift the microbiome and mitigate symptoms of lactose consumption. This review discusses the etiology of lactose intolerance and the efficacy of prebiotic approaches like GOSs and low-dose lactose in symptom management.

Impact of evolution on lifestyle in microbiome

This chapter analyses the interaction between microbiota and humans from an evolutionary point of view. Long-term interactions between gut microbiota and host have been generated as a result of dietary choices through coevolutionary processes, where mutuality of advantage is essential. Likewise, the characteristics of the intestinal environment have made it possible to describe different intrahost evolutionary mechanisms affecting microbiota. For its part, the intestinal microbiota has been of great importance in the evolution of mammals, allowing the diversification of dietary niches, phenotypic plasticity and the selection of host phenotypes. Although the origin of the human intestinal microbial community is still not known with certainty, mother-offspring transmission plays a key role, and it seems that transmissibility between individuals in adulthood also has important implications. Finally, it should be noted that certain aspects inherent to modern lifestyle, including refined diets, antibiotic intake, exposure to air pollutants, microplastics, and stress, could negatively affect the diversity and composition of our gut microbiota. This chapter aims to combine current knowledge to provide a comprehensive view of the interaction between microbiota and humans throughout evolution.

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.

Human adaptations to diet: Biological and cultural coevolution

Modern humans evolved in Africa some 200,000 years ago, and since then, human populations have expanded and diversified to occupy a broad range of habitats and use different subsistence modes. This has resulted in different adaptations, such as differential responses to diseases and different abilities to digest or tolerate certain foods. The shift from a subsistence strategy based on hunting and gathering during the Palaeolithic to a lifestyle based on the consumption of domesticated animals and plants in the Neolithic can be considered one of the most important dietary transitions of Homo sapiens. In this text, we review four examples of gene-culture coevolution: (i) the persistence of the enzyme lactase after weaning, which allows the digestion of milk in adulthood, related to the emergence of dairy farming during the Neolithic; (ii) the population differences in alcohol susceptibility, in particular the ethanol intolerance of Asian populations due to the increased accumulation of the toxic acetaldehyde, related to the spread of rice domestication; (iii) the maintenance of gluten intolerance (celiac disease) with the subsequent reduced fitness of its sufferers, related to the emergence of agriculture and (iv) the considerable variation in the biosynthetic pathway of long-chain polyunsaturated fatty acids in native populations with extreme diets.

Allele age estimators designed for whole genome datasets show only a modest decrease in accuracy when applied to whole exome datasets

Personalized genomics in the healthcare system is becoming increasingly accessible as the costs of sequencing decreases. With the increase in number of genomes, larger numbers of rare variants are being discovered and much work is being done to identify their functional impacts in relation to disease phenotypes. One way to characterize these variants is to estimate the time the mutation entered the population. However, allele age estimators such as Relate, Genealogical Estimator of Variant Age, and time of coalescence, were developed based on the assumption that datasets include the entire genome. We examined the performance of each of these estimators on simulated exome data under a neutral constant population size model and found that each provides usable estimates of allele age from whole-exome datasets. To test the robustness of these methods, analyses were undertaken to simulate data under a population expansion model and background selection. Relate performs the best amongst all three estimators with Pearson coefficients of 0.64 and 0.68 (neutral constant and expansion population model) with a 17 percent and 15 percent drop in accuracy between whole genome and whole exome estimations. Of the three estimators, Relate is best able to parallelize to yield quick results with little resources, however even Relate is only able to scale to thousands of samples making it unable to match the hundreds of thousands of samples being currently released. While more work is needed to expand the capabilities of current methods of estimating allele age, these methods estimate the age of mutations with a modest decrease in performance.

Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

The effect of divergent and parallel selection on the genomic landscape of divergence

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Three linked variants have opposing regulatory effects on isovaleryl-CoA dehydrogenase gene expression

While genome-wide association studies (GWAS) and positive selection scans identify genomic loci driving human phenotypic diversity, functional validation is required to discover the variant(s) responsible. We dissected the IVD gene locus-which encodes the isovaleryl-CoA dehydrogenase enzyme-implicated by selection statistics, multiple GWAS, and clinical genetics as important to function and fitness. We combined luciferase assays, CRISPR/Cas9 genome-editing, massively parallel reporter assays (MPRA), and a deletion tiling MPRA strategy across regulatory loci. We identified three regulatory variants, including an indel, that may underpin GWAS signals for pulmonary fibrosis and testosterone, and that are linked on a positively selected haplotype in the Japanese population. These regulatory variants exhibit synergistic and opposing effects on IVD expression experimentally. Alleles at these variants lie on a haplotype tagged by the variant most strongly associated with IVD expression and metabolites, but with no functional evidence itself. This work demonstrates how comprehensive functional investigation and multiple technologies are needed to discover the true genetic drivers of phenotypic diversity.

Chromatin activity identifies differential gene regulation across human ancestries

BACKGROUND

Current evidence suggests that cis-regulatory elements controlling gene expression may be the predominant target of natural selection in humans and other species. Detecting selection acting on these elements is critical to understanding evolution but remains challenging because we do not know which mutations will affect gene regulation.

RESULTS

To address this, we devise an approach to search for lineage-specific selection on three critical steps in transcriptional regulation: chromatin activity, transcription factor binding, and chromosomal looping. Applying this approach to lymphoblastoid cells from 831 individuals of either European or African descent, we find strong signals of differential chromatin activity linked to gene expression differences between ancestries in numerous contexts, but no evidence of functional differences in chromosomal looping. Moreover, we show that enhancers rather than promoters display the strongest signs of selection associated with sites of differential transcription factor binding.

CONCLUSIONS

Overall, our study indicates that some cis-regulatory adaptation may be more easily detected at the level of chromatin than DNA sequence. This work provides a vast resource of genomic interaction data from diverse human populations and establishes a novel selection test that will benefit future study of regulatory evolution in humans and other species.

Natural and artificial selection of multiple alleles revealed through genomic analyses

Genome-to-phenome research in agriculture aims to improve crops through in silico predictions. Genome-wide association study (GWAS) is potent in identifying genomic loci that underlie important traits. As a statistical method, increasing the sample quantity, data quality, or diversity of the GWAS dataset positively impacts GWAS power. For more precise breeding, concrete candidate genes with exact functional variants must be discovered. Many post-GWAS methods have been developed to narrow down the associated genomic regions and, ideally, to predict candidate genes and causative mutations (CMs). Historical natural selection and breeding-related artificial selection both act to change the frequencies of different alleles of genes that control phenotypes. With higher diversity and more extensive GWAS datasets, there is an increased chance of multiple alleles with independent CMs in a single causal gene. This can be caused by the presence of samples from geographically isolated regions that arose during natural or artificial selection. This simple fact is a complicating factor in GWAS-driven discoveries. Currently, none of the existing association methods address this issue and need to identify multiple alleles and, more specifically, the actual CMs. Therefore, we developed a tool that computes a score for a combination of variant positions in a single candidate gene and, based on the highest score, identifies the best number and combination of CMs. The tool is publicly available as a Python package on GitHub, and we further created a web-based Multiple Alleles discovery (MADis) tool that supports soybean and is hosted in SoyKB (https://soykb.org/SoybeanMADisTool/). We tested and validated the algorithm and presented the utilization of MADis in a pod pigmentation L1 gene case study with multiple CMs from natural or artificial selection. Finally, we identified a candidate gene for the pod color L2 locus and predicted the existence of multiple alleles that potentially cause loss of pod pigmentation. In this work, we show how a genomic analysis can be employed to explore the natural and artificial selection of multiple alleles and, thus, improve and accelerate crop breeding in agriculture.

Human Nutrition Research in the Data Era: Results of 11 Reports on the Effects of a Multiple-Micronutrient-Intervention Study

Large datasets have been used in molecular and genetic research for decades, but only a few studies have included nutrition and lifestyle factors. Our team conducted an n-of-1 intervention with 12 vitamins and five minerals in 9- to 13-year-old Brazilian children and teens with poor healthy-eating indices. A unique feature of the experimental design was the inclusion of a replication arm. Twenty-six types of data were acquired including clinical measures, whole-genome mapping, whole-exome sequencing, and proteomic and a variety of metabolomic measurements over two years. A goal of this study was to use these diverse data sets to discover previously undetected physiological effects associated with a poor diet that include a more complete micronutrient composition. We summarize the key findings of 11 reports from this study that (i) found that LDL and total cholesterol and fasting glucose decreased in the population after the intervention but with inter-individual variation; (ii) associated a polygenic risk score that predicted baseline vitamin B12 levels; (iii) identified metabotypes linking diet intake, genetic makeup, and metabolic physiology; (iv) found multiple biomarkers for nutrient and food groups; and (v) discovered metabolites and proteins that are associated with DNA damage. This summary also highlights the limitations and lessons in analyzing diverse omic data.

Single nucleotide polymorphisms cumulating to genetic variation for fertility in crossbred (Bos taurus × Bos indicus) bull spermatozoa

Spermatozoa from high-fertile (HF) and low-fertile (LF) breeding bulls were subjected to high-throughput next-generation sequencing to identify important Single nucleotide polymorphisms (SNPs) and novel variants associated with fertility. A total of 77,038 genome-wide SNPs were identified, among which, 10,788 were novel variants. A total of 42,290 and 34,748 variants were recorded with 6115 and 4673 novel variants in in HF and LF bulls, respectively. Higher number of SNPs were identified in HF compared to LF bulls. GO analysis of filtered genes with significant variations in HF bulls indicated their involvement in oxidative phosphorylation and metabolic pathways. GO analysis of filtered genes with significant variation in LF bulls revealed their involvement in Ca2++ ion binding, structural constituent of ribosome, and biological processes like translation and ribosomal small subunit assembly. The study identified SNPs in candidate genes including TPT1, BOLA-DRA, CD74, RPS17, RPS28, RPS29, RPL14, RPL13, and RPS27A, which are linked to sperm functionality, survival, oxidative stress, and bull fertility. The identified SNPs could be used in selection of bulls for high fertility and the variation in these genes could be established as an explanation for the fertility differences in bulls upon validation in large number of bulls.

Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Comparing measured dietary variation within and between tropical hunter-gatherer groups to the Paleo Diet

BACKGROUND

Although human diets varied considerably before the spread of agriculture, public perceptions of preagricultural diets have been strongly influenced by the Paleo Diet, which prescribes percentage calorie ranges of 19-35% protein, 22-40% carbohydrate, and 28-47% fat, and prohibits foods with added sugar, dairy, grains, most starchy tubers, and legumes. However, the empirical basis for Paleolithic nutrition remains unclear, with some of its assumptions challenged by the archaeological record and theoretical first principles.

OBJECTIVES

We assessed the variation in diets among tropical hunter-gatherers, including the effect of collection methods on implied macronutrient percentages.

METHODS

We analyzed data on animal food, plant food, and honey consumption by weight and kcal from 15 high-quality published ethnographic studies representing 11 recent tropical hunter-gatherer groups. We used Bayesian analyses to perform inference and included data collection methods and environmental variables as predictors in our models.

RESULTS

Our analyses reveal high levels of variation in animal versus plant foods consumed and in corresponding percentages of protein, fat, and carbohydrates. In addition, studies that weighed food items consumed in and out of camp and across seasons and years reported higher consumption of animal foods, which varied with annual mean temperature.

CONCLUSIONS

The ethnographic evidence from tropical foragers refutes the concept of circumscribed macronutrient ranges modeling preagricultural diets.

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.

Hunter-gatherer diets and activity as a model for health promotion: Challenges, responses, and confirmations

Beginning in 1985, we and others presented estimates of hunter-gatherer (and ultimately ancestral) diet and physical activity, hoping to provide a model for health promotion. The Hunter-Gatherer Model was designed to offset the apparent mismatch between our genes and the current Western-type lifestyle, a mismatch that arguably affects prevalence of many chronic degenerative diseases. The effort has always been controversial and subject to both scientific and popular critiques. The present article (1) addresses eight such challenges, presenting for each how the model has been modified in response, or how the criticism can be rebutted; (2) reviews new epidemiological and experimental evidence (including especially randomized controlled clinical trials); and (3) shows how official recommendations put forth by governments and health authorities have converged toward the model. Such convergence suggests that evolutionary anthropology can make significant contributions to human health.

Association of Estimated Daily Lactose Consumption, Lactase Persistence Genotype (rs4988235), and Gut Microbiota in Healthy Adults in the United States

BACKGROUND

Lactase persistence (LP) is a heritable trait in which lactose can be digested throughout adulthood. Lactase nonpersistent (LNP) individuals who consume lactose may experience microbial adaptations in response to undigested lactose.

OBJECTIVES

The objective of the study was to estimate lactose from foods reported in the Automated Self-Administered 24-Hour Dietary Assessment Tool (ASA24) and determine the interaction between lactose consumption, LP genotype, and gut microbiome in an observational cross-sectional study of healthy adults in the United States (US).

METHODS

Average daily lactose consumption was estimated for 279 healthy US adults, genotyped for the lactase gene -13910G>A polymorphism (rs4988235) by matching ASA24-reported foods to foods in the Nutrition Coordinating Center Food and Nutrient Database. Analysis of covariance was used to identify whether the A genotype (LP) influenced lactose and total dairy consumption, with total energy intake and weight as covariates. The 16S rRNA V4/V5 region, amplified from bacterial DNA extracted from each frozen stool sample, was sequenced using Illumina MiSeq (300 bp paired-end) and analyzed using Quantitative Insights Into Microbial Ecology (QIIME)2 (version 2019.10). Differential abundances of bacterial taxa were analyzed using DESeq2 likelihood ratio tests.

RESULTS

Across a diverse set of ethnicities, LP subjects consumed more lactose than LNP subjects. Lactobacillaceae abundance was highest in LNP subjects who consumed more than 12.46 g/d (upper tercile). Within Caucasians and Hispanics, family Lachnospiraceae was significantly enriched in the gut microbiota of LNP individuals consuming the upper tercile of lactose across both sexes.

CONCLUSIONS

Elevated lactose consumption in individuals with the LNP genotype is associated with increased abundance of family Lactobacillaceae and Lachnospriaceae, taxa that contain multiple genera capable of utilizing lactose. This trial was registered on clinicaltrials.gov as NCT02367287.

Milk Transmission of Mammalian Retroviruses

The transmission of viruses from one host to another typically occurs through horizontal or vertical pathways. The horizontal pathways include transmission amongst individuals, usually through bodily fluids or excretions, while vertical transmission transpires from mother to their offspring, either during pregnancy, childbirth, or breastfeeding. While there are more than 200 human pathogenic viruses to date, only a small number of them are known to be transmitted via breast milk, including cytomegalovirus (CMV), human immunodeficiency virus type 1 (HIV-1), and human T cell lymphotropic virus type 1 (HTLV-1), the latter two belonging to the family Retroviridae. Breast milk transmission is a common characteristic among mammalian retroviruses, but there is a lack of reports summarizing our knowledge regarding this route of transmission of mammalian retroviruses. Here, we provide an overview of the transmission of mammalian exogenous retroviruses with a focus on Orthoretrovirinae, and we highlight whether they have been described or suspected to be transmitted through breast milk, covering various species. We also elaborate on the production and composition of breast milk and discuss potential entry sites of exogenous mammalian retroviruses during oral transmission.

Uncovering Footprints of Natural Selection Through Spectral Analysis of Genomic Summary Statistics

Natural selection leaves a spatial pattern along the genome, with a haplotype distribution distortion near the selected locus that fades with distance. Evaluating the spatial signal of a population-genetic summary statistic across the genome allows for patterns of natural selection to be distinguished from neutrality. Considering the genomic spatial distribution of multiple summary statistics is expected to aid in uncovering subtle signatures of selection. In recent years, numerous methods have been devised that consider genomic spatial distributions across summary statistics, utilizing both classical machine learning and deep learning architectures. However, better predictions may be attainable by improving the way in which features are extracted from these summary statistics. We apply wavelet transform, multitaper spectral analysis, and S-transform to summary statistic arrays to achieve this goal. Each analysis method converts one-dimensional summary statistic arrays to two-dimensional images of spectral analysis, allowing simultaneous temporal and spectral assessment. We feed these images into convolutional neural networks and consider combining models using ensemble stacking. Our modeling framework achieves high accuracy and power across a diverse set of evolutionary settings, including population size changes and test sets of varying sweep strength, softness, and timing. A scan of central European whole-genome sequences recapitulated well-established sweep candidates and predicted novel cancer-associated genes as sweeps with high support. Given that this modeling framework is also robust to missing genomic segments, we believe that it will represent a welcome addition to the population-genomic toolkit for learning about adaptive processes from genomic data.

High-Altitude Andean H194R HIF2A Allele Is a Hypomorphic Allele

For over 10,000 years, Andeans have resided at high altitude where the partial pressure of oxygen challenges human survival. Recent studies have provided evidence for positive selection acting in Andeans on the HIF2A (also known as EPAS1) locus, which encodes for a central transcription factor of the hypoxia-inducible factor pathway. However, the precise mechanism by which this allele might lead to altitude-adaptive phenotypes, if any, is unknown. By analyzing whole genome sequencing data from 46 high-coverage Peruvian Andean genomes, we confirm evidence for positive selection acting on HIF2A and a unique pattern of variation surrounding the Andean-specific single nucleotide variant (SNV), rs570553380, which encodes for an H194R amino acid substitution in HIF-2α. Genotyping the Andean-associated SNV rs570553380 in a group of 299 Peruvian Andeans from Cerro de Pasco, Peru (4,338 m), reveals a positive association with increased fraction of exhaled nitric oxide, a marker of nitric oxide biosynthesis. In vitro assays show that the H194R mutation impairs binding of HIF-2α to its heterodimeric partner, aryl hydrocarbon receptor nuclear translocator. A knockin mouse model bearing the H194R mutation in the Hif2a gene displays decreased levels of hypoxia-induced pulmonary Endothelin-1 transcripts and protection against hypoxia-induced pulmonary hypertension. We conclude the Andean H194R HIF2A allele is a hypomorphic (partial loss of function) allele.

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.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

Feedback circuits are numerous in embryonic gene regulatory networks and offer a stabilizing influence on evolution of those networks

The developmental gene regulatory networks (dGRNs) of two sea urchin species, Lytechinus variegatus (Lv) and Strongylocentrotus purpuratus (Sp), have remained remarkably similar despite about 50 million years since a common ancestor. Hundreds of parallel experimental perturbations of transcription factors with similar outcomes support this conclusion. A recent scRNA-seq analysis suggested that the earliest expression of several genes within the dGRNs differs between Lv and Sp. Here, we present a careful reanalysis of the dGRNs in these two species, paying close attention to timing of first expression. We find that initial expression of genes critical for cell fate specification occurs during several compressed time periods in both species. Previously unrecognized feedback circuits are inferred from the temporally corrected dGRNs. Although many of these feedbacks differ in location within the respective GRNs, the overall number is similar between species. We identify several prominent differences in timing of first expression for key developmental regulatory genes; comparison with a third species indicates that these heterochronies likely originated in an unbiased manner with respect to embryonic cell lineage and evolutionary branch. Together, these results suggest that interactions can evolve even within highly conserved dGRNs and that feedback circuits may buffer the effects of heterochronies in the expression of key regulatory genes.

Pilot study of correlation of selected genetic factors with cribra orbitalia in individuals from a medieval population from Slovakia

OBJECTIVE

The aim of this study is to investigate the potential genetic etiology of cribra orbitalia noted on human skeletal remains.

MATERIALS

We obtained and analyzed ancient DNA of 43 individuals with cribra orbitalia. The analyzed set represented medieval individuals from two cemeteries in western Slovakia, Castle Devín (11th-12th century AD) and Cífer-Pác (8th-9th century AD).

METHODS

We performed a sequence analysis of 5 variants in 3 genes associated with anemia (HBB, G6PD, PKLR), which are the most common pathogenic variants in present day of European populations, and one variant MCM6:c.1917 + 326 C>T (rs4988235) associated with lactose intolerance.

RESULTS

DNA variants associated with anemia were not found in the samples. The allele frequency of MCM6:c.1917 + 326 C was 0.875. This frequency is higher but not statistically significant in individuals displaying cribra orbitalia compared to individuals without the lesion.

SIGNIFICANCE

This study seeks to expand our knowledge of the etiology of cribra orbitalia by exploring the potential association between the lesion and the presence of alleles linked to hereditary anemias and lactose intolerance.

LIMITATIONS

A relatively small set of individuals were analyzed, so an unequivocal conclusion cannot be drawn. Hence, although it is unlikely, a genetic form of anemia caused by rare variants cannot be ruled out.

SUGGESTIONS FOR FURTHER RESEARCH

Genetic research based on larger sample sizes and in more diverse geographical regions.

DNA polymorphisms associated with lactase persistence, self-perceived symptoms of lactose intolerance, milk and dairy consumption, and ancestry, in the Uruguayan population

Uruguay has one of the highest per capita milk intakes worldwide, even with a limited supply of lactose-free products; furthermore, the admixed nature of its population is well known, and various frequencies of lactase persistence (LP) are observed in the source populations. We aimed to contribute to the understanding of the relation between allelic variants associated with LP, milk consumption, digestive symptoms, and genetic ancestry in the Uruguayan population. Samples of saliva or peripheral blood were collected from 190 unrelated individuals from two regions of Uruguay, genotypes for polymorphic sites in a fragment within the LCT enhancer were determined and allelic frequencies calculated in all of them. Data were collected on frequency of milk and dairy consumption and self-reported symptoms in a subsample of 153 individuals. Biparental and maternal ancestry was determined by analyzing individual ancestry markers and mitochondrial DNA. Twenty-nine percentage of individuals reported symptoms attributed to the ingestion of fresh milk, with abdominal pain, bloating and flatulence being the most frequent. European LP-associated allele T-13910 showed a frequency of 33%, while other LP-associated alleles like G-13915 and T-14011 were observed in very low frequencies. Associations between self-reported symptoms, fresh milk intake, and C/T-13910 genotype were statistically significant. No evidence of association between genetic ancestry and C/T-13910 was found, although individuals carrying one T-13910 allele appeared to have more European ancestry. In conclusion, the main polymorphism capable of predicting lactose intolerance in Uruguayans is C/T-13910, although more studies are required to unravel the relation between genotype and lactase activity, especially in heterozygotes.

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.

Gene-by-environment interactions are pervasive among natural genetic variants

Gene-by-environment (GxE) interactions, in which a genetic variant's phenotypic effect is condition specific, are fundamental for understanding fitness landscapes and evolution but have been difficult to identify at the single-nucleotide level. Although many condition-specific quantitative trait loci (QTLs) have been mapped, these typically contain numerous inconsequential variants in linkage, precluding understanding of the causal GxE variants. Here, we introduce BARcoded Cas9 retron precise parallel editing via homology (CRISPEY-BAR), a high-throughput precision genome editing strategy, and use it to map GxE interactions of naturally occurring genetic polymorphisms impacting yeast growth. We identified hundreds of GxE variants within condition-specific QTLs, revealing unexpected genetic complexity. Moreover, we found that 93.7% of non-neutral natural variants within ergosterol biosynthesis pathway genes showed GxE interactions, including many impacting antifungal drug resistance through diverse molecular mechanisms. In sum, our results suggest an extremely complex, context-dependent fitness landscape characterized by pervasive GxE interactions while also demonstrating massively parallel genome editing as an effective means for investigating this complexity.

Tensor decomposition based feature extraction and classification to detect natural selection from genomic data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under non-convex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data while preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx , which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Uncovering the gastrointestinal passage, intestinal epithelial cellular uptake and AGO2 loading of milk miRNAs in neonates using xenomiRs as tracers

BACKGROUND

Human breast milk has a high microRNA (miRNA) content. It remains unknown whether and how milk miRNAs might affect intestinal gene regulation and homeostasis of the developing microbiome after initiation of enteral nutrition. However, this requires that relevant milk miRNA amounts survive gastrointestinal passage, are taken up by cells, and become available to the RNA interference (RNAi) machinery. It seems important to dissect the fate of these miRNAs after oral ingestion and gastrointestinal passage.

OBJECTIVE

Our goal was to analyze the potential transmissibility of milk miRNAs via the gastrointestinal system in neonate humans and a porcine model in vivo to contribute to the discussion whether milk miRNAs could influence gene regulation in neonates and thus might vertically transmit developmental relevant signals.

DESIGN

We performed cross-species profiling of miRNAs via deep-sequencing and utilized dietary xenobiotic taxon-specific milk miRNA (xenomiRs) as tracers in human and porcine neonates, followed by functional studies in primary human fetal intestinal epithelial cells (HIEC-6) using Ad5-mediated miRNA-gene transfer.

RESULTS

Mammals share many milk miRNAs yet exhibit taxon-specific miRNA fingerprints. We traced bovine-specific miRNAs from formula-nutrition in human preterm stool and 9 days after onset of enteral feeding in intestinal cells of preterm piglets. Thereafter, several xenomiRs accumulated in the intestinal cells. Moreover, few hours after introducing enteral feeding in preterm piglets with supplemented reporter miRNAs (cel-miR-39-5p/-3p), we observed their enrichment in blood serum and in AGO2-immunocomplexes from intestinal biopsies.

CONCLUSIONS

Milk-derived miRNAs survived gastrointestinal passage in human and porcine neonates. Bovine-specific miRNAs accumulated in intestinal cells of preterm piglets after enteral feeding with bovine colostrum/formula. In piglets, colostrum supplementation with cel-miR-39-5p/-3p resulted in increased blood levels of cel-miR-39-3p and argonaute RISC catalytic component 2 (AGO2) loading in intestinal cells. This suggests the possibility of vertical transmission of miRNA signaling from milk through the neonatal digestive tract.

Human genetic history on the Tibetan Plateau in the past 5100 years

Using genome-wide data of 89 ancient individuals dated to 5100 to 100 years before the present (B.P.) from 29 sites across the Tibetan Plateau, we found plateau-specific ancestry across plateau populations, with substantial genetic structure indicating high differentiation before 2500 B.P. Northeastern plateau populations rapidly showed admixture associated with millet farmers by 4700 B.P. in the Gonghe Basin. High genetic similarity on the southern and southwestern plateau showed population expansion along the Yarlung Tsangpo River since 3400 years ago. Central and southeastern plateau populations revealed extensive genetic admixture within the plateau historically, with substantial ancestry related to that found in southern and southwestern plateau populations. Over the past ~700 years, substantial gene flow from lowland East Asia further shaped the genetic landscape of present-day plateau populations. The high-altitude adaptive EPAS1 allele was found in plateau populations as early as in a 5100-year-old individual and showed a sharp increase over the past 2800 years.