The Mitogenome of the Subarctic Octocoral Alcyonium digitatum Reveals a Putative tRNAPro Gene Nested within MutS.

We sequenced and analyzed the complete mitogenome of a Norwegian isolate of the octocoral Alcyonium digitatum using the Ion Torrent sequencing technology. The 18,790 bp circular mitochondrial genome was found to harbor the same set of 17 genes, which encode 14 protein subunits, two structural ribosomal RNAs and one tRNA, as reported in other octocorals. In addition, we detected a new tRNAPro-like gene sequence nested within the MutS protein coding region. This putative tRNA gene feature appears to be conserved among the octocorals but has not been reported previously. The A. digitatum mitogenome was also shown to harbor an optional gene (ORFA) that encodes a putative protein of 191 amino acids with unknown function. A mitogenome-based phylogenetic analysis, presented as a maximum likelihood tree, showed that A. digitatum clustered with high statistical confidence with two other Alcyonium species endemic to the Mediterranean Sea and the Southeast Pacific Ocean.

Blood Transcriptome Analysis of Septic Patients Reveals a Long Non-Coding Alu-RNA in the Complement C5a Receptor 1 Gene.

Many severe inflammation conditions are complement-dependent with the complement component C5a-C5aR1 axis as an important driver. At the RNA level, the blood transcriptome undergoes programmed expression of coding and long non-coding RNAs to combat invading microorganisms. Understanding the expression of long non-coding RNAs containing Alu elements in inflammation is important for reconstructing cell fate trajectories leading to severe disease. We have assembled a pipeline for computation mining of new Alu-containing long non-coding RNAs by intersecting immune genes with known Alu coordinates in the human genome. By applying the pipeline to patient bulk RNA-seq data with sepsis, we found immune genes containing 48 Alu insertion as robust candidates for further study. Interestingly, 1 of the 48 candidates was located within the complement system receptor gene C5aR1 and holds promise as a target for RNA therapeutics.

The Mitochondrial Genome of the Sea Anemone Stichodactyla haddoni Reveals Catalytic Introns, Insertion-Like Element, and Unexpected Phylogeny.

A hallmark of sea anemone mitochondrial genomes (mitogenomes) is the presence of complex catalytic group I introns. Here, we report the complete mitogenome and corresponding transcriptome of the carpet sea anemone Stichodactyla haddoni (family Stichodactylidae). The mitogenome is vertebrate-like in size, organization, and gene content. Two mitochondrial genes encoding NADH dehydrogenase subunit 5 (ND5) and cytochrome c oxidase subunit I (COI) are interrupted with complex group I introns, and one of the introns (ND5-717) harbors two conventional mitochondrial genes (ND1 and ND3) within its sequence. All the mitochondrial genes, including the group I introns, are expressed at the RNA level. Nonconventional and optional mitochondrial genes are present in the mitogenome of S. haddoni. One of these gene codes for a COI-884 intron homing endonuclease and is organized in-frame with the upstream COI exon. The insertion-like orfA is expressed as RNA and translocated in the mitogenome as compared with other sea anemones. Phylogenetic analyses based on complete nucleotide and derived protein sequences indicate that S. haddoni is embedded within the family Actiniidae, a finding that challenges current taxonomy.

The shift from early to late types of ribosomes in zebrafish development involves changes at a subset of rRNA 2'-O-Me sites.

During zebrafish development, an early type of rRNA is gradually replaced by a late type that is substantially different in sequence. We applied RiboMeth-seq to rRNA from developmental stages for profiling of 2'-O-Me, to learn if changes in methylation pattern were a component of the shift. We compiled a catalog of 2'-O-Me sites and cognate box C/D guide RNAs comprising 98 high-confidence sites, including 10 sites that were not known from other vertebrates, one of which was specific to late-type rRNA. We identified a subset of sites that changed in methylation status during development and found that some of these could be explained by availability of their cognate SNORDs. Sites that changed during development were enriched in the novel sites revealed in zebrafish. We propose that the early type of rRNA is a specialized form and that its structure and ribose methylation pattern may be an adaptation to features of development, including translation of specific maternal mRNAs.

The mitochondrial transcriptome of the anglerfish Lophius piscatorius.

OBJECTIVE: Analyze key features of the anglerfish Lophius piscatorius mitochondrial transcriptome based on high-throughput total RNA sequencing. RESULTS: We determined the complete mitochondrial DNA and corresponding transcriptome sequences of L. piscatorius. Key features include highly abundant mitochondrial ribosomal RNAs (10-100 times that of mRNAs), and that cytochrome oxidase mRNAs appeared > 5 times more abundant than both NADH dehydrogenase and ATPase mRNAs. Unusual for a vertebrate mitochondrial mRNA, the polyadenylated COI mRNA was found to harbor a 75 nucleotide 3' untranslated region. The mitochondrial genome expressed several non-canonical genes, including the long noncoding RNAs lncCR-H, lncCR-L and lncCOI. Whereas lncCR-H and lncCR-L mapped to opposite strands in a non-overlapping organization within the control region, lncCOI appeared novel among vertebrates. We found lncCOI to be a highly abundant mitochondrial RNA in antisense to the COI mRNA. Finally, we present the coding potential of a humanin-like peptide within the large subunit ribosomal RNA.

Complete loss of the MHC II pathway in an anglerfish, Lophius piscatorius.

Genome studies in fish provide evidence for the adaptability of the vertebrate immune system, revealing alternative immune strategies. The reported absence of the major compatibility complex (MHC) class II pathway components in certain species of pipefish (genus Syngnathus) and cod-like fishes (order Gadiformes) is of particular interest. The MHC II pathway is responsible for immunization and defence against extracellular threats through the presentation of exogenous peptides to T helper cells. Here, we demonstrate the absence of all genes encoding MHC II components (CD4, CD74 A/B, and both classical and non-classical MHC II α/ß) in the genome of an anglerfish, Lophius piscatorius, indicating loss of the MHC II pathway. By contrast, it has previously been reported that another anglerfish, Antennarius striatus, retains all MHC II genes, placing the loss of MHC II in the Lophius clade to their most recent common ancestor. In the three taxa where MHC II loss has occurred, the gene loss has been restricted to four or five core MHC II components, suggesting that, in teleosts, only these genes have functions that are restricted to the MHC II pathway.

A mitochondrial long noncoding RNA in atlantic cod harbors complex heteroplasmic tandem repeat motifs.

A heteroplasmic tandem repeat (HTR) array occupies 100 to 300 bp of the mitochondrial DNA control region in the Atlantic cod, and recently we noted that the repeat appeared integrated in a polyadenylated mitochondrial long noncoding RNA. Here we provide a more detailed analysis of the mitochondrial HTR in the mitochondrial genome of 134 Atlantic cod specimens. We report all specimens to harbor mitochondrial HTRs in the control region, and identified 26 distinct variants among the 402 repeat motifs assessed. Whereas most specimens contained HTR profiles of 2-5 copies consisting of the same 40-bp motif, 22 specimens showed compound HTR arrays of at least two types of motifs present in the same mitochondrial DNA molecule. We found HTR profiles to be highly conserved between different tissue types of a single individual, and strictly maternally inherited in a mating experiment between parental Atlantic cod expressing different HTR profiles and array motifs. We conclude that mitochondrial heteroplasmy in the control region is very common in Atlantic cod, and results in length heterogenity of the long noncoding RNA lncCR-H.

Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair.

Low-level mitochondrial heteroplasmy is a common phenomenon in both normal and cancer cells. Here, we investigate the link between low-level heteroplasmy and mitogenome mutations in a human breast cancer matched cell line by high-throughput sequencing. We identified 23 heteroplasmic sites, of which 15 were common between normal cells (Hs578Bst) and cancer cells (Hs578T). Most sites were clustered within the highly conserved Complex IV and ribosomal RNA genes. Two heteroplasmic variants in normal cells were found as fixed mutations in cancer cells. This indicates a positive selection of these variants in cancer cells. RNA-Seq analysis identified upregulated L-strand specific transcripts in cancer cells, which include three mitochondrial long non-coding RNA molecules. We hypothesize that this is due to two cancer cell-specific mutations in the control region.

Mitochondrial genome variation of Atlantic cod.

OBJECTIVE: The objective of this study was to analyse intraspecific sequence variation of Atlantic cod mitochondrial DNA, based on a comprehensive collection of completely sequenced mitochondrial genomes. RESULTS: We determined the complete mitochondrial DNA sequence of 124 cod specimens from the eastern and western part of the species' distribution range in the North Atlantic Ocean. All specimens harboured a unique mitochondrial DNA haplotype. Nine hundred and fifty-two polymorphic sites were identified, including 109 non-synonymous sites within protein coding regions. Eighteen variable sites were identified as indels, exclusively distributed in structural RNA genes and non-coding regions. Phylogeographic analyses based on 156 available cod mitochondrial genomes did not reveal a clear structure. There was a lack of mitochondrial genetic differentiation between two ecotypes of cod in the eastern North Atlantic, but eastern and western cod were differentiated and mitochondrial genome diversity was higher in the eastern than the western Atlantic, suggesting deviating population histories. The geographic distribution of mitochondrial genome variation seems to be governed by demographic processes and gene flow among ecotypes that are otherwise characterized by localized genomic divergence associated with chromosomal inversions.

Elucidating the Small Regulatory RNA Repertoire of the Sea Anemone Anemonia viridis Based on Whole Genome and Small RNA Sequencing.

Cnidarians harbor a variety of small regulatory RNAs that include microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), but detailed information is limited. Here, we report the identification and expression of novel miRNAs and putative piRNAs, as well as their genomic loci, in the symbiotic sea anemone Anemonia viridis. We generated a draft assembly of the A. viridis genome with putative size of 313 Mb that appeared to be composed of about 36% repeats, including known transposable elements. We detected approximately equal fractions of DNA transposons and retrotransposons. Deep sequencing of small RNA libraries constructed from A. viridis adults sampled at a natural CO2 gradient off Vulcano Island, Italy, identified 70 distinct miRNAs. Eight were homologous to previously reported miRNAs in cnidarians, whereas 62 appeared novel. Nine miRNAs were recognized as differentially expressed along the natural seawater pH gradient. We found a highly abundant and diverse population of piRNAs, with a substantial fraction showing ping-pong signatures. We identified nearly 22% putative piRNAs potentially targeting transposable elements within the A. viridis genome. The A. viridis genome appeared similar in size to that of other hexacorals with a very high divergence of transposable elements resembling that of the sea anemone genus Exaiptasia. The genome encodes and expresses a high number of small regulatory RNAs, which include novel miRNAs and piRNAs. Differentially expressed small RNAs along the seawater pH gradient indicated regulatory gene responses to environmental stressors.

Distinct Small RNA Signatures in Extracellular Vesicles Derived from Breast Cancer Cell Lines.

Breast cancer is a heterogeneous disease, and different subtypes of breast cancer show distinct cellular morphology, gene expression, metabolism, motility, proliferation, and metastatic potential. Understanding the molecular features responsible for this heterogeneity is important for correct diagnosis and better treatment strategies. Extracellular vesicles (EVs) and their associated molecules have gained much attention as players in intercellular communication, ability to precondition specific organs for metastatic invasion, and for their potential role as circulating cancer biomarkers. EVs are released from the cells and contain proteins, DNA, and long and small RNA species. Here we show by high-throughput small RNA-sequencing that EVs from nine different breast cancer cell lines share common characteristics in terms of small RNA content that are distinct from their originating cells. Most strikingly, a highly abundant small RNA molecule derived from the nuclear 28S rRNA is vastly enriched in EVs. The miRNA profiles in EVs correlate with the cellular miRNA expression pattern, but with a few exceptions that includes miR-21. This cancer-associated miRNA is retained in breast cancer cell lines. Finally, we report that EVs from breast cancer cell lines cluster together based on their small RNA signature when compared to EVs derived from other cancer cell lines. Altogether, our data demonstrate that breast cancer cell lines manifest a specific small RNA signature in their released EVs. This opens up for further evaluation of EVs as breast cancer biomarkers.

An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA.

BACKGROUND: Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNA(Thr) and tRNA(Pro) genes ("T-P spacers") have been observed in several taxa, including gadiform species, but information about their biological roles and putative functions is still lacking. RESULTS: Sequence characterization of the complete European hake Merluccius merluccius mitogenome identified a complex T-P spacer ranging in size from 223-532 bp. Further analyses of 32 gadiform species, representing 8 families and 28 genera, revealed the evolutionary preserved presence of T-P spacers across all taxa. Molecular complexity of the T-P spacers was found to be coherent with the phylogenetic relationships, supporting a common ancestral origin and gain of function during codfish evolution. Intraspecific variation of T-P spacer sequences was assessed in 225 Atlantic cod specimens and revealed 26 haplotypes. Pyrosequencing data representing the mito-transcriptome poly (A) fraction in Atlantic cod identified an abundant H-strand specific long noncoding RNA of about 375 nt. The T-P spacer corresponded to the 5' part of this transcript, which terminated within the control region in a tail-to-tail configuration with the L-strand specific transcript (the 7S RNA). CONCLUSIONS: The T-P spacer is inferred to be evolutionary preserved in gadiform mitogenomes due to gain of function through a long noncoding RNA. We suggest that the T-P spacer adds stability to the H-strand specific long noncoding RNA by forming stable hairpin structures and additional protein binding sites.

Performance comparison of digital microRNA profiling technologies applied on human breast cancer cell lines.

MicroRNA profiling represents an important first-step in deducting individual RNA-based regulatory function in a cell, tissue, or at a specific developmental stage. Currently there are several different platforms to choose from in order to make the initial miRNA profiles. In this study we investigate recently developed digital microRNA high-throughput technologies. Four different platforms were compared including next generation SOLiD ligation sequencing and Illumina HiSeq sequencing, hybridization-based NanoString nCounter, and miRCURY locked nucleic acid RT-qPCR. For all four technologies, full microRNA profiles were generated from human cell lines that represent noninvasive and invasive tumorigenic breast cancer. This study reports the correlation between platforms, as well as a more extensive analysis of the accuracy and sensitivity of data generated when using different platforms and important consideration when verifying results by the use of additional technologies. We found all the platforms to be highly capable for microRNA analysis. Furthermore, the two NGS platforms and RT-qPCR all have equally high sensitivity, and the fold change accuracy is independent of individual miRNA concentration for NGS and RT-qPCR. Based on these findings we propose new guidelines and considerations when performing microRNA profiling.

Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals.

Expression and processing of mitochondrial gene transcripts are fundamental to mitochondrial function, but information from early vertebrates like teleost fishes is essentially lacking. We have analyzed mitogenome sequences of ten codfishes (family Gadidae), and provide complete sequences from three new species (Saithe, Pollack and Blue whiting). Characterization of the mitochondrial mRNAs in Saithe and Atlantic cod identified a set of ten poly(A) transcripts, and six UAA stop codons are generated by posttranscriptional polyadenylation. Structural assessment of poly(A) sites is consistent with an RNaseP cleavage activity 5' of tRNA acceptor-like stems. COI, ND5 and ND6 mRNAs were found to harbor 3' UTRs with antisense potential extending into neighboring gene regions. While the 3' UTR of COI mRNA is complementary to the tRNA(Ser UCN) and highly similar to that detected in human mitochondria, the ND5 and ND6 3' UTRs appear more heterogenic. Deep sequencing confirms expression of all mitochondrial mRNAs and rRNAs, and provides information about the precise 5' ends in mature transcripts. Our study supports an overall evolutionary conservation in mitochondrial RNA processing events among vertebrates, but reveals some unique 5' and 3' end characteristics in codfish mRNAs with implications to antisense regulation of gene expression.

RNA deep sequencing of the Atlantic cod transcriptome.

The Atlantic cod (Gadus morhua) is an emerging aquaculture species. Efforts to develop and characterize its genomic recourses, including draft-grade genome sequencing, have been initiated by the research community. The transcriptome represents the whole complement of RNA transcripts in cells and tissues and reflects the expressed genes at various life stages, tissue types, physiological states, and environmental conditions. We are investigating the Atlantic cod transcriptome by Roche 454, Illumina GA, and ABI SOLiD deep sequencing platforms and corresponding bioinformatics. Both embryonic developmental stages and adult tissues are studied. Here we summarize our recent progress in the analyses of nuclear and mitochondrial polyA mRNAs, non-protein-coding intermediate RNAs, and regulatory microRNAs.

Large-scale sequence analyses of Atlantic cod.

The Atlantic cod (Gadus morhua) is a key species in the North Atlantic ecosystem and commercial fisheries, with increasing aquacultural production in several countries. A Norwegian effort to sequence the complete 0.9Gbp genome by the 454 pyrosequencing technology has been initiated and is in progress. Here we review recent progress in large-scale sequence analyses of the nuclear genome, the mitochondrial genome and genome-wide microRNA identification in the Atlantic cod. The nuclear genome will be de novo sequenced with 25 times oversampling. A total of 120 mitochondrial genomes, sampled from several locations in the North Atlantic, are being completely sequenced by Sanger technology in a high-throughput pipeline. These sequences will be included in a new database for maternal marker reference of Atlantic cod diversity. High-throughput 454 sequencing, as well as Evolutionary Image Array (EvoArray) informatics, is used to investigate the complete set of expressed microRNAs and corresponding mRNA targets in various developmental stages and tissues. Information about microRNA profiles will be essential in the understanding of transcriptome complexity and regulation. Finally, developments and perspectives of Atlantic cod aquaculture are discussed in the light of next-generation high-throughput sequence technologies.

Influence of oligoguluronates on alginate gelation, kinetics, and polymer organization.

Structural polysaccharides of the alginate family form gels in aqueous Ca2+-containing solutions by lateral association of chain segments. The effect of adding oligomers of alpha-l-guluronic acid (G blocks) to gelling solutions of alginate was investigated using rheology and atomic force microscopy (AFM). Ca-alginate gels were prepared by in situ release of Ca2+. The gel strength increased with increasing level of calcium saturation of the alginate and decreased with increasing amount of free G blocks. The presence of free G blocks also led to an increased gelation time. The gel point and fractal dimensionalities of the gels were determined based on the rheological characterization. Without added free G blocks the fractal dimension of the gels increased from df = 2.14 to df = 2.46 when increasing [Ca2+] from 10 to 20 mM. This increase was suggested to arise from an increased junction zone multiplicity induced by the increased concentration of calcium ions. In the presence of free G blocks (G block/alginate = 1/1) the fractal dimension increased from 2.14 to 2.29 at 10 mM Ca2+, whereas there was no significant change associated with addition of G blocks at 20 mM Ca2+. These observations indicate that free G blocks are involved in calcium-mediated bonds formed between guluronic acid sequences within the polymeric alginates. Thus, the added oligoguluronate competes with the alginate chains for the calcium ions. The gels and pregel situations close to the gel point were also studied using AFM. The AFM topographs indicated that in situations of low calcium saturation microgels a few hundred nanometers in diameter develop in solution. In situations of higher calcium saturation lateral association of a number of alginate chains are occurring, giving ordered fiber-like structures. These results show that G blocks can be used as modulators of gelation kinetics as well as local network structure formation and equilibrium properties in alginate gels.