Homozygosity mapping in an Irish ALS case–control cohort describes local demographic phenomena and points towards potential recessive risk loci

Associations of genome-wide and regional autozygosity with 96 complex traits in old order Amish

BACKGROUND

Autozygosity, the proportion of the genome that is homozygous by descent, has been associated with variation in multiple health-related traits impacting evolutionary fitness. Autozygosity (FROH) is typically measured from runs of homozygosity (ROHs) that arise when identical-by-descent (IBD) haplotypes are inherited from each parent. Population isolates with a small set of common founders have elevated autozygosity relative to outbred populations.

METHODS

In this study, we examined whether degree of autozygosity was associated with variation in 96 cardiometabolic traits among 7221 Old Order Amish individuals residing in Lancaster County, PA. We estimated the average length of an ROH segment to be 6350 KB, with each individual having on average 17.2 segments 1.5 KB or larger. Measurements of genome-wide and regional FROH were used as the primary predictors of trait variation in association analysis.

RESULTS

In genome-wide FROH analysis, we did not identify any associations that withstood Bonferroni-correction (p = 0.0005). However, on regional FROH analysis, we identified associations exceeding genome-wide thresholds for two traits: serum bilirubin levels, which were significantly associated with a region on chromosome 2 localized to a region surrounding UGT1A10 (p = 1 × 10- 43), and HbA1c levels, which were significantly associated with a region on chromosome 8 localized near CHRNB3 (p = 8 × 10- 10).

CONCLUSIONS

These analyses highlight the potential value of autozygosity mapping in founder populations.

Genome-wide analysis of runs of homozygosity in Pakistani controls with no history of speech or language-related developmental phenotypes

BACKGROUND

Runs of homozygosity (ROHs) analysis of controls provide a convenient resource to minimize the association of false positive results of disease-associated ROHs and genetic variants for simple and complex disorders in individuals from the same population. Evidence for the value of ROHs to speech or language-related traits is restricted due to the absence of population-matched behaviourally defined controls and limited family-based studies.

AIM

This study aims to identify common ROHs in the Pakistani population, focussing on the total length and frequency of ROHs of variable sizes, shared ROHs, and their genomic distribution.

SUBJECTS AND METHODS

We performed homozygosity analysis (in PLINK) of 86 individuals (39 males, 47 females) with no history of speech or language-related phenotypes (controls) who had been genotyped with the Illumina Infinium QC Array-24.

RESULTS

ROHs of 1-<4 megabases (Mb) were frequent in unrelated individuals. We observed ROHs over 20 Mb among six individuals. Over 30 percent of the identified ROHs were shared among several individuals, indicating consanguinity's effect on the Pakistani population.

CONCLUSION

Our findings serve as a foundation for family-based genetic studies of consanguineous families with speech or language-related disorders to ultimately narrow the homozygosity regions of interest to identify pathogenic variants.

GARLIC: Genomic Autozygosity Regions Likelihood-based Inference and Classification

Summary

Runs of homozygosity (ROH) are important genomic features that manifest when identical-by-descent haplotypes are inherited from parents. Their length distributions and genomic locations are informative about population history and they are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. Here, we present software implementing a model-based method ( Pemberton et al., 2012 ) for inferring ROH in genome-wide SNP datasets that incorporates population-specific parameters and a genotyping error rate as well as provides a length-based classification module to identify biologically interesting classes of ROH. Using simulations, we evaluate the performance of this method.

Availability and Implementation

GARLIC is written in C ++. Source code and pre-compiled binaries (Windows, OSX and Linux) are hosted on GitHub ( https://github.com/szpiech/garlic ) under the GNU General Public License version 3.

Contact

zachary.szpiech@ucsf.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetics of Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating, uniformly lethal degenerative disorder of motor neurons that overlaps clinically with frontotemporal dementia (FTD). Investigations of the 10% of ALS cases that are transmitted as dominant traits have revealed numerous gene mutations and variants that either cause these disorders or influence their clinical phenotype. The evolving understanding of the genetic architecture of ALS has illuminated broad themes in the molecular pathophysiology of both familial and sporadic ALS and FTD. These central themes encompass disturbances of protein homeostasis, alterations in the biology of RNA binding proteins, and defects in cytoskeletal dynamics, as well as numerous downstream pathophysiological events. Together, these findings from ALS genetics provide new insight into therapies that target genetically distinct subsets of ALS and FTD.

Relationship between Deleterious Variation, Genomic Autozygosity, and Disease Risk: Insights from The 1000 Genomes Project

Genomic regions of autozygosity (ROAs) represent segments of individual genomes that are homozygous for haplotypes inherited identical-by-descent (IBD) from a common ancestor. ROAs are nonuniformly distributed across the genome, and increased ROA levels are a reported risk factor for numerous complex diseases. Previously, we hypothesized that long ROAs are enriched for deleterious homozygotes as a result of young haplotypes with recent deleterious mutations-relatively untouched by purifying selection-being paired IBD as a consequence of recent parental relatedness, a pattern supported by ROA and whole-exome sequence data on 27 individuals. Here, we significantly bolster support for our hypothesis and expand upon our original analyses using ROA and whole-genome sequence data on 2,436 individuals from The 1000 Genomes Project. Considering CADD deleteriousness scores, we reaffirm our previous observation that long ROAs are enriched for damaging homozygotes worldwide. We show that strongly damaging homozygotes experience greater enrichment than weaker damaging homozygotes, while overall enrichment varies appreciably among populations. Mendelian disease genes and those encoding FDA-approved drug targets have significantly increased rates of gain in damaging homozygotes with increasing ROA coverage relative to all other genes. In genes implicated in eight complex phenotypes for which ROA levels have been identified as a risk factor, rates of gain in damaging homozygotes vary across phenotypes and populations but frequently differ significantly from non-disease genes. These findings highlight the potential confounding effects of population background in the assessment of associations between ROA levels and complex disease risk, which might underlie reported inconsistencies in ROA-phenotype associations.

Insular Celtic population structure and genomic footprints of migration

Previous studies of the genetic landscape of Ireland have suggested homogeneity, with population substructure undetectable using single-marker methods. Here we have harnessed the haplotype-based method fineSTRUCTURE in an Irish genome-wide SNP dataset, identifying 23 discrete genetic clusters which segregate with geographical provenance. Cluster diversity is pronounced in the west of Ireland but reduced in the east where older structure has been eroded by historical migrations. Accordingly, when populations from the neighbouring island of Britain are included, a west-east cline of Celtic-British ancestry is revealed along with a particularly striking correlation between haplotypes and geography across both islands. A strong relationship is revealed between subsets of Northern Irish and Scottish populations, where discordant genetic and geographic affinities reflect major migrations in recent centuries. Additionally, Irish genetic proximity of all Scottish samples likely reflects older strata of communication across the narrowest inter-island crossing. Using GLOBETROTTER we detected Irish admixture signals from Britain and Europe and estimated dates for events consistent with the historical migrations of the Norse-Vikings, the Anglo-Normans and the British Plantations. The influence of the former is greater than previously estimated from Y chromosome haplotypes. In all, we paint a new picture of the genetic landscape of Ireland, revealing structure which should be considered in the design of studies examining rare genetic variation and its association with traits.

A recent genetic study of the UK (People of the British Isles; PoBI) expanded our understanding of population history of the islands, using newly-developed, powerful techniques that harness the rich information embedded in chunks of genetic code called haplotypes. These methods revealed subtle regional diversity across the UK, and, using genetic data alone, timed key migration events into southeast England and Orkney. We have extended these methods to Ireland, identifying regional differences in genetics across the island that adhere to geography at a resolution not previously reported. Our study reveals relative western diversity and eastern homogeneity in Ireland owing to a history of settlement concentrated on the east coast and longstanding Celtic diversity in the west. We show that Irish Celtic diversity enriches the findings of PoBI; haplotypes mirror geography across Britain and Ireland, with relic Celtic populations contributing greatly to haplotypic diversity. Finally, we used genetic information to date migrations into Ireland from Europe and Britain consistent with historical records of Viking and Norman invasions, demonstrating the signatures of these migrations the on modern Irish genome. Our findings demonstrate that genetic structure exists in even small isolated populations, which has important implications for population-based genetic association studies.

Runs of homozygosity: windows into population history and trait architecture

Long runs of homozygosity (ROH) arise when identical haplotypes are inherited from each parent and thus a long tract of genotypes is homozygous. Cousin marriage or inbreeding gives rise to such autozygosity; however, genome-wide data reveal that ROH are universally common in human genomes even among outbred individuals. The number and length of ROH reflect individual demographic history, while the homozygosity burden can be used to investigate the genetic architecture of complex disease. We discuss how to identify ROH in genome-wide microarray and sequence data, their distribution in human populations and their application to the understanding of inbreeding depression and disease risk.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.

Environmental insults: critical triggers for amyotrophic lateral sclerosis

Background

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by a rapid loss of lower and upper motor neurons. As a complex disease, the ageing process and complicated gene-environment interactions are involved in the majority of cases.

Main body

Significant advances have been made in unravelling the genetic susceptibility to ALS with massively parallel sequencing technologies, while environmental insults remain a suspected but largely unexplored source of risk. Several studies applying the strategy of Mendelian randomisation have strengthened the link between environmental insults and ALS, but none so far has proved conclusive. We propose a new ALS model which links the current knowledge of genetic factors, ageing and environmental insults. This model provides a mechanism as to how ALS is initiated, with environmental insults playing a critical role.

Conclusion

The available evidence has suggested that inherited defect(s) could cause mitochondrial dysfunction, which would establish the primary susceptibility to ALS. Further study of the underlying mechanism may shed light on ALS pathogenesis. Environmental insults are a critical trigger for ALS, particularly in the aged individuals with other toxicant susceptible genes. The identification of ALS triggers could lead to preventive strategies for those individuals at risk.

Genetics of consanguinity and inbreeding in health and disease

CONTEXT

Inbreeding increases the level of homozygotes for autosomal recessive disorders and is the major objective in clinical studies. The prevalence of consanguinity and the degree of inbreeding vary from one population to another depending on ethnicity, religion, culture and geography. Global epidemiological studies have revealed that consanguineous unions have been significantly associated with increased susceptibility to various forms of inherited diseases.

OBJECTIVE

The study aimed to determine the role of consanguinity in human health and to highlight the associated risks for various diseases or disorders.

METHODS

PubMed and Google Scholar search engines were used to explore the published literature on consanguinity and its associated risks using the key words "consanguinity", "prevalence", "inbreeding depression", "coefficient of inbreeding", "child health", "mortality", "human health", "homozygosity" and "complex diseases" in different combinations. The studies were screened for eligibility on the basis of their epidemiological relevance.

RESULTS

This comprehensive assessment highlights the deleterious consequences in populations with a higher prevalence of consanguinity among different countries worldwide.

CONCLUSIONS

To avoid the inbreeding load there is the need to improve socioeconomic and educational status and to increase public awareness of reproductive health and anticipated deleterious effects. Pre-marital and pre-conception counselling of consanguineous populations should be an integral part of health policy to train people and make people aware of its harmful consequences. Furthermore, runs of homozygosity (ROH) and whole-exome sequencing (WES) are useful tools in exploring new genomic signatures for the cause of inbreeding depression.

Association of Long Runs of Homozygosity With Alzheimer Disease Among African American Individuals

IMPORTANCE

Mutations in known causal Alzheimer disease (AD) genes account for only 1% to 3% of patients and almost all are dominantly inherited. Recessive inheritance of complex phenotypes can be linked to long (>1-megabase [Mb]) runs of homozygosity (ROHs) detectable by single-nucleotide polymorphism (SNP) arrays.

OBJECTIVE

To evaluate the association between ROHs and AD in an African American population known to have a risk for AD up to 3 times higher than white individuals.

DESIGN, SETTING, AND PARTICIPANTS

Case-control study of a large African American data set previously genotyped on different genome-wide SNP arrays conducted from December 2013 to January 2015. Global and locus-based ROH measurements were analyzed using raw or imputed genotype data. We studied the raw genotypes from 2 case-control subsets grouped based on SNP array: Alzheimer's Disease Genetics Consortium data set (871 cases and 1620 control individuals) and Chicago Health and Aging Project-Indianapolis Ibadan Dementia Study data set (279 cases and 1367 control individuals). We then examined the entire data set using imputed genotypes from 1917 cases and 3858 control individuals.

MAIN OUTCOMES AND MEASURES

The ROHs larger than 1 Mb, 2 Mb, or 3 Mb were investigated separately for global burden evaluation, consensus regions, and gene-based analyses.

RESULTS

The African American cohort had a low degree of inbreeding (F ~ 0.006). In the Alzheimer's Disease Genetics Consortium data set, we detected a significantly higher proportion of cases with ROHs greater than 2 Mb (P = .004) or greater than 3 Mb (P = .02), as well as a significant 114-kilobase consensus region on chr4q31.3 (empirical P value 2 = .04; ROHs >2 Mb). In the Chicago Health and Aging Project-Indianapolis Ibadan Dementia Study data set, we identified a significant 202-kilobase consensus region on Chr15q24.1 (empirical P value 2 = .02; ROHs >1 Mb) and a cluster of 13 significant genes on Chr3p21.31 (empirical P value 2 = .03; ROHs >3 Mb). A total of 43 of 49 nominally significant genes common for both data sets also mapped to Chr3p21.31. Analyses of imputed SNP data from the entire data set confirmed the association of AD with global ROH measurements (12.38 ROHs >1 Mb in cases vs 12.11 in controls; 2.986 Mb average size of ROHs >2 Mb in cases vs 2.889 Mb in controls; and 22% of cases with ROHs >3 Mb vs 19% of controls) and a gene-cluster on Chr3p21.31 (empirical P value 2 = .006-.04; ROHs >3 Mb). Also, we detected a significant association between AD and CLDN17 (empirical P value 2 = .01; ROHs >1 Mb), encoding a protein from the Claudin family, members of which were previously suggested as AD biomarkers.

CONCLUSIONS AND RELEVANCE

To our knowledge, we discovered the first evidence of increased burden of ROHs among patients with AD from an outbred African American population, which could reflect either the cumulative effect of multiple ROHs to AD or the contribution of specific loci harboring recessive mutations and risk haplotypes in a subset of patients. Sequencing is required to uncover AD variants in these individuals.