Runs of homozygosity implicate autozygosity as a schizophrenia risk factor

Discovery of runs-of-homozygosity diplotype clusters and their associations with diseases in UK Biobank

Runs-of-homozygosity (ROH) segments, contiguous homozygous regions in a genome were traditionally linked to families and inbred populations. However, a growing literature suggests that ROHs are ubiquitous in outbred populations. Still, most existing genetic studies of ROH in populations are limited to aggregated ROH content across the genome, which does not offer the resolution for mapping causal loci. This limitation is mainly due to a lack of methods for the efficient identification of shared ROH diplotypes. Here, we present a new method, ROH-DICE (runs-of-homozygous diplotype cluster enumerator), to find large ROH diplotype clusters, sufficiently long ROHs shared by a sufficient number of individuals, in large cohorts. ROH-DICE identified over 1 million ROH diplotypes that span over 100 single nucleotide polymorphisms (SNPs) and are shared by more than 100 UK Biobank participants. Moreover, we found significant associations of clustered ROH diplotypes across the genome with various self-reported diseases, with the strongest associations found between the extended human leukocyte antigen (HLA) region and autoimmune disorders. We found an association between a diplotype covering the homeostatic iron regulator (HFE) gene and hemochromatosis, even though the well-known causal SNP was not directly genotyped or imputed. Using a genome-wide scan, we identified a putative association between carriers of an ROH diplotype in chromosome 4 and an increase in mortality among COVID-19 patients (p-value = 1.82 × 10-11). In summary, our ROH-DICE method, by calling out large ROH diplotypes in a large outbred population, enables further population genetics into the demographic history of large populations. More importantly, our method enables a new genome-wide mapping approach for finding disease-causing loci with multi-marker recessive effects at a population scale.

Genotypic and phenotypic consequences of domestication in dogs

Runs of homozygosity (ROH) are genomic regions that arise when two copies of an identical ancestral haplotype are inherited from parents with a recent common ancestor. In this study, we performed a novel comprehensive analysis to infer genetic diversity among dogs and quantified the association between ROH and non-disease phenotypes. We found distinct patterns of genetic diversity across clades of breed dogs and elevated levels of long ROH, compared to non- domesticated dogs. These high levels of F ROH (inbreeding coefficient) are a consequence of recent inbreeding among domesticated dogs during breed establishment. We identified statistically significant associations between F ROH and height, weight, lifespan, muscled, white head, white chest, furnish, and length of fur. After correcting for population structure, we identified more than 45 genes across the three examined quantitative traits that exceeded the threshold for suggestive significance, indicating significant polygenic inheritance for the complex quantitative phenotypes in dogs.

Influence of autozygosity on common disease risk across the phenotypic spectrum

Autozygosity is associated with rare Mendelian disorders and clinically relevant quantitative traits. We investigated associations between the fraction of the genome in runs of homozygosity (FROH) and common diseases in Genes & Health (n = 23,978 British South Asians), UK Biobank (n = 397,184), and 23andMe. We show that restricting analysis to offspring of first cousins is an effective way of reducing confounding due to social/environmental correlates of FROH. Within this group in G&H+UK Biobank, we found experiment-wide significant associations between FROH and twelve common diseases. We replicated associations with type 2 diabetes (T2D) and post-traumatic stress disorder via within-sibling analysis in 23andMe (median n = 480,282). We estimated that autozygosity due to consanguinity accounts for 5%-18% of T2D cases among British Pakistanis. Our work highlights the possibility of widespread non-additive genetic effects on common diseases and has important implications for global populations with high rates of consanguinity.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

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

BACKGROUND

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

AIM

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

SUBJECTS AND METHODS

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

RESULTS

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

CONCLUSION

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

Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization

Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations.

Exploring the Relationships Between Autozygosity, Educational Attainment, and Cognitive Ability in a Contemporary, Trans-Ancestral American Sample

Previous studies have found significant associations between estimated autozygosity - the proportion of an individual's genome contained in homozygous segments due to distant inbreeding - and multiple traits, including educational attainment (EA) and cognitive ability. In one study, estimated autozygosity showed a stronger association with parental EA than the subject's own EA. This was likely driven by parental EA's association with mobility: more educated parents tended to migrate further from their hometown, and because of the strong correlation between ancestry and geography in the Netherlands, these individuals chose partners farther from their ancestry and therefore more different from them genetically. We examined the associations between estimated autozygosity, cognitive ability, and parental EA in a contemporary sub-sample of adolescents from the Adolescent Brain Cognitive Development Study℠ (ABCD Study®) (analytic N = 6,504). We found a negative association between autozygosity and child cognitive ability consistent with previous studies, while the associations between autozygosity and parental EA were in the expected direction of effect (with greater levels of autozygosity being associated with lower EA) but the effect sizes were significantly weaker than those estimated in previous work. We also found a lower mean level of autozygosity in the ABCD sample compared to previous autozygosity studies, which may reflect overall decreasing levels of autozygosity over generations. Variation in spousal similarities in ancestral background in the ABCD study compared to other studies may explain the pattern of associations between estimated autozygosity, EA, and cognitive ability in the current study.

Schizophrenia: the new etiological synthesis

Schizophrenia has been an evolutionary paradox: it has high heritability, but it is associated with decreased reproductive success. The causal genetic variants underlying schizophrenia are thought to be under weak negative selection. To unravel this paradox, many evolutionary explanations have been suggested for schizophrenia. We critically discuss the constellation of evolutionary hypotheses for schizophrenia, highlighting the lack of empirical support for most existing evolutionary hypotheses-with the exception of the relatively well supported evolutionary mismatch hypothesis. It posits that evolutionarily novel features of contemporary environments, such as chronic stress, low-grade systemic inflammation, and gut dysbiosis, increase susceptibility to schizophrenia. Environmental factors such as microbial infections (e.g., Toxoplasma gondii) can better predict the onset of schizophrenia than polygenic risk scores. However, researchers have not been able to explain why only a small minority of infected people develop schizophrenia. The new etiological synthesis of schizophrenia indicates that an interaction between host genotype, microbe infection, and chronic stress causes schizophrenia, with neuroinflammation and gut dysbiosis mediating this etiological pathway. Instead of just alleviating symptoms with drugs, the parasite x genotype x stress model emphasizes that schizophrenia treatment should focus on detecting and treating possible underlying microbial infection(s), neuroinflammation, gut dysbiosis, and chronic stress.

Homozygous might be hemizygous: CRISPR/Cas9 editing in iPSCs results in detrimental on-target defects that escape standard quality controls

The ability to precisely edit the genome of human induced pluripotent stem cell (iPSC) lines using CRISPR/Cas9 has enabled the development of cellular models that can address genotype to phenotype relationships. While genome editing is becoming an essential tool in iPSC-based disease modeling studies, there is no established quality control workflow for edited cells. Moreover, large on-target deletions and insertions that occur through DNA repair mechanisms have recently been uncovered in CRISPR/Cas9-edited loci. Yet the frequency of these events in human iPSCs remains unclear, as they can be difficult to detect. We examined 27 iPSC clones generated after targeting 9 loci and found that 33% had acquired large, on-target genomic defects, including insertions and loss of heterozygosity. Critically, all defects had escaped standard PCR and Sanger sequencing analysis. We describe a cost-efficient quality control strategy that successfully identified all edited clones with detrimental on-target events and could facilitate the integrity of iPSC-based studies.

Genetic analysis of four consanguineous multiplex families with inflammatory bowel disease

BACKGROUND

Family studies support a genetic predisposition to inflammatory bowel diseases (IBD), but known genetic variants only partially explain the disease heritability. Families with multiple affected individuals potentially harbour rare and high-impact causal variants. Long regions of homozygosity due to recent inbreeding may increase the risk of individuals bearing homozygous loss-of-function variants. This study aimed to identify rare and homozygous genetic variants contributing to IBD.

METHODS

Four families with known consanguinity and multiple cases of IBD were recruited. In a family-specific analysis, we utilised homozygosity mapping complemented by whole-exome sequencing.

RESULTS

We detected a single region of homozygosity shared by Crohn's disease cases from a family of Druze ancestry, spanning 2.6 Mb containing the NOD2 gene. Whole-exome sequencing did not identify any potentially damaging variants within the region, suggesting that non-coding variation may be involved. In addition, affected individuals in the families harboured several rare and potentially damaging homozygous variants in genes with a role in autophagy and innate immunity including LRRK1, WHAMM, DENND3, and C5.

CONCLUSION

This study examined the potential contribution of rare, high-impact homozygous variants in consanguineous families with IBD. While the analysis was not designed to achieve statistical significance, our findings highlight genes or loci that warrant further research. Non-coding variants affecting NOD2 may be of importance in Druze patients with Crohn's disease.

Parental consanguinity among patients with schizophrenia in a rural community of South India: A clinical and genetic investigation

BACKGROUND

Studies from certain regions of the world indicate that consanguineous marriages are a risk factor for the development of schizophrenia in offspring. However the evidence is inconsistent partly due to methodological limitation of which hospital based recruitment contributing to significant bias. The studies from the Indian subcontinent, is scarce, where rates of consanguinity is high.

METHODS

The schizophrenia patients living in a geographically defined rural south Indian community and randomly selected controls dwelling in the same community sharing sociocultural, economic and lifestyle factors were recruited. They were assessed for parental consanguinity using the clinical interviews as well as DNA-based estimates. The latter was conducted by calculating the coefficient of inbreeding 'f'. A participant was considered to have consanguineous parentage if his/her parents shared a common ancestor no more remote than a great-great-grandparent, corresponding to DNA-based estimates of 'f' ≥ 0.0156.

RESULTS

The rates of parental consanguinity assessed by clinical interview were comparable in both groups (Cases: 10.71 %, Controls: 7.25 %; χ² = 0.493, p = 0.4825). However, DNA-based rates of parental consanguinity showed that 'f' was significantly higher among cases than controls (Mann-Whitney U = 11315.5; p = 0.022). Seventy-five cases (62.5 %) and 108 control participants (48.6 %) had 'f' ≥ 0.0156 (χ² = 6.008; p = 0.014). The results were consistent across different quality control measures.

CONCLUSION

Schizophrenia is associated with higher parental consanguinity, suggesting a role for multiple recessive risk alleles in its etiology. Replication in future studies in diverse settings would add further strength to this.

Genomic partitioning of inbreeding depression in humans

Across species, offspring of related individuals often exhibit significant reduction in fitness-related traits, known as inbreeding depression (ID), yet the genetic and molecular basis for ID remains elusive. Here, we develop a method to quantify enrichment of ID within specific genomic annotations and apply it to human data. We analyzed the phenomes and genomes of ∼350,000 unrelated participants of the UK Biobank and found, on average of over 11 traits, significant enrichment of ID within genomic regions with high recombination rates (>21-fold; p < 10^-5), with conserved function across species (>19-fold; p < 10^-4), and within regulatory elements such as DNase I hypersensitive sites (∼5-fold; p = 8.9 × 10^-7). We also quantified enrichment of ID within trait-associated regions and found suggestive evidence that genomic regions contributing to additive genetic variance in the population are enriched for ID signal. We find strong correlations between functional enrichment of SNP-based heritability and that of ID (r = 0.8, standard error: 0.1). These findings provide empirical evidence that ID is most likely due to many partially recessive deleterious alleles in low linkage disequilibrium regions of the genome. Our study suggests that functional characterization of ID may further elucidate the genetic architectures and biological mechanisms underlying complex traits and diseases.

Long tracks of homozygosity predict the severity of alcohol use disorders in an American Indian population

Runs of homozygosity (ROH) arise when an individual inherits two copies of the same haplotype segment. While ROH are ubiquitous across human populations, Native populations-with shared parental ancestry arising from isolation and endogamy-can carry a substantial enrichment for ROH. We have been investigating genetic and environmental risk factors for alcohol use disorders (AUD) in a group of American Indians (AI) who have higher rates of AUD than the general U. S. population. Here we explore whether ROH might be associated with incidence and severity of AUD in this admixed AI population (n = 742) that live on geographically contiguous reservations, using low-coverage whole genome sequences. We have found that the genomic regions in the ROH that were identified in this population had significantly elevated American Indian heritage compared with the rest of the genome. Increased ROH abundance and ROH burden are likely risk factors for AUD severity in this AI population, especially in those diagnosed with severe and moderate AUD. The association between ROH and AUD was mostly driven by ROH of moderate lengths between 1 and 2 Mb. An ROH island on chromosome 1p32.3 and a rare ROH pool on chromosome 3p12.3 were found to be significantly associated with AUD severity. They contain genes involved in lipid metabolism, oxidative stress and inflammatory responses; and OSBPL9 was found to reside on the consensus part of the ROH island. These data demonstrate that ROH are associated with risk for AUD severity in this AI population.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

The Eating Disorders Genetics Initiative (EDGI): study protocol

BACKGROUND

The Eating Disorders Genetics Initiative (EDGI) is an international investigation exploring the role of genes and environment in anorexia nervosa, bulimia nervosa, and binge-eating disorder.

METHODS

A total of 14,500 individuals with eating disorders and 1500 controls will be included from the United States (US), Australia (AU), New Zealand (NZ), and Denmark (DK). In the US, AU, and NZ, participants will complete comprehensive online phenotyping and will submit a saliva sample for genotyping. In DK, individuals with eating disorders will be identified by the National Patient Register, and genotyping will occur using bloodspots archived from birth. A genome-wide association study will be conducted within EDGI and via meta-analysis with other data from the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED).

DISCUSSION

EDGI represents the largest genetic study of eating disorders ever to be conducted and is designed to rapidly advance the study of the genetics of the three major eating disorders (anorexia nervosa, bulimia nervosa, and binge-eating disorder). We will explicate the genetic architecture of eating disorders relative to each other and to other psychiatric and metabolic disorders and traits. Our goal is for EDGI to deliver "actionable" findings that can be transformed into clinically meaningful insights.

TRIAL REGISTRATION

EDGI is a registered clinical trial: clinicaltrials.gov NCT04378101 .

Genome-Wide Detection of Runs of Homozygosity in Laiwu Pigs Revealed by Sequencing Data

Laiwu pigs, distinguished by their high intramuscular fat of 7-9%, is an indigenous pig breed of China, and recent studies also found that Laiwu pigs showed high resistance to Porcine circovirus type 2. However, with the introduction of commercial varieties, the population of Laiwu pigs has declined, and some lineages have disappeared, which could result in inbreeding. Runs of homozygosity (ROH) can be used as a good measure of individual inbreeding status and is also normally used to detect selection signatures so as to map the candidate genes associated with economically important traits. In this study, we used data from Genotyping by Genome Reducing and Sequencing to investigate the number, length, coverage, and distribution patterns of ROH in 93 Chinese Laiwu pigs and identified genomic regions with a high ROH frequency. The average inbreeding coefficient calculated by pedigree was 0.021, whereas that estimated by all detected ROH segments was 0.133. Covering 13.4% of the whole genome, a total of 7,508 ROH segments longer than 1 Mb were detected, whose average length was 3.76 Mb, and short segments (1-5 Mb) dominated. For individuals, the coverage was in the range between 0.56 and 36.86%. For chromosomes, SSC6 had the largest number (n = 688), and the number of ROH in SSC12 was the lowest (n = 215). Thirteen ROH islands were detected in our study, and 86 genes were found within those regions. Some of these genes were correlated with economically important traits, such as meat quality (ECI1, LRP12, NDUFA4L2, GIL1, and LYZ), immunity capacity (IL23A, STAT2, STAT6, TBK1, IFNG, and ITH2), production (DCSTAMP, RDH16, and GDF11), and reproduction (ODF1 and CDK2). A total of six significant Gene Ontology terms and nine significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, most of which were correlated with disease resistance and biosynthesis processes, and one KEGG pathway was related to lipid metabolism. In addition, we aligned all of the ROH islands to the pig quantitative trait loci (QTL) database and finally found eight QTL related to the intramuscular fat trait. These results may help us understand the characteristics of Laiwu pigs and provide insight for future breeding strategies.

The impact of identity by descent on fitness and disease in dogs

Dogs and humans have coexisted together for thousands of years, but it was not until the Victorian Era that humans practiced selective breeding to produce the modern standards we see today. Strong artificial selection during the breed formation period has simplified the genetic architecture of complex traits and caused an enrichment of identity-by-descent (IBD) segments in the dog genome. This study demonstrates the value of IBD segments and utilizes them to infer the recent demography of canids, predict case-control status for complex traits, locate regions of the genome potentially linked to inbreeding depression, and to identify understudied breeds where there is potential to discover new disease-associated variants.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.

Runs of homozygosity are associated with staging of periodontitis in isolated populations

Periodontitis is a common inflammatory disease characterized by a complex etiology, which is the result of a combination of genetic and environmental factors. Genetic variants linked to the periodontitis disease were already investigated, however, little was known regarding the severity of this disease. Recently, long runs of homozygosity (ROH) were associated with several multifactorial diseases. Therefore, in our work, we tried to assess the role of ROH and periodontitis status. We found an association between the excess of homozygosity owing to ROH and staging of periodontitis. More in detail, the total amount of homozygosity owing to ROH is positively associated with an increased severity of periodontitis (P = 0.0001). Regression tree analysis showed the impact of ROH burden in discriminating individuals with mild periodontitis stages I and II and periodontitis stages III and IV (P < 0.001). Furthermore, ROH mapping highlights several regions associated with a severe status of periodontitis (odds ratio > 1). Among them, we found a total of 33 genes. Interestingly, some of these genes were previously associated with granulocyte or platelet measures, both linked to the onset and the progression of periodontal disease. Our results suggest the not only single variants association test could help to risk assessment but even individual genomic features; furthermore, our ROH mapping highlighted the possible role of multiple genes in periodontal development.

Ethnic Identity and Genome Wide Runs of Homozygosity

It is long known that inbreeding increases the detrimental effects of recessive sequence variants in "Runs of Homozygosity" (ROHs). However, although the phenotypic association of ROH has been investigated for a variety of traits, the statistical power of the results often remains limited as a sufficiently high number of cases are available for only a restricted number of traits. In the present study, we aim to analyze the association of runs of homozygosity with the trait "in-group ethnic favoritism". This analysis assumes that if ethnic identity is important for an individual, that individual may tend to marry more frequently within their own group and therefore ROH are expected to increase. We hypothesize that an attitude preferring one's own ethnic group may be associated with a stronger tendency of inbreeding and, as a result, with more and longer ROHs. Accordingly, we investigated the association between the attitude to someone's own ethnicity and ROH, using the Wisconsin Longitudinal data (WLS, total N ~ 9000) as discovery data set and the Brisbane Twin data as replication data set (N ~ 8000). We find that both the number as well as the total length of homozygous segments are significantly positively associated with "in-group ethnic favoritism", independent of the method used for ROH calculation.

Runs of homozygosity and analysis of inbreeding depression

Pedigree information was traditionally used to assess inbreeding. The availability of high-density marker panels provides an alternative to assess inbreeding, particularly in the presence of incomplete and error-prone pedigrees. Assessment of autozygosity across chromosomal segments using runs of homozygosity (ROH) has emerged as a valuable tool to estimate inbreeding due to its general flexibility and ability to quantify the chromosomal contribution to genome-wide inbreeding. Unfortunately, the identification of ROH segments is sensitive to the parameters used during the search process. These parameters are heuristically set, leading to significant variation in the results. The minimum length required to identify an ROH segment has major effects on the estimation of inbreeding and inbreeding depression, yet it is arbitrarily set. To overcome this limitation, a search algorithm to approximate mutation enrichment was developed to determine the minimum length of ROH segments. It consists of finding genome segments with significant effect differences in trait means between animals with high and low burdens of autozygous intervals with a specific length. The minimum length could be determined heuristically as the smallest interval at which a significant signal is detected. The proposed method was tested in an inbred Hereford cattle population genotyped for 30,220 SNPs. Phenotypes recorded for six traits were used for the approximation of mutation loads. The estimated minimum length was around 1 Mb for yearling weight (YW) and average daily gain (ADG) and 4 Mb for birth weight and weaning weight. These trait-specific thresholds estimated using the proposed method could be attributed to a trait-dependent effect of homozygosity. The detection of significant inbreeding effects was well aligned with the estimated thresholds, especially for YW and ADG. Although highly deleterious alleles are expected to be more frequent in recent inbreeding (long ROH), short ROH segments (<5 Mb) could contain a large number of less deleterious mutations with substantial joint effects on some traits (YW and ADG). Our results highlight the importance of accurate estimation of the ROH-based inbreeding and the necessity to consider a trait-specific minimum length threshold for the identification of ROH segments in inbreeding depression analyses. These thresholds could be determined using the proposed method provided the availability of phenotypic information.

Recent Consanguinity and Outbred Autozygosity Are Associated With Increased Risk of Late-Onset Alzheimer's Disease

Prior work in late-onset Alzheimer's disease (LOAD) has resulted in discrepant findings as to whether recent consanguinity and outbred autozygosity are associated with LOAD risk. In the current study, we tested the association between consanguinity and outbred autozygosity with LOAD in the largest such analysis to date, in which 20 LOAD GWAS datasets were retrieved through public databases. Our analyses were restricted to eight distinct ethnic groups: African-Caribbean, Ashkenazi-Jewish European, European-Caribbean, French-Canadian, Finnish European, North-Western European, South-Eastern European, and Yoruba African for a total of 21,492 unrelated subjects (11,196 LOAD and 10,296 controls). Recent consanguinity determination was performed using FSuite v1.0.3, according to subjects' ancestral background. The level of autozygosity in the outbred population was assessed by calculating inbreeding estimates based on the proportion (FROH) and the number (NROH) of runs of homozygosity (ROHs). We analyzed all eight ethnic groups using a fixed-effect meta-analysis, which showed a significant association of recent consanguinity with LOAD (N = 21,481; OR = 1.262, P = 3.6 × 10-4), independently of APOE ∗4 (N = 21,468, OR = 1.237, P = 0.002), and years of education (N = 9,257; OR = 1.274, P = 0.020). Autozygosity in the outbred population was also associated with an increased risk of LOAD, both for F ROH (N = 20,237; OR = 1.204, P = 0.030) and N ROH metrics (N = 20,237; OR = 1.019, P = 0.006), independently of APOE ∗4 [(F ROH, N = 20,225; OR = 1.222, P = 0.029) (N ROH, N = 20,225; OR = 1.019, P = 0.007)]. By leveraging the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data, we determined that LOAD subjects do not show an enrichment of rare, risk-enhancing minor homozygote variants compared to the control population. A two-stage recessive GWAS using ADSP data from 201 consanguineous subjects in the discovery phase followed by validation in 10,469 subjects led to the identification of RPH3AL p.A303V (rs117190076) as a rare minor homozygote variant increasing the risk of LOAD [discovery: Genotype Relative Risk (GRR) = 46, P = 2.16 × 10-6; validation: GRR = 1.9, P = 8.0 × 10-4]. These results confirm that recent consanguinity and autozygosity in the outbred population increase risk for LOAD. Subsequent work, with increased samples sizes of consanguineous subjects, should accelerate the discovery of non-additive genetic effects in LOAD.

Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight

Almond [Prunus dulcis Miller (D.A. Webb)] is the main tree nut species worldwide. Here, genotyping-by-sequencing (GBS) was applied to 149 almond cultivars from the ex situ collections of the Italian Council for Agricultural Research (CREA) and the Spanish National Research Council (CSIC), leading to the detection of 93,119 single-nucleotide polymorphisms (SNPs). The study of population structure outlined four distinct genetic groups and highlighted diversification between the Mediterranean and Californian gene pools. Data on SNP diversity and runs of homozygosity (ROHs) allowed the definition of kinship, inbreeding, and linkage disequilibrium (LD) decay in almond cultivated germplasm. Four-year phenotypic observations, gathered on 98 cultivars of the CREA collection, were used to perform a genome-wide association study (GWAS) and, for the first time in a crop species, homozygosity mapping (HM), resulting in the identification of genomic associations with nut, shell, and seed weight. Both GWAS and HM suggested that loci controlling nut and seed weight are mostly independent. Overall, this study provides insights on the almond cultivation history and delivers information of major interest for almond genetics and breeding. In a broader perspective, our results encourage the use of ROHs in crop science to estimate inbreeding, choose parental combinations minimizing the risk of inbreeding depression, and identify genomic footprints of selection for specific traits.

Discovery of runs-of-homozygosity diplotype clusters and their associations with diseases in UK Biobank

Runs of homozygosity (ROH) segments, contiguous homozygous regions in a genome were traditionally linked to families and inbred populations. However, a growing literature suggests that ROHs are ubiquitous in outbred populations. Still, most existing genetic studies of ROH in populations are limited to aggregated ROH content across the genome, which does not offer the resolution for mapping causal loci. This limitation is mainly due to a lack of methods for efficient identification of shared ROH diplotypes. Here, we present a new method, ROH-DICE, to find large ROH diplotype clusters, sufficiently long ROHs shared by a sufficient number of individuals, in large cohorts. ROH-DICE identified over 1 million ROH diplotypes that span over 100 SNPs and shared by more than 100 UK Biobank participants. Moreover, we found significant associations of clustered ROH diplotypes across the genome with various self-reported diseases, with the strongest associations found between the extended HLA region and autoimmune disorders. We found an association between a diplotype covering the HFE gene and haemochromatosis, even though the well-known causal SNP was not directly genotyped nor imputed. Using genome-wide scan, we identified a putative association between carriers of an ROH diplotype in chromosome 4 and an increase of mortality among COVID-19 patients. In summary, our ROH-DICE method, by calling out large ROH diplotypes in a large outbred population, enables further population genetics into the demographic history of large populations. More importantly, our method enables a new genome-wide mapping approach for finding disease-causing loci with multi-marker recessive effects at population scale.

RG, Bauck S. Genomic mating as sustainable breeding for Chinese indigenous Ningxiang pigs

An important economic reason for the loss of local breeds is that they tend to be less productive, and hence having less market value than commercial breeds. Nevertheless, local breeds often have irreplaceable values, genetically and sociologically. In the breeding programs with local breeds, it is crucial to balance the selection for genetic gain and the maintaining of genetic diversity. These two objectives are often conflicting, and finding the optimal point of the trade-off has been a challenge for breeders. Genomic selection (GS) provides a revolutionary tool for the genetic improvement of farm animals. At the same time, it can increase inbreeding and produce a more rapid depletion of genetic variability of the selected traits in future generations. Optimum-contribution selection (OCS) represents an approach to maximize genetic gain while constraining inbreeding within a targeted range. In the present study, 515 Ningxiang pigs were genotyped with the Illumina Porcine SNP60 array or the GeneSeek Genomic Profiler Porcine 50K array. The Ningxiang pigs were found to be highly inbred at the genomic level. Average locus-wise inbreeding coefficients were 0.41 and 0.37 for the two SNP arrays used, whereas genomic inbreeding coefficients based on runs of homozygosity were 0.24 and 0.25, respectively. Simulated phenotypic data were used to assess the utility of genomic OCS (GOCS) in comparison with GS without inbreeding control. GOCS was conducted under two scenarios, selecting sires only (GOCS_S) or selecting sires and dams (GOCS_SD), while kinships were constrained on selected parents. The genetic gain for average daily body weight gain (ADG) per generation was between 18.99 and 20.55 g with GOCS_S, and between 23.20 and 28.92 with GOCS_SD, and it varied from 25.38 to 48.38 g under GS without controlling inbreeding. While the rate of genetic gain per generation obtained using GS was substantially larger than that obtained by the two scenarios of genomic OCS in the beginning generations of selection, the difference in the genetic gain of ADG between GS and GOCS reduced quickly in latter generations. At generation ten, the difference in the realized rates of genetic gain between GS and GOCS_SD diminished and ended up with even a slightly higher genetic gain with GOCS_SD, due to the rapid loss of genetic variance with GS and fixation of causative genes. The rate of inbreeding was mostly maintained below 5% per generation with genomic OCS, whereas it increased to between 10.5% and 15.3% per generation with GS. Therefore, genomic OCS appears to be a sustainable strategy for the genetic improvement of local breeds such as Ningxiang pigs, but keeping mind that a variety of GOCS methods exist and the optimal forms remain to be exploited further.

Genetic diversity in Iranian indigenous sheep vis-à-vis selected exogenous sheep breeds and wild mouflon

The heterogeneity of climate and different agro-ecological conditions in Iran have resulted in development of 27 indigenous sheep breeds. Wild Asiatic mouflon (Ovis orientalis) is believed to be the ancestor of Iranian sheep. Evaluation of genetic diversity and population structure within and among domestic breeds has important implications for animal breeding programs and genetic resources management. Based on 50K SNP genotype data, we studied the genetic diversity of five indigenous Iranian sheep breeds: Afshari (n = 37), Moghani (n = 34), Qezel (n = 35), Zel (n = 46) and Lori-Bakhtiari (n = 46), and Asiatic mouflon (n = 8) sampled from Iran. Furthermore, genetic diversity and the breed admixture of Iranian sheep were assessed on a larger geographic scale using a reference panel comprising: three indigenous Afghan breeds - Arabi (n = 15), Balouchi (n = 15) and Gadik (n = 15); three indigenous breeds from Turkey and Cyprus - Cyprus Fat Tail (n = 30), Karakas (n = 18) and Norduz (n = 20); and three commercial European breeds - Suffolk (n = 19), Comisana (n = 24) and Engadine Red Sheep (n = 24). The results revealed that the investigated breeds are divided into five genetically distinct clusters according to their geographic origin. Afshari was closest to the local mouflon population and showed signs of mouflon admixture. Qezel was identified as a hybrid sheep breed. Much evidence supported the Afghan breeds being identical. Inbreeding values, which were estimated based on ROHs, were highest for Suffolk (F_ROH  = 0.0544) and lowest for Balouchi (F_ROH  = 0.0078). In conclusion, analysis of selected breeds from neighboring countries along with Asiatic mouflon gave a deeper insight into the evolutionary history and origin of Iranian sheep with important implications for future breed management.

A case of Usher syndrome type IIA caused by a rare USH2A homozygous frameshift variant with maternal uniparental disomy (UPD) in a Chinese family

Usher syndrome encompasses a group of genetically and clinically heterogeneous autosomal recessive disorders with hearing deficiencies and retinitis pigmentosa. The mechanisms underlying the Usher syndrome are highly variable. In the present study, a Chinese family with Usher syndrome was recruited. Whole exome sequencing (WES), Sanger sequencing, homozygosity mapping, short tandem repeat (STR) analysis and segregation analysis were performed. Functional domains of the pathogenic variant for USH2A were analysed. We identified a homozygous frameshift variant c.99_100insT (p.Arg34Serfs*41) in the USH2A gene in the proband that showed discordant segregation in the father. Further homozygosity mapping and STR analysis identified an unusual homozygous variant of proband that originated from maternal uniparental disomy (UPD). The p.Arg34Serfs*41 variant produced a predicted truncated protein that removes all functional domains of USH2A. The variant was not included in the 1000 Human Genomes Project database, ExAC database, HGMD or gnomAD database, but was included in the ClinVar databases as pathogenic. Although USH2A is an autosomal recessive disease, the effects of UPD should be informed in genetic counselling since the recurrence risk of an affected child is greatly reduced when the disease is due to the UPD mechanism. To test potential patients, WES, combined with STR analysis and homozygosity mapping, provides an accurate and useful strategy for genetic diagnosis. In summary, our discoveries can help further the understanding of the molecular pathogenesis of Usher syndrome type IIA to advance the prevention, diagnosis and therapy for this disorder.

Age dependent association of inbreeding with risk for schizophrenia in Egypt

BACKGROUND

Self-reported consanguinity is associated with risk for schizophrenia (SZ) in several inbred populations, but estimates using DNA-based coefficients of inbreeding are unavailable. Further, it is not known whether recessively inherited risk mutations can be identified through homozygosity by descent (HBD) mapping.

METHODS

We studied self-reported and DNA-based estimates of inbreeding among Egyptian patients with SZ (n = 421, DSM IV criteria) and adult controls without psychosis (n = 301), who were evaluated using semi-structured diagnostic interview schedules and genotyped using the Illumina Infinium PsychArray. Following quality control checks, coefficients of inbreeding (F) and regions of homozygosity (ROH) were estimated using PLINK software for HBD analysis. Exome sequencing was conducted in selected cases.

RESULTS

Inbreeding was associated with schizophrenia based on self-reported consanguinity (χ2 = 4.506, 1 df, p = 0.034) and DNA-based estimates for inbreeding (F); the latter with a significant F × age interaction (β = 32.34, p = 0.0047). The association was most notable among patients older than age 40 years. Eleven ROH were over-represented in cases on chromosomes 1, 3, 6, 11, and 14; all but one region is novel for schizophrenia risk. Exome sequencing identified six recessively-acting genes in ROH with loss-of-function variants; one of which causes primary hereditary microcephaly.

CONCLUSIONS

We propose consanguinity as an age-dependent risk factor for SZ in Egypt. HBD mapping is feasible for SZ in adequately powered samples.

AUDACITY: A comprehensive approach for the detection and classification of Runs of Homozygosity in medical and population genomics

Runs of Homozygosity (RoHs) are popular among geneticists as the footprint of demographic processes, evolutionary forces and inbreeding in shaping our genome, and are known to confer risk of Mendelian and complex diseases. Notwithstanding growing interest in their study, there is unmet need for reliable and rapid methods for genomic analyses in large data sets. AUDACITY is a tool integrating novel RoH detection algorithm and autozygosity prediction score for prioritization of mutation-surrounding regions. It processes data in VCF file format, and outperforms existing methods in identifying RoHs of any size. Simulations and analysis of real exomes/genomes show its potential to foster future RoH studies in medical and population genomics.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

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 novel homozygous mutation in GAD1 gene described in a schizophrenic patient impairs activity and dimerization of GAD67 enzyme

Recently, by whole exome sequencing of schizophrenia (SCZ) patients, we identified a subject that was homozygous for a novel missense substitution (c.391 A > G) in the glutamate acid decarboxylase 1 (GAD1) gene. GAD1 encodes for GAD67 enzyme, catalyzing the production of gamma-aminobutyric acid (GABA) from L-glutamic acid. Here, we studied the impact of this mutation on GAD67 activity, dimerization and subcellular localization. Biochemical assay revealed that c.391 A > G reduces GAD67 enzymatic activity by ~30%, probably due to the impaired homodimerization of homozygous mutants as highlighted by proximity ligation assays. The mutational screening of 120 genes of the "GABAergic system" in a cohort of 4,225 SCZ cases and 5,834 controls (dbGaP: phs000473.v1.p2), did not identify other cases that were homozygous for ultra-rare variants in GAD1, but highlighted an increased frequency of cases that were homozygous for rare variants in genes of the GABA system (SCZ: 0.14% vs. Controls: 0.00%; p-value = 0.0055). In conclusion, this study demonstrates the functional impact of c.391 A > G variant and its biological effect makes it a good candidate as risk variant for SCZ. This study also supports an involvement of ultra-rare variants in GABAergic genes in the etiopathogenesis of SCZ.

Recent advances in conservation and population genomics data analysis

New computational methods and next‐generation sequencing (NGS) approaches have enabled the use of thousands or hundreds of thousands of genetic markers to address previously intractable questions. The methods and massive marker sets present both new data analysis challenges and opportunities to visualize, understand, and apply population and conservation genomic data in novel ways. The large scale and complexity of NGS data also increases the expertise and effort required to thoroughly and thoughtfully analyze and interpret data. To aid in this endeavor, a recent workshop entitled “Population Genomic Data Analysis,” also known as “ConGen 2017,” was held at the University of Montana. The ConGen workshop brought 15 instructors together with knowledge in a wide range of topics including NGS data filtering, genome assembly, genomic monitoring of effective population size, migration modeling, detecting adaptive genomic variation, genomewide association analysis, inbreeding depression, and landscape genomics. Here, we summarize the major themes of the workshop and the important take‐home points that were offered to students throughout. We emphasize increasing participation by women in population and conservation genomics as a vital step for the advancement of science. Some important themes that emerged during the workshop included the need for data visualization and its importance in finding problematic data, the effects of data filtering choices on downstream population genomic analyses, the increasing availability of whole‐genome sequencing, and the new challenges it presents. Our goal here is to help motivate and educate a worldwide audience to improve population genomic data analysis and interpretation, and thereby advance the contribution of genomics to molecular ecology, evolutionary biology, and especially to the conservation of biodiversity.

Runs of homozygosity, copy number variation, and risk for depression and suicidal behavior in an Arab Bedouin kindred

OBJECTIVES

Inbreeding increases the probability of homozygosity of deleterious alleles. Inbreeding and runs of homozygosity (ROH) are associated with an increased risk for disease phenotypes, including schizophrenia and other psychiatric disorders. The effects of inbreeding, ROH, homozygous deletions, and other copy number variations (CNVs) on risk for depression and suicide attempt (SA) were quantified in an Arab Bedouin Kindred.

METHODS

We carried out genetic analyses of 439 individuals from an Arab kindred with high rates of depression and suicidal behavior. We obtained complete ascertainment of SAs and first-degree relatives of individuals who have attempted or died by suicide.

RESULTS

We found extensive regions of ROH. On average, 5% of the genome is covered by ROH for these individuals, two-fold higher than ROH rates for individuals from populations of European ancestry. Inbreeding and total length of ROH were not associated with risk for depression or attempt. For CNVs, an increased number of duplications more than 500 kb was associated with an increased risk for attempt (odds ratio: 2.9; P=0.01; 95% confidence interval: 1.3-6.6). Although not significant after correction for multiple testing, the risk for SA appears to increase with copy number for a CNV on chromosome 9p24.1. This possibility is intriguing because the CNV covers GLDC, which encodes glycine dehydrogenase that binds to glycine, a co-agonist at N-methyl-D-aspartate glutamate receptors, and is involved in glutamatergic neurotransmission.

CONCLUSION

Our findings add to the growing evidence of genetic risk factors that act pleiotropically to increase the risk for several neuropsychiatric disorders, including depression and SA, irrespective of ancestry.

Relationships between estimated autozygosity and complex traits in the UK Biobank

Inbreeding increases the risk of certain Mendelian disorders in humans but may also reduce fitness through its effects on complex traits and diseases. Such inbreeding depression is thought to occur due to increased homozygosity at causal variants that are recessive with respect to fitness. Until recently it has been difficult to amass large enough sample sizes to investigate the effects of inbreeding depression on complex traits using genome-wide single nucleotide polymorphism (SNP) data in population-based samples. Further, it is difficult to infer causation in analyses that relate degree of inbreeding to complex traits because confounding variables (e.g., education) may influence both the likelihood for parents to outbreed and offspring trait values. The present study used runs of homozygosity in genome-wide SNP data in up to 400,000 individuals in the UK Biobank to estimate the proportion of the autosome that exists in autozygous tracts—stretches of the genome which are identical due to a shared common ancestor. After multiple testing corrections and controlling for possible sociodemographic confounders, we found significant relationships in the predicted direction between estimated autozygosity and three of the 26 traits we investigated: age at first sexual intercourse, fluid intelligence, and forced expiratory volume in 1 second. Our findings corroborate those of several published studies. These results may imply that these traits have been associated with Darwinian fitness over evolutionary time. However, some of the autozygosity-trait relationships were attenuated after controlling for background sociodemographic characteristics, suggesting that alternative explanations for these associations have not been eliminated. Care needs to be taken in the design and interpretation of ROH studies in order to glean reliable information about the genetic architecture and evolutionary history of complex traits.

Inbreeding is well known to increase the risk of rare, monogenic diseases, and there has been some evidence that it also affects complex traits, such as cognition and educational attainment. However, difficulties can arise when inferring causation in these types of analyses because of the potential for confounding variables (e.g., socioeconomic status) to bias the observed relationships between distant inbreeding and complex traits. In this investigation, we used single-nucleotide polymorphism data in a very large (N > 400,000) sample of seemingly outbred individuals to quantify the degree to which distant inbreeding is associated with 26 complex traits. We found robust evidence that distant inbreeding is inversely associated with fluid intelligence and a measure of lung function, and is positively associated with age at first sex, while other trait associations with inbreeding were attenuated after controlling for background sociodemographic characteristics. Our findings are consistent with evolutionary predictions that fluid intelligence, lung function, and age at first sex have been under selection pressures over time; however, they also suggest that confounding variables must be accounted for in order to reliably interpret results from these types of analyses.

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.

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

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

Artificial selection increased body weight but induced increase of runs of homozygosity in Hanwoo cattle

Artificial selection has been demonstrated to have a rapid and significant effect on the phenotype and genome of an organism. However, most previous studies on artificial selection have focused solely on genomic sequences modified by artificial selection or genomic sequences associated with a specific trait. In this study, we generated whole genome sequencing data of 126 cattle under artificial selection, and 24,973,862 single nucleotide variants to investigate the relationship among artificial selection, genomic sequences and trait. Using runs of homozygosity detected by the variants, we showed increase of inbreeding for decades, and at the same time demonstrated a little influence of recent inbreeding on body weight. Also, we could identify ~0.2 Mb runs of homozygosity segment which may be created by recent artificial selection. This approach may aid in development of genetic markers directly influenced by artificial selection, and provide insight into the process of artificial selection.

Assessing runs of Homozygosity: a comparison of SNP Array and whole genome sequence low coverage data

BACKGROUND

Runs of Homozygosity (ROH) are genomic regions where identical haplotypes are inherited from each parent. Since their first detection due to technological advances in the late 1990s, ROHs have been shedding light on human population history and deciphering the genetic basis of monogenic and complex traits and diseases. ROH studies have predominantly exploited SNP array data, but are gradually moving to whole genome sequence (WGS) data as it becomes available. WGS data, covering more genetic variability, can add value to ROH studies, but require additional considerations during analysis.

RESULTS

Using SNP array and low coverage WGS data from 1885 individuals from 20 world populations, our aims were to compare ROH from the two datasets and to establish software conditions to get comparable results, thus providing guidelines for combining disparate datasets in joint ROH analyses. By allowing heterozygous SNPs per window, using the PLINK homozygosity function and non-parametric analysis, we were able to obtain non-significant differences in number ROH, mean ROH size and total sum of ROH between data sets using the different technologies for almost all populations.

CONCLUSIONS

By allowing 3 heterozygous SNPs per ROH when dealing with WGS low coverage data, it is possible to establish meaningful comparisons between data using SNP array and WGS low coverage technologies.

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.

A novel relationship for schizophrenia, bipolar and major depressive disorder Part 3: Evidence from chromosome 3 high density association screen

Familial clustering of schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD) was systematically reported (Aukes et al, Genet Med 2012, 14, 338-341) and convergent evidence from genetics, symptomatology, and psychopharmacology imply that there are intrinsic connections between these three major psychiatric disorders, for example, any two or even three of these disorders could co-exist in some families. A total of 60, 838 single-nucleotide polymorphisms (SNPs) on chromosome 3 were genotyped by Affymetrix Genome-Wide Human SNP array 6.0 on 119 SCZ, 253 BPD (type-I), 177 MDD patients and 1,000 controls. The population of Shandong province was formed in 14 century and believed that it belongs to homogenous population. Associated SNPs were systematically revealed and outstanding susceptibility genes (CADPS, GRM7,KALRN, LSAMP, NLGN1, PRICKLE2, ROBO2) were identified. Unexpectedly, flanking genes for the associated SNPs distinctive for BPD and/or MDD were replicated in an enlarged cohort of 986 SCZ patients. The evidence from this chromosome 3 analysis supports the notion that both of bipolar and MDD might be subtypes of schizophrenia rather than independent disease entity. Also, a similar finding was detected on chromosome 5, 6, 7, and 8 (Chen et al. Am J Transl Res 2017;9 (5):2473-2491; Curr Mol Med 2016;16(9):840-854; Behav Brain Res 2015;293:241-251; Mol Neurobiol 2016. doi: 10.1007/s12035-016-0102-1). Furthermore, PRICKLE2 play an important role in the pathogenesis of three major psychoses in this population.

Consanguinity Rates Predict Long Runs of Homozygosity in Jewish Populations

OBJECTIVES

Recent studies have highlighted the potential of analyses of genomic sharing to produce insight into the demographic processes affecting human populations. We study runs of homozygosity (ROH) in 18 Jewish populations, examining these groups in relation to 123 non-Jewish populations sampled worldwide.

METHODS

By sorting ROH into 3 length classes (short, intermediate, and long), we evaluate the impact of demographic processes on genomic patterns in Jewish populations.

RESULTS

We find that the portion of the genome appearing in long ROH - the length class most directly related to recent consanguinity - closely accords with data gathered from interviews during the 1950s on frequencies of consanguineous unions in various Jewish groups.

CONCLUSION

The high correlation between 1950s consanguinity levels and coverage by long ROH explains differences across populations in ROH patterns. The dissection of ROH into length classes and the comparison to consanguinity data assist in understanding a number of additional phenomena, including similarities of Jewish populations to Middle Eastern, European, and Central and South Asian non-Jewish populations in short ROH patterns, relative lengths of identity-by-descent tracts in different Jewish groups, and the "population isolate" status of the Ashkenazi Jews.

Detection and quantification of inbreeding depression for complex traits from SNP data

Quantifying the effects of inbreeding is critical to characterizing the genetic architecture of complex traits. This study highlights through theory and simulations the strengths and shortcomings of three SNP-based inbreeding measures commonly used to estimate inbreeding depression (ID). We demonstrate that heterogeneity in linkage disequilibrium (LD) between causal variants and SNPs biases ID estimates, and we develop an approach to correct this bias using LD and minor allele frequency stratified inference (LDMS). We quantified ID in 25 traits measured in [Formula: see text] participants of the UK Biobank, using LDMS, and confirmed previously published ID for 4 traits. We find unique evidence of ID for handgrip strength, waist/hip ratio, and visual and auditory acuity (ID between -2.3 and -5.2 phenotypic SDs for complete inbreeding; [Formula: see text]). Our results illustrate that a careful choice of the measure of inbreeding combined with LDMS stratification improves both detection and quantification of ID using SNP data.

Genome-wide mapping and estimation of inbreeding depression of semen quality traits in a cattle population

Inbreeding depression is known to affect quantitative traits such as male fertility and sperm quality, but the genetic basis for these associations is poorly understood. Most studies have been limited to examining how pedigree- or marker-derived genome-wide autozygosity is associated with quantitative phenotypes. In this study, we analyzed possible associations of genetic features of inbreeding depression with percentage of live spermatozoa and total number of spermatozoa in 19,720 ejaculates obtained from 554 Austrian Fleckvieh bulls during routine artificial insemination programs. Genome-wide inbreeding depression was estimated and genomic regions contributing to inbreeding depression were mapped. Inbreeding depression did affect total number of spermatozoa, and such depression was predicted by pedigree-based inbreeding levels and genome-wide inbreeding levels based on runs of homozygosity (ROH). Genome-wide inbreeding depression did not seem to affect percentage of live spermatozoa. A model incorporating genetic effects of the bull, environmental factors, and additive genetic and ROH status effects of individual single-nucleotide polymorphisms revealed genomic regions significantly associated with ROH status for total number of spermatozoa (4 regions) or percentage of live spermatozoa (5 regions). All but one region contains genes related to spermatogenesis and sperm morphology. These genomic regions contain genes affecting sperm morphogenesis and efficacy. The results highlight that next-generation sequencing may help explain some of the genetic factors contributing to inbreeding depression of sperm quality traits in Fleckvieh bulls.

Structural and genomic variation in preterm birth

BACKGROUND

Runs of homozygosity (ROH) are consecutive homozygous genotypes, which may result from population inbreeding or consanguineous marriages. ROH enhance the expression of recessive traits.

METHODS

We mapped ROH in a case control study of women delivering at term compared with women delivering at or before 34 wk gestation. Gene sets known to be important in risk of preterm birth were examined for their overlap with identified ROH segments.

RESULTS

While we found no evidence of increased burden of ROH or copy number variations in mothers delivering at or before 34 wk compared with term, we identified 424 genome-wide 50 kb segments with significant difference in abundance of overlapping ROH segments in cases vs. controls, P < 0.05. These regions overlap 199 known genes. We found preterm birth associated genes (CXCR4, MYLK, PAK1) and genes shown to have an evolutionary link to preterm (CXCR4, PPP3CB, C6orf57, DUSP13, and SLC25A45) with significant differences in abundance of overlapping ROH blocks in cases vs. controls, P < 0.001.

CONCLUSION

We conclude, while we found no significant burden of ROH, we did identify genomic regions with significantly greater abundance of ROH blocks in women delivering preterm that overlapped genes known to be involved in preterm birth.

Genomics advances the study of inbreeding depression in the wild

Inbreeding depression (reduced fitness of individuals with related parents) has long been a major focus of ecology, evolution, and conservation biology. Despite decades of research, we still have a limited understanding of the strength, underlying genetic mechanisms, and demographic consequences of inbreeding depression in the wild. Studying inbreeding depression in natural populations has been hampered by the inability to precisely measure individual inbreeding. Fortunately, the rapidly increasing availability of high-throughput sequencing data means it is now feasible to measure the inbreeding of any individual with high precision. Here, we review how genomic data are advancing our understanding of inbreeding depression in the wild. Recent results show that individual inbreeding and inbreeding depression can be measured more precisely with genomic data than via traditional pedigree analysis. Additionally, the availability of genomic data has made it possible to pinpoint loci with large effects contributing to inbreeding depression in wild populations, although this will continue to be a challenging task in many study systems due to low statistical power. Now that reliably measuring individual inbreeding is no longer a limitation, a major focus of future studies should be to more accurately quantify effects of inbreeding depression on population growth and viability.