A simple method to control over-alignment in the MAFFT multiple sequence alignment program

Factors influencing estimates of HIV-1 infection timing using BEAST

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3–9 days) and gag (IQR: 5–9 days), whilst the genome placed it at a median of 10 days (IQR: 4–19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

Molecular dating using phylogenetics allows us to estimate the date of an infection from time-stamped within-host sequences alone. There are large datasets of HIV-1 sequences, but genome and gene analyses are not often performed in parallel and rarely with the possibility to compare results against a known narrow window of infection. We showed that all but the longest genes are near-clonal in acute infection, with little information for dating purposes. For infections with single founders, we estimated the eclipse phase—the time between HIV-1 exposure and the first positive diagnostic test—to last between one and two weeks using env, gag, pol and near full-length genomes. This approach could be used to narrow the date of suspected infection in ongoing clinical trials for the prevention of HIV-1 infection.

Genome and time-of-day transcriptome of Wolffia australiana link morphological minimization with gene loss and less growth control

Rootless plants in the genus Wolffia are some of the fastest growing known plants on Earth. Wolffia have a reduced body plan, primarily multiplying through a budding type of asexual reproduction. Here, we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas, which has the smallest genome size in the genus at 357 Mb and has a reduced set of predicted protein-coding genes at about 15,000. Comparison between multiple high-quality draft genome sequences from W. australiana clones confirmed loss of several hundred genes that are highly conserved among flowering plants, including genes involved in root developmental and light signaling pathways. Wolffia has also lost most of the conserved nucleotide-binding leucine-rich repeat (NLR) genes that are known to be involved in innate immunity, as well as those involved in terpene biosynthesis, while having a significant overrepresentation of genes in the sphingolipid pathways that may signify an alternative defense system. Diurnal expression analysis revealed that only 13% of Wolffia genes are expressed in a time-of-day (TOD) fashion, which is less than the typical ∼40% found in several model plants under the same condition. In contrast to the model plants Arabidopsis and rice, many of the pathways associated with multicellular and developmental processes are not under TOD control in W. australiana, where genes that cycle the conditions tested predominantly have carbon processing and chloroplast-related functions. The Wolffia genome and TOD expression data set thus provide insight into the interplay between a streamlined plant body plan and optimized growth.

Pleistocene aridification underlies the evolutionary history of the Caribbean endemic, insular, giant Consolea (Opuntioideae)

PREMISE

The Caribbean islands are in the top five biodiversity hotspots on the planet; however, the biogeographic history of the seasonally dry tropical forest (SDTF) there is poorly studied. Consolea consists of nine species of dioecious, hummingbird-pollinated tree cacti endemic to the West Indies, which form a conspicuous element of the SDTF. Several species are threatened by anthropogenic disturbance, disease, sea-level rise, and invasive species and are of conservation concern. However, no comprehensive phylogeny yet exists for the clade.

METHODS

We reconstructed the phylogeny of Consolea, sampling all species using plastomic data to determine relationships, understand the evolution of key morphological characters, and test their biogeographic history. We estimated divergence times to determine the role climate change may have played in shaping the current diversity of the clade.

RESULTS

Consolea appears to have evolved very recently during the latter part of the Pleistocene on Cuba/Hispaniola likely from a South American ancestor and, from there, moved into the Bahamas, Jamaica, Puerto Rico, Florida, and the Lesser Antilles. The tree growth form is a synapomorphy of Consolea and likely aided in the establishment and diversification of the clade.

CONCLUSIONS

Pleistocene aridification associated with glaciation likely played a role in shaping the current diversity of Consolea, and insular gigantism may have been a key innovation leading to the success of these species to invade the often-dense SDTF. This in-situ Caribbean radiation provides a window into the generation of species diversity and the complexity of the SDTF community within the Antilles.

Unexpected mitochondrial lineage diversity within the genus Alonella Sars, 1862 (Crustacea: Cladocera) across the Northern Hemisphere

Representatives of the genus Alonella Sars (Crustacea: Cladocera: Chydorinae) belong to the smallest known water fleas. Although species of Alonella are widely distributed and often abundant in acidic and mountain water bodies, their diversity is poorly studied. Morphological and genetic approaches have been complicated by the minute size of these microcrustaceans. As a result, taxonomists have avoided revising these species. Here, we present genetic data on Alonella species diversity across the Northern Hemisphere with particular attention to the A. excisa species complex. We analyzed 82 16S rRNA sequences (all newly obtained), and 78 COI sequences (39 were newly obtained). The results revealed at least twelve divergent phylogenetic lineages, possible cryptic species, of Alonella, with different distribution patterns. As expected, the potential species diversity of this genus is significantly higher than traditionally accepted. The A. excisa complex is represented by nine divergent clades in the Northern Hemisphere, some of them have relatively broad distribution ranges and others are more locally distributed. Our results provide a genetic background for subsequent morphological analyses, formal descriptions of Alonella species and detailed phylogeographical studies.

Lasiodiplodia syzygii sp. nov. (Botryosphaeriaceae) causing post-harvest water-soaked brown lesions on Syzygium samarangense in Chiang Rai, Thailand

Background

Syzygium samarangense (Wax apple) is an important tropical fruit tree with high economic and nutrient value and is widely planted in the tropics or subtropics of Asia. Post-harvest water-soaked brown lesions were observed on mature fruits of ornamental wax apples in Chiang Rai Province, Thailand. A fungus with morphological characters, similar to Lasiodiplodia, was consistently isolated from symptomatic fruits. Phylogenetic analyses, based on ITS, LSU, TEF1-a and tub2, revealed that our isolates were closely related to, but phylogenetically distinct from, Lasiodiplodia rubropurpurea.

New information

Morphological comparisons indicated that pycnidia and conidiogenous cells of our strains were significantly larger than L. rubropurpurea. Comparisons of base-pair differences in the four loci confirmed that the species from wax apple was distinct from L. rubropurpurea and a new species, L. syzygii sp. nov., is introduced to accommodate it. Pathogenicity tests confirmed the newly-introduced species as the pathogen of this post-harvest water-soaked brown lesion disease on wax apples.

One new species and one new record of Zasmidium in China

Background

Two hyphomycetous species were collected from leaves of Smilaxchina (Liliales, Smilacaceae) and Cremastraappendiculata (Asparagales, Orchidaceae). ITS barcoding indicated that they belong to the genus Zasmidium.

New information

Morphological data in combination with molecular phylogenetic analyses based on ITS, LSU and rpb2 confirmed that our Chinese strains represented a new species, Zasmidiumliboense and a new record of Z.citri-griseum.

Phylogenetics of Leptocereus (Cactaceae) on Hispaniola: clarifying species limits in the L. weingartianus complex and a new species from the Sierra de Bahoruco

The Antillean genus Leptocereus represents an in-situ radiation among the Greater and Lesser Antilles of 19 currently recognized species. Extensive fieldwork carried out in the Dominican Republic over recent years has revealed that the species limits of Leptocereus of Hispaniola are more complex than previously thought. There are four currently recognized species that occur on the island, L. demissus, L. paniculatus, L. undulosus and L. weingartianus. We evaluate species limits in this group based on DNA sequence data and phylogenetic analysis, morphological characters and a survey of herbarium specimens from across Hispaniola. Based on our analyses, it is clear that at least five species occur on the island of Hispaniola, with the new species from Sierra de Bahoruco, L. velozianus, described here. We provide an identification key, distribution maps and photographic plates for all species on Hispaniola based on our own fieldwork and the study of herbarium specimens. The description of yet another species of Leptocereus on Hispaniola reiterates the importance of the poorly studied, but yet biodiverse, seasonally dry tropical forest in the Antilles.

Fast and Accurate Multiple Sequence Alignment with MSAProbs-MPI

Multiple sequence alignment (MSA) is a central step in many bioinformatics and computational biology analyses. Although there exist many methods to perform MSA, most of them fail when dealing with large datasets due to their high computational cost. MSAProbs-MPI is a publicly available tool ( http://msaprobs.sourceforge.net ) that provides highly accurate results in relatively short runtime thanks to exploiting the hardware resources of multicore clusters. In this chapter, I explain the statistical and biological concepts employed in MSAProbs-MPI to complete the alignments, as well as the high-performance computing techniques used to accelerate it. Moreover, I provide some hints about the configuration parameters that should be used to guarantee high-performance executions.

Multiple Recombination Events and Strong Purifying Selection at the Origin of SARS-CoV-2 Spike Glycoprotein Increased Correlated Dynamic Movements

Our evolutionary and structural analyses revealed that the severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) spike gene is a complex mosaic resulting from several recombination events. Additionally, the fixation of variants has mainly been driven by purifying selection, suggesting the presence of conserved structural features. Our dynamic simulations identified two main long-range covariant dynamic movements of the novel glycoprotein, and showed that, as a result of the evolutionary duality, they are preserved. The first movement involves the receptor binding domain with the N-terminal domain and the C-terminal domain 2 and is maintained across human, bat and pangolin coronaviruses. The second is a complex network of long-range dynamics specific to SARS-CoV-2 involving the novel PRRA and the conserved KR*SF cleavage sites, as well as conserved segments in C-terminal domain 3. These movements, essential for host cell binding, are maintained by hinges conserved across human, bat, and pangolin coronaviruses glycoproteins. The hinges, located around Threonine 333 and Proline 527 within the N-terminal domain and C-terminal domain 2, represent candidate targets for the future development of novel pan-coronavirus inhibitors. In summary, we show that while recombination created a new configuration that increased the covariant dynamic movements of the SARS-CoV-2 glycoprotein, negative selection preserved its inter-domain structure throughout evolution in different hosts and inter-species transmissions.

Genomic characterization of SARS-CoV-2 in Egypt

Introduction

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread throughout the globe, causing a pandemic. In Egypt over 115,000 individuals were infected so far.

Objective

In the present study, the objective is to perform a complete genome sequence of SAR-CoV2 isolated from Egyptian coronavirus disease (COVID-19) patients.

Methods

Nasopharyngeal swabs were collected from 61 COVID-19 patients who attended at National Cancer Institute, Kasr Al-Aini Hospital and the army hospital. Viral RNA was extracted and whole genomic sequencing was conducted using Next Generation Sequencing.

Results

In all cases, the sequenced virus has at least 99% identity to the reference Wuhan 1. The sequence analysis showed 204 distinct genome variations including 114 missense mutations, 72 synonymous mutations, 1 disruptive in-frame deletion, 7 downstream gene mutations, 6 upstream gene mutations, 3 frame-shift deletions, and 1 in-frame deletion. The most dominant clades were G/GH/GR/O and the dominant type is B.

Conclusion

The whole genomic sequence of SARS-CoV2 showed 204 variations in the genomes of the Egyptian isolates, where the Asp614Gly (D614G) substitution is the most common among the samples (60/61). So far, there were no strikingly variations specific to the Egyptian population, at least for this set of samples.

Assembly and Annotation of the Nuclear and Organellar Genomes of a Dwarf Coconut (Chowghat Green Dwarf) Possessing Enhanced Disease Resistance

Coconut (Cocos nucifera L.), an important source of vegetable oil, nutraceuticals, functional foods, and housing materials, provides raw materials for a repertoire of industries engaged in the manufacture of cosmetics, soaps, detergents, paints, varnishes, and emulsifiers, among other products. The palm plays a vital role in maintaining and promoting the sustainability of farming systems of the fragile ecosystems of islands and coastal regions of the tropics. In this study, we present the genome of a dwarf coconut variety "Chowghat Green Dwarf" (CGD) from India, possessing enhanced resistance to root (wilt) disease. Utilizing short reads from the Illumina HiSeq 4000 platform and long reads from the Pacific Biosciences RSII platform, we have assembled the draft genome assembly of 1.93 Gb. The genome is distributed over 26,855 scaffolds, with ∼81.56% of the assembled genome present in scaffolds of lengths longer than 50 kb. About 77.29% of the genome was composed of transposable elements and repeats. Gene prediction yielded 51,953 genes, which upon stringent filtering, based on Annotation Edit Distance, resulted in 13,707 genes, which coded for 11,181 proteins. Among these, we gathered transcript level evidence for a total of 6828 predicted genes based on the RNA-Seq data from different coconut tissues, since they presented assembled transcripts within the genome annotation coordinates. A total of 112 nucleotide-binding and leucine-rich repeat loci, belonging to six classes, were detected. We have also undertaken the assembly and annotation of the CGD chloroplast and mitochondrial genomes. The availability of the dwarf coconut genome shall prove invaluable for deducing the origin of dwarf coconut cultivars, dissection of genes controlling plant habit and fruit color, and accelerated breeding for improved agronomic traits.

Community composition of arctic root-associated fungi mirrors host plant phylogeny

The number of plant species regarded as non-mycorrhizal increases at higher latitudes, and several plant species in the High-Arctic Archipelago Svalbard have been reported as non-mycorrhizal. We used the rRNA ITS2 and 18S gene markers to survey which fungi, as well as other micro-eukaryotes, were associated with roots of 31 arctic plant species not usually regarded as mycorrhizal in Svalbard. We assessed to what degree the root-associated fungi showed any host preference and whether the phylogeny of the plant hosts may mirror the composition of root-associated fungi. Fungal communities were largely structured according to host plant identity and to a less extent by environmental factors. We observed a positive relationship between the phylogenetic distance of host plants and the distance of fungal community composition between samples, indicating that the evolutionary history of the host plants plays a major role for which fungi colonize the plant roots. In contrast to the ITS2 marker, the 18S rRNA gene marker showed that chytrid fungi were prevalently associated with plant roots, together with a wide spectrum of amoeba-like protists and nematodes. Our study confirms that arbuscular mycorrhizal (AM) fungi are present also in arctic environments in low abundance.

Dickeya poaceiphila sp. nov., a plant-pathogenic bacterium isolated from sugar cane (Saccharum officinarum)

The genus Dickeya is an important group of plant pathogens that currently comprises 10 recognized species. Although most Dickeya isolates originated from infected cultivated plants, they are also isolated from water. The genomic sequence of the Australian strain NCPPB 569^T clearly established its separation from the previously characterized Dickeya species. The average nucleotide identity and digital DNA-DNA hybridization values obtained by comparing strain NCPPB 569^T with strains of characterized Dickeya species were lower than 87 and 32 %, respectively, supporting the delineation of a new species. The name Dickeya poaceiphila sp. nov. is proposed for this taxon with the type strain NCPPB 569^T (=CFBP 8731^T). Two other strains isolated in Australia, CFBP 1537 and CFBP 2040, also belong to this species. Phenotypic and genomic comparisons enabled the identification of traits distinguishing D. poaceiphila isolates from strains of other Dickeya species.

Identification and Distribution of Novel Cressdnaviruses and Circular molecules in Four Penguin Species in South Georgia and the Antarctic Peninsula

There is growing interest in uncovering the viral diversity present in wild animal species. The remote Antarctic region is home to a wealth of uncovered microbial diversity, some of which is associated with its megafauna, including penguin species, the dominant avian biota. Penguins interface with a number of other biota in their roles as marine mesopredators and several species overlap in their ranges and habitats. To characterize the circular single-stranded viruses related to those in the phylum Cressdnaviricota from these environmental sentinel species, cloacal swabs (n = 95) were obtained from King Penguins in South Georgia, and congeneric Adélie Penguins, Chinstrap Penguins, and Gentoo Penguins across the South Shetland Islands and Antarctic Peninsula. Using a combination of high-throughput sequencing, abutting primers-based PCR recovery of circular genomic elements, cloning, and Sanger sequencing, we detected 97 novel sequences comprising 40 ssDNA viral genomes and 57 viral-like circular molecules from 45 individual penguins. We present their detection patterns, with Chinstrap Penguins harboring the highest number of new sequences. The novel Antarctic viruses identified appear to be host-specific, while one circular molecule was shared between sympatric Chinstrap and Gentoo Penguins. We also report viral genotype sharing between three adult-chick pairs, one in each Pygoscelid species. Sequence similarity network approaches coupled with Maximum likelihood phylogenies of the clusters indicate the 40 novel viral genomes do not fall within any known viral families and likely fall within the recently established phylum Cressdnaviricota based on their replication-associated protein sequences. Similarly, 83 capsid protein sequences encoded by the viruses or viral-like circular molecules identified in this study do not cluster with any of those encoded by classified viral groups. Further research is warranted to expand knowledge of the Antarctic virome and would help elucidate the importance of viral-like molecules in vertebrate host evolution.

Genotypic characterization and molecular evolution of avian reovirus in poultry flocks from Brazil

Avian reovirus (ARV) is one of the main causes of infectious arthritis/tenosynovitis and malabsorption syndrome (MAS) in poultry. ARVs have been disseminated in Brazilian poultry flocks in the last years. This study aimed to genotype ARVs and to evaluate the molecular evolution of the more frequent ARV lineages detected in Brazilian poultry-producing farms. A total of 100 poultry flocks with clinical signs of tenosynovitis/MAS, from all Brazilian poultry-producing regions were positive for ARV by PCR. Seventeen bird tissues were submitted to cell culture and ARV RNA detection/genotyping by two PCRs. The phylogenetic classification was based on σC gene alignment using a dataset with other Brazilian and worldwide ARVs sequences. ARVs were specifically detected by both PCRs from the 17 cell cultures, and σC gene partial fragments were sequenced. All these sequences were aligned with a total of 451 ARV σC gene data available in GenBank. Phylogenetic analysis demonstrated five well-defined clusters that were classified into lineages I, II, III, IV, and V. Three lineages could be further divided into sub-lineages: I (I vaccine, Ia, Ib), II (IIa, IIb, IIc) and IV (IVa and IVb). Brazilian ARVs were from four lineages/sub-lineages: Ib (48.2%), IIb (22.2%), III (3.7%) and V (25.9%). The Bayesian analysis demonstrated that the most frequent sub-lineage Ib emerged in the world around 1968 and it was introduced into Brazil in 2010, with increasing spread soon after. In conclusion, four different ARV lineages are circulating in Brazilian poultry flocks, all associated with clinical diseases. RESEARCH HIGHLIGHTS One-hundred ARV-positive flocks were detected in all main poultry-producing regions from Brazil. A large dataset of 468 S1 sequences was constructed and divided ARVs into five lineages. Four lineages/sub-lineages (Ib, IIb, III and V) were detected in commercial poultry flocks from Brazil. Brazilian lineages shared a low identity with the commercial vaccine lineage (I vaccine). Sub-lineage Ib emerged around 1968 and was introduced into Brazil in 2010.

Structure of the bacterial ribosome at 2 Å resolution

Using cryo-electron microscopy (cryo-EM), we determined the structure of the Escherichia coli 70S ribosome with a global resolution of 2.0 Å. The maps reveal unambiguous positioning of protein and RNA residues, their detailed chemical interactions, and chemical modifications. Notable features include the first examples of isopeptide and thioamide backbone substitutions in ribosomal proteins, the former likely conserved in all domains of life. The maps also reveal extensive solvation of the small (30S) ribosomal subunit, and interactions with A-site and P-site tRNAs, mRNA, and the antibiotic paromomycin. The maps and models of the bacterial ribosome presented here now allow a deeper phylogenetic analysis of ribosomal components including structural conservation to the level of solvation. The high quality of the maps should enable future structural analyses of the chemical basis for translation and aid the development of robust tools for cryo-EM structure modeling and refinement.

Structure of the bacterial ribosome at 2 Å resolution

Using cryo-electron microscopy (cryo-EM), we determined the structure of the Escherichia coli 70S ribosome with a global resolution of 2.0 Å. The maps reveal unambiguous positioning of protein and RNA residues, their detailed chemical interactions, and chemical modifications. Notable features include the first examples of isopeptide and thioamide backbone substitutions in ribosomal proteins, the former likely conserved in all domains of life. The maps also reveal extensive solvation of the small (30S) ribosomal subunit, and interactions with A-site and P-site tRNAs, mRNA, and the antibiotic paromomycin. The maps and models of the bacterial ribosome presented here now allow a deeper phylogenetic analysis of ribosomal components including structural conservation to the level of solvation. The high quality of the maps should enable future structural analyses of the chemical basis for translation and aid the development of robust tools for cryo-EM structure modeling and refinement.

Entoloma sequestratum, a new species from northern Thailand, and a worldwide key to sequestrate taxa of Entoloma (Entolomataceae)

Based on our study of the morphology and genetics of sporocarps collected in the mountains of northern Thailand, we herein describe Entoloma sequestratum as a new sequestrate member of the Entolomotaceae. This serves as the first report of a sequestrate member of the genus from Thailand. In addition, we provide a worldwide key to all of the described sequestrate members of the genus.

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.

Intragenic Conflict in Phylogenomic Data Sets

Most phylogenetic analyses assume that a single evolutionary history underlies one gene. However, both biological processes and errors can cause intragenic conflict. The extent to which this conflict is present in empirical data sets is not well documented, but if common, could have far-reaching implications for phylogenetic analyses. We examined several large phylogenomic data sets from diverse taxa using a fast and simple method to identify well-supported intragenic conflict. We found conflict to be highly variable between data sets, from 1% to >92% of genes investigated. We analyzed four exemplar genes in detail and analyzed simulated data under several scenarios. Our results suggest that alignment error may be one major source of conflict, but other conflicts remain unexplained and may represent biological signal or other errors. Whether as part of data analysis pipelines or to explore biologically processes, analyses of within-gene phylogenetic signal should become common.

Fragmented mitochondrial genomes evolved in opposite directions between closely related macaque louse Pedicinus obtusus and colobus louse Pedicinus badii

We report for the first time the fragmented mitochondrial (mt) genomes of two Pedicinus species: Pedicinus obtusus and Pedicinus badii, and compared them with the lice of humans and chimpanzees. Despite being congeneric, the two monkey lice are distinct from each other in mt karyotype. The variation in mt karyotype between the two Pedicinus lice is the most pronounced among the congeneric species of sucking lice observed to date and is attributable to the opposite directions between them in mt karyotype evolution. Two of the inferred ancestral mt minichromosomes of the higher primate lice merged as one in the macaque louse whereas one of the ancestral minichromosomes split into two in the colobus louse after these two species diverged from their most recent common ancestor. Our results showed that mt genome fragmentation was a two-way process in the higher primate lice, and minichromosome merger was more common than previously thought.

Viral diversity in oral cavity from Sapajus nigritus by metagenomic analyses

Sapajus nigritus are non-human primates which are widespread in South America. They are omnivores and live in troops of up to 40 individuals. The oral cavity is one of the main entry routes for microorganisms, including viruses. Our study proposed the identification of viral sequences from oral swabs collected in a group of capuchin monkeys (n = 5) living in a public park in a fragment of Mata Atlantica in South Brazil. Samples were submitted to nucleic acid extraction and enrichment, which was followed by the construction of libraries. After high-throughput sequencing and contig assembly, we used a pipeline to identify 11 viral families, which are Herpesviridae, Parvoviridae, Papillomaviridae, Polyomaviridae, Caulimoviridae, Iridoviridae, Astroviridae, Poxviridae, and Baculoviridae, in addition to two complete viral genomes of Anelloviridae and Genomoviridae. Some of these viruses were closely related to known viruses, while other fragments are more distantly related, with 50% of identity or less to the currently available virus sequences in databases. In addition to host-related viruses, insect and small vertebrate-related viruses were also found, as well as plant-related viruses, bringing insights about their diet. In conclusion, this viral metagenomic analysis reveals, for the first time, the profile of viruses in the oral cavity of wild, free ranging capuchin monkeys.

Molecular phylogenetic analysis and morphological reassessments of thief ants identify a new potential case of biological invasions

Species delimitation offered by DNA-based approaches can provide important insights into the natural history and diversity of species, but the cogency of such processes is limited without multigene phylogenies. Recent attempts to barcode various Solenopsidini ant taxa (Hymenoptera: Formicidae: Myrmicinae), including the thief ant Solenopsis saudiensis Sharaf & Aldawood, 2011 described from the Kingdom of Saudi Arabia (KSA), were precipitated by the unexpected existence of a closely related species, the Nearctic S. abdita Thompson, 1989 within the S. molesta species complex native to Florida. This finding left the species status of the former uncertain. Here, we investigated the taxonomy and phylogeny of these two species to determine whether or not S. abdita represents a new global tramp species. We inferred a phylogeny of the two species using DNA sequence data from four nuclear genes (Abd-A, EF1α-F1, EF1α-F2, and Wingless) and one mitochondrial gene (COI) sampled from populations in Florida, Guatemala, Hawaii, and Saudi Arabia. Both species clustered into one distinct and robust clade. The taxonomy of S. saudiensis was re‐examined using morphometrics. A reassessment of the morphological characters used to diagnose the worker and queen castes were consistent with molecular evidence. Based on combined morphological and molecular evidences S. saudiensis is declared as a junior synonym of S. abdita (syn. nov.). In addition, our findings indicate that S. abdita is a novel global tramp species which has a far wider distribution than previously thought and has established itself in many new habitats and different geographic realms.

Groundwater Elusimicrobia are metabolically diverse compared to gut microbiome Elusimicrobia and some have a novel nitrogenase paralog

Currently described members of Elusimicrobia, a relatively recently defined phylum, are animal-associated and rely on fermentation. However, free-living Elusimicrobia have been detected in sediments, soils and groundwater, raising questions regarding their metabolic capacities and evolutionary relationship to animal-associated species. Here, we analyzed 94 draft-quality, non-redundant genomes, including 30 newly reconstructed genomes, from diverse animal-associated and natural environments. Genomes group into 12 clades, 10 of which previously lacked reference genomes. Groundwater-associated Elusimicrobia are predicted to be capable of heterotrophic or autotrophic lifestyles, reliant on oxygen or nitrate/nitrite-dependent respiration, or a variety of organic compounds and Rhodobacter nitrogen fixation (Rnf) complex-dependent acetogenesis with hydrogen and carbon dioxide as the substrates. Genomes from two clades of groundwater-associated Elusimicrobia often encode a new group of nitrogenase paralogs that co-occur with an extensive suite of radical S-Adenosylmethionine (SAM) proteins. We identified similar genomic loci in genomes of bacteria from the Gracilibacteria phylum and the Myxococcales order and predict that the gene clusters reduce a tetrapyrrole, possibly to form a novel cofactor. The animal-associated Elusimicrobia clades nest phylogenetically within two free-living-associated clades. Thus, we propose an evolutionary trajectory in which some Elusimicrobia adapted to animal-associated lifestyles from free-living species via genome reduction.

Characterization of the Role of a Non-GPCR Membrane-Bound CFEM Protein in the Pathogenicity and Germination of Botrytis cinerea

The necrotrophic fungus Botrytis cinerea, is considered a major cause of postharvest losses in a wide range of crops. The common fungal extracellular membrane protein (CFEM), containing a conserved eight-cysteine pattern, was found exclusively in fungi. Previous studies in phytopathogenic fungi have demonstrated the role of membrane-bound and secreted CFEM-containing proteins in different aspects of fungal virulence. However, non-G protein-coupled receptor (non-GPCR) membrane CFEM proteins have not been studied yet in phytopathogenic fungi. In the present study, we have identified a non-GPCR membrane-bound CFEM-containing protein, Bcin07g03260, in the B. cinerea genome, and generated deletion mutants, ΔCFEM-Bcin07g03260, to study its potential role in physiology and virulence. Three independent ΔCFEM-Bcin07g03260 mutants showed significantly reduced progression of a necrotic lesion on tomato (Solanum lycopersicum) leaves. Further analysis of the mutants revealed significant reduction (approximately 20–30%) in conidial germination and consequent germ tube elongation compared with the WT. Our data complements a previous study of secreted ΔCFEM1 mutants of B. cinerea that showed reduced progression of necrotic lesions on leaves, without effect on germination. Considering various functions identified for CFEM proteins in fungal virulence, our work illustrates a potential new role for a non-GPCR membrane CFEM in pathogenic fungi to control virulence in the fungus B. cinerea.

Bat-borne polyomaviruses in Europe reveal an evolutionary history of intrahost divergence with horseshoe bats distributed across the African and Eurasian continents

Polyomaviruses (PyVs) are small, circular dsDNA viruses carried by diverse vertebrates, including bats. Although previous studies have reported several horseshoe bat PyVs collected in Zambia and China, it is still unclear how PyVs evolved in this group of widely dispersed mammals. Horseshoe bats (genus Rhinolophus) are distributed across the Old World and are natural reservoirs of numerous pathogenic viruses. Herein, non-invasive bat samples from European horseshoe bat species were collected in Hungary for PyV identification and novel PyVs with complete genomes were successfully recovered from two different European horseshoe bat species. Genomic and phylogenetic analysis of the Hungarian horseshoe bat PyVs supported their classification into the genera Alphapolyomavirus and Betapolyomavirus. Notably, despite the significant geographical distances between the corresponding sampling locations, Hungarian PyVs exhibited high genetic relatedness with previously described Zambian and Chinese horseshoe bat PyVs, and phylogenetically clustered with these viruses in each PyV genus. Correlation and virus-host relationship analysis suggested that these PyVs co-diverged with their European, African and Asian horseshoe bat hosts distributed on different continents during their evolutionary history. Additionally, assessment of selective pressures over the major capsid protein (VP1) of horseshoe bat PyVs showed sites under positive selection located in motifs exposed to the exterior of the capsid. In summary, our findings revealed a pattern of stable intrahost divergence of horseshoe bat PyVs with their mammalian hosts on the African and Eurasian continents over evolutionary time.

The complete chloroplast genome of Caryota obtusa, an endangered and economically important species

Caryota obtusa is an endangered and economically important species of the Arecaceae. The complete chloroplast genome sequence of this species is a circular molecule of 159,882 bp in size, including a pair of inverted repeats with length of 27,271 bp, separated by a large single-copy (87,645 bp) region and a small single-copy region (17,695 bp). In total, there are 131 genes, encoding 79 protein-coding genes, 40 tRNAs, and 10 rRNA genes, in which 123 genes, 69 CDSs, 37 tRNAs, and 10 rRNAs are unique, respectively. Phylogenetic inference confirmed the monophyly of the Caryota genus and its delimitation in subfamily Coryphoideae.

The complete chloroplast genome sequence of Picea schrenkiana (Pinaceae)

Picea is a phylogenetically complicate genus with great economic and ecological values. Here, we determined the whole complete chloroplast genome of Picea schrenkiana to provide genomic information for phylogenetic analysis of the genus. The plastome of P. schrenkiana is 124,060 bp in size and contains 114 genes, including 74 protein-coding genes, 36 tRNA genes, and four rRNA genes. The overall GC content is 38.7%. Unlike the typical plastome with a conserved quadripartite structure, loss of inverted repeat regions is found in the chloroplast genome. The phylogenetic tree shows that monophyly of P. schrenkiana is well supported.

Genetic Characterization of Influenza A Viruses in Japanese Swine in 2015 to 2019

Understanding the current status of influenza A viruses of swine (IAVs-S) and their evolution at the farm level is important for controlling these pathogens. Efforts to monitor IAVs-S during 2015 to 2019 yielded H1N1, H1N2, and H3N2 viruses. H1 genes in Japanese swine formed a unique clade in the classical swine H1 lineage of 1A.1, and H3 genes originating from 1999–2000 human seasonal influenza viruses appear to have become established among Japanese swine. A(H1N1)pdm09-derived H1 genes became introduced repeatedly and reassorted with endemic IAVs-S, resulting in various combinations of surface and internal genes among pig populations in Japan. At the farm level, multiple introductions of IAVs-S with phylogenetically distinct HA sequences occurred, or IAVs-S derived from a single introduction have persisted for at least 3 years with only a single mutation at the antigenic site of the HA protein. Continued monitoring of IAVs-S is necessary to update and maximize control strategies.

To assess the current status of influenza A viruses of swine (IAVs-S) throughout Japan and to investigate how these viruses persisted and evolve on pig farms, we genetically characterized IAVs-S isolated during 2015 to 2019. Nasal swab samples collected through active surveillance and lung tissue samples collected for diagnosis yielded 424 IAVs-S, comprising 78 H1N1, 331 H1N2, and 15 H3N2 viruses, from farms in 21 sampled prefectures in Japan. Phylogenetic analyses of surface genes revealed that the 1A.1 classical swine H1 lineage has evolved uniquely since the late 1970s among pig populations in Japan. During 2015 to 2019, A(H1N1)pdm09 viruses repeatedly became introduced into farms and reassorted with endemic H1N2 and H3N2 IAVs-S. H3N2 IAVs-S isolated during 2015 to 2019 formed a clade that originated from 1999–2000 human seasonal influenza viruses; this situation differs from previous reports, in which H3N2 IAVs-S derived from human seasonal influenza viruses were transmitted sporadically from humans to swine but then disappeared without becoming established within the pig population. At farms where IAVs-S were frequently isolated for at least 3 years, multiple introductions of IAVs-S with phylogenetically distinct hemagglutinin (HA) genes occurred. In addition, at one farm, IAVs-S derived from a single introduction persisted for at least 3 years and carried no mutations at the deduced antigenic sites of the hemagglutinin protein, except for one at the antigenic site (Sa). Our results extend our understanding regarding the status of IAVs-S currently circulating in Japan and how they genetically evolve at the farm level.

IMPORTANCE Understanding the current status of influenza A viruses of swine (IAVs-S) and their evolution at the farm level is important for controlling these pathogens. Efforts to monitor IAVs-S during 2015 to 2019 yielded H1N1, H1N2, and H3N2 viruses. H1 genes in Japanese swine formed a unique clade in the classical swine H1 lineage of 1A.1, and H3 genes originating from 1999–2000 human seasonal influenza viruses appear to have become established among Japanese swine. A(H1N1)pdm09-derived H1 genes became introduced repeatedly and reassorted with endemic IAVs-S, resulting in various combinations of surface and internal genes among pig populations in Japan. At the farm level, multiple introductions of IAVs-S with phylogenetically distinct HA sequences occurred, or IAVs-S derived from a single introduction have persisted for at least 3 years with only a single mutation at the antigenic site of the HA protein. Continued monitoring of IAVs-S is necessary to update and maximize control strategies.

Yaravirus: A novel 80-nm virus infecting Acanthamoeba castellanii

Here we report the discovery of Yaravirus, a lineage of amoebal virus with a puzzling origin and evolution. Yaravirus presents 80-nm-sized particles and a 44,924-bp dsDNA genome encoding for 74 predicted proteins. Yaravirus genome annotation showed that none of its genes matched with sequences of known organisms at the nucleotide level; at the amino acid level, six predicted proteins had distant matches in the nr database. Complimentary prediction of three-dimensional structures indicated possible function of 17 proteins in total. Furthermore, we were not able to retrieve viral genomes closely related to Yaravirus in 8,535 publicly available metagenomes spanning diverse habitats around the globe. The Yaravirus genome also contained six types of tRNAs that did not match commonly used codons. Proteomics revealed that Yaravirus particles contain 26 viral proteins, one of which potentially representing a divergent major capsid protein (MCP) with a predicted double jelly-roll domain. Structure-guided phylogeny of MCP suggests that Yaravirus groups together with the MCPs of Pleurochrysis endemic viruses. Yaravirus expands our knowledge of the diversity of DNA viruses. The phylogenetic distance between Yaravirus and all other viruses highlights our still preliminary assessment of the genomic diversity of eukaryotic viruses, reinforcing the need for the isolation of new viruses of protists.

Diversity of MSDIN family members in amanitin-producing mushrooms and the phylogeny of the MSDIN and prolyl oligopeptidase genes

Background

Amanitin-producing mushrooms, mainly distributed in the genera Amanita, Galerina and Lepiota, possess MSDIN gene family for the biosynthesis of many cyclopeptides catalysed by prolyl oligopeptidase (POP). Recently, transcriptome sequencing has proven to be an efficient way to mine MSDIN and POP genes in these lethal mushrooms. Thus far, only A. palloides and A. bisporigera from North America and A. exitialis and A. rimosa from Asia have been studied based on transcriptome analysis. However, the MSDIN and POP genes of many amanitin-producing mushrooms in China remain unstudied; hence, the transcriptomes of these speices deserve to be analysed.

Results

In this study, the MSDIN and POP genes from ten Amanita species, two Galerina species and Lepiota venenata were studied and the phylogenetic relationships of their MSDIN and POP genes were analysed. Through transcriptome sequencing and PCR cloning, 19 POP genes and 151 MSDIN genes predicted to encode 98 non-duplicated cyclopeptides, including α-amanitin, β-amanitin, phallacidin, phalloidin and 94 unknown peptides, were found in these species. Phylogenetic analysis showed that (1) MSDIN genes generally clustered depending on the taxonomy of the genus, while Amanita MSDIN genes clustered depending on the chemical substance; and (2) the POPA genes of Amanita, Galerina and Lepiota clustered and were separated into three different groups, but the POPB genes of the three distinct genera were clustered in a highly supported monophyletic group.

Conclusions

These results indicate that lethal Amanita species have the genetic capacity to produce numerous cyclopeptides, most of which are unknown, while lethal Galerina and Lepiota species seem to only have the genetic capacity to produce α-amanitin. Additionally, the POPB phylogeny of Amanita, Galerina and Lepiota conflicts with the taxonomic status of the three genera, suggesting that underlying horizontal gene transfer has occurred among these three genera.

ZapE/Afg1 interacts with Oxa1 and its depletion causes a multifaceted phenotype

ZapE/Afg1 is a component of the inner cell membrane of some eubacteria and the inner mitochondrial membrane of eukaryotes. This protein is involved in FtsZ-dependent division of eubacteria. In the yeast and human mitochondrion, ZapE/Afg1 likely interacts with Oxa1 and facilitates the degradation of mitochondrion-encoded subunits of respiratory complexes. Furthermore, the depletion of ZapE increases resistance to apoptosis, decreases oxidative stress tolerance, and impacts mitochondrial protein homeostasis. It remains unclear whether ZapE is a multifunctional protein, or whether some of the described effects are just secondary phenotypes. Here, we have analyzed the functions of ZapE in Trypanosoma brucei, a parasitic protist, and an important model organism. Using a newly developed proximity-dependent biotinylation approach (BioID2), we have identified the inner mitochondrial membrane insertase Oxa1 among three putative interacting partners of ZapE, which is present in two paralogs. RNAi-mediated depletion of both ZapE paralogs likely affected the function of respiratory complexes I and IV. Consistently, we show that the distribution of mitochondrial ZapE is restricted only to organisms with Oxa1, respiratory complexes, and a mitochondrial genome. We propose that the evolutionarily conserved interaction of ZapE with Oxa1, which is required for proper insertion of many inner mitochondrial membrane proteins, is behind the multifaceted phenotype caused by the ablation of ZapE.

Phylogenetic analysis of SARS-CoV-2 genomes in Turkey

COVID-19 has effectively spread worldwide. As of May 2020, Turkey is among the top ten countries with the most cases. A comprehensive genomic characterization of the virus isolates in Turkey is yet to be carried out. Here, we built a phylogenetic tree with globally obtained 15,277 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes. We identified the subtypes based on the phylogenetic clustering in comparison with the previously annotated classifications. We performed a phylogenetic analysis of the first 30 SARS-CoV-2 genomes isolated and sequenced in Turkey. We suggest that the first introduction of the virus to the country is earlier than the first reported case of infection. Virus genomes isolated from Turkey are dispersed among most types in the phylogenetic tree. We find 2 of the seventeen subclusters enriched with the isolates of Turkey, which likely have spread expansively in the country. Finally, we traced virus genomes based on their phylogenetic placements. This analysis suggested multiple independent international introductions of the virus and revealed a hub for the inland transmission. We released a web application to track the global and interprovincial virus spread of the isolates from Turkey in comparison to thousands of genomes worldwide.

A new species group from the Daphnia curvirostris species complex (Cladocera: Anomopoda) from the eastern Palaearctic: taxonomy, phylogeny and phylogeography

The eastern Palaearctic is a centre of diversity for freshwater cladocerans (Crustacea), but little is known about the evolution and taxonomy of this diversity. Daphnia curvirostris is a Holarctic species complex that has most of its diversity in the eastern Palaearctic. We examined the phylogeography, rates of evolution and taxonomic status for each clade of the D. curvirostris complex using morphological and genetic evidence from four genes. The cybertaxonomical and morphological evidence supported an eastern Palaearctic clade, with at least four species (described here as the Daphnia korovchinskyi sp. nov. group) having diagnostic morphological characters. We also detected convergent morphological characters in the D. curvirostris complex that provided information about species boundaries. Two of the new species (Daphnia koreana sp. nov. and Daphnia ishidai sp. nov.) are known from single ponds and are threatened by human activity. Divergence time estimates suggested an ancient origin (12–28 Mya) for the D. korovchinskyi group, but these estimates are complicated by the small number of calibration points.

The search of sequence variants using a constrained protein evolution simulation approach

Protein engineering or candidate therapeutic peptide optimization are processes in which the identification of relevant sequence variants is critical. Starting from one amino-acid sequence, the choice of the substitutions must meet the objective of not disrupting the structure of the protein, not impacting the main functional properties of the starting entity, while also meeting the condition to enhance some expected property such as thermal stability, resistance to degradation, … Here, we introduce a new approach of sequence evolution that focuses on the objective of not disrupting the structure of the initial protein by embedding a point to point control on the preservation of the local structure at each position in the sequence. For 6 mini-proteins, we find that, starting from a single sequence, our simple approach intrinsically contains information about site-specific rate heterogeneity of substitution, and that it is able to reproduce sequence diversity as can be observed in the sequences available in the Uniref repository. We show that our approach is able to provide information about positions not to substitute and about substitutions not to perform at a given position to maintain structure integrity. Overall, our results demonstrate that point to point preservation of the local structure along a sequence is an important determinant of sequence evolution.

Bacterial Secondary Metabolite Biosynthetic Potential in Soil Varies with Phylum, Depth, and Vegetation Type

Bacteria isolated from soils are major sources of specialized metabolites, including antibiotics and other compounds with clinical value that likely shape interactions among microbial community members and impact biogeochemical cycles. Yet, isolated lineages represent a small fraction of all soil bacterial diversity. It remains unclear how the production of specialized metabolites varies across the phylogenetic diversity of bacterial species in soils and whether the genetic potential for production of these metabolites differs with soil depth and vegetation type within a geographic region. We sampled soils and saprolite from three sites in a northern California Critical Zone Observatory with various vegetation and bedrock characteristics and reconstructed 1,334 metagenome-assembled genomes containing diverse biosynthetic gene clusters (BGCs) for secondary metabolite production. We obtained genomes for prolific producers of secondary metabolites, including novel groups within the Actinobacteria, Chloroflexi, and candidate phylum "Candidatus Dormibacteraeota." Surprisingly, one genome of a candidate phyla radiation (CPR) bacterium coded for a ribosomally synthesized linear azole/azoline-containing peptide, a capacity we found in other publicly available CPR bacterial genomes. Overall, bacteria with higher biosynthetic potential were enriched in shallow soils and grassland soils, with patterns of abundance of BGC type varying by taxonomy.IMPORTANCE Microbes produce specialized compounds to compete or communicate with one another and their environment. Some of these compounds, such as antibiotics, are also useful in medicine and biotechnology. Historically, most antibiotics have come from soil bacteria which can be isolated and grown in the lab. Though the vast majority of soil bacteria cannot be isolated, we can extract their genetic information and search it for genes which produce these specialized compounds. These understudied soil bacteria offer a wealth of potential for the discovery of new and important microbial products. Here, we identified the ability to produce these specialized compounds in diverse and novel bacteria in a range of soil environments. This information will be useful to other researchers who wish to isolate certain products. Beyond their use to humans, understanding the distribution and function of microbial products is key to understanding microbial communities and their effects on biogeochemical cycles.

Conserved motifs in the CRISPR leader sequence control spacer acquisition levels in Type I-D CRISPR-Cas systems

Integrating short DNA fragments at the correct leader-repeat junction is key to successful CRISPR-Cas memory formation. The Cas1-2 proteins are responsible to carry out this process. However, the CRISPR adaptation process additionally requires a DNA element adjacent to the CRISPR array, called leader, to facilitate efficient localization of the correct integration site. In this work, we introduced the core CRISPR adaptation genes cas1 and cas2 from the Type I-D CRISPR-Cas system of Synechocystis sp. 6803 into Escherichia coli and assessed spacer integration efficiency. Truncation of the leader resulted in a significant reduction of spacer acquisition levels and revealed the importance of different conserved regions for CRISPR adaptation rates. We found three conserved sequence motifs in the leader of I-D CRISPR arrays that each affected spacer acquisition rates, including an integrase anchoring site. Our findings support the model in which the leader sequence is an integral part of type I-D adaptation in Synechocystis sp. acting as a localization signal for the adaptation complex to drive CRISPR adaptation at the first repeat of the CRISPR array.

Comparative analysis of mitochondrial genomes among the subfamily Sarcophaginae (Diptera: Sarcophagidae) and phylogenetic implications

The subfamily Sarcophaginae is extremely diverse in morphology, habit and geographical distribution, and usually considered to be of significant ecological, medical, and forensic significance. In the present study, 18 mitochondrial genomes (mitogenomes) of sarcophagid flies were first obtained. The rearrangement and orientation of genes were identical with that of ancestral insects. The degrees of compositional heterogeneity in the datasets were extremely low. Furthermore, 13 protein-coding genes were evolving under purifying selection. The phylogenic relationship of the genus-group taxa Boettcheria + (Sarcophaga + (Peckia + (Ravinia + Oxysarcodexia))) was strongly supported. Four subgenera were recovered as monophyletic (Liopygia, Liosarcophaga, Pierretia, Heteronychia) in addition to Parasarcophaga as polyphyletic. The sister-relationships between S. dux and S. aegyptiaca, S. pingi and S. kawayuensis were recovered, respectively. Moreover, the molecular phylogenetic relationships among the subgenera Helicophagella, Kozlovea, Kramerea, Pandelleisca, Phallocheira, Pseudothyrsocnema, Sinonipponia and Seniorwhitea were rarely put forward prior to this study. This study provides insight into the population genetics, molecular biology, and phylogeny for the subfamily Sarcophaginae, especially for the subgeneric classification of Sarcophaga. However, compared with the enormous species diversity of flesh flies, the available mitogenomes are still limited for recovering the phylogeny of Sarcophaginae.

A Snapshot of SARS-CoV-2 Genome Availability up to April 2020 and its Implications: Data Analysis

BACKGROUND

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been growing exponentially, affecting over 4 million people and causing enormous distress to economies and societies worldwide. A plethora of analyses based on viral sequences has already been published both in scientific journals and through non-peer-reviewed channels to investigate the genetic heterogeneity and spatiotemporal dissemination of SARS-CoV-2. However, a systematic investigation of phylogenetic information and sampling bias in the available data is lacking. Although the number of available genome sequences of SARS-CoV-2 is growing daily and the sequences show increasing phylogenetic information, country-specific data still present severe limitations and should be interpreted with caution.

OBJECTIVE

The objective of this study was to determine the quality of the currently available SARS-CoV-2 full genome data in terms of sampling bias as well as phylogenetic and temporal signals to inform and guide the scientific community.

METHODS

We used maximum likelihood-based methods to assess the presence of sufficient information for robust phylogenetic and phylogeographic studies in several SARS-CoV-2 sequence alignments assembled from GISAID (Global Initiative on Sharing All Influenza Data) data released between March and April 2020.

RESULTS

Although the number of high-quality full genomes is growing daily, and sequence data released in April 2020 contain sufficient phylogenetic information to allow reliable inference of phylogenetic relationships, country-specific SARS-CoV-2 data sets still present severe limitations.

CONCLUSIONS

At the present time, studies assessing within-country spread or transmission clusters should be considered preliminary or hypothesis-generating at best. Hence, current reports should be interpreted with caution, and concerted efforts should continue to increase the number and quality of sequences required for robust tracing of the epidemic.

Identification and Mobilization of a Cryptic Antibiotic Biosynthesis Gene Locus from a Human-Pathogenic Nocardia Isolate

The genus Nocardia contains >50 human pathogenic species that cause a range of illnesses from skin and soft tissue infections to lung and brain infections. However, despite their membership in the most prominent family of secondary metabolite producers (the Actinomycetes), the ability of Nocardia species, especially those that cause human infections, to produce secondary metabolites has not been as well studied. Using genome mining, we have investigated cryptic secondary metabolite biosynthesis gene clusters from Nocardia species and identified a conserved locus within human pathogenic strains of Nocardia brasiliensis and Nocardia vulneris. Direct capture and heterologous expression in a Streptomyces host activated the biosynthetic locus, revealing it to be the source of the brasiliquinones, benz[a]anthraquinone antibiotics whose biosynthetic pathway has remained hidden for over two decades, until now. Our findings highlight these hitherto neglected human pathogenic Nocardia as a source of diverse and important natural products.

Under-the-Radar Dengue Virus Infections in Natural Populations of Aedes aegypti Mosquitoes

The incidence of locally acquired dengue infections increased during the last decade in the United States, compelling a sustained research effort concerning the dengue mosquito vector, Aedes aegypti, and its microbiome, which has been shown to influence virus transmission success. We examined the "metavirome" of four populations of Aedes aegypti mosquitoes collected in 2016 to 2017 in Manatee County, FL. Unexpectedly, we discovered that dengue virus serotype 4 (DENV4) was circulating in these mosquito populations, representing the first documented case of such a phenomenon in the absence of a local DENV4 human case in this county over a 2-year period. We confirmed that all of the mosquito populations carried the same DENV4 strain, assembled its full genome, validated infection orthogonally by reverse transcriptase PCR, traced the virus origin, estimated the time period of its introduction to the Caribbean region, and explored the viral genetic signatures and mosquito-specific virome associations that potentially mediated DENV4 persistence in mosquitoes. We discuss the significance of prolonged maintenance of the DENV4 infections in A. aegypti that occurred in the absence of a DENV4 human index case in Manatee County with respect to the inability of current surveillance paradigms to detect mosquito vector infections prior to a potential local outbreak.IMPORTANCE Since 1999, dengue outbreaks in the continental United States involving local transmission have occurred only episodically and only in Florida and Texas. In Florida, these episodes appear to be coincident with increased introductions of dengue virus into the region through human travel and migration from countries where the disease is endemic. To date, the U.S. public health response to dengue outbreaks has been largely reactive, and implementation of comprehensive arbovirus surveillance in advance of predictable transmission seasons, which would enable proactive preventative efforts, remains unsupported. The significance of our finding is that it is the first documented report of DENV4 transmission to and maintenance within a local mosquito vector population in the continental United States in the absence of a human case during two consecutive years. Our data suggest that molecular surveillance of mosquito populations in high-risk, high-tourism areas of the United States may enable proactive, targeted vector control before potential arbovirus outbreaks.

Genomic and phylogenetic analysis of two guinea pig adenovirus strains recovered from archival lung tissue

Next generation sequencing was used to determine the whole genome sequence for two different strains of guinea pig adenovirus (GPAdV) detected in association with outbreaks of pneumonia in Australia in 1996, and in Germany in 1997 using total DNA extracted from infected archival frozen lung tissue as a template. The length of the determined genomic sequences was 37,031 bp and 37,070 bp, respectively. The nucleotide composition showed a relatively high content of guanine + cytosine (G + C) of 62 %. The 99.6 % nucleotide identity between the two sequenced viruses suggests that they may represent variants of the same genotype. The GPAdV genome exhibits the genomic features of a typical mastadenovirus with at least 32 open reading frames identified. Five novel open reading frames were found at the right end of the genomic sequence. One of them maps to the predicted E3 region and encodes a putative CR1 protein, two map to the E4 region, and two map to the l strand of L1 and L3, respectively. Our phylogenetic analysis of whole genome sequences showed that among the mammalian AdV species described to date, GPAdV is most closely related to MAdV-2 The characterization of this mastadenovirus species offers an opportunity to develop a new small animal model to study mammalian adenovirus pathogenesis.

Substrate recognition and ATPase activity of the E. coli cysteine/cystine ABC transporter YecSC-FliY

Sulfur is essential for biological processes such as amino acid biogenesis, iron-sulfur cluster formation, and redox homeostasis. To acquire sulfur-containing compounds from the environment, bacteria have evolved high-affinity uptake systems, predominant among which is the ABC transporter family. Theses membrane-embedded enzymes use the energy of ATP hydrolysis for transmembrane transport of a wide range of biomolecules against concentration gradients. Three distinct bacterial ABC import systems of sulfur-containing compounds have been identified, but the molecular details of their transport mechanism remain poorly characterized. Here we provide results from a biochemical analysis of the purified Escherichia coli YecSC-FliY cysteine/cystine import system. We found that the substrate-binding protein FliY binds l-cystine, l-cysteine, and d-cysteine with micromolar affinities. However, binding of the l- and d-enantiomers induced different conformational changes of FliY, where the l- enantiomer-substrate-binding protein complex interacted more efficiently with the YecSC transporter. YecSC had low basal ATPase activity that was moderately stimulated by apo FliY, more strongly by d-cysteine-bound FliY, and maximally by l-cysteine- or l-cystine-bound FliY. However, at high FliY concentrations, YecSC reached maximal ATPase rates independent of the presence or nature of the substrate. These results suggest that FliY exists in a conformational equilibrium between an open, unliganded form that does not bind to the YecSC transporter and closed, unliganded and closed, liganded forms that bind this transporter with variable affinities but equally stimulate its ATPase activity. These findings differ from previous observations for similar ABC transporters, highlighting the extent of mechanistic diversity in this large protein family.

Early eukaryotic origins and metazoan elaboration of MAPR family proteins

The membrane-associated progesterone receptor (MAPR) family consists of heme-binding proteins containing a cytochrome b₅ (cytb₅) domain characterized by the presence of a MAPR-specific interhelical insert region (MIHIR) between helices 3 and 4 of the canonical cytb5-domain fold. Animals possess three MAPR genes (PGRMC-like, Neuferricin and Neudesin). Here we show that all three animal MAPR genes were already present in the common ancestor of the opisthokonts (comprising animals and fungi as well as related single-celled taxa). All three MAPR genes acquired extensions C-terminal to the cytb₅ domain, either before or with the evolution of animals. The archetypical MAPR protein, progesterone receptor membrane component 1 (PGRMC1), contains phosphorylated tyrosines Y139 and Y180. The combination of Y139/Y180 appeared in the common ancestor of cnidarians and bilaterians, along with an early embryological organizer and synapsed neurons, and is strongly conserved in all bilaterian animals. A predicted protein interaction motif in the PGRMC1 MIHIR is potentially regulated by Y139 phosphorylation. A multilayered model of animal MAPR function acquisition includes some pre-metazoan functions (e.g., heme binding and cytochrome P450 interactions) and some acquired animal-specific functions that involve regulation of strongly conserved protein interaction motifs acquired by animals (Metazoa). This study provides a conceptual framework for future studies, against which especially PGRMC1's multiple functions can perhaps be stratified and functionally dissected.

New morphological and molecular data for Xystretrum solidum (Gorgoderidae, Gorgoderinae) from Sphoeroides testudineus (Tetraodontiformes, Tetraodontidae) in Mexican waters

Adults of trematodes in the genus Xystretrum Linton, 1910 (Gorgoderidae, Gorgoderinae) are parasites found exclusively in the urinary bladders of tetraodontiform fishes. However, limited and unclear morphological data were used to describe the type species, X. solidum Linton, 1910. Here, we present the first detailed morphological information for a member of Xystretrum. Morphological characters were described using light and scanning electron microscopy (SEM) of Xystretrum specimens from Sphoeroides testudineus (Linnaeus) (Tetraodontiformes, Tetraodontidae), collected at six localities off the northern Yucatan Peninsula coast of the Gulf of Mexico. We also compared sequence fragments of the 28S (region D1–D3) ribosomal DNA and mitochondrial Cytochrome c oxidase subunit 1 (COI) gene with those available for other gorgoderine taxa. We assigned these Xystretrum specimens to X. solidum, despite the incompleteness of published descriptions. The data provide a foundation for future work to validate the identities of X. solidum, X. papillosum Linton, 1910 and X. pulchrum (Travassos, 1920) with new collections from the type localities and hosts. Comparisons of 28S and COI regions described here also provide an opportunity to evaluate the monophyletic status of Xystretrum.

Pseudobactrodesmium (Dactylosporaceae, Eurotiomycetes, Fungi) a Novel Lignicolous Genus

During our ongoing surveys of fungi on submerged wood in the Greater Mekong Subregion, we collected two new species similar to Bactrodesmium longisporum. Pseudobactrodesmium gen. nov. is introduced to accommodate the new species, P. aquaticum, P. chiangmaiensis and B. longisporum is transferred to this genus. Fasciculate conidiophores, enteroblastic conidiogenous cells and subulate to fusiform, phragmoseptate conidia with a tapering apical cell and sheath characterize the genus. Pseudobactrodesmium aquaticum has longer conidia than P. chiangmaiensis. The placement of Pseudobactrodesmium in Dactylosporaceae (Eurotiomycetes) is a novel finding based on analyses of combined LSU, SSU, ITS and RPB2 sequence data. Our study reveals that Pseudobactrodesmium is likely to be a speciose genus with different species in streams around the world.

Evolution of endemic and sylvatic lineages of dengue virus

Recent disease outbreaks have raised awareness of tropical pathogens, especially mosquito-borne viruses. Dengue virus (DENV) is a widely studied mammalian pathogen transmitted by various species of mosquito in the genus Aedes, especially Aedes aegypti and Aedes albopictus. The prevailing view of the research community is that endemic viral lineages that cause epidemics of DENV in humans have emerged over time from sylvatic viral lineages, which persist in wild, non-human primates. These notions have been examined by researchers through phylogenetic analyses of the envelope gene (E) from the four serotypes of DENV (serotypes DENV-1 to DENV-4). In these previous reports, researchers used visual inspection of a phylogeny in order to assert that sylvatic lineages lead to endemic clades. In making this assertion, these researchers also reasserted the model of periodic sylvatic to endemic disease outbreaks. Since that study, there has been a significant increase in data both in terms of metadata (e.g., place and host of isolation) and genetic sequences of DENV. Here, we re-examine the model of sylvatic to endemic shifts in viral lineages through a phylogenetic tree search and character evolution study of metadata on the tree. We built a dataset of nucleotide sequences for 188 isolates of DENV that have metadata on sylvatic or endemic sampling along with three orthologous sequences from West Nile virus as the outgroup for the phylogenetic analysis. In contrast to previous research, we find that there are several shifts from endemic to sylvatic lineages as well as sylvatic to endemic lineages, indicating a much more dynamic model of evolution. We propose that a model that allows oscillation between sylvatic and endemic hosts better captures the dynamics of DENV transmission.

CER16 Inhibits Post-Transcriptional Gene Silencing of CER3 to Regulate Alkane Biosynthesis

The aerial surfaces of land plants have a protective layer of cuticular wax. Alkanes are common components of these waxes, and their abundance is affected by a range of stresses. The CER16 protein has been implicated in alkane biosynthesis in the cuticular wax of Arabidopsis (Arabidopsis thaliana). Here, we identified two new mutant alleles of CER16 in Arabidopsis resulting in production of less wax with dramatically fewer alkanes than the wild type. Map-based cloning with genetic analysis revealed that the cer16 phenotype was caused by complete loss of AT5G44150, encoding a protein with no known domains or motifs. Comparative transcriptomic analysis revealed that transcripts of CER3, previously shown to play a principal role in alkane production, were markedly reduced in the cer16 mutants. To define the relationship between CER3 and CER16, we transformed the full CER3 gene into a cer16 mutant. Transgenic CER3 expression was silenced, and levels of small interfering RNAs targeting CER3 were significantly increased. Mutating two major components of the RNA-silencing machinery in a cer16 genetic background restored CER3 transcript levels to wild-type levels, with the stems restored to wild-type glaucousness. We suggest that CER16 deficiency induces post-transcriptional gene silencing of both endogenous and exogenous expression of CER3.