Spirochaete genome identified in red abalone sample represents a novel genus Candidatus Haliotispira gen. nov. within the order Spirochaetales

A fully assembled spirochaete genome was identified as a contaminating scaffold in our red abalone (Haliotis rufescens) genome assembly. In this paper, we describe the analysis of this bacterial genome. The assembled spirochaete genome is 3.25 Mb in size with 48.5 mol% G+C content. The proteomes of 38 species were compared with the spirochaete genome and it was discovered to form an independent branch within the family Spirochaetaceae on the phylogenetic tree. The comparison of 16S rRNA sequences and average nucleotide identity scores between the spirochaete genome with known species of different families in Spirochaetia indicate that it is an unknown species. Further, the percentage of conserved proteins compared to neighbouring taxa confirm that it does not belong to a known genus within Spirochaetaceae. We propose the name Candidatus Haliotispira prima gen. nov., sp. nov. based on its taxonomic placement and origin. We also tested for the presence of this species in different species of abalone and found that it is also present in white abalone (Haliotis sorenseni). In addition, we highlight the need for better classification of taxa within the class Spirochaetia.

polishCLR: A Nextflow Workflow for Polishing PacBio CLR Genome Assemblies

Long-read sequencing has revolutionized genome assembly, yielding highly contiguous, chromosome-level contigs. However, assemblies from some third generation long read technologies, such as Pacific Biosciences (PacBio) continuous long reads (CLR), have a high error rate. Such errors can be corrected with short reads through a process called polishing. Although best practices for polishing non-model de novo genome assemblies were recently described by the Vertebrate Genome Project (VGP) Assembly community, there is a need for a publicly available, reproducible workflow that can be easily implemented and run on a conventional high performance computing environment. Here, we describe polishCLR (https://github.com/isugifNF/polishCLR), a reproducible Nextflow workflow that implements best practices for polishing assemblies made from CLR data. PolishCLR can be initiated from several input options that extend best practices to suboptimal cases. It also provides re-entry points throughout several key processes, including identifying duplicate haplotypes in purge_dups, allowing a break for scaffolding if data are available, and throughout multiple rounds of polishing and evaluation with Arrow and FreeBayes. PolishCLR is containerized and publicly available for the greater assembly community as a tool to complete assemblies from existing, error-prone long-read data.

Differential Transcriptomic Profiles Following Stimulation with Lipopolysaccharide in Intestinal Organoids from Dogs with Inflammatory Bowel Disease and Intestinal Mast Cell Tumor

Simple Summary

Lipopolysaccharide (LPS) derived from intestinal bacteria is linked to long-lasting inflammation that contributes to the development of intestinal cancer. While much research has been performed on the interplay between LPS and intestinal immune cells, little is known about how LPS influences intestinal epithelial cell structure and function. In this study, we investigated the effects of LPS on the proliferation and function of genes in intestinal organoids derived from dogs with gastrointestinal diseases, including inflammatory bowel disease (IBD) and intestinal mast cell tumor. The goal of this study was to evaluate how LPS affects signaling pathways in intestinal epithelial cells to influence development of a pro-tumor-like environment. Using an ex vivo model system, LPS incubation of organoids activated cancer-causing genes and accelerated the formation of IBD organoids derived from the small and large intestines. In brief, the crosstalk that occurs between the LPS/TLR4 signal transduction pathway and several different metabolic pathways, including primary bile acid biosynthesis and secretion, peroxisome, renin-angiotensin system, glutathione metabolism, and arachidonic acid pathways, may play a prominent role in the development of chronic intestinal inflammation and intestinal cancer.

Abstract

Lipopolysaccharide (LPS) is associated with chronic intestinal inflammation and promotes intestinal cancer progression in the gut. While the interplay between LPS and intestinal immune cells has been well-characterized, little is known about LPS and the intestinal epithelium interactions. In this study, we explored the differential effects of LPS on proliferation and the transcriptome in 3D enteroids/colonoids obtained from dogs with naturally occurring gastrointestinal (GI) diseases including inflammatory bowel disease (IBD) and intestinal mast cell tumor. The study objective was to analyze the LPS-induced modulation of signaling pathways involving the intestinal epithelia and contributing to colorectal cancer development in the context of an inflammatory (IBD) or a tumor microenvironment. While LPS incubation resulted in a pro-cancer gene expression pattern and stimulated proliferation of IBD enteroids and colonoids, downregulation of several cancer-associated genes such as Gpatch4, SLC7A1, ATP13A2, and TEX45 was also observed in tumor enteroids. Genes participating in porphyrin metabolism (CP), nucleocytoplasmic transport (EEF1A1), arachidonic acid, and glutathione metabolism (GPX1) exhibited a similar pattern of altered expression between IBD enteroids and IBD colonoids following LPS stimulation. In contrast, genes involved in anion transport, transcription and translation, apoptotic processes, and regulation of adaptive immune responses showed the opposite expression patterns between IBD enteroids and colonoids following LPS treatment. In brief, the crosstalk between LPS/TLR4 signal transduction pathway and several metabolic pathways such as primary bile acid biosynthesis and secretion, peroxisome, renin–angiotensin system, glutathione metabolism, and arachidonic acid pathways may be important in driving chronic intestinal inflammation and intestinal carcinogenesis.

Systematic in vitro specificity profiling reveals nicking defects in natural and engineered CRISPR-Cas9 variants

Cas9 is an RNA-guided endonuclease in the bacterial CRISPR-Cas immune system and a popular tool for genome editing. The commonly used Streptococcus pyogenes Cas9 (SpCas9) is relatively non-specific and prone to off-target genome editing. Other Cas9 orthologs and engineered variants of SpCas9 have been reported to be more specific. However, previous studies have focused on specificity of double-strand break (DSB) or indel formation, potentially overlooking alternative cleavage activities of these Cas9 variants. In this study, we employed in vitro cleavage assays of target libraries coupled with high-throughput sequencing to systematically compare cleavage activities and specificities of two natural Cas9 variants (SpCas9 and Staphylococcus aureus Cas9) and three engineered SpCas9 variants (SpCas9 HF1, HypaCas9 and HiFi Cas9). We observed that all Cas9s tested could cleave target sequences with up to five mismatches. However, the rate of cleavage of both on-target and off-target sequences varied based on target sequence and Cas9 variant. In addition, SaCas9 and engineered SpCas9 variants nick targets with multiple mismatches but have a defect in generating a DSB, while SpCas9 creates DSBs at these targets. Overall, these differences in cleavage rates and DSB formation may contribute to varied specificities observed in genome editing studies.

The thick aleurone1 Gene Encodes a NOT1 Subunit of the CCR4-NOT Complex and Regulates Cell Patterning in Endosperm

Maize (Zea mays) thick aleurone1 (thk1-R) mutants form multiple aleurone layers in the endosperm and have arrested embryogenesis. Prior studies suggest that thk1 functions downstream of defective kernel1 (dek1) in a regulatory pathway that controls aleurone cell fate and other endosperm traits. The original thk1-R mutant contained an ∼2-Mb multigene deletion, which precluded identification of the causal gene. Here, ethyl methanesulfonate mutagenesis produced additional alleles, and RNA sequencing from developing endosperm was used to identify a candidate gene based on differential expression compared with the wild-type progenitor. Gene editing confirmed the gene identity by producing mutant alleles that failed to complement existing thk1 mutants and that produced multiple-aleurone homozygous phenotypes. Thk1 encodes a homolog of NEGATIVE ON TATA-LESS1, a protein that acts as a scaffold for the CARBON CATABOLITE REPRESSION4-NEGATIVE ON TATA-LESS complex. This complex is highly conserved and essential in all eukaryotes for regulating a wide array of gene expression and cellular activities. Maize also harbors a duplicate locus, thick aleurone-like1, which likely accounts for the ability of thk1 mutants to form viable cells. Transcriptomic analysis indicated that THK1 regulates activities involving cell division, signaling, differentiation, and metabolism. Identification of thk1 provides an important new component of the DEK1 regulatory system that patterns cell fate in endosperm.

Detecting and correcting misclassified sequences in the large-scale public databases

Motivation

As the cost of sequencing decreases, the amount of data being deposited into public repositories is increasing rapidly. Public databases rely on the user to provide metadata for each submission that is prone to user error. Unfortunately, most public databases, such as non-redundant (NR), rely on user input and do not have methods for identifying errors in the provided metadata, leading to the potential for error propagation. Previous research on a small subset of the NR database analyzed misclassification based on sequence similarity. To the best of our knowledge, the amount of misclassification in the entire database has not been quantified. We propose a heuristic method to detect potentially misclassified taxonomic assignments in the NR database. We applied a curation technique and quality control to find the most probable taxonomic assignment. Our method incorporates provenance and frequency of each annotation from manually and computationally created databases and clustering information at 95% similarity.

Results

We found more than two million potentially taxonomically misclassified proteins in the NR database. Using simulated data, we show a high precision of 97% and a recall of 87% for detecting taxonomically misclassified proteins. The proposed approach and findings could also be applied to other databases.

Availability and implementation

Source code, dataset, documentation, Jupyter notebooks and Docker container are available at https://github.com/boalang/nr.

Supplementary information

Supplementary data are available at Bioinformatics online.

CRISPR-Cas12a has widespread off-target and dsDNA-nicking effects

Cas12a (Cpf1) is an RNA-guided endonuclease in the bacterial type V-A CRISPR-Cas anti-phage immune system that can be repurposed for genome editing. Cas12a can bind and cut dsDNA targets with high specificity in vivo, making it an ideal candidate for expanding the arsenal of enzymes used in precise genome editing. However, this reported high specificity contradicts Cas12a's natural role as an immune effector against rapidly evolving phages. Here, we employed high-throughput in vitro cleavage assays to determine and compare the native cleavage specificities and activities of three different natural Cas12a orthologs (FnCas12a, LbCas12a, and AsCas12a). Surprisingly, we observed pervasive sequence-specific nicking of randomized target libraries, with strong nicking of DNA sequences containing up to four mismatches in the Cas12a-targeted DNA-RNA hybrid sequences. We also found that these nicking and cleavage activities depend on mismatch type and position and vary with Cas12a ortholog and CRISPR RNA sequence. Our analysis further revealed robust nonspecific nicking of dsDNA when Cas12a is activated by binding to a target DNA. Together, our findings reveal that Cas12a has multiple nicking activities against dsDNA substrates and that these activities vary among different Cas12a orthologs.

Oral Treatments With Probiotics and Live Salmonella Vaccine Induce Unique Changes in Gut Neurochemicals and Microbiome in Chickens

Cross-talk between the gut microbiota and neurochemicals affects health and well-being of animals. However, little is known about this interaction in chickens despite their importance in food production. Probiotics and live Salmonella vaccines are microbial products commonly given orally to layer pullets to improve health and ensure food safety. This study's objective was to determine how these oral treatments, individually or in combination, would impact the gut environment of chickens. White Leghorn chicks were either non-treated (CON) or orally given probiotics (PRO), a recombinant attenuated Salmonella vaccine (RASV; VAX), or both (P+V). Birds were fed with probiotics daily beginning at 1-day-old and orally immunized with RASV at 4-days-old and boosted 2 weeks post-primary vaccination. At 5 weeks, ceca content, ceca tissues, and small intestinal scrapings (SISs) were collected from ten birds/group post-euthanasia for analyses. Catecholamine, but not serotonergic, metabolism was affected by treatments. Dopamine metabolism, indicated by L-DOPA and DOPAC levels, were increased in P+V birds versus CON and PRO birds. Based on 16S sequencing, beta diversity was more similar among vaccinated birds versus birds given probiotics, suggesting live Salmonella vaccination has a major selective pressure on microbial diversity. Abundances of Akkermansia muciniphila and Enterobacteriaceae positively correlated with levels of tyrosine and norepinephrine, respectively. Both enumeration and 16S sequencing, determined that PRO exhibited the greatest levels of Enterobacteriaceae in the ceca and feces, which was associated with greater IgA production against E. coli virulence factors as tested by ELISA. In summary, we demonstrate that using probiotics alone versus in combination with a live vaccine has major implications in catecholamine production and the microbiota of layer pullets. Additionally, unique correlations between changes in some neurochemicals and specific bacteria have been shown.

Shared data science infrastructure for genomics data

BACKGROUND

Creating a scalable computational infrastructure to analyze the wealth of information contained in data repositories is difficult due to significant barriers in organizing, extracting and analyzing relevant data. Shared data science infrastructures like Boag is needed to efficiently process and parse data contained in large data repositories. The main features of Boag are inspired from existing languages for data intensive computing and can easily integrate data from biological data repositories.

RESULTS

As a proof of concept, Boa for genomics, Boag, has been implemented to analyze RefSeq's 153,848 annotation (GFF) and assembly (FASTA) file metadata. Boag provides a massive improvement from existing solutions like Python and MongoDB, by utilizing a domain-specific language that uses Hadoop infrastructure for a smaller storage footprint that scales well and requires fewer lines of code. We execute scripts through Boag to answer questions about the genomes in RefSeq. We identify the largest and smallest genomes deposited, explore exon frequencies for assemblies after 2016, identify the most commonly used bacterial genome assembly program, and address how animal genome assemblies have improved since 2016. Boag databases provide a significant reduction in required storage of the raw data and a significant speed up in its ability to query large datasets due to automated parallelization and distribution of Hadoop infrastructure during computations.

CONCLUSIONS

In order to keep pace with our ability to produce biological data, innovative methods are required. The Shared Data Science Infrastructure, Boag, provides researchers a greater access to researchers to efficiently explore data in new ways. We demonstrate the potential of a the domain specific language Boag using the RefSeq database to explore how deposited genome assemblies and annotations are changing over time. This is a small example of how Boag could be used with large biological datasets.

The genome of the soybean cyst nematode (Heterodera glycines) reveals complex patterns of duplications involved in the evolution of parasitism genes

BACKGROUND

Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock.

RESULTS

Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons.

CONCLUSIONS

This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.

An Annotated Genome for Haliotis rufescens (Red Abalone) and Resequenced Green, Pink, Pinto, Black, and White Abalone Species

Abalone are one of the few marine taxa where aquaculture production dominates the global market as a result of increasing demand and declining natural stocks from overexploitation and disease. To better understand abalone biology, aid in conservation efforts for endangered abalone species, and gain insight into sustainable aquaculture, we created a draft genome of the red abalone (Haliotis rufescens). The approach to this genome draft included initial assembly using raw Illumina and PacBio sequencing data with MaSuRCA, before scaffolding using sequencing data generated from Chicago library preparations with HiRise2. This assembly approach resulted in 8,371 scaffolds and total length of 1.498 Gb; the N50 was 1.895 Mb, and the longest scaffold was 13.2 Mb. Gene models were predicted, using MAKER2, from RNA-Seq data and all related expressed sequence tags and proteins from NCBI; this resulted in 57,785 genes with an average length of 8,255 bp. In addition, single nucleotide polymorphisms were called on Illumina short-sequencing reads from five other eastern Pacific abalone species: the green (H. fulgens), pink (H. corrugata), pinto (H. kamtschatkana), black (H. cracherodii), and white (H. sorenseni) abalone. Phylogenetic relationships largely follow patterns detected by previous studies based on 1,784,991 high-quality single nucleotide polymorphisms. Among the six abalone species examined, the endangered white abalone appears to harbor the lowest levels of heterozygosity. This draft genome assembly and the sequencing data provide a foundation for genome-enabled aquaculture improvement for red abalone, and for genome-guided conservation efforts for the other five species and, in particular, for the endangered white and black abalone.

SCNBase: a genomics portal for the soybean cyst nematode (Heterodera glycines)

Soybean is an important worldwide crop, and farmers continue to experience significant yield loss due to the soybean cyst nematode (SCN), Heterodera glycines. This soil-borne roundworm parasite is rated the most important pathogen problem in soybean production. The infective nematodes enter into complex interactions with their host plant by inducing the development of specialized plant feeding cells that provide the parasites with nourishment. Addressing the SCN problem will require the development of genomic resources and a global collaboration of scientists to analyze and use these resources. SCNBase.org was designed as a collaborative hub for the SCN genome. All data and analyses are downloadable and can be analyzed with three integrated genomic tools: JBrowse, Feature Search and BLAST. At the time of this writing, a number of genomic and transcriptomic data sets are already available, with 43 JBrowse tracks and 21 category pages describing SCN genomic analyses on gene predictions, transcriptome and read alignments, effector-like genes, expansion and contraction of genomic repeats, orthology and synteny with related nematode species, Single Nucleotide Polymorphism (SNPs) from 15 SCN populations and novel splice sites. Standard functional gene annotations were supplemented with orthologous gene annotations using a comparison to nine related plant-parasitic nematodes, thereby enabling functional annotations for 85% of genes. These annotations led to a greater grasp on the SCN effectorome, which include over 3324 putative effector genes. By designing SCNBase as a hub, future research findings and genomic resources can easily be uploaded and made available for use by others with minimal needs for further curation. By providing these resources to nematode research community, scientists will be empowered to develop novel, more effective SCN management tools.

Vascular Endothelial Growth Factor A and Leptin Expression Associated with Ectopic Proliferation and Retinal Dysplasia in Zebrafish Optic Pathway Tumors

In the central nervous system injury induces cellular reprogramming and progenitor proliferation, but the molecular mechanisms that limit regeneration and prevent tumorigenesis are not completely understood. We previously described a zebrafish optic pathway tumor model in which transgenic Tg(flk1:RFP)is18/+ adults develop nonmalignant retinal tumors. Key pathways driving injury-induced glial reprogramming and regeneration contributed to tumor formation. In this study, we examine a time course of proliferation and present new analyses of the Tg(flk1:RFP)is18/+ dysplastic retina and tumor transcriptomes. Retinal dysplasia was first detected in 3-month-old adults, but was not limited to a specific stem cell or progenitor niche. Pathway analyses suggested a decrease in cellular respiration and increased expression of components of Hif1-α, VEGF, mTOR, NFκβ, and multiple interleukin pathways are associated with early retinal dysplasia. Hif-α targets VEGFA (vegfab) and Leptin (lepb) were both highly upregulated in dysplastic retina; however, each showed distinct expression patterns in neurons and glia, respectively. Phospho-S6 immunolabeling indicated that mTOR signaling is activated in multiple cell populations in wild-type retina and in the dysplastic retina and advanced tumor. Our results suggest that multiple pathways may contribute to the continuous proliferation of retinal progenitors and tumor growth in this optic pathway tumor model. Further investigation of these signaling pathways may yield insight into potential mechanisms to control the proliferative response during regeneration in the nervous system.

Genome, transcriptome and methylome sequencing of a primitively eusocial wasp reveal a greatly reduced DNA methylation system in a social insect

Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect Polistes dominula, accompanied by analysis of caste-related transcriptome and methylome sequence data for adult queens and workers. Polistes dominula possesses a fairly typical hymenopteran genome, but shows very low genomewide GC content and some evidence of reduced genome size. We found numerous caste-related differences in gene expression, with evidence that both conserved and novel genes are related to caste differences. Most strikingly, these -omics data reveal a major reduction in one of the major epigenetic mechanisms that has been previously suggested to be important for caste differences in social insects: DNA methylation. Along with a conspicuous loss of a key gene associated with environmentally responsive DNA methylation (the de novo DNA methyltransferase Dnmt3), these wasps have greatly reduced genomewide methylation to almost zero. In addition to providing a valuable resource for comparative analysis of social insect evolution, our integrative -omics data for this important behavioural and evolutionary model system call into question the general importance of DNA methylation in caste differences and evolution in social insects.

Transcript profiles in longissimus dorsi muscle and subcutaneous adipose tissue: a comparison of pigs with different postweaning growth rates

Although most pigs recover rapidly from stresses associated with the transition of weaning, a portion of the population lags behind their contemporaries in growth performance. The underlying biological and molecular mechanisms involved in postweaning differences in growth performance are poorly understood. The objective of this experiment was to use transcriptional profiling of skeletal muscle and adipose tissue to develop a better understanding of the metabolic basis for poor weaned-pig transition. A total of 1,054 pigs was reared in commercial conditions and weighed at birth, weaning, and 3 wk postweaning. Transition ADG (tADG) was calculated as the ADG for the 3-wk period postweaning. Nine pigs from both the lowest 10th percentile (low tADG) and the 60th to 70th percentile (high tADG) were harvested at 3 wk postweaning. Differential expression analysis was conducted in longissimus dorsi muscle (LM) and subcutaneous adipose tissue using RNA-Seq methodology. In LM, 768 transcripts were differentially expressed (DE), 327 with higher expression in low tADG and 441 with higher expression in high tADG pigs (q < 0.10). Expression patterns measured in LM by RNA-Seq were verified in 30 of 32 transcripts using quantitative PCR. No DE transcripts were identified in adipose tissue. To identify biological functions potentially underlying the effects of tADG on skeletal muscle metabolism and physiology, functional annotation analysis of the DE transcripts was conducted using DAVID and Pathway Studio analytic tools. The group of DE genes with lower expression in LM of low tADG pigs was enriched in genes with functions related to muscle contraction, glucose metabolism, cytoskeleton organization, muscle development, and response to hormone stimulus (enrichment score > 1.3). The list of DE genes with higher expression in low tADG LM was enriched in genes with functions related to protein catabolism (enrichment score > 1.3). Analysis of known gene-gene interactions identified possible regulators of these differences in gene expression in LM of high and low tADG pigs; these include forkhead box O1 (FOXO1), growth hormone (GH1), and the glucocorticoid receptor (NR3C1). Differences in gene expression between poor transitioning pigs and their contemporaries indicate a shift to decreased protein synthesis, increased protein degradation, and reduced glucose metabolism in the LM of low tADG pigs.

Transcriptome profiling of soybean (Glycine max) roots challenged with pathogenic and non-pathogenic isolates of Fusarium oxysporum

BACKGROUND

Fusarium oxysporum is one of the most common fungal pathogens causing soybean root rot and seedling blight in U.S.A. In a recent study, significant variation in aggressiveness was observed among isolates of F. oxysporum collected from roots in Iowa, ranging from highly pathogenic to weakly or non-pathogenic isolates.

RESULTS

We used RNA-seq analysis to investigate the molecular aspects of the interactions of a partially resistant soybean genotype with non-pathogenic/pathogenic isolates of F. oxysporum at 72 and 96 h post inoculation (hpi). Markedly different gene expression profiles were observed in response to the two isolates. A peak of highly differentially expressed genes (HDEGs) was triggered at 72 hpi in soybean roots and the number of HDEGs was about eight times higher in response to the pathogenic isolate compared to the non-pathogenic one (1,659 vs. 203 HDEGs, respectively). Furthermore, the magnitude of induction was much greater in response to the pathogenic isolate. This response included a stronger activation of defense-related genes, transcription factors, and genes involved in ethylene biosynthesis, secondary and sugar metabolism.

CONCLUSIONS

The obtained data provide an important insight into the transcriptional responses of soybean-F. oxysporum interactions and illustrate the more drastic changes in the host transcriptome in response to the pathogenic isolate. These results may be useful in the developing new methods of broadening resistance of soybean to F. oxysporum, including the over-expression of key soybean genes.

Microbiome analyses of pacific white shrimp (Litopenaeus vannamei) collected from disparate geographical locations

In this study, the tail muscle microbiota of pacific white shrimp (Litopenaeus vannamei) sourced from five countries across Central and South America and Southeast Asia were determined and compared. The genomic DNA was sequenced at around 10 × coverage for each geographical location and was assembled de novo for comparative analysis. The assembled sequences for all the lines were classified based on their similarity to the sequences in the public database. We found that there is high correlation among the microbiota of shrimp from disparate regions, as well as the presence of some DNA from bacteria known to cause food poisoning in humans. Sequencing data has been deposited at NCBI-SRA database and can be found under the BioProject ID PRJNA282154.

CRISPR interference and priming varies with individual spacer sequences

CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR associated) systems allow bacteria to adapt to infection by acquiring 'spacer' sequences from invader DNA into genomic CRISPR loci. Cas proteins use RNAs derived from these loci to target cognate sequences for destruction through CRISPR interference. Mutations in the protospacer adjacent motif (PAM) and seed regions block interference but promote rapid 'primed' adaptation. Here, we use multiple spacer sequences to reexamine the PAM and seed sequence requirements for interference and priming in the Escherichia coli Type I-E CRISPR-Cas system. Surprisingly, CRISPR interference is far more tolerant of mutations in the seed and the PAM than previously reported, and this mutational tolerance, as well as priming activity, is highly dependent on spacer sequence. We identify a large number of functional PAMs that can promote interference, priming or both activities, depending on the associated spacer sequence. Functional PAMs are preferentially acquired during unprimed 'naïve' adaptation, leading to a rapid priming response following infection. Our results provide numerous insights into the importance of both spacer and target sequences for interference and priming, and reveal that priming is a major pathway for adaptation during initial infection.

Transcriptional analysis of phloem-associated cells of potato

Background

Numerous signal molecules, including proteins and mRNAs, are transported through the architecture of plants via the vascular system. As the connection between leaves and other organs, the petiole and stem are especially important in their transport function, which is carried out by the phloem and xylem, especially by the sieve elements in the phloem system. The phloem is an important conduit for transporting photosynthate and signal molecules like metabolites, proteins, small RNAs, and full-length mRNAs. Phloem sap has been used as an unadulterated source to profile phloem proteins and RNAs, but unfortunately, pure phloem sap cannot be obtained in most plant species.

Results

Here we make use of laser capture microdissection (LCM) and RNA-seq for an in-depth transcriptional profile of phloem-associated cells of both petioles and stems of potato. To expedite our analysis, we have taken advantage of the potato genome that has recently been fully sequenced and annotated. Out of the 27 k transcripts assembled that we identified, approximately 15 k were present in phloem-associated cells of petiole and stem with greater than ten reads. Among these genes, roughly 10 k are affected by photoperiod. Several RNAs from this day length-regulated group are also abundant in phloem cells of petioles and encode for proteins involved in signaling or transcriptional control. Approximately 22 % of the transcripts in phloem cells contained at least one binding motif for Pumilio, Nova, or polypyrimidine tract-binding proteins in their downstream sequences. Highlighting the predominance of binding processes identified in the gene ontology analysis of active genes from phloem cells, 78 % of the 464 RNA-binding proteins present in the potato genome were detected in our phloem transcriptome.

Conclusions

As a reasonable alternative when phloem sap collection is not possible, LCM can be used to isolate RNA from specific cell types, and along with RNA-seq, provides practical access to expression profiles of phloem tissue. The combination of these techniques provides a useful approach to the study of phloem and a comprehensive picture of the mechanisms associated with long-distance signaling. The data presented here provide valuable insights into potentially novel phloem-mobile mRNAs and phloem-associated RNA-binding proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1844-2) contains supplementary material, which is available to authorized users.

A comparative analysis of methylome profiles of Campylobacter jejuni sheep abortion isolate and gastroenteric strains using PacBio data

Campylobacter jejuni is a leading cause of human gastrointestinal disease and small ruminant abortions in the United States. The recent emergence of a highly virulent, tetracycline-resistant C. jejuni subsp. jejuni sheep abortion clone (clone SA) in the United States, and that strain's association with human disease, has resulted in a heightened awareness of the zoonotic potential of this organism. Pacific Biosciences' Single Molecule, Real-Time sequencing technology was used to explore the variation in the genome-wide methylation patterns of the abortifacient clone SA (IA3902) and phenotypically distinct gastrointestinal-specific C. jejuni strains (NCTC 11168 and 81-176). Several notable differences were discovered that distinguished the methylome of IA3902 from that of 11168 and 81-176: identification of motifs novel to IA3902, genome-specific hypo- and hypermethylated regions, strain level variability in genes methylated, and differences in the types of methylation motifs present in each strain. These observations suggest a possible role of methylation in the contrasting disease presentations of these three C. jejuni strains. In addition, the methylation profiles between IA3902 and a luxS mutant were explored to determine if variations in methylation patterns could be identified that might explain the role of LuxS-dependent methyl recycling in IA3902 abortifacient potential.

Molecular and cellular characterization of a zebrafish optic pathway tumor line implicates glia-derived progenitors in tumorigenesis

In this study we describe the molecular and cellular characterization of a zebrafish mutant that develops tumors in the optic pathway. Heterozygous Tg(flk1:RFP)is18 transgenic adults develop tumors of the retina, optic nerve and optic tract. Molecular and genetic mapping demonstrate the tumor phenotype is linked to a high copy number transgene array integrated in the lincRNA gene lincRNAis18/Zv9_00007276 on chromosome 3. TALENs were used to isolate a 147kb deletion allele that removes exons 2–5 of the lincRNAis18 gene. Deletion allele homozygotes are viable and do not develop tumors, indicating loss of function of the lincRNAis18 locus is not the trigger for tumor onset. Optic pathway tumors in the Tg(flk1:RFP)is18 mutant occur with a penetrance of 80–100% by 1 year of age. The retinal tumors are highly vascularized and composed of rosettes of various sizes embedded in a fibrous matrix. Immunohistochemical analysis showed increased expression of the glial markers GFAP and BLBP throughout retinal tumors and in dysplastic optic nerve. We performed transcriptome analysis of pre-tumorous retina and retinal tumor tissue and found changes in gene expression signatures of radial glia and astrocytes (slc1a3), activated glia (atf3, blbp, apoeb), proliferating neural progenitors (foxd3, nestin, cdh2, her9/hes1), and glioma markers (S100β, vim). The transcriptome also revealed activation of cAMP, Stat3 and Wnt signal transduction pathways. qRT-PCR confirmed >10-fold overexpression of the Wnt pathway components hbegfa, ascl1a, and insm1a. Together the data indicate Müller glia and/or astrocyte-derived progenitors could contribute to the zebrafish Tg(flk1:RFP)is18 optic pathway tumors.

Identification of candidate genes underlying an iron efficiency quantitative trait locus in soybean

Prevalent on calcareous soils in the United States and abroad, iron deficiency is among the most common and severe nutritional stresses in plants. In soybean (Glycine max) commercial plantings, the identification and use of iron-efficient genotypes has proven to be the best form of managing this soil-related plant stress. Previous studies conducted in soybean identified a significant iron efficiency quantitative trait locus (QTL) explaining more than 70% of the phenotypic variation for the trait. In this research, we identified candidate genes underlying this QTL through molecular breeding, mapping, and transcriptome sequencing. Introgression mapping was performed using two related near-isogenic lines in which a region located on soybean chromosome 3 required for iron efficiency was identified. The region corresponds to the previously reported iron efficiency QTL. The location was further confirmed through QTL mapping conducted in this study. Transcriptome sequencing and quantitative real-time-polymerase chain reaction identified two genes encoding transcription factors within the region that were significantly induced in soybean roots under iron stress. The two induced transcription factors were identified as homologs of the subgroup lb basic helix-loop-helix (bHLH) genes that are known to regulate the strategy I response in Arabidopsis (Arabidopsis thaliana). Resequencing of these differentially expressed genes unveiled a significant deletion within a predicted dimerization domain. We hypothesize that this deletion disrupts the Fe-DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT)/bHLH heterodimer that has been shown to induce known iron acquisition genes.

The Medicago genome provides insight into the evolution of rhizobial symbioses

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

Changes in twelve homoeologous genomic regions in soybean following three rounds of polyploidy

With the advent of high-throughput sequencing, the availability of genomic sequence for comparative genomics is increasing exponentially. Numerous completed plant genome sequences enable characterization of patterns of the retention and evolution of genes within gene families due to multiple polyploidy events, gene loss and fractionation, and differential evolutionary pressures over time and across different gene families. In this report, we trace the changes that have occurred in 12 surviving homoeologous genomic regions from three rounds of polyploidy that contributed to the current Glycine max genome: a genome triplication before the origin of the rosids (~130 to 240 million years ago), a genome duplication early in the legumes (~58 million years ago), and a duplication in the Glycine lineage (~13 million years ago). Patterns of gene retention following the genome triplication event generally support predictions of the Gene Balance Hypothesis. Finally, we find that genes in networks with a high level of connectivity are more strongly conserved than those with low connectivity and that the enrichment of these highly connected genes in the 12 highly conserved homoeologous segments may in part explain their retention over more than 100 million years and repeated polyploidy events.

Gene expression patterns are correlated with genomic and genic structure in soybean

Studies have indicated that exon and intron size and intergenic distance are correlated with gene expression levels and expression breadth. Previous reports on these correlations in plants and animals have been conflicting. In this study, next-generation sequence data, which has been shown to be more sensitive than previous expression profiling technologies, were generated and analyzed from 14 tissues. Our results revealed a novel dichotomy. At the low expression level, an increase in expression breadth correlated with an increase in transcript size because of an increase in the number of exons and introns. No significant changes in intron or exon sizes were noted. Conversely, genes expressed at the intermediate to high expression levels displayed a decrease in transcript size as their expression breadth increased. This was due to smaller exons, with no significant change in the number of exons. Taking advantage of the known gene space of soybean, we evaluated the positioning of genes and found significant clustering of similarly expressed genes. Identifying the correlations between the physical parameters of individual genes could lead to uncovering the role of regulation owing to nucleotide composition, which might have potential impacts in discerning the role of the noncoding regions.

An integrative approach to genomic introgression mapping

Near-isogenic lines (NILs) are valuable genetic resources for many crop species, including soybean (Glycine max). The development of new molecular platforms promises to accelerate the mapping of genetic introgressions in these materials. Here, we compare some existing and emerging methodologies for genetic introgression mapping: single-feature polymorphism analysis, Illumina GoldenGate single nucleotide polymorphism (SNP) genotyping, and de novo SNP discovery via RNA-Seq analysis of next-generation sequence data. We used these methods to map the introgressed regions in an iron-inefficient soybean NIL and found that the three mapping approaches are complementary when utilized in combination. The comparative RNA-Seq approach offers several additional advantages, including the greatest mapping resolution, marker depth, and de novo marker utility for downstream fine-mapping analysis. We applied the comparative RNA-Seq method to map genetic introgressions in an additional pair of NILs exhibiting differential seed protein content. Furthermore, we attempted to optimize the comparative RNA-Seq approach by assessing the impact of sequence depth, SNP identification methodology, and post hoc analyses on SNP discovery rates. We conclude that the comparative RNA-Seq approach can be optimized with sufficient sampling and by utilizing a post hoc correction accounting for gene density variation that controls for false discoveries.

RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome

BACKGROUND

Next generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.

RESULTS

The RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.

CONCLUSIONS

This RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.