Search

Molecular subtyping is essential to infer tumor aggressiveness and predict prognosis. In practice, tumor profiling requires in-depth knowledge of bioinformatics tools involved in the processing and analysis of the generated data. Additionally, data incompatibility (e.g., microarray vs. RNA sequencing data) and technical and uncharacterized biological variance between training and test data can pose challenges in classifying individual samples. In this article, we provide a roadmap for implementing bioinformatics frameworks for molecular profiling of human cancers in a clinical diagnostic setting. We describe a framework for integrating several methods for quality control, normalization, batch correction, classification, and reporting and a use case of the framework in breast cancer.

Correlation between rankings of pre- and post-FDR corrected variable importances.

Oral tongue squamous cell carcinomas (OTSCC) are a homogenous group of aggressive tumors in the head and neck region that spread early to lymph nodes and have a higher incidence of regional failure. In addition, there is a rising incidence of oral tongue cancer in younger populations. Studies on functional DNA methylation changes linked with altered gene expression are critical for understanding the mechanisms underlying tumor development and metastasis. Such studies also provide important insight into biomarkers linked with viral infection, tumor metastasis, and patient survival in OTSCC. Therefore, we performed genome-wide methylation analysis of tumors (N = 52) and correlated altered methylation with differential gene expression. The minimal tumor-specific DNA 5-methylcytosine signature identified genes near 16 different differentially methylated regions, which were validated using genomic data from The Cancer Genome Atlas cohort. In our cohort, hypermethylation of MIR10B was significantly associated with the differential expression of its target genes NR4A3 and BCL2L11 (P = 0.0125 and P = 0.014, respectively), which was inversely correlated with disease-free survival (P = 9E−15 and P = 2E−15, respectively) in patients. Finally, differential methylation in FUT3, TRIM5, TSPAN7, MAP3K8, RPS6KA2, SLC9A9, and NPAS3 genes was found to be predictive of certain clinical and epidemiologic parameters.Implications: This study reveals a functional minimal methylation profile in oral tongue tumors with associated risk habits, clinical, and epidemiologic outcomes. In addition, NR4A3 downregulation and correlation with patient survival suggests a potential target for therapeutic intervention in oral tongue tumors. Data from the current study are deposited in the NCBI Geo database (accession number GSE75540). Mol Cancer Res; 14(9); 805–19. ©2016 AACR.

All the supplementary figures

List of ambiguously mapped probe indices and the number of multiply mapped positions in hg19 reference.

All the supplementary text

Multi-tiered data represented in Circos plot

Clinical information for all patients used in the study

Measurement of biological analytes, characterizing flavor in fruits, is a cumbersome, expensive and time-consuming process. Fruits with higher concentration of analytes have greater commercial or nutritional values. Here, we tested a deep learning-based framework with fruit images to predict the class (sweet or sour and high or low) of analytes using images from two types of trees in a single and multi-analyte mode. We used fruit images from kinnow (n= 3,451), an edible hybrid mandarin and neem (n= 1,045), a tree with agrochemical and pharmaceutical properties. We measured sweetness in kinnows and five secondary metabolites in neem fruits (azadirachtin or A, deacetyl-salannin or D, salannin or S, nimbin or N and nimbolide or E) using a refractometer and high-performance liquid chromatography, respectively. We trained the models for 300 epochs, before and after hyper-parameters’ evolution, using 300 generations with 50 epochs/generation, estimated the best models and evaluated their performance on 10% of independent images. The validation F1score and test accuracies were 0.79 and 0.77, and 82.55% and 60.8%, respectively for kinnow and neem A analyte. A multi-analyte model enhanced the neem A model’s prediction to ‘high’ class when the D:N:S’s combined class predictions were high:low:high and to ‘low’ class when D:N’s combined class predictions were low:high respectively. The test accuracy increased further to ~70% with a 10-fold cross-validation error of 0.257 across ten randomly split train:validation:test sets proving the potential of a multi-analyte model to enhance the prediction accuracy, especially when the numbers of images are limiting.

Molecular adaptations are prevalent in carnivorous plants in response to habitat and environmental stress. We used the chloroplast genome and characterized the specific adaptations in the photosystem genes and their interacting partners in Nepenthes khasiana, a carnivorous pitcher plant. When compared with the carnivorous and non-carnivorous groups across Caryophyllales, Lamiales, Poales, Ericales, and Oxalidales, we found Nepenthes-specific changes in psaA, psaB, psaC and psaH. Of these, only a single amino acid change each, G147 in the protein psaA and R40 in the protein ndhD, impacted the three-dimensional structural conformation of the corresponding proteins. Modeling the interaction between the psaA and the ndhD proteins identified group-specific changes between the models between Nepentheceae versus others. The least distance between the structure-impacting residues of psaA and ndhD was 25.9 Å for Nepenthes and 19.4 Å for non-Nepenthes models. Given that the chloroplast ndh and photosystem I subunits form a large super-complex with the light-harvesting carrier proteins from the nucleus to mediate cyclic electron transport, our observations may indicate specific adaptations in the cyclic electron transport arm of the photosynthetic machinery in the Nepenthes species.

Background The incidence of oral tongue squamous cell carcinoma (OTSCC) is increasing among younger birth cohorts. The etiology of early-onset OTSCC (diagnosed before the age of 50 years) and cancer driver genes remain largely unknown. Methods The Sequencing Consortium of Oral Tongue Cancer was established through the pooling of somatic mutation data of oral tongue cancer specimens (n = 227 [107 early-onset cases]) from 7 studies and The Cancer Genome Atlas. Somatic mutations at microsatellite loci and Catalog of Somatic Mutations in Cancer mutation signatures were identified. Cancer driver genes were identified with the MutSigCV and WITER algorithms. Mutation comparisons between early- and typical-onset OTSCC were evaluated via linear regression with adjustments for patient-related factors. Results Two novel driver genes (ATXN1 and CDC42EP1) and 5 previously reported driver genes (TP53, CDKN2A, CASP8, NOTCH1, and FAT1) were identified. Six recurrent mutations were identified, with 4 occurring in TP53. Early-onset OTSCC had significantly fewer nonsilent mutations even after adjustments for tobacco use. No associations of microsatellite locus mutations and mutation signatures with the age of OTSCC onset were observed. Conclusions This international, multicenter consortium is the largest study to characterize the somatic mutational landscape of OTSCC and the first to suggest differences by age of onset. This study validates multiple previously identified OTSCC driver genes and proposes 2 novel cancer driver genes. In analyses by age, early-onset OTSCC had a significantly smaller somatic mutational burden that was not explained by differences in tobacco use. Lay summary This study identifies 7 specific areas in the human genetic code that could be responsible for promoting the development of tongue cancer. Tongue cancer in young patients (under the age of 50 years) has fewer overall changes to the genetic code in comparison with tongue cancer in older patients, but the authors do not think that this is due to differences in smoking rates between the 2 groups. The cause of increasing cases of tongue cancer in young patients remains unclear.

Biofoundry is a place where biomanufacturing meets automation. The highly modular structure of a biofoundry helps accelerate the design–build–test–learn workflow to deliver products fast and in a streamlined fashion. In this perspective, we describe our efforts to build Biofoundry India, where we see the facility add a substantial value in supporting research, innovation and entrepreneurship. We describe three key areas of our focus, harnessing the potential of non-expressing parts of the sequenced genomes, using deep learning in pathway reconstruction and synthesising enzymes and metabolites. Toward the end, we describe specific challenges in building such facility in India and the path to mitigate some of those working with the other biofoundries worldwide.

This is the open peer reviewers comments and recommendations regarding the submitted GigaScience article and/or dataset.

Development of rice cultivars bearing numerous spikelets by breeding approach to increase the yearly production of rice to approximately 800 million metric tons to feed the ever increasing population of the world accompanies poor grain filling in the inferior spikelets preventing achievement of the yield potential. As the initial stages of caryopses development are of much importance for grain filling, spatio-temporal expressions of the miRNAs were studied during these periods in the spikelets of a compact-panicle rice cultivar, Oryza sativa cv. Mahalaxmi, bearing numerous spikelets per panicle to understand the reason of poor grain filling at the level of the initial biochemical events. Differential expression of several known miRNAs between the superior and inferior spikelets suggested great difference in metabolism related to grain filling in the spikelets based on their spatial location on compact panicle. Expressions of five known and four novel miRNAs were validated by Northern. Their targets included the enzymes directly involved in starch biosynthesis like sucrose synthase, starch synthase and pullulanase, besides others. Spatio-temporal expression studies of these miRNAs in the spikelets of Mahalaxmi revealed a pattern of mostly a greater expression in the inferior spikelets compared with the superior ones concomitant with an inverse expression of the target genes, which was not observed in the lax-panicle cultivar Upahar. The study thus revealed that the grain filling in rice is greatly regulated by miRNAs, and these miRNAs or their target genes could be considered for biotechnological interventions for improving grain filling in the rice cultivars of interest.

Purpose Accurate detection of human papillomavirus (HPV) in oral cavity squamous cell carcinoma (OSCC) is essential to understanding the role of HPV in disease prognosis and management of patients. We used different analytes and methods to understand the true prevalence of HPV in a cohort of patients with OSCC with different molecular backgrounds, and we correlated HPV data with patient survival. Methods We integrated data from multiple analytes (HPV DNA, HPV RNA, and p16), assays (immunohistochemistry, polymerase chain reaction [PCR], quantitative PCR [qPCR], and digital PCR), and molecular changes (somatic mutations and DNA methylation) from 153 patients with OSCC to correlate p16 expression, HPV DNA, and HPV RNA with HPV incidence and patient survival. Results High prevalence (33% to 58%) of HPV16/18 DNA did not correlate with the presence of transcriptionally active viral genomes (15%) in tumors. Eighteen percent of the tumors were p16 positive and only 6% were both HPV DNA and HPV RNA positive. Most tumors with relatively high copy number HPV DNA and/or HPV RNA, but not with HPV DNA alone (irrespective of copy number), were wild-type for TP53 and CASP8 genes. In our study, p16 protein, HPV DNA, and HPV RNA, either alone or in combination, did not correlate with patient survival. Nine HPV-associated genes stratified the virus-positive from the virus-negative tumor group with high confidence ( P < .008) when HPV DNA copy number and/or HPV RNA were considered to define HPV positivity, and not HPV DNA alone, irrespective of copy number ( P < .2). Conclusion In OSCC, the presence of both HPV RNA and p16 is rare. HPV DNA alone is not an accurate measure of HPV positivity and therefore may not be informative. HPV DNA, HPV RNA, and p16 do not correlate with patients’ outcome.

Tumor suppression by the extracts of Azadirachta indica (neem) works via anti-proliferation, cell cycle arrest, and apoptosis, demonstrated previously using cancer cell lines and live animal models. However, very little is known about the molecular targets and pathways that the neem extracts and the associated compounds act through. Here, we address this using a genome-wide functional pooled shRNA screen on head and neck squamous cell carcinoma cell line treated with crude neem leaf extracts, known for their anti-tumorigenic activity. By analyzing differences in global clonal sizes of the shRNA-infected cells cultured under no treatment and treatment with neem leaf extract conditions, assayed using next-generation sequencing, we found 225 genes affected the cancer cell growth in the shRNA-infected cells treated with neem extract. Pathway enrichment analyses of whole-genome gene expression data from cells temporally treated with neem extract revealed important roles played by the TGF-β pathway and HSF-1-related gene network. Our results indicate that neem extract simultaneously affects various important molecular signaling pathways in head and neck cancer cells, some of which may be therapeutic targets for this devastating tumor.

BackgroundCAFE MOCHA(Clinical Association of Functionally Established MOlecular CHAnges) is an integrated GUI-driven computational and statistical framework to discover molecular signatures linked to a specific clinical attribute in a cancer type. We testedCAFE MOCHAin head and neck squamous cell carcinoma (HNSCC) for discovering a signature linked to distant metastasis and recurrence (MR) in 517 tumors from TCGA and validated the signature in 18 tumors from an independent cohort.MethodsThe platform integrates mutations and indels, gene expression, DNA methylation and copy number variations to discover a classifier first, predict an incoming tumour for the same by pulling defined class variables into a single framework that incorporates a coordinate geometry-based algorithm, called Complete Specificity Margin Based Clustering (CSMBC) with 100% specificity.CAFE MOCHAclassifies an incoming tumour sample using either a matched normal or a built-in database of normal tissues. The application is packed and deployed using theinstall4jmulti-platform installer.ResultsWe testedCAFE MOCHAto discover a signature for distant metastasis and recurrence in HNSCC. The signature MR44 in HNSCC yielded 80% sensitivity and 100% specificity in the discovery stage and 100% sensitivity and 100% specificity in the validation stage.ConclusionsCAFE MOCHAis a cancer type- and clinical attribute-agnostic computational and statistical framework to discover integrated molecular signature for a specific clinical attribute.CAFE MOCHAis available in GitHub (https://github.com/binaypanda/CAFEMOCHA).

Oral tongue squamous cell carcinomas (OTSCC) are a homogenous group of aggressive tumors in the head and neck region that spread early to lymph nodes and have a higher incidence of regional failure. In addition, there is a rising incidence of oral tongue cancer in younger populations. Studies on functional DNA methylation changes linked with altered gene expression are critical for understanding the mechanisms underlying tumor development and metastasis. Such studies also provide important insight into biomarkers linked with viral infection, tumor metastasis, and patient survival in OTSCC. Therefore, we performed genome-wide methylation analysis of tumors (N = 52) and correlated altered methylation with differential gene expression. The minimal tumor-specific DNA 5-methylcytosine signature identified genes near 16 different differentially methylated regions, which were validated using genomic data from The Cancer Genome Atlas cohort. In our cohort, hypermethylation of MIR10B was significantly associated with the differential expression of its target genes NR4A3 and BCL2L11 (P = 0.0125 and P = 0.014, respectively), which was inversely correlated with disease-free survival (P = 9E−15 and P = 2E−15, respectively) in patients. Finally, differential methylation in FUT3, TRIM5, TSPAN7, MAP3K8, RPS6KA2, SLC9A9, and NPAS3 genes was found to be predictive of certain clinical and epidemiologic parameters. Implications: This study reveals a functional minimal methylation profile in oral tongue tumors with associated risk habits, clinical, and epidemiologic outcomes. In addition, NR4A3 downregulation and correlation with patient survival suggests a potential target for therapeutic intervention in oral tongue tumors. Data from the current study are deposited in the NCBI Geo database (accession number GSE75540). Mol Cancer Res; 14(9); 805–19. ©2016 AACR.

BackgroundSelection of the right reference gene(s) is crucial in the analysis and interpretation of gene expression data. In head and neck cancer, studies evaluating the efficacy of internal reference genes are rare. Here, we present data for a minimal set of candidates as internal control genes for gene expression studies in head and neck cancer.MethodsWe analyzed data from multiple sources (in house whole-genome gene expression microarrays, n=21; TCGA RNA-seq, n=42, and published gene expression studies in head and neck tumors from literature) to come up with a set of genes (discovery set) for their stable expression across tumor and normal tissues. We then performed independent validation of their expression using qPCR in 14 tumor:normal pairs. Genes in the discovery set were ranked using four different algorithms (BestKeeper, geNorm, NormFinder, and comparative delta Ct) and a web-based comparative tool, RefFinder, for their stability and variance in expression across tissues.ResultsOur analyses resulted in 18 genes (discovery set) that had lowest variance and high level of expression across tumor and normal samples. Independent experimental validation and analyses with multiple tools resulted in top ranked five genes (RPL30, RPL27, PSMC5, OAZ1 and MTCH1) out of which, RPL30 (60S ribosomal protein L30) and RPL27 (60S ribosomal protein L27), performed best and were abundantly expressed across tumor and normal tissues.ConclusionsRPL30 and RPL27 are stably expressed in HNSCC and should be used as internal control genes in gene expression in head and neck tumors studies.

Background Poor filling of grains in the basal spikelets of large size panicles bearing numerous spikelets has been a major limitation in attempts to increase the rice production to feed the world’s increasing population. Considering that biotechnological intervention could play important role in overcoming this limitation, the role of cytokinin in grain filling was investigated based on the information on cell proliferating potential of the hormone and reports of its high accumulation in immature seeds. Results A comparative study considering two rice varieties differing in panicle compactness, lax-panicle Upahar and compact-panicle OR-1918, revealed significant difference in grain filling, cytokinin oxidase (CKX) activity and expression, and expression of cell cycle regulators and cytokinin signaling components between the basal and apical spikelets of OR-1918, but not of Upahar. Exogenous application of cytokinin (6-Benzylaminopurine, BAP) to OR-1918 improved grain filling significantly, and this was accompanied by a significant decrease in expression and activity of CKX, particularly in the basal spikelets where the activity of CKX was significantly higher than that in the apical spikelets. Cytokinin application also resulted in significant increase in expression of cell cycle regulators like cyclin dependent kinases and cyclins in the basal spikelets that might be facilitating cell division in the endosperm cells by promoting G1/S phase and G2/M phase transition leading to improvement in grain filling. Expression studies of type-A response regulator (RR) component of cytokinin signaling indicated possible role of OsRR3, OsRR4 and OsRR6 as repressors of CKX expression, much needed for an increased accumulation of CK in cells. Furthermore, the observed effect of BAP might not be solely because of it, but also because of induced synthesis of trans-zeatin (tZ) and N6-(Δ2-isopentenyl)adenine (iP), as reflected from accumulation of tZR (tZ riboside) and iPR (iP riboside), and significantly enhanced expression of an isopentenyl transferase (IPT) isoform. Conclusion The results suggested that seed-specific overexpression of OsRR4 and OsRR6, and more importantly of IPT9 could be an effective biotechnological intervention towards improving the CK level of the developing caryopses leading to enhanced grain filling in rice cultivars bearing large panicles with numerous spikelets, and thereby increasing their yield potential. Electronic supplementary material The online version of this article (10.1186/s12870-018-1279-4) contains supplementary material, which is available to authorized users.

Availability of snake genome sequences has opened up exciting areas of research on comparative genomics and gene diversity. One of the challenges in studying snake genomes is the acquisition of biological material from live animals, especially from the venomous ones. Additionally, in certain countries, Government permission is required to handle live snakes, making the process cumbersome and time-consuming. Here, we report comparative sequence analyses of toxin gene homologs from Russells viper (Daboia russelii) using whole-genome sequencing data obtained from the shed skin. When compared with the major venom proteins in Russells viper studied previously, we found 45-100% sequence similarity between the venom proteins and their skin homologs. Additionally, comparative analyses of 20 toxin gene family homologs provided evidence of unique sequence motifs in nerve growth factor (NGF), platelet derived growth factor (PDGF), Kunitz/Bovine pancreatic trypsin inhibitor (Kunitz BPTI), cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins (CAP) and cysteine-rich secretory protein (CRISP). We identified V11 and T35 in the NGF domain; F23 and A29 in the PDGF domain; N69, K2 and A5 in the CAP domain; and Q17 in the CRISP domain to be responsible for differences in the largest pockets across the protein domain structures in crotalines, viperines and elapids from the in silico structure-based analysis. Similarly, residues F10, Y11 and E20 appear to play an important role in the protein structures across the kunitz protein domain of viperids and elapids. Our study sheds light on the usefulness of studying venom protein homologs from skin, their unique features and evolution in vipers. Data deposition: Russells viper sequence data is deposited in the NCBI SRA database under the accession number SRR5506741 and sequences for the individual venom-associated gene homologs to GenBank (accession numbers in Table S1).

Although a common cause of disease, copy number variants (CNVs) have not routinely been identified from next-generation sequencing (NGS) data in a clinical context. This study aimed to examine the sensitivity and specificity of a widely used software package, ExomeDepth, to identify CNVs from targeted NGS data sets. We benchmarked the accuracy of CNV detection using ExomeDepth v1.1.6 applied to targeted NGS data sets by comparison to CNV events detected through whole-genome sequencing for 25 individuals and determined the sensitivity and specificity of ExomeDepth applied to these targeted NGS data sets to be 100% and 99.8%, respectively. To define quality assurance metrics for CNV surveillance through ExomeDepth, we undertook simulation of single-exon (n=1000) and multiple-exon heterozygous deletion events (n=1749), determining a sensitivity of 97% (n=2749). We identified that the extent of sequencing coverage, the inter- and intra-sample variability in the depth of sequencing coverage and the composition of analysis regions are all important determinants of successful CNV surveillance through ExomeDepth. We then applied these quality assurance metrics during CNV surveillance for 140 individuals across 12 distinct clinical areas, encompassing over 500 potential rare disease diagnoses. All 140 individuals lacked molecular diagnoses after routine clinical NGS testing, and by application of ExomeDepth, we identified 17 CNVs contributing to the cause of a Mendelian disorder. Our findings support the integration of CNV detection using ExomeDepth v1.1.6 with routine targeted NGS diagnostic services for Mendelian disorders. Implementation of this strategy increases diagnostic yields and enhances clinical care.

Epigenetics is a field that encompasses chemical modifications of DNA and histone proteins, both of which alter gene expression without changing the underlying nucleotide sequence. DNA methylation and modifications of histone tails have been studied in detail and are now known to be global gene regulatory mechanisms. An analogous post-transcriptional modification is chemical modification of specific nucleotides in RNA. Study of RNA modifications is a nascent field as yet, and the significance of these marks in controlling cell growth and differentiation is just beginning to be appreciated. The addition of a methyl group to adenosine (N-methyl-6-adenosine) or m6A is the most abundant modification in mammalian mRNAs. Though identified four decades ago, interest in this particular modification was set off by the discovery that the obesity gene FTO was an RNA demethylase. Since then, many studies have investigated m6A modification in different species. In this review, we summarize the current literature and hypotheses about the presence and function of this ubiquitous RNA modification in mammals, viruses, yeast and plants in terms of the consensus sequence and the methyltransferase/demethylation machinery identified thus far. We discuss its potential role in regulating molecular and physiological processes in each of these organisms, especially its role in RNA splicing, RNA degradation and development. We also enlist the methodologies developed so far, both locus-specific and transcriptome-wide, to study this modification. Lastly, we discuss whether m6A alterations have consequences on modulating disease aetiology, and speculate about its potential role in cancer.

Grain amaranths, edible C4 dicots, produce pseudo-cereals high in lysine. Lysine being one of the most limiting essential amino acids in cereals and C4 photosynthesis being one of the most sought-after phenotypes in protein-rich legume crops, the genome of one of the grain amaranths is likely to play a critical role in crop research. We have sequenced the genome and transcriptome of Amaranthus hypochondriacus, a diploid (2n = 32) belonging to the order Caryophyllales with an estimated genome size of 466 Mb. Of the 411 linkage single-nucleotide polymorphisms (SNPs) reported for grain amaranths, 355 SNPs (86%) are represented in the scaffolds and 74% of the 8.6 billion bases of the sequenced transcriptome map to the genomic scaffolds. The genome of A. hypochondriacus, codes for at least 24,829 proteins, shares the paleohexaploidy event with species under the superorders Rosids and Asterids, harbours 1 SNP in 1,000 bases, and contains 13.76% of repeat elements. Annotation of all the genes in the lysine biosynthetic pathway using comparative genomics and expression analysis offers insights into the high-lysine phenotype. As the first grain species under Caryophyllales and the first C4 dicot genome reported, the work presented here will be beneficial in improving crops and in expanding our understanding of angiosperm evolution.

RNA sequencing (RNA-seq) is a powerful technology for identification of novel transcripts (coding, non-coding and splice variants), understanding of transcript structures and estimation of gene and/or allelic expression. There are specific challenges that biologists face in determining the number of replicates to use, total number of sequencing reads to generate for detecting marginally differentially expressed transcripts and the number of lanes in a sequencing flow cell to use for the production of right amount of information. Although past studies attempted answering some of these questions, there is a lack of accessible and biologist-friendly mobile applications to answer these questions. Keeping this in mind, we have developed RNAtor, a mobile application for Android platforms, to aid biologists in correctly designing their RNA-seq experiments. The recommendations from RNAtor are based on simulations and real data.Availability and ImplementationThe Android version of RNAtor is available on Google Play Store and the code from GitHub (https://github.com/binaypanda/RNAtor).

Sample processing protocols that enable compatible recovery of differentially expressed (DE) transcripts and proteins are necessary for integration of the multiomics data applied in the analysis of tumors. In this pilot study, we compared two different isolation methods for extracting RNA and protein from laryngo-pharyngeal tumor tissues and the corresponding adjacent normal sections. In Method 1, RNA and protein were isolated from a single tissue section sequentially and in Method 2, the extraction was carried out using two different sections and two independent and parallel protocols for RNA and protein. RNA and protein from both methods were subjected to RNA-seq and iTRAQ-based LC-MS/MS analysis, respectively. Analysis of data revealed that a higher number of differentially expressed transcripts and proteins were concordant in their regulation trends in Method 1 as compared to Method 2. Cross method comparison of concordant entities revealed that RNA and protein extraction from the same tissue section (Method 1) recovered more concordant entities that are missed in the other extraction method (Method 2) indicating heterogeneity in distribution of these entities in different tissue sections. Method 1 could thus be the method of choice for integrated analysis of transcriptome and proteome data. This article is protected by copyright. All rights reserved

Background Inherited retinal diseases (IRDs) are a clinically and genetically heterogeneous set of disorders, for which diagnostic second-generation sequencing (next-generation sequencing, NGS) services have been developed worldwide. Methods We present the molecular findings of 537 individuals referred to a 105-gene diagnostic NGS test for IRDs. We assess the diagnostic yield, the spectrum of clinical referrals, the variant analysis burden and the genetic heterogeneity of IRD. We retrospectively analyse disease-causing variants, including an assessment of variant frequency in Exome Aggregation Consortium (ExAC). Results Individuals were referred from 10 clinically distinct classifications of IRD. Of the 4542 variants clinically analysed, we have reported 402 mutations as a cause or a potential cause of disease in 62 of the 105 genes surveyed. These variants account or likely account for the clinical diagnosis of IRD in 51% of the 537 referred individuals. 144 potentially disease-causing mutations were identified as novel at the time of clinical analysis, and we further demonstrate the segregation of known disease-causing variants among individuals with IRD. We show that clinically analysed variants indicated as rare in dbSNP and the Exome Variant Server remain rare in ExAC, and that genes discovered as a cause of IRD in the post-NGS era are rare causes of IRD in a population of clinically surveyed individuals. Conclusions Our findings illustrate the continued powerful utility of custom-gene panel diagnostic NGS tests for IRD in the clinic, but suggest clear future avenues for increasing diagnostic yields.

The advent of high-throughput technology has revolutionized biological sciences in the last two decades enabling experiments on the whole genome scale. Data from such large-scale experiments are interpreted at system’s level to understand the interplay among genome, transcriptome, epigenome, proteome, metabolome, and regulome.

Head and neck squamous cell carcinoma (HNSCC) are a heterogeneous group of cancers affecting multiple subsites, including oral cavity. Oral or anterior tongue tumors (OTSCC) are an aggressive group of squamous cell carcinomas, characterized by their early spread to lymph nodes and higher rate of regional failure compared to other oral cavity cancers. There is a rise in the incidence of oral tongue cancer among younger population (

Neem (Azadirachta indica A. Juss.), an evergreen tree of the Meliaceae family, is known for its medicinal, cosmetic, pesticidal and insecticidal properties. We had previously sequenced and published the draft genome of the plant, using mainly short read sequencing data. In this report, we present an improved genome assembly generated using additional short reads from Illumina and long reads from Pacific Biosciences SMRT sequencer. We assembled short reads and error corrected long reads using Platanus, an assembler designed to perform well for heterozygous genomes. The updated genome assembly (v2.0) yielded 3- and 3.5-fold increase in N50 and N75, respectively; 2.6-fold decrease in the total number of scaffolds; 1.25-fold increase in the number of valid transcriptome alignments; 13.4-fold less mis-assembly and 1.85-fold increase in the percentage repeat, over the earlier assembly (v1.0). The current assembly also maps better to the genes known to be involved in the terpenoid biosynthesis pathway. Together, the data represents an improved assembly of the A. indica genome. The raw data described in this manuscript are submitted to the NCBI Short Read Archive under the accession numbers SRX1074131, SRX1074132, SRX1074133, and SRX1074134 (SRP013453).

Oral tongue squamous cell carcinomas (OTSCC) are a homogeneous group of tumors characterized by aggressive behavior, early spread to lymph nodes and a higher rate of regional failure. Additionally, the incidence of OTSCC among younger population (

High throughput molecular profiling and integrated data analysis with tumor tissues require overcoming challenges like tumor heterogeneity and tissue paucity. This study is an attempt to understand and optimize various steps during tissue processing and in establishing pipelines essential for integrated analysis. Towards this effort, we subjected laryngo-pharyngeal primary tumors and the corresponding adjacent normal tissues (n=2) to two RNA and protein isolation methods, one wherein RNA and protein were isolated from the same tissue sequentially (Method 1) and second, wherein the extraction was carried out using two independent methods (Method 2). RNA and protein from both methods were subjected to RNA-seq and iTRAQ based LC-MS/MS analysis. Transcript and peptide identification and quantification was followed by both individual -ome and integrated data analysis. As a result of this analysis, we identified a higher number of total, as well as differentially expressed (DE) transcripts (1329 vs 1134) and proteins (799 vs 408) with fold change ≥ 2.0, in Method 1. Among these, 173 and 86 entities were identified by both transcriptome and proteome analysis in Method 1 and 2, respectively, with higher concordance in the regulation trends observed in the former. The significant cancer related pathways enriched with the individual DE transcript or protein data were similar in both the methods. However, the entities mapping to them were different, allowing enhanced view of the pathways identified after integration of the data and subsequent mapping. The concordant DE transcripts and proteins also revealed key molecules of the pathways with important roles in cancer development. This study thus demonstrates that sequential extraction of the RNA and proteins from the same tissue allows for better profiling of differentially expressed entities and a more accurate integrated data analysis.

Breeding programs for increasing spikelet number in rice have resulted in compactness of the panicle, accompanied by poor grain filling in inferior spikelets. Although the inefficient utilization of assimilate has been indicated as responsible for this poor grain filling, the underlying cause remains elusive. The current study utilized the suppression subtractive hybridization technique to identify 57 and 79 genes that overexpressed in the superior and inferior spikelets (with respect to each other), respectively, of the compact-panicle rice cultivar Mahalaxmi. Functional categorization of these differentially expressed genes revealed a marked metabolic difference between the spikelets according to their spatial location on the panicle. The expression of genes encoding seed storage proteins was dominant in inferior spikelets, whereas genes encoding regulatory proteins, such as serine-threonine kinase, zinc finger protein and E3 ligase, were highly expressed in superior spikelets. The expression patterns of these genes in the inferior and superior spikelets of Mahalaxmi were similar to those observed in another compact-panicle cultivar, OR-1918, but differed from those obtained in two lax-panicle cultivars, Upahar and Lalat. The results first suggest that the regulatory proteins abundantly expressed in the superior spikelets of compact-panicle cultivars and in both the superior and inferior spikelets of lax-panicle cultivars but poorly expressed in the inferior spikelets of compact-panicle cultivars promote grain filling. Second, the high expression of seed-storage proteins observed in the inferior spikelets of compact-panicle cultivars appears to inhibit the grain filling process. Third, the low expression of enzymes of the Krebs cycle in inferior spikelets compared with superior spikelets of compact-panicle cultivars is bound to lead to poor ATP generation in the former and consequently limit starch biosynthesis, an ATP-consuming process, resulting in poor grain filling.