Long homozygous chromosomal segments in reference families from the centre d'Etude du polymorphisme humain

Runs of Homozygosity: Association with Coronary Artery Disease and Gene Expression in Monocytes and Macrophages

Runs of homozygosity (ROHs) are recognized signature of recessive inheritance. Contributions of ROHs to the genetic architecture of coronary artery disease and regulation of gene expression in cells relevant to atherosclerosis are not known. Our combined analysis of 24,320 individuals from 11 populations of white European ethnicity showed an association between coronary artery disease and both the count and the size of ROHs. Individuals with coronary artery disease had approximately 0.63 (95% CI: 0.4-0.8) excess of ROHs when compared to coronary-artery-disease-free control subjects (p = 1.49 × 10(-9)). The average total length of ROHs was approximately 1,046.92 (95% CI: 634.4-1,459.5) kb greater in individuals with coronary artery disease than control subjects (p = 6.61 × 10(-7)). None of the identified individual ROHs was associated with coronary artery disease after correction for multiple testing. However, in aggregate burden analysis, ROHs favoring increased risk of coronary artery disease were much more common than those showing the opposite direction of association with coronary artery disease (p = 2.69 × 10(-33)). Individual ROHs showed significant associations with monocyte and macrophage expression of genes in their close proximity-subjects with several individual ROHs showed significant differences in the expression of 44 mRNAs in monocytes and 17 mRNAs in macrophages when compared to subjects without those ROHs. This study provides evidence for an excess of homozygosity in coronary artery disease in outbred populations and suggest the potential biological relevance of ROHs in cells of importance to the pathogenesis of atherosclerosis.

Old gene, new phenotype: mutations in heparan sulfate synthesis enzyme, EXT2 leads to seizure and developmental disorder, no exostoses

BACKGROUND

Heparan sulfate proteoglycans are vital components of the extracellular matrix and are essential for cellular homeostasis. Many genes are involved in modulating heparan sulfate synthesis, and when these genes are mutated, they can give rise to early-onset developmental disorders affecting multiple body systems. Herein, we describe a consanguineous family of four sibs with a novel disorder, which we designate as seizures-scoliosis-macrocephaly syndrome, characterised by seizures, intellectual disability, hypotonia, scoliosis, macrocephaly, hypertelorism and renal dysfunction.

METHODS

Our application of autozygosity mapping and whole-exome sequencing allowed us to identify mutations in the patients. To confirm the autosomal-recessive mode of inheritance, all available family members were genotyped. We also studied the effect of these mutations on protein expression and function in patient cells and using an in vitro system.

RESULTS

We identified two homozygous mutations p.Met87Arg and p.Arg95 Cys in exostosin 2, EXT2, a ubiquitously expressed gene that encodes a glycosyltransferase required for heparan sulfate synthesis. In patient cells, we observed diminished EXT2 expression and function. We also performed an in vitro assay to determine which mutation has a larger effect on protein expression and observed reduced EXT2 expression in constructs expressing either one of the mutations but a greater reduction when both residues were mutated.

CONCLUSIONS

In short, we have unravelled the genetic basis of a new recessive disorder, seizures-scoliosis-macrocephaly syndrome. Our results have implicated a well-characterised gene in a new developmental disorder and have further illustrated the spectrum of phenotypes that can arise due to errors in glycosylation.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Estimation of inbreeding using pedigree, 50k SNP chip genotypes and full sequence data in three cattle breeds

Background

Levels of inbreeding in cattle populations have increased in the past due to the use of a limited number of bulls for artificial insemination. High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. Various estimators based on different sources, e.g., pedigree or genomic data, have been used to estimate inbreeding coefficients in cattle populations. However, the comparative advantage of using full sequence data to assess inbreeding is unknown. We used pedigree and genomic data at different densities from 50k to full sequence variants to compare how different methods performed for the estimation of inbreeding levels in three different cattle breeds.

Results

Five different estimates for inbreeding were calculated and compared in this study: pedigree based inbreeding coefficient (F_PED); run of homozygosity (ROH)-based inbreeding coefficients (F_ROH); genomic relationship matrix (GRM)-based inbreeding coefficients (F_GRM); inbreeding coefficients based on excess of homozygosity (F_HOM) and correlation of uniting gametes (F_UNI). Estimates using ROH provided the direct estimated levels of autozygosity in the current populations and are free effects of allele frequencies and incomplete pedigrees which may increase in inaccuracy in estimation of inbreeding. The highest correlations were observed between F_ROH estimated from the full sequence variants and the F_ROH estimated from 50k SNP (single nucleotide polymorphism) genotypes. The estimator based on the correlation between uniting gametes (F_UNI) using full genome sequences was also strongly correlated with F_ROH detected from sequence data.

Conclusions

Estimates based on ROH directly reflected levels of homozygosity and were not influenced by allele frequencies, unlike the three other estimates evaluated (F_GRM, F_HOM and F_UNI), which depended on estimated allele frequencies. F_PED suffered from limited pedigree depth. Marker density affects ROH estimation. Detecting ROH based on 50k chip data was observed to give estimates similar to ROH from sequence data. In the absence of full sequence data ROH based on 50k can be used to access homozygosity levels in individuals. However, genotypes denser than 50k are required to accurately detect short ROH that are most likely identical by descent (IBD).

Recessive nephrocerebellar syndrome on the Galloway-Mowat syndrome spectrum is caused by homozygous protein-truncating mutations of WDR73

Galloway-Mowat syndrome (GMS) is a neurodevelopmental disorder characterized by microcephaly, cerebellar hypoplasia, nephrosis, and profound intellectual disability. Jinks et al. extend the GMS spectrum by identifying a novel nephrocerebellar syndrome with selective striatal cholinergic interneuron loss and complete lateral geniculate nucleus delamination, caused by a frameshift mutation in WDR73.

We describe a novel nephrocerebellar syndrome on the Galloway-Mowat syndrome spectrum among 30 children (ages 1.0 to 28 years) from diverse Amish demes. Children with nephrocerebellar syndrome had progressive microcephaly, visual impairment, stagnant psychomotor development, abnormal extrapyramidal movements and nephrosis. Fourteen died between ages 2.7 and 28 years, typically from renal failure. Post-mortem studies revealed (i) micrencephaly without polymicrogyria or heterotopia; (ii) atrophic cerebellar hemispheres with stunted folia, profound granule cell depletion, Bergmann gliosis, and signs of Purkinje cell deafferentation; (iii) selective striatal cholinergic interneuron loss; and (iv) optic atrophy with delamination of the lateral geniculate nuclei. Renal tissue showed focal and segmental glomerulosclerosis and extensive effacement and microvillus transformation of podocyte foot processes. Nephrocerebellar syndrome mapped to 700 kb on chromosome 15, which contained a single novel homozygous frameshift variant (WDR73 c.888delT; p.Phe296Leufs*26). WDR73 protein is expressed in human cerebral cortex, hippocampus, and cultured embryonic kidney cells. It is concentrated at mitotic microtubules and interacts with α-, β-, and γ-tubulin, heat shock proteins 70 and 90 (HSP-70; HSP-90), and the carbamoyl phosphate synthetase 2/aspartate transcarbamylase/dihydroorotase multi-enzyme complex. Recombinant WDR73 p.Phe296Leufs*26 and p.Arg256Profs*18 proteins are truncated, unstable, and show increased interaction with α- and β-tubulin and HSP-70/HSP-90. Fibroblasts from patients homozygous for WDR73 p.Phe296Leufs*26 proliferate poorly in primary culture and senesce early. Our data suggest that in humans, WDR73 interacts with mitotic microtubules to regulate cell cycle progression, proliferation and survival in brain and kidney. We extend the Galloway-Mowat syndrome spectrum with the first description of diencephalic and striatal neuropathology.

Identifying Human Genome-Wide CNV, LOH and UPD by Targeted Sequencing of Selected Regions

Copy-number variations (CNV), loss of heterozygosity (LOH), and uniparental disomy (UPD) are large genomic aberrations leading to many common inherited diseases, cancers, and other complex diseases. An integrated tool to identify these aberrations is essential in understanding diseases and in designing clinical interventions. Previous discovery methods based on whole-genome sequencing (WGS) require very high depth of coverage on the whole genome scale, and are cost-wise inefficient. Another approach, whole exome genome sequencing (WEGS), is limited to discovering variations within exons. Thus, we are lacking efficient methods to detect genomic aberrations on the whole genome scale using next-generation sequencing technology. Here we present a method to identify genome-wide CNV, LOH and UPD for the human genome via selectively sequencing a small portion of genome termed Selected Target Regions (SeTRs). In our experiments, the SeTRs are covered by 99.73%~99.95% with sufficient depth. Our developed bioinformatics pipeline calls genome-wide CNVs with high confidence, revealing 8 credible events of LOH and 3 UPD events larger than 5M from 15 individual samples. We demonstrate that genome-wide CNV, LOH and UPD can be detected using a cost-effective SeTRs sequencing approach, and that LOH and UPD can be identified using just a sample grouping technique, without using a matched sample or familial information.

Characterizing homozygosity across United States, New Zealand and Australian Jersey cow and bull populations

Background

Dairy cattle breeding objectives are in general similar across countries, but environment and management conditions may vary, giving rise to slightly different selection pressures applied to a given trait. This potentially leads to different selection pressures to loci across the genome that, if large enough, may give rise to differential regions with high levels of homozygosity. The objective of this study was to characterize differences and similarities in the location and frequency of homozygosity related measures of Jersey dairy cows and bulls from the United States (US), Australia (AU) and New Zealand (NZ).

Results

The populations consisted of a subset of genotyped Jersey cows born in US (n = 1047) and AU (n = 886) and Jersey bulls progeny tested from the US (n = 736), AU (n = 306) and NZ (n = 768). Differences and similarities across populations were characterized using a principal component analysis (PCA) and a run of homozygosity (ROH) statistic (ROH45), which counts the frequency of a single nucleotide polymorphism (SNP) being in a ROH of at least 45 SNP. Regions that exhibited high frequencies of ROH45 and those that had significantly different ROH45 frequencies between populations were investigated for their association with milk yield traits. Within sex, the PCA revealed slight differentiation between the populations, with the greatest occurring between the US and NZ bulls. Regions with high levels of ROH45 for all populations were detected on BTA3 and BTA7 while several other regions differed in ROH45 frequency across populations, the largest number occurring for the US and NZ bull contrast. In addition, multiple regions with different ROH45 frequencies across populations were found to be associated with milk yield traits.

Conclusion

Multiple regions exhibited differential ROH45 across AU, NZ and US cow and bull populations, an interpretation is that locations of the genome are undergoing differential directional selection. Two regions on BTA3 and BTA7 had high ROH45 frequencies across all populations and will be investigated further to determine the gene(s) undergoing directional selection.

Electronic supplementary material

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

Extensive genome-wide autozygosity in the population isolates of Daghestan

Isolated populations are valuable resources for mapping disease genes, as inbreeding increases genome-wide homozygosity and enhances the ability to map disease alleles on a genetically uniform background within a relatively homogenous environment. The populations of Daghestan are thought to have resided in the Caucasus Mountains for hundreds of generations and are characterized by a high prevalence of certain complex diseases. To explore the extent to which their unique population history led to increased levels of inbreeding, we genotyped >550 000 autosomal single-nucleotide polymorphisms (SNPs) in a set of 14 population isolates speaking Nakh-Daghestanian (ND) languages. The ND-speaking populations showed greatly elevated coefficients of inbreeding, very high numbers and long lengths of Runs of Homozygosity, and elevated linkage disequilibrium compared with surrounding groups from the Caucasus, the Near East, Europe, Central and South Asia. These results are consistent with the hypothesis that most ND-speaking groups descend from a common ancestral population that fragmented into a series of genetic isolates in the Daghestanian highlands. They have subsequently maintained a long-term small effective population size as a result of constant inbreeding and very low levels of gene flow. Given these findings, Daghestanian population isolates are likely to be useful for mapping genes associated with complex diseases.

Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data

BACKGROUND

The current availability of genotypes for very large numbers of single nucleotide polymorphisms (SNPs) is leading to more accurate estimates of inbreeding coefficients and more detailed approaches for detecting inbreeding depression. In the present study, genome-wide information was used to detect inbreeding depression for two reproductive traits (total number of piglets born and number of piglets born alive) in an ancient strain of Iberian pigs (the Guadyerbas strain) that is currently under serious danger of extinction.

METHODS

A total of 109 sows with phenotypic records were genotyped with the PorcineSNP60 BeadChip v1. Inbreeding depression was estimated using a bivariate animal model in which the inbreeding coefficient was included as a covariate. We used two different measures of genomic inbreeding to perform the analyses: inbreeding estimated on a SNP-by-SNP basis and inbreeding estimated from runs of homozygosity. We also performed the analyses using pedigree-based inbreeding.

RESULTS

Significant inbreeding depression was detected for both traits using all three measures of inbreeding. Genome-wide information allowed us to identify one region on chromosome 13 associated with inbreeding depression. This region spans from 27 to 54 Mb and overlaps with a previously detected quantitative trait locus and includes the inter-alpha-trypsin inhibitor gene cluster that is involved with embryo implantation.

CONCLUSIONS

Our results highlight the value of high-density SNP genotyping for providing new insights on where genes causing inbreeding depression are located in the genome. Genomic measures of inbreeding obtained on a SNP-by-SNP basis or those based on the presence/absence of runs of homozygosity represent a suitable alternative to pedigree-based measures to detect inbreeding depression, and a useful tool for mapping studies. To our knowledge, this is the first study in domesticated animals using the SNP-by-SNP inbreeding coefficient to map specific regions within chromosomes associated with inbreeding depression.

Quantifying realized inbreeding in wild and captive animal populations

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression-namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient FPed obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0-0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (FPed=0.013<propIBD=0.064), as expected given that pedigree founders were already related. We suggest that this SNP-based technique is generally useful for quantifying inbreeding at the individual or population level, and we show analytically that it can capture inbreeding loops that reach back up to a few hundred generations.

Effect of inbreeding and endogamy on occlusal traits in human isolates

OBJECTIVE

To discuss the genetic basis of occlusal traits through analysis of the effects of inbreeding in a subdivided isolated community.

SUBJECTS AND METHODS

The sample comprised dental casts of 224 children, aged 7-14 years, from 15 villages of the Island of Hvar, Croatia.

MAIN OUTCOME MEASURES

Studied traits were Angle class, overjet, vertical bite, overbite, and crowding/spacing.

DESIGN

Children with complete grandparental endogamy (all four grandparents born in the village of residence of the examinee) were compared to children with incomplete grandparental endogamy. In addition, children resident in the group of villages with a high prevalence of inbreeding were compared to children resident in the groups of villages with moderate and low prevalence of inbreeding.

RESULTS

In both designs, inbreeding seemed to increase the mean values of overjet, overbite, and vertical bite, while it had little or no effect on crowding/spacing. Angle classes were correlated to inbreeding at the individual level, but this was not supported at the population level. The effects were stronger in the subsample with bilaterally concordant Angle classes.

CONCLUSION

The observed inbreeding effects imply that the genetic basis of some occlusal traits is polygenic and, in considerable part, influenced by rare and recessive genetic variants.

Consanguinity and rare mutations outside of MCCC genes underlie nonspecific phenotypes of MCCD

Purpose

3-Methylcrotonyl-CoA carboxylase deficiency (MCCD) is an autosomal recessive disorder of leucine catabolism that has a highly variable clinical phenotype, ranging from acute metabolic acidosis to nonspecific symptoms such as developmental delay, failure to thrive, hemiparesis, muscular hypotonia, and multiple sclerosis. Implementation of newborn screening for MCCD has resulted in broadening the range of phenotypic expression to include asymptomatic adults. The purpose of this study was to identify factors underlying the varying phenotypes of MCCD.

Methods

We performed exome sequencing on DNA from 33 cases and 108 healthy controls. We examined these data for associations between either MCC mutational status, genetic ancestry, or consanguinity and the absence or presence/specificity of clinical symptoms in MCCD cases.

Results

We determined that individuals with nonspecific clinical phenotypes are highly inbred compared with cases that are asymptomatic and healthy controls. For 5 of these 10 individuals, we discovered a homozygous damaging mutation in a disease gene that is likely to underlie their nonspecific clinical phenotypes previously attributed to MCCD.

Conclusion

Our study shows that nonspecific phenotypes attributed to MCCD are associated with consanguinity and are likely not due to mutations in the MCC enzyme but result from rare homozygous mutations in other disease genes.

A novel LIPE nonsense mutation found using exome sequencing in siblings with late-onset familial partial lipodystrophy

BACKGROUND

Familial lipodystrophies are rare inherited disorders associated with redistribution of body fat and development of dyslipidemia, insulin resistance, and diabetes. We previously reported 2 siblings with unusual late-onset familial partial lipodystrophy in whom heretofore known causative genes had been excluded. We hypothesized they had a mutation in a novel lipodystrophy gene.

METHODS

Our approach centred on whole exome sequencing of the patients' DNA, together with genetic linkage analysis and a bioinformatic prioritization analysis. All candidate variants were assessed in silico and available family members were genotyped to assess segregation of mutations.

RESULTS

Our prioritization algorithm led us to a novel homozygous nonsense variant, namely p.Ala507fsTer563 in the hormone sensitive lipase gene encoding, an enzyme that is differentially expressed in adipocytes and steroidogenic tissues. Pathogenicity of the mutation was supported in bioinformatic analyses and variant cosegregation within the family.

CONCLUSIONS

We have identified a novel nonsense mutation in hormone sensitive lipase gene, which likely explains the lipodystrophy phenotype observed in these patients.

Identity-by-descent approaches identify regions of importance for genetic susceptibility to hereditary esophageal squamous cell carcinoma

Worldwide, the highest prevalence of esophageal cancer (EC) occurs in Northern China. High-density SNP arrays allow identification of identity-by-descent (IBD) segments in genomic DNAs representative of shared common ancestral regions. We utilized IBD approaches to map susceptibility loci associated with low-penetrance SNPs in high-risk Henan hereditary esophageal squamous cell carcinoma (ESCC) patients. Affymetrix GeneChip Human mapping SNP array IBD analysis was performed in 32 Henan family history-positive (FH+) ESCC patients, 18 Henan healthy unrelated individuals, and 45 Chinese individuals from a CHB HapMap dataset using PLink (scoring IBD segments individually) and Beagle (scoring of shared IBD segments among case/case vs. control/control pairs) software. Both analyses identified longer IBD segment lengths associated with FH+ ESCC compared to controls. However, there was no strong evidence for a genetic founder effect. Pairing IBD analysis with BEAGLE identified 8 critical IBD segments residing at 2q32.1-q32.2, 3p22.3-p22.2, 4q21.1-q21.21, 7p22.2, 8q23.2-q23.3, 10q23.33-q24.1, 14q24.3 and 16q11.2-q12.1, which were more significantly shared among case/case compared to control/control. The shared IBD segments in FH+ ESCC samples with no overlap with control/CHB Hapmap may encompass potential cancer susceptibility loci. Selected targeted genes, PLCE1, GPT2, SIAH1 and CYP2C-18, residing within the IBD segments at 10q23.33-q24.1 and 16q11.2-q12.1, had statistically significant differential expression in primary ESCC tissues and are likely involved in ESCC carcinogenesis. The importance of these IBD segments to the etiology and development of ESCC in high-risk areas requires further study with expanded sample sizes. This is the first report employing the pairing IBD approach for elucidation of the genetic basis of hereditary ESCC in Henan by applying high throughput SNP array analysis.

ILDR1: Novel mutation and a rare cause of congenital deafness in the Saudi Arabian population

Hearing impairment is the common human sensorineural disorder and is a genetically heterogeneous phenotype for which more than 100 genomic loci have been mapped so far. ILDR1 located on chromosome 3q13.33, encodes a putative transmembrane receptor containing an immunoglobulin-like domain. We used a combination of autozygosity mapping and candidate gene sequencing to identify a novel mutation in ILDR1, as a causative gene for autosomal-recessive non-syndromic hearing loss (arNSHL) in a consanguineous Saudi family with three affected children. Autozygosity mapping identified a shared region between the affected individuals encompassing ILDR1 on chromosome 3q13.12-3q22.1. Sequencing revealed homozygous 9 base pair duplication, resulting in an in-frame duplication of three amino acids p.(Asn109_Pro111dup). The mutation was segregating with the disease phenotype and is predicted to be pathogenic by SIFT and PROVEAN. The identified mutation is located in the immunoglobulin-type domain of the ILDR1 protein. In silico analysis using I-TASSER server and PyMOL offers the first predictions on the structural and functional consequences of this mutation. To our knowledge, this is the first ILDR1 mutation identified in a Saudi family. Identification of ILDR1 mutation in only one of 100 Saudi familial and sporadic individuals with hearing loss suggests that this mutation is unique to this family and that ILDR1 should be considered as a rare cause of congenital deafness among Saudi Arabian population. Our data also confirms the evidence for ILDR1 allelic heterogeneity and expands the number of familial arNSHL-associated ILDR1 gene mutations.

Allele frequency changes due to hitch-hiking in genomic selection programs

BACKGROUND

Genomic selection makes it possible to reduce pedigree-based inbreeding over best linear unbiased prediction (BLUP) by increasing emphasis on own rather than family information. However, pedigree inbreeding might not accurately reflect loss of genetic variation and the true level of inbreeding due to changes in allele frequencies and hitch-hiking. This study aimed at understanding the impact of using long-term genomic selection on changes in allele frequencies, genetic variation and level of inbreeding.

METHODS

Selection was performed in simulated scenarios with a population of 400 animals for 25 consecutive generations. Six genetic models were considered with different heritabilities and numbers of QTL (quantitative trait loci) affecting the trait. Four selection criteria were used, including selection on own phenotype and on estimated breeding values (EBV) derived using phenotype-BLUP, genomic BLUP and Bayesian Lasso. Changes in allele frequencies at QTL, markers and linked neutral loci were investigated for the different selection criteria and different scenarios, along with the loss of favourable alleles and the rate of inbreeding measured by pedigree and runs of homozygosity.

RESULTS

For each selection criterion, hitch-hiking in the vicinity of the QTL appeared more extensive when accuracy of selection was higher and the number of QTL was lower. When inbreeding was measured by pedigree information, selection on genomic BLUP EBV resulted in lower levels of inbreeding than selection on phenotype BLUP EBV, but this did not always apply when inbreeding was measured by runs of homozygosity. Compared to genomic BLUP, selection on EBV from Bayesian Lasso led to less genetic drift, reduced loss of favourable alleles and more effectively controlled the rate of both pedigree and genomic inbreeding in all simulated scenarios. In addition, selection on EBV from Bayesian Lasso showed a higher selection differential for mendelian sampling terms than selection on genomic BLUP EBV.

CONCLUSIONS

Neutral variation can be shaped to a great extent by the hitch-hiking effects associated with selection, rather than just by genetic drift. When implementing long-term genomic selection, strategies for genomic control of inbreeding are essential, due to a considerable hitch-hiking effect, regardless of the method that is used for prediction of EBV.

A homozygous SLITRK6 nonsense mutation is associated with progressive auditory neuropathy in humans

OBJECTIVES/HYPOTHESIS

SLITRK family proteins control neurite outgrowth and regulate synaptic development. In mice, Slitrk6 plays a role in the survival and innervation of sensory neurons in the inner ear, vestibular apparatus, and retina, and also influences axial eye length. We provide the first detailed description of the auditory phenotype in humans with recessive SLITRK6 deficiency.

STUDY DESIGN

Prospective observational case study.

METHODS

Nine closely related Amish subjects from an endogamous Amish community of Pennsylvania underwent audiologic and vestibular testing. Single nucleotide polymorphism microarrays were used to map the chromosome locus, and Sanger sequencing or high-resolution melt analysis were used to confirm the allelic variant.

RESULTS

All nine subjects were homozygous for a novel nonsense variant of SLITRK6 (c.1240C>T, p.Gln414Ter). Adult patients had high myopia. The 4 oldest SLITRK6 c.1240C>T homozygotes had absent ipsilateral middle ear muscle reflexes (MEMRs). Distortion product otoacoustic emissions (DPOAEs) were absent in all ears tested and the cochlear microphonic (CM) was increased in amplitude and duration in young patients and absent in the two oldest subjects. Auditory brainstem responses (ABRs) were dys-synchronised bilaterally with no reproducible waves I, III, or V at high intensities. Hearing loss and speech reception thresholds deteriorated symmetrically with age, which resulted in severe-to-profound hearing impairment by early adulthood. Vestibular evoked myogenic potentials were normal in three ears and absent in one.

CONCLUSION

Homozygous SLITRK6 c.1240C>T (p.Gln414Ter) nonsense mutations are associated with high myopia, cochlear dysfunction attributed to outer hair cell disease, and progressive auditory neuropathy.

Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors

BACKGROUND

Runs of homozygosity are long, uninterrupted stretches of homozygous genotypes that enable reliable estimation of levels of inbreeding (i.e., autozygosity) based on high-throughput, chip-based single nucleotide polymorphism (SNP) genotypes. While the theoretical definition of runs of homozygosity is straightforward, their empirical identification depends on the type of SNP chip used to obtain the data and on a number of factors, including the number of heterozygous calls allowed to account for genotyping errors. We analyzed how SNP chip density and genotyping errors affect estimates of autozygosity based on runs of homozygosity in three cattle populations, using genotype data from an SNP chip with 777,972 SNPs and a 50 k chip.

RESULTS

Data from the 50 k chip led to overestimation of the number of runs of homozygosity that are shorter than 4 Mb, since the analysis could not identify heterozygous SNPs that were present on the denser chip. Conversely, data from the denser chip led to underestimation of the number of runs of homozygosity that were longer than 8 Mb, unless the presence of a small number of heterozygous SNP genotypes was allowed within a run of homozygosity.

CONCLUSIONS

We have shown that SNP chip density and genotyping errors introduce patterns of bias in the estimation of autozygosity based on runs of homozygosity. SNP chips with 50,000 to 60,000 markers are frequently available for livestock species and their information leads to a conservative prediction of autozygosity from runs of homozygosity longer than 4 Mb. Not allowing heterozygous SNP genotypes to be present in a homozygosity run, as has been advocated for human populations, is not adequate for livestock populations because they have much higher levels of autozygosity and therefore longer runs of homozygosity. When allowing a small number of heterozygous calls, current software does not differentiate between situations where these calls are adjacent and therefore indicative of an actual break of the run versus those where they are scattered across the length of the homozygous segment. Simple graphical tests that are used in this paper are a current, yet tedious solution.

Long runs of homozygosity are enriched for deleterious variation

Exome sequencing offers the potential to study the population-genomic variables that underlie patterns of deleterious variation. Runs of homozygosity (ROH) are long stretches of consecutive homozygous genotypes probably reflecting segments shared identically by descent as the result of processes such as consanguinity, population size reduction, and natural selection. The relationship between ROH and patterns of predicted deleterious variation can provide insight into the way in which these processes contribute to the maintenance of deleterious variants. Here, we use exome sequencing to examine ROH in relation to the distribution of deleterious variation in 27 individuals of varying levels of apparent inbreeding from 6 human populations. A significantly greater fraction of all genome-wide predicted damaging homozygotes fall in ROH than would be expected from the corresponding fraction of nondamaging homozygotes in ROH (p < 0.001). This pattern is strongest for long ROH (p < 0.05). ROH, and especially long ROH, harbor disproportionately more deleterious homozygotes than would be expected on the basis of the total ROH coverage of the genome and the genomic distribution of nondamaging homozygotes. The results accord with a hypothesis that recent inbreeding, which generates long ROH, enables rare deleterious variants to exist in homozygous form. Thus, just as inbreeding can elevate the occurrence of rare recessive diseases that represent homozygotes for strongly deleterious mutations, inbreeding magnifies the occurrence of mildly deleterious variants as well.

A homozygous missense mutation in HERC2 associated with global developmental delay and autism spectrum disorder

We studied a unique phenotype of cognitive delay, autistic behavior, and gait instability segregating in three separate sibships. We initiated genome-wide mapping in two sibships using Affymetrix 10K SNP Mapping Arrays and identified a homozygous 8.2 Mb region on chromosome 15 common to five affected children. We used exome sequencing of two affected children to assess coding sequence variants within the mapped interval. Four novel homozygous exome variants were shared between the two patients; however, only two variants localized to the mapped interval on chromosome 15. A third sibship in an Ohio Amish deme narrowed the mapped interval to 2.6 Mb and excluded one of the two novel homozygous exome variants. The remaining variant, a missense change in HERC2 (c.1781C>T, p.Pro594Leu), occurs in a highly conserved proline residue within an RCC1-like functional domain. Functional studies of truncated HERC2 in adherent retinal pigment epithelium cells suggest that the p.Pro594Leu variant induces protein aggregation and leads to decreased HERC2 abundance. The phenotypic correlation with the mouse Herc1 and Herc2 mutants as well as the phenotypic overlap with Angelman syndrome provide further evidence that pathogenic changes in HERC2 are associated with nonsyndromic intellectual disability, autism, and gait disturbance. Hum Mutat 33:1639-1646, 2012. © 2012 Wiley Periodicals, Inc.

Unexpected relationships and inbreeding in HapMap phase III populations

Correct annotation of the genetic relationships between samples is essential for population genomic studies, which could be biased by errors or omissions. To this end, we used identity-by-state (IBS) and identity-by-descent (IBD) methods to assess genetic relatedness of individuals within HapMap phase III data. We analyzed data from 1,397 individuals across 11 ethnic populations. Our results support previous studies (Pemberton et al., 2010; Kyriazopoulou-Panagiotopoulou et al., 2011) assessing unknown relatedness present within this population. Additionally, we present evidence for 1,657 novel pairwise relationships across 9 populations. Surprisingly, significant Cotterman's coefficients of relatedness K1 (IBD1) values were detected between pairs of known parents. Furthermore, significant K2 (IBD2) values were detected in 32 previously annotated parent-child relationships. Consistent with a hypothesis of inbreeding, regions of homozygosity (ROH) were identified in the offspring of related parents, of which a subset overlapped those reported in previous studies (Gibson et al. 2010; Johnson et al. 2011). In total, we inferred 28 inbred individuals with ROH that overlapped areas of relatedness between the parents and/or IBD2 sharing at a different genomic locus between a child and a parent. Finally, 8 previously annotated parent-child relationships had unexpected K0 (IBD0) values (resulting from a chromosomal abnormality or genotype error), and 10 previously annotated second-degree relationships along with 38 other novel pairwise relationships had unexpected IBD2 (indicating two separate paths of recent ancestry). These newly described types of relatedness may impact the outcome of previous studies and should inform the design of future studies relying on the HapMap Phase III resource.

Length distributions of identity by descent reveal fine-scale demographic history

Data-driven studies of identity by descent (IBD) were recently enabled by high-resolution genomic data from large cohorts and scalable algorithms for IBD detection. Yet, haplotype sharing currently represents an underutilized source of information for population-genetics research. We present analytical results on the relationship between haplotype sharing across purportedly unrelated individuals and a population's demographic history. We express the distribution of IBD sharing across pairs of individuals for segments of arbitrary length as a function of the population's demography, and we derive an inference procedure to reconstruct such demographic history. The accuracy of the proposed reconstruction methodology was extensively tested on simulated data. We applied this methodology to two densely typed data sets: 500 Ashkenazi Jewish (AJ) individuals and 56 Kenyan Maasai (MKK) individuals (HapMap 3 data set). Reconstructing the demographic history of the AJ cohort, we recovered two subsequent population expansions, separated by a severe founder event, consistent with previous analysis of lower-throughput genetic data and historical accounts of AJ history. In the MKK cohort, high levels of cryptic relatedness were detected. The spectrum of IBD sharing is consistent with a demographic model in which several small-sized demes intermix through high migration rates and result in enrichment of shared long-range haplotypes. This scenario of historically structured demographies might explain the unexpected abundance of runs of homozygosity within several populations.

Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations

Using genome-wide SNP data, we calculated genomic inbreeding coefficients (FROH  > 1  Mb , FROH  > 2 Mb , FROH  > 8 Mb and FROH  > 16 Mb ) derived from runs of homozygosity (ROH) of different lengths (>1, >2, >8 and > 16 Mb) as well as from levels of homozygosity (FHOM ). We compared these values of inbreeding coefficients with those calculated from pedigrees (FPED ) of 1422 bulls comprising Brown Swiss (304), Fleckvieh (502), Norwegian Red (499) and Tyrol Grey (117) cattle breeds. For all four breeds, population inbreeding levels estimated by the genomic inbreeding coefficients FROH  > 8 Mb and FROH  > 16 Mb were similar to the levels estimated from pedigrees. The lowest values were obtained for Fleckvieh (FPED  = 0.014, FROH  > 8 Mb  = 0.019 and FROH  > 16 Mb  = 0.008); the highest, for Brown Swiss (FPED  = 0.048, FROH  > 8 Mb  = 0.074 and FROH  > 16 Mb  = 0.037). In contrast, inbreeding estimates based on the genomic coefficients FROH  > 1 Mb and FROH  > 2 Mb were considerably higher than pedigree-derived estimates. Standard deviations of genomic inbreeding coefficients were, on average, 1.3-1.7-fold higher than those obtained from pedigrees. Pearson correlations between genomic and pedigree inbreeding coefficients ranged from 0.50 to 0.62 in Norwegian Red (lowest correlations) and from 0.64 to 0.72 in Tyrol Grey (highest correlations). We conclude that the proportion of the genome present in ROH provides a good indication of inbreeding levels and that analysis based on ROH length can indicate the relative amounts of autozygosity due to recent and remote ancestors.

Genomic patterns of homozygosity in worldwide human populations

Genome-wide patterns of homozygosity runs and their variation across individuals provide a valuable and often untapped resource for studying human genetic diversity and evolutionary history. Using genotype data at 577,489 autosomal SNPs, we employed a likelihood-based approach to identify runs of homozygosity (ROH) in 1,839 individuals representing 64 worldwide populations, classifying them by length into three classes-short, intermediate, and long-with a model-based clustering algorithm. For each class, the number and total length of ROH per individual show considerable variation across individuals and populations. The total lengths of short and intermediate ROH per individual increase with the distance of a population from East Africa, in agreement with similar patterns previously observed for locus-wise homozygosity and linkage disequilibrium. By contrast, total lengths of long ROH show large interindividual variations that probably reflect recent inbreeding patterns, with higher values occurring more often in populations with known high frequencies of consanguineous unions. Across the genome, distributions of ROH are not uniform, and they have distinctive continental patterns. ROH frequencies across the genome are correlated with local genomic variables such as recombination rate, as well as with signals of recent positive selection. In addition, long ROH are more frequent in genomic regions harboring genes associated with autosomal-dominant diseases than in regions not implicated in Mendelian diseases. These results provide insight into the way in which homozygosity patterns are produced, and they generate baseline homozygosity patterns that can be used to aid homozygosity mapping of genes associated with recessive diseases.

Regions of homozygosity identified by SNP microarray analysis aid in the diagnosis of autosomal recessive disease and incidentally detect parental blood relationships

PURPOSE

The purpose of this study was to document the ability of single-nucleotide polymorphism microarray to identify copy-neutral regions of homozygosity, demonstrate clinical utility of regions of homozygosity, and discuss ethical/legal implications when regions of homozygosity are associated with a parental blood relationship.

METHODS

Study data were compiled from consecutive samples sent to our clinical laboratory over a 3-year period. A cytogenetics database identified patients with at least two regions of homozygosity >10 Mb on two separate chromosomes. A chart review was conducted on patients who met the criteria.

RESULTS

Of 3,217 single-nucleotide polymorphism microarrays, 59 (1.8%) patients met inclusion criteria. The percentage of homozygosity ranged from 0.9 to 30.1%, indicating parental relationships from distant to first-degree relatives. First-degree kinship was suspected in the parents of at least 11 patients with regions of homozygosity covering >21.3% of their autosome. In four patients from two families, homozygosity mapping discovered a candidate gene that was sequenced to identify a clinically significant mutation.

CONCLUSION

This study demonstrates clinical utility in the identification of regions of homozygosity, as these regions may aid in diagnosis of the patient. This study establishes the need for careful reporting, thorough pretest counseling, and careful electronic documentation, as microarray has the capability of detecting previously unknown/unreported relationships.

Evidence for additional FREM1 heterogeneity in Manitoba oculotrichoanal syndrome

PURPOSE

Manitoba Oculotrichoanal (MOTA) syndrome is an autosomal recessive disorder present in First Nations families that is characterized by ocular (cryptophthalmos), facial, and genital anomalies. At the commencement of this study, its genetic basis was undefined.

METHODS

Homozygosity analysis was employed to map the causative locus using DNA samples from four probands of Cree ancestry. After single nucleotide polymorphism (SNP) genotyping, data were analyzed and exported to PLINK to identify regions identical by descent (IBD) and common to the probands. Candidate genes within and adjacent to the IBD interval were sequenced to identify pathogenic variants, with analyses of potential deletions or duplications undertaken using the B-allele frequency and log(2) ratio of SNP signal intensity.

RESULTS

Although no shared IBD region >1 Mb was evident on preliminary analysis, adjusting the criteria to permit the detection of smaller homozygous IBD regions revealed one 330 Kb segment on chromosome 9p22.3 present in all 4 probands. This interval comprising 152 SNPs, lies 16 Kb downstream of FRAS1-related extracellular matrix protein 1 (FREM1), and no copy number variations were detected either in the IBD region or FREM1. Subsequent sequencing of both genes in the IBD region, followed by FREM1, did not reveal any mutations.

CONCLUSIONS

This study illustrates the utility of studying geographically isolated populations to identify genomic regions responsible for disease through analysis of small numbers of affected individuals. The location of the IBD region 16 kb from FREM1 suggests the phenotype in these patients is attributable to a variant outside of FREM1, potentially in a regulatory element, whose identification may prove tractable to next generation sequencing. In the context of recent identification of FREM1 coding mutations in a proportion of MOTA cases, characterization of such additional variants offers scope both to enhance understanding of FREM1's role in cranio-facial biology and may facilitate genetic counselling in populations with high prevalences of MOTA to reduce the incidence of this disorder.

Human genetics and genomics a decade after the release of the draft sequence of the human genome

Substantial progress has been made in human genetics and genomics research over the past ten years since the publication of the draft sequence of the human genome in 2001. Findings emanating directly from the Human Genome Project, together with those from follow-on studies, have had an enormous impact on our understanding of the architecture and function of the human genome. Major developments have been made in cataloguing genetic variation, the International HapMap Project, and with respect to advances in genotyping technologies. These developments are vital for the emergence of genome-wide association studies in the investigation of complex diseases and traits. In parallel, the advent of high-throughput sequencing technologies has ushered in the 'personal genome sequencing' era for both normal and cancer genomes, and made possible large-scale genome sequencing studies such as the 1000 Genomes Project and the International Cancer Genome Consortium. The high-throughput sequencing and sequence-capture technologies are also providing new opportunities to study Mendelian disorders through exome sequencing and whole-genome sequencing. This paper reviews these major developments in human genetics and genomics over the past decade.

A genome-wide homozygosity association study identifies runs of homozygosity associated with rheumatoid arthritis in the human major histocompatibility complex

Rheumatoid arthritis (RA) is a chronic inflammatory disorder with a polygenic mode of inheritance. This study examined the hypothesis that runs of homozygosity (ROHs) play a recessive-acting role in the underlying RA genetic mechanism and identified RA-associated ROHs. Ours is the first genome-wide homozygosity association study for RA and characterized the ROH patterns associated with RA in the genomes of 2,000 RA patients and 3,000 normal controls of the Wellcome Trust Case Control Consortium. Genome scans consistently pinpointed two regions within the human major histocompatibility complex region containing RA-associated ROHs. The first region is from 32,451,664 bp to 32,846,093 bp (−log10(p)>22.6591). RA-susceptibility genes, such as HLA-DRB1, are contained in this region. The second region ranges from 32,933,485 bp to 33,585,118 bp (−log10(p)>8.3644) and contains other HLA-DPA1 and HLA-DPB1 genes. These two regions are physically close but are located in different blocks of linkage disequilibrium, and ∼40% of the RA patients' genomes carry these ROHs in the two regions. By analyzing homozygote intensities, an ROH that is anchored by the single nucleotide polymorphism rs2027852 and flanked by HLA-DRB6 and HLA-DRB1 was found associated with increased risk for RA. The presence of this risky ROH provides a 62% accuracy to predict RA disease status. An independent genomic dataset from 868 RA patients and 1,194 control subjects of the North American Rheumatoid Arthritis Consortium successfully validated the results obtained using the Wellcome Trust Case Control Consortium data. In conclusion, this genome-wide homozygosity association study provides an alternative to allelic association mapping for the identification of recessive variants responsible for RA. The identified RA-associated ROHs uncover recessive components and missing heritability associated with RA and other autoimmune diseases.

Consanguinity and the risk of congenital heart disease

Consanguineous unions have been associated with an increased susceptibility to various forms of inherited disease. Although consanguinity is known to contribute to recessive diseases, the potential role of consanguinity in certain common birth defects is less clear, particularly since the disease pathophysiology may involve genetic and environmental/epigenetic factors. In this study, we ask whether consanguinity affects one of the most common birth defects, congenital heart disease, and identify areas for further research into these birth defects, since consanguinity may now impact health on a near-global basis. A systematic review of consanguinity in congenital heart disease was performed, focusing on non-syndromic disease, with the methodologies and results from studies of different ethnic populations compared. The risks for congenital heart disease have been assessed and summarized collectively and by individual lesion. The majority of studies support the view that consanguinity increases the prevalence of congenital heart disease, however, the study designs differed dramatically. Only a few (n = 3) population-based studies that controlled for potential sociodemographic confounding were identified, and data on individual cardiac lesions were limited by case numbers. Overall the results suggest that the risk for congenital heart disease is increased in consanguineous unions in the studied populations, principally at first-cousin level and closer, a factor that should be considered in empiric risk estimates in genetic counseling. However, for more precise risk estimates a better understanding of the underlying disease factors is needed.

Primary immunodeficiencies in highly consanguineous North African populations

The study of inbred populations has contributed remarkably to the description of new autosomal recessive primary immunodeficiencies (PIDs). Here, we examine the pattern of PIDs in North African populations and assess the impact of highly prevalent consanguinity. This review reports on the current status of pediatricians' awareness of PIDs in Egypt, Morocco, and Tunisia, where awareness of PIDs is relatively recent. The phenotypic distribution of PIDs is reported and compared among the three countries and with other populations. Data analysis reveals a prevalence of autosomal recessive forms and a peculiar distribution of major PID categories, particularly more combined immunodeficiencies than antibody disorders. In these endogamous communities, molecular diagnosis is critical to developing a genetic-based preventive approach. The organization of diagnosis and care services in these resource-limited settings faces many obstacles. Autosomal recessive PIDs are overrepresented; thus, it is critical to continue investigation of these diseases in order to better understand the underlying mechanisms and to improve patient care.

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.

Genetic mapping and exome sequencing identify variants associated with five novel diseases

The Clinic for Special Children (CSC) has integrated biochemical and molecular methods into a rural pediatric practice serving Old Order Amish and Mennonite (Plain) children. Among the Plain people, we have used single nucleotide polymorphism (SNP) microarrays to genetically map recessive disorders to large autozygous haplotype blocks (mean = 4.4 Mb) that contain many genes (mean = 79). For some, uninformative mapping or large gene lists preclude disease-gene identification by Sanger sequencing. Seven such conditions were selected for exome sequencing at the Broad Institute; all had been previously mapped at the CSC using low density SNP microarrays coupled with autozygosity and linkage analyses. Using between 1 and 5 patient samples per disorder, we identified sequence variants in the known disease-causing genes SLC6A3 and FLVCR1, and present evidence to strongly support the pathogenicity of variants identified in TUBGCP6, BRAT1, SNIP1, CRADD, and HARS. Our results reveal the power of coupling new genotyping technologies to population-specific genetic knowledge and robust clinical data.

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.

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."

Identification of possible genetic alterations in the breast cancer cell line MCF-7 using high-density SNP genotyping microarray

Context:

Cancer cell lines are used extensively in various research. Knowledge of genetic alterations in these lines is important for understanding mechanisms underlying their biology. However, since paired normal tissues are usually unavailable for comparison, precisely determining genetic alterations in cancer cell lines is difficult. To address this issue, a highly efficient and reliable method is developed. Aims: Establishing a highly efficient and reliable experimental system for genetic profiling of cell lines.

Materials and Methods:

A widely used breast cancer cell line, MCF-7, was genetically profiled with 4,396 single nucleotide polymorphisms (SNPs) spanning 11 whole chromosomes and two other small regions using a newly developed high-throughput multiplex genotyping approach.

Results:

The fractions of homozygous SNPs in MCF-7 (13.3%) were significantly lower than those in the control cell line and in 24 normal human individuals (25.1% and 27.4%, respectively). Homozygous SNPs in MCF-7 were found in clusters. The sizes of these clusters were significantly larger than the expected based on random allelic combination. Fourteen such regions were found on chromosomes 1p, 1q, 2q, 6q, 13, 15q, 16q, 17q and 18p in MCF-7 and two in the small regions.

Conclusions:

These results are generally concordant with those obtained using different approaches but are better in defining their chromosomal positions. The used approach provides a reliable way to detecting possible genetic alterations in cancer cell lines without paired normal tissues.

KinSNP software for homozygosity mapping of disease genes using SNP microarrays

Consanguineous families affected with a recessive genetic disease caused by homozygotisation of a mutation offer a unique advantage for positional cloning of rare diseases. Homozygosity mapping of patient genotypes is a powerful technique for the identification of the genomic locus harbouring the causing mutation. This strategy relies on the observation that in these patients a large region spanning the disease locus is also homozygous with high probability. The high marker density in single nucleotide polymorphism (SNP) arrays is extremely advantageous for homozygosity mapping.

We present KinSNP, a user-friendly software tool for homozygosity mapping using SNP arrays. The software searches for stretches of SNPs which are homozygous to the same allele in all ascertained sick individuals. User-specified parameters control the number of allowed genotyping 'errors' within homozygous blocks. Candidate disease regions are then reported in a detailed, coloured Excel file, along with genotypes of family members and healthy controls. An interactive genome browser has been included which shows homozygous blocks, individual genotypes, genes and further annotations along the chromosomes, with zooming and scrolling capabilities. The software has been used to identify the location of a mutated gene causing insensitivity to pain in a large Bedouin family. KinSNP is freely available from http://bioinfo.bgu.ac.il/bsu/software/kinSNP.

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.

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.

Reconciling the analysis of IBD and IBS in complex trait studies

Identity by descent (IBD) is a fundamental concept in genetics and refers to alleles that are descended from a common ancestor in a base population. Identity by state (IBS) simply refers to alleles that are the same, irrespective of whether they are inherited from a recent ancestor. In modern applications, IBD relationships are estimated from genetic markers in individuals without any known relationship. This can lead to erroneous inference because a consistent base population is not used. We argue that the purpose of most IBD calculations is to predict IBS at unobserved loci. Recognizing this aim leads to better methods to estimating IBD with benefits in mapping genes, estimating genetic variance and predicting inbreeding depression.

A novel autosomal recessive malformation syndrome associated with developmental delay and distinctive facies maps to 16ptel in the Hutterite population

The Hutterites are a genetically isolated Anabaptist group living on the North American prairies; their population numbers over 40,000, the majority of whom are descendants of 89 founders. An autosomal recessive developmental disorder was identified in four patients from two consanguineous Hutterite families. To our knowledge the clinical presentation is unique and undescribed. The patients have distinctive facial features, congenital malformations of the heart and genitourinary system, head circumference at the 2nd centile and developmental delay. The facial features include tall forehead with high anterior hairline, deep-set eyes with short, upslanted palpebral fissures, long nose with low-hanging columella, and thick vermilion of the upper and lower lip. Karyotype and baseline metabolic studies were normal. An identity-by-descent mapping approach was used to localize the gene for this disorder. The patients were genotyped using Affymetrix GeneChip Human Mapping 10 K (Xba 2.0) and 50 K (Xba 240) Arrays which identified a single 5.5 Mb homozygous region at chromosome 16p13.3. To confirm and refine the boundaries of this region, microsatellite markers were used to genotype the patients, their parents, and the available unaffected siblings. The disease locus was refined to a region of 5.1 Mb containing 173 known or predicted genes. No other recessive disorders with similar clinical features are currently mapped to this region. The coding regions of over fifteen genes, prioritized by microarray expression analysis and information available in public databases, have been sequenced, but no potential pathogenic mutations have been identified. The identification of the gene for this syndrome will provide new insights into development and learning.

Homozygosity mapping in patients with cone-rod dystrophy: novel mutations and clinical characterizations

PURPOSE

To determine the genetic defect and to describe the clinical characteristics in a cohort of mainly nonconsanguineous cone-rod dystrophy (CRD) patients.

METHODS

One hundred thirty-nine patients with diagnosed CRD were recruited. Ninety of them were screened for known mutations in ABCA4, and those carrying one or two mutations were excluded from further research. Genome-wide homozygosity mapping was performed in the remaining 108. Known genes associated with autosomal recessive retinal dystrophies located within a homozygous region were screened for mutations. Patients in whom a mutation was detected underwent further ophthalmic examination.

RESULTS

Homozygous sequence variants were identified in eight CRD families, six of which were nonconsanguineous. The variants were detected in the following six genes: ABCA4, CABP4, CERKL, EYS, KCNV2, and PROM1. Patients carrying mutations in ABCA4, CERKL, and PROM1 had typical CRD symptoms, but a variety of retinal appearances on funduscopy, optical coherence tomography, and autofluorescence imaging.

CONCLUSIONS

Homozygosity mapping led to the identification of new mutations in consanguineous and nonconsanguineous patients with retinal dystrophy. Detailed clinical characterization revealed a variety of retinal appearances, ranging from nearly normal to extensive retinal remodeling, retinal thinning, and debris accumulation. Although CRD was initially diagnosed in all patients, the molecular findings led to a reappraisal of the diagnosis in patients carrying mutations in EYS, CABP4, and KCNV2.

An integrated Bayesian analysis of LOH and copy number data

Background

Cancer and other disorders are due to genomic lesions. SNP-microarrays are able to measure simultaneously both genotype and copy number (CN) at several Single Nucleotide Polymorphisms (SNPs) along the genome. CN is defined as the number of DNA copies, and the normal is two, since we have two copies of each chromosome. The genotype of a SNP is the status given by the nucleotides (alleles) which are present on the two copies of DNA. It is defined homozygous or heterozygous if the two alleles are the same or if they differ, respectively. Loss of heterozygosity (LOH) is the loss of the heterozygous status due to genomic events.

Combining CN and LOH data, it is possible to better identify different types of genomic aberrations. For example, a long sequence of homozygous SNPs might be caused by either the physical loss of one copy or a uniparental disomy event (UPD), i.e. each SNP has two identical nucleotides both derived from only one parent. In this situation, the knowledge of the CN can help in distinguishing between these two events.

Results

To better identify genomic aberrations, we propose a method (called gBPCR) which infers the type of aberration occurred, taking into account all the possible influence in the microarray detection of the homozygosity status of the SNPs, resulting from an altered CN level. Namely, we model the distributions of the detected genotype, given a specific genomic alteration and we estimate the parameters involved on public reference datasets. The estimation is performed similarly to the modified Bayesian Piecewise Constant Regression, but with improved estimators for the detection of the breakpoints.

Using artificial and real data, we evaluate the quality of the estimation of gBPCR and we also show that it outperforms other well-known methods for LOH estimation.

Conclusions

We propose a method (gBPCR) for the estimation of both LOH and CN aberrations, improving their estimation by integrating both types of data and accounting for their relationships. Moreover, gBPCR performed very well in comparison with other methods for LOH estimation and the estimated CN lesions on real data have been validated with another technique.

Genome-wide homozygosity signatures and childhood acute lymphoblastic leukemia risk

Recent studies have reported that regions of homozygosity (ROH) in the genome are detectable in outbred populations and can be associated with an increased risk of malignancy. To examine whether homozygosity is associated with an increased risk of developing childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL), we analyzed 824 ALL cases and 2398 controls genotyped for 292 200 tagging SNPs. Across the genome, cumulative distribution of ROH was not significantly different between cases and controls. Four common ROH at 10p11.2-10q11.21, 1p31.1, 19p13.2-3, and 20q11.1-23 were, however, associated with ALL risk at P less than .01 (including 1 ROH to which the erythropoietin receptor [EPOR] gene maps, P = .005) but were nonsignificant after adjusting for multiple testing. Our findings make it unlikely that levels of measured homozygosity, caused by autozygosity, uniparental isodisomy, or hemizygosity, play a major role in defining BCP-ALL risk in predominantly outbred populations.

Genomic and geographic distribution of SNP-defined runs of homozygosity in Europeans

The availability of high-density panels of genetic polymorphisms has led to the discovery of extended regions of apparent autozygosity in the human genome. At the genotype level, these regions present as sizeable stretches, or 'runs', of homozygosity (ROH). Here, we investigated both the genomic and the geographic distribution of ROHs in a large European sample of individuals originating from 23 subpopulations. The genomic ROH distribution was found to be characterized by a pattern of highly significant non-uniformity that was virtually identical in all subpopulations studied. Some 77 chromosomal regions contained ROHs at considerable frequency, thereby forming 'ROH islands' that were not explicable by high linkage disequilibrium alone. At the geographic level, the number and cumulative length of ROHs followed a prominent South to North gradient in agreement with expectations from European population history. The individual ROH length, in contrast, showed only minor and unsystematic geographic variation. While our findings are thus consistent with a larger effective population size in Southern than in Northern Europe, combined with a higher historic population density and mobility, they also indicate that the patterns of meiotic recombination in humans must have been very similar throughout the continent. Extending previous reports of a strong correlation between geography and identity-by-state, our data show that the genomic identity-by-descent patterns of Europeans are also clinal. As a consequence, the planning, design and interpretation of ROH-based genetic studies must take sample origin into account in order for such studies to be sensible and valid.