Runs of homozygosity in European populations

Runs of homozygosity implicate autozygosity as a schizophrenia risk factor

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ~17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time.

Cryptic distant relatives are common in both isolated and cosmopolitan genetic samples

Although a few hundred single nucleotide polymorphisms (SNPs) suffice to infer close familial relationships, high density genome-wide SNP data make possible the inference of more distant relationships such as 2^nd to 9^th cousinships. In order to characterize the relationship between genetic similarity and degree of kinship given a timeframe of 100–300 years, we analyzed the sharing of DNA inferred to be identical by descent (IBD) in a subset of individuals from the 23andMe customer database (n = 22,757) and from the Human Genome Diversity Panel (HGDP-CEPH, n = 952). With data from 121 populations, we show that the average amount of DNA shared IBD in most ethnolinguistically-defined populations, for example Native American groups, Finns and Ashkenazi Jews, differs from continentally-defined populations by several orders of magnitude. Via extensive pedigree-based simulations, we determined bounds for predicted degrees of relationship given the amount of genomic IBD sharing in both endogamous and ‘unrelated’ population samples. Using these bounds as a guide, we detected tens of thousands of 2^nd to 9^th degree cousin pairs within a heterogenous set of 5,000 Europeans. The ubiquity of distant relatives, detected via IBD segments, in both ethnolinguistic populations and in large ‘unrelated’ populations samples has important implications for genetic genealogy, forensics and genotype/phenotype mapping studies.

Turkish population structure and genetic ancestry reveal relatedness among Eurasian populations

Turkey has experienced major population movements. Population structure and genetic relatedness of samples from three regions of Turkey, using over 500,000 SNP genotypes, were compared together with Human Genome Diversity Panel (HGDP) data. To obtain a more representative sampling from Central Asia, Kyrgyz samples (Bishkek, Kyrgyzstan) were genotyped and analysed. Principal component (PC) analysis reveals a significant overlap between Turks and Middle Easterners and a relationship with Europeans and South and Central Asians; however, the Turkish genetic structure is unique. FRAPPE, STRUCTURE, and phylogenetic analyses support the PC analysis depending upon the number of parental ancestry components chosen. For example, supervised STRUCTURE (K=3) illustrates a genetic ancestry for the Turks of 45% Middle Eastern (95% CI, 42-49), 40% European (95% CI, 36-44) and 15% Central Asian (95% CI, 13-16), whereas at K=4 the genetic ancestry of the Turks was 38% European (95% CI, 35-42), 35% Middle Eastern (95% CI, 33-38), 18% South Asian (95% CI, 16-19) and 9% Central Asian (95% CI, 7-11). PC analysis and FRAPPE/STRUCTURE results from three regions in Turkey (Aydin, Istanbul and Kayseri) were superimposed, without clear subpopulation structure, suggesting sample homogeneity. Thus, this study demonstrates admixture of Turkish people reflecting the population migration patterns.

Homozygosity mapping identifies a bile acid biosynthetic defect in an adult with cirrhosis of unknown etiology

The most common inborn error of bile acid metabolism is 3β-hydroxy-Δ(5)-C(27)-steroid oxidoreductase (3β-HSD) deficiency, a disorder that usually presents in early childhood with hepatic dysfunction. Timely diagnosis of this disorder is crucial because it can be effectively treated with primary bile acid replacement. Here we describe a 24-year-old woman from Iran with cirrhosis of unknown etiology. Her sister and a first cousin died of cirrhosis (ages 19 and 6 years) and another 32-year-old first cousin had a self-limited liver disorder in childhood that resolved at age 9 years. The family history suggested that the affected family members were homozygous for a mutant allele inherited identical-by-descent. A genome-wide analysis of 2.4 million single nucleotide polymorphisms was performed to identify regions of homozygosity that were present in the proband and the 32-year-old first cousin, but not in a healthy relative. One of these regions contained the gene encoding 3β-HSD (HSD3B7). Sequence analysis of HSD3B7 revealed that the proband and her 32-year-old cousin were homozygous for a frameshift mutation (c.45_46del AG, p.T15Tfsx27) in exon 1. The diagnosis of 3β-HSD deficiency was confirmed by documenting high levels of 3β-hydroxy-Δ(5) bile acids in the serum of the proband and the 32-year-old first cousin using mass spectrometry. To our knowledge, the 32-year-old relative in this family represents the oldest asymptomatic patient with this disorder.

CONCLUSION

This study highlights the clinical utility of homozygosity mapping in diagnosing autosomal recessive metabolic disorders. This family illustrates the wide variation in expressivity that occurs in 3β-HSD deficiency and underscores the need to consider a bile acid synthetic defect as a possible cause of liver disease in adults.

Homozygosity mapping on a single patient: identification of homozygous regions of recent common ancestry by using population data

Homozygosity mapping has played an important role in detecting recessive mutations using families of consanguineous marriages. However, detection of regions identical and homozygosity by descent (HBD) when family data are not available, or when relationships are unknown, is still a challenge. Making use of population data from high-density SNP genotyping may allow detection of regions HBD from recent common founders in singleton patients without genealogy information. We report a novel algorithm that detects such regions by estimating the population haplotype frequencies (HF) for an entire homozygous region. We also developed a simulation method to evaluate the probability of HBD and linkage to disease for a homozygous region by examining the best regions in unaffected controls from the host population. The method can be applied to diseases of Mendelian inheritance but can also be extended to complex diseases to detect rare founder mutations that affect a very small number of patients using either multiplex families or sporadic cases. Testing of the method on both real cases (singleton affected) and simulated data demonstrated its superb sensitivity and robustness under genetic heterogeneity. Hum Mutat 32:345–353, 2011. © 2011 Wiley-Liss, Inc.

Genome features of "Dark-fly", a Drosophila line reared long-term in a dark environment

Organisms are remarkably adapted to diverse environments by specialized metabolisms, morphology, or behaviors. To address the molecular mechanisms underlying environmental adaptation, we have utilized a Drosophila melanogaster line, termed “Dark-fly”, which has been maintained in constant dark conditions for 57 years (1400 generations). We found that Dark-fly exhibited higher fecundity in dark than in light conditions, indicating that Dark-fly possesses some traits advantageous in darkness. Using next-generation sequencing technology, we determined the whole genome sequence of Dark-fly and identified approximately 220,000 single nucleotide polymorphisms (SNPs) and 4,700 insertions or deletions (InDels) in the Dark-fly genome compared to the genome of the Oregon-R-S strain, a control strain. 1.8% of SNPs were classified as non-synonymous SNPs (nsSNPs: i.e., they alter the amino acid sequence of gene products). Among them, we detected 28 nonsense mutations (i.e., they produce a stop codon in the protein sequence) in the Dark-fly genome. These included genes encoding an olfactory receptor and a light receptor. We also searched runs of homozygosity (ROH) regions as putative regions selected during the population history, and found 21 ROH regions in the Dark-fly genome. We identified 241 genes carrying nsSNPs or InDels in the ROH regions. These include a cluster of alpha-esterase genes that are involved in detoxification processes. Furthermore, analysis of structural variants in the Dark-fly genome showed the deletion of a gene related to fatty acid metabolism. Our results revealed unique features of the Dark-fly genome and provided a list of potential candidate genes involved in environmental adaptation.

Cooperative genome-wide analysis shows increased homozygosity in early onset Parkinson's disease

Parkinson's disease (PD) occurs in both familial and sporadic forms, and both monogenic and complex genetic factors have been identified. Early onset PD (EOPD) is particularly associated with autosomal recessive (AR) mutations, and three genes, PARK2, PARK7 and PINK1, have been found to carry mutations leading to AR disease. Since mutations in these genes account for less than 10% of EOPD patients, we hypothesized that further recessive genetic factors are involved in this disorder, which may appear in extended runs of homozygosity.

We carried out genome wide SNP genotyping to look for extended runs of homozygosity (ROHs) in 1,445 EOPD cases and 6,987 controls. Logistic regression analyses showed an increased level of genomic homozygosity in EOPD cases compared to controls. These differences are larger for ROH of 9 Mb and above, where there is a more than three-fold increase in the proportion of cases carrying a ROH. These differences are not explained by occult recessive mutations at existing loci. Controlling for genome wide homozygosity in logistic regression analyses increased the differences between cases and controls, indicating that in EOPD cases ROHs do not simply relate to genome wide measures of inbreeding. Homozygosity at a locus on chromosome19p13.3 was identified as being more common in EOPD cases as compared to controls. Sequencing analysis of genes and predicted transcripts within this locus failed to identify a novel mutation causing EOPD in our cohort.

There is an increased rate of genome wide homozygosity in EOPD, as measured by an increase in ROHs. These ROHs are a signature of inbreeding and do not necessarily harbour disease-causing genetic variants. Although there might be other regions of interest apart from chromosome 19p13.3, we lack the power to detect them with this analysis.

Genome-wide analysis of epistasis in body mass index using multiple human populations

We surveyed gene–gene interactions (epistasis) in human body mass index (BMI) in four European populations (n<1200) via exhaustive pair-wise genome scans where interactions were computed as F ratios by testing a linear regression model fitting two single-nucleotide polymorphisms (SNPs) with interactions against the one without. Before the association tests, BMI was corrected for sex and age, normalised and adjusted for relatedness. Neither single SNPs nor SNP interactions were genome-wide significant in either cohort based on the consensus threshold (P=5.0E−08) and a Bonferroni corrected threshold (P=1.1E−12), respectively. Next we compared sub genome-wide significant SNP interactions (P<5.0E−08) across cohorts to identify common epistatic signals, where SNPs were annotated to genes to test for gene ontology (GO) enrichment. Among the epistatic genes contributing to the commonly enriched GO terms, 19 were shared across study cohorts of which 15 are previously published genome-wide association loci, including CDH13 (cadherin 13) associated with height and SORCS2 (sortilin-related VPS10 domain containing receptor 2) associated with circulating insulin-like growth factor 1 and binding protein 3. Interactions between the 19 shared epistatic genes and those involving BMI candidate loci (P<5.0E−08) were tested across cohorts and found eight replicated at the SNP level (P<0.05) in at least one cohort, which were further tested and showed limited replication in a separate European population (n>5000). We conclude that genome-wide analysis of epistasis in multiple populations is an effective approach to provide new insights into the genetic regulation of BMI but requires additional efforts to confirm the findings.

Association of medication with the human plasma N-glycome

Glycosylation is highly variable depending on many environmental factors. Using our fully quantitative high-throughput normal phase hydrophilic interaction liquid chromatography platform we have identified glycosylation changes associated with medication in the plasma N-glycome from three different population cohorts: ORCADES from the Orkney Islands in Scotland and CROATIA-Vis and CROATIA-Korcula from the Croatian islands of Vis and Korcula. Associations between glycosylation and the use of hormones (oral contraceptives, hormone replacement therapy), nonsteroidal anti-inflammatory drugs (aspirin and other NSAIDs), oral steroids (prednisolone) and steroid inhalers (beclomethasone) were investigated. Significant differences associated with usage of oral contraceptives were found with increased core-fucosylated biantennary glycans. Decreases in core-fucosylated biantennary glycans, core-fucosylated triantennary glycans with outer-arm fucose, and high mannosylated glycans were associated with the use of anti-inflammatory drugs. All of the changes in glycosylation were independent of blood group status. In conclusion, hormones and anti-inflammatory medication were associated with changes in glycosylation, possibly as a result of the modulatory effect of these drugs on the inflammatory response. In general, cancer is associated with inflammation, and many glycoproteins in the plasma are acute phase related to the host response. These preliminary data indicate the importance of correcting the levels of glycans used as biomarkers for the effects of medication.

Consanguinity in Centre d'Étude du Polymorphisme Humain (CEPH) pedigrees

A set of Centre d'Étude du Polymorphisme Humain (CEPH) cell lines serves as a large reference collection that has been widely used as a benchmark for allele frequencies in the analysis of genetic variants, to create linkage maps of the human genome, to study the genetics of gene expression, to provide samples to the HapMap and 1000 Genomes projects, and for a variety of other applications. An explicit feature of the CEPH collection is that these multigenerational families represent reference panels of known relatedness, consisting mostly of three-generation pedigrees with large sibships, two parents, and grandparents. We applied identity-by-state (IBS) and identity-by-descent (IBD) methods to high-density genotype data from 186 CEPH individuals in 13 families. We identified unexpected relatedness between nominally unrelated grandparents both within and between pedigrees. For one pair, the estimated Cotterman coefficient of relatedness k1 exceeded 0.2, consistent with one-eighth sharing (eg, first-cousins). Unexpectedly, significant IBD2 values were discovered in both second-degree and parent-child relationships. These were accompanied by regions of homozygosity in the offspring, which corresponded to blocks lacking IBS0 in purportedly unrelated parents, consistent with inbreeding. Our findings support and extend a 1999 report, based on the use of short tandem-repeat polymorphisms, that several CEPH families had regions of homozygosity consistent with autozygosity. We benchmarked our IBD approach (called kcoeff) against both RELPAIR and PREST software packages. Our findings may affect the interpretation of previous studies and the design of future studies that rely on the CEPH resource.

Ancient human DNA

The contribution of palaeogenetic data to the study of various aspects of hominin biology and evolution has been significant, and has the potential to increase substantially with the widespread implementation of next generation sequencing techniques. Here we discuss the present state-of-the-art of ancient human DNA analysis and the characteristics of hominin aDNA that make sequence validation particularly complex. A brief overview of the development of anthropological palaeogenetic analysis is given to illustrate the technical challenges motivating recent technological advancements.

Genotyping techniques to address diversity in tumors

Array-based genotyping platforms have during recent years been established as a valuable tool for the characterization of genomic alterations in cancer. The analysis of tumor samples, however, presents challenges for data analysis and interpretation. For example, tumor samples are often admixed with nonaberrant cells that define the tumor microenvironment, such as infiltrating lymphocytes and fibroblasts, or vasculature. Furthermore, tumors often comprise subclones harboring divergent aberrations that are acquired subsequent to the tumor-initiating event. The combined analysis of both genotype and copy number status obtained by array-based genotyping platforms provide opportunities to address these challenges. In this chapter, we present the basic principles for current array-based genotyping platforms and how they can be used to infer genotype and copy number for acquired genomic alterations. We describe how these techniques can be used to resolve tumor ploidy, normal cell admixture, and subclonality. We also exemplify how genotyping techniques can be applied in tumor studies to elucidate the hierarchy among tumor clones, and thus, provide means to study clonal expansion and tumor evolution.

Estimating the allele frequency of autosomal recessive disorders through mutational records and consanguinity: the Homozygosity Index (HI)

In principle mutational records make it possible to estimate frequencies of disease alleles (q) for autosomal recessive disorders using a novel approach based on the calculation of the Homozygosity Index (HI), i.e., the proportion of homozygous patients, which is complementary to the proportion of compound heterozygous patients P(CH). In other words, the rarer the disorder, the higher will be the HI and the lower will be the P(CH). To test this hypothesis we used mutational records of individuals affected with Familial Mediterranean Fever (FMF) and Phenylketonuria (PKU), born to either consanguineous or apparently unrelated parents from six population samples of the Mediterranean region. Despite the unavailability of precise values of the inbreeding coefficient for the general population, which are needed in the case of apparently unrelated parents, our estimates of q are very similar to those of previous descriptive epidemiological studies. Finally, we inferred from simulation studies that the minimum sample size needed to use this approach is 25 patients either with unrelated or first cousin parents. These results show that the HI can be used to produce a ranking order of allele frequencies of autosomal recessive disorders, especially in populations with high rates of consanguineous marriages.

Chromosome 3 anomalies investigated by genome wide SNP analysis of benign, low malignant potential and low grade ovarian serous tumours

Ovarian carcinomas exhibit extensive heterogeneity, and their etiology remains unknown. Histological and genetic evidence has led to the proposal that low grade ovarian serous carcinomas (LGOSC) have a different etiology than high grade carcinomas (HGOSC), arising from serous tumours of low malignant potential (LMP). Common regions of chromosome (chr) 3 loss have been observed in all types of serous ovarian tumours, including benign, suggesting that these regions contain genes important in the development of all ovarian serous carcinomas. A high-density genome-wide genotyping bead array technology, which assayed >600,000 markers, was applied to a panel of serous benign and LMP tumours and a small set of LGOSC, to characterize somatic events associated with the most indolent forms of ovarian disease. The genomic patterns inferred were related to TP53, KRAS and BRAF mutations. An increasing frequency of genomic anomalies was observed with pathology of disease: 3/22 (13.6%) benign cases, 40/53 (75.5%) LMP cases and 10/11 (90.9%) LGOSC cases. Low frequencies of chr3 anomalies occurred in all tumour types. Runs of homozygosity were most commonly observed on chr3, with the 3p12-p11 candidate tumour suppressor region the most frequently homozygous region in the genome. An LMP harboured a homozygous deletion on chr6 which created a GOPC-ROS1 fusion gene, previously reported as oncogenic in other cancer types. Somatic TP53, KRAS and BRAF mutations were not observed in benign tumours. KRAS-mutation positive LMP cases displayed significantly more chromosomal aberrations than BRAF-mutation positive or KRAS and BRAF mutation negative cases. Gain of 12p, which harbours the KRAS gene, was particularly evident. A pathology review reclassified all TP53-mutation positive LGOSC cases, some of which acquired a HGOSC status. Taken together, our results support the view that LGOSC could arise from serous benign and LMP tumours, but does not exclude the possibility that HGOSC may derive from LMP tumours.

No evidence that extended tracts of homozygosity are associated with Alzheimer's disease

We sought to investigate the contribution of extended runs of homozygosity in a genome-wide association dataset of 1,955 Alzheimer's disease cases and 955 elderly screened controls genotyped for 529,205 autosomal single nucleotide polymorphisms. Tracts of homozygosity may mark regions inherited from a common ancestor and could reflect disease loci if observed more frequently in cases than controls. We found no excess of homozygous tracts in Alzheimer's disease cases compared to controls and no individual run of homozygosity showed association to Alzheimer's disease.

Regions of homozygosity in three Southeast Asian populations

The genomes of outbred populations were first shown in 2006 to contain regions of homozygosity (ROHs) of several megabases. Further studies have also investigated the characteristics of ROHs in healthy individuals in various populations but there are no studies on Singapore populations to date. This study aims to identify and investigate the characteristics of ROHs in three Singapore populations. A total of 268 samples (96 Chinese, 89 Malays and 83 Indians) are genotyped on Illumina Human 1 M Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We use the PennCNV algorithm to detect ROHs. We report an abundance of ROHs (≥500 kb), with an average of more than one hundred regions per individual. On average, the Indian population has the lowest number of ROHs and smallest total length of ROHs per individual compared with the Chinese and Malay populations. We further investigate the relationship between the occurrence of ROHs and haplotype frequency, regional linkage disequilibrium (LD) and positive selection. Based on the results of this data set, we find that the frequency of occurrence of ROHs is positively associated with haplotype frequency and regional LD. The majority of regions detected for recent positive selection and regions with differential LD between populations overlap with the ROH loci. When we consider both the location of the ROHs and the allelic form of the ROHs, we are able to separate the populations by principal component analysis, demonstrating that ROHs contain information on population structure and the demographic history of a population.

Molecular genetic analysis of retinitis pigmentosa in Indonesia using genome-wide homozygosity mapping

PURPOSE

Retinitis pigmentosa (RP) is a clinically and genetically heterogeneous retinal disorder. Despite tremendous knowledge about the genes involved in RP, little is known about the genetic causes of RP in Indonesia. Here, we aim to identify the molecular genetic causes underlying RP in a small cohort of Indonesian patients, using genome-wide homozygosity mapping.

METHODS

DNA samples from affected and healthy individuals from 14 Indonesian families segregating autosomal recessive, X-linked, or isolated RP were collected. Homozygosity mapping was conducted using Illumina 6k or Affymetrix 5.0 single nucleotide polymorphism (SNP) arrays. Known autosomal recessive RP (arRP) genes residing in homozygous regions and X-linked RP genes were sequenced for mutations.

RESULTS

In ten out of the 14 families, homozygous regions were identified that contained genes known to be involved in the pathogenesis of RP. Sequence analysis of these genes revealed seven novel homozygous mutations in ATP-binding cassette, sub-family A, member 4 (ABCA4), crumbs homolog 1 (CRB1), eyes shut homolog (Drosophila) (EYS), c-mer proto-oncogene tyrosine kinase (MERTK), nuclear receptor subfamily 2, group E, member 3 (NR2E3) and phosphodiesterase 6A, cGMP-specific, rod, alpha (PDE6A), all segregating in the respective families. No mutations were identified in the X-linked genes retinitis pigmentosa GTPase regulator (RPGR) and retinitis pigmentosa 2 (X-linked recessive; RP2).

CONCLUSIONS

Homozygosity mapping is a powerful tool to identify the genetic defects underlying RP in the Indonesian population. Compared to studies involving patients from other populations, the same genes appear to be implicated in the etiology of recessive RP in Indonesia, although all mutations that were discovered are novel and as such may be unique for this population.

Autozygome decoded

Consanguineous unions permit the "reunion" of ancestral chromosomal segments in a pattern referred to as "autozygosity," which is essentially a special form of homozygosity. This has long been exploited as a gene mapping tool because it is relatively easy to track a recessive mutation by the surrounding pattern of homozygous markers. The recent advent of single nucleotide polymorphism microarrays has rapidly replaced the historical use of microsatellites for this purpose. In this review, the author discusses other exciting opportunities offered by this unique arrangement of the human genome that range from pure clinical to functional genomic applications. The emerging field of whole genome sequencing promises to unlock much of the potential of the "autozygome."

Identity-by-descent-based phasing and imputation in founder populations using graphical models

Accurate knowledge of haplotypes, the combination of alleles co-residing on a single copy of a chromosome, enables powerful gene mapping and sequence imputation methods. Since humans are diploid, haplotypes must be derived from genotypes by a phasing process. In this study, we present a new computational model for haplotype phasing based on pairwise sharing of haplotypes inferred to be Identical-By-Descent (IBD). We apply the Bayesian network based model in a new phasing algorithm, called systematic long-range phasing (SLRP), that can capitalize on the close genetic relationships in isolated founder populations, and show with simulated and real genome-wide genotype data that SLRP substantially reduces the rate of phasing errors compared to previous phasing algorithms. Furthermore, the method accurately identifies regions of IBD, enabling linkage-like studies without pedigrees, and can be used to impute most genotypes with very low error rate.

Assessing the influence of consanguinity on congenital heart disease

Numerous articles have been published linking consanguineous marriage to an elevated prevalence of congenital heart disease, with ventricular septal defects and atrial septal defects the most commonly cited disorders. While initially persuasive, on closer examination many of these studies have fundamental shortcomings in their design and in the recruitment of study subjects and controls. Improved matching of cases and controls, to include recognition of the long-established community boundaries within which most marriages are contracted, and the assessment of consanguinity within specific levels and types of marital union would improve and help to focus the study outcomes. At the same time, major discrepancies between studies in their reported prevalence and types of congenital heart disease suggest an urgent need for greater standardization in the classification and reporting of these disorders.

Detecting autozygosity through runs of homozygosity: a comparison of three autozygosity detection algorithms

Background

A central aim for studying runs of homozygosity (ROHs) in genome-wide SNP data is to detect the effects of autozygosity (stretches of the two homologous chromosomes within the same individual that are identical by descent) on phenotypes. However, it is unknown which current ROH detection program, and which set of parameters within a given program, is optimal for differentiating ROHs that are truly autozygous from ROHs that are homozygous at the marker level but vary at unmeasured variants between the markers.

Method

We simulated 120 Mb of sequence data in order to know the true state of autozygosity. We then extracted common variants from this sequence to mimic the properties of SNP platforms and performed ROH analyses using three popular ROH detection programs, PLINK, GERMLINE, and BEAGLE. We varied detection thresholds for each program (e.g., prior probabilities, lengths of ROHs) to understand their effects on detecting known autozygosity.

Results

Within the optimal thresholds for each program, PLINK outperformed GERMLINE and BEAGLE in detecting autozygosity from distant common ancestors. PLINK's sliding window algorithm worked best when using SNP data pruned for linkage disequilibrium (LD).

Conclusion

Our results provide both general and specific recommendations for maximizing autozygosity detection in genome-wide SNP data, and should apply equally well to research on whole-genome autozygosity burden or to research on whether specific autozygous regions are predictive using association mapping methods.

Genome-wide homozygosity and multiple sclerosis in Orkney and Shetland Islanders

There is strong evidence for both genetic and environmental risk factors comprising the aetiology of multiple sclerosis (MS). While much progress has been made in recent years in identifying common genetic variants using genome-wide association studies, alternative approaches have remained relatively neglected. The prevalence of MS in Orkney and Shetland is among the highest in the world. Previous studies have suggested that a higher degree of parental relatedness in these isolated communities may contribute to the high rates of MS, indicating that recessive effects have an important role in MS aetiology. The Northern Isles Multiple Sclerosis (NIMS) study investigated the potential role of genome-wide homozygosity in MS risk by genotyping 88 MS patients, 89 controls matched by age, sex and ancestry, and a further 89 controls matched for sex and ancestry, but passed the majority of lifetime risk of developing MS (>70 years of age). Three participants were removed on the basis of pedigree-genomic anomalies (n=263). Three measures of genome-wide homozygosity were generated for each individual, and association with MS was assessed using logistic regression models. No effect of genome-wide homozygosity was detected, indicating that inbreeding and consanguinity are not risk factors for MS in this population.

SNP array-based whole genome homozygosity mapping as the first step to a molecular diagnosis in patients with Charcot-Marie-Tooth disease

Considerable non-allelic heterogeneity for autosomal recessively inherited Charcot-Marie-Tooth (ARCMT) disease has challenged molecular testing and often requires a large amount of work in terms of DNA sequencing and data interpretation or remains unpractical. This study tested the value of SNP array-based whole-genome homozygosity mapping as a first step in the molecular genetic diagnosis of sporadic or ARCMT in patients from inbred families or outbred populations with the ancestors originating from the same geographic area. Using 10 K 2.0 and 250 K Nsp Affymetrix SNP arrays, 15 (63%) of 24 CMT patients received an accurate genetic diagnosis. We used our Java-based script eHoPASA CMT—easy Homozygosity Profiling of SNP arrays for CMT patients to display the location of homozygous regions and their extent of marker count and base-pairs throughout the whole genome. CMT4C was the most common genetic subtype with mutations detected in SH3TC2, one (p.E632Kfs13X) appearing to be a novel founder mutation. A sporadic patient with severe CMT was homozygous for the c.250G > C (p.G84R) HSPB1 mutation which has previously been reported to cause autosomal dominant dHMN. Two distantly related CMT1 patients with early disease onset were found to carry a novel homozygous mutation in MFN2 (p.N131S). We conclude that SNP array-based homozygosity mapping is a fast, powerful, and economic tool to guide molecular genetic testing in ARCMT and in selected sporadic CMT patients.

EX-HOM (EXome HOMozygosity): a proof of principle

OBJECTIVE

We provide the proof of principle that exome sequencing of only two affected siblings born to first-cousin parents is capable of directly identifying a single candidate gene for an autosomal recessive disorder. This strategy, which we call EX-HOM (EXome HOMozygosity), combines in a single step the capacity of exome sequencing to identify all the coding variants present in a genome with the property of homozygosity mapping to limit the search for candidate genes to specific chromosomal regions.

METHODS

We sequenced the exomes of two siblings born to first-cousin parents affected with dysmyelinating leukodystrophy and spastic paraparesis caused by a mutation in FA2H. We used exome sequencing data to identify homozygous regions shared by the two affected siblings (EX-HOM regions), compared them with the regions of maximum LOD score obtained with SNP genotyping, and selected the candidate variants within.

RESULTS

We identified regions of shared homozygosity (>1 Mb) accounting for about 290 Mb, containing only 3 candidate variants. Among these, the FA2H mutation remained the only plausible one.

CONCLUSION

In single consanguineous pedigrees with a few affected sibs, EX-HOM can be a one-step approach to identify the candidate genetic defect, bypassing obstacles such as genetic heterogeneity and the need for large pedigrees.

Copy number variation across European populations

Genome analysis provides a powerful approach to test for evidence of genetic variation within and between geographical regions and local populations. Copy number variants which comprise insertions, deletions and duplications of genomic sequence provide one such convenient and informative source. Here, we investigate copy number variants from genome wide scans of single nucleotide polymorphisms in three European population isolates, the island of Vis in Croatia, the islands of Orkney in Scotland and the South Tyrol in Italy. We show that whereas the overall copy number variant frequencies are similar between populations, their distribution is highly specific to the population of origin, a finding which is supported by evidence for increased kinship correlation for specific copy number variants within populations.

Sequencing of high-complexity DNA pools for identification of nucleotide and structural variants in regions associated with complex traits

We have used targeted genomic sequencing of high-complexity DNA pools based on long-range PCR and deep DNA sequencing by the SOLiD technology. The method was used for sequencing of 286 kb from four chromosomal regions with quantitative trait loci (QTL) influencing blood plasma lipid and uric acid levels in DNA pools of 500 individuals from each of five European populations. The method shows very good precision in estimating allele frequencies as compared with individual genotyping of SNPs (r(2) = 0.95, P < 10(-16)). Validation shows that the method is able to identify novel SNPs and estimate their frequency in high-complexity DNA pools. In our five populations, 17% of all SNPs and 61% of structural variants are not available in the public databases. A large fraction of the novel variants show a limited geographic distribution, with 62% of the novel SNPs and 59% of novel structural variants being detected in only one of the populations. The large number of population-specific novel SNPs underscores the need for comprehensive sequencing of local populations in order to identify the causal variants of human traits.

Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data

Inbreeding depression, which refers to reduced fitness among offspring of related parents, has traditionally been studied using pedigrees. In practice, pedigree information is difficult to obtain, potentially unreliable, and rarely assessed for inbreeding arising from common ancestors who lived more than a few generations ago. Recently, there has been excitement about using SNP data to estimate inbreeding (F) arising from distant common ancestors in apparently “outbred” populations. Statistical power to detect inbreeding depression using SNP data depends on the actual variation in inbreeding in a population, the accuracy of detecting that with marker data, the effect size, and the sample size. No one has yet investigated what variation in F is expected in SNP data as a function of population size, and it is unclear which estimate of F is optimal for detecting inbreeding depression. In the present study, we use theory, simulated genetic data, and real genetic data to find the optimal estimate of F, to quantify the likely variation in F in populations of various sizes, and to estimate the power to detect inbreeding depression. We find that F estimated from runs of homozygosity (F_roh), which reflects shared ancestry of genetic haplotypes, retains variation in even large populations (e.g., SD = 0.5% when N_e = 10,000) and is likely to be the most powerful method of detecting inbreeding effects from among several alternative estimates of F. However, large samples (e.g., 12,000–65,000) will be required to detect inbreeding depression for likely effect sizes, and so studies using F_roh to date have probably been underpowered.

Population-genetic comparison of the Sorbian isolate population in Germany with the German KORA population using genome-wide SNP arrays

Background

The Sorbs are an ethnic minority in Germany with putative genetic isolation, making the population interesting for disease mapping. A sample of N = 977 Sorbs is currently analysed in several genome-wide meta-analyses. Since genetic differences between populations are a major confounding factor in genetic meta-analyses, we compare the Sorbs with the German outbred population of the KORA F3 study (N = 1644) and other publically available European HapMap populations by population genetic means. We also aim to separate effects of over-sampling of families in the Sorbs sample from effects of genetic isolation and compare the power of genetic association studies between the samples.

Results

The degree of relatedness was significantly higher in the Sorbs. Principal components analysis revealed a west to east clustering of KORA individuals born in Germany, KORA individuals born in Poland or Czech Republic, Half-Sorbs (less than four Sorbian grandparents) and Full-Sorbs. The Sorbs cluster is nearest to the cluster of KORA individuals born in Poland. The number of rare SNPs is significantly higher in the Sorbs sample. FST between KORA and Sorbs is an order of magnitude higher than between different regions in Germany. Compared to the other populations, Sorbs show a higher proportion of individuals with runs of homozygosity between 2.5 Mb and 5 Mb. Linkage disequilibrium (LD) at longer range is also slightly increased but this has no effect on the power of association studies.

Oversampling of families in the Sorbs sample causes detectable bias regarding higher FST values and higher LD but the effect is an order of magnitude smaller than the observed differences between KORA and Sorbs. Relatedness in the Sorbs also influenced the power of uncorrected association analyses.

Conclusions

Sorbs show signs of genetic isolation which cannot be explained by over-sampling of relatives, but the effects are moderate in size. The Slavonic origin of the Sorbs is still genetically detectable.

Regarding LD structure, a clear advantage for genome-wide association studies cannot be deduced. The significant amount of cryptic relatedness in the Sorbs sample results in inflated variances of Beta-estimators which should be considered in genetic association analyses.

Runs of homozygosity do not influence survival to old age

Runs of homozygosity (ROH) are extended tracts of adjacent homozygous single nucleotide polymorphisms (SNPs) that are more common in unrelated individuals than previously thought. It has been proposed that estimating ROH on a genome-wide level, by making use of the genome-wide single nucleotide polymorphism (SNP) data, will enable to indentify recessive variants underlying complex traits. Here, we examined ROH larger than 1.5 Mb individually and in combination for association with survival in 5974 participants of the Rotterdam Study. In addition, we assessed the role of overall homozygosity, expressed as a percentage of the autosomal genome that is in ROH longer than 1.5 Mb, on survival during a mean follow-up period of 12 years. None of these measures of homozygosity was associated with survival to old age.

A population-based study of copy number variants and regions of homozygosity in healthy Swedish individuals

The abundance of copy number variants (CNVs) and regions of homozygosity (ROHs) have been well documented in previous studies. In addition, their roles in complex diseases and traits have since been increasingly appreciated. However, only a limited amount of CNV and ROH data is currently available for the Swedish population. We conducted a population-based study to detect and characterize CNVs and ROHs in 87 randomly selected healthy Swedish individuals using the Affymetrix SNP Array 6.0. More than 600 CNV loci were detected in the population using two different CNV-detection algorithms (PennCNV and Birdsuite). A total of 196 loci were consistently identified by both algorithms, suggesting their reliability. Numerous disease-associated and pharmacogenetics-related genes were found to be overlapping with common CNV loci such as CFHR1/R3, LCE3B/3C, UGT2B17 and GSTT1. Correlation analysis between copy number polymorphisms (CNPs) and genome-wide association studies-identified single-nucleotide polymorphisms also indicates the potential roles of several CNPs as causal variants for diseases and traits such as body mass index, Crohn's disease and multiple sclerosis. In addition, we also identified a total of 14 815 ROHs 500 kb or 2814 ROHs 1M in the Swedish individuals with an average of 170 and 32 regions detected per individual respectively. Approximately 141 Mb or 4.92% of the genome is homozygous in each individual of the Swedish population. This is the first population-based study to investigate the population characteristics of CNVs and ROHs in the Swedish population. This study found many CNV loci that warrant further investigation, and also highlighted the abundance and importance of investigating ROHs for their associations with complex diseases and traits.

Genome-wide analysis of genetic susceptibility to language impairment in an isolated Chilean population

Specific language impairment (SLI) is an unexpected deficit in the acquisition of language skills and affects between 5 and 8% of pre-school children. Despite its prevalence and high heritability, our understanding of the aetiology of this disorder is only emerging. In this paper, we apply genome-wide techniques to investigate an isolated Chilean population who exhibit an increased frequency of SLI. Loss of heterozygosity (LOH) mapping and parametric and non-parametric linkage analyses indicate that complex genetic factors are likely to underlie susceptibility to SLI in this population. Across all analyses performed, the most consistently implicated locus was on chromosome 7q. This locus achieved highly significant linkage under all three non-parametric models (max NPL = 6.73, P = 4.0 × 10(-11)). In addition, it yielded a HLOD of 1.24 in the recessive parametric linkage analyses and contained a segment that was homozygous in two affected individuals. Further, investigation of this region identified a two-SNP haplotype that occurs at an increased frequency in language-impaired individuals (P = 0.008). We hypothesise that the linkage regions identified here, in particular that on chromosome 7, may contain variants that underlie the high prevalence of SLI observed in this isolated population and may be of relevance to other populations affected by language impairments.

Glycomics meets lipidomics--associations of N-glycans with classical lipids, glycerophospholipids, and sphingolipids in three European populations

Recently, high-throughput technologies have been made available which allow the measurement of a broad spectrum of glycomics and lipidomics parameters in many samples. The aim of this study was to apply these methods and investigate associations between 46 glycan and 183 lipid traits measured in blood of 2041 Europeans from three different local populations (Croatia - VIS cohort; Sweden - NSPHS cohort; Great Britain - ORCADES cohort). N-glycans have been analyzed with High Performance Liquid Chromatography (HPLC) and lipids with Electrospray Ionization Tandem Mass Spectrometry (ESI-MS/MS) covering sterol lipids, glycerolipids, glycerophospholipids and sphingolipids in eight subclasses. Overall, 8418 associations were calculated using linear mixed effect models adjusted for pedigree, sex, age and multiple testing. We found 330 significant correlations in VIS. Pearson's correlation coefficient r ranged from -0.27 to 0.34 with corresponding p-values between 1.45 × 10(-19) and 4.83 × 10(-6), indicating statistical significance. A total of 71 correlations in VIS could be replicated in NSPHS (r = [-0.19; 0.35], p = [4.16 × 10(-18); 9.38 × 10(-5)]) and 31 correlations in VIS were also found in ORCADES (r = [-0.20; 0.24], p = [2.69 × 10(-10); 7.55 × 10(-5)]). However, in total only 10 correlations between a subset of triantennary glycans and unsaturated phosphatidylcholine, saturated ceramide, and sphingomyelin lipids in VIS (r = [0.18; 0.34], p = [2.98 × 10(-21); 1.69 × 10(-06)]) could be replicated in both NSPHS and ORCADES. In summary, the results show strong and consistent associations between certain glycans and lipids in all populations, but also population-specific correlations which may be caused by environmental and genetic differences. These associations point towards potential interactive metabolic pathways.

Identification of novel schizophrenia loci by homozygosity mapping using DNA microarray analysis

The recent development of high-resolution DNA microarrays, in which hundreds of thousands of single nucleotide polymorphisms (SNPs) are genotyped, enables the rapid identification of susceptibility genes for complex diseases. Clusters of these SNPs may show runs of homozygosity (ROHs) that can be analyzed for association with disease. An analysis of patients whose parents were first cousins enables the search for autozygous segments in their offspring. Here, using the Affymetrix® Genome-Wide Human SNP Array 5.0 to determine ROHs, we genotyped 9 individuals with schizophrenia (SCZ) whose parents were first cousins. We identified overlapping ROHs on chromosomes 1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16, 17, 19, 20, and 21 in at least 3 individuals. Only the locus on chromosome 5 has been reported previously. The ROHs on chromosome 5q23.3–q31.1 include the candidate genes histidine triad nucleotide binding protein 1 (HINT1) and acyl-CoA synthetase long-chain family member 6 (ACSL6). Other overlapping ROHs may contain novel rare recessive variants that affect SCZ specifically in our samples, given the highly heterozygous nature of SCZ. Analysis of patients whose parents are first cousins may provide new insights for the genetic analysis of psychiatric diseases.

Does inbreeding affect N-glycosylation of human plasma proteins?

Inbreeding depression and heterosis are the two ends of phenotypic changes defined by the genome-wide homozygosity. The aim of this study was to investigate the association of genetic marker-based homozygosity estimates with 46 N-glycan features measured in human plasma. The study was based on a total of 2,341 subjects, originating from three isolated island communities in Croatia (Vis and Korcula islands) and Scotland (Orkney Islands). Inbreeding estimates were associated with an increase in tetrantennary and tetrasialylated glycans, and a decrease in digalactosylated glycans (P < 0.001). The strength of this association was proportional to the mean cohort-based inbreeding coefficient. Increase in tetraantennary glycans is known to be associated with various tumours and their association with inbreeding might be one of the mechanisms underlying the increased prevalence of tumours reported in some human isolated populations. Further studies are thus merited in order to confirm the association of inbreeding with changes in glycan profiles in other plant and animal populations, thus attempting to establish if glycosylation could indeed be involved in mediation of some phenotypic changes described in inbred and outbred organisms.

SNP-array based whole genome homozygosity mapping: a quick and powerful tool to achieve an accurate diagnosis in LGMD2 patients

A large number of novel disease genes have been identified by homozygosity mapping and the positional candidate approach. In this study we used single nucleotide polymorphism (SNP) array-based, whole genome homozygosity mapping as the first step to a molecular diagnosis in the highly heterogeneous muscle disease, limb girdle muscular dystrophy (LGMD). In a consanguineous family, both affected siblings showed homozygous blocks on chromosome 15 corresponding to the LGMD2A locus. Direct sequencing of CAPN3, encoding calpain-3, identified a homozygous deletion c.483delG (p.Ile162SerfsX17). In a sporadic LGMD patient complete absence of caveolin-3 on Western blot was observed. However, a mutation in CAV3 could not be detected. Homozygosity mapping revealed a large homozygous block at the LGMD2I locus, and direct sequencing of FKRP encoding fukutin-related-protein detected the common homozygous c.826 C>A (p.Leu276Ile) mutation. Subsequent re-examination of this patient's muscle biopsy showed aberrant α-dystroglycan glycosylation. In summary, we show that whole-genome homozygosity mapping using low cost SNP arrays provides a fast and non-invasive method to identify disease-causing mutations in sporadic patients or sibs from consanguineous families in LGMD2. Furthermore, this is the first study describing that in addition to PTRF, encoding polymerase I and transcript release factor, FKRP mutations may cause secondary caveolin-3 deficiency.

UPD detection using homozygosity profiling with a SNP genotyping microarray

Single nucleotide polymorphism (SNP) based chromosome microarrays provide both a high-density whole genome analysis of copy number and genotype. In the past 21 months we have analyzed over 13,000 samples primarily referred for developmental delay using the Affymetrix SNP/CN 6.0 version array platform. In addition to copy number, we have focused on the relative distribution of allele homozygosity (HZ) throughout the genome to confirm a strong association of uniparental disomy (UPD) with regions of isoallelism found in most confirmed cases of UPD. We sought to determine whether a long contiguous stretch of HZ (LCSH) greater than a threshold value found only in a single chromosome would correlate with UPD of that chromosome. Nine confirmed UPD cases were retrospectively analyzed with the array in the study, each showing the anticipated LCSH with the smallest 13.5 Mb in length. This length is well above the average longest run of HZ in a set of control patients and was then set as the prospective threshold for reporting possible UPD correlation. Ninety-two cases qualified at that threshold, 46 of those had molecular UPD testing and 29 were positive. Including retrospective cases, 16 showed complete HZ across the chromosome, consistent with total isoUPD. The average size LCSH in the 19 cases that were not completely HZ was 46.3 Mb with a range of 13.5-127.8 Mb. Three patients showed only segmental UPD. Both the size and location of the LCSH are relevant to correlation with UPD. Further studies will continue to delineate an optimal threshold for LCSH/UPD correlation.

hzAnalyzer: detection, quantification, and visualization of contiguous homozygosity in high-density genotyping datasets

The analysis of contiguous homozygosity (runs of homozygous loci) in human genotyping datasets is critical in the search for causal disease variants in monogenic disorders, studies of population history and the identification of targets of natural selection. Here, we report methods for extracting homozygous segments from high-density genotyping datasets, quantifying their local genomic structure, identifying outstanding regions within the genome and visualizing results for comparative analysis between population samples.

A combined long-range phasing and long haplotype imputation method to impute phase for SNP genotypes

BACKGROUND

Knowing the phase of marker genotype data can be useful in genome-wide association studies, because it makes it possible to use analysis frameworks that account for identity by descent or parent of origin of alleles and it can lead to a large increase in data quantities via genotype or sequence imputation. Long-range phasing and haplotype library imputation constitute a fast and accurate method to impute phase for SNP data.

METHODS

A long-range phasing and haplotype library imputation algorithm was developed. It combines information from surrogate parents and long haplotypes to resolve phase in a manner that is not dependent on the family structure of a dataset or on the presence of pedigree information.

RESULTS

The algorithm performed well in both simulated and real livestock and human datasets in terms of both phasing accuracy and computation efficiency. The percentage of alleles that could be phased in both simulated and real datasets of varying size generally exceeded 98% while the percentage of alleles incorrectly phased in simulated data was generally less than 0.5%. The accuracy of phasing was affected by dataset size, with lower accuracy for dataset sizes less than 1000, but was not affected by effective population size, family data structure, presence or absence of pedigree information, and SNP density. The method was computationally fast. In comparison to a commonly used statistical method (fastPHASE), the current method made about 8% less phasing mistakes and ran about 26 times faster for a small dataset. For larger datasets, the differences in computational time are expected to be even greater. A computer program implementing these methods has been made available.

CONCLUSIONS

The algorithm and software developed in this study make feasible the routine phasing of high-density SNP chips in large datasets.

Consanguinity around the world: what do the genomic data of the HGDP-CEPH diversity panel tell us?

Inbreeding coefficients and consanguineous mating types are usually inferred from population surveys or pedigree studies. Here, we present a method to estimate them from dense genome-wide single-nucleotide polymorphism genotypes and apply it to 940 unrelated individuals from the Human Genome Diversity Panel (HGDP-CEPH). Inbreeding is observed in almost all populations of the panel, and the highest inbreeding levels and frequencies of inbred individuals are found in populations of the Middle East, Central South Asia and the Americas. In these regions, first cousin (1C) marriages are the most frequent, but we also observed marriages between double first cousins (2 × 1C) and between avuncular (AV) pairs. Interestingly, if 2 × 1C marriages are preferred to AV marriages in Central South Asia and the Middle East, the contrary is found in the Americas. There are thus some regional trends but there are also some important differences between populations within a region. Individual results can be found on the CEPH website at ftp://ftp.cephb.fr/hgdp_hbd/.

Identification of genes with allelic imbalance on 6p associated with nasopharyngeal carcinoma in southern Chinese

Nasopharyngeal carcinoma (NPC) is a malignancy of epithelial origin. The etiology of NPC is complex and includes multiple genetic and environmental factors. We employed case-control analysis to study the association of chromosome 6p regions with NPC. In total, 360 subjects and 360 healthy controls were included, and 233 single nucleotide polymorphisms (SNPs) on 6p were examined. Significant single-marker associations were found for SNPs rs2267633 (p = 4.49 × 10(-5)), rs2076483 (most significant, p = 3.36 × 10(-5)), and rs29230 (p=1.43 × 10(-4)). The highly associated genes were the gamma-amino butyric acid B receptor 1 (GABBR1), human leukocyte antigen (HLA-A), and HLA complex group 9 (HCG9). Haplotypic associations were found for haplotypes AAA (located within GABBR1, p-value  = 6.46 × 10(-5)) and TT (located within HLA-A, p = 0.0014). Further investigation of the homozygous genotype frequencies between cases and controls suggested that micro-deletion regions occur in GABBR1 and neural precursor cell expressed developmentally down-regulated 9 (NEDD9). Quantitative real-time polymerase chain reaction (qPCR) using 11 pairs of NPC biopsy samples confirmed the significant decline in GABBR1 and NEDD9 mRNA expression in the cancer tissues compared to the adjacent non-tumor tissue (p<0.05). Our study demonstrates that multiple chromosome 6p susceptibility loci contribute to the risk of NPC, possibly though GABBR1 and NEDD9 loss of function.

Genomic and genealogical investigation of the French Canadian founder population structure

Characterizing the genetic structure of worldwide populations is important for understanding human history and is essential to the design and analysis of genetic epidemiological studies. In this study, we examined genetic structure and distant relatedness and their effect on the extent of linkage disequilibrium (LD) and homozygosity in the founder population of Quebec (Canada). In the French Canadian founder population, such analysis can be performed using both genomic and genealogical data. We investigated genetic differences, extent of LD, and homozygosity in 140 individuals from seven sub-populations of Quebec characterized by different demographic histories reflecting complex founder events. Genetic findings from genome-wide single nucleotide polymorphism data were correlated with genealogical information on each of these sub-populations. Our genomic data showed significant population structure and relatedness present in the contemporary Quebec population, also reflected in LD and homozygosity levels. Our extended genealogical data corroborated these findings and indicated that this structure is consistent with the settlement patterns involving several founder events. This provides an independent and complementary validation of genomic-based studies of population structure. Combined genomic and genealogical data in the Quebec founder population provide insights into the effects of the interplay of two important sources of bias in genetic epidemiological studies, unrecognized genetic structure and cryptic relatedness.

Regions of homozygosity and their impact on complex diseases and traits

Regions of homozygosity (ROHs) are more abundant in the human genome than previously thought. These regions are without heterozygosity, i.e. all the genetic variations within the regions have two identical alleles. At present there are no standardized criteria for defining the ROHs resulting in the different studies using their own criteria in the analysis of homozygosity. Compared to the era of genotyping microsatellite markers, the advent of high-density single nucleotide polymorphism genotyping arrays has provided an unparalleled opportunity to comprehensively detect these regions in the whole genome in different populations. Several studies have identified ROHs which were associated with complex phenotypes such as schizophrenia, late-onset of Alzheimer's disease and height. Collectively, these studies have conclusively shown the abundance of ROHs larger than 1 Mb in outbred populations. The homozygosity association approach holds great promise in identifying genetic susceptibility loci harboring recessive variants for complex diseases and traits.

Genomic runs of homozygosity record population history and consanguinity

The human genome is characterised by many runs of homozygous genotypes, where identical haplotypes were inherited from each parent. The length of each run is determined partly by the number of generations since the common ancestor: offspring of cousin marriages have long runs of homozygosity (ROH), while the numerous shorter tracts relate to shared ancestry tens and hundreds of generations ago. Human populations have experienced a wide range of demographic histories and hold diverse cultural attitudes to consanguinity. In a global population dataset, genome-wide analysis of long and shorter ROH allows categorisation of the mainly indigenous populations sampled here into four major groups in which the majority of the population are inferred to have: (a) recent parental relatedness (south and west Asians); (b) shared parental ancestry arising hundreds to thousands of years ago through long term isolation and restricted effective population size (N_e), but little recent inbreeding (Oceanians); (c) both ancient and recent parental relatedness (Native Americans); and (d) only the background level of shared ancestry relating to continental N_e (predominantly urban Europeans and East Asians; lowest of all in sub-Saharan African agriculturalists), and the occasional cryptically inbred individual. Moreover, individuals can be positioned along axes representing this demographic historic space. Long runs of homozygosity are therefore a globally widespread and under-appreciated characteristic of our genomes, which record past consanguinity and population isolation and provide a distinctive record of the demographic history of an individual's ancestors. Individual ROH measures will also allow quantification of the disease risk arising from polygenic recessive effects.

Population structure and genome-wide patterns of variation in Ireland and Britain

Located off the northwestern coast of the European mainland, Britain and Ireland were among the last regions of Europe to be colonized by modern humans after the last glacial maximum. Further, the geographical location of Britain, and in particular of Ireland, is such that the impact of historical migration has been minimal. Genetic diversity studies applying the Y chromosome and mitochondrial systems have indicated reduced diversity and an increased population structure across Britain and Ireland relative to the European mainland. Such characteristics would have implications for genetic mapping studies of complex disease. We set out to further our understanding of the genetic architecture of the region from the perspective of (i) population structure, (ii) linkage disequilibrium (LD), (iii) homozygosity and (iv) haplotype diversity (HD). Analysis was conducted on 3654 individuals from Ireland, Britain (with regional sampling in Scotland), Bulgaria, Portugal, Sweden and the Utah HapMap collection. Our results indicate a subtle but clear genetic structure across Britain and Ireland, although levels of structure were reduced in comparison with average cross-European structure. We observed slightly elevated levels of LD and homozygosity in the Irish population compared with neighbouring European populations. We also report on a cline of HD across Europe with greatest levels in southern populations and lowest levels in Ireland and Scotland. These results are consistent with our understanding of the population history of Europe and promote Ireland and Scotland as relatively homogenous resources for genetic mapping of rare variants.