Retrospective Analysis of Cancer-Specific Gene Expression Panel for Thyroid Fine Needle Aspiration Specimens.

Background: While molecular testing is a promising strategy for preoperative assessment of cytologically indeterminate thyroid nodules, thyroid fine needle aspiration biopsy (FNA) presents unique challenges for molecular assays, including contaminating peripheral blood mononuclear cells (PBMC) and variable numbers of evaluable epithelial thyroid cells. Moreover, the newly recognized entity, noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP), has added an additional challenge to the currently available molecular diagnostic platforms. New diagnostic tools are still needed to correctly distinguish benign and malignant thyroid nodules preoperatively.Methods: Twenty-two transcript splice variants from 12 genes we previously identified as discriminating benign from malignant thyroid nodules were characterized in 80 frozen thyroid tumors from 8 histological subtypes. Isoforms detectable in PBMC were excluded, and the 5 most discriminating isoforms were further validated by real-time quantitative PCR (qPCR) on intraoperative FNA samples from 59 malignant tumors, 55 benign nodules, and 23 NIFTP samples. The qPCR threshold cycle values for each transcript were normalized to the thyrocyte-specific thyroid peroxidase isoform 1 (TPO1) and z-transformed. Receiver operating characteristic (ROC) analyses of the composite transcript scores were used to evaluate classification of thyroid FNAs by the 5-gene isoform expression panel.Results: A molecular signature was developed by combining expression levels of specific isoforms of CDH3, FNDC4, HMGA2, KLK7, and PLAG1. FNAs containing at least 12-36 thyrocytes were sufficient for this assay. The 5-gene composite score achieved an area under the ROC curve (AUC) of 0.86 for distinguishing malignant from benign nodules, with a specificity of 91%, sensitivity of 75%, negative predictive value of 91% and positive predictive value of 74%.Conclusion: Our newly developed 5-gene isoform expression panel distinguishes benign from malignant thyroid tumors and, may help distinguish benign from malignant thyroid nodules in the context of the new NIFTP subtype.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.

Transcriptome analysis reveals significant differences between primary plasma cell leukemia and multiple myeloma even when sharing a similar genetic background

Primary plasma cell leukemia (pPCL) is a highly aggressive plasma cell dyscrasia characterised by short remissions and very poor survival. Although the 17p deletion is associated with poor outcome and extramedullary disease in MM, its presence does not confer the degree of aggressiveness observed in pPCL. The comprehensive exploration of isoform expression and RNA splicing events may provide novel information about biological differences between the two diseases. Transcriptomic studies were carried out in nine newly diagnosed pPCL and ten MM samples, all of which harbored the 17p deletion. Unsupervised cluster analysis clearly distinguished pPCL from MM samples. In total 3584 genes and 20033 isoforms were found to be deregulated between pPCL and MM. There were 2727 significantly deregulated isoforms of non-differentially expressed genes. Strangely enough, significant differences were observed in the expression of spliceosomal machinery components between pPCL and MM, in respect of the gene, isoform and the alternative splicing events expression. In summary, transcriptome analysis revealed significant differences in the relative abundance of isoforms between pPCL and MM, even when they both had the 17p deletion. The mRNA processing pathway including RNA splicing machinery emerged as one of the most remarkable mechanisms underlying the biological differences between the two entities.

Transcript isoform sequencing reveals widespread promoter-proximal transcriptional termination

SUMMARYHigher organisms achieve optimal gene expression by tightly regulating the transcriptional activity of RNA Polymerase II (RNAPII) along DNA sequences of genes1. RNAPII density across genomes is typically highest where two key choices for transcription occur: near transcription start sites (TSSs) and polyadenylation sites (PASs) at the beginning and end of genes, respectively2,3. Alternative TSSs and PASs amplify the number of transcript isoforms from genes4, but how alternative TSSs connect to variable PASs is unresolved from common transcriptomics methods. Here, we define TSS/PAS pairs for individual transcripts in Arabidopsis thaliana using an improved Transcript Isoform sequencing (TIF-seq) protocol and find on average over four different isoforms corresponding to variable TSS/PAS pairs per expressed gene. While intragenic initiation represents a large source of regulated isoform diversity, we discover that ∼ 14% of expressed genes generate relatively unstable short promoter-proximal RNAs (sppRNAs) from nascent transcript cleavage and polyadenylation shortly after initiation. The location of sppRNAs coincides with increased RNAPII density, indicating these large pools of promoter-stalled RNAPII across genomes are often engaged in transcriptional termination. RNAPII elongation factors progress transcription beyond sites of sppRNA formation, demonstrating RNAPII density near promoters represents a checkpoint for early transcriptional termination that governs full-length gene isoform expression.

Recovery of gene haplotypes from a metagenome

Elucidation of population-level diversity of microbiomes is a significant step towards a complete understanding of the evolutionary, ecological and functional importance of microbial communities. Characterizing this diversity requires the recovery of the exact DNA sequence (haplotype) of each gene isoform from every individual present in the community. To address this, we present Hansel and Gretel: a freely-available data structure and algorithm, providing a software package that reconstructs the most likely haplotypes from metagenomes. We demonstrate recovery of haplotypes from short-read Illumina data for a bovine rumen microbiome, and verify our predictions are 100% accurate with long-read PacBio CCS sequencing. We show that Gretel’s haplotypes can be analyzed to determine a significant difference in mutation rates between core and accessory gene families in an ovine rumen microbiome. All tools, documentation and data for evaluation are open source and available via our repository: https://github.com/samstudio8/gretel