A positively selected FBN1 missense variant reduces height in Peruvian individuals

Recent selection and introgression facilitated high-altitude adaptation in cattle

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Three-Dimensional Morphometric Analysis of Anterior Cerebral Circulation Aneurysms

This article assesses the association between anterior circulation morphometry and the presence of intracranial aneurysm using three-dimensional rotational angiography (3DRA). A retrospective analysis at a Peruvian academic medical center between December 2018 and February 2020 identified 206 patients with unruptured intracranial aneurysms and matched controls who underwent 3DRA. Angiographic images were obtained per standard of care, and measurements of the vasculature were performed using 3DRA vascular automated software. A total of 163 aneurysms and 43 control angiograms were evaluated. Women represented 82.5% of the cases and the mean age was 55.9 years (standard deviation ± 14.2). In multivariate analysis, five specific features were found to be statistically significant predictors for presence of an anterior circulation aneurysm: female sex (odds ratio [OR] = 2.71; p  = 0.048), C-shape of the middle cerebral artery (MCA) (OR = 2.73; p  = 0.018), distal internal carotid artery (ICA) diameter (OR = 3.42; p  = 0.012), ICA bifurcation angle (OR = 1.02; p  = 0.036), and length of the carotid siphon (OR = 1.08; p  = 0.047). Features detected on 3DRA suggest morphological characteristics of the ICA and MCA may be predictive for intracranial aneurysm. Our findings build from prior reports by demonstrating five specific patient and imaging features associated with anterior circulation aneurysms. While 3DRA is the standard of care in many settings, medical centers with resource limitations may not have access to this technique. The demographic and morphological features identified in our study may have correlates that if detected on contrast computed tomography or magnetic resonance imaging studies, may be used to help screen for a higher level of care in select patients.

A Neanderthal/Denisovan GLI3 variant contributes to anatomical variations in mice

Changes in genomic structures underlie phenotypic diversification in organisms. Amino acid-changing mutations affect pleiotropic functions of proteins, although little is known about how mutated proteins are adapted in existing developmental programs. Here we investigate the biological effects of a variant of the GLI3 transcription factor (GLI3R1537C) carried in Neanderthals and Denisovans, which are extinct hominins close to modern humans. R1537C does not compromise protein stability or GLI3 activator-dependent transcriptional activities. In contrast, R1537C affects the regulation of downstream target genes associated with developmental processes. Furthermore, genome-edited mice carrying the Neanderthal/Denisovan GLI3 mutation exhibited various alterations in skeletal morphology. Our data suggest that an extinct hominin-type GLI3 contributes to species-specific anatomical variations, which were tolerated by relaxed constraint in developmental programs during human evolution.

Extremely sparse models of linkage disequilibrium in ancestrally diverse association studies

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using large correlation matrices, but this approach is inefficient, especially in ancestrally diverse studies. In the present study, we introduce LD graphical models (LDGMs), which are an extremely sparse and efficient representation of LD. LDGMs are derived from genome-wide genealogies; statistical relationships among alleles in the LDGM correspond to genealogical relationships among haplotypes. We published LDGMs and ancestry-specific LDGM precision matrices for 18 million common variants (minor allele frequency >1%) in five ancestry groups, validated their accuracy and demonstrated order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implemented an extremely fast multiancestry polygenic prediction method, BLUPx-ldgm, which performs better than a similar method based on the reference LD correlation matrix. LDGMs will enable sophisticated methods that scale to ancestrally diverse genetic association data across millions of variants and individuals.

Multi-ancestry genome-wide meta-analysis of 56,241 individuals identifies LRRC4C, LHX5-AS1 and nominates ancestry-specific loci PTPRK , GRB14 , and KIAA0825 as novel risk loci for Alzheimer's disease: the Alzheimer's Disease Genetics Consortium

Limited ancestral diversity has impaired our ability to detect risk variants more prevalent in non-European ancestry groups in genome-wide association studies (GWAS). We constructed and analyzed a multi-ancestry GWAS dataset in the Alzheimer's Disease (AD) Genetics Consortium (ADGC) to test for novel shared and ancestry-specific AD susceptibility loci and evaluate underlying genetic architecture in 37,382 non-Hispanic White (NHW), 6,728 African American, 8,899 Hispanic (HIS), and 3,232 East Asian individuals, performing within-ancestry fixed-effects meta-analysis followed by a cross-ancestry random-effects meta-analysis. We identified 13 loci with cross-ancestry associations including known loci at/near CR1 , BIN1 , TREM2 , CD2AP , PTK2B , CLU , SHARPIN , MS4A6A , PICALM , ABCA7 , APOE and two novel loci not previously reported at 11p12 ( LRRC4C ) and 12q24.13 ( LHX5-AS1 ). Reflecting the power of diverse ancestry in GWAS, we observed the SHARPIN locus using 7.1% the sample size of the original discovering single-ancestry GWAS (n=788,989). We additionally identified three GWS ancestry-specific loci at/near ( PTPRK ( P =2.4×10 -8 ) and GRB14 ( P =1.7×10 -8 ) in HIS), and KIAA0825 ( P =2.9×10 -8 in NHW). Pathway analysis implicated multiple amyloid regulation pathways (strongest with P adjusted =1.6×10 -4 ) and the classical complement pathway ( P adjusted =1.3×10 -3 ). Genes at/near our novel loci have known roles in neuronal development ( LRRC4C, LHX5-AS1 , and PTPRK ) and insulin receptor activity regulation ( GRB14 ). These findings provide compelling support for using traditionally-underrepresented populations for gene discovery, even with smaller sample sizes.

Identification of two variants in PAX3 and FBN1 in a Chinese family with Waardenburg and Marfan syndrome via whole exome sequencing

Both Warrensburg (WS) and Marfan syndrome (MFS) can impair the vision. Here, we recruited a Chinese family consisting of two WS affected individuals (II:1 and III:3) and five MFS affected individuals( I:1, II:2, III:1, III:2, and III:5) as well as one suspected MFS individual (II:4). Using whole exome sequencing (WES) and subsequent PCR-Sanger sequencing, we identified one novel heterozygous variant NM_000438 (PAX3) c.208 T > C, (p.Cys70Arg) from individuals with WS and one previous reported variant NM_000138 (FBN1) c.2740 T > A, (p.Cys914Ser) from individuals with MFS and co-segregated with the diseases. Real-time PCR and Western blot assay showed that, compared to their wild-type, both mRNAs and proteins of  PAX3 and FBN1 mutants reduced in HKE293T cells. Together, our study identified two disease-causing variants in a same Chinese family with WS and MFS, and confirmed their damaged effects on their genes' expression. Therefore, those findings expand the mutation spectrum of PAX3 and provide a new perspective for the potential therapy.

New insights from GWAS on BMI-related growth traits in a longitudinal cohort of admixed children with Native American and European ancestry

Body-mass index (BMI) is a hallmark of adiposity. In contrast with adulthood, the genetic architecture of BMI during childhood is poorly understood. The few genome-wide association studies (GWAS) on children have been performed almost exclusively in Europeans and at single ages. We performed cross-sectional and longitudinal GWAS for BMI-related traits on 904 admixed children with mostly Mapuche Native American and European ancestries. We found regulatory variants of the immune gene HLA-DQB3 strongly associated with BMI at 1.5 - 2.5 years old. A variant in the sex-determining gene DMRT1 was associated with the age at adiposity rebound (Age-AR) in girls (P = 9.8 × 10 - 9 ). BMI was significantly higher in Mapuche than in Europeans between 5.5 and 16.5 years old. Finally, Age-AR was significantly lower (P = 0.004 ) by 1.94 years and BMI at AR was significantly higher (P = 0.04 ) by 1.2 kg/ m 2 , in Mapuche children compared with Europeans.

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Molecular mechanisms underlying hematophagia revealed by comparative analyses of leech genomes

BACKGROUND

Leeches have been used in traditional Chinese medicine since prehistoric times to treat a spectrum of ailments, but very little is known about their physiological, genetic, and evolutionary characteristics.

FINDINGS

We sequenced and assembled chromosome-level genomes of 3 leech species (bloodsucking Hirudo nipponia and Hirudinaria manillensis and nonbloodsucking Whitmania pigra). The dynamic population histories and genome-wide expression patterns of the 2 bloodsucking leech species were found to be similar. A combined analysis of the genomic and transcriptional data revealed that the bloodsucking leeches have a presumably enhanced auditory sense for prey location in relatively deep fresh water. The copy number of genes related to anticoagulation, analgesia, and anti-inflammation increased in the bloodsucking leeches, and their gene expressions responded dynamically to the bloodsucking process. Furthermore, the expanded FBN1 gene family may help in rapid body swelling of leeches after bloodsucking, and the expanded GLB3 gene family may be associated with long-term storage of prey blood in a leech's body.

CONCLUSIONS

The high-quality reference genomes and comprehensive datasets obtained in this study may facilitate innovations in the artificial culture and strain optimization of leeches.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

The Association between Altitude and Waist-Height Ratio in Peruvian Adults: A Cross-Sectional Data Analysis of a Population-Based Survey

To evaluate the association between altitude and cardiometabolic risk calculated with the weight-height ratio (WHtR) in the Peruvian adult population via the cross-sectional data analysis of the Peruvian Demographic and Health Survey 2021. A total of 26,117 adults from 18 to 64 years of age were included in the analysis. The dependent variable was cardiometabolic risk, defined as "Yes" if the WHtR was ≥0.5 and "No" if the WHtR was <0.5. Exposure was altitude of residence categorized as: <1500 meters above sea level (masl); 1500 to 2499 masl; 2500 to 3499 masl; and ≥3500 masl. Crude and adjusted Poisson regression models were used to calculate prevalence ratios (PR) with 95% confidence intervals (CI). The mean WHtR in the population was 0.59 (standard deviation: 0.08), and 87.6% (95% CI: 86.9-88.2) were classified as at risk. After adjusting for sex, age, education level, well-being index, and area of residence, living at altitudes between 2500 and 3499 masl (aPR: 0.98; 95% CI: 0.96-1.00) and ≥3500 masl (aPR: 0.95; 95% CI: 0.93-0.97) were associated with lower cardiometabolic risk in comparison with living at <1500 masl. An inverse association was identified between living at a higher altitude and the proportion of cardiometabolic risk in the Peruvian adult population. However, at least 8 out of 10 people were identified as at risk in all categories of altitude.

The schizophrenia-associated missense variant rs13107325 regulates dendritic spine density

The missense variant rs13107325 (C/T, p.Ala391Thr) in SLC39A8 consistently showed robust association with schizophrenia in recent genome-wide association studies (GWASs), suggesting the potential pathogenicity of this non-synonymous risk variant. Nevertheless, how this missense variant confers schizophrenia risk remains unknown. Here we constructed a knock-in mouse model (by introducing a threonine at the 393th amino acid of mouse SLC39A8 (SLC39A8-p.393T), which corresponds to rs13107325 (p.Ala391Thr) of human SLC39A8) to explore the potential roles and biological effects of this missense variant in schizophrenia pathogenesis. We assessed multiple phenotypes and traits (associated with rs13107325) of the knock-in mice, including body and brain weight, concentrations of metal ions (including cadmium, zinc, manganese, and iron) transported by SLC39A8, blood lipids, proliferation and migration of neural stem cells (NSCs), cortical development, behaviors and cognition, transcriptome, dendritic spine density, and synaptic transmission. Many of the tested phenotypes did not show differences in SLC39A8-p.393T knock-in and wild-type mice. However, we found that zinc concentration in brain and blood of SLC39A8-p.393T knock-in mice was dysregulated compared with wild-types, validating the functionality of rs13107325. Further analysis indicated that cortical dendritic spine density of the SLC39A8-p.393T knock-in mice was significantly decreased compared with wild-types, indicating the important role of SLC39A8-p.393T in dendritic spine morphogenesis. These results indicated that SLC39A8-p.393T knock-in resulted in decreased dendritic spine density, thus mimicking the dendritic spine pathology observed in schizophrenia. Our study indicates that rs13107325 might confer schizophrenia risk by regulating zinc concentration and dendritic spine density, a featured characteristic that was frequently reported to be decreased in schizophrenia.

Importance of Including Non-European Populations in Large Human Genetic Studies to Enhance Precision Medicine

One goal of genomic medicine is to uncover an individual's genetic risk for disease, which generally requires data connecting genotype to phenotype, as done in genome-wide association studies (GWAS). While there may be clinical promise to employing prediction tools such as polygenic risk scores (PRS), it currently stands that individuals of non-European ancestry may not reap the benefits of genomic medicine because of underrepresentation in large-scale genetics studies. Here, we discuss why this inequity poses a problem for genomic medicine and the reasons for the low transferability of PRS across populations. We also survey the ancestry representation of published GWAS and investigate how estimates of ancestry diversity in GWASparticipants might be biased. We highlight the importance of expanding genetic research in Africa, one of the most underrepresented regions in human genomics research, and discuss issues of ethics, resources, and technology for equitable advancement of genomic medicine.

The sequences of 150,119 genomes in the UK Biobank

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data^1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank³. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

Using the UK Biobank as a global reference of worldwide populations: application to measuring ancestry diversity from GWAS summary statistics

MOTIVATION

Measuring genetic diversity is an important problem because increasing genetic diversity is a key to making new genetic discoveries, while also being a major source of confounding to be aware of in genetics studies.

RESULTS

Using the UK Biobank data, a prospective cohort study with deep genetic and phenotypic data collected on almost 500 000 individuals from across the UK, we carefully define 21 distinct ancestry groups from all four corners of the world. These ancestry groups can serve as a global reference of worldwide populations, with a handful of applications. Here, we develop a method that uses allele frequencies and principal components derived from these ancestry groups to effectively measure ancestry proportions from allele frequencies of any genetic dataset.

AVAILABILITY AND IMPLEMENTATION

This method is implemented in function snp_ancestry_summary of R package bigsnpr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Fully exploiting SNP arrays: a systematic review on the tools to extract underlying genomic structure

Single nucleotide polymorphisms (SNPs) are the most abundant type of genomic variation and the most accessible to genotype in large cohorts. However, they individually explain a small proportion of phenotypic differences between individuals. Ancestry, collective SNP effects, structural variants, somatic mutations or even differences in historic recombination can potentially explain a high percentage of genomic divergence. These genetic differences can be infrequent or laborious to characterize; however, many of them leave distinctive marks on the SNPs across the genome allowing their study in large population samples. Consequently, several methods have been developed over the last decade to detect and analyze different genomic structures using SNP arrays, to complement genome-wide association studies and determine the contribution of these structures to explain the phenotypic differences between individuals. We present an up-to-date collection of available bioinformatics tools that can be used to extract relevant genomic information from SNP array data including population structure and ancestry; polygenic risk scores; identity-by-descent fragments; linkage disequilibrium; heritability and structural variants such as inversions, copy number variants, genetic mosaicisms and recombination histories. From a systematic review of recently published applications of the methods, we describe the main characteristics of R packages, command-line tools and desktop applications, both free and commercial, to help make the most of a large amount of publicly available SNP data.

Populations, Traits, and Their Spatial Structure in Humans

The spatial distribution of genetic variants is jointly determined by geography, past demographic processes, natural selection, and its interplay with environmental variation. A fraction of these genetic variants are "causal alleles" that affect the manifestation of a complex trait. The effect exerted by these causal alleles on complex traits can be independent or dependent on the environment. Understanding the evolutionary processes that shape the spatial structure of causal alleles is key to comprehend the spatial distribution of complex traits. Natural selection, past population size changes, range expansions, consanguinity, assortative mating, archaic introgression, admixture, and the environment can alter the frequencies, effect sizes, and heterozygosities of causal alleles. This provides a genetic axis along which complex traits can vary. However, complex traits also vary along biogeographical and sociocultural axes which are often correlated with genetic axes in complex ways. The purpose of this review is to consider these genetic and environmental axes in concert and examine the ways they can help us decipher the variation in complex traits that is visible in humans today. This initiative necessarily implies a discussion of populations, traits, the ability to infer and interpret "genetic" components of complex traits, and how these have been impacted by adaptive events. In this review, we provide a history-aware discussion on these topics using both the recent and more distant past of our academic discipline and its relevant contexts.

The Opportunities and Challenges of Integrating Population Histories Into Genetic Studies for Diverse Populations: A Motivating Example From Native Hawaiians

There is a well-recognized need to include diverse populations in genetic studies, but several obstacles continue to be prohibitive, including (but are not limited to) the difficulty of recruiting individuals from diverse populations in large numbers and the lack of representation in available genomic references. These obstacles notwithstanding, studying multiple diverse populations would provide informative, population-specific insights. Using Native Hawaiians as an example of an understudied population with a unique evolutionary history, I will argue that by developing key genomic resources and integrating evolutionary thinking into genetic epidemiology, we will have the opportunity to efficiently advance our knowledge of the genetic risk factors, ameliorate health disparity, and improve healthcare in this underserved population.

Prenatal Noninvasive Trio-WES in a Case of Pregnancy-Related Liver Disorder

Liver disease in pregnancy may present as an acute condition related to the gestational period, characterized by pruritus, jaundice, and abnormal liver function. The disease may be misdiagnosed with other liver diseases, some of which may have consequences for fetal health. It is therefore advisable to implement rapid diagnostic strategies to provide information for the management of pregnancy in these conditions. We report the case of a healthy woman with a twin pregnancy from homologous in vitro fertilization (IVF), who in the third trimester presented jaundice and malaise. Biochemical investigations and liver hyperechogenicity raised the suspicion of acute fatty liver disease of pregnancy (AFLP). Non-invasive prenatal whole-exome sequencing (WES) in the trio identified the Phe305Ile heterozygous variant in the ATP8B1 gene. Considering the twin pregnancy, the percentage of the variant versus the wild allele was of 31%, suggesting heterozygosity present in the mother alone. This analysis showed that the mother was affected by benign recurrent intrahepatic cholestasis of pregnancy (ICP1: # 147480) and indicated the opportunity to anticipate childbirth to avoid worsening of the mother’s health. WES after the birth of the twins confirmed the molecular data.

A sex-specific evolutionary interaction between ADCY9 and CETP

Pharmacogenomic studies have revealed associations between rs1967309 in the adenylyl cyclase type 9 (ADCY9) gene and clinical responses to the cholesteryl ester transfer protein (CETP) modulator dalcetrapib, however, the mechanism behind this interaction is still unknown. Here, we characterized selective signals at the locus associated with the pharmacogenomic response in human populations and we show that rs1967309 region exhibits signatures of positive selection in several human populations. Furthermore, we identified a variant in CETP, rs158477, which is in long-range linkage disequilibrium with rs1967309 in the Peruvian population. The signal is mainly seen in males, a sex-specific result that is replicated in the LIMAA cohort of over 3400 Peruvians. Analyses of RNA-seq data further suggest an epistatic interaction on CETP expression levels between the two SNPs in multiple tissues, which also differs between males and females. We also detected interaction effects of the two SNPs with sex on cardiovascular phenotypes in the UK Biobank, in line with the sex-specific genotype associations found in Peruvians at these loci. We propose that ADCY9 and CETP coevolved during recent human evolution due to sex-specific selection, which points toward a biological link between dalcetrapib’s pharmacogene ADCY9 and its therapeutic target CETP.

The distribution of common-variant effect sizes

The genetic effect-size distribution of a disease describes the number of risk variants, the range of their effect sizes and sample sizes that will be required to discover them. Accurate estimation has been a challenge. Here I propose Fourier Mixture Regression (FMR), validating that it accurately estimates real and simulated effect-size distributions. Applied to summary statistics for ten diseases (average [Formula: see text]), FMR estimates that 100,000-1,000,000 cases will be required for genome-wide significant SNPs to explain 50% of SNP heritability. In such large studies, genome-wide significance becomes increasingly conservative, and less stringent thresholds achieve high true positive rates if confounding is controlled. Across traits, polygenicity varies, but the range of their effect sizes is similar. Compared with effect sizes in the top 10% of heritability, including most discovered thus far, those in the bottom 10-50% are orders of magnitude smaller and more numerous, spanning a large fraction of the genome.

Population and comparative genetics of thermotolerance divergence between yeast species

Many familiar traits in the natural world-from lions' manes to the longevity of bristlecone pine trees-arose in the distant past, and have long since fixed in their respective species. A key challenge in evolutionary genetics is to figure out how and why species-defining traits have come to be. We used the thermotolerance growth advantage of the yeast Saccharomyces cerevisiae over its sister species Saccharomyces paradoxus as a model for addressing these questions. Analyzing loci at which the S. cerevisiae allele promotes thermotolerance, we detected robust evidence for positive selection, including amino acid divergence between the species and conservation within S. cerevisiae populations. Because such signatures were particularly strong at the chromosome segregation gene ESP1, we used this locus as a case study for focused mechanistic follow-up. Experiments revealed that, in culture at high temperature, the S. paradoxus ESP1 allele conferred a qualitative defect in biomass accumulation and cell division relative to the S. cerevisiae allele. Only genetic divergence in the ESP1 coding region mattered phenotypically, with no functional impact detectable from the promoter. Our data support a model in which an ancient ancestor of S. cerevisiae, under selection to boost viability at high temperature, acquired amino acid variants at ESP1 and many other loci, which have been constrained since then. Complex adaptations of this type hold promise as a paradigm for interspecies genetics, especially in deeply diverged traits that may have taken millions of years to evolve.

The distinct morphological phenotypes of Southeast Asian aborigines are shaped by novel mechanisms for adaptation to tropical rainforests

Southeast Asian aborigines, the hunter-gatherer populations living in tropical rainforests, exhibit distinct morphological phenotypes, including short stature, dark skin, curly hair and a wide and snub nose. The underlying genetic architecture and evolutionary mechanism of these phenotypes remain a long-term mystery. We conducted whole genome deep sequencing of 81 Cambodian aborigines from eight ethnic groups. Through a genome-wide scan of selective sweeps, we discovered key genes harboring Cambodian-enriched mutations that may contribute to their phenotypes, including two hair morphogenesis genes (TCHH and TCHHL1), one nasal morphology gene (PAX3) and a set of genes (such as ENTPD1-AS1) associated with short stature. The identified new genes and novel mutations suggest an independent origin of the distinct phenotypes in Cambodian aborigines through parallel evolution, refuting the long-standing argument on the common ancestry of these phenotypes among the worldwide rainforest hunter-gatherers. Notably, our discovery reveals that various types of molecular mechanisms, including antisense transcription and epigenetic regulation, contribute to human morphogenesis, providing novel insights into the genetics of human environmental adaptation.

A worldwide map of swine short tandem repeats and their associations with evolutionary and environmental adaptations

Background

Short tandem repeats (STRs) are genetic markers with a greater mutation rate than single nucleotide polymorphisms (SNPs) and are widely used in genetic studies and forensics. However, most studies in pigs have focused only on SNPs or on a limited number of STRs.

Results

This study screened 394 deep-sequenced genomes from 22 domesticated pig breeds/populations worldwide, wild boars from both Europe and Asia, and numerous outgroup Suidaes, and identified a set of 878,967 polymorphic STRs (pSTRs), which represents the largest repository of pSTRs in pigs to date. We found multiple lines of evidence that pSTRs in coding regions were affected by purifying selection. The enrichment of trinucleotide pSTRs in coding sequences (CDS), 5′UTR and H3K4me3 regions suggests that trinucleotide STRs serve as important components in the exons and promoters of the corresponding genes. We demonstrated that, compared to SNPs, pSTRs provide comparable or even greater accuracy in determining the breed identity of individuals. We identified pSTRs that showed significant population differentiation between domestic pigs and wild boars in Asia and Europe. We also observed that some pSTRs were significantly associated with environmental variables, such as average annual temperature or altitude of the originating sites of Chinese indigenous breeds, among which we identified loss-of-function and/or expanded STRs overlapping with genes such as AHR, LAS1L and PDK1. Finally, our results revealed that several pSTRs show stronger signals in domestic pig—wild boar differentiation or association with the analysed environmental variables than the flanking SNPs within a 100-kb window.

Conclusions

This study provides a genome-wide high-density map of pSTRs in diverse pig populations based on genome sequencing data, enabling a more comprehensive characterization of their roles in evolutionary and environmental adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00631-4.

Admixture mapping identifies African and Amerindigenous local ancestry loci associated with fetal growth

Fetal growth is an important determinant of cardiometabolic disease risk during childhood and adulthood. The genetic architecture of fetal growth remains largely understudied in ancestrally diverse populations. We conducted genome-wide admixture mapping scan and analysis of genetic ancestry among Hispanic American, African American, European American, and Asian American pregnant women to identify genetic loci associated with fetal growth measures across 13-40 weeks gestation. Fetal growth measures were associated with genome-wide average African, European, Amerindigenous and East Asian ancestry proportions (P ranged from10^-3 to 4.8 × 10^-2). Admixture mapping analysis identified ten African ancestry loci and three Amerindigenous ancestry loci significantly associated with fetal growth measures at Bonferroni-corrected levels of significance (P ranged from 2.18 × 10^-8 to 3.71 × 10^-6). At the chr2q23.3-24.2 locus in which higher African ancestry was associated with long bone (femur and humerus) lengths, the T allele of rs13030825 (GALNT13) was associated with longer humerus length in African Americans (β = 0.44, P = 6.25 × 10^-6 at week 27; β = 0.39, P = 7.72 × 10^-5 at week 40). The rs13030825 SNP accounted for most of the admixture association at the chr2q23.3-24.2 locus and has substantial allele frequency difference between African and European reference samples (F_ST = 0.55, P = 0.03). Regulatory annotation shows that rs13030825 overlaps with the serum response factor (SRF) transcription factor previously implicated in postnatal bone development of mice. Overall, we identified ancestry-related maternal genetic loci that influence fetal growth, shedding light on molecular pathways that regulate fetal growth and potential effects on health across the lifespan.Clinical trials registration ClinicalTrials.gov, NCT00912132.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

Cut-off points of anthropometric markers associated with hypertension and diabetes in Peru: Demographic and Health Survey 2018

OBJECTIVE

To determine the optimal anthropometric cut-off points for predicting the likelihood ratios of hypertension and diabetes in the Peruvian population.

DESIGN

A cross-sectional study was performed to establish cut-off values for body mass index (BMI), waist circumference (WC), waist:height ratio (WHtR) and Conicity index (C-index) associated with increased risk of hypertension and diabetes. Youden's index (YIndex), area under the curve (AUC), sensitivity and specificity were calculated.

SETTING

Peruvian households.

PARTICIPANTS

Peruvian population over the age of 18 years.

RESULTS

A total of 31 553 subjects were included, 57 % being women. Among the women, 53·06 % belonged to the 25- to 44-year-old age group [mean age: 41·66 in men and 40·02 in women]. The mean BMI, WHtR and C-index values were higher in women 27·49, 0·61, 1·30, respectively, while the mean WC value was higher in men 92·12 cm (sd ± 11·28). The best predictors of hypertension in men were the WHtR (AUC = 0·64) and the C-index (AUC = 0·64) with an optimal cut-off point of 0·57 (YIndex = 0·284) and 1·301 (YIndex = 0·284), respectively. Women showed an AUC of 0·63 and 0·61 in the WHtR and C-index, respectively, with an optimal cut-off of 0·61 (YIndex = 0·236) and 1·323 (YIndex = 0·225). The best predictor for diabetes was the C-index: with an AUC = 0·67 and an optimal cut-off of 1·337 (YIndex = 0·346) for men and an AUC = 0·66 and optimal cut-off of 1·313 (YIndex = 0·319) for women.

CONCLUSIONS

Our findings show that in Peruvian adults, the WHtR and the C-index have the strongest association with hypertension in both sexes. Likewise, the C-index had the strongest association with diabetes.

Acromelic dysplasias: how rare musculoskeletal disorders reveal biological functions of extracellular matrix proteins

Acromelic dysplasias are a group of rare musculoskeletal disorders that collectively present with short stature, pseudomuscular build, stiff joints, and tight skin. Acromelic dysplasias are caused by mutations in genes (FBN1, ADAMTSL2, ADAMTS10, ADAMTS17, LTBP2, and LTBP3) that encode secreted extracellular matrix proteins, and in SMAD4, an intracellular coregulator of transforming growth factor-β (TGF-β) signaling. The shared musculoskeletal presentations in acromelic dysplasias suggest that these proteins cooperate in a biological pathway, but also fulfill distinct roles in specific tissues that are affected in individual disorders of the acromelic dysplasia group. In addition, most of the affected proteins directly interact with fibrillin microfibrils in the extracellular matrix and have been linked to the regulation of TGF-β signaling. Together with recently developed knockout mouse models targeting the affected genes, novel insights into molecular mechanisms of how these proteins regulate musculoskeletal development and homeostasis have emerged. Here, we summarize the current knowledge highlighting pathogenic mechanisms of the different disorders that compose acromelic dysplasias and provide an overview of the emerging biological roles of the individual proteins that are compromised. Finally, we develop a conceptual model of how these proteins may interact and form an "acromelic dysplasia complex" on fibrillin microfibrils in connective tissues of the musculoskeletal system.