Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Bos indicus introgression into (peri‐)alpine cattle breeds – evidence from the analysis of bovine whey protein variants

The major bovine whey proteins, α‐lactalbumin (α‐LA) and β‐lactoglobulin (β‐LG), exhibit breed‐specific genetic variation. The aim of this study was to identify possible new protein variants and determine the distribution of variants across a variety of 18 taurine and indicine cattle breeds applying a DNA‐based sequencing approach. To this end, the open reading frames of the respective genes (LALBA and LGB) were sequenced in 476 animals. Within the LALBA gene, a previously unknown synonymous and a previously undesignated non‐synonymous nucleotide exchange were identified. Furthermore, two known α‐LA variants (A and B) and four known β‐LG variants (A, B, C and W) were determined. The occurrence of typical indicine variants in some taurine cattle breeds, such as Suisse Eringer, German Hinterwälder and Hungarian Grey Steppe, further supports the hypothesis of ancient Bos indicus introgression into (peri‐)alpine cattle breeds.     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

Genome‐wide differential DNA methylation in tropically adapted Creole cattle and their Iberian ancestors

Enhancing climate resilience and sustainable production for animals in harsh environments are important goals for the livestock industry given the predicted impacts of climate change. Rapid adaptation to extreme climatic conditions has already been imposed on livestock species, including those exported after Columbus's arrival in the Americas. We compared the methylomes of two Creole cattle breeds living in tropical environments with their putative Spanish ancestors to understand the epigenetic mechanisms underlying rapid adaptation of a domestic species to a new and more physiologically challenging environment. Reduced representation bisulfite sequencing was used to assess differences in methylation in Creole and Spanish samples and revealed 334 differentially methylated regions using high stringency parameters (P‐value <0.01, ≥4 CpGs within a distance of 200 bp, mean methylation difference >25%) annotated to 263 unique features. Gene ontology analysis revealed candidate genes involved in tropical adaptation processes, including genes differentially hyper‐ or hypomethylated above 80% in Creole samples displaying biological functions related to immune response (IRF6, PTGDR, FAM19A5, PGLYRP1), nervous system (GBX2, NKX2‐8, RPGR), energy management (BTD), heat resistance (CYB561) and skin and coat attributes (LGR6). Our results entail that major environmental changes imposed on Creole cattle have had an impact on their methylomes measurable today, which affects genes implicated in important pathways for adaptation. Although further work is needed, this first characterization of methylation patterns driven by profound environmental change provides a valuable pointer for the identification of biomarkers of resilience for improved cattle performance and welfare under predicted climatic change models.     

γ‐glutamyl transpeptidase 1 specifically suppresses green‐light avoidance via GABAA receptors in Drosophila

Drosophila larvae innately show light avoidance behavior. Compared with robust blue‐light avoidance, larvae exhibit relatively weaker green‐light responses. In our previous screening for genes involved in larval light avoidance, compared with control w¹¹¹⁸ larvae, larvae with γ‐glutamyl transpeptidase 1 (Ggt‐1) knockdown or Ggt‐1 mutation were found to exhibit higher percentage of green‐light avoidance which was mediated by Rhodopsin6 (Rh6) photoreceptors. However, their responses to blue light did not change significantly. By adjusting the expression level of Ggt‐1 in different tissues, we found that Ggt‐1 in malpighian tubules was both necessary and sufficient for green‐light avoidance. Our results showed that glutamate levels were lower in Ggt‐1 null mutants compared with controls. Feeding Ggt‐1 null mutants glutamate can normalize green‐light avoidance, indicating that high glutamate concentrations suppressed larval green‐light avoidance. However, rather than directly, glutamate affected green‐light avoidance indirectly through GABA, the level of which was also lower in Ggt‐1 mutants compared with controls. Mutants in glutamate decarboxylase 1, which encodes GABA synthase, and knockdown lines of the GABA_A receptor, both exhibit elevated levels of green‐light avoidance. Thus, our results elucidate the neurobiological mechanisms mediating green‐light avoidance, which was inhibited in wild‐type larvae.

     

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D‐loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups—T1a, T2, T3, T4 and T5—were identified in Bos taurus, whereas Bos indicus harbored two haplogroups—I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan‐Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

A lycopene β‐cyclase/lycopene ε‐cyclase/light‐harvesting complex‐fusion protein from the green alga Ostreococcus lucimarinus can be modified to produce α‐carotene and β‐carotene at different ratios

Biosynthesis of asymmetric carotenoids such as α‐carotene and lutein in plants and green algae involves the two enzymes lycopene β‐cyclase (LCYB) and lycopene ε‐cyclase (LCYE). The two cyclases are closely related and probably resulted from an ancient gene duplication. While in most plants investigated so far the two cyclases are encoded by separate genes, prasinophyte algae of the order Mamiellales contain a single gene encoding a fusion protein comprised of LCYB, LCYE and a C‐terminal light‐harvesting complex (LHC) domain. Here we show that the lycopene cyclase fusion protein from Ostreococcus lucimarinus catalyzed the simultaneous formation of α‐carotene and β‐carotene when heterologously expressed in Escherichia coli. The stoichiometry of the two products in E. coli could be altered by gradual truncation of the C‐terminus, suggesting that the LHC domain may be involved in modulating the relative activities of the two cyclase domains in the algae. Partial deletions of the linker region between the cyclase domains or replacement of one or both cyclase domains with the corresponding cyclases from the green alga Chlamydomonas reinhardtii resulted in pronounced shifts of the α‐carotene‐to‐β‐carotene ratio, indicating that both the relative activities of the cyclase domains and the overall structure of the fusion protein have a strong impact on the product stoichiometry. The possibility to tune the product ratio of the lycopene cyclase fusion protein from Mamiellales renders it useful for the biotechnological production of the asymmetric carotenoids α‐carotene or lutein in bacteria or fungi.     

PMC,a potent hydrophilic α‐tocopherol derivative,inhibits NF‐κB activation via PP2A but not IκBα‐dependent signals in vascular smooth muscle cells

The hydrophilic α‐tocopherol derivative, 2,2,5,7,8‐pentamethyl‐6‐hydroxychromane (PMC), is a promising alternative to vitamin E in clinical applications. Critical vascular inflammation leads to vascular dysfunction and vascular diseases, including atherosclerosis, hypertension and abdominal aortic aneurysms. In this study, we investigated the mechanisms of the inhibitory effects of PMC in vascular smooth muscle cells (VSMCs) exposed to pro‐inflammatory stimuli, lipopolysaccharide (LPS) combined with interferon (IFN)‐γ. Treatment of LPS/IFN‐γ‐stimulated VSMCs with PMC suppressed the expression of inducible nitric oxide synthase (iNOS) and matrix metalloproteinase‐9 in a concentration‐dependent manner. A reduction in LPS/IFN‐γ‐induced nuclear factor (NF)‐κB activation was also observed in PMC‐treated VSMCs. The translocation and phosphorylation of p65, protein phosphatase 2A (PP2A) inactivation and the formation of reactive oxygen species (ROS) were significantly inhibited by PMC in LPS/IFN‐γ‐activated VSMCs. However, neither IκBα degradation nor IκB kinase (IKK) or ribosomal s6 kinase‐1 phosphorylation was affected by PMC under these conditions. Both treatments with okadaic acid, a PP2A‐selective inhibitor, and transfection with PP2A siRNA markedly reversed the PMC‐mediated inhibition of iNOS expression, NF‐κB‐promoter activity and p65 phosphorylation. Immunoprecipitation analysis of the cellular extracts of LPS/IFN‐γ‐stimulated VSMCs revealed that p65 colocalizes with PP2A. In addition, p65 phosphorylation and PP2A inactivation were induced in VSMCs by treatment with H₂O₂, but neither IκBα degradation nor IKK phosphorylation was observed. These results collectively indicate that the PMC‐mediated inhibition of NF‐κB activity in LPS/IFN‐γ‐stimulated VSMCs occurs through the ROS‐PP2A‐p65 signalling cascade, an IKK‐IκBα‐independent mechanism. Therapeutic interventions using PMC may therefore be beneficial for the treatment of vascular inflammatory diseases.     

Fine‐mapping the POLL locus in Brahman cattle yields the diagnostic marker CSAFG29

The POLL locus has been mapped to the centromeric region of bovine chromosome 1 (BTA1) in both taurine breeds and taurine–indicine crosses in an interval of approximately 1 Mb. It has not yet been mapped in pure‐bred zebu cattle. Despite several efforts, neither causative mutations in candidate genes nor a singular diagnostic DNA marker has been identified. In this study, we genotyped a total of 68 Brahman cattle and 20 Hereford cattle informative for the POLL locus for 33 DNA microsatellites, 16 of which we identified de novo from the bovine genome sequence, mapping the POLL locus to the region of the genes IFNAR2 and SYNJ1. The 303‐bp allele of the new microsatellite, CSAFG29, showed strong association with the POLL allele. We then genotyped 855 Brahman cattle for CSAFG29 and confirmed the association between the 303‐bp allele and POLL. To determine whether the same association was found in taurine breeds, we genotyped 334 animals of the Angus, Hereford and Limousin breeds and 376 animals of the Brangus, Droughtmaster and Santa Gertrudis composite taurine–zebu breeds. The association between the 303‐bp allele and POLL was confirmed in these breeds; however, an additional allele (305 bp) was also associated but not fully predictive of POLL. Across the data, CSAFG29 was in sufficient linkage disequilibrium to the POLL allele in Australian Brahman cattle that it could potentially be used as a diagnostic marker in that breed, but this may not be the case in other breeds. Further, we provide confirmatory evidence that the scur phenotype generally occurs in animals that are heterozygous for the POLL allele.     

Genome‐wide association analysis for quantitative trait loci influencing Warner–Bratzler shear force in five taurine cattle breeds

We performed a genome‐wide association study for Warner–Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single‐nucleotide polymorphisms (SNPs) within μ‐calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within‐ and across‐breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across‐breed analysis were moderately correlated (0.31–0.66) with those from the individual within‐breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across‐breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within‐breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across‐breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome‐wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine‐map the CAPN1 causal mutation to a 4581‐bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.     

LIGHT/IFN‐γ triggers β cells apoptosis via NF‐κB/Bcl2‐dependent mitochondrial pathway

LIGHT recruits and activates naive T cells in the islets at the onset of diabetes. IFN‐γ secreted by activated T lymphocytes is involved in beta cell apoptosis. However, whether LIGHT sensitizes IFNγ‐induced beta cells destruction remains unclear. In this study, we used the murine beta cell line MIN6 and primary islet cells as models for investigating the underlying cellular mechanisms involved in LIGHT/IFNγ – induced pancreatic beta cell destruction. LIGHT and IFN‐γ synergistically reduced MIN6 and primary islet cells viability; decreased cell viability was due to apoptosis, as demonstrated by a significant increase in Annexin V⁺ cell percentage, detected by flow cytometry. In addition to marked increases in cytochrome c release and NF‐κB activation, the combination of LIGHT and IFN‐γ caused an obvious decrease in expression of the anti‐apoptotic proteins Bcl‐2 and Bcl‐xL, but an increase in expression of the pro‐apoptotic proteins Bak and Bax in MIN6 cells. Accordingly, LIGHT deficiency led to a decrease in NF‐κB activation and Bak expression, and peri‐insulitis in non‐obese diabetes mice. Inhibition of NF‐κB activation with the specific NF‐κB inhibitor, PDTC (pyrrolidine dithiocarbamate), reversed Bcl‐xL down‐regulation and Bax up‐regulation, and led to a significant increase in LIGHT‐ and IFN‐γ‐treated cell viability. Moreover, cleaved caspase‐9, ‐3, and PARP (poly (ADP‐ribose) polymerase) were observed after LIGHT and IFN‐γ treatment. Pretreatment with caspase inhibitors remarkably attenuated LIGHT‐ and IFNγ‐induced cell apoptosis. Taken together, our results indicate that LIGHT signalling pathway combined with IFN‐γ induces beta cells apoptosis via an NF‐κB/Bcl2‐dependent mitochondrial pathway.     

A genome‐wide scan for selection signatures in Nellore cattle

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777 000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty‐one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high‐density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.     

The Importance of Toll‐like Receptors in NF‐κB Signaling Pathway Activation by Helicobacter pylori Infection and the Regulators of this Response

Helicobacter pylori (H. pylori) is a common pathogenic bacterium in the stomach that infects almost half of the population worldwide and is closely related to gastric diseases and some extragastric diseases, including iron‐deficiency anemia and idiopathic thrombocytopenic purpura. Both the Maastricht IV/Florence consensus report and the Kyoto global consensus report have proposed the eradication of H. pylori to prevent gastric cancer as H.pylori has been shown to be a major cause of gastric carcinogenesis. The interactions between H. pylori and host receptors induce the release of the proinflammatory cytokines by activating proinflammatory signaling pathways such as nuclear factor kappa B (NF‐κB), which plays a central role in inflammation, immune response, and carcinogenesis. Among these receptors, Toll‐like receptors (TLRs) are classical pattern recognition receptors in the recognition of H. pylori and the mediation of the host inflammatory and immune responses to H. pylori. TLR polymorphisms also contribute to the clinical consequences of H. pylori infection. In this review, we focus on the functions of TLRs in the NF‐κB signaling pathway activated by H. pylori, the regulators modulating this response, and the functions of TLR polymorphisms in H.pylori‐related diseases.