Transcriptional profiling and miRNA-dependent regulatory network analysis of longissimus dorsi muscle during prenatal and adult stages in two distinct pig breeds.

MicroRNAs (miRNAs) and mRNAs establish a complex regulatory network influencing diverse biological pathways including muscle development and growth. Elucidating miRNA-dependent regulatory networks involved in muscle development could provide additional insights into muscle traits largely predefined during prenatal development. The present study aimed to determine differentially expressed transcripts and functional miRNA-mRNA relationships associated with different stages of skeletal muscle development in two pig breeds, German Landrace and Pietrain, distinct in muscle characteristics. A comparative transcriptional profiling of longissimus dorsi muscle tissues from fetuses at 35, 63 and 91 days post-conception as well as adult pigs (180 days postnatum) was performed using the Affymetrix GeneChip porcine genome microarray. Differential expression patterns were identified to be associated with muscularly developmental stages and breed types. The integration of miRNA expression data and ingenuity pathways analysis (ipa) pathway analysis revealed several miRNA-dependent regulatory networks related to muscle growth and development. The present results provide insights into muscle biology for further improvement of porcine meat quality.

Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model.

Traits related to fatness, important as economic factors in pork production, are associated with serious diseases in humans. Genetical genomics is a useful approach for studying the effects of genetic variation at the molecular level in biological systems. Here we applied a whole-genome association analysis to hepatic gene expression traits, focusing on transcripts with expression levels that correlated with fatness traits in a porcine model. A total of 150 crossbred pigs [Pietrain × (German Large White × German Landrace)] were studied for transcript levels in the liver. The 24K Affymetrix expression microarrays and 60K Illumina single nucleotide polymorphism (SNP) chips were used for genotyping. A total of 663 genes, whose expression significantly correlated with the trait "fat area," were analyzed for enrichment of functional annotation groups as defined in the Ingenuity Pathways Knowledge Base (IPKB). Genes involved in metabolism of various macromolecules and nutrients as well as functions related to dynamic cellular processes correlated with fatness traits. Regions affecting the transcription levels of these genes were mapped and revealed 4,727 expression quantitative trait loci (eQTL) at P < 10⁻5, including 448 cis-eQTL. In this study, genome-wide association analysis of trait-correlated expression was successfully used in a porcine model to display molecular networks and list genes relevant to fatness traits.

Association of ZYX polymorphisms with carcass and meat quality traits in commercial pigs.

Zyxin (ZYX) is one of the proteins in focal adhesions along the actin fibers playing a role in actin organization and signal transduction. By radiation hybrid and genetic mapping we assigned ZYX to porcine chromosome 18 in the area of quantitative traits loci for carcass and meat quality and muscle fiber traits and hence considered ZYX a functional positional candidate gene. Analysis of a newly detected SNPs (c.+279 C>T, c.+399 A>G, c.+522 A>G) in pigs from different commercial breeds (Pietrain [Pi], German Landrace [LR], German Large White x German Landrace [F1] and PiF1) revealed a significant association with carcass traits (including: side- and backfat thickness, loin weight and carcass lean content) and meat quality traits (including: pH, color and drip loss). However, the lack of consistent association across all pig populations in this study indicates that the association of the SNPs may be depending on causal mutations in linkage disequilibrium and/or interactions with other loci.

Expression quantitative trait loci analysis of genes in porcine muscle by quantitative real-time RT-PCR compared to microarray data.

Genetic analysis of transcriptional profiling is a promising approach for identifying biological pathways and dissecting the genetics of complex traits. Here, we report on expression quantitative trait loci (eQTL) that were estimated from the quantitative real-time RT-PCR data of 276 F(2) animals and compared with eQTL identified using 74 microarrays. In total, 13 genes were selected that showed trait-dependent expression in microarray experiments and exhibited 21 eQTL. Real-time RT-PCR and microarray data revealed seven cis eQTL in total, of which one was only detected by real-time RT-PCR, one was only detected by microarray analysis, three were consistently found in overlapping intervals and two were in neighbouring intervals on the same chromosome; whereas no trans eQTL was confirmed. We demonstrate that cis regulation is a stable characteristic of individual transcripts. Consequently, a global microarray eQTL analysis of a limited number of samples can be used for exploring functional and regulatory gene networks and scanning for cis eQTL, whereas the subsequent analysis of a subset of likely cis-regulated genes by real-time RT-PCR in a larger number of samples is relevant to narrow down a QTL region by targeting these positional candidate genes. In fact, when modelling SNPs of six genes as fixed effects in the eQTL analysis, eQTL peaks were shifted downwards, experimentally confirming the impact of the respective polymorphic genes, although these SNPs were not located in the regulatory sequence and these shifts occur as a result of linkage disequilibrium in the F(2) population.

The spatial expression pattern of antimicrobial peptides across the healthy bovine udder.

Antimicrobial peptides are key molecules in local host defense. With the aim to better understand the possible involvement of these peptides in the prevention of bovine mastitis, we determined, for the first time to our knowledge, the spatial pattern of constitutive expression of 5 bovine beta-defensins and bovine psoriasin (S100A7) across 5 localizations of the bovine mammary gland applying a quantitative real-time PCR approach. We observed 3 different expression patterns in the healthy udder: 1) constitutive expression of the lingual and tracheal antimicrobial peptides (LAP and TAP), as well that of bovine neutrophil beta-defensins 4 and 10 (BNBD4 and BNBD10), is essentially restricted to the mammary lymph node; 2) bovine beta-defensin 1 (DEFB1) is mainly expressed in the cisternal epithelium and the Rosette of Fürstenberg; 3) strong constitutive mRNA expression of the calcium-binding protein S100A7, which is also known as psoriasin and which has been reported to be highly active against Escherichia coli, was detected in the streak canal. These results indicate a crucial role of S100A7 in the early-stage prevention of coliform mastitis, and the analyzed beta-defensins might be regarded as inducible weapons against already invading pathogens.

The utility of focus group interviews to capture dietary consumption data in the distant past: dairy consumption in Kazakhstan villages 50 years ago.

From 1949 to 1962, residents of several villages in Kazakhstan received substantial doses of radiation to the thyroid gland resulting from nuclear tests conducted at the Semipalatinsk Nuclear Test Site. The primary source of radiation was internal from an intake of radioactive iodine by consumption of contaminated dairy products. A previous research study of childhood exposure and thyroid disease in this region gathered limited data on study participants' dairy intake at the time of the fallout for the purpose of estimating past radiation doses. As many participants were too young at the time of the nuclear tests to recall dietary consumption and existing sources of archival data are limited, it was necessary to interview parents and other village residents who cared for children during this time - older adults ranging in age from 75 to 90 years. Results from 11 focus group interviews conducted in 2007 with 82 women from 4 villages in Kazakhstan yielded group-level estimates of age-, gender-, ethnicity- and village-specific dairy consumption patterns in rural Kazakhstan during the 1950s. Children typically consumed cow's milk with limited consumption of mare, goat and sheep milk; and consumed dairy products such as sour milk (airan), soft cottage cheese (tvorog) and fermented mare milk (koumiss) with the greatest amounts of koumiss reported at ages 15-21 years. The consumption patterns differed by age, and between Kazakh and Russian children, which should lead to different estimates of radiation exposure to the thyroid. This study showed the utility of focus groups to obtain quantitative estimates for dietary intake in the distant past.

Analysis of key molecules of the innate immune system in mammary epithelial cells isolated from marker-assisted and conventionally selected cattle.

Mastitis is the most prevalent infectious disease in dairy herds. Breeding programs considering mastitis susceptibility were adopted as approaches to improve udder health status. In recent decades, conventional selection criteria based on phenotypic characteristics such as somatic cell score in milk have been widely used to select animals. Recently, approaches to incorporate molecular information have become feasible because of the detection of quantitative trait loci (QTL) affecting mastitis resistance. The aims of the study were to explore molecular mechanisms underlying mastitis resistance and the genetic mechanisms underlying a QTL on Bos taurus chromosome 18 found to influence udder health. Primary cell cultures of mammary epithelial cells from heifers that were selected for high or low susceptibility to mastitis were established. Selection based on estimated pedigree breeding value or on the basis of marker-assisted selection using QTL information was implemented. The mRNA expression of 10 key molecules of the innate immune system was measured using quantitative real-time PCR after 1, 6, and 24 h of challenge with heat-inactivated mastitis pathogens (Escherichia coli and Staphylococcus aureus) and expression levels in the high and low susceptibility groups were compared according to selection criteria. In the marker-assisted selection groups, mRNA expression in cells isolated from less-susceptible animals was significantly elevated for toll-like receptor 2, tumor necrosis factor-alpha, IL-1beta, IL-6, IL-8, RANTES (regulated upon activation, normal t-cell expressed and secreted), complement factor C3, and lactoferrin. In the estimated pedigree breeding value groups, mRNA expression was significantly elevated only for V-rel reticuloendotheliosis viral oncogene homolog A, IL-1 beta, and RANTES. These observations provide first insights into genetically determined divergent reactions to pathogens in the bovine mammary gland and indicate that the application of QTL information could be a successful tool for the selection of animals resistant to mastitis.

Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2.

Principal component analysis of traits related to carcass and meat properties were combined with microarray expression data for the identification of functional networks of genes and biological processes taking place during the conversion of muscle to meat. Principal components (PCs) with high loadings of meat quality traits were derived from phenotypic data of 572 animals of a porcine crossbreed population. Microarray data of 74 M. longissimus dorsi samples were correlated with PC datasets. Lists of significantly correlated genes were analyzed for enrichment of functional annotation groups as defined in the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases as well as the Ingenuity Pathways Analysis library. Ubiquitination, phosphorylation, mitochondrion dysfunction, actin, integrin, platelet-derived growth factor, epidermal growth factor, vascular endothelial growth factor, and Ca signaling pathways are correlated with meat quality. Among the significantly trait-associated genes, CAPZB, ANKRD1, and CTBP2 are promoted as candidate genes for meat quality that provide a link between the highlighted pathways. Knowledge of the relevant biological processes and the differential expression of members of the pathway will provide tools that are predictive for traits related to meat quality and that may also be diagnostic for many muscle defects or damages including muscle atrophy, dystrophy, and hypoxia.

Staphylococcus aureus and Escherichia coli cause deviating expression profiles of cytokines and lactoferrin messenger ribonucleic acid in mammary epithelial cells.

Pathogens invading the mammary gland cause a complex signaling network that activates the early immune defense and leads to an outcome of inflammation symptoms. To examine the importance of mammary epithelial cells in these regulations and interactions resulting in a pathogen-related course of mastitis, we characterized the mRNA expression profile of key molecules of the innate immune system by quantitative real-time PCR. Mammary gland epithelial cells isolated on d 42 of lactation from 28 first-lactation Holstein dairy cows were cultured separately under standardized conditions and treated for 1, 6, and 24 h with heat-inactivated gram-positive (Staphylococcus aureus) and gram-negative (Escherichia coli) bacteria. Both pathogens increased mRNA expression patterns of proteins involved in pathogen recognition such as Toll-like receptors and nuclear factor-kappa B, whereas gram-negatives acted as a stronger stimulus. Furthermore, this could be confirmed by the expression profile of the proinflammatory cytokines tumor necrosis factor alpha, IL-1 beta, IL-6, and chemokines such as IL-8 and RANTES (regulated upon activation, normal T-cell expressed and secreted). Remarkably, at a low level of mRNA expression after 1 h of treatment these cytokines and chemokines were expressed at a significantly higher level in Staphyloccocus aureus than in Escherichia coli affected cells. Lactoferrin showed a deviating expression pattern to pathogen stimulation (i.e., at the 1-h measuring point Escherichia coli induced a higher mRNA expression, whereas the highest level was reached after 24 h of stimulation with Staphylococcus aureus). Complement factor 3 was the only measured factor that responded equally to both microorganisms. Our data emphasize the role of mammary epithelial cells in the immune defense of the udder and confirm their contribution to pathogen-related different courses of mastitis.

YAC/BAC contig spanning the MHC class III region of cattle.

A contig of the class III region of the bovine major histocompatibility complex (MHC) was established from bacterial and yeast artificial chromosomes using PCR and BAC-end sequencing. The marker content of individual clones was determined by gene and BAC-end specific PCR, and the location of genes and BAC-ends was confirmed analyzing somatic hybrid cells. A comparative analysis indicated that the content and order of MHC class III genes is strongly conserved between cattle and other mammalian species. Fluorescence in situ hybridization localized the bovine class III region to BTA23q21-->q22. The results show that the collection of sequenced BAC-ends is a powerful resource for generating high-resolution comparative chromosome maps.

Trait-associated expressed hepatic and intestine genes in cattle of different metabolic type--putative functional candidates for nutrient utilization.

The present study aimed at identifying bovine hepatic and intestinal DNA sequences expressed breed specifically as potential functional candidate genes for nutrient transformation. Transcript levels of 29 expressed sequence tags (ESTs) were analysed comparatively in the liver and intestine of growing Charolais and German Holstein bulls by real-time RT-PCR. In previous studies, these ESTs were characterized as differentially displayed in mRNA differential display of cows varying in metabolic type and harbouring single nucleotide polymorphisms. Breed-specific gene expression levels indicate significantly increased hepatic metabolic activity in Charolais and increased intestinal metabolic activity in German Holstein bulls. Transcript levels of six functional genes measured in liver (NDUFB8, NACA, UAP1, SAH) and intestine (FUS/TLS, APOC3), respectively, support this assumption. The observed coincidence of metabolic type-specific expressed ESTs with variant ESTs showing breed-specific allele distribution points to functional genetic variants located in the vicinity of the analysed sequences. In addition, location of most of the breed specifically expressed ESTs within chromosome regions known to be affecting carcass and growth traits in cattle supports the putative candidate gene character of the ESTs identified.

Mapping of quantitative trait loci for lactation persistency traits in German Holstein dairy cattle.

A whole genome scan to map quantitative trait loci (QTL) for persistency of milk yield (PMY), persistency of fat yield (PFY), persistency of protein yield (PPY) and persistency of milk energy yield (PEY) was performed in a granddaughter design in the German Holstein dairy cattle population. The analysis included 16 paternal half-sib families with a total of 872 bulls. The analysis was carried out for the first lactation and for the first three lactations combined using univariate weighted multimarker regression. Controlling the false discovery rate across traits and data sets at a level of 0.15 and treating the four persistency traits as different traits revealed 27 significant QTL. A total of 12 chromosomes showed significant QTL effects on a chromosomewise basis. The DGAT1 effect was highly significant for PPY and protein yield. A haplotype analysis using results of previous studies of the same design revealed a co-segregation of various persistency QTL and QTL affecting health traits like dystocia and stillbirth and functional traits like non-return rate 90 and somatic cell score.

The DGAT1 K232A mutation is not solely responsible for the milk production quantitative trait locus on the bovine chromosome 14.

The gene, acyl-CoA:diacylglycerol acyltransferase1 (DGAT1), was recently identified as the one underlying the quantitative trait locus (QTL) for milk production traits in the centromeric region of the bovine chromosome 14. Until now, 2 alleles, the lysine variant (increasing fat yield, fat and protein percentage) and the alanine variant (increasing protein and milk yield), were postulated at DGAT1. This study investigated whether the diallelic DGAT1 polymorphism is responsible for all the genetic variation at the centromeric region of this chromosome for milk, fat, and protein yield and fat and protein percentage. A statistical model was applied to a granddaughter design to analyze 16 German Holstein families. The model included the diallelic DGAT1 effect and the QTL transition probability estimated for each chromosomal position by a multiple marker approach. Because the regression coefficient of this probability was corrected for the diallelic DGAT1 polymorphism, it represented a putative conditional QTL effect. The effect of the DGAT1 gene was always highly significant. The conditional QTL effect was significant genomewise for fat percentage at the proximal end of the chromosome and for protein percentage at a more distal chromosomal region. Additional chromosomewise significance was found for fat and protein yield. Our results suggest an additional source of genetic variance on this chromosome for these traits; either one or more additional alleles segregating at DGAT1 that were not previously detected, a second quantitative trait locus affecting these traits, or both.

Mapping of QTL for Body Conformation and Behavior in Cattle.

Genome scans for quantitative trait loci (QTL) in farm animals have concentrated on primary production and health traits, and information on QTL for other important traits is rare. We performed a whole genome scan in a granddaughter design to detect QTL affecting body conformation and behavior in dairy cattle. The analysis included 16 paternal half-sib families of the Holstein breed with 872 sons and 264 genetic markers. The markers were distributed across all 29 autosomes and the pseudoautosomal region of the sex chromosomes with average intervals of 13.9 cM and covering an estimated 3155.5 cM. All families were analyzed jointly for 22 traits using multimarker regression and significance thresholds determined empirically by permutation. QTL that exceeded the experiment-wise significance threshold (5% level) were detected on chromosome 6 for foot angle, teat placement, and udder depth, and on chromosome 29 for temperament. QTL approaching experiment-wise significance (10% level) were located on chromosome 6 for general quality of feet and legs and general quality of udder, on chromosome 13 for teat length, on chromosome 23 for general quality of feet and legs, and on chromosome 29 for milking speed. An additional 51 QTL significant at the 5% chromosome-wise level were distributed over 21 chromosomes. This study provides the first evidence for QTL involved in behavior of dairy cattle and identifies QTL for udder conformation on chromosome 6 that could form the basis of recently reported QTL for clinical mastitis.

Clarifications on breakpoints in HSAX and BTAX by comparative mapping of F9, HPRT, and XIST in cattle.

The coagulation factor IX gene (F9), the hypoxanthine phosphoribosyl transferase 1 gene (HPRT1), and the X-inactive specific transcript gene (XIST) were physically assigned in cattle to analyze chromosomal breakpoints on BTAX recently identified by radiation hybrid (RH) mapping experiments. Whereas the FISH assignment of XIST indicates a similar location on the q-arm of the human and cattle X chromosomes, the locus of HPRT1 supported the assumption of a chromosome rearrangement between the distal half of the q-arm of HSAX and the p-arm of BTAX identified by RH mapping. F9 previously located on the q-arm of BTAX was assigned to the p-arm of BTAX using RH mapping and FISH. The suggested new position of F9 close to HPRT1 supports the homology between HSAXq and BTAXp. The F9 locus corresponds with the gene order found in the homologous human chromosome segment. XIST was assigned on BTAXq23, HPRT1 and F9 were mapped to BTAXp22, and the verification of the location of F9 in a 5000 rad cattle-hamster whole genome radiation hybrid panel linked the gene to markers URB10 and HPRT1.

Polymorphism of the bovine CSN1S1 promoter: linkage mapping, intragenic haplotypes, and effects on milk production traits.

The bovine CSN1S1 5' flanking region (CSN1S1-5') was screened for polymorphisms in different cattle breeds. Single-strand conformation polymorphisms (SSCP) and sequence analyses revealed four alleles (1-4), two of them being new allelic forms (3 and 4). Sequences were deposited in GenBank with accession numbers AF549499-502. In alleles 1 and 4, potential transcription factor binding sites are altered by the mutations. Using SSCP analysis, all four alleles were identified in German Holsteins. Six intragenic haplo-types comprising CSN1S1-5' (alleles 1, 2, 3, 4) and exon 17 (CSN1S1*B and C) genotypes were found. Linkage mapping using half-sib families from the German QTL project positioned CSN1S1 between the markers FBN14 and CSN3, with 5.6 cM distance between CSN1S1 and CSN3. Variance analysis, using family and CSN1S1 promoter genotypes as fixed effects, of breeding values and deregressed proofs for milk production traits (milk, fat, and protein yield and also fat and protein percentage) revealed significant effects on protein percentage when all families and genotypes were considered. Contrast calculations assigned a highly significant effect to genotype 24, which was associated with highest LS-means for protein percentage breeding values. As CSN1S1 is one of the main caseins in milk, this could be an effect of mutations in regulatory elements in the promoter region. An effect on milk yield breeding values was indicated for genotype 12, but is probably caused by a linked locus.

Quantitative trait loci mapping of functional traits in the German Holstein cattle population.

A whole-genome scan to detect quantitative trait loci (QTL) for functional traits was performed in the German Holstein cattle population. For this purpose, 263 genetic markers across all autosomes and the pseudoautosomal region of the sex chromosomes were genotyped in 16 granddaughter-design families with 872 sons. The traits investigated were deregressed breedingvalues for maternal and direct effects on dystocia (DYSm, DYSd) and stillbirth (STIm, STId) as well as maternal and paternal effects on nonreturn rates of 90 d (NR90m, NR90p). Furthermore, deregressed breeding values for functional herd life (FHL) and daughter yield deviation for somatic cell count (SCC) were investigated. Weighted multimarker regression analyses across families and permutation tests were applied for the detection of QTL and the calculation of statistical significance. A ten percent genomewise significant QTL was localized for DYSm on chromosome 8 and for SCC on chromosome 18. A further 24 putative QTL exceeding the 5% chromosomewise threshold were detected. On chromosomes 7, 8, 10, 18, and X/Yps, coincidence of QTL for several traits was observed. Our results suggest that loci with influence on udder health may also contribute to genetic variance of longevity. Prior to implementation of these QTL in marker assisted selection programs for functional traits, information about direct and correlated effects of these QTL as well as fine mapping of their chromosomal positions is required.

Mapping of the bovine blood group systems J, N', R', and Z show evidence for oligo-genetic inheritance.

Genes determining the bovine erythrocyte antigens were mapped by linkage analysis. In total 9591 genotypes of 20 grandsire families with 1074 sires from a grand-daughter design were elucidated for the genes determining the erythrocyte antigens EAA, EAB, EAC, EAF, EAJ, EAL, EAM, EAN', EAR', EAS, EAT', and EAZ according to standard paternity testing procedures in the blood typing laboratories. Linkage analyses were performed with 248 microsatellite markers, eight SSCP markers and four polymorphic proteins and enzymes covering the 29 autosomes and the pseudoautosomal region of the sex chromosomes. The number of informative meioses for the blood group systems ranged from 76 to 947. Blood group systems EAM and EAT' were non-informative. Most of the erythrocyte antigen loci showed significant linkage to a single chromosome and were mapped unequivocally. The genes determining erythrocyte antigen EAA, EAB, EAC, EAL, and EAS were mapped to chromosomes 15, 12, 18, 3, and 21, respectively. Lod-score values ranged from 11.43 to 107.83. Moreover, the EAF system could be mapped to chromosome 17. However, the EAN' system previously known as part of the EAF system could be mapped to chromosome 5. In addition, the blood group systems EAJ, the new EAN', EAR', and EAZ, showed significant linkage to microsatellite markers on various chromosomes and also to other blood groups. The appearance of a single blood group system might be therefore either dependent on the existence of other blood group systems or because of an interaction between different loci on various chromosomes as is known in humans and in pigs.