Genomic characterization of equine influenza A subtype H3N8 viruses by long read sequencing and functional analyses of the PB1-F2 virulence factor of A/equine/Paris/1/2018.

Equine influenza virus (EIV) remains a threat to horses, despite the availability of vaccines. Strategies to monitor the virus and prevent potential vaccine failure revolve around serological assays, RT-qPCR amplification, and sequencing the viral hemagglutinin (HA) and neuraminidase (NA) genes. These approaches overlook the contribution of other viral proteins in driving virulence. This study assesses the potential of long-read nanopore sequencing for fast and precise sequencing of circulating equine influenza viruses. Therefore, two French Florida Clade 1 strains, including the one circulating in winter 2018-2019 exhibiting more pronounced pathogenicity than usual, as well as the two currently OIE-recommended vaccine strains, were sequenced. Our results demonstrated the reliability of this sequencing method in generating accurate sequences. Sequence analysis of HA revealed a subtle antigenic drift in the French EIV strains, with specific substitutions, such as T163I in A/equine/Paris/1/2018 and the N188T mutation in post-2015 strains; both substitutions were in antigenic site B. Antigenic site E exhibited modifications in post-2018 strains, with the N63D substitution. Segment 2 sequencing also revealed that the A/equine/Paris/1/2018 strain encodes a longer variant of the PB1-F2 protein when compared to other Florida clade 1 strains (90 amino acids long versus 81 amino acids long). Further biological and biochemistry assays demonstrated that this PB1-F2 variant has enhanced abilities to abolish the mitochondrial membrane potential ΔΨm and permeabilize synthetic membranes. Altogether, our results highlight the interest in rapidly characterizing the complete genome of circulating strains with next-generation sequencing technologies to adapt vaccines and identify specific virulence markers of EIV.

Effect of sainfoin (Onobrychis viciifolia) on cyathostomin eggs excretion, larval development, larval community structure and efficacy of ivermectin treatment in horses.

Alternative strategies to chemical anthelmintics are needed for the sustainable control of equine strongylids. Bioactive forages like sainfoin (Onobrychis viciifolia) could contribute to reducing drug use, with the first hints of in vitro activity against cyathostomin free-living stages observed in the past. We analysed the effect of a sainfoin-rich diet on cyathostomin population and the efficacy of oral ivermectin treatment. Two groups of 10 naturally infected horses were enrolled in a 78-day experimental trial. Following a 1-week adaptation period, they were either fed with dehydrated sainfoin pellets (70% of their diet dry matter) or with alfalfa pellets (control group) for 21-days. No difference was found between the average fecal egg counts (FECs) of the two groups, but a significantly lower increase in larval development rate was observed for the sainfoin group, at the end of the trial. Quantification of cyathostomin species abundances with an ITS-2-based metabarcoding approach revealed that the sainfoin diet did not affect the nemabiome structure compared to the control diet. Following oral ivermectin treatment of all horses on day 21, the drug concentration was lower in horses fed with sainfoin, and cyathostomin eggs reappeared earlier in that group. Our results demonstrated that short-term consumption of a sainfoin-rich diet does not decrease cyathostomin FEC but seems to slightly reduce larval development. Consumption of dehydrated sainfoin pellets also negatively affected ivermectin pharmacokinetics, underscoring the need to monitor horse feeding regimes when assessing ivermectin efficacy in the field.

A new method for long-read sequencing of animal mitochondrial genomes: application to the identification of equine mitochondrial DNA variants.

BACKGROUND: Mitochondrial DNA is remarkably polymorphic. This is why animal geneticists survey mitochondrial genomes variations for fundamental and applied purposes. We present here an approach to sequence whole mitochondrial genomes using nanopore long-read sequencing. Our method relies on the selective elimination of nuclear DNA using an exonuclease treatment and on the amplification of circular mitochondrial DNA using a multiple displacement amplification step. RESULTS: We optimized each preparative step to obtain a 100 million-fold enrichment of horse mitochondrial DNA relative to nuclear DNA. We sequenced these amplified mitochondrial DNA using nanopore sequencing technology and obtained mitochondrial DNA reads that represented up to half of the sequencing output. The sequence reads were 2.3 kb of mean length and provided an even coverage of the mitochondrial genome. Long-reads spanning half or more of the whole mtDNA provided a coverage that varied between 118X and 488X. We evaluated SNPs identified using these long-reads by Sanger sequencing as ground truth and found a precision of 100.0%; a recall of 93.1% and a F1-score of 0.964 using the Twilight horse mtDNA reference. The choice of the mtDNA reference impacted variant calling efficiency with F1-scores varying between 0.947 and 0.964. CONCLUSIONS: Our method to amplify mtDNA and to sequence it using the nanopore technology is usable for mitochondrial DNA variant analysis. With minor modifications, this approach could easily be applied to other large circular DNA molecules.

Multi-species annotation of transcriptome and chromatin structure in domesticated animals.

BACKGROUND: Comparative genomics studies are central in identifying the coding and non-coding elements associated with complex traits, and the functional annotation of genomes is a critical step to decipher the genotype-to-phenotype relationships in livestock animals. As part of the Functional Annotation of Animal Genomes (FAANG) action, the FR-AgENCODE project aimed to create reference functional maps of domesticated animals by profiling the landscape of transcription (RNA-seq), chromatin accessibility (ATAC-seq) and conformation (Hi-C) in species representing ruminants (cattle, goat), monogastrics (pig) and birds (chicken), using three target samples related to metabolism (liver) and immunity (CD4+ and CD8+ T cells). RESULTS: RNA-seq assays considerably extended the available catalog of annotated transcripts and identified differentially expressed genes with unknown function, including new syntenic lncRNAs. ATAC-seq highlighted an enrichment for transcription factor binding sites in differentially accessible regions of the chromatin. Comparative analyses revealed a core set of conserved regulatory regions across species. Topologically associating domains (TADs) and epigenetic A/B compartments annotated from Hi-C data were consistent with RNA-seq and ATAC-seq data. Multi-species comparisons showed that conserved TAD boundaries had stronger insulation properties than species-specific ones and that the genomic distribution of orthologous genes in A/B compartments was significantly conserved across species. CONCLUSIONS: We report the first multi-species and multi-assay genome annotation results obtained by a FAANG project. Beyond the generation of reference annotations and the confirmation of previous findings on model animals, the integrative analysis of data from multiple assays and species sheds a new light on the multi-scale selective pressure shaping genome organization from birds to mammals. Overall, these results emphasize the value of FAANG for research on domesticated animals and reinforces the importance of future meta-analyses of the reference datasets being generated by this community on different species.

A 2.5-kilobase deletion containing a cluster of nine microRNAs in the latency-associated-transcript locus of the pseudorabies virus affects the host response of porcine trigeminal ganglia during established latency.

UNLABELLED: The alphaherpesvirus pseudorabies virus (PrV) establishes latency primarily in neurons of trigeminal ganglia when only the transcription of the latency-associated transcript (LAT) locus is detected. Eleven microRNAs (miRNAs) cluster within the LAT, suggesting a role in establishment and/or maintenance of latency. We generated a mutant (M) PrV deleted of nine miRNA genes which displayed properties that were almost identical to those of the parental PrV wild type (WT) during propagation in vitro. Fifteen pigs were experimentally infected with either WT or M virus or were mock infected. Similar levels of virus excretion and host antibody response were observed in all infected animals. At 62 days postinfection, trigeminal ganglia were excised and profiled by deep sequencing and quantitative RT-PCR. Latency was established in all infected animals without evidence of viral reactivation, demonstrating that miRNAs are not essential for this process. Lower levels of the large latency transcript (LLT) were found in ganglia infected by M PrV than in those infected by WT PrV. All PrV miRNAs were expressed, with highest expression observed for prv-miR-LLT1, prv-miR-LLT2 (in WT ganglia), and prv-miR-LLT10 (in both WT and M ganglia). No evidence of differentially expressed porcine miRNAs was found. Fifty-four porcine genes were differentially expressed between WT, M, and control ganglia. Both viruses triggered a strong host immune response, but in M ganglia gene upregulation was prevalent. Pathway analyses indicated that several biofunctions, including those related to cell-mediated immune response and the migration of dendritic cells, were impaired in M ganglia. These findings are consistent with a function of the LAT locus in the modulation of host response for maintaining a latent state. IMPORTANCE: This study provides a thorough reference on the establishment of latency by PrV in its natural host, the pig. Our results corroborate the evidence obtained from the study of several LAT mutants of other alphaherpesviruses encoding miRNAs from their LAT regions. Neither PrV miRNA expression nor high LLT expression levels are essential to achieve latency in trigeminal ganglia. Once latency is established by PrV, the only remarkable differences are found in the pattern of host response. This indicates that, as in herpes simplex virus, LAT functions as an immune evasion locus.

Validation of novel reference genes for RT-qPCR studies of gene expression in Xenopus tropicalis during embryonic and post-embryonic development.

BACKGROUND: Accurate interpretation of transcriptome profiling by quantitative PCR requires the establishment of species-specific standards. However, the selection of reference genes for assessing RNA expression profiles in Xenopus laevis and Xenopus tropicalis was mostly based on historical reasons and they often only reflect the traditions of a laboratory. RESULTS: We investigated the expression stability of 10 genes (dicer1, drosha, eef1a1, elavl3, gsc, h4, odc1, rpl8, smn2, tbp), 8 of which are commonly used as internal controls in published RT-qPCR experiments. We defined specific primer pairs and evaluated their suitability as reference genes by performing RT-qPCR expression profiling in Xenopus tropicalis. Gene expression stability was assayed in a set of 15 developmental stages from the egg to the froglet, and in dissected embryos. CONCLUSIONS: Overall, we determined a set of qualified reference genes for distinct developmental periods. We recommend the use of dicer1, drosha, eef1a1, and smn2 from early embryonic development up to the end of metamorphosis. During early embryogenesis drosha, eef1a1, smn2 are suitable. For the whole post-embryonic development and for metamorphic stages including pro-metamorphosis and metamorphic climax, we recommend the use of drosha and smn2. These reference genes should prove their usefulness for data comparison across studies.

In vitro maturation of oocytes alters gene expression and signaling pathways in bovine cumulus cells.

In vitro maturation (IVM) of immature oocytes is widely used in assisted reproduction technologies in cattle, and is increasingly used to treat human infertility. The development competence of IVM oocytes, however, is lower than preovulatory, in vivo-matured oocytes. During maturation, cumulus cells (CC) are metabolically coupled with an oocyte and support the acquisition of its developmental potential. Our objective was to identify genes and pathways that were affected by IVM in bovine CC. Microarray transcriptomic analysis of CC enclosing in vitro- or in vivo-mature oocytes revealed 472 differentially expressed genes, including 28% related to apoptosis, correlating with twofold higher cell death after IVM than in vivo, as detected by TUNEL. Genes overexpressed after IVM were significantly enriched in functions involved in cell movement, focal adhesion, extracellular matrix function, and TGF-beta signaling, whereas under-expressed genes were enriched in regulating gene expression, energy metabolism, stress response, and MAP kinases pathway functions. Differential expression of 15 genes, including PAG11 (increased) and TXNIP (decreased), which were never detected in CC before, was validated by real-time RT-PCR. Moreover, protein quantification confirmed the lower abundance of glutathione S-transferase A1 and prostaglandin G/H synthase 2, and the higher abundance of hyaluronan synthase 2 and SMAD4, a member of TGF-beta pathway, in CC after IVM. Phosphorylation levels of SMAD2, MAPK3/1, and MAPK14, but not MAPK8, were higher after IVM that in vivo. In conclusion, IVM provokes the hyper-activation of TGF-beta and MAPK signaling components, modifies gene expression, leads to increased apoptosis in CC, and thus affects oocyte quality.

Species interactions, stability, and resilience of the gut microbiota - Helminth assemblage in horses.

The nature and strength of interactions entertained among helminths and their host gut microbiota remain largely unexplored. Using 40 naturally infected Welsh ponies, we tracked the gut microbiota-cyathostomin temporal dynamics and stability before and following anthelmintic treatment and the associated host blood transcriptomic response. High shedders harbored 14 species of cyathostomins, dominated by Cylicocyclus nassatus. They exhibited a highly diverse and temporal dynamic gut microbiota, with butyrate-producing Clostridia likely driving the ecosystem steadiness and host tolerance toward cyathostomins infection. However, anthelmintic administration sharply bent the microbial community. It disrupted the ecosystem stability and the time-dependent network of interactions, affecting longer term microbial resilience. These observations highlight how anthelmintic treatments alter the triangular relationship of parasite, host, and gut microbiota and open new perspectives for adding nutritional intervention to current parasite management strategies.

Comparison of two molecular barcodes for the study of equine strongylid communities with amplicon sequencing.

Basic knowledge on the biology and epidemiology of equine strongylid species still needs to be improved to contribute to the design of better parasite control strategies. Nemabiome metabarcoding is a convenient tool to quantify and identify species in bulk samples that could overcome the hurdle that cyathostomin morphological identification represents. To date, this approach has relied on the internal transcribed spacer 2 (ITS-2) of the ribosomal RNA gene, with a limited investigation of its predictive performance for cyathostomin communities. Using DNA pools of single cyathostomin worms, this study aimed to provide the first elements to compare performances of the ITS-2 and a cytochrome c oxidase subunit I (COI) barcode newly developed in this study. Barcode predictive abilities were compared across various mock community compositions of two, five and 11 individuals from distinct species. The amplification bias of each barcode was estimated. Results were also compared between various types of biological samples, i.e., eggs, infective larvae or adults. Bioinformatic parameters were chosen to yield the closest representation of the cyathostomin community for each barcode, underscoring the need for communities of known composition for metabarcoding purposes. Overall, the proposed COI barcode was suboptimal relative to the ITS-2 rDNA region, because of PCR amplification biases, reduced sensitivity and higher divergence from the expected community composition. Metabarcoding yielded consistent community composition across the three sample types. However, imperfect correlations were found between relative abundances from infective larvae and other life-stages for Cylicostephanus species using the ITS-2 barcode. While the results remain limited by the considered biological material, they suggest that additional improvements are needed for both the ITS-2 and COI barcodes.

Characterization of a novel Xenopus tropicalis cell line as a model for in vitro studies.

Cell lines are useful tools to facilitate in vitro studies of many biological and molecular processes. We describe a new permanent fibroblast-type cell line obtained from disaggregated Xenopus tropicalis limb bud. The cell line population doubling time was ~24 h. Its karyotype was genetically stable with a chromosome number of 2n = 21 and a chromosome 10 trisomy. These cells could be readily transfected and expressed transgenes faithfully. We obtained stable transformants using transposon-based gene transfer technology. These cells responded to thyroid hormone and thus can provide a complementary research tool to study thyroid hormone signaling events. In conclusion, this cell line baptized "Speedy" should prove useful to couple in vitro and in vivo biological studies in the X. tropicalis frog model.

Influenza A virus protein PB1-F2 exacerbates IFN-beta expression of human respiratory epithelial cells.

The PB1-F2 protein of the influenza A virus (IAV) contributes to viral pathogenesis by a mechanism that is not well understood. PB1-F2 was shown to modulate apoptosis and to be targeted by the CD8(+) T cell response. In this study, we examined the downstream effects of PB1-F2 protein during IAV infection by measuring expression of the cellular genes in response to infection with wild-type WSN/33 and PB1-F2 knockout viruses in human lung epithelial cells. Wild-type virus infection resulted in a significant induction of genes involved in innate immunity. Knocking out the PB1-F2 gene strongly decreased the magnitude of expression of cellular genes implicated in antiviral response and MHC class I Ag presentation, suggesting that PB1-F2 exacerbates innate immune response. Biological network analysis revealed the IFN pathway as a link between PB1-F2 and deregulated genes. Using quantitative RT-PCR and IFN-ß gene reporter assay, we determined that PB1-F2 mediates an upregulation of IFN-ß expression that is dependent on NF-κB but not on AP-1 and IFN regulatory factor-3 transcription factors. Recombinant viruses knocked out for the PB1-F2 and/or the nonstructural viral protein 1 (the viral antagonist of the IFN response) genes provide further evidence that PB1-F2 increases IFN-ß expression and that nonstructural viral protein 1 strongly antagonizes the effect of PB1-F2 on the innate response. Finally, we compared the effect of PB1-F2 variants taken from several IAV strains on IFN-ß expression and found that PB1-F2-mediated IFN-ß induction is significantly influenced by its amino acid sequence, demonstrating its importance in the host cell response triggered by IAV infection.

Antisense oligonucleotides targeting ORF1b block replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).

The ongoing COVID-19 pandemic continues to pose a need for new and efficient therapeutic strategies. We explored antisense therapy using oligonucleotides targeting the severe acute respiratory syndrome coronavirus (SARS-CoV-2) genome. We predicted in silico four antisense oligonucleotides (ASO gapmers with 100% PTO linkages and LNA modifications at their 5' and 3'ends) targeting viral regions ORF1a, ORF1b, N and the 5'UTR of the SARS-CoV-2 genome. Efficiency of ASOs was tested by transfection in human ACE2-expressing HEK-293T cells and monkey VeroE6/TMPRSS2 cells infected with SARS-CoV-2. The ORF1b-targeting ASO was the most efficient, with a 71% reduction in the number of viral genome copies. N- and 5'UTR-targeting ASOs also significantly reduced viral replication by 55 and 63%, respectively, compared to non-related control ASO (ASO-C). Viral titration revealed a significant decrease in SARS-CoV-2 multiplication both in culture media and in cells. These results show that anti-ORF1b ASO can specifically reduce SARS-CoV-2 genome replication in vitro in two different cell infection models. The present study presents proof-of concept of antisense oligonucleotide technology as a promising therapeutic strategy for COVID-19.

A New Method for Sequencing the Mitochondrial Genome by Using Long Read Technology.

We describe a protocol to prepare a multiplexed mtDNA library from a blood sample for performing a long read sequencing of the mitochondrial genome. All steps are carefully described to get a high enrichment of mtDNA relative to total DNA extracted from the blood sample. The obtained mutiplexed library allows the production of long sequence mtDNA reads up to 16.5 kbp with a quality enabling variant-calling by using a portable sequencer (MinION, Oxford Nanopore Technologies).

Gut microbiota resilience in horse athletes following holidays out to pasture.

Elite horse athletes that live in individual boxes and train and compete for hours experience long-term physical and mental stress that compromises animal welfare and alters the gut microbiota. We therefore assessed if a temporary period out to pasture with conspecifics could improve animal welfare and in turn, favorably affect intestinal microbiota composition. A total of 27 athletes were monitored before and after a period of 1.5 months out to pasture, and their fecal microbiota and behavior profiles were compared to those of 18 horses kept in individual boxes. The overall diversity and microbiota composition of pasture and control individuals were temporally similar, suggesting resilience to environmental challenges. However, pasture exposure induced an increase in Ruminococcus and Coprococcus that lasted 1-month after the return to individual boxes, which may have promoted beneficial effects on health and welfare. Associations between the gut microbiota composition and behavior indicating poor welfare were established. Furthermore, withdrawn behavior was associated with the relative abundances of Lachnospiraceae AC2044 group and Clostridiales family XIII. Both accommodate a large part of butyrate-producing bacterial genera. While we cannot infer causality within this study, arguably, these findings suggest that management practices maintained over a longer period of time may moderate the behavior link to the gut ecosystem beyond its resilience potential.

Tissue Resources for the Functional Annotation of Animal Genomes.

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

Priming for welfare: gut microbiota is associated with equitation conditions and behavior in horse athletes.

We simultaneously measured the fecal microbiota and multiple environmental and host-related variables in a cohort of 185 healthy horses reared in similar conditions during a period of eight months. The pattern of rare bacteria varied from host to host and was largely different between two time points. Among a suite of variables examined, equitation factors were highly associated with the gut microbiota variability, evoking a relationship between gut microbiota and high levels of physical and mental stressors. Behavioral indicators that pointed toward a compromised welfare state (e.g. stereotypies, hypervigilance and aggressiveness) were also associated with the gut microbiota, reinforcing the notion for the existence of the microbiota-gut-brain axis. These observations were consistent with the microbiability of behaviour traits (> 15%), illustrating the importance of gut microbial composition to animal behaviour. As more elite athletes suffer from stress, targeting the microbiota offers a new opportunity to investigate the bidirectional interactions within the brain gut microbiota axis.

The miRNA-targeted transcriptome of porcine alveolar macrophages upon infection with Porcine Reproductive and Respiratory Syndrome Virus.

Host miRNAs are known to modulate the cell response to virus infections. We characterized the miRNA-targeted transcriptome of porcine alveolar macrophages (PAMs) at early times after infection with a subtype 1.1 strain of PRRSV (Porcine Reproductive and Respiratory Syndrome Virus). We performed the immunoprecipitation of RISC (RNA-induced Silencing Complex) followed by microarray analysis of the RISC-bound miRNA targets (RIP-Chip) to evaluate the relative enrichment or depletion of expressed genes in RISC. The miRNA-mediated regulation occurred early after PRRSV infection and decreased fast (1,241 and 141 RISC-bound genes at 7 h and 10 h post-infection, respectively); it affected several cell functions with evidence of miRNA buffering of upregulated interferon-related genes. Eight miRNAs were highly enriched in RISC of both control and infected cells with no evidence of differential expression. Although miR-335-5p was the miRNA with most predicted targets among enriched RISC-bound genes, no effects on surface markers, cytokine expression and PRRSV replication were detected upon miR-335-5p mimics of primary PAMs. Our results do not point to specific miRNA-driven mechanisms regulating the early response to infection with this PRRSV 1.1 strain and indicate that the miRNome expressed by steady-state PAMs reacts promptly to counterbalance PRRSV infection by a pervasive modulation of host functions.

Impact of the genetic background on the composition of the chicken plasma MiRNome in response to a stress.

Circulating extra-cellular microRNAs (miRNAs) have emerged as promising minimally invasive markers in human medicine. We evaluated miRNAs isolated from total plasma as biomarker candidates of a response to an abiotic stress (feed deprivation) in a livestock species. Two chicken lines selected for high (R+) and low (R-) residual feed intake were chosen as an experimental model because of their extreme divergence in feed intake and energy metabolism. Adult R+ and R- cocks were sampled after 16 hours of feed deprivation and again four hours after re-feeding. More than 292 million sequence reads were generated by small RNA-seq of total plasma RNA. A total of 649 mature miRNAs were identified; after quality filtering, 148 miRNAs were retained for further analyses. We identified 23 and 19 differentially abundant miRNAs between feeding conditions and between lines respectively, with only two miRNAs identified in both comparisons. We validated a panel of six differentially abundant miRNAs by RT-qPCR on a larger number of plasma samples and checked their response to feed deprivation in liver. Finally, we evaluated the conservation and tissue distribution of differentially abundant miRNAs in plasma across a variety of red jungle fowl tissues. We show that the chicken plasma miRNome reacts promptly to the alteration of the animal physiological condition driven by a feed deprivation stress. The plasma content of stress-responsive miRNAs is strongly influenced by the genetic background, with differences reflecting the phenotypic divergence acquired through long-term selection, as evidenced by the profiles of conserved miRNAs with a regulatory role in energy metabolism (gga-miR-204, gga-miR-let-7f-5p and gga-miR-122-5p). These results reinforce the emerging view in human medicine that even small genetic differences can have a considerable impact on the resolution of biomarker studies, and provide support for the emerging interest in miRNAs as potential novel and minimally invasive biomarkers for livestock species.

Early bovine embryos regulate oviduct epithelial cell gene expression during in vitro co-culture.

In mammals, the oviduct may participate to the regulation of early embryo development. In vitro co-culture of early bovine embryos with bovine oviduct epithelial cells (BOEC) has been largely used to mimic the maternal environment. However, the mechanisms of BOEC action have not been clearly elucidated yet. The aim of this study was to determine the response of BOEC cultures to the presence of developing bovine embryos. A 21,581-element bovine oligonucleotide array was used compare the gene expression profiles of confluent BOEC cultured for 8 days with or without embryos. This study revealed 34 differentially expressed genes (DEG). Of these 34 genes, IFI6, ISG15, MX1, IFI27, IFI44, RSAD2, IFITM1, EPSTI1, USP18, IFIT5, and STAT1 expression increased to the greatest extent due to the presence of embryos with a major impact on antiviral and immune response. Among the mRNAs at least 25 are already described as induced by interferons. In addition, transcript levels of new candidate genes involved in the regulation of transcription, modulation of the maternal immune system and endometrial remodeling were found to be increased. We selected 7 genes and confirmed their differential expression by quantitative RT-PCR. The immunofluorescence imaging of cellular localization of STAT1 protein in BOEC showed a nuclear translocation in the presence of embryos, suggesting the activation of interferon signaling pathway. This first systematic study of BOEC transcriptome changes in response to the presence of embryos in cattle provides some evidences that these cells are able to adapt their transcriptomic profile in response to embryo signaling.