Potential Protein-coding Genes Research Articles

Telomere-to-telomere genome assembly of Oldenlandia diffusa.

We report the complete telomere-to-telomere genome assembly of Oldenlandia diffusa which renowned in traditional Chinese medicine, comprising 16 chromosomes and spanning 499.7Mb. The assembly showcases 28 telomeres and minimal gaps, with a total of only five. Repeat sequences constitute 46.41% of the genome, and 49,701 potential protein-coding genes have been predicted. Compared with O. corymbosa, O. diffusa exhibits chromosome duplication and fusion events, diverging 20.34 million years ago. Additionally, a total of 11 clusters of terpene synthase have been identified. The comprehensive genome sequence, gene catalog, and terpene synthase clusters of O. diffusa detailed in this study will significantly contribute to advancing research in this species' genetic, genomic, and pharmacological aspects.

Telomere-to-telomere genome assembly of matsutake (Tricholoma matsutake)

Here, we report the first telomere-to-telomere genome assembly of matsutake (Tricholoma matsutake), which consists of 13 sequences (spanning 161.0Mb) and a 76kb circular mitochondrial genome. All the 13 sequences were supported with telomeric repeats at the ends. GC-rich regions are located at the middle of the sequences and are enriched with long interspersed nuclear elements (LINEs). Repetitive sequences including long-terminal repeats (LTRs) and LINEs occupy 71.6% of the genome. A total of 21,887 potential protein-coding genes were predicted. The genomic data reported in this study served not only matsutake gene sequences but also genome structures and intergenic sequences. The information gained would be a great reference for exploring the genetics, genomics, and evolutionary study of matsutake in the future, and ultimately facilitate the conservation of this vulnerable genetic resource.

Decontamination and Annotation of the Draft Genome Sequence of the Oomycete Lagenidium giganteum ARSEF 373

Scaffolds of a previously published Lagenidium giganteum ARSEF 373 genome assembly found at GenBank were filtered to remove contaminating sequences. Genome annotation of the 437 scaffolds (total length, 56.2 MB; GC content, 58.8%) with confirmed L. giganteum sequences identified 13,069 potential protein-coding genes, encoding at least 737 predicted secreted proteins and >100 putative translocated effectors.

Nuclear Genome Sequence and Gene Expression of an Intracellular Fungal Endophyte Stimulating the Growth of Cranberry Plants.

Ericaceae thrive in poor soil, which we postulate is facilitated by microbes living inside those plants. Here, we investigate the growth stimulation of the American cranberry (Vaccinium macrocarpon) by one of its fungal endosymbionts, EC4. We show that the symbiont resides inside the epidermal root cells of the host but extends into the rhizosphere via its hyphae. Morphological classification of this fungus is ambiguous, but phylogenetic inference based on 28S rRNA identifies EC4 as a Codinaeella species (Chaetosphaeriaceae, Sordariomycetes, Ascomycetes). We sequenced the genome and transcriptome of EC4, providing the first 'Omics' information of a Chaetosphaeriaceae fungus. The 55.3-Mbp nuclear genome contains 17,582 potential protein-coding genes, of which nearly 500 have the capacity to promote plant growth. For comparing gene sets involved in biofertilization, we annotated the published genome assembly of the plant-growth-promoting Trichoderma hamatum. The number of proteins involved in phosphate transport and solubilization is similar in the two fungi. In contrast, EC4 has ~50% more genes associated with ammonium, nitrate/nitrite transport, and phytohormone synthesis. The expression of 36 presumed plant-growth-promoting EC4 genes is stimulated when the fungus is in contact with the plant. Thus, Omics and in-plantae tests make EC4 a promising candidate for cranberry biofertilization on nutrient-poor soils.

Characterization of the Complete Chloroplast Genome and Phylogenetic Implications of Euonymus microcarpus (Oliv.) Sprague

Euonymus microcarpus (Oliv.) Sprague, is a species of evergreen shrub of the genus Euonymus, family Celastraceae. Here, we extracted the genomic DNA from the leaves of E. microcarpus and constructed a paired-end library. The chloroplast genome of E. microcarpus was generated with the high-throughput sequencing by the illumina Hiseq X Ten platform and de novo assembly. The chloroplast genome had a quadripartite structure, containing a long single copy region with a size of 85,386 bp and a short single copy region with a size of 18,456 bp, separated by two inverted repeat regions of 26,850 bp. The chloroplast genome contained 133 genes identified in total, including 87 potential protein-coding genes, 38 transfer RNA genes, and eight ribosomal RNA genes. A total of 282 simple sequence repeats and 63 long repeats were found. Furthermore, the phylogenetic relationships inferred that E. microcarpus is sister to E. japonicus and E. schensianus. A comparison of the structure of the chloroplast genomes of eight Euonymus species suggests a nucleotide variability of the junction sites and a higher divergence of non-coding regions, compared to the coding regions. The original findings of the study serves as a good reference for chloroplast genome assembly and a valuable foundation for the genetic diversity and evolution of E. microcarpus.

Anatilimnocola floriformis sp. nov., a novel member of the family Pirellulaceae from a boreal lake, and emended description of the genus Anatilimnocola

Planctomycetes of the family Pirellulaceae are commonly addressed as budding aquatic bacteria with a complex lifestyle. Although this family is well represented by cultured and taxonomically characterized isolates, nearly all of them were obtained from brackish or marine habitats. The examples of described freshwater Pirellulaceae planctomycetes are limited to two species only, Pirellula staley and 'Anatilimnocola aggregata'. In this study, we characterized a novel freshwater planctomycete of the genus 'Anatilimnocola', strain PX40T, which was isolated from a boreal eutrophic lake. Strain PX40T was represented by budding, unpigmented, ellipsoidal to pear-shaped cells, which often occurred in characteristic flower-like rosettes. Cells were covered by bundles of fimbriae; crateriform-like structures were localized on a reproductive cell pole only. These planctomycetes were obligately aerobic, heterotrophic bacteria that utilized various sugars and some polysaccharides, and were highly sensitive to NaCl. Growth occurred in the pH range 5.0-7.5 (with an optimum at pH 6.5-7.0), and at temperatures between 15 and 30°C (with an optimum at 22-25°C). The major fatty acids of strain PX40T were C18:1ω9c, C16:0, and 16:1ω7c; cells also contained a wide variety of hydroxy- and dihydroxy-fatty acids and a C31:9 alkene. The major intact polar lipids were diacylglyceryl-(N,N,N)-trimethylhomoserines. The 16S rRNA gene sequence of strain PX40T displayed 96.6% similarity to that of 'Anatilimnocola aggregata' ETA_A8T. The genome of strain PX40T was 8.93Mb in size and contained one copy of rRNA operon, 76 tRNA genes and 7092 potential protein-coding genes. The DNA G+C content was 57.8%. The ANI value between strain PX40T and 'Anatilimnocola aggregata' ETA_A8T was 78.3%, suggesting that these planctomycetes represent distinct species. We, therefore, propose a novel species of the genus 'Anatilimnocola', 'A. floriformis' sp. nov., with strain PX40T (= KCTC 92369T = VKM B-3621T = UQM 41463T) as the type strain.

Construction of a Novel circRNA/miRNA/mRNA Regulatory Network to Explore the Potential Pathogenesis of Wilson's Disease.

Studies show that non-coding RNAs, especially microRNAs (miRNAs) and circular RNAs (circRNAs), and protein-coding genes are involved in the pathophysiology of multi-organ damage caused by Wilson’s disease (WD). However, circRNA expression profiles and their role in initiation and progression of WD kidney injury remain largely unclear at present. Here, we explored potential critical protein-coding genes, miRNAs, and circRNAs, as well as identify competitive endogenous RNAs (ceRNAs) in a WD mouse model by high-throughput sequencing. We investigated the expression profiles of circRNAs, miRNAs, and protein-coding genes, and identified 32 DEcircRs, 45 DEmiRs, and 1623 DEPs. Identified DEcircRs, DEmiRs, and DEPs were used to construct a ceRNA network, which consisted of 15 DEcircRNAs (four upregulated and 11 downregulated), 18 DEmiRNAs (14 upregulated and four downregulated), and 352 DEmRNAs (205 upregulated and 147 downregulated). Further experiments proved that mmu_circ_0001333 and mmu_circ_0000355 acted as sponges of miR-92b-5p, miR-107-3p, and miR-187-3p to regulate the expression of genes including Smad9, Mapk10, and Aldh3a2, which may participate in WD-related kidney injury. Taken together, this study identified the circRNA/miRNA/mRNA network involved in kidney failure in WD, which may serve as a potential biomarker for the pathogenesis of WD.

Populus trichocarpa (Black Cottonwood) As a Model for Plant Genomic Studies: A Review

Numerous plants have been the subject of recent research in the pharmacological, cosmetic, and agro-alimentary domains due to their chemical composition and multiple therapeutic capabilities. Populus trichocarpa is one of the most common trees found in deciduous forests (Salicaceae family). The current study examines Populus trichocarpa as a model plant for plant genomics research, as well as the most recent findings on phytochemical composition and medicinal potential. More than 45,000 potential protein-coding genes were discovered. In the Populus genome, a whole-genome duplication event was discovered, with approximately 8,000 pairs of duplicated genes surviving. Furthermore, the reproductive biology of Populus provides new opportunities and challenges in the study and analysis of natural genetic and phenotypic variation. In the present review, we endeavour to describe and compile the available knowledge on Populus trichocarpa as a model plant for genomic investigations and to bring that material up to date of Populus trichocarpa's phytochemical and medicinal properties.

Complete genome sequence of the cellulolytic planctomycete Telmatocola sphagniphila SP2T and characterization of the first cellulolytic enzyme from planctomycetes

Planctomycetes of the family Gemmataceae are strictly aerobic chemo-organotrophs that display a number of hydrolytic capabilities. A member of this family, Telmatocola sphagniphila SP2T, is the first described planctomycete with experimentally proven ability for growth on cellulose. In this study, the complete genome sequence of strain SP2T was obtained and the genome-encoded determinants of its cellulolytic potential were analyzed. The T. sphagniphila SP2T genome was 6.59 Mb in size and contained over 5200 potential protein-coding genes. The search for enzymes that could be potentially involved in cellulose degradation identified a putative cellulase that contained a domain from the GH44 family of glycoside hydrolases. Homologous enzymes were also revealed in the genomes of two other Gemmataceae planctomycetes, Zavarzinella formosa A10T and Tuwongella immobilis MBLW1T. The gene encoding this predicted cellulase in strain SP2T was expressed in E. coli and the hydrolytic activity of the recombinant enzyme was confirmed in tests with carboxymethyl cellulose but not with crystalline cellulose, xylan, mannan or laminarin. This is the first experimentally characterized cellulolytic enzyme from planctomycetes.

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, 'Kawazu-zakura' and 'Atami-zakura'.

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.

Genomic region associated with pod color variation in pea (Pisum sativum).

Pea (Pisum sativum) was chosen as the research material by Gregor Mendel to discover the laws of inheritance. Out of seven traits studied by Mendel, genes controlling three traits including pod shape, pod color, and flower position have not been identified to date. With the aim of identifying the genomic region controlling pod color, we determined the genome sequence of a pea line with yellow pods. Genome sequence reads obtained using a Nanopore sequencing technology were assembled into 117,981 contigs (3.3 Gb), with an N50 value of 51.2 kb. A total of 531,242 potential protein-coding genes were predicted, of which 519,349 (2.8 Gb) were located within repetitive sequences (2.8 Gb). The assembled sequences were ordered using a reference as a guide to build pseudomolecules. Subsequent genetic and association analyses led to the identification of a genomic region that controls pea pod color. DNA sequences at this genomic location and transcriptome profiles of green and yellow pod lines were analyzed, and genes encoding 3' exoribonucleases were selected as potential candidates controlling pod color. The results presented in this study are expected to accelerate pan-genome studies in pea and facilitate the identification of the gene controlling one of the traits studied by Mendel.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility.

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1Mb of the S2 -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S6a -haplotype of P. inflata and the SN -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

Integrative analysis provides multi-omics evidence for the pathogenesis of placenta percreta.

Pernicious placenta previa with placenta percreta (PP) is a catastrophic condition during pregnancy. However, the underlying pathogenesis remains unclear. In the present study, the placental tissues of normal cases and PP tissues of pernicious placenta previa cases were collected to determine the expression profile of protein‐coding genes, miRNAs, and lncRNAs through sequencing. Weighted gene co‐expression network analysis (WGCNA), accompanied by miRNA target prediction and correlation analysis, were employed to select potential hub protein‐coding genes and lncRNAs. The expression levels of selected protein‐coding genes, Wnt5A and MAPK13, were determined by quantitative PCR and immunohistochemical staining, and lncRNA PTCHD1‐AS and PAPPA‐AS1 expression levels were determined by quantitative PCR and fluorescence in situ hybridization. The results indicated that 790 protein‐coding genes, 382 miRNAs, and 541 lncRNAs were dysregulated in PP tissues, compared with normal tissues. WGCNA identified coding genes in the module (ME) black and ME turquoise modules that may be involved in the pathogenesis of PP. The selected potential hub protein‐coding genes, Wnt5A and MAPK13, were down‐regulated in PP tissues, and their expression levels were positively correlated with the expression levels of PTCHD1‐AS and PAPPA‐AS1. Further analysis demonstrated that PTCHD1‐AS and PAPPA‐AS1 regulated Wnt5A and MAPK13 expression by interacting with specific miRNAs. Collectively, our results provided multi‐omics data to better understand the pathogenesis of PP and help identify predictive biomarkers and therapeutic targets for PP.

Ongoing Coevolution of Wolbachia and a Widespread Invasive Ant, Anoplolepis gracilipes.

While Wolbachia are commonly found among arthropods, intraspecific infection rates can vary substantially across the geographic populations. Here we report nearly 100% prevalence of Wolbachia in the global populations of the yellow crazy ant, Anoplolepis gracilipes. To understand coevolutionary history between Wolbachia and A. gracilipes, we identified single nucleotide polymorphisms (SNPs) in Wolbachia from the ant across 12 geographical regions and compared the phylogeny of SNP-based Wolbachia to patterns of the ant’s mitochondrial DNA (mtDNA) variation. Our results revealed a strong concordance between phylogenies of Wolbachia and host mtDNA, providing immediate evidence of co-divergence. Among eight identified SNP loci separating the genetic clusters of Wolbachia, seven loci are located in potential protein-coding genes, three of which being non-synonymous SNPs that may influence gene functions. We found a Wolbachia hypothetical protein gene with signature of positive selection. These findings jointly allow us to characterize Wolbachia-ant coevolution and also raise a question about mechanism(s) underlying maintenance of high prevalence of Wolbachia during the colonization of this invasive ant.

Frigoriglobus tundricola gen. nov., sp. nov., a psychrotolerant cellulolytic planctomycete of the family Gemmataceae from a littoral tundra wetland

Planctomycetes of the family Gemmataceae are characterized by large genome sizes and cosmopolitan distribution in freshwater and terrestrial environments but their ecological functions remain poorly understood. In this study, we characterized a novel representative of this family, strain PL17T, which was isolated from a littoral tundra wetland and was capable of growth on xylan and cellulose. Cells of this isolate were represented by pink-pigmented spheres that multiplied by budding and occurred singly or in short chains and aggregates. Strain PL17T was obligately aerobic, mildly acidophilic chemoorganotrophic bacterium, which displayed good tolerance of low temperatures. The major fatty acids were C18:0, C16:1ω5, and βOH-C16:1; the major polar lipid was trimethylornithine. The genome of strain PL17T consisted of a 9.83 Mb chromosome and a 24.69kb plasmid. The G+C contents of the chromosomal and plasmid DNA were 67.4 and 62.3mol%, respectively. Over 8900 potential protein-coding genes were identified in the genome including a putative cellulase that contains a domain from the GH5 family of glycoside hydrolases. The genome of strain PL17T contained one linked and one unlinked rRNA operons with 16S rRNA gene sequences displaying 94.5% similarity to that in Gemmata obscuriglobus UQM2246T. Based on the results of comparative phenotypic, chemotaxonomic and phylogenomic analyses, we propose to classify strain PL17T (= CECT 9407T=VKM B-3467T) as representing a novel genus and species of the family Gemmataceae, Frigoriglobus tundricola gen. nov., sp. nov.

Thriving in Wetlands: Ecophysiology of the Spiral-Shaped Methanotroph Methylospira mobilis as Revealed by the Complete Genome Sequence.

Candidatus Methylospira mobilis is a recently described spiral-shaped, micro-aerobic methanotroph, which inhabits northern freshwater wetlands and sediments. Due to difficulties of cultivation, it could not be obtained in a pure culture for a long time. Here, we report on the successful isolation of strain Shm1, the first axenic culture of this unique methanotroph. The complete genome sequence obtained for strain Shm1 was 4.7 Mb in size and contained over 4800 potential protein-coding genes. The array of genes encoding C1 metabolic capabilities in strain Shm1 was highly similar to that in the closely related non-motile, moderately thermophilic methanotroph Methylococcus capsulatus Bath. The genomes of both methanotrophs encoded both low- and high-affinity oxidases, which allow their survival in a wide range of oxygen concentrations. The repertoire of signal transduction systems encoded in the genome of strain Shm1, however, by far exceeded that in Methylococcus capsulatus Bath but was comparable to those in other motile gammaproteobacterial methanotrophs. The complete set of motility genes, the presence of both the molybdenum–iron and vanadium-iron nitrogenases, as well as a large number of insertion sequences were also among the features, which define environmental adaptation of Methylospira mobilis to water-saturated, micro-oxic, heterogeneous habitats depleted in available nitrogen.

Limnoglobus roseus gen. nov., sp. nov., a novel freshwater planctomycete with a giant genome from the family Gemmataceae.

The family Gemmataceae accommodates aerobic, chemoorganotrophic planctomycetes, which inhabit various freshwater ecosystems, wetlands and soils. Here, we describe a novel member of this family, strain PX52T, which was isolated from a boreal eutrophic lake in Northern Russia. This isolate formed pink-pigmented colonies and was represented by spherical cells that occurred singly, in pairs or aggregates and multiplied by budding. Daughter cells were highly motile. PX52T was an obligate aerobic chemoorganotroph, which utilized various sugars and some heteropolysaccharides. Growth occurred at pH 5.0–7.5 (optimum pH 6.5) and at temperatures between 10 and 30 °C (optimum 20–25 °C). The major fatty acids were C18 : 1ɷ7c, C18 : 0 and βOH-C16:0; the major intact polar lipid was trimethylornithine, and the quinone was MK-6. The complete genome of PX52T was 9.38 Mb in size and contained nearly 8000 potential protein-coding genes. Among those were genes encoding a wide repertoire of carbohydrate-active enzymes (CAZymes) including 33 glycoside hydrolases (GH) and 87 glycosyltransferases (GT) affiliated with 17 and 12 CAZy families, respectively. DNA G+C content was 65.6 mol%. PX52T displayed only 86.0–89.8 % 16S rRNA gene sequence similarity to taxonomically described Gemmataceae planctomycetes and differed from them by a number of phenotypic characteristics and by fatty acid composition. We, therefore, propose to classify it as representing a novel genus and species, Limnoglobus roseus gen. nov., sp. nov. The type strain is strain PX52T (=KCTC 72397T=VKM B-3275T).

Hybrid de novo whole-genome assembly and annotation of the model tapeworm Hymenolepis diminuta

Despite the use of Hymenolepis diminuta as a model organism in experimental parasitology, a full genome description has not yet been published. Here we present a hybrid de novo genome assembly based on complementary sequencing technologies and methods. The combination of Illumina paired-end, Illumina mate-pair and Oxford Nanopore Technology reads greatly improved the assembly of the H. diminuta genome. Our results indicate that the hybrid sequencing approach is the method of choice for obtaining high-quality data. The final genome assembly is 177 Mbp with contig N50 size of 75 kbp and a scaffold N50 size of 2.3 Mbp. We obtained one of the most complete cestode genome assemblies and annotated 15,169 potential protein-coding genes. The obtained data may help explain cestode gene function and better clarify the evolution of its gene families, and thus the adaptive features evolved during millennia of co-evolution with their hosts.

Dysregulated Mirnas after Uniform Treatment Predict Outcome of Newly-Diagnosed Multiple Myeloma

Multiple myeloma (MM) is a malignancy of the plasma cell in which clonal plasma cells infiltrate the bone marrow. Although increasingly effective treatment has significantly improved management of the disease, most MM patients ultimately relapse. Previous studies suggest that miRNA expression patterns may function as biomarkers for diagnosis, subtyping, and prognosis of the disease. However, the role of miRNAs in MM relapse and prognosis is not sufficiently understood. Therefore, in this study, we sought to understand the biological and prognostic role of dysregulated miRNAs between newly diagnosed and relapsed MM. We performed RNA-seq as well as microarray-based miRNA profiling on 319 MM samples from newly diagnosed MM patients enrolled in the 2009 IFM/DFCI clinical trial. This included 40 samples from patients with paired samples from diagnosis and at relapse while we had additional 239 samples collected at the time of diagnosis. Using the 40 diagnosis-relapse paired samples, we identified 48 miRNAs that are differentially expressed between diagnosis and relapse (FDR < 0.15). Of these 48 miRNAs, we found five miRNAs (miR-4484, miR-663b, miR-3687, miR-670-5p, miR-4417) whose expression level at diagnosis, expression level at relapse, and log2FC from diagnosis to relapse were significantly associated with overall survival or progression-free survival among the 40 patients with paired samples (p < 0.05). We next investigated the prognostic role of these five miRNAs through unsupervised hierarchical clustering of 239 independent, newly diagnosed MM samples. The two groups identified by the unsupervised hierarchical clustering differed significantly in both overall survival (log rank test, p = 0.009) and progression-free survival (log rank test, p = 0.012), suggesting that these five miRNAs are predictive of disease prognosis. We then performed an integrative analysis of the miRNA and mRNA expression data in our cohort of newly diagnosed samples and identified the potential protein-coding target genes of these five miRNAs. We further filtered the target gene list using TargetScan which predicts the biological targets of miRNAs by searching for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA. After filtering the target genes, we found that each miRNA targets 126 mRNAs [range 26 - 258] and key MM genes, including MAF, BCL2, 14-3-3 and MYEOV, were among the top candidates. The target genes of these select miRNAs are involved in such biological pathways as NF-kB signaling, antigen processing and presentation of endogenous peptide antigen via MHC class I via ER pathway, and positive regulation of cell death and apoptosis. In conclusion, this integrative analysis of miRNA and mRNA expression data from uniformly treated MM samples offers a unique perspective on the role of miRNAs in MM biology and prognosis. The study identified five miRNAs that target key MM genes, play a role in the biological processes associated with MM progression, and carry prognostic information predictive of overall and progression-free survival with prospect for future therapeutic application. Disclosures Richardson: Amgen: Membership on an entity's Board of Directors or advisory committees; Janssen: Membership on an entity's Board of Directors or advisory committees; Bristol-Myers Squibb: Research Funding; Takeda: Membership on an entity's Board of Directors or advisory committees, Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Membership on an entity's Board of Directors or advisory committees, Research Funding; Karyopharm: Membership on an entity's Board of Directors or advisory committees; Sanofi: Membership on an entity's Board of Directors or advisory committees. Moreau:AbbVie: Consultancy, Honoraria; Amgen: Consultancy, Honoraria; Takeda: Consultancy, Honoraria; Celgene: Consultancy, Honoraria; Janssen: Consultancy, Honoraria. Anderson:Sanofi-Aventis: Other: Advisory Board; Bristol-Myers Squibb: Other: Scientific Founder; Oncopep: Other: Scientific Founder; Amgen: Consultancy, Speakers Bureau; Janssen: Consultancy, Speakers Bureau; Takeda: Consultancy, Speakers Bureau; Celgene: Consultancy, Speakers Bureau. Munshi:Celgene: Consultancy; Amgen: Consultancy; Adaptive: Consultancy; Abbvie: Consultancy; Janssen: Consultancy; Takeda: Consultancy; Oncopep: Consultancy.

ChemGenome2.1: An Ab Initio Gene Prediction Software.

Gene prediction, also known as gene identification, gene finding, gene recognition, or gene discovery, is among one of the important problems of molecular biology and is receiving increasing attention due to the advent of large-scale genome sequencing projects. We designed an ab initio model (called ChemGenome) for gene prediction in prokaryotic genomes based on physicochemical characteristics of codons. In this chapter, we present the methodology of the latest version of this model ChemGenome2.1 (CG2.1). The first module of the protocol builds a three-dimensional vector from three calculated quantities for each codon-the double-helical trinucleotide base pairing energy, the base pair stacking energy, and an index of the propensity of a codon for protein-nucleic acid interactions. As this three-dimensional vector moves along any genome, the net orientation of the resultant vector should differ significantly for gene and non-genic regions to make a distinction feasible. The predicted putative protein-coding genes from above parameters are passed through a second module of the protocol which reduces the number of false positives by utilizing a filter based on stereochemical properties of protein sequences. The chemical properties of amino acid side chains taken into consideration are the presence of sp3 hybridized γ carbon atom, hydrogen bond donor ability, short/absence of δ carbon and linearity of the side chains/non-occurrence of bi-dentate forks with terminal hydrogen atoms in the side chain. The final prediction of the potential protein-coding genes is based on the frequency of occurrence of amino acids in the predicted protein sequences and their deviation from the frequency values of Swissprot protein sequences, both at monomer and tripeptide levels. The final screening is based on Z-score. Though CG2.1 is a gene finding tool for prokaryotes, considering the underlying similarity in the chemical and physical properties of DNA among prokaryotes and eukaryotes, we attempted to evaluate its applicability for gene finding in the lower eukaryotes. The results give a hope that the concept of gene finding based on physicochemical model of codons is a viable idea for eukaryotes as well, though, undoubtedly, improvements are needed.