Estimating genomic relationships of metafounders across and within breeds using maximum likelihood, pseudo-expectation-maximization maximum likelihood and increase of relationships.

BACKGROUND: The theory of "metafounders" proposes a unified framework for relationships across base populations within breeds (e.g. unknown parent groups), and base populations across breeds (crosses) together with a sensible compatibility with genomic relationships. Considering metafounders might be advantageous in pedigree best linear unbiased prediction (BLUP) or single-step genomic BLUP. Existing methods to estimate relationships across metafounders Γ are not well adapted to highly unbalanced data, genotyped individuals far from base populations, or many unknown parent groups (within breed per year of birth). METHODS: We derive likelihood methods to estimate Γ . For a single metafounder, summary statistics of pedigree and genomic relationships allow deriving a cubic equation with the real root being the maximum likelihood (ML) estimate of Γ . This equation is tested with Lacaune sheep data. For several metafounders, we split the first derivative of the complete likelihood in a term related to Γ , and a second term related to Mendelian sampling variances. Approximating the first derivative by its first term results in a pseudo-EM algorithm that iteratively updates the estimate of Γ by the corresponding block of the H-matrix. The method extends to complex situations with groups defined by year of birth, modelling the increase of Γ using estimates of the rate of increase of inbreeding ( Δ F ), resulting in an expanded Γ and in a pseudo-EM+ Δ F algorithm. We compare these methods with the generalized least squares (GLS) method using simulated data: complex crosses of two breeds in equal or unsymmetrical proportions; and in two breeds, with 10 groups per year of birth within breed. We simulate genotyping in all generations or in the last ones. RESULTS: For a single metafounder, the ML estimates of the Lacaune data corresponded to the maximum. For simulated data, when genotypes were spread across all generations, both GLS and pseudo-EM(+ Δ F ) methods were accurate. With genotypes only available in the most recent generations, the GLS method was biased, whereas the pseudo-EM(+ Δ F ) approach yielded more accurate and unbiased estimates. CONCLUSIONS: We derived ML, pseudo-EM and pseudo-EM+ Δ F methods to estimate Γ in many realistic settings. Estimates are accurate in real and simulated data and have a low computational cost.

High-throughput DNA sequence analysis elucidates novel insight into the genetic basis of adaptation in local sheep.

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries).

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Asikli Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Asikli Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Asikli Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Polymorphism of growth hormone (GH) gene and its association with performance and body conformation of Harnali sheep.

The present study was carried out to study the polymorphism in the GH gene and its association with various performance and body conformation traits, viz., birth weight (B-WT), weaning weight (W-WT), six-month body weight (6 M-WT), one-year body weight (Y-WT), annual greasy fleece weight (AGFW), body length (BL), body height (BH), heart girth (HG) and paunch girth (PG) in 138 Harnali sheep. PCR-RFLP was performed to identify polymorphism in the targeted region of the GH gene. The PCR product of 422 bp size of the GH gene was amplified encompassing partial exon 2 and inton 3 in Harnali sheep. The PCR product was digested with HaeIII restriction enzyme for the detection of Single nucleotide polymorphism (SNP). The digested products revealed the presence of two genotypes, i.e. AA and AB in the studied population. A > G mutation (A781G) was observed in our resource population. The AA genotype was found to be the predominant genotype (0.62). Chi square value revealed that resource population was not under Hardy-Weinberg equilibrium with respect to target locus. Period of birth was found to have significant effect on W-WT, Y-WT, BL, BH and PG. Sex of animal was found to have significant (P < 0.05) effect on W-WT and highly significant (P < 0.01) effect on 6 M-WT, Y-WT and AGFW in Harnali sheep. The effect of genotype was found to be significant (P < 0.05) on annual greasy fleece weight. AB genotype was found to be associated with higher annual greasy fleece weight and can be used as a potential candidate marker in selection criteria for improving greasy fleece weight in Harnali sheep.

Whole-genome resequencing of Chinese indigenous sheep provides insight into the genetic basis underlying climate adaptation.

BACKGROUND: Chinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations. RESULTS: Based on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability. CONCLUSIONS: Analysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

Identification of the FGB gene polymorphism and analysis of its association with fat deposition traits in Hu sheep.

As a crucial economic trait, fat deposition is directly related to carcass quality and feed efficiency in sheep. The purpose of this study was to investigate the polymorphisms of the FGB gene related to fat deposition and detect the expression features of the FGB gene in different adipose tissues of sheep by using Sanger sequencing, MassARRAY® SNP technique, and quantitative real-time PCR (qRT-PCR). Results showed that in the intron region of the FGB gene, a SNP g. 3378953 A > T has been identified, and significant association was found between perirenal fat weight, perirenal fat relative weight, mesenteric fat weight, and mesenteric fat relative weight (P < 0.05). Moreover, qRT-PCR analysis showed that FGB was expressed in all three adipose tissues, and FGB gene expression level in the AA genotype was significantly lower than that in the AT or TT genotypes (P < 0.05). Therefore, the FGB gene can be used as a candidate gene to reduce fat deposition in Hu sheep breeding, and the selection of the AA genotype in Hu sheep in production practice is more conducive to improving production efficiency.

Complete Mitochondrial Genome for Lucilia cuprina dorsalis (Diptera: Calliphoridae) from the Northern Territory, Australia.

The Australian sheep blowfly, Lucilia cuprina dorsalis, is a major sheep ectoparasite causing subcutaneous myiasis (flystrike), which can lead to reduced livestock productivity and, in severe instances, death of the affected animals. It is also a primary colonizer of carrion, an efficient pollinator, and used in maggot debridement therapy and forensic investigations. In this study, we report the complete mitochondrial (mt) genome of L. c. dorsalis from the Northern Territory (NT), Australia, where sheep are prohibited animals, unlike the rest of Australia. The mt genome is 15,943 bp in length, comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and a non-coding control region. The gene order of the current mt genome is consistent with the previously published L. cuprina mt genomes. Nucleotide composition revealed an AT bias, accounting for 77.5% of total mt genome nucleotides. Phylogenetic analyses of 56 species/taxa of dipterans indicated that L. c. dorsalis and L. sericata are the closest among all sibling species of the genus Lucilia, which helps to explain species evolution within the family Luciliinae. This study provides the first complete mt genome sequence for L. c. dorsalis derived from the NT, Australia to facilitate species identification and the examination of the evolutionary history of these blowflies.

Genetic variant of the sheep E2F8 gene and its associations with litter size.

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Assessing the performance of Moghani crossbred lambs derived from different mating systems with Texel and Booroola sheep.

In our ongoing project, which focuses on the introgression of Booroola/FecB gene and the myostatin (MSTN) gene into purebred Moghani sheep, we assessed the performance of second-generation Moghani crossbreds such as second crossbreds (F2) and initial backcross generation (BC1). These crossbreds were generated through different mating systems, including in-breeding, outcrossing, first paternal backcrossing (PBC1), and first maternal backcrossing (MBC1). Notably, F2 strains exhibited lean tail, woolly fleece and a higher percentage of white coat color compared to BC1. The impact of mating systems and birth types on pre-weaning survival rates was found to be statistically significant (P < 0.0001), with singleton offspring resulting from paternal backcross showing a particularly substantial effect. The F2 crossbred lambs carrying the Booroola gene did not show a statistically significant difference in survivability compared to those carrying the MSTN gene, implying the Booroola prolificacy gene had no significant impact on survival outcomes. However, the occurrence of multiple births had a significant negative impact on lamb survival (P < 0.0001). The PBC1 sheep strains, specifically Texel Tamlet ram strains carrying the MSTN mutation, exhibited superior growth rates compared to others (P < 0.05). Interestingly, the MSTN mutation in the homozygous variant genotype significantly impacts growth rate before weaning compared to other genotypes and pure Moghani sheep (P < 0.05). In conclusion, this study objectively underscores the pivotal role of genetic factors, specifically through strategic mating systems like paternal backcrossing, in enhancing desired traits and growth rates in Moghani sheep, thereby contributing valuable insights to the field of sheep breeding programs.

A working model for the formation of Robertsonian chromosomes.

Robertsonian chromosomes form by fusion of two chromosomes that have centromeres located near their ends, known as acrocentric or telocentric chromosomes. This fusion creates a new metacentric chromosome and is a major mechanism of karyotype evolution and speciation. Robertsonian chromosomes are common in nature and were first described in grasshoppers by the zoologist W. R. B. Robertson more than 100 years ago. They have since been observed in many species, including catfish, sheep, butterflies, bats, bovids, rodents and humans, and are the most common chromosomal change in mammals. Robertsonian translocations are particularly rampant in the house mouse, Mus musculus domesticus, where they exhibit meiotic drive and create reproductive isolation. Recent progress has been made in understanding how Robertsonian chromosomes form in the human genome, highlighting some of the fundamental principles of how and why these types of fusion events occur so frequently. Consequences of these fusions include infertility and Down's syndrome. In this Hypothesis, I postulate that the conditions that allow these fusions to form are threefold: (1) sequence homology on non-homologous chromosomes, often in the form of repetitive DNA; (2) recombination initiation during meiosis; and (3) physical proximity of the homologous sequences in three-dimensional space. This Hypothesis highlights the latest progress in understanding human Robertsonian translocations within the context of the broader literature on Robertsonian chromosomes.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network.

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Network visualization of genes involved in skeletal muscle myogenesis in livestock animals.

BACKGROUND: Muscle growth post-birth relies on muscle fiber number and size. Myofibre number, metabolic and contractile capacities are established pre-birth during prenatal myogenesis. The aim of this study was to identify genes involved in skeletal muscle development in cattle, sheep, and pigs - livestock. RESULTS: The cattle analysis showed significant differences in 5043 genes during the 135-280 dpc period. In sheep, 444 genes differed significantly during the 70-120 dpc period. Pigs had 905 significantly different genes for the 63-91 dpc period.The biological processes and KEGG pathway enrichment results in each species individually indicated that DEGs in cattle were significantly enriched in regulation of cell proliferation, cell division, focal adhesion, ECM-receptor interaction, and signaling pathways (PI3K-Akt, PPAR, MAPK, AMPK, Ras, Rap1); in sheep - positive regulation of fibroblast proliferation, negative regulation of endothelial cell proliferation, focal adhesion, ECM-receptor interaction, insulin resistance, and signaling pathways (PI3K-Akt, HIF-1, prolactin, Rap1, PPAR); in pigs - regulation of striated muscle tissue development, collagen fibril organization, positive regulation of insulin secretion, focal adhesion, ECM-receptor interaction, and signaling pathways (PPAR, FoxO, HIF-1, AMPK). Among the DEGs common for studied animal species, 45 common genes were identified. Based on these, a protein-protein interaction network was created and three significant modules critical for skeletal muscle myogenesis were found, with the most significant module A containing four recognized hub genes - EGFR, VEGFA, CDH1, and CAV1. Using the miRWALK and TF2DNA databases, miRNAs (bta-miR-2374 and bta-miR-744) and transcription factors (CEBPB, KLF15, RELA, ZNF143, ZBTB48, and REST) associated with hub genes were detected. Analysis of GO term and KEGG pathways showed that such processes are related to myogenesis and associated with module A: positive regulation of MAP kinase activity, vascular endothelial growth factor receptor, insulin-like growth factor binding, focal adhesion, and signaling pathways (PI3K-Akt, HIF-1, Rap1, Ras, MAPK). CONCLUSIONS: The identified genes, common to the prenatal developmental period of skeletal muscle in livestock, are critical for later muscle development, including its growth by hypertrophy. They regulate valuable economic characteristics. Enhancing and breeding animals according to the recognized genes seems essential for breeders to achieve superior gains in high-quality muscle mass.

First detection of Ehrlichia chaffeensis, Ehrlichia canis, and Anaplasma ovis in Rhipicephalus bursa ticks collected from sheep, Turkey.

Anaplasmosis and ehrlichiosis are important tick-borne rickettsial diseases of medical and veterinary importance that cause economic losses in livestock. In this study, the prevalence of Anaplasma ovis, Ehrlichia canis and Ehrlichia chaffeensis was investigated in ticks collected from sheep in various farms in Van province, which is located in the Eastern Anatolian Region of Turkey. The ticks used in this study were collected by random sampling in 26 family farm business in 13 districts of Van province. A total of 688 ticks were collected from 88 sheep and 88 tick pools were created. All ticks identified morphologically as Rhipicephalus bursa. Phylogenetic analysis of Chaperonin and 16S rRNA gene sequences confirmed A. ovis, E. canis and E. chaffeensis in this study. Of the 88 tick pools tested, 28.41% (25/88) were positive for at least one pathogen. Anaplasma DNA was detected in five of the 88 pools (5.68%), E. canis DNA was detected in 19 of the 88 pools (21.59%), and E. chaffeensis DNA was detected in one of the 88 pools (1.14%) of R. bursa ticks. To our knowledge, this is the first report describing the presence of A. ovis, E. canis, and E. chaffeensis in R. bursa ticks collected from sheep in Turkey. Further studies are needed to investigate other co-infections in sheep in Turkey.

MiR-23b and miR-133 Cotarget TGFß2/NOTCH1 in Sheep Dermal Fibroblasts, Affecting Hair Follicle Development.

Wool is produced and controlled by hair follicles (HFs). However, little is known about the mechanisms involved in HF development and regulation. Sheep dermal fibroblasts (SDFs) play a key role in the initial stage of HF development. Analyzing the molecular mechanism that regulates early HF development in superfine wool sheep is of great importance for better understanding the HF morphogenesis process and for the breeding of fine wool sheep. Here, we show that two microRNAs (miRNAs) affect the development of HFs by targeting two genes that are expressed by SDFs. Meanwhile, the overexpression and inhibition of oar-miR-23b and oar-miR-133 in SDFs cells and cell proliferation, apoptosis, and migration were further detected using a CCK-8 assay, an Annexin V-FITC assay, a Transwell assay, and flow cytometry. We found that oar-miR-23b, oar-miR-133, and their cotarget genes TGFß2 and NOTCH1 were differentially expressed during the six stages of HF development in superfine wool sheep. Oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs and promoted the apoptosis of SDFs through TGFß2 and NOTCH1. oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs by jointly targeting TGFß2 and NOTCH1, thereby inhibiting the development of superfine wool HFs. Our research provides a molecular marker that can be used to guide the breeding of ultrafine wool sheep.

Structural and sequence analysis of the RPO30 gene of sheeppox and goatpox viruses from India.

Sheeppox and goatpox are transboundary viral diseases of sheep and goats that cause significant economic losses to small and marginal farmers worldwide, including India. Members of the genus Capripoxvirus (CaPV), namely Sheeppox virus (SPPV), Goatpox virus (GTPV), and Lumpy skin disease virus (LSDV), are antigenically similar, and species differentiation can only be accomplished using molecular approaches. The present study aimed to understand the molecular epidemiology and host specificity of SPPV and GTPV circulating in India through sequencing and structural analysis of the RNA polymerase subunit-30 kDa (RPO30) gene. A total of 29 field isolates from sheep (n = 19) and goats (n = 10) belonging to different geographical regions of India during the period: Year 2015 to 2023, were analyzed based on the sequence and structure of the full-length RPO30 gene/protein. Phylogenetically, all the CaPV isolates were separated into three major clusters: SPPV, GTPV, and LSDV. Multiple sequence alignment revealed a highly conserved RPO30 gene, with a stretch of 21 nucleotide deletion in all SPPV isolates. Additionally, the RPO30 gene of the Indian SPPV and GTPV isolates possessed several species-specific conserved signature residues/motifs that could act as genotyping markers. Secondary structure analysis of the RPO30 protein showed four α-helices, two loops, and three turns, similar to that of the E4L protein of vaccinia virus (VACV). All the isolates in the present study exhibited host preferences across different states of India. Therefore, in order to protect vulnerable small ruminants from poxviral infections, it is recommended to take into consideration a homologous vaccination strategy.

Genotyping of ticks: first molecular report of Hyalomma asiaticum and molecular detection of tick-borne bacteria in ticks and blood from Khyber Pakhtunkhwa, Pakistan.

Multiple ticks (Acari: Ixodoidea) carrying Rickettsiales bacteria have significant importance for both human and animal health. Thus, the purpose of this work was to genetically analyze tick species and their associated Rickettsiales bacteria in animal hosts. In order to achieve these objectives, various animals (including camels, cattle, goats, sheep, dogs, and mice) were inspected in four districts (Mardan, Peshawar, Kohat, and Karak) of Khyber Pakhtunkhwa to collect ticks, while blood samples were collected from all the symptomatic and asymptomatic cattle in all four districts. A total of 234 ticks were obtained from 86 out of 143 (60.14%) host animals, which were morphologically identified as Rhipicephalus turanicus, Rhipicephalus microplus, Haemaphysalis cornupunctata, and Hyalomma asiaticum. Among these, their representative ticks (126/234, 53.85%) were processed for molecular confirmation using cytochrome c oxidase (cox1) gene. Obtained cox1 sequences of four different tick species showed 99.72%-100% maximum identity with their corresponding species reported from Pakistan, China, India, and Kazakhstan and clustered phylogenetically. This study presented the first genetic report of Hy. asiaticum ticks in Pakistan. Moreover, genetically confirmed tick species were molecularly analyzed by PCR for detection of Rickettsiales DNA using partial fragments of 16S rDNA, 190-kDa outer membrane protein A (ompA), and 120-kDa outer membrane protein B (ompB) genes. In addition, blood samples were analyzed to identify Rickettsiales bacteria using the aforementioned genes. Rickettsiales bacteria were found in 24/126 (19.05%) ticks and 4/16 (25.00%) in symptomatic cattle's blood. The obtained ompA and ompB sequences from Hy. asiaticum ticks showed 99.73%-99.87% with Candidatus Rickettsia shennongii and unidentified Rickettsia sp., whereas the obtained 16S rDNA sequences from cattle's blood and ticks (Hae. cornupunctata) showed 99.67% highest identity with Anaplasma phagocytophilum. The 16S rDNA sequence of Rickettsiales DNA from Rh. turanicus ticks showed 100% identity with Ehrlichia canis and unidentified Ehrlichia sp. Obtained sequences of Rickettsiales bacteria were grouped along with their respective species in phylogenetic trees, which were previously reported in Greece, Cuba, Iraq, Turkey, Pakistan, South Korea, and China (mainland and Taiwan). This extensive study explores the wide range of damaging ticks and their corresponding tick-borne bacteria in the area, suggesting a possible danger to both livestock and human communities.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn.

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

Fitting of Growth Curves and Estimation of Genetic Relationship between Growth Parameters of Qianhua Mutton Merino.

Qianhua Mutton Merino is a dual-purpose (meat and wool) breed of sheep that has been newly developed in China. In this study, we assessed the growth and development of the Qianhua Mutton Merino sheep breed under house feeding conditions by measuring the body weight and chest circumference of 2300 rams and ewes of this breed aged 0-24 months. Based on the fitting results of three nonlinear growth models, namely Logistic, Gompertz, and von Bertalanffy, in Qianhua Mutton Merino, we selected the von Bertalanffy model because of its highest fitting degree among all models (R2 > 0.977). The significant analysis of the combined fixation of each sheep body's weight and bust took place (A: mature body weight, B: adjustment parameter, K: instant relative growth rate). The results revealed that parameters A, B, and K of body weight and chest circumference have high heritability and thus could be used as target traits for genetic improvement. Moreover, the correlation strength among A, B, and K suggested that these parameters can be used as a reference to adjust the genetic parameters in the growth model to genetically improve the body size of Qianhua Mutton Merino during breeding.