Long-read sequencing for fast and robust identification of correct genome-edited alleles: PCR-based and Cas9 capture methods.

BACKGROUND: Recent developments in CRISPR/Cas9 genome-editing tools have facilitated the introduction of precise alleles, including genetic intervals spanning several kilobases, directly into the embryo. However, the introduction of donor templates, via homology directed repair, can be erroneous or incomplete and these techniques often produce mosaic founder animals. Thus, newly generated alleles must be verified at the sequence level across the targeted locus. Screening for the presence of the desired mutant allele using traditional sequencing methods can be challenging due to the size of the interval to be sequenced, together with the mosaic nature of founders. METHODOLOGY/PRINCIPAL FINDINGS: In order to help disentangle the genetic complexity of these animals, we tested the application of Oxford Nanopore Technologies long-read sequencing at the targeted locus and found that the achievable depth of sequencing is sufficient to offset the sequencing error rate associated with the technology used to validate targeted regions of interest. We have assembled an analysis workflow that facilitates interrogating the entire length of a targeted segment in a single read, to confirm that the intended mutant sequence is present in both heterozygous animals and mosaic founders. We used this workflow to compare the output of PCR-based and Cas9 capture-based targeted sequencing for validation of edited alleles. CONCLUSION: Targeted long-read sequencing supports in-depth characterisation of all experimental models that aim to produce knock-in or conditional alleles, including those that contain a mix of genome-edited alleles. PCR- or Cas9 capture-based modalities bring different advantages to the analysis.

Characterisation and prion transmission study in mice with genetic reduction of sporadic Creutzfeldt-Jakob disease risk gene Stx6.

Sporadic Creutzfeldt-Jakob disease (sCJD), the most common human prion disease, is thought to occur when the cellular prion protein (PrPC) spontaneously misfolds and assembles into prion fibrils, culminating in fatal neurodegeneration. In a genome-wide association study of sCJD, we recently identified risk variants in and around the gene STX6, with evidence to suggest a causal increase of STX6 expression in disease-relevant brain regions. STX6 encodes syntaxin-6, a SNARE protein primarily involved in early endosome to trans-Golgi network retrograde transport. Here we developed and characterised a mouse model with genetic depletion of Stx6 and investigated a causal role of Stx6 expression in mouse prion disease through a classical prion transmission study, assessing the impact of homozygous and heterozygous syntaxin-6 knockout on disease incubation periods and prion-related neuropathology. Following inoculation with RML prions, incubation periods in Stx6-/- and Stx6+/- mice differed by 12 days relative to wildtype. Similarly, in Stx6-/- mice, disease incubation periods following inoculation with ME7 prions also differed by 12 days. Histopathological analysis revealed a modest increase in astrogliosis in ME7-inoculated Stx6-/- animals and a variable effect of Stx6 expression on microglia activation, however no differences in neuronal loss, spongiform change or PrP deposition were observed at endpoint. Importantly, Stx6-/- mice are viable and fertile with no gross impairments on a range of neurological, biochemical, histological and skeletal structure tests. Our results provide some support for a pathological role of Stx6 expression in prion disease, which warrants further investigation in the context of prion disease but also other neurodegenerative diseases considering syntaxin-6 appears to have pleiotropic risk effects in progressive supranuclear palsy and Alzheimer's disease.

Application of long single-stranded DNA donors in genome editing: generation and validation of mouse mutants.

BACKGROUND: Recent advances in clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) genome editing have led to the use of long single-stranded DNA (lssDNA) molecules for generating conditional mutations. However, there is still limited available data on the efficiency and reliability of this method. RESULTS: We generated conditional mouse alleles using lssDNA donor templates and performed extensive characterization of the resulting mutations. We observed that the use of lssDNA molecules as donors efficiently yielded founders bearing the conditional allele, with seven out of nine projects giving rise to modified alleles. However, rearranged alleles including nucleotide changes, indels, local rearrangements and additional integrations were also frequently generated by this method. Specifically, we found that alleles containing unexpected point mutations were found in three of the nine projects analyzed. Alleles originating from illegitimate repairs or partial integration of the donor were detected in eight projects. Furthermore, additional integrations of donor molecules were identified in four out of the seven projects analyzed by copy counting. This highlighted the requirement for a thorough allele validation by polymerase chain reaction, sequencing and copy counting of the mice generated through this method. We also demonstrated the feasibility of using lssDNA donors to generate thus far problematic point mutations distant from active CRISPR cutting sites by targeting two distinct genes (Gckr and Rims1). We propose a strategy to perform extensive quality control and validation of both types of mouse models generated using lssDNA donors. CONCLUSION: lssDNA donors reproducibly generate conditional alleles and can be used to introduce point mutations away from CRISPR/Cas9 cutting sites in mice. However, our work demonstrates that thorough quality control of new models is essential prior to reliably experimenting with mice generated by this method. These advances in genome editing techniques shift the challenge of mutagenesis from generation to the validation of new mutant models.

Synthesis of poly(meth)acrylates with thioether and tertiary sulfonium groups by ARGET ATRP and their use as siRNA delivery agents.

The field of RNA interference depends on the development of safe and efficient carriers for short interfering ribonucleic acid (siRNA) delivery. Conventional cationic monomers for siRNA delivery have utilized the nitrogen heteroatom to produce cationic charges. Here, we polymerized cationic sulfonium (meth)acrylate by activators regenerated by electron transfer (ARGET) atom transfer radical polymerization (ATRP) to form polymers with narrow molecular weight distributions for siRNA delivery. The tertiary sulfonium species was stable toward dealkylation in water but less stable in the polar aprotic solvent dimethyl sulfoxide. Block copolymers poly(ethylene oxide) with poly(meth)acrylate containing sulfonium moieties were prepared as an siRNA delivery platform. Results suggested block copolymers were biocompatible up to 50 µg/mL in vitro and formed polyplexes with siRNA. Additionally, block copolymers protected siRNAs against endonuclease digestion and facilitated knockdown of glyceraldehyde 3-phosphate dehydrogenase (Gapdh) mRNA expression in murine calvarial preosteoblasts. The versatility, biocompatibility, and cationic nature of these tertiary sulfonium groups are expected to find widespread biological applications.

Molecular brakes regulating mTORC1 activation in skeletal muscle following synergist ablation.

The goal of the current work was to profile positive (mTORC1 activation, autocrine/paracrine growth factors) and negative [AMPK, unfolded protein response (UPR)] pathways that might regulate overload-induced mTORC1 (mTOR complex 1) activation with the hypothesis that a number of negative regulators of mTORC1 will be engaged during a supraphysiological model of hypertrophy. To achieve this, mTORC1-IRS-1/2 signaling, BiP/CHOP/IRE1α, and AMPK activation were determined in rat plantaris muscle following synergist ablation (SA). SA resulted in significant increases in muscle mass of ~4% per day throughout the 21 days of the experiment. The expression of the insulin-like growth factors (IGF) were high throughout the 21st day of overload. However, IGF signaling was limited, since IRS-1 and -2 were undetectable in the overloaded muscle from day 3 to day 9. The decreases in IRS-1/2 protein were paralleled by increases in GRB10 Ser(501/503) and S6K1 Thr(389) phosphorylation, two mTORC1 targets that can destabilize IRS proteins. PKB Ser(473) phosphorylation was higher from 3-6 days, and this was associated with increased TSC2 Thr(939) phosphorylation. The phosphorylation of TSC2 (Thr1345) (an AMPK site) was also elevated, whereas phosphorylation at the other PKB site, Thr(1462), was unchanged at 6 days. In agreement with the phosphorylation of Thr(1345), SA led to activation of AMPKα1 during the initial growth phase, lasting the first 9 days before returning to baseline by day 12. The UPR markers CHOP and BiP were elevated over the first 12 days following ablation, whereas IRE1α levels decreased. These data suggest that during supraphysiological muscle loading at least three potential molecular brakes engage to downregulate mTORC1.

Protein kinase D isoforms are dispensable for integrin-mediated lymphocyte adhesion and homing to lymphoid tissues.

Leukocyte function-associated antigen-1 (LFA-1) and very late antigen-4 (VLA-4) integrins are essential for lymphocyte adhesion, trafficking and effector functions. Protein kinase D (PKD) has previously been implicated in lymphocyte integrin regulation through regulation of Rap1 activity. However, the true role of PKD in integrin regulation in primary lymphocytes has not previously been investigated. The major PKD isoform in lymphocytes is PKD2. Here we employed PKD2-deficient mice, a specific PKD kinase inhibitor, as well as PKD-null DT40 B cells to investigate the role of PKD in integrin regulation in lymphocytes. We report that PKD2-deficient lymphocytes bound normally to integrin ligands in static and shear flow adhesion assays. They also homed normally to lymphoid organs after adoptive transfer into wild-type mice. DT40 B cells devoid of any PKD isoforms and primary lymphocytes pretreated with a specific PKD inhibitor bound normally to integrin ligands, indicating that multiple PKD isoforms do not redundantly regulate lymphocyte integrins. In addition, PKD2-deficient lymphocytes, as well as DT40 cells devoid of any PKD isoforms, could activate Rap1 in response to B-cell receptor ligation or phorbol ester treatment. Together, these results show that the PKD family does not play a critical role in lymphocyte integrin-mediated cell adhesion or lymphocyte trafficking in vivo.

Genotyping Genome-Edited Founders and Subsequent Generation.

Targeted nucleases allow the production of many types of genetic mutations directly in the early embryo. However, the outcome of their activity is a repair event of unpredictable nature, and the founder animals that are produced are generally of a mosaic nature. Here, we present the molecular assays and genotyping strategies that will support the screening of the first generation for potential founders and the validation of positive animals in the subsequent generation, according to the type of mutation generated.

Broad-spectrum XX and XY gonadal dysgenesis in patients with a homozygous L193S variant in PPP2R3C.

CONTEXT: Homozygous and heterozygous variants in PPP2R3C are associated with syndromic 46,XY complete gonadal dysgenesis (Myo-Ectodermo-Gonadal Dysgenesis (MEGD) syndrome), and impaired spermatogenesis, respectively. This study expands the role of PPP2R3C in the aetiology of gonadal dysgenesis (GD). METHOD: We sequenced the PPP2R3C gene in four new patients from three unrelated families. The clinical, laboratory, and molecular characteristics were investigated. We have also determined the requirement for Ppp2r3c in mice (C57BL6/N) using CRISPR/Cas9 genome editing. RESULTS: A homozygous c.578T>C (p.L193S) PPP2R3C variant was identified in one 46,XX girl with primary gonadal insufficiency, two girls with 46,XY complete GD, and one undervirilised boy with 46,XY partial GD. The patients with complete GD had low gonadal and adrenal androgens, low anti-Müllerian hormone, and high follicle-stimulating hormone and luteinizing hormone concentrations. All patients manifested characteristic features of MEGD syndrome. Heterozygous Ppp2r3c knockout mice appeared overtly normal and fertile. Inspection of homozygous embryos at 14.5, 9.5, and 8.5 days post coitum(dpc) revealed evidence of dead embryos. We conclude that loss of function of Ppp2r3c is not compatible with viability in mice and results in embryonic death from 7.5 dpc or earlier. CONCLUSION: Our data indicate the essential roles for PPP2R3C in mouse and human development. Germline homozygous variants in human PPP2R3C are associated with distinctive syndromic GD of varying severity in both 46,XY and 46,XX individuals.

Prolonged activation of S6K1 does not suppress IRS or PI-3 kinase signaling during muscle cell differentiation.

BACKGROUND: Myogenesis in C2C12 cells requires the activation of the PI3K/mTOR signaling pathways. Since mTOR signaling can feedback through S6K1 to inhibit the activation of PI3K, the aim of this work was to assess whether feedback from S6K1 played a role in myogenesis and determine whether siRNA mediated knockdown of S6K1 would lead to an increased rate of myotube formation. RESULTS: S6K1 activity increased in a linear fashion following plating and was more than 3-fold higher after Day 3 of differentiation (subconfluent = 11.09 +/- 3.05, Day 3 = 29.34 +/- 3.58). IRS-1 levels tended to increase upon serum withdrawal but decreased approximately 2-fold (subconfluent = 0.88 +/- 0.10, Day 3 = 0.42 +/- 0.06) 3 days following differentiation whereas IRS-2 protein remained stable. IRS-1 associated p85 was significantly reduced upon serum withdrawal (subconfluent = 0.86 +/- 0.07, Day 0 = 0.31 +/- 0.05), remaining low through day 1. IRS-2 associated p85 decreased following serum withdrawal (subconfluent = 0.96 +/- 0.05, Day 1 = 0.56 +/- 0.08) and remained suppressed up to Day 3 following differentiation (0.56 +/- 0.05). Phospho-tyrosine associated p85 increased significantly from subconfluent to Day 0 and remained elevated throughout differentiation. siRNA directed against S6K1 and S6K2 did not result in changes in IRS-1 levels after either 48 or 96 hrs. Furthermore, neither 48 nor 96 hrs of S6K1 knockdown caused a change in myotube formation. CONCLUSIONS: Even though S6K1 activity increases throughout muscle cell differentiation and IRS-1 levels decrease over this period, siRNA suggests that S6K1 is not mediating the decrease in IRS-1. The decrease in IRS-1/2 associated p85 together with the increase in phospho-tyrosine associated p85 suggests that PI3K associates primarily with scaffolds other than IRS-1/2 during muscle cell differentiation.

mVps34 is activated following high-resistance contractions.

Following resistance exercise in the fasted state, both protein synthesis and degradation in skeletal muscle are increased. The addition of essential amino acids potentiates the synthetic response suggesting that an amino acid sensor, which is involved in both synthesis and degradation, may be activated by resistance exercise. One such candidate protein is the class 3 phosphatidylinositol 3OH-kinase (PI3K) Vps34. To determine whether mammalian Vps34 (mVps34) is modulated by high-resistance contractions, mVps34 and S6K1 (an index of mTORC1) activity were measured in the distal hindlimb muscles of rats 0.5, 3, 6 and 18 h after acute unilateral high-resistance contractions with the contralateral muscles serving as a control. In the lengthening tibialis anterior (TA) muscle, S6K1 (0.5 h = 366.3 +/- 112.08%, 3 h = 124.7 +/- 15.96% and 6 h = 129.2 +/- 0%) and mVps34 (3 h = 68.8 +/- 15.1% and 6 h = 36.0 +/- 8.79%) activity both increased, whereas in the shortening soleus and plantaris (PLN) muscles the increase was significantly lower (PLN S6K1 0.5 h = 33.1 +/- 2.29% and 3 h = 47.0 +/- 6.65%; mVps34 3 h = 24.5 +/- 7.92%). HPLC analysis of the TA demonstrated a 25% increase in intramuscular leucine concentration in rats 1.5 h after exercise. A similar level of leucine added to C2C12 cells in vitro increased mVps34 activity 3.2-fold. These data suggest that, following high-resistance contractions, mVps34 activity is stimulated by an influx of essential amino acids such as leucine and this may prolong mTORC1 signalling and contribute to muscle hypertrophy.

Pre-transplant immune factors may be associated with BK polyomavirus reactivation in kidney transplant recipients.

BK polyomavirus (BKPyV) reactivation in kidney transplant recipients can lead to allograft damage and loss. The elements of the adaptive immune system that are permissive of reactivation and responsible for viral control remain incompletely described. We performed a prospective study evaluating BKPyV-specific T-cell response, humoral response and overall T-cell phenotype beginning pre-transplant through one year post-transplant in 28 patients at two centers. We performed an exploratory analysis of risk factors for the development of viremia and viruria as well as compared the immune response to BKPyV in these groups and those who remained BK negative. 6 patients developed viruria and 3 developed viremia. BKPyV-specific CD8+ T-cells increased post-transplant in viremic and viruric but not BK negative patients. BKPyV-specific CD4+ T-cells increased in viremic, but not viruric or BK negative patients. Anti-BKPyV IgG antibodies increased in viruric and viremic patients but remained unchanged in BK negative patients. Viremic patients had a greater proportion of CD8+ effector cells pre-transplant and at 12 months post-transplant. Viremic patients had fewer CD4+ effector memory cells at 3 months post-transplant. Exploratory analysis demonstrated lower CD4 and higher total CD8 proportions, higher anti-BKPyV antibody titers and the cause of renal failure were associated BKPyV reactivation. In conclusion, low CD4, high CD8 and increased effector CD8 cells were found pre-transplant in patients who became viremic, a phenotype associated with immune senescence. This pre-transplant T-cell senescence phenotype could potentially be used to identify patients at increased risk of BKPyV reactivation.

A genomic approach to reveal novel genes associated with myotube formation in the model teleost, Takifugu rubripes.

Little is known about the transcriptional networks that regulate myotube production in vertebrates. In the present study, we have used a genomic approach to discover novel genes associated with myotube formation in fast muscle of the tiger puffer fish, Takifugu rubripes. The number of fast muscle fibers per myotome increased until 1.2 kg body mass, and subsequent growth was by fiber hypertrophy alone. Forward and reverse subtracted cDNA libraries were prepared from a 180-g (myotube +) and a 3.4-kg (myotube -) fish, and 1,452 expressed sequence tags (ESTs) were obtained. After these ESTs were grouped into nonredundant clusters and housekeeping and structural genes were eliminated, 57 genes were selected and quantitative PCR was used to investigate their expression levels in different tissues from independent groups of myotube(-) and myotube(+) fish acclimated to the same environmental conditions and diet. Eleven novel genes were found to be consistently differentially expressed, but only four showed appropriate tissue-specific expression. These four genes were upregulated 5-25 times in fast muscle of myotube(-) relative to myotube(+) growth stages, while their expression remained unchanged in the other tissues studied. The novel genes identified, which are also present in other vertebrate genomes, may play a role in inhibiting myotube formation in vertebrate muscle.

BK polyomavirus reactivation after reduced-intensity double umbilical cord blood cell transplantation.

Serial serum samples from 27 patients who underwent double umbilical cord blood transplantation (dUCBT) were analyzed for BK polyomavirus (BKPyV) DNA by real-time PCR and BKPyV-specific immune globulin by ELISA. Clinical data were collected on all patients. All pre-transplant sera had detectable anti-BKPyV IgG. Fifteen patients (56%) had detectable serum BKPyV DNA (median 8.9 × 10(4) copies/ml; range 4.1 × 10(3)-7.9 × 10(6) copies/ml) a median of 40 days (range, 27-733 days) after dUCBT, with highest viral loads on Day 100 assessment. The cumulative probability of developing BKPyV viremia by Day 100 was 0.52 (95% CI, 0.33-0.71). Six of 15 patients with BKPyV viremia experienced hemorrhagic cystitis by Day 100. By Day 100, there was a trend towards higher BKPyV viral loads in sera of patients with hemorrhagic cystitis than in those BKPyV viremic patients without hemorrhagic cystitis (p = 0.06). BKPyV viremia was associated with significantly higher anti-BKPyV IgM values at 6 months post-dUCBT (P = 0.003). BKPyV viremia occurs early after dUBCT and is associated with a detectable humoral immune response by 6 months post-dUBCT.

Inhibition of myostatin signaling through Notch activation following acute resistance exercise.

Myostatin is a TGFß family member and negative regulator of muscle size. Due to the complexity of the molecular pathway between myostatin mRNA/protein and changes in transcription, it has been difficult to understand whether myostatin plays a role in resistance exercise-induced skeletal muscle hypertrophy. To circumvent this problem, we determined the expression of a unique myostatin target gene, Mighty, following resistance exercise. Mighty mRNA increased by 6 h (82.9 ± 24.21%) and remained high out to 48 h (56.5 ± 19.67%) after resistance exercise. Further examination of the soleus, plantaris and tibialis anterior muscles showed that the change in Mighty mRNA at 6 h correlated with the increase in muscle size associated with this protocol (R(2) = 0.9996). The increase in Mighty mRNA occurred both independent of Smad2 phosphorylation and in spite of an increase in myostatin mRNA (341.8 ± 147.14% at 3 h). The myostatin inhibitor SKI remained unchanged. However, activated Notch, another potential inhibitor of TGFß signaling, increased immediately following resistance exercise (83 ± 11.2%) and stayed elevated out to 6 h (78 ± 16.6%). Electroportion of the Notch intracellular domain into the tibialis anterior resulted in an increase in Mighty mRNA (63 ± 13.4%) that was equivalent to the canonical Notch target HES-1 (94.4 ± 7.32%). These data suggest that acute resistance exercise decreases myostatin signaling through the activation of the TGFß inhibitor Notch resulting in a decrease in myostatin transcriptional activity that correlates well with muscle hypertrophy.

Glycogen content regulates peroxisome proliferator activated receptor-∂ (PPAR-∂) activity in rat skeletal muscle.

Performing exercise in a glycogen depleted state increases skeletal muscle lipid utilization and the transcription of genes regulating mitochondrial ß-oxidation. Potential candidates for glycogen-mediated metabolic adaptation are the peroxisome proliferator activated receptor (PPAR) coactivator-1α (PGC-1α) and the transcription factor/nuclear receptor PPAR-∂. It was therefore the aim of the present study to examine whether acute exercise with or without glycogen manipulation affects PGC-1α and PPAR-∂ function in rodent skeletal muscle. Twenty female Wistar rats were randomly assigned to 5 experimental groups (n = 4): control [CON]; normal glycogen control [NG-C]; normal glycogen exercise [NG-E]; low glycogen control [LG-C]; and low glycogen exercise [LG-E]). Gastrocnemius (GTN) muscles were collected immediately following exercise and analyzed for glycogen content, PPAR-∂ activity via chromatin immunoprecipitation (ChIP) assays, AMPK α1/α2 kinase activity, and the localization of AMPK and PGC-1α. Exercise reduced muscle glycogen by 47 and 75% relative to CON in the NG-E and LG-E groups, respectively. Exercise that started with low glycogen (LG-E) finished with higher AMPK-α2 activity (147%, p<0.05), nuclear AMPK-α2 and PGC-1α, but no difference in AMPK-α1 activity compared to CON. In addition, PPAR-∂ binding to the CPT1 promoter was significantly increased only in the LG-E group. Finally, cell reporter studies in contracting C2C12 myotubes indicated that PPAR-∂ activity following contraction is sensitive to glucose availability, providing mechanistic insight into the association between PPAR-∂ and glycogen content/substrate availability. The present study is the first to examine PPAR-∂ activity in skeletal muscle in response to an acute bout of endurance exercise. Our data would suggest that a factor associated with muscle contraction and/or glycogen depletion activates PPAR-∂ and initiates AMPK translocation in skeletal muscle in response to exercise.

A limited role for PI(3,4,5)P3 regulation in controlling skeletal muscle mass in response to resistance exercise.

BACKGROUND: Since activation of the PI3K/(protein kinase B; PKB/akt) pathway has been shown to alter muscle mass and growth, the aim of this study was to determine whether resistance exercise increased insulin like growth factor (IGF) I/phosphoinositide 3-kinase (PI3K) signalling and whether altering PI(3,4,5)P(3) metabolism genetically would increase load induced muscle growth. METHODOLOGY/PRINCIPAL FINDINGS: Acute and chronic resistance exercise in wild type and muscle specific PTEN knockout mice were used to address the role of PI(3,4,5)P(3) regulation in the development of skeletal muscle hypertrophy. Acute resistance exercise did not increase either IGF-1 receptor phosphorylation or IRS1/2 associated p85. Since insulin/IGF signalling to PI3K was unchanged, we next sought to determine whether inactivation of PTEN played a role in load-induced muscle growth. Muscle specific knockout of PTEN resulted in small but significant increases in heart (PTEN(+/+) = 5.00+/-0.02 mg/g, PTEN(-/-) = 5.50+/-0.09 mg/g), and TA (PTEN(+/+) = 1.74+/-0.04 mg/g, PTEN(-/-) = 1.89 +/-0.03) muscle mass, while the GTN, SOL, EDL and PLN remain unchanged. Following ablation, hypertrophy of the PLN, SOL or EDL muscles was similar between PTEN(-/-) and PTEN(+/+) animals. Even though there were some changes in overload-induced PKB and S6K1 phosphorylation, 1 hr following acute resistance exercise there was no difference in the phosphorylation state of S6K1 Thr389 between genotypes. CONCLUSIONS/SIGNIFICANCE: These data suggest that physiological loading does not lead to the enhanced activation of the PI3K/PKB/mTORC1 axis and that neither PI3K activation nor PTEN, and by extension PI(3,4,5)P(3) levels, play a significant role in adult skeletal muscle growth.

FoxK1 splice variants show developmental stage-specific plasticity of expression with temperature in the tiger pufferfish.

FoxK1 is a member of the highly conserved forkhead/winged helix (Fox) family of transcription factors and it is known to play a key role in mammalian muscle development and myogenic stem cell function. The tiger pufferfish (Takifugu rubripes) orthologue of mammalian FoxK1 (TFoxK1) has seven exons and is located in a region of conserved synteny between pufferfish and mouse. TFoxK1 is expressed as three alternative transcripts: TFoxK1-alpha, TFoxK1-gamma and TFoxK1-delta. TFoxK1-alpha is the orthologue of mouse FoxK1-alpha, coding for a putative protein of 558 residues that contains the forkhead and forkhead-associated domains typical of Fox proteins and shares 53% global identity with its mammalian homologue. TFoxK1-gamma and TFoxK1-delta arise from intron retention events and these transcripts translate into the same 344-amino acid protein with a truncated forkhead domain. Neither are orthologues of mouse FoxK1-beta. In adult fish, the TFoxK1 splice variants were differentially expressed between fast and slow myotomal muscle, as well as other tissues, and the FoxK1-alpha protein was expressed in myogenic progenitor cells of fast myotomal muscle. During embryonic development, TFoxK1 was transiently expressed in the developing somites, heart, brain and eye. The relative expression of TFoxK1-alpha and the other two alternative transcripts varied with the incubation temperature regime for equivalent embryonic stages and the differences were particularly marked at later developmental stages. The developmental expression pattern of TFoxK1 and its localisation to mononuclear myogenic progenitor cells in adult fast muscle indicate that it may play an essential role in myogenesis in T. rubripes.

Myogenin in model pufferfish species: Comparative genomic analysis and thermal plasticity of expression during early development.

Myogenin (Myog) is a muscle-specific basic helix-loop-helix transcription factor that plays an essential role in the specification and differentiation of myoblasts. The myogenin genes from the tiger pufferfish, Takifugu rubripes, and green-spotted pufferfish, Tetraodon nigroviridis, were cloned and a comparative genomic analysis performed. The gene encoding myogenin is composed of three exons and has a relatively similar genomic structure in T. rubripes, T. nigroviridis and human. Introns 1 and 2 were approximately 2-fold and 8-fold longer respectively in human than pufferfish. Myogenin is located in a 100 kb region of conserved synteny between these organisms, corresponding to chromosome 1 in human, chromosome 11 in T. nigroviridis and scaffold 208 in T. rubripes. Pufferfish myogenin contained a serine-rich region at the carboxyl terminus that is highly conserved amongst teleosts. During embryonic development of T. rubripes, myogenin was expressed in a rostral-caudal gradient in the developing somites and subsequently during the pharyngula period in the pectoral fin bud primordia, jaw muscles and extraocular muscle precursors. In T. rubripes, the time required to form a somite pair during the linear phase of somitogenesis ( identical withsomite-interval) was 122 min, 97 min and 50 min in embryos incubated at 15, 18 and 21 degrees C, respectively. Myogenin mRNA transcripts were quantified using qPCR and normalised to the highest level of expression. Peak myogenin expression occurred later with respect to developmental stage (standardised using somite-intervals) and was over 2-fold higher at 21 degrees C than at either 18 or 15 degrees C. Changes in the relative timing and intensity of myogenin expression are a potential mechanism for explaining thermal plasticity of muscle phenotype in larvae via effects on the differentiation programme.