Search

Abstract Background Indigenous buffaloes, as the important livestock species contributing to economy of the country, are the lifeline of livelihood in Nepal. They are distributed across diverse geographical regions of the country and have adapted to various feeding, breeding, and management conditions. The larger group of these native buffalo breeds are present in narrow and stiff hilly terrains. Their dispersal indicates a possible environmental adaptation mechanism, which is crucial for the conservation of these breeds. Results We utilized whole-genome sequencing (WGS) to investigate the genetic diversity, population structure, and selection signatures of Nepalese indigenous buffaloes. We compared 66 whole-genome sequences with 118 publicly available sequences from six river and five swamp buffalo breeds. Genomic diversity parameters indicated genetic variability level in the Nepalese buffaloes comparable to those of Indian breeds, and population genetic structure revealed distinct geography-mediated genetic differentiation among these breeds. We used locus-specific branch length analysis (LSBL) for genome-wide scan, which revealed a list of potentially selected genes in Lime and Parkote breeds that inhabit the hilly region. A gene ontology (GO) analysis discovered that many GO terms were associated with cardiac function regulation. Furthermore, complementary analyses of local selection signatures, tissue expression profiles, and haplotype differences identified candidate genes, including KCNE1, CSF1R, and PDGFRB, related to the regulation of cardiac and pulmonary functions. Conclusions This study is a comprehensive WGS-based genetic analysis of the native Nepalese buffalo breeds. Our study suggested that the Nepalese “hilly” buffaloes, especially the Lime and Parkote breeds, have undergone some characteristic genetic changes and evolved increased cardiac and pulmonary function for their adaptation to the steep hilly terrains of the country.

Abstract To identify the genetic determinants of domestication and productivity of Asian water buffaloes (Bubalus bubalis), 470 genomes of domesticated river and swamp buffaloes along with their putative ancestors, the wild water buffaloes (Bubalus arnee) are sequenced and integrated. The swamp buffaloes inherit the morphology of the wild buffaloes. In contrast, most river buffaloes are unique in their morphology, but their genomes cluster with the wild buffaloes. The levels of genomic diversity in Italian river and Indonesian swamp buffaloes decrease at opposite extremes of their distribution range. Purifying selection prevented the accumulation of harmful loss‐of‐function variants in the Indonesian buffaloes. Genes that evolved rapidly (e.g., GKAP1) following differential selections in the river and swamp buffaloes are involved in their reproduction. Genes related to milk production (e.g., CSN2) and coat color (e.g., MC1R) underwent strong selections in the dairy river buffaloes via soft and hard selective sweeps, respectively. The selective sweeps and single‐cell RNA‐seq data revealed the luminal cells as the key cell type in response to artificial selection for milk production of the dairy buffaloes. These findings show how artificial selection has been driving the evolutionary divergence and genetic differentiation in morphology and productivity of Asian water buffaloes.

Weaning weight is a key indicator of the early growth performance of cattle. An understanding of the genetic mechanisms underlying weaning weight will help increase the accuracy of selection of breeding animals. In order to identify candidate genes associated with weaning weight in Simmental-Holstein crossbred cattle, this study generated RNA-Sequencing (RNA-seq) data from 86 crossbred calves (37 males and 49 famales) and measured their weaning weight and body size traits (wither height, body length, chest girth, rump width, and rump length). Differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) were performed. A total of 498 differentially expressed genes (DEGs) were identified between the low weaning weight (LWW) group and the high weaning weight (HWW) group. Weaning weight was transcriptionally correlated (FDR < 0.05) with four of the eleven co-expression gene modules. By intersecting DEGs and hub genes of the four modules, we identified a final set of 37 candidate genes enriched in growth, development, or immune-related processes. In addition, one co-expression module was significantly correlated with all the five body size traits (P

Milk yield and its composition show individual variation due to the effects of the environment. Previous studies suggest that meteorological variables exert negative effects on milk yield and composition, especially during summer. This study aimed to examine the effects of meteorological variables on bulk milk composition in the Sardinian sheep production system. In this work, a total of 218,170 records belonging to 4562 dairy sheep farms were merged with the meteorological data provided by 60 meteorological stations located on Sardinia Island (Italy). Milk composition in the late spring and summer recorded during a 5-year period was used to evaluate the impact of climate exposure on bulk milk traits. The milk quality was analyzed using a linear mixed model that included the fixed effects of the year of sampling, the flock size, the temperature humidity index (THI) and the random effect of the flock. The variability of milk composition explained by flock and management ranged from 30 to 64%. The flock size exerted a significant effect on milk composition: large flocks characterized by advanced management and feeding techniques resulted in higher milk quality (e.g., higher protein and fat, lower lactose) compared to traditionally managed small flocks. The impact of THI on milk composition was statistically significant across different milk quality traits (p < 0.001); the effect of thermal stress varied according to the month of lactation. For instance, milk fat content in May increased by +0.4% for THI > 76. In June, no relevant differences were observed, whereas a decrease in fat percentage was observed in July as THI values increased (up to −0.5% for THI > 76). While somatic cell counts remained relatively stable across different conditions, total bacterial count showed greater seasonal variability, peaking during warmer periods. In addition, using factor analysis, we developed a multivariate meteorological index (MMI), which explained 51% of the variance of the original meteorological data. MMI was highly correlated with THI (r = 0.75). The same linear mixed model applied for modeling THI was used to assess the effect of MMI on milk traits. Fat, protein fractions and lactose showed significant variation across MMI classes (p-value < 0.001) in the same direction as those based on THI. Overall, our findings underscore the impact of both flock size and environmental conditions on milk quality, with heat stress and traditional versus modern management practices leading to measurable differences in milk traits.

Background: Hepatitis delta virus (HDV) infection represents the most severe form of viral hepatitis and is a significant global health challenge. Bulevirtide (BLV) is a novel therapeutic treatment that has resulted in variable response rates in HBV/HDV-coinfected patients. We evaluated clinical, virological, and polymorphic factors for the purpose of predicting BLV treatment success. Methods: Thirty HBV/HDV-coinfected patients received BLV monotherapy (2 mg/day) for 24 to 48 weeks. Baseline (BL) serum samples were collected to assess clinical parameters and virological markers (HDV RNA, HBV DNA, HBsAg, HBcrAg, anti-HBc IgG) at treatment weeks 24 (TW24) and 48 (TW48). Additionally, full-genome HDV sequencing and a phylogenetic analysis were performed. Finally, analyses of the HDAg protein sequence and HDV RNA secondary structure were conducted to evaluate potential associations with treatment response. Results: A significant reduction in HDV RNA levels was observed at TW48, with a virological response (HDV RNA undetectable or ≥2 Log decline from BL) achieved by 58% of patients. Median BL levels of anti-HBc IgG were significantly different between virological responders (39.3 COI; interquartile range [IQR] 31.6–47.1) and virological non-responders (244.7 COI; IQR 127.0–299.4) (p = 0.0001). HDV genotype 1e was predominant across the cohort, and no specific HDAg polymorphisms predicted the response. However, secondary structure analysis of HDV RNA revealed that a specific pattern of internal loops in the region 63–100 nucleotides downstream of the editing site may influence treatment response by impacting editing efficacy. Conclusions: This study revealed key factors influencing BLV efficacy in HBV/HDV coinfection. Lower baseline anti-HBc IgG levels strongly correlated with a positive virological response, suggesting that the liver’s inflammatory state affects treatment success. Additionally, the analysis of HDV RNA secondary structure in patients receiving BLV treatment revealed a higher editing efficiency in virological responders, highlighting areas for further research.

Livestock expresses complex traits influenced by several factors. The response of animals to variations in climatic factors, such as increases in temperature, may induce heat stress conditions. In this study, animals living at different temperatures were compared using the genome-wide Wright fixation index (FST). A total of 825 genotypes of Sarda breed ewes were divided into two groups based on the flocks’ average temperature over a 20-year period to compute the FST: 395 and 430 sheep were represented in colder and hotter groups, respectively. After LOWESS regression and CONTROL CHART application, 623 significant markers and 97 selection signatures were found. A total of 280 positional candidate genes were retrieved from a public database. Among these genomic regions, we found 51 annotated genes previously associated with heat stress/tolerance in ruminants (FCGR1A, MDH1, UGP2, MYO1G, and HSPB3), as well as immune response and cellular mechanisms related to how animals cope with thermal stress (RIPK1, SERPINB1, SERPINB9, and PELI1). Moreover, other genes were associated with milk fat (SCD, HERC3, SCFD2, and CHUK), body weight, body fat, and intramuscular fat composition (AGPAT2, ABCD2, MFAP32, YTHDC1, SIRT3, SCD, and RNF121), which might suggest the influence of environmental conditions on the genome of Sarda sheep.

Persistence is a strategy used by many viruses to evade eradication by the immune system, ensuring their permanence and transmission within the host and optimizing viral fitness. During persistence, viruses can trigger various phenomena, including target organ damage, mainly due to an inflammatory state induced by infection, as well as cell proliferation and/or immortalization. In addition to immune evasion and chronic inflammation, factors contributing to viral persistence include low-level viral replication, the accumulation of viral mutants, and, most importantly, maintenance of the viral genome and reliance on viral oncoprotein production. This review focuses on the process of genome integration, which may occur at different stages of infection (e.g., HBV), during the chronic phase of infection (e.g., HPV, EBV), or as an essential part of the viral life cycle, as seen in retroviruses (HIV, HTLV-1). It also explores the close relationship between integration, persistence, and oncogenesis. Several models have been proposed to describe the genome integration process, including non-homologous recombination, looping, and microhomology models. Integration can occur either randomly or at specific genomic sites, often leading to genome destabilization. In some cases, integration results in the loss of genomic regions or impairs the regulation of oncogene and/or oncosuppressor expression, contributing to tumor development.

Viral spillover represents the transmission of pathogen viruses from one species to another that can give rise to an outbreak. It is a critical concept that has gained increasing attention, particularly after the SARS-CoV-2 pandemic. However, the term is often used inaccurately to describe events that do not meet the true definition of spillover. This review aims to clarify the proper use of the term and provides a detailed analysis of the mechanisms driving zoonotic spillover, with a focus on the genetic and environmental factors that enable viruses to adapt to new hosts. Key topics include viral genetic variability in reservoir species, biological barriers to cross-species transmission, and the factors that influence viral adaptation and spread in novel hosts. The review also examines the role of evolutionary processes such as mutation and epistasis, alongside ecological conditions that facilitate the emergence of new pathogens. Ultimately, it underscores the need for more accurate predictive models and improved surveillance to better anticipate and mitigate future spillover events.

The αs2-casein is a phosphoprotein secreted in the milk of most mammals, and it is the most hydrophilic of all caseins. Contrary to genes found in ruminants, in donkeys two different encoding genes for donkey αs2-casein (CSN1S2 I and CSN1S2 II) have been identified. However, unlike in ruminants, the variability at these loci has not been characterized in detail in donkeys until now. In this study, we analyze the transcript profile of the donkey CSN1S2 I and CSN1S2 II genes, and we identify and describe the variability of these loci in the Ragusana and Amiatina breeds reared in Italy. The analysis of the CSN1S2 I Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) products and subsequent sequencing showed, in addition to correctly spliced mRNA, seven other minor mRNAs resulting from differential splicing events involving, in various combinations, entire exons (4, 5, 6, and 11), parts of exons (5′ or 3′ end of exon 17), or the recognition of intronic sequences as an exon (exon 12′). Similarly, the transcription analysis of the CSN1S2 II gene revealed a remarkable variability in splicing events, mainly concerning the alternative insertion of an extra exon 7 (named 7′); the first 33 bp of exon 13; or the alternative skipping of exons 9, 10, 11, 12, and 15, and their combinations. At the mRNA level for CSN1S2 I, seven SNPs were observed, five of which led to amino acid changes: p.T73>A, p.I109>V, p.I130>V, p.I146>T, and p.D217>Y. Similarly, nine SNPs were observed at the CSN1S2 II locus, seven of which are non-synonymous: p.L63>F, p.H70>Q, p.D90>N, p.129A>T, p.H131>Y, p.E144>G, and p.F157>S. In addition, the DNA sequencing of exon 17 and flanking introns of the CSN1S2 I gene revealed a G>A transition at the splice acceptor site of CSN1S2 I exon 17 (FM946022.1:c.375-1G>A), resulting in an allele-specific skipping of the first 15 nucleotides of this exon, which encode the peptide 176NKINQ180, and the recognition of an in-frame cryptic splicing acceptor site: arAACAAAATCAACCAG. A genotyping method based on restriction fragment length polymorphism (XbaI PCR-RFLP) was set up for this SNP. In the total population studied (105 Ragusana and 14 Amiatina donkeys), the A allele had a frequency of 0.2437 with no evidence of deviation from the Hardy–Weinberg equilibrium. This study adds new knowledge regarding the genetic variability of αs2-caseins in donkeys and may contribute significantly to the genetic improvement of milk production for this species.

Background: Human papillomavirus (HPV) infection is linked to several cancers, including anal and oral cancers. The incidence of anal cancer is particularly high among HIV-positive men who have sex with men (MSM). DNA methylation markers have shown promise as biomarkers for identifying precancerous lesions and cancer in HPV-infected individuals. The aim of this study was to investigate the correlation of DNA methylation with HPV infection in oral samples and the correlation of DNA methylation with lesion degree in the anal samples of HIV-positive MSM. Methods: This study investigated DNA methylation in oral and anal samples from HIV-positive MSM at the National Institute for Infectious Diseases (INMI) in Rome, Italy. Exfoliated oral epithelial cells and anal samples were collected and analyzed for 28 HPV genotypes using the Allplex 28 HPV assay. DNA methylation was assessed with the PrecursorM+ kit for oral samples and the AnoGyn kit for anal samples, focusing on the promoter regions of specific genes. Results: The study included 63 participants, with a median age of 49 and a median CD4+ count of 705 cells/µL. The oral samples showed HPV16 as the most common type, with 22% testing positive for DNA methylation. The anal samples exhibited HPV-related methylation changes linked to cytological lesions, with a 30% increase in the observed ddCt ratio. Significant differences were found in both ASCL1 and ZNF582 genes, particularly for HSILvsNILM and HSILvsLSIL lesions. Of the samples with an increased ddCt ratio, 80% were from patients over 35 years old, and multiple HPV infections were common. Conclusions: DNA methylation markers could be valuable in identifying high-risk HPV infections in oral samples and detecting potential precancerous lesions in anal samples. These markers may enhance the early detection and prevention strategies for HPV-related cancers in high-risk populations, with follow-up data indicating potential for monitoring lesion progression.

Therapy against the Hepatitis C virus (HCV) has significantly improved with the introduction of direct-acting antiviral drugs (DAAs), achieving over 95% sustained virological response (SVR). Despite this, the development of an effective anti-HCV vaccine remains a critical challenge due to the low number of patients treated with DAAs and the occurrence of HCV reinfections in high-risk groups. Current vaccine strategies aim to stimulate either B-cell or T-cell responses. Vaccines based on E1 and E2 proteins can elicit broad cross-neutralizing antibodies against all major HCV genotypes, though with varying efficiencies and without full protection against infection. In humans, the neutralizing antibodies induced by such vaccines mainly target the AR3 region, but their levels are generally insufficient for broad neutralization. Various HCV proteins expressed through different viral vectors have been utilized to elicit T cell immune responses, showing sustained expansion of HCV-specific effector memory T cells and improved proliferation and polyfunctionality of memory T cells over time. However, despite these advancements, the frequency and effectiveness of T-cell responses remain limited.

Abstract Cattle and the draught force provided by its skeletal muscle have been integral to agro‐ecosystems of agricultural civilization for millennia. However, relatively little is known about the cattle muscle functional genomics (including protein coding genes, non‐coding RNA, etc.). Circular RNAs (circRNAs), as a new class of non‐coding RNAs, can be effectively translated into detectable peptides, which enlightened us on the importance of circRNAs in cattle muscle physiology function. Here, RNA‐seq, Ribosome profiling (Ribo‐seq), and peptidome data are integrated from cattle skeletal muscle, and detected five encoded peptides from circRNAs. It is further identified and functionally characterize a 907‐amino acids muscle‐specific peptide that is named circNEB‐peptide because derived by the splicing of Nebulin (NEB) gene. This peptide localizes to the nucleus and cytoplasm and directly interacts with SKP1 and TPM1, key factors regulating physiological activities of myoblasts, via ubiquitination and myoblast fusion, respectively. The circNEB‐peptide is found to promote myoblasts proliferation and differentiation in vitro, and induce muscle regeneration in vivo. These findings suggest circNEB‐peptide is an important regulator of skeletal muscle regeneration and underscore the possibility that more encoding polypeptides derived by RNAs currently annotated as non‐coding exist.

ABSTRACT: β-lactoglobulin I (β-LG I) is one of the most important whey proteins in donkey milk. However, to our knowledge, there has been no study focusing on the full nucleotide sequences of this gene (BLG I). Current investigation of donkey BLG I gene is very limited with only 2 variants (A and B) characterized so far at the protein level. Recently, a new β-LG I variant, with a significantly higher mass (+1,915 Da) than known variants has been detected. In this study, we report the whole nucleotide sequence of the BLG I gene from 2 donkeys, whose milk samples are characterized by the β-LG I SDS-PAGE band with a normal electrophoretic mobility (18,514.25 Da, β-LG I B1 form) the first, and by the presence of a unique β-LG I band with a higher electrophoretic mobility (20,428.5 Da, β-LG I D form) the latter. A high genetic variability was found all over the 2 sequenced BLG I alleles. In particular, 16 polymorphic sites were found in introns, one in the 5′ flanking region, 3 SNPs in the 5′ untranslated region and one SNP in the coding region (g.1871G > A) located at the 40th nucleotide of exon 2 and responsible for the AA substitutions p.Asp28 > Asn in the mature protein. Two SNPs (g.920–922CAC > TGT and g.1871G/A) were genotyped in 93 donkeys of 2 Italian breeds (60 Ragusana and 33 Amiatina, respectively) and the overall frequencies of g.920–922CAC and g.1871A were 0.3065 and 0.043, respectively. Only the rare allele g.1871A was observed to be associated with the slower migrating β-LG I. Considering this genetic diversity and those found in the database, it was possible to deduce at least 5 different alleles (BLG I A, B, B1, C, D) responsible for 4 potential β-LG I translations. Among these alleles, B1 and D are those characterized in the present research, with the D allele of real novel identification. Haplotype data analysis suggests an evolutionary pathway of donkey BLG I gene and a possible phylogenetic map is proposed. Analyses of mRNA secondary structure showed relevant changes in the structures, as consequence of the g.1871G > A polymorphism, that might be responsible for the recognition of an alternative initiation site providing an additional signal peptide. The extension of 19 AA sequence to the mature protein, corresponding to the canonical signal peptide with an additional alanine residue, is sufficient to provide the observed molecular weight of the slower migrating β-LG I encoded by the BLG I D allele.

The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes and in DGAT1. Nowadays, similar results have not yet been reached in river buffaloes, despite the availability of advanced genomic technologies and accurate phenotype records. The aim of the present study was to investigate and validate the effect of four single nucleotide polymorphisms (SNP) in the CSN1S1, CSN3, SCD and LPL genes on seven milk traits in a larger buffalo population. These SNPs have previously been reported to be associated with, or affect, dairy traits in smaller populations often belonging to one farm. A total of 800 buffaloes were genotyped. The following traits were individually recorded, monthly, throughout each whole lactation period from 2010 to 2021: daily milk yield (dMY, kg), protein yield (dPY, kg) and fat yield (dFY, kg), fat and protein contents (dFP, % and dPP, %), somatic cell count (SCC, 103 cell/mL) and urea (mg/dL). A total of 15,742 individual milk test day records (2496 lactations) were available for 680 buffalo cows, with 3.6 ± 1.7 parities (from 1 to 13) and an average of 6.1 ± 1.2 test day records per lactation. Three out four SNPs in the CSN1S1, CSN3 and LPL genes were associated with at least one of analyzed traits. In particular, the CSN1S1 (AJ005430:c.578C>T) gave favorable associations with all yield traits (dMY, p = 0.022; dPY, p = 0.014; dFY, p = 0.029) and somatic cell score (SCS, p = 0.032). The CSN3 (HQ677596: c.536C>T) was positively associated with SCS (p = 0.005) and milk urea (p = 0.04). Favorable effects on daily milk yield (dMY, p = 0.028), fat (dFP, p = 0.027) and protein (dPP, p = 0.050) percentages were observed for the LPL. Conversely, the SCD did not show any association with milk traits. This is the first example of a confirmation study carried out in the Mediterranean river buffalo for genes of economic interest in the dairy field, and it represents a very important indication for the preselection of young bulls destined for breeding programs aimed at more sustainable dairy production.

Household buffalo dairy farming is gaining popularity nowadays in Bangladesh because of the outstanding food value of buffalo milk as well as the lower production cost of buffalo compared to cattle. An initiative has recently been taken for the genetic improvement of indigenous dairy buffaloes. The present study was carried out to determine the influence of some environmental factors like age, parity, season of calving, calving interval, dry period on the lactation yield, and lactation curve of indigenous dairy buffaloes of Bangladesh. A total of 384 indigenous dairy buffaloes from the 3rd and 4th parity of seven herds under two different agroecological zones covering four seasons were selected and ear tagged for individual buffalo milk recording. A milk yield of 300 days (MY300d) was calculated following the International Committee for Animal Recording (ICAR) and the data were evaluated using the generalized linear model (GLM). In production traits, the mean of calculated lactation period (CLP), calculated lactation yield (CLY), and milk yield of 300 days (MY300d) of the overall population were 267.28 days, 749.36 kg, and 766.92 kg, respectively, whereas calving interval (CI) and dry period (DP) as reproductive traits were 453.06 days and 185.78 days, respectively. The season of calving, age of buffalo cows, population or herd, agroecological zone, calving interval, and dry period had significant effects on production traits (p < 0.05 to p < 0.001). The season of calving, level of milk production of 300 days, population, and agroecological zone significantly affected the reproduction traits (p < 0.01 to p < 0.001). Parity was found to be non-significant for both types of traits. The average peak yield of test day (TD) milk production was highest at TD4 (4.47 kg, 98th day of lactation). The average MY300d of milk production was the highest in the Lalpur buffalo population (1076.13 kg) and the lowest in the buffalo population of Bhola (592.44 kg). The correlations between milk production traits (CLP, CLY, and MY-300d) and reproduction traits (CI and DP) were highly significant (p < 0.01 to p < 0.001). Positive and high correlation was found within milk traits and reproduction traits, but correlation was negative between milk traits and reproduction traits. Therefore, these non-genetic factors should be considered in the future for any genetic improvement program for indigenous dairy buffaloes in Bangladesh.

Coronaviridae is recognized as one of the most rapidly evolving virus family as a consequence of the high genomic nucleotide substitution rates and recombination. The family comprises a large number of enveloped, positive-sense single-stranded RNA viruses, causing an array of diseases of varying severity in animals and humans. To date, seven human coronaviruses (HCoV) have been identified, namely HCoV-229E, HCoV-NL63, HCoV-OC43 and HCoV-HKU1, which are globally circulating in the human population (seasonal HCoV, sHCoV), and the highly pathogenic SARS-CoV, MERS-CoV and SARS-CoV-2. Seasonal HCoV are estimated to contribute to 15–30% of common cold cases in humans; although diseases are generally self-limiting, sHCoV can sometimes cause severe lower respiratory infections and life-threatening diseases in a subset of patients. No specific treatment is presently available for sHCoV infections. Herein we show that the anti-infective drug nitazoxanide has a potent antiviral activity against three human endemic coronaviruses, the Alpha-coronaviruses HCoV-229E and HCoV-NL63, and the Beta-coronavirus HCoV-OC43 in cell culture with IC50 ranging between 0.05 and 0.15 μg/mL and high selectivity indexes. We found that nitazoxanide does not affect HCoV adsorption, entry or uncoating, but acts at postentry level and interferes with the spike glycoprotein maturation, hampering its terminal glycosylation at an endoglycosidase H-sensitive stage. Altogether the results indicate that nitazoxanide, due to its broad-spectrum anti-coronavirus activity, may represent a readily available useful tool in the treatment of seasonal coronavirus infections.

The next generation sequencing has significantly contributed to clarify the genome structure of many species of zootechnical interest. However, to date, some portions of the genome, especially those linked to a heterogametic nature such as the Y chromosome, are difficult to assemble and many gaps are still present. It is well known that the fluorescence in situ hybridization (FISH) is an excellent tool for identifying genes unequivocably mapped on chromosomes. Therefore, FISH can contribute to the localization of unplaced genome sequences, as well as to correct assembly errors generated by comparative bioinformatics. To this end, it is necessary to have starting points; therefore, in this study, we reviewed the physically mapped genes on the Y chromosome of cattle, buffalo, sheep, goats, pigs, horses and alpacas. A total of 208 loci were currently mapped by FISH. 89 were located in the male-specific region of the Y chromosome (MSY) and 119 were identified in the pseudoautosomal region (PAR). The loci reported in MSY and PAR were respectively: 18 and 25 in Bos taurus, 5 and 7 in Bubalus bubalis, 5 and 24 in Ovis aries, 5 and 19 in Capra hircus, 10 and 16 in Sus scrofa, 46 and 18 in Equus caballus. While in Vicugna pacos only 10 loci are reported in the PAR region. The correct knowledge and assembly of all genome sequences, including those of genes mapped on the Y chromosome, will help to elucidate their biological processes, as well as to discover and exploit potentially epistasis effects useful for selection breeding programs.

ABSTRACT: Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (αS1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

Rift Valley fever is a vector-borne zoonotic disease caused by the Rift Valley fever virus (Phlebovirus genus) listed among the eight pathogens included in the Bluepoint list by the WHO. The transmission is mainly vehicled by Aedes and Culex mosquito species. Symptoms of the disease are varied and non-specific, making clinical diagnosis often challenging, especially in the early stages. Due to the difficulty in distinguishing Rift Valley fever from other viral hemorrhagic fevers, as well as many other diseases that cause fever, an early diagnosis of the infection is important to limit its spread and to provide appropriate care to patients. To date, there is no validated point-of-care diagnostic tool. The virus can only be detected in the blood for a brief period, suggesting that molecular methods alone are not sufficient for case determination. For this, it is preferable to combine both molecular and serological tests. The wide distribution of competent vectors in non-endemic areas, together with global climate change, elicit the spread of RVFV to continents other than Africa, making surveillance activities vital to prevent or to limit the impact of human outbreaks and for a rapid identification of positive cases, making diagnosis a key factor for this achievement.

Abstract Aneuploidy is one of the main causes of fetal and embryonic mortality in mammals. Nonetheless, its incidence in domestic ruminants has been investigated little. Indeed, no incidence data have ever been reported for water buffalo. To establish the incidence of aneuploidy in this species, we analysed in vitro matured metaphase II (MII) oocytes with corresponding first polar bodies (I PB) of the river (2n = 50) and swamp (2n = 48) buffaloes. For the first time, six river type probes (corresponding to chromosomes 1–5 and heterosome X), were tested on swamp buffalo metaphases using Multicolor-Fluorescent In Situ Hybridization (M-FISH) before their use on oocytes MII metaphases. Of the 120 total Cumulus Oocyte Complexes (COCs, 60 for each buffalo type) subjected to in vitro maturation, 104 reached the MII stage and were analysed by M-FISH. Haploid chromosome arrangement and visible I PB were observed in 89 of the oocytes (45 in river and 44 in swamp type). In the river type, the analysis revealed one oocyte was disomic for the chromosome X (2.22%). In the swamp type, one oocyte was found to be nullisomic for chromosome X (2.27%); another was found to be nullisomic for chromosome 5 (2.27%). We also observed one oocyte affected by a premature separation of sister chromatids (PSSC) on the chromosome X (2.27%). In both buffalo types, no abnormalities were detected in other investigated chromosomes. Based on merged data, the overall aneuploidy rate for the species was 3.37%. Oocytes with unreduced chromosomes averaged 1.92% across the two types, with 1.96% in river and 1.88% in swamp. The interspecies comparison between these data and cattle and pig published data revealed substantial difference in both total aneuploidy and diploidy rates. Reducing the negative impact of the meiotic segregation errors on the fertility is key to more sustainable breeding, an efficient embryo transfer industry and ex-situ bio-conservation. In this respect, additional M-FISH studies are needed on oocytes of domestic species using larger sets of probes and/or applying next generation sequencing technologies.

Buffalo is an important livestock in Asia. Here, the authors present a comprehensive metagenomic analysis of the microbial communities present in different sites and compartments along the buffalo digestive tract.

Nipah virus (NiV) is a paramyxovirus responsible for a high mortality rate zoonosis. As a result, it has been included in the list of Blueprint priority pathogens. Bats are the main reservoirs of the virus, and different clinical courses have been described in humans. The Bangladesh strain (NiV-B) is often associated with severe respiratory disease, whereas the Malaysian strain (NiV-M) is often associated with severe encephalitis. An early diagnosis of NiV infection is crucial to limit the outbreak and to provide appropriate care to the patient. Due to high specificity and sensitivity, qRT-PCR is currently considered to be the optimum method in acute NiV infection assessment. Nasal swabs, cerebrospinal fluid, urine, and blood are used for RT-PCR testing. N gene represents the main target used in molecular assays. Different sensitivities have been observed depending on the platform used: real-time PCR showed a sensitivity of about 103 equivalent copies/reaction, SYBRGREEN technology’s sensitivity was about 20 equivalent copies/reaction, and in multiple pathogen card arrays, the lowest limit of detection (LOD) was estimated to be 54 equivalent copies/reaction. An international standard for NiV is yet to be established, making it difficult to compare the sensitivity of the different methods. Serological assays are for the most part used in seroprevalence studies owing to their lower sensitivity in acute infection. Due to the high epidemic and pandemic potential of this virus, the diagnosis of NiV should be included in a more global One Health approach to improve surveillance and preparedness for the benefit of public health. Some steps need to be conducted in the diagnostic field in order to become more efficient in epidemic management, such as development of point-of-care (PoC) assays for the rapid diagnosis of NiV.

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3′-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels.

Simple Summary: In extensive breeding systems, environmental conditions strongly influence animal behavior and production. During the process of evolution, animals tended to adapt their morphology and physiology to environmental conditions, leaving genomic imprints. This adaptation can be traced in the animals' genomes, relating environmental features to genome-wide differentiation metrics. In this study, using Sarda sheep living at different temperatures as a case study, we compared their genomes to highlight traces of thermal tolerance and adaptation. Livestock expresses complex traits influenced by several factors. The response of animals to variations in climatic factors, such as increases in temperature, may induce heat stress conditions. In this study, animals living at different temperatures were compared using the genome-wide Wright fixation index (FST). A total of 825 genotypes of Sarda breed ewes were divided into two groups based on the flocks' average temperature over a 20-year period to compute the FST: 395 and 430 sheep were represented in colder and hotter groups, respectively. After LOWESS regression and CONTROL CHART application, 623 significant markers and 97 selection signatures were found. A total of 280 positional candidate genes were retrieved from a public database. Among these genomic regions, we found 51 annotated genes previously associated with heat stress/tolerance in ruminants (FCGR1A, MDH1, UGP2, MYO1G, and HSPB3), as well as immune response and cellular mechanisms related to how animals cope with thermal stress (RIPK1, SERPINB1, SERPINB9, and PELI1). Moreover, other genes were associated with milk fat (SCD, HERC3, SCFD2, and CHUK), body weight, body fat, and intramuscular fat composition (AGPAT2, ABCD2, MFAP32, YTHDC1, SIRT3, SCD, and RNF121), which might suggest the influence of environmental conditions on the genome of Sarda sheep. [ABSTRACT FROM AUTHOR]

Persistence is a strategy used by many viruses to evade eradication by the immune system, ensuring their permanence and transmission within the host and optimizing viral fitness. During persistence, viruses can trigger various phenomena, including target organ damage, mainly due to an inflammatory state induced by infection, as well as cell proliferation and/or immortalization. In addition to immune evasion and chronic inflammation, factors contributing to viral persistence include low-level viral replication, the accumulation of viral mutants, and, most importantly, maintenance of the viral genome and reliance on viral oncoprotein production. This review focuses on the process of genome integration, which may occur at different stages of infection (e.g., HBV), during the chronic phase of infection (e.g., HPV, EBV), or as an essential part of the viral life cycle, as seen in retroviruses (HIV, HTLV-1). It also explores the close relationship between integration, persistence, and oncogenesis. Several models have been proposed to describe the genome integration process, including non-homologous recombination, looping, and microhomology models. Integration can occur either randomly or at specific genomic sites, often leading to genome destabilization. In some cases, integration results in the loss of genomic regions or impairs the regulation of oncogene and/or oncosuppressor expression, contributing to tumor development. [ABSTRACT FROM AUTHOR]

Background: Indigenous buffaloes, as the important livestock species contributing to economy of the country, are the lifeline of livelihood in Nepal. They are distributed across diverse geographical regions of the country and have adapted to various feeding, breeding, and management conditions. The larger group of these native buffalo breeds are present in narrow and stiff hilly terrains. Their dispersal indicates a possible environmental adaptation mechanism, which is crucial for the conservation of these breeds. Results: We utilized whole-genome sequencing (WGS) to investigate the genetic diversity, population structure, and selection signatures of Nepalese indigenous buffaloes. We compared 66 whole-genome sequences with 118 publicly available sequences from six river and five swamp buffalo breeds. Genomic diversity parameters indicated genetic variability level in the Nepalese buffaloes comparable to those of Indian breeds, and population genetic structure revealed distinct geography-mediated genetic differentiation among these breeds. We used locus-specific branch length analysis (LSBL) for genome-wide scan, which revealed a list of potentially selected genes in Lime and Parkote breeds that inhabit the hilly region. A gene ontology (GO) analysis discovered that many GO terms were associated with cardiac function regulation. Furthermore, complementary analyses of local selection signatures, tissue expression profiles, and haplotype differences identified candidate genes, including KCNE1, CSF1R, and PDGFRB, related to the regulation of cardiac and pulmonary functions. Conclusions: This study is a comprehensive WGS-based genetic analysis of the native Nepalese buffalo breeds. Our study suggested that the Nepalese "hilly" buffaloes, especially the Lime and Parkote breeds, have undergone some characteristic genetic changes and evolved increased cardiac and pulmonary function for their adaptation to the steep hilly terrains of the country. [ABSTRACT FROM AUTHOR]

Viral spillover represents the transmission of pathogen viruses from one species to another that can give rise to an outbreak. It is a critical concept that has gained increasing attention, particularly after the SARS-CoV-2 pandemic. However, the term is often used inaccurately to describe events that do not meet the true definition of spillover. This review aims to clarify the proper use of the term and provides a detailed analysis of the mechanisms driving zoonotic spillover, with a focus on the genetic and environmental factors that enable viruses to adapt to new hosts. Key topics include viral genetic variability in reservoir species, biological barriers to cross-species transmission, and the factors that influence viral adaptation and spread in novel hosts. The review also examines the role of evolutionary processes such as mutation and epistasis, alongside ecological conditions that facilitate the emergence of new pathogens. Ultimately, it underscores the need for more accurate predictive models and improved surveillance to better anticipate and mitigate future spillover events. [ABSTRACT FROM AUTHOR]

Weaning weight is a key indicator of the early growth performance of cattle. An understanding of the genetic mechanisms underlying weaning weight will help increase the accuracy of selection of breeding animals. In order to identify candidate genes associated with weaning weight in Simmental-Holstein crossbred cattle, this study generated RNA-Sequencing (RNA-seq) data from 86 crossbred calves (37 males and 49 famales) and measured their weaning weight and body size traits (wither height, body length, chest girth, rump width, and rump length). Differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) were performed. A total of 498 differentially expressed genes (DEGs) were identified between the low weaning weight (LWW) group and the high weaning weight (HWW) group. Weaning weight was transcriptionally correlated (FDR < 0.05) with four of the eleven co-expression gene modules. By intersecting DEGs and hub genes of the four modules, we identified a final set of 37 candidate genes enriched in growth, development, or immune-related processes. In addition, one co-expression module was significantly correlated with all the five body size traits (P < 0.05), from which MX1 was identified as a key candidate gene through protein-protein interaction (PPI) analysis of hub genes. Further evidence from cattle transcriptome-wide association study analysis (TWAS) and human phenome-wide association study (PheWAS) validated significant associations of CACNA1S, SEMA7A, VCAN, CD101, CD19, and CSF2RB with growth and development traits (P < 0.05). Notably, CACNA1S and CD19 were also associated with typical immune traits such as B cell proliferation, differentiation, and activation. In conclusion, this study reveals new candidate genes significantly associated with weaning weight in Simmental-Holstein crossbred cattle, providing a basis for further exploration of the genetic mechanisms behind growth traits of cattle. [ABSTRACT FROM AUTHOR]

Simple Summary: In recent decades, interest in the use of donkey milk for human nutrition has increased, since it may represent a natural substitute for cow's milk for children affected by milk protein allergies. The functional peculiarities of donkey milk are mainly linked to its casein content comparable to that of human milk. This study provides a thorough analysis of transcript isoforms generated by two αs2-casein-encoding genes (CSN1S2 I and CSN1S2 II) in donkeys and a detection of significant genetic diversity at both loci along with sequence comparisons across species. In particular, a key mutation affecting exon 17 splicing in CSN1S2 I was identified, and a genotyping method was developed. These data represent an important step in the understanding of the expression regulation of these genes in donkeys and a useful tool for the genetic improvement of donkey milk production that fulfils special consumer requirements. The αs2-casein is a phosphoprotein secreted in the milk of most mammals, and it is the most hydrophilic of all caseins. Contrary to genes found in ruminants, in donkeys two different encoding genes for donkey αs2-casein (CSN1S2 I and CSN1S2 II) have been identified. However, unlike in ruminants, the variability at these loci has not been characterized in detail in donkeys until now. In this study, we analyze the transcript profile of the donkey CSN1S2 I and CSN1S2 II genes, and we identify and describe the variability of these loci in the Ragusana and Amiatina breeds reared in Italy. The analysis of the CSN1S2 I Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) products and subsequent sequencing showed, in addition to correctly spliced mRNA, seven other minor mRNAs resulting from differential splicing events involving, in various combinations, entire exons (4, 5, 6, and 11), parts of exons (5′ or 3′ end of exon 17), or the recognition of intronic sequences as an exon (exon 12′). Similarly, the transcription analysis of the CSN1S2 II gene revealed a remarkable variability in splicing events, mainly concerning the alternative insertion of an extra exon 7 (named 7′); the first 33 bp of exon 13; or the alternative skipping of exons 9, 10, 11, 12, and 15, and their combinations. At the mRNA level for CSN1S2 I, seven SNPs were observed, five of which led to amino acid changes: p.T73>A, p.I109>V, p.I130>V, p.I146>T, and p.D217>Y. Similarly, nine SNPs were observed at the CSN1S2 II locus, seven of which are non-synonymous: p.L63>F, p.H70>Q, p.D90>N, p.129A>T, p.H131>Y, p.E144>G, and p.F157>S. In addition, the DNA sequencing of exon 17 and flanking introns of the CSN1S2 I gene revealed a G>A transition at the splice acceptor site of CSN1S2 I exon 17 (FM946022.1:c.375-1G>A), resulting in an allele-specific skipping of the first 15 nucleotides of this exon, which encode the peptide 176NKINQ180, and the recognition of an in-frame cryptic splicing acceptor site: arAACAAAATCAACCAG. A genotyping method based on restriction fragment length polymorphism (XbaI PCR-RFLP) was set up for this SNP. In the total population studied (105 Ragusana and 14 Amiatina donkeys), the A allele had a frequency of 0.2437 with no evidence of deviation from the Hardy–Weinberg equilibrium. This study adds new knowledge regarding the genetic variability of αs2-caseins in donkeys and may contribute significantly to the genetic improvement of milk production for this species. [ABSTRACT FROM AUTHOR]