WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update

Immune response gene coexpression network analysis of Arachis hypogaea infected with Aspergillus flavus

Aspergillus flavus (A. flavus) infection and aflatoxin contamination is a major bottleneck for peanut cultivation and value chain industry. In this study, a transcriptomic network study was conducted by retrieving publically available RNA-seq datasets of resistant and susceptible peanut varieties infected by A. flavus separately to understand the peanut defense mechanism against A. flavus. The gene expression analysis revealed differentially expressed genes (DEGs) in response to the different levels of infection and coexpression network of DEGs deciphered hub genes involved in the immune process in resistant and susceptible varieties. The interplay of resistance conferring genes and cell wall related genes was observed through functional enrichment analysis in response to pathogen infection and identified few key genes such as Protein P21, R genes, Pattern Recognition Receptor genes, Pectinesterases, Laccase and Thaumatin-like protein 1b as candidate genes in imparting immune response against A. flavus.

Extraction and identification of biologically important proteins from the medicinal plant God's crown (Phaleria macrocarpa)

The fruit and leaf of God's crown (Phaleria macrocarpa) have been traditionally used to treat a wide variety of diseases. However, the proteins of this tropical plant are still heavily understudied. Three protein extraction methods; phenol (Phe), trichloroacetic acid (TCA)-acetone-phenol (TCA-A-Phe), and ultrasonic (Ult) were compared on the fruit and leaf of P. macrocarpa. The Phe extraction method showed the highest percentage of recovered protein after the resolubilization process for both leaf (12.24%) and fruit (30.41%) based on protein yields of the leaf (6.15 mg/g) and fruit (36.98 mg/g). Phe and TCA-A-Phe extraction methods gave well-resolved bands over a wide range of molecular weights through sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Following liquid chromatography-tandem mass spectrometry analysis, proteins identified through the Phe extraction method were 30%-35% enzymatic proteins, including oxidoreductases, transferases, hydrolases, lyases, isomerases, and ligases that possess various biological functions. PRACTICAL APPLICATIONS: Every part of God's crown plant is traditionally consumed to treat various illnesses. While plant's benefits are well known and have led to a plethora of health products, the proteome remains mostly unknown. This study compares three protein extraction methods for the leaf and fruit of P. macrocarpa and identifies their proteins thru LC-MS/MS coupled with PEAKS. These method comparisons can be a guide for works on other plants as well. In addition, the proteomics data from this study may shed light on the functional properties of these plant parts and their products.

In silico identification of excretory/secretory proteins and drug targets in monogenean parasites

The Excretory/Secretory (ES) proteins of parasites are involved in invasion and colonization of their hosts. In addition, since ES proteins circulate in the extracellular space, they can be more accessible to drugs than other proteins, which makes ES proteins optimal targets for the development of new and better pharmacological strategies. Monogeneans are a group of parasitic Platyhelminthes that includes some pathogenic species problematic for finfish aquaculture. In the present study, 8297 putative ES proteins from four monogenean species which genomic resources are publicly available were identified and functionally annotated by bioinformatic tools. Additionally, for comparative purposes, ES proteins in other parasitic and free-living platyhelminths were identified. Based on data from the monogenean Gyrodactylus salaris, 15 ES proteins are considered potential drug targets. One of them showed homology to 10 cathepsins with known 3D structure. A docking molecular analysis uncovered that the anthelmintic emodepside shows good affinity to these cathepsins suggesting that emodepside can be experimentally tested as a monogenean's cathepsin inhibitor.

LTR-retrotransposon dynamics in common fig (Ficus carica L.) genome

BACKGROUND

Long Terminal Repeat retrotransposons (LTR-REs) are repetitive DNA sequences that constitute a large part of the genome. The improvement of sequencing technologies and sequence assembling strategies has achieved genome sequences with much greater reliability than those of the past, especially in relation to repetitive DNA sequences.

RESULTS

In this study, we analysed the genome of Ficus carica L., obtained using third generation sequencing technologies and recently released, to characterise the complete complement of full-length LTR-REs to study their dynamics during fig genome evolution. A total of 1867 full-length elements were identified. Those belonging to the Gypsy superfamily were the most abundant; among these, the Chromovirus/Tekay lineage was the most represented. For the Copia superfamily, Ale was the most abundant lineage. Measuring the estimated insertion time of each element showed that, on average, Ivana and Chromovirus/Tekay were the youngest lineages of Copia and Gypsy superfamilies, respectively. Most elements were inactive in transcription, both constitutively and in leaves of plants exposed to an abiotic stress, except for some elements, mostly belonging to the Copia/Ale lineage. A relationship between the inactivity of an element and inactivity of genes lying in close proximity to it was established.

CONCLUSIONS

The data reported in this study provide one of the first sets of information on the genomic dynamics related to LTR-REs in a plant species with highly reliable genome sequence. Fig LTR-REs are highly heterogeneous in abundance and estimated insertion time, and only a few elements are transcriptionally active. In general, the data suggested a direct relationship between estimated insertion time and abundance of an element and an inverse relationship between insertion time (or abundance) and transcription, at least for Copia LTR-REs.

The recombination landscape and multiple QTL mapping in a Solanum tuberosum cv. 'Atlantic'-derived F1 population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

Underpinning the molecular programming attributing heat stress associated thermotolerance in tea (Camellia sinensis (L.) O. Kuntze)

The most daunting issue of global climate change is the deleterious impact of extreme temperatures on tea productivity and quality, which has resulted in a quest among researchers and growers. The current study aims to unravel molecular programming underpinning thermotolerance by characterizing heat tolerance and sensitivity response in 20 tea cultivars. The significantly higher negative influence of heat stress was recorded in a sensitive cultivar with reduced water retention (47%), chlorophyll content (33.79%), oxidation potential (32.48%), and increase in membrane damage (76.4%). Transcriptional profiling of most tolerant and sensitive cultivars identified 78 differentially expressed unigenes with chaperon domains, including low and high molecular weight heat shock protein (HSP) and heat shock transcription factors (HSFs) involved in heat shock response (HSR). Further, predicted transcriptional interactome network revealed their key role in thermotolerance via well-co-ordinated transcriptional regulation of aquaporins, starch metabolism, chlorophyll biosynthesis, calcium, and ethylene mediated plant signaling system. The study identified the key role of HSPs (CsHSP90) in regulating HSR in tea, wherein, structure-based molecular docking revealed the inhibitory role of geldanamycin (GDA) on CsHSP90 by blocking ATP binding site at N-terminal domain of predicted structure. Subsequently, GDA mediated leaf disc inhibitor assay further affirmed enhanced HSR with higher expression of CsHSP17.6, CsHSP70, HSP101, and CsHSFA2 genes in tea. Through the current study, efforts were made to extrapolate a deeper understanding of chaperons mediated regulation of HSR attributing thermotolerance in tea.

Streptomyces brasiliscabiei, a new species causing potato scab in south Brazil

This study aimed to characterize six Streptomyces strains associated with potato scab in south Brazil through polyphasic taxonomy involving morphology, pathogenicity and genetic features. These strains were compared with other potato-scab Streptomyces species mainly S. europaeiscabiei, S. scabiei and S. stelliscabiei. South-Brazilian Streptomyces strains were morphologically distinct from the type strains of S. scabiei (CFBP 4517^T) and their genomospecies S. europaeiscabiei (CFBP 4497 ^T) and S. stelliscabiei (CFBP 4521^T), producing a brown substrate mycelium with red borders and cream-grey brown aerial spores. Red-brown diffusible pigment on YME was also observed. The carbon sources L-Arabinose, D-Fructose, D-Glucose, D-Mannitol, meso-Inositol, Raffinose, Rhamnose, Sucrose, D-Xylose were tested for these strains. All strains were pathogenic causing symptoms of necrosis on radish and several potato cultivars commonly used in potato growing areas in Brazil. In greenhouse conditions, the strains caused scab disease and produced deep-pitted lesions covering large areas of the tuber. These results were correlated with presence of pathogenicity marker genes (txtAB, tomA or nec1) detected by PCR amplifications. In both phylogenetic analyses, 16S rRNA and MLSA, Streptomyces sp. Brazilian strains were closely related to S. europaeiscabiei, S. scabiei and S. stelliscabiei species, but they were allocated in separated branches supported by high bootstrap values and/or with low sequence similarity values. Sequencing of whole genome showed an 10,846,379 bp linear chromosome with high GC content (71.3%) consisting of 9179 putative genes, 3 rRNAs, 89 tRNAs and 1 CRISPRS. The molecular data, including genomic features, associated with morphological, biochemical and pathogenic characteristics warrant that the six Streptomyces Brazilian strains represent a new species associated with potato scab in Brazil, which would be named Streptomyces brasiliscabiei with IBSBF 2867^T as the type strain.

Immunoproteomic analysis of fish ectoparasite, Argulus siamensis antigens

AIM

An immunoproteomic approach was followed to identify immunoreactive antigens of fish ectoparasite, Argulus siamensis with rohu (Labeo rohita) immune sera for screening of potential vaccine candidates.

MATERIALS AND RESULTS

The whole adult Argulus antigen was run in 2D electrophoresis with IEF in 7 cm IPG strips of pH 4-7 and SDS-PAGE with 12% acrylamide concentration. Two parallel gels were run; one was stained with silver stain, and the other was Western blotted to nitrocellulose paper (NCP) and reacted with rohu anti-A siamensis sera. Fourteen protein spots corresponding to the spots developed in NCP were picked from the silver-stained gel and subjected to mass spectrometry in MALDI-TOF/TOF. The MS/MS spectra were analysed in MASCOT software with taxonomy 'other metazoa' and the proteins identified based on similarity with the proteins from heterologous species. The gene ontology analysis revealed a majority of proteins being involved in binding activity in 'molecular function' and belonging to metabolic processes in 'biologic process' categories. The possibility of these proteins as vaccine candidates against A siamensis is discussed in the paper.

CONCLUSION

Three of the identified proteins namely, bromodomain-containing protein, anaphase-promoting complex subunit 5 and elongation factor-2 could possibly serve as vaccine candidates against argulosis in carps.

Root Exudates Alter the Expression of Diverse Metabolic, Transport, Regulatory, and Stress Response Genes in Rhizosphere Pseudomonas

Plants live in association with microorganisms that positively influence plant development, vigor, and fitness in response to pathogens and abiotic stressors. The bulk of the plant microbiome is concentrated belowground at the plant root-soil interface. Plant roots secrete carbon-rich rhizodeposits containing primary and secondary low molecular weight metabolites, lysates, and mucilages. These exudates provide nutrients for soil microorganisms and modulate their affinity to host plants, but molecular details of this process are largely unresolved. We addressed this gap by focusing on the molecular dialog between eight well-characterized beneficial strains of the Pseudomonas fluorescens group and Brachypodium distachyon, a model for economically important food, feed, forage, and biomass crops of the grass family. We collected and analyzed root exudates of B. distachyon and demonstrated the presence of multiple carbohydrates, amino acids, organic acids, and phenolic compounds. The subsequent screening of bacteria by Biolog Phenotype MicroArrays revealed that many of these metabolites provide carbon and energy for the Pseudomonas strains. RNA-seq profiling of bacterial cultures amended with root exudates revealed changes in the expression of genes encoding numerous catabolic and anabolic enzymes, transporters, transcriptional regulators, stress response, and conserved hypothetical proteins. Almost half of the differentially expressed genes mapped to the variable part of the strains’ pangenome, reflecting the importance of the variable gene content in the adaptation of P. fluorescens to the rhizosphere lifestyle. Our results collectively reveal the diversity of cellular pathways and physiological responses underlying the establishment of mutualistic interactions between these beneficial rhizobacteria and their plant hosts.

Genome-wide discovery of G-quadruplexes in barley

G-quadruplexes (G4s) are four-stranded nucleic acid structures with closely spaced guanine bases forming square planar G-quartets. Aberrant formation of G4 structures has been associated with genomic instability. However, most plant species are lacking comprehensive studies of G4 motifs. In this study, genome-wide identification of G4 motifs in barley was performed, followed by a comparison of genomic distribution and molecular functions to other monocot species, such as wheat, maize, and rice. Similar to the reports on human and some plants like wheat, G4 motifs peaked around the 5′ untranslated region (5′ UTR), the first coding domain sequence, and the first intron start sites on antisense strands. Our comparative analyses in human, Arabidopsis, maize, rice, and sorghum demonstrated that the peak points could be erroneously merged into a single peak when large window sizes are used. We also showed that the G4 distributions around genic regions are relatively similar in the species studied, except in the case of Arabidopsis. G4 containing genes in monocots showed conserved molecular functions for transcription initiation and hydrolase activity. Additionally, we provided examples of imperfect G4 motifs.

The protective effect and mechanism of epidermal growth factor on necrotizing enterocolitis in a neonatal rat model

BACKGROUND

Necrotizing enterocolitis (NEC) is the most common acquired gastrointestinal emergency in premature infants. This study aimed to investigate the protective effect and mechanism of epidermal growth factor (EGF) on NEC in a neonatal rat model.

METHODS

We randomly divided 50 newborn SD rats into a control group, NEC group, NEC + 50 ng/mL EGF group, NEC + 500 ng/mL EGF group, and NEC + 1,000 ng/mL EGF group, with 10 cases in each group. The appearance of intestinal tissue, physiological status score, inflammatory factor level, HE staining, and pathological score were used to evaluate the protective effect. A one cm tissue sample from the proximal ileum of the ileocecal area of five rats from the NEC group and the group that showed a significant protective effect were extracted for transcriptome sequencing.

RESULTS

The levels of IL-1β and IL-6 in the intestinal mucosa in the NEC + 500 ng/mL EGF group were significantly lower than those in the NEC + 1,000 ng/mL EGF group (P<0.05). Transcriptome sequencing suggested that EGF effects the intestinal barrier, apoptosis, and inflammation of the NEC intestine.

CONCLUSIONS

We conclude that the oral administration of 500 ng/mL EGF effectively inhibits intestinal inflammation in NEC neonatal rat models, thereby affecting the barrier function of the intestinal tract.

Distinct miRNA Signatures and Networks Discern Fetal from Adult Erythroid Differentiation and Primary from Immortalized Erythroid Cells

MicroRNAs (miRNAs) are small non-coding RNAs crucial for post-transcriptional and translational regulation of cellular and developmental pathways. The study of miRNAs in erythropoiesis elucidates underlying regulatory mechanisms and facilitates related diagnostic and therapy development. Here, we used DNA Nanoball (DNB) small RNA sequencing to comprehensively characterize miRNAs in human erythroid cell cultures. Based on primary human peripheral-blood-derived CD34+ (hCD34+) cells and two influential erythroid cell lines with adult and fetal hemoglobin expression patterns, HUDEP-2 and HUDEP-1, respectively, our study links differential miRNA expression to erythroid differentiation, cell type, and hemoglobin expression profile. Sequencing results validated by reverse-transcription quantitative PCR (RT-qPCR) of selected miRNAs indicate shared differentiation signatures in primary and immortalized cells, characterized by reduced overall miRNA expression and reciprocal expression increases for individual lineage-specific miRNAs in late-stage erythropoiesis. Despite the high similarity of same-stage hCD34+ and HUDEP-2 cells, differential expression of several miRNAs highlighted informative discrepancies between both cell types. Moreover, a comparison between HUDEP-2 and HUDEP-1 cells displayed changes in miRNAs, transcription factors (TFs), target genes, and pathways associated with globin switching. In resulting TF-miRNA co-regulatory networks, major therapeutically relevant regulators of globin expression were targeted by many co-expressed miRNAs, outlining intricate combinatorial miRNA regulation of globin expression in erythroid cells.

Zinc Finger-Homeodomain Transcriptional Factors (ZF-HDs) in Wheat (Triticum aestivum L.): Identification, Evolution, Expression Analysis and Response to Abiotic Stresses

Zinc finger-homeodomain transcriptional factors (ZF-HDs), a kind of plant-specific transcription factor, play important roles in plant growth, development and various stress responses. In this study, the genome-wide analysis of the ZF-HD gene family was performed in wheat. A total of 37 TaZF-HD genes were identified in T. aestivum and classified into six groups. The results of a synteny analysis showed that gene replication events contributed to the expansion of the TaZF-HD gene family. The TaZF-HD paralogous gene pairs with similar chromosomal locations in different subgenomes had similar expression patterns. TaZF-HDs were highly induced under PEG (polyethylene glycol), NaCl and cold stress but not induced under heat stress. Gene ontology (GO) annotation and protein-protein interactions suggested that TaZF-HD proteins may participate in various biological processes of plants. These results increase our understanding of ZF-HD genes and provide robust candidate genes for future functional investigations aimed at crop improvement.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

RNA-seq analysis and gene expression dynamics in the salivary glands of the argasid tick Ornithodoros erraticus along the trophogonic cycle

BACKGROUND

The argasid tick Ornithodoros erraticus is the main vector of tick-borne human relapsing fever (TBRF) and African swine fever (ASF) in the Mediterranean Basin. Tick salivary proteins secreted to the host at the feeding interface play critical roles for tick feeding and may contribute to host infection by tick-borne pathogens; accordingly, these proteins represent interesting antigen targets for the development of vaccines aimed at the control and prevention of tick infestations and tick-borne diseases.

METHODS

To identify these proteins, the transcriptome of the salivary glands of O. erraticus was de novo assembled and the salivary gene expression dynamics assessed throughout the trophogonic cycle using Illumina sequencing. The genes differentially upregulated after feeding were selected and discussed as potential antigen candidates for tick vaccines.

RESULTS

Transcriptome assembly resulted in 22,007 transcripts and 18,961 annotated transcripts, which represent 86.15% of annotation success. Most salivary gene expression took place during the first 7 days after feeding (2088 upregulated transcripts), while only a few genes (122 upregulated transcripts) were differentially expressed from day 7 post-feeding onwards. The protein families more abundantly overrepresented after feeding were lipocalins, acid and basic tail proteins, proteases (particularly metalloproteases), protease inhibitors, secreted phospholipases A2, 5'-nucleotidases/apyrases and heme-binding vitellogenin-like proteins. All of them are functionally related to blood ingestion and regulation of host defensive responses, so they can be interesting candidate protective antigens for vaccines.

CONCLUSIONS

The O. erraticus sialotranscriptome contains thousands of protein coding sequences-many of them belonging to large conserved multigene protein families-and shows a complexity and functional redundancy similar to those observed in the sialomes of other argasid and ixodid tick species. This high functional redundancy emphasises the need for developing multiantigenic tick vaccines to reach full protection. This research provides a set of promising candidate antigens for the development of vaccines for the control of O. erraticus infestations and prevention of tick-borne diseases of public and veterinary health relevance, such as TBRF and ASF. Additionally, this transcriptome constitutes a valuable reference database for proteomics studies of the saliva and salivary glands of O. erraticus.

Microscopic and Transcriptomic Analysis of Pollination Processes in Self-Incompatible Taraxacum koksaghyz

The transition of the Russian dandelion Taraxacum koksaghyz (Asteraceae) to a profitable, alternative crop producing natural rubber and inulin requires the optimization of several agronomic traits, cultivation conditions and harvesting procedures to improve the yield. However, efficient breeding is hindered by the obligatory sexual outcrossing of this species. Several other asters have been investigated to determine the mechanism of self-incompatibility, but the underlying molecular basis remains unclear. We therefore investigated the self-pollination and cross-pollination of two compatible T. koksaghyz varieties (TkMS2 and TkMS3) by microscopy and transcriptomic analysis to shed light on the pollination process. Self-pollination showed typical sporophytic self-incompatibility characteristics, with the rare pollen swelling at the pollen tube apex. In contrast, cross-pollination was characterized by pollen germination and penetration of the stigma by the growing pollen tubes. RNA-Seq was used to profile gene expression in the floret tissue during self-pollination and cross-pollination, and the differentially expressed genes were identified. This revealed three candidates for the early regulation of pollination in T. koksaghyz, which can be used to examine self-incompatibility mechanisms in more detail and to facilitate breeding programs.

Transcriptomic Diversity in the Livers of South African Sardines Participating in the Annual Sardine Run

During austral winter, the southern and eastern coastlines of South Africa witness one of the largest animal migrations on the planet, the KwaZulu-Natal sardine run. Hundreds of millions of temperate sardines, Sardinops sagax, form large shoals that migrate north-east towards the subtropical Indian Ocean. Recent studies have highlighted the role that genetic and environmental factors play in sardine run formation. In the present study, we used massively parallel sequencing to assemble and annotate the first reference transcriptome from the liver cells of South African sardines, and to investigate the functional content and transcriptomic diversity. A total of 1,310,530 transcripts with an N50 of 1578 bp were assembled de novo. Several genes and core biochemical pathways that modulate energy production, energy storage, digestion, secretory processes, immune responses, signaling, regulatory processes, and detoxification were identified. The functional content of the liver transcriptome from six individuals that participated in the 2019 sardine run demonstrated heterogeneous levels of variation. Data presented in the current study provide new insights into the complex function of the liver transcriptome in South African sardines.

Marine associated microbial consortium applied to RBBR textile dye detoxification and decolorization: Combined approach and metatranscriptomic analysis

The combination of different microorganisms and their metabolisms makes the use of microbial consortia in bioremediation processes a useful approach. In this sense, this study aimed at structuring and selecting a marine microbial consortium for Remazol Brilliant Blue R (RBBR) detoxification and decolorization. Experimental design was applied to improve the culture conditions, and metatranscriptomic analysis to understand the enzymatic pathways. A promising consortium composed of Mucor racemosus CBMAI 847, Marasmiellus sp. CBMAI 1062, Bacillus subtilis CBMAI 707, and Dietzia maris CBMAI 705 was selected. This consortium showed 52% of detoxification and 86% of decolorization in the validation assays after seven days of incubation in the presence of 500 ppm of RBBR. Reduction in RBBR color and toxicity were achieved by biosorption and microbial metabolisms. Metatranscriptomic data indicate that the consortium was able to decolorize and breakdown the RBBR molecule using a coordinated action of oxidases, oxygenases, and hydrolases. Epoxide hydrolases and glyoxalases expression could be associated with the decrease in toxicity. The efficiency of this marine microbial consortium suggests their use in bioremediation processes of textile effluents.

Comparative transcriptome analysis infers bulb derived in vitro cultures as a promising source for sipeimine biosynthesis in Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) - High value Himalayan medicinal herb

Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) is a critically endangered medicinal herb of immense importance due to its pharmaceutical bioactive compound, especially sipeimine, used for the treatment of chronic respiratory disorders. However, the industrial demand for sipeimine solely depends on its endangered natural habitat. Therefore; there is an utmost need for its biodiversity conservation as well as for the sustainable utilization of phytochemicals. Plant cell culture and transcriptomics-based molecular bioprospection of key regulatory genes involved in sipeimine biosynthesis as such will play a crucial role in exploring the unexplored traits, that are in supply crisis or nearly in extinction stage. De novo comparative transcriptome sequencing of the bulb (in vivo), callus, and regenerated plantlets (in vitro) resulted in more than 150 million high-quality paired-end clean reads that assembled into final 31,428 transcripts. Functional annotation and unigenes classification with multiple public databases such as KEGG, Refseq, Uniprot, TAIR, GO, and COG, etc. along with chemical structures and functional biocatalytic activity analysis of different steroidal alkaloids facilitated the identification of 30 unigenes specific to sipeimine biosynthesis. Additionally, ABC transporters and TFs like bHLH, MYC, MYB, and WRKY suggests their possible role in metabolite translocation and regulation in vivo as well as in vitro tissues. Differential gene expression and quantitative analysis revealed that the MVA pathway probably the predominant route for 5C intermediate (IPP & DMAPP) biosynthesis. Further, the genes involved in the downstream biosynthesis pathway viz. SQLE, CAS1, SMT1, SMO1, SMO2, SC5DL, DHCR7, DHCR24, CYP710A, 3β-HSD, CYP90D2, and CYP374A6 shown similar expression pattern with RNA-Seq and qRT-PCR findings. The positive correlation between higher expression of proposed biosynthetic pathway genes and relatively higher accumulation of sipeimine in differentiated naturally grown bulb tissues (in vivo), undifferentiated cells (callus), and de-differentiated tissues i.e. regenerated plantlets (in vitro) has been evident from the present study. Comprehensive genomic resources created in F. cirrhosa will provide strong evidence of bulb derived in vitro culture as an alternative promising source for steroidal alkaloids biosynthesis and metabolite upscaling through genetic and metabolic engineering.

Comparison of the Microsatellite Distribution Patterns in the Genomes of Euarchontoglires at the Taxonomic Level

Microsatellite or simple sequence repeat (SSR) instability within genes can induce genetic variation. The SSR signatures remain largely unknown in different clades within Euarchontoglires, one of the most successful mammalian radiations. Here, we conducted a genome-wide characterization of microsatellite distribution patterns at different taxonomic levels in 153 Euarchontoglires genomes. Our results showed that the abundance and density of the SSRs were significantly positively correlated with primate genome size, but no significant relationship with the genome size of rodents was found. Furthermore, a higher level of complexity for perfect SSR (P-SSR) attributes was observed in rodents than in primates. The most frequent type of P-SSR was the mononucleotide P-SSR in the genomes of primates, tree shrews, and colugos, while mononucleotide or dinucleotide motif types were dominant in the genomes of rodents and lagomorphs. Furthermore, (A)n was the most abundant motif in primate genomes, but (A)n, (AC)n, or (AG)n was the most abundant motif in rodent genomes which even varied within the same genus. The GC content and the repeat copy numbers of P-SSRs varied in different species when compared at different taxonomic levels, reflecting underlying differences in SSR mutation processes. Notably, the CDSs containing P-SSRs were categorized by functions and pathways using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes annotations, highlighting their roles in transcription regulation. Generally, this work will aid future studies of the functional roles of the taxonomic features of microsatellites during the evolution of mammals in Euarchontoglires.

Transcriptome analysis reveals differentially expressed MYB transcription factors associated with silicon response in wheat

Silicon plays a vital role in plant growth. However, molecular mechanisms in response to silicon have not previously been studied in wheat. In this study, we used RNA-seq technology to identify differentially expressed genes (DEGs) in wheat seedlings treated with silicon. Results showed that many wheat genes responded to silicon treatment, including 3057 DEGs, of which 6.25% (191/3057) were predicted transcription factors (TFs). Approximately 14.67% (28 out of 191) of the differentially expressed TFs belonged to the MYB TF family. Gene ontology (GO) enrichment showed that the highly enriched DEGs were responsible for secondary biosynthetic processes. According to KEGG pathway analysis, the DEGs were related to chaperones and folding catalysts, phenylpropanoid biosynthesis, and protein processing in the endoplasmic reticulum. Moreover, 411 R2R3-MYB TFs were identified in the wheat genome, all of which were classified into 15 groups and accordingly named S1–S15. Among them, 28 were down-regulated under silicon treatment. This study revealed the essential role of MYB TFs in the silicon response mechanism of plants, and provides important genetic resources for breeding silicon-tolerant wheat.

Identifications of immune-responsive genes for adaptative traits by comparative transcriptome analysis of spleen tissue from Kazakh and Suffolk sheep

Aridity and heat are significant environmental stressors that affect sheep adaptation and adaptability, thus influencing immunity, growth, reproduction, production performance, and profitability. The aim of this study was to profile mRNA expression levels in the spleen of indigenous Kazakh sheep breed for comparative analysis with the exotic Suffolk breed. Spleen histomorphology was observed in indigenous Kazakh sheep and exotic Suffolk sheep raised in Xinjiang China. Transcriptome sequencing of spleen tissue from the two breeds were performed via Illumina high-throughput sequencing technology and validated by RT-qPCR. Blood cytokine and IgG levels differed between the two breeds and IgG and IL-1β were significantly higher in Kazakh sheep than in Suffolk sheep (p < 0.05), though spleen tissue morphology was the same. A total of 52.04 Gb clean reads were obtained and the clean reads were assembled into 67,271 unigenes using bioinformatics analysis. Profiling analysis of differential gene expression showed that 1158 differentially expressed genes were found when comparing Suffolk with Kazakh sheep, including 246 up-regulated genes and 912 down-regulated genes. Utilizing gene ontology annotation and pathway analysis, 21 immune- responsive genes were identified as spleen-specific genes associated with adaptive traits and were significantly enriched in hematopoietic cell lineage, natural killer cell-mediated cytotoxicity, complement and coagulation cascades, and in the intestinal immune network for IgA production. Four pathways and up-regulated genes associated with immune responses in indigenous sheep played indispensable and promoting roles in arid and hot environments. Overall, this study provides valuable transcriptome data on the immunological mechanisms related to adaptive traits in indigenous and exotic sheep and offers a foundation for research into adaptive evolution.

Transcriptomic analysis reveals the regulation network of BmKrüppel homolog 1 in the oocyte development of Bombyx mori

Krüppel homolog 1 (Kr-h1), a zinc finger transcription factor, is involved in the metamorphosis and adult reproduction of insects. However, the role of Kr-h1 in reproduction of holometabolic insects remains to be elucidated. The regulation network of Kr-h1-associated genes in the reproduction in Bombyx mori was investigated in this study. The higher expression level of BmKr-h1 in the ovaries was detected during the late pupal stage and adults. RNA interference (RNAi)-mediated depletion of BmKr-h1 in the female at day 6 of pupae resulted in abnormal oocytes at 48 h post-double-stranded RNA treatment, which showed less yolk protein deposition and partially transparent chorion. RNA-seq and subsequent differentially expressed transcripts analysis showed that knockdown of BmKr-h1 caused a decrease in the expression of 2882 genes and an increase in the expression of 2565 genes in the oocytes at day 8 of pupae. Totally, 27 genes coding for transcription factors were down-regulated, while six genes coding for other transcription factors were up-regulated. BmKr-h1 bound to the Kr-h1 binding site of the transcription factors AP-1 (activating protein-1) and FOXG1 to increase their messenger RNA transcripts in the BmN cells, respectively. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses of that positively co-expressed with AP-1 and FOXG1 transcripts showed mainly enrichment in the metabolic-related pathways, the nutrient absorption and the yolk protein absorption processes. These data suggested that BmKr-h1 might directly regulate the metabolic-related pathways, the nutrient absorption and the yolk protein absorption processes or probably through AP-1 and /or FOXG1 to regulate oocyte development.

Comparative transcriptomics and host-specific parasite gene expression profiles inform on drivers of proliferative kidney disease

The myxozoan parasite, Tetracapsuloidesbryosalmonae has a two-host life cycle alternating between freshwater bryozoans and salmonid fish. Infected fish can develop Proliferative Kidney Disease, characterised by a gross lymphoid-driven kidney pathology in wild and farmed salmonids. To facilitate an in-depth understanding of T.bryosalmonae-host interactions, we have used a two-host parasite transcriptome sequencing approach in generating two parasite transcriptome assemblies; the first derived from parasite spore sacs isolated from infected bryozoans and the second from infected fish kidney tissues. This approach was adopted to minimize host contamination in the absence of a complete T.bryosalmonae genome. Parasite contigs common to both infected hosts (the intersect transcriptome; 7362 contigs) were typically AT-rich (60–75% AT). 5432 contigs within the intersect were annotated. 1930 unannotated contigs encoded for unknown transcripts. We have focused on transcripts encoding proteins involved in; nutrient acquisition, host–parasite interactions, development, cell-to-cell communication and proteins of unknown function, establishing their potential importance in each host by RT-qPCR. Host-specific expression profiles were evident, particularly in transcripts encoding proteases and proteins involved in lipid metabolism, cell adhesion, and development. We confirm for the first time the presence of homeobox proteins and a frizzled homologue in myxozoan parasites. The novel insights into myxozoan biology that this study reveals will help to focus research in developing future disease control strategies.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

RNA-Seq Highlights Molecular Events Associated With Impaired Pollen-Pistil Interactions Following Short-Term Heat Stress in Brassica napus

The global climate change is leading to increased frequency of heatwaves with crops getting exposed to extreme temperature events. Such temperature spikes during the reproductive stage of plant development can harm crop fertility and productivity. Here we report the response of short-term heat stress events on the pollen and pistil tissues in a commercially grown cultivar of Brassica napus. Our data reveals that short-term temperature spikes not only affect pollen fitness but also impair the ability of the pistil to support pollen germination and pollen tube growth and that the heat stress sensitivity of pistil can have severe consequences for seed set and yield. Comparative transcriptome profiling of non-stressed and heat-stressed (40°C for 30 min) pollen and pistil (stigma + style) highlighted the underlying cellular mechanisms involved in heat stress response in these reproductive tissues. In pollen, cell wall organization and cellular transport-related genes possibly regulate pollen fitness under heat stress while the heat stress-induced repression of transcription factor encoding transcripts is a feature of the pistil response. Overall, high temperature altered the expression of genes involved in protein processing, regulation of transcription, pollen-pistil interactions, and misregulation of cellular organization, transport, and metabolism. Our results show that short episodes of high-temperature exposure in B. napus modulate key regulatory pathways disrupted reproductive processes, ultimately translating to yield loss. Further investigations on the genes and networks identified in the present study pave a way toward genetic improvement of the thermotolerance and reproductive performance of B. napus varieties.

Temperature modulates immune gene expression in mosquitoes during arbovirus infection

The principal vector of dengue, Zika and chikungunya viruses is the mosquito Aedes aegypti, with its ability to transmit pathogens influenced by ambient temperature. We use chikungunya virus (CHIKV) to understand how the mosquito transcriptome responds to arbovirus infection at different ambient temperatures. We exposed CHIKV-infected mosquitoes to 18, 28 and 32°C, and found that higher temperature correlated with higher virus levels, particularly at 3 days post infection, but lower temperature resulted in reduced virus levels. RNAseq analysis indicated significantly altered gene expression levels in CHIKV infection. The highest number of significantly differentially expressed genes was observed at 28°C, with a more muted effect at the other temperatures. At the higher temperature, the expression of many classical immune genes, including Dicer-2, was not substantially altered in response to CHIKV. The upregulation of Toll, IMD and JAK-STAT pathways was only observed at 28°C. Functional annotations suggested that genes in immune response and metabolic pathways related to energy supply and DNA replication were involved in temperature-dependent changes. Time post infection also led to substantially different gene expression profiles, and this varied with temperature. In conclusion, temperature significantly modulates mosquito gene expression in response to infection, potentially leading to impairment of immune defences at higher temperatures.

Landscape of meiotic crossovers in Hericium erinaceus

Meiotic crossover shows marked interspecific and intraspecific variation, and knowledge about the molecular mechanism of crossover variation remains limited. Herein, we described the genome-wide scanning of crossover in one mushroom-forming fungus Hericium erinaceus. Utilizing the whole-genome single-nucleotide polymorphism (SNP) data-sets of a 127 F₁ haploid progeny, we localized a total of 1316 crossover events and found that they were more likely to occur in the genic than intergenic regions. More than 30 % of the crossovers were concentrated in 59 crossover hotspots that were preferentially located close to chromosome ends. We then examined the genomic features around crossover hotspots. Results showed that the crossover hotspots were associated with increased gene density and guanine-cytosine (GC) content. An 8-bp GC-rich motif (GCGTCAGC) was found to be significantly enriched in these hotspots. The presence of mating-type loci affected the crossover at local scale rather than the overall crossover number. In order to dissect the genetic mechanisms shaping crossover variation, we then conducted quantitative trait locus (QTL) mapping for the total crossovers (TCO) and the crossover events that solely occurred within hotspots (HCO). Genome-wide QTL scanning identified four TCO-QTLs and two HCO-QTLs, which all located within or next to the crossover-hotspots. Crossover variations were shaped by multiple small-effect loci, with individual QTL contributing 6.9 %-11.7 % of variation. A few recombination pathway genes, including Spo11, Msh5, and Smc5 were found to be co-localized with the mapped crossover QTLs. Taken together, findings of this study offer insights into the crossover distribution and genetic factors conferring crossover variation in H. erinaceus, and advance our understandings for meiotic recombination in mushroom-forming fungi.

A phased genome based on single sperm sequencing reveals crossover pattern and complex relatedness in tea plants

For diploid organisms that are highly heterozygous, a phased haploid genome can greatly aid in functional genomic, population genetic and breeding studies. Based on the genome sequencing of 135 single sperm cells of the elite tea cultivar 'Fudingdabai', we herein phased the genome of Camellia sinensis, one of the most popular beverage crops worldwide. High-resolution genetic and recombination maps of Fudingdabai were constructed, which revealed that crossover (CO) positions were frequently located in the 5' and 3' ends of annotated genes, while CO distributions across the genome were random. The low CO frequency in tea can be explained by strong CO interference, and CO simulation revealed the proportion of interference insensitive CO ranged from 5.2% to 11.7%. We furthermore developed a method to infer the relatedness between tea accessions and detected complex kinship and genetic signatures of 106 tea accessions. Among them, 59 accessions were closely related with Fudingdabai and 31 of them were first-degree relatives. We additionally identified genes displaying allele specific expression patterns between the two haplotypes of Fudingdabai and genes displaying significantly differential expression levels between Fudingdabai and other haplotypes. These results lay the foundation for further investigation of genetic and epigenetic factors underpinning the regulation of gene expression and provide insights into the evolution of tea plants as well as a valuable genetic resource for future breeding efforts.

Genome-wide identification, characterization and transcriptional profiling of NHX-type (Na+/H+) antiporters under salinity stress in soybean

This study was aimed at the genome-wide identification, a comprehensive in silico characterization of NHX genes from soybean (Glycine max L.) and their tissue-specific expression under varied levels (0-200 mM NaCl) of salinity stress. A total of nine putative NHX genes were identified from soybean. The phylogenetic analysis confirmed a total of five sub-groups and GmNHXs were distributed in three of them. Bioinformatics analyses confirmed all GmNHXs as ion transporters in nature, and all were localized on the vacuolar membrane. Several cis-acting regulatory elements involved in hormonal signal-responsiveness and abiotic stress including salinity responses were identified in the promoter regions of GmNHXs. Amiloride, which is a known Na⁺/H⁺ exchanger activity inhibitor, binding motifs were observed in all the GmNHXs. Furthermore, the identified GmNHXs were predicted-targets of 75 different miRNA candidates. To gain an insight into the functional divergence of GmNHX transporters, qRT-PCR based gene expression analysis was done in control and salt-treated root, stem and leaf tissues of two contrasting Indian soybean varieties MAUS-47 (tolerant) and Gujosoya-2 (sensitive). The gene up-regulation was tissue-specific and varied amongst the soybean varieties, with higher induction in tolerant variety. Maximum induction was observed in GmNHX2 in root tissues of MAUS-47 at 200 mM NaCl stress. Overall, identified GmNHXs may be explored further as potential gene candidates for soybean improvement.

Comparative genomic analysis of the principal Cryptosporidium species that infect humans

Cryptosporidium parasites are ubiquitous and can infect a broad range of vertebrates and are considered the most frequent protozoa associated with waterborne parasitic outbreaks. The intestine is the target of three of the species most frequently found in humans: C. hominis, C. parvum, and. C. meleagridis. Despite the recent advance in genome sequencing projects for this apicomplexan, a broad genomic comparison including the three species most prevalent in humans have not been published so far. In this work, we downloaded raw NGS data, assembled it under normalized conditions, and compared 23 publicly available genomes of C. hominis, C. parvum, and C. meleagridis. Although few genomes showed highly fragmented assemblies, most of them had less than 500 scaffolds and mean coverage that ranged between 35X and 511X. Synonymous single nucleotide variants were the most common in C. hominis and C. meleagridis, while in C. parvum, they accounted for around 50% of the SNV observed. Furthermore, deleterious nucleotide substitutions common to all three species were more common in genes associated with DNA repair, recombination, and chromosome-associated proteins. Indel events were observed in the 23 studied isolates that spanned up to 500 bases. The highest number of deletions was observed in C. meleagridis, followed by C. hominis, with more than 60 species-specific deletions found in some isolates of these two species. Although several genes with indel events have been partially annotated, most of them remain to encode uncharacterized proteins.

Differential gene expression reveals mechanisms related to habitat divergence between hybridizing orchids from the Neotropical coastal plains

BACKGROUND

Closely related hybridizing species are ideal systems for identifying genomic regions underlying adaptive divergence. Although gene expression plays a central role in determining ecologically-based phenotypic differences, few studies have inferred the role of gene expression for adaptive divergence in Neotropical systems. In this study, we conduct genome-wide expression analysis alongside soil elemental analysis in sympatric and allopatric populations of Epidendrum fulgens and E. puniceoluteum (Orchidaceae), which occur in contrasting adjacent habitats in the Neotropical coastal plains.

RESULTS

These species were highly differentiated by their gene expression profiles, as determined by 18-21% of transcripts. Gene ontology (GO) terms associated with reproductive processes were enriched according to comparisons between species in both allopatric and sympatric populations. Species showed differential expression in genes linked to salt and waterlogging tolerance according to comparisons between species in sympatry, and biological processes related to environmental stimulus appeared as representative among those transcripts associated with edaphic characteristics in each sympatric zone. Hybrids, in their turn, were well differentiated from E. fulgens, but exhibited a similar gene expression profile to flooding-tolerant E. puniceolutem. When compared with parental species, hybrids showed no transcripts with additive pattern of expression and increased expression for almost all transgressive transcripts.

CONCLUSIONS

This study sheds light on general mechanisms promoting ecological differentiation and assortative mating, and suggests candidate genes, such as those encoding catalase and calcium-dependent protein kinase, underling adaptation to harsh edaphic conditions in the Neotropical coastal plains. Moreover, it demonstrates that differential gene expression plays a central role in determining ecologically-based phenotypic differences among co-occurring species and their hybrids.

Comprehensive Transcriptome of the Maize Stalk Borer, Busseola fusca, from Multiple Tissue Types, Developmental Stages, and Parasitoid Wasp Exposures

Busseola fusca (Fuller) (Lepidoptera: Noctuidae), the maize stalk borer, is a widespread crop pest in sub-Saharan Africa that has been the focus of biological research and intensive management strategies. Here, we present a comprehensive annotated transcriptome of B. fusca (originally collected in the Western Province of Kenya) based on ten pooled libraries including a wide array of developmental stages, tissue types, and exposures to parasitoid wasps. Parasitoid wasps have been used as a form of biocontrol to try and reduce crop losses with variable success, in part due to differential infectivities and immune responses among wasps and hosts. We identified a number of loci of interest for pest management, including genes potentially involved in chemoreception, immunity, and response to insecticides. The comprehensive sampling design used expands our current understanding of the transcriptome of this species and deepens the list of potential target genes for future crop loss mitigation, in addition to highlighting candidate loci for differential expression and functional genetic analyses in this important pest species.

A salivary chitinase of Varroa destructor influences host immunity and mite's survival

Varroa destructor is an ectoparasite of honey bees and an active disease vector, which represents one of the most severe threats for the beekeeping industry. This parasitic mite feeds on the host’s body fluids through a wound in the cuticle, which allows food uptake by the mother mite and its progeny, offering a potential route of entrance for infecting microorganisms. Mite feeding is associated with saliva injection, whose role is still largely unknown. Here we try to fill this gap by identifying putative host regulation factors present in the saliva of V. destructor and performing a functional analysis for one of them, a chitinase (Vd-CHIsal) phylogenetically related to chitinases present in parasitic and predatory arthropods, which shows a specific and very high level of expression in the mite’s salivary glands. Vd-CHIsal is essential for effective mite feeding and survival, since it is apparently involved both in maintaining the feeding wound open and in preventing host infection by opportunistic pathogens. Our results show the important role in the modulation of mite-honey bee interactions exerted by a host regulation factor shared by different evolutionary lineages of parasitic arthropods. We predict that the functional characterization of Varroa sialome will provide new background knowledge on parasitism evolution in arthropods and the opportunity to develop new bioinspired strategies for mite control based on the disruption of their complex interactions with a living food source.

Varroa destructor is a parasitic mite of honey bees and a major driver of honey bee colony losses. The feeding mites inject a salivary blend of poorly known molecules, which regulate host physiology. Here, we have identified in silico putative host regulation factors occurring in Varroa saliva and characterized the functional role of a highly expressed chitinase, which is conserved across different evolutionary lineages of parasitic arthropods. This enzyme influences host immune response and mite’s survival. An in-depth functional analysis of Varroa saliva will shed light on parasitism evolution in arthropods and will pave the way towards the development of new bioinspired strategies for mite control.

A Novel Family of Acinetobacter Mega-Plasmids Are Disseminating Multi-Drug Resistance Across the Globe While Acquiring Location-Specific Accessory Genes

Acinetobacter species are emerging as major nosocomial pathogens, aided by their ability to acquire resistance to all classes of antibiotics. A key factor leading to their multi-drug resistance phenotypes is the acquisition of a wide variety of mobile genetic elements, particularly large conjugative plasmids. Here, we characterize a family of 21 multi-drug resistance mega-plasmids in 11 different Acinetobacter species isolated from various locations across the globe. The plasmid family exhibits a highly dynamic and diverse accessory genome, including 221 antibiotic resistance genes (ARGs) that confer resistance to 13 classes of antibiotics. We show that plasmids isolated within the same geographic region are often evolutionarily divergent members of this family based on their core-genome, yet they exhibit a more similar accessory genome. Individual plasmids, therefore, can disseminate to different locations around the globe, where they then appear to acquire diverse sets of accessory genes from their local surroundings. Further, we show that plasmids from several geographic regions were enriched with location-specific functional traits. Together, our findings show that these mega-plasmids can transmit across species boundaries, have the capacity for global dissemination, can accumulate a diverse suite of location-specific accessory genes, and can confer multi-drug resistance phenotypes of significant concern for human health. We therefore highlight this previously undescribed plasmid family as a serious threat to healthcare systems worldwide. These findings also add to the growing concern that mega-plasmids are key disseminators of antibiotic resistance and require global surveillance.

Metabolic flux and transcriptome analyses provide insights into the mechanism underlying zinc sulfate improved β-1,3-D-glucan production by Aureobasidium pullulans

The effects of zinc sulfate at various concentrations on β-1,3-D-glucan (β-glucan) and pullulan production were investigated in flasks, and 0.1 g/L zinc sulfate was found to be the optimum concentration favoring increased β-glucan production. When batch culture of Aureobasidium pullulans CCTCC M 2012259 with 0.1 g/L zinc sulfate was carried out, the maximum dry biomass decreased by 16.9% while β-glucan production significantly increased by 120.5%, compared to results obtained from the control without zinc sulfate addition. To reveal the mechanism underlying zinc sulfate improved β-glucan production, both metabolic flux analysis and RNA-seq analysis were performed. The results indicated that zinc sulfate decreased carbon flux towards biomass formation and ATP supply, down-regulated genes associated with membrane part and cellular components organization, leading to a decrease in dry cell weight. However, zinc sulfate increased metabolic flux towards β-glucan biosynthesis, up-regulated genes related to glycan biosynthesis and nucleotide metabolism, resulting in improved β-glucan production. This study provides insights into the changes in the metabolism of A. pullulans in response to zinc sulfate, and can serve as a valuable reference of genetic information for improving the production of polysaccharides through metabolic engineering.

Integrated transcriptome and proteome analysis provides insight into chilling-induced dormancy breaking in Chimonanthus praecox

Chilling has a critical role in the growth and development of perennial plants. The chilling requirement (CR) for dormancy breaking largely depends on the species. However, global warming is expected to negatively affect chilling accumulation and dormancy release in a wide range of perennial plants. Here, we used Chimonanthus praecox as a model to investigate the CR for dormancy breaking under natural and artificial conditions. We determined the minimum CR (570 chill units, CU) needed for chilling-induced dormancy breaking and analyzed the transcriptomes and proteomes of flowering and non-flowering flower buds (FBs, anther and ovary differentiation completed) with different CRs. The concentrations of ABA and GA3 in the FBs were also determined using HPLC. The results indicate that chilling induced an upregulation of ABA levels and significant downregulation of SHORT VEGETATIVE PHASE (SVP) and FLOWERING LOCUS T (FT) homologs at the transcript level in FBs when the accumulated CR reached 570 CU (IB570) compared to FBs in November (FB.Nov, CK) and nF16 (non-flowering FBs after treatment at 16 °C for -300 CU), which suggested that dormancy breaking of FBs could be regulated by the ABA-mediated SVP-FT module. Overexpression in Arabidopsis was used to confirm the function of candidate genes, and early flowering was induced in 35S::CpFT1 transgenic lines. Our data provide insight into the minimum CR (570 CU) needed for chilling-induced dormancy breaking and its underlying regulatory mechanism in C. praecox, which provides a new tool for the artificial regulation of flowering time and a rich gene resource for controlling chilling-induced blooming.

Sequencing and analysis of gerbera daisy leaf transcriptomes reveal disease resistance and susceptibility genes differentially expressed and associated with powdery mildew resistance

BACKGROUND

RNA sequencing has been widely used to profile genome-wide gene expression and identify candidate genes controlling disease resistance and other important traits in plants. Gerbera daisy is one of the most important flowers in the global floricultural trade, and powdery mildew (PM) is the most important disease of gerbera. Genetic improvement of gerbera PM resistance has become a crucial goal in gerbera breeding. A better understanding of the genetic control of gerbera resistance to PM can expedite the development of PM-resistant cultivars.

RESULTS

The objectives of this study were to identify gerbera genotypes with contrasting phenotypes in PM resistance and sequence and analyze their leaf transcriptomes to identify disease resistance and susceptibility genes differentially expressed and associated with PM resistance. An additional objective was to identify SNPs and SSRs for use in future genetic studies. We identified two gerbera genotypes, UFGE 4033 and 06-245-03, that were resistant and susceptible to PM, respectively. De novo assembly of their leaf transcriptomes using four complementary pipelines resulted in 145,348 transcripts with a N50 of 1124 bp, of which 67,312 transcripts contained open reading frames and 48,268 were expressed in both genotypes. A total of 494 transcripts were likely involved in disease resistance, and 17 and 24 transcripts were up- and down-regulated, respectively, in UFGE 4033 compared to 06-245-03. These gerbera disease resistance transcripts were most similar to the NBS-LRR class of plant resistance genes conferring resistance to various pathogens in plants. Four disease susceptibility transcripts (MLO-like) were expressed only or highly expressed in 06-245-03, offering excellent candidate targets for gene editing for PM resistance in gerbera. A total of 449,897 SNPs and 19,393 SSRs were revealed in the gerbera transcriptomes, which can be a valuable resource for developing new molecular markers.

CONCLUSION

This study represents the first transcriptomic analysis of gerbera PM resistance, a highly important yet complex trait in a globally important floral crop. The differentially expressed disease resistance and susceptibility transcripts identified provide excellent targets for development of molecular markers and genetic maps, cloning of disease resistance genes, or targeted mutagenesis of disease susceptibility genes for PM resistance in gerbera.

Comparative transcriptome analysis implied a ZEP paralog was a key gene involved in carotenoid accumulation in yellow-fleshed sweetpotato

The mechanisms of carotenoid accumulation in yellow-fleshed sweetpotato cultivars are unclear. In this study, we compared the transcriptome profiles of a yellow-fleshed cultivar, Beniharuka (BH) and two of its spontaneous white-fleshed mutants (WH2 and WH3) to reveal the genes involved in yellow flesh. As a result of RNA sequencing, a total of 185 differentially expressed genes (DEGs) were commonly detected in WH2 and WH3 compared to BH. Of these genes, 85 DEGs and 100 DEGs were commonly upregulated and downregulated in WH2 and WH3 compared to BH, respectively. g1103.t1, a paralog of zeaxanthin epoxidase (ZEP), was only DEG common to WH2 and WH3 among 38 genes considered to be involved in carotenoid biosynthesis in storage roots. The expression level of g1103.t1 was also considerably lower in five white-fleshed cultivars than in five yellow-fleshed cultivars. Analysis of carotenoid composition in the storage roots showed that the epoxidised carotenoids were drastically reduced in both WH2 and WH3. Therefore, we propose that the ZEP paralog, g1103.t1, may be involved in carotenoid accumulation through the epoxidation of β-carotene and β-cryptoxanthin in sweetpotato.

Transcriptome analyses reveals the dynamic nature of oil accumulation during seed development of Plukenetia volubilis L

Sacha inchi (Plukenetia volubilis L.) is a shrub native to Amazon rainforests that’s of commercial interest as its seeds contain 35–60% edible oil (dry weight). This oil is one of the healthiest vegetable oils due to its high polyunsaturated fatty acid content and favourable ratio of omega-6 to omega-3 fatty acids. De novo transcriptome assembly and comparative analyses were performed on sacha inchi seeds from five stages of seed development in order to identifying genes associated with oil accumulation and fatty acid production. Of 30,189 unigenes that could be annotated in public databases, 20,446 were differentially expressed unigenes. A total of 14 KEGG pathways related to lipid metabolism were found, and 86 unigenes encoding enzymes involved in α-linolenic acid (ALA) biosynthesis were obtained including five unigenes encoding FATA (Unigene0008403), SAD (Unigene0012943), DHLAT (Unigene0014324), α-CT (Unigene0022151) and KAS II (Unigene0024371) that were significantly up-regulated in the final stage of seed development. A total of 66 unigenes encoding key enzymes involved in the synthesis of triacylglycerols (TAGs) were found, along with seven unigenes encoding PDCT (Unigene0000909), LPCAT (Unigene0007846), Oleosin3 (Unigene0010027), PDAT1 (Unigene0016056), GPDH (Unigene0022660), FAD2 (Unigene0037808) and FAD3 (Unigene0044238); these also proved to be up-regulated in the final stage of seed development.

Turnip mosaic virus in oilseed rape activates networks of sRNA-mediated interactions between viral and host genomes

Virus-induced plant diseases in cultivated plants cause important damages in yield. Although the mechanisms of virus infection are intensely studied at the cell biology level, only little is known about the molecular dialog between the invading virus and the host genome. Here we describe a combinatorial genome-wide approach to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus (TuMV) infection sites in Brassica napus leaves. We show that the induction of host-encoded, virus-activated small interfering RNAs (vasiRNAs) observed in virus-infected tissues is accompanied by site-specific cleavage events on both viral and host RNAs that recalls the activity of small RNA-induced silencing complexes (RISC). Cleavage events also involve virus-derived siRNA (vsiRNA)–directed cleavage of target host transcripts as well as cleavage of viral RNA by both host vasiRNAs and vsiRNAs. Furthermore, certain coding genes act as virus-activated regulatory hubs to produce vasiRNAs for the targeting of other host genes. The observations draw an advanced model of plant-virus interactions and provide insights into the complex regulatory networking at the plant-virus interface within cells undergoing early stages of infection.

Pitzalis et al. use replicative RNAseq, small RNA (sRNA)seq, and parallel analysis of RNA ends (PARE)seq analysis to identify networks of sRNAs-guided post-transcriptional regulation within local Turnip mosaic virus infection sites. This study provides insights into the complex regulatory networking at the plantvirus interface within cells undergoing early stages of infection.

Proteomics of extracellular vesicles produced by Granulicatella adiacens, which causes infective endocarditis

When oral bacteria accidentally enter the bloodstream due to transient tissue damage during dental procedures, they have the potential to attach to the endocardium or an equivalent surface of an indwelling prosthesis and cause infection. Many bacterial species produce extracellular vesicles (EVs) as part of normal physiology, but also use it as a virulence strategy. In this study, it was hypothesized that Granulicatella adiacens produce EVs that possibly help it in virulence. Therefore, the objectives were to isolate and characterize EVs produced by G. adiacens and to investigate its immune-stimulatory effects. The reference strain G. adiacens CCUG 27809 was cultured on chocolate blood agar for 2 days. From subsequent broth culture, the EVs were isolated using differential centrifugation and filtration protocol and then observed using scanning electron microscopy. Proteins in the vesicle preparation were identified by nano LC-ESI-MS/MS. The EVs proteome was analyzed and characterized using different bioinformatics tools. The immune-stimulatory effect of the EVs was studied via ELISA quantification of IL-8, IL-1β and CCL5, major proinflammatory cytokines, produced from stimulated human PBMCs. It was revealed that G. adiacens produced EVs, ranging in diameter from 30 to 250 nm. Overall, G. adiacens EVs contained 112 proteins. The proteome consists of several ribosomal proteins, DNA associated proteins, binding proteins, and metabolic enzymes. It was also shown that these EVs carry putative virulence factors including moonlighting proteins. These EVs were able to induce the production of IL-8, IL-1β and CCL5 from human PBMCs. Further functional characterization of the G. adiacens EVs may provide new insights into virulence mechanisms of this important but less studied oral bacterial species.

Proteomics and Lipidomics Investigations to Decipher the Behavior of Willaertia magna C2c Maky According to Different Culture Modes

Willaertia magna C2c Maky is a free-living amoeba that has demonstrated its ability to inhibit the intracellular multiplication of some Legionella pneumophila strains, which are pathogenic bacteria inhabiting the aquatic environment. The Amoeba, an industry involved in the treatment of microbiological risk in the water and plant protection sectors, has developed a natural biocide based on the property of W. magna to manage the proliferation of the pathogen in cooling towers. In axenic liquid medium, amoebas are usually cultivated in adhesion on culture flask. However, we implemented a liquid culture in suspension using bioreactors in order to produce large quantities of W. magna. In order to investigate the culture condition effects on W. magna, we conducted a study based on microscopic, proteomics and lipidomics analyzes. According to the culture condition, amoeba exhibited two different phenotypes. The differential proteomics study showed that amoebas seemed to promote the lipid metabolism pathway in suspension culture, whereas we observed an upregulation of the carbohydrate pathway in adherent culture. Furthermore, we observed an over-regulation of proteins related to the cytoskeleton for W. magna cells grown in adhesion. Regarding the lipid analysis, suspension and adhesion cell growth showed comparable lipid class compositions. However, the differential lipid analysis revealed differences that confirmed cell phenotype differences observed by microscopy and predicted by proteomics. Overall, this study provides us with a better insight into the biology and molecular processes of W. magna in different culture lifestyles.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

Lipidomic analyses reveal enhanced lipolysis in planthoppers feeding on resistant host plants

The brown planthopper (BPH) (Nilaparvata lugens Stål) is a highly destructive pest that seriously damages rice (Oryza sativa L.) and causes severe yield losses. To better understand the physiological and metabolic mechanisms through which BPHs respond to resistant rice, we combined mass-spectrometry-based lipidomics with transcriptomic analysis and gene knockdown techniques to compare the lipidomes of BPHs feeding on either of the two resistant (NIL-Bph6 and NIL-Bph9) plants or a wild-type, BPH susceptible (9311) plant. Insects that were fed on resistant rice transformed triglyceride (TG) to phosphatidylcholine (PC) and digalactosyldiacylglycerol (DGDG), with these lipid classes showing significant alterations in fatty acid composition. Moreover, the insects that were fed on resistant rice were characterized by prominent expression changes in genes involved in lipid metabolism processes. Knockdown of the NlBmm gene, which encodes a lipase that regulates the mobilization of lipid reserves, significantly increased TG content and feeding performance of BPHs on resistant plants relative to dsGFP-injected BPHs. Our study provides the first detailed description of lipid changes in BPHs fed on resistant and susceptible rice genotypes. Results from BPHs fed on resistant rice plants reveal that these insects can accelerate TG mobilization to provide energy for cell proliferation, body maintenance, growth and oviposition.

Airborne fungus-induced biosynthesis of anthocyanins in Arabidopsis thaliana via jasmonic acid and salicylic acid signaling

Anthocyanins are plant-specific pigments, the biosynthesis of which is stimulated by pathogen infection in several plant species. A. thaliana seedlings injected with airborne fungi can accumulate a high content of anthocyanins. The mechanism involved in fungus-induced anthocyanin accumulation in plants has not been fully described. In this study, the fungus Penicillium corylophilum (P. corylophilum), isolated from an Arabidopsis culture chamber, triggered jasmonic acid (JA), salicylic acid (SA), and anthocyanin accumulation in A. thaliana. Inhibitors of JA and SA biosynthesis suppressed the anthocyanin accumulation induced by P. corylophilum. The anthocyanin content was minimal in both the null mutant of JA-receptor coi1 and the null mutant of SA-receptor npr1 under P. corylophilum stimulation. The results indicate that JA and SA signaling mediated fungus-induced anthocyanin biosynthesis in A. thaliana. P. corylophilum led to different levels of anthocyanin generation in null mutants for MYB75, bHLH, EGL3, and GL3 transcription factors and WD40 protein, demonstrating that multiple MYB-bHLH-WD40 transcription factor complexes participated in fungus-induced anthocyanin accumulation in A. thaliana. The present study will help further elucidate the mechanism of plant resistance to pathogen infection.

Reference transcriptomes and comparative analyses of six species in the threatened rosewood genus Dalbergia

Dalbergia is a pantropical genus with more than 250 species, many of which are highly threatened due to overexploitation for their rosewood timber, along with general deforestation. Many Dalbergia species have received international attention for conservation, but the lack of genomic resources for Dalbergia hinders evolutionary studies and conservation applications, which are important for adaptive management. This study produced the first reference transcriptomes for 6 Dalbergia species with different geographical origins and predicted ~ 32 to 49 K unique genes. We showed the utility of these transcriptomes by phylogenomic analyses with other Fabaceae species, estimating the divergence time of extant Dalbergia species to ~ 14.78 MYA. We detected over-representation in 13 Pfam terms including HSP, ALDH and ubiquitin families in Dalbergia. We also compared the gene families of geographically co-occurring D. cochinchinensis and D. oliveri and observed that more genes underwent positive selection and there were more diverged disease resistance proteins in the more widely distributed D. oliveri, consistent with reports that it occupies a wider ecological niche and has higher genetic diversity. We anticipate that the reference transcriptomes will facilitate future population genomics and gene-environment association studies on Dalbergia, as well as contributing to the genomic database where plants, particularly threatened ones, are currently underrepresented.

Transcriptomic signatures of cold adaptation and heat stress in the winter ant (Prenolepis imparis)

Climate change is a serious threat to biodiversity; it is therefore important to understand how animals will react to this stress. Ectotherms, such as ants, are especially sensitive to the climate as the environmental temperature influences myriad aspects of their biology, from optimal foraging time to developmental rate. In this study, we conducted an RNA-seq analysis to identify stress-induced genes in the winter ant (Prenolepis imparis). We quantified gene expression during heat and cold stress relative to a control temperature. From each of our conditions, we sequenced the transcriptome of three individuals. Our de novo assembly included 13,324 contigs that were annotated against the nr and SwissProt databases. We performed gene ontology and enrichment analyses to gain insight into the physiological processes involved in the stress response. We identified a total of 643 differentially expressed genes across both treatments. Of these, only seven genes were differentially expressed in the cold-stressed ants, which could indicate that the temperature we chose for trials did not induce a strong stress response, perhaps due to the cold adaptations of this species. Conversely, we found a strong response to heat: 426 upregulated genes and 210 downregulated genes. Of these, ten were expressed at a greater than ten-fold change relative to the control. The transcripts we could identify included those encoding for protein folding genes, heat shock proteins, histones, and Ca2+ ion transport. One of these transcripts, hsc70-4L was found to be under positive selection. We also characterized the functional categories of differentially expressed genes. These candidate genes may be functionally conserved and relevant for related species that will deal with rapid climate change.

Enrichment and identification of differentially expressed genes in hepatocellular carcinoma stem-like cells

Cancer stem cells are considered to be tumor-initiating cells. To explain the initiation or progression of hepatocellular carcinoma (HCC), we previously established a culture system that may enrich hepatic cancer stem-like cells (HCSCs). However, the regulatory mechanisms by which HCSCs acquire stem cell properties remain unclear. In the present study, three pairs of HCSCs and case-matched human HCC cells were analyzed by high-throughput screening, and novel biomarkers and pathways for the regulation of HCSCs were identified. The results led to the identification and stratification of 406 differentially expressed genes (DEGs), among which 73 GO terms were found to be significantly associated with DEGs in HCSCs, and only complement and coagulation cascade pathways were identified during the development of HCSCs. By combining the results of the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, it was revealed that 7 genes were downregulated in the complement and coagulation cascade pathways, and 7 miRNAs were predicted to target several downregulated genes involved in these pathways. The results may contribute toward hepatic cancer stem cell studies and novel drug research for HCC treatment.