Accounting for age of onset and family history improves power in genome-wide association studies.

Genome-wide association studies (GWASs) have revolutionized human genetics, allowing researchers to identify thousands of disease-related genes and possible drug targets. However, case-control status does not account for the fact that not all controls may have lived through their period of risk for the disorder of interest. This can be quantified by examining the age-of-onset distribution and the age of the controls or the age of onset for cases. The age-of-onset distribution may also depend on information such as sex and birth year. In addition, family history is not routinely included in the assessment of control status. Here, we present LT-FH++, an extension of the liability threshold model conditioned on family history (LT-FH), which jointly accounts for age of onset and sex as well as family history. Using simulations, we show that, when family history and the age-of-onset distribution are available, the proposed approach yields statistically significant power gains over LT-FH and large power gains over genome-wide association study by proxy (GWAX). We applied our method to four psychiatric disorders available in the iPSYCH data and to mortality in the UK Biobank and found 20 genome-wide significant associations with LT-FH++, compared to ten for LT-FH and eight for a standard case-control GWAS. As more genetic data with linked electronic health records become available to researchers, we expect methods that account for additional health information, such as LT-FH++, to become even more beneficial.

Persistent Gene Flow Suggests an Absence of Reproductive Isolation in an African Antelope Speciation Model.

African antelope diversity is a globally unique vestige of a much richer world-wide Pleistocene megafauna. Despite this, the evolutionary processes leading to the prolific radiation of African antelopes are not well understood. Here, we sequenced 145 whole genomes from both subspecies of the waterbuck (Kobus ellipsiprymnus), an African antelope believed to be in the process of speciation. We investigated genetic structure and population divergence and found evidence of a mid-Pleistocene separation on either side of the eastern Great Rift Valley, consistent with vicariance caused by a rain shadow along the so-called 'Kingdon's Line'. However, we also found pervasive evidence of both recent and widespread historical gene flow across the Rift Valley barrier. By inferring the genome-wide landscape of variation among subspecies, we found 14 genomic regions of elevated differentiation, including a locus that may be related to each subspecies' distinctive coat pigmentation pattern. We investigated these regions as candidate speciation islands. However, we observed no significant reduction in gene flow in these regions, nor any indications of selection against hybrids. Altogether, these results suggest a pattern whereby climatically driven vicariance is the most important process driving the African antelope radiation, and suggest that reproductive isolation may not set in until very late in the divergence process. This has a significant impact on taxonomic inference, as many taxa will be in a gray area of ambiguous systematic status, possibly explaining why it has been hard to achieve consensus regarding the species status of many African antelopes. Our analyses demonstrate how population genetics based on low-depth whole genome sequencing can provide new insights that can help resolve how far lineages have gone along the path to speciation.

A comprehensive map of genetic relationships among diagnostic categories based on 48.6 million relative pairs from the Danish genealogy.

For more than half a century, Denmark has maintained population-wide demographic, health care, and socioeconomic registers that provide detailed information on the interaction between all residents and the extensive national social services system. We leverage this resource to reconstruct the genealogy of the entire nation based on all individuals legally residing in Denmark since 1968. We cross-reference 6,691,426 individuals with nationwide health care registers to estimate heritability and genetic correlations of 10 broad diagnostic categories involving all major organs and systems. Heritability estimates for mental disorders were consistently the highest across demographic cohorts (average h2 = 0.406, 95% CI = [0.403, 0.408]), whereas estimates for cancers were the lowest (average h2 = 0.130, 95% CI = [0.125, 0.134]). The average genetic correlation of each of the 10 diagnostic categories with the other nine was highest for gastrointestinal conditions (average rg = 0.567, 95% CI = [0.566, 0.567]) and lowest for urogenital conditions (average rg = 0.386, 95% CI = [0.385, 0.388]). Mental, pulmonary, gastrointestinal, and neurological conditions had similar genetic correlation profiles.

Worldwide distribution of genetic factors related to severity of COVID-19 infection.

BACKGROUND: Genome-wide association studies of COVID-19 severity have been carried out mostly on European or East Asian populations with small representation of other world regions. Here we explore the worldwide distribution and linkage disequilibrium (LD) patterns of genetic variants previously associated with COVID-19 severity. METHODS: We followed up the results of a large Spanish genome-wide meta-analysis on 26 populations from the 1000 Genomes Project by calculating allele frequencies and LD scores of the nine most significant SNPs. We also used the entire set of summary statistics to compute polygenic risk scores (PRSs) and carried out comparisons at the population and continental level. RESULTS: We observed the strongest differences among continental regions for the five top SNPs in chromosome 3. European, American, and South Asian populations showed similar LD patterns. Average PRSs in South Asian and American populations were consistently higher than those observed in Europeans. While PRS distributions were similar among South Asians, the American populations showed striking differences among them. CONCLUSIONS: Considering the caveats of PRS transferability across ethnicities, our analysis showed that American populations present the highest genetic risk score, hence potentially higher propensity, for COVID-19 severity. Independent validation is warranted with additional summary statistics and phenotype data.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Estimating narrow-sense heritability using family data from admixed populations.

Estimating total narrow-sense heritability in admixed populations remains an open question. In this work, we used extensive simulations to evaluate existing linear mixed-model frameworks for estimating total narrow-sense heritability in two population-based cohorts from Greenland, and compared the results with data from unadmixed individuals from Denmark. When our analysis focused on Greenlandic sib pairs, and under the assumption that shared environment among siblings has a negligible effect, the model with two relationship matrices, one capturing identity by descent and one capturing identity by state, returned heritability estimates close to the true simulated value, while using each of the two matrices alone led to downward biases. When phenotypes correlated with ancestry, heritability estimates were inflated. Based on these observations, we propose a PCA-based adjustment that recovers the true simulated heritability. We use this knowledge to estimate the heritability of ten quantitative traits from the two Greenlandic cohorts, and report differences such as lower heritability for height in Greenlanders compared with Europeans. In conclusion, narrow-sense heritability in admixed populations is best estimated when using a mixture of genetic relationship matrices on individuals with at least one first-degree relative included in the sample.

Sympathetic Nerve Injury in Thyroid Cancer.

The double innervation of the thyroid comes from the sympathetic and parasympathetic nervous system. Injury rates during surgery are at 30% but can be minimized by upwardly preparing the thyroid vessels at the level of thyroid capsule. Several factors have been accused of increasing the risk of injury including age and tumor size. Our aim was to investigate of there is indeed any possible correlations between these factors and a possible increase in injury rates following thyroidectomy. Seven studies were included in the meta-analysis. Statistical correlation was observed for a positive relationship between injury of the sympathetic nerve and thyroid malignancy surgery (p 2 = 74%) No statistical correlations were observed for a negative or positive relationship between injury of the sympathetic nerve and tumor size. There was also no statistically significant value observed for the correlation of the patients' age with the risk of sympathetic nerve injury (p = 0.388). Lack of significant correlation reported could be due to the small number of studies and great heterogeneity between them.

Heritability of bone mineral density in a multivariate family-based study.

There is evidence for a genetic contribution to bone mineral density (BMD×). Different loci affecting BMD have been identified by diverse linkage and genome-wide association studies. We studied the heritability of and the correlations among six densitometric phenotypes and four bone mass/fracture phenotypes. For this purpose, we used a family-based study of the genetics of osteoporosis, the Genetic Analysis of Osteoporosis Project. The primary aim of our study was to examine the roles of genetic and environmental factors in determining osteoporosis-related phenotypes. The project consisted of 11 extended families from Spain. All of them were selected through a proband with osteoporosis. BMD was measured using dual-energy X-ray absorptiometry. The proportion of variance of BMD attributable to significant covariates ranged from 25% (for femoral neck BMD) to 48% (for whole-body total BMD). The vast majority of the densitometric phenotypes had highly significant heritability, ranging from 0.252 (whole-body total BMD) to 0.537 (trochanteric BMD) after correcting for covariate effects. All of the densitometric phenotypes showed high and significant genetic correlations (from -0.772 to -1.000) with a low bone mass/osteopenia condition (Affected 3). Our findings provide additional evidence on the heritability of BMD and a strong genetic correlation between BMD and bone mass/fracture phenotypes in a Spanish population. Our results emphasize the importance of detecting genetic risk factors and the benefit of early diagnosis and especially therapeutic and preventive strategies.

Quantifying the relative importance of genetics and environment on the comorbidity between mental and cardiometabolic disorders using 17 million Scandinavians.

Mental disorders are leading causes of disability and premature death worldwide, partly due to high comorbidity with cardiometabolic disorders. Reasons for this comorbidity are still poorly understood. We leverage nation-wide health records and near-complete genealogies of Denmark and Sweden (n = 17 million) to reveal the genetic and environmental contributions underlying the observed comorbidity between six mental disorders and 15 cardiometabolic disorders. Genetic factors contributed about 50% to the comorbidity of schizophrenia, affective disorders, and autism spectrum disorder with cardiometabolic disorders, whereas the comorbidity of attention-deficit/hyperactivity disorder and anorexia with cardiometabolic disorders was mainly or fully driven by environmental factors. In this work we provide causal insight to guide clinical and scientific initiatives directed at achieving mechanistic understanding as well as preventing and alleviating the consequences of these disorders.

Quantifying the Relative Importance of Genetics and Environment on the Comorbidity between Mental- and Cardiometabolic Disorders: A Comprehensive Analysis of National Register Data from 17 million Scandinavians.

Mental disorders (MDs) are leading causes of disability and premature death worldwide, partly due to high comorbidity with cardiometabolic disorders (CMDs). Reasons for this comorbidity are still poorly understood. We leverage nation-wide health records and complete genealogies of Denmark and Sweden (n=17 million) to reveal the genetic and environmental contributions underlying the observed comorbidity between six MDs and 14 CMDs. Genetic factors contributed about 50% to the comorbidity of schizophrenia, affective disorders, and autism spectrum disorder with CMDs, whereas the comorbidity of attention-deficit/hyperactivity disorder and anorexia with CMDs was mainly or fully driven by environmental factors. These findings provide causal insight to guide clinical and scientific initiatives directed at achieving mechanistic understanding as well as preventing and alleviating the consequences of these disorders.

Spatial principal component analysis points at global genetic structure in the Western Mediterranean.

Many studies have tackled the existence of a genetic barrier in the Strait of Gibraltar between Iberian and North African populations, often with controversial conclusions. Here, we address this issue using a collection of Western Mediterranean populations and two dimensionality reduction methods: principal component analysis (PCA) and spatial PCA (sPCA). Our four different data sets consisted of (i) 16 polymorphic Alu insertions in 12 populations; (ii) 35 single-nucleotide polymorphisms in 13 populations; (iii) 13 short tandem repeats in 11 populations; and (iv) all 64 markers in 9 populations. In all PCA plots, South European and North African samples were visually distinguishable along the first PC. Several smaller clusters were also identifiable, especially on the African side of our geographical setting. sPCA indicated a single global structure for each of the marker sets and no local structures. These results are more compatible with a clinal distribution of allele frequencies rather than with abrupt changes, suggesting that isolation-by-distance, rather than a barrier to gene flow, is a more likely mechanism of genetic differentiation in the Western Mediterranean. An alternative/complementary explanation is progressive introgression from North African to Southwestern European populations.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses.

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains.

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention.

Apolipoprotein E/C1/C4/C2 gene cluster diversity in two native Andean populations: Aymaras and Quechuas.

The APOE/C1/C4/C2 gene cluster presents high relevance in lipid metabolism and, therefore, has important epidemiological implications. Here, we study for the first time the variation patterns of 25 polymorphisms (10 short tandem repeats, STRs, and 15 single nucleotide polymorphismas, SNPs) in two native Andean samples from Bolivia (45 Aymaras and 45 Quechuas) as well as one European sample (n = 41) as external reference. We estimated diversity parameters, linkage disequilibrium patterns, population structure, and possible selective effects. In general, diversity was low and could be partly attributed to selection (probably due to its physiological importance), since the APOE/C1/C4/C2 region was highly conserved compared to the flanking genes in both Bolivians and Europeans. Moreover, the lower gene diversity in Bolivians compared to Europeans for some markers might indicate different demographic histories. Regarding the APOE isoforms, in addition to ɛ3 (94%) and ɛ4 (5%), isoform ɛ2 (1%) was also detected in Bolivians. In relation to previous hypotheses, our results support that genetic drift or founder effects rather than selection for increased cholesterol absorption are the main factors that have shaped the distribution of APOE isoforms observed in South America.

Cross-trait assortative mating is widespread and inflates genetic correlation estimates.

The observation of genetic correlations between disparate human traits has been interpreted as evidence of widespread pleiotropy. Here, we introduce cross-trait assortative mating (xAM) as an alternative explanation. We observe that xAM affects many phenotypes and that phenotypic cross-mate correlation estimates are strongly associated with genetic correlation estimates (R2=74%). We demonstrate that existing xAM plausibly accounts for substantial fractions of genetic correlation estimates and that previously reported genetic correlation estimates between some pairs of psychiatric disorders are congruent with xAM alone. Finally, we provide evidence for a history of xAM at the genetic level using cross-trait even/odd chromosome polygenic score correlations. Together, our results demonstrate that previous reports have likely overestimated the true genetic similarity between many phenotypes.

Fine-scale population structure in five rural populations from the Spanish Eastern Pyrenees using high-coverage whole-genome sequence data.

The area of the Spanish Pyrenees is particularly interesting for studying the demographic dynamics of European rural areas given its orography, the main traditional rural condition of its population and the reported higher patterns of consanguinity of the region. Previous genetic studies suggest a gradient of genetic continuity of the area in the West to East axis. However, it has been shown that micro-population substructure can be detected when considering high-quality NGS data and using spatial explicit methods. In this work, we have analyzed the genome of 30 individuals sequenced at 40× from five different valleys in the Spanish Eastern Pyrenees (SEP) separated by less than 140 km along a west to east axis. Using haplotype-based methods and spatial analyses, we have been able to detect micro-population substructure within SEP not seen in previous studies. Linkage disequilibrium and autozygosity analyses suggest that the SEP populations show diverse demographic histories. In agreement with these results, demographic modeling by means of ABC-DL identify heterogeneity in their effective population sizes despite of their close geographic proximity, and suggests that the population substructure within SEP could have appeared around 2500 years ago. Overall, these results suggest that each rural population of the Pyrenees could represent a unique entity.

The Mediterranean Sea as a barrier to gene flow: evidence from variation in and around the F7 and F12 genomic regions.

BACKGROUND: The Mediterranean has a long history of interactions among different peoples. In this study, we investigate the genetic relationships among thirteen population samples from the broader Mediterranean region together with three other groups from the Ivory Coast and Bolivia with a particular focus on the genetic structure between North Africa and South Europe. Analyses were carried out on a diverse set of neutral and functional polymorphisms located in and around the coagulation factor VII and XII genomic regions (F7 and F12). RESULTS: Principal component analysis revealed a significant clustering of the Mediterranean samples into North African and South European groups consistent with the results from the hierarchical AMOVA, which showed a low but significant differentiation between groups from the two shores. For the same range of geographic distances, populations from each side of the Mediterranean were found to differ genetically more than populations within the same side. To further investigate this differentiation, we carried out haplotype analyses, which provided partial evidence that sub-Saharan gene flow was higher towards North Africa than South Europe. CONCLUSIONS: As there is no consensus between the two genomic regions regarding gene flow through the Sahara, it is hard to reach a solid conclusion about its role in the differentiation between the two Mediterranean shores and more data are necessary to reach a definite conclusion. However our data suggest that the Mediterranean Sea was at least partially a barrier to gene flow between the two shores.

Different evolutionary histories of the coagulation factor VII gene in human populations?

Immoderate blood clotting constitutes a risk factor for cardiovascular disease in modern industrialised societies, but is believed to have conferred a survival advantage, i.e. faster recovery from bleeding, on our ancestors. Here, we investigate the evolutionary history of the Coagulation Factor VII gene (F7) by analysing five cardiovascular-risk-associated mutations from the F7 promoter and nine neutral polymorphisms (six SNPs and three microsatellites) from the flanking region in 16 populations from the broader Mediterranean region, South Saharan Africa and Bolivia (687 individuals in total). Population differentiation and selection tests were performed and linkage disequilibrium patterns were investigated. In all samples, no linkage disequilibrium between adjacent F7 promoter mutations -402 and -401 was observed. No selection signals were detected in any of the samples from the broader Mediterranean region and South Saharan Africa, while some of the data suggested a potential signal of positive selection for the F7 promoter in the Native American samples from Bolivia. In conclusion, our data suggest, although do not prove, different evolutionary histories in the F7 promoter region between Mediterraneans and Amerindians.

Mixed origin of the current Tunisian population from the analysis of Alu and Alu/STR compound systems.

During successive historical periods, Tunisia has been a crossroads of multiple civilizations and their corresponding key population movements. The aim of this study was to provide genetic information relating to the mixed origin of the Tunisian population, and to analyze its genetic relationship with other North African and Mediterranean populations. A set of 16 Alu and 3 Alu/STR compound systems has been analyzed in 268 autochthonous Tunisians from the north-center and the south of the country. Our two sampled populations showed no significant differentiation from one another in any of the three Alu/STR compound systems, whereas the analysis of the 16 Alu markers revealed a significant genetic differentiation between them. A sub-Saharan component shown by the three Alu/STR combinations is more noticeable in our north-center sample than in that of the south. The presence of two Alu/STR combinations specific to North African ancestral populations also suggests that the ancient Berber component is relatively more substantial in the north and center regions than in the south. Our Tunisian samples cluster together with other Berber samples from Morocco and Algeria, underpinning the genetic similarity among North Africans regardless of their current linguistic status (Berber or Arabic).