Exploring gene expression changes in the amphioxus gill after poly(I:C) challenge using digital expression profiling.

Amphioxus, a cephalochordate, is a key model animal for studying the evolution of vertebrate immunity. Recently, studies have revealed that microRNA (miRNA) expression profiles change significantly in the amphioxus gill after immune stimulation, but it remains largely unknown how gene expression responds to immune stress. Elucidating gene expression changes in the amphioxus gill will provide a deeper understanding of the evolution of gill immunity in vertebrates. Here, we used high-throughput RNA sequencing technology (RNA-seq) to conduct tag-based digital gene expression profiling (DGE) analyses of the gills of control Branchiostoma belcheri and of those exposed to the viral mimic, poly(I:C) (pIC). Six libraries were created for the control and treatment groups including three biological replicates per group. A total of 1999 differently expressed genes (DEGs) were obtained, with 571 and 1428 DEGs showing up- or down-regulation, respectively, in the treatment group. Enrichment analysis of gene ontology (GO) terms and pathways revealed that the DEGs were primarily related to immune and defense response, apoptosis, human disease, cancer, protein metabolism, enzyme activity, and regulatory processes. In addition, eight DEGs were randomly selected to validate the RNA-seq data using real-time quantitative PCR (qRT-PCR), and the results confirmed the accuracy of the RNA-seq approach. Next, we screened eight key responding genes to examine the dynamic changes in expression levels at different time points in more detail. The results indicated that expressions of TRADD, MARCH, RNF31, NF-κb, CYP450, TNFRSF6B, IFI and LECT1 were induced to participate in the antiviral response against pIC. This study provides a valuable resource for understanding the role of the amphioxus gill in antiviral immunity and the evolution of gill immunity in vertebrates.

Genome-wide gene expression analysis of amphioxus (Branchiostoma belcheri) following lipopolysaccharide challenge using strand-specific RNA-seq.

Amphioxus is the closest living proxy for exploring the evolutionary origin of the immune system in vertebrates. To understand the immune responses of amphioxus to lipopolysaccharide (LPS), 5 ribosomal RNA (rRNA)-depleted libraries of amphioxus were constructed, including one control (0 h) library and 4 treatment libraries at 6, 12, 24, and 48 h post-injection (hpi) with LPS. The transcriptome of Branchiostoma belcheri was analyzed using strand-specific RNA sequencing technology (RNA-seq). A total of 6161, 6665, 7969, and 6447 differentially expressed genes (DEGs) were detected at 6, 12, 24, and 48 hpi, respectively, compared with expression levels at 0 h. We identified amphioxus genes active during the acute-phase response to LPS at different time points after stimulation. Moreover, to better visualize the resolution phase of the immune process during immune response, we identified 6057 and 5235 DEGs at 48 hpi by comparing with 6 and 24 hpi, respectively. Through real-time quantitative PCR (qRT-PCR) analysis of 12 selected DEGs, we demonstrated the accuracy of the RNA-seq data in this study. Functional enrichment analysis of DEGs demonstrated that most terms were related to defense and immune responses, disease and infection, cell apoptosis, and metabolism and catalysis. Subsequently, we identified 1330, 485, 670, 911, and 1624 time-specific genes (TSGs) at 0, 6, 12, 24, and 48 hpi. Time-specific terms at each of 5 time points were primarily involved in development, immune signaling, signal transduction, DNA repair and stability, and metabolism and catalysis, respectively. As this is the first study to report the transcriptome of an organism with primitive immunity following LPS challenge at multiple time points, it provides gene expression information for further research into the evolution of immunity in vertebrates.

Genome-wide comparison of the protein-coding repertoire reveals fast evolution of immune-related genes in cephalochordates and Osteichthyes superclass.

Amphioxus is used to investigate the origin and evolution of vertebrates. To better understand the characteristics of genome evolution from cephalochordates to Osteichthyes, we conducted a genome-wide pairwise comparison of protein-coding genes within amphioxus (a comparable group) and parallel analyses within Osteichthyes (two comparable groups). A batch of fast-evolving genes in each comparable group was identified. Of these genes, the most fast-evolving genes (top 20) were scrutinized, most of which were involved in immune system. An analysis of the fast-evolving genes showed that they were enriched into gene ontology (GO) terms and pathways primarily involved in immune-related functions. Similarly, this phenomenon was detected within Osteichthyes, and more well-known and abundant GO terms and pathways involving innate immunity were found in Osteichthyes than in cephalochordates. Next, we measured the expression responses of four genes belonging to metabolism or energy production-related pathways to lipopolysaccharide challenge in the muscle, intestine or skin of B. belcheri; three of these genes (HMGCL, CYBS and MDH2) showed innate immune responses. Additionally, some genes involved in adaptive immunity showed fast evolution in Osteichthyes, such as those involving "intestinal immune network for IgA production" or "T-cell receptor signaling pathway". In this study, the fast evolution of immune-related genes in amphioxus and Osteichthyes was determined, providing insights into the evolution of immune-related genes in chordates.

Transcriptome-wide analysis of immune-responsive microRNAs against poly (I:C) challenge in Branchiostoma belcheri by deep sequencing and bioinformatics.

Amphioxus is a key experimental animal for studying the evolution of vertebrate immune system. However, we still do not know about the roles of microRNAs (miRNAs) under viral stress in amphioxus. In this study, we sequenced six small RNA libraries (three biological replicates were included in the treatments challenged by the viral mimic, poly (I:C) (pIC) and control groups, respectively) from Branchiostoma belcheri. A total of 151 known miRNAs, 197 new miRNAs (named novel_mir, including nine conserved miRNAs) were identified by deep sequencing from the six libraries. We primarily focused on differentially expressed miRNAs (DEMs) after pIC challenge. Next, we screened a total of 77 DEMs, including 27 down- and 50 up-regulated DEMs in response to pIC challenge. Furthermore, we used real-time quantitative PCR (qRT-PCR) to verify the expression levels of 10 randomly selected DEMs. Target genes likely regulated by DEMs were predicted, and functional enrichment analyses of these targets were performed using bioinformatics approach. MiRNA targets of DEMs are primarily involved in immune response, diseases, cancer and regulation process, and could be largely linked to 14 immune-related signaling pathways, including NF-kappa B, NOD-like receptor, RIG-I-like receptor and endocytosis. The present study for the first time explores key regulatory roles of miRNAs in the innate antiviral immune response in amphioxus, and will provide insight into the molecular basis of antiviral immunity and evolution of immune-related miRNAs.

Genome-wide gene expression analysis in the amphioxus, Branchiostoma belcheri after poly (I: C) challenge using strand-specific RNA-seq.

The gene expression associated with immune response to bacteria/bacterial mimic has been extensively analyzed in amphioxus, but remains largely unknown about how gene are involved in the immune response to viral invasion at expression level. Here, we analyze the rRNA-depleted transcriptomes of Branchiostoma belcheri using strand-specific RNA-seq in response to the viral mimic, poly (I:C) (pIC). A total of 5,317 differentially expressed genes were detected at treatment group by comparing with control. The gene with the most significant expression changes (top 15) after pIC challenge and 7 immune-related categories involving 58 differently expressed genes were scrutinized. By functional enrichment analysis of differently expressed genes, gene ontology terms involving response to stress and stimulus, apoptosis, catabolic and metabolic processes and enzyme activity were overrepresented, and several pathways related to immune signaling, immune response, cancer, apoptosis, viral disease, metabolism were activated after pIC injection. A positive correlation between the qRT-PCR and strand-specific RNA-seq data confirmed the accuracy of the RNA-seq results. Additionally, the expression of genes encoding NLRC5, CASP1, CASP6, CYP450, CAT, and MDA5 were induced in B. belcheri under pIC challenge. Our experiments provide insight into the immune response of amphioxus to pIC and valuable gene expression information for studying the evolution of antiviral immunity in vertebrates.

Selection of reliable reference genes for normalization of quantitative RT-PCR from different developmental stages and tissues in amphioxus.

Amphioxus is a closest living proxy to the ancestor of cephalochordates with vertebrates, and key animal for novel understanding in the evolutionary origin of vertebrate body plan, genome, tissues and immune system. Reliable analyses using quantitative real-time PCR (qRT-PCR) for answering these scientific questions is heavily dependent on reliable reference genes (RGs). In this study, we evaluated stability of thirteen candidate RGs in qRT-PCR for different developmental stages and tissues of amphioxus by four independent (geNorm, NormFinder, BestKeeper and deltaCt) and one comparative algorithms (RefFinder). The results showed that the top two stable RGs were the following: (1) S20 and 18 S in thirteen developmental stages, (2) EF1A and ACT in seven normal tissues, (3) S20 and L13 in both intestine and hepatic caecum challenged with lipopolysaccharide (LPS), and (4) S20 and EF1A in gill challenged with LPS. The expression profiles of two target genes (EYA and HHEX) in thirteen developmental stages were used to confirm the reliability of chosen RGs. This study identified optimal RGs that can be used to accurately measure gene expression under these conditions, which will benefit evolutionary and functional genomics studies in amphioxus.

[Mapping of five genes resistant to SMV strains in soybean].

Soybean mosaic virus (SMV) is one of the most prevalent pathogens that impact the soybean world widely. Previous reports showed that most of the resistances were controlled by one pair of dominant genes. In this study, a soybean RIL population NJRIKY derived from Kefeng 1 x Nannong 1138-2 was used to study the inheritance of resistance to five SMV strains (Sa,Sc-8, Sc-9, N1, and N3) and mapping of resistant genes. Kefeng 1 is resistant to all five SMV strains while Nannong 1138-2 is susceptible to all five SMV strains. The parents and RIL populations were planted in green house and five SMV strains were inoculated on different populations. The results showed that each ratio of the number of resistant families to that of susceptible families was consistent with 1:1 for the five strains. This indicated that the resistance to each of the five strains was controlled by one dominant gene, respectively. RFLP and SSR markers were used to analyze the RIL population, Mapmaker/Exp 3.0b was used to study the linkage between markers and the resistant genes. Through linkage analysis, Rsa was found linked Rn1, Rn3 and Rsc9 with 21.4 cM, 23.5 cM and 35.3 cM, Rsc8 was found to be linked only Rn1 with 35.8 cM. Multi locus analysis showed that the order and intervals of the five resistance genes were Rsc8-35.8 cM-Rn1-10.3 cM-Rn3-21.5 cM- Rsa-35.8 cM-Rsc9. According to the result of RFLP and SSR analysis, a genetic map was constructed which consisted of 256 markers that covered 3050.9 cM and converged into 22 linkage groups. The five resistance genes were mapped on the linkage group N8-D1b + W. The RFLP markers A691T, K4771, LC5T were found linked to the resistant genes Rn1 and Rn3 with distances of 15.04 cM, 17.82 cM, 15.37 cM, 16.14 cM, 17.82 cM, 16.58 cM.

Raman spectra of a Lower Cambrian ctenophore embryo from southwestern Shaanxi, China.

The Early Cambrian (approximately 540 million years old) Meishucun fossil assemblage of Ningqiang County (Shaanxi Province), China, contains the oldest complex skeletonized organisms known in the geological record. We here report the finding in this assemblage of an exquisitely preserved late-stage embryo of a ctenophore ("comb jelly"), its fine structure documented by confocal laser scanning microscopy and shown by Raman spectroscopy to be composed of carbonaceous kerogen permineralized in apatite. In its spheroidal morphology, the presence of eight comb rows and the absence of tentacles, this embryo resembles an adult ctenophore (Maotianoascus octonarius) known from the immediately younger Chengjiang fauna of Yunnan, China. The oldest ctenophore and the only embryonic comb jelly known from the fossil record, this exceptionally well preserved specimen provides important clues about the early evolution of the phylum Ctenophora and of metazoans in general.

Phosphatized polar lobe-forming embryos from the Precambrian of southwest China.

In developing embryos of some extant spiralian animals, polar lobe formation is one of the symmetry-breaking mechanisms for segregation of maternal cytoplasmic substances to certain blastomeres and not others. Polar lobe formation leads to unique early cleavage morphologies that include trilobed, J-shaped, and five-lobed structures. Fossil embryos similar to modern lobeforming embryos are recognized from the Precambrian Doushantuo Formation phosphates, Weng'an, Guizhou Province, China. These embryos are abundant and form a developmental sequence comparable to different developing stages observed in lobe-forming embryos of extant spiralians. These data imply that lobe formation is an evolutionarily ancient process of embryonic specification.

The first tunicate from the Early Cambrian of South China.

Here we report the discovery of eight specimens of an Early Cambrian fossil tunicate Shankouclava near Kunming (South China). The tunicate identity of this organism is supported by the presence of a large and perforated branchial basket, a sac-like peri-pharyngeal atrium, an oral siphon with apparent oral tentacles at the basal end of the siphonal chamber, perhaps a dorsal atrial pore, and an elongated endostyle on the mid-ventral floor of the pharynx. As in most modern tunicates, the gut is simple and U-shaped, and is connected with posterior end of the pharynx at one end and with an atrial siphon at the other, anal end. Shankouclava differs from Cheungkongella, which was previously called a tunicate. Based on new, more complete "Cheungkongella" specimens that show branching tentacles, this form may be a lophophorate, and in any case is not a tunicate.