Genome Sequence Information Research Articles

Complete Genome Sequence of Bradyrhizobium sp. ORS285, a Photosynthetic Strain Able To Establish Nod Factor-Dependent or Nod Factor-Independent Symbiosis with Aeschynomene Legumes.

ABSTRACTHere, we report the complete genome sequence of Bradyrhizobium sp. strain ORS285, which is able to nodulate Aeschynomene legumes using two distinct strategies that differ in the requirement of Nod factors. The genome sequence information of this strain will help understanding of the different mechanisms of interaction of rhizobia with legumes.

Detection of Verticillium species in Swedish soils using real-time PCR

Verticillium species are soilborne plant pathogens, responsible for big yield losses worldwide. Here, we report improved procedures to generate DNA from Verticillium species imbedded in farm soils. Using new genomic sequence information, primers for V. dahliae, V. albo-atrum, V. tricorpus, and V. longisporum were designed. In a survey of 429 samples from intensively farmed soil of two Swedish regions, only V. dahliae and V. longisporum were identified. A bias towards V. longisporum (40%) was seen in the south, whereas V. dahliae was more frequent in the western region (19%). Analyses of soil and leaf samples from 20 sugar beet fields, where foliar wilting had been observed, revealed V. dahliae DNA in all leaf and soil samples and V. longisporum in 18 soil samples, illustrating host choice and longevity of the V. longisporum microsclerotia. This study demonstrates the applicability of new molecular diagnostic tools that are important for growers of variable crops.

Complete Genome Sequence of Burkholderia stabilis FERMP-21014.

ABSTRACTCholesterol esterase (EC 3.1.1.13) was identified in a bacterium, Burkholderia stabilis strain FERMP-21014. Here, we report the complete genome sequence of B. stabilis FERMP-21014, which has been used in the commercial production of cholesterol esterase. The genome sequence information may be useful for improving production levels of cholesterol esterase.

Retrotransposon distribution and copy number variation in gymnosperm genomes

Retrotransposable elements (REs) and related sequences form a large proportion of conifer genomes. During genome evolution, some RE sequences are degraded or eliminated, but some are evolutionarily stable, and can be identified even in distantly related species. Use of genome sequence information from loblolly pine (Pinus taeda) enables investigation of divergent non-coding RE sequences in other pine and conifer species, including Scots pine (Pinus sylvestris). Non-specific inter-retrotransposon amplified polymorphism technique (IRAP) as well as the amplification polymorphism of 12 RE families were investigated in 80 gymnosperm species. The obtained results were compared with phylogenetic relationships among gymnosperms. Investigation of distantly related gymnosperm species reveals persistent RE sequences, such as IFG and Pineywoods, distributed among a wide range of plant lineages. RE sequence divergence was observed, reflecting periods of inactivity and degradation during speciation of pine lineages, as demonstrated by the delineation of the main pine subgenera. Intraspecific variation of 10 RE copy numbers (CN) between Scots pine individuals ranged from 8.9 to 26.6% of the overall mean estimates. CN analyses were performed in 16 additional gymnosperm species. The analysed pine species contained a similar complement of RE families; however, CN and genome occupation proportions differ. A decrease in RE CN estimates can reflect sequence divergence, associated with independent transposition events. Transposition of some REs can be induced by stress conditions; therefore, even distantly related species inhabiting extreme environments could have similar patterns or distribution of these elements.

Abstract 2590: Analysis of APOBEC3A and APOBEC3B mutational signatures using next-generation sequencing data from cancer cell lines

Abstract The APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) gene family of cytidine deaminases includes evolutionarily conserved genes that play important roles in DNA repair and mRNA editing. Activity of at least two APOBEC family members, APOBEC3A and APOBEC3B, can lead to kataegis, a mutagenic process in cancer cells that generates clusters of closely spaced, single strand specific C->T DNA substitutions. APOBEC mutagenesis has a characteristic signature, most commonly represented by the 5’-Tp(C->T)pW-3’ sequence motif, with additional substitutions also reported. This hypermutation signature and high mRNA expression of APOBEC3A and APOBEC3B have been associated with several cancer types. Most previous studies of APOBEC signatures have examined tumor sequence data from clinical samples, for which limited or no information about drug response was available. We investigated the presence of the mutational signature and mRNA expression patterns of the APOBEC3A and APOBEC3B genes in extensively characterized cell lines, in order to identify those cell lines that carry mutations generated by kataegis, with the aim of establishing associations between the APOBEC mutational signature, individual cancer types, and the patterns of sensitivity to antitumor agents. For this purpose, we analyzed whole exome sequencing (WES) data and mRNA expression of the APOBEC3A and APOBEC3B genes in two resources with extensive drug response data: the NCI-60 cell line panel, which includes 59 human cancer cell lines representing 9 cancer types and drug response information for thousands of anticancer agents, and the Cancer Cell Line Encyclopedia (CCLE), which provides WES, whole genome, and RNA-seq sequence information on hundreds of cancer cell lines and drug response data to over 200 agents. We analyzed WES data of 325 CCLE cell lines and 59 NCI-60 cell lines, with variants identified using GATK pipeline and Varscan2 software. The variants in each cell line were filtered to remove common polymorphisms in dbSNP and 1000 Genome Project databases. We searched the discovered variants for the presence of APOBEC signatures, 5’-Tp(C->AGT)pN-3’, 5’-Tp(C->AGT)pD-3’, and 5’-Tp(C->AGT)pW-3’ in closely spaced (1000 and 10,000 bp) windows that appeared on the same DNA strand. We will discuss the use of optimal filters for detecting APOBEC mutational signatures and will present the analyses of associations between APOBEC signatures, mutational load of the tumor cell lines, APOBEC gene expression, and chemosensitivity to treatment. These results contribute to additional characterization of available cell lines by providing information about specific mutational signatures in different categories of cancer. Our findings may assist with identifying antitumor agents that would be appropriate for treatment of cancer cells with specific signature patterns generated by APOBEC mutagenesis. Citation Format: Suleyman Vural, Julia Krushkal, Richard Simon. Analysis of APOBEC3A and APOBEC3B mutational signatures using next-generation sequencing data from cancer cell lines [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2590. doi:10.1158/1538-7445.AM2017-2590

Draft Genome Report of Bacillus altitudinis SORB11, Isolated from the Indian Sector of the Southern Ocean.

ABSTRACTHere, we present the draft genome sequence of Bacillus altitudinis SORB11, which is tolerant to UV radiation. The strain was isolated from the Indian sector of the Southern Ocean at a depth of 3.8 km. The genome sequence information reported here for B. altitudinis SORB11 gives the basis of its UV resistance mechanism and provides data for further comparative studies with other bacteria resistant to UV radiation.

Human and mouse microarrays-guided expression analysis of membrane protein trafficking-related genes in MDCK cells, a canine epithelial model for apical and basolateral differential protein targeting

MDCK cells are widely used to study the differential targeting of membrane transporters to apical and basolateral membrane but its canine origin limited the commercial tools available for the analysis of protein trafficking machinery. Because apical and basolateral membranes are only found in differentiated epithelial cells, genes critical for differential targeting may be specifically up-regulated upon MDCK cell differentiation. To search for these genes, a cross-species screening strategy was used. We first analyzed the human microarray data for protein trafficking-related genes that were up-regulated in colon carcinoma Caco2 cells upon differentiation. The results of mouse 44K gene expression microarray analysis were then used to extract additional candidate genes that showed higher expression in normal colon epithelium compared to primary embryonic fibroblasts. Finally, NCBI genomic sequence information was used to design RT-PCR primers for 13 candidate and 10 negative control genes and used to analyze MDCK cells at 2, 13 and 17 days after seeding. To determine whether the gene up-regulation was specific in epithelial differentiation, we also performed RT-PCR on rat non-differentiating intestinal IEC-6 cells and mouse C2C12 cells, a differentiating myoblast model. Of the 13 candidate genes, 3 genes, SDCBP2, KIF12, KIF27, met all criteria of specific up-regulation in differentiated MDCK cells. In addition, KIF13A showed up-regulation in differentiated MDCK and C2C12 cells but not in IEC-6 cells cultured for the same duration. The functions of these genes need to be analyzed in the future. This cross-species screening strategy may be useful for other non-human, non-rodent cell models.

Unveiling the Hybrid Genome Structure of Escherichia coli RR1 (HB101 RecA+).

There have been extensive genome sequencing studies for Escherichia coli strains, particularly for pathogenic isolates, because fast determination of pathogenic potential and/or drug resistance and their propagation routes is crucial. For laboratory E. coli strains, however, genome sequence information is limited except for several well-known strains. We determined the complete genome sequence of laboratory E. coli strain RR1 (HB101 RecA+), which has long been used as a general cloning host. A hybrid genome sequence of K-12 MG1655 and B BL21(DE3) was constructed based on the initial mapping of Illumina HiSeq reads to each reference, and iterative rounds of read mapping, variant detection, and consensus extraction were carried out. Finally, PCR and Sanger sequencing-based finishing were applied to resolve non-single nucleotide variant regions with aberrant read depths and breakpoints, most of them resulting from prophages and insertion sequence transpositions that are not present in the reference genome sequence. We found that 96.9% of the RR1 genome is derived from K-12, and identified exact crossover junctions between K-12 and B genomic fragments. However, because RR1 has experienced a series of genetic manipulations since branching from the common ancestor, it has a set of mutations different from those found in K-12 MG1655. As well as identifying all known genotypes of RR1 on the basis of genomic context, we found novel mutations. Our results extend current knowledge of the genotype of RR1 and its relatives, and provide insights into the pedigree, genomic background, and physiology of common laboratory strains.

Genome-Guided Insights into the Plant Growth Promotion Capabilities of the Physiologically Versatile Bacillus aryabhattai Strain AB211.

Bacillus aryabhattai AB211 is a plant growth promoting, Gram-positive firmicute, isolated from the rhizosphere of tea (Camellia sinensis), one of the oldest perennial crops and a major non-alcoholic beverage widely consumed all over the world. The whole genome of B. aryabhattai AB211 was sequenced, annotated and evaluated with special focus on genomic elements related to plant microbe interaction. It’s genome sequence reveals the presence of a 5,403,026 bp chromosome. A total of 5226 putative protein-coding sequences, 16 rRNA, 120 tRNA, 8 ncRNAs, 58 non-protein coding genes, and 11 prophage regions were identified. Genome sequence comparisons between strain AB211 and other related environmental strains of B. aryabhattai, identified about 3558 genes conserved among all B. aryabhattai genomes analyzed. Most of the common genes involved in plant growth promotion activities were found to be present within core genes of all the genomes used for comparison, illustrating possible common plant growth promoting traits shared among all the strains of B. aryabhattai. Besides the core genes, some genes were exclusively identified in the genome of strain AB211. Functional annotation of the genes predicted in the strain AB211 revealed the presence of genes responsible for mineral phosphate solubilization, siderophores, acetoin, butanediol, exopolysaccharides, flagella biosynthesis, surface attachment/biofilm formation, and indole acetic acid production, most of which were experimentally verified in the present study. Genome analysis and experimental evidence suggested that AB211 has robust central carbohydrate metabolism implying that this bacterium can efficiently utilize the root exudates and other organic materials as an energy source. Genes for the production of peroxidases, catalases, and superoxide dismutases, that confer resistance to oxidative stresses in plants were identified in AB211 genome. Besides these, genes for heat shock tolerance, cold shock tolerance, glycine-betaine production, and antibiotic/heavy metal resistance that enable bacteria to survive biotic/abiotic stress were also identified. Based on the genome sequence information and experimental evidence as presented in this study, strain AB211 appears to be metabolically diverse and exhibits tremendous potential as a plant growth promoting bacterium.

TILLING, high-resolution melting (HRM), and next-generation sequencing (NGS) techniques in plant mutation breeding

Induced mutations have been used effectively for plant improvement. Physical and chemical mutagens induce a high frequency of genome variation. Recently, developed screening methods have allowed the detection of single nucleotide polymorphisms (SNPs) and the identification of traits that are difficult to identify at the molecular level by conventional breeding. With the assistance of reverse genetic techniques, sequence variation information can be linked to traits to investigate gene function. Targeting induced local lesions in genomes (TILLING) is a high-throughput technique to identify single nucleotide mutations in a specific region of a gene of interest with a powerful detection method resulted from chemical-induced mutagenesis. The main advantage of TILLING as a reverse genetics strategy is that it can be applied to any species, regardless of genome size and ploidy level. However, TILLING requires laborious and time-consuming steps, and a lack of complete genome sequence information for many crop species has slowed the development of suitable TILLING targets. Another method, high-resolution melting (HRM), which has assisted TILLING in mutation detection, is faster, simpler and less expensive with non-enzymatic screening system. Currently, the sequencing of crop genomes has completely changed our vision and interpretation of genome organization and evolution. Impressive progress in next-generation sequencing (NGS) technologies has paved the way for the detection and exploitation of genetic variation in a given DNA or RNA molecule. This review discusses the applications of TILLING in combination with HRM and NGS technologies for screening of induced mutations and discovering SNPs in mutation breeding programs.

Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation.

Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.

Fine-mapping of genes determining extrafusal fiber properties in murine soleus muscle.

Muscle fiber cross-sectional area (CSA) and proportion of different fiber types are important determinants of muscle function and overall metabolism. Genetic variation plays a substantial role in phenotypic variation of these traits; however, the underlying genes remain poorly understood. This study aimed to map quantitative trait loci (QTL) affecting differences in soleus muscle fiber traits between the LG/J and SM/J mouse strains. Fiber number, CSA, and proportion of oxidative type I fibers were assessed in the soleus of 334 genotyped female and male mice of the F34 generation of advanced intercross lines (AIL) derived from the LG/J and SM/J strains. To increase the QTL detection power, these data were combined with 94 soleus samples from the F2 intercross of the same strains. Transcriptome of the soleus muscle of LG/J and SM/J females was analyzed by microarray. Genome-wide association analysis mapped four QTL (genome-wide P < 0.05) affecting the properties of muscle fibers to chromosome 2, 3, 4, and 11. A 1.5-LOD QTL support interval ranged between 2.36 and 4.67 Mb. On the basis of the genomic sequence information and functional and transcriptome data, we identified candidate genes for each of these QTL. The combination of analyses in F2 and F34 AIL populations with transcriptome and genomic sequence data in the parental strains is an effective strategy for refining QTL and nomination of the candidate genes.

DIRECT ESTIMATION OF GENETIC LIABILITY FOR MENTAL ILLNESS USING ENDOPHENOTYPES

The concept of continuous liability underlies most inferential procedures that are used in complex disease genetics. Statistical genetic methods often attempt to take dichotomous affection status information and project it onto the unobservable liability scale. For substance abuse or any other disorder, this inference about individual genetic liability is often coarse unless additional information can be incorporated. In this talk, I propose a model that leads to improved estimation of individual genetic liability for substance abuse. It incorporates quantitative endophenotypes that are both heritable and genetically correlated with the diagnosis of substance abuse. Family-based designs are optimal for inferences regarding genetic liability since substantial additional information can be drawn from correlated information known about the liability distribution in related individuals. After reviewing the statistical framework required for the multivariate polygenic architecture of affection status and endophenotypes, I develop a best linear unbiased prediction scheme for estimating the genetic liability underlying each endophenotype. Based upon the observed genetic covariance matrix, a multivariate informed single index of genetic liability can be estimated. This index is theoretically largely free of environmental influence and should better represent the true genetic liability for each individual. The resulting liability index can be used as the core phenotype for subsequent analyses searching for specific causal genes and functional genetic variants of relevance to psychiatric disease. As an example, I utilize data from 1,628 Mexican American individuals who are members of large extended pedigrees from the Genetics of Brain Structure and Function Study. Estimates of individual genetic liability for major depressive disorder will be used as examples employing endophenotypes derived from cognitive performance and MRI-based measures of neuroanatomy and neurophysiology. Once estimated, I show how the liability estimates can then be used to search for likely functional non-synonymous deleterious variants that influence addictive behavior by exploiting the extensive whole genome sequence information that we have for this study. The proposed estimation of disease liability using additional information from endophenotypes represents a novel way to explicitly incorporate endophenotypic information into analyses while retaining a direct focus on the genetic basis of the disease itself.

Advancements in Microarray Utility for Detection and Tracking of Foodborne Microbes in the Genomic Era

Outbreaks caused by foodborne microbes pose serious public health and food safety concerns worldwide. There is a great demand for rapid, sensitive and high-throughput methods to detect and track these pathogens in food, water and other environments. Recent advances in DNA genomic technology have enabled high-throughput analyses of strains by capturing total genomic content of strains and with concomitant comparative phylogenies. Microarrays are particularly adept for distilling large amounts of genomic DNA sequence information such as the gene(s) or genetic traits of hundreds of foodborne isolates in a single experiment. Hence, over the past two decades, microarray technology has advanced tremendously due to accessibility to thousands of complete and draft microbial genomes and this progress has led to the design and manufacturing of newer microarrays which can now identify gene sequence variations down to a single nucleotide polymorphism. DNA microarray remains a useful tool for rapid and refined genomic analysis of foodborne microbes. In this review, we will primarily focus our discussion on pathogen detection, serotype identification and tracking the genetic diversity and source of contamination of respective foodborne strains with our first-hand experience in using this technology.

Probing the Role of Cyclic di-GMP Signaling Systems in Disease Using Chinese Radish.

The determination of the genome sequences of pathogenic bacteria has facilitated functional analyses that aim to understand the molecular basis of virulence. In particular, genome sequence information of the pathogen Xanthomonas campestris pathovar campestris has allowed researchers to identify and functionally analyze the role of intracellular signaling involving cyclic di-GMP in black rot disease of crucifers. Here, we describe leaf clipping and spraying methods for testing the virulence of wild type and derived mutants of X. campestris in Chinese radish. These methods address different facets of the disease cycle, which requires the ability to survive epiphytically before entry into the plant and growth and systemic spread within the xylem.

Functional analysis of SH3 domain containing ring finger 2 during the myogenic differentiation of quail myoblast cells.

ObjectiveOwing to the public availability of complete genome sequences, including avian species, massive bioinformatics analyses may be conducted for computational gene prediction and the identification of gene regulatory networks through various informatics tools. However, to evaluate the biofunctional activity of a predicted target gene, in vivo and in vitro functional genomic analyses should be a prerequisite.MethodsDue to a lack of quail genomic sequence information, we first identified the partial genomic structure and sequences of the quail SH3 domain containing ring finger 2 (SH3RF2) gene. Subsequently, SH3RF2 was knocked out using clustered regularly interspaced short palindromic repeat/Cas9 technology and single cell-derived SH3RF2 mutant sublines were established to study the biofunctional activity of SH3RF2 in quail myoblast (QM7) cells during muscle differentiation.ResultsThrough a T7 endonuclease I assay and genotyping analysis, we established an SH3RF2 knockout (KO) QM7#4 subline with 61 and 155 nucleotide deletion mutations in SH3RF2. After the induction of myotube differentiation, the expression profiles were analyzed and compared between regular QM7 and SH3RF2 KO QM7#4 cells by global RNA sequencing and bioinformatics analysis.ConclusionWe did not detect any statistically significant role of SH3RF2 during myotube differentiation in QM7 myoblast cells. However, additional experiments are necessary to examine the biofunctional activity of SH3RF2 in cell proliferation and muscle growth.

Cornerstones of CRISPR-Cas in drug discovery and therapy.

The recent development of CRISPR-Cas systems as easily accessible and programmable tools for genome editing and regulation is spurring a revolution in biology. Paired with the rapid expansion of reference and personalized genomic sequence information, technologies based on CRISPR-Cas are enabling nearly unlimited genetic manipulation, even in previously difficult contexts, including human cells. Although much attention has focused on the potential of CRISPR-Cas to cure Mendelian diseases, the technology also holds promise to transform the development of therapies to treat complex heritable and somatic disorders. In this Review, we discuss how CRISPR-Cas can affect the next generation of drugs by accelerating the identification and validation of high-value targets, uncovering high-confidence biomarkers and developing differentiated breakthrough therapies. We focus on the promises, pitfalls and hurdles of this revolutionary gene-editing technology, discuss key aspects of different CRISPR-Cas screening platforms and offer our perspectives on the best practices in genome engineering.

Pathogenesis of microbial keratitis.

Microbial keratitis is a sight-threatening ocular infection caused by bacteria, fungi, and protist pathogens. Epithelial defects and injuries are key predisposing factors making the eye susceptible to corneal pathogens. Among bacterial pathogens, the most common agents responsible for keratitis include Staphylococcus aureus, Pseudomonas aeruginosa, Streptococcus pneumonia and Serratia species. Fungal agents of corneal infections include both filamentous as well as yeast, including Fusarium, Aspergillus, Phaeohyphomycetes, Curvularia, Paecilomyces, Scedosporium and Candida species, while in protists, Acanthamoeba spp. are responsible for causing ocular disease. Clinical features include redness, pain, tearing, blur vision and inflammation but symptoms vary depending on the causative agent. The underlying molecular mechanisms associated with microbial pathogenesis include virulence factors as well as the host factors that aid in the progression of keratitis, resulting in damage to the ocular tissue. The treatment therefore should focus not only on the elimination of the culprit but also on the neutralization of virulence factors to minimize the damage, in addition to repairing the damaged tissue. A complete understanding of the pathogenesis of microbial keratitis will lead to the rational development of therapeutic interventions. This is a timely review of our current understanding of the advances made in this field in a comprehensible manner. Coupled with the recently available genome sequence information and high throughput genomics technology, and the availability of innovative approaches, this will stimulate interest in this field.

A Guide to the PLAZA 3.0 Plant Comparative Genomic Database.

PLAZA 3.0 is an online resource for comparative genomics and offers a versatile platform to study gene functions and gene families or to analyze genome organization and evolution in the green plant lineage. Starting from genome sequence information for over 35 plant species, precomputed comparative genomic data sets cover homologous gene families, multiple sequence alignments, phylogenetic trees, and genomic colinearity information within and between species. Complementary functional data sets, a Workbench, and interactive visualization tools are available through a user-friendly web interface, making PLAZA an excellent starting point to translate sequence or omics data sets into biological knowledge. PLAZA is available at http://bioinformatics.psb.ugent.be/plaza/ .

Genomic Database Searching.

The availability of reference genome sequences for virtually all species under active research has revolutionized biology. Analyses of genomic variations in many organisms have provided insights into phenotypic traits, evolution and disease, and are transforming medicine. All genomic data from publicly funded projects are freely available in Internet-based databases, for download or searching via genome browsers such as Ensembl, Vega, NCBI's Map Viewer, and the UCSC Genome Browser. These online tools generate interactive graphical outputs of relevant chromosomal regions, showing genes, transcripts, and other genomic landmarks, and epigenetic features mapped by projects such as ENCODE.This chapter provides a broad overview of the major genomic databases and browsers, and describes various approaches and the latest resources for searching them. Methods are provided for identifying genomic locus and sequence information using gene names or codes, identifiers for DNA and RNA molecules and proteins; also from karyotype bands, chromosomal coordinates, sequences, motifs, and matrix-based patterns. Approaches are also described for batch retrieval of genomic information, performing more complex queries, and analyzing larger sets of experimental data, for example from next-generation sequencing projects.