Development of Transgenic Minipigs with Expression of Antimorphic Human Cryptochrome 1

Minipigs have become important biomedical models for human ailments due to similarities in organ anatomy, physiology, and circadian rhythms relative to humans. The homeostasis of circadian rhythms in both central and peripheral tissues is pivotal for numerous biological processes. Hence, biological rhythm disorders may contribute to the onset of cancers and metabolic disorders including obesity and type II diabetes, amongst others. A tight regulation of circadian clock effectors ensures a rhythmic expression profile of output genes which, depending on cell type, constitute about 3–20% of the transcribed mammalian genome. Central to this system is the negative regulator protein Cryptochrome 1 (CRY1) of which the dysfunction or absence has been linked to the pathogenesis of rhythm disorders. In this study, we generated transgenic Bama-minipigs featuring expression of the Cys414-Ala antimorphic human Cryptochrome 1 mutant (hCRY1AP). Using transgenic donor fibroblasts as nuclear donors, the method of handmade cloning (HMC) was used to produce reconstructed embryos, subsequently transferred to surrogate sows. A total of 23 viable piglets were delivered. All were transgenic and seemingly healthy. However, two pigs with high transgene expression succumbed during the first two months. Molecular analyzes in epidermal fibroblasts demonstrated disturbances to the expression profile of core circadian clock genes and elevated expression of the proinflammatory cytokines IL-6 and TNF-α, known to be risk factors in cancer and metabolic disorders.


Introduction
Upholding an entrained biorhythm is pivotal for the timing of metabolic processes (reviewed in [1,2]). Hence, a highly developed hierarchical system of circadian clocks beginning with a master clock residing in the suprachiasmatic nucleus (SCN) of the anterior hypothalamus coordinates rhythmic cell behavior in all peripheral metabolic tissues [3,4]. Cryptochrome (CRYs) proteins belong to a class of light-sensitive flavoproteins that play a core role in the molecular pathways underlying circadian rhythms in mammalian cells (reviewed in [1]). The heterodimeric transcription factor composed of the circadian locomotor output cycles kaput (CLOCK) and the brain and muscle aryl hydrocarbon receptor nuclear translator (ARNT)-like protein 1 (BMAL1) controls the expression of E-box containing genes including those coding for CRY and periodic circadian (PER) proteins. In turn, CRY/PER heterodimers inhibit CLOCK/BMAL1 transactivation activity through a negative transcriptional/translational feed-back loop (TTFL). Timely post-translational modifications of CRY and PER, such as phosphorylation by casein kinase I, target them for proteosomal degradation thereby releasing inhibition of CLOCK/ BMAL1 in the auto-regulatory cycle [5,6].
Regulation of several physiological systems including the endocrine, rest-activity cycle and metabolic system, is heavily intertwined with the circadian clock. In fact, 3-20% of all peripherally active genes are expressed with a periodicity of 24hours and many of these are involved in metabolic processes [1,7,8,9]. Naturally, recent evidence suggests that disturbance of circadian rhythms increases the risk of metabolic disorders such as obesity and type II diabetes [10,11]. Moreover, mice with the sole modification of overexpressing an autologous but mutated CRY1 share key characteristics with type II diabetes mellitus patients [12].
Minipigs have become important biomedical models for human ailments due to the well described anatomical and physiological resemblance to humans. Hence, minipigs as an animal model for at least forty human ailments including cardiovascular and metabolic disorders have been described [13,14,15,16]. Creation of transgenic animals by somatic cell nuclear transfer (SCNT) has been described for a variety of animals including sheep [17], mouse [18,19], cow [20,21] and pig [22,23,24]. In brief, the nuclei of genetically manipulated somatic donor cells are introduced into enucleated oocytes. Subsequently, the reconstructed embryos are implanted in the uterus of recipient sows. Recently a simplified procedure has been devised, termed handmade cloning (HMC) [25,26], which has been applied in a number of cases including the generation of genetically engineered pigs [27]. Hence, minipigs expressing a mutated APP gene associated with early onset of Alzheimer's disease [28] have been produced by this method. More recently minipigs ectopically expressing human integrins as a model for skin inflammation [29], minipigs expressing a human PCSK9 gain-of-function mutant, as a model for early onset atherosclerosis [30], or minipigs expressing the nematode fat-1 gene for improved nutritional value of pigs [31] have also been produced.
Here we describe HMC-based generation of transgenic Bamaminipigs harboring an antimorphic hCRY1 mutation (Cys414-Ala), which previously has been shown to entail disturbance of the circadian rhythm but also phenotypes reminiscent of type II diabetes, such as hyperglycemia and polydipsia [32]. Twenty-three alive and healthy piglets were delivered and all animals presented peripheral exogenous Cry1 expression. Moreover, altered expression patterns of circadian clock and proinflammatory genes were evident in cultured fibroblasts, suggesting a disturbance of peripheral oscillators.

Ethics statement
Animal experiment procedures were approved by the life ethics and biological safety review committee of BGI-Research.

Plasmids and hCry1 AP vector construction
The CMV promoter in pcDNA3.1 (Invitrogen, Paisley, UK) was released by BglII/NheI and exchanged with the CMV early enhancer and chicken beta-actin (CAG) promoter PCR amplified from pEGFP-N1 (Clontech, CA, USA). The new construct was denoted pCAG-Neo. A multiple cloning site was, subsequently, inserted in EcoRI and NotI relaxed pCAG-Neo. The first intron of the chicken beta-actin gene was inserted into NheI/EcoRI. Human Cry1 (hCry1) cDNA was purchased from FulenGen (ID:M0114, FulenGen, Guangzhou, GD, China). The Cys414-Ala mutation was introduced by site-directed mutagenesis (Quik-Change, Stratagene, Santa Clara, CA, USA) employing 59-CAGTTTTTTCACTGCTATGCCCCTGTTGGTTTTGGT-AGG-39 and 59-CCTACCAAAACCAACAGGGGCATAG-CAGTGAAAAAACTG-39 as forward and reverse PCR primers, respectively. The PCR product was digested with DpnI and inserted into the pMD18-T plasmid (TaKaRa, Shiga, Japan) and subsequently transformed into competent DH5a cells. The purified plasmid was digested with AgeI/NotI to release the hCry1 AP CDS fragment, which was then inserted into EcoRI/ NotI digested pCAG-intron-Neo producing the final construct designated pCAG-intron-hCry1 AP .SV40-Neo.

Handmade cloning and embryo transfer
The specific HMC procedure used for pig cloning has previously been described [25,31]. Briefly, cumulus-oocyte complexes (COCs) collected from porcine ovaries were washed and incubated in 4-well plates for in vitro maturation (IVM) at 38.5uC in 5% CO 2 humidified atmosphere. After 41-42 h cumulus cells were separated from the matured oocytes by hyaluronidase treatment. Following partial digestion of the zonae pellucidae, nuclei were removed by bisection. Each cytoplast was transferred into a fusion chamber containing fusion medium and fused to a single transgenic fibroblast with a single 100 V direct current (DC) impulse of 2.0 kV/cm for 9 ms. After one hour, each cytoplast-fibroblast pair was fused with another cytoplast in activation medium using a DC impulse of 0.8 kV/cm for 80 ms. Hereafter, the fused cells were incubated for 4-6 h in porcine zygote medium 3 (PZM-3) supplemented with 5 mg/ml cytochalasin B and 10 mg/ml cyclohexinmide at 38.5uC, 5%CO 2 , 5% O 2 , and 90% N 2 at maximum humidity. Reconstructed embryos developed into transgenic blastocysts during an additional 6 days of culture in PZM-3 medium and were eventually transferred surgically into the uterine horns of recipient sows [33]. Pregnancies were diagnosed by ultrasonography on day 28 post-transfer and monitored every 2 weeks afterwards [25]. Twenty-three live, transgenic minipigs were eventually obtained.

Identification of hCry1 in the transgenic animals
Genomic DNA was phenol-chloroform extracted from transfected primary fetal fibroblasts or tail clips from the 23 transgenic minipigs and a non-transgenic control (#321-1 wt ) at seven days of age. Transgene identification was carried out in PCR reactions of 20 ng genomic DNA, 20pmol of each primer (cACTB intron forward 59-TTCATACCTCTTATCTTCCTCCCA-39 and CRY1 reverse 59-CTTCCACTGCTGCTACAACCTG-39) and 0.5 unit of rTaq polymerase (Takara). PCR conditions were as follows; 5 min at 95uC; 35 cycles of 30 sec. at 95uC, 30 sec. at 60uC, 40 sec. at 72uC; and final extension at 72uC for 5 min. Finally, 4 mL PCR product was electrophoresed on a 1.0% agarose gel.
For relative measurements of hCry1 AP mRNA levels total RNA was Trizol extracted from transfected primary fetal fibroblasts as well as from tail clips of the 23 transgenic minipigs and the nontransgenic control (#321-1 wt ). Subsequently, 1 mg RNA was in one step DNase treated and reverse transcribed using the RevertAid First Strand cDNA Synthesis Kit (Takara) according to manufacturer's instructions. From the total cDNA, 1 mL was used for quantification using an ABI 7500 Real-Time PCR machine (ABI, CA, USA), Power SYBR Green (Takara) and hCry1 specific primers ( Table 1). Conditions were as follows: 30 sec. at 95uC; 35 cycles of 5 sec. at 95uC, 34 sec at 60uC, and 40 sec. at 72uC. GAPDH was used as internal control ( Table 1).

Transcription analysis of cytokine and circadian clock genes
Epidermal fibroblasts from three transgenic pigs (#175-3 AP , #175-5 AP , and #377-3 AP ) and a non-transgenic pig (#321-1 wt ) were expanded from ear-derived skin biopsies. The established cell cultures were seeded in 12-well plates and maintained until 100% confluency at which the cells were synchronized with 50% FBS in DMEM for 2 hours and successively recovered in serum-free medium for 6 hours. At the given time points the cells were harvested and total RNA Trizol extracted. Quantitative RT-PCR for the gene-transcripts of porcine Per2, Cry1, Clock, Bmal1, Interleukin-6 (IL-6), Tumor-necrosis-factor alpha (TNF-a), and the transgenic human Cry1 ( Table 1) was performed as described above.

Statistical analysis
P-values were calculated by a two-tailed Student's t-test or a one-sided ANOVA where appropriate to test the null hypothesis of no difference between the compared groups. The assumption of equal variances was tested by F-test; if significantly different Welch's correction was applied. The assumption of equal SD was tested by Bartlett-test; if significantly different the non-parametric Dunn's Multiple Comparisons Test was used. Assumption of normality was tested using the Kolmogorov and Smirnov method. In all statistical analyses, p-values ,0.05 were considered significant.

Establishment of transgenic donor cells for HMC
With the objective to selectively express an antimorphic human CRY1 (hCRY1) mutant in cultured cells, we firstly constructed a vector comprising two modules; (i) a CAG-intron-hCry1 wt expression cassette composed of the CAG promoter, the first intron of chicken beta-actin (cACTB) gene and the wt hCRY1 gene, and (ii) a downstream SV40-Neo selection cassette composed of the SV40 promoter and the neomycin (Neo) resistance gene. Hereafter, the hCry1 Cys414-Ala antimorphic mutation (hCry1 AP ) was introduced by site-directed mutagenesis generating the final vector pCAG-intron-hCry1 AP .SV40-Neo ( Figure 1A). Cys414 is part of a well conserved CP dipeptide motif across CRY-proteins and phyla, indicating its importance for CRY1 protein functionality. Based on the solved crystal-structure of murine CRY2, the CP motif appears situated within a co-factor binding pocket regulating CRY1 stability through the interaction with the ubiquitin ligase complex SCF FBXL3 , flavin adenine dinucleotide (FAD) and PERs [34]. Although both the catalytic function with PER proteins and repression of the CLOCK/BMAL1 complex appears unaffected by the Cys414-Ala point mutation [12,32], circadian rhythms and cellular function (in e.g. pancreatic betacells) are disturbed. Thus, suggesting that Cry1 AP affects kinetic and dynamic properties or influences non-canonical circadian pathways.
In order to generate hCRY1 AP expressing donor cells, we stably transfected Bama porcine fetal fibroblasts (BM-PFFs) with pCAGintron-hCry1 AP .SV40-Neo and applied a selection pressure with G-418 for 10 days. Subsequently, 25 individual G-418 resistant colonies were isolated and expanded. Genomic DNA extracted from the 25 clones was subjected to PCR analysis with Cry1-Neo specific primers revealing 12 positive clones (data not shown). Subsequently, RT-qPCR on total mRNA extracted from the 12 positive clones specific for hCry1 transcripts demonstrated that hCry1 was expressed in all clones at varying intensities ( Figure 1B). Thus, relative to hCry1 AP clone # 35, the three clones designated # 6, # 12 and # 15 exhibited the highest expression levels, with hCry1 mRNA levels 17, 10 and 8 fold above that of clone # 35, respectively. These high-expressing clones were chosen for HMC. Importantly, the hCry1 transcript was undetectable in un-transfected BM-PFFs.

Generation of transgenic cloned pigs by HMC
The selected hCry1 AP clones (# 6, # 12 and # 15) were individually applied as donor cells for HMC. The blastocyst rate was 45.5%, and a total of 691 reconstituted embryos were surgically transferred to the uteri of six naturally cycling recipient sows, with 90-135 reconstituted embryos transferred per recipient ( Table 2). Of the six recipient sows, five became pregnant from which four went to term delivering a total of 32 piglets of which six were stillborn and two were mummified. Of the living 24 piglets one presented with malformations and was euthanized. Further two succumbed within the first months. The surviving 21 piglets remain healthy and gain weight according to the growth curve of the herd (Figure 2A).

Identification of hCry1 AP expression in the twenty-three transgenic minipigs
All cloned minipigs were evaluated for transgene integration and expression. Hence, transgenic status was demonstrated in all 23 cloned animals by PCR amplification of hCry1 on tail clip extracted genomic DNA obtained from the cloned minipigs at the age of one week. A PCR fragment with the expected size of 746 bp appeared in the lanes for each of the cloned minipigs but not for an age-matched non-transgenic (NT) control ( Figure 2B, upper). Assessment of transgene expression was achieved by conducting RT-PCR on total RNA extracted from the tail clips.
Again an easily detectable band corresponding to hCry1 cDNA was observed for all transgenic minipigs but not for the NT control ( Figure 2B, lower). Analysis of the expression level of Cry1 showed that the relative expression level in the transgenic minipigs varied up to 7-fold compared to that in the lowest expressing transgenic animal (#351-4 AP ) ( Figure 2C). A general higher expression level was detected in animals originating from donor cell # 6 (sow #377 and #351) suggesting a correlation between the expressions observed in the donor cells and in the transgenic minipigs. Interestingly, the expression appears to be more heterogeneous in piglets arising from donor cell # 6 compared to the expression pattern among piglets originating from donor cell # 12 (sow #175 and #208). This suggests that non-inherited positional effects work on the transgene -especially in animals originating from donor cell # 6. Pig #377-1 AP and pig #208-2 AP were the animals displaying the highest relative hCry1 AP expression. Notably, they perished at the age of about 35 days, hinting a possible concentration dependent toxicity.

Altered oscillation of circadian gene expression in the transgenic minipigs
Clock-gene expression levels inherently oscillate in well defined phases with time in intact peripheral tissue. In cell culture this feature is lost but can, however, be reinitiated and synchronized through serum-shock. Hence, in ex vivo fibroblasts this leads to a cell-autonomous oscillation in clock gene expression for at least 48 hours [35,36]. Furthermore, functional CRY1 is required for a sustained circadian rhythm in dissociated murine fibroblasts [37]. In order to evaluate if the introduced point mutation in hCry1 would exert disruption to the cell-autonomous oscillation, we transfected primary porcine fibroblasts with a construct expressing the functional hCRY1 protein (pCAG-intron-hCry1.SV40-Neo) or the antimorphic variant (pCAG-intron-hCry1 AP .SV40-Neo). Following antibiotic selection and RT-qPCR analysis, clones expressing the functional (clone # 3) or antimorphic hCry1 (clone # 3) at comparable intensity and displaying similar growth rates were selected and further expanded ( Figure 3A). At confluency the cells were serum shocked for 2 h after which they were returned to serum-free medium for an additional 6 h. Subsequently, cells were harvested for RT-qPCR analysis every 4 hours through 48 hours.
The mRNA level of pPer2 in the hCry1 transfected cells appears to have a first zenith at Zeitgeber time (ZT) 16 h and a second at ZT 36 h ( Figure 3B). On the other hand, the zeniths of pCry1 mRNA expression are slightly offset occurring at ZT 20 h and 44 h ( Figure 3C). This slight offset but synchronous oscillation of Per2 and Cry1 has previously been documented in SCN as well as peripheral tissue [38]. Interestingly, the rhythmic pattern is only partially maintained in the hCry1 AP transfected cells and, moreover, displays oscillations of greater magnitude, suggesting that the regulation of pPer2 and pCry1 is unbalanced in the presence of the antimorphic hCRY1 AP .
With the intent to examine circadian gene expression patterns in NT as compared to transgenic minipigs we cultured epidermal fibroblasts from transgenic minipigs (#175-3 AP , #175-5 AP and #377-3 AP ) and from a NT control (#321-1 wt ). At confluency the cells were serum shocked for 2 h and then allowed to recover in serum-free medium for another 6 h. Afterwards the cells were harvested for RT-qPCR analysis every 4 hours through 52 hours. Firstly, we validated hCry1 AP mRNA expression in all the three cell populations. After an initial increase in expression level in pig #175-5 AP , all mRNA levels (relative to pig #175-3 AP ) remained within a stable range with only minor yet rhythmic fluctuations. Importantly, no call above background was observed in the control ( Figure 4A). Secondly, we examined porcine Cry1, Per2, Clock and Bmal1 mRNA levels in the cultured cells following synchronization. Given the nature of the feedback regulatory effect of CRY1/PER2 on CLOCK/BMAL1 the corresponding two sets of genes are expected to have peak expression levels occurring in  Figure S1). Notably, the oscillations of Per2 and Cry1 appeared slightly displaced and with a circadian rhythm similar to the one observed in the transfected porcine fibroblasts. Conversely pBaml1 transcript levels displayed a first nadir around ZT 12 and 36, thus, being inversely correlated to pPer2 and, to a lesser extent, pCry1 expression. On the other hand, pClock mRNA levels presented a much weaker circadian pattern although a weak zenith is detectable a ZT 20 and 44. However, this seems to be consistent with previous findings in murine peripheral tissues [39]. In the transgenic animals the expression curves of hCry1, pPer2, pCry1, pClock, and pBmal1 displayed key differences and similarities compared to the NT control (Figure 4 A-E). Oscillation in the mRNA levels of hCry1 follows that of pCry1, indicating that hCry1 is susceptible to the same post-transcriptional regulation otherwise governing the levels of pCry1. Supporting this notion, is the observation that the 39 untranscribed region (UTR) of mouse Cry1 mRNA contains a destabilizing element to which heterogeneous nuclear ribonucleoprotein D (hnRNP D) binds facilitating mRNA degradation [40]. As in the NT control, expression levels of pPer2, pClock and pBaml1 in the transgenic animals follow a rhythmic pattern. However, in the latter the amplitude is increased and the oscillation period decreased. Hence, the ZT for the first zenith is comparable between the control and the transgenic animals. However, cells from the latter seem to complete two cycles within 24 h compared to a single cycle in the NT control.
Group means and analysis of variance between group means at each time point is summarized in Table S1 and S2.
Looking specifically at the acrophases of pCry1 and pPer2 within the first 24 hours, there is a high concordance between the cells derived from the transgenic and NT animals ( Figure S2). Hence for pCry1 the acrophase is at 18 4 h) for pPer2. Taken together, the redundancy in the TTFLs seemingly ensures rhythmicity in the transgenic fibroblasts but with an altered circadian period (t) and amplitude. A slight intra-individual variation in t seems evident but falls in line with previous reports on circadian rhythms in human fibroblasts [41].

Increased activation of proinflammatory markers in the hCry1 AP transgenic minipigs
Next we wanted to determine whether the altered circadian rhythms and/or the expression of the antimorphic hCry1 in the transgenic animals would lead to activation of proinflammatory cytokines. The level of IL-6 and TNF-a mRNA was quantified in cDNA originating from serum shocked hCry1 AP transgenic fibroblasts (derived from pig #377-1 AP , #208-4 AP , and #175-5 AP ) and NT fibroblasts (derived from pig #321-1 wt ). An initial analysis revealed a sharp and linear increase in the mRNA level of IL-6 in the transgenic animals reaching a 7-10 fold increment at ZT 14 compared to ZT 0. Notably, this observation was absent in the NT control ( Figure 5A). To confirm these findings we repeated the experiment with epidermal fibroblasts derived from pig #175-5 AP and a NT control (#321-1 wt ) with measurements every fourth hour through 28 hours ( Figure 5B). The more detailed time course confirmed the initial steep increase in IL-6 mRNA levels until ZT 12 after which the level dropped off reaching a nadir at ZT 20. Comparing the expression levels of IL-6 at nadir and zenith, minipig #175-5 AP presents a level about 3fold higher relative to the NT control. Notably, however, the IL-6 mRNA level in both the transgenic and NT animal revolted with a t of approximately 12 hours. In a parallel experiment assessing TNF-a mRNA levels in the same pool of cells derived from pig #175-5 AP and the NT control showed significant increase in TNFa expression in the transgenic minipig, that is, reaching levels 3-9 fold higher relative to the levels observed in the control ( Figure 5C). This lends further support to the notion of proinflammatory induction in the transgenic cells. Interestingly, a complete circadian cycle was completed in 16 hours in cells from the NT control animal whereas in only took 8 hours for a complete cycle in the cells derived from the transgenic animal.

Circadian oscillation of body temperature
In mammals, body temperature (T b ) shows circadian oscillation. Therefore, recording T b fluctuation is considered an effective way of determining circadian rhythm perpetuations [42]. Body temperature was recorded every second hour through 4 days in five transgenic minipigs (#175-5 AP , #208-4 AP , #351-3 AP , #351-6 AP , and #377-2 AP ) and three NT minipigs (#321-1 wt , #321-3 wt , and #C-28 wt ) ( Figure 6A). Animals in both groups displayed normal T b curves with a steady increase during daylight and a steady decrease during nighttime. Three of the five hCry1 AP and two of the three NT pigs showed rhythmic T b oscillations with t  Table S3). Among the remaining three animals, especially the transgenic minipig #208-4 AP exhibited an abnormal T b curve ( Figure 6B). In summary, these results may imply that expression of the antimorphic mutant hCry1 AP potentially leads to disturbance of the circadian rhythm of T b in the transgenic minipigs.

Discussion
Minipigs have become an attractive mammalian model for research in human ailments including metabolic, cardio-vascular and degenerative disorders. Small body sizes as well as physiological and anatomical resemblance to humans has raised the interest in minipigs and so basic as well as clinical studies have already been conducted. Furthermore, the complete sequencing of the pig genome has further spawned a general interest in genetically modified porcine disease-models [43]. In this report we describe the generation of Bama-minipigs holding a genomically integrated transcription cassette from which a heterologous and antimorphic CRY1 is expressed. The transgene is under the control of a CAG promoter which has shown almost ubiquitous activity in all tissues including the brain and skin [44,45]. Expression of the transgene Figure 3. Expression of the antimorphic but not the wt hCry1 instigate altered expression profiles of pPer2 and pCry1. Primary porcine fibroblasts were stably transfected with pCAG-intron-hCry1.SV40-Neo or pCAG-intron-hCry1 AP .SV40-Neo and relative expression were assessed by RT-qPCR with exon-exon spanning hCry1 primers. (B-C) Serum-shocked cells were harvested every fourth hour through 48 hours. Total RNA was extracted and used for RT-qPCR with exon-exon spanning primers targeting pPer2 or pCry1. GAPDH normalized data relative to hCry1 clone # 3 is depicted as a function of time. The experiment is performed in triplicate and data are presented as mean values with smoothened curves. Grey dashed lines indicate the first and second zenith of pPer2 and pCry1 mRNA expression in hCry1 containing cells. doi:10.1371/journal.pone.0076098.g003 Figure 4. Altered oscillation patterns of key clock transcripts in fibroblasts from transgenic minipigs. Fibroblasts were expanded from skin-biopsies obtained from three transgenic minipigs (#175-3 AP , #175-5 AP , and #377-3 AP ) and one NT control (#321-1 wt ). The cells were serumshocked after which total RNA was extracted every fourth hour through 52 hours. Quantitative RT-PCR was performed with exon-exon spanning primers targeting (A) hCry1 (B) pPer2 (C) pCry1 (D) pClock and (E) pBmal1 normalized to endogenous ACTB. Expression is relative to minipig #175-3 AP in (A) and to the NT control in (B-E) and is depicted as a function of time. The experiment is performed in triplicate and data are presented as mean values with smoothened curves. Grey dashed lines indicate the first, second and third zenith of pPer2 mRNA expression in cells from minipig #175-5 AP . doi:10.1371/journal.pone.0076098.g004 in the 23 live-born, seemingly healthy, cloned animals was verified by PCR and RT-qPCR on DNA and RNA, respectively, extracted from skin biopsies. Curiously, despite originating from only two different donor-cell populations the expression pattern in the transgenic minipigs appeared rather heterogeneous. As the cycle of the external cues (light and feed) are highly controlled it is unlikely that the circadian clocks among the pigs are out of synchronization. However, stochastic events, genomic rearrangements or epigenetic changes during embryogenesis could explain the variation in expression.
Next we wanted to assess to what extent the introduced antimorphic hCry1 AP would be expressed in cultured fibroblasts from the pigs and whether this would alter the autonomous peripheral oscillations of other key clock components. It is well known that several clock effectors including Per2 and Cry1 contain 3'UTRs prone for mRNA degradation [40,46]. Interestingly, the mRNA level of hCry1 AP itself oscillated with a periodicity comparable to pCry1, suggesting that the transgenic Cry1 is subject to the same post-transcriptional regulation as the endogenous Cry1. Furthermore, there appears to be consistency between expression levels of hCry1 AP and the clock genes, hence, minipig #175-5 AP presented the highest relative expression in all assays.
Whereas peripheral oscillations of core clock gene transcripts demonstrated an approximately 24 hours periodicity in epidermal fibroblasts from the NT control pig, the same gene transcripts oscillated with a frequency of about 12 hours in cells from the three transgenic animals tested (#175-3 AP , #175-5 AP , and #377-3 AP ). Importantly, previous reports demonstrated coherence between peripheral oscillations in vitro and in vivo. Furthermore, deviations in the circadian rhythm of clock gene expression were found to be retained in vitro [47,48]. These results, thus, appear to be in line with the abnormal circadian behavior observed in mouse Cry1 AP transgenic mice [32]. This indicates that the regulatory function of CRY1 cannot be rescued by CRY2 or a second regulatory loop, involving the orphan nuclear receptors Rev-Erba and RORa.
Oscillations in the mRNA levels of the clock driving genes pClock and pBaml1 are quantitatively small compared to those of the regulatory arm, pCry1 and pPer2. Importantly, however, there is an anti-phase correlation between the former and pPer2, suggesting that regulation of the cell-autonomous circadian clock is maintained but with a t of only 12 hours. Notably, oscillation of Figure 6. Discrete circadian oscillation of body temperature in hCry1 AP transgenic compared to NT control minipigs. (A) Body temperature was measured using an infrared thermometer in five transgenic minipigs (#175-5 AP , #377-2 AP , #351-3 AP , #351-6 AP , and #208-4 AP ) and three NT control animals (#321-1 wt , #321-2 wt , and #321-3 wt ) every second hour over a time course of 4 days. Mean body temperature is plotted as a function of time 6 standard deviation. (B) Body temperature fluctuations of minipig #208-4 AP compared to the mean from the three NT control minipigs (#321-1 wt , #321-2 wt , and #321-3 wt ) obtained as described above. doi:10.1371/journal.pone.0076098.g006 pPer2 is both very consistent and robust across all transgenic fibroblast cultures. This supports not only the notion of pPer2 being the primary intrinsic Zeitgeber of rhythmicity in peripheral tissue [49,50] but also substantiating the t alternating effect of hCRY1 AP . Furthermore, the expression profiles of Cry1 and pPer2 are synchronized but slightly staggered. This lag in time has previously been described in both the SCN and peripheral tissue [39,51], supporting the concept of PER2 induced stabilization of CRY1 by prevention of ubiquitylation and proteosomal degradation [34,52]. Together with the likewise adopted model of CRY1 induced stabilization of PER2 through similar mechanisms [52,53], we propose a working model in which the hCRY1 AP protein enjoys increased stability leading to augmented pPER2 stability which downstream leads to increased stability of endogenous CRY1. Further supporting this view is the marked difference in expression amplitude between hCry1 and hCry1 AP transfected porcine fibroblasts.
Sleep-loss studies have identified physiological connection between the circadian clock and the immune system [54]. As a chronic state of metainflammation is known to be a substantial risk factor in metabolic disorders and cancer [55,56,57] and because the absence of CRY in knockout mice has been shown to entail increased expression levels of proinflammatory cytokines including TNF-a and IL-6 [58], we sought to determine if the antimorphic variant potentially could lead to a higher steady-state of inflammation in the transgenic animals. A considerable difference was apparent with average levels of IL-6, 4-fold (p = 0.047), and TNF-a, 5-fold (p = 0.0025), higher in fibroblasts from minipig #175-5 AP relative to the control. Notably, the oscillation of IL-6 in pig #175-5 AP displayed a periodicity of about 12 hours which is coherent with the expression profile of hCry1 AP . However, IL-6 in the NT control followed the same pattern although pCry1 had a t of about 24 hours.
CRY proteins have shown to bind adenylyl cyclase (ADCY) leading to a decrease in cyclic AMP production thereby preventing the activation of NF-kB and downstream expression of IL-6 -with the opposite effect taken place in the absence of CRY protein. Our data suggests that hCRY1 AP competes with endogenous CRY proteins for binding to ADCY but that it lacks repressor activity.
Body temperature (T b ) control is a very important projection emanating from the SCN as it plays a role in the resetting of SCN itself but also as an auxiliary Zeitgeber of peripheral oscillators [51,59]. Core T b displays a clear circadian rhythm in mammals with fluctuations of 1-4uC [60]. We therefore wanted to assess whether the presence of hCRY1 AP protein would interfere with T b control in the transgenic animals as previously observed in Cry1 2/ 2 /Cry2 2/2 double knockout mice [61]. The average t and amplitude of the T b rhythm was unchanged (t<23 hours, p#0.005) in the transgenic pigs compared to the controls under normal light-dark and feeding conditions. However, there were a few outliers and especially pig #208-4 AP presented an arrhythmic T b oscillation, suggesting, that the central circadian clock, in at least this animal, is disturbed even under normal entrainment. It has previously been implied that CRY1 is dispensable for persistent SCN rhythmicity [37]. However, it is conceivable that the antimorphic hCRY1 mutant interferes with the circadian clock in a non-canonical fashion.
In summary we have generated 23 Bama-minipigs transgenic for a tissue unspecific antimorphic human Cry1. We have shown stable expression of the transgene in the skin of all cloned minipigs and demonstrated that this gives rise to altered circadian patterns of key peripheral oscillators as well as amplified induction of proinflammatory cytokines. Further analyses are currently undertaken and we anticipate that such porcine models will contribute to a better understanding of the link between the pathogenesis of circadian and metabolic disorders. Figure S1 Circadian rhythmicity in expression patterns of key circadian regulatory genes. Fibroblasts expanded from skin-biopsies obtained from a non-transgenic minipig (#321-1) were serum-shocked after which total RNA was extracted every fourth hour through 52 hours. Quantitative RT-PCR was performed with exon-exon spanning primers targeting hCRY1, pPer2, pCry1, pClock, and pBmal1 normalized to endogenous ACTB. Expression is relative to the first measurement time point (ZT 0) and depicted as a function of time. The experiment is performed in triplicate and data are presented as means. Grey