Architecture Of Complex Traits Research Articles

Evaluation of heritability partitioning approaches in livestock populations

BackgroundHeritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium.ResultsHere, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively).ConclusionsHeritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used.

Joint modeling of human cortical structure: Genetic correlation network and composite-trait genetic correlation

Heritability and genetic covariance/correlation quantify the marginal and shared genetic effects across traits. They offer insights on the genetic architecture of complex traits and diseases. To explore how genetic variations contribute to brain function variations, we estimated heritability and genetic correlation across cortical thickness, surface area, and volume of 33 anatomically predefined regions in left and right hemispheres, using summary statistics of genome-wide association analyses of 31,968 participants in the UK Biobank. To characterize the relationships between these regions of interest, we constructed a genetic network for these regions using recursive two-way cut-offs in similarity matrices defined by genetic correlations. The inferred genetic network matches the brain lobe mapping more closely than the network inferred from phenotypic similarities. We further studied the associations between the genetic network for brain regions and 30 complex traits through a novel composite-linkage disequilibrium score regression method. We identified seven significant pairs, which offer insights on the genetic basis for regions of interest mediated by cortical measures.

Residual network improves the prediction accuracy of genomic selection.

Genetic improvement of complex traits in animal and plant breeding depends on the efficient and accurate estimation of breeding values. Deep learning methods have been shown to be not superior over traditional genomic selection (GS) methods, partially due to the degradation problem (i.e. with the increase of the model depth, the performance of the deeper model deteriorates). Since the deep learning method residual network (ResNet) is designed to solve gradient degradation, we examined its performance and factors related to its prediction accuracy in GS. Here we compared the prediction accuracy of conventional genomic best linear unbiased prediction, Bayesian methods (BayesA, BayesB, BayesC, and Bayesian Lasso), and two deep learning methods, convolutional neural network and ResNet, on three datasets (wheat, simulated and real pig data). ResNet outperformed other methods in both Pearson's correlation coefficient (PCC) and mean squared error (MSE) on the wheat and simulated data. For the pig backfat depth trait, ResNet still had the lowest MSE, whereas Bayesian Lasso had the highest PCC. We further clustered the pig data into four groups and, on one separated group, ResNet had the highest prediction accuracy (both PCC and MSE). Transfer learning was adopted and capable of enhancing the performance of both convolutional neural network and ResNet. Taken together, our findings indicate that ResNet could improve GS prediction accuracy, affected potentially by factors such as the genetic architecture of complex traits, data volume, and heterogeneity.

Integrating meta‐analysis of genome‐wide association study with Pig Genotype‐Tissue Expression resources uncovers the genetic architecture for age at first farrowing in pigs

AbstractAge at first farrowing (AFF) is a reproductive trait with low heritability and high importance in the pig industry. To enhance the statistical power of genome‐wide association study (GWAS) and further explore the genetic nature of AFF, we first conducted GWAS meta‐analysis using three Yorkshire populations, and then integrated the Pig Genotype‐Tissue Expression (PigGTEx) resources to interpret their potential regulatory mechanism. Additionally, we compared the AFF in pig with the age at first birth (AFB) of human using GWAS summary statistics. We detected 18 independent variants in GWAS meta‐analysis and 8 genes in gene‐based association analysis significantly associated with AFF. By integrating the PigGTEx resource, we conducted transcriptome‐wide association study (TWAS) and colocalization analysis on 34 pig tissues. In TWAS, we detected 18 significant gene‐tissue pairs, such as DCAF6 in uterus and CREG1 in blood. In colocalization, we found 111 potential candidate tissue‐gene pairs, such as GJD4 and LYPLAL1. We found that the homologous gene, CHST10, might be the potential candidate gene between humans in AFB and pigs in AFF. In conclusion, integrating GWAS meta‐analysis and PigGTEx resources is a meaningful way to decipher the genetic architecture of complex traits. We found that DCAF6, CREG1, GJD4, and LYPLAL1 are candidate genes with high reliability for AFF in swine. The comparative analysis showed that CHST10 might play a potentially critical role in AFB/AFF across human and pigs.

MPH: fast REML for large-scale genome partitioning of quantitative genetic variation.

Genome partitioning of quantitative genetic variation is useful for dissecting the genetic architecture of complex traits. However, existing methods, such as Haseman-Elston regression and linkage disequilibrium score regression, often face limitations when handling extensive farm animal datasets, as demonstrated in this study. To overcome this challenge, we present MPH, a novel software tool designed for efficient genome partitioning analyses using restricted maximum likelihood. The computational efficiency of MPH primarily stems from two key factors: the utilization of stochastic trace estimators and the comprehensive implementation of parallel computation. Evaluations with simulated and real datasets demonstrate that MPH achieves comparable accuracy and significantly enhances convergence, speed, and memory efficiency compared to widely used tools like GCTA and LDAK. These advancements facilitate large-scale, comprehensive analyses of complex genetic architectures in farm animals. The MPH software is available at https://jiang18.github.io/mph/.

Evolutionary history and biological adaptation of Han Chinese people on the Mongolian Plateau

Complex demographic processes and natural selection pressures are critical to resolving patterns of the molecular genetic basis of adaptive traits or complex diseases. Recent ancient genome data allow us to trace how key traits evolved in different human populations over time, connecting ancient population history with disease susceptibility in western Eurasians. To fill this gap in eastern Eurasians and provide deep insights into the evolutionary history of the population-specific genetic basis of complex biological traits, we explored the evolutionary basis of adaptive traits in one integrative modern and ancient genomic database, including 225 out of 5583 genomes first reported here. We comprehensively characterized the demographic history and biological adaptation of Han Chinese individuals on the Mongolian Plateau based on the allele frequency spectrum and haplotype-resolved fragments. We found strong genetic homogeneity among geographically different Han populations from Inner Mongolia (IMH). We reconstructed their admixture models and demographic events, revealing that IMH had a close genetic relationship with ancient millet farmers and obtained additional gene flow from Altaic-speaking populations. The enrichment of selected candidate genes suggested that essential metabolism-related genes promoted the rapid adaptation of eastern Eurasians to ancient environmental shifts and dietary changes during agricultural innovations. Evolutionary trajectory reconstruction of metabolism-related genes of methylenetetrahydrofolate reductase (MTHFR) and fatty acid desaturase 1 (FADS1) suggested that ancient dietary shifts during the Neolithic transition period promoted the differentiated metabolic rate of folate and fatty acid. We revealed the polygenicity of biological adaptation traits and the pleiotropy of adaptive genes, indicating that complex recent polygenic adaptations, genetic interactions, and genotype–phenotype correlations have contributed to the genetic architecture of complex traits in eastern Eurasians.

Differentiated genomic footprints suggest isolation and long-distance migration of Hmong-Mien populations

BackgroundThe underrepresentation of Hmong-Mien (HM) people in Asian genomic studies has hindered our comprehensive understanding of the full landscape of their evolutionary history and complex trait architecture. South China is a multi-ethnic region and indigenously settled by ethnolinguistically diverse HM, Austroasiatic (AA), Tai-Kadai (TK), Austronesian (AN), and Sino-Tibetan (ST) people, which is regarded as East Asia’s initial cradle of biodiversity. However, previous fragmented genetic studies have only presented a fraction of the landscape of genetic diversity in this region, especially the lack of haplotype-based genomic resources. The deep characterization of demographic history and natural-selection-relevant genetic architecture of HM people was necessary.ResultsWe reported one HM-specific genomic resource and comprehensively explored the fine-scale genetic structure and adaptative features inferred from the genome-wide SNP data of 440 HM individuals from 33 ethnolinguistic populations, including previously unreported She. We identified solid genetic differentiation between HM people and Han Chinese at 7.64‒15.86 years ago (kya) and split events between southern Chinese inland (Miao/Yao) and coastal (She) HM people in the middle Bronze Age period and the latter obtained more gene flow from Ancient Northern East Asians. Multiple admixture models further confirmed that extensive gene flow from surrounding ST, TK, and AN people entangled in forming the gene pool of Chinese coastal HM people. Genetic findings of isolated shared unique ancestral components based on the sharing alleles and haplotypes deconstructed that HM people from the Yungui Plateau carried the breadth of previously unknown genomic diversity. We identified a direct and recent genetic connection between Chinese inland and Southeast Asian HM people as they shared the most extended identity-by-descent fragments, supporting the long-distance migration hypothesis. Uniparental phylogenetic topology and network-based phylogenetic relationship reconstruction found ancient uniparental founding lineages in southwestern HM people. Finally, the population-specific biological adaptation study identified the shared and differentiated natural selection signatures among inland and coastal HM people associated with physical features and immune functions. The allele frequency spectrum of cancer susceptibility alleles and pharmacogenomic genes showed significant differences between HM and northern Chinese people.ConclusionsOur extensive genetic evidence combined with the historical documents supported the view that ancient HM people originated from the Yungui regions associated with ancient “Three-Miao tribes” descended from the ancient Daxi-Qujialing-Shijiahe people. Then, some have recently migrated rapidly to Southeast Asia, and some have migrated eastward and mixed respectively with Southeast Asian indigenes, Liangzhu-related coastal ancient populations, and incoming southward ST people. Generally, complex population migration, admixture, and adaptation history contributed to the complicated patterns of population structure of geographically diverse HM people.

Multi-GPA-Tree: Statistical approach for pleiotropy informed and functional annotation tree guided prioritization of GWAS results.

Genome-wide association studies (GWAS) have successfully identified over two hundred thousand genotype-trait associations. Yet some challenges remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), most with small or moderate effect sizes, making them difficult to detect. Second, many complex traits share a common genetic basis due to 'pleiotropy' and and though few methods consider it, leveraging pleiotropy can improve statistical power to detect genotype-trait associations with weaker effect sizes. Third, currently available statistical methods are limited in explaining the functional mechanisms through which genetic variants are associated with specific or multiple traits. We propose multi-GPA-Tree to address these challenges. The multi-GPA-Tree approach can identify risk SNPs associated with single as well as multiple traits while also identifying the combinations of functional annotations that can explain the mechanisms through which risk-associated SNPs are linked with the traits. First, we implemented simulation studies to evaluate the proposed multi-GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that multi-GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs for multiple traits. Second, we applied multi-GPA-Tree to a systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), and to a Crohn's disease (CD) and ulcertive colitis (UC) GWAS, and functional annotation data including GenoSkyline and GenoSkylinePlus. Our results demonstrate that multi-GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits.

Chemical-tag-based semi-annotated metabolomics facilitates gene identification and specialized metabolic pathway elucidation in wheat.

The importance of metabolite modification and species-specific metabolic pathways has long been recognized. However, linking the chemical structure of metabolites to gene function in order to explore the genetic and biochemical basis of metabolism has not yet been reported in wheat (Triticum aestivum). Here, we profiled metabolic fragment enrichment in wheat leaves and consequently applied chemical-tag-based semi-annotated metabolomics in a genome-wide association study in accessions of wheat. The studies revealed that all 1,483 quantified metabolites have at least one known functional group whose modification is tailored in an enzyme-catalyzed manner and eventually allows efficient candidate gene mining. A Triticeae crop-specific flavonoid pathway and its underlying metabolic gene cluster were elucidated in further functional studies. Additionally, upon overexpressing the major effect gene of the cluster TraesCS2B01G460000 (TaOMT24), the pathway was reconstructed in rice (Oryza sativa), which lacks this pathway. The reported workflow represents an efficient and unbiased approach for gene mining using forward genetics in hexaploid wheat. The resultant candidate gene list contains vast molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets and will ultimately aid in achieving wheat crop improvement.

Sequence-based GWAS meta-analyses for beef production traits

BackgroundCombining the results of within-population genome-wide association studies (GWAS) based on whole-genome sequences into a single meta-analysis (MA) is an accurate and powerful method for identifying variants associated with complex traits. As part of the H2020 BovReg project, we performed sequence-level MA for beef production traits. Five partners from France, Switzerland, Germany, and Canada contributed summary statistics from sequence-based GWAS conducted with 54,782 animals from 15 purebred or crossbred populations. We combined the summary statistics for four growth, nine morphology, and 15 carcass traits into 16 MA, using both fixed effects and z-score methods.ResultsThe fixed-effects method was generally more informative to provide indication on potentially causal variants, although we combined substantially different traits in each MA. In comparison with within-population GWAS, this approach highlighted (i) a larger number of quantitative trait loci (QTL), (ii) QTL more frequently located in genomic regions known for their effects on growth and meat/carcass traits, (iii) a smaller number of genomic variants within the QTL, and (iv) candidate variants that were more frequently located in genes. MA pinpointed variants in genes, including MSTN, LCORL, and PLAG1 that have been previously associated with morphology and carcass traits. We also identified dozens of other variants located in genes associated with growth and carcass traits, or with a function that may be related to meat production (e.g., HS6ST1, HERC2, WDR75, COL3A1, SLIT2, MED28, and ANKAR). Some of these variants overlapped with expression or splicing QTL reported in the cattle Genotype-Tissue Expression atlas (CattleGTEx) and could therefore regulate gene expression.ConclusionsBy identifying candidate genes and potential causal variants associated with beef production traits in cattle, MA demonstrates great potential for investigating the biological mechanisms underlying these traits. As a complement to within-population GWAS, this approach can provide deeper insights into the genetic architecture of complex traits in beef cattle.

Mexican Biobank advances population and medical genomics of diverse ancestries.

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico ata resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

Environmental context of phenotypic plasticity in flowering time in sorghum and rice.

Phenotypic plasticity is an important topic in biology and evolution. However, how to generate broadly applicable insights from individual studies remains a challenge. Here, with flowering time observed from a large geographical region for sorghum and rice genetic populations, we examine the consistency of parameter estimation for reaction norms of genotypes across different subsets of environments and searched for potential strategies to inform the study design. Both sample size and environmental mean range of the subset affected the consistency. The subset with either a large range of environmental mean or a large sample size resulted in genetic parameters consistent with the overall pattern. Furthermore, high accuracy through genomic prediction was obtained for reaction norm parameters of untested genotypes using models built from tested genotypes under the subsets of environments with either a large range or a large sample size. With 1428 and 1674 simulated settings, our analyses suggested that the distribution of environmental index values of a site should be considered in designing experiments. Overall, we showed that environmental context was critical, and considerations should be given to better cover the intended range of the environmental variable. Our findings have implications for the genetic architecture of complex traits, plant-environment interaction, and climate adaptation.

Statistical methods for gene-environment interaction analysis.

Most human complex phenotypes result from multiple genetic and environmental factors and their interactions. Understanding the mechanisms by which genetic and environmental factors interact offers valuable insights into the genetic architecture of complex traits and holds great potential for advancing precision medicine. The emergence of large population biobanks has led to the development of numerous statistical methods aiming at identifying gene-environment interactions (G × E). In this review, we present state-of-the-art statistical methodologies for G × E analysis. We will survey a spectrum of approaches for single-variant G × E mapping, followed by various techniques for polygenic G × E analysis. We conclude this review with a discussion on the future directions and challenges in G × E research.

Development of a 20 K SNP array for the leopard coral grouper, Plectropomus leopardus

Single nucleotide polymorphism (SNP) genotyping array is a powerful tool for providing information on polymorphic variations of species of interest. These genomic variations can be used to investigate the genetic architecture of complex traits, population structures, and genomic selection of breeding programs. The aquaculture of the leopard coral grouper (Plectropomus leopardus), is rapidly developing in South China and Southeast Asia. However, the lack of a high-efficient genotyping platform has hampered the progress of both genetic studies and selective breeding practices of this fish. Here we developed and evaluated a 20 K SNP array namely “DongXin-I” for P. leopardus. Using the genome re-sequencing data of 730 leopard coral groupers, we identified over 3.67 million SNPs, from which 20,737 loci, with even distribution across the genome, were finally involved in the SNP array. Evaluation of the array performance with 232 fish showed a high overall SNP call rates (> 99%) and an average minor allele frequency (MAF) of 0.312. Comparison of genotyping results of 87 fish from both the “DongXin-I” array and genome resequencing demonstrated a high genotyping repeatability and accurate performance of the genotyping array. Furthermore, using the gentyping results generated by the “DongXin-I” array, we revealed the genetic structure of breeding populaitons currently used the in Hainan Province. To summarize, the “DongXin-I” array should be a useful tool for the genetic studies and genomic selection breeding in the aquaculture industry of the coral leopard grouper.

Projecting genetic associations through gene expression patterns highlights disease etiology and drug mechanisms

Genes act in concert with each other in specific contexts to perform their functions. Determining how these genes influence complex traits requires a mechanistic understanding of expression regulation across different conditions. It has been shown that this insight is critical for developing new therapies. Transcriptome-wide association studies have helped uncover the role of individual genes in disease-relevant mechanisms. However, modern models of the architecture of complex traits predict that gene-gene interactions play a crucial role in disease origin and progression. Here we introduce PhenoPLIER, a computational approach that maps gene-trait associations and pharmacological perturbation data into a common latent representation for a joint analysis. This representation is based on modules of genes with similar expression patterns across the same conditions. We observe that diseases are significantly associated with gene modules expressed in relevant cell types, and our approach is accurate in predicting known drug-disease pairs and inferring mechanisms of action. Furthermore, using a CRISPR screen to analyze lipid regulation, we find that functionally important players lack associations but are prioritized in trait-associated modules by PhenoPLIER. By incorporating groups of co-expressed genes, PhenoPLIER can contextualize genetic associations and reveal potential targets missed by single-gene strategies.

UNRAVELING THE GENETIC ARCHITECTURE OF COMPLEX TRAITS IN PLANTS

Complex traits are those whose inheritance does not follow simple and predictable patterns. They are not governed by a single locus, instead, they are determined by several loci and are influenced by the environment. Most of the traits with agronomic interest and economic importance such as resistance to biotic and abiotic stress, and yield, among others, are quantitative traits and their study is based on dissecting the underlying genetic architecture, the number of loci responsible for the variance of a quantitative trait, the relevant contribution made by each locus and their interaction with the environment. This review provides the most relevant conceptual bases for the study of the genetic architecture of complex quantitative traits in plants. The methodologies that allow identifying the loci and candidate genes that govern this type of traits are described, such as QTL mapping by linkage and association mapping. In addition, the incorporation of these loci in phenotype prediction strategies such as marker-assisted selection and genomic selection, exhibits the benefits and limitations of these approaches. Finally, the challenges and perspectives facing the study of the genetic architecture of complex traits in plants are discussed.

Gene-based burden scores identify rare variant associations for 28 blood biomarkers

BackgroundA relevant part of the genetic architecture of complex traits is still unknown; despite the discovery of many disease-associated common variants. Polygenic risk score (PRS) models are based on the evaluation of the additive effects attributable to common variants and have been successfully implemented to assess the genetic susceptibility for many phenotypes. In contrast, burden tests are often used to identify an enrichment of rare deleterious variants in specific genes. Both kinds of genetic contributions are typically analyzed independently. Many studies suggest that complex phenotypes are influenced by both low effect common variants and high effect rare deleterious variants. The aim of this paper is to integrate the effect of both common and rare functional variants for a more comprehensive genetic risk modeling.MethodsWe developed a framework combining gene-based scores based on the enrichment of rare functionally relevant variants with genome-wide PRS based on common variants for association analysis and prediction models. We applied our framework on UK Biobank dataset with genotyping and exome data and considered 28 blood biomarkers levels as target phenotypes. For each biomarker, an association analysis was performed on full cohort using gene-based scores (GBS). The cohort was then split into 3 subsets for PRS construction and feature selection, predictive model training, and independent evaluation, respectively. Prediction models were generated including either PRS, GBS or both (combined).ResultsAssociation analyses of the cohort were able to detect significant genes that were previously known to be associated with different biomarkers. Interestingly, the analyses also revealed heterogeneous effect sizes and directionality highlighting the complexity of the blood biomarkers regulation. However, the combined models for many biomarkers show little or no improvement in prediction accuracy compared to the PRS models.ConclusionThis study shows that rare variants play an important role in the genetic architecture of complex multifactorial traits such as blood biomarkers. However, while rare deleterious variants play a strong role at an individual level, our results indicate that classical common variant based PRS might be more informative to predict the genetic susceptibility at the population level.

Trumpet plots: visualizing the relationship between allele frequency and effect size in genetic association studies.

Recent advances in genome-wide association and sequencing studies have shown that the genetic architecture of complex traits and diseases involves a combination of rare and common genetic variants distributed throughout the genome. One way to better understand this architecture is to visualize genetic associations across a wide range of allele frequencies. However, there is currently no standardized or consistent graphical representation for effectively illustrating these results. Here we propose a standardized approach for visualizing the effect size of risk variants across the allele frequency spectrum. The proposed plots have a distinctive trumpet shape: with the majority of variants having high frequency and small effects, and a small number of variants having lower frequency and larger effects. To demonstrate the utility of trumpet plots in illustrating the relationship between the number of variants, their frequency, and the magnitude of their effects in shaping the genetic architecture of complex traits and diseases, we generated trumpet plots for more than one hundred traits in the UK Biobank. To facilitate their broader use, we developed an R package, 'TrumpetPlots' (available at the Comprehensive R Archive Network) and R Shiny application, 'Shiny Trumpets' (available at https://juditgg.shinyapps.io/shinytrumpets/) that allows users to explore these results and submit their own data.

Characterizing the polygenic architecture of complex traits in populations of East Asian and European descent

To investigate the polygenicity of complex traits in populations of East Asian (EAS) and European (EUR) descents, we leveraged genome-wide data from Biobank Japan, UK Biobank, and FinnGen cohorts. Specifically, we analyzed up to 215 outcomes related to 18 health domains, assessing their polygenic architecture via descriptive statistics, such as the proportion of susceptibility SNPs per trait (πc). While we did not observe EAS–EUR differences in the overall distribution of polygenicity parameters across the phenotypes investigated, there were ancestry-specific patterns in the polygenicity differences between health domains. In EAS, pairwise comparisons across health domains showed enrichment for πc differences related to hematological and metabolic traits (hematological fold-enrichment = 4.45, p = 2.15 × 10–7; metabolic fold-enrichment = 4.05, p = 4.01 × 10–6). For both categories, the proportion of susceptibility SNPs was lower than that observed for several other health domains (EAS-hematological median πc = 0.15%, EAS-metabolic median πc = 0.18%) with the strongest πc difference with respect to respiratory traits (EAS-respiratory median πc = 0.50%; hematological-p = 2.26 × 10–3; metabolic-p = 3.48 × 10–3). In EUR, pairwise comparisons showed multiple πc differences related to the endocrine category (fold-enrichment = 5.83, p = 4.76 × 10–6), where these traits showed a low proportion of susceptibility SNPs (EUR-endocrine median πc = 0.01%) with the strongest difference with respect to psychiatric phenotypes (EUR-psychiatric median πc = 0.50%; p = 1.19 × 10–4). Simulating sample sizes of 1,000,000 and 5,000,000 individuals, we also showed that ancestry-specific polygenicity patterns translate into differences across health domains in the genetic variance explained by susceptibility SNPs projected to be genome-wide significant (e.g., EAS hematological-neoplasm p = 2.18 × 10–4; EUR endocrine-gastrointestinal p = 6.80 × 10–4). These findings highlight that traits related to the same health domains may present ancestry-specific variability in their polygenicity.

GWLD: an R package for genome-wide linkage disequilibrium analysis.

Linkage disequilibrium (LD) analysis is fundamental to the investigation of the genetic architecture of complex traits (e.g. human disease, animal and plant breeding) and population structure and evolution dynamics. However, until now, studies primarily focus on LD status between genetic variants located on the same chromosome. Moreover, genome (re)sequencing produces unprecedented numbers of genetic variants, and fast LD computation becomes a challenge. Here, we have developed GWLD, a parallelized and generalized tool designed for the rapid genome-wide calculation of LD values, including conventional D/D', r2, and (reduced) mutual information (MI and RMI) measures. LD between genetic variants within and across chromosomes can be rapidly computed and visualized in either an R package or a standalone C++ software package. To evaluate the accuracy and speed of LD calculation, we conducted comparisons using 4 real datasets. Interchromosomal LD patterns observed potentially reflect levels of selection intensity across different species. Both versions of GWLD, the R package (https://github.com/Rong-Zh/GWLD/GWLD-R) and the standalone C++ software (https://github.com/Rong-Zh/GWLD/GWLD-C++), are freely available on GitHub.