NCBI BLAST+ integrated into Galaxy

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

The Families of Non-LTR Transposable Elements within Neritimorpha and Other Gastropoda

Repeated sequences, especially transposable elements (TEs), are known to be abundant in some members of the important invertebrate class Gastropoda. TEs that do not have long terminal repeated sequences (non-LTR TEs) are frequently the most abundant type but have not been well characterised in any gastropod. Despite this, sequences in draft gastropod genomes are often described as non-LTR TEs, but without identification to family type. This study was conducted to characterise non-LTR TEs in neritimorph snails, using genomic skimming surveys of three species and the recently published draft genome of Theodoxus fluviatilis. Multiple families of non-LTR TEs from the I, Jockey, L1, R2 and RTE superfamilies were found, although there were notably few representatives of the first of these, which is nevertheless abundant in other Gastropoda. Phylogenetic analyses of amino acid sequences of the reverse transcriptase domain from the elements ORF2 regions found considerable interspersion of representatives of the four neritimorph taxa within non-LTR families and sub-families. In contrast, phylogenetic analyses of sequences from the elements' ORF1 region resolved the representatives from individual species as monophyletic. However, using either region, members of the two species of the Neritidae were closely related, suggesting their potential for investigation of phyletic evolution at the family level.

Orco-dependent survival of odorant receptor neurons in ants

Olfaction is essential for complex social behavior in insects. To discriminate complex social cues, ants evolved an expanded number of odorant receptor (Or) genes. Mutations in the obligate odorant co-receptor gene orco lead to the loss of ~80% of the antennal lobe glomeruli in the jumping ant Harpegnathos saltator. However, the cellular mechanism remains unclear. Here, we demonstrate massive apoptosis of odorant receptor neurons (ORNs) in the mid to late stages of pupal development, possibly due to ER stress in the absence of Orco. Further bulk and single-nucleus transcriptome analysis shows that, although most orco-expressing ORNs die in orco mutants, a small proportion of them survive: They express ionotropic receptor (Ir) genes that form IR complexes. In addition, we found that some Or genes are expressed in mechanosensory neurons and nonneuronal cells, possibly due to leaky regulation from nearby non-Or genes. Our findings provide a comprehensive overview of ORN development and Or expression in H. saltator.

A Case of Food-Borne Salmonellosis in a Corn Snake (Pantherophis guttatus) after a Feeder Mouse Meal

Reptiles are usually asymptomatic carriers of Salmonella, with the manifestation of typical clinical signs of acute forms in adult and non-immunocompromised animals being considered exceptions. In the present case, an adult male corn snake (Pantherophis guttatus) was found dead due to septic shock 48 h after consuming a feeder mouse purchased online. The snake's tissue samples and faeces were cultured for bacteria isolation. Microbiological examinations of the snake and mouse livers revealed the presence of Salmonella enterica subsp. enterica serovar Midway. A whole-genome analysis of these two isolates showed a high correlation between them: they belonged to the strain type ST-357 for the classic MLST scheme and to the strain type ST 171322 for the cgMLST scheme. Also, a virulence gene analysis revealed the presence of stdB and STM3026 genes. This report conveys a case of food-borne salmonellosis in a pet snake, transmitted from a feeder mouse, likely responsible for the snake's death due to septic shock. It highlights the relevance of feeder mice as a source of Salmonella infections in snakes and the associated risks to human health.

Virome analyses of Amblyomma cajennense and Rhipicephalus microplus ticks collected in Colombia

Tick-borne viruses (TBV) have gained public health relevance in recent years due to the recognition of human-associated fatal cases and the increase in tick-borne disease and transmission. However, many tick species have not been studied for their potential to transmit pathogenic viruses, especially those found in Latin America. To gain better understanding of the tick virome, we conducted targeted amplification using broadly-reactive consensus-degenerate pan-viral targeting viruses from the genera Flavivirus, Bandavirus, Uukuvirus, and Orthonairovirus genus. Additionally, we conducted unbiased metagenomic analyses to investigate the presence of viral RNA sequences in Amblyomma cajennense, A. patinoi and Rhipicephalus microplus ticks collected from a horse slaughter plant in Medellín, Colombia. While no viral products were detected by PCR, results of the metagenomic analyses revealed the presence of viral genomes belonging to the genera Phlebovirus, Bandavirus, and Uukuvirus, including Lihan Tick Virus (LTV), which was previously reported in Rhipicephalus microplus from Colombia. Overall, the results emphasized the enormous utility of the next-generation sequencing in identifying virus genetic diversity presents in ticks and other species of vectors and reservoirs.

A Review of Web-Based Metagenomics Platforms for Analysing Next-Generation Sequence Data

Metagenomics has now evolved as a promising technology for understanding the microbial population in the environment. By metagenomics, a number of extreme and complex environment has been explored for their microbial population. Using this technology, researchers have brought out novel genes and their potential characteristics, which have robust applications in food, pharmaceutical, scientific research, and other biotechnological fields. A sequencing platform can provide a sequence of microbial populations in any given environment. The sequence needs to be analysed computationally to derive meaningful information. It is presumed that only bioinformaticians with extensive computational skills can process the sequencing data till the downstream end. However, numerous open-source software and online servers are available to analyse the metagenomic data developed for a biologist with less computational skills. This review is focused on bioinformatics tools such as Galaxy, CSI-NGS portal, ANASTASIA and SHAMAN, EBI- metagenomics, IDseq, and MG-RAST for analysing metagenomic data.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant

BACKGROUND

Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported.

RESULTS

In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted.

CONCLUSIONS

This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Integrated Proteomics and Metabolomics of Safflower Petal Wilting and Seed Development

Safflower (Carthamus tinctorius L.) is an ancient oilseed crop of interest due to its diversity of end-use industrial and food products. Proteomic and metabolomic profiling of its organs during seed development, which can provide further insights on seed quality attributes to assist in variety and product development, has not yet been undertaken. In this study, an integrated proteome and metabolic analysis have shown a high complexity of lipophilic proteins and metabolites differentially expressed across organs and tissues during seed development and petal wilting. We demonstrated that these approaches successfully discriminated safflower reproductive organs and developmental stages with the identification of 2179 unique compounds and 3043 peptides matching 724 unique proteins. A comparison between cotyledon and husk tissues revealed the complementarity of using both technologies, with husks mostly featuring metabolites (99%), while cotyledons predominantly yielded peptides (90%). This provided a more complete picture of mechanisms discriminating the seed envelope from what it protected. Furthermore, we showed distinct molecular signatures of petal wilting and colour transition, seed growth, and maturation. We revealed the molecular makeup shift occurring during petal colour transition and wilting, as well as the importance of benzenoids, phenylpropanoids, flavonoids, and pigments. Finally, our study emphasizes that the biochemical mechanisms implicated in the growing and maturing of safflower seeds are complex and far-reaching, as evidenced by AraCyc, PaintOmics, and MetaboAnalyst mapping capabilities. This study provides a new resource for functional knowledge of safflower seed and potentially further enables the precision development of novel products and safflower varieties with biotechnology and molecular farming applications.

A combined phylogenetic strategy illuminates the evolution of Goniodorididae nudibranchs (Mollusca, Gastropoda, Heterobranchia)

Goniodorididae is a family of small dorid nudibranchs distributed worldwide that feed on entoprocts, ascidians, and bryozoans. The evolutionary relationships between its taxa have been uncertain due to the limited taxa available for phylogenetic analyses; some genera being paraphyletic. The family includes a remarkable number of synonymized genera in which the species richness is unequally distributed, while some genera have dozens of species others are monospecific. Some clades are very uniform morphologically while others are considered highly variable. To increase backbone phylogenetic resolution a target enrichment approach of ultra-conserved elements was aimed at representative Goniodorididae species for the first time. Additionally, we increase species representation by including mitochondrial markers cytochrome c oxidase subunit I and ribosomal RNA 16S as well as nuclear Histone 3 and ribosomal RNA 18S from 109 Goniodorididae species, out of approximately 160 currently valid species. Maximum likelihood and Bayesian inference analyses were performed to infer the phylogeny of the family. As a result, two subfamilies and eleven genera were elucidated. The synonymized genera Bermudella, Cargoa, and Ceratodoris are here resurrected and a new genus, Naisdoris gen. nov., is described. The clades included taxa with shared prey preference, showing that trophic behavior could have driven species evolution and morphological uniqueness within the family Goniodorididae.

Growth, Enzymatic, and Transcriptomic Analysis of xyr1 Deletion Reveals a Major Regulator of Plant Biomass-Degrading Enzymes in Trichoderma harzianum

The regulation of plant biomass degradation by fungi is critical to the carbon cycle, and applications in bioproducts and biocontrol. Trichoderma harzianum is an important plant biomass degrader, enzyme producer, and biocontrol agent, but few putative major transcriptional regulators have been deleted in this species. The T. harzianum ortholog of the transcriptional activator XYR1/XlnR/XLR-1 was deleted, and the mutant strains were analyzed through growth profiling, enzymatic activities, and transcriptomics on cellulose. From plate cultures, the Δxyr1 mutant had reduced growth on D-xylose, xylan, and cellulose, and from shake-flask cultures with cellulose, the Δxyr1 mutant had ~90% lower β-glucosidase activity, and no detectable β-xylosidase or cellulase activity. The comparison of the transcriptomes from 18 h shake-flask cultures on D-fructose, without a carbon source, and cellulose, showed major effects of XYR1 deletion whereby the Δxyr1 mutant on cellulose was transcriptionally most similar to the cultures without a carbon source. The cellulose induced 43 plant biomass-degrading CAZymes including xylanases as well as cellulases, and most of these had massively lower expression in the Δxyr1 mutant. The expression of a subset of carbon catabolic enzymes, other transcription factors, and sugar transporters was also lower in the Δxyr1 mutant on cellulose. In summary, T. harzianum XYR1 is the master regulator of cellulases and xylanases, as well as regulating carbon catabolic enzymes.

Probiotic Potential and Safety Assessment of Lactiplantibacillus plantarum cqf-43 and Whole-Genome Sequence Analysis

This study reports the whole-genome sequence of Lactiplantibacillus plantarum cqf-43 isolated from healthy sow feces. Based on genomic analysis, we performed a comprehensive safety assessment of strain cqf-43, using both in vitro and in vivo experiments, and explored its probiotic potential. The total genome length measures 3,169,201 bp, boasting a GC content of 44.59%. Through phylogenetic analyses, leveraging both 16S rRNA gene and whole-genome sequences, we confidently categorize strain cqf-43 as a member of Lactiplantibacillus. Genome annotation using Prokka unveiled a total of 3141 genes, encompassing 2990 protein-coding sequences, 71 tRNAs, 16 rRNAs, and 1 tmRNA. Functional annotations derived from COG and KEGG databases highlighted a significant abundance of genes related to metabolism, with a notable emphasis on carbohydrate utilization. The genome also revealed the presence of prophage regions and CRISPR-Cas regions while lacking virulence and toxin genes. Screening for antibiotic resistance genes via the CARD database yielded no detectable transferable resistance genes, effectively eliminating the potential for harmful gene transfer. It is worth highlighting that the virulence factors identified via the VFDB database primarily contribute to bolstering pathogen resilience in hostile environments. This characteristic is particularly advantageous for probiotics. Furthermore, the genome is devoid of menacing genes such as hemolysin, gelatinase, and biogenic amine-producing genes. Our investigation also unveiled the presence of three unannotated secondary metabolite biosynthetic gene clusters, as detected by the online tool antiSMASH, suggesting a great deal of unknown potential for this strain. Rigorous in vitro experiments confirmed tolerance of strain cqf-43 in the intestinal environment, its antimicrobial efficacy, sensitivity to antibiotics, absence of hemolysis and gelatinase activity, and its inability to produce biogenic amines. In addition, a 28-day oral toxicity test showed that the strain cqf-43 did not pose a health hazard in mice, further establishing it as a safe strain.

Generation of de novo miRNAs from template switching during DNA replication

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.

A metatranscriptomic evaluation of viruses in field-collected bed bugs

Cimex lectularius, known as the common bed bug, is a widespread hematophagous human ectoparasite and urban pest that is not known to be a vector of any human infectious disease agents. However, few studies in the era of molecular biology have profiled the microorganisms harbored by field populations of bed bugs. The objective of this study was to examine the viruses present in a large sampling of common bed bugs and related bat bugs (Cimex pipistrelle). RNA sequencing was undertaken on an international sampling of > 500 field-collected bugs, and multiple workflows were used to assemble contigs and query these against reference nucleotide databases to identify viral genomes. Shuangao bed bug virus 2, an uncharacterized rhabdovirus previously discovered in Cimex hemipterus from China, was found in several bed bug pools from the USA and Europe, as well as in C. pipistrelle, suggesting that this virus is common among bed bug populations. In addition, Shuangao bed bug virus 1 was detected in a bed bug pool from China, and sequences matching Enterobacteria phage P7 were found in all bed bug pools, indicating the ubiquitous presence of phage-derived elements in the genome of the bed bug or its enterobacterial symbiont. However, viral diversity was low in bed bugs in our study, as no other viral genomes were detected with significant coverage. These results provide evidence against frequent virus infection in bed bugs. Nonetheless, our investigation had several important limitations, and additional studies should be conducted to better understand the prevalence and composition of viruses in bed bugs. Most notably, our study largely focused on insects from urban areas in industrialized nations, thus likely missing infrequent virus infections and those that could occur in rural or tropical environments or developing nations.

Neolamarckia cadamba hosts a putative novel deltapartitivirus: a revelation by transcriptome data-mining

Neolamarckia cadamba (Roxb.) Bosser is a medicinally important, fast-growing, timber-yielding tree species. In the present study, the virome of N. cadamba was explored using the publicly available N. cadamba transcriptome datasets and a putative novel virus, tentatively named as cadamba cryptic virus 1 (CdbCV1), was identified. CdbCV1 contained two genome segments, each coding for a single protein. CdbCV1 RNA1 (1564 nt) encoded for an RNA dependent RNA polymerase (RdRp) protein while CdbCV1 RNA2 (1492 nt) encoded for a coat protein (CP). Phylogenetic and sequence similarity analyses revealed the relatedness of CdbCV1 to pepper cryptic virus 1 and pittosporum cryptic virus 1. Based on the species demarcation criteria, genome organization and phylogeny, CdbCV1 can be regarded a new member of the genus Deltapartitivirus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-023-00845-8.

GhMYC1374 regulates the cotton defense response to cotton aphids by mediating the production of flavonoids and free gossypol

Myelocytomatosis (MYC) transcription factors (TFs) in plants are well-known regulators of plant defense against herbivores. However, the role and mechanism of MYC TFs in cotton (Gossypium hirsutum L.) defense against cotton aphids (Aphis gossypii Glover) remain still elusive. Herein, on the basis of aphid-induced cotton transcriptome analysis, GhMYC1374, a cotton MYC2-like TF that was highly induced by cotton aphid attack, has been identified that confers cotton aphid resistance in cotton. GhMYC1374 was an intranuclear transcription factor with three domains: bHLH-MYC_N, RBR and bHLH_AtAIB_like. GhMYC1374 was induced under cotton aphid feeding, exogenous methyl jasmonate (MeJA) and salicylic acid (SA) treatments. GhMYC1374 transient overexpression in cotton plants enhanced cotton aphid-resistance, while GhMYC1374 silence through VIGS (virus induced gene silencing) decreased cotton aphid-resistance. GhMYC1374 transient overexpression of in cotton plants activated the phenylpropane pathway and promoted the synthesis of flavonoids, and resistance to thus enhanced the cotton resistance against aphids. In contrast, GhMYC1374 silence inhibited the biosynthesis of flavonoids. In addition, GhMYC1374 also positively activated the expression of the biosynthetic genes of free gossypol, leading to the high content of free gossypol. Taken together, our results suggest that GhMYC1374 is involved in the cotton defense response against cotton aphids by regulating the biosynthesis of flavonoids and free gossypol.

A rabies-related lyssavirus from a Nycticeinops schlieffeni bat with neurological signs, South Africa

We report the coding-complete sequence of a lyssavirus, provisionally designated Phala bat lyssavirus (PBLV), characterized using a metagenomics approach. PBLV was identified in a Nycticeinops schlieffeni bat that exhibited neurological signs and died within 24 hours of admission to a wildlife rehabilitation center in Phalaborwa, South Africa.

Characterization and expression analysis of SnRK2, PYL, and ABF/ AREB/ ABI5 gene families in sweet potato

Abscisic acid (ABA) signaling in plants is essential to several aspects of plant development, such as tolerance to environmental stresses and growth. ABA signaling is also important for storage organ formation in crops, such as sweet potato. However, the repertoire of I. batatas ABA signaling gene families has not yet been fully characterized, so that it is unclear which members of these families are necessary for tuberization. Therefore, genome-wide identification of the sweet potato ABF/ AREB/ ABI5, SnRK2, and PYL gene families was performed, along with phylogenetic, motif, cis-regulatory element (CRE), and expression analyses. Nine ABF, eight SnRK2, and eleven PYL gene family members were identified, and there was high sequence conservation among these proteins that were revealed by phylogenetic and motif analyses. The promoter sequences of these genes had multiple CREs that were involved in hormone responses and stress responses. In silico and qRT-PCR expression analyses revealed that these genes were expressed in various tissues and that IbABF3, IbABF4, IbDPBF3, IbDPBF4, IbPYL4, IbSnRK2.1, and IbSnRK2.2 were significantly expressed during storage root development. These results are an important reference that can be used for functional validation studies to better understand how ABA signaling elicits storage root formation at the molecular level.

Characterisation of neonatal Staphylococcus capitis NRCS-A isolates compared with non NRCS-A Staphylococcus capitis from neonates and adults

Staphylococcus capitis is a frequent cause of late-onset sepsis in neonates admitted to Neonatal Intensive Care Units (NICU). One clone of S. capitis, NRCS-A has been isolated from NICUs globally although the reasons for the global success of this clone are not well understood.We analysed a collection of S. capitis colonising babies admitted to two NICUs, one in the UK and one in Germany as well as corresponding pathological clinical isolates. Genome analysis identified a population structure of three groups; non-NRCS-A isolates, NRCS-A isolates, and a group of 'proto NRCS-A' - isolates closely related to NRCS-A but not associated with neonatal infection. All bloodstream isolates belonged to the NRCS-A group and were indistinguishable from strains carried on the skin or in the gut. NRCS-A isolates showed increased tolerance to chlorhexidine and antibiotics relative to the other S. capitis as well as enhanced ability to grow at higher pH values. Analysis of the pangenome of 138 isolates identified characteristic nsr and tarJ genes in both the NRCS-A and proto groups. A CRISPR-cas system was only seen in NRCS-A isolates which also showed enrichment of genes for metal acquisition and transport.We found evidence for transmission of S. capitis NRCS-A within NICU, with related isolates shared between babies and multiple acquisitions by some babies. Our data show NRCS-A strains commonly colonise uninfected babies in NICU representing a potential reservoir for potential infection. This work provides more evidence that adaptation to survive in the gut and on skin facilitates spread of NRCS-A, and that metal acquisition and tolerance may be important to the biology of NRCS-A. Understanding how NRCS-A survives in NICUs can help develop infection control procedures against this clone.

A genomic platform for surveillance and antigen discovery in Plasmodium spp. using long-read amplicon sequencing

Many vaccine candidate proteins in the malaria parasite Plasmodium falciparum are under strong immunological pressure and confer antigenic diversity. We present a sequencing and data analysis platform for the genomic surveillance of the insertion or deletion (indel)-rich antigens merozoite surface protein 1 (MSP1), MSP2, glutamate-rich protein (GLURP), and CSP from P. falciparum using long-read circular consensus sequencing (CCS) in multiclonal malaria isolates. Our platform uses 40 PCR primers per gene to asymmetrically barcode and identify multiclonal infections in pools of up to 384 samples. With msp2, we validated the method using 235 mock infections combining 10 synthetic variants at different concentrations and infection complexities. We applied this strategy to P. falciparum isolates from a longitudinal cohort in Tanzania. Finally, we constructed an analysis pipeline that streamlines the processing and interpretation of epidemiological and antigenic diversity data from demultiplexed FASTQ files. This platform can be easily adapted to other polymorphic antigens of interest in Plasmodium or any other human pathogen.

Climate has contributed to population diversification of Daphnia galeata across Eurasia

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

Characterization and comparative analysis of the Escherichia marmotae M-12 isolate from bank vole (Myodes glareolus)

The Escherichia marmotae is a bacterium of the Enterobacterales order, which was first isolated from the Himalayan marmot (Marmota himalayana). Recently E. marmotae has been shown to cause severe infections in humans. Wild animals were suggested to be a natural reservoir of this bacterium. The present study describes the first case of E. marmotae isolation from an apparently healthy wild bank vole (Myodes glareolus). Phenotype, as well as genotype-based techniques, were applied to characterize E. marmotae M-12 isolate. E. marmotae M-12 had the capsule-positive phenotype, high adhesion to human erythrocytes and HEp-2 cells as well as a low invasion into HEp-2 cells. E. marmotae M-12 was avirulent in mice. The phylogenomic analyses of E. marmotae showed dispersed phylogenetic structure among isolates of different origins. Virulome analysis of M-12 isolate revealed the presence of the following factors: siderophores, heme uptake systems, capsule synthesis, curli and type I fimbriae, flagella proteins, OmpA porin, etc. Comparative virulome analysis among available E. marmotae genomes revealed the presence of capsule K1 genes mostly in pathogenic isolates and OmpA porin presence among all strains. We assume that the K1 capsule and OmpA porin play a key role in the virulence of E. marmotae. Pathogenesis of the latter might be similar to extraintestinal pathogenic E. coli.

De novo transcriptome analysis of Dysoxylum binectariferum to unravel the biosynthesis of pharmaceutically relevant specialized metabolites

The tropical tree, D. binectariferum, is a prominent source of chromone alkaloid rohitukine, which is used in the semi-syntheses of anticancer molecules such as flavopiridol and P-276-00. The biosynthetic pathway of rohitukine or its derivatives is currently unknown in plants. Here, we explored chromone alkaloid biosynthesis in D. binectariferum through targeted transcriptome sequencing. Illumina sequencing of leaves and roots of a year-old D. binectariferum seedling generated, 42.43 and 38.74 million paired-end short reads, respectively. Quality filtering and de novo assembly of the transcriptome generated 274,970 contigs and 126,788 unigenes with an N50 contig length of 1560 bp. The assembly generated 117,619 translated unigene protein sequences and 51,598 non-redundant sequences. Nearly 80% of these non-redundant sequences were annotated to publicly available protein and nucleotide databases, suggesting the completeness and effectiveness of the transcriptome assembly. Using the assembly, we identified a chalcone synthase (CHS) and three type III polyketide synthases (PKS-III; non-CHS type) that are likely to be involved in the biosynthesis of chromone ring/noreugenin moiety of rohitukine. We also identified key enzymes like lysine decarboxylase in the piperidine pathway that make the piperidine moiety of rohitukine. Besides these, the upstream enzymes in flavonoid biosynthesis like phenylalanine ammonia-lyase (PAL), trans-cinnamate 4-hydroxylase (C4H),4-coumarate-CoA ligase (4CL), and chalcone isomerase (CHI) have also been identified. Also, terpene synthases that are likely to be involved in the biosynthesis of various terpenoid scaffolds have been identified. Together, the D. binectariferum transcriptome resource forms a basis for further exploration of biosynthetic pathways of these valuable compounds through functional validation of the candidate genes and metabolic engineering in heterologous hosts. Additionally, the transcriptome dataset generated will serve as an important resource for research on functional genomics and enzyme discovery in D. binectariferum and comparative analysis with other Meliaceae family members.

m6A Profile Dynamics Indicates Regulation of Oyster Development by m6A-RNA Epitranscriptomes

The N6-methylation of RNA adenosines (N6-methyladenosine, m6A) is an important regulator of gene expression with critical implications in vertebrate and insect development. However, the developmental significance of epitranscriptomes in lophotrochozoan organisms remains unknown. Using methylated RNA immunoprecipitation sequencing (MeRIP-seq), we generated transcriptome-wide m6A-RNA methylomes covering the entire development of the oyster from oocytes to juveniles. Oyster RNA classes display specific m6A signatures, with messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) exhibiting distinct profiles and being highly methylated compared to transposable element (TE) transcripts. Epitranscriptomes are dynamic and correspond to the chronological steps of development (cleavage, gastrulation, organogenesis, and metamorphosis), with minimal mRNA and lncRNA methylation at the morula stage followed by a global increase. mRNA m6A levels are correlated with transcript levels, and shifts in methylation profiles correspond to expression kinetics. Differentially methylated transcripts cluster according to embryo-larval stages and bear the corresponding developmental functions (cell division, signal transduction, morphogenesis, and cell differentiation). The m6A level of TE transcripts is also regulated and peaks during the gastrulation. We demonstrate that m6A-RNA methylomes are dynamic and associated with gene expression regulation during oyster development. The putative epitranscriptome implication in the cleavage, maternal-to-zygotic transition, and cell differentiation in a lophotrochozoan model brings new insights into the control and evolution of developmental processes.

Discovery of putative long non-coding RNAs expressed in the eyes of Astyanax mexicanus (Actinopterygii: Characidae)

Astyanax mexicanus is a well-known model species, that has two morphotypes, cavefish, from subterranean rivers and surface fish, from surface rivers. They are morphologically distinct due to many troglomorphic traits in the cavefish, such as the absence of eyes. Most studies on A. mexicanus are focused on eye development and protein-coding genes involved in the process. However, lncRNAs did not get the same attention and very little is known about them. This study aimed to fill this knowledge gap, identifying, describing, classifying, and annotating lncRNAs expressed in the embryo's eye tissue of cavefish and surface fish. To do so, we constructed a concise workflow to assemble and evaluate transcriptomes, annotate protein-coding genes, ncRNAs families, predict the coding potential, identify putative lncRNAs, map them and predict interactions. This approach resulted in the identification of 33,069 and 19,493 putative lncRNAs respectively mapped in cavefish and surface fish. Thousands of these lncRNAs were annotated and identified as conserved in human and several species of fish. Hundreds of them were validated in silico, through ESTs. We identified lncRNAs associated with genes related to eye development. This is the case of a few lncRNAs associated with sox2, which we suggest being isomorphs of the SOX2-OT, a lncRNA that can regulate the expression of sox2. This work is one of the first studies to focus on the description of lncRNAs in A. mexicanus, highlighting several lncRNA targets and opening an important precedent for future studies focusing on lncRNAs expressed in A. mexicanus.

Potato psyllids mount distinct gut responses against two different 'Candidatus Liberibacter solanacearum' haplotypes

'Candidatus Liberibacter solanacearum' (Lso) is a bacterial pathogen infecting several crops and causing damaging diseases. Several Lso haplotypes have been identified. Among the seven haplotypes present in North America, LsoA and LsoB are transmitted by the potato psyllid, Bactericera cockerelli (Šulc), in a circulative and persistent manner. The gut, which is the first organ pathogen encounters, could be a barrier for Lso transmission. However, the molecular interactions between Lso and the psyllid vector at the gut interface remain largely unknown. In this study, we investigated the global transcriptional responses of the adult psyllid gut upon infection with two Lso haplotypes (LsoA and LsoB) using Illumina sequencing. The results showed that each haplotype triggers a unique transcriptional response, with most of the distinct genes elicited by the highly virulent LsoB. The differentially expressed genes were mainly associated with digestion and metabolism, stress response, immunity, detoxification as well as cell proliferation and epithelium renewal. Importantly, distinct immune pathways were triggered by LsoA and LsoB in the gut of the potato psyllid. The information in this study will provide an understanding of the molecular basis of the interactions between the potato psyllid gut and Lso, which may lead to the discovery of novel molecular targets for the control of these pathogens.

Comparitive Analysis of the Chloroplast Genomes of Three Houpoea Plants

The genus Houpoea belongs to the family Magnoliaceae, and the species in this genus have important medicinal values. However, the investigation of the correlation between the evolution of the genus and its phylogeny has been severely hampered by the unknown range of species within the genus and the paucity of research on its chloroplast genome. Thus, we selected three species of Houpoea: Houpoea officinalis var officinalis (OO), Houpoea officinalis var. biloba (OB), and Houpoea rostrata (R). With lengths of 160,153 bp (OO), 160,011 bp (OB), and 160,070 bp (R), respectively, the whole chloroplast genomes (CPGs) of these three Houpoea plants were acquired via Illumina sequencing technology, and the findings were annotated and evaluated. These three chloroplast genomes were revealed by the annotation findings to be typical tetrads. A total of 131, 132, and 120 different genes were annotated. The CPGs of the three species had 52, 47, and 56 repeat sequences, which were primarily found in the ycf2 gene. A useful tool for identifying species is the approximately 170 simple sequence repeats (SSRs) that have been found. The border area of the reverse repetition region (IR) was studied, and it was shown that across the three Houpoea plants, it is highly conservative, with only changes between H. rostrata and the other two plants observed. Numerous highly variable areas (rps3-rps19, rpl32-trnL, ycf1, ccsA, etc.) have the potential to serve as the barcode label for Houpoea, according to an examination of mVISTA and nucleotide diversity (Pi). Phylogenetic relation indicates that Houpoea is a monophyletic taxon, and its genus range and systematic position are consistent with the Magnoliaceae system of Sima Yongkang-Lu Shugang, including five species and varieties of H. officinalis var. officinalis, H. rostrata, H. officinalis var. biloba, Houpoea obovate, and Houpoea tripetala, which evolved and differentiated from the ancestors of Houpoea to the present Houpoea in the above order. This study provides valuable information on the genus Houpoea, enriches the CPG information on Houpoea genus, and provides genetic resources for the further classification of and phylogenetic research on Houpoea.

Rapid defense mechanism suppression during viral- oomycete disease complex formation

Combined infection of the host plant with pathogens involving different parasitic lifestyles may result in synergistic effects that intensify disease symptoms. Understanding the molecular dynamics during concurrent infection provides essential insight into the host response. The transcriptomic pattern of cucumber plants infected with a necrotrophic pathogen, Pythium spinosum, and a biotrophic pathogen, Cucumber green mottle mosaic virus (CGMMV) was studied at different time points, under regimes of single and co-infection. Analysis of CGMMV infection alone revealed a mild influence on host gene expression at the stem base, while the infection by P. spinosum is associated with drastic changes in gene expression. Comparing P. spinosum as a single infecting pathogen with a later co-infection by CGMMV revealed a rapid host response as early as 24 hours post-CGMMV inoculation with a sharp downregulation of genes related to the host defense mechanism against the necrotrophic pathogen. Suppression of the defense mechanism of co-infected plants was followed by severe stress, including 30% plants mortality and an increase of the P. spinosum hyphae. The first evidence of defense recovery against the necrotrophic pathogen only occurred 13 days post-viral infection. These results support the hypothesis that the viral infection of the Pythium pre-infected plants subverted the host defense system and changed the equilibrium obtained with P. spinosum. It also implies a time window in which the plants are most susceptible to P. spinosum after CGMMV infection.

Identification and Profiling of a Novel Bombyx mori latent virus Variant Acutely Infecting Helicoverpa armigera and Trichoplusia ni

Insect cell expression systems are increasingly being used in the medical industry to develop vaccines against diseases such as COVID-19. However, viral infections are common in these systems, making it necessary to thoroughly characterize the viruses present. One such virus is Bombyx mori latent virus (BmLV), which is known to be specific to Bombyx mori and to have low pathogenicity. However, there has been little research on the tropism and virulence of BmLV. In this study, we examined the genomic diversity of BmLV and identified a variant that persistently infects Trichoplusia ni-derived High Five cells. We also assessed the pathogenicity of this variant and its effects on host responses using both in vivo and in vitro systems. Our results showed that this BmLV variant causes acute infections with strong cytopathic effects in both systems. Furthermore, we characterized the RNAi-based immune response in the T. ni cell line and in Helicoverpa armigera animals by assessing the regulation of RNAi-related genes and profiling the generated viral small RNAs. Overall, our findings shed light on the prevalence and infectious properties of BmLV. We also discuss the potential impact of virus genomic diversity on experimental outcomes, which can help interpret past and future research results.

Phylogenomics and morphological evolution of the mega-diverse genus Artemisia (Asteraceae: Anthemideae): implications for its circumscription and infrageneric taxonomy

BACKGROUND AND AIMS

Artemisia is a mega-diverse genus consisting of ~400 species. Despite its medicinal importance and ecological significance, a well-resolved phylogeny for global Artemisia, a natural generic delimitation and infrageneric taxonomy remain missing, owing to the obstructions from limited taxon sampling and insufficient information on DNA markers. Its morphological characters, such as capitulum, life form and leaf, show marked variations and are widely used in its infrageneric taxonomy. However, their evolution within Artemisia is poorly understood. Here, we aimed to reconstruct a well-resolved phylogeny for global Artemisia via a phylogenomic approach, to infer the evolutionary patterns of its key morphological characters and to update its circumscription and infrageneric taxonomy.

METHODS

We sampled 228 species (258 samples) of Artemisia and its allies from both fresh and herbarium collections, covering all the subgenera and its main geographical areas, and conducted a phylogenomic analysis based on nuclear single nucleotide polymorphisms (SNPs) obtained from genome skimming data. Based on the phylogenetic framework, we inferred the possible evolutionary patterns of six key morphological characters widely used in its previous taxonomy.

KEY RESULTS

The genus Kaschgaria was revealed to be nested in Artemisia with strong support. A well-resolved phylogeny of Artemisia consisting of eight highly supported clades was recovered, two of which were identified for the first time. Most of the previously recognized subgenera were not supported as monophyletic. Evolutionary inferences based on the six morphological characters showed that different states of these characters originated independently more than once.

CONCLUSIONS

The circumscription of Artemisia is enlarged to include the genus Kaschgaria. The morphological characters traditionally used for the infrageneric taxonomy of Artemisia do not match the new phylogenetic tree. They experienced a more complex evolutionary history than previously thought. We propose a revised infrageneric taxonomy of the newly circumscribed Artemisia, with eight recognized subgenera to accommodate the new results.

Transcriptomic analysis of cave, surface, and hybrid samples of the isopod Asellus aquaticus and identification of chromosomal location of candidate genes for cave phenotype evolution

BACKGROUND

Transcriptomic methods can be used to elucidate genes and pathways responsible for phenotypic differences between populations. Asellus aquaticus is a freshwater isopod crustacean with surface- and cave-dwelling ecomorphs that differ greatly in multiple phenotypes including pigmentation and eye size. Multiple genetic resources have been generated for this species, but the genes and pathways responsible for cave-specific characteristics have not yet been identified. Our goal was to generate transcriptomic resources in tandem with taking advantage of the species' ability to interbreed and generate hybrid individuals.

RESULTS

We generated transcriptomes of the Rakov Škocjan surface population and the Rak Channel of Planina Cave population that combined Illumina short-read assemblies and PacBio Iso-seq long-read sequences. We investigated differential expression at two different embryonic time points as well as allele-specific expression of F₁ hybrids between cave and surface individuals. RNAseq of F₂ hybrids, as well as genotyping of a backcross, allowed for positional information of multiple candidate genes from the differential expression and allele-specific analyses.

CONCLUSIONS

As expected, genes involved in phototransduction and ommochrome synthesis were under-expressed in the cave samples as compared to the surface samples. Allele-specific expression analysis of F₁ hybrids identified genes with cave-biased (cave allele has higher mRNA levels than the surface allele) and surface-biased expression (surface allele has higher mRNA levels than the cave allele). RNAseq of F₂ hybrids allowed for multiple genes to be placed to previously mapped genomic regions responsible for eye and pigmentation phenotypes. In the future, these transcriptomic resources will guide prioritization of candidates for functional analysis.

Identification of a putative novel cholivirus in the transcriptome of Gymnema sylvestre R. Br

Gymnema sylvestre is a tropical climber species that is widely used in traditional medicine since ages. In the present study, the transcriptome datasets of G. sylvestre available in public domain were screened for the presence of novel plant viral sequences and a putative novel virus tentatively named as Gymnema sylvestre virus 1 (GysV1) was identified. Coding-complete genome segments of GysV1 that are 6.35 kb (RNA1) and 3.98 kb (RNA2) long possessed a single large open reading frame coding for a polyprotein. BLASTp, sequence identity and phylogenetic analyses revealed the relatedness of GysV1 to the members of the subgenus Cholivirus (genus Sadwavirus; family Secoviridae; order Picornavirales). Based on the species demarcation criteria of the family Secoviridae, GysV1 can be regarded as a new cholivirus member.

Broadening the host range and genetic diversity of waikaviruses

Waikaviruses are monopartite, positive sense, single-stranded RNA viruses that cause economically important plant diseases. Despite their importance, waikaviruses are poorly understood and only ten members are currently recognized. The present study on Sequence Read Archive (SRA)-based data-driven virus discovery (DDVD) identified 22 putative new waikaviruses, nearly doubling the number of known waikaviruses, in SRA libraries of diverse plant species, from ferns to trees. Besides, a highly divergent secoviral sequence with distinct genome features was identified in a wheat transcriptome. Other significant findings of the study include identification of a new waikavirus in a library derived from diseased water chestnut sample wherein a caulimovirus was reported, prediction of coiled-coils in hypothetical protein region of waikaviral polyprotein alignment and phylogenetic clustering of tree-infecting waikaviruses. The study not only reiterates the importance of DDVD in unveiling hitherto hidden viral sequences in plant SRA libraries but also deepens our understanding of waikaviral diversity.

Structural genome variants of Pseudomonas aeruginosa clone C and PA14 strains

Plasticity of Pseudomonas aeruginosa chromosomes is mainly driven by an extended accessory genome that is shaped by insertion and deletion events. Further modification of the genome composition can be induced by chromosomal inversion events which lead to relocation of genes in the affected genomic DNA segments, modify the otherwise highly conserved core genome synteny and could even alter the location of the replication terminus. Although the genome of the first sequenced strain, PAO1, displayed such a large genomic inversion, knowledge on such recombination events in the P. aeruginosa population is limited. Several large inversions had been discovered in the late 1990s in cystic fibrosis isolates of the major clonal lineage C by physical genome mapping, and subsequent work on these examples led to the characterization of the DNA at the recombination breakpoints and a presumed recombination mechanism. Since then, the topic was barely addressed in spite of the compilation of thousands of P. aeruginosa genome sequences that are deposited in databases. Due to the use of second-generation sequencing, genome contig assembly had usually followed synteny blueprints provided by the existing reference genome sequences. Inversion detection was not feasible by these approaches, as the respective read lengths did not allow reliable resolution of sequence repeats that are typically found at the borders of inverted segments. In this study, we applied PacBio and MinION long-read sequencing to isolates of the mentioned clone C collection. Confirmation of inversions predicted from the physical mapping data demonstrated that unbiased sequence assembly of such read datasets allows the detection of genomic inversions and the resolution of the recombination breakpoint regions. Additional long-read sequencing of representatives of the other major clonal lineage, PA14, revealed large inversions in several isolates, from cystic fibrosis origin as well as from other sources. These findings indicated that inversion events are not restricted to strains from chronic infection background, but could be widespread in the P. aeruginosa population and contribute to genome plasticity. Moreover, the monitored examples emphasized the role of small mobile DNA units, such as IS elements or transposons, and accessory DNA elements in the inversion-related recombination processes.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Barley RIC157, a potential RACB scaffold protein, is involved in susceptibility to powdery mildew

CRIB motif-containing barley RIC157 is a novel ROP scaffold protein that interacts directly with barley RACB, promotes susceptibility to fungal penetration, and colocalizes with RACB at the haustorial neck. Successful obligate pathogens benefit from host cellular processes. For the biotrophic ascomycete fungus Blumeria hordei (Bh) it has been shown that barley RACB, a small monomeric G-protein (ROP, Rho of plants), is required for full susceptibility to fungal penetration. The susceptibility function of RACB probably lies in its role in cell polarity, which may be co-opted by the pathogen for invasive ingrowth of its haustorium. However, how RACB supports fungal penetration success and which other host proteins coordinate this process is incompletely understood. RIC (ROP-Interactive and CRIB-(Cdc42/Rac Interactive Binding) motif-containing) proteins are considered scaffold proteins which can interact directly with ROPs via a conserved CRIB motif. Here we describe a previously uncharacterized barley RIC protein, RIC157, which can interact directly with RACB in planta. We show that, in the presence of constitutively activated RACB, RIC157 shows a localization at the cell periphery/plasma membrane, whereas it otherwise localizes to the cytoplasm. RIC157 appears to mutually stabilize the plasma membrane localization of the activated ROP. During fungal infection, RIC157 and RACB colocalize at the penetration site, particularly at the haustorial neck. Additionally, transiently overexpressed RIC157 renders barley epidermal cells more susceptible to fungal penetration. We discuss that RIC157 may promote fungal penetration into barley epidermal cells by operating probably downstream of activated RACB.

Changes in annual transcriptome dynamics of a clone of Japanese cedar (Cryptomeria japonica D. Don) planted under different climate conditions

Environmental responses are critical for plant growth and survival under different climate conditions. To elucidate the underlying biological mechanisms of environmental responses in Japanese cedar (Cryptomeria japonica D. Don), the annual transcriptome dynamics of common clonal trees (Godai1) planted at three different climate sites (Yamagata, Ibaraki, and Kumamoto Prefectures) were analyzed using microarrays. Both principal component analysis (PCA) and hierarchical clustering of the microarray data indicated the transition to dormant transcriptome status occurred earlier and the transition to active growth status later in the colder region. Interestingly, PCA also indicated that the transcriptomes of trees grown under three different conditions were similar during the growth period (June to September), whereas the transcriptomes differed between sites during the dormant period (January to March). In between-site comparisons, analyses of the annual expression profiles of genes for sites 'Yamagata vs. Kumamoto', 'Yamagata vs. Ibaraki', and 'Ibaraki vs. Kumamoto' identified 1,473, 1,137, and 925 targets exhibiting significantly different expression patterns, respectively. The total of 2,505 targets that exhibited significantly different expression patterns in all three comparisons may play important roles in enabling cuttings to adapt to local environmental conditions. Partial least-squares regression analysis and Pearson correlation coefficient analysis revealed that air temperature and day length were the dominant factors controlling the expression levels of these targets. GO and Pfam enrichment analyses indicated that these targets include genes that may contribute to environmental adaptation, such as genes related to stress and abiotic stimulus responses. This study provided fundamental information regarding transcripts that may play an important role in adaptation to environmental conditions at different planting sites.

S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence

Self-incompatibility in Prunus species is governed by a single locus consisting of two highly multi-allelic and tightly linked genes, one coding for an F-box protein-i.e., SFB in Prunus- controlling the pollen specificity and one coding for an S-RNase gene controlling the pistil specificity. Genotyping the allelic combination in a fruit tree species is an essential procedure both for cross-based breeding and for establishing pollination requirements. Gel-based PCR techniques using primer pairs designed from conserved regions and spanning polymorphic intronic regions are traditionally used for this task. However, with the great advance of massive sequencing techniques and the lowering of sequencing costs, new genotyping-by-sequencing procedures are emerging. The alignment of resequenced individuals to reference genomes, commonly used for polymorphism detection, yields little or no coverage in the S-locus region due to high polymorphism between different alleles within the same species, and cannot be used for this purpose. Using the available sequences of Japanese plum S-loci concatenated in a rosary-like structure as synthetic reference sequence, we describe a procedure to accurately genotype resequenced individuals that allowed the analysis of the S-genotype in 88 Japanese plum cultivars, 74 of them are reported for the first time. In addition to unraveling two new S-alleles from published reference genomes, we identified at least two S-alleles in 74 cultivars. According to their S-allele composition, they were assigned to 22 incompatibility groups, including nine new incompatibility groups reported here for the first time (XXVII-XXXV).

Characterisation of the Upper Respiratory Tract Virome of Feedlot Cattle and Its Association with Bovine Respiratory Disease

Bovine respiratory disease (BRD) is a major health problem within the global cattle industry. This disease has a complex aetiology, with viruses playing an integral role. In this study, metagenomics was used to sequence viral nucleic acids in the nasal swabs of BRD-affected cattle. The viruses detected included those that are well known for their association with BRD in Australia (bovine viral diarrhoea virus 1), as well as viruses known to be present but not fully characterised (bovine coronavirus) and viruses that have not been reported in BRD-affected cattle in Australia (bovine rhinitis, bovine influenza D, and bovine nidovirus). The nasal swabs from a case-control study were subsequently tested for 10 viruses, and the presence of at least one virus was found to be significantly associated with BRD. Some of the more recently detected viruses had inconsistent associations with BRD. Full genome sequences for bovine coronavirus, a virus increasingly associated with BRD, and bovine nidovirus were completed. Both viruses belong to the Coronaviridae family, which are frequently associated with disease in mammals. This study has provided greater insights into the viral pathogens associated with BRD and highlighted the need for further studies to more precisely elucidate the roles viruses play in BRD.

Application of Cloning-Free Genome Engineering to Escherichia coli

The propagation of foreign DNA in Escherichia coli is central to molecular biology. Recent advances have dramatically expanded the ability to engineer (bacterial) cells; however, most of these techniques remain time-consuming. The aim of the present work was to explore the possibility to use the cloning-free genome editing (CFGE) approach, proposed by Döhlemann and coworkers (2016), for E. coli genetics, and to deepen the knowledge about the homologous recombination mechanism. The E. coli auxotrophic mutant strains FB182 (hisF892) and FB181 (hisI903) were transformed with the circularized wild-type E. coli (i) hisF gene and hisF gene fragments of decreasing length, and (ii) hisIE gene, respectively. His⁺ clones were selected based on their ability to grow in the absence of histidine, and their hisF/hisIE gene sequences were characterized. CFGE method allowed the recombination of wild-type his genes (or fragments of them) within the mutated chromosomal copy, with a different recombination frequency based on the fragment length, and the generation of clones with a variable number of in tandem his genes copies. Data obtained pave the way to further evolutionary studies concerning the homologous recombination mechanism and the fate of in tandem duplicated genes.

A Live Cell Protein Complementation Assay for ORFeome-Wide Probing of Human HOX Interactomes

Biological pathways rely on the formation of intricate protein interaction networks called interactomes. Getting a comprehensive map of interactomes implies the development of tools that allow one to capture transient and low-affinity protein-protein interactions (PPIs) in live conditions. Here we presented an experimental strategy: the Cell-PCA (cell-based protein complementation assay), which was based on bimolecular fluorescence complementation (BiFC) for ORFeome-wide screening of proteins that interact with different bait proteins in the same live cell context, by combining high-throughput sequencing method. The specificity and sensitivity of the Cell-PCA was established by using a wild-type and a single-amino-acid-mutated HOXA9 protein, and the approach was subsequently applied to seven additional human HOX proteins. These proof-of-concept experiments revealed novel molecular properties of HOX interactomes and led to the identification of a novel cofactor of HOXB13 that promoted its proliferative activity in a cancer cell context. Taken together, our work demonstrated that the Cell-PCA was pertinent for revealing and, importantly, comparing the interactomes of different or highly related bait proteins in the same cell context.

DivBrowse-interactive visualization and exploratory data analysis of variant call matrices

BACKGROUND

The sequencing of whole genomes is becoming increasingly affordable. In this context, large-scale sequencing projects are generating ever larger datasets of species-specific genomic diversity. As a consequence, more and more genomic data need to be made easily accessible and analyzable to the scientific community.

FINDINGS

We present DivBrowse, a web application for interactive visualization and exploratory analysis of genomic diversity data stored in Variant Call Format (VCF) files of any size. By seamlessly combining BLAST as an entry point together with interactive data analysis features such as principal component analysis in one graphical user interface, DivBrowse provides a novel and unique set of exploratory data analysis capabilities for genomic biodiversity datasets. The capability to integrate DivBrowse into existing web applications supports interoperability between different web applications. Built-in interactive computation of principal component analysis allows users to perform ad hoc analysis of the population structure based on specific genetic elements such as genes and exons. Data interoperability is supported by the ability to export genomic diversity data in VCF and General Feature Format 3 files.

CONCLUSION

DivBrowse offers a novel approach for interactive visualization and analysis of genomic diversity data and optionally also gene annotation data by including features like interactive calculation of variant frequencies and principal component analysis. The use of established standard file formats for data input supports interoperability and seamless deployment of application instances based on the data output of established bioinformatics pipelines.

Pectobacterium carotovorum Phage vB_PcaM_P7_Pc Is a New Member of the Genus Certrevirus

Pectobacterium carotovorum is an economically important phytopathogen and has been identified as the major causative agent of bacterial soft rot in carrots. Control of this phytopathogen is vital to minimizing carrot harvest losses. As fully efficient control measures to successfully avoid the disease are unavailable, the phage-mediated biocontrol of the pathogen has recently gained scientific attention. In this study, we present a comprehensive characterization of the P. carotovorum phage vB_PcaM_P7_Pc (abbreviated as P7_Pc) that was isolated from infected carrot samples with characteristic soft rot symptoms, which were obtained from storage facilities at market places in Gampaha District, Sri Lanka. P7_Pc is a myovirus, and it exhibits growth characteristics of an exclusively lytic life cycle. It showed visible lysis against four of the tested P. carotovorum strains and one Pectobacterium aroidearum strain. This phage also showed a longer latent period (125 min) than other related phages; however, this did not affect its high phage titter (>1010 PFU/mL). The final assembled genome of P7_Pc is 147,299 bp in length with a G+C content of 50.34%. Of the 298 predicted open reading frames (ORFs) of the genome of P7_Pc, putative functions were assigned to 53 ORFs. Seven tRNA-coding genes were predicted in the genome, while the genome lacked any major genes coding for lysogeny-related products, confirming its virulent nature. The P7_Pc genome shares 96.12% and 95.74% average nucleotide identities with Cronobacter phages CR8 and PBES02, respectively. Phylogenetic and phylogenomic analyses of the genome revealed that P7_Pc clusters well within the clade with the members representing the genus Certrevirus. Currently, there are only 4 characterized Pectobacterium phages (P. atrosepticum phages phiTE and CB7 and Pectobacterium phages DU_PP_I and DU_PP_IV) that are classified under the genus, making the phage P7_Pc the first reported member of the genus isolated using the host bacterium P. carotovorum. The results of this study provide a detailed characterization of the phage P7_Pc, enabling its careful classification into the genus Certrevirus. The knowledge gathered on the phage based on the shared biology of the genus will further aid in the future selection of phage P7_Pc as a biocontrol agent. IMPORTANCE Bacterial soft rot disease, caused by Pectobacterium spp., can lead to significant losses in carrot yields. As current control measures involving the use of chemicals or antibiotics are not recommended in many countries, bacteriophage-mediated biocontrol strategies are being explored for the successful control of these phytopathogens. The successful implementation of such biocontrol strategies relies heavily upon the proper understanding of the growth characteristics and genomic properties of the phage. Further, the selection of taxonomically different phages for the formulation of phage cocktails in biocontrol applications is critical to combat potential bacterial resistance development. This study was conducted to carefully characterize and resolve the phylogenetic placement of the P. carotovorum phage vB_PcaM_P7_Pc by using its biological and genomic properties. Phage P7_Pc has a myovirus morphotype with an exclusively lytic life cycle, and the absence of genes related to lysogeny, toxin production, and antibiotic resistance in its genome confirmed its suitability to be used in environmental applications. Furthermore, P7_Pc is classified under the genus Certrevirus, making it the first reported phage of the genus of the host species, P. carotovorum.

Transcriptome-wide modulation by Sargassum vulgare and Acanthophora spicifera extracts results in a prime-triggered plant signalling cascade in tomato and sweet pepper

Seaweed extracts (SWEs) are becoming integrated into crop production systems due to their multiple beneficial effects including growth promotion and induction of defence mechanisms. However, the comprehensive molecular mechanisms of these effects are yet to be elucidated. The current study investigated the transcriptomic changes induced by SWEs derived from Sargassum vulgare and Acanthophora spicifera on tomato and sweet pepper plants. Tomato and sweet pepper plants were subjected to foliar treatment with alkaline extracts prepared from the above seaweeds. Transcriptome changes in the plants were assessed 72 h after treatments using RNA sequencing. The treated plants were also analysed for defence enzyme activities, nutrient composition and phytohormonal profiles. The results showed the significant enrichment of genes associated with several growth and defence processes including photosynthesis, carbon and nitrogen metabolism, plant hormone signal transduction, plant-pathogen interaction, secondary metabolite metabolism, MAPK signalling and amino acid biosynthesis. Activities of defence enzymes were also significantly increased in SWE-treated plants. Plant nutrient profiling showed significant increases in calcium, potassium, nitrogen, sulphur, boron, copper, iron, manganese, zinc and phosphorous levels in SWE-treated plants. Furthermore, the levels of auxins, cytokinins and gibberellins were also significantly increased in the treated plants. The severity of bacterial leaf spot and early blight incidence in plants treated with SWE was significantly reduced, in addition to other effects like an increase in chlorophyll content, plant growth, and fruit yield. The results demonstrated the complex effect of S. vulgare and A. spicifera extracts on the plants' transcriptome and provided evidence of a strong role of these extracts in increasing plant growth responses while priming the plants against pathogenic attack simultaneously. The current study contributes to the understanding of the molecular mechanisms of SWEs in plants and helps their usage as a viable organic input for sustainable crop production.

In silico analysis of key regulatory networks related to microfibril angle in Populus trichocarpa Hook

Dissection of regulatory network that control wood structure is highly challenging in functional genomics. Nevertheless, due to the availability of genomic, transcriptomic and proteomic sequences, a large amount of information is available for use in achieving this goal. MicroRNAs, which compose a class of small non-coding RNA molecules that inhibit protein translation by targeting mRNA cleavage sites and thus regulate a wide variety of developmental and physiological processes in plants, are important parts of this regulatory network. These findings and the availability of sequence information have made it possible to carry out an in silico analysis to predict and annotate miRNAs and their target genes associated with an important factor affecting wood rigidity, microfibril angle (MFA), throughout the Populus trichocarpa Hook. genome. Our computational approach revealed miRNAs and their targets via ESTs, sequences putatively associated with microfibril angle. In total, 250 miRNAs were identified as RNA molecules with roles in the silencing and post-transcriptional regulation of the expression of nine genes. We found SHY2, IAA4 (ATAUX2–11), BZIP60, AP2, MYB15, ABI3, MYB17, LAF1 and MYB28 as important nodes in a network with possible role in MFA determination. Other co-expressed genes putatively involved in this regulatory system were also identified by construction of a co-expression network. The candidate genes from this study may help unravel the regulatory networks putatively linked to microfibril angle.

First report of the whole genome of Moraxella bovoculi genotype 1 from India and comparative genomics of Moraxella bovoculi to identify genotype-specific markers

Moraxella bovoculi has been isolated frequently from cattle with Infectious bovine keratoconjunctivitis (IBK). Two diverse genotypes of M. bovoculi, 1 and 2 were identified based on whole genome sequence analysis. It is essential to discriminate between the two genotypes to frame prevention and control measures. The whole genome of M. bovoculi TN7 was sequenced and compared to other M. bovoculi strains available in the NCBI database. M. bovoculi TN7 was found to be genotype 1, had an RTX toxin operon and pilA gene that are the known virulence factors in related Moraxella sp., but lacked antimicrobial resistance genes. M. bovoculi was found to have an open pangenome with 4051 (75.31%) accessory genes, and the addition of each new genome adds 18 genes to the pangenome. Comparison of pilin protein amino acid sequences revealed three new sequence types. Furthermore, the presence of linx, nagL, swrC and mdtA genes was found to be genotype 1 specific, whereas hyaD, garR, gbsA, yhdG, gabT, iclR, higB2, hmuU, hmuT and hemS were found only in genotype 2. Polymerase Chain Reaction (PCR) primers were designed and evaluated on strain TN7 plus seven additional strains accessible to us that had not been whole genome sequenced. This initial evaluation of the designed primers for the linX and hyaD genes produced the expected banding patterns on PCR gels for genotypes 1 and 2, respectively, among the 8 strains. The genotype-specific genes identified in this study can be used as markers for accurate diagnosis of genotype 1 isolates and this can aid in the development of autogenous or other molecular vaccines for treatment of infectious bovine keratoconjunctivitis (IBK) in resource-limited research settings.

Analysis of public domain plant transcriptomes expands the phylogenetic diversity of the family Secoviridae

Secoviruses are mono-/bipartite plant-infecting, icosahedral RNA viruses that incite economically important diseases in plants. In the present study, nine secoviruses tentatively named as Ananas comosus secovirus (AcSV), Artocarpus altilis secovirus (AaSV), Boehmeria nivea secovirus (BnSV), Gynostemma pentaphyllum secovirus (GpSV), Orobanche cernua secovirus (OcSV), Paris polyphylla secovirus 1 (PpSV1), Paris polyphylla secovirus 2 (PpSV2), Rhododendron delavayi secovirus (RdSV), and Yucca gloriosa secovirus (YgSV) were identified by probing publicly available transcriptomes of eight plant species. Coding-complete genome/genome segments of all the identified viruses encoding a polyprotein were recovered. Two of the nine identified viruses-AcSV and GpSV were discovered in few of the small RNA libraries of respective plant species. Putative cleavage sites were predicted in polyproteins encoded by AcSV, GpSV, PpSV2 and YgSV genome segments. Phylogenetic and sequence identity analyses revealed that AcSV, GpSV and YgSV, PpSV1 and RdSV putatively belong to the genera- Sadwavirus (sub genus: Cholivirus), Fabavirus, Nepovirus and Waikavirus, respectively, while AaSV, BnSV, and PpSV2 may represent a distinct group of viruses within the family Secoviridae as they could not conclusively be assigned to a single genus.

Development of dim-light vision in the nocturnal reef fish family Holocentridae I: retinal gene expression

Developmental changes to the visual systems of animals are often associated with ecological shifts. Reef fishes experience a change in habitat between larval life in the shallow open ocean to juvenile and adult life on the reef. Some species also change their lifestyle over this period and become nocturnal. While these ecological transitions are well documented, little is known about the ontogeny of nocturnal reef fish vision. Here, we used transcriptomics to investigate visual development in 12 representative species from both subfamilies, Holocentrinae (squirrelfishes) and Myripristinae (soldierfishes), in the nocturnal coral reef fish family, Holocentridae. Results revealed that the visual systems of holocentrids are initially well adapted to photopic conditions with pre-settlement larvae having high levels of cone opsin gene expression and a broad cone opsin gene repertoire (8 genes). At reef settlement, holocentrids started to invest more in their scotopic visual system, and compared with adults, showed upregulation of genes involved in cell differentiation/proliferation. By adulthood, holocentrids had well developed scotopic vision with high levels of rod opsin gene expression, reduced cone opsin gene expression and repertoire (1–4 genes) and upregulated phototransduction genes. Finally, although the two subfamilies shared similar ecologies across development, their visual systems diverged after settlement, with Myripristinae investing more in scotopic vision than Holocentrinae. Hence, both ecology and phylogeny are likely to determine the development of the holocentrid visual system.

Summary: Coral reef fishes in the family Holocentridae remodel their retina at the molecular level to adapt to a nocturnal lifestyle during development.

A Chromosome Level Genome Assembly of a Winter Turnip Rape (Brassica rapa L.) to Explore the Genetic Basis of Cold Tolerance

Winter rapeseed (Brassica rapa L.) is an important overwintering oilseed crop that is widely planted in northwest China and suffers chronic low temperatures in winter. So the cold stress becomes one of the major constraints that limit its production. The currently existing genomes limit the understanding of the cold-tolerant genetic basis of rapeseed. Here we assembled a high-quality long-read genome of B. rapa "Longyou-7" cultivar, which has a cold-tolerant phenotype, and constructed a graph-based pan-genome to detect the structural variations within homologs of currently reported cold-tolerant related genes in the "Longyou-7" genome, which provides an additional elucidation of the cold-tolerant genetic basis of "Longyou-7" cultivar and promotes the development of cold-tolerant breeding in B. rapa.

Role of Staphylococcus agnetis and Staphylococcus hyicus in the Pathogenesis of Buffalo Fly Skin Lesions in Cattle

Buffalo flies (Haematobia irritansexigua) are hematophagous ectoparasites of cattle causing production and welfare impacts in northern Australian herds. Skin lesions associated with buffalo fly infestation and Stephanofilaria nematode infection are manifested as focal dermatitis or ulcerated areas, most commonly on the medial canthus of the eye, along the lateral and ventral neck, and on the abdomen of cattle. For closely related horn flies (Haematobia irritans irritans), Staphylococcus aureus has been suggested as a contributing factor in the development of lesions. To investigate the potential role of bacterial infection in the pathogenesis of buffalo fly lesions, swabs were taken from lesions and normal skin, and bacteria were also isolated from surface washings of buffalo flies and surface-sterilized homogenized flies. Bacterial identification was conducted by matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) and strain typing by repetitive sequence-based PCR (rep-PCR) and DNA sequencing to determine species similarity and virulence factors. Of 50 bacterial isolates collected from lesions, 38 were identified as Staphylococcus agnetis and 12 as Staphylococcus hyicus, whereas four isolates from normal skin were S. hyicus and one was Mammaliicoccus sciuri. Of the Staphylococcus isolates isolated from buffalo flies, five were identified as S. agnetis and three as S. hyicus. Fifty percent of the buffalo fly isolates had rep-PCR genotypic patterns identical to those of the lesion isolates. Genome sequencing of 16 S. agnetis and four S. hyicus isolates revealed closely similar virulence factor profiles, with all isolates possessing exfoliative toxin A and C genes. The findings from this study suggest the involvement of S. agnetis and S. hyicus in buffalo fly lesion pathogenesis. This should be taken into account in the development of effective treatment and control strategies for lesions.

IMPORTANCE Skin lesions in cattle associated with feeding by Haematobia fly species are a significant welfare issue in Australia, North and South America, and Europe. The development of these lesions has been attributed to a number of causal factors, but the exact etiology and pathogenesis were unclear. This study characterized Staphylococcus agnetis and Staphylococcus hyicus strains from cattle skin lesions and in vector flies and demonstrated their role in the pathogenesis of these lesions. These findings will aid the development of targeted and more effective treatment and control strategies for lesions associated with fly infestation in cattle.