Runs of homozygosity identify a recessive locus 12q21.31 for human adult height

Genetic ancestry plays a central role in population pharmacogenomics

Recent studies have pointed out the essential role of genetic ancestry in population pharmacogenetics. In this study, we analyzed the whole-genome sequencing data from The 1000 Genomes Project (Phase 3) and the pharmacogenetic information from Drug Bank, PharmGKB, PharmaADME, and Biotransformation. Here we show that ancestry-informative markers are enriched in pharmacogenetic loci, suggesting that trans-ancestry differentiation must be carefully considered in population pharmacogenetics studies. Ancestry-informative pharmacogenetic loci are located in both protein-coding and non-protein-coding regions, illustrating that a whole-genome analysis is necessary for an unbiased examination over pharmacogenetic loci. Finally, those ancestry-informative pharmacogenetic loci that target multiple drugs are often a functional variant, which reflects their importance in biological functions and pathways. In summary, we develop an efficient algorithm for an ultrahigh-dimensional principal component analysis. We create genetic catalogs of ancestry-informative markers and genes. We explore pharmacogenetic patterns and establish a high-accuracy prediction panel of genetic ancestry. Moreover, we construct a genetic ancestry pharmacogenomic database Genetic Ancestry PhD (http://hcyang.stat.sinica.edu.tw/databases/genetic_ancestry_phd/).

Hsin-Chou Yang et al. examine population structure in several genomic databases and identify that pharmacogenetic loci are enriched for markers of genetic ancestry. Their results suggest that genetic ancestry must be carefully considered in population pharmacogenetics studies.

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

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

A Perception on Genome-Wide Genetic Analysis of Metabolic Traits in Arab Populations

Despite dedicated nation-wide efforts to raise awareness against the harmful effects of fast-food consumption and sedentary lifestyle, the Arab population continues to struggle with an increased risk for metabolic disorders. Unlike the European population, the Arab population lacks well-established genetic risk determinants for metabolic disorders, and the transferability of established risk loci to this population has not been satisfactorily demonstrated. The most recent findings have identified over 240 genetic risk loci (with ~400 independent association signals) for type 2 diabetes, but thus far only 25 risk loci (ADAMTS9, ALX4, BCL11A, CDKAL1, CDKN2A/B, COL8A1, DUSP9, FTO, GCK, GNPDA2, HMG20A, HNF1A, HNF1B, HNF4A, IGF2BP2, JAZF1, KCNJ11, KCNQ1, MC4R, PPARγ, SLC30A8, TCF7L2, TFAP2B, TP53INP1, and WFS1) have been replicated in Arab populations. To our knowledge, large-scale population- or family-based association studies are non-existent in this region. Recently, we conducted genome-wide association studies on Arab individuals from Kuwait to delineate the genetic determinants for quantitative traits associated with anthropometry, lipid profile, insulin resistance, and blood pressure levels. Although these studies led to the identification of novel recessive variants, they failed to reproduce the established loci. However, they provided insights into the genetic architecture of the population, the applicability of genetic models based on recessive mode of inheritance, the presence of genetic signatures of inbreeding due to the practice of consanguinity, and the pleiotropic effects of rare disorders on complex metabolic disorders. This perspective presents analysis strategies and study designs for identifying genetic risk variants associated with diabetes and related traits in Arab populations.

Runs of homozygosity associate with decreased risks of lung cancer in never-smoking East Asian females

Although genome-wide association studies (GWASs) have identified some risk single-nucleotide polymorphisms in East Asian never-smoking females, the unexplained missing heritability is still required to be investigated. Runs of homozygosity (ROHs) are thought to be a type of genetic variation acting on human complex traits and diseases. We detected ROHs in 8,881 East Asian never-smoking women. The summed ROHs were used to fit a logistic regression model which noteworthily revealed a significant association between ROHs and the decreased risk of lung cancer (P < 0.05). We identified 4 common ROHs regions located at 2p22.1, which were significantly associated with decreased risk of lung cancer (P = 2.00 × 10^-4 - 1.35 × 10^-4). Functional annotation was conducted to investigate the regulatory function of ROHs. The common ROHs were overlapped with potential regulatory elements, such as active epigenome elements and chromatin states in lung-derived cell lines. SOS1 and ARHGEF33 were significantly up-regulated as the putative target genes of the identified ROHs in lung cancer samples according to the analysis of differently expressed genes. Our results suggest that ROHs could act as recessive contributing factors and regulatory elements to influence the risk of lung cancer in never-smoking East Asian females.

Comparison of multi-lineage differentiation of hiPSCs reveals novel miRNAs that regulate lineage specification

MicroRNAs (miRNAs) are known to be crucial players in governing the differentiation of human induced pluripotent stem cells (hiPSCs). Despite their utter importance, identifying key lineage specifiers among the myriads of expressed miRNAs remains challenging. We believe that the current practice in mining miRNA specifiers via delineating dynamic fold-changes only is inadequate. Our study, therefore, provides evidence to pronounce "lineage specificity" as another important attribute to qualify for these lineage specifiers. Adopted hiPSCs were differentiated into representative lineages (hepatic, nephric and neuronal) over all three germ layers whilst the depicted miRNA expression changes compiled into an integrated atlas. We demonstrated inter-lineage analysis shall aid in the identification of key miRNAs with lineage-specificity, while these shortlisted candidates were collectively known as "lineage-specific miRNAs". Subsequently, we followed through the fold-changes along differentiation via computational analysis to identify miR-192 and miR-372-3p, respectively, as representative candidate key miRNAs for the hepatic and nephric lineages. Indeed, functional characterization validated that miR-192 and miR-372-3p regulate lineage differentiation via modulation of the expressions of lineage-specific genes. In summary, our presented miRNA atlas is a resourceful ore for the mining of key miRNAs responsible for lineage specification.

GARLIC: Genomic Autozygosity Regions Likelihood-based Inference and Classification

Summary

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

Availability and Implementation

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

Contact

zachary.szpiech@ucsf.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Runs of homozygosity: windows into population history and trait architecture

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

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

Phenotypic Consequences of a Spontaneous Loss of Heterozygosity in a Common Laboratory Strain of Candida albicans

By testing the susceptibility to DNA damaging agents of several Candida albicans mutant strains derived from the commonly used laboratory strain, CAI4, we uncovered sensitivity to methyl methanesulfonate (MMS) in CAI4 and its derivatives, but not in CAF2-1. This sensitivity is not a result of URA3 disruption because the phenotype was not restored after URA3 reintroduction. Rather, we found that homozygosis of a short region of chromosome 3R (Chr3R), which is naturally heterozygous in the MMS-resistant-related strains CAF4-2 and CAF2-1, confers MMS sensitivity and modulates growth polarization in response to MMS. Furthermore, induction of homozygosity in this region in CAF2-1 or CAF4-2 resulted in MMS sensitivity. We identified 11 genes by SNP/comparative genomic hybridization containing only the a alleles in all the MMS-sensitive strains. Four candidate genes, SNF5, POL1, orf19.5854.1, and MBP1, were analyzed by generating hemizygous configurations in CAF2-1 and CAF4-2 for each allele of all four genes. Only hemizygous MBP1a/mbp1b::SAT1-FLIP strains became MMS sensitive, indicating that MBP1a in the homo- or hemizygosis state was sufficient to account for the MMS-sensitive phenotype. In yeast, Mbp1 regulates G1/S genes involved in DNA repair. A second region of homozygosis on Chr2L increased MMS sensitivity in CAI4 (Chr3R homozygous) but not CAF4-2 (Chr3R heterozygous). This is the first example of sign epistasis in C. albicans.

Evidence of Inbreeding in Hodgkin Lymphoma

Genome-wide association studies (GWASs) have identified several, mainly co-dominantly acting, single-nucleotide polymorphisms (SNPs) associated with Hodgkin lymphoma (HL). We searched for recessively acting disease loci by performing an analysis of runs of homozygosity (ROH) based on windows of homozygous SNP-blocks and by calculating genomic inbreeding coefficients on a SNP-wise basis. We used data from a previous GWAS with 906 cases and 1217 controls from a population with a long history of no matings between relatives. Ten recurrent ROHs were identified among 25 055 ROHs across all individuals but their association with HL was not genome-wide significant. All recurrent ROHs showed significant evidence for natural selection. As a novel finding genomic inbreeding among cases was significantly higher than among controls (P = 2.11*10-14) even after correcting for covariates. Higher inbreeding among the cases was mainly based on a group of individuals with a higher average length of ROHs per person. This result suggests a correlation of higher levels of inbreeding with higher cancer incidence and might reflect the existence of recessive alleles causing HL. Genomic inbreeding may result in a higher expression of deleterious recessive genes within a population.

Genetics and psychotic disorders: A fresh look at consanguinity

Consanguineous unions refer to marriages between related individuals who share a common ancestor. These unions are still commonplace in certain regions of the world such as the southern coast of the Mediterranean, throughout the Middle East and South-East Asia. According to available data, couples of second cousins or closer and their offspring currently represent 10.4% of the world's population, thus resulting in increased frequencies of autosomal recessive disorders. Furthermore, consanguinity may be implicated in the increased frequency of multifactorial pathologies such as mental disorders. The few existing epidemiological studies in consanguineous and/or geographically isolated populations confirm that there is a significant association between consanguinity and mental disorders and a higher risk of schizophrenia or bipolar disorders among offspring from consanguineous couples. There exists a strong and complex genetic component in the predisposition to psychotic disorders that has been confirmed in numerous studies. However, the genetic basis of these disorders remains poorly understood. GWAS studies (Genome Wide Association Studies) over the past 10 years have identified a few weak associations, thus refuting the "common diseases-common variants" hypothesis. A model implicating numerous rare variants has been supported by the recent discovery of CNVs (Copy Number Variants) and their statistically significant association with psychiatric disorders such as schizophrenia, bipolar disorders and autism. The study of consanguineous families may contribute to identifying rare variants in homogenous populations who have conserved certain alleles. Major developments in molecular biology techniques would facilitate these studies as well as contributing to identifying major genes. These results emphasize the need for genetic counseling in high-risk communities and the importance of implementing preventive actions and raising awareness concerning the risk of consanguineous unions.

Inbreeding and homozygosity in breast cancer survival

Genome-wide association studies (GWASs) help to understand the effects of single nucleotide polymorphisms (SNPs) on breast cancer (BC) progression and survival. We performed multiple analyses on data from a previously conducted GWAS for the influence of individual SNPs, runs of homozygosity (ROHs) and inbreeding on BC survival. (I.) The association of individual SNPs indicated no differences in the proportions of homozygous individuals among short-time survivors (STSs) and long-time survivors (LTSs). (II.) The analysis revealed differences among the populations for the number of ROHs per person and the total and average length of ROHs per person and among LTSs and STSs for the number of ROHs per person. (III.) Common ROHs at particular genomic positions were nominally more frequent among LTSs than in STSs. Common ROHs showed significant evidence for natural selection (iHS, Tajima's D, Fay-Wu's H). Most regions could be linked to genes related to BC progression or treatment. (IV.) Results were supported by a higher level of inbreeding among LTSs. Our results showed that an increased level of homozygosity may result in a preference of individuals during BC treatment. Although common ROHs were short, variants within ROHs might favor survival of BC and may function in a recessive manner.

Genome-wide Survey of Runs of Homozygosity Identifies Recessive Loci for Bone Mineral Density in Caucasian and Chinese Populations

Runs of homozygosity (ROHs), in which both parental alleles are identical, have been proposed to have recessive effects on multiple human complex diseases. Osteoporosis is a common complex disease characterized by low bone mineral density (BMD), which is highly heritable. And recessive loci that contribute to BMD variations have been identified. In this study, we performed a genome-wide ROHs association study using our SNP array data from three GWAS samples including 4,900 subjects from Caucasian and Chinese populations. Significant results were further subjected to replication in 3,747 additional subjects. ROHs associated with BMD were also tested for associations with osteoporotic fractures in a GWAS fracture sample. Combining results from all the samples, we identified 697 autosomal regions with ROHs. Among these, we detected genome-wide significant associations between BMD and 6 ROHs, including ROH1q31.3, 1p31.1, 3q26.1, 11q12.1, 21q22.1 and 15q22.3 (combined P = 6.29 × 10(-5)-3.17 × 10(-8)). Especially, ROH1p31.1 was found to be associated with increased risk of osteoporotic hip fractures (odds ratio [OR] = 3.71, P = 0.032). To investigate the functional relevance of the identified ROHs, we performed cis-expression quantitative trait locus (eQTL) analysis in lymphoblast cell lines. Three ROHs, including ROH1p31.1, 11q12.1, and 15q22.3, were found to be significantly associated with mRNA expression levels of their nearby genes (PeQTL < 0.05). In summary, our findings reveal that ROHs could play as recessive-acting determinants contributing to the pathogenesis of osteoporosis. Further molecular and functional studies are needed to explore and clarify the potential mechanism.

Estimation of Recent and Ancient Inbreeding in a Small Endogamous Tunisian Community Through Genomic Runs of Homozygosity

Runs of homozygosity (ROHs) are extended genomic regions of homozygous genotypes that record populations' mating patterns in the past. We performed microarray genotyping on 15 individuals from a small isolated Tunisian community. We estimated the individual and population genome-wide level of homozygosity from data on ROH above 0.5 Mb in length. We found a high average number of ROH per individual (48.2). The smallest ROH category (0.5-1.49 Mb) represents 0.93% of the whole genome, while medium-size (1.5-4.99 Mb) and long-size ROH (≥5 Mb) cover 1.18% and 0.95%, respectively. We found that genealogical individual inbreeding coefficients (Fped ) based on three- to four-generation pedigrees are not reliable indicators of the current proportion of genome-wide homozygosity inferred from ROH (FROH ) either for 0.5 or 1.5 Mb ROH length thresholds, while identity-by-descent sharing is a function of shared coancestry. This study emphasizes the effect of reproductive isolation and a prolonged practice of consanguinity that limits the genetic heterogeneity. It also provides evidence of both recent and ancient parental relatedness contribution to the current level of genome-wide homozygosity in the studied population. These findings may be useful for evaluation of long-term effects of inbreeding on human health and for future applications of ROHs in identifying recessive susceptibility genes.

Runs of homozygosity and a cluster of vulvar cancer in young Australian Aboriginal women

OBJECTIVE

A cluster of vulvar cancer exists in young Aboriginal women living in remote communities in Arnhem Land, Australia. A genetic case-control study was undertaken involving 30 cases of invasive vulvar cancer and its precursor lesion, high-grade vulvar intraepithelial neoplasia (VIN), and 61 controls, matched for age and community of residence. It was hypothesized that this small, isolated population may exhibit increased autozygosity, implicating recessive effects as a possible mechanism for increased susceptibility to vulvar cancer.

METHODS

Genotyping data from saliva samples were used to identify runs of homozygosity (ROH) in order to calculate estimates of genome-wide homozygosity.

RESULTS

No evidence of an effect of genome-wide homozygosity on vulvar cancer and VIN in East Arnhem women was found, nor was any individual ROH found to be significantly associated with case status. This study found further evidence supporting an association between previous diagnosis of CIN and diagnosis of vulvar cancer or VIN, but found no association with any other medical history variable.

CONCLUSIONS

These findings do not eliminate the possibility of genetic risk factors being involved in this cancer cluster, but rather suggest that alternative analytical strategies and genetic models should be explored.

Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity

Background

Hanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo.

Results

The Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. ‘Regulation of actin filament length’ was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus.

Conclusions

We identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.

Could inbred cases identified in GWAS data succeed in detecting rare recessive variants where affected sib-pairs have failed?

To detect fully penetrant rare recessive variants that could constitute Mendelian subentities of complex diseases, we propose a novel strategy, the HBD-GWAS strategy, which can be applied to genome-wide association study (GWAS) data. This strategy first involves the identification of inbred individuals among cases using the genome-wide SNP data and then focuses on these inbred affected individuals and searches for genomic regions of shared homozygosity by descent that could harbor rare recessive disease-causing variants. In this second step, analogous to homozygosity mapping, a heterogeneity lod-score, HFLOD, is computed to quantify the evidence of linkage provided by the data. In this paper, we evaluate this strategy theoretically under different scenarios and compare its performances with those of linkage analysis using affected sib-pair (ASP) data. If cases affected by these Mendelian subentities are not enriched in the sample of cases, the HBD-GWAS strategy has almost no power to detect them, unless they explain an important part of the disease prevalence. The HBD-GWAS strategy outperforms the ASP linkage strategy only in a very limited number of situations where there exists a strong allelic heterogeneity. When several rare recessive variants within the same gene are involved, the ASP design indeed often fails to detect the gene, whereas, by focusing on inbred individuals using the HBD-GWAS strategy, the gene might be detected provided very large samples of cases are available.

Contribution of human growth hormone-releasing hormone receptor (GHRHR) gene sequence variation to isolated severe growth hormone deficiency (ISGHD) and normal adult height

OBJECTIVE

Molecular causes of isolated severe growth hormone deficiency (ISGHD) in several genes have been established. The aim of this study was to analyse the contribution of growth hormone-releasing hormone receptor (GHRHR) gene sequence variation to GH deficiency in a series of prepubertal ISGHD patients and to normal adult height.

DESIGN, SUBJECTS AND MEASUREMENTS

A systematic GHRHR gene sequence analysis was performed in 69 ISGHD patients and 60 normal adult height controls (NAHC). Four GHRHR single-nucleotide polymorphisms (SNPs) were genotyped in 248 additional NAHC. An analysis was performed on individual SNPs and combined genotype associations with diagnosis in ISGHD patients and with height-SDS in NAHC.

RESULTS

Twenty-one SNPs were found. P3, P13, P15 and P20 had not been previously described. Patients and controls shared 12 SNPs (P1, P2, P4-P11, P16 and P21). Significantly different frequencies of the heterozygous genotype and alternate allele were detected in P9 (exon 4, rs4988498) and P12 (intron 6, rs35609199); P9 heterozygous genotype frequencies were similar in patients and the shortest control group (heights between -2 and -1 SDS) and significantly different in controls (heights between -1 and +2 SDS). GHRHR P9 together with 4 GH1 SNP genotypes contributed to 6·2% of height-SDS variation in the entire 308 NAHC.

CONCLUSIONS

This study established the GHRHR gene sequence variation map in ISGHD patients and NAHC. No evidence of GHRHR mutation contribution to ISGHD was found in this population, although P9 and P12 SNP frequencies were significantly different between ISGHD and NAHC. Thus, the gene sequence may contribute to normal adult height, as demonstrated in NAHC.

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.

A novel approach of homozygous haplotype sharing identifies candidate genes in autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association 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.

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

Background

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

Method

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

Results

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

Conclusion

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

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

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

Runs of homozygosity do not influence survival to old age

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

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

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

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.