NCBI BLAST+ integrated into Galaxy

Metagenomic analysis of the gut microbiota in piglets either challenged or not with enterotoxigenic Escherichia coli reveals beneficial effects of probiotics on microbiome composition, resistome, digestive function and oxidative stress responses

This study used metagenomic analysis to investigate the gut microbiota and resistome in piglets that were or were not challenged with enterotoxigenic Escherichia coli (ETEC) and had or had not received dietary supplementation with microencapsulated probiotics. The 72 piglets belonged to six groups that were either non-ETEC challenged (groups 1–3) or ETEC challenged (receiving 5ml of 10⁹ CFU/ml pathogenic ETEC strain L3.2 one week following weaning at three weeks of age: groups 4–6). On five occasions at 2, 5, 8, 11, and 14 days of piglet age, groups 2 and 5 were supplemented with 10⁹ CFU/ml of multi-strain probiotics (Lactiplantibacillus plantarum strains 22F and 25F, and Pediococcus acidilactici 72N) while group 4 received 10⁹ CFU/ml of P. acidilactici 72N. Group 3 received 300mg/kg chlortetracycline in the weaner diet to mimic commercial conditions. Rectal faecal samples were obtained for metagenomic and resistome analysis at 2 days of age, and at 12 hours and 14 days after the timing of post-weaning challenge with ETEC. The piglets were all euthanized at 42 days of age. The piglets in groups 2 and 5 were enriched with several desirable microbial families, including Lactobacillaceae, Lachnospiraceae and Ruminococcaceae, while piglets in group 3 had increases in members of the Bacteroidaceae family and exhibited an increase in tetW and tetQ genes. Group 5 had less copper and multi-biocide resistance. Mobile genetic elements IncQ1 and IncX4 were the most prevalent replicons in antibiotic-fed piglets. Only groups 6 and 3 had the integrase gene (intl) class 2 and 3 detected, respectively. The insertion sequence (IS) 1380 was prevalent in group 3. IS3 and IS30, which are connected to dietary intake, were overrepresented in group 5. Furthermore, only group 5 showed genes associated with detoxification, with enrichment of genes associated with oxidative stress, glucose metabolism, and amino acid metabolism compared to the other groups. Overall, metagenomic analysis showed that employing a multi-strain probiotic could transform the gut microbiota, reduce the resistome, and boost genes associated with food metabolism.

Enabling genomic island prediction and comparison in multiple genomes to investigate bacterial evolution and outbreaks

Outbreaks of virulent and/or drug-resistant bacteria have a significant impact on human health and major economic consequences. Genomic islands (GIs; defined as clusters of genes of probable horizontal origin) are of high interest because they disproportionately encode virulence factors, some antimicrobial-resistance (AMR) genes, and other adaptations of medical or environmental interest. While microbial genome sequencing has become rapid and inexpensive, current computational methods for GI analysis are not amenable for rapid, accurate, user-friendly and scalable comparative analysis of sets of related genomes. To help fill this gap, we have developed IslandCompare, an open-source computational pipeline for GI prediction and comparison across several to hundreds of bacterial genomes. A dynamic and interactive visualization strategy displays a bacterial core-genome phylogeny, with bacterial genomes linearly displayed at the phylogenetic tree leaves. Genomes are overlaid with GI predictions and AMR determinants from the Comprehensive Antibiotic Resistance Database (CARD), and regions of similarity between the genomes are also displayed. GI predictions are performed using Sigi-HMM and IslandPath-DIMOB, the two most precise GI prediction tools based on nucleotide composition biases, as well as a novel blast-based consistency step to improve cross-genome prediction consistency. GIs across genomes sharing sequence similarity are grouped into clusters, further aiding comparative analysis and visualization of acquisition and loss of mobile GIs in specific sub-clades. IslandCompare is an open-source software that is containerized for local use, plus available via a user-friendly, web-based interface to allow direct use by bioinformaticians, biologists and clinicians (at https://islandcompare.ca).

The predictive potential of different molecular markers linked to amikacin susceptibility phenotypes in Pseudomonas aeruginosa

Informed antibiotic prescription offers a practical solution to antibiotic resistance problem. With the increasing affordability of different sequencing technologies, molecular-based resistance prediction would direct proper antibiotic selection and preserve available agents. Amikacin is a broad-spectrum aminoglycoside exhibiting higher clinical efficacy and less resistance rates in Ps. aeruginosa due to its structural nature and its ability to achieve higher serum concentrations at lower therapeutic doses. This study examines the predictive potential of molecular markers underlying amikacin susceptibility phenotypes in order to provide improved diagnostic panels. Using a predictive model, genes and variants underlying amikacin resistance have been statistically and functionally explored in a large comprehensive and diverse set of Ps. aeruginosa completely sequenced genomes. Different genes and variants have been examined for their predictive potential and functional correlation to amikacin susceptibility phenotypes. Three predictive sets of molecular markers have been identified and can be used in a complementary manner, offering promising molecular diagnostics. armR, nalC, nalD, mexR, mexZ, ampR, rmtD, nalDSer32Asn, fusA1Y552C, fusA1D588G, arnAA170T, and arnDG206C have been identified as the best amikacin resistance predictors in Ps. aeruginosa while faoAT385A, nuoGA890T, nuoGA574T, lptAT55A, lptAR62S, pstBR87C, gidBE126G, gidBQ28K, amgSE108Q, and rplYQ41L have been identified as the best amikacin susceptibility predictors. Combining different measures of predictive performance together with further functional analysis can help design new and more informative molecular diagnostic panels. This would greatly inform and direct point of care diagnosis and prescription, which would consequently preserve amikacin functionality and usefulness.

Colonization of the live biotherapeutic product VE303 and modulation of the microbiota and metabolites in healthy volunteers

Manipulation of the gut microbiota via fecal microbiota transplantation (FMT) has shown clinical promise in diseases such as recurrent Clostridioides difficile infection (rCDI). However, the variable nature of this approach makes it challenging to describe the relationship between fecal strain colonization, corresponding microbiota changes, and clinical efficacy. Live biotherapeutic products (LBPs) consisting of defined consortia of clonal bacterial isolates have been proposed as an alternative therapeutic class because of their promising preclinical results and safety profile. We describe VE303, an LBP comprising 8 commensal Clostridia strains under development for rCDI, and its early clinical development in healthy volunteers (HVs). In a phase 1a/b study in HVs, VE303 is determined to be safe and well-tolerated at all doses tested. VE303 strains optimally colonize HVs if dosed over multiple days after vancomycin pretreatment. VE303 promotes the establishment of a microbiota community known to provide colonization resistance.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

The core and accessory Hfq interactomes across Pseudomonas aeruginosa lineages

The major RNA-binding protein Hfq interacts with mRNAs, either alone or together with regulatory small noncoding RNAs (sRNAs), affecting mRNA translation and degradation in bacteria. However, studies tend to focus on single reference strains and assume that the findings may apply to the entire species, despite the important intra-species genetic diversity known to exist. Here, we use RIP-seq to identify Hfq-interacting RNAs in three strains representing the major phylogenetic lineages of Pseudomonas aeruginosa. We find that most interactions are in fact not conserved among the different strains. We identify growth phase-specific and strain-specific Hfq targets, including previously undescribed sRNAs. Strain-specific interactions are due to different accessory gene sets, RNA abundances, or potential context- or sequence- dependent regulatory mechanisms. The accessory Hfq interactome includes most mRNAs encoding Type III Secretion System (T3SS) components and secreted toxins in two strains, as well as a cluster of CRISPR guide RNAs in one strain. Conserved Hfq targets include the global virulence regulator Vfr and metabolic pathways involved in the transition from fast to slow growth. Furthermore, we use rGRIL-seq to show that RhlS, a quorum sensing sRNA, activates Vfr translation, thus revealing a link between quorum sensing and virulence regulation. Overall, our work highlights the important intra-species diversity in post-transcriptional regulatory networks in Pseudomonas aeruginosa.

The APAF1_C/WD40 repeat domain-encoding gene from the sea lettuce Ulva mutabilis sheds light on the evolution of NB-ARC domain-containing proteins in green plants

We advance Ulva's genetic tractability and highlight its value as a model organism by characterizing its APAF1_C/WD40 domain-encoding gene, which belongs to a family that bears homology to R genes. The multicellular chlorophyte alga Ulva mutabilis (Ulvophyceae, Ulvales) is native to coastal ecosystems worldwide and attracts both high socio-economic and scientific interest. To further understand the genetic mechanisms that guide its biology, we present a protocol, based on adapter ligation-mediated PCR, for retrieving flanking sequences in U. mutabilis vector-insertion mutants. In the created insertional library, we identified a null mutant with an insertion in an apoptotic protease activating factor 1 helical domain (APAF1_C)/WD40 repeat domain-encoding gene. Protein domain architecture analysis combined with phylogenetic analysis revealed that this gene is a member of a subfamily that arose early in the evolution of green plants (Viridiplantae) through the acquisition of a gene that also encoded N-terminal nucleotide-binding adaptor shared by APAF-1, certain R-gene products and CED-4 (NB-ARC) and winged helix-like (WH-like) DNA-binding domains. Although phenotypic analysis revealed no mutant phenotype, gene expression levels in control plants correlated to the presence of bacterial symbionts, which U. mutabilis requires for proper morphogenesis. In addition, our analysis led to the discovery of a putative Ulva nucleotide-binding site and leucine-rich repeat (NBS-LRR) Resistance protein (R-protein), and we discuss how the emergence of these R proteins in green plants may be linked to the evolution of the APAF1_C/WD40 protein subfamily.

Assembly and Annotation of Escherichia coli Bacteriophage U115

Here, we present the annotated genome sequence of Escherichia coli bacteriophage U115, a T4-like bacteriophage. Phage U115 has a genome length of 166,986 bp and has 286 predicted genes.

Delineation of ISEcp1 and IS26-Mediated Plasmid Fusion Processes by MinION Single-Molecule Long-Read Sequencing

We recently reported the recovery of a novel IncI1 type conjugative helper plasmid which could target mobile genetic elements (MGE) located in non-conjugative plasmid and form a fusion conjugative plasmid to mediate the horizontal transfer of the non-conjugative plasmid. In this study, interactions between the helper plasmid pSa42-91k and two common MGEs, ISEcp1 and IS15DI, which were cloned into a pBackZero-T vector, were monitored during the conjugation process to depict the molecular mechanisms underlying the plasmid fusion process mediated by insertion sequence (IS) elements. The MinION single-molecule long-read sequencing technology can dynamically reveal the plasmid recombination events and produce valuable information on genetic polymorphism and plasmid heterogeneity in different multidrug resistance (MDR) encoding bacteria. Such data would facilitate the development of new strategies to control evolution and dissemination of MDR plasmids.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.)

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Protocol for using NoBadWordsCombiner to merge and minimize "bad words" from BLAST hits against multiple eukaryotic gene annotation databases

Annotating protein-coding genes can be challenging, especially when searching for the best hits against multiple functional databases. This is partly because of "bad words" appearing as top hits, such as hypothetical or uncharacterized proteins. To help alleviate some of these issues, we designed a bioinformatics tool called NoBadWordsCombiner, which efficiently merges the hits from various databases, strengthening gene definitions by minimizing functional descriptions containing "bad words." Unlike other available tools, NoBadWordsCombiner is user friendly, but it does require users to have some general bioinformatics skills, including a basic understanding of the BLAST package and dash shell in Linux/Unix environments. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2021a).

Iocasia fonsfrigidae NS-1 gen. nov., sp. nov., a Novel Deep-Sea Bacterium Possessing Diverse Carbohydrate Metabolic Pathways

Resolving metabolisms of deep-sea microorganisms is crucial for understanding ocean energy cycling. Here, a strictly anaerobic, Gram-negative strain NS-1 was isolated from the deep-sea cold seep in the South China Sea. Phylogenetic analysis based on 16S rRNA gene sequence indicated that strain NS-1 was most closely related to the type strain Halocella cellulosilytica DSM 7362^T (with 92.52% similarity). A combination of phylogenetic, genomic, and physiological traits with strain NS-1, was proposed to be representative of a novel genus in the family Halanaerobiaceae, for which Iocasia fonsfrigidae NS-1 was named. It is noteworthy that I. fonsfrigidae NS-1 could metabolize multiple carbohydrates including xylan, alginate, starch, and lignin, and thereby produce diverse fermentation products such as hydrogen, lactate, butyrate, and ethanol. The expressions of the key genes responsible for carbohydrate degradation as well as the production of the above small molecular substrates when strain NS-1 cultured under different conditions, were further analyzed by transcriptomic methods. We thus predicted that part of the ecological role of Iocasia sp. is likely in the fermentation of products from the degradation of diverse carbohydrates to produce hydrogen as well as other small molecules, which are in turn utilized by other members of cold seep microbes.

The Chemosensory Transcriptome of a Diving Beetle

Insects astoundingly dominate Earth’s land ecosystems and have a huge impact on human life. Almost every aspect of their life relies upon their highly efficient and adaptable chemosensory system. In the air, most chemical signals that are detected at long range are hydrophobic molecules, which insects detect using proteins encoded by multigenic families that emerged following land colonization by insect ancestors, namely the odorant-binding proteins (OBPs) and the odorant receptors (ORs). However, land-to-freshwater transitions occurred in many lineages within the insect tree of life. Whether chemosensory gene repertoires of aquatic insects remained essentially unchanged or underwent more or less drastic modifications to cope with physico-chemical constraints associated with life underwater remains virtually unknown. To address this issue, we sequenced and analyzed the transcriptome of chemosensory organs of the diving beetle Rhantus suturalis (Coleoptera, Dytiscidae). A reference transcriptome was assembled de novo using reads from five RNA-seq libraries (male and female antennae, male and female palps, and wing muscle). It contained 47,570 non-redundant unigenes encoding proteins of more than 50 amino acids. Within this reference transcriptome, we annotated sequences coding 53 OBPs, 48 ORs, 73 gustatory receptors (GRs), and 53 ionotropic receptors (IRs). Phylogenetic analyses notably revealed a large OBP gene expansion (35 paralogs in R. suturalis) as well as a more modest OR gene expansion (9 paralogs in R. suturalis) that may be specific to diving beetles. Interestingly, these duplicated genes tend to be expressed in palps rather than in antennae, suggesting a possible adaptation with respect to the land-to-water transition. This work provides a strong basis for further evolutionary and functional studies that will elucidate how insect chemosensory systems adapted to life underwater.

Determination of the two-component systems regulatory network reveals core and accessory regulations across Pseudomonas aeruginosa lineages

Pseudomonas aeruginosa possesses one of the most complex bacterial regulatory networks, which largely contributes to its success as a pathogen. However, most of its transcription factors (TFs) are still uncharacterized and the potential intra-species variability in regulatory networks has been mostly ignored so far. Here, we used DAP-seq to map the genome-wide binding sites of all 55 DNA-binding two-component systems (TCSs) response regulators (RRs) across the three major P. aeruginosa lineages. The resulting networks encompass about 40% of all genes in each strain and contain numerous new regulatory interactions across most major physiological processes. Strikingly, about half of the detected targets are specific to only one or two strains, revealing a previously unknown large functional diversity of TFs within a single species. Three main mechanisms were found to drive this diversity, including differences in accessory genome content, as exemplified by the strain-specific plasmid in IHMA87 outlier strain which harbors numerous binding sites of conserved chromosomally-encoded RRs. Additionally, most RRs display potential auto-regulation or RR-RR cross-regulation, bringing to light the vast complexity of this network. Overall, we provide the first complete delineation of the TCSs regulatory network in P. aeruginosa that will represent an important resource for future studies on this pathogen.

Laniakea@ReCaS: exploring the potential of customisable Galaxy on-demand instances as a cloud-based service

BACKGROUND

Improving the availability and usability of data and analytical tools is a critical precondition for further advancing modern biological and biomedical research. For instance, one of the many ramifications of the COVID-19 global pandemic has been to make even more evident the importance of having bioinformatics tools and data readily actionable by researchers through convenient access points and supported by adequate IT infrastructures. One of the most successful efforts in improving the availability and usability of bioinformatics tools and data is represented by the Galaxy workflow manager and its thriving community. In 2020 we introduced Laniakea, a software platform conceived to streamline the configuration and deployment of "on-demand" Galaxy instances over the cloud. By facilitating the set-up and configuration of Galaxy web servers, Laniakea provides researchers with a powerful and highly customisable platform for executing complex bioinformatics analyses. The system can be accessed through a dedicated and user-friendly web interface that allows the Galaxy web server's initial configuration and deployment.

RESULTS

"Laniakea@ReCaS", the first instance of a Laniakea-based service, is managed by ELIXIR-IT and was officially launched in February 2020, after about one year of development and testing that involved several users. Researchers can request access to Laniakea@ReCaS through an open-ended call for use-cases. Ten project proposals have been accepted since then, totalling 18 Galaxy on-demand virtual servers that employ ~ 100 CPUs, ~ 250 GB of RAM and ~ 5 TB of storage and serve several different communities and purposes. Herein, we present eight use cases demonstrating the versatility of the platform.

CONCLUSIONS

During this first year of activity, the Laniakea-based service emerged as a flexible platform that facilitated the rapid development of bioinformatics tools, the efficient delivery of training activities, and the provision of public bioinformatics services in different settings, including food safety and clinical research. Laniakea@ReCaS provides a proof of concept of how enabling access to appropriate, reliable IT resources and ready-to-use bioinformatics tools can considerably streamline researchers' work.

The treasure inside barley seeds: microbial diversity and plant beneficial bacteria

BACKGROUND

Bacteria associated with plants can enhance the plants' growth and resistance against phytopathogens. Today, growers aim to reduce the use of mineral fertilizers and pesticides. Since phytopathogens cause severe yield losses in crop production systems, biological alternatives gain more attention. Plant and also seed endophytes have the potential to influence the plant, especially seed-borne bacteria may express their beneficiary impact at initial plant developmental stages. In the current study, we assessed the endophytic seed microbiome of seven genetically diverse barley accessions by 16S rRNA gene amplicon sequencing and verified the in vitro plant beneficial potential of isolated seed endophytes. Furthermore, we investigated the impact of the barley genotype and its seed microbiome on the rhizosphere microbiome at an early growth stage by 16S rRNA gene amplicon sequencing.

RESULTS

The plant genotype displayed a significant impact on the microbiota in both barley seed and rhizosphere. Consequently, the microbial alpha- and beta-diversity of the endophytic seed microbiome was highly influenced by the genotype. Interestingly, no correlation was observed between the endophytic seed microbiome and the single nucleotide polymorphisms of the seven genotypes. Unclassified members of Enterobacteriaceae were by far most dominant. Other abundant genera in the seed microbiome belonged to Curtobacterium, Paenibacillus, Pantoea, Sanguibacter and Saccharibacillus. Endophytes isolated from barley seeds were affiliated to dominant genera of the core seed microbiome, based on their 16S rRNA gene sequence. Most of these endophytic isolates produced in vitro plant beneficial secondary metabolites known to induce plant resistance.

CONCLUSION

Although barley accessions representing high genetic diversity displayed a genotype-dependent endophytic seed microbiome, a core seed microbiome with high relative abundances was identified. Endophytic isolates were affiliated to members of the core seed microbiome and many of them showed plant beneficial properties. We propose therefore that new breeding strategies should consider genotypes with high abundance of beneficial microbes.

A Novel In Silico Benchmarked Pipeline Capable of Complete Protein Analysis: A Possible Tool for Potential Drug Discovery

Simple Summary

Protein interactions govern the majority of an organism’s biological processes. Therefore, to fully understand the functionality of an organism, we must know how proteins work at a molecular level. This study assembled a protocol that enables scientists to construct a protein’s tertiary structure easily and subsequently to investigate its mechanism and function. Each step involved in prediction, validation, and functional analysis of a protein is crucial to obtain an accurate result. We have dubbed this the trifecta analysis. It was clear early in our research that no single study in the literature had previously encompassed the complete trifecta analysis. In particular, studies that recommend free, open-source tools that have been benchmarked for each step are lacking. The present study ensures that predictions are accurate and validated and will greatly benefit new and experienced scientists alike in obtaining a strong understanding of the trifecta analysis, resulting in a domino effect that could lead to drug development.

Abstract

Current in silico proteomics require the trifecta analysis, namely, prediction, validation, and functional assessment of a modeled protein. The main drawback of this endeavor is the lack of a single protocol that utilizes a proper set of benchmarked open-source tools to predict a protein’s structure and function accurately. The present study rectifies this drawback through the design and development of such a protocol. The protocol begins with the characterization of a novel coding sequence to identify the expressed protein. It then recognizes and isolates evolutionarily conserved sequence motifs through phylogenetics. The next step is to predict the protein’s secondary structure, followed by the prediction, refinement, and validation of its three-dimensional tertiary structure. These steps enable the functional analysis of the macromolecule through protein docking, which facilitates the identification of the protein’s active site. Each of these steps is crucial for the complete characterization of the protein under study. We have dubbed this process the trifecta analysis. In this study, we have proven the effectiveness of our protocol using the cystatin C and AChE proteins. Beginning with just their sequences, we have characterized both proteins’ structures and functions, including identifying the cystatin C protein’s seven-residue active site and the AChE protein’s active-site gorge via protein–protein and protein–ligand docking, respectively. This process will greatly benefit new and experienced scientists alike in obtaining a strong understanding of the trifecta analysis, resulting in a domino effect that could expand drug development.

Immune Response in Mice Immunized with Chimeric H1 Antigens

Identification of a universal influenza vaccine candidate has remained a global challenge for both humans and animals. This study describes an approach that uses consensus sequence building to generate chimeric HAs (cHAs): two resultant H1 HA-based chimeras comprising of conserved sequences (within several areas spanning the head and stalk regions) of H1 and H5 or H9 HAs. These cHAs expressed in Drosophila cells (S2) were used to immunize mice. All immunized mice were protected from an infectious H1 virus challenge. Seroconverted mice sera to the H1 cHAs inhibited both the challenge virus and an H5 virus isolate by haemagglutination inhibition (HI) assay. These findings further emphasize that cHAs induce cross-reactive antibodies against conserved areas of both head and stalk regions of the seasonal influenza A (H1N1) pdm09 virus' HA and holds potential for further development of a universal influenza vaccine.

Assembly and Annotation of the Complete Genome Sequence of T4-Like Bacteriophage 132

Here, we describe the complete sequence of bacteriophage 132, a T4-like Escherichia coli phage. Phage 132 has a genome of 166,922-bp length, with 286 predicted genes.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Global population structure of Haemophilus influenzae serotype a (Hia) and emergence of invasive Hia disease: capsule switching or capsule replacement?

The population structure of Hia was examined by interrogation of the H. influenzae MLST website. There were 196 entries of Hia with 55 sequence types (STs) identified (as of September 3, 2020). BURST analysis clustered related STs into four complexes with ST-23, ST-4, ST-21 and ST-62 identified as their ancestral STs. The majority of Hia entries (73.4%) and STs (65.5%) were identified as clonal division I (ST-23 and the ST-4 complexes). Only 43 (21.9%) entries and 14 STs (25.5%) were identified as clonal division II (ST-62 and ST-21 complexes). Current data suggested most invasive Hia belonged to clonal division I and the ST-23 complex while most clonal division II Hia were respiratory isolates with the exception of ST-62 which was common among invasive Hia in the U.S. southwest. Comparison of the capsule bexABCD genes from clonal divisions I and II strains showed sequence diversity with variations following the pattern of clonal divisions. Evidence from the literature and the current study suggests recent emergence of invasive Hia might be related to capsule replacement subsequent to the implementation of the Hib conjugate vaccine and possibly exacerbated by other conjugate vaccines that may have altered the microbial flora of the human respiratory tract.

The Histidine Biosynthetic Genes in the Superphylum Bacteroidota-Rhodothermota-Balneolota-Chlorobiota: Insights into the Evolution of Gene Structure and Organization

One of the most studied metabolic routes is the biosynthesis of histidine, especially in enterobacteria where a single compact operon composed of eight adjacent genes encodes the complete set of biosynthetic enzymes. It is still not clear how his genes were organized in the genome of the last universal common ancestor community. The aim of this work was to analyze the structure, organization, phylogenetic distribution, and degree of horizontal gene transfer (HGT) of his genes in the Bacteroidota-Rhodothermota-Balneolota-Chlorobiota superphylum, a group of phylogenetically close bacteria with different surviving strategies. The analysis of the large variety of his gene structures and organizations revealed different scenarios with genes organized in more or less compact—heterogeneous or homogeneous—operons, in suboperons, or in regulons. The organization of his genes in the extant members of the superphylum suggests that in the common ancestor of this group, genes were scattered throughout the chromosome and that different forces have driven the assembly of his genes in compact operons. Gene fusion events and/or paralog formation, HGT of single genes or entire operons between strains of the same or different taxonomic groups, and other molecular rearrangements shaped the his gene structure in this superphylum.

In Silico Identification and Functional Characterization of Conserved miRNAs in the Genome of Cryptosporidium parvum

Cryptosporidium parvum, a predominant causal agent of a fatal zoonotic protozoan diarrhoeal disease called cryptosporidiosis, bears a worldwide public health concern for childhood mortality and poses a key threat to the dairy and water industries. MicroRNAs (miRNAs), small but powerful posttranscriptional gene silencing RNA molecules, regulate a variety of molecular, biological, and cellular processes in animals and plants. As to the present date, there is a paucity of information regarding miRNAs of C. parvum; hence, this study was used to identify miRNAs in the organism using a comprehensible expressed sequence tag-based homology search approach consisting of a series of computational screening process from the identification of putative miRNA candidates to the functional annotation of the important gene targets in C. parvum. The results revealed a conserved miRNA that targeted 487 genes in the model organism (Drosophila melanogaster) and 85 genes in C. parvum, of which 11 genes had direct involvements in several crucial virulence factors such as environmental oocyst protection, excystation, locomotion, adhesion, invasion, stress protection, intracellular growth, and survival. Besides, 20 genes showed their association with various major pathways dedicated for the ribosomal biosynthesis, DNA repair, transportation, protein production, gene expression, cell cycle, cell proliferation, development, immune response, differentiation, and nutrient metabolism of the organism in the host. Thus, this study provides a strong evidence of great impact of identified miRNA on the biology, virulence, and pathogenesis of C. parvum. Furthermore, the study suggests that the detected miRNA could be a potential epigenomic tool for controlling the protozoon through silencing those virulent and pathway-related target genes.

Genomic characterization of 16S rRNA methyltransferase-producing Escherichia coli isolates from the Parisian area, France

BACKGROUND

The resistance to all aminoglycosides (AGs) conferred by 16S rRNA methyltransferase enzymes (16S-RMTases) is a major public health concern.

OBJECTIVES

To characterize the resistance genotype, its genetic environment and plasmid support, and the phylogenetic relatedness of 16S-RMTase-producing Escherichia coli from France.

METHODS

We screened 137 E. coli isolates resistant to all clinically relevant AGs from nine Parisian hospitals for 16S-RMTases. WGS was performed on clinical isolates with high-level AG resistance (MIC ≥256 mg/L) and their transformants.

RESULTS

Thirty of the 137 AG-resistant E. coli produced 16S-RMTases: 11 ArmA, 18 RmtB and 1 RmtC. The 16S-RMTase producers were also resistant to third-generation cephalosporins (90% due to a blaCTX-M gene), co-trimoxazole, fluoroquinolones and carbapenems (blaNDM and blaVIM genes) in 97%, 83%, 70% and 10% of cases, respectively. Phylogenomic diversity was high in ArmA producers, with 10 different STs, but a similar genetic environment, with the Tn1548 transposon carried by a plasmid closely related to pCTX-M-3 in 6/11 isolates. Conversely, RmtB producers belonged to 12 STs, the most frequent being ST405 and ST complex (STc) 10 (four and four isolates, respectively). The rmtB gene was carried by IncF plasmids in 10 isolates and was found in different genetic environments. The rmtC gene was carried by the pNDM-US plasmid.

CONCLUSIONS

ArmA and RmtB are the predominant 16S-RMTases in France, but their spread follows two different patterns: (i) dissemination of a conserved genetic support carrying armA in E. coli with high levels of genomic diversity; and (ii) various genetic environments surrounding rmtB in clonally related E. coli.

Assessing nucleic acid binding activity of four dinoflagellate cold shock domain proteins from Symbiodinium kawagutii and Lingulodinium polyedra

Background

Dinoflagellates have a generally large number of genes but only a small percentage of these are annotated as transcription factors. Cold shock domain (CSD) containing proteins (CSPs) account for roughly 60% of these. CSDs are not prevalent in other eukaryotic lineages, perhaps suggesting a lineage-specific expansion of this type of transcription factors in dinoflagellates, but there is little experimental data to support a role for dinoflagellate CSPs as transcription factors. Here we evaluate the hypothesis that dinoflagellate CSPs can act as transcription factors by binding double-stranded DNA in a sequence dependent manner.

Results

We find that both electrophoretic mobility shift assay (EMSA) competition experiments and selection and amplification binding (SAAB) assays indicate binding is not sequence specific for four different CSPs from two dinoflagellate species. Competition experiments indicate all four CSPs bind to RNA better than double-stranded DNA.

Conclusion

Dinoflagellate CSPs do not share the nucleic acid binding properties expected for them to function as bone fide transcription factors. We conclude the transcription factor complement of dinoflagellates is even smaller than previously thought suggesting that dinoflagellates have a reduced dependance on transcriptional control compared to other eukaryotes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00368-4.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2^T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2^T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2^T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and ¹³C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2^T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2^T using ¹³C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific ¹³C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

Annotation and Molecular Characterisation of the TaIRO3 and TaHRZ Iron Homeostasis Genes in Bread Wheat (Triticum aestivum L.)

Effective maintenance of plant iron (Fe) homoeostasis relies on a network of transcription factors (TFs) that respond to environmental conditions and regulate Fe uptake, translocation, and storage. The iron-related transcription factor 3 (IRO3), as well as haemerythrin motif-containing really interesting new gene (RING) protein and zinc finger protein (HRZ), are major regulators of Fe homeostasis in diploid species like Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa L.), but remain uncharacterised in hexaploid bread wheat (Triticum aestivum L.). In this study, we have identified, annotated, and characterised three TaIRO3 homoeologs and six TaHRZ1 and TaHRZ2 homoeologs in the bread wheat genome. Protein analysis revealed that TaIRO3 and TaHRZ proteins contain functionally conserved domains for DNA-binding, dimerisation, Fe binding, or polyubiquitination, and phylogenetic analysis revealed clustering of TaIRO3 and TaHRZ proteins with other monocot IRO3 and HRZ proteins, respectively. Quantitative reverse-transcription PCR analysis revealed that all TaIRO3 and TaHRZ homoeologs have unique tissue expression profiles and are upregulated in shoot tissues in response to Fe deficiency. After 24 h of Fe deficiency, the expression of TaHRZ homoeologs was upregulated, while the expression of TaIRO3 homoeologs was unchanged, suggesting that TaHRZ functions upstream of TaIRO3 in the wheat Fe homeostasis TF network.

A capture-based assay for detection and characterization of transposon polymorphisms in maize

Transposons can create allelic diversity that affects gene expression and phenotypic diversity. The detection of transposon polymorphisms at a genome-wide scale across a large population is difficult. Here we developed a targeted sequencing approach to monitor transposon polymorphisms of interest. This approach can interrogate the presence or absence of transposons reliably across various genotypes. Using this approach, we genotyped a set of 965 transposon-related presence/absence polymorphisms in a diverse panel of 16 maize (Zea mays L.) inbred lines that are representative of the major maize breeding groups. About 70% of the selected regions can be effectively assayed in each genotype. The consistency between the capture-based assay and PCR-based assay are 98.6% based on analysis of 24 randomly selected transposon polymorphisms. By integrating the transposon polymorphisms data with gene expression data, ∼18% of the assayed transposon polymorphisms were found to be associated with variable gene expression levels. A detailed analysis of 18 polymorphisms in a larger association panel confirmed the effects of 10 polymorphisms, with one of them having stronger association with expression than nearby SNP markers. The effects of seven polymorphisms were tested using a luciferase-based expression assay, and one was confirmed. Together, this study demonstrates that the targeted sequencing assay is an effective way to explore transposon function in a high-throughput manner.

Complete Genome Sequence of Escherichia coli Bacteriophage U136B

We report the genome sequence of bacteriophage U136B, which is reliant on the lipopolysaccharide and the antibiotic efflux protein TolC for infection of Escherichia coli and is a useful model for studying trade-offs and trade-ups that shape evolution. Phage U136B has a 49,233-bp genome with 87 predicted genes.

Molecular Physiological Characterization of a High Heat Resistant Spore Forming Bacillus subtilis Food Isolate

Bacterial endospores (spores) are among the most resistant living forms on earth. Spores of Bacillus subtilis A163 show extremely high resistance to wet heat compared to spores of laboratory strains. In this study, we found that spores of B. subtilis A163 were indeed very wet heat resistant and released dipicolinic acid (DPA) very slowly during heat treatment. We also determined the proteome of vegetative cells and spores of B. subtilis A163 and the differences in these proteomes from those of the laboratory strain PY79, spores of which are much less heat resistant. This proteomic characterization identified 2011 proteins in spores and 1901 proteins in vegetative cells of B. subtilis A163. Surprisingly, spore morphogenic protein SpoVM had no homologs in B. subtilis A163. Comparing protein expression between these two strains uncovered 108 proteins that were differentially present in spores and 93 proteins differentially present in cells. In addition, five of the seven proteins on an operon in strain A163, which is thought to be primarily responsible for this strain's spores high heat resistance, were also identified. These findings reveal proteomic differences of the two strains exhibiting different resistance to heat and form a basis for further mechanistic analysis of the high heat resistance of B. subtilis A163 spores.

Chromosome-scale genome assembly of Cucumis hystrix-a wild species interspecifically cross-compatible with cultivated cucumber

Cucumis hystrix Chakr. (2n = 2x = 24) is a wild species that can hybridize with cultivated cucumber (C. sativus L., 2n = 2x = 14), a globally important vegetable crop. However, cucumber breeding is hindered by its narrow genetic base. Therefore, introgression from C. hystrix has been anticipated to bring a breakthrough in cucumber improvement. Here, we report the chromosome-scale assembly of C. hystrix genome (289 Mb). Scaffold N50 reached 14.1 Mb. Over 90% of the sequences were anchored onto 12 chromosomes. A total of 23,864 genes were annotated using a hybrid method. Further, we conducted a comprehensive comparative genomic analysis of cucumber, C. hystrix, and melon (C. melo L., 2n = 2x = 24). Whole-genome comparisons revealed that C. hystrix is phylogenetically closer to cucumber than to melon, providing a molecular basis for the success of its hybridization with cucumber. Moreover, expanded gene families of C. hystrix were significantly enriched in "defense response," and C. hystrix harbored 104 nucleotide-binding site-encoding disease resistance gene analogs. Furthermore, 121 genes were positively selected, and 12 (9.9%) of these were involved in responses to biotic stimuli, which might explain the high disease resistance of C. hystrix. The alignment of whole C. hystrix genome with cucumber genome and self-alignment revealed 45,417 chromosome-specific sequences evenly distributed on C. hystrix chromosomes. Finally, we developed four cucumber-C. hystrix alien addition lines and identified the exact introgressed chromosome using molecular and cytological methods. The assembled C. hystrix genome can serve as a valuable resource for studies on Cucumis evolution and interspecific introgression breeding of cucumber.

Sulfur Metabolites Play Key System-Level Roles in Modulating Denitrification

Nitrate-reducing bacteria (NRB) and sulfate-reducing bacteria (SRB) colonize diverse anoxic environments, including soil subsurface, groundwater, and wastewater. NRB and SRB compete for resources, and their interplay has major implications on the global cycling of nitrogen and sulfur species, with undesirable outcomes in some contexts.

Competition between nitrate-reducing bacteria (NRB) and sulfate-reducing bacteria (SRB) for resources in anoxic environments is generally thought to be governed largely by thermodynamics. It is now recognized that intermediates of nitrogen and sulfur cycling (e.g., hydrogen sulfide, nitrite, etc.) can also directly impact NRB and SRB activities in freshwater, wastewater, and sediment and therefore may play important roles in competitive interactions. Here, through comparative transcriptomic and metabolomic analyses, we have uncovered mechanisms of hydrogen sulfide- and cysteine-mediated inhibition of nitrate respiratory growth for the NRB Intrasporangium calvum C5. Specifically, the systems analysis predicted that cysteine and hydrogen sulfide inhibit growth of I. calvum C5 by disrupting distinct steps across multiple pathways, including branched-chain amino acid (BCAA) biosynthesis, utilization of specific carbon sources, and cofactor metabolism. We have validated these predictions by demonstrating that complementation with BCAAs and specific carbon sources relieves the growth inhibitory effects of cysteine and hydrogen sulfide. We discuss how these mechanistic insights give new context to the interplay and stratification of NRB and SRB in diverse environments.

IMPORTANCE Nitrate-reducing bacteria (NRB) and sulfate-reducing bacteria (SRB) colonize diverse anoxic environments, including soil subsurface, groundwater, and wastewater. NRB and SRB compete for resources, and their interplay has major implications on the global cycling of nitrogen and sulfur species, with undesirable outcomes in some contexts. For instance, the removal of reactive nitrogen species by NRB is desirable for wastewater treatment, but in agricultural soils, NRB can drive the conversion of nitrates from fertilizers into nitrous oxide, a potent greenhouse gas. Similarly, the hydrogen sulfide produced by SRB can help sequester and immobilize toxic heavy metals but is undesirable in oil wells where competition between SRB and NRB has been exploited to suppress hydrogen sulfide production. By characterizing how reduced sulfur compounds inhibit growth and activity of NRB, we have gained systems-level and mechanistic insight into the interplay of these two important groups of organisms and drivers of their stratification in diverse environments.

Transcription Inhibitors with XRE DNA-Binding and Cupin Signal-Sensing Domains Drive Metabolic Diversification in Pseudomonas

Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized.

Transcription factors (TFs) are instrumental in the bacterial response to new environmental conditions. They can act as direct signal sensors and subsequently induce changes in gene expression leading to physiological adaptation. Here, by combining transcriptome sequencing (RNA-seq) and cistrome determination (DAP-seq), we studied a family of eight TFs in Pseudomonas aeruginosa. This family, encompassing TFs with XRE-like DNA-binding and cupin signal-sensing domains, includes the metabolic regulators ErfA, PsdR, and PauR and five so-far-unstudied TFs. The genome-wide delineation of their regulons identified 39 regulatory interactions with genes mostly involved in metabolism. We found that the XRE-cupin TFs are inhibitors of their neighboring genes, forming local, functional units encoding proteins with functions in condition-specific metabolic pathways. Growth phenotypes of isogenic mutants highlighted new roles for PauR and PA0535 in polyamines and arginine metabolism. The phylogenetic analysis of this family of regulators across the bacterial kingdom revealed a wide diversity of such metabolic regulatory modules and identified species with potentially higher metabolic versatility. Numerous genes encoding uncharacterized XRE-cupin TFs were found near metabolism-related genes, illustrating the need of further systematic characterization of transcriptional regulatory networks in order to better understand the mechanisms of bacterial adaptation to new environments.

IMPORTANCE Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized. With the recent advances in high-throughput sequencing methods, it is now possible to fill this knowledge gap and help the understanding of how bacteria adapt and thrive in new environments. By leveraging these methods, we provide an example of a comprehensive analysis of an entire family of transcription factors and bring new insights into metabolic and regulatory adaptation in the Pseudomonas genus.

Cryptosporidium species and subtypes in river water and riverbed sediment using next-generation sequencing

This study uncovered the prevalence, harboured species, and subtype diversity of Cryptosporidium species in river water and its sediment from the Apies River in South Africa. Cryptosporidium spp. concentrations in freshwater and its sediment were determined using Ziehl-Neelsen staining and quantitative Polymerase Chain Reaction (qPCR) techniques. Next-generation sequencing (NGS) targeting the 60 kDa glycoprotein (gp60) gene of Cryptosporidium spp. was performed to reveal the species, subtype families and subtypes harboured in freshwater and its sediment. Although the results revealed that water samples had a higher prevalence (30%) compared with sediment (28%), the number of observable Cryptosporidium spp. oocysts in sediment samples (ranging from 4.90 to 5.81 log₁₀ oocysts per 1 Liter) was higher than that of river water samples (ranging from 4.60 to 5.58 log₁₀ oocysts per 1 L) using Ziehl-Neelsen staining. The 18S ribosomal ribonucleic acid (rRNA) gene copy of Cryptosporidium in riverbed sediments ranged from 6.03 to 7.65 log₁₀, whereas in river water, it was found to be between 4.20 and 6.79 log₁₀. Subtyping results showed that in riverbed sediments, Cryptosporidium parvum accounted for 40.72% of sequences, followed by Cryptosporidium hominis with 23.64%, Cryptosporidium cuniculus with 7.10%, Cryptosporidium meleagridis with 4.44% and the least was Cryptosporidium wrairi with 2.59%. A considerable percentage of reads in riverbed sediment (21.25%) was not assigned to any subtype. River water samples had 45.63% of sequences assigned to C. parvum, followed by 30.32% to C. hominis, 17.99% to C. meleagridis and 5.88% to C. cuniculus. The data obtained are concerning, as Cryptosporidium spp. have intrinsic resistance to water treatment processes and low infectious doses, which can pose a risk to human health due to the various uses of water (for human consumption, leisure, and reuse).

Pseudomonas aeruginosa antitoxin HigA functions as a diverse regulatory factor by recognizing specific pseudopalindromic DNA motifs

Type II toxin-antitoxin (TA) systems modulate many essential cellular processes in prokaryotic organisms. Recent studies indicate certain type II antitoxins also transcriptionally regulate other genes, besides neutralizing toxin activity. Herein, we investigated the diverse transcriptional repression properties of type II TA antitoxin PaHigA from Pseudomonas aeruginosa. Biochemical and functional analyses showed that PaHigA recognized variable pseudopalindromic DNA sequences and repressed expression of multiple genes. Furthermore, we presented high resolution structures of apo-PaHigA, PaHigA-P_higBA and PaHigA-P_pa2440 complex, describing how the rearrangements of the HTH domain accounted for the different DNA-binding patterns among HigA homologues. Moreover, we demonstrated that the N-terminal loop motion of PaHigA was associated with its apo and DNA-bound states, reflecting a switch mechanism regulating HigA antitoxin function. Collectively, this work extends our understanding of how the PaHigB/HigA system regulates multiple metabolic pathways to balance the growth and stress response in P. aeruginosa and could guide further development of anti-TA oriented strategies for pathogen treatment.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

Evaluation on the use of Nanopore sequencing for direct characterization of coronaviruses from respiratory specimens, and a study on emerging missense mutations in partial RdRP gene of SARS-CoV-2

Coronavirus disease 2019 (COVID-19) pandemic has been a catastrophic burden to global healthcare systems. The fast spread of the etiologic agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), highlights the need to identify unknown coronaviruses rapidly for prompt clinical and public health decision making. Moreover, owing to the high mutation rate of RNA viruses, periodic surveillance on emerging variants of key virus components is essential for evaluating the efficacy of antiviral drugs, diagnostic assays and vaccines. These 2 knowledge gaps formed the basis of this study. In the first place, we evaluated the feasibility of characterizing coronaviruses directly from respiratory specimens. We amplified partial RdRP gene, a stable genetic marker of coronaviruses, from a collection of 57 clinical specimens positive for SARS-CoV-2 or other human coronaviruses, and sequenced the amplicons with Nanopore Flongle and MinION, the fastest and the most scalable massively-parallel sequencing platforms to-date. Partial RdRP sequences were successfully amplified and sequenced from 82.46% (47/57) of specimens, ranging from 75 to 100% by virus type, with consensus accuracy of 100% compared with Sanger sequences available (n = 40). In the second part, we further compared 19 SARS-CoV-2 RdRP sequences collected from the first to third waves of COVID-19 outbreak in Hong Kong with 22,173 genomes from GISAID EpiCoV™ database. No single nucleotide variants (SNVs) were found in our sequences, and 125 SNVs were observed from global data, with 56.8% being low-frequency (n = 1-47) missense mutations affecting the rear part of RNA polymerase. Among the 9 SNVs found on 4 conserved domains, the frequency of 15438G > T was highest (n = 34) and was predominantly found in Europe. Our data provided a glimpse into the sequence diversity of a primary antiviral drug and diagnostic target. Further studies are warranted to investigate the significance of these mutations.

Analysis of the small chromosomal Prionium serratum (Cyperid) demonstrates the importance of reliable methods to differentiate between mono- and holocentricity

For a long time, the Cyperid clade (Thurniceae-Juncaceae-Cyperaceae) was considered a group of species possessing holocentromeres exclusively. The basal phylogenetic position of Prionium serratum (Thunb.) Drège (Thurniceae) within Cyperids makes this species an important specimen to understand the centromere evolution within this clade. In contrast to the expectation, the chromosomal distribution of the centromere-specific histone H3 (CENH3), alpha-tubulin and different centromere-associated post-translational histone modifications (H3S10ph, H3S28ph and H2AT120ph) demonstrate a monocentromeric organisation of P. serratum chromosomes. Analysis of the high-copy repeat composition resulted in the identification of two centromere-localised satellite repeats. Hence, monocentricity was the ancestral condition for the Juncaceae-Cyperaceae-Thurniaceae Cyperid clade, and holocentricity in this clade has independently arisen at least twice after differentiation of the three families, once in Juncaceae and the other one in Cyperaceae. In this context, methods suitable for the identification of holocentromeres are discussed.

Genetic examination of the marine bacterium Pseudoalteromonas luteoviolacea and effects of its metamorphosis-inducing factors

Pseudoalteromonas luteoviolacea is a globally distributed marine bacterium that stimulates the metamorphosis of marine animal larvae, an important bacteria-animal interaction that can promote the recruitment of animals to benthic ecosystems. Recently, different P. luteoviolacea isolates have been shown to produce two stimulatory factors that can induce tubeworm and coral metamorphosis; Metamorphosis-Associated Contractile structures (MACs) and tetrabromopyrrole (TBP) respectively. However, it remains unclear what proportion of P. luteoviolacea isolates possess the genes encoding MACs, and what phenotypic effect MACs and TBP have on other larval species. Here, we show that 9 of 19 sequenced P. luteoviolacea genomes genetically encode both MACs and TBP. While P. luteoviolacea biofilms producing MACs stimulate the metamorphosis of the tubeworm Hydroides elegans, TBP biosynthesis genes had no effect under the conditions tested. Although MACs are lethal to larvae of the cnidarian Hydractinia symbiologicarpus, P. luteoviolacea mutants unable to produce MACs are capable of stimulating metamorphosis. Our findings reveal a hidden complexity of interactions between a single bacterial species, the factors it produces and two species of larvae belonging to different phyla.

Tissue-Specific Transcriptome Analysis Reveals Candidate Transcripts Associated with the Process of Programmed B Chromosome Elimination in Aegilops speltoides

Some eukaryotes exhibit dramatic genome size differences between cells of different organs, resulting from programmed elimination of chromosomes. Here, we present the first transcriptome analysis of programmed chromosome elimination using laser capture microdissection (LCM)-based isolation of the central meristematic region of Aegilops speltoides embryos where B chromosome (B) elimination occurs. The comparative RNA-seq analysis of meristematic cells of embryos with (Bplus) and without Bs (B0) allowed the identification of 14,578 transcript isoforms (35% out of 41,615 analyzed transcript isoforms) that are differentially expressed during the elimination of Bs. A total of 2908 annotated unigenes were found to be up-regulated in Bplus condition. These genes are either associated with the process of B chromosome elimination or with the presence of B chromosomes themselves. GO enrichment analysis categorized up-regulated transcript isoforms into 27 overrepresented terms related to the biological process, nine terms of the molecular function aspect and three terms of the cellular component category. A total of 2726 annotated unigenes were down-regulated in Bplus condition. Based on strict filtering criteria, 341 B-unique transcript isoforms could be identified in central meristematic cells, of which 70 were functionally annotated. Beside others, genes associated with chromosome segregation, kinetochore function and spindle checkpoint activity were retrieved as promising candidates involved in the process of B chromosome elimination.

A functional m6 A-RNA methylation pathway in the oyster Crassostrea gigas assumes epitranscriptomic regulation of lophotrochozoan development

N⁶ -methyladenosine (m⁶ A) is a prevalent epitranscriptomic mark in eukaryotic RNA, with crucial roles for mammalian and ecdysozoan development. Indeed, m⁶ A-RNA and the related protein machinery are important for splicing, translation, maternal-to-zygotic transition and cell differentiation. However, to date, the presence of an m⁶ A-RNA pathway remains unknown in more distant animals, questioning the evolution and significance of the epitranscriptomic regulation. Therefore, we investigated the m⁶ A-RNA pathway in the oyster Crassostrea gigas, a lophotrochozoan model whose development was demonstrated under strong epigenetic influence. Using mass spectrometry and dot blot assays, we demonstrated that m⁶ A-RNA is actually present in the oyster and displays variations throughout early oyster development, with the lowest levels at the end of cleavage. In parallel, by in silico analyses, we were able to characterize at the molecular level a complete and conserved putative m⁶ A machinery. The expression levels of the identified putative m⁶ A writers, erasers and readers were strongly regulated across oyster development. Finally, RNA pull-down coupled to LC-MS/MS allowed us to prove the actual presence of readers able to bind m⁶ A-RNA and exhibiting specific developmental patterns. Altogether, our results demonstrate the conservation of a complete m⁶ A-RNA pathway in the oyster and strongly suggest its implication in early developmental processes including MZT. This first demonstration and characterization of an epitranscriptomic regulation in a lophotrochozoan model, potentially involved in the embryogenesis, bring new insights into our understanding of developmental epigenetic processes and their evolution.

Genome-Scale Metabolic Model of Xanthomonas phaseoli pv. manihotis: An Approach to Elucidate Pathogenicity at the Metabolic Level

Xanthomonas phaseoli pv. manihotis (Xpm) is the causal agent of cassava bacterial blight, the most important bacterial disease in this crop. There is a paucity of knowledge about the metabolism of Xanthomonas and its relevance in the pathogenic process, with the exception of the elucidation of the xanthan biosynthesis route. Here we report the reconstruction of the genome-scale model of Xpm metabolism and the insights it provides into plant-pathogen interactions. The model, iXpm1556, displayed 1,556 reactions, 1,527 compounds, and 890 genes. Metabolic maps of central amino acid and carbohydrate metabolism, as well as xanthan biosynthesis of Xpm, were reconstructed using Escher (https://escher.github.io/) to guide the curation process and for further analyses. The model was constrained using the RNA-seq data of a mutant of Xpm for quorum sensing (QS), and these data were used to construct context-specific models (CSMs) of the metabolism of the two strains (wild type and QS mutant). The CSMs and flux balance analysis were used to get insights into pathogenicity, xanthan biosynthesis, and QS mechanisms. Between the CSMs, 653 reactions were shared; unique reactions belong to purine, pyrimidine, and amino acid metabolism. Alternative objective functions were used to demonstrate a trade-off between xanthan biosynthesis and growth and the re-allocation of resources in the process of biosynthesis. Important features altered by QS included carbohydrate metabolism, NAD(P)+ balance, and fatty acid elongation. In this work, we modeled the xanthan biosynthesis and the QS process and their impact on the metabolism of the bacterium. This model will be useful for researchers studying host-pathogen interactions and will provide insights into the mechanisms of infection used by this and other Xanthomonas species.

A high-throughput genomic screen identifies a role for the plasmid-borne type II secretion system of Escherichia coli O157:H7 (Sakai) in plant-microbe interactions

Shiga-toxigenic Escherichia coli (STEC) is often transmitted into food via fresh produce plants, where it can cause disease. To identify early interaction factors for STEC on spinach, a high-throughput positive-selection system was used. A bacterial artificial chromosome (BAC) clone library for isolate Sakai was screened in four successive rounds of short-term (2 h) interaction with spinach roots, and enriched loci identified by microarray. A Bayesian hierarchical model produced 115 CDS credible candidates, comprising seven contiguous genomic regions. Of the two candidate regions selected for functional assessment, the pO157 plasmid-encoded type two secretion system (T2SS) promoted interactions, while a chaperone-usher fimbrial gene cluster (loc6) did not. The T2SS promoted bacterial binding to spinach and appeared to involve the EtpD secretin protein. Furthermore, the T2SS genes, etpD and etpC, were expressed at a plant-relevant temperature of 18 °C, and etpD was expressed in planta by E. coli Sakai on spinach plants.

A novel bacterial thiosulfate oxidation pathway provides a new clue about the formation of zero-valent sulfur in deep sea

Zero-valent sulfur (ZVS) has been shown to be a major sulfur intermediate in the deep-sea cold seep of the South China Sea based on our previous work, however, the microbial contribution to the formation of ZVS in cold seep has remained unclear. Here, we describe a novel thiosulfate oxidation pathway discovered in the deep-sea cold seep bacterium Erythrobacter flavus 21–3, which provides a new clue about the formation of ZVS. Electronic microscopy, energy-dispersive, and Raman spectra were used to confirm that E. flavus 21–3 effectively converts thiosulfate to ZVS. We next used a combined proteomic and genetic method to identify thiosulfate dehydrogenase (TsdA) and thiosulfohydrolase (SoxB) playing key roles in the conversion of thiosulfate to ZVS. Stoichiometric results of different sulfur intermediates further clarify the function of TsdA in converting thiosulfate to tetrathionate (⁻O₃S–S–S–SO₃⁻), SoxB in liberating sulfone from tetrathionate to form ZVS and sulfur dioxygenases (SdoA/SdoB) in oxidizing ZVS to sulfite under some conditions. Notably, homologs of TsdA, SoxB, and SdoA/SdoB widely exist across the bacteria including in Erythrobacter species derived from different environments. This strongly indicates that this novel thiosulfate oxidation pathway might be frequently used by microbes and plays an important role in the biogeochemical sulfur cycle in nature.

Pathogenomes of Atypical Non-shigatoxigenic Escherichia coli NSF/SF O157:H7/NM: Comprehensive Phylogenomic Analysis Using Closed Genomes

The toxigenic conversion of Escherichia coli strains by Shiga toxin-converting (Stx) bacteriophages were prominent and recurring events in the stepwise evolution of enterohemorrhagic E. coli (EHEC) O157:H7 from an enteropathogenic (EPEC) O55:H7 ancestor. Atypical, attenuated isolates have been described for both non-sorbitol fermenting (NSF) O157:H7 and SF O157:NM serotypes, which are distinguished by the absence of Stx, the characteristic virulence hallmark of Stx-producing E. coli (STEC). Such atypical isolates either never acquired Stx-phages or may have secondarily lost stx during the course of infection, isolation, or routine subculture; the latter are commonly referred to as LST (Lost Shiga Toxin)-isolates. In this study we analyzed the genomes of 15 NSF O157:H7 and SF O157:NM strains from North America, Europe, and Asia that are characterized by the absence of stx, the virulence hallmark of STEC. The individual genomic basis of the Stx (-) phenotype has remained largely undetermined as the majority of STEC genomes in public genome repositories were generated using short read technology and are in draft stage, posing a major obstacle for the high-resolution whole genome sequence typing (WGST). The application of LRT (long-read technology) sequencing provided us with closed genomes, which proved critical to put the atypical non-shigatoxigenic NSF O157:H7 and SF O157:NM strains into the phylogenomic context of the stepwise evolutionary model. Availability of closed chromosomes for representative Stx (-) NSF O157:H7 and SF O157:NM strains allowed to describe the genomic basis and individual evolutionary trajectories underlying the absence of Stx at high accuracy and resolution. The ability of LRT to recover and accurately assemble plasmids revealed a strong correlation between the strains' featured plasmid genotype and chromosomally inferred clade, which suggests the coevolution of the chromosome and accessory plasmids. The identified ancestral traits in the pSFO157 plasmid of NSF O157:H7 strain LSU-61 provided additional evidence for its intermediate status. Taken together, these observations highlight the utility of LRTs for advancing our understanding of EHEC O157:H7/NM pathogenome evolution. Insights into the genomic and phenotypic plasticity of STEC on a lineage- and genome-wide scale are foundational to improve and inform risk assessment, biosurveillance, and prevention strategies.

Identification and RNAi Profile of a Novel Iflavirus Infecting Senegalese Aedes vexans arabiensis Mosquitoes

The inland floodwater mosquito Aedes vexans (Meigen, 1830) is a competent vector of numerous arthropod-borne viruses such as Rift Valley fever virus (Phenuiviridae) and Zika virus (Flaviviridae). Aedes vexans spp. have widespread Afrotropical distribution and are common European cosmopolitan mosquitoes. We examined the virome of Ae. vexans arabiensis samples from Barkédji village, Senegal, with small RNA sequencing, bioinformatic analysis, and RT-PCR screening. We identified a novel 9494 nt iflavirus (Picornaviridae) designated here as Aedes vexans iflavirus (AvIFV). Annotation of the AvIFV genome reveals a 2782 amino acid polyprotein with iflavirus protein domain architecture and typical iflavirus 5' internal ribosomal entry site and 3' poly-A tail. Aedes vexans iflavirus is most closely related to a partial virus sequence from Venturia canescens (a parasitoid wasp) with 56.77% pairwise amino acid identity. Analysis of AvIFV-derived small RNAs suggests that AvIFV is targeted by the exogenous RNA interference pathway but not the PIWI-interacting RNA response, as ~60% of AvIFV reads corresponded to 21 nt Dicer-2 virus-derived small RNAs and the 24-29 nt AvIFV read population did not exhibit a "ping-pong" signature. The RT-PCR screens of archival and current (circa 2011-2020) Ae. vexans arabiensis laboratory samples and wild-caught mosquitoes from Barkédji suggest that AvIFV is ubiquitous in these mosquitoes. Further, we screened wild-caught European Ae. vexans samples from Germany, the United Kingdom, Italy, and Sweden, all of which tested negative for AvIFV RNA. This report provides insight into the diversity of commensal Aedes viruses and the host RNAi response towards iflaviruses.

The proinflammatory cytokines TNF-α and IL-6 in lumpfish (Cyclopterus lumpus L.) -identification, molecular characterization, phylogeny and gene expression analyses

The proinflammatory cytokines TNF-α and IL-6 are important mediators of inflammatory reactions and orchestrators of the immune system in vertebrate. In this study, we have identified TNF-α and IL-6 in lumpfish, molecular characterized them at mRNA and gene level, performed homology modelling and measured their gene expression in different tissues and upon in vitro stimulation. A comprehensive phylogenetic analysis of TNF-α teleost sequences give novel insight into the TNF -α biology. Interestingly, we identified two isoforms of luIL-6. In normal tissue and leukocyte, the level of luTNF-α transcripts was higher than luIL-6. The expression pattern were parallel, except for brain, eye and gonad, and they displayed a similar induction pattern upon exposure to PAMPs, being most highly upregulated by flagellin. This is the first in-depth characterization of TNF and IL-6 in lumpfish. In recent years, lumpfish has become an important species for the aquaculture industry and establishment of qPCR-assays of luTNF-α and luIL-6 provide a valuable tool to measure effect of immune modulation, such as vaccination, microbiological disease and physiological trials. Lumpfish is also interesting for comparative studies as it represent a phylogenetic group that is poorly described immunologically.

Design of Bacterial Strain-Specific qPCR Assays Using NGS Data and Publicly Available Resources and Its Application to Track Biocontrol Strains

Biological control is emerging as a feasible alternative to chemical pesticides in agriculture. Measuring the microbial biocontrol agent (mBCA) populations in the environment is essential for an accurate environmental and health risk assessment and for optimizing the usage of an mBCA-based plant protection product. We hereby show a workflow to obtain a large number of qPCR markers suitable for robust strain-specific quantification. The workflow starts from whole genome sequencing data and consists of four stages: (i) identifying the strain-specific sequences, (ii) designing specific primer/probe sets for qPCR, and (iii) empirically verifying the performance of the assays. The first two stages involve exclusively computer work, but they are intended for researchers with little or no bioinformatic background: Only a knowledge of the BLAST suite tools and work with spreadsheets are required; a familiarity with the Galaxy environment and next-generation sequencing concepts are strongly advised. All bioinformatic work can be implemented using publicly available resources and a regular desktop computer (no matter the operating system) connected to the Internet. The workflow was tested with five bacterial strains from four different genera under development as mBCAs and yielded thousands of candidate markers and a triplex qPCR assay for each candidate mBCA. The qPCR assays were successfully tested in soils of different natures, water from different sources, and with samples from different plant tissues. The mBCA detection limits and population dynamics in the different matrices are similar to those in qPCR assays designed by other means. In summary, a new accessible, cost-effective, and robust workflow to obtain a large number of strain-specific qPCR markers is presented.

Chromosome-Level Assembly of Drosophila bifasciata Reveals Important Karyotypic Transition of the X Chromosome

The Drosophila obscura species group is one of the most studied clades of Drosophila and harbors multiple distinct karyotypes. Here we present a de novo genome assembly and annotation of D. bifasciata, a species which represents an important subgroup for which no high-quality chromosome-level genome assembly currently exists. We combined long-read sequencing (Nanopore) and Hi-C scaffolding to achieve a highly contiguous genome assembly approximately 193 Mb in size, with repetitive elements constituting 30.1% of the total length. Drosophila bifasciata harbors four large metacentric chromosomes and the small dot, and our assembly contains each chromosome in a single scaffold, including the highly repetitive pericentromeres, which were largely composed of Jockey and Gypsy transposable elements. We annotated a total of 12,821 protein-coding genes and comparisons of synteny with D. athabasca orthologs show that the large metacentric pericentromeric regions of multiple chromosomes are conserved between these species. Importantly, Muller A (X chromosome) was found to be metacentric in D. bifasciata and the pericentromeric region appears homologous to the pericentromeric region of the fused Muller A-AD (XL and XR) of pseudoobscura/affinis subgroup species. Our finding suggests a metacentric ancestral X fused to a telocentric Muller D and created the large neo-X (Muller A-AD) chromosome ∼15 MYA. We also confirm the fusion of Muller C and D in D. bifasciata and show that it likely involved a centromere-centromere fusion.

Cooler temperatures slow the repair of DNA damage in tadpoles exposed to ultraviolet radiation: Implications for amphibian declines at high altitude

Ultraviolet B radiation (UVBR) damages the DNA of exposed cells, causing dimers to form between adjacent pyrimidine nucleotides. These dimers block DNA replication, causing mutations and apoptosis. Most organisms utilize biochemical or biophysical DNA repair strategies to restore DNA structure; however, as with most biological reactions, these processes are likely to be thermally sensitive. Tadpoles exposed to elevated UVBR at low environmental temperatures have significantly higher rates of mortality and developmental deformities compared with tadpoles exposed to the same levels of UVBR at higher environmental temperatures. We hypothesized that low environmental temperatures impair the primary enzymatic (photolyase) DNA repair pathway in amphibians, leading to the accumulation of DNA damage. To test this hypothesis, we compared DNA repair rates and photolyase gene expression patterns in Limnodynastes peronii. Tadpoles were acutely exposed to UVBR for 1 hr at either 20 or 30°C, and we measured DNA damage and photolyase expression levels at intervals following this exposure. Temperature had a significant effect on the rate of DNA repair, with repair at 30°C occurring twice as fast as repair at 20°C. Photolyase gene expression (6-4 PP and CPD) was significantly upregulated by UVBR exposure, with expression levels increasing within 6 hr of UVBR exposure. CPD expression levels were not significantly affected by temperature, but 6-4 PP expression was significantly higher in tadpoles in the 30°C treatment within 12 hr of UVBR exposure. These data support the hypothesis that DNA repair rates are thermally sensitive in tadpoles and may explain why enigmatic amphibian declines are higher in montane regions where UVBR levels are naturally elevated and environmental temperatures are lower.