Events of alternative splicing in head and neck cancer via RNA sequencing - an update.

BACKGROUND: Alternative splicing (AS) is a regulatory mechanism used to create many forms of mature messengers RNAs (mRNAs) from the same gene. Sequencing of RNA (RNA-Seq) is an advanced technology, which has been utilized by different studies to find AS mechanisms in head and neck cancer (HNC). Hitherto, there is no available review that could inform us of the major findings from these studies. Hence, we aim to perform a systematic literature search following PRISMA guidelines to study AS events in HNC identified through RNA-Seq studies. RESULTS: A total of five records were identified that utilized RNA-Seq data for identifying AS events in HNC. Five software was used in these records to identify AS events. Two genes influenced by AS i.e. MLL3 and RPS9 were found to be common in 4 out of 5 records. Likewise, 38 genes were identified to be similar in at least 3 records. CONCLUSIONS: Alternative splicing in HNC is a multifaceted regulatory mechanism of gene expression. It can be studied via RNA-Seq using different bioinformatics tools. Genes MLL3, as well as RPS9, were repeatedly found to be associated with HNC, however needs further functional validation.

WGSSAT: A High-Throughput Computational Pipeline for Mining and Annotation of SSR Markers From Whole Genomes.

Mining and characterization of Simple Sequence Repeat (SSR) markers from whole genomes provide valuable information about biological significance of SSR distribution and also facilitate development of markers for genetic analysis. Whole genome sequencing (WGS)-SSR Annotation Tool (WGSSAT) is a graphical user interface pipeline developed using Java Netbeans and Perl scripts which facilitates in simplifying the process of SSR mining and characterization. WGSSAT takes input in FASTA format and automates the prediction of genes, noncoding RNA (ncRNA), core genes, repeats and SSRs from whole genomes followed by mapping of the predicted SSRs onto a genome (classified according to genes, ncRNA, repeats, exonic, intronic, and core gene region) along with primer identification and mining of cross-species markers. The program also generates a detailed statistical report along with visualization of mapped SSRs, genes, core genes, and RNAs. The features of WGSSAT were demonstrated using Takifugu rubripes data. This yielded a total of 139 057 SSR, out of which 113 703 SSR primer pairs were uniquely amplified in silico onto a T. rubripes (fugu) genome. Out of 113 703 mined SSRs, 81 463 were from coding region (including 4286 exonic and 77 177 intronic), 7 from RNA, 267 from core genes of fugu, whereas 105 641 SSR and 601 SSR primer pairs were uniquely mapped onto the medaka genome. WGSSAT is tested under Ubuntu Linux. The source code, documentation, user manual, example dataset and scripts are available online at https://sourceforge.net/projects/wgssat-nbfgr.

Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes.

Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.

Pedicled Flap vs. Free Flap Reconstruction in Head and Neck Cancers: Clinical Outcome Analysis from a Single Surgical Team.

Head and neck cancers are one of the leading causes of morbidity and mortality in Indian population. Reconstruction and restoration of function are also of paramount importance in these patients. The aim of this study was to report outcomes for free flaps and pedicle flaps in patients with oral cavity cancers. A retrospective analysis of 628 patients with oral cavity cancers who underwent reconstruction with either free flaps or pedicle flaps during 2014-2020 was done. The median age of the cohort was 49 years. The free flap reconstruction was performed in 481 (76%) and pedicle flap in 147 (24%) patients. Among free and pedicle flaps, 27 (5.6%) and 3 (2.1%) respectively had major flap complications and 25 (5.1%) and 14 (9.9%) respectively had minor complications. CCI score > 4 was associated with higher events (p = 0.02) in free flap group. The outcomes of free flaps are similar in comparison to pedicle flaps in patients with oral cavity cancers. The higher CCI score is significantly associated with increased flap-related complications for the free flap group.

The genome of walking catfish Clarias magur (Hamilton, 1822) unveils the genetic basis that may have facilitated the development of environmental and terrestrial adaptation systems in air-breathing catfishes.

The walking catfish Clarias magur (Hamilton, 1822) (magur) is an important catfish species inhabiting the Indian subcontinent. It is considered as a highly nutritious food fish and has the capability to walk to some distance, and survive a considerable period without water. Assembly, scaffolding and several rounds of iterations resulted in 3,484 scaffolds covering ∼94% of estimated genome with 9.88 Mb largest scaffold, and N50 1.31 Mb. The genome possessed 23,748 predicted protein encoding genes with annotation of 19,279 orthologous genes. A total of 166 orthologous groups represented by 222 genes were found to be unique for this species. The Computational Analysis of gene Family Evolution (CAFE) analysis revealed expansion of 207 gene families and 100 gene families have rapidly evolved. Genes specific to important environmental and terrestrial adaptation, viz. urea cycle, vision, locomotion, olfactory and vomeronasal receptors, immune system, anti-microbial properties, mucus, thermoregulation, osmoregulation, air-breathing, detoxification, etc. were identified and critically analysed. The analysis clearly indicated that C. magur genome possessed several unique and duplicate genes similar to that of terrestrial or amphibians' counterparts in comparison to other teleostean species. The genome information will be useful in conservation genetics, not only for this species but will also be very helpful in such studies in other catfishes.

In Silico Mining of Conserved miRNAs of Indian Catfish Clarias batrachus (Linnaeus, 1758) from the Contigs, ESTs, and BAC End Sequences.

MicroRNAs are small non-coding RNAs which play significant role in RNA interference. The present work deals with mining of the conserved miRNA and their target genes from the contigs, ESTs, and BAC end sequences of commercially important catfish, Clarias batrachus, from India. A total of 138, 1 and 1 conserved pre-miRNA sequences, were mined from the contigs, ESTs, and BAC end sequences, respectively. The analysis of families of the conserved pre-miRNA revealed conservation of the fish-specific family mir-430 and other important families, such as mir-455, let-7, mir-133, and mir-137. The mir-455 is involved in hypoxia signaling, let-7 family represents potential anti-tumor molecules involved in human cancer therapy, whereas mir-133 and mir-137 have high therapeutic potentials. Using an alternate computational in silico approach, mining of mature miRNAs resulted in identification of 210 mature miRNAs from contigs, 1 from EST, and 2 each from forward as well as reverse BAC end sequences. Target prediction of these putative miRNAs resulted in the identification of 66,758 and 18,747 target genes in C. batrachus and Danio rerio, respectively. Functional annotation of these miRNAs indicated their involvement in diverse biological functions. The findings of the present study can serve as a valuable resource for further functional genomics studies in C. batrachus.

In silico genome wide mining of conserved and novel miRNAs in the brain and pineal gland of Danio rerio using small RNA sequencing data.

MicroRNAs (miRNAs) are small, non-coding RNA molecules that bind to the mRNA of the target genes and regulate the expression of the gene at the post-transcriptional level. Zebrafish is an economically important freshwater fish species globally considered as a good predictive model for studying human diseases and development. The present study focused on uncovering known as well as novel miRNAs, target prediction of the novel miRNAs and the differential expression of the known miRNA using the small RNA sequencing data of the brain and pineal gland (dark and light treatments) obtained from NCBI SRA. A total of 165, 151 and 145 known zebrafish miRNAs were found in the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Chromosomes 4 and 5 of zebrafish reference assembly GRCz10 were found to contain maximum number of miR genes. The miR-181a and miR-182 were found to be highly expressed in terms of number of reads in the brain and pineal gland, respectively. Other ncRNAs, such as tRNA, rRNA and snoRNA, were curated against Rfam. Using GRCz10 as reference, the subsequent bioinformatic analyses identified 25, 19 and 9 novel miRNAs from the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Targets of the novel miRNAs were identified, based on sequence complementarity between miRNAs and mRNA, by searching for antisense hits in the 3'-UTR of reference RNA sequences of the zebrafish. The discovery of novel miRNAs and their targets in the zebrafish genome can be a valuable scientific resource for further functional studies not only in zebrafish but also in other economically important fishes.

Development and characterization of genic SSR markers from low depth genome sequence of Clarias batrachus (magur).

Indian magur (Clarias batrachus) is an important freshwater catfish, which is listed as endangered under A3cde+4acde ver. 3.1 categories by the IUCN (2015) due to decreasing population trend. Microsatellites or short sequence repeats (SSRs) tagged to genes have been utilized as gene marker. In the present study, 31,814 SSRs of C. batrachus (magur) were identified using microsatellite identification tool programme from the next generation sequencing data generated on Roche 454 and Ion Torrent platforms. A bioinformatics pipeline, with stringent criteria resulted in selection of 1672 microsatellite loci falling in the genic region. Initially, a total of 30 loci were selected for primer development; and of these 14 were successfully amplified and five were found to be polymorphic in 30 individuals of C. batrachus (magur). The observed as well as expected heterozygosity ranged from 0.038 to 0.526 and 0.434 to 0.784, respectively, and the number of observed alleles ranged from three to five. The study reported the application of next generation sequencing technologies for rapid development of microsatellite loci in Indian catfish species, C. batrachus (magur).

Low-depth shotgun sequencing resolves complete mitochondrial genome sequence of Labeo rohita.

Labeo rohita, popularly known as rohu, is a widely cultured species in whole Indian subcontinent. In the present study, we used in-silico approach to resolve complete mitochondrial genome of rohu. Low-depth shotgun sequencing using Roche 454 GS FLX (Branford, Connecticut, USA) followed by de novo assembly in CLC Genomics Workbench version 7.0.4 (Aarhus, Denmark) revealed the complete mitogenome of L. rohita to be 16 606 bp long (accession No. KR185963). It comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 putative control region. The gene order and organization are similar to most vertebrates. The mitogenome in the present investigation has 99% similarity with that of previously reported mitogenomes of rohu and this is also evident from the phylogenetic study using maximum-likelihood (ML) tree method. This study was done to determine the feasibility, accuracy and reliability of low-depth sequence data obtained from NGS platform as compared to the Sanger sequencing. Thus, NGS technology has proven to be competent and a rapid in-silico alternative to resolve the complete mitochondrial genome sequence, thereby reducing labors and time.

Assembly and variation analyses of Clarias batrachus mitogenome retrieved from WGS data and its phylogenetic relationship with other catfishes.

Whole genome sequencing (WGS) using next generation sequencing technologies paves the way to sequence the mitochondrial genomes with greater ease and lesser time. Here, we used the WGS data of Clarias batrachus, generated from Roche 454 and Ion Torrent sequencing platforms, to assemble the complete mitogenome using both de novo and reference based approaches. Both the methods yielded almost similar results and the best assembled mitogenome was of 16,510 bp size (GenBank Acc. No. KM259918). The mitogenome annotation resulted in 13 coding genes, 22 tRNA genes, 2 rRNA genes and one control region, and the gene order was found to be identical with other catfishes. Variation analyses between assembled and the reference (GenBank Acc. No. NC_023923) mitogenome revealed 51 variations. The phylogenetic analysis of coding DNA sequences and tRNA supports the monophyly of catfishes. Two SSRs were identified in C. batrachus mitogenome, out of which one was unique to this species. Based on the relative rate of gene evolution, protein coding mitochondrial genes were found to evolve at a much faster pace than the d-loop, which in turn are followed by the rRNAs; the tRNAs showed wide variability in the rate of sequence evolution, and on average evolve the slowest. Among the coding genes, ND2 evolves most rapidly. The variations present in the coding regions of the mitogenome and their comparative analyses with other catfish species may be useful in species conservation and management programs.