Island populations of large vertebrates have experienced higher extinction rates than mainland populations over long timescales due to demographic stochasticity, genetic drift, and inbreeding. While being more susceptible to extinction and as such potentially targeted for conservation interventions such as genetic rescue, small-island populations can experience relatively less anthropogenic habitat degradation than those on larger islands. Here, we determine the consequences and conservation implications of long-term isolation and recent human activities on genetic diversity of island populations of two forest-dependent mammals endemic to the Wallacea archipelago: the anoa (Bubalus spp.) and babirusa (Babyrousa spp.). Using genomic analyses and habitat suitability models, we show that, compared to closely related species, populations on mainland Sulawesi exhibit low heterozygosity, high inbreeding, a high proportion of deleterious alleles, and experience a high rate of anthropogenic disturbance. In contrast, populations on smaller islands occupy higher-quality habitats, possess fewer deleterious mutations despite exhibiting lower heterozygosity and higher inbreeding. Site frequency spectra indicate that these patterns reflect stronger, long-term purging in smaller-island populations. Our results thus suggest that conservation efforts should focus on protecting small-island high-quality habitats and avoiding translocations from mainland populations. This study highlights the crucial role of small offshore islands for the long-term survival of Wallacea's iconic and indigenous mammals in the face of development on the mainland.

The origin of Arctic charr populations in the lakes of the Italian Alps is not well understood. While some studies have suggested that they are postglacial relics, others have proposed that these populations are the result of intentional stocking efforts dating back to the sixteenth century. Subsequent introductions of Arctic charr to these lakes have made it difficult to untangle the evolutionary history of these Alpine populations. In this study, we examined the distribution of genetic variation among Arctic charr collected throughout their Northern and Southern European ranges at the beginning of the twenty-first century, using nuclear (amplified fragment length polymorphisms) and mitochondrial (control region and cytochrome oxidase I) loci. These analyses revealed the impact of restocking activities, which have resulted in admixture and hybridisation with Brook charr (Salvelinus fontinalis) in Italian alpine lakes and provides clues on the native or pseudo-native origin of Trentino-Alto Adige populations. The lack of detailed historical information, however, makes it difficult to disentangle the postglacial history of the species and to determine how much of the current diversity pattern can be ascribed to the consequences of Pleistocene events or anthropogenic activities.

Sensory beef quality is an important parameter for consumers, but it is difficult to understand what determines quality because of the large variation due to breed, type of animal, sex, and farm management. To address these challenges, meat samples from Longissimus thoracis (LT) were collected 24 h after slaughter from a total of 436 young bulls representing 15 cattle breeds. The samples were analysed for physiological, muscular, and sensory characteristics to evaluate beef quality. Ten groups of variables were identified by hierarchical cluster analysis of individual variables and Principal Component Analysis (PCA). The variables described distinct characteristics, including Physiological traits (animal maturity, growth rate, muscle mass), Sensory traits (tenderness, juiciness, flavor, meat colour), and some Muscular characteristics such as lipid content, maturation enzymes and oxidative metabolism, which are likely to discriminate among breeds. The analyses showed that muscular characteristics are positively associated with beef sensory characteristics whereas physiological characteristics are negatively associated with sensory quality. Dairy breeds produced beef rich in lipids with a strong flavour, French breeds showed fast growth rate and highly flavoured meat, and Italian breeds were characterised by good muscular development. In contrast, British breeds have oxidative muscles and produce beef with a strong flavour. Danish breeds have an intermediate score for the sensorial and muscular characteristics.

In order to identify the selection signatures of the water buffalos across the world, the genomic information of 165 buffalos which belonged to 15 genetic groups of buffaloes was used. The genomic information was obtained from Dryad (doi:10.5061/dryad.h0cc7). The quality control and data filtration were performed using PLINK1.9 software. The genetic clustering and the population structure was examined using the GenABEL and Admixture1.23 software's, respectively. The results of principal component analysis showed that the examined populations could be classified into 4 separate categories. The results of population structure analysis confirmed the results of principal components analysis. The signatures of selection were searched with the help of iHS statistics using the ReHH software. Moreover, the unbiased FST (θ) estimator was calculated using the Plink1.9 software. The 25 and 24 genomic regions, which passed the unbiased FST and iHS statistics thresholds, were identified as selection cues, respectively. Selected regions were aligned on the bovine genome and 411 genes related to selected regions were identified. Of all the identified genes, 53 genes related to olfactory receptors (OR), 51 genes somehow involved in cell membrane structure and animal immunity against pathogens including initiate and regulate the immune response. The identified QTLs related to detected regions, were associated with milk production, milk somatic cells, fertility, ion disease, calving and growth. There is an acceptable consistent between the milk and fat production genes and the related identified QTLs.

The presence of horns in domestic ruminants, such as cattle, sheep and goats, has financial and welfare implications. The genetic interactions that lead to horn development are not known. Hornless, or polled, cattle occur naturally. The known causative DNA variants (Celtic, Friesian, Mongolian and Guarani) are in intergenic regions on bovine chromosome 1, but their functions are not known. It is thought that horns may be derived from cranial neural crest stem cells and the POLLED variants disrupt the migration or proliferation of these cells. Relaxin family peptide receptor 2 (RXFP2) is more highly expressed in developing horns in cattle compared to nearby skin and has been shown to play a role in horn development in sheep. However, the role of RXFP2 in horn formation is not understood. Histological analyses of cranial tissues from homozygous horned and polled cattle fetuses at day 58 of development was carried out to determine the differences in the structure of the horn bud region. Condensed cells were only observed in the horn bud mesenchyme of horned fetuses and could be the progenitor horn cells. The distribution of neural crest markers (SOX10 and NGFR) and RXFP2 between horned and polled tissues by immunohistochemistry was also analysed. However, SOX10 and NGFR were not detected in the condensed cells, and therefore, these cells are either not derived from the neural crest, or have differentiated and no longer express neural crest markers. SOX10 and NGFR were detected in the peripheral nerves, while RXFP2 was detected in peripheral nerves and in the horn bud epidermis. Previous research has shown that RXFP2 variants are associated with horn phenotypes in cattle an sheep. Therefore, the RXFP2 variants may affect the development of the epidermis or peripheral nerves in the horn bud.

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_UOA_SB_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.

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.

Breeding animals able to effectively respond to stress could be a long-term, sustainable, and affordable strategy to improve resilience and welfare in livestock systems. In the present study, the concentrations of 29 plasma biomarkers were used as candidate endophenotypes for metabolic stress response in single-SNP, gene- and haplotype-based GWAS using 739 healthy lactating Italian Holstein cows and 88,271 variants. Significant genetic associations were found in all the 3 GWAS approaches for plasma gamma-glutamyl transferase concentration on BTA17, for paraoxonase on BTA4, and for alkaline phosphatase and zinc on BTA2. On these chromosomes, single-SNP and gene-based chromosome-wide association studies were performed, confirming GWAS findings. The signals identified for paraoxonase, gamma-glutamyl transferase, and alkaline phosphatase were in proximity to the genes coding for them. The heritability of these 4 biomarkers ranged from moderate to high (from 0.39 to 0.54). Plasma biomarkers are known to undergo large changes in concentration during metabolic stress in the transition period, with an interindividual variability in the rate of change and recovery time. Genetics may account in part for these differences. To assess this, we studied a subset of 139 periparturient cows homozygous at 3 SNPs known to be respectively associated with concentration of plasma ceruloplasmin, paraoxonase, and gamma-glutamyl transferase. We compared the immune-metabolic profile measured in plasma at -7, +5, and +30 d relative to calving between groups of opposite homozygotes. A significant effect of the genotype was found on paraoxonase and gamma-glutamyl transferase plasma concentration at all the 3 time points. No evidence for genotype effect was detected for ceruloplasmin. Understanding the genetic control underlying metabolic stress response may suggest new approaches to foster resilience in dairy cows.

Sperm motility is directly related to the ability of sperm to move through the female reproductive tract to reach the ovum. Sperm motility is a complex trait that is influenced by environmental and genetic factors and is associated with male fertility, oocyte penetration rate, and reproductive success of cattle. In this study we carried out a GWAS in Italian Holstein bulls to identify candidate regions and genes associated with variations in progressive and total motility (PM and TM, respectively). After quality control, the final data set consisted of 5,960 records from 949 bulls having semen collected in 10 artificial insemination stations and genotyped at 412,737 SNPs (call rate >95%; minor allele frequency >5%). (Co)variance components were estimated using single trait mixed models, and associations between SNPs and phenotypes were assessed using a genomic BLUP approach. Ten windows that explained the greatest percentage of genetic variance were located on Bos taurus autosomes 1, 2, 4, 6, 7, 23, and 26 for TM and Bos taurus autosomes 1, 2, 4, 6, 8, 16, 23, and 26 for PM. A total of 150 genes for TM and 72 genes for PM were identified within these genomic regions. Gene Ontology enrichment analyses identified significant Gene Ontology terms involved with energy homeostasis, membrane functions, sperm-egg interactions, protection against oxidative stress, olfactory receptors, and immune system. There was significant enrichment of quantitative trait loci for fertility, calving ease, immune response, feed intake, and carcass weight within the candidate windows. These results contribute to understanding the architecture of the genetic control of sperm motility and may aid in the development of strategies to identify subfertile bulls and improve reproductive success.

Bovine colostrum contains a high concentration of immune-related microRNAs (miRNAs) that are packaged in exosomes and are very stable. In this study, 5 immune-related miRNAs (miR-142-5p, miR-150, miR-155, miR-181a, and miR-223) were quantified in dam blood, colostrum, and calf blood using reverse transcription quantitative PCR. Their levels in calf blood after colostrum ingestion were investigated to assess whether miRNAs are transferred from the dam to newborn calves. Three groups of Holstein-Friesian bull calves were bottle-fed 2 L of colostrum or milk from different sources twice per day. The group A calves received colostrum from their own dam and the group B calves were fed foster dam colostrum. Each pair of group A and group B calves were fed identical colostrum from the same milking of the corresponding group A dam for 3 d and then bulk tank milk for 7 d after birth. Group C calves were fed only 2L of "pooled colostrum" from multiple dams d 0 to 4 postpartum, and then fed bulk tank milk thereafter for 7 d after birth. The groups were fed colostrum from different sources and different amounts to assess possible miRNA absorption from the colostrum. All miRNAs were at the highest level in colostrum at d 0 and then decreased rapidly after d 1. The level of miR-150 had the largest decrease from 489 × 106 copies/µL (d 0) to 78 × 106 copies/µL (d 1). MicroRNA-223 and miR-155 were the most abundant in both colostrum and milk. Dam colostrum had significantly higher levels of miR-142-5p, miR-155, and miR-181a than the bulk tank milk. However, only the miR-155 concentration was significantly higher in the dam colostrum than in the pooled colostrum. The concentrations of miRNAs in the colostrum were less than in the cow blood (100- to 1,000-fold less). There was no significant correlation between the level of miRNAs in the dam blood and their colostrum, suggesting that miRNA is synthesized locally by the mammary gland rather than being transferred from the blood. MicroRNA-223 had the highest level in both calf and cow blood compared with the other 4 immune-related miRNAs. Calves were born with high levels of immune-related miRNAs in their blood, and there were no significant differences in miRNA levels between the 3 calf groups at birth or after they were fed different colostrum. This suggests that these miRNAs were not transferred from the colostrum to the newborn calves.

Horns, a form of headgear carried by Bovidae, have ethical and economic implications for ruminant production species such as cattle and goats. Hornless (polled) individuals are preferred. In cattle, four genetic variants (Celtic, Friesian, Mongolian and Guarani) are associated with the polled phenotype, which are clustered in a 300-kb region on chromosome 1. As the variants are intergenic, the functional effect is unknown. The aim of this study was to determine if the POLLED variants affect chromatin structure or disrupt enhancers using publicly available data. Topologically associating domains (TADs) were analyzed using Angus- and Brahman-specific Hi-C reads from lung tissue of an Angus (Celtic allele) cross Brahman (horned) fetus. Predicted bovine enhancers and chromatin immunoprecipitation sequencing peaks for histone modifications associated with enhancers (H3K27ac and H3K4me1) were mapped to the POLLED region. TADs analyzed from Angus- and Brahman-specific Hi-C reads were the same, therefore, the Celtic variant does not appear to affect this level of chromatin structure. The Celtic variant is located in a different TAD from the Friesian, Mongolian, and Guarani variants. Predicted enhancers and histone modifications overlapped with the Guarani and Friesian variants but not the Celtic or Mongolian variants. This study provides insight into the mechanisms of the POLLED variants for disrupting horn development. These results should be validated using data produced from the horn bud region of horned and polled bovine fetuses.

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.

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.

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.

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 peri-conceptional period (42 days in total: -14/+28 from mating), would impact CpG methylation in muscle tissue (Longissimus dorsi) of adult offspring (11.5 months old). Reduced Representation Bisulfite Sequencing (RRBS), 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 (GO) 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.

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.

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.

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.