شناسایی نشانه های انتخاب مرتبط با آترزی دستگاه گوارش در گوساله های هلشتاین

ضایعات انسدادی دستگاه گوارش از مهمترین مشکلات مادرزادی هستند که نتیجه آن، مرگ و میر گوساله ها در مدت کوتاهی بعد از تولد است. ژژنوم شایع ترین مکان بعد از مری برای انسداد کامل یا آترزی است که فقدان مادرزادی یا انسداد کامل قسمتی از حفره روده است. شناخت زودرس آترزی برای پیشگیری از عوارض بعدی آن ضروری است. هدف پژوهش حاضر، شناسایی نشانه های انتخاب با استفاده از آماره نااریب تتا مرتبط با آترزی دستگاه گوارش در گوساله های شیری هلشتاین است. به این منظور، مجموع 466 حیوان برای 777962 جایگاه نشانگری چندشکلی تک نوکلئوتیدی (snp) با استفاده از تراشه های illumina 777k bovinehd  تعیین ژنوتیپ شدند. پس از کنترل کیفیت داده های اولیه در نهایت، 704242 نشانگر snp و 466 راس دام (375 راس شاهد و 91 راس موردی) وارد تجزیه های بعدی شدند. با در نظر گرفتن صدک 9/99 کل ارزش های تتا، هشت ناحیه ژنومی روی کروموزوم‌های 7، 12، 13، 21، 22، 23 (دو نقطه) و 29 شناسایی شدند. بررسی ژن های گزارش شده در این مناطق نشان داد که در داخل یا مجاورت این نواحی، ژن های csf2، siah3، tmem14a و skiv2l قرار داشتند، که با رشد و توسعه جنین، طول روده کوچک، مرگ سلولی و انواع سرطان مرتبط هستند. همچنین بررسی بیوانفورماتیکی این مناطق نشان داد ژن های موجود در این مناطق با جایگاه های ژنی کنترل کننده صفات کمّی مرتبط با حساسیت به بیماری و مسیرهای هستی شناسی مرتبط با تنظیم تفرق و تمایز سلولی، فعالیت پروتئین های کینازی و فعالیت سوخت و ساز سلولی هستند. نتایج این تحقیق می تواند منبع اطلاعاتی ارزشمندی در زمینه شناسایی مناطق ژنومی و در نتیجه ژن های مرتبط با آترزی در گاوهای هلشتاین فراهم آورد.

identification of selective signatures associated with gastrointestinal atresia in holstein calves

introduction: obstructive gastrointestinal (gi) malformations are one of the most important congenital problems resulting in calf mortality within a few days of birth. the most common site for atresia, after the esophagus, is the jejunum. jejunum atresia is the congenital absence or complete blockage of a part of the jejunum lumen. early detection of intestinal obstruction is essential to prevent further complications. intestinal atresia is an underdiagnosed congenital defect in cattle. it results in complete occlusion of the intestinal lumen and, unless surgically corrected, results in death or euthanasia of the affected calf. there is limited information on the incidence of this condition or risk factors, including predisposing alleles, associated with the defect. atresia is a well-known congenital defect of the gastrointestinal system in calves and investigations into the etiology of this condition are warranted. domestication and selection have significantly changed the behavioral and phenotypic traits in modern domestic animals. the selection of animals by humans left detectable signatures on the genome of modern dairy cattle. the identification of these signals can help us to improve the genetic characteristics of economically important traits in goats. over the last decade, interest in the detection of genes or genomic regions that are targeted by selection has been growing. identifying signatures of selection can provide valuable insights about the genes or genomic regions that are or have been under selection pressure, which in turn leads to a better understanding of genotype-phenotype relationships. this study aimed to identify the selection signatures using the unbiased theta method associated with gastrointestinal atresia in holstein dairy calves.materials and methods: for calves with intestinal atresia, muscle tissue (>1 g) was collected from the latissimus dorsi muscle postmortem, and submerged in rna later solution. dna samples from 91 atresia cases and 377 control animals were genotyped using the illumina 777k bovinehd beadchip (illumina inc). the work described here is a case–control association study. single nucleotide polymorphism (snp) missing 5% of data, with maf of <1% and hardy–weinberg equilibrium p-values <10−6 were removed. the genotyping efficiency for samples was also verified, and samples with more than 5% missing data were removed. grouping was done to infer selection signatures based on fst statistic. the bioinformatics investigations were carried out using the ensembl database for bovine genes (assembly ars-ucd1.2), to identify potential candidate genes which already have been reported in/or surrounding genomic regions containing the peak of absolute extreme fst values. the regions corresponding to the upper and lower 0.01% of positive and negative obtained fst scores were considered regions under selection. genes within a 500-kb span of the start and end of the qtl were identified using ensembl 108 on the ars-ucd1.2 bovine genome assembly implemented in biomart. then, using the panther database, the general biological function of the genes was checked. at this stage, it is assumed that genes that belong to a functional class can be considered as a group of genes that have some specific and common characteristics, and the quantitative trait loci (qtls) in the selected region were extracted using the animalgenome database, and the genes were compared with other researches. genecards (http://www.genecards.org) and uniprotkb (http://www.uniprot.org) databases were also used to interpret the function of the obtained genes.results and discussion: with a 99.90 percentile threshold of the obtained theta (θ) values, eight genomic regions on chromosomes 7, 12, 13, 21, 22, 23 (two regions), and 29 in the holstein calves were identified. further investigation using bioinformatics tools showed these genomic regions overlapped with the genes (csf2, siah3, tmem14a, and skiv2l) associated with embryonic development, small intestine length, apoptosis, and several tumors. the population used in our study is small, owing to the challenge of collecting a substantial amount of blood on calves on commercial herds having received the diagnosis of gastrointestinal atresia and ready to be culled. diagnosis and culling of gastrointestinal atresia animals are ineffective preventive measures.  further work is required to identify which farm-specific or management risk factors contribute to the incidence of intestinal atresia.conclusions: the results of this study may provide an important source to facilitate the identification of genomic regions and then, the genes affecting gastrointestinal atresia in claves. however, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high-density genotyping.materials and methods: for calves with intestinal atresia, muscle tissue (>1 g) was collected from the latissimus dorsi muscle postmortem, submerged in rna later solution. dna samples from 91 atresia cases and 377 control animals were genotyped using the illumina 777k bovinehd beadchip (illumina inc). the work described here is a case–control association study. snp missing 5% of data, with maf of <1% and hardy–weinberg equilibrium p-values <10−6 were removed. the genotyping efficiency for samples was also verified, and samples with more than 5% missing data were removed. grouping was done to infer selection signatures based on fst statistic. the bioinformatics investigations were carried out using the ensembl database (cunningham et al., 2022) for bovine genes (assembly ars-ucd1.2), to identify potential candidate genes which already have been reported in/or surrounding genomic regions containing the peak of absolute extreme fst values. the regions corresponding to the upper and lower 0.01% of positive and negative obtained fst scores were considered as regions under selection. genes within a 500-kb span of the start and end of the qtl were identified using ensembl 108 on the ars-ucd1.2 bovine genome assembly implemented in biomart. then, using the panther database, the general biological function of the genes was checked. at this stage, it is assumed that genes that belong to a functional class can be considered as a group of genes that have some specific and common characteristics, and the qtls in the selected region were extracted using the animalgenome database, and the genes were compared with other researches. genecards (http://www.genecards.org) and uniprotkb (http://www.uniprot.org) databases were also used to interpret the function of the obtained genes.results and discussion with 99.90 percentile threshold of the obtained theta (θ) values, eight genomic regions on chromosomes 7, 12, 13, 21, 22, 23 (2 regions), and 29 in holstein calves breed were identified. further investigation using bioinformatics tools showed these genomic regions overlapped with the genes (csf2, siah3, tmem14a, skiv2l) associated with embryonic development, small intestine length, apoptosis, and several tumours. the population used in our study is small, owing to the challenge of collecting a substantial amount of blood on calves on commercial herds having received the diagnosis of gastrointestinal atresia and ready to be culled. diagnosis and culling of gastrointestinal atresia animals are ineffective preventive measures. further work is required to identify which farm-specific or management risk factors contribute to the incidence of intestinal atresia.conclusions: in conclusion, the results of this study may provide an important source to facilitate the identification of genomic regions and then, the genes affecting gastrointestinal atresia in claves. however, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping