Search

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses-calculation of fixation index (FST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals-were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4-25.1 Mbp and on BTA 16 at 42.5-43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74-5.36 Mbp, BTA 15 at 17.80-18.77 Mbp, and BTA 17 at 45.59-45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40-81.69 Mbp, BTA 15 at 16.04-16.62 Mbp, and BTA 18 at 0.19-1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds., Competing Interests: The authors have declared that no competing interests exist.

Reindeer (Rangifer tarandus L. 1758) are an essential element of the Russian Far North, providing a significant source of nutrition for the representatives of 18 ethnicities. The species has wild and domestic forms, which are in constant interaction. The aim of our study was to characterize the genetic structure of domestic and wild reindeer populations, using a genome-wide bovine genotyping array (BovineHD BeadChip). The wild reindeer samples were obtained from the western Taymyr Peninsula population and from the taiga and tundra populations in the Sakha Republic (Yakutia). The domestic populations included the Evenk, Even, and Chukotka-Khargin breeds of Yakutia and the Nenets breed from the Nenets Autonomous district and Murmansk region. The level of genetic diversity was higher for the wild population. Analyzing Neighbor-Net tree, multidimensional scaling, and Structure results, we observed strong genetic population structure and clear differentiation between domestic and wild populations. All regional populations of domestic reindeer were clearly separated, while wild reindeer showed similar genetic backgrounds. Nevertheless, we found contrasting patterns in the genetic structure of the tundra and taiga reindeer, in accordance with their morphological and ecological differences. Thus, our study revealed a clear genetic differentiation between domestic and wild reindeer populations. It provides novel insights into the genetic diversity and structure of reindeer populations, to support resource utilization and aid in the development of genetic improvement strategies and conservation programs for this species., Competing Interests: The authors have declared that no competing interests exist.

Feed conversion efficiency is a measure of how well an animal converts feed into live weight and it is typically expressed as feed conversion ratio (FCR). FCR and related traits like growth rate (e.g. days to 110 kg-D110) are of high interest for animal breeders, farmers and society due to implications on animal performance, feeding costs and environmental sustainability. The objective of this study was to identify genomic regions associated with FCR and D110 in pigs. A total of 952 terminal line boars, showing an individual variation in FCR, were genotyped using 60K SNP-Chips. Markers were tested for associations with estimated breeding values (EBV) for FCR and D110. For FCR, the largest number of associated SNPs was located on chromosomes 4 (30 SNPs), 1 (25 SNPs), X (15 SNPs) and 6 (12 SNPs). The most prominent genomic regions for D110 were identified on chromosomes 15 (10 SNPs), 1 and 4 (both 9 SNPs). The most significantly associated SNPs for FCR and D110 mapped 129.8 Kb from METTL11B (chromosome 4) and 32Kb from MBD5 (chromosome 15), respectively. A list of positional genes, closest to significantly associated SNPs, was used to identify enriched pathways and biological functions related to the QTL for both traits. A number of candidate genes were significantly overrepresented in pathways of immune cell trafficking, lymphoid tissue structure, organ morphology, endocrine system function, lipid metabolism, and energy production. After resequencing the coding region of selected positional and functional candidate genes, six SNPs were genotyped in a subset of boars. SNPs in PRKDC, SELL, NR2E1 and AKRIC3 showed significant associations with EBVs for FCR/D110. The study revealed a number of chromosomal regions and candidate genes affecting FCR/D110 and pointed to corresponding biological pathways related to lipid metabolism, olfactory reception, and also immunological status.

The glucocorticoid receptor (GR) is a central player in the neuroendocrine stress response; it mediates feedback regulation of the hypothalamus-pituitary-adrenal (HPA) axis and physiological actions of glucocorticoids in the periphery. Despite intensive investigations of GR in the context of receptor-ligand interaction, only recently the first naturally occurring gain-of-function substitution, Ala610Val, of the ligand binding domain was identified in mammals. We showed that this mutation underlies a major quantitative trait locus for HPA axis activity in pigs, reducing cortisol production by about 40-50 percent. To unravel the molecular mechanisms behind this gain of function, receptor-ligand interactions were evaluated in silico, in vitro and in vivo. In accordance with previously observed phenotypic effects, the mutant Val610 GR showed significantly increased activation in response to glucocorticoid and non-glucocorticoid steroids, and, as revealed by GR-binding studies in vitro and in pituitary glands, enhanced ligand binding. Concordantly, the protein structure prediction depicted reduced binding distances between the receptor and ligand, and altered interactions in the ligand binding pocket. Consequently, the Ala610Val substitution opens up new structural information for the design of potent GR ligands and to examine effects of the enhanced GR responsiveness to glucocorticoids on the entire organism., Competing Interests: The authors have declared that no competing interests exist.

Variations in egg incubation temperatures can have acute or long-term effects on gene transcription in avian species. Altered gene expression may, in turn, affect muscle traits in poultry and indirectly influence commercial production. To determine how changes in eggshell temperature affect gene expression, incubation temperatures were varied [36.8°C (low), 37.8°C (control), 38.8°C (high)] at specific time periods reflecting two stages of myogenesis [embryonic days (ED) 7-10 and 10-13]. Gene expression was compared between interventions and matching controls by microarrays in broiler breast muscle at ED10 or ED13 and post-hatch at day 35. Early (ED7-10) high incubation temperature (H10ΔC) resulted in 1370 differentially expressed genes (DEGs) in embryos. Ingenuity pathway analysis revealed temporary activation of cell maintenance, organismal development, and survival ability genes, but these effects were not maintained in adults. Late high incubation temperature (ED10-13) (H13ΔC) had slightly negative impacts on development of cellular components in embryos, but a cumulative effect was observed in adults, in which tissue development and nutrition metabolism were affected. Early low incubation temperature (L10ΔC) produced 368 DEGs, most of which were down-regulated and involved in differentiation and formation of muscle cells. In adults, this treatment down-regulated pathways of transcriptional processes, but up-regulated cell proliferation. Late low temperature incubation (L13ΔC) produced 795 DEGs in embryos, and activated organismal survival and post-transcriptional regulation pathways. In adults this treatment activated cellular and organ development, nutrition and small molecule activity, and survival rate, but deactivated size of body and muscle cells. Thermal interventions during incubation initiate immediate and delayed transcriptional responses that are specific for timing and direction of treatment. Interestingly, the transcriptional response to transiently decreased incubation temperature, which did not affect the phenotypes, prompts compensatory effects reflecting resilience. In contrast, higher incubation temperature triggers gene expression and has long-term effects on the phenotype. These mechanisms of considerable phenotypic plasticity contribute to the biodiversity and broaden the basis for managing poultry populations., Competing Interests: The authors have declared that no competing interests exist.

Haematological traits are important traits that show associations with immune and metabolic status, as well as diseases in humans and animals. Mapping genome regions that affect the blood cell traits can contribute to the identification of genomic features useable as biomarkers for immune, disease and metabolic status. A genome-wide association study (GWAS) was conducted using PorcineSNP60 BeadChips. Single-marker and Bayesian multi-marker approaches were integrated to identify genomic regions and corresponding genes overlapping for both methods. GWAS was performed for haematological traits of 591 German Landrace pig. Heritability estimates for haematological traits were medium to high. In total 252 single SNPs associated with 12 haematological traits were identified (NegLog10 of p-value > 5). The Bayesian multi-marker approach revealed 102 QTL regions across the genome, indicated by 1-Mb windows with contribution to additive genetic variance above 0.5%. The integration of both methods resulted in 24 overlapping QTL regions. This study identified overlapping QTL regions from single- and multi-marker approaches for haematological traits. Identifying candidate genes that affect blood cell traits provides the first step towards the understanding of the molecular basis of haematological phenotypes.

The pig shows genetical and physiological resemblance to human, which predestines it as an experimental animal model especially for mucosal physiology. Therefore, the intestinal epithelial cell lines 1 and J2 (IPEC-1, IPEC-J2)--spontaneously immortalised cell lines from the porcine intestine--are important tools for studying intestinal function. A microarray (GeneChip Porcine Genome Array) was performed to compare the genome wide gene expression of IPECs. Different significantly up-regulated pathways were identified, like "lysosome", "pathways in cancer", "regulation of actin cytoskeleton" and "oxidative phosphorylation" in IPEC-J2 in comparison to IPEC-1. On the other hand, "spliceosome", "ribosome", "RNA-degradation" and "tight junction" are significantly down-regulated pathways in IPEC-J2 in comparison to IPEC-1. Examined pathways were followed up by functional analyses. ATP-, oxygen, glucose and lactate-measurement provide evidence for up-regulation of oxidative phosphorylation in IPEC-J2. These cells seem to be more active in their metabolism than IPEC-1 cells due to a significant higher ATP-content as well as a higher O2- and glucose-consumption. The down-regulated pathway "ribosome" was followed up by measurement of RNA- and protein content. In summary, IPEC-J2 is a morphologically and functionally more differentiated cell line in comparison to IPEC-1. In addition, IPEC-J2 cells are a preferential tool for in vitro studies with the focus on metabolism.

In our previous study, we identified an miRNA regulatory network involved in energy metabolism in porcine muscle. To better understand the involvement of miRNAs in cellular ATP production and energy metabolism, here we used C2C12 myoblasts, in which ATP levels increase during differentiation, to identify miRNAs modulating these processes. ATP level, miRNA and mRNA microarray expression profiles during C2C12 differentiation into myotubes were assessed. The results suggest 14 miRNAs (miR-423-3p, miR-17, miR-130b, miR-301a/b, miR-345, miR-15a, miR-16a, miR-128, miR-615, miR-1968, miR-1a/b, and miR-194) as cellular ATP regulators targeting genes involved in mitochondrial energy metabolism (Cox4i2, Cox6a2, Ndufb7, Ndufs4, Ndufs5, and Ndufv1) during C2C12 differentiation. Among these, miR-423-3p showed a high inverse correlation with increasing ATP levels. Besides having implications in promoting cell growth and cell cycle progression, its function in cellular ATP regulation is yet unknown. Therefore, miR-423-3p was selected and validated for the function together with its potential target, Cox6a2. Overexpression of miR-423-3p in C2C12 myogenic differentiation lead to decreased cellular ATP level and decreased expression of Cox6a2 compared to the negative control. These results suggest miR-423-3p as a novel regulator of ATP/energy metabolism by targeting Cox6a2.

Understanding the genetic contributions behind skeletal muscle composition and metabolism is of great interest in medicine and agriculture. Attempts to dissect these complex traits combine genome-wide genotyping, expression data analyses and network analyses. Weighted gene co-expression network analysis (WGCNA) groups genes into modules based on patterns of co-expression, which can be linked to phenotypes by correlation analysis of trait values and the module eigengenes, i.e. the first principal component of a given module. Network hub genes and regulators of the genes in the modules are likely to play an important role in the emergence of respective traits. In order to detect common regulators of genes in modules showing association with meat quality traits, we identified eQTL for each of these genes, including the highly connected hub genes. Additionally, the module eigengene values were used for association analyses in order to derive a joint eQTL for the respective module. Thereby major sites of orchestrated regulation of genes within trait-associated modules were detected as hotspots of eQTL of many genes of a module and of its eigengene. These sites harbor likely common regulators of genes in the modules. We exemplarily showed the consistent impact of candidate common regulators on the expression of members of respective modules by RNAi knockdown experiments. In fact, Cxcr7 was identified and validated as a regulator of genes in a module, which is involved in the function of defense response in muscle cells. Zfp36l2 was confirmed as a regulator of genes of a module related to cell death or apoptosis pathways. The integration of eQTL in module networks enabled to interpret the differentially-regulated genes from a systems perspective. By integrating genome-wide genomic and transcriptomic data, employing co-expression and eQTL analyses, the study revealed likely regulators that are involved in the fine-tuning and synchronization of genes with trait-associated expression.

Defining phenotypes according to molecular features would promote the knowledge of functional traits like behaviour in both human and animal research. Beside physiological states or environmental factors, an innate predisposition of individual coping strategies was discussed, including the proactive and reactive pattern. According to backtest reactivity, animals assigned as high-resisting (proactive) and low-resisting (reactive) were immune challenged with tetanus toxoid in a time course experiment. Using the Affymetrix platform and qPCR, individual coping characteristics were reflected as gene expression signatures in porcine peripheral blood mononuclear cells (PBMC) at naïve state (day 0) and in response to the model antigen (day 14, day 28, and day 140). Further, the blood cell count was analysed at all stages. On the transcriptional level, processes acting on cell communication, vasculogenesis, and blood coagulation were highlighted in high-resisting animals at naïve state (day 0), temporarily blurred due to immune challenge (day 14) but subsequently restored and intensified (day 28). Notably, similar amounts of white and red blood cells, platelets and haematocrit between high-resisting and low-resisting samples suggest coping-specific expression patterns rather than alterations in blood cell distribution. Taken together, the gene expression patterns indicate that proactive pigs might favour molecular pathways enabling an effective strategy for defence and recovery. This corroborates the previously suggested belief, that proactive animals are prone to an increased number of injuries as an evolutionary inherited mechanism. In contrast to previous assumptions, coping-specific immunity in pigs lacks inherited shifts between cellular and humoral immune responses.

The aim of the present study was to characterize in vivo genome-wide transcriptional responses to immune stimulation in order to get insight into the resulting changes of allocation of resources. Vaccination with tetanus toxoid was used as a model for a mixed Th1 and Th2 immune response in pig. Expression profiles of PBMCs (peripheral blood mononuclear cells) before and at 12 time points over a period of four weeks after initial and booster vaccination at day 14 were studied by use of Affymetrix GeneChip microarrays and Ingenuity Pathway Analysis (IPA). The transcriptome data in total comprised more than 5000 genes with different transcript abundances (DE-genes). Within the single time stages the numbers of DE-genes were between several hundred and more than 1000. Ingenuity Pathway Analysis mainly revealed canonical pathways of cellular immune response and cytokine signaling as well as a broad range of processes in cellular and organismal growth, proliferation and development, cell signaling, biosynthesis and metabolism. Significant changes in the expression profiles of PBMCs already occurred very early after immune stimulation. At two hours after the first vaccination 679 DE-genes corresponding to 110 canonical pathways of cytokine signaling, cellular immune response and other multiple cellular functions were found. Immune competence and global disease resistance are heritable but difficult to measure and to address by breeding. Besides QTL mapping of immune traits gene expression profiling facilitates the detection of functional gene networks and thus functional candidate genes.

The improvement of meat quality and production traits has high priority in the pork industry. Many of these traits show a low to moderate heritability and are difficult and expensive to measure. Their improvement by targeted breeding programs is challenging and requires knowledge of the genetic and molecular background. For this study we genotyped 192 artificial insemination boars of a commercial line derived from the Swiss Large White breed using the PorcineSNP60 BeadChip with 62,163 evenly spaced SNPs across the pig genome. We obtained 26 estimated breeding values (EBVs) for various traits including exterior, meat quality, reproduction, and production. The subsequent genome-wide association analysis allowed us to identify four QTL with suggestive significance for three of these traits (p-values ranging from 4.99×10⁻⁶ to 2.73×10⁻⁵). Single QTL for the EBVs pH one hour post mortem (pH1) and carcass length were on pig chromosome (SSC) 14 and SSC 2, respectively. Two QTL for the EBV rear view hind legs were on SSC 10 and SSC 16.

Glucocorticoids produced in the adrenal cortex under the control of the hypothalamic-pituitary axis play a vital role in the maintenance of basal and stress-related homeostasis and influence health and well-being. To identify loci affecting regulation of the hypothalamic-pituitary-adrenal (HPA) axis in the pig we performed a genome-wide association study for two parameters of acute and long-term adrenal activity: plasma cortisol level and adrenal weight. We detected a major quantitative trait locus at the position of the glucocorticoid receptor gene (NR3C1) - a key regulator of HPA axis activity. To determine the causal variant(s), we resequenced the coding region of NR3C1 and found three missense single nucleotide polymorphisms (SNPs). SNP c.1829C>T, leading to a p.Ala610Val substitution in the ligand binding domain, showed large (about 0.6× and 1.2× phenotypic standard deviations for cortisol level and adrenal weight, respectively), and highly significant (2.1E-39≤log10(1/p)≤1.7E+0) negative effects on both traits. We were able to replicate the association in three commercial pig populations with different breed origins. We analyzed effects of the p.Ala610Val substitution on glucocorticoid-induced transcriptional activity of porcine glucocorticoid receptor (GR) in vitro and determined that the substitution introduced by SNP c.1829C>T increased sensitivity of GR by about two-fold. Finally, we found that non-coding polymorphisms in linkage disequilibrium with SNP c.1829C>T have only a minor effect on the expression of NR3C1 in tissues related to the HPA axis. Our findings provide compelling evidence that SNP c.1829C>T in porcine NR3C1 is a gain-of-function mutation with a major effect on the activity of the adrenal gland. Pigs carrying this SNP could provide a new animal model to study neurobiological and physiological consequences of genetically based GR hypersensitivity and adrenal hypofunction.

The intestinal epithelial cell layer represents the border between the luminal and systemic side of the gut. The decision between absorption and exclusion of substances is the quintessential function of the gut and varies along the gut axis. Consequently, potentially toxic substances may reach the basolateral domain of the epithelial cell layer via blood stream. The mycotoxin deoxynivalenol (DON) is a Fusarium derived secondary metabolite known to enter the blood stream and displaying a striking toxicity on the basolateral side of polarised epithelial cell layers in vitro. Here we analysed potential mechanisms of apical and basolateral DON toxicity reflected in the gene expression. We used the jejunum-derived, polarised intestinal porcine epithelial cell line IPEC-J2 as an in vitro cell culture model. Luminal and systemic DON challenge of the epithelial cell layer was mimicked by a DON application from the apical or basolateral compartment of membrane inserts for 72 h. We compared the genome-wide gene expression of untreated and DON-treated IPEC-J2 cells with the GeneChip® Porcine Genome Array of Affymetrix. Low basolateral DON (200 ng/mL) application triggered 10 times more gene transcripts in comparison to the corresponding apical application (2539 versus 267) despite the intactness of the challenged cell layer as measured by transepithelial electrical resistance. Analysis of the regulated genes by bioinformatic resource DAVID identified several groups of biochemical pathways modulated by concentration and orientation of DON application. Selected genes representing pathways of the cellular metabolism, information processing and structural design were analysed in detail by quantitative PCR. Our findings clearly show that apical and basolateral challenge of epithelial cell layers trigger different gene response profiles paralleled with a higher susceptibility towards basolateral challenge. The evaluation of toxicological potentials of mycotoxins should take this difference in gene regulation dependent on route of application into account.

Background: In various animal models pregnancy diets have been shown to affect offspring phenotype. Indeed, the underlying programming of development is associated with modulations in birth weight, body composition, and continual diet-dependent modifications of offspring metabolism until adulthood, producing the hypothesis that the offspring's transcriptome is permanently altered depending on maternal diet., Methodology/principal Findings: To assess alterations of the offspring's transcriptome due to gestational protein supply, German Landrace sows were fed isoenergetic diets containing protein levels of either 30% (high protein--HP) or 12% (adequate protein--AP) throughout their pregnancy. Offspring muscle tissue (M. longissimus dorsi) was collected at 94 days post conception (dpc), and 1, 28, and 188 days post natum (dpn) for use with Affymetrix GeneChip Porcine Genome Arrays and subsequent statistical and Ingenuity pathway analyses. Numerous transcripts were found to have altered abundance at 94 dpc and 1 dpn; at 28 dpn no transcripts were altered, and at 188 dpn only a few transcripts showed a different abundance between diet groups. However, when assessing transcriptional changes across developmental time points, marked differences were obvious among the dietary groups. Depending on the gestational dietary exposure, short- and long-term effects were observed for mRNA expression of genes related to cell cycle regulation, energy metabolism, growth factor signaling pathways, and nucleic acid metabolism. In particular, the abundance of transcripts related to cell cycle remained divergent among the groups during development., Conclusion: Expression analysis indicates that maternal protein supply induced programming of the offspring's genome; early postnatal compensation of the slight growth retardation obvious at birth in HP piglets resulted, as did a permanently different developmental alteration and responsiveness to the common environment of the transcriptome. The transcriptome modulations are interpreted as the molecular equivalent of developmental plasticity of the offspring that necessitates adaptation and maintenance of the organismal phenotype.

Embryonic implantation to establish a pregnancy is a complex process that requires appropriate communication between the embryo and the maternal endometrium. Inadequate uterine receptivity may contribute to the majority of implantation failures. To provide a comprehensive inventory of genes and functional networks that represent the maternal input of the embryo-maternal cross-talk, a longitudinal, holistic study of the endometrial transcriptome in relation to the days of estrous and to the receptivity of the endometrium was performed in bovine. At day 3 of estrous, genes related to cell communication and mitochondrial energy metabolism were differentially expressed among high- and low-receptive endometria (HR, LR); at day 7, transcripts functioning in immune and inflammatory pathways, oxidative stress, and angiogenesis had different abundances. Additionally, temporal transcriptional changes between days 3 and 7 differed considerably among HR and LR. Further, several transcription factors were predicted as relevant for receptivity because they were either differentially expressed among HR and LR animals or are known to be associated with genes we detected to have differential expression. Finally, global DNA methylation varied according to the interaction of receptivity group and day of estrous, and a divergent trend, which correlated with abundance of DNMT1 transcript, was observed in LR and HR along the estrous cycle days. The study revealed that, even in early estrous, transcripts related to cell communication and response to exogenous stimuli, vascularization, and energy supply show divergent expression and longitudinal temporal regulation in HR and LR. Key components of these molecular pathways are known to be dependent on ovarian hormones that promote uterine receptivity.

In rodent models and in humans the impact of gestational diets on the offspring's phenotype was shown experimentally and epidemiologically. The underlying programming of fetal development was shown to be associated with an increased risk of degenerative diseases in adulthood, including the metabolic syndrome. There are clues that diet-dependent modifications of the metabolism during fetal life can persist until adulthood. This leads to the hypothesis that the offspring's transcriptomes show short-term and long-term changes depending on the maternal diet. To this end pregnant German landrace gilts were fed either a high protein diet (HP, 30% CP) or an adequate protein diet (AP, 12% CP) throughout pregnancy. Hepatic transcriptome profiles of the offspring were analyzed at prenatal (94 dpc) and postnatal stages (1, 28, 188 dpn). Depending on the gestational dietary exposure, mRNA expression levels of genes related to energy metabolism, N-metabolism, growth factor signaling pathways, lipid metabolism, nucleic acid metabolism and stress/immune response were affected either in a short-term or in a long-term manner. Gene expression profiles at fetal stage 94 dpc were almost unchanged between the diets. The gestational HP diet affected the hepatic expression profiles at prenatal and postnatal stages. The effects encompassed a modulation of the genome in terms of an altered responsiveness of energy and nutrient sensing pathways. Differential expression of genes related to energy production and nutrient utilization contribute to the maintenance of development and growth performance within physiological norms, however the modulation of these pathways may be accompanied by a predisposition for metabolic disturbances up to adult stages.