An extensive evaluation of read trimming effects on Illumina NGS data analysis

Newly Developed MAGIC Population Allows Identification of Strong Associations and Candidate Genes for Anthocyanin Pigmentation in Eggplant

Multi-parent advanced generation inter-cross (MAGIC) populations facilitate the genetic dissection of complex quantitative traits in plants and are valuable breeding materials. We report the development of the first eggplant MAGIC population (S3 Magic EGGplant InCanum, S3MEGGIC; 8-way), constituted by the 420 S3 individuals developed from the intercrossing of seven cultivated eggplant (Solanum melongena) and one wild relative (S. incanum) parents. The S3MEGGIC recombinant population was genotyped with the eggplant 5k probes SPET platform and phenotyped for anthocyanin presence in vegetative plant tissues (PA) and fruit epidermis (FA), and for the light-insensitive anthocyanic pigmentation under the calyx (PUC). The 7,724 filtered high-confidence single-nucleotide polymorphisms (SNPs) confirmed a low residual heterozygosity (6.87%), a lack of genetic structure in the S3MEGGIC population, and no differentiation among subpopulations carrying a cultivated or wild cytoplasm. Inference of haplotype blocks of the nuclear genome revealed an unbalanced representation of the founder genomes, suggesting a cryptic selection in favour or against specific parental genomes. Genome-wide association study (GWAS) analysis for PA, FA, and PUC detected strong associations with two myeloblastosis (MYB) genes similar to MYB113 involved in the anthocyanin biosynthesis pathway, and with a COP1 gene which encodes for a photo-regulatory protein and may be responsible for the PUC trait. Evidence was found of a duplication of an ancestral MYB113 gene with a translocation from chromosome 10 to chromosome 1 compared with the tomato genome. Parental genotypes for the three genes were in agreement with the identification of the candidate genes performed in the S3MEGGIC population. Our new eggplant MAGIC population is the largest recombinant population in eggplant and is a powerful tool for eggplant genetics and breeding studies.

HiFiAdapterFilt, a memory efficient read processing pipeline, prevents occurrence of adapter sequence in PacBio HiFi reads and their negative impacts on genome assembly

BACKGROUND

Pacific Biosciences HiFi read technology is currently the industry standard for high accuracy long-read sequencing that has been widely adopted by large sequencing and assembly initiatives for generation of de novo assemblies in non-model organisms. Though adapter contamination filtering is routine in traditional short-read analysis pipelines, it has not been widely adopted for HiFi workflows.

RESULTS

Analysis of 55 publicly available HiFi datasets revealed that a read-sanitation step to remove sequence artifacts derived from PacBio library preparation from read pools is necessary as adapter sequences can be erroneously integrated into assemblies.

CONCLUSIONS

Here we describe the nature of adapter contaminated reads, their consequences in assembly, and present HiFiAdapterFilt, a simple and memory efficient solution for removing adapter contaminated reads prior to assembly.

Trimming and Validation of Illumina Short Reads Using Trimmomatic, Trinity Assembly, and Assessment of RNA-Seq Data

Next-generation sequencing (NGS) technologies can generate billions of reads in a single sequencing run. However, with such high-throughput comes quality issues which have to be addressed before undertaking downstream analysis. Quality control on short reads is usually performed at default settings due to a lack of in-depth understanding of a particular software's parameters and their effect if changed on the output. Here we demonstrate how to optimize read trimming using Trimmomatic. We highlight the benefits of trimming by comparing the quality of transcripts assembled using trimmed and untrimmed reads.

Oral microbiome research - A Beginner's glossary

Oral microbiome plays a key role in the etiology of oral diseases and is linked to many diseases in other parts of the body as well. This makes the oral microbiome an area of interest for researchers globally. A meticulous planning of the research project is the first and most crucial step while conducting an oral microbiome study. For beginners in this field, it is essential to be familiar with the terminologies used in oral microbiome research for a better understanding. The purpose of this article is to familiarize new researchers to the frequently used terms for the field of oral microbiome research.

Transcriptomic and metabolic analyses revealed the modulatory effect of vernalization on glucosinolate metabolism in radish (Raphanus sativus L.)

Vernalization is the process by which long-term cold like winter triggers transition to flowering in plants. Many biennial and perennial plants including Brassicaceae family plants require vernalization for floral transition. Not only floral transition, but dynamic physiological and metabolic changes might also take place during vernalization. However, vernalization-mediated metabolic change is merely investigated so far. One of secondary metabolites found in Brassiceceae family plants is glucosinolates (GSLs). GSLs provides defense against pathogens and herbivores attack in plants and also exhibits inhibitory activity against human cancer cell. Profiles of GSLs are highly modulated by different environmental stresses in Brassciaceae family plants. To grasp the effect of vernalization on GSLs metabolic dynamics in radish (Raphanus sativus L.), we performed transcriptomic and metabolic analysis during vernalization in radish. Through transcriptome analysis, we found many GSLs metabolic genes were significantly down-regulated by vernalization in radish plants. Ultra-High Performance Liquid Chromatography analysis also revealed that GSLs compounds were substantially reduced in vernalized radish samples compared to non-vernalized radish samples. Furthermore, we found that repressive histone modification (i.e. H3K27me3) is involved in the modulation of GSLs metabolism via epigenetic suppression of Glucoraphasatin Synthase 1 (GRS1) during vernalization in radish. This study revealed that GSLs metabolism is modulated by vernalization, suggestive of a newly identified target of vernalization in radish.

Evaluating whole-genome sequencing quality metrics for enteric pathogen outbreaks

Background

Whole genome sequencing (WGS) has gained increasing importance in responses to enteric bacterial outbreaks. Common analysis procedures for WGS, single nucleotide polymorphisms (SNPs) and genome assembly, are highly dependent upon WGS data quality.

Methods

Raw, unprocessed WGS reads from Escherichia coli, Salmonella enterica, and Shigella sonnei outbreak clusters were characterized for four quality metrics: PHRED score, read length, library insert size, and ambiguous nucleotide composition. PHRED scores were strongly correlated with improved SNPs analysis results in E. coli and S. enterica clusters.

Results

Assembly quality showed only moderate correlations with PHRED scores and library insert size, and then only for Salmonella. To improve SNP analyses and assemblies, we compared seven read-healing pipelines to improve these four quality metrics and to see how well they improved SNP analysis and genome assembly. The most effective read healing pipelines for SNPs analysis incorporated quality-based trimming, fixed-width trimming, or both. The Lyve-SET SNPs pipeline showed a more marked improvement than the CFSAN SNP Pipeline, but the latter performed better on raw, unhealed reads. For genome assembly, SPAdes enabled significant improvements in healed E. coli reads only, while Skesa yielded no significant improvements on healed reads.

Conclusions

PHRED scores will continue to be a crucial quality metric albeit not of equal impact across all types of analyses for all enteric bacteria. While trimming-based read healing performed well for SNPs analyses, different read healing approaches are likely needed for genome assembly or other, emerging WGS analysis methodologies.

Extensive Genome Exploration of Clostridium botulinum Group III Field Strains

In animals, botulism is commonly sustained by botulinum neurotoxin C, D or their mosaic variants, which are produced by anaerobic bacteria included in Clostridium botulinum group III. In this study, a WGS has been applied to a large collection of C. botulinum group III field strains in order to expand the knowledge on these BoNT-producing Clostridia and to evaluate the potentiality of this method for epidemiological investigations. Sixty field strains were submitted to WGS, and the results were analyzed with respect to epidemiological information and compared to published sequences. The strains were isolated from biological or environmental samples collected in animal botulism outbreaks which occurred in Italy from 2007 to 2016. The new sequenced strains belonged to subspecific groups, some of which were already defined, while others were newly characterized, peculiar to Italian strains and contained genomic features not yet observed. This included, in particular, two new flicC types (VI and VII) and new plasmids which widen the known plasmidome of the species. The extensive genome exploration shown in this study improves the C. botulinum and related species classification scheme, enriching it with new strains of rare genotypes and permitting the highest grade of discrimination among strains for forensic and epidemiological applications.

Genome-wide ultraconserved elements resolve phylogenetic relationships and biogeographic history among Neotropical leaf-nosed bats in the genus Anoura (Phyllostomidae)

AnouraGray, 1838 are Neotropical nectarivorous bats and the most speciose genus within the phyllostomid subfamily Glossophaginae. However, Anoura species limits remain debated, and phylogenetic relationships remain poorly known, because previous studies used limited Anoura taxon sampling or focused primarily on higher-level relationships. Here, we conduct the first phylogenomic study of Anoura by analyzing 2039 genome-wide ultraconserved elements (UCEs) sequenced for 42 individuals from 8 Anoura species/lineages plus two outgroups. Overall, our results based on UCEs resolved relationships in the genus and supported (1) the monophyly of small-bodied Anoura species (previously genus Lonchoglossa); (2) monotypic status of A. caudifer; and (3) nested positions of "A. carishina", A. caudifer aequatoris, and A. geoffroyi peruana specimens within A. latidens, A. caudifer and A. geoffroyi, respectively (suggesting that these taxa are not distinct species). Additionally, (4) phylogenetic networks allowing reticulate edges did not explain gene tree discordance better than the species tree (without introgression), indicating that a coalescent model accounting for discordance solely through incomplete lineage sorting fit our data well. Sensitivity analyses indicated that our species tree results were not adversely affected by varying taxon sampling across loci. Tree calibration and Bayesian coalescent analyses dated the onset of diversification within Anoura to around ∼ 6-9 million years ago in the Miocene, with extant species diverging mainly within the past ∼ 4 million years. We inferred a historical biogeographical scenario for Anoura of parapatric speciation fragmenting the range of a wide-ranging ancestral lineage centered in the Central to Northern Andes, along with Pliocene-Pleistocene dispersal or founder event speciation in Amazonia and the Brazilian Atlantic forest during the last ∼ 2.5 million years.

Microbial Plankton Community Structure and Function Responses to Vitamin B12 and B1 Amendments in an Upwelling System

B vitamins are essential cofactors for practically all living organisms on Earth and are produced by a selection of microorganisms. An imbalance between high demand and limited production, in concert with abiotic processes, may explain the low availability of these vitamins in marine systems. Natural microbial communities from surface shelf water in the productive area off northwestern Spain were enclosed in mesocosms in winter, spring, and summer 2016. In order to explore the impact of B-vitamin availability on microbial community composition (16S and 18S rRNA gene sequence analysis) and bacterial function (metatranscriptomics analysis) in different seasons, enrichment experiments were conducted with seawater from the mesocosms. Our findings revealed that significant increases in phytoplankton or prokaryote biomass associated with vitamin B₁₂ and/or B₁ amendments were not accompanied by significant changes in community composition, suggesting that most of the microbial taxa benefited from the external B-vitamin supply. Metatranscriptome analysis suggested that many bacteria were potential consumers of vitamins B₁₂ and B₁, although the relative abundance of reads related to synthesis was ca. 3.6-fold higher than that related to uptake. Alteromonadales and Oceanospirillales accounted for important portions of vitamin B₁ and B₁₂ synthesis gene transcription, despite accounting for only minor portions of the bacterial community. Flavobacteriales appeared to be involved mostly in vitamin B₁₂ and B₁ uptake, and Pelagibacterales expressed genes involved in vitamin B₁ uptake. Interestingly, the relative expression of vitamin B₁₂ and B₁ synthesis genes among bacteria strongly increased upon inorganic nutrient amendment. Collectively, these findings suggest that upwelling events intermittently occurring during spring and summer in productive ecosystems may ensure an adequate production of these cofactors to sustain high levels of phytoplankton growth and biomass.

IMPORTANCE B vitamins are essential growth factors for practically all living organisms on Earth that are produced by a selection of microorganisms. An imbalance between high demand and limited production may explain the low concentration of these compounds in marine systems. In order to explore the impact of B-vitamin availability on bacteria and algae in the coastal waters off northwestern Spain, six experiments were conducted with natural surface water enclosed in winter, spring, and summer. Our findings revealed that increases in phytoplankton or bacterial growth associated with B₁₂ and/or B₁ amendments were not accompanied by significant changes in community composition, suggesting that most microorganisms benefited from the B-vitamin supply. Our analyses confirmed the role of many bacteria as consumers of vitamins B₁₂ and B₁, although the relative abundance of genes related to synthesis was ca. 3.6-fold higher than that related to uptake. Interestingly, prokaryote expression of B₁₂ and B₁ synthesis genes strongly increased when inorganic nutrients were added. Collectively, these findings suggest that upwelling of cold and nutrient-rich waters occurring during spring and summer in this coastal area may ensure an adequate production of B vitamins to sustain high levels of algae growth and biomass.

A consensus-based ensemble approach to improve transcriptome assembly

BACKGROUND

Systems-level analyses, such as differential gene expression analysis, co-expression analysis, and metabolic pathway reconstruction, depend on the accuracy of the transcriptome. Multiple tools exist to perform transcriptome assembly from RNAseq data. However, assembling high quality transcriptomes is still not a trivial problem. This is especially the case for non-model organisms where adequate reference genomes are often not available. Different methods produce different transcriptome models and there is no easy way to determine which are more accurate. Furthermore, having alternative-splicing events exacerbates such difficult assembly problems. While benchmarking transcriptome assemblies is critical, this is also not trivial due to the general lack of true reference transcriptomes.

RESULTS

In this study, we first provide a pipeline to generate a set of the simulated benchmark transcriptome and corresponding RNAseq data. Using the simulated benchmarking datasets, we compared the performance of various transcriptome assembly approaches including both de novo and genome-guided methods. The results showed that the assembly performance deteriorates significantly when alternative transcripts (isoforms) exist or for genome-guided methods when the reference is not available from the same genome. To improve the transcriptome assembly performance, leveraging the overlapping predictions between different assemblies, we present a new consensus-based ensemble transcriptome assembly approach, ConSemble.

CONCLUSIONS

Without using a reference genome, ConSemble using four de novo assemblers achieved an accuracy up to twice as high as any de novo assemblers we compared. When a reference genome is available, ConSemble using four genome-guided assemblies removed many incorrectly assembled contigs with minimal impact on correctly assembled contigs, achieving higher precision and accuracy than individual genome-guided methods. Furthermore, ConSemble using de novo assemblers matched or exceeded the best performing genome-guided assemblers even when the transcriptomes included isoforms. We thus demonstrated that the ConSemble consensus strategy both for de novo and genome-guided assemblers can improve transcriptome assembly. The RNAseq simulation pipeline, the benchmark transcriptome datasets, and the script to perform the ConSemble assembly are all freely available from: http://bioinfolab.unl.edu/emlab/consemble/ .

Bifidobacteria Strain Typing by Fourier Transform Infrared Spectroscopy

Fourier transform infrared (FTIR) spectroscopy, a technology traditionally used in chemistry to determine the molecular composition of a wide range of sample types, has gained growing interest in microbial typing. It is based on the different vibrational modes of the covalent bonds between atoms of a given sample, as bacterial cells, induced by the absorption of infrared radiation. This technique has been largely used for the study of pathogenic species, especially in the clinical field, and has been proposed also for the typing at different subspecies levels. The high throughput, speed, low cost, and simplicity make FTIR spectroscopy an attractive technique also for industrial applications, in particular, for probiotics. The aim of this study was to compare FTIR spectroscopy with established genotyping methods, pulsed-field gel electrophoresis (PFGE), whole-genome sequencing (WGS), and multilocus sequence typing (MLST), in order to highlight the FTIR spectroscopy potential discriminatory power at strain level. Our study focused on bifidobacteria, an important group of intestinal commensals generally recognized as probiotics. For their properties in promoting and maintaining health, bifidobacteria are largely marketed by the pharmaceutical, food, and dairy industries. Strains belonging to Bifidobacterium longum subsp. longum and Bifidobacterium animalis subsp. lactis were taken into consideration together with some additional type strains. For B. longum subsp. longum, it was possible to discriminate the strains with all the methods used. Although two isolates were shown to be strictly phylogenetically related, constituting a unique cluster, based on PFGE, WGS, and MLST, no clustering was observed with FTIR. For B. animalis subsp. lactis group, PFGE, WGS, and MLST were non-discriminatory, and only one strain was easily distinguished. On the other hand, FTIR discriminated all the isolates one by one, and no clustering was observed. According to these results, FTIR analysis is not only equivalent to PFGE, WGS, and MLST, but also for some strains, in particular, for B. animalis subsp. lactis group, more informative, being able to differentiate strains not discernible with the other two methods based on phenotypic variations likely deriving from certain genetic changes. Fourier transform infrared spectroscopy has highlighted the possibility of using the cell surface as a kind of barcode making tracing strains possible, representing an important aspect in probiotic applications. Furthermore, this work constitutes the first investigation on bifidobacterial strain typing using FTIR spectroscopy.

Paraphyly and low levels of genetic divergence in morphologically distinct taxa: revision of the Pseudoanthidium scapulare complex of carder bees (Apoidea: Megachilidae: Anthidiini)

The Palaearctic complex of anthidiine bees closely related to Pseudoanthidium scapulare has long been a source of unresolved taxonomic and systematic issues. Until now, the number of species in the complex and their geographical distributions were largely unclear, thus complicating the compilation of accurate species checklists and hindering conservation efforts. In order to address these issues, we use morphology and mitochondrial cytochrome c oxidase subunit I (COI) sequences, combined with a thorough examination of the relevant literature and type material, to delimit species within this complex, assign names to species and clarify geographical ranges. An unexpected result was that a certain number of morphologically distinct taxa exhibited low levels of genetic divergence at the COI locus, resulting in species paraphyly. A set of ultra-conserved elements (UCEs) was also sequenced in order to further investigate relationships among these taxa. One morphologically distinct species was also paraphyletic using UCE data, hinting at recent species divergences and genetic exchange at zones of contact between morphologically well-differentiated taxa. The results of our study reveal the presence of ten species in this complex, including a previously overlooked species for western continental Europe. A complete diagnosis of the males and females of these species is provided, as are maps detailing the geographic distributions of each. An illustrated identification key to the males and females of each species is presented. Two new species are described, Pseudoanthidium kaspareki sp. nov. and P. rozeni sp. nov. New synonymy is established for several species and Pseudoanthidium palestinicum and P. tropicum are raised to species level. The new combination, Icteranthidium floripetum comb. nov. is also established. Lectotypes are designated for the following species: Anthidium eversmanni, A. floripetum, A. frontale, A. karakalense, A. nanum and A. reptans. Previously unpublished lectotype designations are published here for A. sinuatum and A. tenellum.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

Analysis of Four New Enterococcus faecalis Phages and Modeling of a Hyaluronidase Catalytic Domain from Saphexavirus

Background: Phage therapy (PT), as a method to treat bacterial infections, needs identification of bacteriophages targeting specific pathogenic host. Enterococcus faecalis, a Gram-positive coccus resident in the human gastrointestinal tract, may become pathogenic in hospitalized patients showing acquired resistance to vancomycin and thus representing a possible target for PT. Materials and Methods: We isolated four phages that infect E. faecalis and characterized them by host range screening, transmission electron microscopy, and genome sequencing. We also identified and three-dimensional modeled a new hyaluronidase enzyme. Results: The four phages belong to Siphoviridae family: three Efquatrovirus (namely vB_EfaS_TV51, vB_EfaS_TV54, and vB_EfaS_TV217) and one Saphexavirus (vB_EfaS_TV16). All of them are compatible with lytic cycle. vB_EfaS_TV16 moreover presents a gene encoding for a hyaluronidase enzyme. Conclusions: The identified phages show features suggesting their useful application in PT, particularly the Saphexavirus that may be of enhanced relevance in PT because of its potential biofilm-digestion capability.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L

BACKGROUND

Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst.

RESULTS

We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech.

CONCLUSIONS

Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

Effects of Dihydrotanshinone I on Proliferation and Invasiveness of Paclitaxel-Resistant Anaplastic Thyroid Cancer Cells

ATC is a very rare, but extremely aggressive form of thyroid malignancy, responsible for the highest mortality rate registered for thyroid cancer. In patients without known genetic aberrations, the current treatment is still represented by palliative surgery and systemic mono- or combined chemotherapy, which is often not fully effective for the appearance of drug resistance. Comprehension of the mechanisms involved in the development of the resistance is therefore an urgent issue to suggest novel therapeutic approaches for this very aggressive malignancy. In this study, we created a model of anaplastic thyroid cancer (ATC) cells resistant to paclitaxel and investigated the characteristics of these cells by analyzing the profile of gene expression and comparing it with that of paclitaxel-sensitive original ATC cell lines. In addition, we evaluated the effects of Dihydrotanshinone I (DHT) on the viability and invasiveness of paclitaxel-resistant cells. ATC paclitaxel-resistant cells highlighted an overexpression of ABCB1 and a hyper-activation of the NF-κB compared to sensitive cells. DHT treatment resulted in a reduction of viability and clonogenic ability of resistant cells. Moreover, DHT induces a decrement of NF-κB activity in SW1736-PTX and 8505C-PTX cells. In conclusion, to the best of our knowledge, the results of the present study are the first to demonstrate the antitumor effects of DHT on ATC cells resistant to Paclitaxel in vitro.

Dihydrotanshinone exerts antitumor effects and improves the effects of cisplatin in anaplastic thyroid cancer cells

Anaplastic thyroid cancer (ATC) is the most aggressive type of thyroid cancer and is responsible for 20‑50% of thyroid cancer‑associated deaths. The absence of response to conventional treatments makes the search for novel therapeutics a clinical challenge. In the present study, the effects of 15,16‑dihydrotanshinone I (DHT), a tanshinone extracted from Salvia miltiorrhiza Bunge (Danshen), which has previously been shown to possess anticancer activity, were examined in two human ATC cell lines. DHT significantly reduced cell viability, which was coupled with an increase in apoptosis. DHT administration also reduced the colony‑forming ability and proliferation of these cells in soft agar and downregulated the expression of epithelial‑to‑mesenchymal transition‑related genes. In addition, DHT significantly reduced MAD2 expression, a target of HuR with a relevant role in ATC. Finally, cotreatment with cisplatin and DHT has a greater effect on cell viability than each compound alone. In conclusion, to the best of our knowledge, the present study is the first to demonstrate that DHT exerts antitumor effects on ATC cells by reducing MAD2 expression levels. Moreover, a synergistic effect of DHT with cisplatin was shown. Further in vivo studies are required to assess this phytochemical compound as a potential adjuvant for the treatment of ATC.

Data on the first functionally-annotated de novo transcriptome assembly for North American flying squirrels (genus Glaucomys)

We report the first functionally-annotated de novo transcriptome assembly for North American flying squirrels (genus Glaucomys). RNA was extracted from tissue samples obtained from two northern flying squirrels and two southern flying squirrels sampled from Ontario, Canada, and sequenced on an Illumina paired-end sequencing platform. We reconstructed 702,228 Glaucomys transcripts using 193,323,120 sequence read pairs and captured sequence homologies, protein domains, and gene function classifications. Introgressive hybridization between northern (Glaucomys sabrinus) and southern flying squirrels (G. volans) has been observed in some areas of North America. However, existing molecular markers lack the resolution to discriminate late-generation introgressants and describe the extent to which hybridization influences the Glaucomys gene pool. These genomic resources can increase the resolution of molecular techniques used to examine the dynamics of the Glaucomys hybrid zone.

Transcriptome and proteome of the corm, leaf and flower of Hypoxis hemerocallidea (African potato)

The corm of Hypoxis hemerocallidea, commonly known as the African potato, is used in traditional medicine to treat several medical conditions such as urinary infections, benign prostate hyperplasia, inflammatory conditions and testicular tumours. The metabolites contributing to the medicinal properties of H. hemerocallidea have been identified in several studies and, more recently, the active terpenoids of the plant were profiled. However, the biosynthetic pathways and the enzymes involved in the production of the terpene metabolites in H. hemerocallidea have not been characterised at a transcriptomic or proteomic level. In this study, total RNA extracted from the corm, leaf and flower tissues of H. hemerocallidea was sequenced on the Illumina HiSeq 2500 platform. A total of 143,549 transcripts were assembled de novo using Trinity and 107,131 transcripts were functionally annotated using the nr, GO, COG, KEGG and SWISS-PROT databases. Additionally, the proteome of the three tissues were sequenced using LC-MS/MS, revealing aspects of secondary metabolism and serving as data validation for the transcriptome. Functional annotation led to the identification of numerous terpene synthases such as nerolidol synthase, germacrene D synthase, and cycloartenol synthase amongst others. Annotations also revealed a transcript encoding the terpene synthase phytoalexin momilactone A synthase. Differential expression analysis using edgeR identified 946 transcripts differentially expressed between the three tissues and revealed that the leaf upregulates linalool synthase compared to the corm and the flower tissues. The transcriptome as well as the proteome of Hypoxis hemerocallidea presented here provide a foundation for future research.

Responses of Coastal Marine Microbiomes Exposed to Anthropogenic Dissolved Organic Carbon

Coastal seawaters receive thousands of organic pollutants. However, we have little understanding of the response of microbiomes to this pool of anthropogenic dissolved organic carbon (ADOC). In this study, coastal microbial communities were challenged with ADOC at environmentally relevant concentrations. Experiments were performed at two Mediterranean sites with different impact by pollutants and nutrients: off the Barcelona harbor ("BCN"), and at the Blanes Bay ("BL"). ADOC additions stimulated prokaryotic leucine incorporation rates at both sites, indicating the use of ADOC as growth substrate. The percentage of "membrane-compromised" cells increased with increasing ADOC, indicating concurrent toxic effects of ADOC. Metagenomic analysis of the BCN community challenged with ADOC showed a significant growth of Methylophaga and other gammaproteobacterial taxa belonging to the rare biosphere. Gene expression profiles showed a taxon-dependent response, with significantly enrichments of transcripts from SAR11 and Glaciecola spp. in BCN and BL, respectively. Further, the relative abundance of transposon-related genes (in BCN) and transcripts (in BL) correlated with the number of differentially abundant genes (in BCN) and transcripts (in BLA), suggesting that microbial responses to pollution may be related to pre-exposure to pollutants, with transposons playing a role in adaptation to ADOC. Our results point to a taxon-specific response to low concentrations of ADOC that impact the functionality, structure and plasticity of the communities in coastal seawaters. This work contributes to address the influence of pollutants on microbiomes and their perturbation to ecosystem services and ocean health.

A Large-Scale and Serverless Computational Approach for Improving Quality of NGS Data Supporting Big Multi-Omics Data Analyses

Various types of analyses performed over multi-omics data are driven today by next-generation sequencing (NGS) techniques that produce large volumes of DNA/RNA sequences. Although many tools allow for parallel processing of NGS data in a Big Data distributed environment, they do not facilitate the improvement of the quality of NGS data for a large scale in a simple declarative manner. Meanwhile, large sequencing projects and routine DNA/RNA sequencing associated with molecular profiling of diseases for personalized treatment require both good quality data and appropriate infrastructure for efficient storing and processing of the data. To solve the problems, we adapt the concept of Data Lake for storing and processing big NGS data. We also propose a dedicated library that allows cleaning the DNA/RNA sequences obtained with single-read and paired-end sequencing techniques. To accommodate the growth of NGS data, our solution is largely scalable on the Cloud and may rapidly and flexibly adjust to the amount of data that should be processed. Moreover, to simplify the utilization of the data cleaning methods and implementation of other phases of data analysis workflows, our library extends the declarative U-SQL query language providing a set of capabilities for data extraction, processing, and storing. The results of our experiments prove that the whole solution supports requirements for ample storage and highly parallel, scalable processing that accompanies NGS-based multi-omics data analyses.

Genome wide association study of agronomic and seed traits in a world collection of proso millet (Panicum miliaceum L.)

BACKGROUND

The climate crisis threatens sustainability of crop production worldwide. Crop diversification may enhance food security while reducing the negative impacts of climate change. Proso millet (Panicum milaceum L.) is a minor cereal crop which holds potential for diversification and adaptation to different environmental conditions. In this study, we assembled a world collection of proso millet consisting of 88 varieties and landraces to investigate its genomic and phenotypic diversity for seed traits, and to identify marker-trait associations (MTA).

RESULTS

Sequencing of restriction-site associated DNA fragments yielded 494 million reads and 2,412 high quality single nucleotide polymorphisms (SNPs). SNPs were used to study the diversity in the collection and perform a genome wide association study (GWAS). A genotypic diversity analysis separated accessions originating in Western Europe, Eastern Asia and Americas from accessions sampled in Southern Asia, Western Asia, and Africa. A Bayesian structure analysis reported four cryptic genetic groups, showing that landraces accessions had a significant level of admixture and that most of the improved proso millet materials clustered separately from landraces. The collection was highly diverse for seed traits, with color varying from white to dark brown and width spanning from 1.8 to 2.6 mm. A GWAS study for seed morphology traits identified 10 MTAs. In addition, we identified three MTAs for agronomic traits that were previously measured on the collection.

CONCLUSION

Using genomics and automated seed phenotyping, we elucidated phylogenetic relationships and seed diversity in a global millet collection. Overall, we identified 13 MTAs for key agronomic and seed traits indicating the presence of alleles with potential for application in proso breeding programs.

Genomic Characterization of hlyF-positive Shiga Toxin-Producing Escherichia coli, Italy and the Netherlands, 2000-2019

Shiga toxin–producing Escherichia coli (STEC) O80:H2 has emerged in Europe as a cause of hemolytic uremic syndrome associated with bacteremia. STEC O80:H2 harbors the mosaic plasmid pR444_A, which combines several virulence genes, including hlyF and antimicrobial resistance genes. pR444_A is found in some extraintestinal pathogenic E. coli (ExPEC) strains. We identified and characterized 53 STEC strains with ExPEC-associated virulence genes isolated in Italy and the Netherlands during 2000–2019. The isolates belong to 2 major populations: 1 belongs to sequence type 301 and harbors diverse stx₂ subtypes, the intimin variant eae-ξ, and pO157-like and pR444_A plasmids; 1 consists of strains belonging to various sequence types, some of which lack the pO157 plasmid, the locus of enterocyte effacement, and the antimicrobial resistance–encoding region. Our results showed that STEC strains harboring ExPEC-associated virulence genes can include multiple serotypes and that the pR444_A plasmid can be acquired and mobilized by STEC strains.

A transcriptomic, metabolomic and cellular approach to the physiological adaptation of tomato fruit to high temperature

High temperatures can negatively influence plant growth and development. Besides yield, the effects of heat stress on fruit quality traits remain poorly characterised. In tomato, insights into how fruits regulate cellular metabolism in response to heat stress could contribute to the development of heat-tolerant varieties, without detrimental effects on quality. In the present study, the changes occurring in wild type tomato fruits after exposure to transient heat stress have been elucidated at the transcriptome, cellular and metabolite level. An impact on fruit quality was evident as nutritional attributes changed in response to heat stress. Fruit carotenogenesis was affected, predominantly at the stage of phytoene formation, although altered desaturation/isomerisation arose during the transient exposure to high temperatures. Plastidial isoprenoid compounds showed subtle alterations in their distribution within chromoplast sub-compartments. Metabolite profiling suggests limited effects on primary/intermediary metabolism but lipid remodelling was evident. The heat-induced molecular signatures included the accumulation of sucrose and triacylglycerols, and a decrease in the degree of membrane lipid unsaturation, which influenced the volatile profile. Collectively, these data provide valuable insights into the underlying biochemical and molecular adaptation of fruit to heat stress and will impact on our ability to develop future climate resilient tomato varieties.

Bioinformatics tools used for whole-genome sequencing analysis of Neisseria gonorrhoeae: a literature review

Whole-genome sequencing (WGS) data are well established for the investigation of gonococcal transmission, antimicrobial resistance prediction, population structure determination and population dynamics. A variety of bioinformatics tools, repositories, services and platforms have been applied to manage and analyze Neisseria gonorrhoeae WGS datasets. This review provides an overview of the various bioinformatics approaches and resources used in 105 published studies (as of 30 April 2021). The challenges in the analysis of N. gonorrhoeae WGS datasets, as well as future bioinformatics requirements, are also discussed.

The probiotic Lactobacillus rhamnosus mimics the dark-driven regulation of appetite markers and melatonin receptors' expression in zebrafish (Danio rerio) larvae: Understanding the role of the gut microbiome

The use of probiotics has been recently considered a novel therapeutic strategy to prevent pathologies such as obesity; however, the specific mechanisms of action by which probiotics exert their beneficial effects on metabolic health remain unclear. The aim of the present study was to investigate the short-term effects of a probiotic Lactobacillus rhamnosus supplementation (PROB) on appetite regulation, growth-related markers, and microbiota diversity in zebrafish (Danio rerio) larvae, compared to a group subjected to a constant darkness photoperiod (DARK), as well as to evaluate the effects of both treatments on melatonin receptors' expression. After a 24 h treatment, both PROB and DARK conditions caused a significant increase in leptin a expression. Moreover, mRNA abundances of leptin b and proopiomelanocortin a were elevated in the PROB group, and DARK showed a similar tendency, supporting a negative regulation of appetite markers by the treatments. Moreover, both PROB and DARK also enhanced the abundances of melatonin receptors transcript (melatonin receptor 1 ba and bb) and protein (melatonin receptor 1) suggesting a potential involvement of melatonin in mediating these effects. Nevertheless, treatments did not exhibit a significant effect on the expression of most of the growth hormone/insulin-like growth factor axis genes evaluated. Finally, only the DARK condition significantly modulated gut microbiota diversity at such short time, altogether highlighting the rapid effects of this probiotic on modulating appetite regulatory and melatonin receptors' expression, without a concomitant variation of gut microbiota.

Lack of parent-of-origin effects in Nasonia jewel wasp: A replication and extension study

In diploid cells, the paternal and maternal alleles are, on average, equally expressed. There are exceptions from this: a small number of genes express the maternal or paternal allele copy exclusively. This phenomenon, known as genomic imprinting, is common among eutherian mammals and some plant species; however, genomic imprinting in species with haplodiploid sex determination is not well characterized. Previous work reported no parent-of-origin effects in the hybrids of closely related haplodiploid Nasonia vitripennis and Nasonia giraulti jewel wasps, suggesting a lack of epigenetic reprogramming during embryogenesis in these species. Here, we replicate the gene expression dataset and observations using different individuals and sequencing technology, as well as reproduce these findings using the previously published RNA sequence data following our data analysis strategy. The major difference from the previous dataset is that they used an introgression strain as one of the parents and we found several loci that resisted introgression in that strain. Our results from both datasets demonstrate a species-of-origin effect, rather than a parent-of-origin effect. We present a reproducible workflow that others may use for replicating the results. Overall, we reproduced the original report of no parent-of-origin effects in the haplodiploid Nasonia using the original data with our new processing and analysis pipeline and replicated these results with our newly generated data.

Prevalence, risk factors and genetic traits of Salmonella Infantis in Dutch broiler flocks

Salmonella Infantis is a poultry-adapted Salmonella enterica serovar that is increasingly reported in broilers and is also regularly identified among human salmonellosis cases. An emerging S. Infantis mega-plasmid (pESI), carrying fitness, virulence and antimicrobial resistance genes, is also increasingly found. We investigated the prevalence, genetic characteristics and risk factors for (pESI-carrying) S. Infantis in broilers. Faecal samples from 379 broiler flocks (in 198 farms with ≥3000 birds) in the Netherlands were tested. A questionnaire about farm characteristics was also administered. Sampling was performed in July 2018-May 2019, three weeks before slaughter. Fourteen flocks (in 10 farms) were S. Infantis-positive, resulting in a 3.7 % flock-level and 5.1 % farm-level prevalence. Based on multi-locus sequence typing (MLST), all isolates belonged to sequence type 32. All but one isolate carried a pESI-like mega-plasmid. Core-genome MLST showed considerable heterogeneity among the isolates, even within the same farm, with a few small clusters detected. The typical pESI-borne multi-resistance pattern to aminoglycosides, sulphonamide and tetracycline (93 %), as well as trimethoprim (71 %), was found. Additionally, resistance to (fluoro)quinolones based on gyrA gene mutations was detected. S. Infantis was found more often in flocks using salinomycin as coccidiostat, where flock thinning was applied or litter quality was poor, whereas employing external cleaning companies, wheat in feed, and vaccination against infectious bronchitis, were protective. Suggestive evidence for vertical transmission from hatcheries was found. A heterogeneous (pESI-carrying) S. Infantis population has established itself in Dutch broiler flocks, calling for further monitoring of its spread and a comprehensive appraisal of control options.

Performance evaluation of pipelines for mapping, variant calling and interval padding, for the analysis of NGS germline panels

Background

Next-generation sequencing (NGS) represents a significant advancement in clinical genetics. However, its use creates several technical, data interpretation and management challenges. It is essential to follow a consistent data analysis pipeline to achieve the highest possible accuracy and avoid false variant calls. Herein, we aimed to compare the performance of twenty-eight combinations of NGS data analysis pipeline compartments, including short-read mapping (BWA-MEM, Bowtie2, Stampy), variant calling (GATK-HaplotypeCaller, GATK-UnifiedGenotyper, SAMtools) and interval padding (null, 50 bp, 100 bp) methods, along with a commercially available pipeline (BWA Enrichment, Illumina®). Fourteen germline DNA samples from breast cancer patients were sequenced using a targeted NGS panel approach and subjected to data analysis.

Results

We highlight that interval padding is required for the accurate detection of intronic variants including spliceogenic pathogenic variants (PVs). In addition, using nearly default parameters, the BWA Enrichment algorithm, failed to detect these spliceogenic PVs and a missense PV in the TP53 gene. We also recommend the BWA-MEM algorithm for sequence alignment, whereas variant calling should be performed using a combination of variant calling algorithms; GATK-HaplotypeCaller and SAMtools for the accurate detection of insertions/deletions and GATK-UnifiedGenotyper for the efficient detection of single nucleotide variant calls.

Conclusions

These findings have important implications towards the identification of clinically actionable variants through panel testing in a clinical laboratory setting, when dedicated bioinformatics personnel might not always be available. The results also reveal the necessity of improving the existing tools and/or at the same time developing new pipelines to generate more reliable and more consistent data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04144-1.

Weighted gene co expression network analysis (WGCNA) with key pathways and hub-genes related to micro RNAs in ischemic stroke

Ischemic stroke (IS) is one of the major causes of death and disability worldwide. However, the specific mechanism of gene interplay and the biological function in IS are not clear. Therefore, more research into IS is necessary. Dataset GSE110993 including 20 ischemic stroke (IS) and 20 control specimens are used to establish both groups and the raw RNA‐seq data were analysed. Weighted gene co‐expression network analysis (WGCNA) was used to screen the key micro‐RNA modules. The centrality of key genes were determined by module membership (mm) and gene significance (GS). The key pathways were identified by enrichment analysis with Kyoto Protocol Gene and Genome Encyclopedia (KEGG), and the key genes were validated by protein‐protein interactions network. Result: Upon investigation, 1185 up‐ and down‐regulated genes were gathered and distributed into three modules in response to their degree of correlation to clinical traits of IS, among which the turquoise module show a trait‐correlation of 0.77. The top 140 genes were further identified by GS and MM. KEGG analysis showed two pathways may evolve in the progress of IS. Discussion: CXCL12 and EIF2a may be important biomarkers for the accurate diagnosis and treatment in IS.

ngsComposer: an automated pipeline for empirically based NGS data quality filtering

Next-generation sequencing (NGS) enables massively parallel acquisition of large-scale omics data; however, objective data quality filtering parameters are lacking. Although a useful metric, evidence reveals that platform-generated Phred values overestimate per-base quality scores. We have developed novel and empirically based algorithms that streamline NGS data quality filtering. The pipeline leverages known sequence motifs to enable empirical estimation of error rates, detection of erroneous base calls and removal of contaminating adapter sequence. The performance of motif-based error detection and quality filtering were further validated with read compression rates as an unbiased metric. Elevated error rates at read ends, where known motifs lie, tracked with propagation of erroneous base calls. Barcode swapping, an inherent problem with pooled libraries, was also effectively mitigated. The ngsComposer pipeline is suitable for various NGS protocols and platforms due to the universal concepts on which the algorithms are based.

Demography, genetic diversity and expansion load in the colonizing species Leontodon longirostris (Asteraceae) throughout its native range

Unravelling the evolutionary processes underlying range expansions is fundamental to understand the distribution of organisms, as well as to predict their future responses to environmental change. Predictions for range expansions include a loss of genetic diversity and an accumulation of deleterious alleles along the expansion axis, which can decrease fitness at the range-front (expansion load). In plants, empirical studies supporting expansion load are scarce, and its effects remain to be tested outside a few model species. Leontodon longirostris is a colonizing Asteraceae with a widespread distribution in the Western Mediterranean, providing a particularly interesting system to gain insight into the factors that can enhance or mitigate expansion load. In this study, we produced a first genome draft for the species, covering 418 Mbp (~53% of the genome). Although incomplete, this draft was suitable to design a targeted sequencing of ~1.5 Mbp in 238 L. longirostris plants from 21 populations distributed along putative colonization routes in the Iberian Peninsula. Inferred demographic history supports a range expansion from southern Iberia around 40,000 years ago, reaching northern Iberia around 25,000 years ago. The expansion was accompanied by a loss of genetic diversity and a significant increase in the proportion of putatively deleterious mutations. However, levels of expansion load in L. longirostris were smaller than those found in other plant species, which can be explained, at least partially, by its high dispersal ability, the self-incompatible mating system, and the fact that the expansion occurred along a strong environmental cline.

A roadmap for the generation of benchmarking resources for antimicrobial resistance detection using next generation sequencing

Next Generation Sequencing technologies significantly impact the field of Antimicrobial Resistance (AMR) detection and monitoring, with immediate uses in diagnosis and risk assessment. For this application and in general, considerable challenges remain in demonstrating sufficient trust to act upon the meaningful information produced from raw data, partly because of the reliance on bioinformatics pipelines, which can produce different results and therefore lead to different interpretations. With the constant evolution of the field, it is difficult to identify, harmonise and recommend specific methods for large-scale implementations over time. In this article, we propose to address this challenge through establishing a transparent, performance-based, evaluation approach to provide flexibility in the bioinformatics tools of choice, while demonstrating proficiency in meeting common performance standards. The approach is two-fold: first, a community-driven effort to establish and maintain "live" (dynamic) benchmarking platforms to provide relevant performance metrics, based on different use-cases, that would evolve together with the AMR field; second, agreed and defined datasets to allow the pipelines' implementation, validation, and quality-control over time. Following previous discussions on the main challenges linked to this approach, we provide concrete recommendations and future steps, related to different aspects of the design of benchmarks, such as the selection and the characteristics of the datasets (quality, choice of pathogens and resistances, etc.), the evaluation criteria of the pipelines, and the way these resources should be deployed in the community.

Caloric Restriction Promotes Immunometabolic Reprogramming Leading to Protection from Tuberculosis

There is a strong relationship between metabolic state and susceptibility to Mycobacterium tuberculosis (MTB) infection, with energy metabolism setting the basis for an exaggerated immuno-inflammatory response, which concurs with MTB pathogenesis. Herein, we show that controlled caloric restriction (CR), not leading to malnutrition, protects susceptible DBA/2 mice against pulmonary MTB infection by reducing bacterial load, lung immunopathology, and generation of foam cells, an MTB reservoir in lung granulomas. Mechanistically, CR induced a metabolic shift toward glycolysis, and decreased both fatty acid oxidation and mTOR activity associated with induction of autophagy in immune cells. An integrated multi-omics approach revealed a specific CR-induced metabolomic, transcriptomic, and proteomic signature leading to reduced lung damage and protective remodeling of lung interstitial tightness able to limit MTB spreading. Our data propose CR as a feasible immunometabolic manipulation to control MTB infection, and this approach offers an unexpected strategy to boost immunity against MTB.

EFECTO DEL FILTRADO DE SECUENCIAS EN EL ENSAMBLADO DEL GENOMA DE Bacillus altitudinis AISLADO DE Ilex paraguariensis

Sin importar el tipo de tecnología aplicada para la secuenciación de un genoma, el filtrado de secuencias es un paso esencial, en el cual aquellas lecturas de baja calidad o parte de estas son eliminadas. En un ensamblado la construcción de un genoma se realiza a partir de la unión de lecturas cortas en cóntigos. Algunos ensambladores miden la relación que existe entre secuencias de una longitud fija (k-mer) que puede verse afectada por la presencia de secuencias de baja calidad. Un enfoque común para evaluar los ensamblados se basa en el análisis del número de cóntigos, la longitud del cóntigo más largo y el valor de N50, definido como la longitud del cóntigo que representa el 50 % de la longitud del conjunto. En este contexto, el presente estudio tuvo como objetivo evaluar el efecto del uso de lecturas crudas y filtradas en los valores de los parámetros de calidad obtenidos en el ensamblado del genoma de la cepa de Bacillus altitudinis19RS3 aislada de Ilex paraguariensis. Se realizó el análisis de calidad de ambos archivos de partida con el softwareFastqC y se filtraron las lecturas con el softwareTrimmomatic. Para el ensamblado se utilizó el softwareSPAdes y para su evaluación la herramienta QUAST. El mejor ensamblado para B. altitudinis19RS3 se obtuvo a partir de las lecturasfiltradas con el valor dek-mer 79, que generó 16 cóntigos mayores a 500 pb con un N50 de 931 914 pb y el cóntigo más largo de 966 271 pb.

Genome Recovery, Functional Profiling, and Taxonomic Classification from Metagenomes

Recovering and annotating bacterial genomes from metagenomes involves a series of complex computational tools that are often difficult to use for researches without a specialistic bioinformatic background. In this chapter we review all the steps that lead from raw reads to a collection of quality-controlled, functionally annotated bacterial genomes and propose a working protocol using state-of-the-art, open source software tools.

Targeted Mutagenesis of the Female-Suppressor SyGI Gene in Tetraploid Kiwifruit by CRISPR/CAS9

Kiwifruit belong to the genus Actinidia with 54 species apparently all functionally dioecious. The sex-determinants of the type XX/XY, with male heterogametic, operate independently of the ploidy level. Recently, the SyGI protein has been described as the suppressor of female development. In the present study, we exploited the CRISPR/Cas9 technology by targeting two different sites in the SyGI gene in order to induce a stable gene knock-out in two tetraploid male accessions of Actinidia chinensis var. chinensis. The two genotypes showed a regenerative efficiency of 58% and 73%, respectively. Despite not yet being able to verify the phenotypic effects on the flower structure, due to the long time required by tissue-cultured kiwifruit plants to flower, we obtained two regenerated lines showing near fixation of a unique modification in their genome, resulting in both cases in the onset of a premature stop codon, which induces the putative gene knock-out. Evaluation of gRNA1 locus for both regenerated plantlets resulted in co-amplification of a minor variant differing from the target region for a single nucleotide. A genomic duplication of the region in proximity of the Y genomic region could be postulated.

Read trimming has minimal effect on bacterial SNP-calling accuracy

Read alignment is the central step of many analytic pipelines that perform variant calling. To reduce error, it is common practice to pre-process raw sequencing reads to remove low-quality bases and residual adapter contamination, a procedure collectively known as ‘trimming’. Trimming is widely assumed to increase the accuracy of variant calling, although there are relatively few systematic evaluations of its effects and no clear consensus on its efficacy. As sequencing datasets increase both in number and size, it is worthwhile reappraising computational operations of ambiguous benefit, particularly when the scope of many analyses now routinely incorporates thousands of samples, increasing the time and cost required. Using a curated set of 17 Gram-negative bacterial genomes, this study initially evaluated the impact of four read-trimming utilities (Atropos, fastp, Trim Galore and Trimmomatic), each used with a range of stringencies, on the accuracy and completeness of three bacterial SNP-calling pipelines. It was found that read trimming made only small, and statistically insignificant, increases in SNP-calling accuracy even when using the highest-performing pre-processor in this study, fastp. To extend these findings, >6500 publicly archived sequencing datasets from Escherichia coli, Mycobacterium tuberculosis and Staphylococcus aureus were re-analysed using a common analytic pipeline. Of the approximately 125 million SNPs and 1.25 million indels called across all samples, the same bases were called in 98.8 and 91.9 % of cases, respectively, irrespective of whether raw reads or trimmed reads were used. Nevertheless, the proportion of mixed calls (i.e. calls where <100 % of the reads support the variant allele; considered a proxy of false positives) was significantly reduced after trimming, which suggests that while trimming rarely alters the set of variant bases, it can affect the proportion of reads supporting each call. It was concluded that read quality- and adapter-trimming add relatively little value to a SNP-calling pipeline and may only be necessary if small differences in the absolute number of SNP calls, or the false call rate, are critical. Broadly similar conclusions can be drawn about the utility of trimming to an indel-calling pipeline. Read trimming remains routinely performed prior to variant calling likely out of concern that doing otherwise would typically have negative consequences. While historically this may have been the case, the data in this study suggests that read trimming is not always a practical necessity.

Inter- and intra-archipelago dynamics of population structure and gene flow in a Polynesian bird

Islands are separated by natural barriers that prevent gene flow between terrestrial populations and promote allopatric diversification. Birds in the South Pacific are an excellent model to explore the interplay between isolation and gene flow due to the region's numerous archipelagos and well-characterized avian communities. The wattled honeyeater complex (Foulehaio spp.) comprises three allopatric species that are widespread and common across Fiji, Tonga, Samoa, and Wallis and Futuna. Here, we explored patterns of diversification within and among these lineages using genomic and morphometric data. We found support for three clades of Foulehaio corresponding to three recognized species. Within F. carunculatus, population genetic analyses identified nine major lineages, most of which were composed of sub-lineages that aligned nearly perfectly to individual island populations. Despite genetic structure and great geographic distance between populations, we found low levels of gene flow between populations in adjacent archipelagos. Additionally, body size of F. carunculatus varied randomly with respect to evolutionary history (as Ernst Mayr predicted), but correlated negatively with island size, consistent with the island rule. Our findings support a hypothesis that widespread taxa can show population structure between immediately adjacent islands, and likely represent many independent lineages loosely connected by gene flow.

Role of DNA methylation in the chromium tolerance of Scenedesmus acutus (Chlorophyceae) and its impact on the sulfate pathway regulation

DNA methylation is a very important epigenetic modification that participates in many biological functions. Although many researches on DNA methylation have been reported in various plant species, few studies have assessed the global DNA methylation pattern in algae. Even more the complex mechanisms by which DNA methylation modulates stress in algae are yet largely unresolved, mainly with respect to heavy metal stress, for which in plants, metal- and species- specific responses were instead evidenced. In this work, we performed a comparative Whole-Genome Bisulfite Sequencing (WGBS) on two strains of the green alga Scenedesmus acutus with different Cr(VI) sensitivity. The pattern of distribution of 5-mC showed significant differences between the two strains concerning both differentially methylated local contexts (CG, CHG and CHH) and Differentially Methylated Regions (DMRs) as well. We also demonstrated that DNA methylation plays an important role in modulating some genes for sulfate uptake/assimilation confirming the involvement of the sulfate pathway in the Cr-tolerance. Our results suggest that DNA methylation may be of particular importance in defining signal specificity associated with Cr-tolerance and in establishing new epigenetic marks which contribute to the adaptation to metal stress and also to transmit the epigenomic traits to the progeny.

Increased biological relevance of transcriptome analyses in human skeletal muscle using a model-specific pipeline

Background

Human skeletal muscle responds to weight-bearing exercise with significant inter-individual differences. Investigation of transcriptome responses could improve our understanding of this variation. However, this requires bioinformatic pipelines to be established and evaluated in study-specific contexts. Skeletal muscle subjected to mechanical stress, such as through resistance training (RT), accumulates RNA due to increased ribosomal biogenesis. When a fixed amount of total-RNA is used for RNA-seq library preparations, mRNA counts are thus assessed in different amounts of tissue, potentially invalidating subsequent conclusions. The purpose of this study was to establish a bioinformatic pipeline specific for analysis of RNA-seq data from skeletal muscles, to explore the effects of different normalization strategies and to identify genes responding to RT in a volume-dependent manner (moderate vs. low volume). To this end, we analyzed RNA-seq data derived from a twelve-week RT intervention, wherein 25 participants performed both low- and moderate-volume leg RT, allocated to the two legs in a randomized manner. Bilateral muscle biopsies were sampled from m. vastus lateralis before and after the intervention, as well as before and after the fifth training session (Week 2).

Result

Bioinformatic tools were selected based on read quality, observed gene counts, methodological variation between paired observations, and correlations between mRNA abundance and protein expression of myosin heavy chain family proteins. Different normalization strategies were compared to account for global changes in RNA to tissue ratio. After accounting for the amounts of muscle tissue used in library preparation, global mRNA expression increased by 43–53%. At Week 2, this was accompanied by dose-dependent increases for 21 genes in rested-state muscle, most of which were related to the extracellular matrix. In contrast, at Week 12, no readily explainable dose-dependencies were observed. Instead, traditional normalization and non-normalized models resulted in counterintuitive reverse dose-dependency for many genes. Overall, training led to robust transcriptome changes, with the number of differentially expressed genes ranging from 603 to 5110, varying with time point and normalization strategy.

Conclusion

Optimized selection of bioinformatic tools increases the biological relevance of transcriptome analyses from resistance-trained skeletal muscle. Moreover, normalization procedures need to account for global changes in rRNA and mRNA abundance.

Genome-wide transcriptome variation landscape in Ruta chalepensis organs revealed potential genes responsible for rutin biosynthesis

Ruta chalepensis L., most commonly known as 'fringed rue,' is an excellent and valuable bioactive plant that produces a range of complex flavonoids, of which rutin is the major compound present in this plant of great pharmaceutical and medicinal significance. The present study is a pioneering attempt to examine the changes in the transcriptomic landscape of leaf, stem, and root tissues and correlate this with rutin quantity in each tissue in order to identify the candidate genes responsible for rutin biosynthesis and to increase genomic resources in fringed rue. Comparative transcriptome sequencing of leaves, stems and roots were performed using the NovaSeq 6000 platform. The de novo transcriptome assembly generated 254,685 transcripts representing 154,018 genes with GC content of 42.60 % and N50 of 2280 bp. Searching assembled transcripts against UniRef90 and SwissProt databases annotated 79.7 % of them as protein coding. The leaf tissues had the highest rutin content followed by stems and roots. Several differentially expressed genes and transcripts relating to rutin biosynthesis were identified in leaves comparing with roots or stems comparing with roots. All the genes known to be involved in rutin biosynthesis showed up-regulation in leaves as compared with roots. These results were confirmed by gene ontology (GO) and pathway enrichment analyses. Up-regulated genes in leaves as compared with roots enriched GO terms with relation to rutin biosynthesis e.g. action of flavonol synthase, biosynthetic mechanism of malonyl-CoA, and action of monooxygenase. Phylogenetic analysis of the rhamnosyltransferase (RT) gene showed that it was highly homologues with RT sequence from Citrus species and all were located in the same clade. This transcriptomic dataset will serve as an important public resource for future genomics and transcriptomic studies in R. chalepensis and will act as a benchmark for the identification and genetic modification of genes involved in the biosynthesis of secondary metabolites.

Functional characterisation of the transcriptome from leaf tissue of the fluoroacetate-producing plant, Dichapetalum cymosum, in response to mechanical wounding

Dichapetalum cymosum produces the toxic fluorinated metabolite, fluoroacetate, presumably as a defence mechanism. Given the rarity of fluorinated metabolites in nature, the biosynthetic origin and function of fluoroacetate have been of particular interest. However, the mechanism for fluorination in D. cymosum was never elucidated. More importantly, there is a severe lack in knowledge on a genetic level for fluorometabolite-producing plants, impeding research on the subject. Here, we report on the first transcriptome for D. cymosum and investigate the wound response for insights into fluorometabolite production. Mechanical wounding studies were performed and libraries of the unwounded (control) and wounded (30 and 60 min post wounding) plant were sequenced using the Illumina HiSeq platform. A combined reference assembly generated 77,845 transcripts. Using the SwissProt, TrEMBL, GO, eggNOG, KEGG, Pfam, EC and PlantTFDB databases, a 69% annotation rate was achieved. Differential expression analysis revealed the regulation of 364 genes in response to wounding. The wound responses in D. cymosum included key mechanisms relating to signalling cascades, phytohormone regulation, transcription factors and defence-related secondary metabolites. However, the role of fluoroacetate in inducible wound responses remains unclear. Bacterial fluorinases were searched against the D. cymosum transcriptome but transcripts with homology were not detected suggesting the presence of a potentially different fluorinating enzyme in plants. Nevertheless, the transcriptome produced in this study significantly increases genetic resources available for D. cymosum and will assist with future research into fluorometabolite-producing plants.

Currently Applied Molecular Assays for Identifying ESR1 Mutations in Patients with Advanced Breast Cancer

Approximately 70% of breast cancers, the leading cause of cancer-related mortality worldwide, are positive for the estrogen receptor (ER). Treatment of patients with luminal subtypes is mainly based on endocrine therapy. However, ER positivity is reduced and ESR1 mutations play an important role in resistance to endocrine therapy, leading to advanced breast cancer. Various methodologies for the detection of ESR1 mutations have been developed, and the most commonly used method is next-generation sequencing (NGS)-based assays (50.0%) followed by droplet digital PCR (ddPCR) (45.5%). Regarding the sample type, tissue (50.0%) was more frequently used than plasma (27.3%). However, plasma (46.2%) became the most used method in 2016-2019, in contrast to 2012-2015 (22.2%). In 2016-2019, ddPCR (61.5%), rather than NGS (30.8%), became a more popular method than it was in 2012-2015. The easy accessibility, non-invasiveness, and demonstrated usefulness with high sensitivity of ddPCR using plasma have changed the trends. When using these assays, there should be a comprehensive understanding of the principles, advantages, vulnerability, and precautions for interpretation. In the future, advanced NGS platforms and modified ddPCR will benefit patients by facilitating treatment decisions efficiently based on information regarding ESR1 mutations.

Systematic comparison and assessment of RNA-seq procedures for gene expression quantitative analysis

RNA-seq is currently considered the most powerful, robust and adaptable technique for measuring gene expression and transcription activation at genome-wide level. As the analysis of RNA-seq data is complex, it has prompted a large amount of research on algorithms and methods. This has resulted in a substantial increase in the number of options available at each step of the analysis. Consequently, there is no clear consensus about the most appropriate algorithms and pipelines that should be used to analyse RNA-seq data. In the present study, 192 pipelines using alternative methods were applied to 18 samples from two human cell lines and the performance of the results was evaluated. Raw gene expression signal was quantified by non-parametric statistics to measure precision and accuracy. Differential gene expression performance was estimated by testing 17 differential expression methods. The procedures were validated by qRT-PCR in the same samples. This study weighs up the advantages and disadvantages of the tested algorithms and pipelines providing a comprehensive guide to the different methods and procedures applied to the analysis of RNA-seq data, both for the quantification of the raw expression signal and for the differential gene expression.