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

Unveiling the impact of hypodermis on gene expression for advancing bioprinted full-thickness 3D skin models

3D skin models have been explored as an alternative method to the use of animals in research and development. Usually, human skin equivalents comprise only epidermis or epidermis/dermis layers. Herein, we leverage 3D bioprinting technology to fabricate a full-thickness human skin equivalent with hypodermis (HSEH). The collagen hydrogel-based structure provides a mimetic environment for skin cells to adhere, proliferate and differentiate. The effective incorporation of the hypodermis layer is evidenced by scanning electron microscopy, immunofluorescence, and hematoxylin and eosin staining. The transcriptome results underscore the pivotal role of the hypodermis in orchestrating the genetic expression of a multitude of genes vital for skin functionality, including hydration, development and differentiation. Accordingly, we evidence the paramount significance of full-thickness human skin equivalents with hypodermis layer to provide an accurate in vitro platform for disease modeling and toxicology studies.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

The coevolution of fungus-ant agriculture

Fungus-farming ants cultivate multiple lineages of fungi for food, but, because fungal cultivar relationships are largely unresolved, the history of fungus-ant coevolution remains poorly known. We designed probes targeting >2000 gene regions to generate a dated evolutionary tree for 475 fungi and combined it with a similarly generated tree for 276 ants. We found that fungus-ant agriculture originated ~66 million years ago when the end-of-Cretaceous asteroid impact temporarily interrupted photosynthesis, causing global mass extinctions but favoring the proliferation of fungi. Subsequently, ~27 million years ago, one ancestral fungal cultivar population became domesticated, i.e., obligately mutualistic, when seasonally dry habitats expanded in South America, likely isolating the cultivar population from its free-living, wet forest-dwelling conspecifics. By revealing these and other major transitions in fungus-ant coevolution, our results clarify the historical processes that shaped a model system for nonhuman agriculture.

Anisakis extracellular vesicles elicit immunomodulatory and potentially tumorigenic outcomes on human intestinal organoids

BACKGROUND

Anisakis spp. are zoonotic nematodes causing mild to severe acute and chronic gastrointestinal infections. Chronic anisakiasis can lead to erosive mucosal ulcers, granulomas and inflammation, potential tumorigenic triggers. How Anisakis exerts its pathogenic potential through extracellular vesicles (EVs) and whether third-stage infective larvae may favor a tumorigenic microenvironment remain unclear.

METHODS

Here, we investigated the parasite's tumorigenic and immunomodulatory capabilities using comparative transcriptomics, qRT-PCR and protein analysis with multiplex ELISA on human intestinal organoids exposed to Anisakis EVs. Moreover, EVs were characterized in terms of shape, size and concentration using classic TEM, SEM and NTA analyses and advanced interferometric NTA.

RESULTS

Anisakis EVs showed classic shape features and a median average diameter of around 100 nm, according to NTA and iNTA. Moreover, a refractive index of 5-20% of non-water content suggested their effective biological cargo. After treatment of human intestinal organoids with Anisakis EVs, an overall parasitic strategy based on mitigation of the immune and inflammatory response was observed. Anisakis EVs impacted gene expression of main cytokines, cell cycle regulation and protein products. Seven key genes related to cell cycle regulation and apoptosis were differentially expressed in organoids exposed to EVs. In particular, the downregulation of EPHB2 and LEFTY1 and upregulation of NUPR1 genes known to be associated with colorectal cancer were observed, suggesting their involvement in tumorigenic microenvironment. A statistically significant reduction in specific mediators of inflammation and cell-cycle regulation from the polarized epithelium as IL-33R, CD40 and CEACAM1 from the apical chambers and IL-1B, GM-CSF, IL-15 and IL-23 from both chambers were observed.

CONCLUSIONS

The results here obtained unravel intestinal epithelium response to Anisakis EVs, impacting host's anthelminthic strategies and revealing for the first time to our knowledge the host-parasite interactions in the niche environment of an emerging accidental zoonosis. Use of an innovative EV characterization approach may also be useful for study of other helminth EVs, since the knowledge in this field is very limited.

Identification of the mechanistic basis of nitrogen responsiveness in two contrasting Setaria italica accessions

Nitrogen (N) is a macronutrient limiting crop productivity with varied requirements across species and genotypes. Understanding the mechanistic basis of N responsiveness by comparing contrasting genotypes could inform the development and selection of varieties with lower N demands, or inform agronomic practices to sustain yields with lower N inputs. Given the established role of millets in ensuring climate-resilient food and nutrition security, we investigated the physiological and genetic basis of nitrogen responsiveness in foxtail millet (Setaria italica L.). We had previously identified genotypic variants linked to N responsiveness, and here we dissect the mechanistic basis of the trait by examining the physiological and molecular behaviour of N responsive (NRp-SI58) and non-responsive (NNRp-SI114) accessions at high and low N. Under high N, NRp-SI58 allocates significantly more biomass to nodes, internodes and roots, more N to developing grains, and is more effective at remobilizing flag leaf N compared with NNRp-SI114. Post-anthesis flag leaf gene expression suggests that differences in N induce much higher transcript abundance in NNRp-SI114 than NRp-SI58, a large proportion of which is potentially regulated by APETALA2 (AP2) transcription factors. Overall, the study provides novel insights into the regulation and manipulation of N responsiveness in S. italica.

Assessing benthic macroinvertebrate communities' spatial heterogeneity in Mediterranean transitional waters through eDNA metabarcoding

Transitional waters are important habitats both for biodiversity and ecological functions, providing valuable natural resources and relevant ecosystem services. However, they are highly susceptible to climate changes and anthropogenic pressures responsible for biodiversity losses and require specific biomonitoring programs. Benthic macroinvertebrates are suitable as ecological indicators of transitional waters, being affected by biological, chemical, and physical conditions of the ecosystems about their life cycles and space-use behaviour. The advent of high-throughput sequencing technologies has allowed biodiversity investigations, at the molecular level, in multiple ecosystems and for different ecological guilds. Benthic macroinvertebrate communities' composition has been investigated, at the molecular level, mainly through DNA extracted from sediments in marine and riverine ecosystems. In this work, benthic macroinvertebrate communities are explored through eDNA metabarcoding from water samples in a Mediterranean transitional water ecosystem. This research highlighted the validity of eDNA metabarcoding as an efficient tool for the assessment of benthic macroinvertebrate community structure in transitional waters, unveiling the spatial heterogeneity of benthic macroinvertebrate communities correlated to the measured environmental gradients. The results suggest that peculiar features of transitional water ecosystems, such as shallow waters and limited currents, facilitate the assessment of benthic macroinvertebrate communities through environmental DNA analysis from surface water samples, opening for more rapid and accurate monitoring programs for these valuable ecosystems.

HuR controls glutaminase RNA metabolism

Glutaminase (GLS) is directly related to cell growth and tumor progression, making it a target for cancer treatment. The RNA-binding protein HuR (encoded by the ELAVL1 gene) influences mRNA stability and alternative splicing. Overexpression of ELAVL1 is common in several cancers, including breast cancer. Here we show that HuR regulates GLS mRNA alternative splicing and isoform translation/stability in breast cancer. Elevated ELAVL1 expression correlates with high levels of the glutaminase isoforms C (GAC) and kidney-type (KGA), which are associated with poor patient prognosis. Knocking down ELAVL1 reduces KGA and increases GAC levels, enhances glutamine anaplerosis into the TCA cycle, and drives cells towards glutamine dependence. Furthermore, we show that combining chemical inhibition of GLS with ELAVL1 silencing synergistically decreases breast cancer cell growth and invasion. These findings suggest that dual inhibition of GLS and HuR offers a therapeutic strategy for breast cancer treatment.

Chasing Sequencing Perfection: Marching Toward Higher Accuracy and Lower Costs

Next-generation sequencing (NGS), represented by Illumina platforms, has been an essential cornerstone of basic and applied research. However, the sequencing error rate of 1 per 1000 bp (10-3) represents a serious hurdle for research areas focusing on rare mutations, such as somatic mosaicism or microbe heterogeneity. By examining the high-fidelity sequencing methods developed in the past decade, we summarized three major factors underlying errors and the corresponding 12 strategies mitigating these errors. We then proposed a novel framework to classify 11 preexisting representative methods according to the corresponding combinatory strategies and identified three trends that emerged during methodological developments. We further extended this analysis to eight long-read sequencing methods, emphasizing error reduction strategies. Finally, we suggest two promising future directions that could achieve comparable or even higher accuracy with lower costs in both NGS and long-read sequencing.

miRNA-27a-3p is involved in the plasticity of differentiated hepatocytes

BACKGROUND

Epigenetic mechanisms, including DNA methylation, histone modifications, and chromatin remodeling, are highly involved in the regulation of hepatocyte viability, proliferation, and plasticity. We have previously demonstrated that repression of H3K27 methylation in differentiated hepatic HepaRG cells by treatment with GSK-J4, an inhibitor of JMJD3 and UTX H3K27 demethylase activity, changed their phenotype, inducing differentiated hepatocytes to proliferate. In addition to the epigenetic enzymatic role in the regulation of the retro-differentiation process, emerging evidence indicate that microRNAs (miRNAs) are involved in controlling hepatocyte proliferation during liver regeneration. Hence, the aim of this work is to investigate the impact of H3K27 methylation on miRNAs expression profile and its role in the regulation of the differentiation status of human hepatic progenitors HepaRG cells.

METHODS

A miRNA-sequencing was carried out in differentiated HepaRG cells treated or not with GSK-J4. Target searching and Gene Ontology analysis were performed to identify the molecular processes modulated by differentially expressed miRNAs. The biological functions of selected miRNAs was further investigated by transfection of miRNAs inhibitors or mimics in differentiated HepaRG cells followed by qPCR analysis, albumin ELISA assay, CD49a FACS analysis and EdU staining.

RESULTS

We identified 12 miRNAs modulated by GSK-J4; among these, miR-27a-3p and miR- 423-5p influenced the expression of several proliferation genes in differentiated HepaRG cells. MiR-27a-3p overexpression increased the number of hepatic cells reentering proliferation. Interestingly, both miR-27a-3p and miR-423-5p did not affect the expression levels of genes involved in the differentiation of progenitors HepaRG cells.

CONCLUSIONS

Modulation of H3K27me3 methylation in differentiated HepaRG cells, by GSK-J4 treatment, influenced miRNA' s expression profile pushing liver cells towards a proliferating phenotype. We demonstrated the involvement of miR-27a-3p in reinducing proliferation of differentiated hepatocytes suggesting a potential role in liver plasticity.

Clinical applications of metagenomics next-generation sequencing in infectious diseases

Infectious diseases are a great threat to human health. Rapid and accurate detection of pathogens is important in the diagnosis and treatment of infectious diseases. Metagenomics next-generation sequencing (mNGS) is an unbiased and comprehensive approach for detecting all RNA and DNA in a sample. With the development of sequencing and bioinformatics technologies, mNGS is moving from research to clinical application, which opens a new avenue for pathogen detection. Numerous studies have revealed good potential for the clinical application of mNGS in infectious diseases, especially in difficult-to-detect, rare, and novel pathogens. However, there are several hurdles in the clinical application of mNGS, such as: (1) lack of universal workflow validation and quality assurance; (2) insensitivity to high-host background and low-biomass samples; and (3) lack of standardized instructions for mass data analysis and report interpretation. Therefore, a complete understanding of this new technology will help promote the clinical application of mNGS to infectious diseases. This review briefly introduces the history of next-generation sequencing, mainstream sequencing platforms, and mNGS workflow, and discusses the clinical applications of mNGS to infectious diseases and its advantages and disadvantages.

Identification of a Unique Genomic Region in Sweet Chestnut (Castanea sativa Mill.) That Controls Resistance to Asian Chestnut Gall Wasp Dryocosmus kuriphilus Yasumatsu

The Asian chestnut gall wasp (ACGW) (Hymenoptera Dryocosmus kuriphilus Yasumatsu) is a severe pest of sweet chestnut (Castanea sativa Mill.) with a strong impact on growth and nut production. A comparative field trial in Central Italy, including provenances from Spain, Italy, and Greece, was screened for ACGW infestation over consecutive years. The Greek provenance Hortiatis expressed a high proportion of immune plants and was used to perform a genome-wide association study based on DNA pool sequencing (Pool-GWAS) by comparing two DNA pools from 25 susceptible and 25 resistant plants. DNA pools were sequenced with 50X coverage depth. Sequence reads were aligned to a C. mollissima reference genome and the pools were compared to identify SNPs associated with resistance. Twenty-one significant SNPs were identified and highlighted a small genomic region on pseudochromosome 3 (Chr 3), containing 12 candidate genes of three gene families: Cytochrome P450, UDP-glycosyltransferase, and Rac-like GTP-binding protein. Functional analyses revealed a putative metabolic gene cluster related to saccharide biosynthesis in the genomic regions associated with resistance that could be involved in the production of a toxic metabolite against parasites. The comparison with previous genetic studies confirmed the involvement of Chr 3 in the control of resistance to ACGW.

Genomic prediction of resistance to Hymenoscyphus fraxineus in common ash (Fraxinus excelsior L.) populations

The increase in introduced insect pests and pathogens due to anthropogenic environmental changes has become a major concern for tree species worldwide. Common ash (Fraxinus excelsior L.) is one of such species facing a significant threat from the invasive fungal pathogen Hymenoscyphus fraxineus. Some studies have indicated that the susceptibility of ash to the pathogen is genetically determined, providing some hope for accelerated breeding programs that are aimed at increasing the resistance of ash populations. To address this challenge, we used a genomic selection strategy to identify potential genetic markers that are associated with resistance to the pathogen causing ash dieback. Through genome-wide association studies (GWAS) of 300 common ash individuals from 30 populations across Poland (ddRAD, dataset A), we identified six significant SNP loci with a p-value ≤1 × 10-4 associated with health status. To further evaluate the effectiveness of GWAS markers in predicting health status, we considered two genomic prediction scenarios. Firstly, we conducted cross-validation on dataset A. Secondly, we trained markers on dataset A and tested them on dataset B, which involved whole-genome sequencing of 20 individuals from two populations. Genomic prediction analysis revealed that the top SNPs identified via GWAS exhibited notably higher prediction accuracies compared to randomly selected SNPs, particularly with a larger number of SNPs. Cross-validation analyses using dataset A showcased high genomic prediction accuracy, predicting tree health status with over 90% accuracy across the top SNP sets ranging from 500 to 10,000 SNPs from the GWAS datasets. However, no significant results emerged for health status when the model trained on dataset A was tested on dataset B. Our findings illuminate potential genetic markers associated with resistance to ash dieback, offering support for future breeding programs in Poland aimed at combating ash dieback and bolstering conservation efforts for this invaluable tree species.

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

Negative plant-soil feedback in Arabidopsis thaliana: Disentangling the effects of soil chemistry, microbiome, and extracellular self-DNA

Nutrient deficiency, natural enemies and litter autotoxicity have been proposed as possible mechanisms to explain species-specific negative plant-soil feedback (PSF). Another potential contributor to negative PSF is the plant released extracellular self-DNA during litter decay. In this study, we sought to comprehensively investigate these hypotheses by using Arabidopsis thaliana (L.) Heynh as a model plant in a feedback experiment. The experiment comprised a conditioning phase and a response phase in which the conditioned soils underwent four treatments: (i) addition of activated carbon, (ii) washing with tap water, (iii) sterilization by autoclaving, and (iv) control without any treatment. We evaluated soil chemical properties, microbiota by shotgun sequencing and the amount of A. thaliana extracellular DNA in the differently treated soils. Our results showed that washing and sterilization treatments mitigated the negative PSF effect. While shifts in soil chemical properties were not pronounced, significant changes in soil microbiota were observed, especially after sterilization. Notably, plant biomass was inversely associated with the content of plant self-DNA in the soil. Our results suggest that the negative PSF observed in the conditioned soil was associated to increased amounts of soilborne pathogens and plant self-DNA. However, fungal pathogens were not limited to negative conditions, butalso found in soils enhancing A.thaliana growth. In-depth multivariate analysis highlights that the hypothesis of negative PSF driven solely by pathogens lacks consistency. Instead, we propose a multifactorial explanation for the negative PSF buildup, in which the accumulation of self-DNA weakens the plant's root system, making it more susceptible to pathogens.

A comparative analysis of plastome evolution in autotrophic Piperales

PREMISE

Many plastomes of autotrophic Piperales have been reported to date, describing a variety of differences. Most studies focused only on a few species or a single genus, and extensive, comparative analyses have not been done. Here, we reviewed publicly available plastome reconstructions for autotrophic Piperales, reanalyzed publicly available raw data, and provided new sequence data for all previously missing genera. Comparative plastome genomics of >100 autotrophic Piperales were performed.

METHODS

We performed de novo assemblies to reconstruct the plastomes of newly generated sequence data. We used Sanger sequencing and read mapping to verify the assemblies and to bridge assembly gaps. Furthermore, we reconstructed the phylogenetic relationships as a foundation for comparative plastome genomics.

RESULTS

We identified a plethora of assembly and annotation issues in published plastome data, which, if unattended, will lead to an artificial increase of diversity. We were able to detect patterns of missing and incorrect feature annotation and determined that the inverted repeat (IR) boundaries were the major source for erroneous assembly. Accounting for the aforementioned issues, we discovered relatively stable junctions of the IRs and the small single-copy region (SSC), whereas the majority of plastome variations among Piperales stems from fluctuations of the boundaries of the IR and the large single-copy (LSC) region.

CONCLUSIONS

This study of all available plastomes of autotrophic Piperales, expanded by new data for previously missing genera, highlights the IR-LSC junctions as a potential marker for discrimination of various taxonomic levels. Our data indicates a pseudogene-like status for cemA and ycf15 in various Piperales. Based on a review of published data, we conclude that incorrect IR-SSC boundary identification is the major source for erroneous plastome assembly. We propose a gold standard for assembly and annotation of high-quality plastomes based on de novo assembly methods and appropriate references for gene annotation.

An Integrated Comprehensive Peptidomics and In Silico Analysis of Bioactive Peptide-Rich Milk Fermented by Three Autochthonous Cocci Strains

Bioactive peptides (BPs) are molecules of paramount importance with great potential for the development of functional foods, nutraceuticals or therapeutics for the prevention or treatment of various diseases. A functional BP-rich dairy product was produced by lyophilisation of bovine milk fermented by the autochthonous strains Lactococcus lactis subsp. lactis ZGBP5-51, Enterococcus faecium ZGBP5-52 and Enterococcus faecalis ZGBP5-53 isolated from the same artisanal fresh cheese. The efficiency of the proteolytic system of the implemented strains in the production of BPs was confirmed by a combined high-throughput mass spectrometry (MS)-based peptidome profiling and an in silico approach. First, peptides released by microbial fermentation were identified via a non-targeted peptide analysis (NTA) comprising reversed-phase nano-liquid chromatography (RP nano-LC) coupled with matrix-assisted laser desorption/ionisation-time-of-flight/time-of-flight (MALDI-TOF/TOF) MS, and then quantified by targeted peptide analysis (TA) involving RP ultrahigh-performance LC (RP-UHPLC) coupled with triple-quadrupole MS (QQQ-MS). A combined database and literature search revealed that 10 of the 25 peptides identified in this work have bioactive properties described in the literature. Finally, by combining the output of MS-based peptidome profiling with in silico bioactivity prediction tools, three peptides (75QFLPYPYYAKPA86, 40VAPFPEVFGK49, 117ARHPHPHLSF126), whose bioactive properties have not been previously reported in the literature, were identified as potential BP candidates.

Mozambican Coffea accessions from Ibo and Quirimba Islands: identification and geographical distribution

Mozambique does not have a tradition of farming Coffea arabica or Coffea canephora, the two species that dominate the worldwide coffee market. However, native coffee plants have been growing spontaneously and in some cases cultivated in the Ibo and Quirimba islands in the north of the country and Inhambane province in the south. Historically there has been confusion over the precise taxonomic classification of these indigenous coffee plants, with different botanists identifying the species as C. racemosa, C. zanguebariae or various synonyms of both. The present research aims to clarify the subject and provide new information on these little-described coffee species which may prove valuable as new breeding material for future cultivars, something that is sorely needed to face the present and future challenges of coffee production. Leaf samples were collected from 40 accessions from Ibo Island, Quirimba Island and Inhambane province. The samples were sequenced by whole-genome technology and WGS reads were filtered to identify relevant SNP variants. Diversity among the samples was assessed by PCA, and a phylogenetic tree including several Coffea species was built using additional data available in public databases. Experimental data confirm the presence of C. zanguebariae as the only coffee species present in both Ibo and Quirimba Islands, while it appears that C. racemosa is exclusive to the southern Inhambane province. The present research provides the most detailed analysis so far on the genetic identity of the traditional Mozambican coffee crops. This is the prerequisite for undertaking further scientific studies on these almost unknown coffee species and for starting agronomic development programs for the economic revival of Ibo and Quirimba islands based on coffee cultivation. Furthermore, these species could provide much-needed genetic material for the breeding of new hybrids with the two main commercial coffee species.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Seed bacterial microbiota in post-submergence tolerant and sensitive barley genotypes

Flooding is a predominant abiotic stress for cultivated plants, including barley. This cereal crop shows a large adaptability to different environmental conditions, suggesting the presence of key traits to tolerate adverse conditions. During germination, genetic variations account for dissimilarities in flooding tolerance. However, differences in the seed microbiota may also contribute to tolerance/sensitivity during seedling establishment. This work investigated differences in microbiome among the grains of barley accessions. Two barley phenotypes were compared, each either tolerant or sensitive to a short submergence period followed by a recovery. The study used a metataxonomic analysis based on 16S ribosomal RNA gene sequencing and subsequent functional prediction. Our results support the hypothesis that bacterial microbiota inhabiting the barley seeds are different between sensitive and tolerant barley accessions, which harbour specific bacterial phyla and families. Finally, bacteria detected in tolerant barley accessions show a peculiar functional enrichment that suggests a possible connection with successful germination and seedling establishment.

Analysis of >3400 worldwide eggplant accessions reveals two independent domestication events and multiple migration-diversification routes

Eggplant (Solanum melongena) is an important Solanaceous crop, widely cultivated and consumed in Asia, the Mediterranean basin, and Southeast Europe. Its domestication centers and migration and diversification routes are still a matter of debate. We report the largest georeferenced and genotyped collection to this date for eggplant and its wild relatives, consisting of 3499 accessions from seven worldwide genebanks, originating from 105 countries in five continents. The combination of genotypic and passport data points to the existence of at least two main centers of domestication, in Southeast Asia and the Indian subcontinent, with limited genetic exchange between them. The wild and weedy eggplant ancestor S. insanum shows admixture with domesticated S. melongena, similar to what was described for other fruit-bearing Solanaceous crops such as tomato and pepper and their wild ancestors. After domestication, migration and admixture of eggplant populations from different regions have been less conspicuous with respect to tomato and pepper, thus better preserving 'local' phenotypic characteristics. The data allowed the identification of misclassified and putatively duplicated accessions, facilitating genebank management. All the genetic, phenotypic, and passport data have been deposited in the Open Access G2P-SOL database, and constitute an invaluable resource for understanding the domestication, migration and diversification of this cosmopolitan vegetable.

Draft whole-genome sequences of Bacillus strains of potential probiotic for poultry

Bacillus species have been widely selected and used as probiotics for humans and animals. In this article, we reported draft whole-genome sequences of four Bacillus strains isolated from sourdough and chicken cecum and previously selected as potential probiotics for poultry. These genome sequences will provide a foundation for further characterization and understanding of their probiotic attributes.

Draft whole-genome sequences of competitive exclusion Ligilactobacillus salivarius strains for poultry

Competitive exclusion (CE) bacteria have been used to control the colonization of chickens by major foodborne pathogens. In this article, we report draft whole-genome sequences of three Ligilactobacillus salivarius strains isolated from chicken gastrointestinal tracts and previously selected as CE for poultry. These genome sequences will provide a foundation for further characterization and understanding of their CE attributes.

Microbiome depiction through user-adapted bioinformatic pipelines and parameters

Introduction. The role of the microbiome in health and disease continues to be increasingly recognized. However, there is significant variability in the bioinformatic protocols for analysing genomic data. This, in part, has impeded the potential incorporation of microbiomics into the clinical setting and has challenged interstudy reproducibility. In microbial compositional analysis, there is a growing recognition for the need to move away from a one-size-fits-all approach to data processing.Gap Statement. Few evidence-based recommendations exist for setting parameters of programs that infer microbiota community profiles despite these parameters significantly impacting the accuracy of taxonomic inference.Aim. To compare three commonly used programs (DADA2, QIIME2, and mothur) and optimize them into four user-adapted pipelines for processing paired-end amplicon reads. We aim to increase the accuracy of compositional inference and help standardize microbiomic protocol.Methods. Two key parameters were isolated across four pipelines: filtering sequence reads based on a whole-number error threshold (maxEE) and truncating read ends based on a quality score threshold (QTrim). Closeness of sample inference was then evaluated using a mock community of known composition.Results. We observed that raw genomic data lost were proportionate to how stringently parameters were set. Exactly how much data were lost varied by pipeline. Accuracy of sample inference correlated with increased sequence read retention. Falsely detected taxa and unaccounted for microbial constituents were unique to pipeline and parameter. Implementation of optimized parameter values led to better approximation of the known mock community.Conclusions. Microbial compositions generated based on the 16S rRNA marker gene should be interpreted with caution. To improve microbial community profiling, bioinformatic protocols must be user-adapted. Analysis should be performed with consideration for the select target amplicon, pipelines and parameters used, and taxa of interest.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.

Metabolic and Molecular Rearrangements of Sauvignon Blanc (Vitis vinifera L.) Berries in Response to Foliar Applications of Specific Dry Yeast

Dry yeast extracts (DYE) are applied to vineyards to improve aromatic and secondary metabolic compound content and wine quality; however, systematic information on the underpinning molecular mechanisms is lacking. This work aimed to unravel, through a systematic approach, the metabolic and molecular responses of Sauvignon Blanc berries to DYE treatments. To accomplish this, DYE spraying was performed in a commercial vineyard for two consecutive years. Berries were sampled at several time points after the treatment, and grapes were analyzed for sugars, acidity, free and bound aroma precursors, amino acids, and targeted and untargeted RNA-Seq transcriptional profiles. The results obtained indicated that the DYE treatment did not interfere with the technological ripening parameters of sugars and acidity. Some aroma precursors, including cys-3MH and GSH-3MH, responsible for the typical aromatic nuances of Sauvignon Blanc, were stimulated by the treatment during both vintages. The levels of amino acids and the global RNA-seq transcriptional profiles indicated that DYE spraying upregulated ROS homeostatic and thermotolerance genes, as well as ethylene and jasmonic acid biosynthetic genes, and activated abiotic and biotic stress responses. Overall, the data suggested that the DYE reduced berry oxidative stress through the regulation of specific subsets of metabolic and hormonal pathways.

Development and application of Single Primer Enrichment Technology (SPET) SNP assay for population genomics analysis and candidate gene discovery in lettuce

Single primer enrichment technology (SPET) is a novel high-throughput genotyping method based on short-read sequencing of specific genomic regions harboring polymorphisms. SPET provides an efficient and reproducible method for genotyping target loci, overcoming the limits associated with other reduced representation library sequencing methods that are based on a random sampling of genomic loci. The possibility to sequence regions surrounding a target SNP allows the discovery of thousands of closely linked, novel SNPs. In this work, we report the design and application of the first SPET panel in lettuce, consisting of 41,547 probes spanning the whole genome and designed to target both coding (~96%) and intergenic (~4%) regions. A total of 81,531 SNPs were surveyed in 160 lettuce accessions originating from a total of 10 countries in Europe, America, and Asia and representing 10 horticultural types. Model ancestry population structure clearly separated the cultivated accessions (Lactuca sativa) from accessions of its presumed wild progenitor (L. serriola), revealing a total of six genetic subgroups that reflected a differentiation based on cultivar typology. Phylogenetic relationships and principal component analysis revealed a clustering of butterhead types and a general differentiation between germplasm originating from Western and Eastern Europe. To determine the potentiality of SPET for gene discovery, we performed genome-wide association analysis for main agricultural traits in L. sativa using six models (GLM naive, MLM, MLMM, CMLM, FarmCPU, and BLINK) to compare their strength and power for association detection. Robust associations were detected for seed color on chromosome 7 at 50 Mbp. Colocalization of association signals was found for outer leaf color and leaf anthocyanin content on chromosome 9 at 152 Mbp and on chromosome 5 at 86 Mbp. The association for bolting time was detected with the GLM, BLINK, and FarmCPU models on chromosome 7 at 164 Mbp. Associations were detected in chromosomal regions previously reported to harbor candidate genes for these traits, thus confirming the effectiveness of SPET for GWAS. Our findings illustrated the strength of SPET for discovering thousands of variable sites toward the dissection of the genomic diversity of germplasm collections, thus allowing a better characterization of lettuce collections.

An empirical evaluation of the estimation of inbreeding depression from molecular markers under suboptimal conditions

Inbreeding depression (ID), the reduction in fitness due to inbreeding, is typically measured by the regression of the phenotypic values of individuals for a particular trait on their corresponding inbreeding coefficients (F). While genealogical records can provide these coefficients, they may be unavailable or incomplete, making molecular markers a useful alternative. The power to detect ID and its accuracy depend on the variation of F values of individuals, the sample sizes available, and the accuracy in the estimation of individual fitness traits and F values. In this study, we used Drosophila melanogaster to evaluate the effectiveness of molecular markers in estimating ID under suboptimal conditions. We generated two sets of 100 pairs of unrelated individuals from a large panmictic population and mated them for two generations to produce non-inbred and unrelated individuals (F = 0) and inbred individuals (full-sib progeny; F = 0.25). Using these expected genealogical F values, we calculated inbreeding depression for two fitness-related traits, pupae productivity and competitive fitness. We then sequenced the males from 17 non-inbred pairs and 17 inbred pairs to obtain their genomic inbreeding coefficients and estimate ID for the two traits. The scenario assumed was rather restrictive in terms of estimation of ID because: (1) the individuals belonged to the same generation of a large panmictic population, leading to low variation in individual F coefficients; (2) the sample sizes were small; and (3) the traits measured depended on both males and females while only males were sequenced. Despite the challenging conditions of our study, we found that molecular markers provided estimates of ID that were comparable to those obtained from simple pedigree estimations with larger sample sizes. The results therefore suggest that genomic measures of inbreeding are useful to provide estimates of inbreeding depression even under very challenging scenarios.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Tomato Plant Microbiota under Conventional and Organic Fertilization Regimes in a Soilless Culture System

Tomato is the main vegetable cultivated under soilless culture systems (SCSs); production of organic tomato under SCSs has increased due to consumer demands for healthier and environmentally friendly vegetables. However, organic tomato production under SCSs has been associated with low crop performance and fruit quality defects. These agricultural deficiencies could be linked to alterations in tomato plant microbiota; nonetheless, this issue has not been sufficiently addressed. Thus, the main goal of the present study was to characterize the rhizosphere and phyllosphere of tomato plants cultivated under conventional and organic SCSs. To accomplish this goal, tomato plants grown in commercial greenhouses under conventional or organic SCSs were tested at 8, 26, and 44 weeks after seedling transplantation. Substrate (n = 24), root (n = 24), and fruit (n = 24) composite samples were subjected to DNA extraction and high-throughput 16S rRNA gene sequencing. The present study revealed that the tomato core microbiota was predominantly constituted by Proteobacteria, Actinobacteria, and Firmicutes. Remarkably, six bacterial families, Bacillaceae, Microbacteriaceae, Nocardioidaceae, Pseudomonadaceae, Rhodobacteraceae, and Sphingomonadaceae, were shared among all substrate, rhizosphere, and fruit samples. Importantly, it was shown that plants under organic SCSs undergo a dysbiosis characterized by significant changes in the relative abundance of Bradyrhizobiaceae, Caulobacteraceae, Chitinophagaceae, Enterobacteriaceae, Erythrobacteraceae, Flavobacteriaceae, Nocardioidaceae, Rhodobacteraceae, and Streptomycetaceae. These results suggest that microbial alterations in substrates, roots, and fruits could be potential factors in contributing to the crop performance and fruit quality deficiencies observed in organic SCSs.

The genomic and bioclimatic characterization of Ethiopian barley (Hordeum vulgare L.) unveils challenges and opportunities to adapt to a changing climate

The climate crisis is impacting agroecosystems and threatening food security of millions of smallholder farmers. Understanding the potential for current and future climatic adaptation of local crop agrobiodiversity may guide breeding efforts and support resilience of agriculture. Here, we combine a genomic and climatic characterization of a large collection of traditional barley varieties from Ethiopia, a staple for local smallholder farmers cropping in challenging environments. We find that the genomic diversity of barley landraces can be partially traced back to geographic and environmental diversity of the landscape. We employ a machine learning approach to model Ethiopian barley adaptation to current climate and to identify areas where its existing diversity may not be well adapted in future climate scenarios. We use this information to identify optimal trajectories of assisted migration compensating to detrimental effects of climate change, finding that Ethiopian barley diversity bears opportunities for adaptation to the climate crisis. We then characterize phenology traits in the collection in two common garden experiments in Ethiopia, using genome-wide association approaches to identify genomic loci associated with timing of flowering and maturity of the spike. We combine this information with genotype-environment associations finding that loci involved in flowering time may also explain environmental adaptation. Our data show that integrated genomic, climatic, and phenotypic characterizations of agrobiodiversity may provide breeding with actionable information to improve local adaptation in smallholder farming systems.

Long-lasting impact of chitooligosaccharide application on strigolactone biosynthesis and fungal accommodation promotes arbuscular mycorrhiza in Medicago truncatula

The establishment of arbuscular mycorrhiza (AM) between plants and Glomeromycotina fungi is preceded by the exchange of chemical signals: fungal released Myc-factors, including chitooligosaccharides (CO) and lipo-chitooligosaccharides (LCO), activate plant symbiotic responses, while root-exuded strigolactones stimulate hyphal branching and boost CO release. Furthermore, fungal signaling reinforcement through CO application was shown to promote AM development in Medicago truncatula, but the cellular and molecular bases of this effect remained unclear. Here, we focused on long-term M. truncatula responses to CO treatment, demonstrating its impact on the transcriptome of both mycorrhizal and nonmycorrhizal roots over several weeks and providing an insight into the mechanistic bases of the CO-dependent promotion of AM colonization. CO treatment caused the long-lasting regulation of strigolactone biosynthesis and fungal accommodation-related genes. This was mirrored by an increase in root didehydro-orobanchol content, and the promotion of accommodation responses to AM fungi in root epidermal cells. Lastly, an advanced downregulation of AM symbiosis marker genes was observed at the latest time point in CO-treated plants, in line with an increased number of senescent arbuscules. Overall, CO treatment triggered molecular, metabolic, and cellular responses underpinning a protracted acceleration of AM development.

Recovery and Analysis of Long-Read Metagenome-Assembled Genomes

The development of long-read nucleic acid sequencing is beginning to make very substantive impact on the conduct of metagenome analysis, particularly in relation to the problem of recovering the genomes of member species of complex microbial communities. Here we outline bioinformatics workflows for the recovery and characterization of complete genomes from long-read metagenome data and some complementary procedures for comparison of cognate draft genomes and gene quality obtained from short-read sequencing and long-read sequencing.

The DNA/RNA helicase DHX9 contributes to the transcriptional program of the androgen receptor in prostate cancer

Background

Prostate cancer (PC) is the most commonly diagnosed male malignancy and an important cause of mortality. Androgen deprivation therapy is the first line treatment but, unfortunately, a large part of patients evolves to a castration-resistant stage, for which no effective cure is currently available. The DNA/RNA helicase DHX9 is emerging as an important regulator of cellular processes that are often deregulated in cancer.

Methods

To investigate whether DHX9 modulates PC cell transcriptome we performed RNA-sequencing analyses upon DHX9 silencing in the androgen-responsive cell line LNCaP. Bioinformatics and functional analyses were carried out to elucidate the mechanism of gene expression regulation by DHX9. Data from The Cancer Genome Atlas were mined to evaluate the potential role of DHX9 in PC.

Results

We found that up-regulation of DHX9 correlates with advanced stage and is associated with poor prognosis of PC patients. High-throughput RNA-sequencing analysis revealed that depletion of DHX9 in androgen-sensitive LNCaP cells affects expression of hundreds of genes, which significantly overlap with known targets of the Androgen Receptor (AR). Notably, AR binds to the DHX9 promoter and induces its expression, while Enzalutamide-mediated inhibition of AR activity represses DHX9 expression. Moreover, DHX9 interacts with AR in LNCaP cells and its depletion significantly reduced the recruitment of AR to the promoter region of target genes and the ability of AR to promote their expression in response to 5α-dihydrotestosterone. Consistently, silencing of DXH9 negatively affected androgen-induced PC cell proliferation and migration.

Conclusions

Collectively, our data uncover a new role of DHX9 in the control of the AR transcriptional program and establish the existence of an oncogenic DHX9/AR axis, which may represent a new druggable target to counteract PC progression.

Divergent gene expression responses in two Baltic Sea heterotrophic model bacteria to dinoflagellate dissolved organic matter

Phytoplankton release massive amounts of dissolved organic matter (DOM) into the water column during recurring blooms in coastal waters and inland seas. The released DOM encompasses a complex mixture of both known and unknown compounds, and is a rich nutrient source for heterotrophic bacteria. The metabolic activity of bacteria during and after phytoplankton blooms can hence be expected to reflect the characteristics of the released DOM. We therefore investigated if bacterioplankton could be used as “living sensors” of phytoplankton DOM quantity and/or quality, by applying gene expression analyses to identify bacterial metabolisms induced by DOM. We used transcriptional analysis of two Baltic Sea bacterial isolates (Polaribacter sp. BAL334 [Flavobacteriia] and Brevundimonas sp. BAL450 [Alphaproteobacteria]) growing with DOM from axenic cultures of the dinoflagellate Prorocentrum minimum. We observed pronounced differences between the two bacteria both in growth and the expressed metabolic pathways in cultures exposed to dinoflagellate DOM compared with controls. Differences in metabolic responses between the two isolates were caused both by differences in gene repertoire between them (e.g. in the SEED categories for membrane transport, motility and photoheterotrophy) and the regulation of expression (e.g. fatty acid metabolism), emphasizing the importance of separating the responses of different taxa in analyses of community sequence data. Similarities between the bacteria included substantially increased expression of genes for Ton and Tol transport systems in both isolates, which are commonly associated with uptake of complex organic molecules. Polaribacter sp. BAL334 showed stronger metabolic responses to DOM harvested from exponential than stationary phase dinoflagellates (128 compared to 26 differentially expressed genes), whereas Brevundimonas sp. BAL450 responded more to the DOM from stationary than exponential phase dinoflagellates (33 compared to 6 differentially expressed genes). These findings suggest that shifts in bacterial metabolisms during different phases of phytoplankton blooms can be detected in individual bacterial species and can provide insights into their involvement in DOM transformations.

Species and genetic diversity relationships in benthic macroinvertebrate communities along a salinity gradient

BACKGROUND

Species- and genetic diversity can change in parallel, resulting in a species-genetic diversity correlation (SGDC) and raising the question if the same drivers influence both biological levels of diversity. The SGDC can be either positive or negative, depending on whether the species diversity and the genetic diversity of the measured species respond in the same or opposite way to drivers. Using a traditional species diversity approach together with ultra-conserved elements and high throughput sequencing, we evaluated the SGDCs in benthic macrofauna communities in the Baltic Sea, a geologically young brackish water sea characterised by its steep salinity gradient and low species richness. Assessing SGDCs from six focal marine invertebrate species from different taxonomic groups and with differing life histories and ecological functions on both a spatial and temporal scale gives a more comprehensive insight into the community dynamics of this young ecosystem and the extrinsic factors that might drive the SGDCs.

RESULTS

No significant correlations between species diversity and genetic diversity were found for any of the focal species. However, both negative and positive trends of SGDCs for the individual focal species were observed. When examining the environmental drivers, no common trends between the species were found, even when restricting the analysis to specific taxonomic classes. Additionally, there were no common environmental factors driving the diversity relationships for species sharing the same SGDC trend (positive or negative). Local population dynamics, together with the invasion history of the individual species and their unique adaptation to the distinctive environment of the Baltic Sea, are expected to be of major influence on the outcome of the SGDCs.

CONCLUSIONS

The present results highlight the importance of assessing SGDCs using multiple species, not just a single indicator species. This emphasises a need to pay attention to the ecology and life history of the focal species. This study also provides insight into the large differences in both patterns and drivers of genetic diversity, which is important when including genetic biodiversity in conservation plans. We conclude that the effects of environmental and biological factors and processes that affects diversity patterns at both the community and genetic levels are likely species dependent, even in an environment such as the Baltic Sea with strong environmental gradients.

Role of Bifidobacterium pseudocatenulatum in Degradation and Consumption of Xylan-Derived Carbohydrates

Xylans, a family of xylose-based polysaccharides, are dietary fibers resistant to digestion. They therefore reach the large intestine intact; there, they are utilized by members of the gut microbiota. They are initially broken down by primary degraders that utilize extracellular xylanases to cleave xylan into smaller oligomers. The resulting xylooligosaccharides (XOS) can either be further metabolized directly by primary degraders or cross-feed secondary consumers, including Bifidobacterium. While several Bifidobacterium species have metabolic systems for XOS, most grow poorly on longer-chain XOS and xylan substrates. In this study, we isolated strains of Bifidobacterium pseudocatenulatum and observed that some, including B. pseudocatenulatum ED02, displayed growth on XOS with a high degree of polymerization (DP) and straight-chain xylan, suggesting a primary degrader phenotype that is rare in Bifidobacterium. In silico analyses revealed that only the genomes of these xylan-fermenting (xylan+) strains contained an extracellular GH10 endo-β-1.4 xylanase, a key enzyme for primary degradation of xylan. The presence of an extracellular xylanase was confirmed by the appearance of xylan hydrolysis products in cell-free supernatants. Extracellular xylanolytic activity was only detected in xylan+ strains, as indicated by the production of XOS fragments with a DP of 2 to 6, identified by thin-layer chromatography (TLC) and high-performance liquid chromatography (HPLC). Additionally, in vitro fecal fermentations revealed that strains with a xylan+ phenotype can persist with xylan supplementation. These results indicate that xylan+ B. pseudocatenulatum strains may have a competitive advantage in the complex environment of the gastrointestinal tract, due to their ability to act as primary degraders of xylan through extracellular enzymatic degradation. IMPORTANCE The beneficial health effects of dietary fiber are now well established. Moreover, low fiber consumption is associated with increased risks of metabolic and systemic diseases. This so-called "fiber gap" also has a profound impact on the composition of the gut microbiome, leading to a disrupted or dysbiotic microbiota. Therefore, understanding the mechanisms by which keystone bacterial species in the gut utilize xylans and other dietary fibers may provide a basis for developing strategies to restore gut microbiome function. The results described here provide biochemical and genetic evidence for primary xylan utilization by human-derived Bifidobacterium pseudocatenulatum and show also that cooperative utilization of xylans occurs among other members of this species.

Development of SNP Markers from GWAS for Selecting Seed Coat and Aleurone Layers in Brown Rice (Oryza sativa L.)

Ninety-five percent of the general nutrients in rice are concentrated in the rice bran and germ, and many nutrients such as vitamins, minerals, dietary fiber, and essential fatty acids, as well as antioxidants such as tocopherol, are lost during milling. In this study, we investigated the thickness of seed coat and aleurone layers using a 294 rice core collection, and found candidate genes related to thickness of seed coat and aleurone layers, by performing a genome wide association study (GWAS) analysis using whole genome resequencing data. Two primer pairs that can be used as high-resolution melting (HRM) markers were developed. As a result of genotyping BC₂F₂ individuals derived from a cross between “Samgwang” and “Seolgaeng”, and using corresponding HRM markers, it was possible to finally develop HRM markers for selecting seed coat and aleurone layer thickness. This is expected to be used as basic data for the application of gene editing using CRISPR/Cas9 technology and for establishing a breeding strategy for high eating quality rice using molecular genetic technology.

Review of the Current State of Freely Accessible Web Tools for the Analysis of 16S rRNA Sequencing of the Gut Microbiome

Owing to the emergence and improvement of high-throughput technology and the associated reduction in costs, next-generation sequencing (NGS) technology has made large-scale sampling and sequencing possible. With the large volume of data produced, the processing and downstream analysis of data are important for ensuring meaningful results and interpretation. Problems in data analysis may be encountered if researchers have little experience in using programming languages, especially if they are clinicians and beginners in the field. A strategy for solving this problem involves ensuring easy access to commercial software and tools. Here, we observed the current status of free web-based tools for microbiome analysis that can help users analyze and handle microbiome data effortlessly. We limited our search to freely available web-based tools and identified MicrobiomeAnalyst, Mian, gcMeta, VAMPS, and Microbiome Toolbox. We also highlighted the various analyses that each web tool offers, how users can analyze their data using each web tool, and noted some of their limitations. From the abovementioned list, gcMeta, VAMPS, and Microbiome Toolbox had several issues that made the analysis more difficult. Over time, as more data are generated and accessed, more users will analyze microbiome data. Thus, the availability of free and easily accessible web tools can enable the easy use and analysis of microbiome data, especially for those users with less experience in using command-line interfaces.

A comparative analysis of genomic and phenomic predictions of growth-related traits in 3-way coffee hybrids

Genomic prediction has revolutionized crop breeding despite remaining issues of transferability of models to unseen environmental conditions and environments. Usage of endophenotypes rather than genomic markers leads to the possibility of building phenomic prediction models that can account, in part, for this challenge. Here, we compare and contrast genomic prediction and phenomic prediction models for 3 growth-related traits, namely, leaf count, tree height, and trunk diameter, from 2 coffee 3-way hybrid populations exposed to a series of treatment-inducing environmental conditions. The models are based on 7 different statistical methods built with genomic markers and ChlF data used as predictors. This comparative analysis demonstrates that the best-performing phenomic prediction models show higher predictability than the best genomic prediction models for the considered traits and environments in the vast majority of comparisons within 3-way hybrid populations. In addition, we show that phenomic prediction models are transferrable between conditions but to a lower extent between populations and we conclude that chlorophyll a fluorescence data can serve as alternative predictors in statistical models of coffee hybrid performance. Future directions will explore their combination with other endophenotypes to further improve the prediction of growth-related traits for crops.

Differential Expression of Genes between a Tolerant and a Susceptible Maize Line in Response to a Sugarcane Mosaic Virus Infection

To uncover novel genes associated with the Sugarcane mosaic virus (SCMV) response, we used RNA-Seq data to analyze differentially expressed genes (DEGs) and transcript expression pattern clusters between a tolerant/resistant (CI-RL1) and a susceptible (B73) line, in addition to the F1 progeny (CI-RL1xB73). A Gene Ontology (GO) enrichment of DEGs led us to propose three genes possibly associated with the CI-RL1 response: a heat shock 90-2 protein and two ABC transporters. Through a clustering analysis of the transcript expression patterns (CTEPs), we identified two genes putatively involved in viral systemic spread: the maize homologs to the PIEZO channel (ZmPiezo) and to the Potyvirus VPg Interacting Protein 1 (ZmPVIP1). We also observed the complex behavior of the maize eukaryotic factors ZmeIF4E and Zm-elfa (involved in translation), homologs to eIF4E and eEF1α in A. thaliana. Together, the DEG and CTEPs results lead us to suggest that the tolerant/resistant CI-RL1 response to the SCMV encompasses the action of diverse genes and, for the first time, that maize translation factors are associated with viral interaction.

Maximizing Molecular Data From Low-Quality Fluid-Preserved Specimens in Natural History Collections

Over the past decade, museum genomics studies have focused on obtaining DNA of sufficient quality and quantity for sequencing from fluid-preserved natural history specimens, primarily to be used in systematic studies. While these studies have opened windows to evolutionary and biodiversity knowledge of many species worldwide, published works often focus on the success of these DNA sequencing efforts, which is undoubtedly less common than obtaining minimal or sometimes no DNA or unusable sequence data from specimens in natural history collections. Here, we attempt to obtain and sequence DNA extracts from 115 fresh and 41 degraded samples of homalopsid snakes, as well as from two degraded samples of a poorly known snake, Hydrablabes periops. Hydrablabes has been suggested to belong to at least two different families (Natricidae and Homalopsidae) and with no fresh tissues known to be available, intractable museum specimens currently provide the only opportunity to determine this snake’s taxonomic affinity. Although our aim was to generate a target-capture dataset for these samples, to be included in a broader phylogenetic study, results were less than ideal due to large amounts of missing data, especially using the same downstream methods as with standard, high-quality samples. However, rather than discount results entirely, we used mapping methods with references and pseudoreferences, along with phylogenetic analyses, to maximize any usable molecular data from our sequencing efforts, identify the taxonomic affinity of H. periops, and compare sequencing success between fresh and degraded tissue samples. This resulted in largely complete mitochondrial genomes for five specimens and hundreds to thousands of nuclear loci (ultra-conserved loci, anchored-hybrid enrichment loci, and a variety of loci frequently used in squamate phylogenetic studies) from fluid-preserved snakes, including a specimen of H. periops from the Field Museum of Natural History collection. We combined our H. periops data with previously published genomic and Sanger-sequenced datasets to confirm the familial designation of this taxon, reject previous taxonomic hypotheses, and make biogeographic inferences for Hydrablabes. A second H. periops specimen, despite being seemingly similar for initial raw sequencing results and after being put through the same protocols, resulted in little usable molecular data. We discuss the successes and failures of using different pipelines and methods to maximize the products from these data and provide expectations for others who are looking to use DNA sequencing efforts on specimens that likely have degraded DNA.

Life Science Identifier (Hydrablabes periops)

urn:lsid:zoobank.org:pub:F2AA44 E2-D2EF-4747-972A-652C34C2C09D.

Bacterial Endophytes Contribute to Rice Seedling Establishment Under Submergence

Flooding events caused by severe rains and poor soil drainage can interfere with plant germination and seedling establishment. Rice is one of the cereal crops that has unique germination strategies under flooding. One of these strategies is based on the fast coleoptile elongation in order to reach the water surface and re-establish the contact with the air. Microorganisms can contribute to plant health via plant growth promoters and provide protection from abiotic stresses. To characterise the community composition of the microbiome in rice germination under submergence, a 16S rRNA gene profiling metagenomic analysis was performed of temperate japonica rice varieties Arborio and Lamone seedlings, which showed contrasting responses in terms of coleoptile length when submerged. This analysis showed a distinct microbiota composition of Arborio seeds under submergence, which are characterised by the development of a long coleoptile. To examine the potential function of microbial communities under submergence, culturable bacteria were isolated, identified and tested for plant growth-promoting activities. A subgroup of isolated bacteria showed the capacity to hydrolyse starch and produce indole-related compounds under hypoxia. Selected bacteria were inoculated in seeds to evaluate their effect on rice under submergence, showing a response that is dependent on the rice genotype. Our findings suggest that endophytic bacteria possess plant growth-promoting activities that can substantially contribute to rice seedling establishment under submergence.

Analysis and pharmacological modulation of senescence in human epithelial stem cells

Human epithelial stem cells (ESCs) are characterized by long‐term regenerative properties, much dependent on the tissue of origin and varying during their lifespan. We analysed such variables in cultures of ESCs isolated from the skin, conjunctiva, limbus and oral mucosa of healthy donors and patients affected by ectrodactyly‐ectodermal dysplasia‐clefting syndrome, a rare genetic disorder caused by mutations in the p63 gene. We cultured cells until exhaustion in the presence or in the absence of DAPT (γ‐secretase inhibitor; N‐[N‐(3, 5‐difluorophenacetyl)‐L‐alanyl]‐S‐phenylglycine T‐butyl ester). All cells were able to differentiate in vitro but exhibited variable self‐renewal potential. In particular, cells carrying p63 mutations stopped prematurely, compared with controls. Importantly, administration of DAPT significantly extended the replicative properties of all stem cells under examination. RNA sequencing analysis revealed that distinct sets of genes were up‐ or down‐regulated during their lifetime, thus allowing to identify druggable gene networks and off‐the‐shelf compounds potentially dealing with epithelial stem cell senescence. These data will expand our knowledge on the genetic bases of senescence and potentially pave the way to the pharmacological modulation of ageing in epithelial stem cells.

MycoVarP: Mycobacterium Variant and Drug Resistance Prediction Pipeline for Whole-Genome Sequence Data Analysis

Whole-genome sequencing (WGS) provides a comprehensive tool to analyze the bacterial genomes for genotype–phenotype correlations, diversity of single-nucleotide variant (SNV), and their evolution and transmission. Several online pipelines and standalone tools are available for WGS analysis of Mycobacterium tuberculosis (Mtb) complex (MTBC). While they facilitate the processing of WGS data with minimal user expertise, they are either too general, providing little insights into bacterium-specific issues such as gene variations, INDEL/synonymous/PE-PPE (IDP family), and drug resistance from sample data, or are limited to specific objectives, such as drug resistance. It is understood that drug resistance and lineage-specific issues require an elaborate prioritization of identified variants to choose the best target for subsequent therapeutic intervention. Mycobacterium variant pipeline (MycoVarP) addresses these specific issues with a flexible battery of user-defined and default filters. It provides an end-to-end solution for WGS analysis of Mtb variants from the raw reads and performs two quality checks, viz, before trimming and after alignments of reads to the reference genome. MycoVarP maps the annotated variants to the drug-susceptible (DS) database and removes the false-positive variants, provides lineage identification, and predicts potential drug resistance. We have re-analyzed the WGS data reported by Advani et al. (2019) using MycoVarP and identified some additional variants not reported so far. We conclude that MycoVarP will help in identifying nonsynonymous, true-positive, drug resistance–associated variants more effectively and comprehensively, including those within the IDP of the PE-PPE/PGRS family, than possible from the currently available pipelines.

Differential gene expression and chemical patterns of an intertidal crab inhabiting a polluted port and an adjacent marine protected area

The acquisition of data to safeguard marine protected areas located close to ports is important in order to develop plans that allow effective protection from pollution as well as sustainable development of the port. The area Secche della Meloria is a Marine Protected Area (MPA-MEL) three miles from Livorno Harbour (LH), which is characterized by a long history of pollution. Here we studied the bioaccumulation and transcriptomic patterns of the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) (Crustacea; Brachyura, Grapsidae), inhabiting the two selected sites. Results showed that the two crab populations are significantly different in their chemical composition of trace elements and Polyciclic Aromatic Hydrocarbons (PAHs), and gene expression patterns (1280 DEGs). Enrichment analysis indicated that crabs at LH had the highest stress response genes, and they were associated with higher levels of bioaccumulation detected in body tissues. We are confident that the significant differential gene expression profiles observed between crabs, characterized by significant chemical differences, is associated with responses to contaminant exposure.

Toward a Molecular Approach to Chronotype Assessment

The aim of the present study was to develop a Polygenic Score-based model for molecular chronotype assessment. Questionnaire-based phenotypical chronotype assessment was used as a reference. In total, 54 extremely morning/morning (MM/M; 35 females, 39.7 ± 3.8 years) and 44 extremely evening/evening (EE/E; 20 females, 27.3 ± 7.7 years) individuals donated a buccal DNA sample for genotyping by sequencing of the entire genetic variability of 19 target genes known to be involved in circadian rhythmicity and/or sleep duration. Targeted genotyping was performed using the single primer enrichment technology and a specifically designed panel of 5526 primers. Among 2868 high-quality polymorphisms, a cross-validation approach lead to the identification of 83 chronotype predictive variants, including previously known and also novel chronotype-associated polymorphisms. A large (35 single-nucleotide polymorphisms [SNPs]) and also a small (13 SNPs) panel were obtained, both with an estimated predictive validity of approximately 80%. Potential mechanistic hypotheses for the role of some of the newly identified variants in modulating chronotype are formulated. Once validated in independent populations encompassing the whole range of chronotypes, the identified panels might become useful within the setting of both circadian public health initiatives and precision medicine.

Isoprene enhances leaf cytokinin metabolism and induces early senescence

Isoprene, a major biogenic volatile hydrocarbon of climate-relevance, indisputably mitigates abiotic stresses in emitting plants. However functional relevance of constitutive isoprene emission in unstressed plants remains contested. Isoprene and cytokinins (CKs) are synthesized from a common substrate and pathway in chloroplasts. It was postulated that isoprene emission may affect CK-metabolism. Using transgenic isoprene-emitting (IE) Arabidopsis and isoprene nonemitting (NE) RNA-interference grey poplars (paired with respective NE and IE genotypes), the life of individual IE and NE leaves from emergence to abscission was followed under stress-free conditions. We monitored plant growth rate, aboveground developmental phenotype, modelled leaf photosynthetic energy status, quantified the abundance of leaf CKs, analysed Arabidopsis and poplar leaf transcriptomes by RNA-sequencing in presence and absence of isoprene during leaf senescence. Isoprene emission by unstressed leaves enhanced the abundance of CKs (isopentenyl adenine and its precursor) by > 200%, significantly upregulated genes coding for CK-synthesis, CK-signalling and CK-degradation, hastened plant development, increased chloroplast metabolic rate, altered photosynthetic energy status, induced early leaf senescence in both Arabidopsis and poplar. IE leaves senesced sooner even in decapitated poplars where source-sink relationships and hormone homeostasis were perturbed. Constitutive isoprene emission significantly accelerates CK-led leaf and organismal development and induces early senescence independent of growth constraints. Isoprene emission provides an early-riser evolutionary advantage and shortens lifecycle duration to assist rapid diversification in unstressed emitters.

Importance of clitellar tissue in the regeneration ability of earthworm Eudrilus eugeniae

Among the annelids, earthworms are renowned for their phenomenal ability to regenerate the lost segments. The adult earthworm Eudrilus eugeniae contains 120 segments and the body segments of the earthworm are divided into pre-clitellar, clitellar and post-clitellar segments. The present study denoted that clitellum plays vital role in the successful regeneration of the species. We have performed histological studies to identify among the three skin layers of the earthworm, which cellular layer supports the blastema formation and regeneration of the species. The histological evidences denoted that the proliferation of the longitudinal cell layer at the amputation site is crucial for the successful regeneration of the earthworm and it takes place only in the presence of an intact clitellum. Besides we have performed clitellar transcriptome analysis of the earthworm Eudrilus eugeniae to monitor the key differentially expressed genes and their associated functions and pathways controlling the clitellar tissue changes during both anterior and posterior regeneration of the earthworm. A total of 4707 differentially expressed genes (DEGs) were identified between the control clitellum and clitellum of anterior regenerated earthworms and 4343 DEGs were detected between the control clitellum and clitellum of posterior regenerated earthworms. The functional enrichment analysis confirmed the genes regulating the muscle mass shape and structure were significantly downregulated and the genes associated with response to starvation and anterior-posterior axis specification were significantly upregulated in the clitellar tissue during both anterior and posterior regeneration of the earthworm. The RNA sequencing data of clitellum and the comparative transcriptomic analysis were helpful to understand the complex regeneration process of the earthworm.

Mitochondrial phylogenomics and mitogenome organization in the parasitoid wasp family Braconidae (Hymenoptera: Ichneumonoidea)

BACKGROUND

Mitochondrial (mt) nucleotide sequence data has been by far the most common tool employed to investigate evolutionary relationships. While often considered to be more useful for shallow evolutionary scales, mt genomes have been increasingly shown also to contain valuable phylogenetic information about deep relationships. Further, mt genome organization provides another important source of phylogenetic information and gene reorganizations which are known to be relatively frequent within the insect order Hymenoptera. Here we used a dense taxon sampling comprising 148 mt genomes (132 newly generated) collectively representing members of most of the currently recognised subfamilies of the parasitoid wasp family Braconidae, which is one of the largest radiations of hymenopterans. We employed this data to investigate the evolutionary relationships within the family and to assess the phylogenetic informativeness of previously known and newly discovered mt gene rearrangements.

RESULTS

Most subfamilial relationships and their composition obtained were similar to those recovered in a previous phylogenomic study, such as the restoration of Trachypetinae and the recognition of Apozyginae and Proteropinae as valid braconid subfamilies. We confirmed and detected phylogenetic signal in previously known as well as novel mt gene rearrangements, including mt rearrangements within the cyclostome subfamilies Doryctinae and Rogadinae.

CONCLUSIONS

Our results showed that both the mt genome DNA sequence data and gene organization contain valuable phylogenetic signal to elucidate the evolution within Braconidae at different taxonomic levels. This study serves as a basis for further investigation of mt gene rearrangements at different taxonomic scales within the family.

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.