Disruption of Mitochondrion-To-Nucleus Interaction in Deceased Cloned Piglets.

Most animals produced by somatic cell nuclear transfer (SCNT) are heteroplasmic for mitochondrial DNA (mtDNA). Oxidative phosphorylation (OXPHOS) in clones therefore requires the coordinated expression of genes encoded by the nuclear DNA and the two sources of mitochondria. Such interaction is rarely studied because most clones are generated using slaughterhouse oocytes of unrecorded origin. Here we traced the maternal lineages of seven diseased and five one-month-old live cloned piglets by sequencing their mtDNA. Additionally by using a 13K oligonucleotide microarray, we compared the expression profiles of nuclear and mtDNA-encoded genes that are involved in mitochondrial functions and regulation between the cloned groups and their age-matched controls (n=5 per group). We found that the oocytes used to generate the cloned piglets were of either the Large White or Duroc background, and oocyte genetic background was not related to the clones' survival. Expression profiles of mtDNA-encoded genes in clones and controls showed intermixed clustering patterns without treatment or maternal lineage-dependency. In contrast, clones and controls clustered separately for their global and nuclear DNA-encoded mitochondrial genes in the lungs for both the deceased and live groups. Functional annotation of differentially expressed genes encoded by both nuclear and mtDNA revealed abnormal gene expression in the mitochondrial OXPHOS pathway in deceased clones. Among the nine differentially expressed genes of the OXPHOS pathway, seven were down-regulated in deceased clones compared to controls, suggesting deficiencies in mitochondrial functions. Together, these data demonstrate that the coordination of expression of mitochondrial genes encoded by nuclear and mtDNA is disrupted in the lung of diseased clones.

Introduction

Since the first report of cloning success in mammals by SCNT [1], mammalian oocytes have been recognized as possessing all the necessary components for complete reprogramming [2]. The number of healthy births, however, is still low and most cloned animals suffer from a variety of abnormal phenotypes [1, 3, 4]. Nevertheless, SCNT is still more efficient in reprogramming capacity than the induced pluripotent stem cell technology, which reprograms in the range of 0.1%.

In the last decade, efforts have been made to compare cloned animals to those of conventional breeding in order to further improve cloning efficiency [5-7]. A major area of difference lies in mitochondria which are important in ATP generation and in programmed cell death [8]. Conventionally bred animals are homoplasmic for the maternal mitochondria, while most SCNT animals are heteroplasmic, containing mitochondria from both the recipient oocytes and the somatic donor cells [9-12].

The mitochondrial genome is composed of 37 genes encoding 13 subunits of the respiratory enzyme complexes in the OXPHOS pathway as well as 22 tRNAs and two rRNAs. The majority of the elements in the OXPHOS pathway as well as other proteins involved in mitochondrial regulation and biogenesis (approximately 1,500 proteins) is encoded by nuclear DNA [13, 14]. Therefore, the coordinated expression of the nuclear and mtDNA and the interaction of the protein products are critical in maintaining normal mitochondrial and host cell functions. Mitochondrial dysfunction also feeds back to the host cell by altering expression of numerous nuclear DNA-encoded genes to readjust its metabolic profile. For example, in yeast, retrograde response from mitochondria to the nucleus influences many cellular activities under both normal and pathophysiological conditions [15]. Likewise, mitochondrial OXPHOS is defective in cytoplasmic hybrids containing rat mtDNA in mouse cells due to impairment of coordinated assembly of nuclear- and mtDNA-encoded OXPHOS subunits [16]. In addition, mutations in mtDNA can induce genomic instability as well as tumorigenesis [17].

Oocytes used in cloning farm animals are usually obtained from slaughterhouses [18], and thus animals cloned from the same donor cell line may not be entirely identical due to the different mtDNA. Heteroplasmy may interfere with the development of the cloned embryos and thus cloning outcome due to improper interaction of the donor nuclear DNA and recipient oocyte’s mitochondria in both intra- and, more likely, inter-species nuclear transfer. Several indirect lines of evidence support such a possibility. For example, reconstructed cytoplasmic hybrid embryos from Bos Taurus (Brown Swiss) granulosa cells and oocytes that contained B. taurus A (Simmental), B. taurus B (Simmental), or Bos indicus (Dwarf Zebu) cytoplasm have different developmental potentials in bovine [19]. Additionally, inter-subspecies SCNT using cytoplasts from oocytes of crossbred goats (Saanen male symbol x Boer female symbol descendant) improves embryo-fetal development than cytoplasts from Saanen oocytes [20]. Furthermore, higher efficiencies were observed in homoplasmic cloning in both cattle and sheep [21, 22]. However, whether mtDNA-encoded genes are coordinately expressed with those encoded by the nucleus has not been studied in clones.

In this study, we used DNA sequencing to identify the maternal lineages of oocytes in cloned pigs and subsequently analyzed the expression profiles of global, nuclear DNA- and mtDNA-encoded genes involved in OXPHOS in deceased and live cloned pigs. We found that expression profiles of mtDNA-encoded genes of clones and controls clustered independently from treatment or maternal lineage. Clones and controls clustered separately for their global and nuclear DNA-encoded mitochondrial genes in the lungs for both deceased and live groups. Functional annotation analyses of differentially expressed genes encoded by the nuclear DNA and mtDNA genes revealed defects in the mitochondrial OXPHOS pathway in deceased clones.

Materials and Methods

RNA and DNA isolation

All animal experimentation was conducted at the University of Missouri. All procedures were approved by the University of Missouri Institutional Animal Care and Use Committee. Veterinary staff was consulted on a daily basis regarding the health and care of the cloned piglets, humane sacrifice, and tissue collection.

Frozen lung tissues of cloned piglets by SCNT and of age-matched controls by conventional breeding were provided by Dr. Randall Prather of the University of Missouri-Columbia [4, 23]. Seven clones that died shortly after birth were designated as the “deceased” group, and five one-month old clones as the “live” group. Each cloned group was age-matched with five controls from conventional reproduction. Total RNA was isolated using TRIzol (Invitrogen, Carlsbad, CA) and treated with DNase I (Invitrogen) to remove any possible genomic DNA contamination. Isolated RNA was stored at -80°C until utilization. Total DNA was isolated by back extraction using back extraction buffer (4 M guanidine thiocyanate, 50 mM sodium citrate, 1M Tris-base) from the inter- and organic phases of TRIzol extraction, followed by RNase A (Invitrogen) treatment to remove any possible contamination of RNA. Isolated DNA was stored at -20°C until further analysis.

Mitochondrial DNA (mtDNA) sequencing

A 595-bp fragment of the mitochondrial NADH dehydrogenase 1 (ND1, 957 bp in length) was amplified and sequenced to determine genetic background of oocytes (“maternal”) used for cloning. The following PCR primers were designed based on the Sus scrofa mt DNA (GenBank access number AF034253): forward 5’-CCT ACT GGC CGT AGC ATT-3’, reverse 5’- GAA TCG TGG GTA TGA TGC TC-3’. PCR was performed using the Taq PCR kit (New England Biolab, Ipswich, MA). Briefly, 5 ng of total DNA in a 50 μl reaction mixture (1X standard Taq reaction buffer, 0.2 mM dNTP mixture, 0.75 μM of each primer, and 2.5 U Taq DNA Polymerase) was subjected to an initial denaturation at 94°C for 3 min, followed by 30 cycles of denaturation at 94°C for 20 sec, annealing at 60°C for 20 sec, extension at 72°C for 30 sec, and last extension at 72°C for 5 min. After purification by using ExoSAP-IT (USB, Cleveland, OH), the PCR products were sequenced by using the BigDye Terminator v1.1 Cycle Sequencing kit with ABI PRISM Model 3700 (Applied Biosystems Inc., Foster City, CA). Briefly, 30 μM of purified PCR products in a 20 μl reaction mixture (1X Terminator Ready Reaction Mix, 1X BigDye Sequencing Buffer, and 0.75 μM forward or reverse primer) were subjected to an initial denaturation at 94°C for 5 min, followed by 25 cycles of denaturation at 96°C for 10 sec, annealing at 50°C for 15 sec, and extension at 60°C for 4 min. The reaction mixture was then purified with DyeEx 2.0 (Qiagen, Valencia, CA), and electrophoresis was performed at Biotechnology/Bioservice Center at the University of Connecticut. Sequence similarity was measured by Basic Local Alignment Search Tool (BLAST, http://blast.ncbi.nlm.nih.gov/Blast.cgi). Bit score from the in silico sequence BLAST was used to determine maternal lineage by ND1 sequence, and lineage-specific mtDNA polymorphisms were identified with Sequencher (version 4.1.4, Gene Codes, Ann Arbor, MI) by comparing with pig mtDNA reference sequences.

DNA Microarray

The porcine Array-Ready Oligo Set (version 1.1, Operon Technologies, Alameda, CA) [24] was used in this study. Out of the 13,310 total probes, 566 genes were involved in mitochondrial functions including 534 and 32 genes from nuclear and mitochondrial DNA, respectively. Reference RNA was isolated and pooled from brain, kidney, liver, and lung tissues of male and female pig fetuses [25, 26]. Aminoallyl-linked cDNA was generated from 10 μg of the reference and extracted total RNA by reverse transcription, followed by labeling with Cy3 or Cy5, and hybridized to microarray slides as described [25, 27]. A total of 44 microarrays were used by co-hybridizing reference and extracted RNA with dye-swap. After scanning, local background subtraction was performed in GenePix Pro 6.0 (Molecular Devices, Union City, CA), and probes were flagged as “present” if they had at least 70% of feature pixels with more than two standard deviations (SD) above the background in either the Cy3 or Cy5 channel. Signal intensities of the “present probes” were loaded into GeneSpring 6.1 (Agilent Technologies, Palo Alto, CA). All microarrays were normalized by the Lowess method [28]. Probes were assigned as “informative” when they were present in either the reference or sample in more than 90% of the microarrays with raw expression values > 100, and a SD < 1.4 [29]. Informative probes were used for further analysis. Hierarchical clustering was conducted with three different informative probe sets: 1) all informative probes (global); informative probes involved in mitochondrial functions and encoded by the 2) nuclear DNA and 3) mtDNA. Differentially expressed (DE) genes were defined as those with a fold change ≥ 1.5 and a p-value ≤ 0.05 (Welch t-test adjusted by using Benjamini and Hochberg false discovery rate). Functional annotation was followed with RefSeq IDs from human orthologs of the DE genes using the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov) [30].

Quantitative real-time RT-PCR (qPCR)

Nine genes involved in mitochondrial functions, of which 2 are encoded by mtDNA, were selected for qPCR validation. Using GenBank or Sus Scrofa Gene Index (SsGI) release 12.0, primers were designed to generate PCR products that harbored probe sequences of the microarray (Table 1). Tubulin beta-2 (TBB2) and the reference RNA were used as the internal control and as the calibrator, respectively. qPCR was performed with ABI 7500 Fast Real-Time PCR System (Applied Biosystems). Briefly, 2 ng of cDNA in a 20 μl reaction mixture [1X SYBR Green PCR master mix (Applied Biosystems), 0.3 μM of each primer] was subjected to an initial denaturation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 sec and annealing/extension at 60°C for 1 min. Data were analyzed using ABI prism SDS software (Applied Biosystems), and the relative expression levels of the target genes were calculated using the comparative CT method according to the manufacturer’s instruction. Parametric data were tested for equal variance by using an F-test, and when it was not significant, data were subject to one-sided t-test (p < 0.05).

Table 1

Primer sequences for qPCR.

Gene	Type	Sequence	Tm	Amplicon (bp)	SsGI Release 12.0 accession number
ACSL1	Forward Primer	AGTCTTCACCCTGAATTATTCTCCAT	59	80	TC240600
	Reverse Primer	TTCCGAAGCTCTGGCCTTTT	60
CEBPA	Forward Primer	GGACCACCAGTCTTGTCTGTACTG	59	80	TC292331
	Reverse Primer	TCCCTCTGGGATGGACTGATT	59
COX4I1	Forward Primer	GCAGGATGTTGGCTACCAGAGT	59	92	TC275877
	Reverse Primer	CTTCACGACGCTTCCATGTG	59
COX7C	Forward Primer	GCGCAATGTTGGGACAGAGT	60	79	TC298342
	Reverse Primer	TTCCCTGGACCCTCCTCATAG	59
GLUL	Forward Primer	GGTGGGAAATGAGGTAGGAAAAT	58	81	TC162304
	Reverse Primer	CGAACATAATCAGGGCCTTTAGC	60
NDUFA4	Forward Primer	GGAGGTACTGGAGCAGCACTGTAT	60	74	TC182398
	Reverse Primer	TTCTTCCTGTCCCAACAGACATC	59
NOTCH3	Forward Primer	ACTTCACTGCATTCCAGATGGA	58	111	TC257056
	Reverse Primer	GCCAGTTCCCAAAGGGATTC	59
P38	Forward Primer	GGGCCAAGGTGTCTCCATTT	60	96	TC287045
	Reverse Primer	CCTTCCTTCTCGCTCCAGTTG	60
TBB2	Forward Primer	ACCCGAGGGACCTCTTTATTCA	60	84	TC279724
	Reverse Primer	AGACAGGAAGAAGAACAGATACACACA	59
UQCR	Forward Primer	ATCTTCCTTCTTAAACTTGCCATTG	58	76	TC182529
	Reverse Primer	GCTGGTATGGGCCACTGACT	59

Results

Determination of the “maternal background” of the cloned piglets

All microarray data from this study have been submitted to NCBI under the accession number of GSE68877. All sequences of mtDNA obtained in the current study were of either the Large White or the Duroc breed. For convenience, the ND1 transcription start site was designated as Base #1 in the sequenced fragment. Seven breed-specific ND1 gene polymorphisms between the Large White and Duroc (Table 2; T342C, T369C, T420C, G438A, G459A, C489A, and C712T from the Large White to the Duroc) were found in the 371-bp fragment analyzed. Alignment with mtDNA reference sequences confirmed lineage-specific ND1 gene polymorphisms at the expected base positions as well. None of these are anticipated to cause translational variants.

Table 2

Mitochondrial ND1 gene sequence analysis.

Group		Base comparison									BLAST hit	Matched strain	Bit score
Deceased	Base position	316	342	369	420	438	459	489	528	712
	Reference	C	T a	T a	T a	G a	G a	C a	T	C a	AF486874	Large White
	Control ID
	263–1	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1172
	263–2	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1172
	263–3	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1172
	263–4	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1172
	263–5	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1172
	Clone ID
	628–1	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1179
	629–1	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1179
	635–1	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1179
	635–2	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	637–1	C	T	T	T	G	G	C	C b	C	AF486874	Large White	1179
	642–1	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	671–5	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
Live	Base position	316	342	369	420	438	459	489	528	712
	Reference	C	C a	C a	C a	A a	A a	A a	T	T a	AY337045	Duroc
	control ID
	42–7	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	42–8	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	42–11	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	42–12	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	47–2	C	C	C	C	A	A	A	T	T	AY337045	Duroc	1179
	Clone ID
	43–1	C	T	T	T	G	G	C	T	C	AF486874	Large White	1179
	43–2	C	T	T	T	G	G	C	T	C	AF486874	Large White	1179
	43–3	C	T	T	T	G	G	C	T	C	AF486874	Large White	1179
	44–1	T2	C	C	C	A	A	A	T	T	AY337045	Duroc	1172
	45–1	C	T	T	T	G	G	C	T	C	AF486874	Large White	1179

aStrain-specific DNA polymorphism

bNonsense mutation

We then used the mtDNA sequences of the two breeds to determine the maternal background of the clones and their respective controls. The controls in the deceased and the live groups matched to the Large White (AF486874, bit score = 1,172), and Duroc (AY337045, bit score = 1,179), respectively. For the deceased clones, four gave the best match against the Large White (bit score = 1,179) and the other three against Duroc (bit score = 1,179). Four of the live clones matched to the Large White (bit score = 1,179) and the other to Duroc (bit score = 1,172).

Hierarchical clustering of expression profiles of global, nuclear DNA- and mtDNA-encoded genes involved in mitochondrial functions

Expression profiles of global, nuclear DNA- and mtDNA-encoded genes of mitochondrial functions of the clones were compared with those of their conventionally bred counterparts (Fig 1). In the deceased group, global expression of 9,297 informative probes (69.8%) was largely clustered by treatment (clone vs. control). The correlation coefficient (r) of the controls was 0.925 compared with that of the clones at r = 0.673, demonstrating more variability in the cloned group. Likewise, 534 (4%) nuclear DNA-encoded genes of mitochondrial functions from the informative probes revealed that samples were clustered by treatment with r = 0.924 and 0.638 for controls and clones, respectively. Neither the global nor the nuclear DNA-encoded genes clustered with the maternal background. In contrast, the 32 (0.2%) mtDNA-encoded genes did not cluster by treatment or breed with the lowest correlation at r = 0.473.

Fig 1

Hierarchical clustering of global, nuclear DNA (nDNA)- and mitochondrial DNA (mtDNA)-encoded genes of deceased (top panel) and live clones (bottom panel) vs. age-matched controls.

From left to right, clustering with all informative probes, nuclear DNA- and mitochondrial DNA-encoded probes involved in mitochondrial functions were depicted. r = correlation coefficient, red line = clones, blue line = controls, red animal ID = Duroc, blue animal ID = Large White.

In the live clone group, 9,955 probes (74.8%) were found to be informative for further analysis. As is with the dead group, hierarchical clustering of the informative probes also separated samples by treatment with r = 0.871 and 0.913 for controls and clones, respectively. Similarly, 451 nuclear DNA-encoded genes for mitochondrial functions (3.4%) were clustered by treatment with r = 0.879 and 0.921 for controls and clones, respectively. For the 32 mtDNA-encoded genes (0.2%), most samples separated by treatment with the exception of one clone of the Duroc background (Animal ID: 44–1) and a control of the Large White background (Animal ID: 47–2).

Functional annotation of differentially expressed genes involved in mitochondrial functions

In the deceased group, 46 differentially expressed (DE) genes were identified, among which 43 were encoded by the nuclear DNA and 3 by the mtDNA. Seventeen of the nuclear-encoded genes were up-regulated and 26 down-regulated compared to controls (Table 3). Functional annotation of these DE genes revealed that the following nine involved in the OXPHOS pathway (Fig 2): NADH dehydrogenase (ubiquinone) 1 alpha subcomplex 4 (NDUFA4), ubiquinol-cytochrome c reductase complex III subunit VII (UQCRQ), ubiquinol-cytochrome c reductase (UQCR), ATPase V1 subunit E1 (ATP6V1E1), ATP synthase F1 complex beta polypeptide (ATP5B), and ATP synthase F0 complex subunit C1 (ATP5G1), all encoded by the nuclear DNA, as well as cytochrome c oxidase subunit 1 (COX1), cytochrome c oxidase subunit 2 (COX2), and NADH dehydrogenase 3 (ND3), all encoded by the mtDNA.

Table 3

Differentially expressed mitochondrial genes in deceased cloned piglets.

Encoded from	GenBank ID	TC ID	Fold change c	Gene Symbol	Gene Name
nDNA a	NM_004169	TC182203	3.65	SHMT1	SERINE HYDROXYMETHYLTRANSFERASE 1
	NM_004052	TC182058	2.64	BNIP3	BCL2/ADENOVIRUS E1B 19KDA INTERACTING PROTEIN 3
	NM_002065	TC162304	2.12	GLUL	GLUTAMATE-AMMONIA LIGASE (GLUTAMINE SYNTHETASE)
	NM_000817	TC189704	2.07	GAD1	GLUTAMATE DECARBOXYLASE 1
	NM_000274	TC181897	2.05	OAT	ORNITHINE AMINOTRANSFERASE
	NM_000663	TC186487	1.97	ABAT	4-AMINOBUTYRATE AMINOTRANSFERASE
	NM_002612	TC167221	1.91	PDK4	PYRUVATE DEHYDROGENASE KINASE, ISOZYME 4
	NM_001437	TC169575	1.88	ESR2	ESTROGEN RECEPTOR 2 (ER BETA)
	NM_001006642	TC166867	1.87	SLC25A25	SOLUTE CARRIER FAMILY 25 (MITOCHONDRIAL CARRIER; PHOSPHATE CARRIER), MEMBER 25
	NM_021232	N/A	1.73	PRODH2	PROLINE DEHYDROGENASE (OXIDASE) 2
	NM_001696	TC162504	1.62	ATP6V1E1	ATPASE, H+ TRANSPORTING, LYSOSOMAL 31KDA, V1 SUBUNIT E1
	NM_000349	TC162876	1.61	STAR	STEROIDOGENIC ACUTE REGULATOR
	NM_012094	TC164220	1.57	PRDX5	PEROXIREDOXIN 5
	NM_022036	TC164081	1.56	GPRC5C	G PROTEIN-COUPLED RECEPTOR, FAMILY C, GROUP 5, MEMBER C
	NM_020166	TC186613	1.55	MCCC1	METHYLCROTONOYL-COENZYME A CARBOXYLASE 1 (ALPHA)
	NM_004415	TC188422	1.53	DSP	DESMOPLAKIN
	NM_005763	TC190934	1.52	AASS	AMINOADIPATE-SEMIALDEHYDE SYNTHASE
	NM_002080	TC163473	0.67	GOT2	GLUTAMIC-OXALOACETIC TRANSAMINASE 2, MITOCHONDRIAL (ASPARTATE AMINOTRANSFERASE 2)
	NM_006111	TC185178	0.66	ACAA2	ACETYL-COENZYME A ACYLTRANSFERASE 2 (MITOCHONDRIAL 3-OXOACYL-COENZYME A THIOLASE)
	NM_014297	TC181027	0.66	ETHE1	ETHYLMALONIC ENCEPHALOPATHY 1
	NM_000308	TC182001	0.66	PPGB	PROTECTIVE PROTEIN FOR BETA-GALACTOSIDASE
	NM_000397	TC163629	0.66	CYBB	CYTOCHROME B-245, BETA POLYPEPTIDE
	NM_006793	TC181710	0.65	PRDX3	PEROXIREDOXIN 3
	NM_003679	TC168142	0.65	KMO	KYNURENINE 3-MONOOXYGENASE (KYNURENINE 3-HYDROXYLASE)
	NM_015469	TC188504	0.65	NIPSNAP3A	NIPSNAP HOMOLOG 3A
	NM_000016	TC183520	0.64	ACADM	ACYL-COENZYME A DEHYDROGENASE, C-4 TO C-12 STRAIGHT CHAIN
	NM_001686	TC162934	0.64	ATP5B	ATP SYNTHASE, H+ TRANSPORTING, MITOCHONDRIAL F1 COMPLEX, BETA POLYPEPTIDE
	NM_020312	TC182787	0.64	COQ9	COENZYME Q9 HOMOLOG
	NM_012259	TC168760	0.64	HEY2	HAIRY/ENHANCER-OF-SPLIT RELATED WITH YRPW MOTIF 2
	NM_002489	TC182398	0.64	NDUFA4	NADH DEHYDROGENASE (UBIQUINONE) 1 ALPHA SUBCOMPLEX, 4, 9KDA
	NM_014402	TC164365	0.63	UQCRQ	UBIQUINOL-CYTOCHROME C REDUCTASE, COMPLEX III SUBUNIT VII, 9.5KDA
	NM_002396	TC167924	0.63	ME2	MALIC ENZYME 2, NAD(+)-DEPENDENT, MITOCHONDRIAL
	NM_001985	TC163232	0.62	ETFB	ELECTRON-TRANSFER-FLAVOPROTEIN, BETA POLYPEPTIDE
	NM_014481	TC166660	0.62	APEX2	APEX NUCLEASE (APURINIC/APYRIMIDINIC ENDONUCLEASE) 2
	NM_006830	TC182529	0.62	UQCR	UBIQUINOL-CYTOCHROME C REDUCTASE, 6.4KDA SUBUNIT
	NM_020548	TC182240	0.61	DBI	DIAZEPAM BINDING INHIBITOR (GABA RECEPTOR MODULATOR, ACYL-COENZYME A BINDING PROTEIN)
	NM_012343	TC189901	0.58	NNT	NICOTINAMIDE NUCLEOTIDE TRANSHYDROGENASE
	NM_016497	TC165872	0.55	MRPL51	MITOCHONDRIAL RIBOSOMAL PROTEIN L51
	NM_005175	TC163315	0.55	ATP5G1	ATP SYNTHASE, H+ TRANSPORTING, MITOCHONDRIAL F0 COMPLEX, SUBUNIT C1 (SUBUNIT 9)
	NM_002979	TC163162	0.54	SCP2	STEROL CARRIER PROTEIN 2
	NM_006082	TC181275	0.52	TUBA6	TUBULIN, ALPHA, UBIQUITOUS
	NM_001876	TC164450	0.46	CPT1A	CARNITINE PALMITOYLTRANSFERASE 1A
	NM_001948	TC184267	0.45	DUT	DUTP PYROPHOSPHATASE
mtDNA b	NM_006818	TC168152	0.65	COX1	CYTOCHOROME C OXIDASE SUBUNIT 1
	N/A	TC181301	0.65	COX2	CYTOCHROME C OXIDASE SUBUNIT 2
	N/A	TC190307	2.11	ND3	NADH DEHYDROGENASE 3

anDNA: nuclear DNA

bmtDNA: mitochondrial DNA

cFold change = clone/age-matched control

Fig 2

Dysregulated oxidative phosphorylation pathway of deceased cloned piglets.

Genes in green boxes were down-regulated in clones (fold change > 1.5). Genes in orange box and red box were upregulated in clones by more than 1.5-, and 2.0-fold, respectively. Data were analyzed by using the Fisher’s exact test (p < 9.95 x 10–7, Benjamini and Hochberg FDR = 0.00019).

In the live clones, 35 nuclear DNA-encoded DE genes were found (Table 4). Among them, 24 were up-regulated and 11 were down-regulated compared to controls. Three mtDNA-encoded genes were also identified as differentially expressed. Functional analysis revealed that six nuclear genes were involved in the “FATTY ACID METABOLISM” pathway. These are dodecenoyl-coenzyme A delta isomerase (DCI), enoyl coenzyme A hydratase short chain 1 (ECHS1), glutaryl-coenzyme A dehydrogenase (GCDH), L-3-hydrozyacyl-coenzyme A dehydrogenase short chain (HADHSC), acetyl-coenzyme A acyltransferase 2 (ACAA2), and fatty-acid-coenzyme A ligase long chain 1 (ACSL1). However, no significant functional annotation was enriched with the DE genes encoded by mtDNA.

Table 4

Differentially expressed mitochondrial genes in live cloned piglets.

Encoded from	GenBank ID	TC ID	Fold change c	Gene Symbol	Gene Name
nDNA a	NM_014585	TC165023	2.13	SLC40A1	SOLUTE CARRIER FAMILY 40 (IRON-REGULATED TRANSPORTER), MEMBER 1
	NM_004331	TC184073	1.79	BNIP3L	BCL2/ADENOVIRUS E1B 19KDA INTERACTING PROTEIN 3-LIKE
	NM_170711	TC183283	1.70	DAZAP1	DAZ ASSOCIATED PROTEIN 1
	NM_001752	TC181995	1.69	CAT	CATALASE
	NM_012343	TC189901	1.65	NNT	NICOTINAMIDE NUCLEOTIDE TRANSHYDROGENASE
	NM_001919	TC181889	1.64	DCI	DODECENOYL-COENZYME A DELTA ISOMERASE (3,2 TRANS-ENOYL-COENZYME A ISOMERASE)
	NM_015917	TC161988	1.63	GSTK1	GLUTATHIONE S-TRANSFERASE SUBUNIT 13 HOMOLOG
	NM_002979	TC163162	1.63	SCP2	STEROL CARRIER PROTEIN 2
	NM_000456	TC185923	1.61	SUOX	SULFITE OXIDASE
	NM_005881	TC163821	1.59	BCKDK	BRANCHED CHAIN KETOACID DEHYDROGENASE KINASE
	NM_018838	TC164907	1.58	NDUFA12	NADH DEHYDROGENASE (UBIQUINONE) 1 ALPHA SUBCOMPLEX, 12
	NM_004092	TC181253	1.57	ECHS1	ENOYL COENZYME A HYDRATASE, SHORT CHAIN, 1, MITOCHONDRIAL
	NM_000159	TC182620	1.57	GCDH	GLUTARYL-COENZYME A DEHYDROGENASE
	NM_005327	TC183224	1.56	HADHSC	L-3-HYDROXYACYL-COENZYME A DEHYDROGENASE, SHORT CHAIN
	NM_005984	TC163540	1.56	SLC25A1	SOLUTE CARRIER FAMILY 25 (MITOCHONDRIAL CARRIER; CITRATE TRANSPORTER), MEMBER 1
	NM_000398	TC183476	1.55	CYB5R3	CYTOCHROME B5 REDUCTASE 3
	NM_006111	TC185178	1.54	ACAA2	ACETYL-COENZYME A ACYLTRANSFERASE 2 (MITOCHONDRIAL 3-OXOACYL-COENZYME A THIOLASE)
	NM_004453	TC166131	1.53	ETFDH	ELECTRON-TRANSFERRING-FLAVOPROTEIN DEHYDROGENASE
	NM_001359	TC182911	1.53	DECR1	2,4-DIENOYL COA REDUCTASE 1, MITOCHONDRIAL
	NM_001609	TC167600	1.52	ACADSB	ACYL-COENZYME A DEHYDROGENASE, SHORT/BRANCHED CHAIN
	NM_032549	TC167347	1.51	IMMP2L	IMP2 INNER MITOCHONDRIAL MEMBRANE PEPTIDASE-LIKE
	NM_002858	TC166623	1.51	ABCD3	ATP-BINDING CASSETTE, SUB-FAMILY D (ALD), MEMBER 3
	NM_012181	TC161966	1.50	FKBP8	FK506 BINDING PROTEIN 8, 38KDA
	NM_030791	TC164539	1.50	SGPP1	SPHINGOSINE-1-PHOSPHATE PHOSPHATASE 1
	NM_018947	TC163123	0.65	CYCS	CYTOCHROME C, SOMATIC
	NM_001995	TC163394	0.64	ACSL1	FATTY-ACID-COENZYME A LIGASE, LONG-CHAIN 1
	NM_018188	TC162178	0.63	ATAD3A	ATPASE FAMILY, AAA DOMAIN CONTAINING 3A
	NM_005729	TC182963	0.60	PPIF	PEPTIDYLPROLYL ISOMERASE F (CYCLOPHILIN F)
	NM_031212	TC182541	0.58	SLC25A28	SOLUTE CARRIER FAMILY 25, MEMBER 28
	NM_000032	TC166637	0.58	ALAS2	AMINOLEVULINATE, DELTA-, SYNTHASE 2
	NM_002065	TC165015	0.54	GLUL	GLUTAMATE-AMMONIA LIGASE (GLUTAMINE SYNTHETASE)
	NM_001006642	TC166867	0.48	SLC25A25	SOLUTE CARRIER FAMILY 25 (MITOCHONDRIAL CARRIER; PHOSPHATE CARRIER), MEMBER 25
	NM_199187	TC162826	0.41	KRT18	KERATIN 18
	NM_000636	TC164641	0.25	SOD2	SUPEROXIDE DISMUTASE 2, MITOCHONDRIAL
	NM_002612	N/A	0.20	PDK4	PYRUVATE DEHYDROGENASE KINASE, ISOZYME 4
mtDNA b	N/A	TC176636	1.60	ATP6	ATPASE SUBUNIT 6
	N/A	TC166635	0.66	COX3	CYTOCHROME C OXIDASE SUBUNIT 3
	NM_016357	TC167974	0.62	ND2	NADH DEHYDROGENASE SUBUNIT 2

anDNA: nuclear DNA

bmtDNA: mitochondrial DNA

cFold change = clone/age-matched control

The expression levels of nine genes were validated by qPCR (Fig 3). They are GLUL, CEBPA, NOTCH3 (up-regulated in clones), NDUFA4, UQCR, ACSL1 (down-regulated in clones), and COX4l1, COX7C, P38 (control level of expression in clones). Although slight differences in the fold changes of three genes were observed (GLUL, ACSL1, and P38), overall the qPCR results confirmed the directionality of the microarray.

Fig 3

Comparisons of representative gene expressions from qPCR and microarray.

Black bars indicate mean gene expression data from qPCR, and gray bars from microarray. Error bars mean ± SD. Data were analyzed by one-sided t-test (*p < 0.05, ***p < 0.001).

Discussion

Using sequence variations in the mtDNA-encoded ND1 gene, we first determined the “maternal background” of clones and controls used in this study and subsequently analyzed if expression of genes involved in mitochondrial functions were normally expressed in the deceased and the live clones or if they are related to their maternal lineages. We found that all cloned piglets had either the Large White or the Duroc as the maternal background. These breeds are the most common commercial pig lineages of the world and serve to produce many mixed-lineaged commercial hogs [31]. Although these European domestic pig lineages are closely related in their maternal lines, Large White is supposed to have separated from the European lineages around 58,000 year ago. In addition, the genetic distance between the Large White and the Duroc is 40-fold longer than that among European domestic lineages [32]; therefore, it should be feasible to distinguish between the maternal lineage of the Large White and the Duroc by using mitochondrial ND 1 gene sequence By sequence alignment, we found seven lineage-specific ND1 polymorphisms from the 371-bp target fragment. Interestingly, these polymorphisms were not supposed to induce translational variants. Therefore, these ND1 polymorphisms are good candidates for lineage determination while not affecting function s of the mitochondrial OXPHOS complex I. We did not find maternal lineage-dependent survivability of cloned piglets. Whether oocytes of different maternal breed background critically impact the cloning efficiency or survivability in pigs requires additional study and larger datasets.

The persistence of donor cell-derived mtDNA in cloned animals has been reported in four mammalian species including pigs [9]. The degree of heteroplasmy; however, has been found to be highly variable, ranging from 0 to 57% when donor cell was an embryonic blastomere determined by restriction digestion and Southern blotting [10]. Most somatic cell clones only contain a low fraction (<1%) of the donor cell mtDNA as determined by the more sensitive qPCR method [11]. In general, it is thought that the onset of physiological effects of pathogenic heteroplasmy begins when the ratio of mutant to normal mtDNA exceeds a specific threshold. For example, heteroplasmy-induced mitochondrial dysfunction contributes to pulmonary hypertension associated with cardiac output in mice. Quantification analysis revealed that most organs in these mice have more than 30% of mutant mtDNA [12]. In addition, a very low proportion of mutant heteroplasmy (less than 5%) is of questionable clinical significance since it may be due to laboratory error [33]. In preliminary trials of this study, we determined that a higher than 25% heteroplasmy would be detectable by sequencing the ND1. None of the cloned animals has such a high level of heteroplasmy (data not shown). Therefore, we excluded the possibility of pathogenesis contributed by donor cell-derived mtDNA.

Developmental competence of reconstructed embryos is significantly affected by the cytoplasm of the recipient oocytes. For example, oocytes from the hybrid strain BDF1 produce the highest cloning efficiency in the mouse [34]. Using oocytes of three different breeds and the same donor cell line, both transferable cloned embryos and viable bovine cloned fetuses were found to vary in phenotype based on different oocyte cytoplasm [19]. These cytoplasmic effects continued to affect birth weight, crump-rump length, and femoral length of cloned fetuses. In inter-subspecies SCNT embryos in goats, full term birth rates were highest in intra-species SCNT controls, followed by NT using recipient oocytes of crossbred goats with the lowest birth rate resulting from Boer donor cells into Saanen oocytes [20]. Furthermore, structure and quantity of mitochondria in the recipient oocytes have also been shown to influence the developmental potential of reconstructed bovine embryos [35]. For instance, oocytes with different mtDNA haplotype have a different copy number of mtDNA and ATP content. These oocytes have differential expression levels of COX1 and COX3, key components of the OXPHOS pathway and ATP generation, which is suspected to affect the cloning efficiency. These results shed light on the relevance of the maternal oocyte lineage to generate cloned animals by SCNT. In contrast, the maternal lineages revealed in this study were found to not be critical for embryo-fetal development or at least for viability of cloned piglets. Therefore, it is likely that the recipient oocytes retrieved at least from the Large White and/or the Duroc do not affect the cloning efficiency in pigs as opposed to observations in the mouse or bovine where the source of the oocyte is important.

Despite the lack of correlation of clone pig survival and mitochondrial lineages, gene expression analyses revealed an interruption of the interaction between genes encoded by the nuclear and mtDNA. Lack of expression coordination, especially abnormal expression of mtDNA-encoded genes, can affect nuclear DNA expression at the cellular level. For example, in the respiratory-deficient cell line lacking mtDNA (ρ°cells), the expression levels of numerous nuclear DNA-encoded genes for mitochondrial biogenesis and functions are dysregulated [36]. Kulawiec et al. reported that in ρ°cells, UQCRC1, a nuclear DNA-encoded gene, is down-regulated by as much as 10-fold. Furthermore, in both breast and ovarian carcinomas, UQCRC1 expression is positively correlated with mtDNA-encoded COX2 expression [17]. In this study, expression profiles of mtDNA-encoded genes in the live clone group showed treatment-dependent clustering. However, clustering analysis based on the expression profiles of mtDNA-encoded genes of the deceased clones showed an intermixed pattern with their age-matched control. These results imply that disruption in mitochondria-to-nucleus interaction occurred only in the deceased cloned piglets. Such lack of coordination may affect energy metabolism in the lungs and contribute to the overall deficiency in full development of clones.

Annotation analysis revealed that defects in the OXPHOS pathway could have contributed early death of cloned piglets. Several genes for OXPHOS pathway were differentially expressed in the deceased cloned piglets. For example, COX1 and COX2, two of the three mtDNA-encoded genes involved in OXPHOS complex IV, were significantly down-regulated only in the deceased clones. NDUFA4 and ND3 are complex I subunits encoded from nuclear and mtDNA, respectively. Mitochondrial genes like ND3 are down-regulated under pathological conditions such as glioblastoma [37] and hypoxia [38], and it was accompanied by the reduction of the activities of the respiratory enzyme complexes. Furthermore, down-regulation of these genes leads to loss of complex IV and complex I assemblies and activities [39, 40]. Another down-regulated gene in clones, UQCRQ, encodes a ubiquinone-binding protein (QP-C) as a subunit of complex III. A deletion in this gene decreases cytochrome b content which affects the assembly or maintenance of complex III [41]. ATP5B is a beta subunit of ATP synthase that catalyzes ATP formation using the energy of proton flux through the mitochondrial inner membrane during OXPHOS [42], and it is down-regulated in ulcerative colitis affecting the colon mucosa with other mitochondrial proteins [43]. Multiple copies of subunit C in the transmembrane portion (F0) of the ATP synthase transport rotate protons across the mitochondrial inner membrane to the F1-ATPase. ATP5G1 is one of the subunit C genes, and mature forms of ATP5G1, ATP5G2, and ATP5G3 are identical [44]. In alcohol consuming rats, all of these genes were dysregulated accompanying mitochondrial injury in the pancreas [45]. The abnormal expression of all of the genes above in deceased clones likely resulted in reduced ATP generation and pulmonary functionality.

In conclusion, this study demonstrates that the lineage of recipient oocytes may have little or no effect on the viability of cloned piglets by SCNT. However, cloned pigs may develop abnormal expression of both nuclear and mtDNA-encoded genes especially in the OXPHOS pathway, leading to a disruption of coordination of nuclear and mitochondrion interactions.