Discovery of an Active RAG Transposon Illuminates the Origins of V(D)J Recombination.

Co-option of RAG1 and RAG2 for antigen receptor gene assembly by V(D)J recombination was a crucial event in the evolution of jawed vertebrate adaptive immunity. RAG1/2 are proposed to have arisen from a transposable element, but definitive evidence for this is lacking. Here, we report the discovery of ProtoRAG, a DNA transposon family from lancelets, the most basal extant chordates. A typical ProtoRAG is flanked by 5-bp target site duplications and a pair of terminal inverted repeats (TIRs) resembling V(D)J recombination signal sequences. Between the TIRs reside tail-to-tail-oriented, intron-containing RAG1-like and RAG2-like genes. We demonstrate that ProtoRAG was recently active in the lancelet germline and that the lancelet RAG1/2-like proteins can mediate TIR-dependent transposon excision, host DNA recombination, transposition, and low-efficiency TIR rejoining using reaction mechanisms similar to those used by vertebrate RAGs. We propose that ProtoRAG represents a molecular "living fossil" of the long-sought RAG transposon.

The IL-17 pathway intertwines with neurotrophin and TLR/IL-1R pathways since its domain shuffling origin.

The IL-17 pathway displays remarkably diverse functional modes between different subphyla, classes, and even orders, yet its driving factors remains elusive. Here, we demonstrate that the IL-17 pathway originated through domain shuffling between a Toll-like receptor (TLR)/IL-1R pathway and a neurotrophin-RTK (receptor-tyrosine-kinase) pathway (a Trunk-Torso pathway). Unlike other new pathways that evolve independently, the IL-17 pathway remains intertwined with its donor pathways throughout later evolution. This intertwining not only influenced the gains and losses of domains and components in the pathway but also drove the diversification of the pathway's functional modes among animal lineages. For instance, we reveal that the crustacean female sex hormone, a neurotrophin inducing sex differentiation, could interact with IL-17Rs and thus be classified as true IL-17s. Additionally, the insect prothoracicotropic hormone, a neurotrophin initiating ecdysis in Drosophila by binding to Torso, could bind to IL-17Rs in other insects. Furthermore, IL-17R and TLR/IL-1R pathways maintain crosstalk in amphioxus and zebrafish. Moreover, the loss of the Death domain in the pathway adaptor connection to IκB kinase and stress-activated protein kinase (CIKSs) dramatically reduced their abilities to activate nuclear factor-kappaB (NF-κB) and activator protein 1 (AP-1) in amphioxus and zebrafish. Reinstating this Death domain not only enhanced NF-κB/AP-1 activation but also strengthened anti-bacterial immunity in zebrafish larvae. This could explain why the mammalian IL-17 pathway, whose CIKS also lacks Death, is considered a weak signaling activator, relying on synergies with other pathways. Our findings provide insights into the functional diversity of the IL-17 pathway and unveil evolutionary principles that could govern the pathway and be used to redesign and manipulate it.

Both IRAK3 and IRAK1 Activate the MyD88-TRAF6 Pathway in Zebrafish.

IL-1R-associated kinases (IRAKs) are signal transducers of the TLR/IL-1R-MyD88-TRAF6 pathways. Vertebrates possess two IRAK lineages, IRAK1/2/3 and IRAK4. In mammals, IRAK4/IRAK1 and IRAK4/IRAK2 are pathway enhancers, whereas IRAK3 is a repressor. However, in bony fish, IRAK2 is absent, and it remains elusive how fish IRAK1/3/4 functionally differ from their mammalian counterparts. In this study, we explored this using the zebrafish model. First, we showed that in human 293T cells, zebrafish IRAK1 and IRAK4 were components of the Myddosome (MyD88-IRAK4-IRAK1) complex, with IRAK1 serving as a potent pathway enhancer. Then, we discovered two zebrafish IRAK3 variants: one (IRAK3a) contains an N-terminal Death domain, a middle pseudokinase domain, and a C-terminal TRAF6-binding domain, whereas the other (IRAK3b) lost both the kinase and TRAF6-binding domains. This truncation of IRAK3 variants could be a conserved phenomenon in fish, because it is also observed in trout and grass carp. We proceeded to show that zebrafish IRAK3a acts as a pathway enhancer by binding with MyD88 and TRAF6, but its activity is milder than IRAK1, possibly because it has no kinase activity. Zebrafish IRAK3b, however, plays a sheer negative role, apparently because of its lack of kinase and TRAF6-binding domains. Moreover, zebrafish IRAK3a/3b inhibit the activity of IRAK1/4, not by interacting with IRAK1/4 but possibly by competing for MyD88 and TRAF6. Finally, we have verified the essential activities of zebrafish IRAK1/3a/3b/4 in zebrafish cells and embryos. In summary, to our knowledge, our findings provide new insights into the molecular functions of fish IRAKs and the evolution of the IRAK functional modes in vertebrates.

Differential roles of the fish chitinous membrane in gut barrier immunity and digestive compartments.

The chitin-based peritrophic matrix (PM) is a structure critical for both gut immunity and digestion in invertebrates. PM was traditionally considered lost in all vertebrates, but a PM-like chitinous membrane (CM) has recently been discovered in fishes, which may increase the knowledge on vertebrate gut physiology and structural evolution. Here, we show that in zebrafish, the CM affects ingestion behavior, microbial homeostasis, epithelial renewal, digestion, growth, and longevity. Young mutant fish without CM appear healthy and are able to complete their life cycle normally, but with increasing age they develop gut inflammation, resulting in gut atrophy. Unlike mammals, zebrafish have no visible gel-forming mucin layers to protect their gut epithelia, but at least in young fish, the CM is not a prerequisite for the antibacterial gut immunity. These findings provide new insights into the role of the CM in fish prosperity and its eventual loss in tetrapods. These findings may also help to improve fish health and conservation, as well as to advance the understanding of vertebrate gut physiology and human intestinal diseases.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

A chimeric hairpin DNA aptamer-based biosensor for monitoring the therapeutic drug bevacizumab.

The accurate and rapid detection of specific antibodies in blood is very important for efficient diagnosis and precise treatment. Conventional methods often suffer from time-consuming operations and/or a narrow detection range. In this work, for the rapid determination of bevacizumab in plasma, a series of chimeric hairpin DNA aptamer-based probes were designed by the modification, labeling and theoretical computation of an original aptamer. Then, the dissociation constant of the modified hairpin DNA to bevacizumab was measured and screened using microscale thermophoresis. The best chimeric hairpin DNA aptamer-based probe was then selected, and a one-step platform for the rapid determination of bevacizumab was constructed. This strategy has the advantages of being simple, fast and label-free. Because of the design and screening of the hairpin DNA, as well as the optimization of the concentration and electrochemical parameters, a low detection limit of 0.37 pM (0.054 ng mL-1) with a wide linear range (1 pM-1 µM) was obtained. Finally, the rationally constructed biosensor was successfully applied to the determination of bevacizumab in spiked samples, and it showed good accuracy and precision. This method is expected to truly realize accurate and rapid detection of bevacizumab and provides a new idea for the precise treatment of diseases.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Functional Variation of IL-1R-Associated Kinases in the Conserved MyD88-TRAF6 Pathway during Evolution.

IL-1R-associated kinases (IRAK) are important regulators in the TLR/IL-1R pathways, but their function appears inconsistent between Drosophila, bony fishes, and vertebrates. This causes a difficulty to understand the IRAK functions. As a step to reveal the evolution of IRAKs, in this study, we performed comparative and functional analysis of IRAKs by exploiting the amphioxus, a pivotal taxon connecting invertebrates and vertebrates. Sequence and phylogenetic analysis indicated three major IRAK lineages: IRAK1/2/3 is a vertebrate-specific lineage, IRAK4 is an ancient lineage conserved between invertebrate and vertebrates, and Pelle is another ancient lineage that is preserved in protostomes and invertebrate deuterostomes but lost in vertebrate deuterostomes. Pelle is closer neither to IRAK4 nor to IRAK1/2/3, hence suggesting no clear functional analogs to IRAK1/2/3 in nonvertebrates. Functional analysis showed that both amphioxus IRAK4 and Pelle could suppress NF-κB activation induced by MyD88 and TRAF6, which are unlike mammalian and Drosophila IRAKs, but, surprisingly, similar to bony fish IRAK4. Also unlike Drosophila IRAKs, no interaction was detected between amphioxus IRAK4 and Pelle, although both of them were shown capable of binding MyD88. These findings, together with previous reports, show that unlike other signal transducers in the TLR/IL-1R pathways, such as MyD88 and TRAF6, the functions of IRAKs are highly variable during evolution and very specialized in different major animal taxa. Indeed, we suggest that the functional variability of IRAKs might confer plasticity to the signal transduction of the TLR/IL-1R pathways, which in return helps the species to evolve against the pathogens.

Broad distribution, high diversity and ancient origin of the ApeC-containing proteins.

Apextrin C-terminal (ApeC) is a novel protein domain with unknown functions, although early studies suggest that some ApeC-containing proteins (ACPs) bind to carbohydrates and have a role in development and immunity. Here we investigated the taxonomic distribution, sequence diversification and origin of ACPs in Metazoa. Most ACPs are present in invertebrates from aquatic or moist environments, including cnidarians, mollusks, echinoderms, cephalochordates, flatworms, water bears, nematodes and annelids. However, ACPs are absent in vertebrates and in most arthropod lineages (e.g. insects and crustaceans) except arachnids. ACPs apparently undergo rapid turnover and diversification, hence no orthologs could be found between (sub)phyla. ApeC can function either as a standalone domain or as a partner domain. It has been found to pair up with over ten different domain types in different ACPs. The partner domains are related to immunity, extracellular matrix, protein-protein and protein-carbohydrate interactions. Notably, the domain pair with the widest taxonomic distribution is MACPF/perforin-ApeC, which represent a classic group of ACPs called apextrins. ApeC also frequently pairs up with itself to form dual-ApeC modules in different phyla. Notably, in parasite flatworms, dual-ApeCs are present in 70% of ACPs and all inherited from a common ancestor. The broad distribution of MACPF-ApeC and dual-ApeC suggest their conserved yet unknown functions. We also discovered distant ApeC homologs in bacteria, hence tracing the origin of ApeC back to prokaryotes. Our findings show that ApeC has an ancient origin and is able to function alone or in complex domain architectures, though it is less prevalent than other versatile domains such as immunoglobulin domains and C-type lectin domains. This work provides a foundation for further functional study of this novel domain type.

The family of amphioxus chitin synthases offers insight into the evolution of chitin formation in chordates.

Chitin is a very important and widely-used biopolymer in fungi and lower metazoans, but mysteriously disappears in mammals. Recent studies reveal that at least lower vertebrates have chitin synthases (CS) and use them to synthesize endogenous chitin. Amphioxus, a basal chordate, therefore becomes critical to understand the evolution of CS, as it occupies the transitional position from invertebrates to vertebrates, and is considered as a good proxy to the chordate ancestor. Here, by exploiting multiple genome assemblies, high-depth RNA-seq data and synteny relations, we identify 11-12 CS genes for each amphioxus species. It represents the largest CS gene pool ever found in eukaryotes so far. As comparison, most metazoans have one or two CSs. Amphioxus is the only chordate that has both the very ancient type-I CS family and the more broadly distributed type-II CS family. Specifically, amphioxus has only one type-II CS but 10-11 type-I CSs, which means that amphioxus is the only metazoan with a greatly expanded type-I CS family. Further analysis suggests that the chordate ancestor have at least one type-II CS and an expanded of type-I CS family. We hypothesize that: these ancient CSs are mostly retained in amphioxus; but the whole type-I CS family was lost in urochordates and vertebrates; the type-II CS was later duplicated into two lineages in vertebrates and followed by stochastic losses, till all type-II CSs were eventually lost in birds and mammals. Finally, our expression profiling and preliminary gene knockout analysis suggest that amphioxus CSs could have highly diverse but mildly overlapping functions in various tissues and organs. Taken together, these findings not only provide insights into the evolution of chordate CSs, lay a foundation for further functional study of the chordate CSs. After all, it is mysterious that our chordate ancestor needed so many isoenzymes for chitin formation.

New insights on unspecific peroxygenases: superfamily reclassification and evolution.

BACKGROUND: Unspecific peroxygenases (UPO) (EC 1.11.2.1) represent an intriguing oxidoreductase sub-subclass of heme proteins with peroxygenase and peroxidase activity. With over 300 identified substrates, UPOs catalyze numerous oxidations including 1- or 2- electron oxygenation, selective oxyfunctionalizations, which make them most significant in organic syntheses and potentially attractive as industrial biocatalysts. There are very few UPOs available with distinct properties, notably, MroUPO which shows behavior ranging between UPO and another heme-thiolate peroxidase, called Chloroperoxidase (CPO). It prompted us to search for more UPOs in fungal kingdom which led us to studying their relationship with CPO. RESULTS: In this study, we searched for novel UPOs in more than 800 fungal genomes and found 113 putative UPO-encoding sequences distributed in 35 different fungal species (or strains), amongst which single sequence per species were subjected to phylogeny study along with CPOs. Our phylogenetic study show that the UPOs are distributed in Basidiomycota and Ascomycota phyla of fungi. The sequence analysis helped to classify the UPOs into five distinct subfamilies: classic AaeUPO and four new subfamilies with potential new traits. We have also shown that each of these five subfamilies (supported by) have their own signature motifs. Surprisingly, some of the CPOs appeared to be a type of UPOs indicating that they were previously identified incorrectly. Selection pressure was observed on important motifs in UPOs which could have driven their functional divergence. Furthermore, the sites having different evolutionary rates caused by the functional divergence were also identified on some motifs along with the other relevant amino acid residues. Finally, we predicted critical amino acids responsible for the functional divergence in the UPOs and identified some sequence differences among UPOs, CPOs, and MroUPO to predict it's ranging behavior. CONCLUSION: This study discovers new UPOs, provides a glimpse of their evolution from CPOs, and presents new insight on their functional divergence. We present a new classification of UPOs and shed new light on its phylogenetics. These different UPOs may exhibit a wide range of characteristics and specificities which may help in various fields of synthetic chemistry and industrial biocatalysts, and may as well lead to an advancement towards the understanding of physiological role of UPOs in fungi.

HaploMerger2: rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly.

SUMMARY: De novo assembly is a difficult issue for heterozygous diploid genomes. The advent of high-throughput short-read and long-read sequencing technologies provides both new challenges and potential solutions to the issue. Here, we present HaploMerger2 (HM2), an automated pipeline for rebuilding both haploid sub-assemblies from the polymorphic diploid genome assembly. It is designed to work on pre-existing diploid assemblies, which are typically created by using de novo assemblers. HM2 can process any diploid assemblies, but it is especially suitable for diploid assemblies with high heterozygosity (≥3%), which can be difficult for other tools. This pipeline also implements flexible and sensitive assembly error detection, a hierarchical scaffolding procedure and a reliable gap-closing method for haploid sub-assemblies. Using HM2, we demonstrate that two haploid sub-assemblies reconstructed from a real, highly-polymorphic diploid assembly show greatly improved continuity. AVAILABILITY AND IMPLEMENTATION: Source code, executables and the testing dataset are freely available at https://github.com/mapleforest/HaploMerger2/releases/. CONTACT: hshengf2@mail.sysu.edu.cn. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic Regulation of Tandem 3' Untranslated Regions in Zebrafish Spleen Cells during Immune Response.

Alternative polyadenylation (APA) has been found to be involved in tumorigenesis, development, and cell differentiation, as well as in the activation of several subsets of immune cells in vitro. Whether APA takes place in immune responses in vivo is largely unknown. We profiled the variation in tandem 3' untranslated regions (UTRs) in pathogen-challenged zebrafish and identified hundreds of APA genes with ∼ 10% being immune response genes. The detected immune response APA genes were enriched in TLR signaling, apoptosis, and JAK-STAT signaling pathways. A greater number of microRNA target sites and AU-rich elements were found in the extended 3' UTRs than in the common 3' UTRs of these APA genes. Further analysis suggested that microRNA and AU-rich element-mediated posttranscriptional regulation plays an important role in modulating the expression of APA genes. These results indicate that APA is extensively involved in immune responses in vivo, and it may be a potential new paradigm for immune regulation.

Characterization of Amphioxus IFN Regulatory Factor Family Reveals an Archaic Signaling Framework for Innate Immune Response.

The IFN regulatory factor (IRF) family encodes transcription factors that play important roles in immune defense, stress response, reproduction, development, and carcinogenesis. Although the origin of the IRF family has been dated back to multicellular organisms, invertebrate IRFs differ from vertebrate IRFs in genomic structure and gene synteny, and little is known about their functions. Through comparison of multiple amphioxus genomes, in this study we suggested that amphioxus contains nine IRF members, whose orthologs are supposed to be shared among three amphioxus species. As the orthologs to the vertebrate IRF1 and IRF4 subgroups, Branchiostoma belcheri tsingtauense (bbt)IRF1 and bbtIRF8 bind the IFN-stimulated response element (ISRE) and were upregulated when amphioxus intestinal cells were stimulated with poly(I:C). As amphioxus-specific IRFs, both bbtIRF3 and bbtIRF7 bind ISRE. When activated, they can be phosphorylated by bbtTBK1 and then translocate into nucleus for target gene transcription. As transcriptional repressors, bbtIRF2 and bbtIRF4 can inhibit the transcriptional activities of bbtIRF1, 3, 7, and 8 by competing for the binding of ISRE. Interestingly, amphioxus IRF2, IRF8, and Rel were identified as target genes of bbtIRF1, bbtIRF7, and bbtIRF3, respectively, suggesting a dynamic feedback regulation among amphioxus IRF and NF-κB. Collectively, to our knowledge we present for the first time an archaic IRF signaling framework in a basal chordate, shedding new insights into the origin and evolution of vertebrate IFN-based antiviral networks.

APASdb: a database describing alternative poly(A) sites and selection of heterogeneous cleavage sites downstream of poly(A) signals.

Increasing amounts of genes have been shown to utilize alternative polyadenylation (APA) 3'-processing sites depending on the cell and tissue type and/or physiological and pathological conditions at the time of processing, and the construction of genome-wide database regarding APA is urgently needed for better understanding poly(A) site selection and APA-directed gene expression regulation for a given biology. Here we present a web-accessible database, named APASdb (http://mosas.sysu.edu.cn/utr), which can visualize the precise map and usage quantification of different APA isoforms for all genes. The datasets are deeply profiled by the sequencing alternative polyadenylation sites (SAPAS) method capable of high-throughput sequencing 3'-ends of polyadenylated transcripts. Thus, APASdb details all the heterogeneous cleavage sites downstream of poly(A) signals, and maintains near complete coverage for APA sites, much better than the previous databases using conventional methods. Furthermore, APASdb provides the quantification of a given APA variant among transcripts with different APA sites by computing their corresponding normalized-reads, making our database more useful. In addition, APASdb supports URL-based retrieval, browsing and display of exon-intron structure, poly(A) signals, poly(A) sites location and usage reads, and 3'-untranslated regions (3'-UTRs). Currently, APASdb involves APA in various biological processes and diseases in human, mouse and zebrafish.

Two apextrin-like proteins mediate extracellular and intracellular bacterial recognition in amphioxus.

Animals exploit different germ-line-encoded proteins with various domain structures to detect the signature molecules of pathogenic microbes. These molecules are known as pathogen-associated molecular patterns (PAMPs), and the host proteins that react with PAMPs are called pattern recognition proteins (PRPs). Here, we present a novel type of protein domain structure capable of binding to bacterial peptidoglycan (PGN) and the minimal PGN motif muramyl dipeptide (MDP). This domain is designated as apextrin C-terminal domain (ApeC), and its presence was confirmed in several invertebrate phyla and subphyla. Two apextrin-like proteins (ALP1 and ALP2) were identified in a basal chordate, the Japanese amphioxus Branchiostoma japonicum (bj). bjALP1 is a mucosal effector secreted into the gut lumen to agglutinate the Gram-positive bacterium Staphylococcus aureus via PGN binding. Neutralization of secreted bjALP1 by anti-bjALP1 monoclonal antibodies caused serious damage to the gut epithelium and rapid death of the animals after bacterial infection. bjALP2 is an intracellular PGN sensor that binds to TNF receptor-associated factor 6 (TRAF6) and prevents TRAF6 from self-ubiquitination and hence from NF-κB activation. MDP was found to compete with TRAF6 for bjALP2, which released TRAF6 to activate the NF-κB pathway. BjALP1 and bjALP2 therefore play distinct and complementary functions in amphioxus gut mucosal immunity. In conclusion, discovery of the ApeC domain and the functional analyses of amphioxus ALP1 and ALP2 allowed us to define a previously undocumented type of PRP that is represented across different animal phyla.

Association of OCT4, SOX2, and NANOG expression with oral squamous cell carcinoma progression.

BACKGROUND: OCT4, SOX2, and NANOG are major transcription factors related to stem cell self-renewal and differentiation. The aim of this study was to examine the association of OCT4, SOX2, and NANOG expression levels with the development and prognosis of patients with oral squamous cell carcinoma (OSCC). MATERIALS AND METHODS: Expression levels of OCT4, SOX2, and NANOG were evaluated by immunohistochemistry with tissue microarray slides of 436 OSCC, 362 corresponding tumor-adjacent normal (CTAN) tissues, and 71 normal uvula epithelium tissues. The clinicopathologic and follow-up data of the OSCC patients were recorded. RESULTS: OCT4 expression was significantly higher in normal and CTAN tissues than in tumor tissue (both P < 0.001). SOX2 expression in CTAN tissue was significantly higher than that in normal (P = 0.021) and tumor tissues (P < 0.001). However, NANOG expression was significantly higher in CTAN (P = 0.014) and tumor tissues (P = 0.009) than in normal tissue. Higher OCT4 and SOX2 expressions were associated with earlier AJCC stage (P = 0.002 and P < 0.001), small tumor size (P = 0.017 and P = 0.001), and the absence of lymph node metastasis (P = 0.015 and P = 0.025). Higher levels of SOX2 expression were associated with better disease-specific survival (P = 0.002) even after adjustment for clinicopathologic factors. DISCUSSION: OCT4 and SOX2 are biomarkers of tumorigenesis and early stage OSCC. SOX2 is an independent prognostic factor for OSCC.

Ancestral genetic complexity of arachidonic acid metabolism in Metazoa.

Eicosanoids play an important role in inducing complex and crucial physiological processes in animals. Eicosanoid biosynthesis in animals is widely reported; however, eicosanoid production in invertebrate tissue is remarkably different to vertebrates and in certain respects remains elusive. We, for the first time, compared the orthologs involved in arachidonic acid (AA) metabolism in 14 species of invertebrates and 3 species of vertebrates. Based on parsimony, a complex AA-metabolic system may have existed in the common ancestor of the Metazoa, and then expanded and diversified through invertebrate lineages. A primary vertebrate-like AA-metabolic system via cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome P450 (CYP) pathways was further identified in the basal chordate, amphioxus. The expression profiling of AA-metabolic enzymes and lipidomic analysis of eicosanoid production in the tissues of amphioxus supported our supposition. Thus, we proposed that the ancestral complexity of AA-metabolic network diversified with the different lineages of invertebrates, adapting with the diversity of body plans and ecological opportunity, and arriving at the vertebrate-like pattern in the basal chordate, amphioxus.

Dynamic landscape of tandem 3' UTRs during zebrafish development.

Tandem 3' untranslated regions (UTRs), produced by alternative polyadenylation (APA) in the terminal exon of a gene, could have critical roles in regulating gene networks. Here we profiled tandem poly(A) events on a genome-wide scale during the embryonic development of zebrafish (Danio rerio) using a recently developed SAPAS method. We showed that 43% of the expressed protein-coding genes have tandem 3' UTRs. The average 3' UTR length follows a V-shaped dynamic pattern during early embryogenesis, in which the 3' UTRs are first shortened at zygotic genome activation, and then quickly lengthened during gastrulation. Over 4000 genes are found to switch tandem APA sites, and the distinct functional roles of these genes are indicated by Gene Ontology analysis. Three families of cis-elements, including miR-430 seed, U-rich element, and canonical poly(A) signal, are enriched in 3' UTR-shortened/lengthened genes in a stage-specific manner, suggesting temporal regulation coordinated by APA and trans-acting factors. Our results highlight the regulatory role of tandem 3' UTR control in early embryogenesis and suggest that APA may represent a new epigenetic paradigm of physiological regulations.

HaploMerger: reconstructing allelic relationships for polymorphic diploid genome assemblies.

Whole-genome shotgun assembly has been a long-standing issue for highly polymorphic genomes, and the advent of next-generation sequencing technologies has made the issue more challenging than ever. Here we present an automated pipeline, HaploMerger, for reconstructing allelic relationships in a diploid assembly. HaploMerger combines a LASTZ-ChainNet alignment approach with a novel graph-based structure, which helps to untangle allelic relationships between two haplotypes and guides the subsequent creation of reference haploid assemblies. The pipeline provides flexible parameters and schemes to improve the contiguity, continuity, and completeness of the reference assemblies. We show that HaploMerger produces efficient and accurate results in simulations and has advantages over manual curation when applied to real polymorphic assemblies (e.g., 4%-5% heterozygosity). We also used HaploMerger to analyze the diploid assembly of a single Chinese amphioxus (Branchiostoma belcheri) and compared the resulting haploid assemblies with EST sequences, which revealed that the two haplotypes are not only divergent but also highly complementary to each other. Taken together, we have demonstrated that HaploMerger is an effective tool for analyzing and exploiting polymorphic genome assemblies.