Leveling the Playing Field in Homozygosity Mapping Using Map Distances

Studies of linkage disequilibrium (LD) and its variation in the genome are of central importance for understanding evolutionary history, population structure, and selective sweeps. Extreme forms of the latter may result in runs of homozygosity (ROH). In human gene mapping, long ROHs are the basis for homozygosity mapping (HM) with length measured in terms of Mb (10⁶ base pairs physical distance). LD varies greatly over the human genome so that long ROHs tend to occur preferentially in regions of high LD and ROHs of the same length in different regions are not strictly comparable. Thus, in human gene mapping, LD appears as a confounder that needs to be taken into account in the interpretation of ROHs. The effect of varying LD can be mitigated by working on a scale of centimorgans (cM, genetic distance) instead of Mb. We demonstrate this effect for HapMap 3 data on chromosome 19 and show examples with different ROH lengths depending on whether physical or genetic lengths are used. These results suggest that HM should preferably be done on genetic rather than physical distances.     

Genome-wide analysis of runs of homozygosity identifies new susceptibility regions of lung cancer in Han Chinese

Runs of homozygosity (ROHs) are a class of important but poorly studied genomic variations and may be involved in individual susceptibility to diseases. To better understand ROH and its relationship with lung cancer, we performed a genome-wide ROH analysis of a subset of a previous genome-wide case-control study (1,473 cases and 1,962 controls) in a Han Chinese population. ROHs were classified into two classes, based on lengths, intermediate and long ROHs, to evaluate their association with lung cancer risk using existing genome-wide single nucleotide polymorphism (SNP) data. We found that the overall level of intermediate ROHs was significantly associated with a decreased risk of lung cancer (odds ratio = 0.63; 95% confidence interval: 0.51-0.77; P = 4.78×10⁻⁶ ), while the long ROHs seemed to be a risk factor of lung cancer. We also identified one ROH region at 14q23.1 that was consistently associated with lung cancer risk in the study. These results indicated that ROHs may be a new class of variation which may be associated with lung cancer risk, and genetic variants at 14q23.1 may be involved in the development of lung cancer.     

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.     

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.     

Extended runs of homozygosity at 17q11.2: an association with type‐2 NF1 deletions?

Large deletions in the NF1 gene region at 17q11.2 are caused by nonallelic homologous recombination (NAHR). The recurrent type‐2 NF1 deletions span 1.2 Mb, with breakpoints in the SUZ12 gene and SUZ12P. Type‐2 NF1 deletions occur preferentially during mitosis and are associated with somatic mosaicism. A panel of 16 type‐2 NF1 deletions was used as a model system in which to investigate whether extended homozygosity across 17q11.2 might be associated with somatic deletion. Using SNP arrays, a 3.2 Mb interval encompassing the NF1 deletion region was found to harbor runs of homozygosity (ROHs) in different human populations. However, ROHs ≥500 kb directly flanking the NF1 deletion region on both sides were not found to occur disproportionately in NF1 patients harboring type‐2 deletions compared to controls. Although low allelic diversity in 17q11.2 is unlikely to be a key factor in promoting NAHR‐mediated somatic type‐2 deletions, a specific ROH of 588 kb (roh1), located some 525 kb proximal to the deletion interval, was found to occur more frequently (P=0.012) in the type‐2 deletion patients compared with controls. We postulate that roh1 may act remotely, via an as yet unknown mechanism, to increase the frequency of somatic recombination between the distally duplicated SUZ12 sequences. Hum Mutat 30:1–10, 2010. © 2010 Wiley‐Liss, Inc.     

Estimating Human Inbreeding Coefficients: Comparison of Genealogical and Marker Heterozygosity Approaches

We have used genealogies and genomic polymorphisms to estimate individual inbreeding coefficients (F) in 50 subjects with an expected range (based on recent genealogies) of F from 0.0 to 0.0625. The estimates were based on two approaches, using genotypes respectively from 410 microsatellite markers (410‐STR panel) and from 10,000 SNPs (10K‐SNP panel). The latter was performed in a sub‐sample of 15 individuals. We concluded that for both marker panels measures of inbreeding based on the excess of homozygosity over Hardy‐Weinberg expectation were not closely correlated with 4‐5 generation genealogical F‐values. For the 10K‐SNP panel we found two alternative measures which correlated more closely with F, based respectively on standard errors and on paired homozygosity of nearby SNPs over distances of 2‐4 cM. We propose an empirical method for estimating standard errors and hence individual F‐values, based on the variation between individual autosomes. This method could provide useful estimates of average F‐values for groups of individuals in population‐based studies of the effects of inbreeding/homozygosity on quantitative traits.     

Regions of homozygosity identified by oligonucleotide SNP arrays: evaluating the incidence and clinical utility

Copy neutral segments with allelic homozygosity, also known as regions of homozygosity (ROHs), are frequently identified in cases interrogated by oligonucleotide single-nucleotide polymorphism (oligo-SNP) microarrays. Presence of ROHs may be because of parental relatedness, chromosomal recombination or rearrangements and provides important clues regarding ancestral homozygosity, consanguinity or uniparental disomy. In this study of 14 574 consecutive cases, 832 (6%) were found to harbor one or more ROHs over 10 Mb, of which 651 cases (78%) had multiple ROHs, likely because of identity by descent (IBD), and 181 cases (22%) with ROHs involving a single chromosome. Parental relatedness was predicted to be first degree or closer in 5%, second in 9% and third in 19%. Of the 181 cases, 19 had ROHs for a whole chromosome revealing uniparental isodisomy (isoUPD). In all, 25 cases had significant ROHs involving a single chromosome; 5 cases were molecularly confirmed to have a mixed iso- and heteroUPD15 and 1 case each with segmental UPD9pat and segmental UPD22mat; 17 cases were suspected to have a mixed iso- and heteroUPD including 2 cases with small supernumerary marker and 2 cases with mosaic trisomy. For chromosome 15, 12 (92%) of 13 molecularly studied cases had either Prader–Willi or Angelman syndrome. Autosomal recessive disorders were confirmed in seven of nine cases from eight families because of the finding of suspected gene within a ROH. This study demonstrates that ROHs are much more frequent than previously recognized and often reflect parental relatedness, ascertain autosomal recessive diseases or unravel UPD in many cases.     

Extended tracts of homozygosity in outbred human populations

Long tracts of consecutive homozygous single nucleotide polymorphisms (SNPs) can arise in the genome through a number of mechanisms. These include inbreeding in which an individual inherits chromosomal segments that are identical by descent from each parent. However, recombination and other processes break up chromosomal segments over generations. The longest tracts are therefore to be expected in populations with an appreciable degree of inbreeding. We examined the length, number and distribution of long tracts of homozygosity in the apparently outbred HapMap populations. We observed 1393 tracts exceeding 1 Mb in length among the 209 unrelated HapMap individuals. The longest was an uninterrupted run of 3922 homozygous SNPs spanning 17.9 Mb in a Japanese individual. We find that homozygous tracts are significantly more common in regions with high linkage disequilibrium and low recombination, and the location of tracts is similar across all populations. The Yoruba sample has the fewest long tracts per individual, consistent with a larger number of generations (and hence amount of recombination) since the founding of that population. Our results suggest that multiple-megabase-scale ancestral haplotypes persist in outbred human populations in broad genomic regions which have lower than average recombination rates. We observed three outlying individuals who have exceptionally long and numerous homozygous tracts that are not associated with recombination suppressed areas of the genome. We consider that this reflects a high level of relatedness in their ancestry which is too recent to have been influenced by the local recombination intensity. Possible alternative mechanisms and the implications of long homozygous tracts in the genome are discussed.     

Genome-wide estimates of inbreeding in unrelated individuals and their association with cognitive ability

The consequence of reduced cognitive ability from inbreeding has long been investigated, mainly restricted to cousin–cousin marriages. Molecular genetic techniques now allow us to test the relationship between increased ancestral inbreeding and cognitive ability in a population of traditionally unrelated individuals. In a representative UK sample of 2329 individuals, we used genome-wide SNP data to estimate the percentage of the genome covered by runs of homozygous SNPs (ROH). This was tested for association with general cognitive ability, as well as measures of verbal and non-verbal ability. Further, association was tested between these traits and specific ROH. Burden of ROH was not associated with cognitive ability after correction for multiple testing, although burden of ROH was nominally associated with increased non-verbal cognitive ability (P=0.03). Moreover, although no individual ROH was significantly associated with cognitive ability, there was a significant bias towards increased cognitive ability in carriers of ROH (P=0.002). A potential explanation for these results is increased positive assortative mating in spouses with higher cognitive ability, although we found no evidence in support of this hypothesis in a separate sample. Reduced minor allele frequency across the genome was associated with higher cognitive ability, which could contribute to an apparent increase in ROH. This may reflect minor alleles being more likely to be deleterious.     

Risk of breast and prostate cancer is not associated with increased homozygosity in outbred populations

Regions of restricted genetic heterogeneity due to identity by descent (autozygosity) are known to confer susceptibility to a number of diseases. Regions of germline homozygosity (ROHs) of 1–2 Mb, the result of autozygosity, are detectable at high frequency in outbred populations. Recent studies have reported that ROHs, possibly through exposing recessive disease-causing alleles or alternative mechanisms, are associated with an increased cancer risk. To examine whether homozygosity is associated with breast or prostate cancer risk, we analysed 500K single-nucleotide polymorphism data from two genome-wide association studies conducted by the Cancer Genetics Markers of Susceptibility initiatives (http://cgems.cancer.gov/). Six common ROHs were associated with breast cancer risk and four with prostate cancer (P<0.01). Intriguingly, one of the breast cancer ROHs maps to 6q22.31–6q22.3, a region that has been previously shown to confer breast cancer risk. Although none of the ROHs remained significantly associated with cancer risk after adjustment for multiple testing, a number of ROHs merit further interrogation. However, our findings provide no strong evidence that levels of measured homozygosity, whatever their aetiology (autozygosity, uniparental isodisomy or hemizygosity), confer an increased risk of developing breast or prostate cancer in predominantly outbred populations.     

Association Between Autozygosity and Major Depression: Stratification Due to Religious Assortment

The effects of inbreeding on the health of offspring can be studied by measuring genome-wide autozygosity as the proportion of the genome in runs of homozygosity (F _roh) and relate F _roh to outcomes such as psychiatric phenotypes. To successfully conduct these studies, the main patterns of variation for genome-wide autozygosity between and within populations should be well understood and accounted for. Within population variation was investigated in the Dutch population by comparing autozygosity between religious and non-religious groups. The Netherlands have a history of societal segregation and assortment based on religious affiliation, which may have increased parental relatedness within religious groups. Religion has been associated with several psychiatric phenotypes, such as major depressive disorder (MDD). We investigated whether there is an association between autozygosity and MDD, and the extent to which this association can be explained by religious affiliation. All F _roh analyses included adjustment for ancestry-informative principal components (PCs) and geographic factors. Religious affiliation was significantly associated with autozygosity, showing that F _roh has the ability to capture within population differences that are not captured by ancestry-informative PCs or geographic factors. The non-religious group had significantly lower F _roh values and significantly more MDD cases, leading to a nominally significant negative association between autozygosity and depression. After accounting for religious affiliation, MDD was not associated with F _roh, indicating that the relation between MDD and inbreeding was due to stratification. This study shows how past religious assortment and recent secularization can have genetic consequences in a relatively small country. This warrants accounting for the historical social context and its effects on genetic variation in association studies on psychiatric and other related traits.     

Inbreeding among Caribbean Hispanics from the Dominican Republic and its effects on risk of Alzheimer disease

BackgroundInbreeding can be associated with a modification of disease risk due to excess homozygosity of recessive alleles affecting a wide range of phenotypes. We estimated the inbreeding coefficient in Caribbean Hispanics and examined its effects on risk of late-onset Alzheimer disease.MethodsThe inbreeding coefficient was calculated in 3,392 subjects (1,451 late-onset Alzheimer disease patients and 1,941 age-matched healthy controls) of Caribbean Hispanic ancestry using 177,997 nearly independent single-nucleotide polymorphisms from genome-wide array. The inbreeding coefficient was estimated using the excess homozygosity method with and without adjusting for admixture.ResultsThe average inbreeding coefficient in Caribbean Hispanics without accounting for admixture was F = 0.018 (±0.048), suggesting a mating equivalent to that of second cousins or second cousins once removed. Adjusting for admixture from three parent populations, the average inbreeding coefficient was found to be 0.0034 (±0.019) or close to third-cousin mating. Inbreeding coefficient was a significant predictor of Alzheimer disease when age, sex, and APOE genotype were used as adjusting covariates (P = 0.03).ConclusionThe average inbreeding coefficient of this population is significantly higher than that of the general Caucasian populations in North America. The high rate of inbreeding resulting in increased frequency of recessive variants is advantageous for the identification of rare variants associated with late-onset Alzheimer disease.Genet Med17 8, 639–643.     

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. © 2011 Wiley‐Liss, Inc.     

HDR‐del: A tool based on Hamming distance for prioritizing pathogenic chromosomal deletions in exome sequencing

High‐density oligonucleotide arrays have widely been used to detect pathogenic chromosomal deletions. In addition to high‐density oligonucleotide arrays, programs using whole‐exome sequencing have become available for estimating copy‐number variations using depth of coverage. Here, we propose a new statistical method, HDR‐del, to prioritize pathogenic chromosomal deletions based on Hamming distance in exome sequencing. In vcf (variant call format) files generated from exome sequencing, hemizygous chromosomal deletion regions lack heterozygous variants and lead to apparent long runs of homozygosity (ROH). In our Hamming distance ratio (HDR)‐del approach, we calculate the “difference” in heterozygous status between an affected individual and control individuals using the HDR over all candidate chromosomal deletion regions defined as ROH longer than 1Mbp. Using a suitable test statistic, which is expected to be large for a true pathogenic deletion region, we prioritize candidate chromosomal deletion regions based on this statistic. In our approach, we were able to considerably narrow down true pathogenic chromosomal deletion regions, which were confirmed by high‐density oligonucleotide arrays in four mitochondrial disease patients. Our HDR‐del approach represents an easy method for detecting chromosomal deletions.     

Genetic meta‐analysis of obsessive–compulsive disorder and self‐report compulsive symptoms

We investigated whether obsessive–compulsive (OC) symptoms from a population‐based sample could be analyzed to detect genetic variants influencing obsessive–compulsive disorder (OCD). We performed a genome‐wide association studies (GWAS) on the obsession (rumination and impulsions) and compulsion (checking, washing, and ordering/precision) subscales of an abbreviated version of the Padua Inventory (N = 8,267 with genome‐wide genotyping and phenotyping). The compulsion subscale showed a substantial and significant positive genetic correlation with an OCD case–control GWAS (r _G = 0.61, p = .017) previously published by the Psychiatric Genomics Consortium (PGC‐OCD). The obsession subscale and the total Padua score showed no significant genetic correlations (r _G = ?0.02 and r _G = 0.42, respectively). A meta‐analysis of the compulsive symptoms GWAS with the PGC‐OCD revealed no genome‐wide significant Single‐Nucleotide Polymorphisms (SNPs combined N = 17,992, indicating that the power is still low for individual SNP effects). A gene‐based association analysis, however, yielded two novel genes (WDR7 and ADCK1). The top 250 genes in the gene‐based test also showed a significant increase in enrichment for psychiatric and brain‐expressed genes. S‐Predixcan testing showed that for genes expressed in hippocampus, amygdala, and caudate nucleus significance increased in the meta‐analysis with compulsive symptoms compared to the original PGC‐OCD GWAS. Thus, the inclusion of dimensional symptom data in genome‐wide association on clinical case–control GWAS of OCD may be useful to find genes for OCD if the data are based on quantitative indices of compulsive behavior. SNP‐level power increases were limited, but aggregate, gene‐level analyses showed increased enrichment for brain‐expressed genes related to psychiatric disorders, and increased association with gene expression in brain tissues with known emotional, reward processing, memory, and fear‐formation functions.     

Colorectal cancer risk is not associated with increased levels of homozygosity in Saudi Arabia

PurposeRuns of homozygosity (ROHs) represent a measure of the extent of autozygosity and are correlated with the extent of inbreeding. Recently, it has been suggested that ROHs may contribute to the risk of colorectal cancer (CRC). The high rate of consanguinity and CRC in the Saudi population prompted us to test the role of autozygosity in the CRC risk.MethodsWe compared 48 Saudi CRC patients to 100 ethnically matched controls, processed on the Affymetrix 250K StyI SNP GeneChip platform and analyzed using the plink package.ResultsWe could find no evidence of a significant relationship between autozygosity and CRC risk.ConclusionThe negative results in our study add additional significance to what has been previously reported in literature, as this is the first study to address these questions in an inbred population. Our subgroup analysis of patients with microsatellite unstable–positive tumors as compared with other groups did not significantly change our results. Although these results do not rule out the presence of recessively acting CRC-predisposing genes in a small percentage of patients, which our relatively small sample size could not capture, they suggest that such genes are unlikely to account for the disturbingly high incidence of CRC in our consanguineous population.Genet Med 2012:14(8):720–728     

Interpretation and reporting of large regions of homozygosity and suspected consanguinity/uniparental disomy, 2021 revision: A technical standard of the American College of Medical Genetics and Genomics (ACMG)

Genomic testing, including single-nucleotide variation (formerly single-nucleotide polymorphism)–based chromosomal microarray and exome and genome sequencing, can detect long regions of homozygosity (ROH) within the genome. Genomic testing can also detect possible uniparental disomy (UPD). Platforms that can detect ROH and possible UPD have matured since the initial American College of Medical Genetics and Genomics (ACMG) standard was published in 2013, and the detection of ROH and UPD by these platforms has shown utility in diagnosis of patients with genetic/genomic disorders. The presence of these segments, when distributed across multiple chromosomes, may indicate a familial relationship between the proband’s parents. This technical standard describes the detection of possible consanguinity and UPD by genomic testing, as well as the factors confounding the inference of a specific parental relationship or UPD. Current bioethical and legal issues regarding detection and reporting of consanguinity are also discussed.     

TRAPPC9-related neurodevelopmental disorder: Report of a homozygous deletion in TRAPPC9 due to paternal uniparental isodisomy

TRAPPC9 loss-of-function biallelic variants are associated with an autosomal recessive intellectual disability syndrome (Online Mendelian Inheritance of Man no. 613192), also characterized by microcephaly, hypertelorism, obesity, growth delay, and behavioral differences. Here, we describe an 8-year-old Hispanic female with neurodevelopmental disorder, partial epilepsy, microcephaly, bilateral cleft lip and alveolus, growth delay, and dysmorphic features. She had abnormal myelination, mega cisterna magna, and colpocephaly on brain magnetic resonance imaging (MRI). Microarray showed a single ~146 Mb region of homozygosity (ROH) encompassing all of Chromosome 8, consistent with uniparental isodisomy (UPD). Exome sequencing performed in-house did not identify single nucleotide variants to explain her phenotype. Algorithms developed in-house and further evaluation of BAM files revealed a homozygous deletion overlapping Exon 2 in TRAPPC9 within the ROH. Subsequent del/dup analyses with exon-level oligo array confirmed a likely pathogenic deletion in TRAPPC9 (NM_031466.5): arr[GRCh37] 8q24.3(141460661_141461780)x0. Our case highlights the implications of downstream analyses from UPD/ROH given the increased risk for AR conditions, the strengths of combining orthologous molecular methods to establish a diagnosis and further delineates the TRAPPC9-related phenotype in an individual of Hispanic ancestry.