Relationships between estimated autozygosity and complex traits in the UK Biobank

Genotypic and phenotypic consequences of domestication in dogs

Runs of homozygosity (ROH) are genomic regions that arise when two copies of an identical ancestral haplotype are inherited from parents with a recent common ancestor. In this study, we performed a novel comprehensive analysis to infer genetic diversity among dogs and quantified the association between ROH and non-disease phenotypes. We found distinct patterns of genetic diversity across clades of breed dogs and elevated levels of long ROH, compared to non- domesticated dogs. These high levels of F ROH (inbreeding coefficient) are a consequence of recent inbreeding among domesticated dogs during breed establishment. We identified statistically significant associations between F ROH and height, weight, lifespan, muscled, white head, white chest, furnish, and length of fur. After correcting for population structure, we identified more than 45 genes across the three examined quantitative traits that exceeded the threshold for suggestive significance, indicating significant polygenic inheritance for the complex quantitative phenotypes in dogs.

Long-term health outcomes from inbreeding in a historical Swedish population: longevity, fertility, and impairments

BACKGROUND

Inbreeding, arising from consanguinity between related parents, has been observed to impact the health of individuals, typically attributed to biological factors. Nevertheless, these effects may be influenced by the social and environmental conditions. The prevalence of consanguineous marriages increased in certain parts of Sweden after it became legal in 1844, which offers a unique opportunity to study and understand the effects of inbreeding on health.

AIM

The objective of this study is to explore the potential impact of inbreeding on the longevity, fertility, and impairments of individuals born in the Skellefteå region, Sweden, between 1890 and 1905, with a follow-up period extending until 1950.

SUBJECTS AND METHODS

The level of inbreeding is calculated using micro-level parish register data and related to longevity, fertility, and impairments using regression analysis.

RESULTS

Inbreeding is shown to be associated with longevity, fertility, and impairments. It seems to affect the risk of stillbirth and impairments and male longevity and fertility.

CONCLUSION

Inbreeding seems to have had a detrimental effect on some health outcomes in this historical population under study.

Influence of autozygosity on common disease risk across the phenotypic spectrum

Autozygosity is associated with rare Mendelian disorders and clinically relevant quantitative traits. We investigated associations between the fraction of the genome in runs of homozygosity (FROH) and common diseases in Genes & Health (n = 23,978 British South Asians), UK Biobank (n = 397,184), and 23andMe. We show that restricting analysis to offspring of first cousins is an effective way of reducing confounding due to social/environmental correlates of FROH. Within this group in G&H+UK Biobank, we found experiment-wide significant associations between FROH and twelve common diseases. We replicated associations with type 2 diabetes (T2D) and post-traumatic stress disorder via within-sibling analysis in 23andMe (median n = 480,282). We estimated that autozygosity due to consanguinity accounts for 5%-18% of T2D cases among British Pakistanis. Our work highlights the possibility of widespread non-additive genetic effects on common diseases and has important implications for global populations with high rates of consanguinity.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

The evolution of personality disorders: A review of proposals

Personality disorders (PDs) are currently considered dysfunctions. However, personality differences are older than humanity and are ubiquitous in nature, from insects to higher primates. This suggests that a number of evolutionary mechanisms-other than dysfunctions-may be able to maintain stable behavioral variation in the gene pool. First of all, apparently maladaptive traits may actually improve fitness by enabling better survival or successful mating or reproduction, as exemplified by neuroticism, psychopathy, and narcissism. Furthermore, some PDs may harm important biological goals while facilitating others, or may be globally beneficial or detrimental depending on environmental circumstances or body condition. Alternatively, certain traits may form part of life history strategies: Coordinated suites of morphological, physiological and behavioral characters that optimize fitness through alternative routes and respond to selection as a whole. Still others may be vestigial adaptations that are no longer beneficial in present times. Finally, variation may be adaptative in and by itself, as it reduces competition for finite resources. These and other evolutionary mechanisms are reviewed and illustrated through human and non-human examples. Evolutionary theory is the best-substantiated explanatory framework across the life sciences, and may shed light on the question of why harmful personalities exist at all.

Exploring the Relationships Between Autozygosity, Educational Attainment, and Cognitive Ability in a Contemporary, Trans-Ancestral American Sample

Previous studies have found significant associations between estimated autozygosity - the proportion of an individual's genome contained in homozygous segments due to distant inbreeding - and multiple traits, including educational attainment (EA) and cognitive ability. In one study, estimated autozygosity showed a stronger association with parental EA than the subject's own EA. This was likely driven by parental EA's association with mobility: more educated parents tended to migrate further from their hometown, and because of the strong correlation between ancestry and geography in the Netherlands, these individuals chose partners farther from their ancestry and therefore more different from them genetically. We examined the associations between estimated autozygosity, cognitive ability, and parental EA in a contemporary sub-sample of adolescents from the Adolescent Brain Cognitive Development Study℠ (ABCD Study®) (analytic N = 6,504). We found a negative association between autozygosity and child cognitive ability consistent with previous studies, while the associations between autozygosity and parental EA were in the expected direction of effect (with greater levels of autozygosity being associated with lower EA) but the effect sizes were significantly weaker than those estimated in previous work. We also found a lower mean level of autozygosity in the ABCD sample compared to previous autozygosity studies, which may reflect overall decreasing levels of autozygosity over generations. Variation in spousal similarities in ancestral background in the ABCD study compared to other studies may explain the pattern of associations between estimated autozygosity, EA, and cognitive ability in the current study.

Parent-offspring inference in inbred populations

Genealogical relationships are fundamental components of genetic studies. However, it is often challenging to infer correct and complete pedigrees even when genome-wide information is available. For example, inbreeding can obscure genetic differences between individuals, making it difficult to even distinguish first-degree relatives such as parent-offspring from full siblings. Similarly, genotyping errors can interfere with the detection of genetic similarity between parents and their offspring. Inbreeding is common in natural, domesticated, and experimental populations and genotyping of these populations often has more errors than in human data sets, so efficient methods for building pedigrees under these conditions are necessary. Here, we present a new method for parent-offspring inference in inbred pedigrees called specific parent-offspring relationship estimation (spore). spore is vastly superior to existing pedigree-inference methods at detecting parent-offspring relationships, in particular when inbreeding is high or in the presence of genotyping errors, or both. spore therefore fills an important void in the arsenal of pedigree inference tools.

Limiting distribution of X-chromosomal coalescence times under first-cousin consanguineous mating

By providing additional opportunities for coalescence within families, the presence of consanguineous unions in a population reduces coalescence times relative to non-consanguineous populations. First-cousin consanguinity can take one of six forms differing in the configuration of sexes in the pedigree of the male and female cousins who join in a consanguineous union: patrilateral parallel, patrilateral cross, matrilateral parallel, matrilateral cross, bilateral parallel, and bilateral cross. Considering populations with each of the six types of first-cousin consanguinity individually and a population with a mixture of the four unilateral types, we examine coalescent models of consanguinity. We previously computed, for first-cousin consanguinity models, the mean coalescence time for X-chromosomal loci and the limiting distribution of coalescence times for autosomal loci. Here, we use the separation-of-time-scales approach to obtain the limiting distribution of coalescence times for X-chromosomal loci. This limiting distribution has an instantaneous coalescence probability that depends on the probability that a union is consanguineous; lineages that do not coalesce instantaneously coalesce according to an exponential distribution. We study the effects on the coalescence time distribution of the type of first-cousin consanguinity, showing that patrilateral-parallel and patrilateral-cross consanguinity have no effect on X-chromosomal coalescence time distributions and that matrilateral-parallel consanguinity decreases coalescence times to a greater extent than does matrilateral-cross consanguinity.

Ancestry-related distribution of Runs of homozygosity and functional variants in Qatari population

Background

Describing how genetic history shapes the pattern of medically relevant variants could improve the understanding of how specific loci interact with each other and affect diseases and traits prevalence. The Qatari population is characterized by a complex history of admixture and substructure, and the study of its population genomic features would provide valuable insights into the genetic landscape of functional variants. Here, we analyzed the genomic variation of 186 newly-genotyped healthy individuals from the Qatari peninsula.

Results

We discovered an intricate genetic structure using ancestry related analyses. In particular, the presence of three different clusters, Cluster 1, Cluster 2 and Cluster 3 (with Near Eastern, South Asian and African ancestry, respectively), was detected with an additional fourth one (Cluster 4) with East Asian ancestry. These subpopulations show differences in the distribution of runs of homozygosity (ROH) and admixture events in the past, ranging from 40 to 5 generations ago. This complex genetic history led to a peculiar pattern of functional markers under positive selection, differentiated in shared signals and private signals. Interestingly we found several signatures of shared selection on SNPs in the FADS2 gene, hinting at a possible common evolutionary link to dietary intake. Among the private signals, we found enrichment for markers associated with HDL and LDL for Cluster 1(Near Eastern ancestry) and Cluster 3 (South Asian ancestry) and height and blood traits for Cluster 2 (African ancestry).

The differences in genetic history among these populations also resulted in the different frequency distribution of putative loss of function variants. For example, homozygous carriers for rs2884737, a variant linked to an anticoagulant drug (warfarin) response, are mainly represented by individuals with predominant Bedouin ancestry (risk allele frequency G at 0.48).

Conclusions

We provided a detailed catalogue of the different ancestral pattern in the Qatari population highlighting differences and similarities in the distribution of selected variants and putative loss of functions. Finally, these results would provide useful guidance for assessing genetic risk factors linked to consanguinity and genetic ancestry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01087-1.

British surname origins, population structure and health outcomes-an observational study of hospital admissions

Population structure is a confounder on pathways linking genotypes to health outcomes. This study examines whether the historical, geographical origins of British surnames are associated with health outcomes today. We coded hospital admissions of over 30 million patients in England between 1999 and 2013 to their British surname origin and divided their diagnoses into 125 major disease categories (of which 94 were complete-case). A base population was constructed with patients’ first admission of any kind. Age- and sex-standardised odds ratios were calculated with logistic regression using patients with ubiquitous English surnames such as “Smith” as reference (alpha = .05; Benjamini–Hochberg false discovery rate (FDR) = .05). The results were scanned for “signals”, where a branch of related surname origins all had significantly higher or lower risk. Age- and sex-standardised admission (alpha = .05) was calculated for each signal across area deprivation and surname origin density quintiles. Signals included three branches of English surnames (disorders of teeth and jaw, fractures, upper gastrointestinal disorders). Although the signal with fractures was considered unusual overall, 2 out of the 9 origins in the branch would only be significant at a FDR > .05: OR 0.92 (95% confidence interval 0.86–0.98) and 0.70 (0.55–0.90). The risk was only different in the quintile with the highest density of that group. Differential risk remained when studied across quintiles of area deprivation. The study shows that surname origins are associated with diverse health outcomes and thus act as markers of population structure over and above area deprivation.

Correspondence Between Genomic- and Genealogical/Coalescent-Based Inference of Homozygosity by Descent in Large French-Canadian Genealogies

Research on the genetics of complex traits overwhelmingly focuses on the additive effects of genes. Yet, animal studies have shown that non-additive effects, in particular homozygosity effects, can shape complex traits. Recent investigations in human studies found some significant homozygosity effects. However, most human populations display restricted ranges of homozygosity by descent (HBD), making the identification of homozygosity effects challenging. Founder populations give rise to higher HBD levels. When deep genealogical data are available in a founder population, it is possible to gain information on the time to the most recent common ancestor (MRCA) from whom a chromosomal segment has been transmitted to both parents of an individual and in turn to that individual. This information on the time to MRCA can be combined with the time to MRCA inferred from coalescent models of gene genealogies. HBD can also be estimated from genomic data. The extent to which the genomic HBD measures correspond to the genealogical/coalescent measures has not been documented in founder populations with extensive genealogical data. In this study, we used simulations to relate genomic and genealogical/coalescent HBD measures. We based our simulations on genealogical data from two ongoing studies from the French-Canadian founder population displaying different levels of inbreeding. We simulated single-nucleotide polymorphisms (SNPs) in a 1-Mb genomic segment from a coalescent model in conjunction with the observed genealogical data. We compared genealogical/coalescent HBD to two genomic methods of HBD estimation based on hidden Markov models (HMMs). We found that genomic estimates of HBD correlated well with genealogical/coalescent HBD measures in both study genealogies. We described generation time to coalescence in terms of genomic HBD estimates and found a large variability in generation time captured by genomic HBD when considering each SNP. However, SNPs in longer segments were more likely to capture recent time to coalescence, as expected. Our study suggests that estimating the coalescent gene genealogy from the genomic data to use in conjunction with observed genealogical data could provide valuable information on HBD.

Polygenic burden could explain high rates of affective disorders in a community with restricted founder population

This study investigates if genetic factors could contribute to the high rate of mood disorders reported in a U.S. community known to have a restricted early founder population (confirmed here through runs of homozygosity analysis). Polygenic scores (PGSs) for eight common diseases, disorders, or traits, including psychiatric disorders, were calculated in 274 participants (125 mood disorder cases) who each reported three or four grandparents born in the community. Ancestry-matched controls were selected from the UK Biobank (UKB; three sets of N = 1,822 each). The mean PGSs were significantly higher in the community for major depression PRS (p = 2.1 × 10^-19 , 0.56 SD units), bipolar disorder (p = 2.5 × 10^-15 , 0.56 SD units), and schizophrenia (p = 3.8 × 10^-21 , 0.64 SD units). The PGSs were not significantly different between the community participants and UKB controls for the traits of body mass index, Type 2 diabetes, coronary artery disease, and chronotype. The mean PGSs for height were significantly lower in the community sample compared to controls (-0.21 SD units, p = 1.2 × 10^-5 ). The results are consistent with enrichment of polygenic risk factors for psychiatric disorders in this community.

Genetic Diversity and Identification of Homozygosity-Rich Genomic Regions in Seven Italian Heritage Turkey (Meleagris gallopavo) Breeds

Italian autochthonous turkey breeds are an important reservoir of genetic biodiversity that should be maintained with an in vivo approach. The aim of this study, part of the TuBAvI national project on biodiversity, was to use run of homozygosity (ROH), together with others statistical approaches (e.g., Wright's F-statistics, principal component analysis, ADMIXTURE analysis), to investigate the genomic diversity in several heritage turkey breeds. We performed a genome-wide characterization of ROH-rich regions in seven autochthonous turkey breeds, i.e., Brianzolo (Brzl), Bronzato Comune Italiano (BrCI), Bronzato dei Colli Euganei (CoEu), Parma e Piacenza (PrPc), Nero d'Italia (NeIt), Ermellinato di Rovigo (ErRo) and Romagnolo (Roma). ROHs were detected based on a 650K SNP genotyping. ROH_islands were identified as homozygous ROH regions shared by at least 75% of birds (within breed). Annotation of genes was performed with DAVID. The admixture analyses revealed that six breeds are unique populations while the Roma breed consists in an admixture of founder populations. Effective population size estimated on genomic data shows a numeric contraction. ROH_islands harbour genes that may be interesting for target selection in commercial populations also. Among them the PTGS2 and PLA2G4A genes on chr10 were related to reproduction efficiency. This is the first study mapping genetic variation in autochthonous turkey populations. Breeds were genetically different among them, with the Roma breed proving to be a mixture of the other breeds. The ROH_islands identified harboured genes peculiar to the selection that occurred in heritage breeds. Finally, this study releases previously undisclosed information on existing genetic variation in the turkey species.

Genomic partitioning of inbreeding depression in humans

Across species, offspring of related individuals often exhibit significant reduction in fitness-related traits, known as inbreeding depression (ID), yet the genetic and molecular basis for ID remains elusive. Here, we develop a method to quantify enrichment of ID within specific genomic annotations and apply it to human data. We analyzed the phenomes and genomes of ∼350,000 unrelated participants of the UK Biobank and found, on average of over 11 traits, significant enrichment of ID within genomic regions with high recombination rates (>21-fold; p < 10^-5), with conserved function across species (>19-fold; p < 10^-4), and within regulatory elements such as DNase I hypersensitive sites (∼5-fold; p = 8.9 × 10^-7). We also quantified enrichment of ID within trait-associated regions and found suggestive evidence that genomic regions contributing to additive genetic variance in the population are enriched for ID signal. We find strong correlations between functional enrichment of SNP-based heritability and that of ID (r = 0.8, standard error: 0.1). These findings provide empirical evidence that ID is most likely due to many partially recessive deleterious alleles in low linkage disequilibrium regions of the genome. Our study suggests that functional characterization of ID may further elucidate the genetic architectures and biological mechanisms underlying complex traits and diseases.

Long tracks of homozygosity predict the severity of alcohol use disorders in an American Indian population

Runs of homozygosity (ROH) arise when an individual inherits two copies of the same haplotype segment. While ROH are ubiquitous across human populations, Native populations-with shared parental ancestry arising from isolation and endogamy-can carry a substantial enrichment for ROH. We have been investigating genetic and environmental risk factors for alcohol use disorders (AUD) in a group of American Indians (AI) who have higher rates of AUD than the general U. S. population. Here we explore whether ROH might be associated with incidence and severity of AUD in this admixed AI population (n = 742) that live on geographically contiguous reservations, using low-coverage whole genome sequences. We have found that the genomic regions in the ROH that were identified in this population had significantly elevated American Indian heritage compared with the rest of the genome. Increased ROH abundance and ROH burden are likely risk factors for AUD severity in this AI population, especially in those diagnosed with severe and moderate AUD. The association between ROH and AUD was mostly driven by ROH of moderate lengths between 1 and 2 Mb. An ROH island on chromosome 1p32.3 and a rare ROH pool on chromosome 3p12.3 were found to be significantly associated with AUD severity. They contain genes involved in lipid metabolism, oxidative stress and inflammatory responses; and OSBPL9 was found to reside on the consensus part of the ROH island. These data demonstrate that ROH are associated with risk for AUD severity in this AI population.

Autozygosity mapping and time-to-spontaneous delivery in Norwegian parent-offspring trios

Parental genetic relatedness may lead to adverse health and fitness outcomes in the offspring. However, the degree to which it affects human delivery timing is unknown. We use genotype data from ≃25 000 parent-offspring trios from the Norwegian Mother, Father and Child Cohort Study to optimize runs of homozygosity (ROH) calling by maximizing the correlation between parental genetic relatedness and offspring ROHs. We then estimate the effect of maternal, paternal and fetal autozygosity and that of autozygosity mapping (common segments and gene burden test) on the timing of spontaneous onset of delivery. The correlation between offspring ROH using a variety of parameters and parental genetic relatedness ranged between −0.2 and 0.6, revealing the importance of the minimum number of genetic variants included in an ROH and the use of genetic distance. The optimized compared to predefined parameters showed a ≃45% higher correlation between parental genetic relatedness and offspring ROH. We found no evidence of an effect of maternal, paternal nor fetal overall autozygosity on spontaneous delivery timing. Yet, through autozygosity mapping, we identified three maternal loci TBC1D1, SIGLECs and EDN1 gene regions reducing the median time-to-spontaneous onset of delivery by ≃2–5% (P-value < 2.3 × 10⁻⁶). We also found suggestive evidence of a fetal locus at 3q22.2, near the RYK gene region (P-value = 2.0 × 10⁻⁶). Autozygosity mapping may provide new insights on the genetic determinants of delivery timing beyond traditional genome-wide association studies, but particular and rigorous attention should be given to ROH calling parameter selection.

Autozygosity influences cardiometabolic disease-associated traits in the AWI-Gen sub-Saharan African study

The analysis of the effects of autozygosity, measured as the change of the mean value of a trait among offspring of genetic relatives, reveals the existence of directional dominance or overdominance. In this study we detect evidence of the effect of autozygosity in 4 out of 13 cardiometabolic disease-associated traits using data from more than 10,000 sub-Saharan African individuals recruited from Ghana, Burkina Faso, Kenya and South Africa. The effect of autozygosity on these phenotypes is found to be sex-related, with inbreeding having a significant decreasing effect in men but a significant increasing effect in women for several traits (body mass index, subcutaneous adipose tissue, low-density lipoproteins and total cholesterol levels). Overall, the effect of inbreeding depression is more intense in men. Differential effects of inbreeding depression are also observed between study sites with different night-light intensity used as proxy for urban development. These results suggest a directional dominant genetic component mediated by environmental interactions and sex-specific differences in genetic architecture for these traits in the Africa Wits-INDEPTH partnership for Genomic Studies (AWI-Gen) cohort.

Hybrid Versus Autochthonous Turkey Populations: Homozygous Genomic Regions Occurrences Due to Artificial and Natural Selection

Simple Summary

In this study we investigate the genomic differentiation of traditional Mexican turkey breeds and commercial hybrid strains. The analysis aimed to identify the effects of different types of selection on the birds’ genome structure. Mexican turkeys are characterized by an adaptive selection to their specific original environment; on the other hand, commercial hybrid strains are directionally selected to maximize productive traits and to reduce production costs. The Mexican turkeys were grouped in two geographic subpopulations, while high genomic homogeneity was found in hybrid birds. Traditional breeds and commercial strains are clearly differentiated from a genetic point of view. Inbreeding coefficients were here calculated with different approaches. A clear effect of selection for productive traits was recorded.

Abstract

The Mexican turkey population is considered to be the descendant of the original domesticated wild turkey and it is distinct from hybrid strains obtained by the intense artificial selection activity that has occurred during the last 40 years. In this study 30 Mexican turkeys were genomically compared to 38 commercial hybrids using 327,342 SNP markers in order to elucidate the differences in genome variability resulting from different types of selection, i.e., only adaptive for Mexican turkey, and strongly directional for hybrids. Runs of homozygosity (ROH) were detected and the two inbreeding coefficients (F and F_ROH) based on genomic information were calculated. Principal component and admixture analyses revealed two different clusters for Mexican turkeys (MEX_cl_1 and MEX_cl_2) showing genetic differentiation from hybrids (HYB) (F_ST equal 0.168 and 0.167, respectively). A total of 3602 ROH were found in the genome of the all turkeys populations. ROH resulted mainly short in length and the ROH_island identified in HYB (n = 9), MEX_cl_1 (n = 1), and MEX_cl_2 (n = 2) include annotated genes related to production traits: abdominal fat (percentage and weight) and egg characteristics (egg shell color and yolk weight). F and F_ROH resulted correlated to each other only for Mexican populations. Mexican turkey genomic variability allows us to separate the birds into two subgroups according to the geographical origin of samples, while the genomic homogeneity of hybrid birds reflected the strong directional selection occurring in this population.

Unbalanced Sample Size Introduces Spurious Correlations to Genome-Wide Heterozygosity Analyses

Excess of heterozygosity (H) is a widely used measure of genetic diversity of a population. As high-throughput sequencing and genotyping data become readily available, it has been applied to investigating the associations of genome-wide genetic diversity with human diseases and traits. However, these studies often report contradictory results. In this paper, we present a meta-analysis of five whole-exome studies to examine the association of H scores with Alzheimer's disease. We show that the mean H score of a group is not associated with the disease status, but ot is associated with the sample size. Across all five studies, the group with more samples has a significantly lower H score than the group with fewer samples. To remove potential confounders in empirical data sets, we perform computer simulations to create artificial genomes controlled for the number of polymorphic loci, the sample size, and the allele frequency. Analyses of these simulated data confirm the negative correlation between the sample size and the H score. Furthermore, we find that genomes with a large number of rare variants also have inflated H scores. These biases altogether can lead to spurious associations between genetic diversity and the phenotype of interest. Based on these findings, we advocate that studies shall balance the sample sizes when using genome-wide H scores to assess genetic diversities of different populations, which helps improve the reproducibility of future research.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

[Epigenomic variations manifesting as a loss of heterozygosity affecting imprinted genes represent a molecular mechanism of autism spectrum disorders and intellectual disability in children]

AIM

Long continuous stretches of homozygosity (LCSH) are regularly detected in studies using molecular karyotyping (SNP array). Despite this type of variation being able to provide meaningful data on the parents' kinship, uniparental disomy and chromosome rearrangements, LCSH are rarely considered as a possible epigenetic cause of neurodevelopmental disorders. Despite their direct relationship to imprinting, LCSH in imprinted loci have not been considered in terms of pathogenicity. The present work is aimed at studying LCSH in chromosomal regions containing imprinted genes previously associated with disease in children with idiopathic intellectual disability, autism, congenital malformations and/or epilepsy.

MATERIAL AND METHODS

Five hundred and four patients with autism spectrum disorders and intellectual disability were examined.

RESULTS

LCSH affecting imprinted loci associated with various diseases were identified in 40 (7.9%) individuals. Chromosomal region 7q21.3 was affected in twenty three cases, 15q11.2 in twelve, 11p15.5 in five, 7q32.2 in four. Four patients had 2 LCSH affecting imprinted loci. Besides one LCSH in 7q31.33q32.3 (~4 Mbp) region, all LCSH were 1-1.6 Mbp. Clinically, these cases resembled the corresponding imprinting diseases (e.g. Silver-Russell, Beckwith-Wiedemann, Prader-Willi, Angelman syndromes). Parental kinship was identified in 8 cases (1.59%), which were not affected by LCSH at imprinted loci.

CONCLUSION

The present study shows that LCSH affecting chromosomal regions 7q21.3, 7q32.2, 11p15.5 and 15p11.2 occur in about 7.9% of children with intellectual disability, autism, congenital malformations and/or epilepsy. Consequently, this type of epigenetic mutations is obviously common in a group of children with neurodevelopmental disorders. LCSH less than 2.5-10 Mbp are usually ignored in molecular karyotyping (SNP array) studies and, therefore, an important epigenetic cause of intellectual disability, autism or epilepsy with high probability remains without attention.

Genome-Wide Homozygosity Mapping Reveals Genes Associated With Cognitive Ability in Children From Saudi Arabia

Recent studies of the genetic foundations of cognitive ability rely on large samples (in extreme, hundreds of thousands) of individuals from relatively outbred populations of mostly European ancestry. Hypothesizing that the genetic foundation of cognitive ability depends on the broader population-specific genetic context, we performed a genome-wide association study and homozygosity mapping of cognitive ability estimates obtained through latent variable modeling in a sample of 354 children from a consanguineous population of Saudi Arabia. Approximately half of the sample demonstrated significantly elevated homozygosity levels indicative of inbreeding, and among those with elevated levels, homozygosity was negatively associated with cognitive ability. Further homozygosity mapping identified a specific run, inclusive of the GRIA4 gene, that survived corrections for multiple testing for association with cognitive ability. The results suggest that in a consanguineous population, a notable proportion of the variance in cognitive ability in the normal range in children might be regulated by population-specific mechanisms such as patterns of elevated homozygosity. This observation has implications for the field's understanding of the etiological bases of intelligence and its variability around the world.

Extreme inbreeding in a European ancestry sample from the contemporary UK population

In most human societies, there are taboos and laws banning mating between first- and second-degree relatives, but actual prevalence and effects on health and fitness are poorly quantified. Here, we leverage a large observational study of ~450,000 participants of European ancestry from the UK Biobank (UKB) to quantify extreme inbreeding (EI) and its consequences. We use genotyped SNPs to detect large runs of homozygosity (ROH) and call EI when >10% of an individual's genome comprise ROHs. We estimate a prevalence of EI of ~0.03%, i.e., ~1/3652. EI cases have phenotypic means between 0.3 and 0.7 standard deviation below the population mean for 7 traits, including stature and cognitive ability, consistent with inbreeding depression estimated from individuals with low levels of inbreeding. Our study provides DNA-based quantification of the prevalence of EI in a European ancestry sample from the UK and measures its effects on health and fitness traits.