Adding unaligned sequences into an existing alignment using MAFFT and LAST

Stictaflakusiorum and S.kukwae-two additional new species from the Neotropics (Peltigerales, Peltigeraceae)

Two additional species of Sticta are described as new to science based on material from Bolivia and Peru and supported by phylogenetic analysis of the fungal ITS barcoding marker. The two new species represent lineages within clade I on the global Sticta phylogeny. Stictaflakusiorum Ossowska, B. Moncada & Lücking is a species in the S.humboldtii morphodeme and is characterized by lobes partly to entirely covered with white hairs, also covering the margins of submarginal and laminal apothecia, and the scabrid basal membrane of cyphellae, which is white to yellow, or partly brown, and when yellow K+ purple. The taxon was discovered at a single locality in Bolivia, but it is closely related to a potentially new Sticta species from Peru, which is here left undescribed. The other new species, S.kukwae Ossowska, Magain & Sérus., belongs to the S.weigelii morphodeme. It has lobes with sinuous margins and dark, palmate to corymbose phyllidia. It was collected at several locations in Peru and a single locality in Bolivia.

Genomic characterization of a dog-mediated rabies outbreak in El Pedregal, Arequipa, Peru

BACKGROUND

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru.

METHODS

We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were generated, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins.

RESULTS

Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to the new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil.

CONCLUSION

Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Impact of Alignments on the Accuracy of Protein Subcellular Localization Predictions

Alignments in bioinformatics refer to the arrangement of sequences to identify regions of similarity that can indicate functional, structural, or evolutionary relationships. They are crucial for bioinformaticians as they enable accurate predictions and analyses in various applications, including protein subcellular localization. The predictive model used in this article is based on a deep - convolutional architecture. We tested configurations of Deep N-to-1 convolutional neural networks of various depths and widths during experimentation for the evaluation of better-performing values across a diverse set of eight classes. For without alignment assessment, sequences are encoded using one-hot encoding, converting each character into a numerical representation, which is straightforward for non-numerical data and useful for machine learning models. For with alignments assessment, multiple sequence alignments (MSAs) are created using PSI-BLAST, capturing evolutionary information by calculating frequencies of residues and gaps. The average difference in peak performance between models with alignments and without alignments is approximately 15.82%. The average difference in the highest accuracy achieved with alignments compared with without alignments is approximately 15.16%. Thus, extensive experimentation indicates that higher alignment accuracy implies a more reliable model and improved prediction accuracy, which can be trusted to deliver consistent performance across different layers and classes of subcellular localization predictions. This research provides valuable insights into prediction accuracies with and without alignments, offering bioinformaticians an effective tool for better understanding while potentially reducing the need for extensive experimental validations. The source code and datasets are available at http://distilldeep.ucd.ie/SCL8/.

High resistance barrier and prophylactic protection in preclinical models of SARS-CoV-2 with two siRNA combination

RNA interference is a natural antiviral mechanism that could be harnessed to combat SARS-CoV-2 infection by targeting and destroying the viral RNA. We identified potent lipophilic small interfering RNA (siRNA) conjugates targeting highly conserved regions of SARS-CoV-2 outside of the spike-encoding region capable of achieving ≥3-log viral reduction. Serial passaging studies demonstrated that a two-siRNA combination prevented development of resistance compared to a single siRNA approach. Viral resistance to single siRNA treatment occurred due to emergence of point mutations at critical positions required for siRNA-mediated target binding and cleavage, which led to a loss of siRNA efficacy. With a two-siRNA combination, emergence of mutations within the siRNA binding site was abolished. When delivered intranasally, two-siRNA combination protected Syrian hamsters from weight loss and lung pathology by viral infection upon prophylactic administration but not following onset of infection. Together, the data support potential utility of RNAi as a prophylactic approach with high resistance barrier to counteract SARS-CoV-2 emergent variants and complement vaccination. Most importantly, given that the siRNAs can be rapidly developed from a new pathogen sequence, this strategy has implications as a new type of preventive medicine that may protect against future coronavirus pandemics.

Human Circovirus in Patients with Hepatitis, Hong Kong

Circovirus human is a new viral species that includes the human circovirus (HCirV), which has been linked to hepatitis in immunocompromised persons. We investigated prevalence of HCirV infection in 278 patients with hepatitis and 184 asymptomatic persons using real-time PCR and sequencing assays. HCirV viremia and sequences were found in 8 (2.9%) hepatitis patients and no asymptomatic patients. Alternate causes of hepatitis (hepatitis E and cholangitis) were clearly identifiable in 2 HCirV-infected patients. HCirV could not be ruled out as a contributor to hepatitis in the remaining 6 patients, 4 of whom were immunocompromised. Persistent infections were documented in 3 patients, but only 1 had relapsing hepatitis. One HCirV patient displayed symptoms of an infectious mononucleosis-like syndrome. Isolates clustered with known HCirV strains from France and China. HCirV-derived virus-like particles bound to PLC/PRF/5 and Hep-G2 human hepatoma cells but not to lung epithelial cells, indicating hepatic tropism.

Regulators of aerobic and anaerobic methane oxidation in two pristine temperate peatland types

Despite covering <5% of Earth's terrestrial area, peatlands are crucial for global carbon storage and are hot spots of methane cycling. This study examined the dynamics of aerobic and anaerobic methane oxidation in two undisturbed peatlands: a fen and a spruce swamp forest. Using microcosm incubations, we investigated the effect of ammonium addition, at a level similar to current N pollution processes, on aerobic methane oxidation. Our findings revealed higher methane consumption rates in fen compared to swamp peat, but no effect of ammonium amendment on methane consumption was found. Members of Methylocystis and Methylocella were the predominant methanotrophs in both peatlands. Furthermore, we explored the role of ferric iron and sulfate as electron acceptors for the anaerobic oxidation of methane (AOM). AOM occurred without the addition of an external electron acceptor in the fen, but not in the swamp peat. AOM was stimulated by sulfate and ferric iron addition in the swamp peat and inhibited by ferric iron in the fen. Our findings suggest that aerobic methane oxidizers are not N-limited in these peatlands and that there is an intrinsic potential for AOM in these environments, partially facilitated by ferric iron and sulfate acting as electron acceptors.

Emergence of highly pathogenic avian influenza viruses H5N1 and H5N5 in white-tailed eagles, 2021-2023

Highly pathogenic avian influenza (HPAI) poses a substantial threat to several raptors. Between 2021 and 2023, HPAI viruses (HPAIVs) of the Goose/Guangdong lineage H5 clade 2.3.4.4b became widespread in wild birds in Norway, and H5N1 and H5N5 viruses were detected in 31 white-tailed eagles (Haliaeetus albicilla, WTEs). Post-mortem examinations of four WTEs revealed no macroscopic pathological findings. Microscopic examinations showed the presence of myocardial and splenic necroses and a few lesions in the brain. In situ hybridization revealed the presence of the virus in several organs, suggesting a multisystemic infection. The detection of HPAIV H5N5 in a WTE in February 2022 marked the first recorded occurrence of this subtype in Norway. Since then, the virus has persisted, sporadically being detected in WTEs and other wild bird species. Phylogenetic analyses reveal that at least two distinct incursions of HPAIV H5N1 Eurasian (EA) genotype C affected WTEs, likely introduced by migratory birds from Eurasia and seabirds entering from Western and Central Europe. Some WTE isolates from 2021 to 2022 clustered with those from Canada and Ireland, aligning with the transatlantic spread of H5N1. Others were related to the 2021 mass mortality of great skuas in the UK or outbreaks in seabird populations, including gannets, gulls and terns, during 2022 in the North Sea region. This suggests that the WTEs were likely preying on the affected birds. Our study highlights that WTEs can act as sentinels for some HPAIV strains, but the absence of several known circulating genotypes in WTEs suggests varying pathogenic effects on this species.

Identifying the multiple drivers of cactus diversification

Our understanding of the complexity of forces at play in the rise of major angiosperm lineages remains incomplete. The diversity and heterogeneous distribution of most angiosperm lineages is so extraordinary that it confounds our ability to identify simple drivers of diversification. Using machine learning in combination with phylogenetic modelling, we show that five separate abiotic and biotic variables significantly contribute to the diversification of Cactaceae. We reconstruct a comprehensive phylogeny, build a dataset of 39 abiotic and biotic variables, and predict the variables of central importance, while accounting for potential interactions between those variables. We use state-dependent diversification models to confirm that five abiotic and biotic variables shape diversification in the cactus family. Of highest importance are diurnal air temperature range, soil sand content and plant size, with lesser importance identified in isothermality and geographic range size. Interestingly, each of the estimated optimal conditions for abiotic variables were intermediate, indicating that cactus diversification is promoted by moderate, not extreme, climates. Our results reveal the potential primary drivers of cactus diversification, and the need to account for the complexity underlying the evolution of angiosperm lineages.

Genomic Characterization of a Dog-Mediated Rabies Outbreak in El Pedregal, Arequipa, Peru

Background

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru.

Methods

We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were analyzed, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins.

Results

Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to a new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil.

Conclusion

Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences. Evidence for specific base pairs turns out to be more informative than a position-wise profile for the conservation of the pairing status. A comparison with chemical probing data, furthermore, strongly suggests that phylogenetic base pairing data are more informative than position-specific data on (un)pairedness as obtained from chemical probing experiments. In this context we demonstrate, in addition, that the conversion of signal from probing data into pseudo-energies is possible using thermodynamic structure predictions as a reference instead of known RNA structures.

A giant virus infecting the amoeboflagellate Naegleria

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Additional new species and new records of the genus Sticta (lichenised Ascomycota, lobarioid Peltigeraceae) from Bolivia

Four species of the genus Sticta are described as new from Bolivia, based on morphological examination and phylogenetic analysis of the fungal ITS barcoding marker. Additionally, two species are reported as new to Bolivia (their identification confirmed by molecular data) and one previously reported species is confirmed by molecular data for the first time. Detailed morphological and anatomical descriptions are provided for all new species. Two of the new species, S.isidiolobulata Ossowska, B. Moncada, Lücking & Kukwa and S.madidiensis Ossowska, B. Moncada, Lücking & Kukwa belong to clade I, as defined in previous studies. In contrast, S.montepunkuensis Ossowska, B. Moncada, Lücking & Kukwa and S.macrolobata Ossowska, B. Moncada, Lücking & Kukwa, also described here as new to science, belong to clade III. Stictaisidiolobulata has an irregular to suborbicular thallus of medium size, with isidia developing into spathulate lobules, cyanobacterial photobiont and apothecia with entire to weakly-crenate margins. The large irregular thallus of the cyanobacteria-associated S.macrolobata has broad lobes, apothecia with verrucous to tomentose margins and cyphellae with raised margins, whereas S.madidiensis has a medium-sized, palmate to irregular thallus with a stipe, but without vegetative propagules and apothecia. Stictamontepunkuensis has large and irregular thalli with green algae as photobiont, apothecia with crenate to verrucous margins and urceolate cyphellae with a wide pore and a scabrid basal membrane. Two species, S.beauvoisii Delise and S.riparia Merc.-Díaz are reported as new to Bolivia (the latter also as new to South America) and belong to clade III. Stictatomentosa (Sw.) Ach., species confirmed from Bolivia by molecular data, belongs to clade II. Stictabeauvoisii is characterised by a smooth yellowish-brown upper surface with darker apices and abundant, marginal isidia and a brown lower surface with golden-chocolate brown primary tomentum and sparse, golden-brown rhizines. Stictariparia has a strongly branched thallus, with undulate lobes and abundant, marginal, palmate, grey to dark brown phyllidia and greyish-brown lower surface with the primary tomentum absent towards the margins. Stictatomentosa has palmate, bluish thalli with white cilia and abundant, submarginal apothecia and creamy-white lower surface with a sparse, white primary tomentum.

Five New Species of the Genus Hymenogaster (Hymenogastraceae, Agaricales) from Northern China

In this study, five new species from China, Hymenogaster latisporus, H. minisporus, H. papilliformis, H. perisporius, and H. variabilis, are described and illustrated based on morphological and molecular evidence. Hymenogaster latisporus was distinguished from other species of the genus by the subglobose, broad ellipsoidal, ovoid basidiospores (average = 13.7 μm × 11.6 μm) with sparse verrucose and ridge-like ornamentation (1-1.2 μm high); H. minisporus by the ellipsoidal to broadly ellipsoidal and small basidiospores (average = 11.7 μm × 9.5 μm); H. papilliformis was characterized by the whitish to cream-colored basidiomes, and broadly fusiform to citriform basidiospores with a pronounced apex (2-3 μm, occasionally up to 4 μm high), papillary, distinct warts and ridges, and pronounced appendix (2-3 μm long); H. perisporius by the dirty white to pale yellow basidiomes, broad ellipsoidal to ellipsoidal, and yellow-brown to dark-brown basidiospores with warts and gelatinous perisporium; H. variabilis by the peridium with significant changes in thickness (167-351 μm), and broad ellipsoidal to subglobose basidiospores ornamented with sparse warts and ridges. An ITS/LSU-based phylogenetic analysis supported the erection of the five new species. A key for Hymenogaster species from northern China is provided.

Global and genetic diversity of SARS-CoV-2 in wastewater

The analysis of SARS-CoV-2 in wastewater has enabled us to better understand the spread and evolution of the virus worldwide. To deepen our understanding of its epidemiological and genomic characteristics, we analyzed 10,147 SARS-CoV-2 sequences from 5 continents and 21 countries that were deposited in the GISAID database up until January 31, 2023. Our results revealed over 100 independent lineages of the virus circulating in water samples from March 2020 to January 2023, including variants of interest and concern. We observed four clearly defined periods of global distribution of these variants over time, with one variant being replaced by another. Interestingly, we found that SARS-CoV-2 water-borne sequences from different countries had a close phylogenetic relationship. Additionally, 40 SARS-CoV-2 water-borne sequences from Europe and the USA did not show any phylogenetic relationship with SARS-CoV-2 human sequences. We also identified a significant number of non-synonymous mutations, some of which were detected in previously reported cryptic lineages. Among the countries analyzed, France and the USA showed the highest degree of sequence diversity, while Austria reported the highest number of genomes (6,296). Our study provides valuable information about the epidemiological and genomic diversity of SARS-CoV-2 in wastewater, which can be employed to support public health initiatives and preparedness.

Parvovirus B19 and Human Parvovirus 4 Encode Similar Proteins in a Reading Frame Overlapping the VP1 Capsid Gene

Viruses frequently contain overlapping genes, which encode functionally unrelated proteins from the same DNA or RNA region but in different reading frames. Yet, overlapping genes are often overlooked during genome annotation, in particular in DNA viruses. Here we looked for the presence of overlapping genes likely to encode a functional protein in human parvovirus B19 (genus Erythroparvovirus), using an experimentally validated software, Synplot2. Synplot2 detected an open reading frame, X, conserved in all erythroparvoviruses, which overlaps the VP1 capsid gene and is under highly significant selection pressure. In a related virus, human parvovirus 4 (genus Tetraparvovirus), Synplot2 also detected an open reading frame under highly significant selection pressure, ARF1, which overlaps the VP1 gene and is conserved in all tetraparvoviruses. These findings provide compelling evidence that the X and ARF1 proteins must be expressed and functional. X and ARF1 have the exact same location (they overlap the region of the VP1 gene encoding the phospholipase A2 domain), are both in the same frame (+1) with respect to the VP1 frame, and encode proteins with similar predicted properties, including a central transmembrane region. Further studies will be needed to determine whether they have a common origin and similar function. X and ARF1 are probably translated either from a polycistronic mRNA by a non-canonical mechanism, or from an unmapped monocistronic mRNA. Finally, we also discovered proteins predicted to be expressed from a frame overlapping VP1 in other species related to parvovirus B19: porcine parvovirus 2 (Z protein) and bovine parvovirus 3 (X-like protein).

Molecular evolution of 2022 multi-country outbreak-causing monkeypox virus Clade IIb

The monkeypox virus (Mpoxv) Clade IIb viruses that caused an outbreak in 2017-18 in Nigeria and its genetically related viruses have been detected in many countries and caused multi-country outbreak in 2022. Since the pandemic-causing Mpoxv Clade IIb viruses are closely related to Clade IIa viruses which mostly cause endemic, the Clade IIb Mpoxv might have certain specific genetic variations that are still largely unknown. Here, we have systematically analyzed genetic alterations in different clades of Mpox viruses. The results suggest that the Mpoxv Clade IIb have genetic variations in terms of genomic gaps, frameshift mutations, in-frame nonsense mutations, amino acid tandem repeats, and APOBEC3 mutations. Further, we observed specific genetic variations in the multiple genes specific for Clade I and Clade IIb, and exclusive genetic variations for Clade IIa and Clade IIb. Collectively, findings shed light on the evolution and genetic variations in the outbreak of 2022 causing Mpoxv Clade IIb.

Evolutionary dynamics of the LTR-retrotransposon crapaud in the Podospora anserina species complex and the interaction with repeat-induced point mutations

BACKGROUND

The genome of the filamentous ascomycete Podospora anserina shows a relatively high abundance of retrotransposons compared to other interspersed repeats. The LTR-retrotransposon family crapaud is particularly abundant in the genome, and consists of multiple diverged sequence variations specifically localized in the 5' half of both long terminal repeats (LTRs). P. anserina is part of a recently diverged species-complex, which makes the system ideal to classify the crapaud family based on the observed LTR variation and to study the evolutionary dynamics, such as the diversification and bursts of the elements over recent evolutionary time.

RESULTS

We developed a sequence similarity network approach to classify the crapaud repeats of seven genomes representing the P. anserina species complex into 14 subfamilies. This method does not utilize a consensus sequence, but instead it connects any copies that share enough sequence similarity over a set sequence coverage. Based on phylogenetic analyses, we found that the crapaud repeats likely diversified in the ancestor of the complex and have had activity at different time points for different subfamilies. Furthermore, while we hypothesized that the evolution into multiple subfamilies could have been a direct effect of escaping the genome defense system of repeat induced point mutations, we found this not to be the case.

CONCLUSIONS

Our study contributes to the development of methods to classify transposable elements in fungi, and also highlights the intricate patterns of retrotransposon evolution over short timescales and under high mutational load caused by nucleotide-altering genome defense.

Unraveling viral drug targets: a deep learning-based approach for the identification of potential binding sites

The coronavirus disease 2019 (COVID-19) pandemic has spurred a wide range of approaches to control and combat the disease. However, selecting an effective antiviral drug target remains a time-consuming challenge. Computational methods offer a promising solution by efficiently reducing the number of candidates. In this study, we propose a structure- and deep learning-based approach that identifies vulnerable regions in viral proteins corresponding to drug binding sites. Our approach takes into account the protein dynamics, accessibility and mutability of the binding site and the putative mechanism of action of the drug. We applied this technique to validate drug targeting toward severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein S. Our findings reveal a conformation- and oligomer-specific glycan-free binding site proximal to the receptor binding domain. This site comprises topologically important amino acid residues. Molecular dynamics simulations of Spike in complex with candidate drug molecules bound to the potential binding sites indicate an equilibrium shifted toward the inactive conformation compared with drug-free simulations. Small molecules targeting this binding site have the potential to prevent the closed-to-open conformational transition of Spike, thereby allosterically inhibiting its interaction with human angiotensin-converting enzyme 2 receptor. Using a pseudotyped virus-based assay with a SARS-CoV-2 neutralizing antibody, we identified a set of hit compounds that exhibited inhibition at micromolar concentrations.

Sign of APOBEC editing, purifying selection, frameshift, and in-frame nonsense mutations in the microevolution of lumpy skin disease virus

The lumpy skin disease virus (LSDV), which mostly affects ruminants and causes huge-economic loss, was endemic in Africa, caused outbreaks in the Middle East, and was recently detected in Russia, Serbia, Greece, Bulgaria, Kazakhstan, China, Taiwan, Vietnam, Thailand, and India. However, the role of evolutionary drivers such as codon selection, negative/purifying selection, APOBEC editing, and genetic variations such as frameshift and in-frame nonsense mutations in the LSDVs, which cause outbreaks in cattle in various countries, are still largely unknown. In the present study, a frameshift mutation in LSDV035, LSDV019, LSDV134, and LSDV144 genes and in-frame non-sense mutations in LSDV026, LSDV086, LSDV087, LSDV114, LSDV130, LSDV131, LSDV145, LSDV154, LSDV155, LSDV057, and LSDV081 genes were revealed among different clusters. Based on the available complete genome sequences, the prototype wild-type cluster-1.2.1 virus has been found in other than Africa only in India, the wild-type cluster-1.2.2 virus found in Africa were spread outside Africa, and the recombinant viruses spreading only in Asia and Russia. Although LSD viruses circulating in different countries form a specific cluster, the viruses detected in each specific country are distinguished by frameshift and in-frame nonsense mutations. Furthermore, the present study has brought to light that the selection pressure for codons usage bias is mostly exerted by purifying selection, and this process is possibly caused by APOBEC editing. Overall, the present study sheds light on microevolutions in LSDV, expected to help in future studies towards disturbed ORFs, epidemiological diagnostics, attenuation/vaccine reverts, and predicting the evolutionary direction of LSDVs.

Plastid-localized xanthorhodopsin increases diatom biomass and ecosystem productivity in iron-limited surface oceans

Microbial rhodopsins are photoreceptor proteins that convert light into biological signals or energy. Proteins of the xanthorhodopsin family are common in eukaryotic photosynthetic plankton including diatoms. However, their biological role in these organisms remains elusive. Here we report on a xanthorhodopsin variant (FcR1) isolated from the polar diatom Fragilariopsis cylindrus. Applying a combination of biophysical, biochemical and reverse genetics approaches, we demonstrate that FcR1 is a plastid-localized proton pump which binds the chromophore retinal and is activated by green light. Enhanced growth of a Thalassiora pseudonana gain-of-function mutant expressing FcR1 under iron limitation shows that the xanthorhodopsin proton pump supports growth when chlorophyll-based photosynthesis is iron-limited. The abundance of xanthorhodopsin transcripts in natural diatom communities of the surface oceans is anticorrelated with the availability of dissolved iron. Thus, we propose that these proton pumps convey a fitness advantage in regions where phytoplankton growth is limited by the availability of dissolved iron.

Phylogeny and species diversity of the genus Helvella with emphasis on eighteen new species from China

Helvella is a widespread, frequently encountered fungal group appearing in forests, but the species diversity and molecular phylogeny of Helvella in China remains incompletely understood. In this work, we performed comprehensive phylogenetic analyses using multilocus sequence data. Six datasets were employed, including a five-locus concatenated dataset (ITS, nrLSU, tef1-α, rpb2, hsp), a two-locus concatenated dataset (ITS, nrLSU), and four single-locus datasets (ITS) that were divided based on the four different phylogenetic clades of Helvella recognized in this study. A total of I 946 sequences were used, of which 713 were newly generated, including 170 sequences of ITS, 174 sequences of nrLSU, 131 sequences of tef1-α, 107 sequences of rpb2 and 131 sequences of hsp. The phylogeny based on the five-locus concatenated dataset revealed that Helvellas. str. is monophyletic and four phylogenetic clades are clearly recognized, i.e., Acetabulum clade, Crispa clade, Elastica clade, and Lacunosa clade. A total of 24 lineages or subclades were recognized, II of which were new, the remaining 13 corresponding with previous studies. Chinese Helvella species are distributed in 22 lineages across four clades. Phylogenetic analyses based on the two-locus concatenated dataset and four single-locus datasets confirmed the presence of at least 93 phylogenetic species in China. Among them, 58 are identified as known species, including a species with a newly designated lectotype and epitype, 18 are newly described in this paper, and the remaining 17 taxa are putatively new to science but remain unnamed due to the paucity or absence of ascomatal materials. In addition, the Helvella species previously recorded in China are discussed. A list of 76 confirmed species, including newly proposed species, is provided. The occurrence of H. crispa and H. elastica are not confirmed although both are commonly recorded in China. Citation: Mao N, Xu YY, Zhang YX, Zhou H, Huang XB, Hou CL, Fan L (2023). Phylogeny and species diversity of the genus Helvella with emphasis on eighteen new species from China. Fungal Systematics and Evolution 12: 111-152. doi: 10.3114/fuse.2023.12.08.

Human immunodeficiency virus 1 5'-leader mutations in plasma viruses before and after the development of reverse transcriptase inhibitor-resistance mutations

Human immunodeficiency virus 1 (HIV-1) reverse transcriptase (RT) initiation depends on interaction between viral 5'-leader RNA, RT and host tRNA3Lys. Therefore, we sought to identify co-evolutionary changes between the 5'-leader and RT in viruses developing RT-inhibitor resistance mutations. We sequenced 5'-leader positions 37-356 of paired plasma virus samples from 29 individuals developing the nucleoside RT inhibitor (NRTI)-resistance mutation M184V, 19 developing a non-nucleoside RT inhibitor (NNRTI)-resistance mutation and 32 untreated controls. 5'-Leader variants were defined as positions where ≥20 % of next-generation sequencing (NGS) reads differed from the HXB2 sequence. Emergent mutations were defined as nucleotides undergoing a ≥4-fold change in proportion between baseline and follow-up. Mixtures were defined as positions containing ≥2 nucleotides each present in ≥20 % of NGS reads. Among 80 baseline sequences, 87 positions (27.2 %) contained a variant; 52 contained a mixture. Position 201 was the only position more likely to develop a mutation in the M184V (9/29 vs 0/32; P=0.0006) or NNRTI-resistance (4/19 vs 0/32; P=0.02; Fisher's exact test) groups than the control group. Mixtures at positions 200 and 201 occurred in 45.0 and 28.8 %, respectively, of baseline samples. Because of the high proportion of mixtures at these positions, we analysed 5'-leader mixture frequencies in two additional datasets: five publications reporting 294 dideoxyterminator clonal GenBank sequences from 42 individuals and six National Center for Biotechnology Information (NCBI) BioProjects reporting NGS datasets from 295 individuals. These analyses demonstrated position 200 and 201 mixtures at proportions similar to those in our samples and at frequencies several times higher than at all other 5'-leader positions. Although we did not convincingly document co-evolutionary changes between RT and 5'-leader sequences, we identified a novel phenomenon, wherein positions 200 and 201 immediately downstream of the HIV-1 primer binding site exhibited an extraordinarily high likelihood of containing a nucleotide mixture. Possible explanations for the high mixture rates are that these positions are particularly error-prone or provide a viral fitness advantage.

HIV-1 5'-Leader Mutations in Plasma Viruses Before and After the Development of Reverse Transcriptase Inhibitor-Resistance Mutations

Background

HIV-1 RT initiation depends on interaction between viral 5'-leader RNA, RT, and host tRNA3Lys. We therefore sought to identify co-evolutionary changes between the 5'-leader and RT in viruses developing RT-inhibitor resistance mutations.

Methods

We sequenced 5'-leader positions 37-356 of paired plasma virus samples from 29 individuals developing the NRTI-resistance mutation M184V, 19 developing an NNRTI-resistance mutation, and 32 untreated controls. 5'-leader variants were defined as positions where ≥20% of NGS reads differed from the HXB2 sequence. Emergent mutations were defined as nucleotides undergoing ≥4-fold change in proportion between baseline and follow-up. Mixtures were defined as positions containing ≥2 nucleotides each present in ≥20% of NGS reads.

Results

Among 80 baseline sequences, 87 positions (27.2%) contained a variant; 52 contained a mixture. Position 201 was the only position more likely to develop a mutation in the M184V (9/29 vs. 0/32; p=0.0006) or NNRTI-resistance (4/19 vs. 0/32; p=0.02; Fisher's Exact Test) groups than the control group. Mixtures at positions 200 and 201 occurred in 45.0% and 28.8%, respectively, of baseline samples. Because of the high proportion of mixtures at these positions, we analyzed 5'-leader mixture frequencies in two additional datasets: five publications reporting 294 dideoxyterminator clonal GenBank sequences from 42 individuals and six NCBI BioProjects reporting NGS datasets from 295 individuals. These analyses demonstrated position 200 and 201 mixtures at proportions similar to those in our samples and at frequencies several times higher than at all other 5'-leader positions.

Conclusions

Although we did not convincingly document co-evolutionary changes between RT and 5'-leader sequences, we identified a novel phenomenon, wherein positions 200 and 201, immediately downstream of the HIV-1 primer binding site exhibited an extraordinarily high likelihood of containing a nucleotide mixture. Possible explanations for the high mixture rates are that these positions are particularly error-prone or provide a viral fitness advantage.

Classification of Promoter Sequences from Human Genome

We have developed a new method for promoter sequence classification based on a genetic algorithm and the MAHDS sequence alignment method. We have created four classes of human promoters, combining 17,310 sequences out of the 29,598 present in the EPD database. We searched the human genome for potential promoter sequences (PPSs) using dynamic programming and position weight matrices representing each of the promoter sequence classes. A total of 3,065,317 potential promoter sequences were found. Only 1,241,206 of them were located in unannotated parts of the human genome. Every other PPS found intersected with either true promoters, transposable elements, or interspersed repeats. We found a strong intersection between PPSs and Alu elements as well as transcript start sites. The number of false positive PPSs is estimated to be 3 × 10-8 per nucleotide, which is several orders of magnitude lower than for any other promoter prediction method. The developed method can be used to search for PPSs in various eukaryotic genomes.

On seven undescribed leaf insect species revealed within the recent "Tree of Leaves" (Phasmatodea, Phylliidae)

With the recent advance in molecular phylogenetics focused on the leaf insects (Phasmatodea, Phylliidae), gaps in knowledge are beginning to be filled. Yet, shortcomings are also being highlighted, for instance, the unveiling of numerous undescribed phylliid species. Here, some of these taxa are described, including Phylliumiyadaonsp. nov. from Mindoro Island, Philippines; Phylliumsamarensesp. nov. from Samar Island, Philippines; Phylliumortizisp. nov. from Mindanao Island, Philippines; Pulchriphylliumheraclessp. nov. from Vietnam; Pulchriphylliumdelisleisp. nov. from South Kalimantan, Indonesia; and Pulchriphylliumbhaskaraisp. nov. from Java, Indonesia. Several additional specimens of these species together with a seventh species described herein, Pulchriphylliumanangusp. nov. from southwestern India, were incorporated into a newly constructed phylogenetic tree. Additionally, two taxa that were originally described as species, but in recent decades have been treated as subspecies, are elevated back to species status to reflect their unique morphology and geographic isolation, creating the following new combinations: Pulchriphylliumscythe (Gray, 1843) stat. rev., comb. nov. from Bangladesh and northeastern India, and Pulchriphylliumcrurifolium (Audinet-Serville, 1838) stat. rev., comb. nov. from the Seychelles islands. Lectotype specimens are also designated for Pulchriphylliumscythe (Gray, 1843) stat. rev., comb. nov. and Pulchriphylliumcrurifolium (Audinet-Serville, 1838) stat. rev., comb. nov. from original type material.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

Drug resistance prediction for Mycobacterium tuberculosis with reference graphs

Tuberculosis is a global pandemic disease with a rising burden of antimicrobial resistance. As a result, the World Health Organization (WHO) has a goal of enabling universal access to drug susceptibility testing (DST). Given the slowness of and infrastructure requirements for phenotypic DST, whole-genome sequencing, followed by genotype-based prediction of DST, now provides a route to achieving this. Since a central component of genotypic DST is to detect the presence of any known resistance-causing mutations, a natural approach is to use a reference graph that allows encoding of known variation. We have developed DrPRG (Drug resistance Prediction with Reference Graphs) using the bacterial reference graph method Pandora. First, we outline the construction of a Mycobacterium tuberculosis drug resistance reference graph. The graph is built from a global dataset of isolates with varying drug susceptibility profiles, thus capturing common and rare resistance- and susceptible-associated haplotypes. We benchmark DrPRG against the existing graph-based tool Mykrobe and the haplotype-based approach of TBProfiler using 44 709 and 138 publicly available Illumina and Nanopore samples with associated phenotypes. We find that DrPRG has significantly improved sensitivity and specificity for some drugs compared to these tools, with no significant decreases. It uses significantly less computational memory than both tools, and provides significantly faster runtimes, except when runtime is compared to Mykrobe with Nanopore data. We discover and discuss novel insights into resistance-conferring variation for M. tuberculosis - including deletion of genes katG and pncA - and suggest mutations that may warrant reclassification as associated with resistance.

Speciation across the Earth driven by global cooling in terrestrial orchids

Although climate change has been implicated as a major catalyst of diversification, its effects are thought to be inconsistent and much less pervasive than localized climate or the accumulation of species with time. Focused analyses of highly speciose clades are needed in order to disentangle the consequences of climate change, geography, and time. Here, we show that global cooling shapes the biodiversity of terrestrial orchids. Using a phylogeny of 1,475 species of Orchidoideae, the largest terrestrial orchid subfamily, we find that speciation rate is dependent on historic global cooling, not time, tropical distributions, elevation, variation in chromosome number, or other types of historic climate change. Relative to the gradual accumulation of species with time, models specifying speciation driven by historic global cooling are over 700 times more likely. Evidence ratios estimated for 212 other plant and animal groups reveal that terrestrial orchids represent one of the best-supported cases of temperature-spurred speciation yet reported. Employing >2.5 million georeferenced records, we find that global cooling drove contemporaneous diversification in each of the seven major orchid bioregions of the Earth. With current emphasis on understanding and predicting the immediate impacts of global warming, our study provides a clear case study of the long-term impacts of global climate change on biodiversity.

Cytoplasmic Polyadenylation Is an Ancestral Hallmark of Early Development in Animals

Differential regulation of gene expression has produced the astonishing diversity of life on Earth. Understanding the origin and evolution of mechanistic innovations for control of gene expression is therefore integral to evolutionary and developmental biology. Cytoplasmic polyadenylation is the biochemical extension of polyadenosine at the 3'-end of cytoplasmic mRNAs. This process regulates the translation of specific maternal transcripts and is mediated by the Cytoplasmic Polyadenylation Element-Binding Protein family (CPEBs). Genes that code for CPEBs are amongst a very few that are present in animals but missing in nonanimal lineages. Whether cytoplasmic polyadenylation is present in non-bilaterian animals (i.e., sponges, ctenophores, placozoans, and cnidarians) remains unknown. We have conducted phylogenetic analyses of CPEBs, and our results show that CPEB1 and CPEB2 subfamilies originated in the animal stem lineage. Our assessment of expression in the sea anemone, Nematostella vectensis (Cnidaria), and the comb jelly, Mnemiopsis leidyi (Ctenophora), demonstrates that maternal expression of CPEB1 and the catalytic subunit of the cytoplasmic polyadenylation machinery (GLD2) is an ancient feature that is conserved across animals. Furthermore, our measurements of poly(A)-tail elongation reveal that key targets of cytoplasmic polyadenylation are shared between vertebrates, cnidarians, and ctenophores, indicating that this mechanism orchestrates a regulatory network that is conserved throughout animal evolution. We postulate that cytoplasmic polyadenylation through CPEBs was a fundamental innovation that contributed to animal evolution from unicellular life.

Annelid functional genomics reveal the origins of bilaterian life cycles

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2-11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called 'head larvae' of other bilaterians13-17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

HMMerge: an ensemble method for multiple sequence alignment

Motivation

Despite advances in method development for multiple sequence alignment over the last several decades, the alignment of datasets exhibiting substantial sequence length heterogeneity, especially when the input sequences include very short sequences (either as a result of sequencing technologies or of large deletions during evolution) remains an inadequately solved problem.

Results

We present HMMerge, a method to compute an alignment of datasets exhibiting high sequence length heterogeneity, or to add short sequences into a given 'backbone' alignment. HMMerge builds on the technique from its predecessor alignment methods, UPP and WITCH, which build an ensemble of profile HMMs to represent the backbone alignment and add the remaining sequences into the backbone alignment using the ensemble. HMMerge differs from UPP and WITCH by building a new 'merged' HMM from the ensemble, and then using that merged HMM to align the query sequences. We show that HMMerge is competitive with WITCH, with an advantage over WITCH when adding very short sequences into backbone alignments.

Availability and implementation

HMMerge is freely available at https://github.com/MinhyukPark/HMMerge.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Genome Analysis of 10K SARS-COV-2 Sequences to Identify the Presence of Single-Nucleotide Polymorphisms

A new type of coronavirus was identified in Wuhan, China, in December 2019, which was named SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus-2). The high mutation rate of SARS-CoV2 makes it challenging to develop effective vaccines for all variants. Substitution is the most common type of mutation that occurs in SARS-CoV-2. This research was conducted to identify the genetic variability of mutations in SNP of SARS-CoV-2 and analyse the impact. About 15,000 sequences of SARS-CoV-2 were downloaded from GISAID, which were isolated from 33 different countries around the world from February 2020 to July 2021. Sequence analysis was done using the MAFFT and the Nextclade. The results of this study are expected to help identify conserved regions in SARS-CoV-2 which can be used as probes for the virus identification process and can be used as target areas in vaccine development. Furthermore the results showed that the most common variants were variants 20B, 20A, and 20I (Alpha), with a population percentage of 32.12%, 23.95% and 17.39% of the total population, respectively. Furthermore, SNPs were called in the samples using the SNP-sites and extracted using Excel. Of the 10,107 sequences of SARSCoV-2 studied, 154 SNPs were found with the highest number of SNPs in the spike, nsp3 and nucleocapsid genes. The ratio of the number of mutations to the most extensive sequence length was in the ORF8, ORF7a, and ORF7b genes with respective values of 0.537, 0.474, and 0.419.

Afforestation can lower microbial diversity and functionality in deep soil layers in a semiarid region

Afforestation is an effective approach to rehabilitate degraded ecosystems, but often depletes deep soil moisture. Presently, it is not known how an afforestation-induced decrease in moisture affects soil microbial community and functionality, hindering our ability to understand the sustainability of the rehabilitated ecosystems. To address this issue, we examined the impacts of 20 years of afforestation on soil bacterial community, co-occurrence pattern, and functionalities along vertical profile (0-500 cm depth) in a semiarid region of China's Loess Plateau. We showed that the effects of afforestation with a deep-rooted legume tree on cropland were greater in deep than that of in top layers, resulting in decreased bacterial beta diversity, more responsive bacterial taxa and functional groups, increased homogeneous selection, and decreased network robustness in deep soils (120-500 cm). Organic carbon and nitrogen decomposition rates and multifunctionality also significantly decreased by afforestation, and microbial carbon limitation significantly increased in deep soils. Moreover, changes in microbial community and functionality in deep layer was largely related to changes in soil moisture. Such negative impacts on deep soils should be fully considered for assessing afforestation's eco-environment effects and for the sustainability of ecosystems because deep soils have important influence on forest ecosystems in semiarid and arid climates.

New species of Sticta (lichenised Ascomycota, lobarioid Peltigeraceae) from Bolivia suggest a high level of endemism in the Central Andes

Six species of Sticta are described as new to science on the basis of material from Bolivia and supported by phylogenetic analysis of the fungal ITS barcoding marker. The species were resolved in all three of the clades (I, II, III) widespread and common in the Neotropics, as defined in an earlier study on the genus. Comparison with material from neighbouring countries (i.e. Colombia, Ecuador, Peru) suggests that these new species may be potentially endemic to the Bolivian Yungas ecoregion. For each species, a detailed morphological and anatomical description is given. Stictaamboroensis Ossowska, Kukwa, B. Moncada & Lücking is a medium-sized green-algal species with laminal to submarginal apothecia with hirsute margins and with light to dark brown lower tomentum. Stictaaymara Ossowska, Kukwa, B. Moncada, Flakus, Rodriguez-Flakus & Lücking is a comparatively small cyanobacterial taxon with Nostoc as photobiont, laminal, richly branched, aggregate isidia and a golden to chocolate-brown lower tomentum. The medium-sized, cyanobacterial S.bicellulata Ossowska, Kukwa, B. Moncada & Lücking has cyanobacterial photobiont, bicellular ascospores, apothecia with white to golden-brown hairs on the margins, K+ violet apothecial margin (ring around disc) and epihymenium and a white to dark brown lower tomentum. In contrast, the green-algal species, S.carrascoensis Ossowska, Kukwa, B. Moncada & Lücking is characterised by its large size, apothecia with dark brown hairs on the margins and a yellow medulla. The cyanobacterial S.catharinae Ossowska, B. Moncada, Kukwa, Flakus, Rodriguez-Flakus & Lücking forms stipitate thalli with Nostoc as photobiont, abundant, laminal to submarginal apothecia and a golden-brown lower tomentum. Finally, the cyanobacterial S.pseudoimpressula Ossowska, Kukwa, B. Moncada & Lücking produces laminal apothecia with an orange-yellow line of pruina along the margins which reacts K+ carmine-red. In addition to the six new Bolivian taxa, the cyanobacterial S.narinioana B. Moncada, Ossowska & Lücking is described as new from Colombia and it represents the closely-related sister species of the Bolivian S.aymara; it differs from the latter largely in the marginal instead of laminal isidia.

Detection of human pathogenic bacteria in rectal DNA samples from Zalophus californianus in the Gulf of California, Mexico

Human intrusions into undisturbed wildlife areas greatly contribute to the emergence of infectious diseases. To minimize the impacts of novel emerging infectious diseases (EIDs) on human health, a comprehensive understanding of the microbial species that reside within wildlife species is required. The Gulf of California (GoC) is an example of an undisturbed ecosystem. However, in recent decades, anthropogenic activities within the GoC have increased. Zalophus californianus has been proposed as the main sentinel species in the GoC; hence, an assessment of sea lion bacterial microbiota may reveal hidden risks for human health. We evaluated the presence of potential human pathogenic bacterial species from the gastrointestinal (GI) tracts of wild sea lions through a metabarcoding approach. To comprehensively evaluate this bacterial consortium, we considered the genetic information of six hypervariable regions of 16S rRNA. Potential human pathogenic bacteria were identified down to the species level by integrating the RDP and Pplacer classifier outputs. The combined genetic information from all analyzed regions suggests the presence of at least 44 human pathogenic bacterial species, including Shigella dysenteriae and Bacillus anthracis. Therefore, the risks of EIDs from this area should be not underestimated.

Neogastropod (Mollusca, Gastropoda) phylogeny: A step forward with mitogenomes

The Neogastropoda (Mollusca, Gastropoda) encompass more than 15,000 described species of marine predators, including several model organisms in toxinology, embryology and physiology. However, their phylogenetic relationships remain mostly unresolved and their classification unstable. We took advantage of the many mitogenomes published in GenBank to produce a new molecular phylogeny of the neogastropods. We completed the taxon sampling by using an in-house bioinformatic pipeline to retrieve mitochondrial genes from 13 transcriptomes, corresponding to five families not represented in GenBank, for a final dataset of 113 taxa. Because mitogenomic data are prone to reconstruction artefacts, eight different evolutionary models were applied to reconstruct phylogenetic trees with IQTREE, RAxML and MrBayes. If the over-parametrization of some models produced trees with aberrant internal long branches, the global topology of the trees remained stable over models and softwares, and several relationships were revealed or found supported here for the first time. However, even if our dataset encompasses 60% of the valid families of neogastropods, some key taxa are missing and should be added in the future before proposing a revision of the classification of the neogastropods. Our study also demonstrates that even complex models struggle to satisfactorily handle the evolutionary history of mitogenomes, still leading to long-branch attractions in phylogenetic trees. Other approaches, such as reduced-genome strategies, must be envisaged to fully resolve the neogastropod phylogeny.

Phylogeny and evolution of morphological structures in a highly diverse lineage of fruiting-body-forming amoebae, order Trichiales (Myxomycetes, Amoebozoa)

Early phylogenetic studies refuted most previous assumptions concerning the evolution of the morphological traits in the fruiting bodies of the order Trichiales and did not detect discernible evolutionary patterns, yet they were based on a limited number of species. We infer a new Trichiales phylogeny based on three independently inherited genetic regions (nuclear and mitochondrial), with a fair taxonomic sampling encompassing its broad diversity. Besides, we study the evolutionary history of some key morphological characters. According to the new phylogeny, most fruiting body traits in Trichiales systematics do not represent exclusive synapomorphies or autapomorphies for most monophyletic groups. Instead, the evolution of the features derived from the peridium, stalk, capillitium, and spores showed intricate patterns, and character state transitions occurred rather within- than between clades. Thus, we should consider other evolutionary scenarios instead of assuming the homology of some characters. According to these results, we propose a new classification of Trichiales, including the creation of a new genus, Gulielmina, the resurrection of the family Dictydiaethaliaceae and the genus Ophiotheca, and the proporsal of 13 new combinations for species of the genera Arcyria (1), Hemitrichia (2), Ophiotheca (2), Oligonema (4), Gulielmina (3), and Perichaena (1).

Metagenomic Analysis Using Phylogenetic Placement-A Review of the First Decade

Phylogenetic placement refers to a family of tools and methods to analyze, visualize, and interpret the tsunami of metagenomic sequencing data generated by high-throughput sequencing. Compared to alternative (e. g., similarity-based) methods, it puts metabarcoding sequences into a phylogenetic context using a set of known reference sequences and taking evolutionary history into account. Thereby, one can increase the accuracy of metagenomic surveys and eliminate the requirement for having exact or close matches with existing sequence databases. Phylogenetic placement constitutes a valuable analysis tool per se, but also entails a plethora of downstream tools to interpret its results. A common use case is to analyze species communities obtained from metagenomic sequencing, for example via taxonomic assignment, diversity quantification, sample comparison, and identification of correlations with environmental variables. In this review, we provide an overview over the methods developed during the first 10 years. In particular, the goals of this review are 1) to motivate the usage of phylogenetic placement and illustrate some of its use cases, 2) to outline the full workflow, from raw sequences to publishable figures, including best practices, 3) to introduce the most common tools and methods and their capabilities, 4) to point out common placement pitfalls and misconceptions, 5) to showcase typical placement-based analyses, and how they can help to analyze, visualize, and interpret phylogenetic placement data.

Human-to-dog transmission of SARS-CoV-2, Colombia

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of the current COVID-19 pandemic, has evolved to have a wide range of hosts, including non-human primates, wild and domestic animals. The ACE2 protein has a high level of conservation and is the common receptor invertebrate species for a viral infection to occur; this receptor could give rise to anthroponotic events. This article describes the first event of symptomatic transmission in Latin America from a human to a dog by the B.1.625 lineage of SARS-CoV-2. We found 21 shared mutations in the complete genomes of viral sequences from owners and dogs. Further phylogenetic and molecular analysis showed that 100% co-localization of the clade helps to understand human-animal transmission. Prediction of the Spike protein structure of the sequenced virus and docking analyzes showed that the E484K mutation in the receptor-binding domain (RBD) could contribute to the viral affinity of dACE2. Therefore, close contact between SARS-CoV-2-infected humans and pets should be avoided to prevent the emergence of novel mutations of public health importance from anthroponotic events.

Functional and Phylogenetic Characterization of Bacteria in Bovine Rumen Using Fractionation of Ruminal Fluid

Cattle productivity depends on our ability to fully understand and manipulate the fermentation process of plant material that occurs in the bovine rumen, which ultimately leads to the improvement of animal health and increased productivity with a reduction in environmental impact. An essential step in this direction is the phylogenetic and functional characterization of the microbial species composing the ruminal microbiota. To address this challenge, we separated a ruminal fluid sample by size and density using a sucrose density gradient. We used the full sample and the smallest fraction (5%), allowing the enrichment of bacteria, to assemble metagenome-assembled genomes (MAGs). We obtained a total of 16 bacterial genomes, 15 of these enriched in the smallest fraction of the gradient. According to the recently proposed Genome Taxonomy Database (GTDB) taxonomy, these MAGs belong to Bacteroidota, Firmicutes_A, Firmicutes, Proteobacteria, and Spirochaetota phyla. Fifteen MAGs were novel at the species level and four at the genus level. The functional characterization of these MAGs suggests differences from what is currently known from the genomic potential of well-characterized members from this complex environment. Species of the phyla Bacteroidota and Spirochaetota show the potential for hydrolysis of complex polysaccharides in the plant cell wall and toward the production of B-complex vitamins and protein degradation in the rumen. Conversely, the MAGs belonging to Firmicutes and Alphaproteobacteria showed a reduction in several metabolic pathways; however, they have genes for lactate fermentation and the presence of hydrolases and esterases related to chitin degradation. Our results demonstrate that the separation of the rumen microbial community by size and density reduced the complexity of the ruminal fluid sample and enriched some poorly characterized ruminal bacteria allowing exploration of their genomic potential and their functional role in the rumen ecosystem.

Protist Predation Influences the Temperature Response of Bacterial Communities

Temperature strongly influences microbial community structure and function, in turn contributing to global carbon cycling that can fuel further warming. Recent studies suggest that biotic interactions among microbes may play an important role in determining the temperature responses of these communities. However, how predation regulates these microbiomes under future climates is still poorly understood. Here, we assess whether predation by a key global bacterial consumer-protists-influences the temperature response of the community structure and function of a freshwater microbiome. To do so, we exposed microbial communities to two cosmopolitan protist species-Tetrahymena thermophila and Colpidium sp.-at two different temperatures, in a month-long microcosm experiment. While microbial biomass and respiration increased with temperature due to community shifts, these responses changed over time and in the presence of protists. Protists influenced microbial biomass and respiration rate through direct and indirect effects on bacterial community structure, and predator presence actually reduced microbial respiration at elevated temperature. Indicator species analyses showed that these predator effects were mostly determined by phylum-specific bacterial responses to protist density and cell size. Our study supports previous findings that temperature is an important driver of microbial communities but also demonstrates that the presence of a large predator can mediate these responses to warming.

Developments in Algorithms for Sequence Alignment: A Review

The continuous development of sequencing technologies has enabled researchers to obtain large amounts of biological sequence data, and this has resulted in increasing demands for software that can perform sequence alignment fast and accurately. A number of algorithms and tools for sequence alignment have been designed to meet the various needs of biologists. Here, the ideas that prevail in the research of sequence alignment and some quality estimation methods for multiple sequence alignment tools are summarized.

Application of the MAHDS Method for Multiple Alignment of Highly Diverged Amino Acid Sequences

The aim of this work was to compare the multiple alignment methods MAHDS, T-Coffee, MUSCLE, Clustal Omega, Kalign, MAFFT, and PRANK in their ability to align highly divergent amino acid sequences. To accomplish this, we created test amino acid sequences with an average number of substitutions per amino acid (x) from 0.6 to 5.6, a total of 81 sets. Comparison of the performance of sequence alignments constructed by MAHDS and previously developed algorithms using the CS and Z score criteria and the benchmark alignment database (BAliBASE) indicated that, although the quality of the alignments built with MAHDS was somewhat lower than that of the other algorithms, it was compensated by greater statistical significance. MAHDS could construct statistically significant alignments of artificial sequences with x ≤ 4.8, whereas the other algorithms (T-Coffee, MUSCLE, Clustal Omega, Kalign, MAFFT, and PRANK) could not perform that at x > 2.4. The application of MAHDS to align 21 families of highly diverged proteins (identity < 20%) from Pfam and HOMSTRAD databases showed that it could calculate statistically significant alignments in cases when the other methods failed. Thus, MAHDS could be used to construct statistically significant multiple alignments of highly divergent protein sequences, which accumulated multiple mutations during evolution.

Paxilloboletus gen. nov., a new lamellate bolete genus from tropical Africa

This study presents Paxilloboletus gen. nov., a new lamellate bolete genus represented by two tropical African species, Paxilloboletus africanus sp. nov. and Paxilloboletus latisporus sp. nov. Although the new taxa strongly resemble Paxillus (Paxillaceae), they lack clamp connections and form a separate generic clade within the Boletaceae phylogeny. The new species are lookalikes, morphologically only separable by their spore morphology. Descriptions and illustrations of the new genus and new species are given, as well as comments on ecology, distribution, and morphological differences with other gilled Boletaceae.

Intestinal microbial communities and Holdemanella isolated from HIV+/- men who have sex with men increase frequencies of lamina propria CCR5+ CD4+ T cells

Men who have sex with men (MSM), regardless of HIV infection status, have an intestinal microbiome that is compositionally distinct from men who have sex with women (MSW) and women. We recently showed HIV-negative MSM have elevated levels of intestinal CD4⁺ T cells expressing CCR5, a critical co-receptor for HIV. Whether elevated expression of CCR5 is driven by the altered gut microbiome composition in MSM has not been explored. Here we used in vitro stimulation of gut Lamina Propria Mononuclear Cells (LPMCs) with whole intact microbial cells isolated from stool to demonstrate that fecal bacterial communities (FBCs) from HIV-positive/negative MSM induced higher frequencies of CCR5⁺ CD4⁺ T cells compared to FBCs from HIV-negative MSW and women. To identify potential microbial drivers, we related the frequency of CCR5⁺ CD4⁺ T cells to the abundance of individual microbial taxa in rectal biopsy of HIV-positive/negative MSM and controls, and Holdemanella biformis was strongly associated with increased frequency of CCR5⁺ CD4⁺ T cells. We used in vitro stimulation of gut LPMCs with the type strain of H. biformis, a second strain of H.biformis and an isolate of the closely related Holdemanella porci , cultured from either a HIV-positive or a HIV-negative MSM stool. H. porci elevated the frequency of both CCR5⁺ CD4⁺ T cells and the ratio of TNF-α/IL-10 Genomic comparisons of the 3 Holdemanella isolates revealed unique cell wall and capsular components, which may be responsible for their differences in immunogenicity. These findings describe a novel mechanism potentially linking intestinal dysbiosis in MSM to HIV transmission and mucosal pathogenesis.

An accurate method for identifying recent recombinants from unaligned sequences

Motivation

Recombination is a fundamental process in molecular evolution, and the identification of recombinant sequences is thus of major interest. However, current methods for detecting recombinants are primarily designed for aligned sequences. Thus they struggle with analyses of highly diverse genes, such as the var genes of the malaria parasite Plasmodium falciparum, which are known to diversify primarily through recombination.

Results

We introduce an algorithm to detect recent recombinant sequences from a dataset without a full multiple alignment. Our algorithm can handle thousands of gene-length sequences without the need for a reference panel. We demonstrate the accuracy of our algorithm through extensive numerical simulations; in particular, it maintains its effectiveness in the presence of insertions and deletions. We apply our algorithm to a dataset of 17,335 DBLα types in var genes from Ghana, observing that sequences belonging to the same ups group or domain subclass recombine amongst themselves more frequently, and that non-recombinant DBLα types are more conserved than recombinant ones.

Availability

Source code is freely available at https://github.com/qianfeng2/detREC_program.

Supplementary information

Supplementary data are available at Bioinformatics online.

Species diversity, phylogeny, endemism and geography of the truffle genus Tuber in China based on morphological and molecular data

The genus Tuber (Tuberaceae, Pezizales) is an important fungal group of Ascomycota both economically and ecologically. However, the species diversity, phylogenetic relationships, and geographic distribution of Tuber species in China remains poorly understood, primarily because descriptions of many new species relied heavily on morphological features with molecular data either not sought or ignored. The misapplication of European and North American names further added to confusion regarding the taxonomy of Chinese Tuber species. In this study, we examined more than 1000 specimens from China, and performed a comprehensive phylogenetic analysis for Chinese Tuber species using ITS sequences and multilocus sequence data. To infer the phylogeny of Chinese Tuber spp., 11 molecular datasets were assembled, including a concatenated internal transcribed spacers of the nuc rDNA (ITS), nuc rDNA 28S subunit (LSU), translation elongation factor 1-alpha (tef1-α), and RNA polymerase II subunit (rpb2) dataset as well as 10 ITS datasets (totally including 1435 sequences from 828 collections with 597 newly generated sequences, and 168 sequences from the types of 63 species). Our phylogenetic tree based on a concatenated multilocus dataset revealed that all Chinese Tuber species nested in nine phylogenetic clades (phylogroups), including Aestivum, Excavatum, Latisporum, Macrosporum, Maculatum, Melanosporum, Puberulum, Rufum and Turmericum. Of these, five phylogroups (Macrosporum, Maculatum, Melanosporum, Puberulum and Rufum) are shared across the continents of Asia, Europe and North America; two phylogroups (Aestivum and Excavatum) are shared by Europe and Asia; and the phylogroups Turmericum and Latisporum are endemic only to Asia. Phylogenetic trees based on 10 ITS datasets confirmed the presence of at least 82 phylogenetic species in China. Of these, 53 are identified as known species, including three new records for China, and 25 species are identified as new to science. Of the new species, nine are described and illustrated in this paper, and the others remain unnamed due to the paucity or absence of ascomatal materials. Accordingly, the confirmed, excluded and doubtful Tuber species in China are discussed. Tuber species showed high endemism. Of the 82 phylogenetic species found in China, 68 species occur only in China, six species are also found in other regions in Asia, and only eight species (T. anniae, T. excelsum-reticulatum, T. formosanum, T. maculatum, T. wenchuanense, Tuber sp. CHN-3, Tuber sp. CHN-10 and Tuber sp. CHN-11) are shared with other continents. Most Tuber species have a small and limited distribution in China, but a few, such as T. formosanum and T. parvomurphium, are widely distributed across China. Some phylogenetically closely related species, such as T. liaotongense and T. subglobosum, as well as T. xuanhuaense and T. lijiangense, show a pattern of allopatric distribution.

MAGUS+eHMMs: improved multiple sequence alignment accuracy for fragmentary sequences

SUMMARY

Multiple sequence alignment is an initial step in many bioinformatics pipelines, including phylogeny estimation, protein structure prediction and taxonomic identification of reads produced in amplicon or metagenomic datasets, etc. Yet, alignment estimation is challenging on datasets that exhibit substantial sequence length heterogeneity, and especially when the datasets have fragmentary sequences as a result of including reads or contigs generated by next-generation sequencing technologies. Here, we examine techniques that have been developed to improve alignment estimation when datasets contain substantial numbers of fragmentary sequences. We find that MAGUS, a recently developed MSA method, is fairly robust to fragmentary sequences under many conditions, and that using a two-stage approach where MAGUS is used to align selected 'backbone sequences' and the remaining sequences are added into the alignment using ensembles of Hidden Markov Models further improves alignment accuracy. The combination of MAGUS with the ensemble of eHMMs (i.e. MAGUS+eHMMs) clearly improves on UPP, the previous leading method for aligning datasets with high levels of fragmentation.

AVAILABILITY AND IMPLEMENTATION

UPP is available on https://github.com/smirarab/sepp, and MAGUS is available on https://github.com/vlasmirnov/MAGUS. MAGUS+eHMMs can be performed by running MAGUS to obtain the backbone alignment, and then using the backbone alignment as an input to UPP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.