GenBank

Molecular Phylogeny of the Family Cordulegastridae (Odonata) Worldwide

In this study, we present the first attempt at a molecular phylogenetic analysis of the entire family of Cordulegastridae involving 60% of its known species. Our analysis is in favor of reclassification of the members of the family into four genera: (i) the monophyletic genus Anotogaster Selys, 1854, with the number of known species reduced by three synonymizations; (ii) the genus Cordulegaster Leach in Brewster, 1815 including all members of the boltonii group and, as a preliminary solution, the American species C. virginiae Novelo-Gutiérrez, 2018 and, very tentatively, C. diadema Selys, 1868. The bidentata group forms a genus of its own, for which we restored the name Thecagaster Selys 1854, stat. rev. Cordulegaster pekinensis McLachlan in Selys, 1886, currently considered as Neallogaster pekinensis, was placed by us in Thecagaster as well. The genus Neallogaster Cowley, 1934 needs further investigation involving all remaining species listed in it. The genus Zoraena Kirby, 1890, stat. rev., was recovered to accommodate the remaining American species of Cordulegaster. We synonymized three species of Anotogaster: Anotogaster gregoryi Fraser, 1923 = Anotogaster xanthoptera Lohmann, 1993, syn. nov.; Anotogasterkuchenbeiseri (Förster, 1899) = Anotogaster antehumeralis Lohmann, 1993, syn. nov.; Anotogaster kuchenbeiseri (Förster, 1899) = Anotogaster cornutifrons Lohmann, 1993, syn. nov., based on examination of the existing type specimens. The type of specimens of A. klossi Fraser, 1919 = A. flaveola Lohmann 1993 syn. confirm., were also examined, and their synonymy was confirmed. The isolated populations of A. sieboldii (Selys, 1854) from the archipelagos of Okinawa and Amami Oshima in Japan, respectively, should be regarded as separate species, which will be described elsewhere. Furthermore, we suggest the synonymization of Cordulegaster parvistigma Selys 1873 syn. nov. with Thecagaster brevistigma (Selys 1854) comb. restaur.

Draft genome sequence of the deep-subsurface Ciceribacter sp. strain T2.26MG-112.2, a second Rhizobiaceae isolated from the Iberian Pyrite Belt at 492.6 mbs

T2.26MG-112.2 is a Ciceribacter strain that has been isolated from the deep subsurface of the Iberian Pyrite Belt. We report its draft genome consisting of a chromosome of ≈4.9 Mb and a plasmid of 357 kb. The annotation reveals 4,824 coding sequences, 48 tRNA genes, and 1 rRNA operon.

Omics technology draws a comprehensive heavy metal resistance strategy in bacteria

The rapid industrial revolution significantly increased heavy metal pollution, becoming a major global environmental concern. This pollution is considered as one of the most harmful and toxic threats to all environmental components (air, soil, water, animals, and plants until reaching to human). Therefore, scientists try to find a promising and eco-friendly technique to solve this problem i.e., bacterial bioremediation. Various heavy metal resistance mechanisms were reported. Omics technologies can significantly improve our understanding of heavy metal resistant bacteria and their communities. They are a potent tool for investigating the adaptation processes of microbes in severe conditions. These omics methods provide unique benefits for investigating metabolic alterations, microbial diversity, and mechanisms of resistance of individual strains or communities to harsh conditions. Starting with genome sequencing which provides us with complete and comprehensive insight into the resistance mechanism of heavy metal resistant bacteria. Moreover, genome sequencing facilitates the opportunities to identify specific metal resistance genes, operons, and regulatory elements in the genomes of individual bacteria, understand the genetic mechanisms and variations responsible for heavy metal resistance within and between bacterial species in addition to the transcriptome, proteome that obtain the real expressed genes. Moreover, at the community level, metagenome, meta transcriptome and meta proteome participate in understanding the microbial interactive network potentially novel metabolic pathways, enzymes and gene species can all be found using these methods. This review presents the state of the art and anticipated developments in the use of omics technologies in the investigation of microbes used for heavy metal bioremediation.

TcZC3HTTP, a regulatory element that contributes to Trypanosoma cruzi cell proliferation

Post-transcriptional regulation of gene expression is a critical process for adapting to and surviving Trypanosoma cruzi, a parasite with a complex life cycle. RNA-binding proteins (RBPs) are key players in this regulation, forming ribonucleoprotein complexes (messenger ribonucleoproteins) and RNA granules that control transcript stability, localization, degradation, and translation modulation. Understanding the specific roles of individual RBPs is crucial for unraveling the details of this regulatory network. In this study, we generated null mutants of the TcZC3HTTP gene, a specific RBP in the Trypanosoma family characterized by a C3H zinc finger and a DNAJ domain associated with RNA and protein binding, respectively. Through cell growth assays, we demonstrated that the absence of TcZC3HTTP or the expression of an additional tagged version impacted epimastigote growth, indicating its contribution to cell proliferation. TcZC3HTTP was found to associate with mRNAs involved in cell cycle and division in epimastigotes, while in nutritionally stressed parasites it exhibited associations with mRNAs coding for other RBPs and rRNA. Furthermore, our analysis identified that TcZC3HTTP protein partners were different during normal growth conditions compared to starvation conditions, with the latter showing enrichment of ribosomal proteins and other RBPs. Therefore, this study provides insights into TcZC3HTTP's role in the post-transcriptional regulation of gene expression during normal growth and nutritional stress in T. cruzi, uncovering its versatile functions in different cellular contexts.IMPORTANCEUnderstanding how Trypanosoma cruzi, the causative agent of Chagas disease, regulates gene expression is crucial for developing targeted interventions. In this study, we investigated the role of TcZC3HTTP, an RNA-binding protein, in post-transcriptional regulation. Our findings demonstrate that TcZC3HTTP is relevant for the growth and proliferation of epimastigotes, a stage of the parasite's life cycle. We identified its associations with specific mRNAs involved in cell cycle and division and its interactions with enzymes and other RNA-binding proteins (RBPs) under normal and starvation conditions. These insights shed light on the regulatory network underlying gene expression in T. cruzi and reveal the multifaceted functions of RBPs in this parasite. Such knowledge enhances our understanding of the parasite's biology and opens avenues for developing novel therapeutic strategies targeting post-transcriptional gene regulation in T. cruzi.

Digger: directed annotation of immunoglobulin and T cell receptor V, D, and J gene sequences and assemblies

SUMMARY

Knowledge of immunoglobulin and T cell receptor encoding genes is derived from high-quality genomic sequencing. High-throughput sequencing is delivering large volumes of data, and precise, high-throughput approaches to annotation are needed. Digger is an automated tool that identifies coding and regulatory regions of these genes, with results comparable to those obtained by current expert curational methods.

AVAILABILITY AND IMPLEMENTATION

Digger is published under open source license at https://github.com/williamdlees/Digger and is available as a Python package and a Docker container.

Small extracellular vesicles from irradiated lung epithelial cells promote the activation of fibroblasts in pulmonary fibrosis

BACKGROUND

Alveolar epithelial injury and dysfunction are the risk factors for radiation-induced pulmonary fibrosis (RIPF). However, it is not clear about the relationship between RIPF and the small extracellular vesicles (sEV) secreted by irradiated alveolar epithelial cells. Based on the activation of fibroblasts, this study explored the role of sEV derived from alveolar epithelial cells in RIPF and the potential mechanisms.

METHODS

Transmission electron microscopy (TEM), nanoparticle tracking analysis (NTA), and western blotting were used to characterize sEV. Western blotting was used to detect fibrosis-associated proteins. Cell counts and transwell assays were used to evaluate the proliferation and migration ability of fibroblasts. RT-PCR was used to observe the extracellular matrix (ECM) synthesized by fibroblasts, miRNA changes in the sEV were determined by second-generation sequencing.

RESULTS

TEM, NTA, and western blotting showed the extracellular vesicles with a double-layer membrane structure of approximately 100 nm in diameter. The sEV derived from irradiated A549, HBEC3-KT, and MLE12 cells upregulated FN1 and alpha-SMA proteins expression in fibroblasts and drove the fibroblast to myofibroblast transition, and the sEV from irradiated mouse bronchoalveolar lavage fluid (BALF) affirmed the same results. In addition, the sEV derived from irradiated alveolar epithelial cells significantly increased the migration ability of fibroblasts and the expression of extracellular matrix proteins such as FN1. The results of miRNA sequencing of sEV in BALF of rats with RIPF showed that the metabolic pathway may be important for miRNA to regulate the activation of fibroblasts.

CONCLUSION

The sEV derived from radiated pulmonary epithelial cells promote the activation, migration and extracellular matrix proteins expression of lung fibroblasts; miRNA in sEV may be an important molecular that affects the activation of lung fibroblasts.

Characterization of an EPS-producing bifidobacterial strain based on integration of phenotypic and complete genome sequencing data

Bifidobacterium and Lactobacillus are known to be common members of the human intestinal microbiota, which play important roles in maintaining the homeostasis of host gut microenvironment. Several bifidobacterial and lactobacilli strains have been used as probiotics for health benefits. The exopolysaccharides (EPSs) produced by strains from Bifidobacterium and Lactobacillus are considered as beneficial traits mediating these beneficial effects. In this study, 21 strains belonging to Bifidobacterium and Lactobacillus were isolated from healthy infants' stool and were screened for EPS-producing ability. Among these strains, Bifidobacterium longum XZM1 showed the highest EPS productivity, which was further confirmed and characterized. The complete genome of strain XZM1 was sequenced, which revealed the presence of a gene cluster for EPS production. Furthermore, comparative genome analysis was performed among XZM1 and other strains from B. longum species. Following purification, the molecular weight (Mw) of EPS from XZM1 was determined as 4023 Da (Mw) through gel permeation chromatography. Analysis of the EPS hydrolysates revealed that the EPS was composed of mannose, glucose, galactose, arabinose, and fucose. Additionally, the EPS exhibited higher scavenging abilities toward hydroxyl than 1,1-diphenyl-2-picrylhydrazyl free radical. Overall, these results suggest that XZM1 from B. longum species may be a promising probiotic candidate.

Non-coding RNAs identification and regulatory networks in pathogen-host interaction in the microsporidia congenital infection

BACKGROUND

The interaction networks between coding and non-coding RNAs (ncRNAs) including long non-coding RNA (lncRNA), covalently closed circular RNA (circRNA) and miRNA are significant to elucidate molecular processes of biological activities and interactions between host and pathogen. Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-feeding industry. However, little is known about ncRNAs that take place in the microsporidia congenital infection. Here we conducted whole-transcriptome RNA-Seq analyses to identify ncRNAs and regulatory networks for both N. bombycis and host including silkworm embryos and larvae during the microsporidia congenital infection.

RESULTS

A total of 4,171 mRNAs, 403 lncRNA, 62 circRNAs, and 284 miRNAs encoded by N. bombycis were identified, among which some differentially expressed genes formed cross-talk and are involved in N. bombycis proliferation and infection. For instance, a lncRNA/circRNA competing endogenous RNA (ceRNA) network including 18 lncRNAs, one circRNA, and 20 miRNAs was constructed to describe 14 key parasites genes regulation, such as polar tube protein 3 (PTP3), ricin-B-lectin, spore wall protein 4 (SWP4), and heat shock protein 90 (HSP90). Regarding host silkworm upon N. bombycis congenital infection, a total of 14,889 mRNAs, 3,038 lncRNAs, 19,039 circRNAs, and 3,413 miRNAs were predicted based on silkworm genome with many differentially expressed coding and non-coding genes during distinct developmental stages. Different species of RNAs form interacting network to modulate silkworm biological processes, such as growth, metamorphosis and immune responses. Furthermore, a lncRNA/circRNA ceRNA network consisting of 140 lncRNAs, five circRNA, and seven miRNAs are constructed hypothetically to describe eight key host genes regulation, such as Toll-6, Serpin-6, inducible nitric oxide synthase (iNOS) and Caspase-8. Notably, cross-species analyses indicate that parasite and host miRNAs play a vital role in pathogen-host interaction in the microsporidia congenital infection.

CONCLUSION

This is the first comprehensive pan-transcriptome study inclusive of both N. bombycis and its host silkworm with a specific focus on the microsporidia congenital infection, and show that ncRNA-mediated regulation plays a vital role in the microsporidia congenital infection, which provides a new insight into understanding the basic biology of microsporidia and pathogen-host interaction.

Comparative Analysis of CRISPR-Cas Systems in Pseudomonas Genomes

Pseudomonas is a bacterial genus with some saprophytic species from land and others associated with opportunistic infections in humans and animals. Factors such as pathogenicity or metabolic aspects have been related to CRISPR-Cas, and in silico studies into it have focused more on the clinical and non-environmental setting. This work aimed to perform an in silico analysis of the CRISPR-Cas systems present in Pseudomonas genomes. It analyzed 275 complete genomic sequences of Pseudomonas taken from the NCBI database. CRISPR loci were obtained from CRISPRdb. The genes associated with CRISPR (cas) and CAS proteins, and the origin and diversity of spacer sequences, were identified and compared by BLAST. The presence of self-targeting sequences, PAMs, and the conservation of DRs were visualized using WebLogo 3.6. The CRISPR-like RNA secondary structure prediction was analyzed using RNAFold and MFold. CRISPR structures were identified in 19.6% of Pseudomonas species. In all, 113 typical CRISPR arrays with 18 putative cas were found, as were 2050 spacers, of which 52% showed homology to bacteriophages, 26% to chromosomes, and 22% to plasmids. No potential self-targeting was detected within the CRISPR array. All the found DRs can form thermodynamically stable secondary RNA structures. The comparison of the CRISPR/Cas system can help understand the environmental adaptability of each evolutionary lineage of clinically and environmentally relevant species, providing data support for bacterial typing, traceability, analysis, and exploration of unconventional CRISPR.

Emergence of a novel sublineage, MYMBD21 under SA-2018 lineage of Foot-and-Mouth Disease Virus serotype O in Bangladesh

Foot-and-Mouth Disease (FMD) hinders the growth of the livestock industry in endemic countries like Bangladesh. The management and prevention of FMD are severely impacted by the high mutation rate and subsequent frequent generation of newer genotypes of the causative agent, Foot-and-Mouth Disease Virus (FMDV). The current study was conducted in nine districts of Bangladesh during 2019-21 to characterize the circulating FMDV strains based on the VP1 sequence analysis, the major antigenic recognition site providing serotype specificity and high variability of FMDV. This study detected the first emergence of the SA-2018 lineage in Bangladesh along with the predominance of Ind-2001e (or Ind-2001BD1) sublineage of ME-SA topotype under serotype O during 2019-21. The mutational spectrum, evolutionary divergence analysis and multidimensional plotting confirmed the isolates collected from Mymensingh districts, designated as MYMBD21 as a novel sublineage under the SA-2018 lineage. Analysis of the amino acid sequence revealed several changes in the G-H loop, B-C loop and C-terminal region of VP1, revealing a 12-13% divergence from the existing vaccine strains and a 95% VP1 protein homology, with most of the mutations potentially considerable as vaccine escape mutations, evidenced by three-dimensional structural analysis. This is the first report on the emergence of the SA-2018 lineage of ME-SA topotype of FMDV serotype O in Bangladesh, as well as a possible mutational trend towards the emergence of a distinct sublineage under SA-2018 lineage, which calls for in-depth genome-wide analysis and monitoring of the FMD situation in the country to implement a strategic vaccination and effective FMD control program.

In silico annotation of a hypothetical protein from Listeria monocytogenes EGD-e unfolds a toxin protein of the type II secretion system

The gram-positive bacterium Listeria monocytogenes is an important foodborne intracellular pathogen that is widespread in the environment. The functions of hypothetical proteins (HP) from various pathogenic bacteria have been successfully annotated using a variety of bioinformatics strategies. In this study, a HP Imo0888 (NP_464414.1) from the Listeria monocytogenes EGD-e strain was annotated using several bioinformatics tools. Various techniques, including CELLO, PSORTb, and SOSUIGramN, identified the candidate protein as cytoplasmic. Domain and motif analysis revealed that the target protein is a PemK/MazF-like toxin protein of the type II toxin-antitoxin system (TA) which was consistent with BLASTp analysis. Through secondary structure analysis, we found the random coil to be the most frequent. The Alpha Fold 2 Protein Structure Prediction Database was used to determine the three-dimensional (3D) structure of the HP using the template structure of a type II TA PemK/MazF family toxin protein (DB ID_AFDB: A0A4B9HQB9) with 99.1% sequence identity. Various quality evaluation tools, such as PROCHECK, ERRAT, Verify 3D, and QMEAN were used to validate the 3D structure. Following the YASARA energy minimization method, the target protein's 3D structure became more stable. The active site of the developed 3D structure was determined by the CASTp server. Most pathogens that harbor TA systems create a crucial risk to human health. Our aim to annotate the HP Imo088 found in Listeria could offer a chance to understand bacterial pathogenicity and identify a number of potential targets for drug development.

Mycorrhizal diversity and community composition in co-occurring Cypripedium species

Orchids commonly rely on mycorrhizal fungi to obtain the necessary resources for seed germination and growth. Whereas most photosynthetic orchids typically associate with so-called rhizoctonia fungi to complete their life cycle, there is increasing evidence that other fungi may be involved as well and that the mycorrhizal communities associated with orchids may be more diverse. Coexisting orchid species also tend to associate with different fungi to reduce competition for similar resources and to increase long-term population viability. However, few studies have related the mycorrhizal communities in the rhizosphere to communities found in the roots of closely related coexisting orchid species. In this study, we used high-throughput sequencing to investigate the diversity and community composition of orchid mycorrhizal fungi in the roots and the rhizosphere of four Cypripedium species growing in forests in Northeast China. The results showed that the investigated Cypripedium species associated with a wide variety of fungi including members of Tulasnellaceae, Psathyrellaceae, and Herpotrichiellaceae, whereas members of Russulaceae, Cortinariaceae, Thelephoraceae, and Herpotrichiellaceae showed high abundance in rhizosphere soils. The diversity of fungi detected in the rhizosphere soil was much higher than that in the roots. The observed variation in fungal communities in Cypripedium roots was not related to forest site or orchid species. On the other hand, variation in mycorrhizal communities of rhizosphere soil was significantly related to sampling site. These results indicate that orchid mycorrhizal communities in the rhizosphere display considerable variation among sites and that orchids use only a subset of the locally available fungi. Future studies focusing on the fine-scale spatial distribution of orchid mycorrhizal fungi and more detailed assessments of local environmental conditions will provide novel insights into the mechanisms explaining variation of fungal communities in both orchid roots and the rhizosphere.

C.acnes in the joint, is it all just a false positive?

BACKGROUND

Cutibacterium Acnes (C.acnes) has been linked to several shoulder pathologies. An alternative hypothesis suggests it only occurs in the joint secondary to previous instrumentation. Our hypothesis was patients with previous instrumentation would have C.acnes in their joint if it was in skin.

MATERIALS AND METHODS

Sixty-six patients undergoing arthroscopic shoulder surgery had biopsies taken from the affected joint at the time of surgery, along with control biopsies of subdermal fat. The extended culture results were assessed and correlated to previous intervention.

RESULTS

35% tested positive for C.acnes in their joint. 78% were male. 53% had absence of C.acnes in both skin and joint and 29% had presence in both (p = 0.0001). 15% with previous surgery had C.acnes. 53% with previous injection had C.acnes. 25% of patients with virgin joints had C.acnes. There was no statistical difference in the presence of C.acnes in the joint between those with previous instrumentation and without.

CONCLUSION

The significant factors for joint C.acnes were male sex and the presence of the bacteria in the fat. Previous instrumentation was not correlated with C.acnes in the joint. This raises the question of whether the process of biopsy itself may lead to inoculation of the joint.

Genome-scale signatures of adaptive gene expression changes in an invasive seaweed Gracilaria vermiculophylla

Invasive species can successfully and rapidly colonize new niches and expand ranges via founder effects and enhanced tolerance towards environmental stresses. However, the underpinning molecular mechanisms (i.e., gene expression changes) facilitating rapid adaptation to harsh environments are still poorly understood. The red seaweed Gracilaria vermiculophylla, which is native to the northwest Pacific but invaded North American and European coastal habitats over the last 100 years, provides an excellent model to examine whether enhanced tolerance at the level of gene expression contributed to its invasion success. We collected G. vermiculophylla from its native range in Japan and from two non-native regions along the Delmarva Peninsula (Eastern United States) and in Germany. Thalli were reared in a common garden for 4 months at which time we performed comparative transcriptome (mRNA) and microRNA (miRNA) sequencing. MRNA-expression profiling identified 59 genes that were differently expressed between native and non-native thalli. Of these genes, most were involved in metabolic pathways, including photosynthesis, abiotic stress, and biosynthesis of products and hormones in all four non-native sites. MiRNA-based target-gene correlation analysis in native/non-native pairs revealed that some target genes are positively or negatively regulated via epigenetic mechanisms. Importantly, these genes are mostly associated with metabolism and defence capability (e.g., metal transporter Nramp5, senescence-associated protein, cell wall-associated hydrolase, ycf68 protein and cytochrome P450-like TBP). Thus, our gene expression results indicate that resource reallocation to metabolic processes is most likely a predominant mechanism contributing to the range-wide persistence and adaptation of G. vermiculophylla in the invaded range. This study, therefore, provides molecular insight into the speed and nature of invasion-mediated rapid adaption.

MiMeDB: the Human Microbial Metabolome Database

The Human Microbial Metabolome Database (MiMeDB) (https://mimedb.org) is a comprehensive, multi-omic, microbiome resource that connects: (i) microbes to microbial genomes; (ii) microbial genomes to microbial metabolites; (iii) microbial metabolites to the human exposome and (iv) all of these 'omes' to human health. MiMeDB was established to consolidate the growing body of data connecting the human microbiome and the chemicals it produces to both health and disease. MiMeDB contains detailed taxonomic, microbiological and body-site location data on most known human microbes (bacteria and fungi). This microbial data is linked to extensive genomic and proteomic sequence data that is closely coupled to colourful interactive chromosomal maps. The database also houses detailed information about all the known metabolites generated by these microbes, their structural, chemical and spectral properties, the reactions and enzymes responsible for these metabolites and the primary exposome sources (food, drug, cosmetic, pollutant, etc.) that ultimately lead to the observed microbial metabolites in humans. Additional, extensively referenced data about the known or presumptive health effects, measured biosample concentrations and human protein targets for these compounds is provided. All of this information is housed in richly annotated, highly interactive, visually pleasing database that has been designed to be easy to search, easy to browse and easy to navigate. Currently MiMeDB contains data on 626 health effects or bioactivities, 1904 microbes, 3112 references, 22 054 reactions, 24 254 metabolites or exposure chemicals, 648 861 MS and NMR spectra, 6.4 million genes and 7.6 billion DNA bases. We believe that MiMeDB represents the kind of integrated, multi-omic or systems biology database that is needed to enable comprehensive multi-omic integration.

An integrative re-evaluation of Typhlatya shrimp within the karst aquifer of the Yucatán Peninsula, Mexico

The Yucatán Peninsula, Mexico is a carbonate platform well-known for extensive karst networks of densely stratified aquifer ecosystems. This aquifer supports diverse anchialine fauna, including species of the globally distributed anchialine shrimp genus Typhlatya (Atyidae). Four species (T. campecheae, T. pearsei, T. dzilamensis and T. mitchelli) are endemic to the Peninsula, of which three are federally listed in Mexico. This first integrative evaluation (i.e., molecular, morphological, broad geographic and type locality sampling, and environmental data) of Yucatán Typhlatya reveals considerable species identity conflict in prior phylogenetic assessments, broad species ranges, syntopy within cave systems and five genetic lineages (of which two are new to science). Despite sampling from the type locality of endangered T. campecheae, specimens (and molecular data) were indistinguishable from vulnerable T. pearsei. Ancestral/divergence reconstructions support convergent evolution of a low-salinity ancestor for a post-Paleogene arc Yucatán + Cuba Typhlatya clade within the anchialine Atyidae clade. A secondary adaptation for the coastal-restricted euryhaline (2–37 psu), Typhlatya dzilamensis (unknown conservation status) was identified, while remaining species lineages were low-salinity (< 5 psu) adapted and found within the meteoric lens of inland and coastal caves. This study demonstrates the need for integrative/interdisciplinary approaches when conducting biodiversity assessments in complex and poorly studied aquifers.

Integrated transcriptome and small RNA sequencing in revealing miRNA-mediated regulatory network of floral bud break in Prunus mume

MicroRNAs is one class of small non-coding RNAs that play important roles in plant growth and development. Though miRNAs and their target genes have been widely studied in many plant species, their functional roles in floral bud break and dormancy release in woody perennials is still unclear. In this study, we applied transcriptome and small RNA sequencing together to systematically explore the transcriptional and post-transcriptional regulation of floral bud break in P. mume. Through expression profiling, we identified a few candidate genes and miRNAs during different developmental stage transitions. In total, we characterized 1,553 DEGs associated with endodormancy release and 2,084 DEGs associated with bud flush. Additionally, we identified 48 known miRNAs and 53 novel miRNAs targeting genes enriched in biological processes such as floral organ morphogenesis and hormone signaling transudation. We further validated the regulatory relationship between differentially expressed miRNAs and their target genes combining computational prediction, degradome sequencing, and expression pattern analysis. Finally, we integrated weighted gene co-expression analysis and constructed miRNA-mRNA regulatory networks mediating floral bud flushing competency. In general, our study revealed the miRNA-mediated networks in modulating floral bud break in P. mume. The findings will contribute to the comprehensive understanding of miRNA-mediated regulatory mechanism governing floral bud break and dormancy cycling in wood perennials.

RLEP LAMP for the laboratory confirmation of leprosy: towards a point-of-care test

BACKGROUND

Nucleic acid-based amplification tests (NAAT), above all (q)PCR, have been applied for the detection of Mycobacterium leprae in leprosy cases and household contacts with subclinical infection. However, their application in the field poses a range of technical challenges. Loop-mediated isothermal amplification (LAMP), as a promising point-of-care NAAT does not require sophisticated laboratory equipment, is easy to perform, and is applicable for decentralized diagnosis at the primary health care level. Among a range of gene targets, the M. leprae specific repetitive element RLEP is regarded as highly sensitive and specific for diagnostic applications.  METHODS: Our group developed and validated a dry-reagent-based (DRB) RLEP LAMP, provided product specifications for customization of a ready-to-use kit (intended for commercial production) and compared it against the in-house prototype. The assays were optimized for application on a Genie® III portable fluorometer. For technical validation, 40 "must not detect RLEP" samples derived from RLEP qPCR negative exposed and non-exposed individuals, as well as from patients with other conditions and a set of closely related mycobacterial cultures, were tested together with 25 "must detect RLEP" samples derived from qPCR confirmed leprosy patients. For clinical validation, 150 RLEP qPCR tested samples were analyzed, consisting of the following categories: high-positive samples of multibacillary (MB) leprosy patients (> 10.000 bacilli/extract), medium-positive samples of MB leprosy patients (1.001-10.000 bacilli/extract), low-positive samples of MB leprosy patients (1-1.000 bacilli/extract), endemic controls and healthy non-exposed controls; each n = 30.  RESULTS: Technical validation: both LAMP formats had a limit of detection of 1.000 RLEP copies, i.e. 43-27 bacilli, a sensitivity of 92% (in-house protocol)/100% (ready-to-use protocol) and a specificity of 100%. Reagents were stable for at least 1 year at 22 °C. Clinical validation: Both formats showed a negativity rate of 100% and a positivity rate of 100% for high-positive samples and 93-100% for medium positive samples, together with a positive predictive value of 100% and semi-quantitative results. The positivity rate for low-positive samples was 77% (in-house protocol)/43% (ready-to-use protocol) and differed significantly between both formats.  CONCLUSIONS: The ready-to-use RLEP DRB LAMP assay constitutes an ASSURED test ready for field-based evaluation trials aiming for routine diagnosis of leprosy at the primary health care level.

Vector-Borne Blood Parasites of the Great-Tailed Grackle (Quiscalus mexicanus) in East-Central Texas, USA

Great-tailed grackles (Quiscalus mexicanus) have dramatically expanded into North America over the past century. However, little is known about the blood that parasites they support. Here, for the first time, we document an assemblage of trypanosome, haemosporida, and filarial nematodes co-circulating in invasive great-tailed grackles. Between February and July, 2015, 61 individuals were captured in an urban environment of College Station, Texas. Field microscopy and molecular diagnostics indicate that 52% (24/46) were visually infected with filarioid nematodes, 24% (11/46) with avian trypanosomes, and 73% (n = 44/60) with haemosporida parasites, such as Haemoproteus (Parahaemoproteus) and Plasmodium cathemerium. Overall, 87% of great-tailed grackles were infected with blood parasites. Although 50% of individuals hosted parasites from multiple phylum, no patterns of parasite assembly were observed. Results indicate that great-tailed grackles can support a relatively high level of blood parasitism. However, the consequences for avian health remain to be determined.

PolyDAT: A Generic Data Schema for Polymer Characterization

Polymers are stochastic materials that represent distributions of different molecules. In general, to quantify the distribution, polymer researchers rely on a series of chemical characterizations that each reveal partial information on the distribution. However, in practice, the exact set of characterizations that are carried out, as well as how the characterization data are aggregated and reported, is largely nonstandard across the polymer community. This scenario makes polymer characterization data highly disparate, thereby significantly slowing down the development of polymer informatics. In this work, a proposal on how structural characterization data can be organized is presented. To ensure that the system can apply universally across the entire polymer community, the proposed schema, PolyDAT, is designed to embody a minimal congruent set of vocabulary that is common across different domains. Unlike most chemical schemas, where only data pertinent to the species of interest are included, PolyDAT deploys a multi-species reaction network construct, in which every characterization on relevant species is collected to provide the most comprehensive profile on the polymer species of interest. Instead of maintaining a comprehensive list of available characterization techniques, PolyDAT provides a handful of generic templates, which align closely with experimental conventions and cover most types of common characterization techniques. This allows flexibility for the development and inclusion of new measurement methods. By providing a standard format to digitalize data, PolyDAT serves not only as an extension to BigSMILES that provides the necessary quantitative information but also as a standard channel for researchers to share polymer characterization data.

Data integration enables global biodiversity synthesis

The accessibility of global biodiversity information has surged in the past two decades, notably through widespread funding initiatives for museum specimen digitization and emergence of large-scale public participation in community science. Effective use of these data requires the integration of disconnected datasets, but the scientific impacts of consolidated biodiversity data networks have not yet been quantified. To determine whether data integration enables novel research, we carried out a quantitative text analysis and bibliographic synthesis of >4,000 studies published from 2003 to 2019 that use data mediated by the world's largest biodiversity data network, the Global Biodiversity Information Facility (GBIF). Data available through GBIF increased 12-fold since 2007, a trend matched by global data use with roughly two publications using GBIF-mediated data per day in 2019. Data-use patterns were diverse by authorship, geographic extent, taxonomic group, and dataset type. Despite facilitating global authorship, legacies of colonial science remain. Studies involving species distribution modeling were most prevalent (31% of literature surveyed) but recently shifted in focus from theory to application. Topic prevalence was stable across the 17-y period for some research areas (e.g., macroecology), yet other topics proportionately declined (e.g., taxonomy) or increased (e.g., species interactions, disease). Although centered on biological subfields, GBIF-enabled research extends surprisingly across all major scientific disciplines. Biodiversity data mobilization through global data aggregation has enabled basic and applied research use at temporal, spatial, and taxonomic scales otherwise not possible, launching biodiversity sciences into a new era.

Horizontal gene transfer and recombination analysis of SARS-CoV-2 genes helps discover its close relatives and shed light on its origin

BACKGROUND

The SARS-CoV-2 pandemic is one of  the greatest  global medical and social challenges that have emerged in recent history. Human coronavirus strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein.

RESULTS

We provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-CoV-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215-1425] of gene S and region [534-727] of gene N. Moreover, some statistically significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster.

CONCLUSIONS

The results of our horizontal gene transfer and recombination analysis suggest that SARS-CoV-2 could not only be a chimera virus resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of this dangerous virus.

Draft Genome Sequence of Pseudomonas sp. Strain T2.31D-1, Isolated from a Drilling Core Sample Obtained 414 Meters below Surface in the Iberian Pyrite Belt

We report the draft genome of Pseudomonas sp. strain T2.31D-1, which was isolated from a drilling core sample obtained 414 m below surface in the Iberian Pyrite Belt. The genome consists of a 4.7-Mb chromosome, with 4,428 coding sequences, one rRNA operon, and 59 tRNA genes, and a 31.8-kb plasmid.

We report the draft genome of Pseudomonas sp. strain T2.31D-1, which was isolated from a drilling core sample obtained 414 m below surface in the Iberian Pyrite Belt. The genome consists of a 4.7-Mb chromosome with 4,428 coding sequences, 1 rRNA operon, 59 tRNA genes, and a 31.8-kb plasmid.

Isolation, expression, and characterization of raw starch degrading α-amylase from a marine lake Bacillus megaterium NL3

A land-locked marine lake Kakaban with its significant ecological paramaters provides a unique habitat for bacteria with novel biotechnology potential that uses a diverse array of catalytic agents, including α-amylase. Aiming at the isolation of raw starch degrading α-amylase from marine biodiversity, a gene encoding BmaN2 from a sea anemone associated bacterium Bacillus megaterium NL3 was cloned and expressed in Escherichia coli ArcticExpress (DE3). It comprises an open reading frame of 1,563 nucleotides encoding BmaN2 of 520 amino acids and belongs to the glycoside hydrolase family 13 subfamily 36 (GH13_36). This α-amylase has a maximum activity at pH 6.0 and 60 °C with a specific activity of 28.7 U mg-1. The BmaN2 activity is enhanced strongly by Ca2+ but inhibited by EDTA. BmaN2 also exhibits high catalytic efficiency on soluble starch with k cat /K M value of 14.1 mL mg-1 s-1. Despite no additional starch-binding domain, BmaN2 is able to hydrolyze various raw starches, such as wheat, corn, cassava, potato, rice, sago, and canna, in which granular wheat is the preferred substrate for BmaN2. These characteristics indicate that BmaN2 is a promising raw starch degrading enzyme within the subfamily GH13_36.

Environmental DNA surveys detect distinct metazoan communities across abyssal plains and seamounts in the western Clarion Clipperton Zone

The deep seafloor serves as a reservoir of biodiversity in the global ocean, with >80% of invertebrates at abyssal depths still undescribed. These diverse and remote deep-sea communities are critically under-sampled and increasingly threatened by anthropogenic impacts, including future polymetallic nodule mining. Using a multigene environmental DNA (eDNA) metabarcoding approach, we characterized metazoan communities sampled from sediments, polymetallic nodules and seawater in the western Clarion Clipperton Zone (CCZ) to test the hypotheses that deep seamounts (a) are species richness hotspots in the abyss, (b) have structurally distinct communities in comparison to other deep-sea habitats, and (c) that seafloor particulate organic carbon (POC) flux and polymetallic nodule density are positively correlated with metazoan diversity. eDNA metabarcoding was effective at characterizing distinct biotas known to occur in association with different abyssal substrate types (e.g., nodule- and sediment-specific fauna), with distinct community composition and few taxa shared across substrates. Seamount faunas had higher overall taxonomic richness, and different community composition and biogeography than adjacent abyssal plains, with seamount communities displaying less connectivity between regions than comparable assemblages on the abyssal plains. Across an estimated gradient of low to moderate POC flux, we find lowest taxon richness at the lowest POC flux, as well as an effect of nodule size on community composition. Our results suggest that while abyssal seamounts are important reservoirs of metazoan diversity in the CCZ, given limited taxonomic overlap between seamount and plains fauna, conservation of seamount assemblages will be insufficient to protect biodiversity and ecosystem function in regions targeted for mining.

The Mosaic Architecture of NRPS-PKS in the Arbuscular Mycorrhizal Fungus Gigaspora margarita Shows a Domain With Bacterial Signature

As obligate biotrophic symbionts, arbuscular mycorrhizal fungi (AMF) live in association with most land plants. Among them, Gigaspora margarita has been deeply investigated because of its peculiar features, i.e., the presence of an intracellular microbiota with endobacteria and viruses. The genome sequencing of this fungus revealed the presence of some hybrid non-ribosomal peptide synthases-polyketide synthases (NRPS-PKS) that have been rarely identified in AMF. The aim of this study is to describe the architecture of these NRPS-PKS sequences and to understand whether they are present in other fungal taxa related to G. margarita. A phylogenetic analysis shows that the ketoacyl synthase (KS) domain of one G. margarita NRPS-PKS clusters with prokaryotic sequences. Since horizontal gene transfer (HGT) has often been advocated as a relevant evolutionary mechanism for the spread of secondary metabolite genes, we hypothesized that a similar event could have interested the KS domain of the PKS module. The bacterial endosymbiont of G. margarita, Candidatus Glomeribacter gigasporarum (CaGg), was the first candidate as a donor, since it possesses a large biosynthetic cluster involving an NRPS-PKS. However, bioinformatics analyses do not confirm the hypothesis of a direct HGT from the endobacterium to the fungal host: indeed, endobacterial and fungal sequences show a different evolution and potentially different donors. Lastly, by amplifying a NRPS-PKS conserved fragment and mining the sequenced AMF genomes, we demonstrate that, irrespective of the presence of CaGg, G. margarita, and some other related Gigasporaceae possess such a sequence.

A quick and sensitive diagnostic tool for detection of Maize streak virus

Maize streak virus disease (MSVD), caused by Maize streak virus (MSV; genus Mastrevirus), is one of the most severe and widespread viral diseases that adversely reduces maize yield and threatens food security in Africa. An effective control and management of MSVD requires robust and sensitive diagnostic tests capable of rapid detection of MSV. In this study, a loop-mediated isothermal amplification (LAMP) assay was designed for the specific detection of MSV. This test has shown to be highly specific and reproducible and able to detect MSV in as little as 10 fg/µl of purified genomic DNA obtained from a MSV-infected maize plant, a sensitivity 105 times higher to that obtained with polymerase chain reaction (PCR) in current general use. The high degree of sequence identity between Zambian and other African MSV isolates indicate that this LAMP assay can be used for detecting MSV in maize samples from any region in Africa. Furthermore, this assay can be adopted in minimally equipped laboratories and with potential use in plant clinic laboratories across Africa strengthening diagnostic capacity in countries dealing with MSD.

Analysis of 16S rRNA genes reveals reduced Fusobacterial community diversity when translocating from saliva to GI sites

Fusobacterium nucleatum is a Gram-negative oral commensal anaerobe which has been increasingly implicated in various gastrointestinal (GI) disorders, including inflammatory bowel disease, appendicitis, GI cancers. The oral cavity harbors a diverse group of Fusobacterium, and it is postulated that F. nucleatum in the GI tract originate from the mouth. It is not known, however, if all oral Fusobacterium translocate to the GI sites with equal efficiencies. Therefore, we amplified 16S rRNA genes of F. nucleatum and F. periodonticum, two closely related oral species from matched saliva, gastric aspirates, and colon or ileal pouch aspirates of three patients with inflammatory bowel disease (IBD) and three healthy controls, and saliva alone from seven patients with either active IBD or IBD in remission. The 16S rRNA gene amplicons were cloned, and the DNA sequences determined by Sanger sequencing. The results demonstrate that fusobacterial community composition differs more significantly between the oral and GI sites than between different individuals. The oral communities demonstrate the highest level of variation and have the richest pool of unique sequences, with certain nodes/strains enriched in the GI tract and others diminished during translocation. The gastric and colon/pouch communities exhibit reduced diversity and are more closely related, possibly due to selective pressure in the GI tract. This study elucidates selective translocation of oral fusobacteria to the GI tract. Identification of specific transmissible clones will facilitate risk assessment for developing Fusobacterium-implicated GI disorders.

Potential therapeutic target identification in the novel 2019 coronavirus: insight from homology modeling and blind docking study

Background

The 2019-nCoV which is regarded as a novel coronavirus is a positive-sense single-stranded RNA virus. It is infectious to humans and is the cause of the ongoing coronavirus outbreak which has elicited an emergency in public health and a call for immediate international concern has been linked to it. The coronavirus main proteinase which is also known as the 3C-like protease (3CLpro) is a very important protein in all coronaviruses for the role it plays in the replication of the virus and the proteolytic processing of the viral polyproteins. The resultant cytotoxic effect which is a product of consistent viral replication and proteolytic processing of polyproteins can be greatly reduced through the inhibition of the viral main proteinase activities. This makes the 3C-like protease of the coronavirus a potential and promising target for therapeutic agents against the viral infection.

Results

This study describes the detailed computational process by which the 2019-nCoV main proteinase coding sequence was mapped out from the viral full genome, translated and the resultant amino acid sequence used in modeling the protein 3D structure. Comparative physiochemical studies were carried out on the resultant target protein and its template while selected HIV protease inhibitors were docked against the protein binding sites which contained no co-crystallized ligand.

Conclusion

In line with results from this study which has shown great consistency with other scientific findings on coronaviruses, we recommend the administration of the selected HIV protease inhibitors as first-line therapeutic agents for the treatment of the current coronavirus epidemic.

Draft Genome Sequence of Shewanella sp. Strain T2.3D-1.1, Isolated from 121.8 Meters Deep in the Subsurface of the Iberian Pyrite Belt

Shewanella sp. strain T2.3D-1.1 was isolated from the deep subsurface of the Iberian Pyrite Belt. We report its draft genome sequence, consisting of 49 scaffolds, with a chromosome of ≈4.6 Mb and a 23.8-kb plasmid. The chromosome annotation identified 4,068 coding DNA sequences, 1 rRNA operon, and 108 tRNA genes.

Shewanella sp. strain T2.3D-1.1 was isolated from the deep subsurface of the Iberian Pyrite Belt. We report its draft genome sequence, consisting of 49 scaffolds, with a chromosome of ≈4.6 Mb and a 23.8-kb plasmid. The chromosome annotation identified 4,068 coding DNA sequences, 1 rRNA operon, and 108 tRNA genes.

Artificial Neural Networks Model for Predicting Type 2 Diabetes Mellitus Based on VDR Gene FokI Polymorphism, Lipid Profile and Demographic Data

Type 2 diabetes mellitus (T2DM) is a multifactorial disease associated with many genetic polymorphisms; among them is the FokI polymorphism in the vitamin D receptor (VDR) gene. In this case-control study, samples from 82 T2DM patients and 82 healthy controls were examined to investigate the association of the FokI polymorphism and lipid profile with T2DM in the Jordanian population. DNA was extracted from blood and genotyped for the FokI polymorphism by polymerase chain reaction (PCR) and DNA sequencing. Lipid profile and fasting blood sugar were also measured. There were significant differences in high-density lipoprotein (HDL) cholesterol and triglyceride levels between T2DM and control samples. Frequencies of the FokI polymorphism (CC, CT and TT) were determined in T2DM and control samples and were not significantly different. Furthermore, there was no significant association between the FokI polymorphism and T2DM or lipid profile. A feed-forward neural network (FNN) was used as a computational platform to predict the persons with diabetes based on the FokI polymorphism, lipid profile, gender and age. The accuracy of prediction reached 88% when all parameters were included, 81% when the FokI polymorphism was excluded, and 72% when lipids were only included. This is the first study investigating the association of the VDR gene FokI polymorphism with T2DM in the Jordanian population, and it showed negative association. Diabetes was predicted with high accuracy based on medical data using an FNN. This highlights the great value of incorporating neural network tools into large medical databases and the ability to predict patient susceptibility to diabetes.

Egg parasitoids of Arboridia apicalis (Nawa, 1913) (Hemiptera, Cicadellidae), a leafhopper pest of grapevines in Japan, with description of a new species of Anagrus Haliday, 1833 (Hymenoptera, Mymaridae)

Several species of egg parasitoids (Hymenoptera: Mymaridae and Trichogrammatidae) of the leafhopper pest of grapevines in Japan, Arboridia (Arboridia) apicalis (Nawa) (Hemiptera, Cicadellidae), were reared and identified for the first time. Using a combination of genetic and morphological evidence, Anagrus (Anagrus) arboridiae Triapitsyn & Adachi-Hagimori, sp. nov. (Mymaridae) is described and illustrated from Honshu Island (Shimane Prefecture) and Kyushu Island (Miyazaki Prefecture). It is shown to be different from Anagrus (Anagrus) japonicus Sahad and A.flaviapex Chiappini & Lin, to which it is most similar; the latter species was originally described from China and is newly recorded here from Okinawa Island, Japan. Mitochondrial and nuclear ribosomal DNA sequence data provide clear evidence for the separation of A.arboridiae from A.flaviapex, A.japonicus, and some other members of the Anagrus (Anagrus) atomus (L.) species group. Two other species of Anagrus Haliday, A. (Anagrus) avalae Soyka and A.atomus, are also identified in Japan from eggs of the leafhoppers Edwardsianaishidae (Matsumura) and Eurhadina?betularia Anufriev, respectively. An updated key to females of the Japanese species of Anagrus is given. Oligositapallida Kryger (a new record for Japan), Oligosita sp., and an Aphelinoidea (Aphelinoidea) sp. (Trichogrammatidae) were the other, although much less abundant, apparent egg parasitoids of A.apicalis in Shimane Prefecture, mainly in non-organic vineyards.

A First NGS Investigation Suggests No Association Between Viruses and Canine Cancers

Approximately 10–15% of worldwide human cancers are attributable to viral infection. When operating as carcinogenic elements, viruses may act with various mechanisms, but the most important is represented by viral integration into the host genome, causing chromosome instability, genomic mutations, and aberrations. In canine species, few reports have described an association between viral integration and canine cancers, but more comprehensive studies are needed. The advancement of next-generation sequencing and the cost reduction have resulted in a progressive increasing of sequencing data in veterinary oncology offering an opportunity to study virome in canine cancers. In this study, we have performed viral detection and integration analyses using VirusFinder2 software tool on available whole-genome and whole-exome sequencing data of different canine cancers. Several viral sequences were detected in lymphomas, hemangiosarcomas, melanomas, and osteosarcomas, but no reliable integration sites were identified. Even if with some limitations such as the depth and type of sequencing, a restricted number of available nonhuman genomes software, and a limited knowledge on endogenous retroviruses in the canine genome, results are compelling. However, further experiments are needed, and similarly to feline species, dedicated analysis tools for the identification of viral integration sites in canine cancers are required.

In silico and empirical evaluation of twelve metabarcoding primer sets for insectivorous diet analyses

During the most recent decade, environmental DNA metabarcoding approaches have been both developed and improved to minimize the biological and technical biases in these protocols. However, challenges remain, notably those relating to primer design. In the current study, we comprehensively assessed the performance of ten COI and two 16S primer pairs for eDNA metabarcoding, including novel and previously published primers. We used a combined approach of in silico, in vivo-mock community (33 arthropod taxa from 16 orders), and guano-based analyses to identify primer sets that would maximize arthropod detection and taxonomic identification, successfully identify the predator (bat) species, and minimize the time and financial costs of the experiment. We focused on two insectivorous bat species that live together in mixed colonies: the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy's bat (Myotis emarginatus). We found that primer degeneracy is the main factor that influences arthropod detection in silico and mock community analyses, while amplicon length is critical for the detection of arthropods from degraded DNA samples. Our guano-based results highlight the importance of detecting and identifying both predator and prey, as guano samples can be contaminated by other insectivorous species. Moreover, we demonstrate that amplifying bat DNA does not reduce the primers' capacity to detect arthropods. We therefore recommend the simultaneous identification of predator and prey. Finally, our results suggest that up to one-third of prey occurrences may be unreliable and are probably not of primary interest in diet studies, which may decrease the relevance of combining several primer sets instead of using a single efficient one. In conclusion, this study provides a pragmatic framework for eDNA primer selection with respect to scientific and methodological constraints.

Identification of Known and Novel MicroRNAs in Raspberry Organs Through High-Throughput Sequencing

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in plants by negatively affecting gene expression. Studies on the identification of miRNAs and their functions in various plant species and organs have significantly contributed to plant development research. In the current study, we utilized high-throughput sequencing to detect the miRNAs in the root, stem, and leaf tissues of raspberry (Rubus idaeus). A total of more than 35 million small RNA reads ranging in size from 18 to 35 nucleotides were obtained, with 147 known miRNAs and 542 novel miRNAs identified among the three organs. Sequence verification and the relative expression profiles of the six known miRNAs were investigated by stem-loop quantitative real-time PCR. Furthermore, the potential target genes of the known and novel miRNAs were predicted and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway annotation. Enrichment analysis of the GO-associated biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Moreover, the miRNA target prediction revealed that most of the targets predicted as transcription factor-coding genes are involved in cellular and metabolic processes. This report is the first to identify miRNAs in raspberry. The detected miRNAs were analyzed by cluster analysis according to their expression, which revealed that these conservative miRNAs are necessary for plant functioning. The results add novel miRNAs to the raspberry transcriptome, providing a useful resource for the further elucidation of the functional roles of miRNAs in raspberry growth and development.

Environmental DNA surveys detect distinct metazoan communities across abyssal plains and seamounts in the western Clarion Clipperton Zone

The deep seafloor serves as a reservoir of biodiversity in the global ocean, with >80% of invertebrates at abyssal depths still undescribed. These diverse and remote deep-sea communities are critically under-sampled and increasingly threatened by anthropogenic impacts, including future polymetallic nodule mining. Using a multigene environmental DNA (eDNA) metabarcoding approach, we characterized metazoan communities sampled from sediments, polymetallic nodules and seawater in the western Clarion Clipperton Zone (CCZ) to test the hypotheses that deep seamounts (a) are species richness hotspots in the abyss, (b) have structurally distinct communities in comparison to other deep-sea habitats, and (c) that seafloor particulate organic carbon (POC) flux and polymetallic nodule density are positively correlated with metazoan diversity. eDNA metabarcoding was effective at characterizing distinct biotas known to occur in association with different abyssal substrate types (e.g., nodule- and sediment-specific fauna), with distinct community composition and few taxa shared across substrates. Seamount faunas had higher overall taxonomic richness, and different community composition and biogeography than adjacent abyssal plains, with seamount communities displaying less connectivity between regions than comparable assemblages on the abyssal plains. Across an estimated gradient of low to moderate POC flux, we find lowest taxon richness at the lowest POC flux, as well as an effect of nodule size on community composition. Our results suggest that while abyssal seamounts are important reservoirs of metazoan diversity in the CCZ, given limited taxonomic overlap between seamount and plains fauna, conservation of seamount assemblages will be insufficient to protect biodiversity and ecosystem function in regions targeted for mining.

Species delimitation analyses of NE Atlantic Chaetozone (Annelida, Cirratulidae) reveals hidden diversity among a common and abundant marine annelid

The polychaetes of the family Cirratulidae (Annelida) are common inhabitants in continental shelf benthic environments and considered an important group of organisms in environmental monitoring surveys. The family represents a taxonomic and systematic challenge, as monophyly of genera and evolutionary relationships within the family remain to be explored in a proper phylogenetic framework. Bitentaculate cirratulids, especially the genus Chaetozone, form one of the most species-diverse group of polychaetes worldwide. In this study, we aimed at evaluating the species diversity of the genus Chaetozonein benthic environments in the North East Atlantic by molecular means. We tested whether traditional morphological diagnostic characters are able to discriminate between the species hypothesis after species delimitation analyses, and assessed monophyly of the genera involved. Two DNA markers were sequenced from about 200 specimens belonging to Chaetozone, Aphelochaeta, Dodecaceria, Cirriformia and Cirratulus - the universal mitochondrial barcoding region COI, and the D1-D2 regions of the nuclear 28S rRNA - and analyzed with Bayesian inference, Maximum Likelihood and the species delimitation methods mPTP and GMYC. The first phylogeny of the family Cirratulidae is inferred and the genera Chaetozone, Dodecaceria and Cirratulus are recovered monophyletic. A total of 14 clusters of sequences - corresponding to species of Chaetozone - were found in the study area, and only one of them is here referred to a nominal species, Chaetozone setosa. Our results reveal several species complexes in the genus Chaetozone, that some of these independent lineages are unnamed and undescribed, and that morphological diagnostic features are in most cases unable to discriminate between the most similar species.

Natural Communities of Carotenogenic Chlorophyte Haematococcus lacustris and Bacteria from the White Sea Coastal Rock Ponds

Haematococcus lacustris is a biotechnologically important green unicellular alga producing widely used keta-karotenoid astaxanthin. In natural habitats, it exists in the form of algal-bacterial community, and under laboratory conditions, it is also accompanied by bacteria. The issue of the bacterial composition of industrial algal cultures is widely recognized as important. However, there is a dearth of information about bacterial composition of H. lacustris communities. In current work, we analyze the composition of natural H. lacustris communities from the White Sea coastal temporal rock ponds. For the first time, a 16S rRNA gene-based metagenome of natural H. lacustris bacterial communities has been generated. Main results of its analysis are as follow. Bacterial families Comamonadaceae, Cytophagaceae, Xanthomonadaceae, Acetobacteraceae, Rhodobacteraceae, and Rhodocyclaceae were observed in all studied H. lacustris natural communities. They also contained genera Hydrogenophaga and Cytophaga. Bacteria from the Hydrogenophaga genus were present in H. lacustris cultures after their isolation under the conditions of laboratory cultivation. Similar to other planktonic microalgae, H. lacustris forms a phycosphere around the cells. In this zone, bacteria attached to the algal surface. The contact between H. lacustris and bacteria is maintained even after sample drying. The study provides information about possible members of H. lacustris core microbiome, which can be presented in the industrial and laboratory cultures of the microalga.

ppiGReMLIN: a graph mining based detection of conserved structural arrangements in protein-protein interfaces

Background

Protein-protein interactions (PPIs) are fundamental in many biological processes and understanding these interactions is key for a myriad of applications including drug development, peptide design and identification of drug targets. The biological data deluge demands efficient and scalable methods to characterize and understand protein-protein interfaces. In this paper, we present ppiGReMLIN, a graph based strategy to infer interaction patterns in a set of protein-protein complexes. Our method combines an unsupervised learning strategy with frequent subgraph mining in order to detect conserved structural arrangements (patterns) based on the physicochemical properties of atoms on protein interfaces. To assess the ability of ppiGReMLIN to point out relevant conserved substructures on protein-protein interfaces, we compared our results to experimentally determined patterns that are key for protein-protein interactions in 2 datasets of complexes, Serine-protease and BCL-2.

Results

ppiGReMLIN was able to detect, in an automatic fashion, conserved structural arrangements that represent highly conserved interactions at the specificity binding pocket of trypsin and trypsin-like proteins from Serine-protease dataset. Also, for the BCL-2 dataset, our method pointed out conserved arrangements that include critical residue interactions within the conserved motif LXXXXD, pivotal to the binding specificity of BH3 domains of pro-apoptotic BCL-2 proteins towards apoptotic suppressors. Quantitatively, ppiGReMLIN was able to find all of the most relevant residues described in literature for our datasets, showing precision of at least 69% up to 100% and recall of 100%.

Conclusions

ppiGReMLIN was able to find highly conserved structures on the interfaces of protein-protein complexes, with minimum support value of 60%, in datasets of similar proteins. We showed that the patterns automatically detected on protein interfaces by our method are in agreement with interaction patterns described in the literature.

Macrostylis metallicola spec. nov.-an isopod with geographically clustered genetic variability from a polymetallic-nodule area in the Clarion-Clipperton Fracture Zone

Background

The Clarion-Clipperton Fracture Zone (CCFZ) in the Northeast Central Pacific Ocean is a region of heightened scientific and public interest because of its wealth in manganese nodules. Due to a poor ecological understanding at the abyssal seafloor and limited knowledge of the organisms inhabiting this area, huge efforts in alpha taxonomy are required. To predict and manage potential hazards associated with future mining, taxonomy is an essential first step to grasp fundamental ecosystem traits, such as biogeographic patterns, connectivity, and the potential for post-impact recolonization. Amongst samples from the Global Sea Mineral Resources NV exploration area (EA) in the CCFZ an undescribed species of the isopod crustacean family Macrostylidae was discovered. Previously, it has been reported from two other nearby regions, the Institut Français de Recherche pour l’Exploitation de la Mer and BGR EAs. There it was one of the more widely distributed and abundant species of the benthic macrofauna and exhibited geographically structured populations. It nevertheless remained taxonomically undescribed so far.

Methods

The new species is described by means of integrative taxonomy. Morphologically, macro photography, confocal microscopy, scanning electron microscopy and light microscopy were used to describe the species and to get first insights on its phylogenetic origin. Additionally, mitochondrial DNA markers were used to test the morphological allocation of the two dimorphic sexes and juvenile stages, to analyze geographic patterns of genetic differentiation, and to study intra-and inter-species relationships, also in light of previously published population genetics on this species.

Results

The new species, Macrostylis metallicola spec. nov., is a typical representative of Macrostylidae as recognizable from the fossosoma, prognathous cephalothorax, and styliform uropods. It can be morphologically distinguished from congeners by a combination of character states which include the autapomorphic shape of the first pleopod of the copulatory male. A sexual dimorphism, as expressed by a peculiar sequence of article length-width ratios of the male antennula, indicates a relationship with M. marionaeKniesz, Brandt & Riehl (2018) and M. longipesHansen (1916) amongst other species sharing this dimorphism. Mitochondrial genetic markers point in a similar direction. M. metallicola appears to be amongst the more common and widely distributed components of the benthic macrofauna in this region which may suggest a resilience of this species to future mining activities because of its apparent potential for recolonization of impacted sites from adjacent areas of particular environmental interest. The genetic data, however, show geographic clustering of its genetic variability, pointing towards a limited potential for dispersal. Local extinction of populations could potentially not be compensated quickly and would mean a loss of genetic diversity of this species.

Characterization of Dermotheca sp. Infection in a midwestern state-endangered salamander (Ambystoma platineum) and a co-occurring common species (Ambystoma texanum)

Ichthyosporean parasites (order Dermocystida) can cause morbidity and mortality in amphibians, but their ecology and epidemiology remain understudied. We investigated the prevalence, gross and histologic appearance, and molecular phylogeny of a novel dermocystid in the state-endangered silvery salamander (Ambystoma platineum) and the co-occurring, non-threatened small-mouthed salamander (Ambystoma texanum) from Illinois. Silvery salamanders (N = 610) were sampled at six ephemeral wetlands from 2016 to 2018. Beginning in 2017, 1-3 mm raised, white skin nodules were identified in 24 silvery salamanders and two small-mouthed salamanders from five wetlands (prevalence = 0-11.1%). Skin biopsy histology (N = 4) was consistent with dermocystid sporangia, and necropsies (N = 3) identified infrequent hepatic sporangia. Parasitic 18S rRNA sequences (N = 5) from both salamander species were identical, and phylogenetic analysis revealed a close relationship to Dermotheca viridescens. Dermocystids were not identified in museum specimens from the same wetlands (N = 125) dating back to 1973. This is the first report of Dermotheca sp. affecting caudates in the Midwestern United States. Future research is needed to determine the effects of this pathogen on individual and population health, and to assess whether this organism poses a threat to the conservation of ambystomatid salamanders.

IntaRNAhelix-composing RNA–RNA interactions from stable inter-molecular helices boosts bacterial sRNA target prediction

Efficient computational tools for the identification of putative target RNAs regulated by prokaryotic sRNAs rely on thermodynamic models of RNA secondary structures. While they typically predict RNA–RNA interaction complexes accurately, they yield many highly-ranked false positives in target screens. One obvious source of this low specificity appears to be the disability of current secondary-structure-based models to reflect steric constraints, which nevertheless govern the kinetic formation of RNA–RNA interactions. For example, often — even thermodynamically favorable — extensions of short initial kissing hairpin interactions are kinetically prohibited, since this would require unwinding of intra-molecular helices as well as sterically impossible bending of the interaction helix. Another source is the consideration of instable and thus unlikely subinteractions that enable better scoring of longer interactions. In consequence, the efficient prediction methods that do not consider such effects show a high false positive rate.

To increase the prediction accuracy we devise IntaRNAhelix, a dynamic programming algorithm that length-restricts the runs of consecutive inter-molecular base pairs (perfect canonical stackings), which we hypothesize to implicitly model the steric and kinetic effects. The novel method is implemented by extending the state-of-the-art tool IntaRNA.

Our comprehensive bacterial sRNA target prediction benchmark demonstrates significant improvements of the prediction accuracy and enables more than 40-times faster computations. These results indicate — supporting our hypothesis — that stable helix composition increases the accuracy of interaction prediction models compared to the current state-of-the-art approach.

Molecular Biomarkers and Influential Factors of Denitrification in a Full-Scale Biological Nitrogen Removal Plant

Three denitrifying bacteria, Paracoccus spp., Thauera spp., Pseudomonas-like spp., and two functional genes, nitrate reductase (narG and napA), were studied as potential biomarkers for total nitrogen removal. These bacterial genera and the functional genes showed significant negative correlations with total nitrogen in the effluent (TN_eff). Thauera spp. had the highest correlation (r = −0.793, p < 0.001) with TN_eff, and narG-like and napA genes also showed significant correlations (r = −0.663 and −0.643, respectively), suggesting functional genes have equal validity to 16S rRNA genes in monitoring denitrification performance. The most explanatory variables were a combination of constituents, with temperature emerging as the most important in Pearson’s correlation and redundancy analysis. Thauera spp. had the highest correlation with temperature (r = 0.739) followed closely by Paracoccus spp. (r = 0.705). Denitrification was also significantly affected by pH (r = 0.369), solids retention time (r = −0.377), total nitrogen_in (r = 0.635), and organic matter in the influent (biochemical oxygen demand and chemical oxygen demand; r = 0.320 and 0.522, respectively). Our data verified that major denitrifiers’ 16S rRNA genes and nitrate reductase genes were better biomarkers than the biomass concentration, and any of the biomarkers could track denitrification in real time.

Exploring snake occurrence records: Spatial biases and marginal gains from accessible social media

A species' distribution provides fundamental information on: climatic niche, biogeography, and conservation status. Species distribution models often use occurrence records from biodiversity databases, subject to spatial and taxonomic biases. Deficiencies in occurrence data can lead to incomplete species distribution estimates. We can incorporate other data sources to supplement occurrence datasets. The general public is creating (via GPS-enabled cameras to photograph wildlife) incidental occurrence records that may present an opportunity to improve species distribution models. We investigated (1) occurrence data of a cryptic group of animals: non-marine snakes, in a biodiversity database (Global Biodiversity Information Facility (GBIF)) and determined (2) whether incidental occurrence records extracted from geo-tagged social media images (Flickr) could improve distribution models for 18 tropical snake species. We provide R code to search for and extract data from images using Flickr's API. We show the biodiversity database's 302,386 records disproportionately originate from North America, Europe and Oceania (250,063, 82.7%), with substantial gaps in tropical areas that host the highest snake diversity. North America, Europe and Oceania averaged several hundred records per species; whereas Asia, Africa and South America averaged less than 35 per species. Occurrence density showed similar patterns; Asia, Africa and South America have roughly ten-fold fewer records per 100 km2than other regions. Social media provided 44,687 potential records. However, including them in distribution models only marginally impacted niche estimations; niche overlap indices were consistently over 0.9. Similarly, we show negligible differences in Maxent model performance between models trained using GBIF-only and Flickr-supplemented datasets. Model performance appeared dependent on species, rather than number of occurrences or training dataset. We suggest that for tropical snakes, accessible social media currently fails to deliver appreciable benefits for estimating species distributions; but due to the variation between species and the rapid growth in social media data, may still be worth considering in future contexts.

A reference library for Canadian invertebrates with 1.5 million barcodes, voucher specimens, and DNA samples

The reliable taxonomic identification of organisms through DNA sequence data requires a well parameterized library of curated reference sequences. However, it is estimated that just 15% of described animal species are represented in public sequence repositories. To begin to address this deficiency, we provide DNA barcodes for 1,500,003 animal specimens collected from 23 terrestrial and aquatic ecozones at sites across Canada, a nation that comprises 7% of the planet's land surface. In total, 14 phyla, 43 classes, 163 orders, 1123 families, 6186 genera, and 64,264 Barcode Index Numbers (BINs; a proxy for species) are represented. Species-level taxonomy was available for 38% of the specimens, but higher proportions were assigned to a genus (69.5%) and a family (99.9%). Voucher specimens and DNA extracts are archived at the Centre for Biodiversity Genomics where they are available for further research. The corresponding sequence and taxonomic data can be accessed through the Barcode of Life Data System, GenBank, the Global Biodiversity Information Facility, and the Global Genome Biodiversity Network Data Portal.