Lineage replacement and evolution captured by the United Kingdom Covid Infection Survey

The Office for National Statistics COVID-19 Infection Survey is a large household-based surveillance study based in the United Kingdom. Here, we report on the epidemiological and evolutionary dynamics of SARS-CoV-2 determined by analysing sequenced samples collected up until 13th November 2021. We observed four distinct sweeps or partial-sweeps, by lineages B.1.177, B.1.1.7/Alpha, B.1.617.2/Delta, and finally AY.4.2, a sublineage of B.1.617.2, with each sweeping lineage having a distinct growth advantage compared to their predecessors. Evolution was characterised by steady rates of evolution and increasing diversity within lineages, but with step increases in divergence associated with each sweeping major lineage, leading to a faster overall rate of evolution and fluctuating levels of diversity. These observations highlight the value of viral sequencing integrated into community surveillance studies to monitor the viral epidemiology and evolution of SARS-CoV-2, and potentially other pathogens, particularly as routine PCR testing is phased out or in settings where large-scale sequencing is not feasible.

Genome of the parasitoid wasp Dinocampus coccinellae reveals extensive duplications, accelerated evolution, and independent origins of thelytokous parthenogeny and solitary behavior

Dinocampus coccinellae (Hymenoptera: Braconidae) is a generalist parasitoid wasp that parasitizes >50 species of predatory lady beetles (Coleoptera: Coccinellidae), with thelytokous parthenogeny as its primary mode of reproduction. Here we present the first high quality genome of D. coccinellae using a combination of short and long read sequencing technologies, followed by assembly and scaffolding of chromosomal segments using Chicago+HiC technologies. We also present a first-pass ab initio and a reference-based genome annotation, and resolve timings of divergence and evolution of (1) solitary behavior vs eusociality, (2) arrhenotokous vs thelytokous parthenogenesis, and (3) rates of gene loss and gain among Hymenopteran lineages. Our study finds (1) at least two independent origins of eusociality and solitary behavior among Hymenoptera, (2) two independent origins of thelytokous parthenogenesis from ancestral arrhenotoky, and (3) accelerated rates of gene duplications, loss, and gain along the lineages leading to D. coccinellae. Our work both affirms the ancient divergence of Braconid wasps from ancestral Hymenopterans and accelerated rates of evolution in response to adaptations to novel hosts, including polyDNA viral co-evolution.

Six Newly Sequenced Chloroplast Genomes From Trentepohliales: The Inflated Genomes, Alternative Genetic Code and Dynamic Evolution

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26–36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata, which does have an inverted repeat. The duplicated trnD-GTC, petD, and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales.

KARAKTERISASI BERBASIS MARKA MOLEKULER ITS2 TERHADAP SUB-SPESIES KOMPLEKS Anopheles vagus vagus DAN Anopheles vagus limosus

The presence of intraspecies variations of An. vagus later categorized as the subspecies of An. vagus vagus and An. vagus limosus, could be an obstacle to the identification process, which is an important step for malaria vector’s competence characterization. Based on morphological identification, those subspecies could be distinguished by the presences of pale scales in prehumeral and pale bands in proboscis. The objective of this research was to compare subspecies complexes of An. vagus morphologically and molecularly using Internal Transcribed Spacer 2 (ITS2). Anopheles samples were collected from Bangsring, Banyuwangi. Their phylogenetic tree was constructed by using NJ method based on their ITS2 sequences. BLAST result showed that An. vagus vagus and An. vagus limosus were similar to An. vagus FJ654649.1 from East Java Indonesia and East Timor based on its 99% homology and their molecular distance. The Neighbour Joining (NJ) tree grouped those subspecies in one clade with a boostrap value of 82%. This subspeciation might be due to the different rates of evolution. ITS2 sequences of An. vagus vagus and An. vagus limosus were submitted to GenBank with the accession number of MW314227.1 and MW319822.1, respectively. Kemunculan variasi intraspesies An. vagus yang kemudian dikategorikan sebagai subspesies An. vagus vagus dan An. vagus limosus menjadi kendala dalam proses identifikasi yang merupakan langkah penting dalam menentukan kompetensi vektor malaria. Berdasarkan karakter morfologi, subspesies tersebut dibedakan dengan adanya sisik pucat pada bagian prehumeral dan pita pucat pada probosis. Penelitian ini bertujuan untuk membandingkan subspesies An. vagus secara morfologis dan molekuler menggunakan Internal Transcribed Spacer 2 (ITS2). Nyamuk Anopheles didapatkan dari Bangsring, Banyuwangi. Konstruksi pohon filogeni dilakukan berdasarkan sekuen ITS2 yang dianalisis menggunakan metode NJ. Hasil BLAST menunjukkan, ITS2 An. vagus vagus dan An. vagus limosus memiliki tingkat homologi 99% dan jarak evolusi molekuler terendah dengan An. vagus FJ654649.1 dari Jawa Timur Indonesia dan Timor Timur. Pohon NJ mengelompokkan subspesies tersebut dalam satu klade dengan nilai boostrap 82%. Hal ini dapat terjadi karena perbedaan kecepatan evolusi yang memungkinkan terjadinya subspesiasi. Urutan basa ITS2 dari An. vagus vagus dan An. vagus limosus telah didaftarkan ke GenBank dengan nomor aksesi MW314227.1 dan MW319822.1.

High rates of evolution preceded shifts to sex-biased gene expression in Leucadendron, the most sexually dimorphic angiosperms

Differences between males and females are usually more subtle in dioecious plants than animals, but strong sexual dimorphism has evolved convergently in the South African Cape plant genus Leucadendron. Such sexual dimorphism in leaf size is expected largely to be due to differential gene expression between the sexes. We compared patterns of gene expression in leaves among ten Leucadendron species across the genus. Surprisingly, we found no positive association between sexual dimorphism in morphology and the number or the percentage of sex-biased genes. Sex bias in most sex-biased genes evolved recently and was species-specific. We compared rates of evolutionary change in expression for genes that were sex-biased in one species but unbiased in others and found that sex-biased genes evolved faster in expression than un-biased genes. This greater rate of expression evolution of sex-biased genes, also documented in animals, might suggest the possible role of sexual selection in the evolution of gene expression. However, our comparative analysis clearly indicates that the more rapid rate of expression evolution of sex-biased genes predated the origin of bias, and shifts towards bias were depleted in signatures of adaptation. Our results are thus more consistent with the view that sex bias is simply freer to evolve in genes less subject to constraints in expression level.

Heterogeneity in conservation of multifunctional partner enzymes with meiotic importance, CDK2 kinase and BRCA1 ubiquitin ligase

The evolution of proteins can be accompanied by changes not only to their amino acid sequences, but also their structural and spatial molecular organization. Comparison of the protein conservation within different taxonomic groups (multifunctional, or highly specific) allows to clarify their specificity and the direction of evolution. Two multifunctional enzymes, cyclin-dependent kinase 2 (CDK2) and BRCA1 ubiquitin ligase, that are partners in some mitotic and meiotic processes were investigated in the present work. Two research methods, bioinformatics and immunocytochemical, were combined to examine the conservation levels of the two enzymes. It has been established that CDK2 is a highly conserved protein in different taxonomic lineages of the eukaryotic tree. Immunocytochemically, a conserved CDK2 pattern was revealed in the meiotic autosomes of five rodent species and partially in domestic turkey and clawed frog. Nevertheless, variable CDK2 distribution was detected at the unsynapsed segments of the rodent X chromosomes. BRCA1 was shown to be highly conserved only within certain mammalian taxa. It was also noted that in those rodent nuclei, where BRCA1 specifically binds to antigens, asynaptic regions of sex chromosomes were positive. BRCA1 staining was not always accompanied by specific binding, and a high nonspecificity in the nucleoplasm was observed. Thus, the studies revealed different conservation of the two enzymes at the level of protein structure as well as at the level of chromosome behavior. This suggests variable rates of evolution due to both size and configuration of the protein molecules and their multifunctionality.

Rapid Initial Morphospace Expansion and Delayed Morphological Disparity Peak in the First 100 Million Years of the Archosauromorph Evolutionary Radiation

Adaptive radiations have played a major role in generating modern and deep-time biodiversity. The Triassic radiation of the Archosauromorpha was one of the most spectacular vertebrate radiations, giving rise to many highly ecomorphologically varied lineages—including the dinosaurs, pterosaurs, and stem-crocodylians—that dominated the larger-bodied land fauna for the following 150 Ma, and ultimately gave rise to today’s > 10,000 species of birds and crocodylians. This radiation provides an outstanding testbed for hypotheses relating to adaptive radiations more broadly. Recent studies have started to characterize the tempo and mode of the archosauromorph early adaptive radiation, indicating very high initial rates of evolution, non-competitive niche-filling processes, and previously unrecognized morphological disparity even among non-crown taxa. However, these analyses rested primarily either on discrete characters or on geometric morphometrics of the cranium only, or even failed to fully include phylogenetic information. Here we expand previous 2D geometric morphometric cranial datasets to include new taxa and reconstructions, and create an analogous dataset of the pelvis, thereby allowing comparison of anatomical regions and the transition from “sprawling” to “upright” posture to be examined. We estimated morphological disparity and evolutionary rates through time. All sampled clades showed a delayed disparity peak for sum of variances and average nearest neighbor distances in both the cranium and pelvis, with disparity likely not saturated by the end of the studied time span (Late Jurassic); this contrasts with smaller radiations, but lends weight to similar results for large, ecomorphologically-varied groups. We find lower variations in pelvic than cranial disparity among Triassic-Jurassic archosaurs, which may be related to greater morphofunctional constraints on the pelvis. Contrasting with some previous work, but also confirming some previous findings during adaptive radiations, we find relatively widespread evidence of correlation between sampled diversity and disparity, especially at the largest phylogenetic scales and using average displacement rather than sum of variances as disparity metric; this also demonstrates the importance of comparing disparity metrics, and the importance of phylogenetic scale. Stem and crown archosauromorphs show a morphological diversification of both the cranium and pelvis with higher initial rates (Permian–Middle Triassic and at the base of major clades) followed by lower rates once diversification into niches has occurred (Late Triassic–Jurassic), indicating an “early burst” pattern sensu lato. Our results provide a more detailed and comprehensive picture of the early archosauromorph radiation and have significant bearing on the understanding of deep-time adaptive radiations more broadly, indicating widespread patterns of delayed disparity peaks, initial correlation of diversity and disparity, and evolutionary early bursts.

Evolution of Methyltransferase-Like (METTL) Proteins in Metazoa: A Complex Gene Family Involved in Epitranscriptomic Regulation and Other Epigenetic Processes

The methyltransferase-like (METTL) proteins constitute a family of seven-beta-strand methyltransferases with S-adenosyl methionine-binding domains that modify DNA, RNA, and proteins. Methylation by METTL proteins contributes to the epigenetic, and in the case of RNA modifications, epitranscriptomic regulation of a variety of biological processes. Despite their functional importance, most investigations of the substrates and functions of METTLs within metazoans have been restricted to model vertebrate taxa. In the present work, we explore the evolutionary mechanisms driving the diversification and functional differentiation of 33 individual METTL proteins across Metazoa. Our results show that METTLs are nearly ubiquitous across the animal kingdom, with most having arisen early in metazoan evolution (i.e., occur in basal metazoan phyla). Individual METTL lineages each originated from single independent ancestors, constituting monophyletic clades, which suggests that each METTL was subject to strong selective constraints driving its structural and/or functional specialization. Interestingly, a similar process did not extend to the differentiation of nucleoside-modifying and protein-modifying METTLs (i.e., each METTL type did not form a unique monophyletic clade). The members of these two types of METTLs also exhibited differences in their rates of evolution. Overall, we provide evidence that the long-term evolution of METTL family members was driven by strong purifying selection, which in combination with adaptive selection episodes, led to the functional specialization of individual METTL lineages. This work contributes useful information regarding the evolution of a gene family that fulfills a variety of epigenetic functions, and can have profound influences on molecular processes and phenotypic traits.

Prevalence and Molecular Characterization of Human Bocavirus Detected in Croatian Children with Respiratory Infection

Human bocavirus (HBoV) 1 is considered an important respiratory pathogen, while the role of HBoV2-4 in clinical disease remains somewhat controversial. Since, they are characterized by a rapid evolution, worldwide surveillance of HBoVs’ genetics is necessary. This study explored the prevalence of HBoV genotypes in pediatric patients with respiratory tract infection in Croatia and studied their phylogeny. Using multiplex PCR for 15 respiratory viruses, we investigated 957 respiratory samples of children up to 18 years of age with respiratory tract infection obtained from May 2017 to March 2021 at two different hospitals in Croatia. Amplification of HBoV near-complete genome or three overlapping fragments was performed, sequenced, and their phylogenetic inferences constructed. HBoV was detected in 7.6% children with a median age of 1.36 years. Co-infection was observed in 82.2% samples. Sequencing was successfully performed on 29 HBoV positive samples, and all belonged to HBoV1. Croatian HBoV1 sequences are closely related to strains isolated worldwide, and no phylogenetic grouping based on mono- or co-infection cases or year of isolation was observed. Calculated rates of evolution for HBoV1 were 10−4 and 10−5 substitutions per site and year. Recombination was not detected among sequences from this study.

Genomic Epidemiology and Evolution of Duck Hepatitis A Virus

Duck hepatitis A virus (DHAV), an avian picornavirus, causes high-mortality acute disease in ducklings. Among the three serotypes, DHAV-1 is globally distributed, whereas DHAV-2 and DHAV-3 serotypes are chiefly restricted to Southeast Asia. In this study, we analyzed the genomic evolution of DHAV-1 strains using extant GenBank records and genomic sequences of 10 DHAV-1 strains originating from a large disease outbreak in 2004–2005, in Hungary. Recombination analysis revealed intragenotype recombination within DHAV-1 as well as intergenotype recombination events involving DHAV-1 and DHAV-3 strains. The intergenotype recombination occurred in the VP0 region. Diversifying selection seems to act at sites of certain genomic regions. Calculations estimated slightly lower rates of evolution of DHAV-1 (mean rates for individual protein coding regions, 5.6286 × 10−4 to 1.1147 × 10−3 substitutions per site per year) compared to other picornaviruses. The observed evolutionary mechanisms indicate that whole-genome-based analysis of DHAV strains is needed to better understand the emergence of novel strains and their geographical dispersal.