Next Article in Journal
IGF-I Medium Supplementation Improves Singly Cultured Cat Oocyte Maturation and Embryo Development In Vitro
Previous Article in Journal
A DNA Regulatory Element Haplotype at Zinc Finger Genes Is Associated with Host Resilience to Small Ruminant Lentivirus in Two Sheep Populations
Previous Article in Special Issue
Gene Expression Profiling Reveals New Pathways and Genes Associated with Visna/Maedi Viral Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 7
10.3390/ani11071908
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed

by
Monika Olech
1,*,
Katarzyna Ropka-Molik
2,
Tomasz Szmatoła
2,3,
Katarzyna Piórkowska
2 and
Jacek Kuźmak
1
1
Department of Biochemistry, National Veterinary Research Institute, 24-100 Pulawy, Poland
2
Department of Animal Molecular Biology, National Research Institute of Animal Production, Krakowska 1, 32-083 Balice, Poland
3
Center for Experimental and Innovative Medicine, University of Agriculture in Krakow, Rędzina 1c, 30-248 Krakow, Poland
*
Author to whom correspondence should be addressed.
Animals 2021, 11(7), 1908; https://doi.org/10.3390/ani11071908
Submission received: 10 June 2021 / Accepted: 23 June 2021 / Published: 26 June 2021
(This article belongs to the Special Issue Recent Advances in Small Ruminant Lentiviruses)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The purpose of the study was to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs provirus concentration. A total of 14 SNPs were detected, 6 SNPs in the TLR7 gene locus and 8 SNPs in the TLR8 gene. These SNPs were located in intron, 3′UTR and 5′UTR regions and within the coding sequences leading to the synonymous mutations. Our results revealed that 9 out 14 identified polymorphisms were associated with the SRLVs proviral concentration. This finding supports a role for genetic variations of TLR7 and TLR8 in SRLVs infection.

Abstract

Toll-like receptors (TLRs) 7 and 8 are important in single-stranded viral RNA recognition, so genetic variation of these genes may play a role in SRLVs infection and disease progression. Present study aimed to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs provirus concentration as index of disease progression. A total of 14 SNPs were detected, 6 SNPs in the TLR7 gene locus and 8 SNPs in the TLR8 gene. Nine of the 14 identified polymorphisms, 4 in the TLR7 gene and 5 in TLR8 gene, were significantly associated with the SRLVs proviral concentration. These SNPs were located in 3′UTR, 5′UTR and intron sequences as well as in the coding sequences, but they led to silent changes. Homozygous genotypes of three TLR7 SNPs (synonymous variant 1:50703293, 3′UTR variant 1:50701297 and 5′UTR variant 1:50718645) were observed in goats with lower provirus copy number as well as in seronegative animals. The results obtained in this study suggest that SNPs of TLR7/TLR8 genes may induce differential innate immune response towards SRLVs affecting proviral concentration and thereby disease pathogenesis and progression. These findings support a role for genetic variations of TLR7 and TLR8 in SRLVs infection and warrants further studies on the effect of TLR7/TLR8 polymorphisms on SRLVs infection in different populations.

1. Introduction

Maedi visna virus (MVV) and caprine arthritis encephalitis virus (CAEV) also referred as small ruminant lentiviruses (SRLVs) are two related retroviruses which infect sheep and goats. These viruses infect monocytes, macrophages and dendritic cells and despite immune response cause a lifelong infection which can persist for months in latent or subclinical form. The most prevalent clinical signs as an outcome of SRLVs are associated with arthritis, neurological disorders, mastitis, emaciation and pneumonia [1]. Transmission occurs from infected dams to offspring by colostrum/milk consumption and between adults mainly through respiratory secretion [2]. Since no effective vaccines are available, the current control strategies against SRLVs infections are mostly based on the detection and culling of infected animals [3]. What is more, SRLVs infections have become a worldwide problem bringing considerable financial losses in the small ruminant industry [3,4,5]. In Poland, any SRLVs control programs have been never implemented, and thus, infections with SRLVs are quite common. The overall true prevalence at the flock level reached 33.3% and 71.9% in sheep and goats, respectively [6,7].
In the course of SRLVs infection, their genome is integrated into the host genome in the form of proviral DNA which load is a factor determining disease prediction. It was shown that high proviral load corresponded to higher lesion score [8,9]. On the other hand, the presence of animals with low SRLVs proviral load may suggest the implication of host factors that may restrict and control viral replication. Such animals which are also referred as long term non-progressor show competent humoral immune response in the absence of virus replication leading to lower SRLVs transmission [10,11,12]. This different reactivity to SRLVs infection may suggest different host immune response ability to control this infection.
In SRLVs-infected animals, both innate and adaptive immunity are induced. Moreover, several genes associated with resistance/susceptibility to SRLVs infection and disease outcome have been identified [13,14]. Toll-like receptors (TLRs) are host pattern recognition receptors (PRRs) that play a pivotal role in the innate immune system. PRRs activate the innate immune system through different signaling pathways upon recognition of variety of structurally conserved molecules derived from pathogens known as PAMPs (pathogen-associated molecular patterns). TLRs are type I transmembrane proteins which contain three domains, an ectodomain (ECD) containing leucine-rich repeat (LRR) motifs that mediate PAMPs recognition and a cytoplasmic Toll/interleukin-1 receptor (TIR) domain linked by a single transmembrane (TM) domain [15]. TLR involved in response to viral infection are TLR1-TLR4 and TLR6-TLR9 [16]. In retroviral infections, TLR7 and TLR8, which share a high degree of structural similarity and recognize viral ssRNA (single-stranded RNA) within endosomal compartments, are of particular interest [17]. Activation of these TLRs initiates the MyD88-dependent pathway, mainly in NF-κB activation for inducing the expression of pro-inflammatory cytokines, chemokines, and type I and type III interferons [18,19]. It is well known that cytokines impairing expression is a key point in SRLVs-dependent immunity, pathogenies, and appearance of lesions [20,21]. However, so far, the role of lentivirus-induced TLR signaling has not been widely studied in small ruminants. It was showed that in the course of SRLVs infection, TLR7 and TLR8 become activated, inducing IFN-α, IL-6, TNF-α production and expression of antiviral proteins [21]. Upregulation of TLR7 and TLR8 was noted in naturally SRLVs-infected sheep showing lung lesions [13].
Single nucleotide polymorphisms (SNPs) of genes encoding proteins involved in innate response have aroused much attention, and a number of studies have been conducted to identify SNPs in these genes in different species [22,23,24]. It has been found that genetic variants of TLRs can be associated with different outcome of several diseases [15,25]. Polymorphisms within TLR7 and TLR8 have been linked to susceptibility and progression of different viral diseases including these caused by Human Immunodeficiency virus (HIV-1), Hepatitis C virus (HCV), Chikungunya virus (CHIKV) and Crimean-Congo hemorrhagic fever (CCHF) virus [26,27,28,29]. Mikula et al. [30,31] suggested that mutations in TLR7, TLR8 and TLR9 may play an important role as host factor predisposing sheep for infection with SRLVs. Despite these observations, there is currently a lack of any studies reporting polymorphisms of caprine genes encoding TLR7 and TLR8 and their possible role in SRLVs infection. Herein, this study was conducted to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs infection and provirus concentration. We focused on goats of Carpathian breed because represent a remnant population of an ancient breed; thus, analysis in this population could be important for its conservation.

2. Materials and Methods

2.1. Animals and Sample Preparation

The study was performed in one flock, counting 32 adult goats representing Carpathian breed. All 32 goats were examined in this study. The goats of Carpathian breed were widely found in Carpathian Mountain in Poland in the 19th and 20th century and then became an extinct breed. In 2005, all of these goats were moved to the National Research Institute of Animal Production in Cracow where this flock has been recreated and currently is covered by the genetic resources protection program, supported by the Ministry of Agriculture and Rural Development; however, its risk of extinction is high.
All goats were clinically healthy and maintained at the same environmental and feeding conditions. Age of goats ranged from 2 to 10 years (average ~6 years). Blood samples were taken by jugular venipuncture in EDTA tubes and in second tubes for serum collection, and the blood was collected from all animals at the same day. DNA extraction was performed using Sherlock AX isolation kit (A&A Biotechnology, Gdynia, Poland) and quality of preparation was checked using NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA). Status of goats for SRLVs infection was confirmed by serological testing using ELISA (ID Screen MVV/CAEV Indirect Screening test, IDVet, France), according to the manufacturer’s recommendations. All procedures associated with animal handling and treatments were approved by the Local Ethical Committee on Animal Testing at the University of Life Sciences in Lublin (Poland).

2.2. SNP Identification–Variant Calling

As a first approach, the screening for SNPs identification was performed using already available RNA sequencing (RNA-seq) data obtained from 12 goats from tested flock (data available under GEO GSE168160 accession number). For SNPs identification, the quality of raw data was measured with the use of FastQC software [32] which was followed by trimming procedure that focused on removing adapter content, reads of low quality (phred quality < 20) and low read length (minimal read length set to 35) with the use of Trimmomatic software [33]. Then, the quality of trimmed reads was checked again to ensure the effectiveness of trimming. The next step was mapping procedure which was maintained with the use of Tophat software [34] to Capra Hircus ARS1 genome. Then, the duplicates were marked with the use of MarkDuplicates function of Picard Tools, and finally, the variant calling procedure was utilized with the use of Freebayes software [35] with a minimum coverage threshold set to 5. The filtration of the obtained variants was done with the use of VCFtools software [36] with the following parameters: minimum coverage set to 10; minimum quality set to 30; minimum combined coverage set to 120 and minimum combined quality set to 360. The annotation of the TLR8 and TLR7 gene variants was done with the use of Ensembl Variant Effect Predictor online software [37].

2.3. SNP Genotyping

Based on information on RNA-seq SNPs, for both TLR8 and TLR7 loci, the primers span gene’s regions containing polymorphisms were designed and PCR amplification, Sanger sequencing and SNPs analysis were performed on samples of all 32 goats. The primers span gene’s region containing SNPs was designed using Primer3 (v. 0.4.0) based on the reference sequence shown in Table 1 and Supplementary Figures S1 and S2. The PCR amplification products were obtained using AmpliTaq Gold™ 360 Master Mix (Thermo Fisher Scientific, Waltham, MA, USA) and purified using enzyme mixture-EPPiC (A&A Biotechnology, Gdynia, Poland). The Sanger sequencing was performed using BigDye™ Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific) and BigDye XTerminator™ Purification Kit (Thermo Fisher Scientific) according to the manufacturer’s protocol. The amplicons were capillary sequenced on 3500xL Genetic Analyzer (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA). Data were analyzed using Data Collection Software (Applied Biosystems). All identified variants were submitted to European Variation Archive (EVA) and received the Project ID PRJEB42246.

2.4. Poviral Load Quantification

DNA extracted from peripheral blood leukocytes (PBLs) [38] was quantified by the real time PCRs using Rotor-Gene Q Series ver. 2.0.3 (Qiagen) with primers and probe specifically designed for SRLV A5 subtype, which circulation in this flock was confirmed [39]. Sequence of forward and reverse primers and probe were CA5F (5′ TGGGAGTAGGACAAACAAATCA 3′), CA5R (5′ TGACATAT GCCTTACTGCTCTC 3′) and CA5P (5′ 6-FAM-TCACCCATTGTAGGCATAGCTGCC-BHQ-1 3′), respectively. A reference plasmid encompassing the target gag region was generated by the cloning of a 625 bp fragment into pDrive plasmid used to generate a standard curve based on 10-fold serial dilutions of plasmid DNA from 108 to 10. Amplification was performed in a total volume of 20 μL, according to the following cycling conditions: initial incubation and polymerase activation at 95 °C for 15 min and followed by 45 cycles of 94 °C for 60 s and 60 °C for 60 s. The reaction mixture for each PCR test contained 10 μL 2× QuantiTect Multiplex NoROX PCR buffer (Qiagen, Hilden, Germany), 400 nM of each of the primers, 200 nM of the specific probe, 5 μL of extracted genomic DNA. A non-template control (DEPC H2O) was included in each run. All samples were tested in duplicate, and the results were expressed as a mean copy number of provirus per 500 ng of genomic DNA of each goat.

2.5. Statistical Analysis

To identify association between genotypes and proviral load, a simple linear model approach was performed using t-test. In this analysis SNP was used as a classification variable and proviral load values were used as analysis variables. Quantitative results given as means ± SD presented specific significant differences between animals within given genotype groups. The association between SRLVs DNA proviral load and age of tested goats were estimated using the Kruskal–Wallis H test. For this purpose, the animals were divided into 4 groups: 1 (2–3 years old), 2 (4–6 years old), 3 (7–8 years old) and 4 (9–10 years old). p-value of <0.05 was considered statistically significant. All statistical analyses were performed using SAS Enterprise (Statistical Analysis System, Version 8.02, 2001, Wadowice, Poland).
The genotype distribution of the SNPs of TLR7 and TLR8 were tested for deviation from the Hardy-Weinberg equilibrium (HWE) by means of χ2 analysis using Court lab- HW calculator. p-value of <0.05 were considered statistically significant. Phase program was used to assess frequency of haplotype distribution, and Haploview program [40] was used (default mode) to generate haplotype block structures and calculate LD values between SNPs.

3. Results

3.1. Serology and qPCR

Out of 32 goats used in this study, 29 were found seropositive by ELISA and positive to quantitative polymerase chain reaction (qPCR) confirming the infection with SRLVs. Three goats were negative in both ELISA and qPCR. The average number of proviral copies in positive samples varied from 1 to 263 per 500 ng of genomic DNA. The correlation between proviral load and age of goats was not statistically significant (H = 6.893613; p < 0.0754; Kruskal-Wallis test).

3.2. TLR7 Gene-Detected SNPs, Allele and Haplotype Frequency

The variant calling method allowed to identify four polymorphisms in TLR7 gene locus. Two of them were 3′UTR (untranslated regions) variants (C/T 1:50701297 and T/C 1:50702074), one 5′UTR variant (C/T 1:50718645) and one synonymous polymorphism (T/C 1: 50703293) (Table 2). Additionally, Sanger sequencing performed on specimen collected from all goats revealed two additional SNPs, one in promoter and second in intron 1 regions (G/A 1: 50718760 and C/T 1: 50718466), which was not found in the first approach when whole blood transcriptomes were sequenced (Table 3).
The allele frequency analysis showed that two of all analyzed SNPs (promoter G/A 1:50718760 and 3′UTR T/C 1:50702074) were monomorphic but showed the mutant allele compared to the reference sequence. For the rest four polymorphisms, three genotypes were identified. The genotype and allele distributions of these four SNPs were shown in Table 4.
The Synonymous Variant (SV) T/C; 3′UTR and 5′UTR polymorphisms showed the similar allele distribution. These SNPs were most abundantly present in heterozygous goats (Table 4). Distribution of CC and TT genotype was from 16% to 31% and from 23% to 31%, respectively. For 3′UTR (C/T 1:50701297) and 5′UTR SNPs (C/T 1:50718645), only five animals, including 3 seronegative goats, carried CC genotype. For SV T/C 1:50703293 polymorphism, these goats carried TT genotype. The frequency of intron 1 (C/T 1:50718466) allele showed the predominance of CC genotype (91%), while opposite homozygotes TT and heterozygotes CT were 3% and 6%, respectively. One out 3 seronegative goats carried TT genotype while 2 other goats carried CC genotype. Moreover, the analyzed population was not in Hardy-Weinberg equilibrium according to intron C/T 1:50718466 (Table 4). Haplotype analysis allowed to detect 8 haplotype regions. Three haplotypes represented frequency from 10% to 51% while two haplotypes had frequency lower that 1% (Table 5). The linkage disequilibrium (LD) analysis showed the presence of one LD block and strong LD between four SNPs (Supplementary Materials: Table S1 and Figure S3).

3.3. TLR8 Gene-Detected SNPs, Allele and Haplotype Frequency

In TLR8 gene, using RNA-seq, six polymorphisms were detected from which two were the 3′UTR variants and the rest were the synonymous SNPs (Table 2). Additionally, Sanger sequencing allowed detecting two extra SNPs in the 3′UTR region: (1:50659346 C/A, 1:50659136 A/C) (Table 3).
In TLR8 locus, 3′UTR T/C 1:50659202 SNP was monomorphic but showed the mutant allele compared to the reference sequence. For 1:50664682 A/G SNP, two genotypes were identified. The frequencies of identified AA and GT genotypes were 95% and 5%, respectively. For the rest six polymorphisms, three genotypes were identified. One type of homozygote (AA or CC) was found in most goats, including also seronegative animals while heterozygotes goats accounted for 32–42% (Table 6).
The haplotype analysis revealed the presence of 10 haplotypes from which frequency ranged from 29.7% to 13.30% (Table 5). The linkage disequilibrium (LD) analysis showed the one LD block involving 6 polymorphisms (Supplementary Materials: Table S1 and Figure S3).

3.4. Association between SNPs and Provirus Copy Number

Because only three goats were uninfected with SRLVs, it was impossible to create control group with an equivalent number of animals to those of serologically positive and compare allele and genotype frequencies between infected and uninfected goats. Therefore, only associations between identified SNPs and provirus copy number of SRLVs were estimated for both TLR genes.
It was showed that nine of the identified polymorphisms showed significant association with provirus copy number. For TLR7 gene, four polymorphisms were significantly associated with SRLVs provirus copy number (Figure 1). For TLR7 intron variant 1:50718466, goats with CC genotype had significantly higher proviral load that CT genotype goats (p value < 0.05). For three other SNPs (synonymous variant 1:50703293, 3′UTR variant 1:50701297 and 5′UTR variant 1:50718645), the heterozygote goats CT were characterized by higher provirus copy number (p < 0.05).
For TLR8 gene, five polymorphisms were significantly associated with provirus copy number (Figure 2). For TLR8 1:50666071 synonymous variant, homozygotes GG goats had lower proviral load than heterozygotes AG goats (p < 0.05). For two synonymous TLR8 polymorphisms (1:50664064 and 1:50664208), the goats with the CC genotype showed the lower copy number of SRLVs compared to CT animals. For 3′UTR variant 1:50659346 AC genotype, goats had significantly higher proviral load than CC genotype goats (p value < 0.05). The exact opposite was the case with 1:50659136 SNP, for which homozygotes AA goats had significantly higher SRLVs proviral load than both other genotypes (p < 0.01).

4. Discussion

The present study is the first attempt showing the genetic variability in TLR7 and TLR8 genes in SRLVs infected goats as well as association of SNPs with SRLVs proviral concentration.
We identified 6 polymorphisms in the TLR7 gene locus and 8 SNPs in the TLR8 gene. These SNPs were located in intron, 3′UTR and 5′UTR regions and within the coding sequences leading to the synonymous mutations which do not result in any amino acid change in encoded protein, due to genetic code redundancy. Synonymous substitutions and mutations, affecting noncoding DNA regions, are often considered silent mutations; however, it is not always the case that the mutation is silent. This kind of polymorphisms can produce different effects on gene expression leading to functional differences of various significance. Several recent reports consider such mutations, particularly concerning human diseases [41,42]. Silent mutations may affect translational efficiency and protein folding by changing codons read by tRNAs [43,44] and alter mRNA stability structure or splicing leading to variation in protein expression [45]. SNPs in the promoter regions could affect their activity and regulation producing changes in gene expression levels. Moreover, through different mechanisms, silent mutations may influence gene regulation, differences in mRNA and protein abundance, and proteins’ structure and functionality [45]. Untranslated regions (UTRs), which are localized at both ends of transcript (mRNA), make numerous conformational structures, tridimensional loops and hairpins that interact with numerous proteins and other functional and regulatory compounds like ribosomes or microRNA. UTRs are known to play crucial roles in the post-transcriptional regulation of gene expression. The main role of 5′UTR is controlling of translation efficiency as well as transcript stabilization while 3′UTR is mostly implicated in regulation of transcript stabilization, including modulation of the transport of mRNAs out of the nucleus [46,47]. Furthermore 3′UTR can be microRNA target sites, where microRNAs bind and regulate genetic expression. Therefore, SNPs or mutations in these regions might alter the existing target sites for microRNAs [48]. The importance of UTRs in regulating gene expression is underlined by the finding that mutations that alter the UTR can lead to serious pathology [49].
The role of TLR7 and TLR8 SNPs in SRLVs infection surely requires a comparison of infected and uninfected animals. Unfortunately, only 3 out of 32 tested goats were uninfected, so statistical analysis was not possible. However, it is worth pointing out that goats from tested flock were serologically examined several times over the past few years, and only these 3 goats remained always negative. Serological studies revealed that very young goats presented antibodies against SRLVs while these 3 goats, which have 3, 7 and 10 years old, were uninfected. This clearly indicates an implication of host factors that may restrict and control SRLVs replication. Because it was impossible to create control group with an equivalent number of animals to those of serologically positive and compare allele and genotype frequencies between infected and uninfected goats, we focused on investigating the association between TLR7 and TLR8 mutations and SRLVs proviral concentration, as important index of disease progression. It has been demonstrated that host control of infection with SRLVs, including provirus level, may have a genetic basis [10,50,51,52]. Genetic studies in humans have pointed to a role for TLR SNPs, especially TLR7, in HIV infection and their association with viral loads. The presence of the most frequent TLR7 polymorphisms, TLR7 Gln11Leu (rs179008), was associated with increased viral load and altered CD4 T cell counts during HIV infection [53,54] while polymorphisms in TLR7 (re179010) and TLR8 (rs3764880) may reduce the risk of disease by participating in inhibition of viral load leading to the slower progression of infection [28,55,56]. Our results revealed that 9 out 14 identified polymorphisms were significantly associated with the SRLVs proviral concentration. In particular, homozygous genotypes of three TLR7 SNPs (synonymous variant 1:50703293, 3′UTR variant 1:50701297 and 5′UTR variant 1:50718645) were observed in goats with lower provirus copy number as well as in seronegative animals. However, additional analysis is needed to confirm this association. In this analysis, the effort should be directed to goat’s samples collected at different time points to check the evolution of proviral load over the time.
The direct influence of the SNPs detected in this study on phenotypic features is unknown. Similar to what was showed for SRLVs infection in sheep [57], it can be assumed that during SRLVs infection in goats, TLR7 and TLR8 became activated, inducing production of antiviral cytokine, including IFN5-α, IL-6 and TNF-α production, which are critical in the development of antiviral immune response [21]. Activation of IFN type I pathway leads to expression of genes with IFN-sensitive response elements (ISREs) which fight viruses through different mechanisms. It has been confirmed that SRLVs contain ISRE and that IFN may modulate viral transcription promoting an antiviral state [14,58]. Therefore, we can speculate that SNPs detected in this study may affect the level of TLR7/8 expression, causing the differences in the production of downstream cytokines, like IFN-α, which finally can lead to the differences in SRLVs proviral concentration.
Studies on the association between TLR variants and HIV infection/disease progression suggested that they may be specific for different human populations [59]. Diversity reported between breed’s susceptibility/resistance to SRLVs [50,60] may suggest that SNPs detected in this study may be typical for goats of Carpathian breed. We especially focused on goats of this breed because the Carpathian goat is an ancient breed which in the 19th/20th century was present in Carpathian Mountain in Poland and then became an extinct breed. It emphasizes the need to select animals for higher resistance to infections. Considering the fact that animals belonging to this flock showed close to 100% seropositivity against SRLVs, the identification of genetic markers associated with SRLVs proviral load may be important for further breeding program. Removal of infected animals, especially those with high proviral load can limit the spread of the virus since such animals are highly efficient in shedding the virus [10].
In conclusion, this is the first report showing single nucleotide polymorphisms in genes encoding TLR7/8 in goats and the association between TLR7/8 polymorphisms and SRLVs proviral concentration in goats of Carpathian breed. Limited number of animals tested in this study and lack of possibility to compare infected and non-infected animals are undoubtedly limitations of this work. However, the obtained results suggest that SNPs of TLR7/TLR8 genes may induce differential innate immune response towards SRLVs affecting proviral concentration and thereby disease progression. These findings support a role of genetic variations in TLR7 and TLR8 in the course of infection with SRLVs and warrants further studies on TLR7/TLR8 polymorphisms.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11071908/s1, Figure S1: The exact TLR7 SNPs localization based on transcript and genomic sequence, Figure S2: The exact TLR8 SNPs localization based on transcript and genomic sequence, Figure S3: The visualization of LD and haplotype blocks based on Haploview software for TLR7 and TLR8 genes, Table S1: The linkage disequilibrium (LD) between identified SNPs in TLRs genes.

Author Contributions

Conceptualization, J.K., M.O. and K.R.-M.; methodology, K.R.-M. and K.P.; software, T.S.; formal analysis, M.O. and K.R.-M.; investigation, K.P., K.R.-M. and M.O.; writing—original draft preparation, M.O. and K.R.-M.; writing—review and editing, M.O., K.R.-M., J.K., T.S. and K.P.; visualization, K.R.-M. and M.O.; supervision, J.K. and K.R.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

All procedures associated with animal handling and treatments were approved (No 37/2016) by the Local Ethical Committee on Animal Testing at the University of Life Sciences in Lublin (Poland), and written informed consent was obtained from the owner.

Informed Consent Statement

Not applicable.

Data Availability Statement

RNA-seq data are available under GEO GSE168160 accession number. Rest of data generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Polledo, L.; González, J.; Benavides, J.; Morales, S.; Martínez-Fernández, B.; Delgado, L.; Reina, R.; Glaria, I.; Pérez, V.; Ferreras, M.C.; et al. Patterns of lesion and local host cellular immune response in natural cases of ovine maedi-visna. J. Comp. Pathol. 2012, 147, 1–10. [Google Scholar] [CrossRef]
  2. Minguijón, E.; Reina, R.; Pérez, M.; Polledo, L.; Villoria, M.; Ramírez, H.; Leginagoikoa, I.; Badiola, J.J.; García-Marín, J.F.; De Andrés, D.; et al. Small ruminant lentivirus infections and diseases. Vet. Microbiol. 2015, 181, 75–89. [Google Scholar] [CrossRef] [PubMed]
  3. Kalogianni, A.I.; Bossis, I.; Ekateriniadou, L.V.; Gelasakis, A.I. Etiology, Epizootiology and Control of Maedi-Visna in Dairy Sheep: A Review. Animals 2020, 10, 616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Reina, R.; Berriatua, E.; Luján, L.; Juste, R.; Sánchez, A.; De Andrés, D.; Amorena, B. Prevention strategies against small ruminant lentiviruses: An update. Vet. J. 2009, 182, 31–37. [Google Scholar] [CrossRef] [PubMed]
  5. Pérez, M.; Biescas, E.; De Andrés, X.; Leginagoikoa, I.; Salazar, E.; Berriatua, E.; Reina, R.; Bolea, R.; De Andrés, D.; Juste, R.A.; et al. Visna/maedi virus serology in sheep: Survey, risk factors and implementation of a successful control programme in Aragón (Spain). Vet. J. 2010, 186, 221–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Olech, M.; Osiński, Z.; Kuźmak, J. Bayesian estimation of seroprevalence of small ruminant lentiviruses in sheep from Poland. Prev. Vet. Med. 2017, 147, 66–78. [Google Scholar] [CrossRef] [PubMed]
  7. Kaba, J.; Czopowicz, M.; Ganter, M.; Nowicki, M.; Witkowski, L.; Nowicka, D.; Szaluś-Jordanow, O. Risk factors associated with seropositivity to small ruminant lentiviruses in goat herds. Res. Vet. Sci. 2013, 94, 225–227. [Google Scholar] [CrossRef]
  8. Ravazzolo, A.P.; Nenci, C.; Vogt, H.R.; Waldvogel, A.; Obexer-Ruff, G.; Peterhans, E.; Bertoni, G. Viral load, organ distribution, histopathological lesions, and cytokine mRNA expression in goats infected with a molecular clone of the caprine arthritis encephalitis virus. Virology 2006, 350, 116–127. [Google Scholar] [CrossRef] [PubMed]
  9. Herrmann-Hoesing, L.M.; Noh, S.M.; White, S.N.; Snekvik, K.R.; Truscott, T.; Knowles, D.P. Peripheral ovine progressive pneumonia provirus levels correlate with and predict histological tissue lesion severity in naturally infected sheep. Clin. Vaccine Immunol. 2009, 16, 551–557. [Google Scholar] [CrossRef] [Green Version]
  10. Crespo, H.; Bertolotti, L.; Proffiti, M.; Cascio, P.; Cerruti, F.; Acutis, P.L.; De Andrés, D.; Reina, R.; Rosati, S. Low proviral small ruminant lentivirus load as biomarker of natural restriction in goats. Vet. Microbiol. 2016, 192, 152–162. [Google Scholar] [CrossRef]
  11. Stonos, N.; Wootton, S.K.; Karrow, N. Immunogenetics of small ruminant lentiviral infections. Viruses 2014, 6, 3311–3333. [Google Scholar] [CrossRef] [PubMed]
  12. Cardinaux, L.; Zahno, M.L.; Deubelbeiss, M.; Zanoni, R.; Vogt, H.R.; Bertoni, G. Virological and phylogenetic characterization of attenuated small ruminant lentivirus isolates eluding efficient serological detection. Vet. Microbiol. 2013, 162, 572–581. [Google Scholar] [CrossRef]
  13. Larruskain, A.; Bernales, I.; Luján, L.; De Andrés, D.; Amorena, B.; Jugo, B.M. Expression analysis of 13 ovine immune response candidate genes in Visna/Maedi disease progression. Comp. Immunol. Microbiol. Infect. Dis. 2013, 36, 405–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gomez-Lucia, E.; Barquero, N.; Domenech, A. Maedi-Visna virus: Current perspectives. Vet. Med. 2018, 21, 11–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Skevaki, C.; Pararas, M.; Kostelidou, K.; Tsakris, A.; Routsias, J.G. Single nucleotide polymorphisms of Toll-like receptors and susceptibility to infectious diseases. Clin. Exp. Immunol. 2015, 180, 165–177. [Google Scholar] [CrossRef] [Green Version]
  16. Lester, S.N.; Li, K. Toll-like receptors in antiviral innate immunity. J. Mol. Biol. 2014, 426, 1246–1264. [Google Scholar] [CrossRef]
  17. Heil, F.; Hemmi, H.; Hochrein, H.; Ampenberger, F.; Kirschning, C.; Akira, S.; Lipford, G.; Wagner, H.; Bauer, S. Species-specific recognition of single-stranded RNA via toll-like receptor 7 and 8. Science 2004, 303, 1526–1529. [Google Scholar] [CrossRef] [Green Version]
  18. Akira, S.; Takeda, K. Toll-like receptor signalling. Nat. Rev. Immunol. 2004, 4, 499–511. [Google Scholar] [CrossRef]
  19. De Pablo-Maiso, L.; Doménech, A.; Echeverría, I.; Gómez-Arrebola, C.; De Andrés, D.; Rosati, S.; Gómez-Lucia, E.; Reina, R. Prospects in Innate Immune Responses as Potential Control Strategies against Non-Primate Lentiviruses. Viruses 2018, 10, 435. [Google Scholar] [CrossRef] [Green Version]
  20. Jarczak, J.; Kaba, J.; Reczyńska, D.; Bagnicka, E. Impaired Expression of Cytokines as a Result of Viral Infections with an Emphasis on Small Ruminant Lentivirus Infection in Goats. Viruses 2016, 8, 186. [Google Scholar] [CrossRef]
  21. Blacklaws, B.A. Small ruminant lentiviruses: Immunopathogenesis of visna-maedi and caprine arthritis and encephalitis virus. Comp. Immunol. Microbiol. Infect. Dis. 2012, 35, 259–269. [Google Scholar] [CrossRef]
  22. Fang, X.; Liu, X.; Meng, C.; Fu, Y.; Wang, X.; Li, B.; Tu, F.; Zhao, F.; Ren, S. Breed-linked polymorphisms of porcine toll-like receptor 2 (TLR2) and TLR4 and the primary investigation on their relationship with prevention against Mycoplasma pneumoniae and bacterial LPS challenge. Immunogenetics 2013, 65, 829–834. [Google Scholar] [CrossRef]
  23. Kloch, A.; Wenzel, M.A.; Laetsch, D.R.; Michalski, O.; Bajer, A.; Behnke, J.M.; Welc-Falęciak, R.; Piertney, S.B. Signatures of balancing selection in toll-like receptor (TLRs) genes–novel insights from a free-living rodent. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef] [PubMed]
  24. Novák, K. Functional polymorphisms in Toll-like receptor genes for innate immunity in farm animals. Vet. Immunol. Immunopathol. 2014, 157, 1–11. [Google Scholar] [CrossRef]
  25. Beutler, B.; Jiang, Z.; Georgel, P.; Crozat, K.; Croker, B.; Rutschmann, S.; Du, X.; Hoebe, K. Genetic analysis of host resistance: Toll-like receptor signaling and immunity at large. Annu. Rev. Immunol. 2006, 24, 353–389. [Google Scholar] [CrossRef] [Green Version]
  26. Dutta, S.K.; Tripathi, A. Association of toll-like receptor polymorphisms with susceptibility to chikungunya virus infection. Virology 2017, 511, 207–213. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, C.H.; Eng, H.L.; Lin, K.H.; Chang, C.H.; Hsieh, C.A.; Lin, Y.L.; Lin, T.M. TLR7 and TLR8 gene variations and susceptibility to hepatitis C virus infection. PLoS ONE 2011, 6, e26235. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, T.; Zhu, J.; Su, B.; Cao, L.; Li, Z.; Wei, H.; Huang, X.; Zheng, K.; Li, A.; Chen, N.; et al. Effects of TLR7 Polymorphisms on the Susceptibility and Progression of HIV-1 Infection in Chinese MSM Population. Front. Immunol. 2020, 11, 589010. [Google Scholar] [CrossRef]
  29. Arslan, S.; Engin, A.; Özbilüm, N.; Bakır, M. Toll-like receptor 7 Gln11Leu, c.4-151A/G, and +1817G/T polymorphisms in Crimean Congo hemorrhagic fever. J. Med. Virol. 2015, 87, 1090–1095. [Google Scholar] [CrossRef] [PubMed]
  30. Mikula, I.; Bhide, M.; Pastorekova, S.; Mikula, I. Characterization of ovine TLR7 and TLR8 protein coding regions, detection of mutations and Maedi Visna virus infection. Vet. Immunol. Immunopathol. 2010, 138, 51–59. [Google Scholar] [CrossRef]
  31. Mikula, I., Jr.; Mikula, I., Sr. Characterization of ovine Toll-like receptor 9 protein coding region, comparative analysis, detection of mutations and maedi visna infection. Dev. Comp. Immunol. 2011, 35, 182–192. [Google Scholar] [CrossRef]
  32. De Sena Brandine, G.; Smith, A.D. Falco: High-speed FastQC emulation for quality control of sequencing data. F1000Research 2019, 8, 1874. [Google Scholar] [CrossRef] [Green Version]
  33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  34. Ghosh, S.; Chan, C.K. Analysis of RNA-Seq Data Using TopHat and Cufflinks. Methods Mol. Biol. 2016, 1374, 339–361. [Google Scholar]
  35. Sánchez-Sevilla, J.F.; Vallarino, J.G.; Osorio, S.; Bombarely, A.; Posé, D.; Merchante, C.; Botella, M.A.; Amaya, I.; Valpuesta, V. Gene expression atlas of fruit ripening and transcriptome assembly from RNA-seq data in octoploid strawberry (Fragaria × ananassa). Sci. Rep. 2017, 7, 13737. [Google Scholar] [CrossRef]
  36. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. 1000 Genomes Project Analysis Group. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  37. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [Green Version]
  38. Olech, M.; Valas, S.; Kuźmak, J. Epidemiological survey in single-species flocks from Poland reveals expanded genetic and antigenic diversity of small ruminant lentiviruses. PLoS ONE 2018, 13, e0193892. [Google Scholar] [CrossRef] [Green Version]
  39. Olech, M.; Kuźmak, J. Molecular Characterization of Small Ruminant Lentiviruses of Subtype A5 Detected in Naturally Infected but Clinically Healthy Goats of Carpathian Breed. Pathogens 2020, 9, 992. [Google Scholar] [CrossRef] [PubMed]
  40. Gabriel, S.B.; Schaffner, S.F.; Nguyen, H.; Moore, J.M.; Roy, J.; Blumenstiel, B.; Higgins, J.; DeFelice, M.; Lochner, A.; Faggart, M.; et al. The structure of haplotype blocks in the human genome. Science 2002, 296, 2225–2229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Sauna, Z.E.; Kimchi-Sarfaty, C. Understanding the contribution of synonymous mutations to human disease. Nat. Rev. Genet. 2011, 12, 683–691. [Google Scholar] [CrossRef]
  42. Fåhraeus, R.; Marin, M.; Olivares-Illana, V. Whisper mutations: Cryptic messages within the genetic code. Oncogene 2016, 35, 3753–3759. [Google Scholar] [CrossRef]
  43. Hunt, R.C.; Simhadri, V.L.; Iandoli, M.; Sauna, Z.E.; Kimchi-Sarfaty, C. Exposing synonymous mutations. Trends Genet. 2014, 30, 308–321. [Google Scholar] [CrossRef]
  44. Fernández-Calero, T.; Cabrera-Cabrera, F.; Ehrlich, R.; Marín, M. Silent Polymorphisms: Can the tRNA Population Explain Changes in Protein Properties? Life 2016, 6, 9. [Google Scholar] [CrossRef] [Green Version]
  45. Diederichs, S.; Bartsch, L.; Berkmann, J.C.; Fröse, K.; Heitmann, J.; Hoppe, C.; Iggena, D.; Jazmati, D.; Karschnia, P.; Linsenmeier, M.; et al. The dark matter of the cancer genome: Aberrations in regulatory elements, untranslated regions, splice sites, non-coding RNA and synonymous mutations. EMBO Mol. Med. 2016, 8, 442–457. [Google Scholar] [CrossRef]
  46. Carmody, S.R.; Wente, S.R. mRNA nuclear export at a glance. J. Cell Sci. 2009, 122, 1933–1937. [Google Scholar] [CrossRef] [Green Version]
  47. Pickering, B.M.; Willis, A.E. The implications of structured 5′ untranslated regions on translation and disease. Semin. Cell Dev. Biol. 2005, 16, 39–47. [Google Scholar] [CrossRef]
  48. Georges, M.; Clop, A.; Marcq, F.; Takeda, H.; Pirottin, D.; Hiard, S.; Tordoir, X.; Caiment, F.; Meish, F.; Bibé, B.; et al. Polymorphic microRNA-target interactions: A novel source of phenotypic variation. Cold Spring Harb. Symp. Quant. Biol. 2006, 71, 343–350. [Google Scholar] [CrossRef] [Green Version]
  49. Conne, B.; Stutz, A.; Vassalli, J.D. The 3′ untranslated region of messenger RNA: A molecular ‘hotspot’ for pathology? Nat. Med. 2000, 6, 637–641. [Google Scholar] [CrossRef]
  50. Herrmann-Hoesing, L.M.; White, S.N.; Mousel, M.R.; Lewis, G.S.; Knowles, D.P. Ovine progressive pneumonia provirus levels associate with breed and Ovar-DRB1. Immunogenetics 2008, 60, 749–758. [Google Scholar] [CrossRef]
  51. White, S.N.; Mousel, M.R.; Reynolds, J.O.; Lewis, G.S.; Herrmann-Hoesing, L.M. Common promoter deletion is associated with 3.9-fold differential transcription of ovine CCR5 and reduced proviral level of ovine progressive pneumonia virus. Anim. Genet. 2009, 40, 583–589. [Google Scholar] [CrossRef]
  52. Alshanbari, F.A.; Mousel, M.R.; Reynolds, J.O.; Herrmann-Hoesing, L.M.; Highland, M.A.; Lewis, G.S.; White, S.N. Mutations in Ovis aries TMEM154 are associated with lower small ruminant lentivirus proviral concentration in one sheep flock. Anim. Genet. 2014, 45, 565–571. [Google Scholar] [CrossRef] [PubMed]
  53. Said, E.A.; Al-Yafei, F.; Zadjali, F.; Hasson, S.S.; Al-Balushi, M.S.; Al-Mahruqi, S.; Koh, C.Y.; Al-Naamani, K.; Al-Busaidi, J.Z.; Idris, M.A.; et al. Association of single-nucleotide polymorphisms in TLR7 (Gln11Leu) and TLR9 (1635A/G) with a higher CD4T cell count during HIV infection. Immunol. Lett. 2014, 160, 58–64. [Google Scholar] [CrossRef]
  54. Anokhin, V.V.; Bakhteeva, L.B.; Khasanova, G.R.; Khaiboullina, S.F.; Martynova, E.V.; Tillett, R.L.; Schlauch, K.A.; Lombardi, V.C.; Rizvanov, A.A. Previously Unidentified Single Nucleotide Polymorphisms in HIV/AIDS Cases Associate with Clinical Parameters and Disease Progression. Biomed. Res. Int. 2016, 2016, 2742648. [Google Scholar] [CrossRef]
  55. Browne, E.P. The Role of Toll-Like Receptors in Retroviral Infection. Microorganisms 2020, 8, 1787. [Google Scholar] [CrossRef]
  56. Oh, D.Y.; Taube, S.; Hamouda, O.; Kücherer, C.; Poggensee, G.; Jessen, H.; Eckert, J.K.; Neumann, K.; Storek, A.; Pouliot, M.; et al. A functional toll-like receptor 8 variant is associated with HIV disease restriction. J. Infect. Dis. 2008, 198, 701–709. [Google Scholar] [CrossRef] [Green Version]
  57. Larruskain, A.; Jugo, B.M. Retroviral infections in sheep and goats: Small ruminant lentiviruses and host interaction. Viruses 2013, 5, 2043–2061. [Google Scholar] [CrossRef] [Green Version]
  58. Murphy, B.; Hillman, C.; Castillo, D.; Vapniarsky, N.; Rowe, J. The presence or absence of the gamma-activated site determines IFN gamma-mediated transcriptional activation in CAEV promoters cloned from the mammary gland and joint synovium of a single CAEV-infected goat. Virus Res. 2012, 163, 537–545. [Google Scholar] [CrossRef]
  59. Willie, B.; Hall, N.B.; Stein, C.M.; Jurevic, R.J.; Weinberg, A.; Mehlotra, R.K.; Zimmerman, P.A. Association of Toll-like receptor polymorphisms with HIV status in North Americans. Genes Immun. 2014, 15, 569–577. [Google Scholar] [CrossRef] [Green Version]
  60. Peterhans, E.; Greenland, T.; Badiola, J.; Harkiss, G.; Bertoni, G.; Amorena, B.; Eliaszewicz, M.; Juste, R.A.; Krassnig, R.; Lafont, J.-P.; et al. Routes of transmission and consequences of small ruminant lentiviruses (SRLVs) infection and eradication schemes. Vet. Res. 2004, 35, 257–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Animals 11 01908 g001 550
Figure 1. The association between genotype within TLR7 genes and SRLVs copy number. *: p < 0.05.
Figure 1. The association between genotype within TLR7 genes and SRLVs copy number. *: p < 0.05.
Animals 11 01908 g001
Animals 11 01908 g002a 550Animals 11 01908 g002b 550
Figure 2. The association between genotype within TLR8 genes and SRLVs copy number. *: p < 0.05; **: p < 0.01.
Figure 2. The association between genotype within TLR8 genes and SRLVs copy number. *: p < 0.05; **: p < 0.01.
Animals 11 01908 g002aAnimals 11 01908 g002b
Table
Table 1. Primers used for amplification of TLR7 and TLR8 polymorphisms.
Table 1. Primers used for amplification of TLR7 and TLR8 polymorphisms.
Reference SequenceSequences of Primers (5′-3′)OrientationSizeIdentified SNP Based on RNA-Seq Data
TLR7
ENSCHIG00000021012		GACATGAGGTTGCCCTGATTGF338 bpLWLT01000021.1_50718645_C/T
TATGATTTGGGTCCCTTCCCCR
TCTCTGAGTTTCTTACCTTTGGGA F285 bpLWLT01000021.1_50703293_T/C
GCTGAGGTCCAGATGTCGCR
CTTGATCGGCTCACCCATGCF297 bpLWLT01000021.1_50702074_T/C
AAGTGCAGTTCTTGGTGACTGR
ACCTGCTTAAATATGGCTCTGGF418 bpLWLT01000021.1_50701297_C/T
TGGATGAGGCAACTTTTCTTTGGR
TLR8
ENSCHIG00000016625		GGGCATTTCTCAACCTCAAAF316 bpLWLT01000021.1_50666071_A/G
GTTAGCGAGCTTGGCAGACTR
AATGCGCACTATTTCCGAATF301 bpLWLT01000021.1_50664682_A/G
TCCTGGGCAAGTTAAGGAAGR
GCAGCCTGATACACCTCGATF424 bpLWLT01000021.1_50664208_C/T LWLT01000021.1_50664064_C/T
TGGGATGTGGACAGAGACCTR
CTGATGTTGCTGGGAGTCCTF576 bpLWLT01000021.1_50659559_C/T
LWLT01000021.1_50659202_T/C
TAGCTGAGAGGGGAATTGCCR
F—forward, R—reverse.
Table
Table 2. New polymorphisms within goat TLR7 and TLR8 genes identified using RNA-seq data.
Table 2. New polymorphisms within goat TLR7 and TLR8 genes identified using RNA-seq data.
LocationMutant AlleleVariant TypeGene ReferenceTranscript ReferencecDNA PositionCodons
TLR7
LWLT01000021.1_50701297_C/TLWLT01000021.1:50701297-50701297T3_prime_UTR_variantENSCHIG00000021012ENSCHIT00000031299.14304-
LWLT01000021.1_50702074_T/CLWLT01000021.1:50702074-50702074C3_prime_UTR_variantENSCHIG00000021012ENSCHIT00000031299.13527-
LWLT01000021.1_50703293_T/CLWLT01000021.1:50703293-50703293Csynonymous_variantENSCHIG00000021012ENSCHIT00000031309.12088ggA/ggG
LWLT01000021.1_50718645_C/TLWLT01000021.1:50718645-50718645T5_prime_UTR_variantENSCHIG00000021012ENSCHIT00000031299.123-
TLR8
LWLT01000021.1_50659202_T/CLWLT01000021.1:50659202-50659202C3_prime_UTR_variantENSCHIG00000016625ENSCHIT00000024133.14708-
LWLT01000021.1_50659559_C/TLWLT01000021.1:50659559-50659559T3_prime_UTR_variantENSCHIG00000016625ENSCHIT00000024133.14351-
LWLT01000021.1_50664064_C/TLWLT01000021.1:50664064-50664064Tsynonymous_variantENSCHIG00000016625ENSCHIT00000024133.12460gcG/gcA
LWLT01000021.1_50664208_G/ALWLT01000021.1:50664208-50664208Asynonymous_variantENSCHIG00000016625ENSCHIT00000024133.12316ttC/ttT
LWLT01000021.1_50664682_A/GLWLT01000021.1:50664682-50664682Gsynonymous_variantENSCHIG00000016625ENSCHIT00000024133.11842ctT/ctC
LWLT01000021.1_50666071_A/GLWLT01000021.1:50666071-50666071Gsynonymous_variantENSCHIG00000016625ENSCHIT00000024133.1453aaT/aaC
Table
Table 3. New polymorphisms within goat TLR7 and TLR8 genes identified using Sanger sequencing method.
Table 3. New polymorphisms within goat TLR7 and TLR8 genes identified using Sanger sequencing method.
LocationMutant AlleleVariant TypeGene ReferenceTranscript ReferencecDNA Position
TLR7
LWLT01000021.1_50718760_G/ALWLT01000021.1:50718760-1: 50718760GPromoterENSCHIG00000021012ENSCHIT00000031299.1-
LWLT01000021.1_0718466_C/TLWLT01000021.1:50717866-50717866CIntron 1ENSCHIG00000021012ENSCHIT00000031299.1-
TLR8
LWLT01000021.1_50659346_C/ALWLT01000021.1: 50659346-50659346C3_prime_UTR_variantENSCHIG00000016625ENSCHIT00000024133.14564
LWLT01000021.1_50659136_A/CLWLT01000021.1: 50659136-50659136A3_prime_UTR_variantENSCHIG00000016625ENSCHIT00000024133.14774
Table
Table 4. The genotype distribution for SNPs detected in TLR7 gene.
Table 4. The genotype distribution for SNPs detected in TLR7 gene.
SNPGenotypesAllelesHWE p Value
CCCTTTCT
SV T/C 1:5070329311 (31%)16 (46%)5 (23%)38 (0.54)26 (0.46)ns
CCCTTTCT
3′UTR C/T 1:507012975 (23%)16 (46%)11 (31%)38 (0.54)26 (0.46)ns
CCCTTTCT
5′UTR C/T 1:507186455 (16%)17 (53%)10 (31%)27 (0.42)37 (0.58)ns
CCCTTTCT
INTRON1 C/T 1:5071846629 (91%)2 (6%)1 (3%)58 (0.94)6 (0.06)0.00001
HWE—Hardy-Weinberg equilibrium; ns—not significant.
Table
Table 5. The haplotypes frequency of TLR7 and TLR8 genes.
Table 5. The haplotypes frequency of TLR7 and TLR8 genes.
HaplotypeFrequencyS.E.
TLR7
1T/T/C/T/A/C51%0.003
2C/T/T/C/A/C29%0.004
3C/T/T/C/A/T10%0.004
4T/T/T/C/A/C4%0.002
5T/T/C/C/A/C4%0.0003
6T/T/C/T/G/C1.20%0.000
7T/T/T/C/A/T<1%0.002
8T/T/C/T/A/T<1%0.003
TLR8
1C/C/C/G/A/A/C/A29.70%0.010
2C/T/T/A/A/G/A/A19.90%0.003
3C/C/C/G/A/A/C/C14.80%0.008
4C/T/C/G/A/A/C/A13.40%0.007
5C/T/C/G/A/A/C/C13.30%0.008
6C/T/T/A/A/G/C/A5.70%0.001
7C/T/C/G/G/A/C/A2.70%0.003
8C/C/T/A/A/G/A/A<1%0.002
9C/T/C/G/G/A/C/C<1%0.003
10C/T/T/A/A/G/A/C<1%0.002
SNPs order: TLR7-LWLT01000021.1-50701297; 50702074; 50703293; 50718645; 50718760; 50718466; TLR8–LWLT01000021.1 -50659202; 50659559; 50664064; 50664208; 50664682; 50666071; 50659346; 50659136.
Table
Table 6. The genotype distribution for SNPs detected in TLR8 gene.
Table 6. The genotype distribution for SNPs detected in TLR8 gene.
SNPGenotypesAllelesHWE p Value
AAGAGGAG
SV A/G 1:5066607118 (56%)11 (34%)3 (10%)47 (0.73)17 (0.27)ns
CCCTTTGA
SV G/A 1:50664220818 (56%)11 (34%)3 (10%)47 (0.73)17 (0.27)ns
CCCTTTCT
SV C/T 1:5066406418 (56%)11 (34%)3 (10%)47 (0.73)17 (0.27)ns
AAACCCAC
3′UTR A/C 1:5065913617 (55%)10 (32%)4 (13%)42 (0.71)18 (0.29)ns
CCACAACA
3′UTR C/A 1:5065934617 (55%)13 (42%)1 (3%)47 (0.76)15 (0.24)ns
CCTCTTCT
3″UTR C/T 1:5065955915 (48%)13 (42%)3 (10%)43 (0.69)19 (0.31)ns
HWE—Hardy-Weinberg equilibrium; ns—not significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Olech, M.; Ropka-Molik, K.; Szmatoła, T.; Piórkowska, K.; Kuźmak, J. Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed. Animals 2021, 11, 1908. https://doi.org/10.3390/ani11071908

AMA Style

Olech M, Ropka-Molik K, Szmatoła T, Piórkowska K, Kuźmak J. Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed. Animals. 2021; 11(7):1908. https://doi.org/10.3390/ani11071908

Chicago/Turabian Style

Olech, Monika, Katarzyna Ropka-Molik, Tomasz Szmatoła, Katarzyna Piórkowska, and Jacek Kuźmak. 2021. "Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed" Animals 11, no. 7: 1908. https://doi.org/10.3390/ani11071908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop