Dataset Information

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances.

ABSTRACT: Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our aim was to identify ROH islands influenced by marker coverage gaps or CNV, using Illumina BovineHD BeadChip (777 K) single nucleotide polymorphism (SNP) data for Austrian Brown Swiss, Tyrol Grey and Pinzgauer cattle.ROH were detected using clustering, and ROH islands were determined from population inbreeding levels for each marker. CNV were detected using a multivariate copy number analysis method and a hidden Markov model. SNP coverage gaps were defined as genomic regions with intermarker distances on average longer than 9.24 kb. ROH islands that overlapped CNV regions (CNVR) or SNP coverage gaps were considered as potential artefacts. Permutation tests were used to determine if overlaps between CNVR with copy losses and ROH islands were due to chance. Diversity of the haplotypes in the ROH islands was assessed by haplotype analyses.In Brown Swiss, Tyrol Grey and Pinzgauer, we identified 13, 22, and 24 ROH islands covering 26.6, 389.0 and 35.8 Mb, respectively, and we detected 30, 50 and 71 CNVR derived from CNV by using both algorithms, respectively. Overlaps between ROH islands, CNVR or coverage gaps occurred for 7, 14 and 16 ROH islands, respectively. About 37, 44 and 52% of the ROH islands coverage in Brown Swiss, Tyrol Grey and Pinzgauer, respectively, were affected by copy loss. Intersections between ROH islands and CNVR were small, but significantly larger compared to ROH islands at random locations across the genome, implying an association between ROH islands and CNVR. Haplotype diversity for reliable ROH islands was lower than for ROH islands that intersected with copy loss CNVR.Our findings show that a significant proportion of the ROH islands in the bovine genome are artefacts due to CNV or SNP coverage gaps.

SUBMITTER: Nandolo W

PROVIDER: S-EPMC6106898 | biostudies-other | 2018 Aug

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances.

Nandolo Wilson W Utsunomiya Yuri T YT Mészáros Gábor G Wurzinger Maria M Khayadzadeh Negar N Torrecilha Rafaela B P RBP Mulindwa Henry A HA Gondwe Timothy N TN Waldmann Patrik P Ferenčaković Maja M Garcia José F JF Rosen Benjamin D BD Bickhart Derek D van Tassell Curt P CP Curik Ino I Sölkner Johann J

Genetics, selection, evolution : GSE 20180822 1

<h4>Background</h4>Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our ...[more]

PMID: 30134820

Similar Datasets

Project description:BACKGROUND:Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (FROH), genomic relationship matrix approach (FGRM) and based on the observed versus expected number of homozygous genotypes (FHOM), and from pedigree-based coefficient (FPED). RESULTS:ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH1-2 Mb) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of FROH, FGRM, FHOM, and FPED approaches. FPED estimates ranged from 0.00 to 0.327 and FROH from 0.001 to 0.201. Low to moderate correlations were observed between FPED-FROH and FGRM-FROH, with values ranging from -0.11 to 0.51. Low to high correlations were observed between FROH-FHOM and moderate between FPED-FHOM and FGRM-FHOM. Correlations between FROH from different lengths and FPED gradually increased with ROH length. CONCLUSIONS:Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low FPED-FROH correlations for small segments indicate that FPED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Project description:Simple Summary Wenchang chickens are the only chicken breed listed in the ‘animal genetic resources in China (poultry)’ in Hainan Province and are famous for their excellent meat quality. Protection of this genetic resource may ensure poultry husbandry’s sustainable and successful development. For more effective conservation, development, and utilization of this genetic resource, we investigated the diversity, degree of inbreeding, and runs of homozygosity (ROH) patterns for Wenchang chickens using whole-genome sequencing data. Our analysis showed that the genetic diversity of Wenchang chickens was relatively high. Selection signal analysis of Wenchang chickens based on ROH found some candidate genes that were putatively associated with meat quality traits and stress resistance traits, such as disease resistance and heat tolerance. Abstract Wenchang chickens, a native breed in the Hainan province of China, are famous for their meat quality and adaptability to tropical conditions. For effective management and conservation, in the present study, we systematically investigated the characteristics of genetic variations and runs of homozygosity (ROH) along the genome using re-sequenced whole-genome sequencing data from 235 Wenchang chickens. A total of 16,511,769 single nucleotide polymorphisms (SNPs) and 53,506 ROH segments were identified in all individuals, and the ROH of Wenchang chicken were mainly composed of short segments (0–1 megabases (Mb)). On average, 5.664% of the genome was located in ROH segments across the Wenchang chicken samples. According to several parameters, the genetic diversity of the Wenchang chicken was relatively high. The average inbreeding coefficient of Wenchang chickens based on FHOM, FGRM, and FROH was 0.060 ± 0.014, 0.561 ± 0.020, and 0.0566 ± 0.01, respectively. A total of 19 ROH islands containing 393 genes were detected on 9 different autosomes. Some of these genes were putatively associated with growth performance (AMY1a), stress resistance (THEMIS2, PIK3C2B), meat traits (MBTPS1, DLK1, and EPS8L2), and fat deposition (LANCL2, PPARγ). These findings provide a better understanding of the degree of inbreeding in Wenchang chickens and the hereditary basis of the characteristics shaped under selection. These results are valuable for the future breeding, conservation, and utilization of Wenchang and other chicken breeds.

Project description:BackgroundExome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES.MethodsIn a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array.ResultsThe QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities.ConclusionsOur clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Dataset Information

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances.

Publications

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets