An alternative splicing event in the Pax-3 paired domain identifies the linker region as a key determinant of paired domain DNA-binding activity

The PAX Genes: Roles in Development, Cancer, and Other Diseases

Since their 1986 discovery in Drosophila, Paired box (PAX) genes have been shown to play major roles in the early development of the eye, muscle, skeleton, kidney, and other organs. Consistent with their roles as master regulators of tissue formation, the PAX family members are evolutionarily conserved, regulate large transcriptional networks, and in turn can be regulated by a variety of mechanisms. Losses or mutations in these genes can result in developmental disorders or cancers. The precise mechanisms by which PAX genes control disease pathogenesis are well understood in some cases, but much remains to be explored. A deeper understanding of the biology of these genes, therefore, has the potential to aid in the improvement of disease diagnosis and the development of new treatments.

Post-Transcriptional Regulation in Skeletal Muscle Development, Repair, and Disease

Skeletal muscle formation is a complex process that requires tight spatiotemporal control of key myogenic factors. Emerging evidence suggests that RNA processing is crucial for the regulation of these factors, and that multiple post-transcriptional regulatory pathways work dependently and independently of one another to enable precise control of transcripts throughout muscle development and repair. Moreover, disruption of these pathways is implicated in neuromuscular disease, and the recent development of RNA-mediated therapies shows enormous promise in the treatment of these disorders. We discuss the overlapping post-transcriptional regulatory pathways that mediate muscle development, how these pathways are disrupted in neuromuscular disorders, and advances in RNA-mediated therapies that present a novel approach to the treatment of these diseases.

The Switch from NF-YAl to NF-YAs Isoform Impairs Myotubes Formation

NF-YA, the regulatory subunit of the trimeric transcription factor (TF) NF-Y, is regulated by alternative splicing (AS) generating two major isoforms, “long” (NF-YAl) and “short” (NF-YAs). Muscle cells express NF-YAl. We ablated exon 3 in mouse C2C12 cells by a four-guide CRISPR/Cas9n strategy, obtaining clones expressing exclusively NF-YAs (C2-YAl-KO). C2-YAl-KO cells grow normally, but are unable to differentiate. Myogenin and—to a lesser extent, MyoD— levels are substantially lower in C2-YAl-KO, before and after differentiation. Expression of the fusogenic Myomaker and Myomixer genes, crucial for the early phases of the process, is not induced. Myomaker and Myomixer promoters are bound by MyoD and Myogenin, and Myogenin overexpression induces their expression in C2-YAl-KO. NF-Y inactivation reduces MyoD and Myogenin, but not directly: the Myogenin promoter is CCAAT-less, and the canonical CCAAT of the MyoD promoter is not bound by NF-Y in vivo. We propose that NF-YAl, but not NF-YAs, maintains muscle commitment by indirectly regulating Myogenin and MyoD expression in C2C12 cells. These experiments are the first genetic evidence that the two NF-YA isoforms have functionally distinct roles.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

The expression and function of PAX3 in development and disease

The PAX3 gene encodes a member of the PAX family of transcription factors that is characterized by a highly conserved paired box motif. The PAX3 protein is a transcription factor consisting of an N-terminal DNA binding domain (containing a paired box and homeodomain) and a C-terminal transcriptional activation domain. This protein is expressed during development of skeletal muscle, central nervous system and neural crest derivatives, and regulates expression of target genes that impact on proliferation, survival, differentiation and motility in these lineages. Germline mutations of the murine Pax3 and human PAX3 genes cause deficiencies in these developmental lineages and result in the Splotch phenotype and Waardenburg syndrome, respectively. Somatic genetic rearrangements that juxtapose the PAX3 DNA binding domain to the transcriptional activation domain of other transcription factors deregulate PAX3 function and contribute to the pathogenesis of the soft tissue cancers alveolar rhabdomyosarcoma and biphenotypic sinonasal sarcoma. The wild-type PAX3 protein is also expressed in other cancers related to developmental lineages that normally express this protein and exerts phenotypic effects related to its normal developmental role.

Diversification of the muscle proteome through alternative splicing

Background

Skeletal muscles express a highly specialized proteome that allows the metabolism of energy sources to mediate myofiber contraction. This muscle-specific proteome is partially derived through the muscle-specific transcription of a subset of genes. Surprisingly, RNA sequencing technologies have also revealed a significant role for muscle-specific alternative splicing in generating protein isoforms that give specialized function to the muscle proteome.

Main body

In this review, we discuss the current knowledge with respect to the mechanisms that allow pre-mRNA transcripts to undergo muscle-specific alternative splicing while identifying some of the key trans-acting splicing factors essential to the process. The importance of specific splicing events to specialized muscle function is presented along with examples in which dysregulated splicing contributes to myopathies. Though there is now an appreciation that alternative splicing is a major contributor to proteome diversification, the emergence of improved “targeted” proteomic methodologies for detection of specific protein isoforms will soon allow us to better appreciate the extent to which alternative splicing modifies the activity of proteins (and their ability to interact with other proteins) in the skeletal muscle. In addition, we highlight a continued need to better explore the signaling pathways that contribute to the temporal control of trans-acting splicing factor activity to ensure specific protein isoforms are expressed in the proper cellular context.

Conclusions

An understanding of the signal-dependent and signal-independent events driving muscle-specific alternative splicing has the potential to provide us with novel therapeutic strategies to treat different myopathies.

Electronic supplementary material

The online version of this article (10.1186/s13395-018-0152-3) contains supplementary material, which is available to authorized users.

Alternative Splicing of Transcription Factors Genes in Muscle Physiology and Pathology

Skeletal muscle formation is a multi-step process that is governed by complex networks of transcription factors. The regulation of their functions is in turn multifaceted, including several mechanisms, among them alternative splicing (AS) plays a primary role. On the other hand, altered AS has a role in the pathogenesis of numerous muscular pathologies. Despite these premises, the causal role played by the altered splicing pattern of transcripts encoding myogenic transcription factors in neuromuscular diseases has been neglected so far. In this review, we systematically investigate what has been described about the AS patterns of transcription factors both in the physiology of the skeletal muscle formation process and in neuromuscular diseases, in the hope that this may be useful in re-evaluating the potential role of altered splicing of transcription factors in such diseases.

It's a bit over, is that ok? The subtle surplus from tandem alternative splicing

Tandem alternative splice sites (TASS) form a defined class of alternative splicing and give rise to mRNA insertion/deletion variants with only small size differences. Previous work has confirmed evolutionary conservation of TASS elements while many cases show only low tissue specificity of isoform ratios. We pinpoint stochasticity and noise as important methodological issues for the dissection of TASS isoform patterns. Resolving such uncertainties, a recent report showed regulation in a cell culture system, with shifts of alternative splicing isoform ratios dependent on cell density. This novel type of regulation affects not only multiple TASS isoforms, but also other alternative splicing classes, in a concerted manner. Here, we discuss how specific regulatory network architectures may be realized through the novel regulation type and highlight the role of differential isoform functions as a key step in order to better understand the functional role of TASS.

PAX transcription factors in neural crest development

The nine vertebrate PAX transcription factors (PAX1-PAX9) play essential roles during early development and organogenesis. Pax genes were identified in vertebrates using their homology with the Drosophila melanogaster paired gene DNA-binding domain. PAX1-9 functions are largely conserved throughout vertebrate evolution, in particular during central nervous system and neural crest development. The neural crest is a vertebrate invention, which gives rise to numerous derivatives during organogenesis, including neurons and glia of the peripheral nervous system, craniofacial skeleton and mesenchyme, the heart outflow tract, endocrine and pigment cells. Human and mouse spontaneous mutations as well as experimental analyses have evidenced the critical and diverse functions of PAX factors during neural crest development. Recent studies have highlighted the role of PAX3 and PAX7 in neural crest induction. Additionally, several PAX proteins - PAX1, 3, 7, 9 - regulate cell proliferation, migration and determination in multiple neural crest-derived lineages, such as cardiac, sensory, and enteric neural crest, pigment cells, glia, craniofacial skeleton and teeth, or in organs developing in close relationship with the neural crest such as the thymus and parathyroids. The diverse PAX molecular functions during neural crest formation rely on fine-tuned modulations of their transcriptional transactivation properties. These modulations are generated by multiple means, such as different roles for the various isoforms (formed by alternative splicing), or posttranslational modifications which alter protein-DNA binding, or carefully orchestrated protein-protein interactions with various co-factors which control PAX proteins activity. Understanding these regulations is the key to decipher the versatile roles of PAX transcription factors in neural crest development, differentiation and disease.

Characterization of Pax3 and Pax7 genes and their expression patterns during different development and growth stages of Japanese pufferfish Takifugu rubripes

Pax3 and Pax7 are the regulators and markers of muscle progenitors and satellite cells that contribute to the embryonic development and postembryonic growth of skeletal muscle in vertebrates, as well as to its repair and regeneration. However, information regarding them in vertebrate genome model, torafugu Takifugu rubripes, has remained unknown. Therefore, as an initial step, here we characterized Pax3 and Pax7 from torafugu and investigated their expression patterns during different developmental stages by RT-PCR. In silico analysis with the Fugu genome database (ver. 4.0) yielded two distinct genes each for Pax3 (Pax3a and Pax3b) and Pax7 (Pax7a and Pax7b). The 75th amino acid, glutamine (Gln75), from the N-terminus was replaced by proline in the paired box domain (PD) of Pax3a. One single cDNA clone encoding Pax3a had deletion of Gln75 in PD, suggesting the presence of alternatively spliced variants (Q+/Q-). This was further supported by identification of two adjacent alternative 3' splice acceptor sites which produce Pax3b Q+ (aagCAGGGA) and Q- (aagcagGGA) variants. Interestingly, torafugu Pax7a, but not Pax7b, had an insert encoding five amino acid residues (SGEAS) in a C-terminal region of PD in two out of three cDNA clones. Genomic analysis showed two alternate splice donor sites at exon 4 of Pax7a. In synteny analysis, torafugu Pax3a showed syntenic relationship with the corresponding regions in other teleosts only, whereas Pax3b and Pax7b showed high syntenic relationship with the corresponding regions of both mammals and other teleosts. RT-PCR revealed that expression of Pax3a and Pax3b transcripts was restricted to embryonic stages only, whereas those of Pax7a and Pax7b was continued to be expressed in larvae and importantly those of Pax7a were found in adult skeletal muscles. Therefore, Pax3 appears to be most important for primary myogenesis and Pax7 for secondary myogenesis and growth by hyperplasia in fish. In this regard, the transcripts of torafugu Pax3 and Pax7 genes might be used for further investigation as a marker for identification of muscle precursor cells during different phases of growth, and this ambiguity is the next target of our research.

The duplicated paired box protein 7 (pax7) genes differentially transcribed during Japanese flounder (Paralichthys olivaceus) embryogenesis

PAX are important regulators of developmental processes. PAX7 plays crucial roles in patterning of the dorsal central nervous system (CNS), neural crest (NC), and skeletal muscle. Here, we identified six spliced isoforms of pax7a and one pax7b and characterized their expression patterns. All of flounder Pax7a-1, Pax7a-2, Pax7a-3, and Pax7b contain a conserved paired domain (PD), an octapeptide motif (OP), and a paired type homeodomain (HD). However, the PD of Pax7a-4 and the HD of Pax7a-5 are not intact, and there is no HD in Pax7a-4 and Pax7a-6. pax7a and pax7b show distinct spatiotemporal expression patterns during embryogenesis. Whole-mount in situ hybridization demonstrates that the expression patterns of pax7a and pax7b are overlapping but distinguishable in the dorsal central nervous system. pax7a is expressed in most part of the brain and the neural tube, while pax7b is expressed exclusively in the diencephalon and the midbrain. In addition, pax7a is also expressed in the cranial NC and the trunk NC. RT-PCR results show that there were different expression patterns between the different isoforms. These results indicate subfunction partitioning of the duplicated pax7 genes. The duplicated pax7 may provide additional flexibility in fine-tuning neurogenesis and somitogenesis.

Pushing the annotation of cellular activities to a higher resolution: Predicting functions at the isoform level

In past decades, the experimental determination of protein functions was expensive and time-consuming, so numerous computational methods were developed to speed up and guide the process. However, most of these methods predict protein functions at the gene level and do not consider the fact that protein isoforms (translated from alternatively spliced transcripts), not genes, are the actual function carriers. Now, high-throughput RNA-seq technology is providing unprecedented opportunities to unravel protein functions at the isoform level. In this article, we review recent progress in the high-resolution functional annotations of protein isoforms, focusing on two methods developed by the authors. Both methods can integrate multiple RNA-seq datasets for comprehensively characterizing functions of protein isoforms.

Pax factors in transcription and epigenetic remodelling

The nine Pax transcription factors that constitute the mammalian family of paired domain (PD) factors play key roles in many developmental processes. As DNA binding transcription factors, they exhibit tremendous variability and complexity in their DNA recognition patterns. This is ascribed to the presence of multiple DNA binding structural domains, namely helix-turn-helix (HTH) domains. The PD contains two HTH subdomains and four of the nine Pax factors have an additional HTH domain, the homeodomain (HD). We now review these diverse DNA binding modalities together with their properties as transcriptional activators and repressors. The action of Pax factors on gene expression is also exerted through recruitment of chromatin remodelling complexes that introduce either activating or repressive chromatin marks. Interestingly, the recent demonstration that Pax7 has pioneer activity, the unique property to "open" chromatin, further underlines the mechanistic versatility and the developmental importance of these factors.

nagnag: Identification and quantification of NAGNAG alternative splicing using RNA-Seq data

Regulation of proteome diversity by alternative splicing has been widely demonstrated in plants and animals. NAGNAG splicing, which was recently defined as a tissue specific event, results in the production of two distinct isoforms that are distinguished by three nucleotides (NAG) as a consequence of the intron proximal or distal to the splice site. Since the NAGNAG mechanism is not well characterized, tools for the identification and quantification of NAGNAG splicing events remain under-developed. Here we report nagnag, an R-based NAGNAG splicing detection tool, which accurately identifies and quantifies NAGNAG splicing events using RNA-Seq. Overall, nagnag produces user-friendly visualization reports and highlights differences between the DNA/RNA/protein across the identified isoforms of the reported gene. The package is available on https://sourceforge.net/projects/nagnag/files/; or http://genome.sdau.edu.cn/research/software/nagnag.html.

The olive flounder (Paralichthys olivaceus) Pax3 homologues are highly conserved, encode multiple isoforms and show unique expression patterns

Pax genes encode a highly conserved family of transcription factors that play crucial roles in the formation of tissues and organs during development. Pax3 plays crucial roles in patterning of the dorsal central nervous system (CNS), neural crest and skeletal muscle. Here, we identified two spliced isoforms of Pax3a and three spliced isoforms of Pax3b and characterized their expression patterns. Both of flounder Pax3a-1 and Pax3b-1 contain the conserved paired domain (PD), an octapeptide motif (OP), and a paired type homeodomain (HD). But the PD domain in Pax3a-2 and Pax3b-3 is not intact and there is no HD in Pax3b-2 and Pax3b-3. Pax3a and Pax3b show distinct temporal expression patterns during embryogenesis. Whole-mount in situ hybridization demonstrates that Pax3a and Pax3b are expressed in overlapping patterns in the dorsal central nervous system, with some subtle regional differences between the two genes. In addition, Pax3a is scattered in the somites while Pax3b is specifically expressed in the newly forming somites. RT-PCR results have shown that there were different expression patterns between the different isoforms. These results indicate subfunction partitioning of the duplicated Pax3 genes. The duplicated Pax3 may provide additional flexibility in fine-tuning neurogenesis and somitogenesis.

Computational evidence of NAGNAG alternative splicing in human large intergenic noncoding RNA

NAGNAG alternative splicing plays an essential role in biological processes and represents a highly adaptable system for posttranslational regulation of gene function. NAGNAG alternative splicing impacts a myriad of biological processes. Previous studies of NAGNAG largely focused on messenger RNA. To the best of our knowledge, this is the first study testing the hypothesis that NAGNAG alternative splicing is also operative in large intergenic noncoding RNA (lincRNA). The RNA-seq data sets from recent deep sequencing studies were queried to test our hypothesis. NAGNAG alternative splicing of human lincRNA was identified while querying two independent RNA-seq data sets. Within these datasets, 31 NAGNAG alternative splicing sites were identified in lincRNA. Notably, most exons of lincRNA containing NAGNAG acceptors were longer than those from protein-coding genes. Furthermore, presence of CAG coding appeared to participate in the splice site selection. Finally, expression of the isoforms of NAGNAG lincRNA exhibited tissue specificity. Together, this study improves our understanding of the NAGNAG alternative splicing in lincRNA.

A subtle alternative splicing event gives rise to a widely expressed human RNase k isoform

Subtle alternative splicing leads to the formation of RNA variants lacking or including a small number of nucleotides. To date, the impact of subtle alternative splicing phenomena on protein biosynthesis has been studied in frame-preserving incidents. On the contrary, mRNA isoforms derived from frame-shifting events were poorly studied and generally characterized as non-coding. This work provides evidence for a frame-shifting subtle alternative splicing event which results in the production of a novel protein isoform. We applied a combined molecular approach for the cloning and expression analysis of a human RNase κ transcript (RNase κ-02) which lacks four consecutive bases compared to the previously isolated RNase κ isoform. RNase κ-02 mRNA is expressed in all human cell lines tested end encodes the synthesis of a 134-amino-acid protein by utilizing an alternative initiation codon. The expression of RNase κ-02 in the cytoplasm of human cells was verified by Western blot and immunofluorescence analysis using a specific polyclonal antibody developed on the basis of the amino-acid sequence difference between the two protein isoforms. The results presented here show that subtle changes during mRNA splicing can lead to the expression of significantly altered protein isoforms.

Pax3 isoforms in sensory neurogenesis: expression and function in the ophthalmic trigeminal placode

BACKGROUND

In the trigeminal placode, Pax3 is classified as necessary but not sufficient for sensory neuron differentiation. One hypothesis is that different Pax3 isoforms regulate cellular differentiation uniquely. Pax3 is known to sometimes activate and sometimes repress gene transcription, and its activity can be dependent on the isoforms present. Pax3 isoforms had not previously been characterized in chick sensory neurogenesis.

RESULTS

Reverse transcriptase-polymerase chain reaction (PCR) analysis revealed three well-expressed Pax3 splice variants: full-length (flPax3), Pax3V1, and Pax3V2. Each was characterized for its effect on neurogenesis by misexpression in placodal ectoderm. The differences observed were more apparent under conditions of enhanced neurogenesis (by means of Notch inhibition), where flPax3 and Pax3V1 caused failed differentiation, while Pax3V2 misexpression resembled the neuronal differentiation seen in controls. Quantitative PCR analysis revealed a progressive increase in Pax3 expression, but no significant change in relative isoform expression. Of interest, Notch inhibition led to a significant increase in Pax3 expression.

CONCLUSIONS

We can conclude that: (1) flPax3 and Pax3V1 inhibit neuronal differentiation; (2) Pax3V2 is permissive for neuronal differentiation; (3) while absolute levels change over time, relative splice form expression levels are largely maintained in the trigeminal placode domain; and (4) Pax3 expression generally increases in response to Notch inhibition.

Genome-wide study of NAGNAG alternative splicing in Arabidopsis

NAGNAG alternative splicing is one type of alternative splicing in mammals and plants. There are two opposite arguments regarding the mechanism of this NAGNAG event, i.e. whether splice variation is controllable by the cell or is just biological noise. In this paper, we systematically investigated NAGNAG acceptors in Arabidopsis thaliana using both cDNA/EST and RNA-Seq data. We identified 9,473 NAGNAG motifs, including 529 cDNA/EST-confirmed NAGNAG acceptors. A nomenclature tree for this type of alternative splicing was defined based on the cDNA/EST validation, location in the exon, sequence and expression level. Low expression of some NAGNAG motifs was observed in various tissues or pathogen-infected samples, indicating the existence of background splicing. Tissue-specific or treatment-specific differences in the dynamic profiles suggest that some NAGNAG acceptors are highly regulated.

High-resolution functional annotation of human transcriptome: predicting isoform functions by a novel multiple instance-based label propagation method

Alternative transcript processing is an important mechanism for generating functional diversity in genes. However, little is known about the precise functions of individual isoforms. In fact, proteins (translated from transcript isoforms), not genes, are the function carriers. By integrating multiple human RNA-seq data sets, we carried out the first systematic prediction of isoform functions, enabling high-resolution functional annotation of human transcriptome. Unlike gene function prediction, isoform function prediction faces a unique challenge: the lack of the training data--all known functional annotations are at the gene level. To address this challenge, we modelled the gene-isoform relationships as multiple instance data and developed a novel label propagation method to predict functions. Our method achieved an average area under the receiver operating characteristic curve of 0.67 and assigned functions to 15 572 isoforms. Interestingly, we observed that different functions have different sensitivities to alternative isoform processing, and that the function diversity of isoforms from the same gene is positively correlated with their tissue expression diversity. Finally, we surveyed the literature to validate our predictions for a number of apoptotic genes. Strikingly, for the famous 'TP53' gene, we not only accurately identified the apoptosis regulation function of its five isoforms, but also correctly predicted the precise direction of the regulation.

Pax3 and Zic1 trigger the early neural crest gene regulatory network by the direct activation of multiple key neural crest specifiers

Neural crest development is orchestrated by a complex and still poorly understood gene regulatory network. Premigratory neural crest is induced at the lateral border of the neural plate by the combined action of signaling molecules and transcription factors such as AP2, Gbx2, Pax3 and Zic1. Among them, Pax3 and Zic1 are both necessary and sufficient to trigger a complete neural crest developmental program. However, their gene targets in the neural crest regulatory network remain unknown. Here, through a transcriptome analysis of frog microdissected neural border, we identified an extended gene signature for the premigratory neural crest, and we defined novel potential members of the regulatory network. This signature includes 34 novel genes, as well as 44 known genes expressed at the neural border. Using another microarray analysis which combined Pax3 and Zic1 gain-of-function and protein translation blockade, we uncovered 25 Pax3 and Zic1 direct targets within this signature. We demonstrated that the neural border specifiers Pax3 and Zic1 are direct upstream regulators of neural crest specifiers Snail1/2, Foxd3, Twist1, and Tfap2b. In addition, they may modulate the transcriptional output of multiple signaling pathways involved in neural crest development (Wnt, Retinoic Acid) through the induction of key pathway regulators (Axin2 and Cyp26c1). We also found that Pax3 could maintain its own expression through a positive autoregulatory feedback loop. These hierarchical inductions, feedback loops, and pathway modulations provide novel tools to understand the neural crest induction network.

Arginine68 is an essential residue for the C-terminal cleavage of human Atg8 family proteins

Background

Autophagy is a conserved cellular process that degrades and recycles cytoplasmic components via a lysosomal pathway. The phosphatidylethanolamine (PE)-conjugation of the Atg8 protein plays an important role in the yeast autophagy process. In humans, six Atg8 homologs, including MAP1LC3A, MAP1LC3B, MAP1LC3C (refer to LC3A, LC3B, and LC3C hereafter), GABARAP, GABARAPL1, and GABARAPL2 have been reported. All of them can be conjugated to PE through a ubiquitin-like conjugation system, and be located to autophagosomes.

Results

In this study, we found 3 new alternative splicing isoforms in LC3B, GABARAP, and GABARAPL1, (designated as LC3B-a, GABARAP-a and GABARAPL1-a, respectively). None of them can go through the PE-conjugation process and be located to autophagosomes. Interestingly, compared with LC3B, LC3B-a has a single amino acid (Arg68) deletion due to the NAGNAG alternative splicing in intron 3. Through structural simulations, we found that the C-terminal tail of LC3B-a is less mobile than that of LC3B, thus affecting its C-terminal cleavage by human ATG4 family proteins. Furthermore, we found that Arg68 is an essential residue facilitating the interaction between human Atg8 family proteins and ATG4B by forming a salt bridge with Asp171 of ATG4B. Depletion of this salt bridge reduces autophagosomes formation and autophagic flux under both normal and nutrition starvation conditions.

Conclusions

These results suggest Arg68 is an essential residue for the C-terminal cleavage of Atg8 family proteins during the autophagy process.

A novel missense mutation of the paired box 3 gene in a Turkish family with Waardenburg syndrome type 1

PURPOSE

Screening of mutations in the paired box 3 (PAX3) gene in three generations of a Turkish family with Waardenburg syndrome type 1 (WS1).

METHODS

WS1 was diagnosed in a 13-month-old girl according to the WS Consortium criteria. Detailed family history of the proband revealed eight affected members in three generations. Routine clinical and audiological examination and ophthalmologic evaluation were performed on eight affected and five healthy members of the study family. Dystopia canthorum was detected in all affected patients; however, a brilliant blue iris was present in five patients who also had mild retinal hypopigmentation. Genomic DNA was extracted from the peripheral blood of affected and unaffected individuals in the family as well as 50 unrelated healthy volunteers. All coding exons and adjacent intronic regions of PAX3 were sequenced directly.

RESULTS

A novel missense heterozygous c.788T>G mutation was identified in eight patients. This nucleotide alteration was not found in unaffected members of the study family or in the 50 unrelated control subjects. The mutation causes V263G amino-acid substitution in the homeodomain of the PAX3 protein, which represents the 45(th) residue of helix 3.

CONCLUSIONS

We identified a novel missense c.788T>G mutation in PAX3 in a family with Waardenburg syndrome with intrafamilial phenotypic heterogeneity.

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ∼15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.

Alternate PAX3 and PAX7 C-terminal isoforms in myogenic differentiation and sarcomagenesis

OBJECTIVE

Pax3 and Pax7 are closely related genes that are involved in commitment of cells to a myogenic lineage during skeletal muscle development and regeneration. Several Pax3 and Pax7 transcripts are expressed from the genes, generating different isoforms with potentially distinct DNA binding and transactivation properties. The aim of this study was to investigate the implication of Pax3 and Pax7 C-terminal isoforms during myogenic differentiation and tumorigenesis, since fusions involving these genes are commonly associated with alveolar rhabdomyosarcoma (ARMS).

METHODS

Uncommitted (mouse mesenchymal stem cells, MSCs) and committed (C2C12) myogenic precursor cells were stably transfected with PAX3/FKHR and PAXC7/ FKHR fusion genes. We analysed gene and protein expression comparing the newly generated cells with the parental cells, to determine the functional importance of Pax3 and Pax7 C-terminal isoforms.

RESULTS

We found that the transcript Pax3c was expressed at low levels in undifferentiated C2C12 and MSCs cells, but its expression levels increased considerably at later stages of differentiation. However, expression levels of Pax3d transcript increased only slightly after differentiation. Pax7 transcripts, present before differentiation in committed C2C12 cells, but absent in uncommitted MSCs, increased noticeably in MSCs after differentiation. We also found that the presence of PAX/FKHR fusions prevented both C2C12 and MSC cells from terminal myogenic differentiation and increased the expression of discrete endogenous Pax3/7 transcripts, in particular Pax3d and Pax7B.

CONCLUSIONS

Our results suggest that both Pax3 and Pax7 transcripts are required for commitment of cells to the myogenic lineage, with each transcript having a distinct role. More specifically, the Pax3c isoform may be required for terminal myogenic differentiation whereas the Pax3d isoform may be involved in undifferentiated cell maintenance and/or proliferation.

Evolution of major milk proteins in Mus musculus and Mus spretus mouse species: a genoproteomic analysis

Background

Due to their high level of genotypic and phenotypic variability, Mus spretus strains were introduced in laboratories to investigate the genetic determinism of complex phenotypes including quantitative trait loci. Mus spretus diverged from Mus musculus around 2.5 million years ago and exhibits on average a single nucleotide polymorphism (SNP) in every 100 base pairs when compared with any of the classical laboratory strains. A genoproteomic approach was used to assess polymorphism of the major milk proteins between SEG/Pas and C57BL/6J, two inbred strains of mice representative of Mus spretus and Mus musculus species, respectively.

Results

The milk protein concentration was dramatically reduced in the SEG/Pas strain by comparison with the C57BL/6J strain (34 ± 9 g/L vs. 125 ± 12 g/L, respectively). Nine major proteins were identified in both milks using RP-HPLC, bi-dimensional electrophoresis and MALDI-Tof mass spectrometry. Two caseins (β and α_s1) and the whey acidic protein (WAP), showed distinct chromatographic and electrophoresis behaviours. These differences were partly explained by the occurrence of amino acid substitutions and splicing variants revealed by cDNA sequencing. A total of 34 SNPs were identified in the coding and 3'untranslated regions of the SEG/Pas Csn1s1 (11), Csn2 (7) and Wap (8) genes. In addition, a 3 nucleotide deletion leading to the loss of a serine residue at position 93 was found in the SEG/Pas Wap gene.

Conclusion

SNP frequencies found in three milk protein-encoding genes between Mus spretus and Mus musculus is twice the values previously reported at the whole genome level. However, the protein structure and post-translational modifications seem not to be affected by SNPs characterized in our study. Splicing mechanisms (cryptic splice site usage, exon skipping, error-prone junction sequence), already identified in casein genes from other species, likely explain the existence of multiple α_s1-casein isoforms both in SEG/Pas and C57BL/6J strains. Finally, we propose a possible mechanism by which the hallmark tandem duplication of a 18-nt exon (14 copies) may have occurred in the mouse genome.

Sequence features involved in the mechanism of 3' splice junction wobbling

BACKGROUND

Alternative splicing is an important mechanism mediating the diversified functions of genes in multicellular organisms, and such event occurs in around 40-60% of human genes. Recently, a new splice-junction wobbling mechanism was proposed that subtle modifications exist in mRNA maturation by alternatively choosing at 5'- GTNGT and 3'- NAGNAG, which created single amino acid insertion and deletion isoforms.

RESULTS

By browsing the Alternative Splicing Database information, we observed that most 3' alternative splice site choices occur within six nucleotides of the dominant splice site and the incidence significantly decreases further away from the dominant acceptor site. Although a lower frequency of alternative splicing occurs within the intronic region (alternative splicing at the proximal AG) than in the exonic region (alternative splicing at the distal AG), alternative AG sites located within the intronic region show stronger potential as the acceptor. These observations revealed that the choice of 3' splice sites during 3' splicing junction wobbling could depend on the distance between the duplicated AG and the branch point site (BPS). Further mutagenesis experiments demonstrated that the distance of AG-to-AG and BPS-to-AG can greatly influence 3' splice site selection. Knocking down a known alternative splicing regulator, hSlu7, failed to affect wobble splicing choices.

CONCLUSION

Our results implied that nucleotide distance between proximal and distal AG sites has an important regulatory function. In this study, we showed that occurrence of 3' wobble splicing occurs in a distance-dependent manner and that most of this wobble splicing is probably caused by steric hindrance from a factor bound at the neighboring tandem motif sequence.

Mutational analysis of the eyeless gene and phenotypic rescue reveal that an intact Eyeless protein is necessary for normal eye and brain development in Drosophila

Pax6 genes encode evolutionarily highly conserved transcription factors that are required for eye and brain development. Despite the characterization of mutations in Pax6 homologs in a range of organisms, and despite functional studies, it remains unclear what the relative importance is of the various parts of the Pax6 protein. To address this, we have studied the Drosophila Pax6 homolog eyeless. Specifically, we have generated new eyeless alleles, each with single missense mutations in one of the four domains of the protein. We show that these alleles result in abnormal eye and brain development while maintaining the OK107 eyeless GAL4 activity from which they were derived. We performed in vivo functional rescue experiments by expressing in an eyeless-specific pattern Eyeless proteins in which either the paired domain, the homeodomain, or the C-terminal domain was deleted. Rescue of the eye and brain phenotypes was only observed when full-length Eyeless was expressed, while all deletion constructs failed to rescue. These data, along with the phenotypes observed in the four newly characterized eyeless alleles, demonstrate the requirement for an intact Eyeless protein for normal Drosophila eye and brain development. They also suggest that some endogenous functions may be obscured in ectopic expression experiments.

Molecular determinants and evolutionary dynamics of wobble splicing

Alternative splicing at tandem splice sites (wobble splicing) is widespread in many species, but the mechanisms specifying the tandem sites remain poorly understood. Here, we used synaptotagmin I as a model to analyze the phylogeny of wobble splicing spanning more than 300 My of insect evolution. Phylogenetic analysis indicated that the occurrence of species-specific wobble splicing was related to synonymous variation at tandem splice sites. Further mutagenesis experiments demonstrated that wobble splicing could be lost by artificially induced synonymous point mutations due to destruction of splice acceptor sites. In contrast, wobble splicing could not be correctly restored through mimicking an ancestral tandem acceptor by artificial synonymous mutation in in vivo splicing assays, which suggests that artificial tandem splice sites might be incompatible with normal wobble splicing. Moreover, combining comparative genomics with hybrid minigene analysis revealed that alternative splicing has evolved from the 3′ tandem donor to the 5′ tandem acceptor in Culex pipiens, as a result of an evolutionary shift of cis element sequences from 3′ to 5′ splice sites. These data collectively suggest that the selection of tandem splice sites might not simply be an accident of history but rather in large part the result of coevolution between splice site and cis element sequences as a basis for wobble splicing. An evolutionary model of wobble splicing is proposed.

Pax genes in embryogenesis and oncogenesis

The paired box genes are a family of nine developmental control genes, which in human beings (PAX) and mice (Pax) encode nuclear transcription factors. The temporal and spatial expressions of these highly conserved genes are tightly regulated during foetal development including organogenesis. PAY/Paxgenes are switched off during the terminal differentiation of most structures. Specific mutations within a number of PAX/Pax genes lead to developmental abnormalities in both human beings and mice. Mutation in PAX3 causes Waardenburg syndrome, and craniofacial-deafness-hand syndrome. The Splotch phenotype in mouse exhibits defects in neural crest derivatives such as, pigment cells, sympathetic ganglia and cardiac neural crest-derived structures. The PAX family also plays key roles in several human malignancies. In particular, PAX3 is involved in rhabdomyosarcoma and tumours of neural crest origin, including melanoma and neuroblastoma. This review critically evaluates the roles of PAX/Pax in oncogenesis. It especially highlights recent advances in knowledge of how their genetic alterations directly interfere in the transcriptional networks that regulate cell differentiation, proliferation, migration and survival and may contribute to oncogenesis.

Widespread and subtle: alternative splicing at short-distance tandem sites

Alternative splicing at donor or acceptor sites located just a few nucleotides apart is widespread in many species. It results in subtle changes in the transcripts and often in the encoded proteins. Several of these tandem splice events contribute to the repertoire of functionally different proteins, whereas many are neutral or deleterious. Remarkably, some of the functional events are differentially spliced in tissues or developmental stages, whereas others exhibit constant splicing ratios, indicating that function is not always associated with differential splicing. Stochastic splice site selection seems to play a major role in these processes. Here, we review recent progress in understanding functional and evolutionary aspects as well as the mechanism of splicing at short-distance tandem sites.

Assessing the fraction of short-distance tandem splice sites under purifying selection

Many alternative splice events result in subtle mRNA changes, and most of them occur at short-distance tandem donor and acceptor sites. The splicing mechanism of such tandem sites likely involves the stochastic selection of either splice site. While tandem splice events are frequent, it is unknown how many are functionally important. Here, we use phylogenetic conservation to address this question, focusing on tandems with a distance of 3-9 nucleotides. We show that previous contradicting results on whether alternative or constitutive tandem motifs are more conserved between species can be explained by a statistical paradox (Simpson's paradox). Applying methods that take biases into account, we found higher conservation of alternative tandems in mouse, dog, and even chicken, zebrafish, and Fugu genomes. We estimated a lower bound for the number of alternative sites that are under purifying (negative) selection. While the absolute number of conserved tandem motifs decreases with the evolutionary distance, the fraction under selection increases. Interestingly, a number of frameshifting tandems are under selection, suggesting a role in regulating mRNA and protein levels via nonsense-mediated decay (NMD). An analysis of the intronic flanks shows that purifying selection also acts on the intronic sequence. We propose that stochastic splice site selection can be an advantageous mechanism that allows constant splice variant ratios in situations where a deviation in this ratio is deleterious.

Overlapping alternative donor splice sites in the human genome

Over 50% of donor splice sites in the human genome have a potential alternative donor site at a distance of three to six nucleotides. Conservation of these potential sites is determined by the consensus requirements and by its exonic or intronic location. Several hundred pairs of overlapping sites are confirmed to be alternatively spliced as both sites in a pair are supported by a protein, by a full-length mRNA, or by expressed sequence tags (ESTs) from at least two independent clone libraries. Overlapping sites may clash with consensus requirements. Pairs with a site shift of four nucleotides are the most abundant, despite the frameshift in the protein-coding region that they introduce. The site usage in pairs is usually uneven, and the major site is more frequently conserved in other mammalian genomes. Overlapping alternative donor sites and acceptor sites may have different functional roles: alternative splicing of overlapping acceptor sites leads mainly to microvariations in protein sequences; whereas alternative donor sites often lead to frameshifts and thus either yield major differences in the protein sequence and structure, or generate nonsense-mediated decay-inducing mRNA isoforms likely involved in regulated unproductive splicing pathways.

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.

Phylogenetically widespread alternative splicing at unusual GYNGYN donors

BACKGROUND

Splice donor sites have a highly conserved GT or GC dinucleotide and an extended intronic consensus sequence GTRAGT that reflects the sequence complementarity to the U1 snRNA. Here, we focus on unusual donor sites with the motif GYNGYN (Y stands for C or T; N stands for A, C, G, or T).

RESULTS

While only one GY functions as a splice donor for the majority of these splice sites in human, we provide computational and experimental evidence that 110 (1.3%) allow alternative splicing at both GY donors. The resulting splice forms differ in only three nucleotides, which results mostly in the insertion/deletion of one amino acid. However, we also report the insertion of a stop codon in four cases. Investigating what distinguishes alternatively from not alternatively spliced GYNGYN donors, we found differences in the binding to U1 snRNA, a strong correlation between U1 snRNA binding strength and the preferred donor, over-represented sequence motifs in the adjacent introns, and a higher conservation of the exonic and intronic flanks between human and mouse. Extending our genome-wide analysis to seven other eukaryotic species, we found alternatively spliced GYNGYN donors in all species from mouse to Caenorhabditis elegans and even in Arabidopsis thaliana. Experimental verification of a conserved GTAGTT donor of the STAT3 gene in human and mouse reveals a remarkably similar ratio of alternatively spliced transcripts in both species.

CONCLUSION

In contrast to alternative splicing in general, GYNGYN donors in addition to NAGNAG acceptors enable subtle protein variations.

TassDB: a database of alternative tandem splice sites

Subtle alternative splice events at tandem splice sites are frequent in eukaryotes and substantially increase the complexity of transcriptomes and proteomes. We have developed a relational database, TassDB (TAndem Splice Site DataBase), which stores extensive data about alternative splice events at GYNGYN donors and NAGNAG acceptors. These splice events are of subtle nature since they mostly result in the insertion/deletion of a single amino acid or the substitution of one amino acid by two others. Currently, TassDB contains 114 554 tandem splice sites of eight species, 5209 of which have EST/mRNA evidence for alternative splicing. In addition, human SNPs that affect NAGNAG acceptors are annotated. The database provides a user-friendly interface to search for specific genes or for genes containing tandem splice sites with specific features as well as the possibility to download large datasets. This database should facilitate further experimental studies and large-scale bioinformatics analyses of tandem splice sites. The database is available at .

Characterisation and expression of the paired box protein 7 (Pax7) gene in polymorphic Arctic charr (Salvelinus alpinus)

Arctic charr (Salvelinus alpinus L.) from Lake Thingvallavatn, Iceland occur as four distinct morphs: large benthivorous (LB), dwarf benthivorous (DB), piscivorous (PI) and planktonivorous (PL). The morphs differ with respect to body size, head morphology, growth rate, and life history. The aim of this study was to investigate the paired box protein 7 (Pax7) gene as a candidate for such polymorphisms due to its importance in cranio-facial, skeletal muscle, and central nervous system development. No variation in coding and intronic sequences was found between morphs. We identified 10 alternate Pax7 isoforms with insertions/deletions: a four-residue (GNRT) deletion, a GEASS insertion truncated by the first serine residue (GEAS), and a thirteen-residue insertion (GQYA/TGPEYVYCGT). The latter insertion with a threonine (T) contains a putative casein kinase II (CK-2) phosphorylation site. Pax7 spatial expression patterns were identical in embryos of DB-, LB-, and PL-morphs, and were similar to those described for zebrafish Pax7c, but a difference in temporal expression for segmentation was observed between DB and LB morphs. At the end of segmentation, novel expression was observed in the mandibular region as two bilateral domains. The potential role of multiple alternative splicing of the Pax7 gene for the generation of different Arctic charr morphs is briefly discussed.

Alternative splicing regulation at tandem 3' splice sites

Alternative splicing (AS) constitutes a major mechanism creating protein diversity in humans. Previous bioinformatics studies based on expressed sequence tag and mRNA data have identified many AS events that are conserved between humans and mice. Of these events, ∼25% are related to alternative choices of 3′ and 5′ splice sites. Surprisingly, half of all these events involve 3′ splice sites that are exactly 3 nt apart. These tandem 3′ splice sites result from the presence of the NAGNAG motif at the acceptor splice site, recently reported to be widely spread in the human genome. Although the NAGNAG motif is common in human genes, only a small subset of sites with this motif is confirmed to be involved in AS. We examined the NAGNAG motifs and observed specific features such as high sequence conservation of the motif, high conservation of ∼30 bp at the intronic regions flanking the 3′ splice site and overabundance of cis-regulatory elements, which are characteristic of alternatively spliced tandem acceptor sites and can distinguish them from the constitutive sites in which the proximal NAG splice site is selected. Our findings imply that AS at tandem splice sites and constitutive splicing of the distal NAG are highly regulated.

Frequent occurrence of protein isoforms with or without a single amino acid residue by subtle alternative splicing: the case of Gln in DRPLA affects subcellular localization of the products

Protein isoforms with or without a single amino acid residue make a subtle difference. It has been documented on a few genes that alternative splicing generated such isoforms; however, the fact has attracted little attention. We became aware of a subtle sequence difference in DRPLA, a polyglutamine disease gene for dentatorubral pallidoluysian atrophy. Some reported cDNA sequences lacked 3 nucleotides (nt) (CAG), which were positioned apart from the expandable and polymorphic CAG repeats and also coded for glutamine. We experimentally confirmed that the difference was indeed generated by alternative splicing utilizing two acceptors separated by 3 nt. In DRPLA, the expression ratio of two mRNA isoforms was almost constant among tissues, with the CAG-included form being major. The glutamine-included protein isoform was more predominantly localized in the nucleus. Database searching revealed that alternative splice acceptors, as well as donors, are frequently situated very close to each other. We experimentally confirmed two mRNA isoforms of 3 nt difference in more than 200 cases by RT-PCR and found interesting features associated with this phenomena. Inclusion of 3 nt tends to result in single amino acid inclusion despite the phase of translational frame. The expression ratio sometimes varied extensively among tissues.

Regulation of gene expression by Pax6 in ocular cells: a case of tissue-preferred expression of crystallins in lens

Lens development is an excellent model for genetic and biochemical studies of embryonic induction, cell cycle regulation, cellular differentiation and signal transduction. Differentiation of lens is characterized by lens-preferred expression and accumulation of water-soluble proteins, crystallins. Crystallins are required for light transparency, refraction and maintenance of lens integrity. Here, we review mechanisms of lens-preferred expression of crystallin genes by employing synergism between developmentally regulated DNA-binding transcription factors: Pax6, c-Maf, MafA/L-Maf, MafB, NRL, Sox2, Sox1, RARbeta/RXRbeta, RORalpha, Prox1, Six3, gammaFBP-B and HSF2. These factors are differentially expressed in lens precursor cells, lens epithelium and primary and secondary lens fibers. They exert their function in combination with ubiquitously expressed factors (e.g. AP-1, CREB, pRb, TFIID and USF) and co-activators/chromatin remodeling proteins (e.g. ASC-2 and CBP/p300). A special function belongs to Pax6, a paired domain and homeodomain-containing protein, which is essential for lens formation. Pax6 is expressed in lens progenitor cells before the onset of crystallin expression and it serves as an important regulatory factor required for expression of c-Maf, MafA/L-Maf, Six3, Prox1 and retinoic acid signaling both in lens precursor cells and the developing lens. The roles of these factors are illustrated by promoter studies of mouse alphaA-, alphaB-, gammaF- and guinea pig zeta-crystallins. Pax6 forms functional complexes with a number of transcription factors including the retinoblastoma protein, pRb, MafA, Mitf and Sox2. We present novel data showing that pRb antagonizes Pax6-mediated activation of the alphaA-crystallin promoter likely by inhibiting binding of Pax6 to DNA.

Isolation and expression analysis of a Pax group III gene from the crustacean Cherax destructor

Pax genes encode transcription factors that are critical regulators of key developmental processes in evolutionarily diverse animal phyla. Here we report the first isolation of a Pax gene from a crustacean: a Pax group III gene we have termed CdpaxIII that contains highly conserved DNA-binding domains, the paired domain and homeodomain. CdpaxIII is expressed in the embryo, in adult limb muscle during both quiescence and regeneration, and during the distinct process of epimorphic limb regeneration. Interestingly, CdpaxIII is expressed as two distinct alternate transcripts, one of which is novel in lacking a large portion of its paired domain.

Co-expression of alternatively spliced forms of PAX3, PAX7, PAX3-FKHR and PAX7-FKHR with distinct DNA binding and transactivation properties in rhabdomyosarcoma

PAX3 and PAX7 encode transcription factors implicated in the pathogenesis of rhabdomyosarcoma (RMS), including alveolar RMS in which chromosomal translocations generate PAX3-FKHR and PAX7-FKHR fusions. Previous studies of wild-type PAX3 and PAX7 identified alternative splicing events that modify the paired box and generate 2 isoforms of PAX3 (Q+ and Q-) and 4 isoforms of PAX7 (Q+GL+, Q+GL-, Q-GL+, Q-GL-). In our study, we investigated alternative splicing of the wild-type and fusion forms of PAX3 and PAX7 in alveolar and embryonal RMS and assessed the functional implications. For PAX3 and PAX3-FKHR, the Q+ and Q- isoforms were consistently co-expressed in RMS tumors with slightly higher levels of the Q+ isoform. For PAX7 and PAX7-FKHR, there was a consistent pattern of co-expression of the 4 isoforms in RMS tumors: Q+GL- > Q+GL+ >/= Q-GL- > Q-GL+. DNA binding analysis demonstrated that PAX3 and PAX3-FKHR Q- isoforms exhibit higher affinity than corresponding Q+ isoforms for class I sites and no difference for class II sites. For PAX7 and PAX7-FKHR, the relative affinity was Q-GL- > Q+GL- > Q-GL+ >/= Q+GL+ for class I sites and Q-GL-, Q+GL- > Q-GL+, Q+GL+ for class II sites. Finally, the transcriptional activities of the PAX3-FKHR and PAX7-FKHR isoforms on reporter plasmids varied over a 5-fold and 50-fold range, respectively, in accord with the differences in DNA binding activity. In conclusion, these studies reveal that PAX3, PAX7 and their fusions with FKHR are each expressed in RMS tumors as a consistent mixture of functionally distinct isoforms.

Zebrafish pax8 is required for otic placode induction and plays a redundant role with Pax2 genes in the maintenance of the otic placode

Vertebrate Pax2 and Pax8 proteins are closely related transcription factors hypothesized to regulate early aspects of inner ear development. In zebrafish and mouse, Pax8 expression is the earliest known marker of otic induction, and Pax2 homologs are expressed at slightly later stages of placodal development. Analysis of compound mutants has not been reported. To facilitate analysis of zebrafish pax8, we completed sequencing of the entire gene, including the 5' and 3' UTRs. pax8 transcripts undergo complex alternative splicing to generate at least ten distinct isoforms. Two different subclasses of pax8 splice isoforms encode different translation initiation sites. Antisense morpholinos (MOs) were designed to block translation from both start sites, and four additional MOs were designed to target different exon-intron boundaries to block splicing. Injection of MOs, individually and in various combinations, generated similar phenotypes. Otic induction was impaired, and otic vesicles were small. Regional ear markers were expressed correctly, but hair cell production was significantly reduced. This phenotype was strongly enhanced by simultaneously disrupting either of the co-inducers fgf3 or fgf8, or another early regulator, dlx3b, which is thought to act in a parallel pathway. In contrast, the phenotype caused by disrupting foxi1, which is required for pax8 expression, was not enhanced by simultaneously disrupting pax8. Disrupting pax8, pax2a and pax2b did not further impair otic induction relative to loss of pax8 alone. However, the amount of otic tissue gradually decreased in pax8-pax2a-pax2b-deficient embryos such that no otic tissue was detectable by 24 hours post-fertilization. Loss of otic tissue did not correlate with increased cell death, suggesting that otic cells dedifferentiate or redifferentiate as other cell type(s). These data show that pax8 is initially required for normal otic induction, and subsequently pax8, pax2a and pax2b act redundantly to maintain otic fate.

Cross-talk between the paired domain and the homeodomain of Pax3: DNA binding by each domain causes a structural change in the other domain, supporting interdependence for DNA Binding

The Pax3 protein has two DNA binding domains, a Paired domain (PD) and a paired-type Homeo domain (HD). Although the PD and HD can bind to cognate DNA sequences when expressed individually, genetic and biochemical data indicate that the two domains are functionally interdependent in intact Pax3. The mechanistic basis of this functional interdependence is unknown and was studied by protease sensitivity. Pax3 was modified by the creation of Factor Xa cleavage sites at discrete locations in the PD, the HD, and in the linker segment joining the PD and the HD (Xa172, Xa189, and Xa216) in individual Pax3 mutants. The effect of Factor Xa insertions on protein stability and on DNA binding by the PD and the HD was measured using specific target site sequences. Independent insertions at position 100 in the linker separating the first from the second helix-turn-helix motif of the PD and at position 216 immediately upstream of the HD were found to be readily accessible to Factor Xa cleavage. The effect of DNA binding by the PD or the HD on accessibility of Factor Xa sites inserted in the same or in the other domain was monitored and quantitated for multiple mutants bearing different numbers of Xa sites at each position. In general, DNA binding reduced accessibility of all sites, suggesting a more compact and less solvent-exposed structure of DNA-bound versus DNA-free Pax3. Results of dose response and time course experiments were consistent and showed that DNA binding by the PD not only caused a local structural change in the PD but also caused a conformational change in the HD (P3OPT binding to Xa216 mutants); similarly, DNA binding by the HD also caused a conformational change in the PD (P2 binding to Xa100 mutants). These results provide a structural basis for the functional interdependence of the two DNA binding domains of Pax3.

Functional interactions between alternatively spliced forms of Pax6 in crystallin gene regulation and in haploinsufficiency

Pax6 is essential for development of the eye, olfactory system, brain and pancreas. Haploinsufficiency of Pax6 causes abnormal eye development. Two forms of Pax6 protein, PAX6 and PAX6(5a), differ in a 14 amino acid insertion encoded by an alternatively spliced exon 5a in the N-terminal DNA-binding paired domain (PD), and they are simultaneously expressed. Here, we show that PAX6 and PAX6(5a) together synergistically activate transcription from promoters recognized by Pax6 PD and PD5a, but not by their homeodomain. This synergism promotes activation of transcription by c-Maf and MafA on the alphaB-crystallin promoter, and is required for transcriptional co-activation by RARbeta/RXRbeta and PAX6/PAX6(5a) on the gammaF-crystallin promoter. To determine the role of this synergism in haploinsufficiency, we tested four human missense (G18W, R26G, G64V and R128C) and one nonsense (R317X) mutants, with reporters driven by Pax6 PD consensus binding sites and the alphaB-crystallin promoter. The simultaneous activity of Pax6 proteins [PAX6, mutated PAX6, PAX6(5a) and mutated PAX6(5a)] modeling haploinsufficiency yielded results not predicted by properties of individual PAX6 or PAX6(5a). Taken together, these results indicate that complex ocular phenotypes due to Pax6 haploinsufficiency originate, at least partially, from functional interactions between alternatively spliced PAX6 and PAX6(5a) variants and other factors, e.g. MafA/c-Maf.

Mutational analysis of the Acropora millepora PaxD paired domain highlights the importance of the linker region for DNA binding

Pax transcription factors are found in animals, from simple sponges to insects and vertebrates. The defining feature of Pax proteins is the DNA-binding paired domain (PD), which consists of two helix-turn-helix subdomains, joined with a linker region. Despite high specificity in vivo, the paired domains of different Pax proteins bind similar consensus DNA sequences in vitro. Using bandshift techniques, we show here that the paired domain of the Acropora millepora PaxD protein, which unambiguously belongs to the Pax3/7 group, does not bind to three defined paired domain-binding sites. Domain swapping experiments and site-directed mutagenesis identified two amino acid residues in the linker region of the paired domain as critical to DNA binding; G70 and S71 are highly conserved in Pax proteins, but differ in PaxD (L70 and N71). The PaxD data thus highlight the importance of the linker region, and particularly G70 and S71, in DNA binding by Pax proteins.