Microinjection of antisense oligonucleotides into living mouse testis enables lncRNA function study

Background Long non-coding RNAs (lncRNAs) have been the focus of ongoing research in a diversity of cellular processes. LncRNAs are abundant in mammalian testis, but their biological function remains poorly known. Results Here, we established an antisense oligonucleotides (ASOs)-based targeting approach that can efficiently knock down lncRNA in living mouse testis. We cloned the full-length transcript of lncRNA Tsx (testis-specific X-linked) and defined its testicular localization pattern. Microinjection of ASOs through seminiferous tubules in vivo significantly lowered the Tsx levels in both nucleus and cytoplasm. This effect lasted no less than 10 days, conducive to the generation and maintenance of phenotype. Importantly, ASOs performed better in depleting the nuclear Tsx and sustained longer effect than small interfering RNAs (siRNAs). In addition to the observation of an elevated number of apoptotic germ cells upon ASOs injection, which recapitulates the documented description of Tsx knockout, we also found a specific loss of meiotic spermatocytes despite overall no impact on meiosis and male fertility. Conclusions Our study detailed the characterization of Tsx and illustrates ASOs as an advantageous tool to functionally interrogate lncRNAs in spermatogenesis.

Genome-wide survey and characterization of transcription factors in the silk gland of the silkworm, Bombyx mori

Transcription factors (TFs) are key proteins that modulate gene transcription and thereby lead to changes in the gene expression profile and the subsequent alteration of cellular functions. In the silk gland (SG) of silkworm Bombyx mori, an important silk-producing insect, TFs are of vital importance in the regulation of silk protein synthesis in this organ. However, which TFs exist and express in the SG remains largely unknown. Here, we report the large-scale identification of TFs in the SG based on available full-length transcript sequences and the most recent version of silkworm genome data. In total, 348 candidate TFs were identified by strict filtration and were classified into 56 TF families. Chromosomal distribution, motif composition, and phylogenetic relationship analyses revealed the typical characteristics of these TFs. In addition, the expression patterns of 348 TFs in various tissues of B. mori, especially the SG of fourth-molt (4LM) and day-3 and day-4 fifth-instar (5L3D and 5L4D) larvae, were investigated based on public RNA-seq data and gene microarray data, followed by spatiotemporal verification of TF expression levels by quantitative real-time PCR (qRT-PCR). This report describes the first comprehensive analysis of TFs in the B. mori SG. The results can serve as a baseline for further studies of the roles of TFs in the B. mori SG.

EPCO-14. MULTIFACETED TRANSCRIPTOMIC AND PROTEOMIC ANALYSES IDENTIFIED PUTATIVE ALTERNATIVE SPLICING-DERIVED CELL SURFACE ANTIGENS IN GLIOMA

BACKGROUND To develop effective immunotherapy for gliomas, it is crucial to expand the repertoire of targetable antigens. Recent studies have suggested that alternative splicing (AS), or its deriving tumor-specific junctions (“neojunctions”), could generate cryptic amino acid sequences that can be a source of neoantigens. In this study, we investigated neojunctions based on multifaceted transcriptomic and proteomic analyses, seeking the potential cell surface antigens that may be targeted by CAR. METHODS For screening, we analyzed bulk RNA-sequencing data of TCGA-GBM/LGG with high tumor purity (n = 429) and GTEx normal tissues (n = 9,166). Cohorts of spatially mapped intratumoral samples and longitudinally collected tumors were used to determine clonality and stability of the candidate neojunctions. Nanopore long-read amplicon sequencing was deployed to confirm the full-length transcript sequence. Their protein-level expression was explored by analyzing the Clinical Proteomic Tumor Analysis Consortium (CPTAC)-GBM proteomics dataset. RESULTS In the screening analysis comparing TCGA and GTEx datasets, we identified 218 neojunctions with adequate expression, prevalence, and tumor-specificity. Of these, 12 were predicted to be cell-surface antigens. Eight of the 12, such as BCAN, DLL3, and PTPRZ1, were also observed in multiple cases of another validation dataset. In the analysis of tumors with spatially mapped intratumoral samples, 7 of the 12 were recurrently detected in no less than 50% of the samples in multiple cases. In addition, 5 of the 12 were found to be conserved in primary and recurrent pairs of tumors in multiple cases. Full-length transcript sequencing corroborated our predictions based on short reads, and also demonstrated more complex AS patterns. Finally, CPTAC-GBM proteomics analysis identified one cryptic peptide that substantiated the corresponding transcriptome-based prediction. CONCLUSION: We identified neojunctions with the potential to generate cell-surface antigens. These multifaceted transcriptomic and proteomic analyses provide the rationale to pursue the development of immunotherapy targeting neojunction-derived antigens.

Transcriptome in Combination Proteome Unveils the Phenylpropane Pathway Involved in Garlic (Allium sativum) Greening

Garlic (Allium sativum) is an important vegetable crop that is widely used in cooking and medicine. The greening phenomenon of garlic severely decreases the quality of garlic and hinders garlic processing. To study the mechanism of garlic greening, comprehensive full-length transcript sets were constructed. We detected the differences in greening between Pizhou (PZ) garlic and Laiwu (LW) garlic that were both stored at −2.5°C and protected from light at the same time. The results showed that 60,087 unigenes were respectively annotated to the NR, KEGG, GO, Pfam, eggNOG and Swiss Prot databases, and a total of 30,082 unigenes were annotated. The analysis of differential genes and differential proteins showed that PZ garlic and LW garlic had 923 differentially expressed genes (DEGs), of which 529 genes were up regulated and 394 genes were downregulated. Through KEGG and GO enrichment analysis, it was found that the most significant way of enriching DEGs was the phenylpropane metabolic pathway. Proteomics analysis found that there were 188 differentially expressed proteins (DAPs), 162 up-regulated proteins, and 26 down-regulated proteins between PZ garlic and LW garlic. The content of 10 proteins related to phenylpropanoid biosynthesis in PZ garlic was significantly higher than that of LW garlic. This study explored the mechanisms of garlic greening at a molecular level and further discovered that the formation of garlic green pigment was affected significantly by the phenylpropanoid metabolic pathway. This work provided a theoretical basis for the maintenance of garlic quality during garlic processing and the future development of the garlic processing industries.

Mapping and modeling the genomic basis of differential RNA isoform expression at single-cell resolution with LR-Split-seq

The rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.

SMA-miRs: miR-181a-5p, -324-5p, -451a are overexpressed in spinal muscular atrophy skeletal muscle and serum samples

Background:Spinal muscular atrophy (SMA) is a neuromuscular disorder characterized by the degeneration of the second motor-neuron. The phenotype ranges from very severe to very mild forms. All patients have the homozygous loss of the SMN1 gene and a variable number of SMN2 (generally two-to-four copies), inversely related with the severity. The amazing results of the available treatments have made compelling the need of prognostic biomarkers to predict the progression trajectories of patients. Beside the SMN2 products, few other biomarkers have been evaluated so far, including some miRs.Methods:We performed whole miRNome analysis of muscle samples of patients and controls (14 biopsies and 9 cultures). The levels of muscle differentially expressed miRs were evaluated in serum samples (51 patients and 37 controls) and integrated with SMN2 copies, SMN2-full length transcript levels in blood and age (SMA-score). Results:Over 100 miRs were differentially expressed in SMA muscle; three of them (HSA-miR-181a-5p, -324-5p, -451a; SMA-miRs) were significantly up-regulated in serum of patients. The severity predicted by the SMA-score was related with that of the clinical classification at a correlation coefficient of 0.87 (p<10-5).Conclusions:miRNome analyses suggest the primary involvement of skeletal muscle in SMA pathogenesis; the SMA-miRs are likely actively released in the blood flow, even if their function and target cells require to be elucidated. The accuracy of the SMA-score needs to be verified in replicative studies: if confirmed, its use could be crucial for the routine prognostic assessment, also in pre-symptomatic patients. Funding:Telethon Italia (grant # GGP12116).

Chromosome-scale genome assembly of Eustoma grandiflorum, the first complete genome sequence in family Gentianaceae

Eustoma grandiflorum (Raf.) Shinn., is an annual herbaceous plant native to the southern United States, Mexico, and the Greater Antilles. It has a large flower with a variety of colors and an important flower crop. In this study, we established a chromosome-scale de novo assembly of E. grandiflorum by integrating four genomic and genetic approaches: (1) Pacific Biosciences (PacBio) Sequel deep sequencing, (2) error correction of the assembly by Illumina short reads, (3) scaffolding by chromatin conformation capture sequencing (Hi-C), and (4) genetic linkage maps derived from an F2 mapping population. The 36 pseudomolecules and unplaced 64 scaffolds were created with total length of 1,324.8 Mb. Full-length transcript sequencing was obtained by PacBio Iso-Seq sequencing for gene prediction on the assembled genome, Egra_v1. A total of 36,619 genes were predicted on the genome as high confidence HC) genes. Of the 36,619, 25,936 were annotated functions by ZenAnnotation. Genetic diversity analysis was also performed for nine commercial E. grandiflorum varieties bred in Japan, and 254,205 variants were identified. This is the first report of the construction of reference genome sequences in E. grandiflorum as well as in the family Gentianaceae.

Microinjection of Antisense Oligonucleotides into Living Mouse Testis Enables LncRNA Function Study

BackgroundLong non-coding RNAs (lncRNAs) has been the focus of ongoing research in a diversity of cellular processes. LncRNAs are abundant in mammalian testis, but their biological functions remain poorly known.ResultsHere, we established an antisense oligonucleotides (ASOs)-based targeting approach that can efficiently knock down lncRNA in living mouse testis. We cloned the full-length transcript of lncRNA Tsx (testis-specific X-linked) and defined its testicular localization pattern. Microinjection of ASOs through seminiferous tubules in vivo significantly lowered the Tsx levels in both nucleus and cytoplasm. This effect lasted no less than 10 days, conducive to the generation and maintenance of phenotype. Importantly, ASOs performed better in depleting the nuclear Tsx and sustained longer effect than small interfering RNAs (siRNAs). In addition to the observation of an elevated number of apoptotic germ cells upon ASOs injection, which recapitulates the documented description of Tsx knockout, we also found a specific loss of meiotic spermatocytes despite overall normal process of meiosis.ConclusionsOur study detailed the characterization of Tsx and illustrates ASOs as an advantageous tool to functionally interrogate lncRNAs in spermatogenesis.

Single molecule, full-length transcript sequencing provides insight into the TPS gene family in Paeonia ostii

Background The tree peony (Paeonia section Moutan DC), one of the traditional famous flowers with both ornamental and medicinal value, was widely used in China. Surprisingly little is known about the full-length transcriptome sequencing in tree peony, limiting the research on its gene function and molecular mechanism. The trehalose phosphate phosphatase (TPS) family genes has been found to affect plant growth and development and the function of TPS genes in Paeonia ostii is unknown. Methods In our study, we performed single molecule, full-length transcript sequencing in P. ostii. 10 TPS family members were identified from PacBio sequencing for bioinformatics analysis and transcriptional expression analysis. Results A total of 230,736 reads of insert (ROI) sequences and 114,215 full-Length non-chimeric reads (FLNC) were obtained for further ORFs and transcription factors prediction, SSR analysis and lncRNA identification. NR, Swissprot, GO, COG, KOG, Pfam and KEGG databases were used to obtain annotation information of transcripts. 10 TPS family members were identified with molecular weights between 48.0 to 108.5 kD and isoelectric point between 5.61 to 6.37. Furthermore, we found that TPS family members contain conserved TPP or TPS domain. Based on phylogenetic tree analysis, PoTPS1 protein was highly similar to AtTPS1 protein in Arabidopsis. Finally, we analyzed the expression levels of all TPS genes in P. ostii and found PoTPS5 expressed at the highest level. In conclusion, this study combined the results of the transcriptome to systematically analyze the 10 TPS family members, and sets a framework for further research of this important gene family in development of tree peony.