Fine mapping of loci on BTA8 associated to antibody response to Mycobacterium avium paratuberculosis in cattle

Paratuberculosis (ParaTB) or Johne's disease, caused by Mycobacterium avium subspecies paratuberculosis commonly known as MAP in cattle, is a chronic gastroenteritis characterized by diarrhoea, decreased milk production and ultimately death. MAP is responsible for huge economic losses, particularly in dairy cattle herds. Susceptibility to MAP infection has been found to be heritable with heritability estimates ranging from 0.06 to 0.102. The definition of an infected animal can be based either on the presence of anti-MAP antibodies in the serum, or by direct demonstration of MAP in tissue or faeces by culture or PCR. Several studies have addressed the identification of genetic loci associated with MAP susceptibility. The objective of this study was to refine a locus associated with antibody response to Mycobacterium avium paratuberculosi (MAP). Using a genome- wide association analysis, a single nucleotide polymorphism on Bos taurus autosome BTA8 namely the SNP rs43161947 at posi- tion 35398490 with a p-value of 7.02 e-05, has previously been identified by the authors as associated with MAP infection. Fine mapping of the region was conducted with 100 single nucleotide polymorphisms spanning a region between BTA8: 34422912 and BTA8: 364553881 covering 2 Mega bases (Mb) designed in to cover 1 Mb ahead and after the SNP identified on BTA8. The 2 Mb region on BTA8 was evaluated within a group of 966 Holstein cows collected from routine ParaTB screening in the province of Lodi in Italy, in an area with a high prevalence of ParaTB. Animals were defined as ParaTB positive based on the detection of serum antibodies produced in response to MAP infection using the ID-screen\uae ELISA test (ID VET Montpellier, France). Of the 966 samples, 483 were MAP antibody positive (cases) and 483 MAP antibody negative (MAP negative controls). All animals were female, and cases and MAP negative controls were from the same farm tested on the same day.Using a single marker associ- ation analysis, conducted within the R statistical environment, we identified 3 different QTLs within the 2 Mega base region, under the main QTL on BTA8 associated with antibody response to MAP, in position 34.700.000, 35.800.000 and 36.400.000 bp. This reveals the complexity of the genetic architecture of thetrait and confirms the need to further explore the genome with fine mapping approaches, or by the use of whole genome sequencing to investigate complex traits, such as disease resistance

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range

Adaptive introgression from indicine cattle into white cattle breeds from Central Italy

Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies

Insights in the genome of mycobacterium avium subsp. paratuberculosis by Next Generation Sequencing approaches

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of paratuberculosis - or Johne's disease \u2013 that affects farmed and wild animals worldwide, causing negative economic consequences particularly relevant in the livestock sector of dairy cattle and beef (1). Recent estimates say that more than 50 % of the herds in Europe and North America are infected (3). In Italy, a study conducted in the Lombardy and Veneto regions reveals that about 70 % of dairy herds are infected (2). The disease shows high variability in the progression and symptoms that may be due to the genetic variability of the host, the pathogen, or a combination of the two. Understanding the mechanism responsible of this variability could be of paramount importance for the control of the disease (1, 3). Aim of this work was to study the genomic variability of MAP isolated from dairy cattle from different farms distributed in several Italian regions through the use of Next Generation Sequencing (NGS) techniques. The preliminary results on 15 strains are presented

Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes

Background: Goat populations that are characterized within the AdaptMap project cover a large part of the worldwide distribution of this species and provide the opportunity to assess their diversity at a global scale. We analysed genome-wide 50 K single nucleotide polymorphism (SNP) data from 144 populations to describe the global patterns of molecular variation, compare them to those observed in other livestock species, and identify the drivers that led to the current distribution of goats. Results: A high degree of genetic variability exists among the goat populations studied. Our results highlight a strong partitioning of molecular diversity between and within continents. Three major gene pools correspond to goats from Europe, Africa and West Asia. Dissection of sub-structures disclosed regional gene pools, which reflect the main post-domestication migration routes. We also identified several exchanges, mainly in African populations, and which often involve admixed and cosmopolitan breeds. Extensive gene flow has taken place within specific areas (e.g., south Europe, Morocco and Mali-Burkina Faso-Nigeria), whereas elsewhere isolation due to geographical barriers (e.g., seas or mountains) or human management has decreased local gene flows. Conclusions: After domestication in the Fertile Crescent in the early Neolithic era (ca. 12,000 YBP), domestic goats that already carried differentiated gene pools spread to Europe, Africa and Asia. The spread of these populations determined the major genomic background of the continental populations, which currently have a more marked subdivision than that observed in other ruminant livestock species. Subsequently, further diversification occurred at the regional level due to geographical and reproductive isolation, which was accompanied by additional migrations and/or importations, the traces of which are still detectable today. The effects of breed formation were clearly detected, particularly in Central and North Europe. Overall, our results highlight a remarkable diversity that occurs at the global scale and is locally partitioned and often affected by introgression from cosmopolitan breeds. These findings support the importance of long-term preservation of goat diversity, and provide a useful framework for investigating adaptive introgression, directing genetic improvement and choosing breeding targets

VarGoats project: a dataset of 1159 whole-genome sequences to dissect Capra hircus global diversity

Background: Since their domestication 10,500 years ago, goat populations with distinctive genetic backgrounds have adapted to a broad variety of environments and breeding conditions. The VarGoats project is an international 1000-genome resequencing program designed to understand the consequences of domestication and breeding on the genetic diversity of domestic goats and to elucidate how speciation and hybridization have modeled the genomes of a set of species representative of the genus Capra. Findings: A dataset comprising 652 sequenced goats and 507 public goat sequences, including 35 animals representing eight wild species, has been collected worldwide. We identified 74,274,427 single nucleotide polymorphisms (SNPs) and 13,607,850 insertion-deletions (InDels) by aligning these sequences to the latest version of the goat reference genome (ARS1). A Neighbor-joining tree based on Reynolds genetic distances showed that goats from Africa, Asia and Europe tend to group into independent clusters. Because goat breeds from Oceania and Caribbean (Creole) all derive from imported animals, they are distributed along the tree according to their ancestral geographic origin. Conclusions: We report on an unprecedented international effort to characterize the genome-wide diversity of domestic goats. This large range of sequenced individuals represents a unique opportunity to ascertain how the demographic and selection processes associated with post-domestication history have shaped the diversity of this species. Data generated for the project will also be extremely useful to identify deleterious mutations and polymorphisms with causal effects on complex traits, and thus will contribute to new knowledge that could be used in genomic prediction and genome-wide association studies

Observation of two new Ξb−\Xi_b^- baryon resonances

Two structures are observed close to the kinematic threshold in the $\Xi_b^0 \pi^-$ mass spectrum in a sample of proton-proton collision data, corresponding to an integrated luminosity of 3.0 fb$^{-1}$ recorded by the LHCb experiment. In the quark model, two baryonic resonances with quark content $bds$ are expected in this mass region: the spin-parity $J^P = \frac{1}{2}^+$ and $J^P=\frac{3}{2}^+$ states, denoted $\Xi_b^{\prime -}$ and $\Xi_b^{*-}$. Interpreting the structures as these resonances, we measure the mass differences and the width of the heavier state to be $m(\Xi_b^{\prime -}) - m(\Xi_b^0) - m(\pi^{-}) = 3.653 \pm 0.018 \pm 0.006$ MeV$/c^2$, $m(\Xi_b^{*-}) - m(\Xi_b^0) - m(\pi^{-}) = 23.96 \pm 0.12 \pm 0.06$ MeV$/c^2$, $\Gamma(\Xi_b^{*-}) = 1.65 \pm 0.31 \pm 0.10$ MeV, where the first and second uncertainties are statistical and systematic, respectively. The width of the lighter state is consistent with zero, and we place an upper limit of $\Gamma(\Xi_b^{\prime -}) < 0.08$ MeV at 95% confidence level. Relative production rates of these states are also reported.Comment: 17 pages, 2 figure

Measurement of the Bs0→J/ψηB_{s}^{0} \rightarrow J/\psi \eta lifetime

Using a data set corresponding to an integrated luminosity of $3 fb^{-1}$, collected by the LHCb experiment in $pp$ collisions at centre-of-mass energies of 7 and 8 TeV, the effective lifetime in the $B^0_s \rightarrow J/\psi \eta$ decay mode, $\tau_{\textrm{eff}}$, is measured to be $\tau_{\textrm{eff}} = 1.479 \pm 0.034~\textrm{(stat)} \pm 0.011 ~\textrm{(syst)}$ ps. Assuming $CP$ conservation, $\tau_{\textrm{eff}}$ corresponds to the lifetime of the light $B_s^0$ mass eigenstate. This is the first measurement of the effective lifetime in this decay mode.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2016-017.htm

Measurement of the mass and lifetime of the Ωb−\Omega_b^- baryon

A proton-proton collision data sample, corresponding to an integrated luminosity of 3 fb$^{-1}$ collected by LHCb at $\sqrt{s}=7$ and 8 TeV, is used to reconstruct $63\pm9$ $\Omega_b^-\to\Omega_c^0\pi^-$, $\Omega_c^0\to pK^-K^-\pi^+$ decays. Using the $\Xi_b^-\to\Xi_c^0\pi^-$, $\Xi_c^0\to pK^-K^-\pi^+$ decay mode for calibration, the lifetime ratio and absolute lifetime of the $\Omega_b^-$ baryon are measured to be \begin{align*} \frac{\tau_{\Omega_b^-}}{\tau_{\Xi_b^-}} &= 1.11\pm0.16\pm0.03, \\ \tau_{\Omega_b^-} &= 1.78\pm0.26\pm0.05\pm0.06~{\rm ps}, \end{align*} where the uncertainties are statistical, systematic and from the calibration mode (for $\tau_{\Omega_b^-}$ only). A measurement is also made of the mass difference, $m_{\Omega_b^-}-m_{\Xi_b^-}$, and the corresponding $\Omega_b^-$ mass, which yields \begin{align*} m_{\Omega_b^-}-m_{\Xi_b^-} &= 247.4\pm3.2\pm0.5~{\rm MeV}/c^2, \\ m_{\Omega_b^-} &= 6045.1\pm3.2\pm 0.5\pm0.6~{\rm MeV}/c^2. \end{align*} These results are consistent with previous measurements.Comment: 11 pages, 5 figures, All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2016-008.htm