Unraveling incompatibility between wheat and the fungal pathogen Zymoseptoria tritici through apoplastic proteomics.

BACKGROUND: Hemibiotrophic fungal pathogen Zymoseptoria tritici causes severe foliar disease in wheat. However, current knowledge of molecular mechanisms involved in plant resistance to Z. tritici and Z. tritici virulence factors is far from being complete. The present work investigated the proteome of leaf apoplastic fluid with emphasis on both host wheat and Z. tritici during the compatible and incompatible interactions. RESULTS: The proteomics analysis revealed rapid host responses to the biotrophic growth, including enhanced carbohydrate metabolism, apoplastic defenses and stress, and cell wall reinforcement, might contribute to resistance. Compatibility between the host and the pathogen was associated with inactivated plant apoplastic responses as well as fungal defenses to oxidative stress and perturbation of plant cell wall during the initial biotrophic stage, followed by the strong induction of plant defenses during the necrotrophic stage. To study the role of anti-oxidative stress in Z. tritici pathogenicity in depth, a YAP1 transcription factor regulating antioxidant expression was deleted and showed the contribution to anti-oxidative stress in Z. tritici, but was not required for pathogenicity. This result suggests the functional redundancy of antioxidants in the fungus. CONCLUSIONS: The data demonstrate that incompatibility is probably resulted from the proteome-level activation of host apoplastic defenses as well as fungal incapability to adapt to stress and interfere with host cell at the biotrophic stage of the interaction.

The application of artificial intelligence in EUS.

Artificial intelligence (AI) is an epoch-making technology, among which the 2 most advanced parts are machine learning and deep learning algorithms that have been further developed by machine learning, and it has been partially applied to assist EUS diagnosis. AI-assisted EUS diagnosis has been reported to have great value in the diagnosis of pancreatic tumors and chronic pancreatitis, gastrointestinal stromal tumors, esophageal early cancer, biliary tract, and liver lesions. The application of AI in EUS diagnosis still has some urgent problems to be solved. First, the development of sensitive AI diagnostic tools requires a large amount of high-quality training data. Second, there is overfitting and bias in the current AI algorithms, leading to poor diagnostic reliability. Third, the value of AI still needs to be determined in prospective studies. Fourth, the ethical risks of AI need to be considered and avoided.

Mitochondrial DNA Leakage Promotes Persistent Pancreatic Acinar Cell Injury in Acute Pancreatitis via the cGAS-STING-NF-κB Pathway.

Previous research has shown that the activation of the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway in macrophages can promote severe acute pancreatitis through the release of inflammatory factors. The role of this pathway in pancreatic acinar cells, however, has not been studied, and understanding its mechanism could be crucial. We analysed plasma from 50 acute pancreatitis (AP) patients and 10 healthy donors using digital PCR, which links mitochondrial DNA (mtDNA) levels to the severity of AP. Single-cell sequencing of the pancreas during AP revealed differentially expressed genes and pathways in acinar cells. Experimental studies using mouse and cell models, which included mtDNA staining and quantitative PCR, revealed mtDNA leakage and the activation of STING-related pathways, indicating potential inflammatory mechanisms in AP. In conclusion, our study revealed that the mtDNA-STING-nuclear factor κB(NF-κB) pathway in pancreatic acinar cells could be a novel pathogenic factor in AP.

Single-cell RNA sequencing reveals the process of CA19-9 production and dynamics of the immune microenvironment between CA19-9 (+) and CA19-9 (-) PDAC.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the main types of malignant tumor of the digestive system, and patient prognosis is affected by difficulties in early diagnosis, poor treatment response, and a high postoperative recurrence rate. Carbohydrate antigen 19-9 (CA19-9) has been widely used as a biomarker for the diagnosis and postoperative follow-up of PDAC patients. Nevertheless, the production mechanism and potential role of CA19-9 in PDAC progression have not yet been elucidated. METHODS: We performed single-cell RNA sequencing on six samples pathologically diagnosed as PDAC (three CA19-9-positive and three CA19-9-negative PDAC samples) and two paracarcinoma samples. We also downloaded and integrated PDAC samples (each from three CA19-9-positive and CA19-9-negative patients) from an online database. The dynamics of the proportion and potential function of each cell type were verified through immunofluorescence. Moreover, we built an in vitro coculture cellular model to confirm the potential function of CA19-9. RESULTS: Three subtypes of cancer cells with a high ability to produce CA19-9 were identified by the markers TOP2A , AQP5 , and MUC5AC . CA19-9 production bypass was discovered on antigen-presenting cancer-associated fibroblasts (apCAFs). Importantly, the proportion of immature ficolin-1 positive (FCN1+) macrophages was high in the CA19-9-negative group, and the proportion of mature M2-like macrophages was high in the CA19-9-positive group. High proportions of these two macrophage subtypes were associated with an unfavourable clinical prognosis. Further experiments indicated that CA19-9 could facilitate the transformation of M0 macrophages into M2 macrophages in the tumor microenvironment. CONCLUSIONS: Our study described CA19-9 production at single-cell resolution and the dynamics of the immune atlas in CA19-9-positive and CA19-9-negative PDAC. CA19-9 could promote M2 polarization of macrophage in the pancreatic tumor microenvironment.

Microplastics are detected in human gallstones and have the ability to form large cholesterol-microplastic heteroaggregates.

Ubiquitous pollution due to microplastics through the food chain is a major cause of various deleterious effects on the human health. The aim of this study was to determine the existence of microplastics and the internal mechanism of microplastics as accelerators of cholelithiasis. Gallstones were collected from 16 patients after cholecystectomy, and microplastics in the gallstones were detected through laser direct infrared and pyrolysis gas chromatographymass spectrometry examinations. Mice model of gallstone were constructed with or without different diameters of microplastic (0.5, 5 and 50 µm). The affinity between microplastic and cholesterol or bilirubin was tested by co-culturing and qualified using molecular dynamics simulations. Finally, altered gut microbiota among the groups were identified using 16 s rRNA sequencing. The presence of microplastics in the gallstones of all the patients were confirmed. Microplastic content was significantly higher in younger chololithiasis patients (age<50 years). Mice fed a high-cholesterol diet with microplastic drinks showed more severe chololithiasis. In terms of the mechanism, microplastics showed a higher affinity for cholesterol than for bilirubin. Significant alterations in the gut microbiota have also been identified after microplastic intake in mice. Our study revealed the presence of microplastics in human gallstones, showcasing their potential to aggravate chololithiasis by forming large cholesterol-microplastic heteroaggregates and altering the gut microbiota.

Proteotranscriptomic analyses of the midgut and Malpighian tubules after a sublethal concentration of Cry1Ab exposure on Spodoptera litura.

BACKGROUND: Cry1Ab has emerged as a bio-insecticide to control Spodoptera litura (Lepidoptera: Noctuidae). However, the sublethal effects of Cry1Ab on the physiological changes and molecular level of S. litura have not been well documented. Our aims in this study were to assess the sublethal effect of Cry1Ab on S. litura, including midgut and Malpighian tubules as targets. RESULTS: After sublethal Cry1Ab exposure, distinct histological alterations were mainly observed in the midgut. Furthermore, the results of comparative RNA sequencing and tandem mass tag-based proteomics showed that, in the midgut, most differential expression genes (DEGs) were up-regulated and significantly enriched in the serine protease activity pathway, and up-regulated differential expression proteins (DEPs) were mainly associated with the oxidative phosphorylation pathway, whereas the down-regulated involved in the ribosome pathways. In the Malpighian tubules, DEGs and DEPs were significantly enriched in the ribosome pathway. We proposed that ribosome may act as a universal target in energy metabolism with other pathways via the results of protein-protein interaction analysis. Further, by verification of the mRNA expression of some Cry protein receptor and detoxification genes after Cry1Ab treatment, it was suggested that the ribosomal proteins (RPs) possibly participate in influencing the Bt-resistance of S. litura larvae under sublethal Cry1Ab exposure. CONCLUSION: Under sublethal Cry1Ab exposure, the midgut of S. litura was damaged, and the proteotranscriptomic analysis elucidated that Cry1Ab disrupted the energy homeostasis of larvae. Furthermore, we emphasized the potential role of ribosomes in sublethal Cry1Ab exposure. © 2024 Society of Chemical Industry.

Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in Brassica rapa.

BACKGROUND: The species Brassica rapa (2n=20, AA) is an important vegetable and oilseed crop, and serves as an excellent model for genomic and evolutionary research in Brassica species. With the availability of whole genome sequence of B. rapa, it is essential to further determine the activity of all functional elements of the B. rapa genome and explore the transcriptome on a genome-wide scale. Here, RNA-seq data was employed to provide a genome-wide transcriptional landscape and characterization of the annotated and novel transcripts and alternative splicing events across tissues. RESULTS: RNA-seq reads were generated using the Illumina platform from six different tissues (root, stem, leaf, flower, silique and callus) of the B. rapa accession Chiifu-401-42, the same line used for whole genome sequencing. First, these data detected the widespread transcription of the B. rapa genome, leading to the identification of numerous novel transcripts and definition of 5'/3' UTRs of known genes. Second, 78.8% of the total annotated genes were detected as expressed and 45.8% were constitutively expressed across all tissues. We further defined several groups of genes: housekeeping genes, tissue-specific expressed genes and co-expressed genes across tissues, which will serve as a valuable repository for future crop functional genomics research. Third, alternative splicing (AS) is estimated to occur in more than 29.4% of intron-containing B. rapa genes, and 65% of them were commonly detected in more than two tissues. Interestingly, genes with high rate of AS were over-represented in GO categories relating to transcriptional regulation and signal transduction, suggesting potential importance of AS for playing regulatory role in these genes. Further, we observed that intron retention (IR) is predominant in the AS events and seems to preferentially occurred in genes with short introns. CONCLUSIONS: The high-resolution RNA-seq analysis provides a global transcriptional landscape as a complement to the B. rapa genome sequence, which will advance our understanding of the dynamics and complexity of the B. rapa transcriptome. The atlas of gene expression in different tissues will be useful for accelerating research on functional genomics and genome evolution in Brassica species.

Mechanisms of sex differentiation and sex reversal in hermaphrodite fish as revealed by the Epinephelus coioides genome.

Most grouper species are functional protogynous hermaphrodites, but the genetic basis and the molecular mechanisms underlying the regulation of this unique reproductive strategy remain enigmatic. In this study, we report a high-quality chromosome-level genome assembly of the representative orange-spotted grouper (Epinephelus coioides). No duplication or deletion of sex differentiation-related genes was found in the genome, suggesting that sex development in this grouper may be related to changes in regulatory sequences or environmental factors. Transcriptomic analyses showed that aromatase and retinoic acid are probably critical to promoting ovarian fate determination, and follicle-stimulating hormone triggers the female-to-male sex change. Socially controlled sex-change studies revealed that, in sex-changing fish, the brain's response to the social environment may be mediated by activation of the phototransduction cascade and the melatonin synthesis pathway. In summary, our genomic and experimental results provide novel insights into the molecular mechanisms of sex differentiation and sex change in the protogynous groupers.

Transcriptomic analysis of 'Suli' pear (Pyrus pyrifolia white pear group) buds during the dormancy by RNA-Seq.

BACKGROUND: Bud dormancy is a critical developmental process that allows perennial plants to survive unfavorable environmental conditions. Pear is one of the most important deciduous fruit trees in the world, but the mechanisms regulating bud dormancy in this species are unknown. Because genomic information for pear is currently unavailable, transcriptome and digital gene expression data for this species would be valuable resources to better understand the molecular and biological mechanisms regulating its bud dormancy. RESULTS: We performed de novo transcriptome assembly and digital gene expression (DGE) profiling analyses of 'Suli' pear (Pyrus pyrifolia white pear group) using the Illumina RNA-seq system. RNA-Seq generated approximately 100 M high-quality reads that were assembled into 69,393 unigenes (mean length = 853 bp), including 14,531 clusters and 34,194 singletons. A total of 51,448 (74.1%) unigenes were annotated using public protein databases with a cut-off E-value above 10-5. We mainly compared gene expression levels at four time-points during bud dormancy. Between Nov. 15 and Dec. 15, Dec. 15 and Jan. 15, and Jan. 15 and Feb. 15, 1,978, 1,024, and 3,468 genes were differentially expressed, respectively. Hierarchical clustering analysis arranged 190 significantly differentially-expressed genes into seven groups. Seven genes were randomly selected to confirm their expression levels using quantitative real-time PCR. CONCLUSIONS: The new transcriptomes offer comprehensive sequence and DGE profiling data for a dynamic view of transcriptomic variation during bud dormancy in pear. These data provided a basis for future studies of metabolism during bud dormancy in non-model but economically-important perennial species.

Methods of a New Chronic Pancreatitis and Spontaneous Pancreatic Cancer Mouse Model Using Retrograde Pancreatic Duct Injection of Dibutyltin Dichloride.

The current study aimed to develop a new chronic pancreatitis and spontaneous pancreatic cancer model on C57/BL6 mouse through retrograde pancreatic duct injection of dibutyltin dichloride (DBTC) and explore its basic pathological changes as compared to the previous published chronic pancreatitis model through tail vein injection of DBTC with alcohol drinking. C57/BL6 mice were randomly divided into 3 groups: CG (control group; n = 15), VG (tail vein injection of DBTC (8 mg/kg) with 10% alcohol drinking group; n = 20), and PG (retrograde pancreatic duct injection of DBTC group (1 mg/kg); n = 30). Five mice in each group were sacrificed at a specific time point after the first treatment. The pathological section was observed. The activities of amylase, bilirubin, and hyaluronic acid in serum were determined. The expression of fibronectin, COL1A1, α-SMA, MMP-1, and TIMP-1 in the pancreas was assayed. Severe fibrosis of the pancreas with inflammatory cell infiltration could be observed on day 21 in the PG. In the VG, slight fibrosis of the pancreas with inflammatory cell infiltration was observed on day 28. There were significant differences in serum amylase, bilirubin, and hyaluronic acid levels between the PG and VG. The protein level of COL1A1 and α-SMA significantly increased in the PG. The mRNA expression of TIMP-1 is upregulated and the MMP-1 mRNA level is downregulated in the PG. Finally, typical neoplastic pathological change is significantly obvious in the PG. In conclusion, we established and validated a new chronic pancreatitis (CP) and spontaneous pancreatic cancer mouse model through retrograde injection of DBTC into the pancreatic duct. Previously reported mouse model through tail vein injection of DBTC with alcohol drinking could not cause obvious CP and neoplastic pathological change in mice.

Single-cell sequencing reveals heterogeneity between pancreatic adenosquamous carcinoma and pancreatic ductal adenocarcinoma with prognostic value.

Pancreatic adenosquamous carcinoma (ASPC) is a rare subtype of pancreatic cancer with lethal malignancy, and few studies have focused on the heterogeneity of ASPC. Here, we performed a single-cell sequencing procedure on pancreatic tumor tissue from an ASPC patient and a patient with high-grade intraductal papillary mucinous neoplasm (IPMN). Through the combined analysis of single-cell sequencing data from five pancreatic ductal adenocarcinoma (PDAC) patients, one IPMN patient, and one ASPC patient in a public database, we identified 11 main types of cells, including macrophages, B cells, cancer stem cells, ductal cells, fibroblasts, endo/stellate cells, neutrophils, acinar cells, T cells, natural killer (NK) cells, dendritic cells, and mast cells. Then, the different characteristics and differentiation paths of the immune microenvironment among IPMN, ASPC, and PDAC in macrophages, T cells, and cancer-associated fibroblasts (CAFs) were identified through multiple bioinformatics analyses. Two novel special cancer-associated fibroblasts were identified as nCAFs and imCAFs. Then, cancer cells in duct cells were identified using the infercnv software. Two ASPC-specific subgroups of cancer cells with squamous cell features were identified. Finally, the identified specific CAFs and cancer cells were mapped to TCGA-PAAD cohort through the cibersoftx software. All of these identified subgroups were calculated to have a significant prognostic value in pancreatic cancer patients. These findings will promote the clinical application of single-cell sequencing data of pancreatic cancer and deepen our understanding of ASPC.

Circadian regulation of night feeding and daytime detoxification in a formidable Asian pest Spodoptera litura.

Voracious feeding, trans-continental migration and insecticide resistance make Spodoptera litura among the most difficult Asian agricultural pests to control. Larvae exhibit strong circadian behavior, feeding actively at night and hiding in soil during daytime. The daily pattern of larval metabolism was reversed, with higher transcription levels of genes for digestion (amylase, protease, lipase) and detoxification (CYP450s, GSTs, COEs) in daytime than at night. To investigate the control of these processes, we annotated nine essential clock genes and analyzed their transcription patterns, followed by functional analysis of their coupling using siRNA knockdown of interlocked negative feedback system core and repressor genes (SlituClk, SlituBmal1 and SlituCwo). Based on phase relationships and overexpression in cultured cells the controlling mechanism seems to involve direct coupling of the circadian processes to E-boxes in responding promoters. Additional manipulations involving exposure to the neonicotinoid imidacloprid suggested that insecticide application must be based on chronotoxicological considerations for optimal effectiveness.

Lepidopteran wing scales contain abundant cross-linked film-forming histidine-rich cuticular proteins.

Scales are symbolic characteristic of Lepidoptera; however, nothing is known about the contribution of cuticular proteins (CPs) to the complex patterning of lepidopteran scales. This is because scales are resistant to solubilization, thus hindering molecular studies. Here we succeeded in dissolving developing wing scales from Bombyx mori, allowing analysis of their protein composition. We identified a distinctive class of histidine rich (His-rich) CPs (6%-45%) from developing lepidopteran scales by LC-MS/MS. Functional studies using RNAi revealed CPs with different histidine content play distinct and critical roles in constructing the microstructure of the scale surface. Moreover, we successfully synthesized films in vitro by crosslinking a 45% His-rich CP (BmorCPR152) with laccase2 using N-acetyl- dopamine or N-ß-alanyl-dopamine as the substrate. This molecular study of scales provides fundamental information about how such a fine microstructure is constructed and insights into the potential application of CPs as new biomaterials.

[Efficacy of paliperidone extended-release tablets in the improvement of social functions in schizophrenics: a randomized and controlled study].

OBJECTIVE: To explore the efficacy of paliperidone extended-release tablets in the improvement of social functions in schizophrenics. METHODS: A total of 81 schizophrenics were randomly divided into study group with paliperidone extended-release tablets and control group with risperidone for a 12-week treatment. They were assessed and analyzed by positive and negative symptoms scales (PANSS), social disability screening schedule (SDSS) and treatment emergent symptom scale (TESS) at baseline, 6(th) weekend and 12(th) weekend. RESULTS: In study group, the factors and total scores of PANSS in the 12(th) weekend of treatment [(12.0 ± 2.8), (12.1 ± 3.6), (26.2 ± 5.0), (50.2 ± 8.7)] were all significantly lower than those at baseline [(24.7 ± 5.3), (23.8 ± 3.6), (45.0 ± 2.9), (93.5 ± 6.8)] (t = 9.60-16.78, P < 0.05). In study group, the positive factor, negative factor and total scores of PANSS in the 12(th) weekend of treatment [(12.0 ± 2.8), (12.1 ± 3.6), (50.2 ± 8.7)] were all significantly lower than those in the 6(th) weekend of treatment [(14.2 ± 1.8), (14.6 ± 2.4), (56.5 ± 6.4)] (t = 2.58-4.26, P < 0.05). In the 12(th) weekend of treatment, the factors and total scores of PANSS in study group [(12.0 ± 2.8), (12.1 ± 3.6), (26.2 ± 5.0), (50.2 ± 8.7)] were all significantly lower than those in control group [(16.9 ± 4.9), (18.7 ± 5.3), (32.5 ± 5.1), (68.1 ± 13.0)] (t = -4.28--5.67, P < 0.05). In study group, the total scores of SDSS in the 12(th) weekend of treatment (5.93 ± 2.78) were significantly lower than those at baseline (13.9 ± 3.4) (t = 10.83, P < 0.05). In study group, the total scores of SDSS in the 12(th) weekend of treatment (5.9 ± 2.8) were significantly lower than those in the 6(th) weekend of treatment (7.6 ± 2.9) (t = 5.21, P < 0.05). But there was no significant improvement in control group (t = 1.88, P > 0.05). In the 12(th) weekend of treatment, the total scores of SDSS in study group (5.9 ± 2.8) were significantly lower than those in control group (8.8 ± 2.9) (t = -4.49, P < 0.05). No severe adverse effect was reported in either group. CONCLUSION: Paliperidone extended-release tablets are effective to improve social functions and psychiatric symptoms of schizophrenics.

High-throughput identification of heavy metal binding proteins from the byssus of chinese green mussel (Perna viridis) by combination of transcriptome and proteome sequencing.

The Byssus, which is derived from the foot gland of mussels, has been proved to bind heavy metals effectively, but few studies have focused on the molecular mechanisms behind the accumulation of heavy metals by the byssus. In this study, we integrated high-throughput transcriptome and proteome sequencing to construct a comprehensive protein database for the byssus of Chinese green mussel (Perna viridis), aiming at providing novel insights into the molecular mechanisms by which the byssus binds to heavy metals. Illumina transcriptome sequencing generated a total of 55,670,668 reads. After filtration, we obtained 53,047,718 clean reads and subjected them to de novo assembly using Trinity software. Finally, we annotated 73,264 unigenes and predicted a total of 34,298 protein coding sequences. Moreover, byssal samples were analyzed by proteome sequencing, with the translated protein database from the foot transcriptome as the reference for further prediction of byssal proteins. We eventually determined 187 protein sequences in the byssus, of which 181 proteins are reported for the first time. Interestingly, we observed that many of these byssal proteins are rich in histidine or cysteine residues, which may contribute to the byssal accumulation of heavy metals. Finally, we picked one representative protein, Pvfp-5-1, for recombinant protein synthesis and experimental verification of its efficient binding to cadmium (Cd2+) ions.

Genome-wide annotation and comparative analysis of cuticular protein genes in the noctuid pest Spodoptera litura.

Insect cuticle is considered an adaptable and versatile building material with roles in the construction and function of exoskeleton. Its physical properties are varied, as the biological requirements differ among diverse structures and change during the life cycle of the insect. Although the bulk of cuticle consists basically of cuticular proteins (CPs) associated with chitin, the degree of cuticular sclerotization is an important factor in determining its physical properties. Spodoptera litura, the tobacco cutworm, is an important agricultural pest in Asia. Compared to the domestic silkworm, Bombyx mori, another lepidopteran whose CP genes have been well annotated, S. litura has a shorter life cycle, hides in soil during daytime beginning in the 5th instar and is exposed to soil in the pupal stage without the protection of a cocoon. In order to understand how the CP genes may have been adapted to support the characteristic life style of S. litura, we searched its genome and found 287 putative cuticular proteins that can be classified into 9 CP families (CPR with three groups (RR-1, RR-2, RR-3), CPAP1, CPAP3, CPF, CPFL, CPT, CPG, CPCFC and CPLCA), and a collection of unclassified CPs named CPH. There were also 112 cuticular proteins enriched in Histidine residues with content varying from 6% to 30%, comprising many more His-rich cuticular proteins than B. mori. A phylogenetic analysis between S. litura, M. sexta and B. mori uncovered large expansions of RR-1 and RR-2 CPs, forming large gene clusters in different regions of S. litura chromosome 9. We used RNA-seq analysis to document the expression profiles of CPs in different developmental stages and tissues of S. litura. The comparative genomic analysis of CPs between S. litura and B. mori integrated with the unique behavior and life cycle of the two species offers new insights into their contrasting ecological adaptations.

DNA methylation on N6-adenine in lepidopteran Bombyx mori.

Increasing knowledge of DNA methylation that occurs on the sixth position of adenine (N6-methyladenine, 6 mA) has emerged as a novel epigenetic mark in eukaryotes and plays an important role in regulating gene transcription, DNA replication and repair, transposable activities, and others. Here, we show DNA 6 mA methylation is present in Bombyx mori, a lepidopteran model insect, and identify the 6 mA methyltransferase, METTL4, and 6 mA demethylase, NMAD, which regulate the levels of 6 mA in embryogenesis and cultured cells of B. mori. Importantly, RNAi knockdown of METTL4 and NMAD not only induce cell cycle arrest at G1 phase but also result in defects of chromosome alignments at metaphase. We further demonstrate that 6 mA methylation is widely distributed across the genome of B. mori by 6 mA-Seq and primarily enriched in the regulatory regions as well as gene bodies. Integrated analysis of 6 mA-Seq and RNA-Seq reveals that 6 mA methylation in B. mori is preferentially related with lowly expressed genes and negatively correlated with active gene transcription, which provides a novel regulatory mechanism of DNA 6 mA methylation on target genes. Altogether, these data identify 6 mA methylation in B. mori and demonstrate a crucial role of 6 mA signaling in controlling cell cycle progression.

A reference gene set construction using RNA-seq of multiple tissues of Chinese giant salamander, Andrias davidianus.

BACKGROUND: Chinese giant salamander (CGS) is the largest extant amphibian species in the world. Owing to its evolutionary position and four peculiar phenomenon of life (longevity, starvation tolerance, regenerative ability, and hatch without sunshine), it is an invaluable model species for research. However, lack of genomic resources leads to fewer study progresses in these fields, due to its huge genome of ∼50 GB making it extremely difficult to be assembled. RESULTS: We reported the sequenced transcriptome of more than 20 tissues from adult CGS using Illumina Hiseq 2000 technology, and a total of 93 366 no-redundancy transcripts with a mean length of 1326 bp were obtained. We developed for the first time an efficient pipeline to construct a high-quality reference gene set of CGS and obtained 26 135 coding genes. BUSCO and homologous assessment showed that our assembly captured 70.6% of vertebrate universal single-copy orthologs, and this coding gene set had a higher proportion of completeness CDS with comparable quality of the protein sets of Tibetan frog. CONCLUSIONS: These highest quality data will provide a valuable reference gene set to the subsequent research of CGS. In addition, our strategy of de novo transcriptome assembly and protein identification is applicable to similar studies.