GET_PANGENES: calling pangenes from plant genome alignments confirms presence-absence variation.

Crop pangenomes made from individual cultivar assemblies promise easy access to conserved genes, but genome content variability and inconsistent identifiers hamper their exploration. To address this, we define pangenes, which summarize a species coding potential and link back to original annotations. The protocol get_pangenes performs whole genome alignments (WGA) to call syntenic gene models based on coordinate overlaps. A benchmark with small and large plant genomes shows that pangenes recapitulate phylogeny-based orthologies and produce complete soft-core gene sets. Moreover, WGAs support lift-over and help confirm gene presence-absence variation. Source code and documentation: https://github.com/Ensembl/plant-scripts .

Development of a Low-cost NGS Test for the Evaluation of Thyroid Nodules.

Ultrasound-guided fine needle aspiration cytology (FNAC) is the preferred method of identifying malignancy in palpable thyroid nodules using the Bethesda reporting system. However, in around 30-40% of FNACs (Bethesda categories III, IV, and V), the results are indeterminate and surgery is required to confirm malignancy. Out of those who undergo surgery, only 10-40% of patients in these categories are found to have malignancies, thus proving surgery to be unnecessary for some patients or to be incomplete in others. While molecular testing on thyroid FNAC material is part of the American Thyroid Association (ATA) guidelines in evaluating thyroid nodules, it is currently unavailable in India due to cost constraints. In this study, we prospectively collected FNAC samples from sixty-nine patients who presented with palpable thyroid nodules. We designed a cost-effective next-generation sequencing (NGS) test to query multiple variants in the DNA and RNA isolated from the fine needle aspirate. The identification of oncogenic variants was considered to be indicative of malignancy, and confirmed by surgical histopathology. The panel showed an overall sensitivity of 81.25% and a specificity of 100%, while in the case of Bethesda categories III, IV, and V, the sensitivity was higher (87.5%) and the specificity was established at 100%. The panel could thereby serve as a rule-in test for the diagnosis of thyroid cancer and therefore help identify patients who require surgery, especially in the indeterminate Bethesda categories III, IV, and V.

Determining Cost-Optimal Next-Generation Sequencing Panels for Rare Disease and Pharmacogenomics Testing.

BACKGROUND: Multi-gene panel sequencing using next-generation sequencing (NGS) methods is a key tool for genomic medicine. However, with an estimated 140 000 genomic tests available, current system inefficiencies result in high genetic-testing costs. Reduced testing costs are needed to expand the availability of genomic medicine. One solution to improve efficiency and lower costs is to calculate the most cost-effective set of panels for a typical pattern of test requests. METHODS: We compiled rare diseases, associated genes, point prevalence, and test-order frequencies from a representative laboratory. We then modeled the costs of the relevant steps in the NGS process in detail. Using a simulated annealing-based optimization procedure, we determined panel sets that were more cost-optimal than whole exome sequencing (WES) or clinical exome sequencing (CES). Finally, we repeated this methodology to cost-optimize pharmacogenomics (PGx) testing. RESULTS: For rare disease testing, we show that an optimal choice of 4-6 panels, uniquely covering genes that comprise 95% of the total prevalence of monogenic diseases, saves $257-304 per sample compared with WES, and $66-135 per sample compared with CES. For PGx, we show that the optimal multipanel solution saves $6-7 (27%-40%) over a single panel covering all relevant gene-drug associations. CONCLUSIONS: Laboratories can reduce costs using the proposed method to obtain and run a cost-optimal set of panels for specific test requests. In addition, payers can use this method to inform reimbursement policy.

Differentially Expressed MicroRNAs Link Cellular Physiology to Phenotypic Changes in Rice Under Stress Conditions.

Plant microRNAs (miRNAs) and their target genes have important functional roles in nutrition deficiency and stress response. However, the underlying mechanisms relating relative expression of miRNAs and target mRNAs to morphological adjustments are not well defined. By combining miRNA expression profiles, corresponding target genes and transcription factors that bind to computationally identified over-represented cis-regulatory elements (CREs) common in miRNAs and target gene promoters, we implement a strategy that identifies a set of differentially expressed regulatory interactions which, in turn, relate underlying cellular mechanisms to some of the phenotypic changes observed. Integration of experimentally reported individual interactions with identified regulatory interactions explains how (i) during mineral deficiency osa-miR167 inhibits shoot growth but activates adventitious root growth by influencing free auxin content; (ii) during sulfur deficiency osa-miR394 is involved in adventitious root growth inhibition, sulfur and iron homeostasis, and auxin-mediated regulation of sulfur homeostasis; (iii) osa-miR399 contributes to cross-talk between cytokinin and phosphorus deficiency signaling; and (iv) a feed-forward loop involving the osa-miR166, trihelix and HD-ZIP III transcription factors may regulate leaf senescence during drought. This strategy not only identifies various regulatory interactions connecting phenotypic changes with cellular or molecular events triggered by stress, but also provides a framework to deepen our understanding of stress cellular physiology.

Identification and validation of a virus-inducible ta-siRNA-generating TAS4 locus in tomato.

Trans-acting small interfering RNAs (ta-siRNAs) are a class of endogenous small RNA, associated with posttranscriptional gene silencing. Their biogenesis requires an initial microRNA (miRNA)-mediated cleavage of precursor RNA. Around 20 different ta-siRNA-producing loci (TASs), whose sequences are conserved, are reported in plants. In tomato, two TAS gene families have been identified, which are found to target auxin response factor gene and bacterial spot disease resistance protein Bs4 gene. Using high-throughput computational and experimental approach, we identified a new locus-producing ta-siRNA in tomato. We have also identified the putative miRNA regulating the production of ta-siRNA from this locus. The ta-siRNAs generated from TAS4 were up-regulated upon infection with a DNA virus. The potential targets of ta-siRNAs were predicted to be variety of proteins including MYB transcription factors and cell cycle regulators for some of the ta-siRNAs produced.

miRMOD: a tool for identification and analysis of 5' and 3' miRNA modifications in Next Generation Sequencing small RNA data.

In the past decade, the microRNAs (miRNAs) have emerged to be important regulators of gene expression across various species. Several studies have confirmed different types of post-transcriptional modifications at terminal ends of miRNAs. The reports indicate that miRNA modifications are conserved and functionally significant as it may affect miRNA stability and ability to bind mRNA targets, hence affecting target gene repression. Next Generation Sequencing (NGS) of the small RNA (sRNA) provides an efficient and reliable method to explore miRNA modifications. The need for dedicated software, especially for users with little knowledge of computers, to determine and analyze miRNA modifications in sRNA NGS data, motivated us to develop miRMOD. miRMOD is a user-friendly, Microsoft Windows and Graphical User Interface (GUI) based tool for identification and analysis of 5' and 3' miRNA modifications (non-templated nucleotide additions and trimming) in sRNA NGS data. In addition to identification of miRNA modifications, the tool also predicts and compares the targets of query and modified miRNAs. In order to compare binding affinities for the same target, miRMOD utilizes minimum free energies of the miRNA:target and modified-miRNA:target interactions. Comparisons of the binding energies may guide experimental exploration of miRNA post-transcriptional modifications. The tool is available as a stand-alone package to overcome large data transfer problems commonly faced in web-based high-throughput (HT) sequencing data analysis tools. miRMOD package is freely available at http://bioinfo.icgeb.res.in/miRMOD.

3' and 5' microRNA-end post-biogenesis modifications in plant transcriptomes: Evidences from small RNA next generation sequencing data analysis.

The processing of miRNA from its precursors is a precisely regulated process and after biogenesis, the miRNAs are amenable to different kinds of modifications by the addition or deletion of nucleotides at the terminal ends. However, the mechanism and functions of such modifications are not well studied in plants. In this study, we have specifically analysed the terminal end non-templated miRNA modifications, using NGS data of rice, tomato and Arabidopsis small RNA transcriptomes from different tissues and physiological conditions. Our analysis reveals template independent terminal end modifications in the mature as well as passenger strands of the miRNA duplex. Interestingly, it is also observed that miRNA sequences terminating with a cytosine (C) at the 3' end undergo a higher percentage of 5' end modifications. The terminal end modifications did not correlate with the miRNA abundances and are independent of tissue types, physiological conditions and plant species. Our analysis indicates that the addition of nucleotides at miRNA ends is not influenced by the absence of RNA dependent RNA polymerase 6. Moreover the terminal end modified miRNAs are also observed amongst AGO1 bound small RNAs and have potential to alter target, indicating its important functional role in repression of gene expression.

Prediction and characterization of Tomato leaf curl New Delhi virus (ToLCNDV) responsive novel microRNAs in Solanum lycopersicum.

Tomato leaf curl New Delhi virus (ToLCNDV) infects tomato (Solanum lycopersicum) plants and causes severe crop losses. As the microRNAs (miRNAs) are deregulated during stressful events, such as biotic stress, we wanted to study the effect of ToLCNDV infection on tomato miRNAs. We constructed two libraries, isolating small RNAs (sRNAs) from healthy (HT) and ToLCNDV infected (IT) tomato leaves, and sequenced the library-specific sRNAs using the next generation sequencing (NGS) approach. These data helped predict 112 mature miRNA sequences employing the miRDeep-P program. A substantial number (58) of the sequences were 24-mer in size, which was a bit surprising. Based on the calculation of precision values, 53 novel miRNAs were screened from the predicted sequences. Nineteen of these were chosen for expression analysis; a northern blot analysis showed 15 to be positive. Many of the predicted miRNAs were up-regulated following viral infection. The target genes of the miRNAs were also predicted and the expression analysis of selected transcripts showed a typical inverse relation between the accumulation of target transcripts and the abundance of corresponding miRNAs. Furthermore, the cleavage sites of the target transcripts for three novel miRNAs were mapped, confirming the correct annotation of the miRNA-targets. The sRNA deep sequencing clearly revealed that the virus modulated global miRNA expression in the host. The validated miRNAs (Tom_4; Tom_14; Tom_17; Tom_21; Tom_29; Tom_43) could be valuable tools for understanding the ToLCNDV-tomato interaction, ultimately leading to the development of a virus-resistant tomato plant.