Deep Amplicon Sequencing Reveals Culture Selection of Mycobacterium Tuberculosis by Clinical Samples.

The commonly-used drug susceptibility testing (DST) relies on bacterial culture and faces shortcomings such as long turnaround time and clone/subclone selection. We developed a targeted deep amplification sequencing (DAS) method directly applied to clinical specimens. In this DAS panel, we examined 941 drug-resistant mutations associated with 20 anti-tuberculosis drugs with an initial amount of 4 pg DNA and reduced clinical testing time from 20 days to two days. A prospective study was conducted using 115 clinical specimens mainly with Xpert® Mycobacterium tuberculosis/rifampicin (Xpert MTB/RIF) assay positive to evaluate drug-resistant mutation detection. DAS was performed on culture-free specimens, while culture-dependent isolates were used for phenotypic DST, DAS, and whole-genome sequencing (WGS). For in silico molecular DST, our result based on DAS panel revealed the similar accuracy to three published reports based on WGS. For 82 isolates, application of DAS showed better sensitivity (93.03% vs. 92.16%), specificity (96.10% vs. 95.02%), and accuracy (91.33% vs. 90.62%) than Mykrobe software using WGS. Compared to culture-dependent WGS, culture-free DAS provides a full picture of sequence variation at population level, exhibiting in detail the gain-and-loss variants caused by bacterial culture. Our study performs a systematic verification of the advantages of DAS in clinical applications and comprehensively illustrates the discrepancy in Mycobacterium tuberculosis before and after culture.

A tool to automatically design multiplex PCR primer pairs for specific targets using diverse templates.

Multiplex PCR is an increasingly popular method for identifying species, investigating environmental diversity, and conducting phylogenetic analysis. The complexity and increasing availability of diverse templates necessitate a highly automated approach to design degenerate primer pairs for specific targets with multiple sequences. Existing tools for degenerate primer design suffer from poor maintenance, semi-automation, low adaptability, and low tolerance for gaps. We developed PMPrimer, a Python-based tool for automated design and evaluation of multiplex PCR primer pairs for specific targets using diverse templates. PMPrimer automatically designs optimal multiplex PCR primer pairs using a statistical-based template filter; performs multiple sequence alignment, conserved region identification, and primer design; and evaluates the primers based on template coverage, taxon specificity, and target specificity. PMPrimer identifies conserved regions using Shannon's entropy method, tolerates gaps using a haplotype-based method, and evaluates multiplex PCR primer pairs based on template coverage and taxon specificity. We tested PMPrimer using datasets with diverse levels of conservation, sizes, and applications, including tuf genes of Staphylococci, hsp65 genes of Mycobacteriaceae, and 16S ribosomal RNA genes of Archaea. PMPrimer showed outstanding performance compared with existing tools and experimental validated primers. PMPrimer is available as a Python package at https://github.com/AGIScuipeng/PMPrimer .

The wild allotetraploid sesame genome provides novel insights into evolution and lignan biosynthesis.

INTRODUCTION: The wild tetraploid sesame (Sesamum schinzianum), an ancestral relative of diploid cultivated sesame, grows in the tropical desert of the African Plateau. As a valuable seed resource, wild sesame has several advantageous traits, such as strong environmental adaptability and an extremely high content of sesamolin in its seeds. High-quality genome assembly is essential for a detailed understanding of genome structure, genome evolution and crop improvement. OBJECTIVES: Here, we generated two high-quality chromosome-scale genomes from S. schinzianum and a cultivated diploid elite sesame (Sesamum indicum L.) to investigate the potential genetic basis underlying these traits of wild sesame. METHODS: The long-read data from PacBio Sequel II platform and high-throughput chromosome conformation capture (Hi-C) data were used to construct high-quality sesame genome. Then dissecting the molecular mechanisms of sesame evolution and lignan biosynthesis through comparative genomics and transcriptomics. RESULTS: We found evidence of divergent evolution that involved differences in the number, sequence and expression level of homologous genes between the two sets of subgenomes from allotetraploids in S. schinzianum, all of which might be driven by subfunctionalization after polyploidization. Furthermore, it was found that a great number of genes involved in the stress response have undergone positive selection and resulted from gene family expansion in the wild sesame genome compared with the cultivated sesame genome, which, overall, is associated with adaptative evolution to the environment. We hypothesized that the sole functional member CYP92B14 (SscC22g35272) could be associated with high content of sesamolin in wild sesame seeds. CONCLUSION: This study provides high-quality wild allotetraploid sesame and cultivated sesame genomes, reveals evolutionary features of the allotetraploid genome and provides novel insights into lignan synthesis pathways. Meanwhile, the wild sesame genome will be an important resource to conduct comparative genomic and evolutionary studies and plant improvement programmes.

A Chromosome-level Genome Assembly of Wild Castor Provides New Insights into its Adaptive Evolution in Tropical Desert.

Wild castor grows in the high-altitude tropical desert of the African Plateau, a region known for high ultraviolet radiation, strong light, and extremely dry condition. To investigate the potential genetic basis of adaptation to both highland and tropical deserts, we generated a chromosome-level genome sequence assembly of the wild castor accession WT05, with a genome size of 316 Mb, a scaffold N50 of 31.93 Mb, and a contig N50 of 8.96 Mb, respectively. Compared with cultivated castor and other Euphorbiaceae species, the wild castor exhibits positive selection and gene family expansion for genes involved in DNA repair, photosynthesis, and abiotic stress responses. Genetic variations associated with positive selection were identified in several key genes, such as LIG1, DDB2, and RECG1, involved in nucleotide excision repair. Moreover, a study of genomic diversity among wild and cultivated accessions revealed genomic regions containing selection signatures associated with the adaptation to extreme environments. The identification of the genes and alleles with selection signatures provides insights into the genetic mechanisms underlying the adaptation of wild castor to the high-altitude tropical desert and would facilitate direct improvement of modern castor varieties.

Improved 93-11 Genome and Time-Course Transcriptome Expand Resources for Rice Genomics.

In 2002, the first crop genome was published using the rice cultivar 93-11, which is the progenitor of the first super-hybrid rice. The genome sequence has served as a reference genome for the indica cultivars, but the assembly has not been updated. In this study, we update the 93-11 genome assembly to a gap-less sequence using ultra-depth single molecule real-time (SMRT) reads, Hi-C sequencing, reference-guided, and gap-closing approach. The differences in the genome collinearity and gene content between the 93-11 and the Nipponbare reference genomes confirmed to map the indica cultivar sequencing data to the 93-11 genome, instead of the reference. Furthermore, time-course transcriptome data showed that the expression pattern was consistently correlated with the stages of seed development. Alternative splicing of starch synthesis-related genes and genomic variations of waxy make it a novel resource for targeted breeding. Collectively, the updated high quality 93-11 genome assembly can improve the understanding of the genome structures and functions of Oryza groups in molecular breeding programs.

Sequencing of Chinese castor lines reveals genetic signatures of selection and yield-associated loci.

Oil produced by castor (Ricinus communis) has broad industrial applications. However, knowledge on the genetic diversity, especially genetic alterations that occurred during domestication and subsequent traits selection, of this oil crop is limited. Here, our population genomics analyses show that the Chinese castors have developed a geographic pattern, classified into the southern-, the middle-, and the northern-China groups. We detect a number of candidate genomic loci that are associated with the selection signals during the geographical differentiation and domestication. Using genome-wide association analysis, we identify candidate genes associated with nine agronomically important traits. One of the candidate genes encoding a glycosyltransferase related to cellulose and lignin biosynthesis is associated with both capsule dehiscence and endocarp thickness. We hypothesize that the abundance of cellulose or lignin in endocarp is an important factor for capsule dehiscence. Our results provide foundation for castor breeding and genetic study.

RGAAT: A Reference-based Genome Assembly and Annotation Tool for New Genomes and Upgrade of Known Genomes.

The rapid development of high-throughput sequencing technologies has led to a dramatic decrease in the money and time required for de novo genome sequencing or genome resequencing projects, with new genome sequences constantly released every week. Among such projects, the plethora of updated genome assemblies induces the requirement of version-dependent annotation files and other compatible public dataset for downstream analysis. To handle these tasks in an efficient manner, we developed the reference-based genome assembly and annotation tool (RGAAT), a flexible toolkit for resequencing-based consensus building and annotation update. RGAAT can detect sequence variants with comparable precision, specificity, and sensitivity to GATK and with higher precision and specificity than Freebayes and SAMtools on four DNA-seq datasets tested in this study. RGAAT can also identify sequence variants based on cross-cultivar or cross-version genomic alignments. Unlike GATK and SAMtools/BCFtools, RGAAT builds the consensus sequence by taking into account the true allele frequency. Finally, RGAAT generates a coordinate conversion file between the reference and query genomes using sequence variants and supports annotation file transfer. Compared to the rapid annotation transfer tool (RATT), RGAAT displays better performance characteristics for annotation transfer between different genome assemblies, strains, and species. In addition, RGAAT can be used for genome modification, genome comparison, and coordinate conversion. RGAAT is available at https://sourceforge.net/projects/rgaat/ and https://github.com/wushyer/RGAAT_v2 at no cost.

REDO: RNA Editing Detection in Plant Organelles Based on Variant Calling Results.

RNA editing is a post-transcriptional or cotranscriptional process that changes the sequence of the precursor transcript by substitutions, insertions, or deletions. Almost all of the land plants undergo RNA editing in organelles (plastids and mitochondria). Although several software tools have been developed to identify RNA editing events, there has been a great challenge to distinguish true RNA editing events from genome variation, sequencing errors, and other factors. Here we introduce REDO, a comprehensive application tool for identifying RNA editing events in plant organelles based on variant call format files from RNA-sequencing data. REDO is a suite of Perl scripts that illustrate a bunch of attributes of RNA editing events in figures and tables. REDO can also detect RNA editing events in multiple samples simultaneously and identify the significant differential proportion of RNA editing loci. Comparing with similar tools, such as REDItools, REDO runs faster with higher accuracy, and more specificity at the cost of slightly lower sensitivity. Moreover, REDO annotates each RNA editing site in RNAs, whereas REDItools reports only possible RNA editing sites in genome, which need additional steps to obtain RNA editing profiles for RNAs. Overall, REDO can identify potential RNA editing sites easily and provide several functions such as detailed annotations, statistics, figures, and significantly differential proportion of RNA editing sites among different samples.

DRDB: An Online Date Palm Genomic Resource Database.

Background: Date palm (Phoenix dactylifera L.) is a cultivated woody plant with agricultural and economic importance in many countries around the world. With the advantages of next generation sequencing technologies, genome sequences for many date palm cultivars have been released recently. Short sequence repeat (SSR) and single nucleotide polymorphism (SNP) can be identified from these genomic data, and have been proven to be very useful biomarkers in plant genome analysis and breeding. Results: Here, we first improved the date palm genome assembly using 130X of HiSeq data generated in our lab. Then 246,445 SSRs (214,901 SSRs and 31,544 compound SSRs) were annotated in this genome assembly; among the SSRs, mononucleotide SSRs (58.92%) were the most abundant, followed by di- (29.92%), tri- (8.14%), tetra- (2.47%), penta- (0.36%), and hexa-nucleotide SSRs (0.19%). The high-quality PCR primer pairs were designed for most (174,497; 70.81% out of total) SSRs. We also annotated 6,375,806 SNPs with raw read depth≥3 in 90% cultivars. To further reduce false positive SNPs, we only kept 5,572,650 (87.40% out of total) SNPs with at least 20% cultivars support for downstream analyses. The high-quality PCR primer pairs were also obtained for 4,177,778 (65.53%) SNPs. We reconstructed the phylogenetic relationships among the 62 cultivars using these variants and found that they can be divided into three clusters, namely North Africa, Egypt - Sudan, and Middle East - South Asian, with Egypt - Sudan being the admixture of North Africa and Middle East - South Asian cultivars; we further confirmed these clusters using principal component analysis. Moreover, 34,346 SSRs and 4,177,778 SNPs with PCR primers were assigned to shared cultivars for cultivar classification and diversity analysis. All these SSRs, SNPs and their classification are available in our database, and can be used for cultivar identification, comparison, and molecular breeding. Conclusion:DRDB is a comprehensive genomic resource database of date palm. It can serve as a bioinformatics platform for date palm genomics, genetics, and molecular breeding. DRDB is freely available at http://drdb.big.ac.cn/home.

MicroRNA Expression in Multistage Date Fruit Development.

MicroRNAs (miRNAs) are small endogenous noncoding RNAs. Plant miRNAs are known to play important regulatory roles in homeostasis, stress response, and diverse developmental processes. Here, we describe the identification of conserved miRNAs in date palm (Phoenix dactylifera L.) based on transcriptomic data acquired across multistage fruit development and genome sequences, which include 238 plant conserved miRNAs and 276 novel P. dactylifera-specific miRNAs.

ISVASE: identification of sequence variant associated with splicing event using RNA-seq data.

BACKGROUND: Exon recognition and splicing precisely and efficiently by spliceosome is the key to generate mature mRNAs. About one third or a half of disease-related mutations affect RNA splicing. Software PVAAS has been developed to identify variants associated with aberrant splicing by directly using RNA-seq data. However, it bases on the assumption that annotated splicing site is normal splicing, which is not true in fact. RESULTS: We develop the ISVASE, a tool for specifically identifying sequence variants associated with splicing events (SVASE) by using RNA-seq data. Comparing with PVAAS, our tool has several advantages, such as multi-pass stringent rule-dependent filters and statistical filters, only using split-reads, independent sequence variant identification in each part of splicing (junction), sequence variant detection for both of known and novel splicing event, additional exon-exon junction shift event detection if known splicing events provided, splicing signal evaluation, known DNA mutation and/or RNA editing data supported, higher precision and consistency, and short running time. Using a realistic RNA-seq dataset, we performed a case study to illustrate the functionality and effectiveness of our method. Moreover, the output of SVASEs can be used for downstream analysis such as splicing regulatory element study and sequence variant functional analysis. CONCLUSIONS: ISVASE is useful for researchers interested in sequence variants (DNA mutation and/or RNA editing) associated with splicing events. The package is freely available at https://sourceforge.net/projects/isvase/ .

Complete Sequence and Analysis of Coconut Palm (Cocos nucifera) Mitochondrial Genome.

Coconut (Cocos nucifera L.), a member of the palm family (Arecaceae), is one of the most economically important crops in tropics, serving as an important source of food, drink, fuel, medicine, and construction material. Here we report an assembly of the coconut (C. nucifera, Oman local Tall cultivar) mitochondrial (mt) genome based on next-generation sequencing data. This genome, 678,653bp in length and 45.5% in GC content, encodes 72 proteins, 9 pseudogenes, 23 tRNAs, and 3 ribosomal RNAs. Within the assembly, we find that the chloroplast (cp) derived regions account for 5.07% of the total assembly length, including 13 proteins, 2 pseudogenes, and 11 tRNAs. The mt genome has a relatively large fraction of repeat content (17.26%), including both forward (tandem) and inverted (palindromic) repeats. Sequence variation analysis shows that the Ti/Tv ratio of the mt genome is lower as compared to that of the nuclear genome and neutral expectation. By combining public RNA-Seq data for coconut, we identify 734 RNA editing sites supported by at least two datasets. In summary, our data provides the second complete mt genome sequence in the family Arecaceae, essential for further investigations on mitochondrial biology of seed plants.

Identification and analysis of mouse non-coding RNA using transcriptome data.

Transcripts are expressed spatially and temporally and they are very complicated, precise and specific; however, most studies are focused on protein-coding related genes. Recently, massively parallel cDNA sequencing (RNA-seq) has emerged to be a new and promising tool for transcriptome research, and numbers of non-coding RNAs, especially lincRNAs, have been widely identified and well characterized as important regulators of diverse biological processes. In this study, we used ultra-deep RNA-seq data from 15 mouse tissues to study the diversity and dynamic of non-coding RNAs in mouse. Using our own criteria, we identified totally 16,249 non-coding genes (21,569 non-coding RNAs) in mouse. We annotated these non-coding RNAs by diverse properties and found non-coding RNAs are generally shorter, have fewer exons, express in lower level and are more strikingly tissue-specific compared with protein-coding genes. Moreover, these non-coding RNAs show significant enrichment with transcriptional initiation and elongation signals including histone modifications (H3K4me3, H3K27me3 and H3K36me3), RNAPII binding sites and CAGE tags. The gene set enrichment analysis (GSEA) result revealed several sets of lincRNAs associated with diverse biological processes such as immune effector process, muscle development and sexual reproduction. Taken together, this study provides a more comprehensive annotation of mouse non-coding RNAs and gives an opportunity for future functional and evolutionary study of mouse non-coding RNAs.

Identification and analysis of house-keeping and tissue-specific genes based on RNA-seq data sets across 15 mouse tissues.

Recently, RNA-seq has become widely used technology for transcriptome profiling due to its single-base accuracy and high-throughput speciality. In this study, we applied a computational approach on an integrated RNA-seq dataset across 15 normal mouse tissues, and consequently assigned 8408 house-keeping (HK) genes and 2581 tissue-specific (TS) genes among UCSC RefGene annotation. Apart from some basic genomic features, we also performed expression, function and pathway analysis with clustering, DAVID and Ingenuity Pathway Analysis, indicating the physiological connections (tissues) and diverse biological roles of HK genes (fundamental processes) and TS genes (tissue-corresponding processes). Moreover, we used RT-PCR method to test 18 candidate HK genes and finally identified a novel list of highly stable internal control genes: Ywhae, Ddb 1, Eif4h, etc. In summary, this study provides a new HK gene and TS gene resource for further genetic and evolution research and helps us better understand morphogenesis and biological diversity in mouse.

Profiling microRNA expression during multi-staged date palm (Phoenix dactylifera L.) fruit development.

MicroRNAs (miRNAs) play crucial roles in multiple stages of plant development and regulate gene expression at posttranscriptional and translational levels. In this study, we first identified 238 conserved miRNAs in date palm (Phoenix dactylifera) based on a high-quality genome assembly and defined 78 fruit-development-associated (FDA) miRNAs, whose expression profiles are variable at different fruit development stages. Using experimental data, we subsequently detected 276 novel P. dactylifera-specific FDA miRNAs and predicted their targets. We also revealed that FDA miRNAs function mainly in regulating genes involved in starch/sucrose metabolisms and other carbon metabolic pathways; among them, 221 FDA miRNAs exhibit negative correlation with their corresponding targets, which suggests their direct regulatory roles on mRNA targets. Our data define a comprehensive set of conserved and novel FDA miRNAs along with their expression profiles, which provide a basis for further experimentation in assigning discrete functions of these miRNAs in P. dactylifera fruit development.

Transcriptomic study of the red palm weevil Rhynchophorus ferrugineus embryogenesis.

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-ß, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

A large-scale gene discovery for the red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae).

The red palm weevil (RPW; Rhynchophorus ferrugineus) is a devastating pest of palms, prevalent in the Middle East as well as many other regions of the world. Here, we report a large-scale de novo complementary DNA (cDNA) sequencing effort that acquired ∼5 million reads and assembled them into 26 765 contigs from 12 libraries made from samples of different RPW developmental stages based on the Roche/454 GS FLX platform. We annotated these contigs based on the publically available known insect genes and the Tribolium castaneum genome assembly. We find that over 80% of coding sequences (CDS) from the RPW contigs have high-identity homologs to known proteins with complete CDS. Gene expression analysis shows that the pupa and larval stages have the highest and lowest expression levels, respectively. In addition, we also identified more than 60 000 single nucleotide polymorphisms and 1 200 simple sequence repeat markers. This study provides the first large-scale cDNA dataset for RPW, a much-needed resource for future molecular studies.

Genome sequence of the date palm Phoenix dactylifera L.

Date palm (Phoenix dactylifera L.) is a cultivated woody plant species with agricultural and economic importance. Here we report a genome assembly for an elite variety (Khalas), which is 605.4 Mb in size and covers >90% of the genome (~671 Mb) and >96% of its genes (~41,660 genes). Genomic sequence analysis demonstrates that P. dactylifera experienced a clear genome-wide duplication after either ancient whole genome duplications or massive segmental duplications. Genetic diversity analysis indicates that its stress resistance and sugar metabolism-related genes tend to be enriched in the chromosomal regions where the density of single-nucleotide polymorphisms is relatively low. Using transcriptomic data, we also illustrate the date palm's unique sugar metabolism that underlies fruit development and ripening. Our large-scale genomic and transcriptomic data pave the way for further genomic studies not only on P. dactylifera but also other Arecaceae plants.

Evolutional and functional analysis of a serine protease in Spodoptera litura.

Spodoptera litura is a threatening agricultural insect in tropical and subtropical areas and accounts for tremendous annual crop losses. As seen in virtually all insect species, serine proteases (SPs) are crucial to S. litura. The expression pattern of SPs from the midgut of S. litura was studied through expressed sequence tags (ESTs) analysis. One of SP (SlSP1) was chosen for detailed study, because the expression of the gene was midgut and larvae specific. SlSP1 was conducted as a model of its evolution, structure, and potential binding activity with corresponding substrates. SlSP1 is composed of 255 amino acids including a signal peptide at N-terminal followed by a putative activation peptide and the mature protein along with five putative phosphorylation sites, three disulphide bridges, and two N-glycosylation positions. At least nine conserved motifs were obtained in multiple sequence alignments. Some conserved residues, such as the catalytic triad His84, Asp127, and Ser229 as well as six cysteines at position 66, 82, 194, 211, 223, and 247, were examined. After homology modeling and molecular dynamics simulation, the resultant three-dimensional (3D) structure of SlSP1 was docked with the substrates 2PTC-Arg and 2PTC-Lys, respectively. Molecular Mechanic/Poisson-Boltzmann surface area analysis was applied to anticipate optimal binding mode and crucial active sites of this enzyme. The residues Trp28, Gly187, Aso188, Arg249, Ile250, Lys246, and Lys278 are crucial for the substrate binding and molecule process. This information can be used in logical design of SPs inhibitors. New inhibitors may be a basis for development of a new pest control technology.