A 7872 cDNA microarray and its use in bovine functional genomics

Newly identified interferon tau-responsive Hes family BHLH transcription factor 4 and cytidine/uridine monophosphate kinase 2 genes in peripheral blood granulocytes during early pregnancy in cows

In cattle, interferon-stimulated genes (ISGs) such as ISG15, MX1, MX2, and OAS1 are known as classic ISGs that are highly involved in the implantation process. Various molecules play a crucial role in the mechanisms underlying ISG effects. Although microarray analyses have highlighted the expression of various molecules during the implantation period, these molecules remain incompletely characterized. In the present study, various specifically expressed genes were selected and their characteristics were examined. The microarray data from peripheral blood leukocytes derived from artificially inseminated cows and granulocytes obtained from embryo-transferred cows, respectively, were used to identify new ISG candidates. Seven common genes, including ISG15 and OAS1, were confirmed, but only 4 of the 5 genes were amplified by reverse transcription quantitative polymerase chain reaction. In addition, 3 expressed sequence tags (ESTs) exhibited significantly greater expression in granulocytes from pregnant cows than that observed in bred nonpregnant cows, and the expression in granulocytes increased after interferon-tau stimulation. Sequence alignment revealed similar sequences within 2 ESTs on the Hairy and enhancer of split (Hes) family basic helix-loop-helix transcription factor 4 (HES4) gene. An additional EST was identified as cytidine/uridine monophosphate kinase 2 (CMPK2). In silico analysis facilitated the identification of transcription factor-binding sequences, including an interferon-stimulated response element and interferon regulatory factor-binding sites, within the promoter region of HES4 and CMPK2. These genes may function as new ISGs in the context of implantation and may participate in the coordination of the feto-maternal interface in cows.

Review: Enhancing gastrointestinal health in dairy cows

Due to their high energy requirements, high-yielding dairy cows receive high-grain diets. This commonly jeopardises their gastrointestinal health by causing subacute ruminal acidosis (SARA) and hindgut acidosis. These disorders can disrupt nutrient utilisations, impair the functionalities of gastrointestinal microbiota, and reduce the absorptive and barrier capacities of gastrointestinal epithelia. They can also trigger inflammatory responses. The symptoms of SARA are not only due to a depressed rumen pH. Hence, the diagnosis of this disorder based solely on reticulo-rumen pH values is inaccurate. An accurate diagnosis requires a combination of clinical examinations of cows, including blood, milk, urine and faeces parameters, as well as analyses of herd management and feed quality, including the dietary contents of NDF, starch and physical effective NDF. Grain-induced SARA increases acidity and shifts availabilities of substrates for microorganisms in the reticulo-rumen and hindgut and can result in a dysbiotic microbiota that are characterised by low richness, diversity and functionality. Also, amylolytic microorganisms become more dominant at the expense of proteolytic and fibrolytic ones. Opportunistic microorganisms can take advantage of newly available niches, which, combined with reduced functionalities of epithelia, can contribute to an overall reduction in nutrient utilisation and increasing endotoxins and pathogens in digesta and faeces. The reduced barrier function of epithelia increases translocation of these endotoxins and other immunogenic compounds out of the digestive tract, which may be the cause of inflammations. This needs to be confirmed by determining the toxicity of these compounds. Cows differ in their susceptibility to poor gastrointestinal health, due to variations in genetics, feeding history, diet adaptation, gastrointestinal microbiota, metabolic adaptation, stress and infections. These differences may also offer opportunities for the management of gastrointestinal health. Strategies to prevent SARA include balancing the diet for physical effective fibre, non-fibre carbohydrates and starch, managing the different fractions of non-fibre carbohydrates, and consideration of the type and processing of grain and forage digestibility. Gastrointestinal health disorders due to high grain feeding may be attenuated by a variety of feed supplements and additives, including buffers, antibiotics, probiotics/direct fed microbials and yeast products. However, the efficacy of strategies to prevent these disorders must be improved. This requires a better understanding of the mechanisms through which these strategies affect the functionality of gastrointestinal microbiota and epithelia, and the immunity, inflammation and 'gastrointestinal-health robustness' of cows. More representative models to induce SARA are also needed.

Skeletal muscle and liver gene expression profiles in finishing steers supplemented with Amaize

Our main objective was to evaluate the effects of feeding α-amylase (Amaize, Alltech Inc., Nicholasville, KY, USA) for 140 days on skeletal muscle and liver gene transcription in beef steers. Steers fed Amaize had lower average daily gain (p = .03) and gain:feed ratio (p = .05). No differences (p > .10) in serum metabolites or carcass traits were detected between the two groups but Amaize steers tended (p < .15) to have increased 12th rib fat depth. Microarray analysis of skeletal muscle revealed 21 differentially expressed genes (DEG), where 14 were up-regulated and seven were down-regulated in Amaize-fed steers. The bioinformatics analysis indicated that metabolic pathways involved in fat formation and deposition, stress response, and muscle function were activated, while myogenesis was inhibited in Amaize-fed steers. The quantitative PCR results for liver revealed a decrease (p < .01) in expression of fatty acid binding protein 1 (FABP1) and 3-hydroxybutyrate dehydrogenase 1 (BDH1) with Amaize. Because these genes are key for intracellular fatty acid transport, oxidation and ketone body production, data suggest a reduction in hepatic lipid catabolism. Future work to investigate potential positive effects of Amaize on cellular stress response, muscle function, and liver function in beef cattle appears warranted.

Microarray analysis of differential gene expression profiles in blood cells of naturally BLV-infected and uninfected Holstein-Friesian cows

The aim of the present study was to examine gene expression changes in response to bovine leukemia virus (BLV) infection, in an effort to determine genes that take a part in molecular events leading to persistent lymphocytosis (PL), and to better define genes involved in host response to BLV infection. Using bovine 70-mer oligonucleotide spotted microarrays (BLOPlus) and qRT-PCR validation, we studied global gene expression profiles in blood cells in vivo of 12 naturally BLV-infected Polish Holstein cows, and 12 BLV non-infected controls of the same breed and reared in herds with high BLV seroprevalence. With an arbitrary cut-off value of 1.5-fold change in gene expression, we identified the down-regulation of 212 genes (M value ≤-0.585) and the up-regulation of 158 genes (M value of ≥0.585) at 1% false discovery rate in BLV-positive animals in comparison to the BLV-negative group. The gene set enrichment analysis demonstrated that the differentially expressed (DE) genes could be classified to diverse biological processes, including immune response of host blood cells. Interestingly, our data indicated the potential involvement of the innate immunity, including complement system activation, NK-cell cytotoxicity and TREM-1 signaling, during the BLV-induced pathogenesis. We showed the occurrence of numerous regulatory processes that are targeted by BLV-infection. We also suggest that a complex network of interrelated pathways is disturbed, causing the interruption of the control of B-cell proliferation and programmed cell death.

Systems Analysis of Early Host Gene Expression Provides Clues for Transient Mycobacterium avium ssp avium vs. Persistent Mycobacterium avium ssp paratuberculosis Intestinal Infections

It has long been a quest in ruminants to understand how two very similar mycobacterial species, Mycobacterium avium ssp. paratuberculosis (MAP) and Mycobacterium avium ssp. avium (MAA) lead to either a chronic persistent infection or a rapid-transient infection, respectively. Here, we hypothesized that when the host immune response is activated by MAP or MAA, the outcome of the infection depends on the early activation of signaling molecules and host temporal gene expression. To test our hypothesis, ligated jejuno-ileal loops including Peyer’s patches in neonatal calves were inoculated with PBS, MAP, or MAA. A temporal analysis of the host transcriptome profile was conducted at several times post-infection (0.5, 1, 2, 4, 8 and 12 hours). When comparing the transcriptional responses of calves infected with the MAA versus MAP, discordant patterns of mucosal expression were clearly evident, and the numbers of unique transcripts altered were moderately less for MAA-infected tissue than were mucosal tissues infected with the MAP. To interpret these complex data, changes in the gene expression were further analyzed by dynamic Bayesian analysis. Bayesian network modeling identified mechanistic genes, gene-to-gene relationships, pathways and Gene Ontologies (GO) biological processes that are involved in specific cell activation during infection. MAP and MAA had significant different pathway perturbation at 0.5 and 12 hours post inoculation. Inverse processes were observed between MAP and MAA response for epithelial cell proliferation, negative regulation of chemotaxis, cell-cell adhesion mediated by integrin and regulation of cytokine-mediated signaling. MAP inoculated tissue had significantly lower expression of phagocytosis receptors such as mannose receptor and complement receptors. This study reveals that perturbation of genes and cellular pathways during MAP infection resulted in host evasion by mucosal membrane barrier weakening to access entry in the ileum, inhibition of Ca signaling associated with decreased phagosome-lysosome fusion as well as phagocytosis inhibition, bias toward Th2 cell immune response accompanied by cell recruitment, cell proliferation and cell differentiation; leading to persistent infection. Contrarily, MAA infection was related to cellular responses associated with activation of molecular pathways that release chemicals and cytokines involved with containment of infection and a strong bias toward Th1 immune response, resulting in a transient infection.

Effects of High Hydrostatic Pressure on Expression Profiles of In Vitro Produced Vitrified Bovine Blastocysts

High hydrostatic pressure (HHP) has been used to pre-condition embryos before essential, yet potentially detrimental procedures such as cryopreservation. However, the mechanisms for HHP are poorly understood. We treated bovine blastocysts with three different HHP (40, 60 and 80 MPa) in combination with three recovery periods (0, 1 h, 2 h post HHP). Re-expansion rates were significantly higher at 40 and 60 but lower at 80 MPa after vitrification-warming in the treated groups than controls. Microarray analysis revealed 399 differentially expressed transcripts, representing 254 unique genes, among different groups. Gene ontology analysis indicated that HHP at 40 and 60 MPa promoted embryo competence through down-regulation of genes in cell death and apoptosis, and up-regulation of genes in RNA processing, cellular growth and proliferation. In contrast, 80 MPa up-regulated genes in apoptosis, and down-regulated protein folding and cell cycle-related genes. Moreover, gene expression was also influenced by the length of the recovery time after HHP. The significantly over-represented categories were apoptosis and cell death in the 1 h group, and protein folding, response to unfolded protein and cell cycle in the 2 h group compared to 0 h. Taken together, HHP promotes competence of vitrified bovine blastocysts through modest transcriptional changes.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.

Systems biology analysis of gene expression during in vivo Mycobacterium avium paratuberculosis enteric colonization reveals role for immune tolerance

Survival and persistence of Mycobacterium avium subsp. paratuberculosis (MAP) in the intestinal mucosa is associated with host immune tolerance. However, the initial events during MAP interaction with its host that lead to pathogen survival, granulomatous inflammation, and clinical disease progression are poorly defined. We hypothesize that immune tolerance is initiated upon initial contact of MAP with the intestinal Peyer's patch. To test our hypothesis, ligated ileal loops in neonatal calves were infected with MAP. Intestinal tissue RNAs were collected (0.5, 1, 2, 4, 8 and 12 hrs post-infection), processed, and hybridized to bovine gene expression microarrays. By comparing the gene transcription responses of calves infected with the MAP, informative complex patterns of expression were clearly visible. To interpret these complex data, changes in the gene expression were further analyzed by dynamic Bayesian analysis, and genes were grouped into the specific pathways and gene ontology categories to create a holistic model. This model revealed three different phases of responses: i) early (30 min and 1 hr post-infection), ii) intermediate (2, 4 and 8 hrs post-infection), and iii) late (12 hrs post-infection). We describe here the data that include expression profiles for perturbed pathways, as well as, mechanistic genes (genes predicted to have regulatory influence) that are associated with immune tolerance. In the Early Phase of MAP infection, multiple pathways were initiated in response to MAP invasion via receptor mediated endocytosis and changes in intestinal permeability. During the Intermediate Phase, perturbed pathways involved the inflammatory responses, cytokine-cytokine receptor interaction, and cell-cell signaling. During the Late Phase of infection, gene responses associated with immune tolerance were initiated at the level of T-cell signaling. Our study provides evidence that MAP infection resulted in differentially regulated genes, perturbed pathways and specifically modified mechanistic genes contributing to the colonization of Peyer's patch.

Uncoupled embryonic and extra-embryonic tissues compromise blastocyst development after somatic cell nuclear transfer

Somatic cell nuclear transfer (SCNT) is the most efficient cell reprogramming technique available, especially when working with bovine species. Although SCNT blastocysts performed equally well or better than controls in the weeks following embryo transfer at Day 7, elongation and gastrulation defects were observed prior to implantation. To understand the developmental implications of embryonic/extra-embryonic interactions, the morphological and molecular features of elongating and gastrulating tissues were analysed. At Day 18, 30 SCNT conceptuses were compared to 20 controls (AI and IVP: 10 conceptuses each); one-half of the SCNT conceptuses appeared normal while the other half showed signs of atypical elongation and gastrulation. SCNT was also associated with a high incidence of discordance in embryonic and extra-embryonic patterns, as evidenced by morphological and molecular “uncoupling”. Elongation appeared to be secondarily affected; only 3 of 30 conceptuses had abnormally elongated shapes and there were very few differences in gene expression when they were compared to the controls. However, some of these differences could be linked to defects in microvilli formation or extracellular matrix composition and could thus impact extra-embryonic functions. In contrast to elongation, gastrulation stages included embryonic defects that likely affected the hypoblast, the epiblast, or the early stages of their differentiation. When taking into account SCNT conceptus somatic origin, i.e. the reprogramming efficiency of each bovine ear fibroblast (Low: 0029, Med: 7711, High: 5538), we found that embryonic abnormalities or severe embryonic/extra-embryonic uncoupling were more tightly correlated to embryo loss at implantation than were elongation defects. Alternatively, extra-embryonic differences between SCNT and control conceptuses at Day 18 were related to molecular plasticity (high efficiency/high plasticity) and subsequent pregnancy loss. Finally, because it alters re-differentiation processes in vivo, SCNT reprogramming highlights temporally and spatially restricted interactions among cells and tissues in a unique way.

Role of SPI-1 secreted effectors in acute bovine response to Salmonella enterica Serovar Typhimurium: a systems biology analysis approach

Salmonella enterica Serovar Typhimurium (S. Typhimurium) causes enterocolitis with diarrhea and polymorphonuclear cell (PMN) influx into the intestinal mucosa in humans and calves. The Salmonella Type III Secretion System (T3SS) encoded at Pathogenicity Island I translocates Salmonella effector proteins SipA, SopA, SopB, SopD, and SopE2 into epithelial cells and is required for induction of diarrhea. These effector proteins act together to induce intestinal fluid secretion and transcription of C-X-C chemokines, recruiting PMNs to the infection site. While individual molecular interactions of the effectors with cultured host cells have been characterized, their combined role in intestinal fluid secretion and inflammation is less understood. We hypothesized that comparison of the bovine intestinal mucosal response to wild type Salmonella and a SipA, SopABDE2 effector mutant relative to uninfected bovine ileum would reveal heretofore unidentified diarrhea-associated host cellular pathways. To determine the coordinated effects of these virulence factors, a bovine ligated ileal loop model was used to measure responses to wild type S. Typhimurium (WT) and a ΔsipA, sopABDE2 mutant (MUT) across 12 hours of infection using a bovine microarray. Data were analyzed using standard microarray analysis and a dynamic bayesian network modeling approach (DBN). Both analytical methods confirmed increased expression of immune response genes to Salmonella infection and novel gene expression. Gene expression changes mapped to 219 molecular interaction pathways and 1620 gene ontology groups. Bayesian network modeling identified effects of infection on several interrelated signaling pathways including MAPK, Phosphatidylinositol, mTOR, Calcium, Toll-like Receptor, CCR3, Wnt, TGF-β, and Regulation of Actin Cytoskeleton and Apoptosis that were used to model of host-pathogen interactions. Comparison of WT and MUT demonstrated significantly different patterns of host response at early time points of infection (15 minutes, 30 minutes and one hour) within phosphatidylinositol, CCR3, Wnt, and TGF-β signaling pathways and the regulation of actin cytoskeleton pathway.

Comparative analysis of the early transcriptome of Brucella abortus--infected monocyte-derived macrophages from cattle naturally resistant or susceptible to brucellosis

Brucellosis is a worldwide zoonotic infectious disease that has a significant economic impact on animal production and human public health. We characterized the gene expression profile of B. abortus-infected monocyte-derived macrophages (MDMs) from naïve cattle naturally resistant (R) or susceptible (S) to brucellosis using a cDNA microarray technology. Our data indicate that (1) B. abortus induced a slightly increased genome activation in R MDMs and a down-regulated transcriptome in S MDMs, during the onset of infection, (2) R MDMs had the ability to mount a type 1 immune response against B. abortus infection which was impaired in S cells, and (3) the host cell activity was not altered after 12 h post-B. abortus infection in R MDMs while the cell cycle was largely arrested in infected S MDMs at 12 h p.i. These results contribute to an improved understanding of how host responses may be manipulated to prevent infection by brucellae.

A small set of extra-embryonic genes defines a new landmark for bovine embryo staging

Axis specification in mouse is determined by a sequence of reciprocal interactions between embryonic and extra-embryonic tissues so that a few extra-embryonic genes appear as 'patterning' the embryo. Considering these interactions as essential, but lacking in most mammals the genetically driven approaches used in mouse and the corresponding patterning mutants, we examined whether a molecular signature originating from extra-embryonic tissues could relate to the developmental stage of the embryo proper and predict it. To this end, we have profiled bovine extra-embryonic tissues at peri-implantation stages, when gastrulation and early neurulation occur, and analysed the subsequent expression profiles through the use of predictive methods as previously reported for tumour classification. A set of six genes (CALM1, CPA3, CITED1, DLD, HNRNPDL, and TGFB3), half of which had not been previously associated with any extra-embryonic feature, appeared significantly discriminative and mainly dependent on embryonic tissues for its faithful expression. The predictive value of this set of genes for gastrulation and early neurulation stages, as assessed on naive samples, was remarkably high (93%). In silico connected to the bovine orthologues of the mouse patterning genes, this gene set is proposed as a new trait for embryo staging. As such, this will allow saving the bovine embryo proper for molecular or cellular studies. To us, it offers as well new perspectives for developmental phenotyping and modelling of embryonic/extra-embryonic co-differentiation.

Transferase activity function and system development process are critical in cattle embryo development

Microarray gene expression experiments often consider specific developmental stages, tissue sources, or reproductive technologies. This focus hinders the understanding of the cattle embryo transcriptome. To address this, four microarray experiments encompassing three developmental stages (7, 25, 280 days), two tissue sources (embryonic or extra-embryonic), and two reproductive technologies (artificial insemination or AI and somatic cell nuclear transfer or NT) were combined using two sets of meta-analyses. The first set of meta-analyses uncovered 434 genes differentially expressed between AI and NT (regardless of stage or source) that were not detected by the individual-experiment analyses. The molecular function of transferase activity was enriched among these genes that included ECE2, SLC22A1, and a gene similar to CAMK2D. Gene POLG2 was over-expressed in AI versus NT 7-day embryos and was under-expressed in AI versus NT 25-day embryos. Gene HAND2 was over-expressed in AI versus NT extra-embryonic samples at 280 days yet under-expressed in AI versus NT embryonic samples at 7 days. The second set of meta-analyses uncovered enrichment of system, organ, and anatomical structure development among the genes differentially expressed between 7- and 25-day embryos from either reproductive technology. Genes PRDX1and SLC16A1 were over-expressed in 7- versus 25-day AI embryos and under-expressed in 7- versus 25-day NT embryos. Changes in stage were associated with high number of differentially expressed genes, followed by technology and source. Genes with transferase activity may hold a clue to the differences in efficiency between reproductive technologies.

Functional annotation of novel lineage-specific genes using co-expression and promoter analysis

Background

The diversity of placental architectures within and among mammalian orders is believed to be the result of adaptive evolution. Although, the genetic basis for these differences is unknown, some may arise from rapidly diverging and lineage-specific genes. Previously, we identified 91 novel lineage-specific transcripts (LSTs) from a cow term-placenta cDNA library, which are excellent candidates for adaptive placental functions acquired by the ruminant lineage. The aim of the present study was to infer functions of previously uncharacterized lineage-specific genes (LSGs) using co-expression, promoter, pathway and network analysis.

Results

Clusters of co-expressed genes preferentially expressed in liver, placenta and thymus were found using 49 previously uncharacterized LSTs as seeds. Over-represented composite transcription factor binding sites (TFBS) in promoters of clustered LSGs and known genes were then identified computationally. Functions were inferred for nine previously uncharacterized LSGs using co-expression analysis and pathway analysis tools. Our results predict that these LSGs may function in cell signaling, glycerophospholipid/fatty acid metabolism, protein trafficking, regulatory processes in the nucleus, and processes that initiate parturition and immune system development.

Conclusions

The placenta is a rich source of lineage-specific genes that function in the adaptive evolution of placental architecture and functions. We have shown that co-expression, promoter, and gene network analyses are useful methods to infer functions of LSGs with heretofore unknown functions. Our results indicate that many LSGs are involved in cellular recognition and developmental processes. Furthermore, they provide guidance for experimental approaches to validate the functions of LSGs and to study their evolution.

Differential gene expression in bovine elongated (Day 17) embryos produced by somatic cell nucleus transfer and in vitro fertilization

Somatic cloning in cattle is associated with impaired embryo development, caused by inappropriate epigenetic reprogramming during embryogenesis; however, there is a paucity of data regarding gene expression at the critical elongation and peri-implantation stages. The objective of the present study was to identify genes differentially expressed in bovine cloned embryos at Day 17 of development (Day 0=day of nucleus transfer or IVF). Day 7 blastocysts (Hand Made Cloned or IVP) were transferred to recipient cattle and collected at Day 17. The efficiency of recovery of elongated embryos was similar, however cloned embryos elongated less than IVP embryos (91.8+/-45.8 vs. 174+/-50mm) and fewer had embryonic discs (63 vs. 83%). Qualitative and quantitative PCR detected expression of OCT4, NANOG, IFNtau, EOMES, FGF4, SOX2, and CDX2 in all IVP embryos. In most cloned embryos, NANOG and FGF4 were absent (verified by qPCR); NANOG, EOMES, and FGF4 were underexpressed, whereas IFNtau was overexpressed in cloned embryos. Based on qPCRs, other genes, i.e., SPARC, SNRB1, and CBPP22, were down-regulated in cloned embryos, whereas HSP70 and TDKP1 were overexpressed. In bovine microarrays, 47 genes (3.6%) were deregulated in cloned embryos, including several involved in trophoblast growth and differentiation. In conclusion, we inferred that these data were indicative of incomplete epigenetic reprogramming after cloning; this could lead to aberrant gene expression and subsequently early pregnancy loss. There was an apparent association between incomplete morphological elongation and aberrant reprogramming of a subset of genes critical for early embryonic development.

Endometrium as an early sensor of in vitro embryo manipulation technologies

Implantation is crucial for placental development that will subsequently impact fetal growth and pregnancy success with consequences on postnatal health. We postulated that the pattern of genes expressed by the endometrium when the embryo becomes attached to the mother uterus could account for the final outcome of a pregnancy. As a model, we used the bovine species where the embryo becomes progressively and permanently attached to the endometrium from day 20 of gestation onwards. At that stage, we compared the endometrial genes profiles in the presence of an in vivo fertilized embryo (AI) with the endometrial patterns obtained in the presence of nuclear transfer (SCNT) or in vitro fertilized embryos (IVF), both displaying lower and different potentials for term development. Our data provide evidence that the endometrium can be considered as a biological sensor able to fine-tune its physiology in response to the presence of embryos whose development will become altered much later after the implantation process. Compared with AI, numerous biological functions and several canonical pathways with a major impact on metabolism and immune function were found to be significantly altered in the endometrium of SCNT pregnancies at implantation, whereas the differences were less pronounced with IVF embryos. Determining the limits of the endometrial plasticity at the onset of implantation should bring new insights on the contribution of the maternal environment to the development of an embryo and the success of pregnancy.

Embryonic gene expression profiling using microarray analysis

Microarray technology enables the interrogation of thousands of genes at one time and therefore a systems level of analysis. Recent advances in the amplification of RNA, genome sequencing and annotation, and the lower cost of developing microarrays or purchasing them commercially, have facilitated the analysis of single preimplantation embryos. The present review discusses the components of embryonic expression profiling and examines current research that has used microarrays to study the effects of in vitro production and nuclear transfer.

Biological interpretations of transcriptomic profiles in mammalian oocytes and embryos

The characterization of gene-expression profiles in oocytes and embryos is critical to understand the influence of genetic and environmental factors on preimplantation and fetal development. Numerous gene-expression microarray studies using different platforms and species are offering insights into the biological processes extensively represented among the genes exhibiting differential expression. Major advances on understanding the direct relationship between gene expression and developmental competence are being reported. Integration of information across studies using meta-analysis techniques can increase the precision and accuracy to identify expression profiles associated with embryo development. Gene network and pathway analyses are offering insights into gene interactions and expression profiles of embryos. All these advances are cementing the way toward a comparative and systems approach to understanding the complex processes underlying vertebrate development.

Comparative gene expression analysis of bovine nuclear-transferred embryos with different developmental potential by cDNA microarray and real-time PCR to determine genes that might reflect calf normality

Placental abnormalities are the main factor in the high incidence of somatic cell clone abnormalities. The expression of several trophoblast cell-specific molecules is enhanced during gestational days 7 to 14. To determine the possible genes whose expression patterns might reflect calf normality, we first compared the gene expression profiles on day 15 between in vitro-fertilized (IVF) embryos and two types of somatic cell nuclear-transferred embryos with either a high (FNT) or low (CNT) incidence of neonatal abnormalities using a cDNA microarray containing 16 of 21 placenta-specific genes developed from tissues collected across gestation. To identify significant genes from the screening of day 15 embryos, genes with a less than two-fold difference in expression between IVF and CNT embryos, and those with a greater than two-fold difference between IVF and FNT and between CNT and FNT were considered to contribute to clone abnormalities. These two comparisons revealed 18 down-regulated and 18 upregulated genes of the 1722 genes examined. We then examined the expression levels of 10 genes with known functions in eight-cell and blastocyst-stage embryos by real-time PCR. The mRNA expression pattern of interferon (IFN)-tau, a trophectoderm-related gene, differed between IVF, CNT, and FNT eight-cell embryos; few or none of the IVF or CNT eight-cell embryos expressed IFN-tau mRNA, but all eight-cell FNT embryos expressed IFN-tau. IFN-tau mRNA expression was significantly higher in IVF blastocysts, however, than in nuclear-transferred blastocysts. Average IFN-tau mRNA expression in FNT blastocysts was not different from that in CNT blastocysts, due to one CNT blastocyst with high expression. The precise relation between early expression of IFN-tau mRNA and inferior developmental potential in cloned embryos should be examined further.

Aberrant gene expression patterns in placentomes are associated with phenotypically normal and abnormal cattle cloned by somatic cell nuclear transfer

Transcription profiling of placentomes derived from somatic cell nuclear transfer (SCNT, n = 20), in vitro fertilization (IVF, n = 9), and artificial insemination (AI, n = 9) at or near term development was performed to better understand why SCNT and IVF often result in placental defects, hydrops, and large offspring syndrome (LOS). Multivariate analysis of variance was used to distinguish the effects of SCNT, IVF, and AI on gene expression, taking into account the effects of parturition (term or preterm), sex of fetus, breed of dam, breed of fetus, and pathological finding in the offspring (hydrops, normal, or other abnormalities). Differential expression of 20 physiologically important genes was confirmed with quantitative PCR. The largest effect on placentome gene expression was attributable to whether placentas were collected at term or preterm (i.e., whether the collection was because of disease or to obtain stage-matched controls) followed by placentome source (AI, IVF, or SCNT). Gene expression in SCNT placentomes was dramatically different from AI ( n = 336 genes; 276 >2-fold) and from IVF ( n = 733 genes; 162 >2-fold) placentomes. Functional analysis of differentially expressed genes (DEG) showed that IVF has significant effects on genes associated with cellular metabolism. In contrast, DEG associated with SCNT are involved in multiple pathways, including cell cycle, cell death, and gene expression. Many DEG were shared between the gene lists for IVF and SCNT comparisons, suggesting that common pathways are affected by the embryo culture methods used for IVF and SCNT. However, the many unique gene functions and pathways affected by SCNT suggest that cloned fetuses may be starved and accumulating toxic wastes due to placental insufficiency caused by reprogramming errors. Many of these genes are candidates for hydrops and LOS.

Identification of differentially expressed genes in distinct skeletal muscles in cattle using cDNA microarray

The 788-gene microarray was manufactured using selected elements from three different cDNA libraries in order to identify molecular processes that determine phenotypic characteristics between loin (M. longissimus thoracis) and round (M. semimembranosus) muscles. Microarray analyses identified 24 differentially expressed genes between the two muscles investigated. Five of the genes were verified by quantitative RT-PCR and three of them were mapped on bovine chromosomes using 5,000 rad bovine radiation hybrid (RH) panel. The map locations indicated that they were mapped in the same chromosomal regions where IMF and growth QTLs were located, suggesting that they are most possible positional candidate genes for the traits.

Nutrition-induced ketosis alters metabolic and signaling gene networks in liver of periparturient dairy cows

Dairy cows are highly susceptible after parturition to developing liver lipidosis and ketosis, which are costly diseases to farmers. A bovine microarray platform consisting of 13,257-annotated oligonucleotides was used to study hepatic gene networks underlying nutrition-induced ketosis. On day 5 postpartum, 14 Holstein cows were randomly assigned to ketosis-induction (n = 7) or control (n = 7) groups. Cows in the ketosis-induction group were fed at 50% of day 4 intake until they developed signs of clinical ketosis, and cows in the control group were fed ad libitum throughout the treatment period. Liver was biopsied at 10-14 (ketosis) or 14 days postpartum (controls). Feed restriction increased blood concentrations of nonesterified fatty acids and beta-hydroxybutyrate, but decreased glucose. Liver triacylglycerol concentration also increased. A total of 2,415 genes were altered by ketosis (false discovery rate = 0.05). Ingenuity Pathway Analysis revealed downregulation of genes associated with oxidative phosphorylation, protein ubiquitination, and ubiquinone biosynthesis with ketosis. Other molecular adaptations included upregulation of genes and nuclear receptors associated with cytokine signaling, fatty acid uptake/transport, and fatty acid oxidation. Genes downregulated during ketosis included several associated with cholesterol metabolism, growth hormone signaling, proton transport, and fatty acid desaturation. Feed restriction and ketosis resulted in previously unrecognized alterations in gene network expression underlying key cellular functions and discrete metabolic events. These responses might help explain well-documented physiological adaptations to reduced feed intake in early postpartum cows and, thus, provide molecular targets that might be useful in prevention and treatment of liver lipidosis and ketosis.

Technical Note: Bovine Oviduct and Endometrium Array Version 1: A Tailored Tool for Studying Bovine Endometrium Biology and Pathophysiology

Fertility problems are the main reason for slaughter of high-performance milk cows, because elongated calving intervals result in financial losses for the farmer and retard genetic progress. Genetic improvement of fertility would be of great benefit, but functional traits for effective selection are missing. Recent advances in functional genomics tools like DNA microarrays could be the key to identify gene expression patterns in the endometrium that correlate with maternal fertility. Therefore, a first version of a bovine oviduct and endometrium cDNA array was established that contains a set of 1,440 cDNA clones and long oligonucleotides representing 950 different genes. The major part of these genes has been identified in a series of differential gene expression studies in endometrium (different stages of the estrous cycle, d 18 pregnant vs. nonpregnant) and in oviduct epithelial cells (different stages of the estrous cycle) using a combination of subtracted cDNA libraries and cDNA array hybridization. Furthermore, cDNA clones of genes, which showed no changes in their mRNA levels in the analyzed tissues, were added as controls. Reproducibility of the array hybridization, a comparison with the Affymetrix bovine genome array, and confirmation of differential gene expression with reverse transcription-quantitative PCR is shown. Potential future applications include systematic studies of interactions between metabolic status and functionality of the endometrium to identify genes that could be used for differential diagnosis of fertility problems. Further, endometrium transcriptome profiles may serve as novel traits to improve fertility by genetic selection.

Discovery and characterization of 91 novel transcripts expressed in cattle placenta

Background

Among the eutherian mammals, placental architecture varies to a greater extent than any other tissue. The diversity of placental types, even within a single mammalian order suggests that genes expressed in placenta are under strong Darwinian selection. Thus, the ruminant placenta may be a rich source of genes to explore adaptive evolutionary responses in mammals. The aim of our study was to identify novel transcripts expressed in ruminant placenta, and to characterize them with respect to their expression patterns, organization of coding sequences in the genome, and potential functions.

Results

A combination of bioinformatics, comparative genomics and transcript profiling was used to identify and characterize 91 novel transcripts (NTs) represented in a cattle placenta cDNA library. These NTs have no significant similarity to any non-ferungulate DNA or RNA sequence. Proteins longer than 100 aa were predicted for 29 NTs, and 21 are candidate non-coding RNAs. Eighty-six NTs were found to be expressed in one or more of 18 different tissues, with 39 (42%) showing tissue-preference, including six that were expressed exclusively in placentome. The authenticity of the NTs was confirmed by their alignment to cattle genome sequence, 42 of which showed evidence of mRNA splicing. Analysis of the genomic context where NT genes reside revealed 61 to be in intergenic regions, whereas 30 are within introns of known genes. The genes encoding the NTs were found to be significantly associated with subtelomeric regions.

Conclusion

The 91 lineage-specific transcripts are a useful resource for studying adaptive evolutionary responses of the ruminant placenta. The presence of so many genes encoding NTs in cattle but not primates or rodents suggests that gene loss and gain are important mechanisms of genome evolution in mammals. Furthermore, the clustering of NT genes within subtelomeric regions suggests that such regions are highly dynamic and may foster the birth of novel genes. The sequencing of additional vertebrate genomes with defined phylogenetic relationships will permit the search for lineage-specific genes to take on a more evolutionary context that is required to understand their origins and functions.

Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle

Achieving greater understanding of the genomic influence on milk synthesis in dairy cows represents a daunting challenge. Bovine-specific microarrays have allowed for high-throughput gene expression analysis of the mammary transcriptome. However, real-time PCR (qPCR) still represents the method of choice for accurate expression profiling of small numbers of genes and verification of key microarray relationships. This method is extremely sensitive but requires data normalization to account for analytical errors. Ideally, expression of genes used as internal controls should not be affected by specific treatments or physiological state. Mammary biopsies were collected from five cows each at −15, 1, 15, 30, 60, 120, and 240 days relative to parturition for gene expression profiling. We evaluated expression of nine genes ( RPS9, ACTB, GAPD, GTP, ITGB4BP, MRPL39, RPS23, RPS15, and UXT) that could serve as internal controls in mammary tissue using qPCR. Due to gradual increases in mammary RNA concentration (μg/mg tissue) over lactation, all genes investigated experienced a dilution effect. We used pairwise comparison of expression ratios to analyze the reliability of these genes as internal controls. UXT, RPS9, and RPS15 had the most stable expression ratios across cow and time. We also assessed co-regulation among genes through network analysis. Network analysis suggested co-regulation among most of the genes examined, with MYC playing a central role. Pairwise comparison was suitable for finding appropriate internal controls in mammary gland tissue. Results showed that the geometrical average of UXT, RPS9, and RPS15 expression could be used as internal control for longitudinal mammary gene expression profiling.

Analysis of uteroplacental-specific molecules and their functions during implantation and placentation in the bovine

In cattle, the mechanisms underlying implantation and placental development are still unclear. Synepitheliochorial placentation in cattle is noninvasive, and thus generates limited interest in terms of degradation and remodeling of endometrial tissues. The overall purpose of this study was three-fold: (1) to examine the gene circuitry around the implantation window, (2) to understand development of the placenta during the peri-implantation period by using a uteroplacental cDNA microarray, and (3) to study the roles of molecules involved in endometrial remodeling. Bovine trophoblastic binucleate cell-specific molecules, such as pregnancy-associated glycoproteins (PAGs), placental lactogen (PL), and prolactin-related proteins (PRPs), were markedly expressed in binucleate cells (BNCs) around implantation. The expression of PRP-1 was specific to the caruncular (CAR) area of the gravid uterine horn. Gelatinases (MMP-2 and -9) in association with heparanase may be central to endometrial remodeling. In situ hybridization analyses of PAGs, PRPs, PL, and heparanase suggested that BNCs expressed these molecules simultaneously. Future studies will further investigate the specific roles of these molecules in placentogenesis. The uteroplacental cDNA microarray presented cascades of molecular signatures not only for the endometrium but also for the intricate dialogue at the level of the feto-maternal interface in cattle. Placentome morphogenesis potentially parallels the dynamic multigenic circuitry and regulates the cell cycle in the endometrium. The roles of BNCs and their secreted molecules remain an enigma, particularly with regard to the adhesion process and endometrial remodeling, which is the focus of this study.

Gene expression and maintenance of pregnancy in bovine: roles of trophoblastic binucleate cell-specific molecules

Cell to cell interaction plays a pivotal role in the regulation of placentogenesis and exchange of stage-specific developmental signals between the fetal and maternal units. Specifically, these interactions are paramount for programmed fetal growth, maternal adaptation to pregnancy and coordination of parturition. However, little is known about the precise regulation of placentation and maintenance of gestation in cattle. Therefore, the aim of the present study was to decipher the complex networks of cell communication to gain an insight into the multifaceted developmental process and understand the profound consequences of flawed communication. In the ruminant, the binucleate cell plays a central role in forming the structures and secretions at the fetomaternal interface that are crucial in establishing and maintaining pregnancy. Herein, we summarise differences in the abundance of specific RNA transcripts in the bovine cotyledon and caruncle using global gene expression profiling and further investigate the relationship of mRNA abundance for selected pregnancy-specific genes of interest (identified from microarray studies) that are localised exclusively to the binucleate cell, such as placental lactogen, prolactin-related proteins and pregnancy-associated glycoproteins. The results suggest that a well-orchestrated transcriptional command from binucleate cells is pivotal to the establishment and progression of pregnancy in cattle.

Plane of nutrition prepartum alters hepatic gene expression and function in dairy cows as assessed by longitudinal transcript and metabolic profiling

Liver metabolism and health in dairy cows during the periparturient period are affected by plane of nutrition prepartum. Long-term adaptations in hepatic gene expression are important for complete understanding of liver function. We examined temporal gene expression profiles during the dry period and early lactation in liver of Holstein cows fed moderate dietary energy ad libitum or restricted during the entire dry period using a microarray consisting of 7,872 annotated cattle cDNA inserts and quantitative RT-PCR. We identified 85 genes with expression patterns that were affected by level of energy intake prepartum over time. Restricted energy intake prepartum resulted in more pronounced upregulation of genes with key functions in hepatic fatty acid oxidation ( CPT1A, ADIPOR2), gluconeogenesis ( PC), and cholesterol synthesis ( SC4MOL). Ad libitum feeding upregulated a number of genes associated with liver triacylglycerol synthesis ( DGAT1) and proinflammatory cytokines ( TNFAIP3). Genomic responses to ad libitum feeding were accompanied by increased incorporation of palmitate to esterified products in vitro and increased liver triacylglycerol concentration in vivo. Overall, gene expression profiles due to plane of nutrition prepartum partly explained differences in rates of liver palmitate metabolism, blood serum metabolite concentrations, and liver tissue triacylglycerol concentration. Our data show that moderate overfeeding of energy in the dry period, in the absence of obesity, results in transcriptional changes predisposing cows to fatty liver and perhaps compromising overall liver health during the periparturient period. In this context, controlled energy intake may confer an advantage to the cow by triggering hepatic molecular adaptations well ahead of parturition.

Discovery of eight novel divergent homologs expressed in cattle placenta

Ten divergent homologs were identified using a subtractive bioinformatic analysis of 12,614 cattle placenta expressed sequence tags followed by comparative, evolutionary, and gene expression studies. Among the 10 divergent homologs, 8 have not been identified previously. These were named as follows: cattle cerebrum and skeletal muscle-specific transcript 1 (CSSMST1), cattle intestine-specific transcript 1 (CIST1), hepatitis A virus cellular receptor 1 amino-terminal domain-containing protein (HAVCRNDP), prolactin-related proteins 8, 9, and 11 (PRP8, PRP9, and PRP11, respectively) and secreted and transmembrane protein 1A and 1B (SECTM1A and SECTM1B, respectively). In addition, two previously known divergent genes were identified, trophoblast Kunitz domain protein 1 (TKDP1) and a new splice variant of TKDP4. Nucleotide substitution analysis provided evidence for positive selection in members of the PRP gene family, SECTM1A and SECTM1B. Gene expression profiles, motif predictions, and annotations of homologous sequences indicate immunological and reproductive functions of the divergent homologs. The genes identified in this study are thus of evolutionary and physiological importance and may have a role in placental adaptations.

Emerging technologies for identifying superior dairy cows in New Zealand

The performance of animals is determined by the interaction of their genes with environmental circumstances. Accordingly, animals exhibiting superior performance are not necessarily the animals with the best genes nor are they the best choice of parents. Statistical analyses of production records for repeated traits, e.g. lactation yields and reproductive performance, show that part of the variation in performance among animals in the same herd and year is due to genetic differences, and the remainder is due to so-called residual or environmental factors that are not passed on to offspring. These within-herd environmental factors can be partitioned into a component that affects performance throughout an animal's lifetime, and a part that is unique to each observation. The process of animal evaluation from pedigree and performance records partitions the superiority of each cow into these three components. Reliable assessment of the genetic merit of bulls has required progeny testing, and for cows has required observation of their own individual performance. Selection on the genetic or breeding value component has systematically improved animal performance over recent decades, but has been limited by the age at which assessments of genetic merit are available. Emerging molecular technologies can read DNA sequences or measure RNA expression and have allowed the identification of a number of chromosome regions, and a few specific genes in those regions, that influence economic performance. This information allows better characterisation of the relationships between animals and more accurate predictions of genetic merit in bulls without progeny information and in cows that have yet to produce their own performance record. At some stage, enough genes responsible for variation in performance will be identified to allow faster genetic progress through selection of animals at young ages and therefore more rapid turnover of the generations. Mechanisms that modify gene expression have been identified and these may ultimately allow animals to be selected at an early age for lifetime productivity, accounting for processes that modify gene expression and lead to differences in performance that are not reflected by DNA sequence information. This review describes the status of these emerging technologies and their likely role in the improvement of dairy cattle.

Advances in livestock genomics: Opening the barn door

Genome research in animals used in agriculture has progressed rapidly in recent years, moving from rudimentary genome maps to trait maps to gene discovery. These advances are the result of animal genome projects following closely in the footsteps of the Human Genome Project, which has opened the door to genome research in farm animals. In return, genome research in livestock species is contributing to our understanding of chromosome evolution and to informing the human genome. Enhancement of these contributions plus the much anticipated application of DNA-based tools to animal health and production can be expected as livestock genomics enters its sequencing era.

Global gene expression profiles reveal significant nuclear reprogramming by the blastocyst stage after cloning

Nuclear transfer (NT) has potential applications in agriculture and biomedicine, but the technology is hindered by low efficiency. Global gene expression analysis of clones is important for the comprehensive study of nuclear reprogramming. Here, we compared global gene expression profiles of individual bovine NT blastocysts with their somatic donor cells and fertilized control embryos using cDNA microarray technology. The NT embryos' gene expression profiles were drastically different from those of their donor cells and closely resembled those of the naturally fertilized embryos. Our findings demonstrate that the NT embryos have undergone significant nuclear reprogramming by the blastocyst stage; however, problems may occur during redifferentiation for tissue genesis and organogenesis, and small reprogramming errors may be magnified downstream in development.

Temporal gene expression profiling of liver from periparturient dairy cows reveals complex adaptive mechanisms in hepatic function

Long-term molecular adaptations in liver from high-producing dairy cows are virtually unknown. Liver from five Holstein cows was biopsied at −65, −30, −14, +1, +14, +28, and +49 days relative to parturition for transcript profiling using a microarray consisting of 7,872 annotated cattle cDNA inserts. More than 5,000 cDNA elements represented on the microarray were expressed in liver. From this set we identified 62 differentially expressed genes related to physiological state, with a false discovery rate threshold of P = 0.20. The dominant expression pattern consisted of upregulation from day −30 through day +1, followed by downregulation through day +28. There was a threefold decrease from day −65 through day +14 in expression of IGFBP3, GSTM5, and PDPK1. These genes mediate IGF-I transport, oxidative stress, and glucose homeostasis, respectively. IGFBP3, EIF4B, and GSTM5 mRNA levels were positively correlated with blood serum total protein. Correlation analysis showed positive associations between serum nonesterified fatty acids and mRNA expression for SAA1, CPT1A, ACADVL, and TFAP2A. Transcript levels of ACSL1, PPARA, and TFAP2A were positively correlated with serum β-hydroxybutyrate. Expression patterns for certain genes (e.g., IGFBP3, HNF4A, GPAM) revealed adaptations commencing well ahead of parturition, suggesting they are regulated by factors other than periparturient hormonal environment. Results provide evidence that hepatic inflammatory responses occurring near parturition initiate or augment adipose catabolism. In this context, cytokines, acute-phase proteins, and serum nonesterified fatty acids are key players in periparturient cow metabolism. We propose a model for integrating gene expression, metabolite, and liver composition data to explain physiological events in placenta, adipose, and liver during the periparturient period.

Construction and validation of a Bovine Innate Immune Microarray

BACKGROUND

Microarray transcript profiling has the potential to illuminate the molecular processes that are involved in the responses of cattle to disease challenges. This knowledge may allow the development of strategies that exploit these genes to enhance resistance to disease in an individual or animal population.

RESULTS

The Bovine Innate Immune Microarray developed in this study consists of 1480 characterised genes identified by literature searches, 31 positive and negative control elements and 5376 cDNAs derived from subtracted and normalised libraries. The cDNA libraries were produced from 'challenged' bovine epithelial and leukocyte cells. The microarray was found to have a limit of detection of 1 pg/microg of total RNA and a mean slide-to-slide correlation co-efficient of 0.88. The profiles of differentially expressed genes from Concanavalin A (ConA) stimulated bovine peripheral blood lymphocytes were determined. Three distinct profiles highlighted 19 genes that were rapidly up-regulated within 30 minutes and returned to basal levels by 24 h; 76 genes that were up-regulated between 2-8 hours and sustained high levels of expression until 24 h and 10 genes that were down-regulated. Quantitative real-time RT-PCR on selected genes was used to confirm the results from the microarray analysis. The results indicate that there is a dynamic process involving gene activation and regulatory mechanisms re-establishing homeostasis in the ConA activated lymphocytes. The Bovine Innate Immune Microarray was also used to determine the cross-species hybridisation capabilities of an ovine PBL sample.

CONCLUSION

The Bovine Innate Immune Microarray has been developed which contains a set of well-characterised genes and anonymous cDNAs from a number of different bovine cell types. The microarray can be used to determine the gene expression profiles underlying innate immune responses in cattle and sheep.

The expanding role of microarrays in the investigation of macrophage responses to pathogens

In the last few years, microarray technology has emerged as the method of choice for large-scale gene expression studies. It provides an efficient and rapid method to investigate the entire transcriptome of a cell. No research field has benefited more from microarray technology than the study of the exquisite interplay between pathogens and hosts. Numerous microarray studies have now been published in this field, which have provided insights into the mechanisms of host defence and the tactics employed by pathogens to circumvent these protection strategies. These studies have led to a more comprehensive understanding of the host immune response and identified new avenues of research for potential control strategies against pathogens. In the past, research has concentrated on human and mouse microarrays to investigate host-pathogen interactions, regardless of the host species. This trend is changing with the ever-expanding sequence resources now available for many pathogen and host species, including livestock animals. The use of species-specific microarrays has furthered our understanding of host-pathogen interactions for particular organisms and aided in the annotation of unknown genes. Macrophages play a central role in the host's innate and adaptive immune responses to pathogens. These cells are in the first line of defence and interact with a wide range of pathogens; many of which have evolved strategies to circumvent the macrophage defence mechanisms and survive within these cells. In this report, we review the wealth of studies using microarray technology to investigate the response of macrophages to pathogens. These studies illustrate how microarray technology has expanded our understanding of the dialogue between macrophage and pathogen and provide examples of the benefits and pitfalls of using this technique. Furthermore, we discuss the resources available to use microarray analysis to study the immune response of a non-human, non-rodent species, the cow.