Genomic imbalance modulates transposable element expression in maize.

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

BICOSS: Bayesian iterative conditional stochastic search for GWAS.

BACKGROUND: Single marker analysis (SMA) with linear mixed models for genome wide association studies has uncovered the contribution of genetic variants to many observed phenotypes. However, SMA has weak false discovery control. In addition, when a few variants have large effect sizes, SMA has low statistical power to detect small and medium effect sizes, leading to low recall of true causal single nucleotide polymorphisms (SNPs). RESULTS: We present the Bayesian Iterative Conditional Stochastic Search (BICOSS) method that controls false discovery rate and increases recall of variants with small and medium effect sizes. BICOSS iterates between a screening step and a Bayesian model selection step. A simulation study shows that, when compared to SMA, BICOSS dramatically reduces false discovery rate and allows for smaller effect sizes to be discovered. Finally, two real world applications show the utility and flexibility of BICOSS. CONCLUSIONS: When compared to widely used SMA, BICOSS provides higher recall of true SNPs while dramatically reducing false discovery rate.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize.

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize.

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

Sequence of the supernumerary B chromosome of maize provides insight into its drive mechanism and evolution.

B chromosomes are enigmatic elements in thousands of plant and animal genomes that persist in populations despite being nonessential. They circumvent the laws of Mendelian inheritance but the molecular mechanisms underlying this behavior remain unknown. Here we present the sequence, annotation, and analysis of the maize B chromosome providing insight into its drive mechanism. The sequence assembly reveals detailed locations of the elements involved with the cis and trans functions of its drive mechanism, consisting of nondisjunction at the second pollen mitosis and preferential fertilization of the egg by the B-containing sperm. We identified 758 protein-coding genes in 125.9 Mb of B chromosome sequence, of which at least 88 are expressed. Our results demonstrate that transposable elements in the B chromosome are shared with the standard A chromosome set but multiple lines of evidence fail to detect a syntenic genic region in the A chromosomes, suggesting a distant origin. The current gene content is a result of continuous transfer from the A chromosomal complement over an extended evolutionary time with subsequent degradation but with selection for maintenance of this nonvital chromosome.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize.

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize.

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Magnitude of modulation of gene expression in aneuploid maize depends on the extent of genomic imbalance.

Aneuploidy has profound effects on an organism, typically more so than polyploidy, and the basis of this contrast is not fully understood. A dosage series of the maize long arm of chromosome 1 (1L) was used to compare relative global gene expression in different types and degrees of aneuploidy to gain insights into how the magnitude of genomic imbalance as well as hypoploidy affects global gene expression. While previously available methods require a selective examination of specific genes, RNA sequencing provides a whole-genome view of gene expression in aneuploids. Most studies of global aneuploidy effects have concentrated on individual types of aneuploids because multiple dose aneuploidies of the same genomic region are difficult to produce in most model genetic organisms. The genetic toolkit of maize allows the examination of multiple ploidies and 1-4 doses of chromosome arms. Thus, a detailed examination of expression changes both on the varied chromosome arms and elsewhere in the genome is possible, in both hypoploids and hyperploids, compared with euploid controls. Previous studies observed the inverse trans effect, in which genes not varied in DNA dosage were expressed in a negative relationship to the varied chromosomal region. This response was also the major type of changes found globally in this study. Many genes varied in dosage showed proportional expression changes, though some were seen to be partly or fully dosage compensated. It was also found that the effects of aneuploidy were progressive, with more severe aneuploids producing effects of greater magnitude.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants.

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

An empirical bayesian approach for testing gene expression fold change and its application in detecting global dosage effects.

We are motivated by biological studies intended to understand global gene expression fold change. Biologists have generally adopted a fixed cutoff to determine the significance of fold changes in gene expression studies (e.g. by using an observed fold change equal to two as a fixed threshold). Scientists can also use a t-test or a modified differential expression test to assess the significance of fold changes. However, these methods either fail to take advantage of the high dimensionality of gene expression data or fail to test fold change directly. Our research develops a new empirical Bayesian approach to substantially improve the power and accuracy of fold-change detection. Specifically, we more accurately estimate gene-wise error variation in the log of fold change. We then adopt a t-test with adjusted degrees of freedom for significance assessment. We apply our method to a dosage study in Arabidopsis and a Down syndrome study in humans to illustrate the utility of our approach. We also present a simulation study based on real datasets to demonstrate the accuracy of our method relative to error variance estimation and power in fold-change detection. Our developed R package with a detailed user manual is publicly available on GitHub at https://github.com/cuiyingbeicheng/Foldseq.

Conditions of embryo culture from days 5 to 7 of development alter the DNA methylome of the bovine fetus at day 86 of gestation.

PURPOSE: We tested whether in vitro production (IVP) causes changes in DNA methylation in fetal liver and skeletal muscle and if exposure of cultured embryos to colony-stimulating factor 2 (CSF2) alters DNA methylation. METHODS: Female fetuses were produced by artificial insemination or transfer of an IVP embryo. Embryos were treated from days 5 to 7 after fertilization with CSF2 or vehicle. DNA methylation in fetal liver and skeletal muscle was determined by post-bisulfite adaptor tagging-based sequencing. The degree of DNA methylation for CpG sites in 50-bp windows of the promoter region 500 bp upstream of the transcriptional start site was compared between treatments. RESULTS: For liver, there were 12 genes (6% of those analyzed) in which DNA methylation was affected by treatment, with one 50-bp window per gene affected by treatment. For muscle, the degree of DNA methylation was affected by treatment for 32 windows (19% of the total windows analyzed) representing 28 distinct genes (23% of analyzed genes). For 19 of the 28 genes in muscle, the greatest deviation in DNA methylation was for the CSF2 group. CONCLUSION: Results are consistent with alterations in the methylome being one of the mechanisms by which IVP can result in altered fetal development and postnatal function in the resultant offspring. In addition, results indicate that maternally derived cell-signaling molecules can regulate the pattern of DNA methylation.

Modeling allele-specific expression at the gene and SNP levels simultaneously by a Bayesian logistic mixed regression model.

BACKGROUND: High-throughput sequencing experiments, which can determine allele origins, have been used to assess genome-wide allele-specific expression. Despite the amount of data generated from high-throughput experiments, statistical methods are often too simplistic to understand the complexity of gene expression. Specifically, existing methods do not test allele-specific expression (ASE) of a gene as a whole and variation in ASE within a gene across exons separately and simultaneously. RESULTS: We propose a generalized linear mixed model to close these gaps, incorporating variations due to genes, single nucleotide polymorphisms (SNPs), and biological replicates. To improve reliability of statistical inferences, we assign priors on each effect in the model so that information is shared across genes in the entire genome. We utilize Bayesian model selection to test the hypothesis of ASE for each gene and variations across SNPs within a gene. We apply our method to four tissue types in a bovine study to de novo detect ASE genes in the bovine genome, and uncover intriguing predictions of regulatory ASEs across gene exons and across tissue types. We compared our method to competing approaches through simulation studies that mimicked the real datasets. The R package, BLMRM, that implements our proposed algorithm, is publicly available for download at https://github.com/JingXieMIZZOU/BLMRM . CONCLUSIONS: We will show that the proposed method exhibits improved control of the false discovery rate and improved power over existing methods when SNP variation and biological variation are present. Besides, our method also maintains low computational requirements that allows for whole genome analysis.

Altered microRNA expression profiles in large offspring syndrome and Beckwith-Wiedemann syndrome.

The use of assisted reproductive technologies (ART) can induce a congenital overgrowth condition in humans and ruminants, namely Beckwith-Wiedemann syndrome (BWS) and large offspring syndrome (LOS), respectively. Shared phenotypes and epigenotypes have been found between BWS and LOS. We have observed global misregulation of transcripts in bovine foetuses with LOS. microRNAs (miRNAs) are important post-transcriptional gene expression regulators. We hypothesize that there is miRNA misregulation in LOS and that this misregulation is shared with BWS. In this study, small RNA sequencing was conducted to investigate miRNA expression profiles in bovine and human samples. We detected 407 abundant known miRNAs and predicted 196 putative miRNAs from the bovine sequencing results and identified 505 abundant miRNAs in human tongue. Differentially expressed miRNAs (DE-miRNAs) were identified between control and LOS groups in all tissues analysed as well as between BWS and control human samples. DE-miRNAs were detected from several miRNA clusters including DLK1-DIO3 genomic imprinted cluster in LOS and BWS. DNA hypermethylation was associated with downregulation of miRNAs in the DLK1-DIO3. mRNA targets of the DE-miRNAs were predicted and signalling pathways associated with control of organ size (including the Hippo signalling pathway), cell proliferation, apoptosis, cell survival, cell cycle, and cell adhesion were found to be enriched with these genes. Yes associated protein 1 (YAP1) is the core effector of the Hippo signalling pathway, and increased level of active (non-phosphorylated) YAP1 protein was detected in skeletal muscle of LOS foetuses. Overall, our data provide evidence of miRNA misregulation in LOS and BWS.

Differential gene expression in response to eCry3.1Ab ingestion in an unselected and eCry3.1Ab-selected western corn rootworm (Diabrotica virgifera virgifera LeConte) population.

Diabrotica virgifera virgifera LeConte, the western corn rootworm (WCR) is one of the most destructive pests in the U.S. Corn Belt. Transgenic maize lines expressing various Cry toxins from Bacillus thuringiensis have been adopted as a management strategy. However, resistance to many Bt toxins has occurred. To investigate the mechanisms of Bt resistance we carried out RNA-seq using Illumina sequencing technology on resistant, eCry3.1Ab-selected and susceptible, unselected, whole WCR neonates which fed on seedling maize with and without eCry3.1Ab for 12 and 24 hours. In a parallel experiment RNA-seq experiments were conducted when only the midgut of neonate WCR was evaluated from the same treatments. After de novo transcriptome assembly we identified differentially expressed genes (DEGs). Results from the assemblies and annotation indicate that WCR neonates from the eCry3.1Ab-selected resistant colony expressed a small number of up and down-regulated genes following Bt intoxication. In contrast, unselected susceptible WCR neonates expressed a large number of up and down-regulated transcripts in response to intoxication. Annotation and pathway analysis of DEGs between susceptible and resistant whole WCR and their midgut tissue revealed genes associated with cell membrane, immune response, detoxification, and potential Bt receptors which are likely related to eCry3.1Ab resistance. This research provides a framework to study the toxicology of Bt toxins and mechanism of resistance in WCR, an economically important coleopteran pest species.

Lean maternal hyperglycemia alters offspring lipid metabolism and susceptibility to diet-induced obesity in mice†.

We previously developed a model of gestational diabetes mellitus (GDM) in which dams exhibit glucose intolerance, insulin resistance, and reduced insulin response to glucose challenge only during pregnancy, without accompanying obesity. Here, we aimed to determine how lean gestational glucose intolerance affects offspring risk of metabolic dysfunction. One cohort of offspring was sacrificed at 19 weeks, and one at 31 weeks, with half of the second cohort placed on a high-fat, high-sucrose diet (HFHS) at 23 weeks. Exposure to maternal glucose intolerance increased weights of HFHS-fed offspring. Chow-fed offspring of GDM dams exhibited higher body fat percentages at 4, 12, and 20 weeks of age. At 28 weeks, offspring of GDM dams fed the HFHS but not the chow diet (CD) also had higher body fat percentages than offspring of controls (CON). Exposure to GDM increased the respiratory quotient (Vol CO2/Vol O2) in offspring. Maternal GDM increased adipose mRNA levels of peroxisome proliferator-activated receptor gamma (Pparg) and adiponectin (Adipoq) in 31-week-old CD-fed male offspring, and increased mRNA levels of insulin receptor (Insr) and lipoprotein lipase (Lpl) in 31-week-old male offspring on both diets. In liver at 31 weeks, mRNA levels of peroxisome proliferator-activated receptor alpha (Ppara) were elevated in CD-fed male offspring of GDM dams, and male offspring of GDM dams exhibited higher mRNA levels of Insr on both diets. Neither fasting insulin nor glucose tolerance was affected by exposure to GDM. Our findings show that GDM comprising glucose intolerance only during pregnancy programs increased adiposity in offspring, and suggests increased insulin sensitivity of subcutaneous adipose tissue.

A Bayesian hidden Markov model for detecting differentially methylated regions.

Alterations in DNA methylation have been linked to the development and progression of many diseases. The bisulfite sequencing technique presents methylation profiles at base resolution. Count data on methylated and unmethylated reads provide information on the methylation level at each CpG site. As more bisulfite sequencing data become available, these data are increasingly needed to infer methylation aberrations in diseases. Automated and powerful algorithms also need to be developed to accurately identify differentially methylated regions between treatment groups. This study adopts a Bayesian approach using the hidden Markov model to account for inherent dependence in read count data. Given the expense of sequencing experiments, few replicates are available for each treatment group. A Bayesian approach that borrows information across an entire chromosome improves the reliability of statistical inferences. The proposed hidden Markov model considers location dependence among genomic loci by incorporating correlation structures as a function of genomic distance. An iterative algorithm based on expectation-maximization is designed for parameter estimation. Methylation states are inferred by identifying the optimal sequence of latent states from observations. Real datasets and simulation studies that mimic the real datasets are used to illustrate the reliability and success of the proposed method.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa.

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Detecting differentially expressed genes for syndromes by considering change in mean and dispersion simultaneously.

BACKGROUND: Using next-generation sequencing technology to measure gene expression, an empirically intriguing question concerns the identification of differentially expressed genes across treatment groups. Existing methods aim to identify genes whose mean expressions differ among treatment groups by assuming equal dispersion across all groups. For syndromes, however, various combinations of gene expression alterations can result in the same disease, leading to greater heteroscedasticity in the biological replicates in the disease group compared to the normal group. Traditional methods that only consider changes in the mean will fail to fully analyze gene expression in such a scenario. In addition, sequencing technology is relatively expensive; most labs can only afford a few replicates per treatment group, which poses further challenges to reliably estimating the mean and dispersion under each treatment condition. RESULTS: We designed an empirical Bayes method and a pooled permutation test to simultaneously consider the change in mean and dispersion across treatment groups. We further computed confidence intervals based on Bayes estimates to identify differentially expressed genes that are unique to each disease sample as well as those that are common across all disease samples. We illustrated our method by applying it to gene expression data from a large offspring syndrome experiment, which motivated this study. We compared our method to competing approaches through simulation studies that mimicked the real datasets to demonstrate the effectiveness of our proposed method. CONCLUSIONS: We will show that, compared to popular methods that only aim to find the difference in the mean, our method can capture greater variation in the disease group to effectively identify differentially expressed genes for syndromes.

Porcine Fetal-Derived Fibroblasts Alter Gene Expression and Mitochondria to Compensate for Hypoxic Stress During Culture.

The Warburg effect is characterized by decreased mitochondrial oxidative phosphorylation and increased glycolytic flux in adequate oxygen. The preimplantation embryo has been described to have characteristics of the Warburg effect, including similar changes in gene expression and mitochondria, which are more rudimentary in appearance. We hypothesized hypoxia would facilitate anaerobic glycolysis in fibroblasts thereby promoting gene expression and media metabolite production reflecting the Warburg effect hallmarks in early embryos. Additionally, we speculated that hypoxia would induce a rudimentary small mitochondrial phenotype observed in several cell types evidenced to demonstrate the Warburg effect. While many have examined the role hypoxia plays in pathological conditions, few studies have investigated changes in primary cells which could be used in somatic cell nuclear transfer. We found that cells grown in 1.25% O2 had normal cell viability and more, but smaller mitochondria. Several hypoxia-inducible genes were identified, including seven genes for glycolytic enzymes. In conditioned media from hypoxic cells, the quantities of gluconolactone, cytosine, and uric acid were decreased indicating higher consumption than control cells. These results indicate that fibroblasts alter gene expression and mitochondria to compensate for hypoxic stress and maintain viability. Furthermore, the metabolic changes observed, making them more similar to preimplantation embryos, could be facilitating nuclear reprogramming making these cells more amendable to future use in somatic cell nuclear transfer.

Glutamine supplementation enhances development of in vitro-produced porcine embryos and increases leucine consumption from the medium.

Improper composition of culture medium contributes to reduced viability of in vitro-produced embryos. Glutamine (Gln) is a crucial amino acid for preimplantation embryos as it supports proliferation and is involved in many different biosynthetic pathways. Previous transcriptional profiling revealed several upregulated genes related to Gln transport and metabolism in in vitro-produced porcine blastocysts compared to in vivo-produced counterparts, indicating a potential deficiency in the culture medium. Therefore, the objective of this study was to determine the effects of Gln supplementation on in vitro-produced porcine embryo development, gene expression, and metabolism. Cleaved embryos were selected and cultured in MU2 medium supplemented with 1 mM Gln (control), 3.75 mM Gln (+Gln), 3.75 mM GlutaMAX (+Max), or 3.75 mM alanine (+Ala) until day 6. Embryos cultured with +Gln or +Max had increased development to the blastocyst stage and total number of nuclei compared to the control (P < 0.05). Moreover, expression of misregulated transcripts involved in glutamine and glutamate transport and metabolism was corrected when embryos were cultured with +Gln or +Max. Metabolomics analysis revealed increased production of glutamine and glutamate into the medium by embryos cultured with +Max and increased consumption of leucine by embryos cultured with +Gln or +Max. As an indicator of cellular health, mitochondrial membrane potential was increased when embryos were cultured with +Max which was coincident with decreased apoptosis in these blastocysts. Lastly, two embryo transfers by using embryos cultured with +Max resulted in viable piglets, confirming that this treatment is consistent with in vivo developmental competence.