Impact of Library Preparation on Downstream Analysis and Interpretation of RNA-Seq Data: Comparison between Illumina PolyA and NuGEN Ovation Protocol

Objectives

The sequencing by the PolyA selection is the most common approach for library preparation. With limited amount or degraded RNA, alternative protocols such as the NuGEN have been developed. However, it is not yet clear how the different library preparations affect the downstream analyses of the broad applications of RNA sequencing.

Methods and Materials

Eight human mammary epithelial cell (HMEC) lines with high quality RNA were sequenced by Illumina’s mRNA-Seq PolyA selection and NuGEN ENCORE library preparation. The following analyses and comparisons were conducted: 1) the numbers of genes captured by each protocol; 2) the impact of protocols on differentially expressed gene detection between biological replicates; 3) expressed single nucleotide variant (SNV) detection; 4) non-coding RNAs, particularly lincRNA detection; and 5) intragenic gene expression.

Results

Sequences from the NuGEN protocol had lower (75%) alignment rate than the PolyA (over 90%). The NuGEN protocol detected fewer genes (12–20% less) with a significant portion of reads mapped to non-coding regions. A large number of genes were differentially detected between the two protocols. About 17–20% of the differentially expressed genes between biological replicates were commonly detected between the two protocols. Significantly higher numbers of SNVs (5–6 times) were detected in the NuGEN samples, which were largely from intragenic and intergenic regions. The NuGEN captured fewer exons (25% less) and had higher base level coverage variance. While 6.3% of reads were mapped to intragenic regions in the PolyA samples, the percentages were much higher (20–25%) for the NuGEN samples. The NuGEN protocol did not detect more known non-coding RNAs such as lincRNAs, but targeted small and “novel” lincRNAs.

Conclusion

Different library preparations can have significant impacts on downstream analysis and interpretation of RNA-seq data. The NuGEN provides an alternative for limited or degraded RNA but it has limitations for some RNA-seq applications.     

Molecular dynamics simulation of complex Histones Deacetylase (HDAC) Class II Homo Sapiens with suberoylanilide hydroxamic acid (SAHA) and its derivatives as inhibitors of cervical cancer

Cervical cancer is second most common cancer in woman worldwide. Cervical cancer caused by human papillomavirus (HPV) oncogene. Inhibition of histone deacetylase (HDAC) activity has been known as a potential strategy for cancer therapy. SAHA is an HDAC inhibitor that has been used in cancer therapy but still has side effects. SAHA modification proposed to minimize side effects. Triazole attachment on the chain of SAHA has been known to enhance the inhibition ability of SAHA and less toxic. In this study, it will be carried out with molecular dynamic simulations of SAHA modifications consisting ligand 1a, 2a and, 2c to interact with six HDAC in hydrated conditions. To all six HDAC Class II, performed docking with SAHA and a modified inhibitor. The docking results were then carried out molecular dynamics simulations to determine the inhibitor affinities in hydrated conditions. The molecular dynamic simulations results show better affinities of ligand 2c with HDAC 4, 6, and 7 than SAHA itself, and good affinity was also shown by ligand 2a and 1c on HDAC 5 and 9. The results of this study can be a reference to obtain better inhibitors.     

Cytotoxic T-lymphocyte elicited therapeutic vaccine candidate targeting cancer against MAGE-A11 carcinogenic protein

Immunotherapy is a breakthrough approach for cancer treatment and prevention. By exploiting the fact that cancer cells have overexpression of tumor antigens responsible for its growth and progression, which can be identified and removed by boosting the immune system. In silico techniques have provided efficient ways for developing preventive measures to ward off cancer. Herein, we have designed a potent cytotoxic T-lymphocyte epitope to elicit a desirable immune response against carcinogenic melanoma-associated antigen-A11. Potent epitope was predicted using reliable algorithms and characterized by advanced computational avenue CABS molecular dynamics simulation, for full flexible binding with HLA-A*0201 and androgen receptor to large-scale rearrangements of the complex system. Results showed the potent immunogenic construct (KIIDLVHLL), from top epitopes using five algorithms. Molecular docking analyses showed the strong binding of epitope with HLA-A*0201 and androgen receptor with docking score of −780.6 and −641.06 kcal/mol, respectively. Molecular dynamics simulation analysis revealed strong binding of lead epitope with androgen receptor by involvement of 127 elements through atomic-model study. Full flexibility study showed stable binding of epitope with an average root mean square deviation (RMSD) 2.21 Å and maximum RMSD value of 6.48 Å in optimal cluster density area. The epitope also showed remarkable results with radius of gyration 23.0777 Å, world population coverage of 39.08% by immune epitope database, and transporter associated with antigen processing (TAP) affinity IC₅₀ value of 2039.65 nm. Moreover, in silico cloning approach confirmed the expression and translation capacity of the construct within a suitable expression vector. The present study paves way for a potential immunogenic construct for prevention of cancer.     

Gibberellic Acid-Priming Promotes Fluoride Tolerance in a Susceptible Indica Rice Cultivar by Regulating the Antioxidant and Phytohormone Homeostasis

Journal of Plant Growth Regulation - Excessive utilization of groundwater for anthropogenic purposes has led to severe depletion of the water table, resulting in contamination of fluorides from the...     

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.     

Cross-genera amplification of Cajanus spp. specific SSR markers in Clitoria ternatea (L.) and their application in genetic diversity studies

Clitoria ternatea (L.) is a medicinal leguminous plant and is cultivated to cater the need of herbal industries and asthetic purposes. The unavailability of steady molecular marker impedes the genetic improvement of C. ternatea. In the present study, transferability of 98 pairs of Cajanus spp. specific SSR primers were assessed among 14 genotypes of C. ternatea, varied for their flower color, floral architecture and bio-metabolite (taraxerol and delphinidin) content, and out of them 43 had successfully amplified the fragments. Among them, 36 pairs of primers showed 100% transferability, whereas rest seven varied from 42.86 to 92.85% transferability. The transferable 43 pairs of SSR primers generated 196 alleles across the 14 genotypes and the AMOVA analysis showed moderate genetic variation (55.1%) among the genotypes of C. ternatea, which was also reinforced by Nei’s genetic distance and gene identity estimates derived haplotype matrix. Similarly, both the principal coordinate analysis and dendrogram grouped these 14 genotypes of C. ternatea into two major clusters based on SSR allele distribution and frequency, and the clustering pattern is in accordance with petal color but in contrast to floral architecture. MCheza based outlier analysis revealed 16 alleles for balancing selection, which are putatively involved in the maintenance of genetic polymorphism in C. ternatea. Moreover, the estimates of molecular diversity and bio-metabolite content revealed the possible use of these genotypes in future breeding programme of this species.Electronic supplementary materialThe online version of this article (10.1007/s12298-020-00907-x) contains supplementary material, which is available to authorized users.     

Genetic control of fusion pore expansion in the epidermis of Caenorhabditis elegans

Developmental cell fusion is found in germlines, muscles, bones, placentae, and stem cells. In Caenorhabditis elegans 300 somatic cells fuse during development. Although there is extensive information on the early intermediates of viral-induced and intracellular membrane fusion, little is known about late stages in membrane fusion. To dissect the pathway of cell fusion in C. elegans embryos, we use genetic and kinetic analyses using live-confocal and electron microscopy. We simultaneously monitor the rates of multiple cell fusions in developing embryos and find kinetically distinct stages of initiation and completion of membrane fusion in the epidermis. The stages of cell fusion are differentially blocked or retarded in eff-1 and idf-1 mutants. We generate kinetic cell fusion maps for embryos grown at different temperatures. Different sides of the same cell differ in their fusogenicity: the left and right membrane domains are fusion-incompetent, whereas the anterior and posterior membrane domains fuse with autonomous kinetics in embryos. All but one cell pair can initiate the formation of the largest syncytium. The first cell fusion does not trigger a wave of orderly fusions in either direction. Ultrastructural studies show that epidermal syncytiogenesis require eff-1 activities to initiate and expand membrane merger.     

The karrikin ‘calisthenics’: Can compounds derived from smoke help in stress tolerance?

Karrikins (KARs) are unique butenolides derived as a by‐product of incomplete combustion during wildfire. Some receptive plant species respond to KARs in the form of accelerated germination. These molecules originate from stress to mediate tolerance against different sub‐optimal conditions like oxidative stress, drought and low‐light intensity (shade stress). KARs promote seed germination, seedling establishment and ecological diversity by accelerating the abundance of selective communities of plants. The signaling pathway is closely related, yet unique from strigolactones (SLs). Due to the structural relatedness with SLs, KARs have potential roles in mediating abiotic stress tolerance in plants. In addition, the KAR directly/indirectly interact with crucial phytohormones like abscisic acid, gibberellic acid, auxins and ethylene. This article is a summarized update on KAR research in recent times. The exhaustive discussions would be beneficial for understanding the extraordinary strategy devised by nature to protect plants from stress using a molecule which is itself generated out of stress.     

Targeting Glycinebetaine for Abiotic Stress Tolerance in Crop Plants: Physiological Mechanism,Molecular Interaction and Signaling

In the era of climate change, abiotic stresses (e.g., salinity, drought, extreme temperature, flooding, metal/metalloid(s), UV radiation, ozone, etc.) are considered as one of the most complex environmental constraints that restricts crop production worldwide. Introduction of stress-tolerant crop cultivars is the most auspicious way of surviving this constraint, and to produce these types of tolerant crops. Several bioengineering mechanisms involved in stress signaling are being adopted in this regard. One example of this kind of manipulation is the osmotic adjustment. The quarternary ammonium compound glycinebetaine (GB), also originally referred to as betaine is a methylated glycine derivative. Among the betaines, GB is the most abundant one in plants, which is mostly produced in response to dehydration caused by different abiotic stresses like drought, salinity, and extreme temperature. Glycinebetaine helps in decreased accumulation and detoxification of ROS, thereby restoring photosynthesis and reducing oxidative stress. It takes part in stabilizing membranes and macromolecules. It is also involved in the stabilization and protection of photosynthetic components, such as ribulose-1, 5-bisphosphate carboxylase/oxygenase, photosystem II and quarternary enzyme and protein complex structures under environmental stresses. Glycinebetaine was found to perform in chaperone-induced protein disaggregation. In addition, GB can confer stress tolerance in very low concentrations, and it acts in activating defense responsive genes with stress protection. Recently, field application of GB has also shown protective effects against environmental adversities increasing crop yield and quality. In this review, we will focus on the role of GB in conferring abiotic stress tolerance and the possible ways to engineer GB biosynthesis in plants.     

High‐Performance Hydrogen Evolution by Ru Single Atoms and Nitrided‐Ru Nanoparticles Implanted on N‐Doped Graphitic Sheet

The most efficient electrocatalyst for the hydrogen evolution reaction (HER) is a Pt‐based catalyst, but its high cost and nonperfect efficiency hinder wide‐ranging industrial/technological applications. Here, an electrocatalyst of both ruthenium (Ru) single atoms (SAs) and N‐doped‐graphitic(G_N)‐shell‐covered nitrided‐Ru nanoparticles (NPs) (having a Ru‐N_x shell) embedded on melamine‐derived G_N matrix { 1 : [Ru(SA)+Ru(NP)@RuN_x@G_N]/G_N}, which exhibits superior HER activity in both acidic and basic media, is presented. In 0.5 m H₂SO₄/1 m KOH solutions, 1 shows diminutive “negative overpotentials” (?η = |η| = 10/7 mV at 10 mA cm^?2, lowest ever) and high exchange current densities (4.70/1.96 mA cm^?2). The remarkable HER performance is attributed to the near‐zero free energies for hydrogen adsorption/desorption on Ru(SAs) and the increased conductivity of melamine‐derived G_N sheets by the presence of nitrided‐Ru(NPs). The nitridation process forming nitrided‐Ru(NPs), which are imperfectly covered by a G_N shell, allows superb long‐term operation durability. The catalyst splits water into molecular oxygen and hydrogen at 1.50/1.40 V (in 0.1 m HClO₄/1 m KOH), demonstrating its potential as a ready‐to‐use, highly effective energy device for industrial applications.