Unveiling the Role of Endoplasmic Reticulum Stress Pathways in Canine Demodicosis.

Canine demodicosis is a prevalent skin disease caused by overpopulation of a commensal species of Demodex mite, yet its precise cause remains unknown. Research suggests that T-cell exhaustion, increased immunosuppressive cytokines, induction of regulatory T cells and increased expression of immune checkpoint inhibitors may contribute to its pathogenesis. This study aimed to gain a deeper understanding of the molecular changes occurring in canine demodicosis using mass spectrometry and pathway enrichment analysis. The results indicate that endoplasmic reticulum stress promotes canine demodicosis through regulation of three linked signalling pathways: eIF2, mTOR, and eIF4 and p70S6K. These pathways are involved in the modulation of Toll-like receptors, most notably TLR2, and have been shown to play a role in the pathogenesis of skin diseases in both dogs and humans. Moreover, these pathways are also implicated in the promotion of immunosuppressive M2 phenotype macrophages. Immunohistochemical analysis, utilising common markers of dendritic cells and macrophages, verified the presence of M2 macrophages in canine demodicosis. The proteomic analysis also identified immunological disease, organismal injury and abnormalities and inflammatory response as the most significant underlying diseases and disorders associated with canine demodicosis. This study demonstrates that Demodex mites, through ER stress, unfolded protein response and M2 macrophages contribute to an immunosuppressive microenvironment, thereby assisting in their proliferation.

Selection signatures for local and regional adaptation in Chinese Mongolian horse breeds reveal candidate genes for hoof health.

BACKGROUND: Thousands of years of natural and artificial selection since the domestication of the horse has shaped the distinctive genomes of Chinese Mongolian horse populations. Consequently, genomic signatures of selection can provide insights into the human-mediated selection history of specific traits and evolutionary adaptation to diverse environments. Here, we used genome-wide SNPs from five distinct Chinese Mongolian horse populations to identify genomic regions under selection for the population-specific traits, gait, black coat colour, and hoof quality. Other global breeds were used to identify regional-specific signatures of selection. RESULTS: We first identified the most significant selection peak for the Wushen horse in the region on ECA23 harbouring DMRT3, the major gene for gait. We detected selection signatures encompassing several genes in the Baicha Iron Hoof horse that represent good biological candidates for hoof health, including the CSPG4, PEAK1, EXPH5, WWP2 and HAS3 genes. In addition, an analysis of regional subgroups (Asian compared to European) identified a single locus on ECA3 containing the ZFPM1 gene that is a marker of selection for the major domestication event leading to the DOM2 horse clade. CONCLUSIONS: Genomic variation at these loci in the Baicha Iron Hoof may be leveraged in other horse populations to identify animals with superior hoof health or those at risk of hoof-related pathologies. The overlap between the selection signature in Asian horses with the DOM2 selection peak raises questions about the nature of horse domestication events, which may have involved a prehistoric clade other than DOM2 that has not yet been identified.

Inbreeding depression and durability in the North American Thoroughbred horse.

The proportion of the genome containing runs of homozygosity (ROH) affects production traits in livestock populations. In European and Australasian Thoroughbreds inbreeding, quantified using ROH (FROH ), is associated with the probability of ever racing. Here, we measured FROH using 333 K SNP genotypes from 768 Thoroughbred horses born in North America to evaluate the effect of inbreeding on racing traits in that region. Among North American horses, FROH was not associated (p = 0.518) with the probability of ever racing but was significantly associated with the number of race starts (p = 0.002). Among raced horses, those with a 10% higher FROH than the mean inbreeding coefficient were predicted to have 3.5 fewer race starts compared to horses with a mean inbreeding coefficient. Considering the trend of increasing inbreeding and a decline in the average number of race starts per runner in North America, mitigating inbreeding in the population could positively influence racing durability.

Integrative genomics analysis highlights functionally relevant genes for equine behaviour.

Behavioural plasticity enables horses entering an exercise training programme to adapt with reduced stress. We characterised SNPs associated with behaviour in yearling Thoroughbred horses using genomics analyses for two phenotypes: (1) handler-assessed coping with early training events [coping] (n = 96); and (2) variation in salivary cortisol concentration at the first backing event [cortisol] (n = 34). Using RNA-seq derived gene expression data for amygdala and hippocampus tissues from n = 2 Thoroughbred stallions, we refined the SNPs to those with functional relevance to behaviour by cross-referencing to the 500 most highly expressed genes in each tissue. The SNPs of high significance (q < 0.01) were in proximity to genes (coping - GABARAP, NDM, OAZ1, RPS15A, SPARCL1, VAMP2; cortisol - CEBPA, COA3, DUSP1, HNRNPH1, RACK1) with biological functions in social behaviour, autism spectrum disorder, suicide, stress-induced anxiety and depression, Alzheimer's disease, neurodevelopmental disorders, neuroinflammatory disease, fear-induced behaviours and alcohol and cocaine addiction. The strongest association (q = 0.0002) was with NDN, a gene previously associated with temperament in cattle. This approach highlights functionally relevant genes in the behavioural adaptation of Thoroughbred horses that will contribute to the development of genetic markers to improve racehorse welfare.

Inbreeding depression and the probability of racing in the Thoroughbred horse.

Small effective population sizes and active inbreeding can lead to inbreeding depression due to deleterious recessive mutations exposed in the homozygous state. The Thoroughbred racehorse has low levels of population genetic diversity, but the effects of genomic inbreeding in the population are unknown. Here, we quantified inbreeding based on runs of homozygosity (ROH) using 297 K SNP genotypes from 6128 horses born in Europe and Australia, of which 13.2% were unraced. We show that a 10% increase in inbreeding (FROH) is associated with a 7% lower probability of ever racing. Moreover, a ROH-based genome-wide association study identified a haplotype on ECA14 which, in its homozygous state, is linked to a 32.1% lower predicted probability of ever racing, independent of FROH. The haplotype overlaps a candidate gene, EFNA5, that is highly expressed in cartilage tissue, which when damaged is one of the most common causes of catastrophic musculoskeletal injury in racehorses. Genomics-informed breeding aiming to reduce inbreeding depression and avoid damaging haplotype carrier matings will improve population health and racehorse welfare.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria.

BACKGROUND: Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection. RESULTS: The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool-gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively. CONCLUSIONS: The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics.

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

Refinement of Global Domestic Horse Biogeography Using Historic Landrace Chinese Mongolian Populations.

The Mongolian horse is one of the oldest extant horse populations and although domesticated, most animals are free-ranging and experience minimal human intervention. As an ancient population originating in one of the key domestication centers, the Mongolian horse may play a key role in understanding the origins and recent evolutionary history of horses. Here we describe an analysis of high-density genome-wide single-nucleotide polymorphism (SNP) data in 40 globally dispersed horse populations (n = 895). In particular, we have focused on new results from Chinese Mongolian horses (n = 100) that represent 5 distinct populations. These animals were genotyped for 670K SNPs and the data were analyzed in conjunction with 35K SNP data for 35 distinct breeds. Analyses of these integrated SNP data sets demonstrated that the Chinese Mongolian populations were genetically distinct from other modern horse populations. In addition, compared to other domestic horse breeds, the Chinese Mongolian horse populations exhibited relatively high genomic diversity. These results suggest that, in genetic terms, extant Chinese Mongolian horses may be the most similar modern populations to the animals originally domesticated in this region of Asia. Chinese Mongolian horse populations may therefore retain ancestral genetic variants from the earliest domesticates. Further genomic characterization of these populations in conjunction with archaeogenetic sequence data should be prioritized for understanding recent horse evolution and the domestication process that has led to the wealth of diversity observed in modern global horse breeds.

Equine skeletal muscle adaptations to exercise and training: evidence of differential regulation of autophagosomal and mitochondrial components.

BACKGROUND: A single bout of exercise induces changes in gene expression in skeletal muscle. Regular exercise results in an adaptive response involving changes in muscle architecture and biochemistry, and is an effective way to manage and prevent common human diseases such as obesity, cardiovascular disorders and type II diabetes. However, the biomolecular mechanisms underlying such responses still need to be fully elucidated. Here we performed a transcriptome-wide analysis of skeletal muscle tissue in a large cohort of untrained Thoroughbred horses (n = 51) before and after a bout of high-intensity exercise and again after an extended period of training. We hypothesized that regular high-intensity exercise training primes the transcriptome for the demands of high-intensity exercise. RESULTS: An extensive set of genes was observed to be significantly differentially regulated in response to a single bout of high-intensity exercise in the untrained cohort (3241 genes) and following multiple bouts of high-intensity exercise training over a six-month period (3405 genes). Approximately one-third of these genes (1025) and several biological processes related to energy metabolism were common to both the exercise and training responses. We then developed a novel network-based computational analysis pipeline to test the hypothesis that these transcriptional changes also influence the contextual molecular interactome and its dynamics in response to exercise and training. The contextual network analysis identified several important hub genes, including the autophagosomal-related gene GABARAPL1, and dynamic functional modules, including those enriched for mitochondrial respiratory chain complexes I and V, that were differentially regulated and had their putative interactions 're-wired' in the exercise and/or training responses. CONCLUSION: Here we have generated for the first time, a comprehensive set of genes that are differentially expressed in Thoroughbred skeletal muscle in response to both exercise and training. These data indicate that consecutive bouts of high-intensity exercise result in a priming of the skeletal muscle transcriptome for the demands of the next exercise bout. Furthermore, this may also lead to an extensive 're-wiring' of the molecular interactome in both exercise and training and include key genes and functional modules related to autophagy and the mitochondrion.

Differential gene expression in the endometrium reveals cytoskeletal and immunological genes in lactating dairy cows genetically divergent for fertility traits.

Profitable milk production in dairy cows requires good reproductive performance. Calving interval is a trait used to measure reproductive efficiency. Herein we used a novel lactating Holstein cow model of fertility that displayed genetic and phenotypic divergence in calving interval, a trait used to define reproductive performance using a national breeding index in Ireland. Cows had similar genetic merit for milk production traits, but either very good genetic merit for fertility (Fert+; n=7) or very poor genetic merit for fertility (Fert-; n=6). We tested the hypothesis that Fert+ cows would have a corresponding detectable difference in endometrial gene expression compared with the Fert- cows. To do this, we sequenced the transcriptome of endometrial biopsies collected on Day 7 of the oestrous cycle (non-pregnant). This is an important stage for uterine remodelling and initiation of histotroph secretion. Significant differential expression (false discovery rate-adjusted P<0.1) of 403 genes between Fert+ and Fert- cows was found. A novel network-based functional analysis highlighted 123 genes from three physiologically relevant networks of the endometrium: (1) actin and cytoskeletal components; (2) immune function; and (3) ion transportation. In particular, our results indicate an overall downregulation of inflammation-related genes and an upregulation of multiple ion transporters and gated-voltage channels and cytoskeletal genes in Fert+ cows. These three topics, which are discussed in terms of the uterus and in the context of fertility, provide molecular evidence for an association between gene expression in the uterine environment and genetic merit for fertility in dairy cows.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

GOexpress: an R/Bioconductor package for the identification and visualisation of robust gene ontology signatures through supervised learning of gene expression data.

BACKGROUND: Identification of gene expression profiles that differentiate experimental groups is critical for discovery and analysis of key molecular pathways and also for selection of robust diagnostic or prognostic biomarkers. While integration of differential expression statistics has been used to refine gene set enrichment analyses, such approaches are typically limited to single gene lists resulting from simple two-group comparisons or time-series analyses. In contrast, functional class scoring and machine learning approaches provide powerful alternative methods to leverage molecular measurements for pathway analyses, and to compare continuous and multi-level categorical factors. RESULTS: We introduce GOexpress, a software package for scoring and summarising the capacity of gene ontology features to simultaneously classify samples from multiple experimental groups. GOexpress integrates normalised gene expression data (e.g., from microarray and RNA-seq experiments) and phenotypic information of individual samples with gene ontology annotations to derive a ranking of genes and gene ontology terms using a supervised learning approach. The default random forest algorithm allows interactions between all experimental factors, and competitive scoring of expressed genes to evaluate their relative importance in classifying predefined groups of samples. CONCLUSIONS: GOexpress enables rapid identification and visualisation of ontology-related gene panels that robustly classify groups of samples and supports both categorical (e.g., infection status, treatment) and continuous (e.g., time-series, drug concentrations) experimental factors. The use of standard Bioconductor extension packages and publicly available gene ontology annotations facilitates straightforward integration of GOexpress within existing computational biology pipelines.

Definition of the cattle killer cell Ig-like receptor gene family: comparison with aurochs and human counterparts.

Under selection pressure from pathogens, variable NK cell receptors that recognize polymorphic MHC class I evolved convergently in different species of placental mammal. Unexpectedly, diversified killer cell Ig-like receptors (KIRs) are shared by simian primates, including humans, and cattle, but not by other species. Whereas much is known of human KIR genetics and genomics, knowledge of cattle KIR is limited to nine cDNA sequences. To facilitate comparison of the cattle and human KIR gene families, we determined the genomic location, structure, and sequence of two cattle KIR haplotypes and defined KIR sequences of aurochs, the extinct wild ancestor of domestic cattle. Larger than its human counterpart, the cattle KIR locus evolved through successive duplications of a block containing ancestral KIR3DL and KIR3DX genes that existed before placental mammals. Comparison of two cattle KIR haplotypes and aurochs KIR show the KIR are polymorphic and the gene organization and content appear conserved. Of 18 genes, 8 are functional and 10 were inactivated by point mutation. Selective inactivation of KIR3DL and activating receptor genes leaves a functional cohort of one inhibitory KIR3DL, one activating KIR3DX, and six inhibitory KIR3DX. Functional KIR diversity evolved from KIR3DX in cattle and from KIR3DL in simian primates. Although independently evolved, cattle and human KIR gene families share important function-related properties, indicating that cattle KIR are NK cell receptors for cattle MHC class I. Combinations of KIR and MHC class I are the major genetic factors associated with human disease and merit investigation in cattle.

High-throughput transcriptomics reveals common and strain-specific responses of human macrophages to infection with Mycobacterium abscessus Smooth and Rough variants.

BACKGROUND: Mycobacterium abscessus (MAB) is an emerging pathogen causing pulmonary infections in those with inflammatory lung disorders, such as Cystic Fibrosis (CF), and is associated with the highest fatality rate among rapidly growing mycobacteria (RGM). Phenotypically, MAB manifests as either a Smooth (MAB-S) or a Rough (MAB-R) morphotype, which differ in their levels of cell wall glycopeptidolipids (GPLs) and in their pathogenicity in vivo. As one of the primary immune cells encountered by MAB, we sought to examine the early transcriptional events within macrophages, following infection with both MAB-S or MAB-R. RESULTS: We sampled the transcriptomes (mRNA and miRNA) of THP-1-derived macrophages infected with both MAB-R and MAB-S at 1, 4 and 24 h post-infection (hpi) using RNA-seq. A core set of 606 genes showed consistent expression profiles in response to both morphotypes, corresponding to the early transcriptional response to MAB. The core response is type I Interferon (IFN)-driven, involving the NF-κB and MAPK signaling pathways with concomitant pro-inflammatory cytokine production, and network analysis identified STAT1, EGR1, and SRC as key hub and bottleneck genes. MAB-S elicited a more robust transcriptional response at both the mRNA and miRNA levels, which was reflected in higher cytokine levels in culture supernatants. The transcriptional profiles of macrophages infected with both morphotypes were highly correlated, however, and a direct comparison identified few genes to distinguish them. Most of the induced miRNAs have previously been associated with mycobacterial infection and overall miRNA expression patterns were similarly highly correlated between the morphotypes. CONCLUSIONS: The report here details the first whole transcriptome analysis of the early macrophage response to MAB infection. The overall picture at the early stages of macrophage infection is similar to that of other mycobacteria, reflected in a core type I IFN and pro-inflammatory cytokine response. Large-scale transcriptional differences in the host response to the different MAB morphotypes are not evident in the early stages of infection, however the subset of genes with distinct expression profiles suggest potentially interesting differences in internal trafficking of MAB within macrophages.

Receptor-mediated recognition of mycobacterial pathogens.

Mycobacteria are a genus of bacteria that range from the non-pathogenic Mycobacterium smegmatis to Mycobacterium tuberculosis, the causative agent of tuberculosis in humans. Mycobacteria primarily infect host tissues through inhalation or ingestion. They are phagocytosed by host macrophages and dendritic cells. Here, conserved pathogen-associated molecular patterns (PAMPs) on the surface of mycobacteria are recognized by phagocytic pattern recognition receptors (PRRs). Several families of PRRs have been shown to non-opsonically recognize mycobacterial PAMPs, including membrane-bound C-type lectin receptors, membrane-bound and cytosolic Toll-like receptors and cytosolic NOD-like receptors. Recently, a possible role for intracellular cytosolic PRRs in the recognition of mycobacterial pathogens has been proposed. Here, we discuss currentideas on receptor-mediated recognition of mycobacterial pathogens by macrophages and dendritic cells.

Transcriptomics analysis of the bovine endometrium during the perioestrus period.

During the oestrous cycle, the bovine endometrium undergoes morphological and functional changes, which are regulated by alterations in the levels of oestrogen and progesterone and consequent changes in gene expression. To clarify these changes before and after oestrus, RNA-seq was used to profile the transcriptome of oestrus-synchronized beef heifers. Endometrial samples were collected from 29 animals, which were slaughtered in six groups beginning 12 h after the withdrawal of intravaginal progesterone releasing devices until seven days post-oestrus onset (luteal phase). The groups represented proestrus, early oestrus, metoestrus and early dioestrus (luteal phase). Changes in gene expression were estimated relative to gene expression at oestrus. Ingenuity Pathway Analysis (IPA) was used to identify canonical pathways and functional processes of biological importance. A total of 5,845 differentially expressed genes (DEGs) were identified. The lowest number of DEGs was observed at the 12 h post-oestrus time point, whereas the greatest number was observed at Day 7 post-oestrus onset (luteal phase). A total of 2,748 DEGs at this time point did not overlap with any other time points. Prior to oestrus, Neurological disease and Organismal injury and abnormalities appeared among the top IPA diseases and functions categories, with upregulation of genes involved in neurogenesis. Lipid metabolism was upregulated before oestrus and downregulated at 48h post-oestrus, at which point an upregulation of immune-related pathways was observed. In contrast, in the luteal phase the Lipid metabolism and Small molecule biochemistry pathways were upregulated.

Functional analysis of the Mycobacterium bovis AF2122/97 PhoPR system.

The PhoPR system is a master regulator in Mycobacterium tuberculosis. A key difference between M. tuberculosis and Mycobacterium bovis is a G71I substitution in the M. bovis PhoR orthologue. Functional studies of the M. bovis PhoPR system have generated conflicting findings, with some research suggesting that the M. bovis PhoPR is defective while others indicate it is functional. We sought to revisit the functionality of the M. bovis PhoPR system. To address this, we constructed a phoPR mutant in the reference strain M. bovis AF2122/97. We employed a combination of growth assays and transcriptomics analyses to assess the phenotype of the mutant vs wild type and complemented strains. We found that the M. bovis AF2122/97 ΔphoPR mutant showed a growth defect on solid and liquid media compared to the wild type and complemented strains. The transcriptome of the M. bovis AF2122/97 ΔphoPR mutant was also altered as compared to wild type, including differential expression of genes involved in lipid metabolism and secretion. Our work provides further insight into the activity of PhoPR in M. bovis and underlines the importance of the PhoPR system as a master regulator of gene expression in the Mycobacterium tuberculosis complex.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle.

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Recent selection and introgression facilitated high-altitude adaptation in cattle.

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro.

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.