RNA-Seq analysis uncovers transcriptomic variations between morphologically similar in vivo- and in vitro-derived bovine blastocysts
- Ashley M Driver†1,
- Francisco Peñagaricano†2,
- Wen Huang1,
- Khawaja R Ahmad3,
- Katie S Hackbart1,
- Milo C Wiltbank1 and
- Hasan Khatib2Email author
© Driver et al; licensee BioMed Central Ltd. 2012
Received: 15 November 2011
Accepted: 28 March 2012
Published: 28 March 2012
A valuable tool for both research and industry, in vitro fertilization (IVF) has applications range from gamete selection and preservation of traits to cloning. Although IVF has achieved worldwide use, with approximately 339,685 bovine embryos transferred in 2010 alone, there are still continuing difficulties with efficiency. It is rare to have more than 40% of fertilized in vitro cattle oocytes reach blastocyst stage by day 8 of culture, and pregnancy rates are reported as less than 45% for in vitro produced embryos. To investigate potential influences in-vitro fertilization (IVF) has on embryonic development, this study compares in vivo- and in vitro-derived bovine blastocysts at a similar stage and quality grade (expanded, excellent quality) to determine the degree of transcriptomic variation beyond morphology using RNA-Seq.
A total of 26,906,451 and 38,184,547 fragments were sequenced for in vitro and in vivo embryo pools, respectively. We detected expression for a total of 17,634 genes, with 793 genes showing differential expression between the two embryo populations with false discovery rate (FDR) < 0.05. There were also 395 novel transcribed units found, of which 45 were differentially expressed (FDR < 0.05). In addition, 4,800 genes showed evidence of alternative splicing, with 873 genes displaying differential alternative splicing between the two pools (FDR < 0.05). Using GO enrichment analysis, multiple biological pathways were found to be significantly enriched for differentially expressed genes (FDR < 0.01), including cholesterol and sterol synthesis, system development, and cell differentiation.
Thus, our results support that IVF may influence at the transcriptomic level and that morphology is limited in full characterization of bovine preimplantation embryos.
KeywordsIn vivo embryo In vitro fertilization RNA-Seq Transcriptome Alternative splicing
The use of in vitro fertilization (IVF) in cattle and in other species such as humans has seen rapid increases in the last few decades. In 2010 alone, over 339,685 in vitro-produced bovine embryos were transferred worldwide . In addition 60,190 infants were born in the United States in 2009 from assisted reproductive technologies (ART), 99% of which involved IVF . In addition, the International Committee for Monitoring Assisted Reproductive Technology's (ICMART) estimated that in 2003, approximately 232,000 babies had been born by ART worldwide .
Although the statistics reflect a striking amount of IVF use, they also depict relatively low success rates. It is rare to have more than 40% of fertilized in vitro cattle oocytes reach blastocyst stage, and of those transferred pregnancy rates are less than 45% for in vitro techniques . Similarly in humans, statistics from 2007 show an approximate 32-37% live birth rate from transfers . As such, there is increasing pressure to better understand the embryo's biological framework in order to interpret the seemingly low developmental rates. The reported sensitivity of the embryo to environmental factors, such as pH and ionic stress during the period in which they are manipulated for in-vitro procedures has raised further importance for more in-depth studies .
Morphologically, distinction between in vivo and IVF pre-implantation embryos have been confirmed by characteristics such as differing cytoplasmic color , lipid droplet composition , and mitochondrial content . A comparison of in vitro and in vivo blastocysts via light microscopy showed lower microvilli coverage and lipid content in the in vitro-produced embryos . In addition to morphology, previous studies have ascertained differences metabolically  and at the chromosomal level, with in vitro embryos incurring higher incidences of abnormalities such as mixiploidy .
Differential gene expression has gained increasing attention as a mechanism underlying the phenotypic abnormalities between in vivo and in vitro embryos [9, 12–15]. In addition to cattle, these differences have also been reported in sheep  and mice , suggesting that the IVF process itself may influence the biological framework of the early developing embryo. These studies, however, utilized expression microarrays, which provide limited transcriptomic data. In addition, in vivo embryos were derived from a limited number of super-ovulated donor cows, which presents a confounding variable of potential superovulation effects. In contrast, this study utilizes in vivo embryos from non-superovulated cows to eliminate potential sources of transcriptomic influence to characterize transcriptomic variations due to the in vitro process itself. We hypothesize that in vivo and in vitro embryo populations differ in their biological and transcriptomic characteristics and that the identification of these differences can assist in understanding the molecular mechanisms influencing development of IVF embryos. This study reports a fine scale assessment of the differential gene expression and alternative splicing between bovine in vivo- and in vitro-derived embryos to provide potential transcriptomic characteristics for further investigation.
Deep sequencing of the bovine blastocyst transcriptomes
Summary of sequence read alignments to the reference genome
In vitro embryos
In vivo embryos
Pair end reads
26,906,451 x 2
38,184,547 x 2
Total sequenced fragments
Total mapped fragments
Uniquely mapped fragments
Fragments mapped to autosomes and X chromosome
Fragments mapped to annotated genes
Fragments mapped to annotated exons
Fragments overlapped with annotated introns
Analysis of sequencing reads was done using Tophat software  for alignment with the reference genome (btau4.0). Of the total sequenced fragments, 83% and 85% were mapped to the reference genome for in vitro and in vivo embryos, respectively, and of these, 91% and 92%, respectively, were uniquely mapped to specific regions in the bovine genome (Table 1). Of those uniquely mapped, 19,556,926 and 28,326,794 were mapped to one of the autosomes or X chromosome for in vitro and in vivo embryos, respectively (Table 1).
Identification and analysis of novel transcribed units
Using stringent criteria, our analysis uncovered a total of 395 novel transcribed units. Of these, 45 show differential expression (FDR < 0.05), with 23 transcripts over-expressed in the in vivo embryos and 22 over-expressed in the in vitro embryos. Three transcripts had exclusive expression in the in vivo population while one transcript was exclusive to in vitro-produced embryos. When compared to known cattle ESTs from the UCSC genome browser , we found that 237 novel transcripts were supported by at least five cow ESTs with total coverage greater than 50%. Twenty of these transcripts showed significant differential expression.
Identification and validation of differentially expressed genes
Fold change differences in expression of genes differentially expressed between in vivo and in vitro embryos at FDR < 0.05
Genes over-expressed in in vivo embryos
Genes over-expressed in in vitro embryos
Differential alternative splicing between in vivo and in vitro embryos
Classification of genes with simple alternative splicing events or differential alternative splicing
Alternative splicing event
No. of events
No. of genes
No. of genes that showed differential splicing*
Alternative last exon
Gene ontology analysis
Fold difference and q-values of 11 cholesterol biosynthesis pathway genes
Depth of RNA-Seq analysis
Our current knowledge on dynamic changes in gene expression during embryonic development has for the most part been based on results from microarray technology. Although providing insight, microarrays are based on a priori knowledge and oligonucleotide design and may not detect lowly expressed transcripts or provide complete genome coverage due to lack of probe availability. Also, expression microarrays do not provide information regarding alternative splicing and novel transcripts. Therefore, we utilized RNA-Seq technology in this study to perform an in-depth analysis of the transcriptomic landscape of bovine in vivo- and IVF-derived embryos at the blastocyst stage. This is the first report of a high-resolution snapshot of transcriptomic differences between in vivo and in vitro embryos beyond expression for any species using RNA-Seq technology.
Another strength of RNA-Seq technology is the ability to improve genome annotation by the discovery of novel transcripts. In our study, we identified 307,646 bp covering novel exonic regions from 395 novel transcript units. Of these, 45 transcripts showed differential expression between in vivo and in vitro embryos (FDR < 0.05), with 23 expressed higher in the in vivo embryos and 22 in the in vitro embryos. In addition, three novel transcripts were exclusive to the in vivo embryos and one to the in vitro. Thus, this is evidence that not only are there uncharacterized regions of the bovine genome, but also that these novel regions may serve important biological roles.
Although the two embryo groups possessed similar morphological appearance and would have received a similar rating for developmental stage and embryo quality, they possessed numerous transcriptomic differences. Of the 793 genes showing differential expression, 413 had greater than a 10-fold difference in mRNA concentrations.
In addition, RNA-Seq technology provided valuable information regarding alternative and novel splice variants reflecting more complex mechanisms of RNA regulation. A recent study by Corcoran et al.  found that many genes involved in RNA processing and translation were differentially expressed between in vivo and in vitro blastocyst populations. Our study extends these observations by identifying a number of alternatively spliced genes between in vivo and in vitro populations.
Overall, our study detected 873 genes with alternative splicing, and of these, 782 showed differential splicing without showing overall differential gene expression. This is critical as alternative splicing can influence protein production and biological function. Associations between differential alternative splicing and embryonic development have been reported by Zhang et al. . This study reported differential expression levels of alternative transcripts for heat shock proteins in developmentally arrested in vitro bovine embryos in comparison to developed blastocysts, with the suggestion that differential processing of RNA is associated with regulation of embryonic development. Thus, our study is presenting transcriptomic information that would otherwise not be detected by standard microarray and yet have key biological importance for development.
Beyond individual gene differences, our study also revealed several biological pathways that were significantly enriched or depleted based on substantial differences in embryonic mRNA concentrations. The 23 pathways showing the highest level of significance can be found in Additional file 1.
Of particular interest was the cholesterol biosynthesis pathway, which had 11 genes with increased expression in the in vitro embryos. Differential expression levels of the cholesterol biosynthesis pathway between our embryo populations could be supported by a study conducted by Tint et al. establishing that the mother provides necessary cholesterol to the in vivo embryo throughout the pre-implantation period . Since in vitro embryos lack contact with maternal tissues, and thus maternal sources of cholesterol, there could be increased de novo production of cholesterol.
Of the differential genes listed in Table 4, many have been shown to be critical to the cholesterol biosynthesis process. For example, lanosterol synthase (LSS) catalyzes (S)-2,3 oxidosqualene into the parental compounds for cholesterol and ultimate steroid synthesis in mammals . In addition, the gene HMGCR contributes to the major rate-limiting step for conversion of HMG-CoA into mevalonate . HMGCR is also regulated through a feedback system, which could be altered by a lack of maternal influence and delivery on the embryo . Should the proper levels of these genes be perturbed, it could have significant effects on an organism by altering products of cholesterol (i.e. steroids). Interestingly it should be noted that another pathway related to steroid production (sterol biosynthesis) was also enriched between in vivo and in vitro embryos.
In addition, the cell differentiation pathway was enriched which is intriguing as differentiation is not only critical for the blastocyst formation but to ensure proper implantation and placenta development . This finding supports earlier histological studies showing fewer ICM cells in in vitro-produced bovine embryos compared to in vivo . Detection of this differential pathway in our study provides further evidence for its potential importance, thus helping to further the ongoing characterization of in vivo and in vitro embryos. In conclusion, our results suggest that IVF may influence at a larger scale by altering biological pathways critical for early development.
One limitation to this study is the potential shortening of reads due to the RNA amplification process, which tends to favor the 3' end of the sequence. In addition this may cause a lack of complete information on full length alternative splicing. Due to limited amounts of RNA in individual embryos, two rounds of linear amplification were necessary. Although it remains unclear as to the degree of transcript shortening that occurs, it should be acknowledged as it may prevent construction of full length transcripts as discussed by Huang and Khatib . Nonetheless, amplification provided a means to utilize individual embryos for validation and pool construction.
A second limitation to this study was the lack of RNA-Seq replicates. Nonetheless, in a prior RNA-Seq study by our lab, two pools were utilized and results were consistent with qRT-PCR validations and prior microarray results that used replicates . Thus, even with limited number of samples, RNA-Seq results from this study provide valuable insight into the genomic characteristics of embryo populations. As such, although there are acknowledged limitations, using RNA-Seq approach as a hypothesis-generating study provides numerous opportunities for further investigation to better understand early embryonic development.
Use of IVF and in vitro embryo culture has been found to be a valuable tool for both research and industry with applications ranging from infertility treatment, and gamete selection to cloning. In addition, the bovine embryo model has become an interesting and popular biological model to help in understanding potential problems or optimization procedures for human embryonic practices . As such, results of this study may help elucidate the molecular mechanisms leading to seemingly low success rates of IVF embryos and may help to unlock important factors underlying the early embryo's growth and sensitivity during the pre-implantation period. Overall, our results support the initial hypothesis that in vitro embryos differ from their in vivo counterparts on multiple transcriptomic levels. In addition, this study highlighted the importance of studying RNA processing beyond the levels of expression and the limitations morphology alone has on embryo characterization. Candidate genes and pathways uncovered in this study can be pursued in future studies with the goal of improvement of in vitro fertilization and culture systems and to further our understanding of factors involved in optimal embryonic development.
In vivo embryo collection and grading
In vitro maturation, fertilization, and embryo culture
Ovaries from Holstein cows were obtained from a local abattoir, and from them antral follicles (2-6 mm) were aspirated. Recovered oocytes were then washed with warm Tyrode's albumin lactate pyruvate (TALP)-Hepes and underwent maturation, fertilization, and culture as outlined in detail by Khatib et al.  and briefly here. For fertilization, frozen-thawed semen from one bull used in producing the in vivo embryos was prepared by centrifuging through a discontinuous Percoll gradient (45-90%) to recover motile sperm and adjusted to a final concentration of 1 × 106/ml . After 8 days of culture, embryos were morphologically assessed and those with an IETS grading of 7-1 (expanded blastocyst, excellent quality) were collected individually in RNALater (Ambion) and frozen at -20°C until RNA extraction (Figure 3).
RNA extraction, amplification, and pooling
Total RNA was extracted from individual embryos using RNaqueous Micro (Ambion). Due to limitations in the amount of RNA in each embryo, linear amplification was performed using the MessageAmp II aRNA amplification kit (Ambion). Briefly, cDNA was primed by a T7 promoter tagged poly-dT primer. After the second strand synthesis, the double stranded cDNA was used as template in an in vitro transcription reaction using T7 polymerase with unlabeled NTPs. Two rounds of linear amplification provided a proper concentration of RNA for analysis. In order to minimize biological variance between the in vivo and in vitro embryo groups for the transcriptomic comparisons, five grade 7-1 embryos from each group from the same sire were used. Samples of the amplified RNA from individual embryos were then pooled (n = 5 per pool) and prepped for RNA-sequencing.
Libraries of amplified RNA for each pool were prepared following the Illumina mRNA-Seq protocol. Sequencing libraries were created from 50 ng samples and sequenced with Illumina's HiSeq 2000 at the University of Wisconsin-Madison Biotechnology Center. Libraries were barcoded, multiplexed, and sequenced in one HiSeq 2000 lane. A 'fragment' was defined as a cDNA fragment sequenced from both ends. Approximately 60 million fragments were sequenced for both libraries and de-multiplexed.
Mapping reads to the reference genome
Sequencing reads were mapped to the reference genome (btau4.0) using the software package Tophat (v1.2.0)  implemented on Galaxy . The alignment was first performed in each of the samples independently. Novel splice junctions discovered by Tophat and known splice junctions from the Ensembl annotation were then combined and supplied to Tophat for a second alignment such that the junction database was the same for each sample. This mapping strategy allows a full utilization of the novel junctions identified in both samples. A maximum of two mismatches and a minimum length of 25 bp per segment were allowed. Moreover, reads that mapped equally well to more than 40 genomic locations were discarded.
Assembly of transcripts and estimation of abundance
Cufflinks (v1.0.3) was used to assemble transcript models from RNA-Seq alignments and to estimate transcript and gene expression . This program uses graph theory to find a parsimonious set of transcript models that comply with alignments. In addition, it estimates isoform and gene expression (fragments per kilobase exon per million mapped fragments, or FPKM) by optimizing a likelihood function containing transcript abundances as parameters. In our analysis, abundances of transcripts were upper-quartile normalized and also corrected for sequence bias in order to improve expression estimates .
Identification of novel transcribed units
Novel transcribed units were defined as regions containing multi-exonic transcripts that were: 1) at least 1,000 bp away from known gene boundaries; 2) of length ≥ 250 bp; 3) with an average coverage ≥ 5; 4) with less than 50% of repetitive sequences; and 5) with splice junctions supported by at least two junction alignments. For external support of our novel transcripts, cow EST alignments were downloaded from the UCSC genome browser and compared to the novel transcribed units using BEDTools .
Overall gene expression and GO enrichment analysis
Differential expression of annotated genes and novel transcribed units was tested using Cuffdiff, a program part of the Cufflinks package for testing differential gene expression . In our study, genes with a false discovery rate (FDR) < 0.05 were considered significant. FDR is defined as the expected proportion of false positives among all significant hypotheses . The Gene Ontology (GO) enrichment analysis was performed using the GOseq (v1.4.0) package  that is available in the R language/environment . Importantly, GOseq adjusts for gene selection bias due to difference in gene lengths, which is known to affect the variance of gene expression estimates . Biological pathways with a FDR < 0.01 were considered significant.
Analysis of alternative splicing
Differential alternative splicing of genes that showed more than one isoform was tested also with Cuffdiff using the Shannon-Jensen divergence metric to measure relative abundances between transcripts . A FDR < 0.05 was chosen as the significance threshold for detecting differential alternative splicing between the two treatment groups. Alternative splicing events were classified according to Wang et al.  using a custom perl script.
Validation of RNA-Seq data
Primers used for the validation of gene expression and alternative splicing
Primer sequence 5'→3'
Amplicon size (bp)
Primers for gene expression validation
Primers for alternative splicing validation
This study was supported by USDA Hatch grant WIS-142-PRJ17PH from the University of Wisconsin-Madison. KA was supported by the Higher Education Commission, Pakistan. The authors thank Ricky Monson and John Parrish for assistance with the IVF experiments. The authors also thank Jerry Guenther, Mary Herlihy, and Abdulkadir Keskin for the in vivo embryo collection.
- Stroud B: IETS 2011 Statistics and Data Retrieval Committee Report The year 2010 worldwide statistics of embryo transfer in domestic farm animals. Embryo Transfer Newsletter. 2011, 29: 1-23.
- CDC- Assisted Reproductive Technology: [http://www.cdc.gov/art/]
- Nygren K, Sullivan E, Zegers-Hochschild F, Mansour T, Ishihara O, Adamson D, de Mouzon J: International Committee for Monitoring Assisted Reproductive Technology (ICMART) world report: assisted reproductive technology 2003. Fertil Steril. 2011, 95: 2209-2222.e17. 10.1016/j.fertnstert.2011.03.058.View ArticlePubMed
- Farin PW, Piedrahita JA, Farin CE: Errors in development of fetuses and placentas from in vitro-produced bovine embryos. Theriogenology. 2006, 5: 178-191.View Article
- Gardner DK, Lane M: Ex vivo early embryo development and effects on gene expression and imprinting. Reprod Fertil Dev. 2005, 17: 361-370. 10.1071/RD04103.View ArticlePubMed
- Abe H, Otoi T, Tachikawa S, Yamashita S, Satoh T, Hoshi H: Fine structure of bovine morulae and blastocysts in vivo and in vitro. Anat Embryo. 1999, 199: 519-527. 10.1007/s004290050249.View Article
- Fair T, Lonergan P, Dinnyes A, Cottell DC, Hyttel P, Ward FA, Boland MP: Ultrastructure of bovine blastocysts following cryopreservation: effect of method of blastocyst production. Mol Reprod Dev. 2001, 58: 186-195. 10.1002/1098-2795(200102)58:2<186::AID-MRD8>3.0.CO;2-N.View ArticlePubMed
- Crosier AE, Farin PW, Dykstra MJ, Alexander JE, Farin CE: Ultrastructural morphometry of bovine blastocysts produced in vivo or in vitro. Biol Reprod. 2001, 64: 1375-1385. 10.1095/biolreprod64.5.1375.View ArticlePubMed
- Rizos D, Fair T, Papadopoulos S, Boland MP, Longergan P: Developmental, qualitative, and ultrastructural differences between ovine and bovine embryos produced in vivo or in vitro. Mol Reprod Dev. 2002, 62: 320-327. 10.1002/mrd.10138.View ArticlePubMed
- Khurana NK, Niemann H: Energy metabolism in preimplantation bovine embryos derived in vitro or in vivo. Biol Reprod. 2000, 62: 847-856. 10.1095/biolreprod62.4.847.View ArticlePubMed
- Hyttel P, Vluff D, Laurincik J, Schmidt M, Thomsen PD, Avery B, Callesen H, Rath D, Niemann H, Rosenkranz H, Shellander K, Ochs RL, Greve T: Risks of in-vitro production of cattle and swine embryos: aberrations in chromosome numbers, ribosomal RNA gene activation and perinatal physiology. Hum Reprod. 2000, 5 (Suppl 5): 87-97.View Article
- Wrenzycki C, Herrmann D, Keskintepe L, Martins A, Sirisathien S, Brackett B, Niemann H: Effects of culture system and protein supplementation on mRNA expression in pre-implantation bovine embryos. Hum Reprod. 2001, 16: 893-901. 10.1093/humrep/16.5.893.View ArticlePubMed
- Corcoran D, Fair T, Park S, Rizos D, Patel OV, Smith GW, Coussens PM, Ireland JJ, Boland MP, Evans AC, Lonergan P: Suppressed expression of genes involved in transcription and translation in in vitro compared with in vivo cultured bovine embryos. Reproduction. 2006, 131: 651-660. 10.1530/rep.1.01015.View ArticlePubMed
- Kues WA, Sudheer S, Herrmann D, Carnwath JW, Havlicek V, Besenfelder U, Lehrach H, Adjaye J, Niemann H: Genome-wide expression profiling reveals distinct clusters of transcriptional regulation during bovine preimplantation development in vivo. Proc Natl Acad Sci USA. 2008, 105: 19768-19773. 10.1073/pnas.0805616105.PubMed CentralView ArticlePubMed
- Smith SL, Everts RE, Sung LY, Du F, Page RL, Henderson B, Rodriguez-Zas SL, Nedambale TL, Renard JP, Lewin HA, Yang X, Tian XC: Gene expression of single bovine embryos uncovers significant effects of in vitro maturation, fertilization, and culture. Mol Reprod Dev. 2009, 76: 38-47. 10.1002/mrd.20927.View ArticlePubMed
- Powell K, Rooke JA, McEvoy TG, Ashworth CJ, Robinson JJ, Wilmut I, Young LE, Sinclair KD: Zygote donor nitrogen metabolism and in vitro embryo culture perturbs in utero development and IGF2R expression in ovine fetal tissues. Theriogenology. 2006, 66: 1901-1912. 10.1016/j.theriogenology.2006.05.008.View ArticlePubMed
- Giritharan G, Talbi S, Donjacour A, Di Sebastiano F, Dobson AT, Rinaudo PF: Effect of in vitro fertilization on gene expression and development of mouse preimplantation embryos. Reproduction. 2007, 134: 63-72. 10.1530/REP-06-0247.View ArticlePubMed
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMed
- UCSC Genome Browser. [http://genome.ucsc.edu/]
- Werner T: Next generation sequencing in functional genomics. Brief Bioinform. 2010, 11: 499-511. 10.1093/bib/bbq018.View ArticlePubMed
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mary C, Kinsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.PubMed CentralView ArticlePubMed
- Zhang B, Peñagaricano F, Driver A, Chen H, Khatib H: Differential expression of heat shock protein genes and their splice variants in bovine preimplantation embryos. J Dairy Sci. 2011, 94: 4174-4182. 10.3168/jds.2010-4137.View ArticlePubMed
- Tint GS, Yu H, Shang Q, Xu G, Patel SB: The use of Dhcr7 knockout mouse to accurately determine the origin of fetal sterols. J Lipid Res. 2006, 47: 1535-1541. 10.1194/jlr.M600141-JLR200.PubMed CentralView ArticlePubMed
- Villagra A, Ulloa N, Zhang X, Zhigang Y, Sotomayor E, Seto E: Histone Deacetylase 3 Down-regulates Cholesterol Synthesis through Repression of Lanosterol Synthase Gene Expression. J Bio Chem. 2007, 282: 35457-35470. 10.1074/jbc.M701719200.View Article
- Burg S, Espenshade P: Regulation of HMG-CoA reductase in mammals and yeast. Prog Lipid Res. 2011, 50: 403-410. 10.1016/j.plipres.2011.07.002.PubMed CentralView ArticlePubMed
- Senner CE, Hemberger M: Regulation of early trophoblast differentiation- lessons from the mouse. Placenta. 2010, 31: 944-950. 10.1016/j.placenta.2010.07.013.View ArticlePubMed
- Marquant-Le Guienne B, Gerard M, Solari A, Thibault C: In vitro culture of bovine egg fertilized either in vitro or in vivo. Reprod Nutr Dev. 1989, 29: 559-568. 10.1051/rnd:19890505.View ArticlePubMed
- Huang W, Khatib H: Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq. BMC Genomics. 2010, 11: 711-10.1186/1471-2164-11-711.PubMed CentralView ArticlePubMed
- Menezo Y, Herubel F: Mouse and bovine models for human IVF. Reprod BioMed Online. 2002, 4: 170-175. 10.1016/S1472-6483(10)61936-0.View ArticlePubMed
- Souza AH, Ayres H, Ferreira RM, Wiltbank MC: A new presynchronization system (Double-Ovsynch) increases fertility at first postpartum timed AI in lactating dairy cows. Theriogenology. 2008, 70: 208-215. 10.1016/j.theriogenology.2008.03.014.View ArticlePubMed
- Hackbart KS, Ferreira RM, Dietsche AA, Socha MT, Shaver RD, Wiltbank MC, Fricke PM: Effect of dietary organic zinc, manganese, copper, and cobalt supplementation on milk production, follicular growth, embryo quality, and tissue mineral concentrations in dairy cows. J Anim Sci. 2010, 88: 3856-3870. 10.2527/jas.2010-3055.View ArticlePubMed
- Stringfellow DA: Manual of the International Embryo Transfer Society: a procedural guide and general information for the use of embryo transfer technology, emphasizing sanitary precautions. The Society. 1998
- Khatib H, Monson RL, Schutzkus V, Kohl DM, Rosa GJM, Rutledge JJ: Mutations in the STAT5A Gene Are Associated with Embryonic Survival and Milk Composition in Cattle. J Dairy Sci. 2008, 91: 784-793. 10.3168/jds.2007-0669.View ArticlePubMed
- Parrish JJ, Krogenaes A, Susko-Parrish JL: Effect of bovine sperm separation by either swim-up or Percoll method on success of in vitro fertilization and early embryonic development. Theriogenology. 1995, 44: 859-869. 10.1016/0093-691X(95)00271-9.View ArticlePubMed
- Goecks J, Nekrutenko A, Taylor J: Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010, 11: R86-10.1186/gb-2010-11-8-r86.PubMed CentralView ArticlePubMed
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology. 2010, 28: 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMed
- Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L: Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011, 12: R22-10.1186/gb-2011-12-3-r22.PubMed CentralView ArticlePubMed
- Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.PubMed CentralView ArticlePubMed
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Statist Soc B. 1995, 57: 963-971.
- Young MD, Wakefield MJ, Smyth GK, Oshlack A: Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010, 11: R14-10.1186/gb-2010-11-2-r14.PubMed CentralView ArticlePubMed
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2009, Vienna, Austria: R Foundation for Statistical Computing
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: research0034-PubMed CentralView ArticlePubMed
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.