Analysis of 13 cell types reveals evidence for the expression of numerous novel primate- and tissue-specific microRNAs.

Two decades after the discovery of the first animal microRNA (miRNA), the number of miRNAs in animal genomes remains a vexing question. Here, we report findings from analyzing 1,323 short RNA sequencing samples (RNA-seq) from 13 different human tissue types. Using stringent thresholding criteria, we identified 3,707 statistically significant novel mature miRNAs at a false discovery rate of ≤ 0.05 arising from 3,494 novel precursors; 91.5% of these novel miRNAs were identified independently in 10 or more of the processed samples. Analysis of these novel miRNAs revealed tissue-specific dependencies and a commensurate low Jaccard similarity index in intertissue comparisons. Of these novel miRNAs, 1,657 (45%) were identified in 43 datasets that were generated by cross-linking followed by Argonaute immunoprecipitation and sequencing (Ago CLIP-seq) and represented 3 of the 13 tissues, indicating that these miRNAs are active in the RNA interference pathway. Moreover, experimental investigation through stem-loop PCR of a random collection of newly discovered miRNAs in 12 cell lines representing 5 tissues confirmed their presence and tissue dependence. Among the newly identified miRNAs are many novel miRNA clusters, new members of known miRNA clusters, previously unreported products from uncharacterized arms of miRNA precursors, and previously unrecognized paralogues of functionally important miRNA families (e.g., miR-15/107). Examination of the sequence conservation across vertebrate and invertebrate organisms showed 56.7% of the newly discovered miRNAs to be human-specific whereas the majority (94.4%) are primate lineage-specific. Our findings suggest that the repertoire of human miRNAs is far more extensive than currently represented by public repositories and that there is a significant number of lineage- and/or tissue-specific miRNAs that are uncharacterized.

The human platelet: strong transcriptome correlations among individuals associate weakly with the platelet proteome.

BACKGROUND: For the anucleate platelet it has been unclear how well platelet transcriptomes correlate among different donors or across different RNA profiling platforms, and what the transcriptomes' relationship is with the platelet proteome. We profiled the platelet transcriptome of 10 healthy young males (5 white and 5 black) with no notable clinical history using RNA sequencing and by Affymetrix microarray. RESULTS: We found that the abundance of platelet mRNA transcripts was highly correlated across the 10 individuals, independently of race and of the employed technology. Our RNA-seq data showed that these high inter-individual correlations extend beyond mRNAs to several categories of non-coding RNAs. Pseudogenes represented a notable exception by exhibiting a difference in expression by race. Comparison of our mRNA signatures to a publicly available quantitative platelet proteome showed that most (87.5%) identified platelet proteins had a detectable corresponding mRNA. However, a high number of mRNAs that were present in the transcriptomes of all 10 individuals had no representation in the proteome. Spearman correlations of the relative abundances for those genes represented by both an mRNA and a protein showed a weak (~0.3) connection. Further analysis of the overlapping and non-overlapping platelet mRNAs and proteins identified gene groups corresponding to distinct cellular processes. CONCLUSIONS: The results of our analyses provide novel insights for platelet biology, show only a weak connection between the platelet transcriptome and proteome, and indicate that it is feasible to assemble a platelet mRNA-ome that can serve as a reference for future platelet transcriptomic studies of human health and disease.

The complex transcriptional landscape of the anucleate human platelet.

BACKGROUND: Human blood platelets are essential to maintaining normal hemostasis, and platelet dysfunction often causes bleeding or thrombosis. Estimates of genome-wide platelet RNA expression using microarrays have provided insights to the platelet transcriptome but were limited by the number of known transcripts. The goal of this effort was to deep-sequence RNA from leukocyte-depleted platelets to capture the complex profile of all expressed transcripts. RESULTS: From each of four healthy individuals we generated long RNA (≥40 nucleotides) profiles from total and ribosomal-RNA depleted RNA preparations, as well as short RNA (<40 nucleotides) profiles. Analysis of ~1 billion reads revealed that coding and non-coding platelet transcripts span a very wide dynamic range (≥16 PCR cycles beyond ß-actin), a result we validated through qRT-PCR on many dozens of platelet messenger RNAs. Surprisingly, ribosomal-RNA depletion significantly and adversely affected estimates of the relative abundance of transcripts. Of the known protein-coding loci, ~9,500 are present in human platelets. We observed a strong correlation between mRNAs identified by RNA-seq and microarray for well-expressed mRNAs, but RNASeq identified many more transcripts of lower abundance and permitted discovery of novel transcripts. CONCLUSIONS: Our analyses revealed diverse classes of non-coding RNAs, including: pervasive antisense transcripts to protein-coding loci; numerous, previously unreported and abundant microRNAs; retrotransposons; and thousands of novel un-annotated long and short intronic transcripts, an intriguing finding considering the anucleate nature of platelets. The data are available through a local mirror of the UCSC genome browser and can be accessed at: http://cm.jefferson.edu/platelets_2012/.

Whole-exome sequencing of DNA from peripheral blood mononuclear cells (PBMC) and EBV-transformed lymphocytes from the same donor.

BACKGROUND: The creation of lymphoblastoid cell lines (LCLs) through Epstein-Barr virus (EBV) transformation of B-lymphocytes can result in a valuable biomaterial for cell biology research and a renewable source of DNA. While LCLs have been used extensively in cellular and genetic studies, the process of cell transformation and expansion during culturing may introduce genomic changes that may impact their use and the interpretation of subsequent genetic findings. RESULTS: We performed whole exome sequencing on a tetrad family using DNA derived from peripheral blood mononuclear cells (PBMCs) and LCLs from each individual. We generated over 4.7 GB of mappable sequence to a 125X read coverage per sample. An average of 19,354 genetic variants were identified. Comparison of the two DNA sources from each individual showed an average concordance rate of 95.69%. By lowering the variant calling parameters, the concordance rate between the paired samples increased to 99.82%. Sanger sequencing of a subset of the remaining discordant variants did confirm the presence of de novo mutations arising in LCLs. CONCLUSIONS: By varying software stringency parameters, we identified 99% concordance between DNA sequences derived from the two different sources from the same donors. These results suggest that LCLs are an appropriate representation of the genetic material of the donor and suggest that EBV transformation can result in low-level generation of de novo mutations. Therefore, use of PBMC or early passage EBV-transformed cells is recommended. These findings have broad-reaching implications, as there are thousands of LCLs in public biorepositories and individual laboratories.

Genotyping beta-globin gene mutations on copolymer-coated glass slides with the ligation detection reaction.

BACKGROUND: Methods are needed to analyze small amounts of samples for variation in disease-causing genes. One means is to couple the sensitivity and multiplexing capability of the ligation detection reaction (LDR) with the use of simple glass slides specifically functionalized with a novel polymer coating to enhance sensitivity. METHODS: We developed an array-based genotyping assay based on glass slides coated with copolymer (N,N-dimethylacrylamide, N,N-acryloyloxysuccinimide, and 3-(trimethoxysilyl)propyl methacrylate). The assay consists of an LDR with genomic DNA followed by a universal PCR (U-PCR) of genomic DNA-templated LDR product. The LDR occurs in the presence of 3 primers for each sequence variant under investigation: 2 distinguishing primers (allele specific and perfectly complementary to wild-type and mutant alleles) and 1 common locus-specific primer. The 2 allele-specific primers have different capture sequences for binding different complementary probes on a tag array. The LDR product templated from genomic DNA is made fluorescent during the U-PCR via incorporation of a Cy5-labeled universal primer into all LDR products; detection occurs on the coated glass slides. RESULTS: The assay was designed to detect 7 prevalent mutations in the beta-globin gene (HBB, hemoglobin, beta) in a multiplex format, and signals for the different alleles are detected by their fluorescence. The assay was applied to 40 genomic DNA samples from both control individuals and patients with known beta-thalassemia mutations. Results show good correspondence between the patients' genotypes as assessed by DNA sequence analysis and those generated from the LDR assays. CONCLUSIONS: The developed technology allows accurate identification of sequence variants in a simple, cost-effective way and offers good flexibility for scaling to other applications with different numbers of single-nucleotide polymorphisms or mutations to be detected.

Identification of APC gene mutations in colorectal cancer using universal microarray-based combinatorial sequencing-by-hybridization.

Familial adenomatous polyposis (FAP) is an autosomal dominant inherited form of colorectal cancer, caused mostly by mutations in the APC gene. Due to the wide variety of mutations found and the large size of the APC gene, several methods of mutation detection are used, which can be time consuming and costly. Here we demonstrate a new method of mutation detection in the APC gene using an array-based approach termed combinatorial sequencing-by-hybridization (cSBH). In cSBH, a universal probe set is attached to a support and a second one is in solution. Two-probe ligation occurs when a DNA strand from the target PCR product consecutively anneals to both unlabeled array-bound and solution-phase dye-labeled probe, creating all target-complementary long-labeled probes attached to the surface. A standard array reader scores fluorescent signals at each array position. Cell lines and patient DNA with known APC gene mutations were analyzed using a cSBH-based HyChip trade mark product. Results show that this universal hexamer (6-mer) chip can successfully detect a range of mutations. Results are very robust for a continuous readout of 3.6 kb from a PCR target, with 99.97% accuracy on a single HyChip trade mark slide. cSBH is a fast, cost-efficient method for first stage mutation screening in the APC or any other gene.

Erythroid-induced commitment of K562 cells results in clusters of differentially expressed genes enriched for specific transcription regulatory elements.

Understanding regulation of fetal and embryonic hemoglobin expression is critical, since their expression decreases clinical severity in sickle cell disease and beta-thalassemia. K562 cells, a human erythroleukemia cell line, can differentiate along erythroid or megakaryocytic lineages and serve as a model for regulation of fetal/embryonic globin expression. We used microarray expression profiling to characterize transcriptomes from K562 cells treated for various times with hemin, an inducer of erythroid commitment. Approximately 5,000 genes were expressed irrespective of treatment. Comparative expression analysis (CEA) identified 899 genes as differentially expressed; analysis by the self-organizing map (SOM) algorithm clustered 425 genes into 8 distinct expression patterns, 322 of which were shared by both analyses. Differential expression of a subset of genes was validated by real-time RT-PCR. Analysis of 5'-flanking regions from differentially expressed genes by PAINT v3.0 software showed enrichment in specific transcription regulatory elements (TREs), some localizing to different expression clusters. This finding suggests coordinate regulation of cluster members by specific TREs. Finally, our findings provide new insights into rate-limiting steps in the appearance of heme-containing hemoglobin tetramers in these cells.

Interpretation of fetal hemoglobin only on newborn screening for hemoglobinopathy.

Newborn screening for hemoglobinopathies rarely produces a fetal hemoglobin only result; it is most consistent with beta-thalassemia major, although other diagnoses are possible. The authors describe two unrelated African-American babies born in North Carolina whose newborn screening revealed fetal hemoglobin only. Both had a relatively benign clinical and hematologic picture. Molecular analyses indicated that both children are compound heterozygotes for beta-thalassemia and pancellular (deletional) hereditary persistence of fetal hemoglobin, a rare and apparently benign condition. Accurate interpretation of the fetal hemoglobin only result on newborn screening requires thorough evaluation, including family studies and molecular analysis.