Estimating maximum likelihood phylogenies with PhyML

Investigation of Genetic Relationships Between Hanseniaspora Species Found in Grape Musts Revealed Interspecific Hybrids With Dynamic Genome Structures

Hanseniaspora, a predominant yeast genus of grape musts, includes sister species recently reported as fast evolving. The aim of this study was to investigate the genetic relationships between the four most closely related species, at the population level. A multi-locus sequence typing strategy based on five markers was applied on 107 strains, confirming the clear delineation of species H. uvarum, H. opuntiae, H. guilliermondii, and H. pseudoguilliermondii. Huge variations were observed in the level of intraspecific nucleotide diversity, and differences in heterozygosity between species indicate different life styles. No clear population structure was detected based on geographical or substrate origins. Instead, H. guilliermondii strains clustered into two distinct groups, which may reflect a recent step toward speciation. Interspecific hybrids were detected between H. opuntiae and H. pseudoguilliermondii. Their characterization using flow cytometry, karyotypes and genome sequencing showed different genome structures in different ploidy contexts: allodiploids, allotriploids, and allotetraploids. Subculturing of an allotriploid strain revealed chromosome loss equivalent to one chromosome set, followed by an auto-diploidization event, whereas another auto-diploidized tetraploid showed a segmental duplication. Altogether, these results suggest that Hanseniaspora genomes are not only fast evolving but also highly dynamic.

Longitudinal within-host evolution of HIV Nef-mediated CD4, HLA and SERINC5 downregulation activity: a case study

The HIV accessory protein Nef downregulates the viral entry receptor CD4, the Human Leukocyte Antigen (HLA)-A and -B molecules, the Serine incorporator 5 (SERINC5) protein and other molecules from the infected cell surface, thereby promoting viral infectivity, replication and immune evasion. The nef locus also represents one of the most genetically variable regions in the HIV genome, and nef sequences undergo substantial evolution within a single individual over the course of infection. Few studies however have simultaneously characterized the impact of within-host nef sequence evolution on Nef protein function over prolonged timescales. Here, we isolated 50 unique Nef clones by single-genome amplification over an 11-year period from the plasma of an individual who was largely naïve to antiretroviral treatment during this time. Together, these clones harbored nonsynonymous substitutions at 13% of nef’s codons. We assessed their ability to downregulate cell-surface CD4, HLA and SERINC5 and observed that all three Nef functions declined modestly over time, where the reductions in CD4 and HLA downregulation (an average of 0.6% and 2.0% per year, respectively) achieved statistical significance. The results from this case study support all three Nef activities as being important to maintain throughout untreated HIV infection, but nevertheless suggest that, despite nef’s mutational plasticity, within-host viral evolution can compromise Nef function, albeit modestly, over prolonged periods.

Novel RNA viruses associated with Plasmodium vivax in human malaria and Leucocytozoon parasites in avian disease

Eukaryotes of the genus Plasmodium cause malaria, a parasitic disease responsible for substantial morbidity and mortality in humans. Yet, the nature and abundance of any viruses carried by these divergent eukaryotic parasites is unknown. We investigated the Plasmodium virome by performing a meta-transcriptomic analysis of blood samples taken from patients suffering from malaria and infected with P. vivax, P. falciparum or P. knowlesi. This resulted in the identification of a narnavirus-like sequence, encoding an RNA polymerase and restricted to P. vivax samples, as well as an associated viral segment of unknown function. These data, confirmed by PCR, are indicative of a novel RNA virus that we term Matryoshka RNA virus 1 (MaRNAV-1) to reflect its analogy to a "Russian doll": a virus, infecting a parasite, infecting an animal. Additional screening revealed that MaRNAV-1 was abundant in geographically diverse P. vivax derived from humans and mosquitoes, strongly supporting its association with this parasite, and not in any of the other Plasmodium samples analyzed here nor Anopheles mosquitoes in the absence of Plasmodium. Notably, related bi-segmented narnavirus-like sequences (MaRNAV-2) were retrieved from Australian birds infected with a Leucocytozoon—a genus of eukaryotic parasites that group with Plasmodium in the Apicomplexa subclass hematozoa. Together, these data support the establishment of two new phylogenetically divergent and genomically distinct viral species associated with protists, including the first virus likely infecting Plasmodium parasites. As well as broadening our understanding of the diversity and evolutionary history of the eukaryotic virosphere, the restriction to P. vivax may be of importance in understanding P. vivax-specific biology in humans and mosquitoes, and how viral co-infection might alter host responses at each stage of the P. vivax life-cycle.

While parasites are a major cause of human disease, they can themselves be infected by viruses. We asked whether three of the major malaria-causing parasites in humans—Plasmodium vivax, P. falciparum and P. knowlesi—were also infected by viruses. To this end we performed total RNA-Sequencing (“meta-transcriptomics”) on human blood samples infected with these Plasmodium species. This resulted in the discovery of an abundant bi-segmented virus—Matryoshka RNA virus 1 (MaRNAV-1)—in all P. vivax samples tested (but no other Plasmodium species) that contains a replicase segment related to those of narnaviruses, arguably the simplest type of RNA viruses discovered to date. By screening for MaRNAV-1 in a larger set of Plasmodium species we revealed a strong specificity between this virus and P. vivax, as well as the presence of a related virus—MaRNAV-2—in avian Leucocytozoon hematozoa parasites. This is the first discovery of a Plasmodium-associated virus and will assist in revealing the deep evolutionary history of RNA viruses and our understanding of Plasmodium biology and disease processes.

Insight into diversity of bacteria belonging to the order Rickettsiales in 9 arthropods species collected in Serbia

Rickettsiales bacteria in arthropods play a significant role in both public health and arthropod ecology. However, the extensive genetic diversity of Rickettsiales endosymbionts of arthropods is still to be discovered. In 2016, 515 arthropods belonging to 9 species of four classes (Insecta, Chilopoda, Diplopoda and Arachnida) were collected in Serbia. The presence and genetic diversity of Rickettsiales bacteria were evaluated by characterizing the 16S rRNA (rrs), citrate synthase (gltA) and heat shock protein (groEL) genes. The presence of various Rickettsiales bacteria was identified in the majority of tested arthropod species. The results revealed co-circulation of five recognized Rickettsiales species including Rickettsia, Ehrlichia and Wolbachia, as well as four tentative novel species, including one tentative novel genus named Neowolbachia. These results suggest the remarkable genetic diversity of Rickettsiales bacteria in certain arthropod species in this region. Furthermore, the high prevalence of spotted fever group Rickettsia in Ixodes ricinus ticks highlights the potential public health risk of human Rickettsia infection.

The dynamic proteome of influenza A virus infection identifies M segment splicing as a host range determinant

Pandemic influenza A virus (IAV) outbreaks occur when strains from animal reservoirs acquire the ability to infect and spread among humans. The molecular basis of this species barrier is incompletely understood. Here we combine metabolic pulse labeling and quantitative proteomics to monitor protein synthesis upon infection of human cells with a human- and a bird-adapted IAV strain and observe striking differences in viral protein synthesis. Most importantly, the matrix protein M1 is inefficiently produced by the bird-adapted strain. We show that impaired production of M1 from bird-adapted strains is caused by increased splicing of the M segment RNA to alternative isoforms. Strain-specific M segment splicing is controlled by the 3' splice site and functionally important for permissive infection. In silico and biochemical evidence shows that avian-adapted M segments have evolved different conserved RNA structure features than human-adapted sequences. Thus, we identify M segment RNA splicing as a viral host range determinant.

Bioprospecting desert plant Bacillus endophytic strains for their potential to enhance plant stress tolerance

Plant growth-promoting bacteria (PGPB) are known to increase plant tolerance to several abiotic stresses, specifically those from dry and salty environments. In this study, we examined the endophyte bacterial community of five plant species growing in the Thar desert of Pakistan. Among a total of 368 culturable isolates, 58 Bacillus strains were identified from which the 16 most divergent strains were characterized for salt and heat stress resilience as well as antimicrobial and plant growth-promoting (PGP) activities. When the 16 Bacillus strains were tested on the non-host plant Arabidopsis thaliana, B. cereus PK6-15, B. subtilis PK5-26 and B. circulans PK3-109 significantly enhanced plant growth under salt stress conditions, doubling fresh weight levels when compared to uninoculated plants. B. circulans PK3-15 and PK3-109 did not promote plant growth under normal conditions, but increased plant fresh weight by more than 50% when compared to uninoculated plants under salt stress conditions, suggesting that these salt tolerant Bacillus strains exhibit PGP traits only in the presence of salt. Our data indicate that the collection of 58 plant endophytic Bacillus strains represents an important genomic resource to decipher plant growth promotion at the molecular level.

The genetic basis of adaptive evolution in parasitic environment from the Angiostrongylus cantonensis genome

Angiostrongylus cantonensis (rat lungworm) is the etiological agent of angiostrongyliasis, mainly causing eosinophilic meningitis or meningoencephalitis in human. Although the biology of A. cantonensis is relatively well known, little is understood about the mechanisms of the parasite’s development and survival in definitive hosts, or its adaptation to a broad range of snail intermediate hosts. Here, we generate a high-quality assembly of a well-defined laboratory strain of A. cantonensis from Guangzhou, China, by using Illumina and PacBio sequencing technologies. We undertake comparative analyses with representative helminth genomes and explore transcriptomic data throughout key developmental life-cycles of the parasite. We find that part of retrotransposons and gene families undergo multiple waves of expansions. These include extracellular superoxide dismutase (EC-SOD) and astacin-like proteases which are considered to be associated with invasion and survival of the parasite. Furthermore, these paralogs from different sub-clades based on phylogeny, have different expression patterns in the molluscan and rodent stages, suggesting divergent functions under the different parasitic environment. We also find five candidate convergent signatures in the EC-SOD proteins from flukes and one sub-clade of A. cantonensis. Additionally, genes encoding proteolytic enzymes, involved in host hemoglobin digestion, exhibit expansion in A. cantonensis as well as two other blood-feeding nematodes. Overall, we find several potential adaptive evolutionary signatures in A. cantonensis, and also in some other helminths with similar traits. The genome and transcriptomes provide a useful resource for detailed studies of A. cantonensis-host adaptation and an in-depth understanding of the global-spread of angiostrongyliasis.

Angiostrongylus cantonensis, rat lungworm, is a common pathogen that causes human eosinophilic meningitis via eating contaminated food. Human angiostrongyliasis has been reported globally. This worm has a complex life-cycle, which includes an especially wide range of snails as intermediate hosts, making it more difficult to eradicate. In this study, we sequenced the genome and transcriptome, and performed comparative analyses to study the potential genetics of its biology using short-read and long-read sequencing technologies. We revealed some potential adaptive evolution in the genome, such as the expansion of retrotransposons and gene families encoding antioxidant enzymes, invasion, migration and digestion related proteases. Specifically, we found a potential clue suggesting convergent evolution of EC-SODs in Angiostrongylus and flukes, all of which require snails as intermediate hosts. These results provide an abundant data resource to study the biology and evolution of A. cantonensis and showed some potential targets against A. cantonensis and helminths with similar traits.

Helicobacter pylori in ancient human remains

The bacterium Helicobacter pylori (H. pylori) infects the stomachs of approximately 50% of all humans. With its universal occurrence, high infectivity and virulence properties it is considered as one of the most severe global burdens of modern humankind. It has accompanied humans for many thousands of years, and due to its high genetic variability and vertical transmission, its population genetics reflects the history of human migrations. However, especially complex demographic events such as the colonisation of Europe cannot be resolved with population genetic analysis of modern H. pylori strains alone. This is best exemplified with the reconstruction of the 5300-year-old H. pylori genome of the Iceman, a European Copper Age mummy. Our analysis provided precious insights into the ancestry and evolution of the pathogen and underlined the high complexity of ancient European population history. In this review we will provide an overview on the molecular analysis of H. pylori in mummified human remains that were done so far and we will outline methodological advancements in the field of ancient DNA research that support the reconstruction and authentication of ancient H. pylori genome sequences.

Gene Composition as a Potential Barrier to Large Recombinations in the Bacterial Pathogen Klebsiella pneumoniae

Klebsiella pneumoniae (Kp) is one of the most important nosocomial pathogens worldwide, able to cause multiorgan infections and hospital outbreaks. One of the most widely disseminated lineage of Kp is the clonal group 258 (CG258), which includes the highly resistant "high-risk" sequence types ST258 and ST11. Genomic investigations revealed that very large recombination events have occurred during the emergence of Kp lineages. A striking example is provided by ST258, which has undergone a recombination event that replaced over 1 Mb of the genome with DNA from an unrelated Kp donor. Although several examples of this phenomenon have been documented in Kp and other bacterial species, the significance of these very large recombination events for the emergence of either hypervirulent or resistant clones remains unclear. Here, we present an analysis of 834 Kp genomes that provides data on the frequency of these very large recombination events (defined as those involving >100 kb), their distribution within the genome, and the dynamics of gene flow within the Kp population. We note that very large recombination events occur frequently, and in multiple lineages, and that the majority of recombinational exchanges are clustered within two overlapping genomic regions, which have been involved by recombination events with different frequencies. Our results also indicate that certain lineages are more likely to act as donors to CG258. Furthermore, comparison of gene content in CG258 and non-CG258 strains agrees with this pattern, suggesting that the success of a large recombination depends on gene composition in the exchanged genomic portion.

Identification of Norovirus and Human Parechovirus in Patients With Hand, Foot and Mouth Disease Syndrome

BACKGROUND

Hand, foot and mouth disease (HFMD) is caused mostly by enteroviruses. However, other viral agents also can cause similar syndromes, and hence, the infections they cause are often misdiagnosed clinically. To determine non-enterovirus etiologic agents in HFMD-like cases, we screened enterovirus-negative samples collected from the patients who were clinically diagnosed as HFMD in China.

METHODS

Two hundred enterovirus-negative samples were collected previously in Wenzhou city of Zhejiang province, China. Both high throughput sequencing and RT-PCR were used to screen viral agents. In addition, their clinical features were analyzed.

RESULTS

Norovirus (NoV) and human parechovirus (HPeV) were identified from 22 (11.00%) and 9 (4.50%) samples, respectively. In addition, the complete genome sequences were recovered from 4 NoV-positive samples, and the VP1/3Dpol gene sequences were recovered from 5 HPeV-positive samples. Phylogenetic analyses of the NoV sequences revealed that they were closely related to those circulated in other regions of China. Notably, 4 genotypes of HPeVs, including HPeV-1, HPeV-4, HPeV-5 and HPeV-14, were found, indicating high genetic diversity of the virus. Frequent recombination between various genotypes was also observed in the HPeVs. Although most of the patients presented with the clinical features of HFMD, 4 patients infected with NoV GII.4 and 3 patients infected with HPeV-1 (1) and HPeV-4 (2) were characterized with diarrhea. Finally, tonsillitis, convulsion and granulocytopenia were observed in 1 NoV GII.4 patient, while liver dysfunction was found in 1 NoV GII.17 patient.

CONCLUSIONS

These data reveal the variety of agents in the cases clinically diagnosed as HFMD.

Emergence of Southern Rice Black-Streaked Dwarf Virus in the Centuries-Old Chinese Yuanyang Agrosystem of Rice Landraces

Southern rice black-streaked dwarf virus (SRBSDV), which causes severe disease symptoms in rice (Oriza sativa L.) has been emerging in the last decade throughout northern Vietnam, southern Japan and southern, central and eastern China. Here we attempt to quantify the prevalence of SRBSDV in the Honghe Hani rice terraces system (HHRTS)-a Chinese 1300-year-old traditional rice production system. We first confirm that genetically diverse rice varieties are still being cultivated in the HHRTS and categorize these varieties into three main genetic clusters, including the modern hybrid varieties group (MH), the Hongyang improved modern variety group (HY) and the traditional indica landraces group (TIL). We also show over a 2-year period that SRBSDV remains prevalent in the HHRTS (20.1% prevalence) and that both the TIL (17.9% prevalence) and the MH varieties (5.1% prevalence) were less affected by SRBSDV than were the HY varieties (30.2% prevalence). Collectively we suggest that SRBSDV isolates are freely moving within the HHRTS and that TIL, HY and MH rice genetic clusters are not being preferentially infected by particular SRBSDV lineages. Given that SRBSDV can cause 30-50% rice yield losses, our study emphasizes both the need to better monitor the disease in the HHRTS, and the need to start considering ways to reduce its burden on rice production.

Mode and tempo of human hepatitis virus evolution

Human viral hepatitis, a major cause of morbidity and mortality worldwide, is caused by highly diverse viruses with different genetic, ecological, and pathogenetic features. Technological advances that allow throughput sequencing of viral genomes, as well as the development of computational tools to analyze such genome data, have largely expanded our knowledge on the host range and evolutionary history of human hepatitis viruses. Thus, with the exclusion of hepatitis D virus, close or distant relatives of these human pathogens were identified in a number of domestic and wild mammals. Also, sequences of human viral strains isolated from different geographic locations and over different time-spans have allowed the application of phylogeographic and molecular dating approaches to large viral phylogenies. In this review, we summarize the most recent insights into our understanding of the evolutionary events and ecological contexts that determined the origin and spread of human hepatitis viruses.

A complex evolutionary relationship between HHV-6A and HHV-6B

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6, iciHHV6). We reconstructed the population structure of HHV-6 and we show that HHV-6A genomes diverged less than HHV-6B genomes from the ancestral common HHV-6A/B population. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, exogenous viral and iciHHV-6B populations were similar in terms of ancestry components, with no evident geographic structuring. Most HHV-6B genomes sampled to date derive from viral populations that experienced considerable drift. However, a population of HHV-6 exogenous viruses, currently classified as HHV-6B and sampled in New York state, formed a separate cluster (NY cluster) and harbored a considerable portion of HHV-6A-like ancestry. Recombination detection methods identified these viruses as interspecies recombinants, but phylogenetic reconstruction indicated that the recombination signals are due to shared ancestry. In analogy to iciHHV-6A, NY cluster viruses have high nucleotide diversity and constant population size. We propose that HHV-6A sequences and the NY cluster population diverged from an ancestral HHV-6A-like population. A relatively recent bottleneck of the NY (or a related) population with subsequent expansion originated most HHV-6B genomes currently sampled. Our findings indicate that the distinction between HHV-6A and -6B is not as clear-cut as previously thought. More generally, epidemiological and clinical surveys would benefit from taking HHV-6 genetic diversity into account.

Genome-wide survey and expression analyses of the GRAS gene family in Brassica napus reveals their roles in root development and stress response

Genome-wide identification, classification, expression analyses, and functional characterization of GRAS genes in oil crop, Brassica napus, indicate their importance in root development and stress response. GRAS proteins are a plant-specific transcription factor gene family involved in tissues development and stress response. We classified 87 putative GRAS genes in the Brassica napus genome (BnGRASs) into 13 subfamilies by phylogenetic analysis. The C-terminal GRAS domains of Brassica napus (B. napus) proteins were less conserved among subfamilies, but were conserved within each subfamily. A series of analyses revealed that 89.7% of the BnGRASs did not have intron insertions, and 24 specific-motifs were found at the N-terminal. A highly conserved microRNA 171 (miRNA171) target was observed specifically in the HAM subfamily across land plants. A total of 868 pairs of interaction proteins were predicted, the primary of which were transcription factors involved in transcriptional regulation and signal transduction. Integrated comparative analysis of GRAS genes across 26 species of algae, mosses, ferns, gymnosperms, and angiosperms revealed that this gene family originated in early mosses and was classified into 19 subfamilies, 14 of which may have originated prior to bryophyte evolution. RNA-Seq analysis demonstrated that most BnGRASs were widely expressed in different tissues/organs at different stages in B. napus, and 24 BnGRASs were highly/specifically expressed in roots. Results from a qRT-PCR analysis suggested that two BnGRASs belonging to SCR and LISCL subfamilies potentially have important roles in the stress response of roots.

Persistent Core Populations Shape the Microbiome Throughout the Water Column in the North Pacific Subtropical Gyre

Marine microbial communities are responsible for many important ecosystem processes in the oceans. Their variability across time and depths is well recognized, but mostly at a coarse-grained taxonomic resolution. To gain a deeper perspective on ecological patterns of bacterioplankton diversity in the North Pacific Subtropical Gyre, we characterized bacterioplankton communities throughout the water column at a fine-grained taxonomic level with a focus on temporally persistent (core) populations. Considerable intra-clade microdiversity was evident in virtually every microbial clade examined. While some of the most abundant populations comprised only a small fraction of the intra-clade microdiversity, they formed a temporally persistent core within a more diverse array of less abundant ephemeral populations. The depth-stratified population structure within many phylogenetically disparate clades suggested that ecotypic variation was the rule among most planktonic bacterial and archaeal lineages. Our results suggested that the abundant, persistent core populations comprised the bulk of the biomass within any given clade. As such, we postulate that these core populations are largely responsible for microbially driven ecosystem processes, and so represent ideal targets for elucidating key microbial processes in the open-ocean water column.

Plastome of Quercus xanthoclada and comparison of genomic diversity amongst selected Quercus species using genome skimming

The genus Quercus L. contains several of the most economically important species for timber production in the Northern Hemisphere. It was one of the first genera described, but genetic diversity at a global scale within and amongst oak species remains unclear, despite numerous regional or species-specific assessments. To evaluate global plastid diversity in oaks, we sequenced the complete chloroplast of Quercus xanthoclada and compared its sequence with those available from other main taxonomic groups in Quercus. We quantify genomic divergence amongst oaks and performed a sliding window analysis to detect the most variable regions amongst members of the various clades, as well as divergent regions occurring in specific pairs of species. We identified private and shared SNPs amongst oaks species and sections and stress the need for a large global assessment of genetic diversity in this economically and ecologically important genus.

Phylogenetic Methods Inconsistently Predict the Direction of HIV Transmission Among Heterosexual Pairs in the HPTN 052 Cohort

BACKGROUND

We evaluated use of phylogenetic methods to predict the direction of human immunodeficiency virus (HIV) transmission.

METHODS

For 33 pairs of HIV-infected patients (hereafter, "index patients") and their partners who acquired genetically linked HIV infection during the study, samples were collected from partners and index patients close to the time when the partner seroconverted (hereafter, "SC samples"); for 31 pairs, samples collected from the index patient at an earlier time point (hereafter, "early index samples") were also available. Phylogenies were inferred using env next-generation sequences (1 tree per pair/subtype). The direction of transmission (DoT) predicted from each tree was classified as correct or incorrect on the basis of which sequences (those from the index patient or the partner) were closest to the root. DoT was also assessed using maximum parsimony to infer ancestral node states for 100 bootstrap trees.

RESULTS

DoT was predicted correctly for both single-pair and subtype-specific trees in 22 pairs (67%) by using SC samples and in 23 pairs (74%) by using early index samples. DoT was predicted incorrectly for 4 pairs (15%) by using SC or early index samples. In the bootstrap analysis, DoT was predicted correctly for 18 pairs (55%) by using SC samples and for 24 pairs (73%) by using early index samples. DoT was predicted incorrectly for 7 pairs (21%) by using SC samples and for 4 pairs (13%) by using early index samples.

CONCLUSIONS

Phylogenetic methods based solely on the tree topology of HIV env sequences, particularly without consideration of phylogenetic uncertainty, may be insufficient for determining DoT.

Phylogenetic Analyses of Some Melanoleuca Species (Agaricales, Tricholomataceae) in Northern China, With Descriptions of Two New Species and the Identification of Seven Species as a First Record

Two new species (Melanoleuca galerina and M. subgrammopodia) and seven new recorded species from northern China are described here using morphological and molecular methods. Melanoleuca galerina is mainly characterized by its hygrophanous pileus, decurrent lamellae, fibrous stipe and spores with round warts. Key characteristics of M. subgrammopodia include its discolored pileus, fibrous stipe and urticiform cystidia. The divergence time of Melanoleuca fungi as well as the phylogenetic relationships within this genus were analyzed using DNA sequences of the internal transcribed spacer (ITS) and the nuclear large subunit rDNA (nrLSU) gene fragments. Analyses revealed that morphological identifications and phylogenetic relationships were consistent with the results of divergence time, thereby confirming that M. galerina and M. subgrammopodia are new species.

PLAZA 4.0: an integrative resource for functional, evolutionary and comparative plant genomics

PLAZA (https://bioinformatics.psb.ugent.be/plaza) is a plant-oriented online resource for comparative, evolutionary and functional genomics. The PLAZA platform consists of multiple independent instances focusing on different plant clades, while also providing access to a consistent set of reference species. Each PLAZA instance contains structural and functional gene annotations, gene family data and phylogenetic trees and detailed gene colinearity information. A user-friendly web interface makes the necessary tools and visualizations accessible, specific for each data type. Here we present PLAZA 4.0, the latest iteration of the PLAZA framework. This version consists of two new instances (Dicots 4.0 and Monocots 4.0) providing a large increase in newly available species, and offers access to updated and newly implemented tools and visualizations, helping users with the ever-increasing demands for complex and in-depth analyzes. The total number of species across both instances nearly doubles from 37 species in PLAZA 3.0 to 71 species in PLAZA 4.0, with a much broader coverage of crop species (e.g. wheat, palm oil) and species of evolutionary interest (e.g. spruce, Marchantia). The new PLAZA instances can also be accessed by a programming interface through a RESTful web service, thus allowing bioinformaticians to optimally leverage the power of the PLAZA platform.

Challenges in funding and developing genomic software: roots and remedies

The computer software used for genomic analysis has become a crucial component of the infrastructure for life sciences. However, genomic software is still typically developed in an ad hoc manner, with inadequate funding, and by academic researchers not trained in software development, at substantial costs to the research community. I examine the roots of the incongruity between the importance of and the degree of investment in genomic software, and I suggest several potential remedies for current problems. As genomics continues to grow, new strategies for funding and developing the software that powers the field will become increasingly essential.

Genome Sequencing of the Japanese Eel (Anguilla japonica) for Comparative Genomic Studies on tbx4 and a tbx4 Gene Cluster in Teleost Fishes

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

Genome of the African cassava whitefly Bemisia tabaci and distribution and genetic diversity of cassava-colonizing whiteflies in Africa

The whitefy Bemisia tabaci, a species complex consisting of many morphologically indistinguishable species divided into distinct clades, is one of the most globally important agricultural pests and plant virus vectors. Cassava-colonizing B. tabaci transmits viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). Half of all cassava plants in Africa are affected by these viral diseases, resulting in annual production losses of more than US$ 1 billion. Here we report the draft genome of the cassava whitefly B. tabaci Sub-Saharan Africa - East and Central Africa (SSA-ECA), the super-abundant population that has been associated with the rapid spread of viruses causing the pandemics of CMD and CBSD. The SSA-ECA genome assembled from Illumina short reads has a total size of 513.7 Mb and a scaffold N50 length of 497 kb, and contains 15,084 predicted protein-coding genes. Phylogenetic analysis suggests that SSA-ECA diverged from MEAM1 around 5.26 million years ago. A comprehensive genetic analysis of cassava-colonizing B. tabaci in Africa was also conducted, in which a total of 243 whitefly specimens were collected from 18 countries representing all major cassava-growing regions in the continent and genotyped using NextRAD sequencing. Population genomic analyses confirmed the existence of six major populations linked by gene flow and inferred the distribution patterns of these populations across the African continent. The genome of SSA-ECA and the genetic findings provide valuable resources and guidance to facilitate whitefly research and the development of strategies to control cassava viral diseases spread by whiteflies.

Divergence dating using mixed effects clock modelling: An application to HIV-1

The need to estimate divergence times in evolutionary histories in the presence of various sources of substitution rate variation has stimulated a rich development of relaxed molecular clock models. Viral evolutionary studies frequently adopt an uncorrelated clock model as a generic relaxed molecular clock process, but this may impose considerable estimation bias if discrete rate variation exists among clades or lineages. For HIV-1 group M, rate variation among subtypes has been shown to result in inconsistencies in time to the most recent common ancestor estimation. Although this calls into question the adequacy of available molecular dating methods, no solution to this problem has been offered so far. Here, we investigate the use of mixed effects molecular clock models, which combine both fixed and random effects in the evolutionary rate, to estimate divergence times. Using simulation, we demonstrate that this model outperforms existing molecular clock models in a Bayesian framework for estimating time-measured phylogenies in the presence of mixed sources of rate variation, while also maintaining good performance in simpler scenarios. By analysing a comprehensive HIV-1 group M complete genome data set we confirm considerable rate variation among subtypes that is not adequately modelled by uncorrelated relaxed clock models. The mixed effects clock model can accommodate this rate variation and results in a time to the most recent common ancestor of HIV-1 group M of 1920 (1915-25), which is only slightly earlier than the uncorrelated relaxed clock estimate for the same data set. The use of complete genome data appears to have a more profound impact than the molecular clock model because it reduces the credible intervals by 50 per cent relative to similar estimates based on short envelope gene sequences.

Evolutionary journey and characterisation of a novel pan-gene associated with beer strains of Saccharomyces cerevisiae

The sequencing of over a thousand Saccharomyces cerevisiae genomes revealed a complex pangenome. Over one third of the discovered genes are not present in the S. cerevisiae core genome but instead are often restricted to a subset of yeast isolates and thus may be important for adaptation to specific environmental niches. We refer to these genes as "pan-genes," being part of the pangenome but not the core genome. Here, we describe the evolutionary journey and characterisation of a novel pan-gene, originally named hypothetical (HYPO) open-reading frame. Phylogenetic analysis reveals that HYPO has been predominantly retained in S. cerevisiae strains associated with brewing but has been repeatedly lost in most other fungal species during evolution. There is also evidence that HYPO was horizontally transferred at least once, from S. cerevisiae to Saccharomyces paradoxus. The phylogenetic analysis of HYPO exemplifies the complexity and intricacy of evolutionary trajectories of genes within the S. cerevisiae pangenome. To examine possible functions for Hypo, we overexpressed a HYPO-GFP fusion protein in both S. cerevisiae and Saccharomyces pastorianus. The protein localised to the plasma membrane where it accumulated initially in distinct foci. Time-lapse fluorescent imaging revealed that when cells are grown in wort, Hypo-gfp fluorescence spreads throughout the membrane during cell growth. The overexpression of Hypo-gfp in S. cerevisiae or S. pastorianus strains did not significantly alter cell growth in medium-containing glucose, maltose, maltotriose, or wort at different concentrations.

Phylogenetic position of the Drosophila fima and dentissima lineages, and the status of the D. melanogaster species group

The subgenus Sophophora of Drosophila, which includes D. melanogaster, is an important model for the study of molecular evolution, comparative genomics, and evolutionary developmental biology. Numerous phylogenetic studies have examined species relationships in the well-known melanogaster, obscura, willistoni, and saltans species groups, as well as the relationships among these clades. In contrast, other species groups of Sophophora have been relatively neglected and have not been subjected to molecular phylogenetic analysis. Here, we focus on the endemic African Drosophila fima and dentissima lineages. We find that both these clades fall within the broadly defined melanogaster species group, but are otherwise distantly related to each other. The new phylogeny supports pervasive divergent and convergent evolution of male-specific grasping structures (sex combs). We discuss the implications of these results for defining the boundaries of the melanogaster species group, and weigh the relative merits of "splitting" and "lumping" approaches to the taxonomy of this key model system.

Poly-Saturated Dolichols from Filamentous Fungi Modulate Activity of Dolichol-Dependent Glycosyltransferase and Physical Properties of Membranes

Mono-saturated polyprenols (dolichols) have been found in almost all Eukaryotic cells, however, dolichols containing additional saturated bonds at the ω-end, have been identified in A. fumigatus and A. niger. Here we confirm using an LC-ESI-QTOF-MS analysis, that poly-saturated dolichols are abundant in other filamentous fungi, Trichoderma reesei, A. nidulans and Neurospora crassa, while the yeast Saccharomyces cerevisiae only contains the typical mono-saturated dolichols. We also show, using differential scanning calorimetry (DSC) and fluorescence anisotropy of 1,6-diphenyl-l,3,5-hexatriene (DPH) that the structure of dolichols modulates the properties of membranes and affects the functioning of dolichyl diphosphate mannose synthase (DPMS). The activity of this enzyme from T. reesei and S. cerevisiae was strongly affected by the structure of dolichols. Additionally, the structure of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) model membranes was more strongly disturbed by the poly-saturated dolichols from Trichoderma than by the mono-saturated dolichols from yeast. By comparing the lipidome of filamentous fungi with that from S. cerevisiae, we revealed significant differences in the PC/PE ratio and fatty acids composition. Filamentous fungi differ from S. cerevisiae in the lipid composition of their membranes and the structure of dolichols. The structure of dolichols profoundly affects the functioning of dolichol-dependent enzyme, DPMS.

LTK is an ER-resident receptor tyrosine kinase that regulates secretion

The endoplasmic reticulum (ER) is a major regulator of cellular proteostasis. However, only little is known about signaling molecules resident to this organelle. Centonze et al. identify LTK as the first ER-resident receptor tyrosine kinase and show that it stimulates secretory trafficking out of the ER.

The endoplasmic reticulum (ER) is a key regulator of cellular proteostasis because it controls folding, sorting, and degradation of secretory proteins. Much has been learned about how environmentally triggered signaling pathways regulate ER function, but only little is known about local signaling at the ER. The identification of ER-resident signaling molecules will help gain a deeper understanding of the regulation of ER function and thus of proteostasis. Here, we show that leukocyte tyrosine kinase (LTK) is an ER-resident receptor tyrosine kinase. Depletion of LTK as well as its pharmacologic inhibition reduces the number of ER exit sites and slows ER-to-Golgi transport. Furthermore, we show that LTK interacts with and phosphorylates Sec12. Expression of a phosphoablating mutant of Sec12 reduces the efficiency of ER export. Thus, LTK-to-Sec12 signaling represents the first example of an ER-resident signaling module with the potential to regulate proteostasis.

Clinical Metagenomic Sequencing for Diagnosis of Meningitis and Encephalitis

BACKGROUND

Metagenomic next-generation sequencing (NGS) of cerebrospinal fluid (CSF) has the potential to identify a broad range of pathogens in a single test.

METHODS

In a 1-year, multicenter, prospective study, we investigated the usefulness of metagenomic NGS of CSF for the diagnosis of infectious meningitis and encephalitis in hospitalized patients. All positive tests for pathogens on metagenomic NGS were confirmed by orthogonal laboratory testing. Physician feedback was elicited by teleconferences with a clinical microbial sequencing board and by surveys. Clinical effect was evaluated by retrospective chart review.

RESULTS

We enrolled 204 pediatric and adult patients at eight hospitals. Patients were severely ill: 48.5% had been admitted to the intensive care unit, and the 30-day mortality among all study patients was 11.3%. A total of 58 infections of the nervous system were diagnosed in 57 patients (27.9%). Among these 58 infections, metagenomic NGS identified 13 (22%) that were not identified by clinical testing at the source hospital. Among the remaining 45 infections (78%), metagenomic NGS made concurrent diagnoses in 19. Of the 26 infections not identified by metagenomic NGS, 11 were diagnosed by serologic testing only, 7 were diagnosed from tissue samples other than CSF, and 8 were negative on metagenomic NGS owing to low titers of pathogens in CSF. A total of 8 of 13 diagnoses made solely by metagenomic NGS had a likely clinical effect, with 7 of 13 guiding treatment.

CONCLUSIONS

Routine microbiologic testing is often insufficient to detect all neuroinvasive pathogens. In this study, metagenomic NGS of CSF obtained from patients with meningitis or encephalitis improved diagnosis of neurologic infections and provided actionable information in some cases. (Funded by the National Institutes of Health and others; PDAID ClinicalTrials.gov number, NCT02910037.).

Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity

Over the previous decades, numerous studies focused on how oceanic islands have contributed to determine the phylogenetic relationships and times of origin and diversification of different endemic lineages. The Macaronesian Islands (i.e., Azores, Madeira, Selvagens, Canaries, and Cabo Verde), harbour biotas with exceptionally high levels of endemism. Within the region, the vascular plants and reptiles constitute two of the most important radiations. In this study we compare relevant published phylogenetic data and diversification rates retrieved within Cabo Verde endemic lineages and discuss the importance of choosing appropriate phylogeny-based methods to investigate diversification dynamics on islands. From this selective literature-based review, we summarize the software packages used in Macaronesian studies and discuss their adequacy considering the published data to obtain well-supported phylogenies in the target groups. We further debate the importance of Next Generation Sequencing (NGS), to investigate the evolutionary processes of diversification in the Macaronesian Islands. Analysis of genomic data provides phylogenetic resolution for rapidly evolving species radiations, suggesting a great potential to improve the phylogenetic signal and divergence time estimates in insular lineages. The most important Macaronesian reptile radiations provide good case-studies to compare classical phylogenetic methods with new tools, such as phylogenomics, revealing a high value for research on this hotspot area.

A cluster of the first reported Plasmodium ovale spp. infections in Peru occuring among returning UN peace-keepers, a review of epidemiology, prevention and diagnostic challenges in nonendemic regions

BACKGROUND

Plasmodium ovale curtisi and Plasmodium ovale wallikeri are regarded as less virulent forms of malaria with a geographic distribution including Southeast Asia, Central and West Africa, and is increasingly reported as an infection in returning travellers. A species of malaria that may have delayed or relapsing presentations similar to Plasmodium vivax, the clinical presentation of P. ovale spp. has been described to have prepatent periods of 2 weeks or slightly longer with reports of relapse following primary infection out to 8-9 months. This presentation may be obscured further in the setting of anti-malarial exposure, with report of delayed primary infection out to 4 years. Presented is a cluster of 4 imported P. ovale spp. cases in returning Peruvian military personnel assigned to United Nations peace-keeping operations in the Central African Republic.

CASE PRESENTATION

From January to December 2016, Peruvian peace-keepers were deployed in support of United Nations (UN) operations in the Central African Republic (CAR). While serving abroad, Navy, Army, and Air Force members experienced 223 episodes of Plasmodium falciparum malaria following interruption of prophylaxis with mefloquine. Diagnosis was made using rapid diagnostics tests (RDTs) and/or smear with no coinfections identified. Cases of malaria were treated with locally-procured artemether-lumefantrine. Returning to Peru in January 2017, 200 peace-keepers were screened via thick and thin smear while on weekly mefloquine prophylaxis with only 1 showing nucleic acid within red blood cells consistent with Plasmodium spp. and 11 reporting syndromes of ill-defined somatic complaints. Between a period of 5 days to 11 months post return, 4 cases of P. ovale spp. were diagnosed using smear and polymerase chain reaction (PCR) following febrile complaints. All cases were subsequently treated with chloroquine and primaquine, with cure of clinical disease and documented clearance of parasitaemia.

CONCLUSION

These patients represent the first imported cases in Peru of this species of malaria as well as highlight the challenges in implementing population level prophylaxis in a deployed environment, and the steps for timely diagnosis and management in a non-endemic region where risk of introduction for local transmission exists.

Evidence for both Intermittent and Persistent Compartmentalization of HIV-1 in the Female Genital Tract

HIV-1 has been shown to evolve independently in different anatomical compartments, but studies in the female genital tract have been inconclusive. Here, we examined evidence of compartmentalization using HIV-1 subtype C envelope (Env) glycoprotein genes (gp160) obtained from matched cervicovaginal lavage (CVL) and plasma samples over 2 to 3 years of infection. HIV-1 gp160 amplification from CVL was achieved for only 4 of 18 acutely infected women, and this was associated with the presence of proinflammatory cytokines and/or measurable viremia in the CVL. Maximum likelihood trees and divergence analyses showed that all four individuals had monophyletic compartment-specific clusters of CVL- and/or plasma-derived gp160 sequences at all or some time points. However, two participants (CAP177 and CAP217) had CVL gp160 diversity patterns that differed from those in plasma and showed restricted viral flow from the CVL. Statistical tests of compartmentalization revealed evidence of persistent compartment-specific gp160 evolution in CAP177, while in CAP217 this was intermittent. Lastly, we identified several Env sites that distinguished viruses in these two compartments; for CAP177, amino acid differences arose largely through positive selection, while insertions/deletions were more common in CAP217. In both cases these differences contributed to substantial charge changes spread across the Env. Our data indicate that, in some women, HIV-1 populations within the genital tract can have Env genetic features that differ from those of viruses in plasma, which could impact the sensitivity of viruses in the genital tract to vaginal microbicides and vaccine-elicited antibodies.IMPORTANCE Most HIV-1 infections in sub-Saharan Africa are acquired heterosexually through the genital mucosa. Understanding the properties of viruses replicating in the female genital tract, and whether these properties differ from those of more commonly studied viruses replicating in the blood, is therefore important. Using longitudinal CVL and plasma-derived sequences from four HIV-1 subtype C-infected women, we found fewer viral migrations from the genital tract to plasma than in the opposite direction, suggesting a mucosal sieve effect from the genital tract to the blood compartment. Evidence for both persistent and intermittent compartmentalization between the genital tract and plasma viruses during chronic infection was detected in two of four individuals, perhaps explaining previously conflicting findings. In cases where compartmentalization occurred, comparison of CVL- and plasma-derived HIV sequences indicated that distinct features of viral populations in the CVL may affect the efficacy of microbicides and vaccines designed to provide mucosal immunity.

Exploring the sequence, function, and evolutionary space of protein superfamilies using sequence similarity networks and phylogenetic reconstructions

Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions by enabling the observation and investigation of complex sequence-structure-function and evolutionary relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) and phylogenetic reconstructions to map the functional divergence and evolutionary history of protein superfamilies. We exemplify this approach using the nitroreductase (NTR) flavoenzyme superfamily, demonstrating that SSN investigations can provide a rapid and effective means to classify groups of proteins, expose sequence similarity relationships across the global scale of a protein superfamily, and efficiently support detailed phylogenetic analyses. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes, their evolution, and their associated physiological roles.

Genome sequence of the corn leaf aphid (Rhopalosiphum maidis Fitch)

BACKGROUND

The corn leaf aphid (Rhopalosiphum maidis Fitch) is the most economically damaging aphid pest on maize (Zea mays), one of the world's most important grain crops. In addition to causing direct damage by removing photoassimilates, R. maidis transmits several destructive maize viruses, including maize yellow dwarf virus, barley yellow dwarf virus, sugarcane mosaic virus, and cucumber mosaic virus.

FINDINGS

The genome of a parthenogenetically reproducing R. maidis clone was assembled with a combination of Pacific Biosciences (207-fold coverage) and Illumina (83-fold coverage) sequencing. The 689 assembled contigs, which have an N50 size of 9.0 megabases (Mb) and a low level of heterozygosity, were clustered using Phase Genomics Hi-C interaction maps. Consistent with the commonly observed 2n = 8 karyotype of R. maidis, most of the contigs (473 spanning 321 Mb) were successfully oriented into 4 scaffolds. The genome assembly captured the full length of 95.8% of the core eukaryotic genes, indicating that it is highly complete. Repetitive sequences accounted for 21.2% of the assembly, and a total of 17,629 protein-coding genes were predicted with integrated evidence from ab initio and homology-based gene predictions and transcriptome sequences generated with both Pacific Biosciences and Illumina. An analysis of likely horizontally transferred genes identified 2 from bacteria, 7 from fungi, 2 from protozoa, and 9 from algae. Repeat elements, transposons, and genes encoding likely detoxification enzymes (cytochrome P450s, glutathione S-transferases, carboxylesterases, uridine diphosphate-glucosyltransferases, and ABC transporters) were identified in the genome sequence. Other than Buchnera aphidicola (642,929 base pairs, 602 genes), no endosymbiont bacteria were found in R. maidis.

CONCLUSIONS

A high-quality R. maidis genome was assembled at the chromosome level. This genome sequence will enable further research related to ecological interactions, virus transmission, pesticide resistance, and other aspects of R. maidis biology. It also serves as a valuable resource for comparative investigation of other aphid species.

A whole genome scan of SNP data suggests a lack of abundant hard selective sweeps in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

The pathogenic fungus Sclerotinia sclerotiorum infects over 600 species of plant. It is present in numerous environments throughout the world and causes significant damage to many agricultural crops. Fragmentation and lack of gene flow between populations may lead to population sub-structure. Within discrete recombining populations, positive selection may lead to a ‘selective sweep’. This is characterised by an increase in frequency of a favourable allele leading to reduction in genotypic diversity in a localised genomic region due to the phenomenon of genetic hitchhiking. We aimed to assess whether isolates of S. sclerotiorum from around the world formed genotypic clusters associated with geographical origin and to determine whether signatures of population-specific positive selection could be detected. To do this, we sequenced the genomes of 25 isolates of S. sclerotiorum collected from four different continents–Australia, Africa (north and south), Europe and North America (Canada and the northen United States) and conducted SNP based analyses of population structure and selective sweeps. Among the 25 isolates, there was evidence for two major population clusters. One of these consisted of 11 isolates from Canada, the USA and France (population 1), and the other consisted of nine isolates from Australia and one from Morocco (population 2). The rest of the isolates were genotypic outliers. We found that there was evidence of outcrossing in these two populations based on linkage disequilibrium decay. However, only a single candidate selective sweep was observed, and it was present in population 2. This sweep was close to a Major Facilitator Superfamily transporter gene, and we speculate that this gene may have a role in nutrient uptake from the host. The low abundance of selective sweeps in the S. sclerotiorum genome contrasts the numerous examples in the genomes of other fungal pathogens. This may be a result of its slow rate of evolution and low effective recombination rate due to self-fertilisation and vegetative reproduction.

Genomic profile of Brazilian methicillin-resistant Staphylococcus aureus resembles clones dispersed worldwide

PURPOSE

Comparative genomic analysis of strains may help us to better understand the wide diversity of their genetic profiles. The aim of this study was to analyse the genomic features of the resistome and virulome of Brazilian first methicillin-resistant Staphylococcus aureus (MRSA) isolates and their relationship to other Brazilian and international MRSA strains.

METHODOLOGY

The whole genomes of three MRSA strains previously isolated in Vitória da Conquista were sequenced, assembled, annotated and compared with other MRSA genomes. A phylogenetic tree was constructed and the pan-genome and accessory and core genomes were constructed. The resistomes and virulomes of all strains were identified.Results/Key findings. Phylogenetic analysis of all 49 strains indicated different clones showing high similarity. The pan-genome of the analysed strains consisted of 4484 genes, with 31 % comprising the gene portion of the core genome, 47 % comprising the accessory genome and 22 % being singletons. Most strains showed at least one gene related to virulence factors associated with immune system evasion, followed by enterotoxins. The strains showed multiresistance, with the most recurrent genes conferring resistance to beta-lactams, fluoroquinolones, aminoglycosides and macrolides.

CONCLUSIONS

Our comparative genomic analysis showed that there is no pattern of virulence gene distribution among the clones analysed in the different regions. The Brazilian strains showed similarity with clones from several continents.

HIV Subtype and Nef-Mediated Immune Evasion Function Correlate with Viral Reservoir Size in Early-Treated Individuals

The HIV accessory protein Nef modulates key immune evasion and pathogenic functions, and its encoding gene region exhibits high sequence diversity. Given the recent identification of early HIV-specific adaptive immune responses as novel correlates of HIV reservoir size, we hypothesized that viral factors that facilitate the evasion of such responses-namely, Nef genetic and functional diversity-might also influence reservoir establishment and/or persistence. We isolated baseline plasma HIV RNA-derived nef clones from 30 acute/early-infected individuals who participated in a clinical trial of early combination antiretroviral therapy (cART) (<6 months following infection) and assessed each Nef clone's ability to downregulate CD4 and human leukocyte antigen (HLA) class I in vitro We then explored the relationships between baseline clinical, immunological, and virological characteristics and the HIV reservoir size measured 48 weeks following initiation of suppressive cART (where the reservoir size was quantified in terms of the proviral DNA loads as well as the levels of replication-competent HIV in CD4⁺ T cells). Maximal within-host Nef-mediated downregulation of HLA, but not CD4, correlated positively with post-cART proviral DNA levels (Spearman's R = 0.61, P = 0.0004) and replication-competent reservoir sizes (Spearman's R = 0.36, P = 0.056) in univariable analyses. Furthermore, the Nef-mediated HLA downregulation function was retained in final multivariable models adjusting for established clinical and immunological correlates of reservoir size. Finally, HIV subtype B-infected persons (n = 25) harbored significantly larger viral reservoirs than non-subtype B-infected persons (2 infected with subtype CRF01_AE and 3 infected with subtype G). Our results highlight a potentially important role of viral factors-in particular, HIV subtype and accessory protein function-in modulating viral reservoir establishment and persistence.IMPORTANCE While combination antiretroviral therapies (cART) have transformed HIV infection into a chronic manageable condition, they do not act upon the latent HIV reservoir and are therefore not curative. As HIV cure or remission should be more readily achievable in individuals with smaller HIV reservoirs, achieving a deeper understanding of the clinical, immunological, and virological determinants of reservoir size is critical to eradication efforts. We performed a post hoc analysis of 30 participants of a clinical trial of early cART who had previously been assessed in detail for their clinical, immunological, and reservoir size characteristics. We observed that the HIV subtype and autologous Nef-mediated HLA downregulation function correlated with the viral reservoir size measured approximately 1 year post-cART initiation. Our findings highlight virological characteristics-both genetic and functional-as possible novel determinants of HIV reservoir establishment and persistence.

Positive Selection at Key Residues in the HIV Envelope Distinguishes Broad and Strain-Specific Plasma Neutralizing Antibodies

Millions of people are still being infected with HIV decades after the first recognition of the virus. Currently, no vaccine is able to elicit bNAbs that will prevent infection by global HIV strains. Several studies have implicated HIV Env diversity in the development of breadth. However, Env evolution in individuals who fail to develop breadth despite mounting potent strain-specific neutralizing responses has not been well defined. Using longitudinal neutralization, epitope mapping, and sequence data from 14 participants, we found that overall measures of viral diversity were similar in all donors. However, the number of positively selected sites within Env epitopes was higher in bNAb participants than in strain-specific donors. We further identified common sites that were positively selected as bNAbs developed. These data indicate that while viral diversity is required for breadth, this should be highly targeted to specific residues to shape the elicitation of bNAbs by vaccination.

The development of HIV broadly neutralizing antibodies (bNAbs) has previously been shown to be associated with viral evolution and high levels of genetic diversity in the HIV envelope (Env) glycoprotein. However, few studies have examined Env evolution in those who fail to develop neutralization breadth in order to assess whether bNAbs result from distinct evolutionary pathways. We compared Env evolution in eight HIV-1-infected participants who developed bNAbs to six donors with similar viral loads who did not develop bNAbs over three years of infection. We focused on Env V1V2 and C3V4, as these are major targets for both strain-specific neutralizing antibodies (nAbs) and bNAbs. Overall evolutionary rates (ranging from 9.92 × 10⁻³ to 4.1 × 10⁻² substitutions/site/year) and viral diversity (from 1.1% to 6.5%) across Env, and within targeted epitopes, did not distinguish bNAb donors from non-bNAb donors. However, bNAb participants had more positively selected residues within epitopes than those without bNAbs, and several of these were common among bNAb donors. A comparison of the kinetics of strain-specific nAbs and bNAbs indicated that selection pressure at these residues increased with the onset of breadth. These data suggest that highly targeted viral evolution rather than overall envelope diversity is associated with neutralization breadth. The association of shared positively selected sites with the onset of breadth highlights the importance of diversity at specific positions in these epitopes for bNAb development, with implications for the development of sequential and cocktail immunization strategies.

IMPORTANCE Millions of people are still being infected with HIV decades after the first recognition of the virus. Currently, no vaccine is able to elicit bNAbs that will prevent infection by global HIV strains. Several studies have implicated HIV Env diversity in the development of breadth. However, Env evolution in individuals who fail to develop breadth despite mounting potent strain-specific neutralizing responses has not been well defined. Using longitudinal neutralization, epitope mapping, and sequence data from 14 participants, we found that overall measures of viral diversity were similar in all donors. However, the number of positively selected sites within Env epitopes was higher in bNAb participants than in strain-specific donors. We further identified common sites that were positively selected as bNAbs developed. These data indicate that while viral diversity is required for breadth, this should be highly targeted to specific residues to shape the elicitation of bNAbs by vaccination.

Genomic Characterization and Virulence Potential of Two Fusarium oxysporum Isolates Cultured from the International Space Station

Two isolates of Fusarium oxysporum, ISS-F3 and ISS-F4, were cultured from the dining table on the International Space Station (ISS). Genomic analyses using EF-1α sequences, presence/absence of effector proteins, k-mer comparisons, and single nucleotide polymorphisms indicate that these two strains are genomically different from 65 known sequenced strains. Functional analysis revealed that ISS-F3/F4 had higher relative abundances of polyketide synthase domains than a non-plant-pathogenic soil isolate, used for biocontrol properties (Fo47), and a clinical isolate (FOSC-3a). Putative secondary metabolite analysis indicates that ISS-F3/F4 may produce yet-unreported polyketides and nonribosomal peptides. While genomic analysis showed that these ISS strains are unlikely to be plant pathogens, a virulence assay using an immunocompromised Caenorhabditis elegans model of fusariosis revealed that they were virulent and may represent opportunistic pathogens in animals, including humans. However, its effects on the health of immunocompromised humans warrant further study. IMPORTANCE This is the first study to isolate and characterize F. oxysporum isolates from a built environment, as well as one that has been exposed to space. The characterization and analysis of these two genomes may have important implications for the medical, agricultural, and food industries as well as for the health of the crew who coinhabit the ISS with these strains.

ProtParCon: A Framework for Processing Molecular Data and Identifying Parallel and Convergent Amino Acid Replacements

Studying parallel and convergent amino acid replacements in protein evolution is frequently used to assess adaptive evolution at the molecular level. Identifying parallel and convergent replacements involves multiple steps and computational routines, such as multiple sequence alignment, phylogenetic tree inference, ancestral state reconstruction, topology tests, and simulation of sequence evolution. Here, we present ProtParCon, a Python 3 package that provides a common interface for users to process molecular data and identify parallel and convergent amino acid replacements in orthologous protein sequences. By integrating several widely used programs for computational biology, ProtParCon implements general functions for handling multiple sequence alignment, ancestral-state reconstruction, maximum-likelihood phylogenetic tree inference, and sequence simulation. ProtParCon also contains a built-in pipeline that automates all these sequential steps, and enables quick identification of observed and expected parallel and convergent amino acid replacements under different evolutionary assumptions. The most up-to-date version of ProtParCon, including scripts containing user tutorials, the full API reference and documentation are publicly and freely available under an open source MIT License via GitHub. The latest stable release is also available on PyPI (the Python Package Index).

Inconsistent temporal patterns of genetic variation of HCV among high-risk subjects may impact inference of transmission networks

Hepatitis-C Virus (HCV) sequences are often used to establish networks of people who inject drugs (PWID). However, the degree to which within-host evolutionary dynamics affect those inferences has not been carefully studied. Here, we analyzed 702 longitudinally-sampled HCV E1 sequences from 88 HCV+ people who inject drugs (PWID) in the Baltimore Before and After Acute Study of Hepatitis (BBAASH) cohort. Individuals were tested for HCV RNA over multiple visits to the clinic, and the HCV E1 gene was sequenced for HCV+ samples. Genetic clustering was performed on the full set of sequences using a 3% genetic distance threshold to define epidemiological linkage. Maximum-likelihood (ML) phylogenies were inferred to assess evolutionary relationships. We found 22 clusters containing sequences sampled over five or more years (long-term clusters, LTC), of which 17 had >1 subject. In six of the multi-subject LTC, one subject had a sequence sampled >3 years earlier or later than the next-closest subject in the cluster (time-gap LTC). ML trees showed that, in three of the time-gap LTC, two subjects had identical sequences despite 7-10 years separating the sampling times. In four of the time-gap LTC for whom additional data were available, the subject with the later detected shared variant had both different variants and visits with no detectable HCV RNA (RNA-) prior to the appearance of the shared variant. In the subject with the earlier detection of the shared variant, different variants and RNA- visits were also detected in multiple cases subsequent to appearance of the shared variant. Complex patterns of shared viral variation among PWID reflect on-going re-infection, multiple transmission partners, and/or inconsistent detection of viral variants. Our results suggest that transmission events are currently underestimated by analysis of sequences at a single point in time.

Virulence factors, prevalence and potential transmission of extraintestinal pathogenic Escherichia coli isolated from different sources: recent reports

Extraintestinal pathogenic E. coli (ExPEC) are facultative pathogens that are part of the normal human intestinal flora. The ExPEC group includes uropathogenic E. coli (UPEC), neonatal meningitis E. coli (NMEC), sepsis-associated E. coli (SEPEC), and avian pathogenic E. coli (APEC). Virulence factors (VF) related to the pathogenicity of ExPEC are numerous and have a wide range of activities, from those related to bacteria colonization to those related to virulence, including adhesins, toxins, iron acquisition factors, lipopolysaccharides, polysaccharide capsules, and invasins, which are usually encoded on pathogenicity islands (PAIs), plasmids and other mobile genetic elements. Mechanisms underlying the dynamics of ExPEC transmission and the selection of virulent clones are still poorly understood and require further research. The time shift between colonization of ExPEC and the development of infection remains problematic in the context of establishing the relation between consumption of contaminated food and the appearance of first disease symptoms. What appears to be most difficult is to prove that ExPEC strains cause disease symptoms and to examine the mechanism of transition from the asymptomatic colonization of the intestines to the spreading of the bacteria outside the digestive system. A significant problem for researchers who are trying to ascribe ExPEC transmission to food, people or the environment is to draw the distinction between colonization of ExPEC and infection. Food safety is an important challenge for public health both at the production stage and in the course of its processing and distribution. Examination of the genetic similarity of ExPEC strains will allow to determine their origin from different sources. Many levels of genotyping have been proposed in which the typing of strains, plasmids and genes is compared in order to obtain a more complete picture of this complex problem. The aim of our study was to characterize E. coli strains isolated from humans, animals and food for the presence of bacterial genes encoding virulence factors such as toxins, and iron acquisition systems (siderophores) in the context of an increasing spread of ExPEC infections.

Misidentification of Burkholderia pseudomallei and Other Burkholderia Species From Pediatric Infections in Mexico

Burkholderia pseudomallei and Burkholderia cepacia complex are poorly studied in Mexico. The genotypic analysis of 38 strains isolated from children with pneumonia were identified and showed that both Burkholderia groups were present in patients. From our results, it is plausible to suggest that new species are among the analyzed strains.

Adding function to the genome of African Salmonella Typhimurium ST313 strain D23580

Salmonella Typhimurium sequence type (ST) 313 causes invasive nontyphoidal Salmonella (iNTS) disease in sub-Saharan Africa, targeting susceptible HIV+, malarial, or malnourished individuals. An in-depth genomic comparison between the ST313 isolate D23580 and the well-characterized ST19 isolate 4/74 that causes gastroenteritis across the globe revealed extensive synteny. To understand how the 856 nucleotide variations generated phenotypic differences, we devised a large-scale experimental approach that involved the global gene expression analysis of strains D23580 and 4/74 grown in 16 infection-relevant growth conditions. Comparison of transcriptional patterns identified virulence and metabolic genes that were differentially expressed between D23580 versus 4/74, many of which were validated by proteomics. We also uncovered the S. Typhimurium D23580 and 4/74 genes that showed expression differences during infection of murine macrophages. Our comparative transcriptomic data are presented in a new enhanced version of the Salmonella expression compendium, SalComD23580: http://bioinf.gen.tcd.ie/cgi-bin/salcom_v2.pl. We discovered that the ablation of melibiose utilization was caused by three independent SNP mutations in D23580 that are shared across ST313 lineage 2, suggesting that the ability to catabolize this carbon source has been negatively selected during ST313 evolution. The data revealed a novel, to our knowledge, plasmid maintenance system involving a plasmid-encoded CysS cysteinyl-tRNA synthetase, highlighting the power of large-scale comparative multicondition analyses to pinpoint key phenotypic differences between bacterial pathovariants.

Computational Prediction of De Novo Emerged Protein-Coding Genes

De novo genes, that is, protein-coding genes originating from previously noncoding sequence, have gone from being considered impossibly unlikely to being recognized as an important source of genetic novelty in eukaryotic genomes. It is clear that de novo gene evolution is a rare but consistent feature of eukaryotic genomes, being detected in every genome studied. However, different studies often use different computational methods, and the numbers and identities of the detected genes vary greatly. Here we present a coherent protocol for the computational identification of de novo genes by comparative genomics. The method described uses homology searches, identification of syntenic regions, and ancestral sequence reconstruction to produce high-confidence candidates with robust evidence of de novo emergence. It is designed to be easily applicable given the basic knowledge of bioinformatic tools and scalable so that it can be applied on large and small datasets.

A Not-So-Long Introduction to Computational Molecular Evolution

In this chapter, we give a not-so-long and self-contained introduction to computational molecular evolution. In particular, we present the emergence of the use of likelihood-based methods, review the standard DNA substitution models, and introduce how model choice operates. We also present recent developments in inferring absolute divergence times and rates on a phylogeny, before showing how state-of-the-art models take inspiration from diffusion theory to link population genetics, which traditionally focuses at a taxonomic level below that of the species, and molecular evolution. Although this is not a cookbook chapter, we try and point to popular programs and implementations along the way.

Polymorphism analysis of the apxIA gene of Actinobacillus pleuropneumoniae serovar 5 isolated in swine herds from Brazil

The bacterium Actinobacillus pleuropneumoniae is the etiological agent of Contagious Porcine Pleuropneumonia, a disease responsible for economic losses in the swine industry worldwide. A. pleuropneumoniae is capable of producing proteinaceous exotoxins responsible for inducing hemorrhagic lesions, one of which is ApxI. Few studies have conducted an in-depth evaluation of polymorphisms of the nucleotides that make up the ApxI toxin gene. Here we analyze the polymorphisms of the apxIA gene region of A. pleuropneumoniae serovar 5 isolated from swine in different regions in Brazil and report the results of molecular sequencing and phylogenetic analysis. Analysis of the apxIA gene in 60 isolates revealed the presence of genetic diversity and variability. The polymorphisms in the nucleotide sequences determined the grouping of the Brazilian sequences and five more sequences from the GenBank database into 14 different haplotypes, which formed three main groups and revealed the presence of mutations in the nucleotide sequences. The estimation of selection pressures suggests the occurrence of genetic variations by positive selective pressure on A. pleuropneumoniae in large groups of animals in relatively small spaces. These conditions presumably favor the horizontal dissemination of apxIA gene mutations within bacterial populations with host reservoirs. As a result, the same serovar can demonstrate different antigenic capacities due to mutations in the apxIA gene. These alterations in sequences of the apxIA gene could occur in other areas of countries with intense swine production, which could lead to differences in the pathogenicity and immunogenicity of each serovar and have implications for the clinical status or diagnosis of A. pleuropneumoniae.