Search

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Newcastle Disease Virus (NDV) genotype VII (GVII) is becoming the predominant strain of NDV in the poultry industry. It causes high mortality even in vaccinated chickens with a common NDV genotype II vaccine (GII-vacc). To overcome this, the killed GVII vaccine has been used to prevent NDV outbreaks. However, the debate about vaccine differences remains ongoing. Hence, this study investigated the difference in chickens' responses to the two vaccines at the molecular level. The spleen transcriptomes from vaccinated chickens reveal that GVII-vacc affected the immune response by downregulating neuroinflammation. It also enhanced a synaptogenesis pathway that operates typically in the nervous system, suggesting a mechanism for the neurotrophic effect of this strain. We speculated that the down-regulated immune system regulation correlated with protecting the nervous system from excess leukocytes and cytokine activity. In contrast, GII-vacc inhibited apoptosis by downregulating PERK/ATF4/CHOP as part of the unfolded protein response pathway but did not affect the expression of the same synaptogenesis pathway. Thus, the application of GVII-vacc needs to be considered in countries where GVII is the leading cause of NDV outbreaks. The predicted molecular signatures may also be used in developing new vaccines that trigger specific genes in the immune system in combating NDV outbreaks.

This study delves into the pathogenesis of virulent genotype VII strains of the Newcastle disease virus (NDV), focusing on experimentally infected birds. Predominant and consistent lesions observed include bursal atrophy and extensive depletion of all lymphoid tissues. Immunohistochemistry (IHC) analysis, targeting apoptosis (Caspase-3), necroptosis (MLKL), and NDV markers, indicates that bursal atrophy is linked to a non-apoptotic programmed cell death pathway known as "necroptosis". Repair assisted damage detection (RADD) of the bursa reveal oxidative DNA damage patterns consistent with programmed cell death, aligning with MLKL expression. Contrastingly, in the spleen, our findings suggest that necrosis (non-programmed cell death) predominantly contributes to lymphoid depletion. This conclusion is supported by evidence of karyorrhexis, fibrinous inflammation, RADD analyses, and IHC. Moreover, in addition to being pathogenic in its own right, NDV caused extensive and rapid lymphoid depletion that should be expected to contribute to profound immunosuppression. The elucidation of necroptosis in NDV-infected chickens provides a good rationale to investigate this mechanism in other paramyxoviral diseases such as human measles.

The whole genomic sequence of fowl adenovirus C (FAdV-4) strain FAdV-4/Pasouk, isolated from chickens with hepatitis-hydropericardium syndrome (HHS) from an outbreak in Iran, has been deposited in GenBank under accession number ON652872. Notably, this FAdV-4 isolate exhibited significant genetic similarities to contemporary isolates originating from China, indicating a shared ancestry., Competing Interests: The authors declare no conflict of interest.

Introduction: Sex differences in prenatal growth may contribute to sex-dependent programming effects on postnatal phenotype., Methods: We integrated for the first time phenotypic, histomorphological, clinico-chemical, endocrine and gene expression analyses in a single species, the bovine conceptus at mid-gestation., Results: We demonstrate that by mid-gestation, before the onset of accelerated growth, the female conceptus displays asymmetric lower growth compared to males. Female fetuses were smaller with lower ponderal index and organ weights than males. However, their brain:body weight, brain:liver weight and heart:body weight ratios were higher than in males, indicating brain and heart 'sparing'. The female placenta weighed less and had lower volumes of trophoblast and fetal connective tissue than the male placenta. Female umbilical cord vessel diameters were smaller, and female-specific relationships of body weight and brain:liver weight ratios with cord vessel diameters indicated that the umbilico-placental vascular system creates a growth-limiting environment where blood flow is redistributed to protect brain and heart growth. Clinico-chemical indicators of liver perfusion support this female-specific growth-limiting phenotype, while lower insulin-like growth factor 2 (IGF2) gene expression in brain and heart, and lower circulating IGF2, implicate female-specific modulation of key endocrine mediators by nutrient supply., Conclusion: This mode of female development may increase resilience to environmental perturbations in utero and contribute to sex-bias in programming outcomes including susceptibility to non-communicable diseases., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Estrella, Gatford, Xiang, Javadmanesh, Ghanipoor-Samami, Nattrass, Shuaib, McAllister, Beckman, Thomsen, Clifton, Owens, Roberts, Hiendleder and Kind.)

Background: Most DNA methylation studies have used a single reference genome with little attention paid to the bias introduced due to the reference chosen. Reference genome artifacts and genetic variation, including single nucleotide polymorphisms (SNPs) and structural variants (SVs), can lead to differences in methylation sites (CpGs) between individuals of the same species. We analyzed whole-genome bisulfite sequencing data from the fetal liver of Angus (Bos taurus taurus), Brahman (Bos taurus indicus), and reciprocally crossed samples. Using reference genomes for each breed from the Bovine Pangenome Consortium, we investigated the influence of reference genome choice on the breed and parent-of-origin effects in methylome analyses., Results: Our findings revealed that ∼75% of CpG sites were shared between Angus and Brahman, ∼5% were breed specific, and ∼20% were unresolved. We demonstrated up to ∼2% quantification bias in global methylation when an incorrect reference genome was used. Furthermore, we found that SNPs impacted CpGs 13 times more than other autosomal sites (P < $5 \times {10}^{ - 324}$) and SVs contained 1.18 times (P < $5 \times {10}^{ - 324}$) more CpGs than non-SVs. We found a poor overlap between differentially methylated regions (DMRs) and differentially expressed genes (DEGs) and suggest that DMRs may be impacting enhancers that target these DEGs. DMRs overlapped with imprinted genes, of which 1, DGAT1, which is important for fat metabolism and weight gain, was found in the breed-specific and sire-of-origin comparisons., Conclusions: This work demonstrates the need to consider reference genome effects to explore genetic and epigenetic differences accurately and identify DMRs involved in controlling certain genes., (© The Author(s) 2024. Published by Oxford University Press GigaScience.)

More people in the world depend on water buffalo for their livelihoods than on any other domesticated animals, but its genetics is still not extensively explored. The 1000 Buffalo Genomes Project (1000BGP) provides genetic resources for global buffalo population study and tools to breed more sustainable and productive buffaloes. Here we report the most contiguous swamp buffalo genome assembly (PCC&#95;UOA&#95;SB&#95;1v2) with substantial resolution of telomeric and centromeric repeats, ∼4-fold more contiguous than the existing reference river buffalo assembly and exceeding a recently published male swamp buffalo genome. This assembly was used along with the current reference to align 140 water buffalo short-read sequences and produce a public genetic resource with an average of ∼41 million single nucleotide polymorphisms per swamp and river buffalo genome. Comparison of the swamp and river buffalo sequences showed ∼1.5% genetic differences, and estimated divergence time occurred 3.1 million years ago (95% CI, 2.6-4.9). The open science model employed in the 1000BGP provides a key genomic resource and tools for a species with global economic relevance., (© The Author(s) 2024. Published by Oxford University Press on behalf of GigaScience.)

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman ( Bos taurus indicus ) and Angus ( Bos taurus taurus ) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets ( p -value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 MacPhillamy, Ren, Chen, Hiendleder and Low.)

Background: The presence of Aggregated Lymphoid Nodules Area (ALNA) is a notable anatomical characteristic observed in the abomasum of Bactrian camels. This area is comprised of two separate regions, namely the Reticular Mucosal Folds Region (RMFR) and the Longitudinal Mucosal Folds Region (LMFR). The histological properties of ALNA exhibit significant similarities to those of Peyer's patches (PPs) found in the gastrointestinal system. The functional characteristics of ALNA were examined in relation to mucosal immunity in the gastrointestinal system., Results: We used iTRAQ-based proteomic analysis on twelve Bactrian camels to measure the amount of proteins expressed in ALNA. In the experiment, we sampled the RMFR and LMFR separately from the ALNA and compared their proteomic quantification results with samples from the PPs. A total of 1253 proteins were identified, among which 39 differentially expressed proteins (DEPs) were found between RMFR and PPs, 33 DEPs were found between LMFR and PPs, and 22 DEPs were found between LMFR and RMFR. The proteins FLNA, MYH11, and HSPB1 were chosen for validation using the enzyme-linked immunosorbent assay (ELISA), and the observed expression profiles were found to be in agreement with the results obtained from the iTRAQ study. The InnateDB database was utilized to get data pertaining to immune-associated proteins in ALNA. It was observed that a significant proportion, specifically 76.6%, of these proteins were found to be associated with the same orthogroups as human immune-related genes. These proteins are acknowledged to be associated with a diverse range of functions, encompassing the uptake, processing and presentation of antigens, activation of lymphocytes, the signaling pathways of T-cell and B-cell receptors, and the control of actin polymerization., Conclusions: The experimental results suggest that there are parallels in the immune-related proteins found in ALNA and PPs. Although there are variations in the structures of LMFR and RMFR, the proteins produced in both structures exhibit a high degree of similarity and perform comparable functions in the context of mucosal immune responses., (© 2023. BioMed Central Ltd., part of Springer Nature.)

The Iso-Seq method of full-length cDNA sequencing is suitable to quantify differentially expressed genes (DEGs), transcripts (DETs) and transcript usage (DTU). However, the higher cost of Iso-Seq relative to RNA-seq has limited the comparison of both methods. Transcript abundance estimated by RNA-seq and deep Iso-Seq data for fetal liver from two cattle subspecies were compared to evaluate concordance. Inter-sample correlation of gene- and transcript-level abundance was higher within technology than between technologies. Identification of DEGs between the cattle subspecies depended on sequencing method with only 44 genes identified by both that included 6 novel genes annotated by Iso-Seq. There was a pronounced difference between Iso-Seq and RNA-seq results at transcript-level wherein Iso-Seq revealed several magnitudes more transcript abundance and usage differences between subspecies. Factors influencing DEG identification included size selection during Iso-Seq library preparation, average transcript abundance, multi-mapping of RNA-seq reads to the reference genome, and overlapping coordinates of genes. Some DEGs called by RNA-seq alone appear to be sequence duplication artifacts. Among the 44 DEGs identified by both technologies some play a role in immune system, thyroid function and cell growth. Iso-Seq revealed hidden transcriptional complexity in DEGs, DETs and DTU genes between cattle subspecies previously missed by RNA-seq., (© 2023. The Author(s).)

Mycobacterium avium ssp. paratuberculosis (MAP), causes Johne's disease (JD), or paratuberculosis, a chronic enteritis of ruminants, which in goats is characterized by ileal lesions. The work described here is a case-control association study using the Illumina Caprine SNP50 BeadChip to unravel the genes involved in susceptibility of goats to JD. Goats in herds with a high occurrence of Johne's disease were classified as healthy or infected based on the level of serum antibodies against MAP, and 331 animals were selected for the association study. Goats belonged to the Jonica (157) and Siriana breeds (174). Whole-genome association analysis identified one region suggestive of significance associated with an antibody response to MAP on chromosome 7 (p-value = 1.23 × 10 -5 ). These results provide evidence for genetic loci involved in the antibody response to MAP in goats., (© 2022 Stichting International Foundation for Animal Genetics.)

Aims: The aim of this study was to assess a phage-displayed MilA protein of Myc. bovis in an indirect ELISA for the detection of Myc. bovis antibodies in milk samples., Methods and Results: The desired sequence of milA gene was synthesized and cloned into pCANTAB-F12 phagemid vector. The expression of the MilA on the phage surface was confirmed by Western blotting. The recombinant phage was used in the development of an indirect ELISA to detect Myc. bovis antibodies in milk samples. There was a significant agreement between the results of phage-based ELISA and recombinant GST-MilA ELISA for the detection of Myc. bovis antibodies in milk samples., Conclusions: The inexpensive and convenient phage-based ELISA can be used instead of recombinant protein/peptide ELISA as an initial screening of Myc. bovis-associated mastitis., Significance and Impact of Study: Mastitis associated with Myc. bovis is a continuous and serious problem in the dairy industry. Sero-monitoring of Myc. bovis infection cases are one of the key factors for surveillance of the infections in dairy farms. Despite the existence of some commercially serological assays for Myc. bovis antibodies, they have some limitations regarding their sensitivity and availability. The development of accurate diagnosis tools could contribute to control programmes of Myc. bovis-associated mastitis in the dairy herds., (© 2022 The Authors. Journal of Applied Microbiology published by John Wiley & Sons Ltd on behalf of Society for Applied Microbiology.)

Experimental and epidemiological studies suggest that maternal nutritional status during early pregnancy, including the period around the time of conception, may induce long-lasting epigenetic changes in the offspring. However, this remains largely unexplored in livestock. Therefore, the objective of this study was to evaluate if modification of the maternal diet of sheep (CTR: control; UND: 50% undernutrition) during the periconceptional period (42 d in total: -14/+28 from mating), would impact CpG methylation in muscle tissue (Longissimus dorsi) of adult offspring (11.5 mo old). Reduced representation bisulfite sequencing identified 262 (Edge-R, FDR < 0.05) and 686 (logistic regression, FDR < 0.001) differentially methylated regions (DMRs) between the UND and CTR groups. Gene ontology analysis identified genes related to development, functions of the muscular system, and steroid hormone receptor activity within the DMRs. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the muscle of the resulting offspring, with possible implications for cardiac dysfunction, muscle physiology, and meat production., (© The Author(s) 2022. Published by Oxford University Press on behalf of the American Society of Animal Science. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: The gaur (Bos gaurus) is the largest extant wild bovine species, native to South and Southeast Asia, with unique traits, and is listed as vulnerable by the International Union for Conservation of Nature (IUCN)., Results: We report the first gaur reference genome and identify three biological pathways including lysozyme activity, proton transmembrane transporter activity, and oxygen transport with significant changes in gene copy number in gaur compared to other mammals. These may reflect adaptation to challenges related to climate and nutrition. Comparative analyses with domesticated indicine (Bos indicus) and taurine (Bos taurus) cattle revealed genomic signatures of artificial selection, including the expansion of sperm odorant receptor genes in domesticated cattle, which may have important implications for understanding selection for male fertility., Conclusions: Apart from aiding dissection of economically important traits, the gaur genome will also provide the foundation to conserve the species., (© 2022. The Author(s).)

Newcastle disease virus genotype VII (NDV-GVII) is a highly contagious pathogen responsible for pandemics that have caused devastating economic losses in the poultry industry. Several features in the transcription of NDV mRNA, including differentially expressed genes across the viral genome, are shared with that for other single, non-segmented, negative-strand viruses. Previous studies measuring viral gene expression using northern blotting indicated that the NDV transcription produced non-equimolar levels of viral mRNAs. However, deep high-throughput sequencing of virus-infected tissues can provide a better insight into the patterns of viral transcription. In this report, the transcription pattern of virulent NDV-GVII was analysed using RNA-seq and qRT-PCR. This study revealed the transcriptional profiling of these highly pathogenic NDV-GVII genes: NP:P:M:F:HN:L, in which there was a slight attenuation at the NP:P and HN:L gene boundaries. Our result also provides a fully comprehensive qPCR protocol for measuring viral transcript abundance that may be more convenient for laboratories where accessing RNA-seq is not feasible., (© 2022. The Author(s).)

Linkage disequilibrium (LD) across the genome provides information to identify the genes and variations related to quantitative traits in genome-wide association studies (GWAS) and for the implementation of genomic selection (GS). LD can also be used to evaluate genetic diversity and population structure and reveal genomic regions affected by selection. LD structure and Ne were assessed in a set of 83 water buffaloes, comprising Azeri (AZI), Khuzestani (KHU), and Mazandarani (MAZ) breeds from Iran, Kundi (KUN) and Nili-Ravi (NIL) from Pakistan, Anatolian (ANA) buffalo from Turkey, and buffalo from Egypt (EGY). The values of corrected r 2 (defined as the correlation between two loci) of adjacent SNPs for three pooled Iranian breeds (IRI), ANA, EGY, and two pooled Pakistani breeds (PAK) populations were 0.24, 0.28, 0.27, and 0.22, respectively. The corrected r 2 between SNPs decreased with increasing physical distance from 100 Kb to 1 Mb. The LD values for IRI, ANA, EGY, and PAK populations were 0.16, 0.23, 0.24, and 0.21 for less than 100Kb, respectively, which reduced rapidly to 0.018, 0.042, 0.059, and 0.024, for a distance of 1 Mb. In all the populations, the decay rate was low for distances greater than 2Mb, up to the longest studied distance (15 Mb). The r 2 values for adjacent SNPs in unrelated samples indicated that the Affymetrix Axiom 90 K SNP genomic array was suitable for GWAS and GS in these populations. The persistency of LD phase (PLDP) between populations was assessed, and results showed that PLPD values between the populations were more than 0.9 for distances of less than 100 Kb. The Ne in the recent generations has declined to the extent that breeding plans are urgently required to ensure that these buffalo populations are not at risk of being lost. We found that results are affected by sample size, which could be partially corrected for; however, additional data should be obtained to be confident of the results., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Rahimmadar, Ghaffari, Mokhber and Williams.)

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Newcastle disease virus (NDV) has caused significant outbreaks in South-East Asia, particularly in Indonesia in recent years. Recently emerged genotype VII NDVs (NDV-GVII) have shifted their tropism from gastrointestinal/respiratory tropism to a lymphotropic virus, invading lymphoid organs including spleen and bursa of Fabricius to cause profound lymphoid depletion. In this study, we aimed to identify candidate genes and biological pathways that contribute to the disease caused by this velogenic NDV-GVII. A transcriptomic analysis based on RNA-Seq of spleen was performed in chickens challenged with NDV-GVII and a control group. In total, 6361 genes were differentially expressed that included 3506 up-regulated genes and 2855 down-regulated genes. Real-Time PCR of ten selected genes validated the RNA-Seq results as the correlation between them is 0.98. Functional and network analysis of Differentially Expressed Genes (DEGs) showed altered regulation of ElF2 signalling, mTOR signalling, proliferation of cells of the lymphoid system, signalling by Rho family GTPases and synaptogenesis signalling in spleen. We have also identified modified expression of IFIT5, PI3K, AGT and PLP1 genes in NDV-GVII infected chickens. Our findings in activation of autophagy-mediated cell death, lymphotropic and synaptogenesis signalling pathways provide new insights into the molecular pathogenesis of this newly emerged NDV-GVII., (© 2021. The Author(s).)

Eukaryotic genomes are highly organised within the nucleus of a cell, allowing widely dispersed regulatory elements such as enhancers to interact with gene promoters through physical contacts in three-dimensional space. Recent chromosome conformation capture methodologies such as Hi-C have enabled the analysis of interacting regions of the genome providing a valuable insight into the three-dimensional organisation of the chromatin in the nucleus, including chromosome compartmentalisation and gene expression. Complicating the analysis of Hi-C data, however, is the massive amount of identified interactions, many of which do not directly drive gene function, thus hindering the identification of potentially biologically functional 3D interactions. In this review, we collate and examine the downstream analysis of Hi-C data with particular focus on methods that prioritise potentially functional interactions. We classify three groups of approaches: structural-based discovery methods, e.g. A/B compartments and topologically associated domains, detection of statistically significant chromatin interactions, and the use of epigenomic data integration to narrow down useful interaction information. Careful use of these three approaches is crucial to successfully identifying potentially functional interactions within the genome., (© 2021. The Author(s).)

The ability to assess the welfare of animals is dependent on our ability to accurately determine their emotional (affective) state, with particular emphasis being placed on the identification of positive emotions. The challenge remains that current physiological and behavioral indices are either unable to distinguish between positive and negative emotional states, or they are simply not suitable for a production environment. Therefore, the development of novel measures of animal emotion is a necessity. Here we investigated the efficacy of microRNA (miRNA) in the brain and blood as biomarkers of emotional state in the pig. Female Large White × Landrace pigs ( n = 24) were selected at weaning and trained to perform a judgment bias test (JBT), before being exposed for 5 weeks to either enriched ( n = 12) or barren housing ( n = 12) conditions. Pigs were tested on the JBT once prior to treatment, and immediately following treatment. MiRNA and neurotransmitters were analyzed in blood and brain tissue after euthanasia. Treatment had no effect on the outcomes of the JBT. There was also no effect of treatment on miRNA expression in blood or the brain (FDR p > 0.05). However, pigs exposed to enriched housing had elevated dopamine within the striatum compared to pigs in barren housing ( p = 0.02). The results imply that either (a) miRNAs are not likely to be valid biomarkers of a positive affective state, at least under the type of conditions employed in this study, or (b) that the study design used to modify affective state was not able to create differential affective states, and therefore establish the validity of miRNA as biomarkers.

Research on fatty acids (FA) is important because their intake is related to human health. NIRS can be a useful tool to estimate the FA of beef but due to the high moisture and the high absorbance of water makes it difficult to calibrate the analyses. This work evaluated near-infrared reflectance spectroscopy as a tool to assess the total fatty acid composition and the phospholipid fraction of fatty acids of beef using freeze-dried meat. An average of 22 unrelated pure breed young bulls from 15 European breeds were reared on a common concentrate-based diet. A total of 332 longissimus thoracis steaks were analysed for fatty acid composition and a freeze-dried sample was subjected to near-infrared spectral analysis. 220 samples (67%) were used as a calibration set with the remaining 110 (33%) being used for validation of the models obtained. There was a large variation in the total FA concentration across the animals giving a good data set for the analysis and whilst the coefficient of variation was nearly 68% for the monounsaturated FA it was only 27% for the polyunsaturated fatty acids (PUFA). PLS method was used to develop the prediction models. The models for the phospholipid fraction had a low R 2 p and high standard error, while models for neutral lipid had the best performance, in general. It was not possible to obtain a good prediction of many individual PUFA concentrations being present at low concentrations and less variable than other FA. The best models were developed for Total FA, saturated FA, 9c18:1 and 16:1 with R 2 p greater than 0.76. This study indicates that NIRS is a feasible and useful tool for screening purposes and it has the potential to predict most of the FA of freeze-dried beef.

Background: There are two genetically distinct subspecies of cattle, Bos taurus taurus and Bos taurus indicus, which arose from independent domestication events. The two types of cattle show substantial phenotypic differences, some of which emerge during fetal development and are reflected in birth outcomes, including birth weight. We explored gene expression profiles in the placenta and four fetal tissues at mid-gestation from one taurine (Bos taurus taurus; Angus) and one indicine (Bos taurus indicus; Brahman) breed and their reciprocal crosses., Results: In total 120 samples were analysed from a pure taurine breed, an indicine breed and their reciprocal cross fetuses, which identified 6456 differentially expressed genes (DEGs) between the two pure breeds in at least one fetal tissue of which 110 genes were differentially expressed in all five tissues examined. DEGs shared across tissues were enriched for pathways related to immune and stress response functions. Only the liver had a substantial number of DEGs when reciprocal crossed were compared among which 310 DEGs were found to be in common with DEGs identified between purebred livers; these DEGs were significantly enriched for metabolic process GO terms. Analysis of DEGs across purebred and crossbred tissues suggested an additive expression pattern for most genes, where both paternal and maternal alleles contributed to variation in gene expression levels. However, expression of 5% of DEGs in each tissue was consistent with parent of origin effects, with both paternal and maternal dominance effects identified., Conclusions: These data identify candidate genes potentially driving the tissue-specific differences between these taurine and indicine breeds and provide a biological insight into parental genome effects underlying phenotypic differences in bovine fetal development.

Subclinical mastitis, an economically challenging disease of dairy cattle, is associated with an increased use of antimicrobials which reduces milk quantity and quality. It is more common than clinical mastitis and far more difficult to detect. Recently, much attention has been paid to the development of machine-learning expert systems for early detection of subclinical mastitis from milking features. However, differences between animals within a farm as well as between farms, particularly across multiple years, are major obstacles to the generalisation of machine learning models. Here, for the first time, we integrated scaling by quartiling with classification based on associations in a multi-year study to deal with farm heterogeneity by discovery of multiple patterns towards mastitis. The data were obtained from one farm comprising Holstein Friesian cows in Ongaonga, New Zealand, using an electronic automated monitoring system. The data collection was repeated annually over 3 consecutive years. Some discovered rules, such as when the milking peak flow is low, electrical conductivity (EC) of milk is low, milk lactose is low, milk fat is high, and milk volume is low, the cow has subclinical mastitis, reached high confidence (>70%) in multiple years. On averages, over 3 years, low level of milk lactose and high value of milk EC were part of 93% and 83.8% of all subclinical mastitis detecting rules, offering a reproducible pattern of subclinical mastitis detection. The scaled year-independent combinational rules provide an easy-to-apply and cost-effective machine-learning expert system for early detection of hidden mastitis using milking parameters.

Somatic cell nuclear transfer (SCNT) is a key technology with broad applications that range from production of cloned farm animals to derivation of patient-matched stem cells or production of humanized animal organs for xenotransplantation. However, effects of aberrant epigenetic reprogramming on gene expression compromise cell and organ phenotype, resulting in low success rate of SCNT. Standard SCNT procedures include enucleation of recipient oocytes before the nuclear donor cell is introduced. Enucleation removes not only the spindle apparatus and chromosomes of the oocyte but also the perinuclear, mitochondria rich, ooplasm. Here, we use a Bos taurus SCNT model with in vitro fertilized (IVF) and in vivo conceived controls to demonstrate a ∼50% reduction in mitochondrial DNA (mtDNA) in the liver and skeletal muscle, but not the brain, of SCNT fetuses at day 80 of gestation. In the muscle, we also observed significantly reduced transcript abundances of mtDNA-encoded subunits of the respiratory chain. Importantly, mtDNA content and mtDNA transcript abundances correlate with hepatomegaly and muscle hypertrophy of SCNT fetuses. Expression of selected nuclear-encoded genes pivotal for mtDNA replication was similar to controls, arguing against an indirect epigenetic nuclear reprogramming effect on mtDNA amount. We conclude that mtDNA depletion is a major signature of perturbations after SCNT. We further propose that mitochondrial perturbation in interaction with incomplete nuclear reprogramming drives abnormal epigenetic features and correlated phenotypes, a concept supported by previously reported effects of mtDNA depletion on the epigenome and the pleiotropic phenotypic effects of mtDNA depletion in humans. This provides a novel perspective on the reprogramming process and opens new avenues to improve SCNT protocols for healthy embryo and tissue development., Competing Interests: MR is employed by the company Bayern-Genetik GmbH, Grub, Germany. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Bebbere, Ulbrich, Giller, Zakhartchenko, Reichenbach, Reichenbach, Verma, Wolf, Ledda and Hiendleder.)

Haplotype-resolved genome assemblies are important for understanding how combinations of variants impact phenotypes. To date, these assemblies have been best created with complex protocols, such as cultured cells that contain a single-haplotype (haploid) genome, single cells where haplotypes are separated, or co-sequencing of parental genomes in a trio-based approach. These approaches are impractical in most situations. To address this issue, we present FALCON-Phase, a phasing tool that uses ultra-long-range Hi-C chromatin interaction data to extend phase blocks of partially-phased diploid assembles to chromosome or scaffold scale. FALCON-Phase uses the inherent phasing information in Hi-C reads, skipping variant calling, and reduces the computational complexity of phasing. Our method is validated on three benchmark datasets generated as part of the Vertebrate Genomes Project (VGP), including human, cow, and zebra finch, for which high-quality, fully haplotype-resolved assemblies are available using the trio-based approach. FALCON-Phase is accurate without having parental data and performance is better in samples with higher heterozygosity. For cow and zebra finch the accuracy is 97% compared to 80-91% for human. FALCON-Phase is applicable to any draft assembly that contains long primary contigs and phased associate contigs.

Stress in livestock reduces productivity and is a welfare concern. At a physiological level, stress is associated with the activation of inflammatory responses and increased levels of harmful reactive oxygen species. Biomarkers that are indicative of stress could facilitate the identification of more stress-resilient animals. We examined twenty-one metabolic, immune response, and liver function biomarkers that have been associated with stress in 416 Italian Simmental and 436 Italian Holstein cows which were genotyped for 150K SNPs. Single-SNP and haplotype-based genome-wide association studies were carried out to assess whether the variation in the levels in these biomarkers is under genetic control and to identify the genomic loci involved. Significant associations were found for the plasma levels of ceruloplasmin ( Bos taurus chromosome 1-BTA1), paraoxonase (BTA4) and γ-glutamyl transferase (BTA17) in the individual breed analysis that coincided with the position of the genes coding for these proteins, suggesting that their expression is under cis-regulation. A meta-analysis of both breeds identified additional significant associations with paraoxonase on BTA 16 and 26. Finding genetic associations with variations in the levels of these biomarkers suggests that the selection for high or low levels of expression could be achieved rapidly. Whether the level of expression of the biomarkers correlates with the response to stressful situations has yet to be determined.

Mastitis represents one of the major economic and health threats to the livestock sector associated with reduction in milk quality, loss of production and is a major reason for culling. Somatic cell score (SCS) is used as a criterion in breeding programmes to select cows genetically less susceptible to mastitis. The relevance of SCS as a predictor of udder health and susceptibility to mastitis is still untested in goats. In this study, two lines of French Alpine goats selected for extreme breeding values for somatic cell scores, one line with high SCS (HSCS) and the other with low SCS (LSCS), were used to test the hypothesis that the mammary response and function differed between the lines. The aim of the present study was to investigate differences in the early immune response in caprine mammary gland tissues challenged with Staphylococcus aureus, one of the main pathogens responsible for the intra-mammary infection in small ruminants, using transcriptomic and histopathology analyses. The comparison between HSCS and LSCS goat lines, showed differences in the response at the histological level for inflammation, presence of neutrophils and micro-abscess formation, and at the molecular level in the expression of CXCL8, IL-6, NFKBIZ and IL-1β. CXCL8 and CXCL2 genes, which showed a higher level of expression in the experimentally infected HSCS line. The molecular data and histopathology both suggested that following S. aureus infection, mobilization, recruitment, infiltration, and chemotaxis of neutrophil, leads to a more severe inflammation in the HSCS compared to LSCS animals. Our results represent an initial basis for further studies to unravel the genetic basis of early mastitis inflammatory responses and the selection of dairy animals more resistant to bacterial mastitis., (Copyright © 2021 Elsevier B.V. All rights reserved.)

The swamp buffalo is a domesticated animal commonly found in Southeast Asia. It is a highly valued agricultural animal for smallholders, but the production of this species has unfortunately declined in recent decades due to rising farm mechanization. While swamp buffalo still plays a role in farmland cultivation, this species' purposes has shifted from draft power to meat, milk, and hide production. The current status of swamp buffaloes in Southeast Asia is still understudied compared to its counterparts such as the riverine buffaloes and cattle. This review discusses the background of swamp buffalo, with an emphasis on recent work on this species in Southeast Asia, and associated genetics and genomics work such as cytogenetic studies, phylogeny, domestication and migration, genetic sequences and resources. Recent challenges to realize the potential of this species in the agriculture industry are also discussed. Limited genetic resource for swamp buffalo has called for more genomics work to be done on this species including decoding its genome. As the economy progresses and farm mechanization increases, research and development for swamp buffaloes are focused on enhancing its productivity through understanding the genetics of agriculturally important traits. The use of genomic markers is a powerful tool to efficiently utilize the potential of this animal for food security and animal conservation. Understanding its genetics and retaining and maximizing its adaptability to harsher environments are a strategic move for food security in poorer nations in Southeast Asia in the face of climate change., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Pineda, Flores, Herrera and Low.)

Background: Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types., Results: The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction., Conclusions: Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Major health conditions in sheep contribute to substantial economic losses throughout the sheepmeat supply chain in Australia. A systematic review was undertaken to explore the measurable impact of six conditions: arthritis, sheep measles, pleurisy, pneumonia, grass seeds and rib fractures, on the production of lamb and mutton across the meat value chain. Peer-reviewed scientific literature from three databases and non-peer-reviewed articles and reports from Australian government and non-government websites were searched between 11 and 17 November 2019. Original articles, including studies conducted in Australia and New Zealand, that had measurable impacts on conditions of interest were included. The search yielded 16 articles and reports and were classified as producer impact and/or processor impact studies. Mortalities were quantified for pneumonia and arthritis, with pneumonia having the highest impact for producers. Grass seed infestation resulted in the highest impact on carcase and liveweight losses compared to arthritis and pneumonia. Arthritis had the highest trim weight losses for both lamb and mutton and the highest rate of carcase condemnation. Grass seed was the only condition where other impacts on the processor (chain speed and staff relocation to the boning room) were quantified. Although quantifiable production and processing losses were available for some conditions, this review has highlighted that limited quantifiable data based on scientifically sound research were not available for other conditions. The evidence for some conditions found in this review can be used to target future research activities and to further assist producers in making informed management decisions on prevention and control., (© 2020 Australian Veterinary Association.)

New technologies and analysis methods are enabling genomic structural variants (SVs) to be detected with ever-increasing accuracy, resolution and comprehensiveness. To help translate these methods to routine research and clinical practice, we developed a sequence-resolved benchmark set for identification of both false-negative and false-positive germline large insertions and deletions. To create this benchmark for a broadly consented son in a Personal Genome Project trio with broadly available cells and DNA, the Genome in a Bottle Consortium integrated 19 sequence-resolved variant calling methods from diverse technologies. The final benchmark set contains 12,745 isolated, sequence-resolved insertion (7,281) and deletion (5,464) calls ≥50 base pairs (bp). The Tier 1 benchmark regions, for which any extra calls are putative false positives, cover 2.51 Gbp and 5,262 insertions and 4,095 deletions supported by ≥1 diploid assembly. We demonstrate that the benchmark set reliably identifies false negatives and false positives in high-quality SV callsets from short-, linked- and long-read sequencing and optical mapping.

Despite the considerable advances that have been made to improve dairy food safety, there is rising concern that pasteurization is not sufficient for the destruction of plasmid-mediated antimicrobial resistance (AMR) genes of resistant bacteria and could stimulate bacteria to enter into a viable but nonculturable (VBNC) state. In the current study, we surveyed the prevalence of 1 genomic and 9 plasmid-mediated AMR genes in 100 samples (bulk tank milk and milk filter socks) at the farm level and 152 commercial milk samples (pasteurized and UHT milks) and assessed the VBNC state in dairy bacteria. Results revealed that sul2 was the most prevalent plasmid-mediated gene in milk filter socks (96%), bulk tank milk (48%), pasteurized milk (68%), and UHT (43%) milk; in contrast, mecA was not detected in any sample. Additionally, commercial pasteurization (as currently practiced) failed to decrease the prevalence of the bla TEM-B1 (43%), tetK (30%), and tetA (55%) plasmid-mediated AMR genes; thus, commercial pasteurization may be one of the factors creating the VBNC state in some dairy bacteria. Continued research is necessary to identify bacterial species entering the VBNC state after pasteurization, to assess their potential hazard level and shed more light on the expression and possibility of horizontal gene transfer of those plasmid-mediated AMR genes., (Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

The role of microbiota in gut-brain communication has led to the development of probiotics promoting brain health. Here we report a genomic study of a Lactobacillus fermentum PS150 and its patented bioactive protein, elongation factor Tu (EF-Tu), which is associated with cognitive improvement in rats. The L. fermentum PS150 circular chromosome is 2,238,401 bp and it consists of 2281 genes. Chromosome comparisons with other L. fermentum strains highlighted a cluster of glycosyltransferases as potential candidate probiotic factors besides EF-Tu. Molecular evolutionary analyses on EF-Tu genes (tuf) in 235 bacteria species revealed one to three copies of the gene per genome. Seven tuf pseudogenes were found and three species only possessed pseudogenes, which is an unprecedented finding. Protein variability analysis of EF-Tu showed five highly variable residues (40 K, 41G, 42 L, 44 K, and 46E) on the protein surface, which warrant further investigation regarding their potential roles as binding sites., (Copyright © 2020. Published by Elsevier Inc.)

More people globally depend on the water buffalo than any other domesticated species, and as the most closely related domesticated species to cattle they can provide important insights into the shared evolutionary basis of domestication. Here, we sequence the genomes of 79 water buffalo across seven breeds and compare patterns of between breed selective sweeps with those seen for 294 cattle genomes representing 13 global breeds. The genomic regions under selection between cattle breeds significantly overlap regions linked to stature in human genetic studies, with a disproportionate number of these loci also shown to be under selection between water buffalo breeds. Investigation of potential functional variants in the water buffalo genome identifies a rare example of convergent domestication down to the same mutation having independently occurred and been selected for across domesticated species. Cross-species comparisons of recent selective sweeps can consequently help identify and refine important loci linked to domestication.

Weaning is a stressful time for piglets, often leading to weight loss and is associated with increased morbidity and mortality. A leading cause for these post-weaning problems is enteric dysbiosis and methods to improve piglet health at this crucial developmental stage are needed. This study aimed to determine whether an enteric dysbiosis caused by weaning could be corrected via a faecal microbiota transplantation (FMT) from healthy piglets from a previous wean. Two or four focal piglets per litter were assigned to one of two treatments; FMT two days post weaning ( n = 21; FMT) or a control which received saline two days post weaning ( n = 21; CON). FMT consisted of homogenised donor faeces administered orally at 3 mL/kg. Weaning occurred at 18 days of age and weights and faecal samples were collected on days 18, 20, 24 and 35. 16S rRNA amplicon analysis was used to assess the faecal microbiota of piglets. FMT increased Shannon's diversity post weaning ( p < 0.001) and reduced the scratch score observed at 24 days of age ( p < 0.001). The bacterial populations significantly differed in composition at each taxonomic level. In FMT pigs, significant increases in potentially pathogenic Escherichia coli were observed. However, increases in beneficial bacteria Lactobacillus mucosae and genera Fibrobacteres and Bacteroidetes were also observed in FMT treated animals. To our knowledge, this is the first study to observe a significant effect on piglet faecal microbiota following a single FMT administered post weaning. Therefore, FMT post weaning can potentially alleviate enteric dysbiosis.

Diatoms are the largest group of heterokont algae with more than 100,000 species. As one of the single-celled photosynthetic organisms that inhabit marine, aquatic and terrestrial ecosystems, diatoms contribute ~ 45% of global primary production. Despite their ubiquity and environmental significance, very few diatom plastid genomes (plastomes) have been sequenced and studied. This study explored patterns of nucleotide substitution rates of diatom plastids across the entire suite of plastome protein-coding genes for 40 taxa representing the major clades. The highest substitution rate was lineage-specific within the araphid 2 taxon Astrosyne radiata and radial 2 taxon Proboscia sp. Rate heterogeneity was also evident in different functional classes and individual genes. Similar to land plants, proteins genes involved in photosynthetic metabolism have lower synonymous and nonsynonymous substitutions rates than those involved in transcription and translation. Significant positive correlations were identified between substitution rates and measures of genomic rearrangements, including indels and inversions, which is a similar result to what was found in legume plants. This work advances the understanding of the molecular evolution of diatom plastomes and provides a foundation for future studies.

In Earth's near-surface environments, gold biogeochemical cycling involves gold dissolution and precipitation processes, which are partly attributed to bacteria. These biogeochemical processes as well as abrasion (via physical transport) are known to act upon gold particles, thereby resulting in particle transformation including the development of pure secondary gold and altered morphology, respectively. While previous studies have inferred gold biogeochemical cycling from gold particles obtained from natural environments, little is known about how metal contamination in an environment could impact this cycle. Therefore, this study aims to infer how potentially toxic metal contaminants could affect the structure and chemistry of gold particles and therefore the biogeochemical cycling of gold. In doing so, river sediments and gold particles from the De Kaap Valley, South Africa, were analysed using both microanalytical and molecular techniques. Of the metal contaminants detected in the sediment, mercury can chemically interact with gold particles thereby directly altering particle morphology and "erasing" textural evidence indicative of particle transformation. Other metal contaminants (including mercury) indirectly affect gold cycling by exerting a selective pressure on bacteria living on the surface of gold particles. Particles harbouring gold-tolerant bacteria with diverse metal resistant genes, such as Arthrobacter sp. and Pseudomonas sp., contained nearly two times more secondary gold relative to particles harbouring bacteria with less gold-tolerance. In conclusion, metal contaminants can have a direct or indirect effect on gold biogeochemical cycling in natural environments impacted by anthropogenic activity., Competing Interests: Declaration of competing interest The author declares no competing interest., (Copyright © 2020 Elsevier B.V. All rights reserved.)

The full-genome sequences of strains chicken/Indonesia/Cilebut/010WJ/2015 and chicken/Indonesia/ITA/012WJ/1951, isolated in West Java, Indonesia, in 2015 and 1951, respectively, were examined. Chicken/Indonesia/Cilebut/010WJ/2015 (genotype VII) caused a 2015 disease outbreak in Indonesia, and chicken/Indonesia/ITA/012WJ/1951 (genotype VI) is used as a standard strain for challenge in Newcastle disease virus (NDV) vaccine trials., (Copyright © 2020 Pandarangga et al.)

The genomes of two newly emerged Newcastle disease virus strains, chicken/Indonesia/Mega/001WJ/2013 and chicken/Indonesia/Cimanglid/002WJ/2015, from disease outbreaks in chickens in Indonesia are reported. Phylogenetic analysis of different genotypes of Newcastle disease virus using the F gene coding sequences suggests that these two strains belong to genotype VII.2, in class II of avian paramyxoviruses., (Copyright © 2020 Doan et al.)

Here, we report two genomes of newly emerged strains of Newcastle disease virus (NDV), Chicken/Indonesia/Tangerang/004WJ/14 and Chicken/Indonesia/VD/003WJ/11, from disease outbreaks in chickens in Indonesia. Phylogenetic study results of the fusion (F) protein's gene-coding sequences of different genotypes of NDV revealed that these two strains belong to genotype VII.1 in the class II cluster of avian paramyxoviruses., (Copyright © 2020 Rabiei et al.)

The genotype of an organism is stable throughout its life; however, its epigenome is dynamic and can be altered in response to environmental factors, such as diet. Inheritance of acquired epigenetic modifications by the next generation occurs through the germline, although the precise mechanisms remain to be elucidated. Here, we used a sheep model to evaluate if modification of the maternal diet (CTR; control, UND: undernutrition; FA: undernutrition and folic acid supplementation) during the peri-conceptional period affects the genome-wide methylation status of the gametes of male offspring. Sperm DNA methylation, measured by Reduced Representation Bisulfite Sequencing (RRBS), identified Differentially Methylated Regions (DMR) in offspring that experienced in utero undernutrition, both in UND (244) and FA (240), compared with CTR. Gene ontology (GO) analysis identified DMRs in categories related to sperm function, therefore we investigated whether the fertilizing capacity of the semen from the three groups differed in an in vitro fertilization assay. Spermatozoa from the undernourished groups showed lower motility and sperm chromatin structure abnormalities, represented by a higher percentage of DNA fragmentation and an increased number of immature cells, compared with CTR. While good quality blastocysts were obtained from all three groups, the proportion of embryos reaching the blastocyst stage was reduced in the UND vs CTR, an effect partially rescued by the FA treatment. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the semen of the resulting offspring, the effects of which in next generation remain to be elucidated.

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations.

Antimicrobial use in animals and the potential development of antimicrobial resistance is a global concern. So, non-antimicrobial techniques for animal disease control are needed. This study aimed to determine whether neonatal ceftiofur (CF) treatment affects piglet faecal microbiomes and whether faecal microbiome transplantation (FMT) can correct it. Two focal piglets per sow were assigned to treatments as follows: cffresh ( n = 6) received CF (3 mg/kg intramuscular) at 7 d and fresh FMT at 13 d; cffrozen ( n = 7) received CF at 7 d and frozen FMT at 13 d; CF ( n = 8) received CF at 7 d and no FMT; and no CF ( n = 5) received no CF or FMT. DNA was extracted from faecal samples collected on days 7, 13, and 18 for 16S rRNA amplicon analysis. All faecal blends used for the FMT consisted of pooled donor pig faeces at 1:2 ratio with saline, delivered orally at 3 mL/kg. Alpha and beta diversity metrics increased with age (p < 0.05). However, no effect of antibiotic or FMT treatment was evident in 13 and 18 d old piglets ( p > 0.05). Although no effect of treatment was observed, information regarding microbial membership during lactation was gained., Competing Interests: The authors declare no conflict of interest.

In Hong Kong, there is a cattle population of ~1,200 individuals of uncertain origin and genetic diversity. This population shows heterogeneous morphology, both in body type and pigmentation. Once used as draught animals by the local farmers, they were abandoned around the 1970s due to changes in the economy, and since then have lived as feral populations. To explore the origins of these cattle, we analysed ~50k genotype data of 21 Hong Kong feral cattle, along with data from 703 individuals of 36 cattle populations of European, African taurine, and Asian origin, the wild x domestic hybrid gayal, plus two wild bovine species, gaur and banteng. To reduce the effect of ascertainment bias ~4k loci that are polymorphic in the two wild species were selected for further analysis. The stringent SNP selection we applied resulted in increased heterozygosity across all populations studies, compared with the full panel of SNP, thus reducing the impact of ascertainment bias and facilitating the comparison of divergent breeds of cattle. Our results showed that Hong Kong feral cattle have relatively high levels of genetic distinctiveness, possibly due to the low level of artificial selection, and a likely common ancestry with wild species. We found signs of a putative taurine introgression, probably dating to the import of north European breeds during the British colonialism of Hong Kong. We showed that Hong Kong feral cattle, are distinct from Bos taurus and Bos indicus breeds. Our results highlight the distinctiveness of Hong Kong feral cattle and stress the conservation value of this indigenous breed that is likely to harbour adaptive genetic variation, which is a fundamental livestock resource in the face of climate change and diversifying market demands., Competing Interests: The authors have declared that no competing interests exist.