Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development.

BACKGROUND: Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes. RESULTS: We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7). CONCLUSIONS: Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Interleukin-6 supplementation improves bovine conceptus elongation and transcriptomic indicators of developmental competence†.

A high incidence of pregnancy failures occurs in cattle during the second week of pregnancy as blastocysts transition into an elongated conceptus. This work explored whether interleukin-6 supplementation during in vitro embryo production would improve subsequent conceptus development. Bovine embryos were treated with 0 or 100 ng/mL recombinant bovine interleukin-6 beginning on day 5 post-fertilization. At day 7.5 post-fertilization, blastocysts were transferred into estrus synchronized beef cows (n = 5 recipients/treatment, 10 embryos/recipient). Seven days after transfer (day 14.5), cows were euthanized to harvest reproductive tracts and collect conceptuses. Individual conceptus lengths and stages were recorded before processing for RNA sequencing. Increases in conceptus recovery, length, and the proportion of tubular and filamentous conceptuses were detected in conceptuses derived from interleukin-6-treated embryos. The interleukin-6 treatment generated 591 differentially expressed genes in conceptuses (n = 9-10/treatment). Gene ontology enrichment analyses revealed changes in transcriptional regulation, DNA-binding, and antiviral actions. Only a few differentially expressed genes were associated with extraembryonic development, but several differentially expressed genes were associated with embryonic regulation of transcription, mesoderm and ectoderm development, organogenesis, limb formation, and somatogenesis. To conclude, this work provides evidence that interleukin-6 treatment before embryo transfer promotes pre-implantation conceptus development and gene expression in ways that resemble the generation of a robust conceptus containing favorable abilities to survive this critical period of pregnancy.

Tight gene co-expression in BCB positive cattle oocytes and their surrounding cumulus cells.

BACKGROUND: Cytoplasmic and nuclear maturation of oocytes, as well as interaction with the surrounding cumulus cells, are important features relevant to the acquisition of developmental competence. METHODS: Here, we utilized Brilliant cresyl blue (BCB) to distinguish cattle oocytes with low activity of the enzyme Glucose-6-Phosphate Dehydrogenase, and thus separated fully grown (BCB positive) oocytes from those in the growing phase (BCB negative). We then analyzed the developmental potential of these oocytes, mitochondrial DNA (mtDNA) copy number in single oocytes, and investigated the transcriptome of single oocytes and their surrounding cumulus cells of BCB positive versus BCB negative oocytes. RESULTS: The BCB positive oocytes were twice as likely to produce a blastocyst in vitro compared to BCB- oocytes (P < 0.01). We determined that BCB negative oocytes have 1.3-fold more mtDNA copies than BCB positive oocytes (P = 0.004). There was no differential transcript abundance of genes expressed in oocytes, however, 172 genes were identified in cumulus cells with differential transcript abundance (FDR < 0.05) based on the BCB staining of their oocyte. Co-expression analysis between oocytes and their surrounding cumulus cells revealed a subset of genes whose co-expression in BCB positive oocytes (n = 75) and their surrounding cumulus cells (n = 108) compose a unique profile of the cumulus-oocyte complex. CONCLUSIONS: If oocytes transition from BCB negative to BCB positive, there is a greater likelihood of producing a blastocyst, and a reduction of mtDNA copies, but there is no systematic variation of transcript abundance. Cumulus cells present changes in transcript abundance, which reflects in a dynamic co-expression between the oocyte and cumulus cells.

Fine-tuned adaptation of embryo-endometrium pairs at implantation revealed by transcriptome analyses in Bos taurus.

Interactions between embryo and endometrium at implantation are critical for the progression of pregnancy. These reciprocal actions involve exchange of paracrine signals that govern implantation and placentation. However, it remains unknown how these interactions between the conceptus and the endometrium are coordinated at the level of an individual pregnancy. Under the hypothesis that gene expression in endometrium is dependent on gene expression of extraembryonic tissues and genes expressed in extraembryonic tissues are dependent of genes expressed in the endometrium, we performed an integrative analysis of transcriptome profiles of paired extraembryonic tissue and endometria obtained from cattle (Bos taurus) pregnancies initiated by artificial insemination. We quantified strong dependence (|r| > 0.95, empirical false discovery rate [eFDR] < 0.01) in transcript abundance of genes expressed in the extraembryonic tissues and genes expressed in the endometrium. The profiles of connectivity revealed distinct coexpression patterns of extraembryonic tissues with caruncular and intercaruncular areas of the endometrium. Notably, a subset of highly coexpressed genes between extraembryonic tissue (n = 229) and caruncular areas of the endometrium (n = 218, r > 0.9999, eFDR < 0.001) revealed a blueprint of gene expression specific to each pregnancy. Gene ontology analyses of genes coexpressed between extraembryonic tissue and endometrium revealed significantly enriched modules with critical contribution for implantation and placentation, including "in utero embryonic development," "placenta development," and "regulation of transcription." Coexpressing modules were remarkably specific to caruncular or intercaruncular areas of the endometrium. The quantitative association between genes expressed in extraembryonic tissue and endometrium emphasize a coordinated communication between these two entities in mammals. We provide evidence that implantation in mammalian pregnancy relies on the ability of the extraembryonic tissue and the endometrium to develop a fine-tuned adaptive response characteristic of each pregnancy.

The blueprint of RNA storages relative to oocyte developmental competence in cattle (Bos taurus).

From the time oocytes leave quiescence, there are constant microenvironmental influences contributing to development, thus acquiring developmental competence is not a simple, linear phenomenon. During folliculogenesis, oocytes experience many morphological and cytological changes that contribute toward the acquisition of developmental competence, a process defined by an oocyte's ability to progress through folliculogenesis, be fertilized, undergo cleavage, and develop into an embryo. Many factors, such as ovarian follicle size, cow age, and the morphology of the cumulus-oocyte complex, have been extensively investigated to understand this process. In parallel to aiding in the understanding of oocyte biology, these features have been used to characterize an oocyte's ability to achieve competence. In addition, oocytes undergo intense gene transcription and protein translation to accumulate the maternal stores. When the oocyte is fully grown, most genes are transcriptionally inactive, and the chromatin is densely compacted. More recently, RNA profiling has been used to further define the transcriptional parameters that are associated with oocyte development. Here, focusing on cattle, we provide an overview of the experimental models commonly used to understand the underlying biology related to oocyte developmental competence. We compiled public data and showed that cattle oocytes can express over 15 000 protein-coding genes, suggesting a complex transcriptome landscape. Surprisingly, less than 2% of the expressed genes have been linked to developmental competence. The identification of the gene products that contribute to oocyte development, and understanding their biological function, are a vital component of our quest toward defining oocyte developmental competence at the molecular level.

Massive dysregulation of genes involved in cell signaling and placental development in cloned cattle conceptus and maternal endometrium.

A major unresolved issue in the cloning of mammals by somatic cell nuclear transfer (SCNT) is the mechanism by which the process fails after embryos are transferred to the uterus of recipients before or during the implantation window. We investigated this problem by using RNA sequencing (RNA-seq) to compare the transcriptomes in cattle conceptuses produced by SCNT and artificial insemination (AI) at day (d) 18 (preimplantation) and d 34 (postimplantation) of gestation. In addition, endometrium was profiled to identify the communication pathways that might be affected by the presence of a cloned conceptus, ultimately leading to mortality before or during the implantation window. At d 18, the effects on the transcriptome associated with SCNT were massive, involving more than 5,000 differentially expressed genes (DEGs). Among them are 121 genes that have embryonic lethal phenotypes in mice, cause defects in trophoblast and placental development, and/or affect conceptus survival in mice. In endometria at d 18, <0.4% of expressed genes were affected by the presence of a cloned conceptus, whereas at d 34, ∼36% and <0.7% of genes were differentially expressed in intercaruncular and caruncular tissues, respectively. Functional analysis of DEGs in placental and endometrial tissues suggests a major disruption of signaling between the cloned conceptus and the endometrium, particularly the intercaruncular tissue. Our results support a "bottleneck" model for cloned conceptus survival during the periimplantation period determined by gene expression levels in extraembryonic tissues and the endometrial response to altered signaling from clones.

Functional signaling and gene regulatory networks between the oocyte and the surrounding cumulus cells.

BACKGROUND: The maturation and successful acquisition of developmental competence by an oocyte, the female gamete, during folliculogenesis is highly dependent on molecular interactions with somatic cells. Most of the cellular interactions identified, thus far, are modulated by growth factors, ions or metabolites. We hypothesized that this interaction is also modulated at the transcriptional level, which leads to the formation of gene regulatory networks between the oocyte and cumulus cells. We tested this hypothesis by analyzing transcriptome data from single oocytes and the surrounding cumulus cells collected from antral follicles employing an analytical framework to determine interdependencies at the transcript level. RESULTS: We overlapped our transcriptome data with putative protein-protein interactions and identified hundreds of ligand-receptor pairs that can transduce paracrine signaling between an oocyte and cumulus cells. We determined that 499 ligand-encoding genes expressed in oocytes and cumulus cells are functionally associated with transcription regulation (FDR < 0.05). Ligand-encoding genes with specific expression in oocytes or cumulus cells were enriched for biological functions that are likely associated with the coordinated formation of transzonal projections from cumulus cells that reach the oocyte's membrane. Thousands of gene pairs exhibit significant linear co-expression (absolute correlation > 0.85, FDR < 1.8 × 10- 5) patterns between oocytes and cumulus cells. Hundreds of co-expressing genes showed clustering patterns associated with biological functions (FDR < 0.5) necessary for a coordinated function between the oocyte and cumulus cells during folliculogenesis (i.e. regulation of transcription, translation, apoptosis, cell differentiation and transport). CONCLUSION: Our analyses revealed a complex and functional gene regulatory circuit between the oocyte and surrounding cumulus cells. The regulatory profile of each cumulus-oocyte complex is likely associated with the oocytes' developmental potential to derive an embryo.

Transcriptome profiles in peripheral white blood cells at the time of artificial insemination discriminate beef heifers with different fertility potential.

BACKGROUND: Infertility is a longstanding limitation in livestock production with important economic impact for the cattle industry. Female reproductive traits are polygenic and lowly heritable in nature, thus selection for fertility is challenging. Beef cattle operations leverage estrous synchronization in combination with artificial insemination (AI) to breed heifers and benefit from an early and uniform calving season. A couple of weeks following AI, heifers are exposed to bulls for an opportunity to become pregnant by natural breeding (NB), but they may also not become pregnant during this time period. Focusing on beef heifers, in their first breeding season, we hypothesized that: a- at the time of AI, the transcriptome of peripheral white blood cells (PWBC) differs between heifers that become pregnant to AI and heifers that become pregnant late in the breeding season by NB or do not become pregnant during the breeding season; and b- the ratio of transcript abundance between genes in PWBC classifies heifers according to pregnancy by AI, NB, or failure to become pregnant. RESULTS: We generated RNA-sequencing data from 23 heifers from two locations (A: six AI-pregnant and five NB-pregnant; and B: six AI-pregnant and six non-pregnant). After filtering out lowly expressed genes, we quantified transcript abundance for 12,538 genes. The comparison of gene expression levels between AI-pregnant and NB-pregnant heifers yielded 18 differentially expressed genes (DEGs) (ADAM20, ALDH5A1, ANG, BOLA-DQB, DMBT1, FCER1A, GSTM3, KIR3DL1, LOC107131247, LOC618633, LYZ, MNS1, P2RY12, PPP1R1B, SIGLEC14, TPPP, TTLL1, UGT8, eFDR≤0.02). The comparison of gene expression levels between AI-pregnant and non-pregnant heifers yielded six DEGs (ALAS2, CNKSR3, LOC522763, SAXO2, TAC3, TFF2, eFDR≤0.05). We calculated the ratio of expression levels between all gene pairs and assessed their potential to classify samples according to experimental groups. Considering all samples, relative expression from two gene pairs correctly classified 10 out of 12 AI-pregnant heifers (P = 0.0028) separately from the other 11 heifers (NB-pregnant, or non-pregnant). CONCLUSION: The transcriptome profile in PWBC, at the time of AI, is associated with the fertility potential of beef heifers. Transcript levels of specific genes may be further explored as potential classifiers, and thus selection tools, of heifer fertility.

Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencing.

It remains an open question when and how the first cell fate decision is made in mammals. Using deep single-cell RNA-seq of matched sister blastomeres, we report highly reproducible inter-blastomere differences among 10 2-cell and five 4-cell mouse embryos. Inter-blastomere gene expression differences dominated between-embryo differences and noise, and were sufficient to cluster sister blastomeres into distinct groups. Dozens of protein-coding genes exhibited reproducible bimodal expression in sister blastomeres, which cannot be explained by random fluctuations. The protein expression of one gene out of four of these bimodal genes tested, Gadd45a, exhibited clear inter-blastomeric contrasts. We traced some of the bimodal mRNA expressions to embryonic genome activation, and others to blastomere-specific RNA depletion. Inter-blastomere differences created coexpression gene networks that were much stronger and larger than those that can possibly be created by random noise. The highly correlated gene pairs at the 4-cell stage overlapped with those showing the same directions of differential expression between inner cell mass (ICM) and trophectoderm (TE). These data substantiate the hypothesis of inter-blastomere differences in 2- and 4-cell mouse embryos, and associate these differences with ICM/TE differences.

Time-variant clustering model for understanding cell fate decisions.

Both spatial characteristics and temporal features are often the subjects of concern in physical, social, and biological studies. This work tackles the clustering problems for time course data in which the cluster number and clustering structure change with respect to time, dubbed time-variant clustering. We developed a hierarchical model that simultaneously clusters the objects at every time point and describes the relationships of the clusters between time points. The hidden layer of this model is a generalized form of branching processes. A reversible-jump Markov Chain Monte Carlo method was implemented for model inference, and a feature selection procedure was developed. We applied this method to explore an open question in preimplantation embryonic development. Our analyses using single-cell gene expression data suggested that the earliest cell fate decision could start at the 4-cell stage in mice, earlier than the commonly thought 8- to 16-cell stage. These results together with independent experimental data from single-cell RNA-seq provided support against a prevailing hypothesis in mammalian development.

In vitro maturation alters gene expression in bovine oocytes.

Gene expression profiling of in vivo- and in vitro-matured bovine oocytes can identify transcripts related to the developmental potential of oocytes. Nonetheless, the effects of in vitro culturing oocytes are yet to be fully understood. We tested the effects of in vitro maturation on the transcript profile of oocytes collected from Bos taurus indicus cows. We quantified the expression of 1488 genes in in vivo- and in vitro-matured oocytes. Of these, 51 genes were up-regulated, whereas 56 were down-regulated (≥2-fold) in in vivo-matured oocytes in comparison with in vitro-matured oocytes. Quantitative real-time polymerase chain reaction (PCR) of nine genes confirmed the microarray results of differential expression between in vivo- and in vitro-matured oocytes (EZR, EPN1, PSEN2, FST, IGFBP3, RBBP4, STAT3, FDPS and IRS1). We interrogated the results for enrichment of Gene Ontology categories and overlap with protein-protein interactions. The results revealed that the genes altered by in vitro maturation are mostly related to the regulation of oocyte metabolism. Additionally, analysis of protein-protein interactions uncovered two regulatory networks affected by the in vitro culture system. We propose that the differentially expressed genes are candidates for biomarkers of oocyte competence. In vitro oocyte maturation can affect the abundance of specific transcripts and are likely to deplete the developmental competence.

Messenger RNAs in metaphase II oocytes correlate with successful embryo development to the blastocyst stage.

The mRNAs accumulated in oocytes provide support for embryo development until embryo genomic activation. We hypothesized that the maternal mRNA stock present in bovine oocytes is associated with embryo development until the blastocyst stage. To test our hypothesis, we analyzed the transcriptome of the oocyte and correlated the results with the embryo development. Our goal was to identify genes expressed in the oocyte that correlate with its ability to develop to the blastocyst stage. A fraction of oocyte cytoplasm was biopsied using micro-aspiration and stored for further expression analysis. Oocytes were activated chemically, cultured individually and classified according to their capacity to develop in vitro to the blastocyst stage. Microarray analysis was performed on mRNA extracted from the oocyte cytoplasm fractions and correlated with its ability to develop to the blastocyst stage (good quality oocyte) or arrest at the 8-16-cell stage (bad quality oocyte). The expression of 4320 annotated genes was detected in the fractions of cytoplasm that had been collected from oocytes matured in vitro. Gene ontology classification revealed that enriched gene expression of genes was associated with certain biological processes: 'RNA processing', 'translation' and 'mRNA metabolic process'. Genes that are important to the molecular functions of 'RNA binding' and 'translation factor activity, RNA binding' were also enriched in oocytes. We identified 29 genes with differential expression between the two groups of oocytes compared (good versus bad quality). The content of mRNAs expressed in metaphase II oocytes influences the activation of the embryonic genome and enables further develop to the blastocyst stage.

Protocol for the electroporation of CRISPR-Cas for DNA and RNA targeting in Bos taurus zygotes.

The use of CRISPR-Cas9 ribonucleoproteins has revolutionized manipulation of genomes. Here, we present a protocol for the electroporation of CRISPR-Cas for DNA and RNA targeting in Bos taurus zygotes. First, we describe steps for production and preparation of presumptive zygotes for electroporation. The first electroporation introduces ribonucleoproteins formed by Cas9D10A with two guide RNAs to target DNA, and the second introduces the same ribonucleoprotein complex to target DNA plus Cas13a with one guide RNA to target RNAs. For complete details on the use and execution of this protocol, please refer to Nix et al.1.

Higher abundance of 2-dehydro-d-gluconate in the plasma of sub-fertile or infertile Bos taurus heifers.

Infertility or subfertility impacts approximately 5% and 15% of dairy and beef heifers (Bos taurus), respectively. Heifers that do not produce a calf within an optimum window of time have a significant negative impact on the profitability and sustainability of the cattle industry. Selection of heifers based on their fertility potential remains a challenge yet to be resolved. Here, we tested the hypothesis that heifers of different fertility potential have differing metabolome signatures in their plasma. We obtained blood from Bos taurus heifers at their first artificial insemination and processed the samples to separate the plasma. The heifers were classified based on their reproductive outcome as fertile (pregnant and delivered a calf after their first artificial insemination (AI)) or sub-fertile (Angus heifers: no pregnancy after two AI and exposure to a bull; Holstein heifers: no pregnancy by the third AI). We tested the relative abundance of 140 metabolites obtained from 22 heifers (Angus fertile n = 5, Angus sub-fertile n = 7, Holstein fertile N = 5, Holstein sub-fertile N = 5). The metabolite 2-Dehydro-D-gluconate (C6H10O7) was significantly more abundant in the plasma of sub-fertile heifers in both breeds (1.4-fold, false discovery rate < 0.1). In the context that a small proportion of circulating metabolites in the plasma were quantified in this study, the results show that the metabolomic profile in the blood stream may be associated with heifer fertility potential.

During the development of heifers for cow replacement, producers must invest substantial resources in each animal for over 15 mo. While the use of resources is equivalent across heifers being developed on a farm, a substantial proportion of the animals will not produce a calf (approximately 5% and approximately 15% of dairy and beef heifers, respectively). In this study, we identified one metabolite with higher abundance in the plasma of dairy and beef heifers with low chances of producing a calf by 25 mo of age.

Altered microRNA composition in the uterine lumen fluid in cattle (Bos taurus) pregnancies initiated by artificial insemination or transfer of an in vitro produced embryo.

BACKGROUND: MicroRNAs (miRNAs) are presented in the uterine lumen of many mammals, and in vitro experiments have determined that several miRNAs are important for the regulation of endometrial and trophoblast functions. Our aim was to identify and contrast the miRNAs present in extracellular vesicles (EVs) in the uterine lumen fluid (ULF) at the onset of attachment in cattle pregnancies (gestation d 18) initiated by artificial insemination (AI) or by the transfer of an in vitro-produced blastocyst (IVP-ET). A third group had no conceptus after the transfer of an IVP embryo. RESULTS: The abundance of 263 annotated miRNAs was quantified in the EVs collected from ULF. There was an increase in the transcript abundance of 20 miRNAs in the ULF EVs from the AI pregnant group, while 4 miRNAs had a lower abundance relative to the group not containing a conceptus. Additionally, 4 miRNAs were more abundant in ULF EVs in the AI pregnant group relative to IVP-ET group (bta-mir-17, bta-mir-7-3, MIR7-1, MIR18A). Specific miRNAs in the ULF EVs were co-expressed with messenger RNAs expressed in extra-embryonic tissues and endometrium, including genes that are known to be their targets. CONCLUSIONS: The results provide biological insights into the participation of miRNAs in the regulation of trophoblast proliferation and differentiation, as well as in endometrium receptivity. The knowledge that in vitro cultured embryos can contribute to the altered abundance of specific miRNAs in the uterine lumen can lead to the development of corrective approaches to reduce conceptus losses during the first month of pregnancy in cattle.

Non-Alcoholic Fatty Liver Disease Induced by Feeding Medium-Chain Fatty Acids Upregulates Cholesterol and Lipid Homeostatic Genes in Skeletal Muscle of Neonatal Pigs.

Non-alcoholic fatty liver disease (NAFLD) is a range of disorders characterized by lipid accumulation in hepatocytes. Although this spectrum of disorders is associated with adult obesity, recent evidence suggests that this condition could also occur independently of obesity, even in children. Previously, we reported that pigs fed a formula containing medium-chain fatty acids (MCFAs) developed hepatic steatosis and weighed less than those fed an isocaloric formula containing long-chain fatty acids (LCFAs). Our objective was to determine the association between NAFLD and the skeletal muscle transcriptome in response to energy and lipid intake. Neonatal pigs were fed one of three formulas: a control formula (CONT, n = 6) or one of two isocaloric high-energy formulas containing either long (LCFA, n = 6) or medium (MCFA, n = 6) chain fatty acids. Pigs were fed for 22 d, and tissues were collected. Body weight at 20 and 22 d was greater for LCFA-fed pigs than their CONT or MCFA counterparts (p < 0.005). Longissimus dorsi weight was greater for LCFA compared with MCFA, while CONT was intermediate (p < 0.05). Lean gain and protein deposition were greater for LCFA than for CONT and MCFA groups (p < 0.01). Transcriptomic analysis revealed 36 differentially expressed genes (DEGs) between MCFA and LCFA, 53 DEGs between MCFA and CONT, and 52 DEGs between LCFA and CONT (FDR < 0.2). Feeding formula high in MCFAs resulted in lower body and muscle weights. Transcriptomics data suggest that the reduction in growth was associated with a disruption in cholesterol metabolism in skeletal muscles.

The Ruminant Telomere-to-Telomere (RT2T) Consortium.

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

Changes in WNT signaling-related gene expression associated with development and cloning in bovine extra-embryonic and endometrial tissues during the peri-implantation period.

We determined if somatic cell nuclear transfer (SCNT) cloning is associated with WNT-related gene expression in cattle development, and if the expression of genes in the WNT pathway changes during the peri-implantation period. Extra-embryonic and endometrial tissues were collected at gestation days 18 and 34 (d18, d34). WNT5A, FZD4, FZD5, LRP5, CTNNB1, GNAI2, KDM1A, BCL2L1, and SFRP1 transcripts were localized in extra-embryonic tissue, whereas SFRP1 and DKK1 were localized in the endometrium. There were no differences in the localization of these transcripts in extra-embryonic tissue or endometrium from SCNT or artificial insemination (AI) pregnancies. Expression levels of WNT5A were 11-fold greater in the allantois of SCNT than AI samples. In the trophoblast, expression of WNT5A, FZD5, CTNNB1, and DKK1 increased significantly from d18 to d34, whereas expression of KDM1A and SFRP1 decreased, indicating that implantation is associated with major changes in WNT signaling. SCNT was associated with altered WNT5A expression in trophoblasts, with levels increasing 2.3-fold more in AI than SCNT conceptuses from d18 to d34. In the allantois, expression of WNT5A increased 6.3-fold more in SCNT than AI conceptuses from d18 to d34. Endometrial tissue expression levels of the genes tested did not differ between AI or SCNT pregnancies, although expression of individual genes showed variation across developmental stages. Our results demonstrate that SCNT is associated with altered expression of specific WNT-related genes in extra-embryonic tissue in a time- and tissue-specific manner. The pattern of gene expression in the WNT pathway suggests that noncanonical WNT signal transduction is important for implantation of cattle conceptuses.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing.

BACKGROUND: A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation. RESULTS: In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework. CONCLUSIONS: Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

A multi-omics analysis identifies molecular features associated with fertility in heifers (Bos taurus).

Infertility or subfertility is a critical barrier to sustainable cattle production, including in heifers. The development of heifers that do not produce a calf within an optimum window of time is a critical factor for the profitability and sustainability of the cattle industry. In parallel, heifers are an excellent biomedical model for understanding the underlying etiology of infertility because well-nourished heifers can still be infertile, mostly because of inherent physiological and genetic causes. Using a high-density single nucleotide polymorphism (SNP) chip, we collected genotypic data, which were analyzed using an association analysis in PLINK with Fisher's exact test. We also produced quantitative transcriptome data and proteome data. Transcriptome data were analyzed using the quasi-likelihood test followed by the Wald's test, and the likelihood test and proteome data were analyzed using a generalized mixed model and Student's t-test. We identified two SNPs significantly associated with heifer fertility (rs110918927, chr12: 85648422, P = 6.7 × 10-7; and rs109366560, chr11:37666527, P = 2.6 × 10-5). We identified two genes with differential transcript abundance (eFDR ≤ 0.002) between the two groups (Fertile and Sub-Fertile): Adipocyte Plasma Membrane Associated Protein (APMAP, 1.16 greater abundance in the Fertile group) and Dynein Axonemal Intermediate Chain 7 (DNAI7, 1.23 greater abundance in the Sub-Fertile group). Our analysis revealed that the protein Alpha-ketoglutarate-dependent dioxygenase FTO was more abundant in the plasma collected from Fertile heifers relative to their Sub-Fertile counterparts (FDR < 0.05). Lastly, an integrative analysis of the three datasets identified a series of molecular features (SNPs, gene transcripts, and proteins) that discriminated 21 out of 22 heifers correctly based on their fertility category. Our multi-omics analyses confirm the complex nature of female fertility. Very importantly, our results also highlight differences in the molecular profile of heifers associated with fertility that transcend the constraints of breed-specific genetic background.