Analysis of immune-related microRNAs in cows and newborn calves.

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.

Long read isoform sequencing reveals hidden transcriptional complexity between cattle subspecies.

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.

Genome-wide association study for antibody response to Mycobacterium avium ssp. paratubeculosis in goats.

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.

Gaur genome reveals expansion of sperm odorant receptors in domesticated cattle.

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.

Weaknesses in Experimental Design and Reporting Decrease the Likelihood of Reproducibility and Generalization of Recent Cardiovascular Research.

Recent evidence indicates that many clinical and preclinical studies are not reproducible. Prominent causes include design and implementation issues, low statistical power, unintentional bias, and incomplete reporting in the published literature. The primary goal of this study was to assess the quality of published research in three prominent cardiovascular research journals by examining statistical power and assessing the adherence to augmented ARRIVE guidelines (Animal Research: Reporting of In Vivo Experiments). For unpaired t-tests, the average median power for a 20% and 50% change was 0.27 ± 0.06 and 0.88 ± 0.08, respectively. For analysis of guidelines, 40 categories were assessed with a 0-2 scale. Although many strengths were observed, several key elements that were needed for reproducibility were inadequate, including differentiation of primary and secondary outcomes, power calculations for group size, allocation methods, use of randomization and blinding, checks for normality, reports of attrition, and adverse events of subjects, and assessment of bias. A secondary goal was to examine whether a required checklist improved the quality of reporting; those results indicated that a checklist improved compliance and quality of reporting, but adequacy levels in key categories were still too low. Overall, the findings of this study indicated that the probability for reproducibility of many clinical and preclinical cardiovascular research studies was low because of incomplete reporting, low statistical power, and lack of research practices that decrease experimental bias. Expansion of group sizes to increase power, use of detailed checklists, and closer monitoring for checklist adherence by editors and journals should remediate many of these deficits and increase the likelihood of reproducibility.

Linkage Disequilibrium and Effective Population Size of Buffalo Populations of Iran, Turkey, Pakistan, and Egypt Using a Medium Density SNP Array.

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.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes.

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.

Near-Infrared Reflectance Spectroscopy for Predicting the Phospholipid Fraction and the Total Fatty Acid Composition of Freeze-Dried Beef.

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

Using Omics Approaches in the Discovery of Biomarkers for Early Diagnosis of Johne's Disease in Sheep and Goats.

Johne's disease (JD) is caused by Mycobacterium avium subsp. paratuberculosis (MAP) and is an important and emerging problem in livestock; therefore, its control and prevention is a priority to reduce economic losses and health risks. Most JD research has been carried out on cattle, but interest in the pathogenesis and diagnosis of this disease in sheep and goats is greatest in developing countries. Sheep and goats are also a relevant part of livestock production in Europe and Australia, and these species provide an excellent resource to study and better understand the mechanism of survival of MAP and gain insights into possible approaches to control this disease. This review gives an overview of the literature on paratuberculosis in sheep and goats, highlighting the immunological aspects and the potential for "omics" approaches to identify effective biomarkers for the early detection of infection. As JD has a long incubation period before the disease becomes evident, early diagnosis is important to control the spread of the disease.

Distinctive gene expression patterns and imprinting signatures revealed in reciprocal crosses between cattle sub-species.

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.

Predicting colostrum and calf blood components based on refractometry.

Provision of good quality colostrum is essential for the passive immunity and nutrition of newborn calves. In order to better predict the quality of colostrum and the transfer of passive immunity, the relationships between colostrum components and between calf serum components were examined in this study. Samples of bulk tank milk, colostrum pooled from several cows 0-4 d postpartum, and colostrum collected from individual cows twice daily for 3 d post-partum were compared. With the exception of fat percentage, there were strong correlations between the levels of the components in the pooled colostrum and in the individual cow colostrum collected 0-1 d postpartum. The correlations between total solids as measured by Brix refractometry and total protein, immunoglobulin G (IgG), lactose % and protein % in colostrum within 1 d postpartum and pooled colostrum were 0.92, 0.90, -0.88 and 0.98, respectively. These high correlations enabled these colostrum components to be accurately predicted from Brix % and therefore, the volume of colostrum required to feed neonate calves can be optimised based on Brix refractometry to avoid failure of passive immunity transfer. To assess whether the components obtained from colostrum were correlated in calf blood, newborn calves were separated from their dams before suckling and blood sampled before feeding (day 0), and on days 1 and 7, after receiving colostrum or milk twice a day. The correlations between glucose, total protein, IgG, and gamma-glutamyl transferase (GGT) levels in the calf blood were lower than the correlations observed between the colostrum components. The highest correlation was between serum protein measured by refractometer and serum IgG within one week postpartum. GGT activity was not a good indicator of serum IgG levels. However, serum protein refractometer measurements predicted serum IgG level with high accuracy, providing an on-farm test to determine that calves have received sufficient passive immunity and colostrum components.

Colostrum source and passive immunity transfer in dairy bull calves.

Colostrum is essential for good neonate health; however, it is not known whether different calves absorb the nutrients from colostrum equally well. In this study, the absorption of protein, IgG, and γ-glutamyl transferase was compared in newborn dairy bull calves for 1 wk after feeding colostrum from different sources. Thirty-five Holstein-Friesian bull calves were randomly allocated into 3 groups and fed colostrum within 4 h after birth. Group A calves (n = 12) were bottle fed colostrum from their own dam for 3 d. Colostrum from these group A cows was also used as foster cow colostrum for the group B calves (n = 12), such that each group A and B calf pair received identical colostrum from each milking of the respective group A dam (10% of birth weight per day). The group C calves (n = 11) were fed 1 bottle (2 L) of pooled colostrum and transition milk (referred to as pooled colostrum), as was the standard practice on the dairy farm. The pooled colostrum was collected from the other dairy cows on the farm 0 to 4 d postpartum and stored at 4°C for less than 12 h. Blood was sampled from calves before the first feeding and at 1, 2, 3, and 7 d after birth. Levels of total solids, total protein, and IgG were higher in the dam colostrum than in the pooled colostrum. At birth, there were no differences between the calf groups for any measurements, and all calves had very low IgG levels. After receiving colostrum, the glucose, plasma γ-glutamyl transferase, serum total protein, and IgG concentrations increased significantly in all calves. There were no differences in any blood measurements at any time point between the pairs of group A and group B calves that received colostrum from the same cow except for the IgG concentration 2 d after birth. However, the group A calves had a higher total serum protein level and IgG concentration than the group C calves for all the time points after the first feeding. The group B calves had a higher IgG concentration than the group C calves on d 1, 2, and 7 after birth. Compared with groups A and B, there was no difference in the proportion of calves in group C that failed to have passive immunity transferred adequately based on the IgG threshold (<10 g/L). Thus, the calves receiving identical colostrum from the same cow had the same levels of IgG, and even the pooled colostrum provided sufficient transfer of IgG as the calves were fed within 4 h after birth.

Depth and strain rate-dependent mechanical response of chondrocytes in reserve zone cartilage subjected to compressive loading.

The role of the growth plate reserve zone is not well understood. It has been proposed to serve as a source of stem cells and to produce morphogens that control the alignment of clones in preparation for the transition into the proliferative zone. We hypothesized that if such a role exists, there are likely to be mechanoregulatory stimuli in cellular response through the depth of the reserve zone. A poroelastic multiscale finite element model of bone/growth-plate/bone was developed for examining the reserve zone cell transient response when compressed to 5% of the cartilage thickness at strain rates of 0.18%/s, 5%/s, 50%/s, and 200%/s. Chondrocyte maximum principal strains, height-, width-, and membrane-strains were found to be highly dependent on reserve zone tissue depth and strain rate. Cell-level strains and fluid transmembrane outflow from the cell were influenced by the permeability of the calcified cartilage between subchondral bone plate and reserve zone and by the applied strain rate. Cell strain levels in the lower reserve zone were less sensitive to epiphyseal permeability than in the upper reserve zone. In contrast, the intracellular fluid pressures were relatively uniform with reserve zone tissue depth and less sensitive to epiphyseal permeability. Fluid shear stress, induced by fluid flow over the cell surface, provided mechanoregulatory signals potentially sufficient to stimulate reserve zone chondrocytes near the subchondral bone plate interface. These results suggest that the strain rate and tissue depth dependence of cell-level strains and cell surface fluid shear stress may provide mechanoregulatory cues in the reserve zone.

Extended haplotype-phasing of long-read de novo genome assemblies using Hi-C.

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.

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events.

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.

Comparison of the response of mammary gland tissue from two divergent lines of goat with high and low milk somatic cell scores to an experimental Staphylococcus aureus infection.

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.

Whole genome analysis of water buffalo and global cattle breeds highlights convergent signatures of domestication.

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.

Haplotype-resolved genomes provide insights into structural variation and gene content in Angus and Brahman cattle.

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.

A genetically unique Chinese cattle population shows evidence of common ancestry with wild species when analysed with a reduced ascertainment bias SNP panel.

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.

Combinatorial Discriminant Analysis Applied to RNAseq Data Reveals a Set of 10 Transcripts as Signatures of Exposure of Cattle to Mycobacterium avium subsp. paratuberculosis.

Paratuberculosis or Johne's disease in cattle is a chronic granulomatous gastroenteritis caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP). Paratuberculosis is not treatable; therefore, the early identification and isolation of infected animals is a key point to reduce its incidence. In this paper, we analyse RNAseq experimental data of 5 ELISA-negative cattle exposed to MAP in a positive herd, compared to 5 negative-unexposed controls. The purpose was to find a small set of differentially expressed genes able to discriminate between exposed animals in a preclinical phase from non-exposed controls. Our results identified 10 transcripts that differentiate between ELISA-negative, clinically healthy, and exposed animals belonging to paratuberculosis-positive herds and negative-unexposed animals. Of the 10 transcripts, five (TRPV4, RIC8B, IL5RA, ERF, CDC40) showed significant differential expression between the three groups while the remaining 5 (RDM1, EPHX1, STAU1, TLE1, ASB8) did not show a significant difference in at least one of the pairwise comparisons. When tested in a larger cohort, these findings may contribute to the development of a new diagnostic test for paratuberculosis based on a gene expression signature. Such a diagnostic tool could allow early interventions to reduce the risk of the infection spreading.