Rare and common variant discovery by whole-genome sequencing of 101 Thoroughbred racehorses.
Ontology highlight
ABSTRACT:
SUBMITTER: Tozaki T
PROVIDER: S-EPMC8346562 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Thoroughbreds are some of the most famous racehorses worldwide and are currently animals of high economic value. To understand genomic variability in Thoroughbreds, we identified genome-wide insertions and deletions (INDELs) and obtained their allele frequencies in this study. INDELs were obtained from whole-genome sequencing data of 101 Thoroughbred racehorses by mapping sequence reads to the horse reference genome. By integrating individual data, 1,453,349 and 113,047 INDELs were identified in the autosomal (1-31) and X chromosomes, respectively, while 18 INDELs were identified on the mitochondrial genome, totaling 1,566,414 INDELs. Of those, 779,457 loci (49.8%) were novel INDELs, while 786,957 loci (50.2%) were already registered in Ensembl. The sizes of diallelic INDELs ranged from -286 to +476, and the majority, 717,736 (52.14%) and 220,672 (16.03%), were 1-bp and 2-bp variants, respectively. Numerous INDELs were found to have lower frequencies of alternative (Alt) alleles. Many rare variants with low Alt allele frequencies (<0.5%) were also detected. In addition, 5955 loci were genotyped as having a minor allele frequency of 0.5 and being heterogeneous genotypes in all the horses. While short-read sequencing and its mapping to reference genome is a simple way of detecting variants, fake variants may be detected. Therefore, our data help to identify true variants in Thoroughbred horses. The INDEL database we constructed will provide useful information for genetic studies and industrial applications in Thoroughbred horses, including a gene-editing test for gene-doping control and a parentage test using INDELs for horse registration and identification.
| S-EPMC10048024 | biostudies-literature
Project description:The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of <1.0%. This detailed analysis detected signatures for purifying selection on regulatory elements as well as coding regions. We also catalogue structural variants, including 3.4 million insertions and deletions, and 25,923 genic copy-number variants. The 1KJPN was effective for imputing genotypes of the Japanese population genome wide. These data demonstrate the value of high-coverage sequencing for constructing population-specific variant panels, which covers 99.0% SNVs of minor allele frequency ≥0.1%, and its value for identifying causal rare variants of complex human disease phenotypes in genetic association studies.
| S-EPMC4560751 | biostudies-literature
Project description:The domestication process of plants and animals typically involves intense inbreeding and directional selection for various traits. Here, we genotyped 370 Japanese Thoroughbred horses using the recently developed 670k SNP array and performed various genome-wide analysis also using genotype data of other horse breeds. We identified a number of regions showing interesting patterns of polymorphisms. For instance, the region containing the MC1R locus associated with chestnut coat color may have been targeted by selection for a different mutation much earlier on than the recent selection for chestnut color. We also identified regions that show signatures of selection specific to Thoroughbreds. In addition, we found that intense inbreeding early in the history of the Thoroughbred breed and also before the formation of the breed has a significant impact on the genomic architecture of modern Thoroughbreds. Our study demonstrates that the horse 670k array can be utilized to gain important insight into the domestication process of horses and to understand the genetic basis of the phenotypic diversity in horses.
| S-EPMC6655603 | biostudies-literature
Project description:Essential tremor (ET) is one of the most common movement disorders. The etiology of ET remains largely unexplained. Whole genome sequencing (WGS) is likely to be of value in understanding a large proportion of ET with Mendelian and complex disease inheritance patterns. In ET families with Mendelian inheritance patterns, WGS may lead to gene identification where WES analysis failed to identify the causative single nucleotide variant (SNV) or indel due to incomplete coverage of the entire coding region of the genome, in addition to accurate detection of larger structural variants (SVs) and copy number variants (CNVs). Alternatively, in ET families with complex disease inheritance patterns with gene x gene and gene x environment interactions enrichment of functional rare coding and non-coding variants may explain the heritability of ET. We performed WGS in eight ET families (n = 40 individuals) enrolled in the Family Study of Essential Tremor. The analysis included filtering WGS data based on allele frequency in population databases, rare SNV and indel classification and association testing using the Mixed-Model Kernel Based Adaptive Cluster (MM-KBAC) test. A separate analysis of rare SV and CNVs segregating within ET families was also performed. Prioritization of candidate genes identified within families was performed using phenolyzer. WGS analysis identified candidate genes for ET in 5/8 (62.5%) of the families analyzed. WES analysis in a subset of these families in our previously published study failed to identify candidate genes. In one family, we identified a deleterious and damaging variant (c.1367G>A, p.(Arg456Gln)) in the candidate gene, CACNA1G, which encodes the pore forming subunit of T-type Ca(2+) channels, CaV3.1, and is expressed in various motor pathways and has been previously implicated in neuronal autorhythmicity and ET. Other candidate genes identified include SLIT3 which encodes an axon guidance molecule and in three families, phenolyzer prioritized genes that are associated with hereditary neuropathies (family A, KARS, family B, KIF5A and family F, NTRK1). Functional studies of CACNA1G and SLIT3 suggest a role for these genes in ET disease pathogenesis.
| S-EPMC6690583 | biostudies-literature
Project description:Metabolomics-genome-wide association studies (mGWAS) have uncovered many metabolic quantitative trait loci (mQTLs) influencing human metabolic individuality, though predominantly in European cohorts. By combining whole-exome sequencing with a high-resolution metabolomics profiling for a highly consanguineous Middle Eastern population, we discover 21 common variant and 12 functional rare variant mQTLs, of which 45% are novel altogether. We fine-map 10 common variant mQTLs to new metabolite ratio associations, and 11 common variant mQTLs to putative protein-altering variants. This is the first work to report common and rare variant mQTLs linked to diseases and/or pharmacological targets in a consanguineous Arab cohort, with wide implications for precision medicine in the Middle East.
| S-EPMC5780481 | biostudies-literature
Project description:For many complex traits/diseases, it is believed that rare variants account for some of the missing heritability that cannot be explained by common variants. Sequencing a large number of samples through DNA pooling is a cost-effective strategy to discover rare variants and to investigate their associations with phenotypes. Overlapping pool designs provide further benefit because such approaches can potentially identify variant carriers, which is important for downstream applications of association analysis of rare variants. However, existing algorithms for analysing sequence data from overlapping pools are limited.We propose a complete data analysis framework for overlapping pool designs, with novelties in all three major steps: variant pool and variant locus identification, variant allele frequency estimation and variant sample decoding. The framework can be used in combination with any design matrix. We have investigated its performance based on two different overlapping designs and have compared it with three state-of-the-art methods, by simulating targeted sequencing and by pooling real sequence data. Results on both datasets show that our algorithm has made significant improvements over existing ones. In conclusion, successful discovery of rare variants and identification of variant carriers using overlapping pool strategies critically depend on many steps, from generation of design matrixes to decoding algorithms. The proposed framework in combination with the design matrixes generated based on the Chinese remainder theorem achieves best overall results.Source code of the program, termed VIP for Variant Identification by Pooling, is available at http://cbc.case.edu/VIP.
| S-EPMC3530907 | biostudies-other
Project description: Not available
| S-EPMC7446391 | biostudies-literature
Project description:Whole-genome sequencing (WGS) studies are being widely conducted in order to identify rare variants associated with human diseases and disease-related traits. Classical single-marker association analyses for rare variants have limited power, and variant-set-based analyses are commonly used by researchers for analyzing rare variants. However, existing variant-set-based approaches need to pre-specify genetic regions for analysis; hence, they are not directly applicable to WGS data because of the large number of intergenic and intron regions that consist of a massive number of non-coding variants. The commonly used sliding-window method requires the pre-specification of fixed window sizes, which are often unknown as a priori, are difficult to specify in practice, and are subject to limitations given that the sizes of genetic-association regions are likely to vary across the genome and phenotypes. We propose a computationally efficient and dynamic scan-statistic method (Scan the Genome [SCANG]) for analyzing WGS data; this method flexibly detects the sizes and the locations of rare-variant association regions without the need to specify a prior, fixed window size. The proposed method controls for the genome-wise type I error rate and accounts for the linkage disequilibrium among genetic variants. It allows the detected sizes of rare-variant association regions to vary across the genome. Through extensive simulated studies that consider a wide variety of scenarios, we show that SCANG substantially outperforms several alternative methods for detecting rare-variant-associations while controlling for the genome-wise type I error rates. We illustrate SCANG by analyzing the WGS lipids data from the Atherosclerosis Risk in Communities (ARIC) study.
| S-EPMC6507043 | biostudies-literature
Project description:Training has a great impact on the physiology of an athlete and, like all stressful stimuli, can trigger an innate immune response and inflammation, which is part of a wider coping strategy of the host to restore homeostasis. The Thoroughbred racehorse is a valid animal model to investigate these changes thanks to its homogeneous training and highly selected genetic background. The aim of this study was to investigate modifications of the innate immune response and inflammation in young untrained Thoroughbred racehorses during the first training season through haematological and molecular investigations. Twenty-nine Thoroughbred racehorses were followed during their incremental 3-month sprint exercise schedule. Blood collection was performed at time 0 (T0; before starting the intense training period), 30 days after T0 (T30), and 90 days after T0 (T90). Haematological parameters (red and white blood cells, haemoglobin, and platelets) were evaluated and haematocrit (HCT), mean corpuscular haemoglobin concentration (MCHC), and red cells width distribution + standard deviation (RDW-SD) were calculated. Moreover, via RT-qPCR, we investigated the expression of, Interleukin 1? (IL-1?), Interleukin 4 (IL-4) Interleukin 6 (IL-6), Interleukin 2 (IL-2), Interleukin 3 (IL-3), Interleukin 5 (IL-5) Interleukin 8 (IL-8), Trasformig Growth Factor ? and ? (TGF-?), Tumor necrosis factor ? (TNF-?), and Interferon ? (IFN-?)genes. Main corpuscular volume (MCV) showed a significant (p = 0.008) increase at T90. Main corpuscular haemoglobin (MCH) and haemoglobin concentration (MCHC) values were significantly augmented at both T30 (p < 0.001) and T90 (p < 0.001). Basophils were significant increased at T30 (p = 0.02) and eosinophils were significantly increased at T90 (p = 0.03). Significant differences in gene expression were found for all the genes under study, with the exception of IFN-? and TNF-?. In particular, IL-2 (T30, p = 0.011; T90, p = 0.015), IL-4 (T30, p = 0.009; T90, p < 0.001), and IL-8 (T30, p < 0.001; T90, p < 0.001) genes were significantly upregulated at both T30 and T90 with respect to T0, TGF-? was intensely downregulated at T30 (p < 0.001), IL-5 gene expression was significantly decreased at T90 (p = 0.001), while IL-1? (p = 0.005) and IL-3 (p = 0.001) expression was strongly augmented at the same time. This study highlighted long-term adjustments of O2 transport capability that can be reasonably traced back to exercise adaptation. Moreover, the observed changes of granulocyte numbers and functions and inflammatory cytokine gene expression confirm a major role of the innate immune system in the response to the complex of stressful stimuli experienced during the training period.
| S-EPMC7600081 | biostudies-literature
Project description:Blepharospasm (BSP) is a type of focal dystonia characterized by involuntary orbicularis oculi spasms that are usually bilateral, synchronous, and symmetrical. Despite strong evidence for genetic contributions to BSP, progress in the field has been constrained by small cohorts, incomplete penetrance, and late age of onset. Although several genetic etiologies for dystonia have been identified through whole-exome sequencing (WES), none of these are characteristically associated with BSP as a singular or predominant manifestation. We performed WES on 31 subjects from 21 independent pedigrees with BSP. The strongest candidate sequence variants derived from in silico analyses were confirmed with bidirectional Sanger sequencing and subjected to cosegregation analysis. Cosegregating deleterious variants (GRCH37/hg19) in CACNA1A (NM_001127222.1: c.7261_7262delinsGT, p.Pro2421Val), REEP4 (NM_025232.3: c.109C>T, p.Arg37Trp), TOR2A (NM_130459.3: c.568C>T, p.Arg190Cys), and ATP2A3 (NM_005173.3: c.1966C>T, p.Arg656Cys) were identified in four independent multigenerational pedigrees. Deleterious variants in HS1BP3 (NM_022460.3: c.94C>A, p.Gly32Cys) and GNA14 (NM_004297.3: c.989_990del, p.Thr330ArgfsTer67) were identified in a father and son with segmental cranio-cervical dystonia first manifest as BSP. Deleterious variants in DNAH17, TRPV4, CAPN11, VPS13C, UNC13B, SPTBN4, MYOD1, and MRPL15 were found in two or more independent pedigrees. To our knowledge, none of these genes have previously been associated with isolated BSP, although other CACNA1A mutations have been associated with both positive and negative motor disorders including ataxia, episodic ataxia, hemiplegic migraine, and dystonia. Our WES datasets provide a platform for future studies of BSP genetics which will demand careful consideration of incomplete penetrance, pleiotropy, population stratification, and oligogenic inheritance patterns.
| S-EPMC6081235 | biostudies-literature