Parallelization of MAFFT for large-scale multiple sequence alignments

Divergence time estimation of genus Tribolium by extensive sampling of highly conserved orthologs

Tribolium castaneum, the red flour beetle, is among the most well-studied eukaryotic genetic model organisms. Tribolium often serves as a comparative bridge from highly derived Drosophila traits to other organisms. Simultaneously, as a member of the most diverse order of metazoans, Coleoptera, Tribolium informs us about innovations that accompany hyper diversity. However, understanding the tempo and mode of evolutionary innovation requires well-resolved, time-calibrated phylogenies, which are not available for Tribolium. The most recent effort to understand Tribolium phylogenetics used two mitochondrial and three nuclear markers. The study concluded that the genus may be paraphyletic and reported a broad range for divergence time estimates. Here we employ recent advances in Bayesian methods to estimate the relationships and divergence times among Tribolium castaneum, T. brevicornis, T. confusum, T. freemani, and Gnatocerus cornutus using 1368 orthologs conserved across all five species and an independent substitution rate estimate. We find that the most basal split within Tribolium occurred ~86 Mya [95% HPD 85.90-87.04 Mya] and that the most recent split was between T. freemani and T. castaneum at ~14 Mya [95% HPD 13.55-14.00]. Our results are consistent with broader phylogenetic analyses of insects and suggest that Cenozoic climate changes played a role in the Tribolium diversification.

Pedobacter gandavensis sp. nov., Pedobacter foliorum sp. nov. and Pedobacter planticolens sp. nov., isolated from leaves of Arabidopsis thaliana

Three rod-shaped, non-spore-forming, yellow or pale-yellow pigmented bacteria with distinct MALDI-TOF mass spectra were isolated from the phyllosphere of Arabidopsis thaliana seedlings. Their 16S rRNA gene sequences demonstrated that these isolates belong to the genus Pedobacter. The nearest phylogenetic neighbours of strain LMG 31462^T were Pedobacter steynii DSM 19110^T (98.3 % 16S rRNA sequence similarity) and Pedobacter caeni LMG 22862^T (98.3 %); the nearest phylogenetic neighbours of strain LMG 31463^T were Pedobacter panaciterrae Gsoil 042^T (98.3 %) and Pedobacter nutrimenti DSM 27372^T (98.1 %); and the nearest phylogenetic neighbours of strain LMG 31464^T were Pedobacter boryungensis BR-9^T (99.0 %) and Pedobacter daejeonensis THG-DN3.18^T (98.7 %). Average nucleotide identity analyses between the whole genome sequences of the three strains and of the type strains of their respective nearest-neighbour taxa yielded values well below the species delineation threshold and thus confirmed that the three strains represented a novel Pedobacter species each. An extensive phenotypic comparison and an analysis of whole-cell fatty acid components yielded distinctive phenotypic characteristics for each of these strains. We therefore propose to classify these isolates as three novel species, for which we propose the names Pedobacter gandavensis with LMG 31462^T (=R-74704^T=CECT 30149^T) as the type strain, Pedobacter foliorum with LMG 31463^T (=R-74623^T=CECT 30150^T) as the type strain and Pedobacter planticolens with LMG 31464^T (=R-74626^T=CECT 30151^T) as the type strain.

'Candidatus Sarmatiella mevalonica' endosymbiont of the ciliate Paramecium provides insights on evolutionary plasticity among Rickettsiales

Members of the bacterial order Rickettsiales are obligatorily associated with a wide range of eukaryotic hosts. Their evolutionary trajectories, in particular concerning the origin of shared or differential traits among distant sub-lineages, are still poorly understood. Here, we characterized a novel Rickettsiales bacterium associated with the ciliate Paramecium tredecaurelia and phylogenetically related to the Rickettsia genus. Its genome encodes significant lineage-specific features, chiefly the mevalonate pathway gene repertoire, involved in isoprenoid precursor biosynthesis. Not only this pathway has never been described in Rickettsiales, it also is very rare among bacteria, though typical in eukaryotes, thus likely representing a horizontally acquired trait. The presence of these genes could enable an efficient exploitation of host-derived intermediates for isoprenoid synthesis. Moreover, we hypothesize the reversed reactions could have replaced canonical pathways for producing acetyl-CoA, essential for phospholipid biosynthesis. Additionally, we detected phylogenetically unrelated mevalonate pathway genes in metagenome-derived Rickettsiales sequences, likely indicating evolutionary convergent effects of independent horizontal gene transfer events. Accordingly, convergence, involving both gene acquisitions and losses, is highlighted as a relevant evolutionary phenomenon in Rickettsiales, possibly favoured by plasticity and comparable lifestyles, representing a potentially hidden origin of other more nuanced similarities among sub-lineages.

ResiDB: An automated database manager for sequence data

The amount of publicly available DNA sequence data is drastically increasing, making it a tedious task to create sequence databases necessary for the design of diagnostic assays. The selection of appropriate sequences is especially challenging in genes affected by frequent point mutations such as antibiotic resistance genes. To overcome this issue, we have designed the webtool resiDB, a rapid and user-friendly sequence database manager for bacteria, fungi, viruses, protozoa, invertebrates, plants, archaea, environmental and whole genome shotgun sequence data. It automatically identifies and curates sequence clusters to create custom sequence databases based on user-defined input sequences. A collection of helpful visualization tools gives the user the opportunity to easily access, evaluate, edit, and download the newly created database. Consequently, researchers do no longer have to manually manage sequence data retrieval, deal with hardware limitations, and run multiple independent software tools, each having its own requirements, input and output formats. Our tool was developed within the H2020 project FAPIC aiming to develop a single diagnostic assay targeting all sepsis-relevant pathogens and antibiotic resistance mechanisms. ResiDB is freely accessible to all users through https://residb.ait.ac.at/.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

Integrative Analysis of Whole Genome Sequencing and Phenotypic Resistance Toward Prediction of Trimethoprim-Sulfamethoxazole Resistance in Staphylococcus aureus

As whole genome sequencing is becoming more accessible and affordable for clinical microbiological diagnostics, the reliability of genotypic antimicrobial resistance (AMR) prediction from sequencing data is an important issue to address. Computational AMR prediction can be performed at multiple levels. The first-level approach, such as simple AMR search relies heavily on the quality of the information fed into the database. However, AMR due to mutations are often undetected, since this is not included in the database or poorly documented. Using co-trimoxazole (trimethoprim-sulfamethoxazole) resistance in Staphylococcus aureus, we compared single-level and multi-level analysis to investigate the strengths and weaknesses of both approaches. The results revealed that a single mutation in the AMR gene on the nucleotide level may produce false positive results, which could have been detected if protein sequence analysis would have been performed. For AMR predictions based on chromosomal mutations, such as the folP gene of S. aureus, natural genetic variations should be taken into account to differentiate between variants linked to genetic lineage (MLST) and not over-estimate the potential resistant variants. Our study showed that careful analysis of the whole genome data and additional criterion such as lineage-independent mutations may be useful for identification of mutations leading to phenotypic resistance. Furthermore, the creation of reliable database for point mutations is needed to fully automatized AMR prediction.

Earliest detection to date of SARS-CoV-2 in Florida: Identification together with influenza virus on the main entry door of a university building, February 2020

In February and March, 2020, environmental surface swab samples were collected from the handle of the main entry door of a major university building in Florida, as part of a pilot surveillance project screening for influenza. Samples were taken at the end of regular classroom hours, between the dates of February 1–5 and February 19-March 4, 2020. Influenza A(H1N1)pdm09 virus was isolated from the door handle on four of the 19 days sampled. Both SARS-CoV-2 and A(H1N1)pdm09 virus were detected in a sample collected on February 21, 2020. Based on sequence analysis, the Florida SARS-CoV-2 strain (designated UF-11) was identical to strains being identified in Washington state during the same time period, while the earliest similar sequences were sampled in China/Hubei between Dec 30^th 2019 and Jan 5^th 2020. The first human case of COVID-19 was not officially reported in Florida until March 1^st. In an analysis of sequences from COVID-19 patients in this region of Florida, there was only limited evidence of subsequent dissemination of the UF-11 strain. Identical or highly similar strains, possibly related through a common transmission chain, were detected with increasing frequency in Washington state between end of February and beginning of March. Our data provide further documentation of the rapid early spread of SARS-CoV-2 and underscore the likelihood that closely related strains were cryptically circulating in multiple U.S. communities before the first “official” cases were recognized.

The complete chloroplast genome sequence of the Sechium edule (Jacq.) Swartz. (Cucurbitaceae)

Sechium edule (Jacq.) Swartz is an important vegetable with both food and medicinal values. The complete chloroplast genome sequence of S. edule has been reported in this study. The total genome size is 154,558 bp in length and contains a pair of inverted repeats (IRs) of 19,128 bp, which were separated by large single-copy (LSC) and small single-copy (SSC) of 98,806 and 17,496 bp, respectively. A total of 122 genes were predicted including 78 protein-coding genes, 8 rRNA genes, and 36 tRNA genes. Further, the phylogenetic analysis confirmed that S. edule belongs to the family Cucurbitaceae. The complete chloroplast genome of S. edule would play a significant role in the development of molecular markers for plant phylogenetic and population genetic studies.

Full length genomic sanger sequencing and phylogenetic analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in Nigeria

In an outbreak, effective detection of the aetiological agent(s) involved using molecular techniques is key to efficient diagnosis, early prevention and management of the spread. However, sequencing is necessary for mutation monitoring and tracking of clusters of transmission, development of diagnostics and for vaccines and drug development. Many sequencing methods are fast evolving to reduce test turn-around-time and to increase through-put compared to Sanger sequencing method; however, Sanger sequencing remains the gold standard for clinical research sequencing with its 99.99% accuracy This study sought to generate sequence data of SARS-CoV-2 using Sanger sequencing method and to characterize them for possible site(s) of mutations. About 30 pairs of primers were designed, synthesized, and optimized using endpoint PCR to generate amplicons for the full length of the virus. Cycle sequencing using BigDye Terminator v.3.1 and capillary gel electrophoresis on ABI 3130xl genetic analyser were performed according to the manufacturers’ instructions. The sequence data generated were assembled and analysed for variations using DNASTAR Lasergene 17 SeqMan Ultra. Total length of 29,760bp of SARS-CoV-2 was assembled from the sample analysed and deposited in GenBank with accession number: MT576584. Blast result of the sequence assembly shows a 99.97% identity with the reference sequence. Variations were noticed at positions: nt201, nt2997, nt14368, nt16535, nt20334, and nt28841-28843, which caused amino acid alterations at the S (aa614) and N (aa203-204) regions. The mutations observed at S and N-gene in this study may be indicative of a gradual changes in the genetic coding of the virus hence, the need for active surveillance of the viral genome.

Long rDNA amplicon sequencing of insect-infecting nephridiophagids reveals their affiliation to the Chytridiomycota and a potential to switch between hosts

Nephridiophagids are unicellular eukaryotes that parasitize the Malpighian tubules of numerous insects. Their life cycle comprises multinucleate vegetative plasmodia that divide into oligonucleate and uninucleate cells, and sporogonial plasmodia that form uninucleate spores. Nephridiophagids are poor in morphological characteristics, and although they have been tentatively identified as early-branching fungi based on the SSU rRNA gene sequences of three species, their exact position within the fungal tree of live remained unclear. In this study, we describe two new species of nephridiophagids (Nephridiophaga postici and Nephridiophaga javanicae) from cockroaches. Using long-read sequencing of the nearly complete rDNA operon of numerous further species obtained from cockroaches and earwigs to improve the resolution of the phylogenetic analysis, we found a robust affiliation of nephridiophagids with the Chytridiomycota-a group of zoosporic fungi that comprises parasites of diverse host taxa, such as microphytes, plants, and amphibians. The presence of the same nephridiophagid species in two only distantly related cockroaches indicates that their host specificity is not as strict as generally assumed.

Pseudooceanicola algae sp. nov., isolated from the marine macroalga Fucus spiralis, shows genomic and physiological adaptations for an algae-associated lifestyle

The genus Pseudooceanicola from the alphaproteobacterial Roseobacter group currently includes ten validated species. We herein describe strain Lw-13eT, the first Pseudooceanicola species from marine macroalgae, isolated from the brown alga Fucus spiralis abundant at European and North American coasts. Physiological and pangenome analyses of Lw-13eT showed corresponding adaptive features. Adaptations to the tidal environment include a broad salinity tolerance, degradation of macroalgae-derived substrates (mannitol, mannose, proline), and resistance to several antibiotics and heavy metals. Notably, Lw-13eT can degrade oligomeric alginate via PL15 alginate lyase encoded in a polysaccharide utilization locus (PUL), rarely described for roseobacters to date. Plasmid localization of the PUL strengthens the importance of mobile genetic elements for evolutionary adaptations within the Roseobacter group. PL15 homologs were primarily detected in marine plant-associated metagenomes from coastal environments but not in the open ocean, corroborating its adaptive role in algae-rich habitats. Exceptional is the tolerance of Lw-13eT against the broad-spectrum antibiotic tropodithietic acid, produced by Phaeobacter spp. co-occurring in coastal habitats. Furthermore, Lw-13eT exhibits features resembling terrestrial plant-bacteria associations, i.e. biosynthesis of siderophores, terpenes and volatiles, which may contribute to mutual bacteria-algae interactions. Closest described relative of Lw-13eT is Pseudopuniceibacterium sediminis CY03T with 98.4% 16S rRNA gene sequence similarity. However, protein sequence-based core genome phylogeny and average nucleotide identity indicate affiliation of Lw-13eT with the genus Pseudooceanicola. Based on phylogenetic, physiological and (chemo)taxonomic distinctions, we propose strain Lw-13eT (=DSM 29013T=LMG 30557T) as a novel species with the name Pseudooceanicola algae.

Transmission dynamics of SARS-CoV-2 within-host diversity in two major hospital outbreaks in South Africa

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes acute, highly transmissible respiratory infection in humans and a wide range of animal species. Its rapid global spread has resulted in a major public health emergency, necessitating commensurately rapid research to improve control strategies. In particular, the ability to effectively retrace transmission chains in outbreaks remains a major challenge, partly due to our limited understanding of the virus' underlying evolutionary dynamics within and between hosts. We used high-throughput sequencing whole-genome data coupled with bottleneck analysis to retrace the pathways of viral transmission in two nosocomial outbreaks that were previously characterised by epidemiological and phylogenetic methods. Additionally, we assessed the mutational landscape, selection pressures, and diversity at the within-host level for both outbreaks. Our findings show evidence of within-host selection and transmission of variants between samples. Both bottleneck and diversity analyses highlight within-host and consensus-level variants shared by putative source-recipient pairs in both outbreaks, suggesting that certain within-host variants in these outbreaks may have been transmitted upon infection rather than arising de novo independently within multiple hosts. Overall, our findings demonstrate the utility of combining within-host diversity and bottleneck estimations for elucidating transmission events in SARS-CoV-2 outbreaks, provide insight into the maintenance of viral genetic diversity, provide a list of candidate targets of positive selection for further investigation, and demonstrate that within-host variants can be transferred between patients. Together these results will help in developing strategies to understand the nature of transmission events and curtail the spread of SARS-CoV-2.

The structure, function and evolution of a complete human chromosome 8

The complete assembly of each human chromosome is essential for understanding human biology and evolution^1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the β-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

Independent infections of porcine deltacoronavirus among Haitian children

Coronaviruses have caused three major epidemics since 2003, including the ongoing SARS-CoV-2 pandemic. In each case, the emergence of coronavirus in our species has been associated with zoonotic transmissions from animal reservoirs^1,2, underscoring how prone such pathogens are to spill over and adapt to new species. Among the four recognized genera of the family Coronaviridae, human infections reported so far have been limited to alphacoronaviruses and betacoronaviruses^3-5. Here we identify porcine deltacoronavirus strains in plasma samples of three Haitian children with acute undifferentiated febrile illness. Genomic and evolutionary analyses reveal that human infections were the result of at least two independent zoonoses of distinct viral lineages that acquired the same mutational signature in the genes encoding Nsp15 and the spike glycoprotein. In particular, structural analysis predicts that one of the changes in the spike S1 subunit, which contains the receptor-binding domain, may affect the flexibility of the protein and its binding to the host cell receptor. Our findings highlight the potential for evolutionary change and adaptation leading to human infections by coronaviruses outside of the previously recognized human-associated coronavirus groups, particularly in settings where there may be close human-animal contact.

Genome-Wide Identification of Tannase Genes and Their Function of Wound Response and Astringent Substances Accumulation in Juglandaceae

Tannins are important polyphenol compounds with different component proportions in different plant species. The plants in the Juglandaceae are rich in tannins, including condensed tannins and hydrolyzable tannins. In this study, we identified seven tannase genes (TAs) responsible for the tannin metabolism from walnut, pecan, and Chinese hickory, and three nut tree species in the Juglandaceae, which were divided into two groups. The phylogenetic and sequence analysis showed that TA genes and neighboring clade genes (TA-like genes) had similar sequences compared with other carboxylesterase genes, which may be the origin of TA genes produced by tandem repeat. TA genes also indicated higher expressions in leaf than other tissues and were quickly up-regulated at 3 h after leaf injury. During the development of the seed coat, the expression of the synthesis-related gene GGTs and the hydrolase gene TAs was continuously decreased, resulting in the decrease of tannin content in the dry sample of the seed coat of Chinese hickory. However, due to the reduction in water content during the ripening process, the tannin content in fresh sample increased, so the astringent taste was obvious at the mature stage. In addition, the CcGGTs' expression was higher than CiGGTs in the initiation of development, but CcTAs continued to be down-regulated while CiTA2a and CiTA2b were up-regulated, which may bring about the significant differences in tannin content and astringent taste between Chinese hickory and pecan. These results suggested the crucial role of TAs in wound stress of leaves and astringent ingredient accumulation in seed coats of two nut tree species in the Juglandaceae.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

The complete chloroplast genome sequence of the Dioscorea esculenta (Lour.) Burkill (Dioscoreaceae)

Dioscorea esculenta (Lour.) Burkill is an essential tuber crop with pharmacological effects in the family Dioscoreaceae. The complete chloroplast genome of D. esculenta was determined in this study. The total genome size is 153, 437 bp in length and demonstrates a typical quadripartite structure containing a large single copy (LSC, 83,628 bp) and a small single copy (SSC, 18,893 bp), separated by a pair of inverted repeats (IRa, IRb) of 25,458 bp. The GC content of the complete chloroplast genome sequence is 37.07%. A total of 131 genes were predicted including 86 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic tree analysis of 25 species belonging to the genus Dioscorea indicated that D. esculenta and D. sansibarensis were clustered into one branch.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Transposable elements (TEs) drive genome evolution and are a notable source of pathogenesis, including cancer. While CpG methylation regulates TE activity, the locus-specific methylation landscape of mobile human TEs has to date proven largely inaccessible. Here, we apply new computational tools and long-read nanopore sequencing to directly infer CpG methylation of novel and extant TE insertions in hippocampus, heart, and liver, as well as paired tumor and non-tumor liver. As opposed to an indiscriminate stochastic process, we find pronounced demethylation of young long interspersed element 1 (LINE-1) retrotransposons in cancer, often distinct to the adjacent genome and other TEs. SINE-VNTR-Alu (SVA) retrotransposons, including their internal tandem repeat-associated CpG island, are near-universally methylated. We encounter allele-specific TE methylation and demethylation of aberrantly expressed young LINE-1s in normal tissues. Finally, we recover the complete sequences of tumor-specific LINE-1 insertions and their retrotransposition hallmarks, demonstrating how long-read sequencing can simultaneously survey the epigenome and detect somatic TE mobilization.

Microbial Shifts Following Five Years of Cover Cropping and Tillage Practices in Fertile Agroecosystems

Metagenomics in agricultural research allows for searching for bioindicators of soil health to characterize changes caused by management practices. Cover cropping (CC) improves soil health by mitigating nutrient losses, yet the benefits depend on the tillage system used. Field studies searching for indicator taxa within these systems are scarce and narrow in their scope. Our goal was to identify bioindicators of soil health from microbes that were responsive to CC (three levels) and tillage (chisel tillage, no-till) treatments after five years under field conditions. We used rRNA gene-based analysis via Illumina HiSeq2500 technology with QIIME 2.0 processing to characterize the microbial communities. Our results indicated that CC and tillage differentially changed the relative abundances (RAs) of the copiotrophic and oligotrophic guilds. Corn-soybean rotations with legume-grass CC increased the RA of copiotrophic decomposers more than rotations with grass CC, whereas rotations with only bare fallows favored stress-tolerant oligotrophs, including nitrifiers and denitrifiers. Unlike bacteria, fewer indicator fungi and archaea were detected; fungi were poorly identified, and their responses were inconsistent, while the archaea RA increased under bare fallow treatments. This is primary information that allows for understanding the potential for managing the soil community compositions using cover crops to reduce nutrient losses to the environment.

Understanding the causes of errors in eukaryotic protein-coding gene prediction: a case study of primate proteomes

Background

Recent advances in sequencing technologies have led to an explosion in the number of genomes available, but accurate genome annotation remains a major challenge. The prediction of protein-coding genes in eukaryotic genomes is especially problematic, due to their complex exon–intron structures. Even the best eukaryotic gene prediction algorithms can make serious errors that will significantly affect subsequent analyses.

Results

We first investigated the prevalence of gene prediction errors in a large set of 176,478 proteins from ten primate proteomes available in public databases. Using the well-studied human proteins as a reference, a total of 82,305 potential errors were detected, including 44,001 deletions, 27,289 insertions and 11,015 mismatched segments where part of the correct protein sequence is replaced with an alternative erroneous sequence. We then focused on the mismatched sequence errors that cause particular problems for downstream applications. A detailed characterization allowed us to identify the potential causes for the gene misprediction in approximately half (5446) of these cases. As a proof-of-concept, we also developed a simple method which allowed us to propose improved sequences for 603 primate proteins.

Conclusions

Gene prediction errors in primate proteomes affect up to 50% of the sequences. Major causes of errors include undetermined genome regions, genome sequencing or assembly issues, and limitations in the models used to represent gene exon–intron structures. Nevertheless, existing genome sequences can still be exploited to improve protein sequence quality. Perspectives of the work include the characterization of other types of gene prediction errors, as well as the development of a more comprehensive algorithm for protein sequence error correction.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Comparative Analyses of the Chloroplast Genomes of Patchouli Plants and Their Relatives in Pogostemon (Lamiaceae)

Pogostemon Desf., the largest genus of the tribe Pogostemoneae (Lamiaceae), consists of ca. 80 species distributed mainly from South and Southeast Asia to China. The genus contains many patchouli plants, which are of great economic importance but taxonomically difficult. Therefore, it is necessary to characterize more chloroplast (cp) genomes for infrageneric phylogeny analyses and species identification of Pogostemon, especially for patchouli plants. In this study, we newly generated four cp genomes for three patchouli plants (i.e., Pogostemon plectranthoides Desf., P. septentrionalis C. Y. Wu et Y. C. Huang, and two cultivars of P. cablin (Blanoco) Benth.). Comparison of all samples (including online available cp genomes of P. yatabeanus (Makino) Press and P. stellatus (Lour.) Kuntze) suggested that Pogostemon cp genomes are highly conserved in terms of genome size and gene content, with a typical quadripartite circle structure. Interspecific divergence of cp genomes has been maintained at a relatively low level, though seven divergence hotspot regions were identified by stepwise window analysis. The nucleotide diversity (Pi) value was correlated significantly with gap proportion (indels), but significantly negative with GC content. Our phylogenetic analyses based on 80 protein-coding genes yielded high-resolution backbone topologies for the Lamiaceae and Pogostemon. For the overall mean substitution rates, the synonymous (d_S) and nonsynonymous (d_N) substitution rate values of protein-coding genes varied approximately threefold, while the d_N values among different functional gene groups showed a wider variation range. Overall, the cp genomes of Pogostemon will be useful for phylogenetic reconstruction, species delimitation and identification in the future.

Diversity and evolution of B-family DNA polymerases

B-family DNA polymerases (PolBs) represent the most common replicases. PolB enzymes that require RNA (or DNA) primed templates for DNA synthesis are found in all domains of life and many DNA viruses. Despite extensive research on PolBs, their origins and evolution remain enigmatic. Massive accumulation of new genomic and metagenomic data from diverse habitats as well as availability of new structural information prompted us to conduct a comprehensive analysis of the PolB sequences, structures, domain organizations, taxonomic distribution and co-occurrence in genomes. Based on phylogenetic analysis, we identified a new, widespread group of bacterial PolBs that are more closely related to the catalytically active N-terminal half of the eukaryotic PolEpsilon (PolEpsilonN) than to Escherichia coli Pol II. In Archaea, we characterized six new groups of PolBs. Two of them show close relationships with eukaryotic PolBs, the first one with PolEpsilonN, and the second one with PolAlpha, PolDelta and PolZeta. In addition, structure comparisons suggested common origin of the catalytically inactive C-terminal half of PolEpsilon (PolEpsilonC) and PolAlpha. Finally, in certain archaeal PolBs we discovered C-terminal Zn-binding domains closely related to those of PolAlpha and PolEpsilonC. Collectively, the obtained results allowed us to propose a scenario for the evolution of eukaryotic PolBs.

Alkalibaculum sporogenes sp. nov., isolated from a terrestrial mud volcano and emended description of the genus Alkalibaculum

A novel anaerobic, endospore-forming bacterium (strain M08 DMB^T) was isolated from a terrestrial mud volcano (Taman Peninsula, Russia). Cells of the strain were motile rods 1.3-2.0 µm long and 0.4 µm in diameter. The temperature range for growth was 5-42 °C, with an optimum at 30 °C. The pH range for growth was H 6.5-11.0, with an optimum at pH 8.0. Growth of strain M08 DMB^T was observed at NaCl concentrations of 0-5.0 % (w/v) with an optimum at 1.0 %. Strain M08 DMB^T utilized 3,4-dimethoxybenzoic acid, 2-methoxyphenol, carbon monoxide, glucose, fructose, mannose, xylose and yeast extract. The end product of glucose fermentation was acetate. The DNA G+C content of strain M08 DMB^T was 32.3 mol% (obtained via whole genome sequencing). The closest phylogenetic relative of strain M08 DMB^T was Alkalibaculum bacchi (family Eubacteriaceae, class Clostridia) with 95.17 % 16S rRNA gene sequence similarity. Based on the phenotypic, genotypic and phylogenetic characteristics of the isolate, strain M08 DMB^T is considered to represent a novel species of the genus Alkalibaculum, for which the name Alkalibaculum sporogenes sp. nov. is proposed. The type strain of Alkalibaculum sporogenes is M08 DMB^T (=KCTC 15840^T=VKM B-3387^T).

Reference transcriptomes and comparative analyses of six species in the threatened rosewood genus Dalbergia

Dalbergia is a pantropical genus with more than 250 species, many of which are highly threatened due to overexploitation for their rosewood timber, along with general deforestation. Many Dalbergia species have received international attention for conservation, but the lack of genomic resources for Dalbergia hinders evolutionary studies and conservation applications, which are important for adaptive management. This study produced the first reference transcriptomes for 6 Dalbergia species with different geographical origins and predicted ~ 32 to 49 K unique genes. We showed the utility of these transcriptomes by phylogenomic analyses with other Fabaceae species, estimating the divergence time of extant Dalbergia species to ~ 14.78 MYA. We detected over-representation in 13 Pfam terms including HSP, ALDH and ubiquitin families in Dalbergia. We also compared the gene families of geographically co-occurring D. cochinchinensis and D. oliveri and observed that more genes underwent positive selection and there were more diverged disease resistance proteins in the more widely distributed D. oliveri, consistent with reports that it occupies a wider ecological niche and has higher genetic diversity. We anticipate that the reference transcriptomes will facilitate future population genomics and gene-environment association studies on Dalbergia, as well as contributing to the genomic database where plants, particularly threatened ones, are currently underrepresented.

A database for risk assessment and comparative genomic analysis of foodborne Vibrio parahaemolyticus in China

Vibrio parahaemolyticus is a major foodborne pathogen worldwide. The increasing number of cases of V. parahaemolyticus infections in China indicates an urgent need to evaluate the prevalence and genetic diversity of this pathogenic bacterium. In this paper, we introduce the Foodborne Vibrio parahaemolyticus genome database (FVPGD), the first scientific database of foodborne V. parahaemolyticus distribution and genomic data in China, based on our previous investigations of V. parahaemolyticus contamination in different kinds of food samples across China from 2011 to 2016. The dataset includes records of 2,499 food samples and 643 V. parahaemolyticus strains from supermarkets and marketplaces distributed over 39 cities in China; 268 whole-genome sequences have been deposited in this database. A spatial view on the risk situations of V. parahaemolyticus contamination in different food types is provided. Additionally, the database provides a functional interface of sequence BLAST, core genome multilocus sequence typing, and phylogenetic analysis. The database will become a powerful tool for risk assessment and outbreak investigations of foodborne pathogens in China.

Metagenome-assembled genome binning methods with short reads disproportionately fail for plasmids and genomic Islands

Metagenomic methods enable the simultaneous characterization of microbial communities without time-consuming and bias-inducing culturing. Metagenome-assembled genome (MAG) binning methods aim to reassemble individual genomes from this data. However, the recovery of mobile genetic elements (MGEs), such as plasmids and genomic islands (GIs), by binning has not been well characterized. Given the association of antimicrobial resistance (AMR) genes and virulence factor (VF) genes with MGEs, studying their transmission is a public-health priority. The variable copy number and sequence composition of MGEs makes them potentially problematic for MAG binning methods. To systematically investigate this issue, we simulated a low-complexity metagenome comprising 30 GI-rich and plasmid-containing bacterial genomes. MAGs were then recovered using 12 current prediction pipelines and evaluated. While 82-94 % of chromosomes could be correctly recovered and binned, only 38-44 % of GIs and 1-29 % of plasmid sequences were found. Strikingly, no plasmid-borne VF nor AMR genes were recovered, and only 0-45 % of AMR or VF genes within GIs. We conclude that short-read MAG approaches, without further optimization, are largely ineffective for the analysis of mobile genes, including those of public-health importance, such as AMR and VF genes. We propose that researchers should explore developing methods that optimize for this issue and consider also using unassembled short reads and/or long-read approaches to more fully characterize metagenomic data.

Natranaerofaba carboxydovora gen. nov., sp. nov., an extremely haloalkaliphilic CO-utilizing acetogen from a hypersaline soda lake representing a novel deep phylogenetic lineage in the class 'Natranaerobiia'

An anaerobic enrichment with CO from sediments of hypersaline soda lakes resulted in a methane‐forming binary culture, whereby CO was utilized by a bacterium and not the methanogenic partner. The bacterial isolate ANCO1 forms a deep‐branching phylogenetic lineage at the level of a new family within the class ‘Natranaerobiia’. It is an extreme haloalkaliphilic and moderate thermophilic acetogen utilizing CO, formate, pyruvate and lactate as electron donors and thiosulfate, nitrate (reduced to ammonia) and fumarate as electron acceptors. The genome of ANCO1 encodes a full Wood–Ljungdahl pathway allowing for CO oxidation and acetogenic conversion of pyruvate. A locus encoding Nap nitrate reductase/NrfA ammonifying nitrite reductase is also present. Thiosulfate respiration is encoded by a Phs/Psr‐like operon. The organism obviously relies on Na‐based bioenergetics, since the genome encodes for the Na⁺‐Rnf complex, Na⁺‐F1F0 ATPase and Na⁺‐translocating decarboxylase. Glycine betaine serves as a compatible solute. ANCO1 has an unusual membrane polar lipid composition dominated by diethers, more common among archaea, probably a result of adaptation to multiple extremophilic conditions. Overall, ANCO1 represents a unique example of a triple extremophilic CO‐oxidizing anaerobe and is classified as a novel genus and species Natranaerofaba carboxydovora in a novel family Natranaerofabacea.

Geographical reconstruction of the SARS-CoV-2 outbreak in Lombardy (Italy) during the early phase

The first identification of autochthonous transmission of SARS‐CoV‐2 in Italy was documented by the Laboratory of Clinical Microbiology, Virology and Bioemergencies of L. Sacco Hospital (Milano, Italy) on 20th February 2020 in a 38 years old male patient, who was found positive for pneumonia at the Codogno Hospital. Thereafter Lombardy has reported the highest prevalence of COVID‐19 cases in the country, especially in Milano, Brescia and Bergamo provinces. The aim of this study was to assess the potential presence of different viral clusters belonging to the six main provinces involved in Lombardy COVID‐19 cases in order to highlight peculiar province‐dependent viral characteristics. A phylogenetic analysis was conducted on 20 full length genomes obtained from patients addressing to several Lombard hospitals from February 20th to April 4th, 2020, aligned with 41 Italian viral genome assemblies available on GISAID database as of 30th March, 2020: two main monophyletic clades, containing 8 and 53 isolates, respectively, were identified. Noteworthy, Bergamo isolates mapped inside the small clade harbouring M gene D3G mutation. The molecular clock analysis estimated a cluster divergence approximately one month before the first patient identification, supporting the hypothesis that different SARS‐CoV‐2 strains had spread worldwide at different times, but their presence became evident only in late February along with Italian epidemic emergence. Therefore, this epidemiological reconstruction suggests that virus initial circulation in Lombardy was ascribable to multiple introduction. The phylogenetic reconstruction robustness, however, will be improved when more genomic sequences are available, in order to guarantee a complete epidemiological surveillance.

Northern Italy was the most SARS‐CoV‐2 pandemic interested country area.

A phylogeographical analysis was conducted to investigate virus entry and circulation in Italy.

Two main monophyletic clades were identified, containing 8 and 53 isolates, respectively.

The estimated cluster divergence in mid January supported the hypothesis of different viral strains spreading at different times and multiple inputs Lombardy region.

Uncovering the shell game with barcodes: diversity of meiofaunal Caecidae snails (Truncatelloidea, Caenogastropoda) from Central America

Caecidae is a species-rich family of microsnails with a worldwide distribution. Typical for many groups of gastropods, caecid taxonomy is largely based on overt shell characters. However, identification of species using shell characteristics is problematic due to their rather uniform, tubular shells, the presence of different growth stages, and a high degree of intraspecific variability. In the present study, a first integrative approach to caecid taxonomy is provided using light-microscopic investigation with microsculptural analyses and multi-marker barcoding, in conjunction with molecular species delineation analyses (ABGD, haplotype networks, GMYC, and bPTP). In total 132 specimens of Caecum and Meioceras collected during several sampling trips to Central America were analyzed and delineated into a minimum of 19 species to discuss putative synonyms, and supplement the original descriptions. Molecular phylogenetic analyses suggest Meiocerasnitidum and M.cubitatum should be reclassified as Caecum, and the genus Meioceras might present a junior synonym of Caecum. Meiofaunal caecids morphologically resembling C.glabrum from the Northeast Atlantic are a complex of cryptic species with independent evolutionary origins, likely associated with multiple habitat shifts to the mesopsammic environment. Caecuminvisibile Egger & Jörger, sp. nov. is formally described based on molecular diagnostic characters. This first integrative approach towards the taxonomy of Caecidae increases the known diversity, reveals the need for a reclassification of the genus Caecum and serves as a starting point for a barcoding library of the family, thereby enabling further reliable identifications of these taxonomically challenging microsnails in future studies.

A survey of RNA secondary structural propensity encoded within human herpesvirus genomes: global comparisons and local motifs

There are nine herpesviruses known to infect humans, of which Epstein-Barr virus (EBV) is the most widely distributed (>90% of adults infected). This ubiquitous virus is implicated in a variety of cancers and autoimmune diseases. Previous analyses of the EBV genome revealed numerous regions with evidence of generating unusually stable and conserved RNA secondary structures and led to the discovery of a novel class of EBV non-coding (nc)RNAs: the stable intronic sequence (sis)RNAs. To gain a better understanding of the roles of RNA structure in EBV biology and pathogenicity, we revisit EBV using recently developed tools for genome-wide motif discovery and RNA structural characterization. This corroborated previous results and revealed novel motifs with potential functionality; one of which has been experimentally validated. Additionally, since many herpesviruses increasingly rival the seroprevalence of EBV (VZV, HHV-6 and HHV-7 being the most notable), analyses were expanded to include all sequenced human Herpesvirus RefSeq genomes, allowing for genomic comparisons. In total 10 genomes were analyzed, for EBV (types 1 and 2), HCMV, HHV-6A, HHV-6B, HHV-7, HSV-1, HSV-2, KSHV, and VZV. All resulting data were archived in the RNAStructuromeDB (https://structurome.bb.iastate.edu/herpesvirus) to make them available to a wide array of researchers.

Phylogenetic relationships of Grapsoidea and insights into the higher phylogeny of Brachyuran

Decapoda is one of the most diverse crustacean orders, and has become an important research subject. However, the phylogenetic relationships among the main lineages of Decapoda remain uncertain, especially in the order Brachyura. Herein, we sequenced the whole mitochondrial genome of V. litterata and constructed a phylogenetic tree to understand its phylogenetic relationships with other species. The results showed that the mitochondrial genome of V. litterata was generally similar to mitogenomes of Metazoa reported in the literature, with a size of 16,247 bp, 37 genes, and a control region. Both AT-skew and GC-skew were negative, indicating more abundant Cs and Ts than Gs and As. The gene arrangement of V. litterata is identical to those of Eriocheir hepuensis, Cyclograpsus granulosus, Hemigrapsus sanguineus, Helicana wuana, and Helice tientsinensis but differs from the pancrustacean ground pattern and typical arrangement of Brachyuran crabs. Phylogenetic reconstruction showed that V. litterata belongs to the Varunidae.

Genomic Variability of Phytophthora palmivora Isolates from Different Oil Palm Cultivation Regions in Colombia

Palm oil is the most consumed vegetable oil globally, and Colombia is the largest palm oil producer in South America and fourth worldwide. However, oil palm plantations in Colombia are affected by bud rot disease caused by the oomycete Phytophthora palmivora, leading to significant economic losses. Infection processes by plant pathogens involve the secretion of effector molecules, which alter the functioning or structure of host cells. Current long-read sequencing technologies provide the information needed to produce high-quality genome assemblies, enabling a comprehensive annotation of effectors. Here, we describe the development of genomic resources for P. palmivora, including a high-quality genome assembly based on long and short-read sequencing data, intraspecies variability for 12 isolates from different oil palm cultivation regions in Colombia, and a catalog of over 1,000 candidate effector proteins. A total of 45,416 genes were annotated from the new genome assembled in 2,322 contigs adding to 165.5 Mbp, which represents an improvement of two times more gene models, 33 times better contiguity, and 11 times less fragmentation compared with currently available genomic resources for the species. Analysis of nucleotide evolution in paralogs suggests a recent whole-genome duplication event. Genetic differences were identified among isolates showing variable virulence levels. We expect that these novel genomic resources contribute to the characterization of the species and the understanding of the interaction of P. palmivora with oil palm and could be further exploited as tools for the development of effective strategies for disease control.

Genome Sequence of Trichoderma lixii MUT3171, A Promising Strain for Mycoremediation of PAH-Contaminated Sites

Mono- and polycyclic aromatic hydrocarbons (PAHs) are widespread and recalcitrant pollutants that threaten both environmental and human health. By exploiting the powerful enzymatic machinery of fungi, mycoremediation in contaminated sites aims at removing a wide range of pollutants in a cost-efficient and environmentally friendly manner. Next-generation sequencing (NGS) techniques are powerful tools for understanding the molecular basis of biotransformation of PAHs by selected fungal strains, allowing genome mining to identify genetic features of biotechnological value. Trichoderma lixii MUT3171, isolated from a historically PAH-contaminated soil in Italy, can grow on phenanthrene, as a sole carbon source. Here, we report the draft genome sequence of T. lixii MUT3171 obtained with high-throughput sequencing method. The genome of T. lixii MUT3171 was compared with other 14 Trichoderma genomes, highlighting both shared and unique features that can shed a light on the biotransformation of PAHs. Moreover, the genes potentially involved in the production of important biosurfactants and bioactive molecules have been investigated. The gene repertoire of T. lixii MUT3171 indicates a high degrading potential and provides hints on putative survival strategies in a polluted environment.

Ancestral Absence of Electron Transport Chains in Patescibacteria and DPANN

Recent discoveries suggest that the candidate superphyla Patescibacteria and DPANN constitute a large fraction of the phylogenetic diversity of Bacteria and Archaea. Their small genomes and limited coding potential have been hypothesized to be ancestral adaptations to obligate symbiotic lifestyles. To test this hypothesis, we performed cell-cell association, genomic, and phylogenetic analyses on 4,829 individual cells of Bacteria and Archaea from 46 globally distributed surface and subsurface field samples. This confirmed the ubiquity and abundance of Patescibacteria and DPANN in subsurface environments, the small size of their genomes and cells, and the divergence of their gene content from other Bacteria and Archaea. Our analyses suggest that most Patescibacteria and DPANN in the studied subsurface environments do not form specific physical associations with other microorganisms. These data also suggest that their unusual genomic features and prevalent auxotrophies may be a result of ancestral, minimal cellular energy transduction mechanisms that lack respiration, thus relying solely on fermentation for energy conservation.

Complete chloroplast genome sequence of Acer nikoense (Sapindaceae)

Acer nikoense (Sapindaceae: Acer) is a deciduous tree, belonging to the Ser. Grisea of Sect. Trifoliata. Its complete genome sequence was obtained using genome Illumina pair-end sequencing data. It had a typical quadripartite structure with 155,952 bp in length, consisting of a large single-copy region (85,720 bp) and a small single-copy region (18,072 bp), as well as a pair of inverted repeats (26,080 bp). The total GC content was 37.9%. A total of 113 unique genes were annotated, including 30 tRNAs, 4 rRNAs, and 79 protein-coding genes. The phylogenetic analysis indicated that A. nikoense and A. triflorum were the most closely related.

Phylotranscriptomics in Cucurbitaceae Reveal Multiple Whole-Genome Duplications and Key Morphological and Molecular Innovations

The ability of climbing plants to grow upward along others to reach the canopy for photosynthesis is hypothesized as a key innovation in flowering plants. Most members of the Cucurbitaceae, a family containing ∼1000 species and many important crops, are climbers and have characteristic tendrils and pepo fruits. Here, we present 127 newly sequenced transcriptomes and genomes along with other datasets for a total of 136 cucurbits representing all tribes to establish a robust Cucurbitaceae phylogeny containing eight highly resolved major clades. We analyzed whole-genome duplication, diversification dynamics, and ancestral morphologies, and found that after early genome duplication event(s), a burst of diversification and morphological innovations in flower, fruit, and root characters occurred under the climate optimum in the Early Eocene. Species radiation during the Mid-Eocene Climatic Optimum also coincided with several morphological changes shared by 80% of cucurbits. We found that the cucurbit-specific tendril identity gene TEN originated from a paleo-polyploidization event at the origin of the family. Our results support the hypothesis that cucurbit diversifications were probably driven by increased genetic diversity following polyploidizations and by trait morphological innovations under paleo-climate upheavals. Our study provides a phylogenetic framework and new insights into morphological and genomic changes underlying the adaptive evolution of Cucurbitaceae.

Comparison of eight complete plastid genomes from three moss families Amblystegiaceae, Calliergonaceae and Pylaisiaceae

We sequenced and assembled eight complete plastid genomes from three closely related pleurocarpous moss families: Amblystegium serpens, Campyliadelphus stellatus, Cratoneuron filicinum, Drepanocladus aduncus, and Leptodictyum humile (Amblystegiaceae), Calliergon sarmentosum and Warnstorfia exannulata (Calliergonaceae), and Calliergonella cuspidata (Pylaisiaceae). The newly generated plastid genomes range from 124,256 to 124,819 bp, with two inverted repeat regions (9,624–9,696 bp) separated by a large single-copy region (86,422–86,924 bp) and a small single-copy region (18,430–18,514 bp). All these plastid genomes encode 116 unique genes including 82 protein-coding genes, 30 tRNA genes and four rRNAgenes. The overall GC content is between 28.6%–29.3%. Phylogenetic analysis showed that all Amblystegiaceae species Amblystegium serpens, Campyliadelphus stellatus, Cratoneuron filicinum, Drepanocladus aduncus, Leptodictyum humile, and Sanionia uncinata clustered in one clade, which is sister to the Pylaisiaceae species Calliergonella cuspidata. The two Calliergonaceae species Calliergon sarmentosum and Warnstorfia exannulata form a clade and is sister to Amblystegiaceae and Pylaisiaceae.

Within- and cross-species predictions of plant specialized metabolism genes using transfer learning

Plant specialized metabolites mediate interactions between plants and the environment and have significant agronomical/pharmaceutical value. Most genes involved in specialized metabolism (SM) are unknown because of the large number of metabolites and the challenge in differentiating SM genes from general metabolism (GM) genes. Plant models like Arabidopsis thaliana have extensive, experimentally derived annotations, whereas many non-model species do not. Here we employed a machine learning strategy, transfer learning, where knowledge from A. thaliana is transferred to predict gene functions in cultivated tomato with fewer experimentally annotated genes. The first tomato SM/GM prediction model using only tomato data performs well (F-measure = 0.74, compared with 0.5 for random and 1.0 for perfect predictions), but from manually curating 88 SM/GM genes, we found many mis-predicted entries were likely mis-annotated. When the SM/GM prediction models built with A. thaliana data were used to filter out genes where the A. thaliana-based model predictions disagreed with tomato annotations, the new tomato model trained with filtered data improved significantly (F-measure = 0.92). Our study demonstrates that SM/GM genes can be better predicted by leveraging cross-species information. Additionally, our findings provide an example for transfer learning in genomics where knowledge can be transferred from an information-rich species to an information-poor one.

Gluconobacter dominates the gut microbiome of the Asian palm civet Paradoxurus hermaphroditus that produces kopi luwak

Coffee beans derived from feces of the civet cat are used to brew coffee known as kopi luwak (the Indonesian words for coffee and palm civet, respectively), which is one of the most expensive coffees in the world owing to its limited supply and strong market demand. Recent metabolomics studies have revealed that kopi luwak metabolites differ from metabolites found in other coffee beans. To produce kopi luwak, coffee beans are first eaten by civet cats. It has been proposed that fermentation inside the civet cat digestive tract may contribute to the distinctively smooth flavor of kopi luwak, but the biological basis has not been determined. Therefore, we characterized the microbiome of civet cat feces using 16S rRNA gene sequences to determine the bacterial taxa that may influence fermentation processes related to kopi luwak. Moreover, we compared this fecal microbiome with that of 14 other animals, revealing that Gluconobacter is a genus that is, uniquely found in feces of the civet cat. We also found that Gluconobacter species have a large number of cell motility genes, which may encode flagellar proteins allowing colonization of the civet gut. In addition, genes encoding enzymes involved in the metabolism of hydrogen sulfide and sulfur-containing amino acids were over-represented in Gluconobacter. These genes may contribute to the fermentation of coffee beans in the digestive tract of civet cats.

Monoclonal antibodies for the S2 subunit of spike of SARS-CoV-1 cross-react with the newly-emerged SARS-CoV-2

Background

A novel coronavirus, SARS-CoV-2, which emerged at the end of 2019 and causes COVID-19, has resulted in worldwide human infections. While genetically distinct, SARS-CoV-1, the aetiological agent responsible for an outbreak of severe acute respiratory syndrome (SARS) in 2002–2003, utilises the same host cell receptor as SARS-CoV-2 for entry: angiotensin-converting enzyme 2 (ACE2). Parts of the SARS-CoV-1 spike glycoprotein (S protein), which interacts with ACE2, appear conserved in SARS-CoV-2.

Aim

The cross-reactivity with SARS-CoV-2 of monoclonal antibodies (mAbs) previously generated against the S protein of SARS-CoV-1 was assessed.

Methods

The SARS-CoV-2 S protein sequence was aligned to those of SARS-CoV-1, Middle East respiratory syndrome (MERS) and common-cold coronaviruses. Abilities of mAbs generated against SARS-CoV-1 S protein to bind SARS-CoV-2 or its S protein were tested with SARS-CoV-2 infected cells as well as cells expressing either the full length protein or a fragment of its S2 subunit. Quantitative ELISA was also performed to compare binding of mAbs to recombinant S protein.

Results

An immunogenic domain in the S2 subunit of SARS-CoV-1 S protein is highly conserved in SARS-CoV-2 but not in MERS and human common-cold coronaviruses. Four murine mAbs raised against this immunogenic fragment could recognise SARS-CoV-2 S protein expressed in mammalian cell lines. In particular, mAb 1A9 was demonstrated to detect S protein in SARS-CoV-2-infected cells and is suitable for use in a sandwich ELISA format.

Conclusion

The cross-reactive mAbs may serve as useful tools for SARS-CoV-2 research and for the development of diagnostic assays for COVID-19.

Third-generation sequencing-based mapping and visualization of single nucleotide polymorphism, meiotic recombination, illegitimate mutation and repeat-induced point mutation

Generation of new genetic diversity by crossover (CO) and non-crossover (NCO) is a fundamental process in eukaryotes. Fungi have played critical roles in studying this process because they permit tetrad analysis, which has been used by geneticists for several decades to determine meiotic recombination products. New genetic variations can also be generated in zygotes via illegitimate mutation (IM) and repeat-induced point mutation (RIP). RIP is a genome defense mechanism for preventing harmful expansion of transposable elements or duplicated sequences in filamentous fungi. Although the exact mechanism of RIP is unknown, the C:G to T:A mutations might result from DNA cytosine methylation. A comprehensive approach for understanding the molecular mechanisms underlying these important processes is to perform high-throughput mapping of CO, NCO, RIP and IM in zygotes bearing large numbers of heterozygous variant markers. To this aim, we developed ‘TSETA’, a versatile and user-friendly pipeline that utilizes high-quality and chromosome-level genome sequences involved in a single meiotic event of the industrial workhorse fungus Trichoderma reesei. TSETA not only can be applied to most sexual eukaryotes for genome-wide tetrad analysis, it also outcompetes most currently used methods for calling out single nucleotide polymorphisms between two or more intraspecies strains or isolates.

Chromosome Level Genome Assembly of Andrographis paniculata

Andrographis paniculata (Chinese name: Chuanxinlian) is an annual dicotyledonous medicinal plant widely grown in China and Southeast Asia. The dried plant has a highly acclaimed usage in the traditional Chinese medicine for its antipyretic, anti-inflammatory, and analgesic effects. In order to help delineate the biosynthetic pathways of various secondary metabolites, we report in this study a high-quality reference genome for A. paniculata. With the help of both PacBio single molecule real time sequencing and Illumina sequencing reads for error correction, the A. paniculata genome was assembled into a total size of 284 Mb with a contig N50 size of 5.14 Mb. The contigs were further assembled into 24 pseudo-chromosomes by the Hi-C technique. We also analyzed the gene families (e.g., KSL, and CYP450) whose protein products are essential for synthesizing bioactive compounds in A. paniculata. In conclusion, the high-quality A. paniculata genome assembly builds the foundation for decoding the biosynthetic pathways of various medicinal compounds.

Comparative Genomic Analysis of Rapidly Evolving SARS-CoV-2 Reveals Mosaic Pattern of Phylogeographical Distribution

The outbreak of coronavirus disease 2019 (COVID-19) that started in Wuhan, China, in December 2019 has spread worldwide, emerging as a global pandemic. The severe respiratory pneumonia caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has so far claimed more than 0.38 million lives and has impacted human lives worldwide. However, as the novel SARS-CoV-2 virus displays high transmission rates, the underlying genomic severity is required to be fully understood. We studied the complete genomes of 95 SARS-CoV-2 strains from different geographical regions worldwide to uncover the pattern of the spread of the virus. We show that there is no direct transmission pattern of the virus among neighboring countries, suggesting that its spread is a result of travel of infected humans to different countries. We revealed unique single nucleotide polymorphisms (SNPs) in nonstructural protein 13 (nsp13), nsp14, nsp15, and nsp16 (ORF1b polyproteins) and in the S-protein within 10 viral isolates from the United States. These viral proteins are involved in RNA replication and binding with the human receptors, indicating that the viral variants that are circulating in the population of the United States are different from those circulating in the populations of other countries. In addition, we found an amino acid addition in nsp16 (mRNA cap-1 methyltransferase) of a U.S. isolate (GenBank accession no. MT188341.1) leading to a shift in the amino acid frame from position 2540 onward. Through comparative structural analysis of the wild-type and mutant proteins, we showed that this addition of a phenylalanine residue renders the protein in the mutant less stable, which might affect mRNA cap-1 methyltransferase function. We further analyzed the SARS-CoV-2-human interactome, which revealed that the interferon signaling pathway is targeted by orf1ab during infection and that it also interacts with NF-κB-repressing factor (NKRF), which is a potential regulator of interleukin-8 (IL-8). We propose that targeting this interaction may subsequently improve the health condition of COVID-19 patients. Our analysis also emphasized that SARS-CoV-2 manipulates spliceosome machinery during infection; hence, targeting splicing might affect viral replication. In conclusion, the replicative machinery of SARS-CoV-2 is targeting interferon and the notch signaling pathway along with spliceosome machinery to evade host challenges.IMPORTANCE The COVID-19 pandemic continues to storm the world, with over 6.5 million cases worldwide. The severity of the disease varies with the territories and is mainly influenced by population density and age factor. In this study, we analyzed the transmission pattern of 95 SARS-CoV-2 genomes isolated from 11 different countries. Our study also revealed several nonsynonymous mutations in ORF1b and S-proteins and the impact on their structural stability. Our analysis showed the manipulation of host system by viral proteins through SARS-CoV-2-human protein interactome, which can be useful to understand the impact of virus on human health.

Cultivar-specific transcriptome and pan-transcriptome reconstruction of tetraploid potato

Although the reference genome of Solanum tuberosum Group Phureja double-monoploid (DM) clone is available, knowledge on the genetic diversity of the highly heterozygous tetraploid Group Tuberosum, representing most cultivated varieties, remains largely unexplored. This lack of knowledge hinders further progress in potato research. In conducted investigation, we first merged and manually curated the two existing partially-overlapping DM genome-based gene models, creating a union of genes in Phureja scaffold. Next, we compiled available and newly generated RNA-Seq datasets (cca. 1.5 billion reads) for three tetraploid potato genotypes (cultivar Désirée, cultivar Rywal, and breeding clone PW363) with diverse breeding pedigrees. Short-read transcriptomes were assembled using several de novo assemblers under different settings to test for optimal outcome. For cultivar Rywal, PacBio Iso-Seq full-length transcriptome sequencing was also performed. EvidentialGene redundancy-reducing pipeline complemented with in-house developed scripts was employed to produce accurate and complete cultivar-specific transcriptomes, as well as to attain the pan-transcriptome. The generated transcriptomes and pan-transcriptome represent a valuable resource for potato gene variability exploration, high-throughput omics analyses, and breeding programmes.

Measurement(s)	genome • RNA • sequence_assembly • transcriptome
Technology Type(s)	digital curation • RNA sequencing • sequence assembly process
Factor Type(s)	cultivar
Sample Characteristic - Organism	Solanum tuberosum

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12562733

Description of two new species of the genus Heterochelamon Türkay & Dai, 1997 (Crustacea: Decapoda: Brachyura: Potamidae), from southern China

This study describes two new species of freshwater crab of the genus Heterochelamon Türkay & Dai, 1997 from southern China, H. huidongense from Guangdong Province and H. jinxiuense from Guangxi Zhuang Autonomous Region. The two new species can be differentiated from congeners by characters derived from the shape of the epibranchial tooth, external orbital angle, cheliped proportions and structure of the male first gonopod. The present study brings the number of Heterochelamon species to seven. We used the mitochondrial 16S rRNA gene for a molecular analysis and the results are consistent with the morphological features that support the recognition of two new taxa.

Complete chloroplast genome sequence of Acer sutchuenense subsp. tienchuanenge (Sapindaceae)

Acer sutchuenense subsp. tienchuanenge (Sapindaceae: Acer) is an endangered deciduous arbor species and endemic to China. Being obtained by using genome Illumina pair-end sequencing data, the complete chloroplast genome of A. sutchuenense subsp. tienchuanenge had a typical quadripartite structure, with 156,063 bp long, including a large single-copy (LSC) region of 85,772 bp, a small single-copy (SSC) region of 18,117 bp, and a pair of inverted repeats (IRs) (each 26,087 bp in length). A total of 136 genes were annotated, of which 113 are unique genes, including 30 tRNAs, 4 rRNAs, and 79 protein-coding genes. The overall GC content was 37.9%. The phylogenetic analysis suggested that A. sutchuenense subsp. tienchuanenge was the most closely related to A. griseum and A. triflorum. The complete chloroplast genome of A. sutchuenense subsp. tienchuanenge is valuable for assessment and conservation of genetic resources and further for phylogenetic study of Acer L.

Tailor-made multiple sequence alignments using the PRALINE 2 alignment toolkit

Summary

PRALINE 2 is a toolkit for custom multiple sequence alignment workflows. It can be used to incorporate sequence annotations, such as secondary structure or (DNA) motifs, into the alignment scoring, as well as to customize many other aspects of a progressive multiple alignment workflow.

Availability and implementation

PRALINE 2 is implemented in Python and available as open source software on GitHub: https://github.com/ibivu/PRALINE/.

The complete chloroplast genome of Acanthopanax brachypus (Araliaceae), a famous medicinal plant in China

Acanthopanax brachypus (Araliaceae) is an important medicinal plant originated from China. Here, we reported the first chloroplast genome sequence of A. brachypus. The plastome of A. brachypu is 156,802 bp in length, including a large single-copy (LSC) of 86,742 bp, a small single-copy (SSC) of 18,184 bp, and a pair of inverted repeat regions (IRa and IRb) of 25,938 bp. It contains 113 unique genes consisting of 80 protein-coding genes, 29 tRNA genes, and 4 rRNA genes.