Design and characterization of a high-resolution multiple-SNP capture array by target sequencing for sheep

The efficiency of molecular breeding largely depends on inexpensive genotyping arrays. In this study, we aimed to develop an ovine high-resolution multiple-single-nucleotide polymorphism (SNP) capture array, based on genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology. All the markers were from 40K captured regions, including genes located within selective sweep regions, breed-specific regions, quantitative trait loci (QTL), and the potential functional SNPs on the sheep genome. The results showed that a total of 210K high-quality SNPs were identified in the 40K regions, indicating a high average capture ratio (99.7%) for the target genomic regions. Using genotyped data (n = 317) from liquid chip technology, we further performed genome-wide association studies (GWAS) to detect the genetic loci affecting sheep hair types and teat number. A single significant association signal for hair types was identified on 6.7-7.1 Mb of chromosome 25. The IRF2BP2 gene (chr25: 7,067,974-7,071,785), which is located within this genomic region, has been previously known to be involved in hair/wool traits in sheep. The results further showed a new candidate region around 26.4 Mb of chromosome 13, between the ARHGAP21 and KIAA1217 genes, that was significantly related to teat number in sheep. The haplotype patterns of this region also showed differences in animals with 2, 3, or 4 teats. Advances in using the high-accuracy and low-cost liquid chip are expected to accelerate sheep genomic and breeding studies in the coming years.

Rapid detection of deformed wing virus in honeybee using ultra-rapid qPCR and a DNA-chip

Fast and accurate detection of viral RNA pathogens is important in apiculture. A polymerase chain reaction (PCR)-based detection method has been developed, which is simple, specific, and sensitive. In this study, we rapidly (in 1 min) synthesized cDNA from the RNA of deformed wing virus (DWV)-infected bees (Apis mellifera), and then, within 10 min, amplified the target cDNA by ultra-rapid qPCR. The PCR products were hybridized to a DNA-chip for confirmation of target gene specificity. The results of this study suggest that our method might be a useful tool for detecting DWV, as well as for the diagnosis of RNA virus-mediated diseases on-site.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Genome-wide association and genotype by environment interactions for growth traits in U.S. Gelbvieh cattle

BACKGROUND

Single nucleotide polymorphism (SNP) arrays have facilitated discovery of genetic markers associated with complex traits in domestic cattle; thereby enabling modern breeding and selection programs. Genome-wide association analyses (GWAA) for growth traits were conducted on 10,837 geographically diverse U.S. Gelbvieh cattle using a union set of 856,527 imputed SNPs. Birth weight (BW), weaning weight (WW), and yearling weight (YW) were analyzed using GEMMA and EMMAX (via imputed genotypes). Genotype-by-environment (GxE) interactions were also investigated.

RESULTS

GEMMA and EMMAX produced moderate marker-based heritability estimates that were similar for BW (0.36-0.37, SE = 0.02-0.06), WW (0.27-0.29, SE = 0.01), and YW (0.39-0.41, SE = 0.01-0.02). GWAA using 856K imputed SNPs (GEMMA; EMMAX) revealed common positional candidate genes underlying pleiotropic QTL for Gelbvieh growth traits on BTA6, BTA7, BTA14, and BTA20. The estimated proportion of phenotypic variance explained (PVE) by the lead SNP defining these QTL (EMMAX) was larger and most similar for BW and YW, and smaller for WW. Collectively, GWAAs (GEMMA; EMMAX) produced a highly concordant set of BW, WW, and YW QTL that met a nominal significance level (P ≤ 1e-05), with prioritization of common positional candidate genes; including genes previously associated with stature, feed efficiency, and growth traits (i.e., PLAG1, NCAPG, LCORL, ARRDC3, STC2). Genotype-by-environment QTL were not consistent among traits at the nominal significance threshold (P ≤ 1e-05); although some shared QTL were apparent at less stringent significance thresholds (i.e., P ≤ 2e-05).

CONCLUSIONS

Pleiotropic QTL for growth traits were detected on BTA6, BTA7, BTA14, and BTA20 for U.S. Gelbvieh beef cattle. Seven QTL detected for Gelbvieh growth traits were also recently detected for feed efficiency and growth traits in U.S. Angus, SimAngus, and Hereford cattle. Marker-based heritability estimates and the detection of pleiotropic QTL segregating in multiple breeds support the implementation of multiple-breed genomic selection.

Genomic Divergence in Swedish Warmblood Horses Selected for Equestrian Disciplines

The equestrian sport horse Swedish Warmblood (SWB) originates from versatile cavalry horses. Most modern SWB breeders have specialized their breeding either towards show jumping or dressage disciplines. The aim of this study was to explore the genomic structure of SWB horses to evaluate the presence of genomic subpopulations, and to search for signatures of selection in subgroups of SWB with high or low breeding values (EBVs) for show jumping. We analyzed high density genotype information from 380 SWB horses born in the period 2010-2011, and used Principal Coordinates Analysis and Discriminant Analysis of Principal Components to detect population stratification. Fixation index and Cross Population Extended Haplotype Homozygosity scores were used to scan the genome for potential signatures of selection. In accordance with current breeding practice, this study highlights the development of two separate breed subpopulations with putative signatures of selection in eleven chromosomes. These regions involve genes with known function in, e.g., mentality, endogenous reward system, development of connective tissues and muscles, motor control, body growth and development. This study shows genetic divergence, due to specialization towards different disciplines in SWB horses. This latter evidence can be of interest for SWB and other horse studbooks encountering specialized breeding.

Family-specific genotype arrays increase the accuracy of pedigree-based imputation at very low marker densities

BACKGROUND

In this paper, we evaluate the performance of using family-specific low-density genotype arrays to increase the accuracy of pedigree-based imputation. Genotype imputation is a widely used tool that decreases the costs of genotyping a population by genotyping the majority of individuals on a low-density array and using statistical regularities between the low-density and high-density individuals to fill in the missing genotypes. Previous work on population-based imputation has found that it is possible to increase the accuracy of imputation by maximizing the number of informative markers on an array. In the context of pedigree-based imputation, where the informativeness of a marker depends only on the genotypes of an individual's parents, it may be beneficial to select the markers on each low-density array on a family-by-family basis.

RESULTS

In this paper, we examined four family-specific low-density marker selection strategies and evaluated their performance in the context of a real pig breeding dataset. We found that family-specific or sire-specific arrays could increase imputation accuracy by 0.11 at one marker per chromosome, by 0.027 at 25 markers per chromosome and by 0.007 at 100 markers per chromosome.

CONCLUSIONS

These results suggest that there may be room to use family-specific genotyping for very-low-density arrays particularly if a given sire or sire-dam pairing have a large number of offspring.

MicroRNA profile of immune response in gills of zebrafish (Danio rerio) upon Staphylococcus aureus infection

The gills of fish are large mucosal surfaces that are very important portals for pathogen entry. Investigations have shown that microRNAs (miRNAs) are key regulators of immune response to bacterial infections in the gills of fish; however, how miRNA expression changes in response to infection by Gram-positive bacteria remains largely unknown. To further investigate the immunological role of miRNAs in fish gills under pathogen stress induced by Gram-positive bacterial infection, this study investigated Staphylococcus aureus (SA)-induced changes in the miRNAs levels in gills of adult zebrafish (Danio rerio). miRNA microarrays were used to analyze expression profiles of known miRNA in the gills of zebrafish in response to SA infection and compared these to uninfected control fish. A total of 30 differentially expressed miRNAs (DEMs) were identified. Target genes likely regulated by DEMs were predicted, and functional enrichment analyses were performed. The results indicated that DEM targets were primarily involved in innate immune processes, apoptosis, defense responses, and antibacterial responses. Pathways involving bacterial infection, innate immunity, metabolic process, disease, and apoptosis were mediated by DEMs. Furthermore, real-time quantitative PCR experiments for nine key SA-responsive DEMs that regulated the "SA infection" pathway validated the accuracy of microarray results. Dynamic variations in gene expression were surveyed in detail for these key SA-responsive DEMs for PBS control and at 6, 12, 24, and 48 h after SA challenge in detail. This study provides novel insight into the mechanisms underlying the miRNA regulation during the SA-induced immune response in zebrafish gills, and provides basic knowledge on the innate immune response against Gram-positive bacterial infection in bony fish.

Mammary Gland Transcriptome and Proteome Modifications by Nutrient Restriction in Early Lactation Holstein Cows Challenged with Intra-Mammary Lipopolysaccharide

: The objective is to study the effects of nutrient restrictions, which induce a metabolic imbalance on the inflammatory response of the mammary gland in early lactation cows. The aim is to decipher the molecular mechanisms involved, by comparing a control, with a restriction group, a transcriptome and proteome, after an intra-mammary lipopolysaccharide challenge. Multi-parous cows were either allowed ad libitum intake of a lactation diet (n = 8), or a ration containing low nutrient density (n = 8; 48% barley straw and dry matter basis) for four days starting at 24 ± 3 days in milk. Three days after the initiation of their treatments, one healthy rear mammary quarter of 12 lactating cows was challenged with 50 µg of lipopolysaccharide (LPS). Transcriptomic and proteomic analyses were performed on mammary biopsies obtained 24 h after the LPS challenge, using bovine 44K microarrays, and nano-LC-MS/MS, respectively. Restriction-induced deficits in energy, led to a marked negative energy balance (41 versus 97 ± 15% of Net Energy for Lactation (NEL) requirements) and metabolic imbalance. A microarray analyses identified 25 differentially expressed genes in response to restriction, suggesting that restriction had modified mammary metabolism, specifically β-oxidation process. Proteomic analyses identified 53 differentially expressed proteins, which suggests that the modification of protein synthesis from mRNA splicing to folding. Under-nutrition influenced mammary gland expression of the genes involved in metabolism, thereby increasing β-oxidation and altering protein synthesis, which may affect the response to inflammation.

New host shift from human to cows within Staphylococcus aureus involved in bovine mastitis and nasal carriage of animal's caretakers

Staphylococcus aureus is a commensal and pathogen of both humans and bovines. While the epidemiology of both groups has been extensively studied individually, little is known about the potential zoonotic transfer from animal strains to human being and vice versa. To determine the S. aureus prevalence of bovine mastitis in Algeria and the zoonotic transfer of strains to human beings, mastitis milk samples were collected, and professionals in a close contact with bovines were nasal swabbed. S. aureus isolates were all characterized by methicillin resistance and spa-typing. DNA microarrays analysis was performed on a subset of strains in order to detect other virulence factors, including toxins, and to assign the isolates to theirs MLST clonal complexes. Overall, 116/222 (52.3%) cows suffered from mastitis, whose 38.8% (45/116) infected with S. aureus. Human nasal carriage was of 38% (49/129), with only 4 MRSA carriers (3.1%). A higher diversity of spa-types was observed in human (35/50) than in bovine (18/67) isolates, with a predominance of clonal complexes CC97 and CC22 in bovines. The typical animal clone CC97 was occasionally detected in human beings. Conversely, the CC22 S. aureus clone largely switched from humans to bovines. Our study highlights the potential dynamics of animal and human S. aureus strains in the farm environment in Algeria, which may represent a health threat in both populations.

Accuracy of genotype imputation in Labrador Retrievers

The dog is a valuable model species for the genetic analysis of complex traits, and the use of genotype imputation in dogs will be an important tool for future studies. It is of particular interest to analyse the effect of factors like single nucleotide polymorphism (SNP) density of genotyping arrays and relatedness between dogs on imputation accuracy due to the acknowledged genetic and pedigree structure of dog breeds. In this study, we simulated different genotyping strategies based on data from 1179 Labrador Retriever dogs. The study involved 5826 SNPs on chromosome 1 representing the high density (HighD) array; the low-density (LowD) array was simulated by masking different proportions of SNPs on the HighD array. The correlations between true and imputed genotypes for a realistic masking level of 87.5% ranged from 0.92 to 0.97, depending on the scenario used. A correlation of 0.92 was found for a likely scenario (10% of dogs genotyped using HighD, 87.5% of HighD SNPs masked in the LowD array), which indicates that genotype imputation in Labrador Retrievers can be a valuable tool to reduce experimental costs while increasing sample size. Furthermore, we show that genotype imputation can be performed successfully even without pedigree information and with low relatedness between dogs in the reference and validation sets. Based on these results, the impact of genotype imputation was evaluated in a genome-wide association analysis and genomic prediction in Labrador Retrievers.

Frequency and distribution of 152 genetic disease variants in over 100,000 mixed breed and purebred dogs

Knowledge on the genetic epidemiology of disorders in the dog population has implications for both veterinary medicine and sustainable breeding. Limited data on frequencies of genetic disease variants across breeds exists, and the disease heritage of mixed breed dogs remains poorly explored to date. Advances in genetic screening technologies now enable comprehensive investigations of the canine disease heritage, and generate health-related big data that can be turned into action. We pursued population screening of genetic variants implicated in Mendelian disorders in the largest canine study sample examined to date by examining over 83,000 mixed breed and 18,000 purebred dogs representing 330 breeds for 152 known variants using a custom-designed beadchip microarray. We further announce the creation of MyBreedData (www.mybreeddata.com), an online updated inherited disorder prevalence resource with its foundation in the generated data. We identified the most prevalent, and rare, disease susceptibility variants across the general dog population while providing the first extensive snapshot of the mixed breed disease heritage. Approximately two in five dogs carried at least one copy of a tested disease variant. Most disease variants are shared by both mixed breeds and purebreds, while breed- or line-specificity of others is strongly suggested. Mixed breed dogs were more likely to carry a common recessive disease, whereas purebreds were more likely to be genetically affected with one, providing DNA-based evidence for hybrid vigor. We discovered genetic presence of 22 disease variants in at least one additional breed in which they were previously undescribed. Some mutations likely manifest similarly independently of breed background; however, we emphasize the need for follow up investigations in each case and provide a suggested validation protocol for broader consideration. In conclusion, our study provides unique insight into genetic epidemiology of canine disease risk variants, and their relevance for veterinary medicine, breeding programs and animal welfare.

Like any human, dogs may suffer from or pass on a variety of inherited disorders. Knowledge of how likely a typical dog is to carry an inherited disorder in its genome, and which disorders are the most common and relevant ones across dog breeds, is valuable for both veterinary care and breeding of healthy dogs. We have explored the largest global dog study sample collected to date, consisting of more than 100,000 mixed breed and purebred dogs, to advance research on this subject. We found that mixed breed dogs and purebred dogs potentially suffer from many of the same inherited disorders, and that around two in five dogs carried at least one of the conditions that we screened for. A dog carrying an inherited disorder is not a “bad dog”–but we humans responsible for breeding selections do need to make sustainable decisions avoiding inbreeding, i.e. mating of dogs that are close relatives. The disease prevalence information we generated during this study is made available online (www.mybreeddata.com), as a free tool for breed and kennel clubs, breeders, as well as the veterinary and scientific community.

Genomewide association mapping and pathway analysis of meat tenderness in Polled Nellore cattle

Brazil is one of the world's largest beef exporters, although the product has a low price due to quality issues. The meat exported by Brazil is considered medium and low quality by international buyers, mainly due to lack of tenderness. The predominant Zebu breeds (80% Nellore) are known for producing tougher beef than taurine breeds. Nonetheless, some studies have shown that there is substantial genetic variability for tenderness within the Nellore breed, although it is a difficult trait to improve by conventional selection methods. Therefore, the aim of this study was to perform a genomewide association study (GWAS) and a gene set enrichment analysis to identify genomic regions and biologically relevant pathways associated with meat tenderness in Polled Nellore cattle. Data consisted of Warner-Bratzler shear force values of LM from 427 Polled Nellore animals divided into 3 experimental slaughters (years 2005, 2008, and 2010). The animals were genotyped with either the Illumina BovineHD BeadChip (777k, on 61 samples) or the GGP Indicus HD chip (77k, on 366 samples). Single nucleotide polymorphisms were excluded when the call rate was <90%, the Hardy-Weinberg proportions -value was <1% (Fisher exact test, Bonferroni adjusted), and the minor allele frequency was <1%. Imputation from the GGP Indicus HD chip to the Illumina BovineHD BeadChip was performed using the FImput program. Genomewide association analysis was performed using the Efficient Mixed Model Association eXpedited (EMMAx) and the population parameters previously determined (P3D) methods. The GWAS was complemented with a gene set enrichment analysis performed using the FatiGO procedure. Significant markers ( < 0.0001) explaining a larger proportion of variation than other significant SNPs were located on chromosomes 3, 13, 17, 20, 21, and 25, indicating QTL associated with meat tenderness throughout the genome. Additionally, gene set analysis identified 22 Gene Ontology functional terms and 2 InterPro entries that showed significant enrichment of genes associated with tenderness. The functional categories included protein tyrosine and serine/threonine kinase activity, calcium ion binding, lipid metabolic process, and growth factors, among others. These results help to elucidate the genetic architecture and metabolic pathways underlying this trait, which is of extreme economic and social importance to Brazil, because Nellore is the dominant beef cattle breed in the country.

Multiplex RT-PCR and Automated Microarray for Detection of Eight Bovine Viruses

Microarrays can be a useful tool for pathogen detection as it allow for simultaneous interrogation of the presence of a large number of genetic sequences in a sample. However, conventional microarrays require extensive manual handling and multiple pieces of equipment for printing probes, hybridization, washing and signal detection. In this study, a reverse transcription (RT)-PCR with an accompanying novel automated microarray for simultaneous detection of eight viruses that affect cattle [vesicular stomatitis virus (VSV), bovine viral diarrhoea virus type 1 and type 2, bovine herpesvirus 1, bluetongue virus, malignant catarrhal fever virus, rinderpest virus (RPV) and parapox viruses] is described. The assay accurately identified a panel of 37 strains of the target viruses and identified a mixed infection. No non-specific reactions were observed with a panel of 23 non-target viruses associated with livestock. Vesicular stomatitis virus was detected as early as 2 days post-inoculation in oral swabs from experimentally infected animals. The limit of detection of the microarray assay was as low as 1 TCID50 /ml for RPV. The novel microarray platform automates the entire post-PCR steps of the assay and integrates electrophoretic-driven capture probe printing in a single user-friendly instrument that allows array layout and assay configuration to be user-customized on-site.

Rapid Communication: MicroRNA co-expression network reveals apoptosis in the reproductive tract during molting in laying hens

The aim of this study was to determine the regulatory mechanisms of molting and recrudescence via studying the micro-RNA (miRNA) expression in the oviduct of laying hens. We performed a cDNA microarray analysis in the magnum tissue from the oviduct to identify the whole miRNA profiles through the molting and recrudescence periods. A total of 35 laying hens (47-wk-old) were divided into 7 groups (0 d: a control group; 6 and 12 d: 2 molting-period groups fed on a high-zinc diet; and 20, 25, 30, and 35 d: 4 recrudescence-period groups fed on a normal diet after a 12-d period on a high-zinc diet). An miRNA co-expression network (miRCN) was generated using the differentially expressed miRNA (DEM) according to the entire data integration. The significantly co-expressed miRNA ( = 111) were highly differentially expressed from 12 to 20 d, which was a transition period between molting and recrudescence, while their expression patterns were contrary to the estrogen changes. The targets of highly connected miRNA ( = 12) indicated the significant biological pathways and gene ontology (GO) terms, such as MAPK and Wnt signaling and magnesium-ion binding, which are associated with apoptotic activities. These results suggest that the miRNA of the miRCN might play a role in the apoptotic progression of the reproductive tract during molting.

Genome wide host gene expression analysis in mice experimentally infected with Pasteurella multocida

Pasteurella multocida causes acute septicemic and respiratory diseases, including haemorrhagic septicaemia, in cattle and buffalo with case fatality of 100%. In the present study, mice were infected with P. multocida (1.6 × 10³ cfu, intraperitoneal) to evaluate host gene expression profile at early and late stages of infection using high throughput microarray transcriptome analyses. Several differentially expressed genes (DEGs) at both the time points were identified in P.multocida infected spleen, liver and lungs. Functional annotation of these DEGs showed enrichment of key pathways such as TLR, NF-κB, MAPK, TNF, JAK-STAT and NOD like receptor signaling pathways. Several DEGs overlapped across different KEGG pathways indicating a crosstalk between them. The predicted protein—protein interaction among these DEGs suggested, that the recognition of P. multocida LPS or outer membrane components by TLR4 and CD14, results in intracellular signaling via MyD88, IRAKs and/or TRAF6 leading to activation of NFκB and MAPK pathways and associated cytokines.

Sexual dimorphisms of mRNA and miRNA in human/murine heart disease

BACKGROUND

Sexual dimorphisms are well recognized in various cardiac diseases such as ischemic cardiomyopathy (ICM), hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Thorough understanding of the underlying genetic programs is crucial to optimize treatment strategies specified for each gender. By performing meta-analysis and microarray analysis, we sought to comprehensively characterize the sexual dimorphisms in the healthy and diseased heart at the level of both mRNA and miRNA transcriptome.

RESULTS

Existing mRNA microarray data of both mouse and human heart were integrated, identifying dozens/ hundreds of sexually dimorphic genes in healthy heart, ICM, HCM, and DCM. These sexually dimorphic genes overrepresented gene ontologies (GOs) important for cardiac homeostasis. Further, microarray of miRNA, isolated from mouse sham left ventricle (LV) (n = 6 & n = 5 for male & female) and chronic MI LV (n = 19 & n = 19) and from human normal LV (n = 6 & n = 6) and ICM LV (n = 4 & n = 5), was conducted. This revealed that 13 mouse miRNAs are sexually dimorphic in MI and 6 in normal heart. In human, 3 miRNAs were sexually dimorphic in ICM and 15 in normal heart. These data revealed miRNA-mRNA networks that operate in a sexually-biased fashion.

CONCLUSIONS

mRNA and miRNA transcriptome of normal and disease heart show significant sex differences, which might impact the cardiac homeostasis. Together this study provides the first comprehensive picture of the genome-wide program underlying the heart sexual dimorphisms, laying the foundation for gender specific treatment strategies.

Linkage disequilibrium, persistence of phase, and effective population size in Spanish local beef cattle breeds assessed through a high-density single nucleotide polymorphism chip

Linkage disequilibrium (LD) and persistence of phase are fundamental approaches for exploring the genetic basis of economically important traits in cattle, including the identification of QTL for genomic selection and the estimation of effective population size () to determine the size of the training populations. In this study, we have used the Illumina BovineHD chip in 168 trios of 7 Spanish beef cattle breeds to obtain an overview of the magnitude of LD and the persistence of LD phase through the physical distance between markers. Also, we estimated the time of divergence based on the persistence of the LD phase and calculated past from LD estimates using different alternatives to define the recombination rate. Estimates of average (as a measure of LD) for adjacent markers were close to 0.52 in the 7 breeds and decreased with the distance between markers, although in long distances, some LD still remained (0.07 and 0.05 for markers 200 kb and 1 Mb apart, respectively). A panel with a lower boundary of 38,000 SNP would be necessary to launch a successful within-breed genomic selection program. Persistence of phase, measured as the pairwise correlations between estimates of in 2 breeds at short distances (10 kb), was in the 0.89 to 0.94 range and decreased from 0.33 to 0.52 to a range of 0.01 to 0.08 when marker distance increased from 200 kb to 1 Mb, respectively. The magnitude of the persistence of phase between the Spanish beef breeds was similar to those found in dairy breeds. For across-breed genomic selection, the size of the SNP panels must be in the range of 50,000 to 83,000 SNP. Estimates of past showed values ranging from 26 to 31 for 1 generation ago in all breeds. The divergence among breeds occurred between 129 and 207 generations ago. The results of this study are relevant for the future implementation of within- and across-breed genomic selection programs in the Spanish beef cattle populations. Our results suggest that a reduced subset of the SNP panel would be enough to achieve an adequate precision of the genomic predictions.

Accuracy of genomic breeding values for meat tenderness in Polled Nellore cattle

Zebu () cattle, mostly of the Nellore breed, comprise more than 80% of the beef cattle in Brazil, given their tolerance of the tropical climate and high resistance to ectoparasites. Despite their advantages for production in tropical environments, zebu cattle tend to produce tougher meat than Bos taurus breeds. Traditional genetic selection to improve meat tenderness is constrained by the difficulty and cost of phenotypic evaluation for meat quality. Therefore, genomic selection may be the best strategy to improve meat quality traits. This study was performed to compare the accuracies of different Bayesian regression models in predicting molecular breeding values for meat tenderness in Polled Nellore cattle. The data set was composed of Warner-Bratzler shear force (WBSF) of longissimus muscle from 205, 141, and 81 animals slaughtered in 2005, 2010, and 2012, respectively, which were selected and mated so as to create extreme segregation for WBSF. The animals were genotyped with either the Illumina BovineHD (HD; 777,000 from 90 samples) chip or the GeneSeek Genomic Profiler (GGP Indicus HD; 77,000 from 337 samples). The quality controls of SNP were Hard-Weinberg Proportion -value ≥ 0.1%, minor allele frequency > 1%, and call rate > 90%. The FImpute program was used for imputation from the GGP Indicus HD chip to the HD chip. The effect of each SNP was estimated using ridge regression, least absolute shrinkage and selection operator (LASSO), Bayes A, Bayes B, and Bayes Cπ methods. Different numbers of SNP were used, with 1, 2, 3, 4, 5, 7, 10, 20, 40, 60, 80, or 100% of the markers preselected based on their significance test (-value from genomewide association studies [GWAS]) or randomly sampled. The prediction accuracy was assessed by the correlation between genomic breeding value and the observed WBSF phenotype, using a leave-one-out cross-validation methodology. The prediction accuracies using all markers were all very similar for all models, ranging from 0.22 (Bayes Cπ) to 0.25 (Bayes B). When preselecting SNP based on GWAS results, the highest correlation (0.27) between WBSF and the genomic breeding value was achieved using the Bayesian LASSO model with 15,030 (3%) markers. Although this study used relatively few animals, the design of the segregating population ensured wide genetic variability for meat tenderness, which was important to achieve acceptable accuracy of genomic prediction. Although all models showed similar levels of prediction accuracy, some small advantages were observed with the Bayes B approach when higher numbers of markers were preselected based on their -values resulting from a GWAS analysis.

Differential gene expression profiling of endometrium during the mid-luteal phase of the estrous cycle between a repeat breeder (RB) and non-RB cows

BACKGROUND

Repeat breeding directly affects reproductive efficiency in cattle due to an increase in services per conception and calving interval. This study aimed to investigate whether changes in endometrial gene expression profile are involved in repeat breeding in cows. Differential gene expression profiles of the endometrium were investigated during the mid-luteal phase of the estrous cycle between repeat breeder (RB) and non-RB cows using microarray analysis.

METHODS

The caruncular (CAR) and intercaruncular (ICAR) endometrium of both ipsilateral and contralateral uterine horns to the corpus luteum were collected from RB (inseminated at least three times but not pregnant) and non-RB cows on Day 15 of the estrous cycle (4 cows/group). Global gene expression profiles of these endometrial samples were analyzed with a 15 K custom-made oligo-microarray for cattle. Immunohistochemistry was performed to investigate the cellular localization of proteins of three identified transcripts in the endometrium.

RESULTS

Microarray analysis revealed that 405 and 397 genes were differentially expressed in the CAR and ICAR of the ipsilateral uterine horn of RB, respectively when compared with non-RB cows. In the contralateral uterine horn, 443 and 257 differentially expressed genes were identified in the CAR and ICAR of RB, respectively when compared with non-RB cows. Gene ontology analysis revealed that genes involved in development and morphogenesis were mainly up-regulated in the CAR of RB cows. In the ICAR of both the ipsilateral and contralateral uterine horns, genes related to the metabolic process were predominantly enriched in the RB cows when compared with non-RB cows. In the analysis of the whole uterus (combining the data above four endometrial compartments), RB cows showed up-regulation of 37 genes including PRSS2, GSTA3 and PIPOX and down-regulation of 39 genes including CHGA, KRT35 and THBS4 when compared with non-RB cows. Immunohistochemistry revealed that CHGA, GSTA3 and PRSS2 proteins were localized in luminal and glandular epithelial cells and stroma of the endometrium.

CONCLUSION

The present study showed that endometrial gene expression profiles are different between RB and non-RB cows. The identified candidate endometrial genes and functions in each endometrial compartment may contribute to bovine reproductive performance.

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds

Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.

Temperature-dependent regulation of gene expression in Japanese flounder Paralichthys olivaceus kidney after Edwardsiella tarda formalin-killed cells

Temperature affects the activities of the immune system and the susceptibility of fish to pathogens. To investigate the modulation of temperature on immune related gene expression in formalin-killed cells (FKC) of Edwardsiella tarda-injected Japanese flounder Paralichthys olivaceus, fish reared at 15 or 22 °C were injected with FKC of E. tarda. The up-regulation of immune related genes was detected in FKC-injected fish at both temperatures by qPCR. The mRNA expression of IFNγ was highly up-regulated at 6 h post injection (hpi) in FKC-injected fish at 15 °C, whereas at 22 °C, strong up-regulation of the gene was detected at 3 hpi The mRNA expression level of IRF1 was detected from 3 hpi to day 14 post injection in fish reared at 15 °C, but the gene was up-regulated from 3 to 6 hpi in fish reared at 22 °C. Comprehensive gene expression profiling showed that immune related genes are differentially expressed between 15 and 22 °C. Genes involved in the IFNγ signaling pathway were up-regulated at 22 °C but not at 15 °C. These results demonstrate that gene(s) involved in IFNγ signaling pathway in Japanese flounder stimulated with FKC of E. tarda are regulated by temperature.

Gene transcription profiling in wild and laboratory-exposed eels: Effect of captivity and in situ chronic exposure to pollution

Aquatic ecosystems are subjected to a variety of man-induced stressors but also vary spatially and temporally due to variation in natural factors. In such complex environments, it remains difficult to detect, dissociate and evaluate the effects of contaminants in wild organisms. In this context, the aim of this study was to test whether the hepatic transcriptome profile of fish may be used to detect in situ exposure to a particular contaminant. Transcriptomic profiles from laboratory-exposed and wild eels sampled along a contamination gradient were compared. During laboratory experiments, fish were exposed during 45days to different pollutants (Hg, PCBs, OCPs or Cd) or natural factors (temperature, salinity or low food supply) at levels close to those found in the sampling sites. A strong difference was observed between the transcriptomic profiles obtained from wild and laboratory-exposed animals (whatever the sites or experimental conditions), suggesting a general stress induced by captivity in the laboratory. Among the biological functions that were up-regulated in laboratory eels in comparison to wild eels, histone modification was the most represented. This finding suggests that laboratory conditions could affect the epigenome of fish and thus modulate the transcriptional responses developed by fish in response to pollutant exposure. Among experimental conditions, only the transcription profiles of laboratory animals exposed to cold temperature were correlated with those obtained from wild fish, and more significantly with fish from contaminated sites. Common regulated genes were mainly involved in cell differentiation and liver development, suggesting that stem/progenitor liver cells could be involved in the adaptive response developed by fish chronically exposed to pollutant mixtures.

Ontogenetic onset of immune-relevant genes in the common sole (Solea solea)

Fish are free-living organisms since initial stages of development and are exposed to numerous pathogens before their lymphoid organs have matured and adaptive immunity has developed. Susceptibility to diseases and juvenile mortality represent key critical factors for aquaculture. In this context, the characterization of the appearance kinetics of the immune system key members will be useful in understanding the ability of a particular species in generating immune protection against invading pathogens at different developmental stages. The present study characterized, for the first time, the transcriptional onset of un-explored relevant genes of both innate and adaptive immune system during the Solea solea ontogenesis. Gene expression profiles of immune relevant genes was investigated, by means of DNA microarray, in ten developmental stages, from hatching (1 day post-hatching, dph) to accomplishment of the juvenile form (33 dph). The obtained results revealed that transcripts encoding relevant members of innate immune repertoire, such as lysozyme, AMPs (hepcidin, β-defensin), PPRs and complement components are generally characterized by high expression levels at first stages (i.e. hatch and first feeding) indicating protection from environmental pathogens even at early development. Transcription of adaptive immune genes (i.e. Class I and class II MHC, TCRs) differs from that of the innate immune system. Their onset coincides with metamorphosis and larvae-to-juvenile transition, and likely overlaps with the appearance and maturation of the main lymphoid organs. Finally, data collected suggest that at the end of metamorphosis S. solea cell-mediated immune system hasn't still undergone full maturation.

Differentially expressed genes after viral haemorrhagic septicaemia virus infection in olive flounder (Paralichthys olivaceus)

A strain of viral haemorrhagic septicaemia virus (VHSV) was isolated from cultured olive flounder (Paralichthys olivaceus) during epizootics in South Korean. This strain showed high mortality to olive flounder in in vivo challenge experiment. The complete genomic RNA sequences were determined and phylogenetic analysis of the amino acid sequences of glycoprotein revealed that this isolate was grouped into genotype IVa of genus Novirhabdovirus. Expression profile of genes in olive flounder was analyzed at day 1 and day3 after infection with this VHSV isolate by using cDNA microarray containing olive flounder 13K cDNA clones. Microarray analysis revealed 785 up-regulated genes and 641 down-regulated genes by at least two-fold in virus-infected fish compared to healthy control groups. Among 785 up-regulated genes, we identified seven immune response-associated genes, including the interferon (IFN)-induced 56-kDa protein (IFI56), suppressor of cytokine signaling 1 (SOCS1), interleukin 8 (IL-8), cluster of differentiation 83 (CD83), α-globin (HBA), VHSV-induced protein-6 (VHSV6), and cluster of differentiation antigen 9 (CD9). Our results confirm previous reports that even virulent strain of VHSV induces expression of genes involved in protective immunity against VHSV.

Gene expression profiling in melanised sites of Atlantic salmon fillets

Black spots, which deteriorate quality of Atlantic salmon fillets represent a significant problem for commercial aquaculture. These areas are characterized with accumulation of melanomacrophages, occasional formation of granulomas and substitution of skeletal muscle with connective tissue. A number of possible causative agents have been suggested including vaccination and infection with piscine reovirus (PRV). We report transcriptome profiling of melanised foci with oligonucleotide DNA microarrays. Analyses revealed a multitude of differentially expressed genes associated with melanogenesis, metabolic changes and formation of scar. The immune profile was characterized with inflammation, preferential activation of classical complement pathway, MHCII and helper T cells combined with strong B cells responses and massive induction of immunoglobulins; innate antiviral responses were relatively weak in sharp contrast to PRV-caused heart and skeletal muscle inflammation and other viral infections. A panel of immune genes with specific activation in dark spots was found, most up-regulated were CD209-like lectin (44-fold) and prostaglandin reductase (11-fold). Further, RNA sequencing was performed on the same material to search for the presence of putative pathogens. Transcripts of prokaryotic rRNA with exclusive or preferential location in black spots were found. Results suggest mild chronic inflammation initiated with trauma, bacterial or viral infection followed by sustained immune responses to opportunistic microorganisms as a realistic scenario of dark spots formation.

Gene Expression and Antimicrobial Activity of Bovine Macrophages in Response to Mycobacterium avium subsp. paratuberculosis

We evaluated gene expression and antimicrobial responses of bovine monocyte—derived macrophages incubated with Mycobacterium avium subsp. paratuberculosis (M. a. ptb), the causative agent of Johne's disease. Gene expression was evaluated by the use of human noncompetitive high-density oligonucleotide microarrays. Bovine messenger RNA hybridized with 14.2–18.2% of the 12,600 oligonucleotide probe sets. When macrophages incubated with M. a. ptb were compared with nonactivated control macrophages, macrophages activated by addition of interferon-γ and lipopolysaccharide, and macrophages incubated with Mycobacterium avium subspecies avium (M. a. a), 47, 79, and 27 genes, respectively, were differentially expressed. Differential expression of six of these genes was confirmed using reverse transcriptase polymerase chain reaction. Several functional assays were performed to evaluate the potential relevance of differentially expressed genes to host defense. Macrophages phagocytizing M. a. a had a greater capacity to kill the organisms and to acidify phagosomes and a greater degree of apoptosis than did macrophages incubated with M. a. ptb. The results of these studies indicate that multiple genes and metabolic pathways are differentially expressed by macrophages ingesting mycobacterial organisms. Although the intracellular fate of mycobacterial organisms appears to be dependent on a complex interaction between macrophage and organism, phagosome acidification and apoptosis may play central roles in organism survival.

Characterization of 26 Isolates of Staphylococcus Aureus, Predominantly from Dairy Sheep, Using Four Different Techniques of Molecular Epidemiology

Little information is available regarding the molecular epidemiology of Staphylococcus aureus–induced mastitis in dairy sheep. In this study, 4 different typing techniques were compared in typing 26 S. aureus isolates, predominantly from cases of subclinical mastitis in dairy ewes. The 4 techniques were pulsed-field gel electrophoresis (PFGE), restriction fragment length polymorphism (RFLP) on 2 genes (coagulase and clumping factor B), randomly amplified polymorphic DNA–polymerase chain reaction (PCR) (RAPD-PCR), and multilocus sequence typing (MLST). On the basis of discriminatory power as the key parameter of typing systems, MLST and PFGE were found to be the most powerful techniques. The MLST and PFGE could contribute to epidemiological surveillance and evaluation of mastitis control programs, by documenting prevalence and dissemination of endemic clones in infected populations. The results of this study show that a single clone of S. aureus is widely distributed in infected ewe mammary glands.

Microarray analysis of gene expression in olive flounder liver infected with viral haemorrhagic septicaemia virus (VHSV)

The most fatal viral pathogen in olive flounder Paralichthys olivaceus, is viral hemorrhagic septicemia virus, which afflicts over 48 species of freshwater and marine fish. Here, we performed gene expression profiling on transcripts isolated from VHSV-infected olive flounder livers using a 13 K cDNA microarray chip. A total of 1832 and 1647 genes were upregulated and down-regulated over two-fold, respectively, after infection. A variety of immune-related genes showing significant changes in gene expression were identified in upregulated genes through gene ontology annotation. These genes were grouped into categories such as antibacterial peptide, antigen-recognition and adhesion molecules, apoptosis, cytokine-related pathway, immune system, stress response, and transcription factor and regulatory factors. To verify the cDNA microarray data, we performed quantitative real-time PCR, and the results were similar to the microarray data. In conclusion, these results may be useful for the identification of specific genes or for the diagnosis of VHSV infection in flounder.

Success of interventions in mastitis problems with Staphylococcus aureus after the introduction of an automatic milking system

Staphylococcus aureus (S. aureus) is often the cause of mastitis problems in dairy herds and causes great economic losses. In this study, isolates from a dairy herd with a known S. aureus mastitis problem were examined by means of molecular methods (spa typing, PFGE, and DNA microarray) to investigate their epidemiological relationship and the success of intervention measures. The investigated dairy farm has a herd size of 60 cows and uses a fully automated milking system for milk production. A S. aureus strain, which contaminated the automated milking system and was subsequently spread among the herd through the latter, was suspected to be the origin of the mastitis problem within the herd. Thanks to the applied molecular methods, the common origin of the S. aureus isolates from the collected milk and swab samples could be shown. By culling chronically infected cows, optimising dry cow management and ensuring reliable intermediate cluster disinfection, the bulk milk somatic cell count improved.

Differential Gene Expression in Liver, Gill, and Olfactory Rosettes of Coho Salmon (Oncorhynchus kisutch) After Acclimation to Salinity

Most Pacific salmonids undergo smoltification and transition from freshwater to saltwater, making various adjustments in metabolism, catabolism, osmotic, and ion regulation. The molecular mechanisms underlying this transition are largely unknown. In the present study, we acclimated coho salmon (Oncorhynchus kisutch) to four different salinities and assessed gene expression through microarray analysis of gills, liver, and olfactory rosettes. Gills are involved in osmotic regulation, liver plays a role in energetics, and olfactory rosettes are involved in behavior. Between all salinity treatments, liver had the highest number of differentially expressed genes at 1616, gills had 1074, and olfactory rosettes had 924, using a 1.5-fold cutoff and a false discovery rate of 0.5. Higher responsiveness of liver to metabolic changes after salinity acclimation to provide energy for other osmoregulatory tissues such as the gills may explain the differences in number of differentially expressed genes. Differentially expressed genes were tissue- and salinity-dependent. There were no known genes differentially expressed that were common to all salinity treatments and all tissues. Gene ontology term analysis revealed biological processes, molecular functions, and cellular components that were significantly affected by salinity, a majority of which were tissue-dependent. For liver, oxygen binding and transport terms were highlighted. For gills, muscle, and cytoskeleton-related terms predominated and for olfactory rosettes, immune response-related genes were accentuated. Interaction networks were examined in combination with GO terms and determined similarities between tissues for potential osmosensors, signal transduction cascades, and transcription factors.

Pathogen surveillance in wild bottlenose dolphins Tursiops truncatus

The number and prevalence of diseases is rapidly increasing in the marine ecosystem. Although there is an increase in the number of marine diseases observed world-wide, current understanding of the pathogens associated with marine mammals is limited. An important need exists to develop and apply platforms for rapid detection and characterization of pathogenic agents to assess, prevent and respond to disease outbreaks. In this study, a broad-spectrum molecular detection technology capable of detecting all sequenced microbial organisms, the Lawrence Livermore Microbial Detection Array, was used to assess the microbial agents that could be associated with wild Atlantic dolphins. Blowhole, gastric, and fecal samples from 8 bottlenose dolphins were collected in Charleston, SC, as part of the dolphin assessment effort. The array detected various microbial agents from the dolphin samples. Clostridium perfringens was most prevalent in the samples surveyed using the microarray. This pathogen was also detected using microbiological culture techniques. Additionally, Campylobacter sp., Staphylococcus sp., Erwinia amylovora, Helicobacter pylori, and Frankia sp. were also detected in more than one dolphin using the microarray, but not in culture. This study provides the first survey of pathogens associated with 3 tissue types in dolphins using a broad-spectrum microbial detection microarray and expands insight on the microbial community profile in dolphins.

Transcriptome analyses of Atlantic salmon (Salmo salar L.) erythrocytes infected with piscine orthoreovirus (PRV)

Heart and skeletal muscle inflammation (HSMI) is a widespread disease of farmed Atlantic salmon (Salmo salar L.) and is associated with piscine orthoreovirus (PRV) infection. PRV is detectable in blood long before development of pathology in cardiac- and skeletal muscle appear, and erythrocytes have been identified as important target cells for the virus. The effects of PRV infection on cellular processes of erythrocytes are not known, but haemolytic anemia or systemic lysis of erythrocytes does not seem to occur, even with high virus loads in erythrocytes. In this study, gene expression profiling performed with high-density oligonucleotide microarray showed that PRV infection of erythrocytes induced a large panel of virus responsive genes. These involved interferon-regulated antiviral genes, as well as genes involved in antigen presentation via MHC class I. PRV infection also stimulated negative immune regulators. In contrast, a large number of immune genes expressed prior to infection were down-regulated. Moderate reduction of expression was also found for many genes encoding components of cytoskeleton and myofiber, proteins involved in metabolism, ion exchange, cell-cell interactions as well as growth factors and regulators of differentiation. PRV did not affect expression of genes involved in heme biosynthesis, gas exchange or erythrocyte-specific markers, but some regulators of erythropoiesis showed decreased transcription levels. These results indicate that PRV infection activates innate antiviral immunity in salmon erythrocytes, but suppresses other gene expression programs. Gene expression profiles suggest major phenotypic changes in PRV infected erythrocytes, but the functional consequences remain to be explored.

Temperature-dependent regulation of gene expression in poly (I:C)-treated Japanese flounder, Paralichthys olivaceus

Gene expression profiling of poly (I:C)-treated Japanese flounder, Paralichthys olivaceus, under different temperatures was investigated using microarray analysis. The response was analyzed in spleen tissue at 3 and 24 h post injection (hpi) at 15 °C and 25 °C. A large number of genes in fish treated with poly (I:C) at 25 °C were expressed at 3 hpi, whereas the expression profiles at 24 hpi appeared to be similar to those of the controls. Cluster analysis of the different expression profiles showed three distinct groups of up-regulated genes in fish reared at 15 °C. These were early (3 hpi), early-to-late (3 and 24 hpi), and late (24 hpi) up-regulated genes. These genes included type I IFN-related genes and inflammatory genes. Among the up-regulated genes, most of the type I IFN-related genes played early-to-late- and late-responding genes at 15 °C but early-responding genes at 25 °C. Thus, several up-regulated genes in these groups from the microarray result were further verified by qPCR. These results indicate that the type I IFN gene expressions of P. olivaceus treated with poly (I:C) can be regulated in a temperature-dependent manner.

Evaluation of the impact of camelina oil-containing diets on the expression of genes involved in the innate anti-viral immune response in Atlantic cod (Gadus morhua)

To improve sustainability of aquaculture, especially for carnivorous species like Atlantic cod, replacement of fish oil-based diets with vegetable oil-based diets has been studied. The use of vegetable oil in fish feeds can significantly change the fatty acid composition of fish tissues, and given the importance of fatty acids in inflammation and immunity, this change could potentially impact the immune response and health of the fish. The oilseed Camelina sativa is a promising source for this vegetable oil, because of the high oil content of its seeds (40%), a higher n-3 fatty acid content than most other oilseeds, and a high amount of γ-tocopherol. This study aims to investigate the effect of the replacement of dietary fish oil with oil from Camelina sativa on the immune response of Atlantic cod, as measured by the gene expression in spleen. Juvenile cod were fed on a fish oil-based diet (FO) or one of two diets in which camelina oil replaced 40% or 80% of fish oil (40CO and 80CO respectively) for 67 days, after which they were injected with either the viral mimic polyriboinosinic polyribocytidylic acid (pIC), or phosphate-buffered saline (PBS) as a control. Microarray analysis was used to determine the effect of the diet on the basal spleen transcriptome (pre-injection), and on the response to pIC (24 h post-injection). No marked differences in the spleen transcriptome were found between the three diets, either before or after injection with pIC. All fish, regardless of diet, showed a strong anti-viral response 24 h after pIC injection, with more than 500 genes having a significant difference of expression of 2-fold or higher compared to the PBS-injected fish for the FO, 40CO and 80CO diets. Gene Ontology annotation analysis of the three pIC-responsive gene lists indicated they were highly similar, and that the term 'immune system process' was significantly enriched in the pIC-responsive gene lists for all three diets. QPCR analysis for 5 genes with a known function in the anti-viral innate immune response (LGP2, STAT1, IRF1, ISG15 and viperin) showed modestly (smaller than 2-fold) up-regulated basal expression of LGP2, IRF1 and STAT1 in fish fed 40CO compared to the other diets. After pIC injection, all 5 genes were significantly and strongly up-regulated in pIC-injected fish compared to PBS-injected fish, but no significant differences were found between any of the diets. In conclusion, replacement of up to 80% of fish oil with camelina oil in Atlantic cod diets does not have a strong effect on basal spleen gene expression. Atlantic cod fed on camelina oil-containing diets are capable of mounting a strong anti-viral immune response, which is comparable to that in cod fed with a fish oil diet.

Old world monkeys and new age science: the evolution of nonhuman primate systems virology

Nonhuman primate (NHP) biomedical models are critical to our understanding of human health and disease, yet we are still in the early stages of developing sufficient tools to support primate genomic research that allow us to better understand the basis of phenotypic traits in NHP models of disease. A mere 7 years ago, the limited NHP transcriptome profiling that was being performed was done using complementary DNA arrays based on human genome sequences, and the lack of NHP genomic information and immunologic reagents precluded the use of NHPs in functional genomic studies. Since then, significant strides have been made in developing genomics capabilities for NHP research, from the rhesus macaque genome sequencing project to the construction of the first macaque-specific high-density oligonucleotide microarray, paving the way for further resource development and additional primate sequencing projects. Complete published draft genome sequences are now available for the chimpanzee ( Chimpanzee Sequencing Analysis Consortium 2005), bonobo ( Prufer et al. 2012), gorilla ( Scally et al. 2012), and baboon ( Ensembl.org 2013), along with the recently completed draft genomes for the cynomolgus macaque and Chinese rhesus macaque. Against this backdrop of both expanding sequence data and the early application of sequence-derived DNA microarrays tools, we will contextualize the development of these community resources and their application to infectious disease research through a literature review of NHP models of acquired immune deficiency syndrome and models of respiratory virus infection. In particular, we will review the use of -omics approaches in studies of simian immunodeficiency virus and respiratory virus pathogenesis and vaccine development, emphasizing the acute and innate responses and the relationship of these to the course of disease and to the evolution of adaptive immunity.

Transcriptome profiling reveals insight into distinct immune responses to Aeromonas salmonicida in gill of two rainbow trout strains

The fish gills represent a crucial organ for the communication with the aquatic environment. Transcriptional changes in gills of two hatchery rainbow trout strains in response to injection with the potent pathogen Aeromonas salmonicida were detected by global gene expression profiling using a 4×44K oligonucleotide microarray. Emphasis was placed on "day 3 postinfection" representing a decisive time point for the resolution of inflammation. The comparison of features and pathways differentially regulated in branchial tissues revealed that the local breeding strain BORN and imported American rainbow trout apply common and specific immune strategies. In gills of infected BORN trout, we observed a dynamic regulation of genes controlling NF-κB pathways and the induction of factors promoting the development of myeloid cells, whereas an increased expression of lysozyme and immunoglobulin genes was obvious in gills of infected import trout. In order to prove the relevance of the array-predicted candidates as well as well-known immune genes for gill immunity, a subsequent in vitro experiment was conducted. Altogether, we uncovered dynamic but moderate changes in the expression of a broad range of immune-relevant features implying the gill's involvement in pathogen defense strategies.

Pyroptosis and adaptive immunity mechanisms are promptly engendered in mesenteric lymph-nodes during pig infections with Salmonella enterica serovar Typhimurium

In this study, we explored the transcriptional response and the morphological changes occurring in porcine mesenteric lymph-nodes (MLN) along a time course of 1, 2 and 6 days post infection (dpi) with Salmonella Typhimurium. Additionally, we analysed the expression of some Salmonella effectors in tissue to complete our view of the processes triggered in these organs upon infection. The results indicate that besides dampening apoptosis, swine take advantage of the flagellin and prgJ expression by Salmonella Typhimuriun to induce pyroptosis in MLN, preventing bacterial dissemination. Furthermore, cross-presentation of Salmonella antigens was inferred as a mechanism that results in a rapid clearance of pathogen by cytotoxic T cells. In summary, although the Salmonella Typhimurium strain employed in this study was able to express some of its major virulence effectors in porcine MLN, a combination of early innate and adaptive immunity mechanisms might overcome virulence strategies employed by the pathogen, enabling the host to protect itself against bacterial spread beyond gut-associated lymph-nodes. Interestingly, we deduced that clathrin-mediated endocytosis could contribute to mechanisms of pathogen virulence and/or host defence in MLN of Salmonella infected swine. Taken together, our results are useful for a better understanding of the critical protective mechanisms against Salmonella that occur in porcine MLN to prevent the spread of infection beyond the intestine.

Expression of inflammation-mediated cluster of genes as a new marker of canine mammary malignancy

Because canine mammary tumours constitute a serious clinical problem and there are no good prognostic markers (only histopathological variables are used), the aim of the presented study was to find new malignancy markers as well as to identify intracellular pathways and biological processes characteristic for canine mammary malignancy. We compared gene expression of the most malignant mammary tumours (poorly differentiated cancers of the 3rd grade of malignancy) with less malignant tumours (well differentiated cancers of the 1st grade of malignancy). The results of our study indicated that in dogs the number of tumour-infiltrating myeloid cells or expression of myeloid-specific antigens by cancer cells is related to the cancer progression and may constitute a new marker of malignancy, however further studies in this field are required.

High mortality of juvenile gilthead sea bream (Sparus aurata) from photobacteriosis is associated with alternative macrophage activation and anti-inflammatory response: results of gene expression profiling of early responses in the head kidney

The halophilic bacterium Photobacterium damselae subsp. piscicida (Phdp) represents a substantial health problem for several fish species in aquaculture. Bacteria that reside free and inside phagocytes cause acute and chronic forms of photobacteriosis. Infections of juveniles rapidly kill up to 90-100% fish. Factors underlying failure of the immune protection against bacteria remain largely unknown. The reported study used a transcriptomic approach to address this issue. Juvenile sea breams (0.5 g) were challenged by immersion in salt water containing 2.89 × 10(8) CFU of a virulent Phdp and the head kidney was sampled after 24- and 48-h. Analyses were performed using the second version of a 44 k oligonucleotide DNA microarray that represents 19,734 sea bream unique transcripts and covers diverse immune pathways. Expression changes of selected immune genes were validated with qPCR. Results suggested rapid recognition of the pathogen, as testified by up-regulation of lectins and antibacterial proteins (bactericidal permeability-increasing protein lectins, lysozyme, intracellular and extracellular proteases), chemokines and chemokine receptors. Increased expression of proteins involved in iron and heme metabolism also could be a response against bacteria that are dependent on iron. However, negative regulators of immune/inflammatory response were preponderant among the up-regulated genes. A remarkable finding was the increased expression of IL-10 in concert with up-regulation of arginase I and II and proteins of the polyamine biosynthesis pathway that diverts the arginine flux from the production of reactive nitrogen species. Such expression changes are characteristic for alternatively activated macrophages that do not develop acute inflammatory responses. Immune suppression can be induced by the host to reduce tissue damages or by the pathogen to evade host response.

Expression profile of immune-related genes in Lates calcarifer infected by Cryptocaryon irritans

Cryptocaryon irritans causes Cyptocaryonosis or white spot disease in a wide range of marine fish including Lates calcarifer (Asian seabass). However, the immune response of this fish to the parasite is still poorly understood. In this study, quantitative polymerase chain reaction (qPCR) was performed to assess the expression profile of immune-related genes in L. calcarifer infected by C. irritans. A total of 21 immune-related genes encoding various functions in the fish immune system were utilized for the qPCR analysis. The experiment was initiated with the infection of juvenile fish by exposure to theronts from 200 C. irritans cysts, and non-infected juvenile fish were used as controls. Spleen, liver, gills and kidney tissues were harvested at three days post-infection from control and infected fish. In addition, organs were also harvested on day-10 post-infection from fish that had been allowed to recover from day-4 up to day-10 post-infection. L. calcarifer exhibited pathological changes on day-3 post-infection with the characteristic presence of white spots on the entire fish body, excessive mucus production and formation of a flap over the fish eye. High quality total RNA was extracted from all tissues and qPCR was performed. The qPCR analysis on the cohort of 21 immune-related genes of the various organs harvested on day-3 post-infection demonstrated that most genes were induced significantly (p < 0.05) in all tissues, particularly liver (11/21 genes) and kidney (11/21). The expression profile demonstrated that induction of the MHC Class IIα gene was the highest compared to the other genes followed by serum amyloid A, CC chemokine and hepcidin-2 precursor genes. In fish that were allowed to recover from the C. irritans infection (10 days post-infection), expression of the immune-related genes was down-regulated to levels similar to the control fish. These results provide insights into the interaction between C. irritans and L. calcarifer and suggest that the innate immune system plays an important role in early defence against parasite infection allowing the fish to eventually recover from the infection.

Effects of chronic cortisol administration on global expression of GR and the liver transcriptome in Sparus aurata

The present work was designed to assess the effects of artificially increased high plasma cortisol levels induced by slow-release cortisol implants on the mRNA abundance of the glucocorticoid receptor (GR) in different organs of Sparus aurata (Gilthead sea bream), as well as to evaluate global transcriptional changes in the liver, using the Aquagenomics S. aurata oligo-nucleotide microarray technology. For that purpose, groups of fish were intraperitoneally injected with implants containing two different concentrations of cortisol (50 or 200 μg/g body weight). Blood and organs were sampled after 7 and 14 days of cortisol implantation. Only fish with 200 μg/g implants exhibited a significant rise in plasma cortisol. Thus, we evaluated the expression of the GR in different organs in these fish 7 and 14 days post-implantation. GR mRNA abundance was upregulated in head kidney and heart of fish at both sampling times. In liver and muscle, GR mRNA abundance was upregulated after 14 days, whereas in gills, the GR mRNA transcript was upregulated earlier, at day 7. These results suggest that increased plasma cortisol induced by a slow-release implant of cortisol mimics the overall effects of stress and affects the expression of GR mRNA in a time- and organ-specific manner. Data obtained with the Aquagenomics S. aurata oligo-nucleotide microarray allowed the identification of a total of 491 cortisol-responsive transcripts and highlight the strong intensity of transcriptional modulation in liver of fish implanted with cortisol after 7 days, in contrast to that observed at day 14. Transcriptional remodeling highlighted a significant activity in carbohydrate metabolism mainly in the gluconeogenic pathway linked to downregulation of inflammatory and immune response processes in implanted fish.

Early response of gene expression in the distal intestine of Atlantic salmon (Salmo salar L.) during the development of soybean meal induced enteritis

Plant products in general and soybeans in particular can challenge the function and health of the intestinal tract. Salmonids develop an intestinal inflammation when fed diets containing soybean meal (SBM) and certain other legume ingredients. In the present study a 44K oligonucleotide salmonid microarray, qPCR and histology were used to investigate early response mechanisms in the distal intestine of Atlantic salmon (Salmo salar L.) during the first week of oral exposure to a diet containing 20% extracted SBM. The distal intestine transcriptome was profiled on days 1, 2, 3, 5 and 7 and compared to a control group fed fishmeal as the sole protein source. Histological evaluation of the distal intestine revealed the first signs of inflammation on day 5. The most prominent gene expression changes were seen on days 3 and 5. Up-regulation in immune-related genes was observed during the first 5 days, including GTPase IMAP family members, NF-kB-related genes and regulators of T cell and B cell function. Many functional genes involved in lipid metabolism, proteolysis, transport, metabolism and detoxification were initially up-regulated on days 1-3, possibly as an attempt by the tissue to compensate for the initiating immune response. Cell repair and extracellular matrix remodeling genes were up-regulated (heparanase, collagenase) on days 3 and 5. Down regulation of genes related to endocytosis, exocytosis, detoxification, transporters and metabolic processes from day 5 indicated initiation of dysfunction of digestive and metabolic functions that may occur as a result of inflammation or as a response to the introduction of soybean meal in the diet. This is the first study conducting transcriptomic profiling to characterize early responses during the development of SBMIE. Switching Atlantic salmon from a fishmeal to a 20% SBM diet resulted in rapid changes to the intestinal transcriptome, indicating an immune reaction with subsequent impaired epithelial barrier function and other vital intestinal functions.

Virulence and antibiotic resistance of Escherichia coli isolated from rooks

With regard to antibiotic resistance studies in various model animals in the urban environment, the presented study focused on the rook, many behavioural and ecological aspects of which are important from an epidemiological point of view. A total of 130 Escherichia coli strains isolated from rook faeces during a two-year period (2011-2012) were investigated for antibiotic resistance and virulence. Resistance to ampicillin (60%) and streptomycin (40%) were the most frequent, followed by resistance to fluoroquinolones (ciprofloxacin-22% and enrofloxacin-24%), tetracycline (18%), cotrimoxazol (17%) and florfenicol (14%). Ceftiofur resistance occured in 10.7% of strains and cefquinom resistance in 1.5% of strains. Twenty-five E.coli strains with a higher level of MICs of cephalosporins (over 2mg/L of ceftazidime and ceftriaxon) and fluoroquinolones were selected for detection of betalactamase genes (CTX-M, CMY), plasmid-mediated quinolone resistance qnrS, integrase 1, and for APEC (avian pathogenic E.coli) virulence factors (iutA, cvaC, iss, tsh, ibeA, papC, kpsII). Genes of CTX-M1, CMY-2, integrase 1, papC, cvaC, iutA were detected in one strain of E.coli, and qnrS, integrase 1, iss, cvaC, tsh were detected in another E.coli. DNA microarray revealed the absence of verotoxin and enterotoxin genes and pathogenicity islands. The results show that rooks can serve as a reservoir of antibiotic-resistant E. coli with avian pathogenic virulence factors for the human population, and potentially transmit such E.coli over long distances.

Cryptocaryon irritans infection induces the acute phase response in Lates calcarifer: a transcriptomic perspective

Cryptocaryoniasis (also known as marine white spot disease) is mediated by Cryptocaryon irritans. This obligate ectoparasitic protozoan infects virtually all marine teleosts, which includes Lates calcarifer, a highly valuable aquaculture species. Little is known about L. calcarifer-C. irritans interactions. This study was undertaken to gain an informative snapshot of the L. calcarifer transcriptomic response over the course of C. irritans infection. An in-house fabricated cDNA microarray slides containing 3872 probes from L. calcarifer liver and spleen cDNA libraries were used as a tool to investigate the response of L. calcarifer to C. irritans infection. Juvenile fish were infected with parasites for four days, and total RNA was extracted from liver tissue, which was harvested daily. We compared the transcriptomes of C. irritans-infected liver to uninfected liver over an infection period of four days; the comparison was used to identify the genes with altered expression levels in response to C. irritans infection. The greatest number of infection-modulated genes was recorded at 2 and 3 days post-infection. These genes were mainly associated with the immune response and were associated in particular with the acute phase response. Acute phase proteins such as hepcidin, C-type lectin and serum amyloid A are among the highly modulated genes. Our results indicate that an induced acute phase response in L. calcarifer toward C. irritans infection is similar to the responses observed in bacterial infections of teleosts. This response demonstrates the importance of first line defenses in teleost innate immune responses against ectoparasite infection.

Transcriptomic analysis of the host response to early stage salmonid alphavirus (SAV-1) infection in Atlantic salmon Salmo salar L

Salmon pancreas disease, caused by salmonid alphavirus (SAV) of the family Togaviridae, is an economically important disease affecting farmed Atlantic salmon (Salmo salar L.) in Scotland, Norway, and Ireland. The virus causes characteristic lesions in the pancreas, heart, kidney and skeletal muscle of infected fish. The mechanisms responsible for the pathology and the immune responses elicited in infected Atlantic salmon are not fully understood. A microarray-based study was therefore performed to evaluate the host transcriptomic response during the early stages of an experimentally-induced SAV-1 infection. Atlantic salmon parr were injected intra-peritoneally with viral cell culture supernatant or cell culture supernatant without virus. RNA, extracted from head kidney sampled from infected and control fish at 1, 3 and 5 days post-injection (d.p.i.), was interrogated with the 17 k TRAITS/SGP cDNA microarray. The greatest number of significantly differentially expressed genes was recorded at 3 d.p.i., mainly associated with immune and defence mechanisms, including genes involved in interferon I pathways and Major Histocompatibility Complex Class I and II responses. Genes associated with apoptosis and cellular stress were also found to be differentially expressed between infected and uninfected individuals, as were genes involved in inhibiting viral attachment and replication. The microarray results were validated by follow-on analysis of eight genes by real-time PCR. The findings of the study reflect mechanisms used by the host to protect itself during the early stages of SAV-1 infection. In particular, there was evidence of rapid induction of interferon-mediated responses similar to those seen during mammalian alphavirus infections, and also early involvement of an adaptive immune response. This study provides essential knowledge to assist in the development of effective control and management strategies for SAV-1 infection.

Pan-genomic analysis of bovine monocyte-derived macrophage gene expression in response to in vitro infection with Mycobacterium avium subspecies paratuberculosis

Mycobacterium avium subspecies paratuberculosis is the causative agent of Johne's disease, an intestinal disease of ruminants with major economic consequences. Infectious bacilli are phagocytosed by host macrophages upon exposure where they persist, resulting in lengthy subclinical phases of infection that can lead to immunopathology and disease dissemination. Consequently, analysis of the macrophage transcriptome in response to M. avium subsp. paratuberculosis infection can provide valuable insights into the molecular mechanisms that underlie Johne's disease. Here, we investigate pan-genomic gene expression in bovine monocyte-derived macrophages (MDM) purified from seven age-matched females, in response to in vitro infection with M. avium subsp. paratuberculosis (multiplicity of infection 2:1) at intervals of 2 hours, 6 hours and 24 hours post-infection (hpi). Differentially expressed genes were identified by comparing the transcriptomes of the infected MDM to the non-infected control MDM at each time point (adjusted P-value threshold ≤ 0.10). 1050 differentially expressed unique genes were identified 2 hpi, with 974 and 78 differentially expressed unique genes detected 6 and 24 hpi, respectively. Furthermore, in the infected MDM the number of upregulated genes exceeded the number of downregulated genes at each time point, with the fold-change in expression for the upregulated genes markedly higher than that for the downregulated genes. Inspection and systems biology analysis of the differentially expressed genes revealed an enrichment of genes involved in the inflammatory response, cell signalling pathways and apoptosis. The transcriptional changes associated with cellular signalling and the inflammatory response may reflect different immuno-modulatory mechanisms that underlie host-pathogen interactions during infection.

Effect of ovarian stimulation on oocyte gene expression in cattle

The objective was to analyze the impact of follicle stimulating hormone (FSH, ovarian stimulation) on the transcriptome of in vivo bovine oocytes three times around the luteinizing hormone (LH) surge. In vivo bovine oocytes were collected 2 h pre-LH surge, 6 h post-LH surge, and 22 h post-LH surge in both naturally ovulating and superovulated animals. To assess potential changes in gene levels, samples were hybridized using a custom bovine microarray. Two series of hybridizations were performed: the first comparing natural vs. stimulated cycles, the second according to time of collection. Among the potential candidates, 13 genes were selected according to their degree of differential expression and their potential link to oocyte competence. Measurements of their relative mRNA levels was made using QPCR. Gene candidates BTG4 (P = 0.0006), PTTG1 (P = 0.0027), PAPOLA (P = 0.0245), and LEO1 (P = 0.0393) had higher mRNA levels in oocytes treated with FSH for all collection times when compared to oocytes produced through the natural cycle. Among our selected candidates, only one gene, GDF9 (P = 0.0261), was present at a higher level in oocytes collected at -2 h and 6 h than 22 h post-LH for all treatments, regardless of the presence of FSH. Although the number of genes influenced by ovarian stimulation seemed low, the observed differences occurred at a time of minimal transcriptional activity and supported the potential impact on the future embryo. These impacts could have been epigenetic in nature, as embryo quality was not reported to be different from stimulated animals.

Chromosome analyses in dogs

Cytogenetics is the study of normal and abnormal chromosomes. Every species is characterized by a given number of chromosomes that can be recognized by their specific shape. The chromosomes are arranged according to standard classification schemes for the respective species. While pre- and postnatal chromosome analyses investigate the constitutional karyotype, tumor cytogenetics is focused on the detection of clonal acquired, tumor-associated chromosome aberrations. Cytogenetic investigations in dogs are of great value especially for breeders dealing with fertility problems within their pedigrees, for veterinarians and last but not least for the dog owners. Dogs and humans share a variety of genetic diseases, including cancer. Thus, the dog has become an increasingly important model for genetic diseases. However, cytogenetic analyses of canine cells are complicated by the complex karyotype of the dog. Only just 15 years ago, a standard classification scheme for the complete canine karyotype was established. For chromosome analyses of canine cells the same steps of chromosome preparation are used as in human cytogenetics. There are few reports about cytogenetic changes in non-neoplastic cells, involving predominantly the sex chromosomes. Cytogenetic analyses of different entities of canine tumors revealed that, comparable to human tumors, tumors of the dog are often characterized by clonal chromosome aberrations, which might be used as diagnostic and prognostic markers. The integration of modern techniques (molecular genetic approaches, adaptive computer programs) will facilitate and complete conventional cytogenetic studies. However, conventional cytogenetics is still non-replaceable.