Cloning, sequencing, and characterization of a heat- and alkali-stable type I pullulanase from Anaerobranca gottschalkii

The gene encoding a type I pullulanase was identified from the genome sequence of the anaerobic thermoalkaliphilic bacterium Anaerobranca gottschalkii. In addition, the homologous gene was isolated from a gene library of Anaerobranca horikoshii and sequenced. The proteins encoded by these two genes showed 39% amino acid sequence identity to the pullulanases from the thermophilic anaerobic bacteria Fervidobacterium pennivorans and Thermotoga maritima. The pullulanase gene from A. gottschalkii (encoding 865 amino acids with a predicted molecular mass of 98 kDa) was cloned and expressed in Escherichia coli strain BL21(DE3) so that the protein did not have the signal peptide. Accordingly, the molecular mass of the purified recombinant pullulanase (rPulAg) was 96 kDa. Pullulan hydrolysis activity was optimal at pH 8.0 and 70 degrees C, and under these physicochemical conditions the half-life of rPulAg was 22 h. By using an alternative expression strategy in E. coli Tuner(DE3)(pLysS), the pullulanase gene from A. gottschalkii, including its signal peptide-encoding sequence, was cloned. In this case, the purified recombinant enzyme was a truncated 70-kDa form (rPulAg'). The N-terminal sequence of purified rPulAg' was found 252 amino acids downstream from the start site, presumably indicating that there was alternative translation initiation or N-terminal protease cleavage by E. coli. Interestingly, most of the physicochemical properties of rPulAg' were identical to those of rPulAg. Both enzymes degraded pullulan via an endo-type mechanism, yielding maltotriose as the final product, and hydrolytic activity was also detected with amylopectin, starch, beta-limited dextrins, and glycogen but not with amylose. This substrate specificity is typical of type I pullulanases. rPulAg was inhibited by cyclodextrins, whereas addition of mono- or bivalent cations did not have a stimulating effect. In addition, rPulAg' was stable in the presence of 0.5% sodium dodecyl sulfate, 20% Tween, and 50% Triton X-100. The pullulanase from A. gottschalkii is the first thermoalkalistable type I pullulanase that has been described.     

AmyA, an alpha-amylase with beta-cyclodextrin-forming activity, and AmyB from the thermoalkaliphilic organism Anaerobranca gottschalkii: two alpha-amylases adapted to their different cellular localizations

Two alpha-amylase genes from the thermophilic alkaliphile Anaerobranca gottschalkii were cloned, and the corresponding enzymes, AmyA and AmyB, were investigated after purification of the recombinant proteins. Based on their amino acid sequences, AmyA is proposed to be a lipoprotein with extracellular localization and thus is exposed to the alkaline milieu, while AmyB apparently represents a cytoplasmic enzyme. The amino acid sequences of both enzymes bear high similarity to those of GHF13 proteins. The different cellular localizations of AmyA and AmyB are reflected in their physicochemical properties. The alkaline pH optimum (pH 8), as well as the broad pH range, of AmyA activity (more than 50% activity between pH 6 and pH 9.5) mirrors the conditions that are encountered by an extracellular enzyme exposed to the medium of A. gottschalkii, which grows between pH 6 and pH 10.5. AmyB, on the other hand, has a narrow pH range with a slightly acidic pH optimum at 6 to 6.5, which is presumably close to the pH in the cytoplasm. Also, the intracellular AmyB is less tolerant of high temperatures than the extracellular AmyA. While AmyA has a half-life of 48 h at 70 degrees C, AmyB has a half-life of only about 10 min at that temperature, perhaps due to the lack of stabilizing constituents of the cytoplasm. AmyA and AmyB were very similar with respect to their substrate specificity profiles, clearly preferring amylose over amylopectin, pullulan, and glycogen. Both enzymes also hydrolyzed alpha-, beta-, and gamma-cyclodextrin. Very interestingly, AmyA, but not AmyB, displayed high transglycosylation activity on maltooligosaccharides and also had significant beta-cyclodextrin glycosyltransferase (CGTase) activity. CGTase activity has not been reported for typical alpha-amylases before. The mechanism of cyclodextrin formation by AmyA is unknown.     

Capturing single-copy nuclear genes,organellar genomes,and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae

With the decreasing cost and availability of many newly developed bioinformatics pipelines, next-generation sequencing (NGS) has revolutionized plant systematics in recent years. Genome skimming has been widely used to obtain high-copy fractions of the genomes, including plastomes, mitochondrial DNA (mtDNA), and nuclear ribosomal DNA (nrDNA). In this study, through simulations, we evaluated the optimal (minimum) sequencing depth and performance for recovering single-copy nuclear genes (SCNs) from genome skimming data, by subsampling genome resequencing data and generating 10 data sets with different sequencing coverage in silico. We tested the performance of four data sets (plastome, nrDNA, mtDNA, and SCNs) obtained from genome skimming based on phylogenetic analyses of the Vitis clade at the genus level and Vitaceae at the family level, respectively. Our results showed that optimal minimum sequencing depth for high-quality SCNs assembly via genome skimming was about 10× coverage. Without the steps of synthesizing baits and enrichment experiments, coupled with incredibly low sequencing costs, we showcase that deep genome skimming (DGS) is as effective for capturing large data sets of SCNs as the widely used Hyb-Seq approach, in addition to capturing plastomes, mtDNA, and entire nrDNA repeats. DGS may serve as an efficient and economical alternative and may be superior to the popular target enrichment/Hyb-Seq approach.     

How can we deliver the large plant genomes? Strategies and perspectives

The first sequenced plant genome, from the small mustard plant Arabidopsis thaliana, was published at the end of 2000. The sequencing of the rice genome is well under way. The sizes of plant genomes vary by a factor of up to 1000, and many important crop plants have genomes that are several times larger than the human genome. To gain insight into the gene toolbox of plant species, numerous large-scale EST sequencing projects have been launched successfully, and analysis procedures are constantly being refined to add maximum value to the sequence data. In addition, an alternative approach to exclude repetitive noncoding DNA and to enrich sequence libraries for gene-containing genomic regions has been developed. This strategy has the potential to deliver information about both genes and regulatory regions outside the transcribed regions.     

Gene identification using rice genome sequences

Identifying useful gene(s) is one of the most important objectives of plant geneticists. Various strategies can be used, which are based on the characteristics of plant reproduction and available technology. Rice is the first model crop whose whole genome sequence has been reported. In addition, information on the whole genome sequences of two important rice subspecies (japonica and indica rice) is also available. Rice is a self-pollinating crop and methods of artificial crossing are relatively easy to perform; such methods enable the production of numerous seeds for genetic analyses. Based on these features, a map-based cloning (i.e., positional cloning) strategy has been successfully applied over the last decade to identify rice genes. Recently, advanced next-generation sequencing (NGS) technology was used to ascertain the genome sequences of individual plants, opening up a new strategy for gene identification. This strategy has been used successfully to identify the genes responsible for certain qualitative traits in rice. However, to identify the gene(s) involved in a quantitative trait, a map-based cloning strategy is still required after quantitative trait loci analysis using NGS technology. In this review, we discuss both map-based cloning (which is still the primary strategy used to identify rice genes) and NGS-based strategies.     

De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an edible root vegetable crop that is cultivated worldwide and whose genome has been sequenced. Here we report the complete nucleotide sequence of the radish cultivar WK10039 chloroplast (cp) genome, along with a de novo assembly strategy using whole genome shotgun sequence reads obtained by next generation sequencing. The radish cp genome is 153,368 bp in length and has a typical quadripartite structure, composed of a pair of inverted repeat regions (26,217 bp each), a large single copy region (83,170 bp), and a small single copy region (17,764 bp). The radish cp genome contains 87 predicted protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Sequence analysis revealed the presence of 91 simple sequence repeats (SSRs) in the radish cp genome.     

Complete genome sequence of Paenibacillus sp. strain JDR-2

Paenibacillus sp. strain JDR-2, an aggressively xylanolytic bacterium isolated from sweetgum (Liquidambar styraciflua) wood, is able to efficiently depolymerize, assimilate and metabolize 4-O-methylglucuronoxylan, the predominant structural component of hardwood hemicelluloses. A basis for this capability was first supported by the identification of genes and characterization of encoded enzymes and has been further defined by the sequencing and annotation of the complete genome, which we describe. In addition to genes implicated in the utilization of β-1,4-xylan, genes have also been identified for the utilization of other hemicellulosic polysaccharides. The genome of Paenibacillus sp. JDR-2 contains 7,184,930 bp in a single replicon with 6,288 protein-coding and 122 RNA genes. Uniquely prominent are 874 genes encoding proteins involved in carbohydrate transport and metabolism. The prevalence and organization of these genes support a metabolic potential for bioprocessing of hemicellulose fractions derived from lignocellulosic resources.     

Sequencing XMET genes to promote genotype-guided risk assessment and precision medicine

High-throughput next generation sequencing(NGS) is a shotgun approach applied in a parallel fashion by which the genome is fragmented and sequenced through small pieces and then analyzed either by aligning to a known reference genome or by de novo assembly without reference genome. This technology has led researchers to conduct an explosion of sequencing related projects in multidisciplinary fields of science. However, due to the limitations of sequencing-based chemistry, length of sequencing reads and the complexity of genes, it is difficult to determine the sequences of some portions of the human genome, leaving gaps in genomic data that frustrate further analysis. Particularly, some complex genes are difficult to be accurately sequenced or mapped because they contain high GC-content and/or low complexity regions, and complicated pseudogenes, such as the genes encoding xenobiotic metabolizing enzymes and transporters(XMETs). The genetic variants in XMET genes are critical to predicate interindividual variability in drug efficacy, drug safety and susceptibility to environmental toxicity. We summarized and discussed challenges, wet-lab methods, and bioinformatics algorithms in sequencing "complex" XMET genes, which may provide insightful information in the application of NGS technology for implementation in toxicogenomics and pharmacogenomics.     

Targeted analysis of orthologous phytochrome A regions of the sorghum,maize, and rice genomes using comparative gene-island sequencing

A "gene-island" sequencing strategy has been developed that expedites the targeted acquisition of orthologous gene sequences from related species for comparative genome analysis. A 152-kb bacterial artificial chromosome (BAC) clone from sorghum (Sorghum bicolor) encoding phytochrome A (PHYA) was fully sequenced, revealing 16 open reading frames with a gene density similar to many regions of the rice (Oryza sativa) genome. The sequences of genes in the orthologous region of the maize (Zea mays) and rice genomes were obtained using the gene-island sequencing method. BAC clones containing the orthologous maize and rice PHYA genes were identified, sheared, subcloned, and probed with the sorghum PHYA-containing BAC DNA. Sequence analysis revealed that approximately 75% of the cross-hybridizing subclones contained sequences orthologous to those within the sorghum PHYA BAC and less than 25% contained repetitive and/or BAC vector DNA sequences. The complete sequence of four genes, including up to 1 kb of their promoter regions, was identified in the maize PHYA BAC. Nine orthologous gene sequences were identified in the rice PHYA BAC. Sequence comparison of the orthologous sorghum and maize genes aided in the identification of exons and conserved regulatory sequences flanking each open reading frame. Within genomic regions where micro-colinearity of genes is absolutely conserved, gene-island sequencing is a particularly useful tool for comparative analysis of genomes between related species.     

Comparative Genome Analysis of Wheat Blue Dwarf Phytoplasma,an Obligate Pathogen That Causes Wheat Blue Dwarf Disease in China

Wheat blue dwarf (WBD) disease is an important disease that has caused heavy losses in wheat production in northwestern China. This disease is caused by WBD phytoplasma, which is transmitted by Psammotettix striatus. Until now, no genome information about WBD phytoplasma has been published, seriously restricting research on this obligate pathogen. In this paper, we report a new sequencing and assembling strategy for phytoplasma genome projects. This strategy involves differential centrifugation, pulsed-field gel electrophoresis, whole genome amplification, shotgun sequencing, de novo assembly, screening of contigs from phytoplasma and the connection of phytoplasma contigs. Using this scheme, the WBD phytoplasma draft genome was obtained. It was comprised of six contigs with a total size of 611,462 bp, covering ∼94% of the chromosome. Five-hundred-twenty-five protein-coding genes, two operons for rRNA genes and 32 tRNA genes were identified. Comparative genome analyses between WBD phytoplasma and other phytoplasmas were subsequently carried out. The results showed that extensive arrangements and inversions existed among the WBD, OY-M and AY-WB phytoplasma genomes. Most protein-coding genes in WBD phytoplasma were found to be homologous to genes from other phytoplasmas; only 22 WBD-specific genes were identified. KEGG pathway analysis indicated that WBD phytoplasma had strongly reduced metabolic capabilities. However, 46 transporters were identified, which were involved with dipeptides/oligopeptides, spermidine/putrescine, cobalt and Mn/Zn transport, and so on. A total of 37 secreted proteins were encoded in the WBD phytoplasma chromosome and plasmids. Of these, three secreted proteins were similar to the reported phytoplasma virulence factors TENGU, SAP11 and SAP54. In addition, WBD phytoplasma possessed several proteins that were predicted to play a role in its adaptation to diverse environments. These results will provide clues for research on the pathogenic mechanisms of WBD phytoplasma and will also provide a perspective about the genome sequencing of other phytoplasmas and obligate organisms.     

TILLING - a shortcut in functional genomics

Recent advances in large-scale genome sequencing projects have opened up new possibilities for the application of conventional mutation techniques in not only forward but also reverse genetics strategies. TILLING (Targeting Induced Local Lesions IN Genomes) was developed a decade ago as an alternative to insertional mutagenesis. It takes advantage of classical mutagenesis, sequence availability and high-throughput screening for nucleotide polymorphisms in a targeted sequence. The main advantage of TILLING as a reverse genetics strategy is that it can be applied to any species, regardless of its genome size and ploidy level. The TILLING protocol provides a high frequency of point mutations distributed randomly in the genome. The great mutagenic potential of chemical agents to generate a high rate of nucleotide substitutions has been proven by the high density of mutations reported for TILLING populations in various plant species. For most of them, the analysis of several genes revealed 1 mutation/200–500 kb screened and much higher densities were observed for polyploid species, such as wheat. High-throughput TILLING permits the rapid and low-cost discovery of new alleles that are induced in plants. Several research centres have established a TILLING public service for various plant species. The recent trends in TILLING procedures rely on the diversification of bioinformatic tools, new methods of mutation detection, including mismatch-specific and sensitive endonucleases, but also various alternatives for LI-COR screening and single nucleotide polymorphism (SNP) discovery using next-generation sequencing technologies. The TILLING strategy has found numerous applications in functional genomics. Additionally, wide applications of this throughput method in basic and applied research have already been implemented through modifications of the original TILLING strategy, such as Ecotilling or Deletion TILLING.     

Chemical transmission in the sea anemone Nematostella vectensis: A genomic perspective

The sequencing of the starlet sea anemone (Nematostella vectensis) genome provides opportunities to investigate the function and evolution of genes associated with chemical neurotransmission and hormonal signaling. This is of particular interest because sea anemones are anthozoans, the phylogenetically basal cnidarians least changed from the common ancestors of cnidarians and bilaterian animals, and because cnidarians are considered the most basal metazoans possessing a nervous system. This analysis of the genome has yielded 20 orthologues of enzymes and nicotinic receptors associated with cholinergic function, an even larger number of genes encoding enzymes, receptors and transporters for glutamatergic (28) and GABAergic (34) transmission, and two orthologues of purinergic receptors. Numerous genes encoding enzymes (14), receptors (60) and transporters (5) for aminergic transmission were identified, along with four adenosine-like receptors and one nitric oxide synthase. Diverse neuropeptide and hormone families are also represented, mostly with genes encoding prepropeptides and receptors related to varying closeness to RFamide (17) and tachykinin (14), but also galanin (8), gonadotropin-releasing hormones and vasopressin/oxytocin (5), melanocortins (11), insulin-like peptides (5), glycoprotein hormones (7), and uniquely cnidarian peptide families (44). Surprisingly, no muscarinic acetylcholine receptors were identified and a large number of melatonin-related, but not serotonin, orthologues were found. Phylogenetic tree construction and inspection of multiple sequence alignments reveal how evolutionarily and functionally distant chemical transmitter-related proteins are from those of higher metazoans.     

Design of efficient simplified genomic DNA and bisulfite sequencing in large plant populations

The next generation sequencing enables generation of high resolution and high throughput data for structure sequence of any genome at a fast declining cost. This opens opportunity for population based genetic and genomic analyses. In many applications, whole genome sequencing or re-sequencing is unnecessary or prohibited by budget limits. The Reduced Representation Genome Sequencing (RRGS), which sequences only a small proportion of the genome of interest, has been proposed to deal with the situations. Several forms of RRGS are proposed and implemented in the literature. When applied to plant or crop species, the current RRGS protocols shared a key drawback that a significantly high proportion (up to 60%) of sequence reads to be generated may be of non-genomic origin but attributed to chloroplast DNA or rRNA genes, leaving an exceptional low efficiency of the sequencing experiment. We recommended and discussed here the design of optimized simplified genomic DNA and bisulfite sequencing strategies, which may greatly improves efficiency of the sequencing experiments by bringing down the presentation of the undesirable sequencing reads to less than 10% in the whole sequence reads. The optimized RAD-seq and RRBS-seq methods are potentially useful for sequence variant screening and genotyping in large plant/crop populations.     

全基因组测序及其在遗传性疾病研究及诊断中的应用

最近,随着测序成本的不断降低,数据分析策略的不断提升,全基因组测序（whole-genome sequencing,WGS）已经在癌症、孟德尔遗传病、复杂疾病的致病基因检测中得到了一定运用,并逐步走向了临床诊断。全基因组测序不但可以检测编码区和非编码区的点突变（SNVs）和插入缺失（InDels）,还可以在全基因组范围内检测拷贝数变异（copy number variation,CNV）以及结构变异（structure variation,SV）。本文详细地介绍了全基因组测序的标准生物信息分析流程与方法,及其在疾病研究、临床诊断中的应用,并对全基因组测序在医学遗传学中的应用与研究进展,以及数据分析方面面临的挑战进行了概述。     

Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome

The application of a new gene-based strategy for sequencing the wheat mitochondrial genome shows its structure to be a 452528 bp circular molecule, and provides nucleotide-level evidence of intra-molecular recombination. Single, reciprocal and double recombinant products, and the nucleotide sequences of the repeats that mediate their formation have been identified. The genome has 55 genes with exons, including 35 protein-coding, 3 rRNA and 17 tRNA genes. Nucleotide sequences of seven wheat genes have been determined here for the first time. Nine genes have an exon–intron structure. Gene amplification responsible for the production of multicopy mitochondrial genes, in general, is species-specific, suggesting the recent origin of these genes. About 16, 17, 15, 3.0 and 0.2% of wheat mitochondrial DNA (mtDNA) may be of genic (including introns), open reading frame, repetitive sequence, chloroplast and retro-element origin, respectively. The gene order of the wheat mitochondrial gene map shows little synteny to the rice and maize maps, indicative that thorough gene shuffling occurred during speciation. Almost all unique mtDNA sequences of wheat, as compared with rice and maize mtDNAs, are redundant DNA. Features of the gene-based strategy are discussed, and a mechanistic model of mitochondrial gene amplification is proposed.     

The malaria genome sequencing project: complete sequence of Plasmodium falciparum chromosome 2

An international consortium has been formed to sequence the entire genome of the human malaria parasite Plasmodium falciparum. We sequenced chromosome 2 of clone 3D7 using a shotgun sequencing strategy. Chromosome 2 is 947 kb in length, has a base composition of 80.2% A + T, and contains 210 predicted genes. In comparison to the Saccharomyces cerevisiae genome, chromosome 2 has a lower gene density, a greater proportion of genes containing introns, and nearly twice as many proteins containing predicted non-globular domains. A group of putative surface proteins was identified, rifins, which are encoded by a gene family comprising up to 7% of the protein-encoding gene in the genome. The rifins exhibit considerable sequence diversity and may play an important role in antigenic variation. Sixteen genes encoded on chromosome 2 showed signs of a plastid or mitochondrial origin, including several genes involved in fatty acid biosynthesis. Completion of the chromosome 2 sequence demonstrated that the A + T-rich genome of P. falciparum can be sequenced by the shotgun approach. Within 2-3 years, the sequence of almost all P. falciparum genes will have been determined, paving the way for genetic, biochemical, and immunological research aimed at developing new drugs and vaccines against malaria.     

利用木质纤维素类生物质产微生物絮凝剂Pseudomonas boreopolis GO2的全基因组测序及比较基因组学分析

【目的】Pseudomonas boreopolis GO2可以利用木质纤维素类生物质为唯一碳源发酵产微生物絮凝剂。解析菌株GO2的全基因组特征可为利用木质纤维素类生物质定向合成多糖型微生物絮凝剂提供分子基础。【方法】利用Illumina NovaSeq测序平台对菌株GO2进行测序,用SMRT等软件进行基因组组装、系统发育分析、基因预测和功能注释,并与4株近缘模式株进行了比较基因组分析。【结果】菌株GO2基因组大小为4 498 896 bp,GC含量为69.5%,共编码3 906个基因。菌株GO2与Pseudomonas boreopolis JCM 13306的16S r RNA基因相似性、平均核苷酸一致性(average nucleotide identity, ANI)、DNA-DNA杂交(DNA-DNA hybridization, DDH)值最高,分别为99.93%、98.36%和88.00%,将菌株GO2命名为Pseudomonas boreopolis GO2。比较基因组分析发现,GO2与4个近缘模式菌株共有2 348个直系同源核心基因,主要参与碳水化合物代谢、氨基酸代谢...     

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene‐rich regions. Gene‐enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene‐enrichment strategy, we have compared assemblies using methyl‐filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA‐seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single‐nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.