Divergent Lineages Research Articles

Patterns of morphological diversification are influenced by dietary evolution in a highly species-rich lizard radiation

The diversification of lineages is facilitated or constrained by the simultaneous evolution of multiple components of the phenotype that interact with each other during the course of speciation. When evolutionary radiations are adaptive, lineages proliferate via the emergence of multiple phenotypic optima that underlie diversification of species across multiple ecological niches. When radiations are non-adaptive, lineage proliferation unfolds constrained by similar (or nearly identical) correlations among traits that keep phenotypic and ecological diversity across newly emerging species within a single optimum. Nature offers very few opportunities where both types of diversification occur between closely related and highly diverse lineages. The Liolaemidae family of South American lizards offers unique such opportunities given two speciose lineages that have rapidly proliferated via adaptive (Liolaemus) and non-adaptive (Phymaturus) radiations. We analyze body shape in lizards in association with type of diet (herbivory, omnivory or carnivory). In these lizards, diet types have been suggested to be linked to body size. Our results confirm this hypothesis, with three body size optima tightly linked to all three diet types when radiation is both adaptive and non-adaptive. Diet reconstruction along their evolutionary history showed that the common ancestor of Liolaemidae was likely omnivorous, which is matched by ancestral reconstruction of body size. Phylogenetic PCA revealed that herbivorous species generally have more differentiated body shape than insectivores and omnivores. Herbivorous species have evolved larger heads, shorter hindlimbs and a small difference between forelimb and hindlimb length. In contrast, omnivores and insectivores have smaller heads and longer hindlimbs. Collectively, trophic niche plays an important role in defining body shape and size across species within lineages, and the patterns of trait–ecology correlations remain consistent when lineages have diversified via adaptive and non-adaptive radiation.

High prevalence of haemosporidian parasites in Eurasian jays

Avian haemosporidians are vector-borne parasites, infecting a great variety of birds. The order Passeriformes has the highest average infection probability; nevertheless, some common species of Passeriformes have been rather poorly studied. We investigated haemosporidians in one such species, the Eurasian jay Garrulus glandarius (Corvidae), from a forest population in Hesse, Central Germany. All individuals were infected with at least one haemosporidian genus (overall prevalence: 100%). The most common infection pattern was a mixed Haemoproteus and Leucocytozoon infection, whereas no Plasmodium infection was detected. Results on lineage diversity indicate a rather pronounced host-specificity of Haemoproteus and Leucocytozoon lineages infecting birds of the family Corvidae.

Short-wavelength-sensitive 1 (SWS1) opsin gene duplications and parallel visual pigment tuning support ultraviolet communication in damselfishes (Pomacentridae).

Damselfishes (Pomacentridae) are one of the most behaviourally diverse, colourful and species-rich reef fish families. One remarkable characteristic of damselfishes is their communication in ultraviolet (UV) light. Not only are they sensitive to UV, they are also prone to have UV-reflective colours and patterns enabling social signalling. Using more than 50 species, we aimed to uncover the evolutionary history of UV colour and UV vision in damselfishes. All damselfishes had UV-transmitting lenses, expressed the UV-sensitive SWS1 opsin gene, and most displayed UV-reflective patterns and colours. We find evidence for several tuning events across the radiation, and while SWS1 gene duplications are generally very rare among teleosts, our phylogenetic reconstructions uncovered two independent duplication events: one close to the base of the most species-rich clade in the subfamily Pomacentrinae, and one in a single Chromis species. Using amino acid comparisons, we found that known spectral tuning sites were altered several times in parallel across the damselfish radiation (through sequence change and duplication followed by sequence change), causing repeated shifts in peak spectral absorbance of around 10 nm. Pomacentrinae damselfishes expressed either one or both copies of SWS1, likely to further finetune UV-signal detection and differentiation. This highly advanced and modified UV vision among damselfishes, in particular the duplication of SWS1 among Pomacentrinae, might be seen as a key evolutionary innovation that facilitated the evolution of the exuberant variety of UV-reflectance traits and the diversification of this coral reef fish lineage.

The Sequential Direct and Indirect Effects of Mountain Uplift, Climatic Niche, and Floral Trait Evolution on Diversification Dynamics in an Andean Plant Clade.

Why and how organismal lineages radiate is commonly studied through either assessing abiotic factors (biogeography, geomorphological processes, and climate) or biotic factors (traits and interactions). Despite increasing awareness that both abiotic and biotic processes may have important joint effects on diversification dynamics, few attempts have been made to quantify the relative importance and timing of these factors, and their potentially interlinked direct and indirect effects, on lineage diversification. We here combine assessments of historical biogeography, geomorphology, climatic niche, vegetative, and floral trait evolution to test whether these factors jointly, or in isolation, explain diversification dynamics of a Neotropical plant clade (Merianieae, Melastomataceae). After estimating ancestral areas and the changes in niche and trait disparity over time, we employ Phylogenetic Path Analyses as a synthesis tool to test eleven hypotheses on the individual direct and indirect effects of these factors on diversification rates. We find strongest support for interlinked effects of colonization of the uplifting Andes during the mid-Miocene and rapid abiotic climatic niche evolution in explaining a burst in diversification rate in Merianieae. Within Andean habitats, later increases in floral disparity allowed for the exploitation of wider pollination niches (i.e., shifts from bee to vertebrate pollinators), but did not affect diversification rates. Our approach of including both vegetative and floral trait evolution, rare in assessments of plant diversification in general, highlights that the evolution of woody habit and larger flowers preceded the colonization of the Andes, but was likely critical in enabling the rapid radiation in montane environments. Overall, and in concert with the idea that ecological opportunity is a key element of evolutionary radiations, our results suggest that a combination of rapid niche evolution and trait shifts was critical for the exploitation of newly available niche space in the Andes in the mid-Miocene. Further, our results emphasize the importance of incorporating both abiotic and biotic factors into the same analytical framework if we aim to quantify the relative and interlinked effects of these processes on diversification.

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk.

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

A preliminary taxonomic revision of the girii clade of South Asian Cnemaspis Strauch, 1887 (Squamata: Gekkonidae) with the description of four new species from southern Maharashtra, India

The girii clade is the most diverse clade of South Asian Cnemaspis in the Northern Western Ghats, including ten valid species. Recent collections of four divergent lineages from southern Maharashtra prompted us to examine the type material of all species in the girii clade, which led to the discovery of numerous errors and discrepancies in the literature. We take this opportunity to provide an updated and expanded morphological dataset including diagnoses and comparisons for all species of the girii clade and to build an almost complete species-level phylogeny including the first ND2 sequences of C. rajgadensis and C. uttaraghati, as well as the four new species described herein (but excluding C. fortis). We describe the four divergent mitochondrial lineages as new species using morphological data, and provide a key to species of the girii clade. Finally, we list some inconsistences in previous publications and data for the girii clade.

Holarctic Lineages Cannot Inform Diversity and Evolution in the Neotropics – the barklice family Psocidae as a case study

Abstract Despite tropical species comprising nearly 60% of Psocidae species, previous studies examining the Psocidae phylogeny have undersampled tropical diversity (<40% species in trees). Here we discuss the systematics of the Psocidae based on the most comprehensive species-level sampling of the Psocidae. We infer the phylogenetic position of 43 previously unsampled Neotropical species in the Psocidae phylogeny. We find that Neotropical psocids are generally not closely related to morphologically similar taxa in the Holarctic region. Consequently, the monophyletic status for the major groups within Psocidae (subfamilies and tribes) is recovered only when Holarctic groups are sampled (7–10 of 11 higher-level groups are monophyletic) but violated when Neotropical species are also sampled in the tree (1 of 11 higher-level groups are monophyletic). Our study pinpoints at the downfalls of simply extending taxonomic knowledge from lineages of a certain area (i.e. Holarctic) to inform diversity and evolution of lineages in other regions (i.e. Neotropics).

A new species of terrestrial toad of the Rhinellafestae group (Anura, Bufonidae) from the highlands of the Central Cordillera of the Andes of Colombia.

The genus Rhinella (Bufonidae) comprises 92 species of Neotropical toads. In Colombia, Rhinella is represented by 22 recognized species, of which nine belong to the Rhinellafestae group. Over the past decade, there has been increasing evidence of cryptic diversity within this group, particularly in the context of Andean forms. Specimens of Rhinella collected in high Andean forests on both slopes of the Central Cordillera in Colombia belong to an undescribed species, Rhinellakumandaysp. nov. Genetic analyses using the mitochondrial 16S rRNA gene indicated that the individuals belong to the festae species group. However, they can be distinguished from other closely related species such as Rhinellaparaguas and Rhinellatenrec by a combination of morphological traits including the presence of tarsal fold, a moderate body size, and substantial genetic divergence in the 16S rRNA gene (> 5%). Through this integrative approach, the specimens from the Central Cordillera of Colombia are considered an evolutionary divergent lineage that is sister to R.paraguas, and described as a new species. Rhinellakumandaysp. nov. is restricted to the Central Cordillera of Colombia inhabiting both slopes in the departments of Caldas and Tolima, in an elevational range between 2420 and 3758 m. With the recognition of this new species, the genus Rhinella now comprises 93 species with 23 of them found in Colombia, and ten species endemic to the country.

Abstract 1761: Comprehensive profiling of the FAT3-associated gene signature in deciphering immunophenotypes of gastrointestinal cancers: Analysis of institutional cohort and TCGA dataset

Abstract Gastrointestinal (GI) cancers are among the most prevalent cancers affecting the US population. Anatomically, GI cancer includes cancers of the organs in the digestive tract, from the esophagus to the rectum. Advances in molecular oncology have started to transform the therapeutic landscape and offer tremendous promise for patients across diverse lineages. However, the equitable benefit of this approach has been hampered due to tumor heterogeneity and the lack of accurate, tumor-agnostic biomarkers with prognosis and predictive utility. Thus, there is a need to identify novel biomarkers that affect gastrointestinal malignancies to improve the management of cancer patients. This study analyzed a single institution's whole genome sequencing dataset to explore the genomic variation in GI cancers. Between 2020 and 2022, 302 AU patients underwent molecular profiling, including whole-exome sequencing (WES) profiling at Caris Life Sciences. WES is an NGS assay that analyzes the DNA sequences of all protein-coding exons in the genome, representing approximately 1-2% of the human genome with coverage of over 22,000 genes. We analyzed the TCGA dataset for comparative analysis, which included cancer data from >1800 GI cancer patients. In addition, survival and network analyses were performed to identify a 6-gene FAT3-related signature for GI cancers. In our cohort, we identified FAT3 (Fat atypical cadherin 3) as the most frequently mutated gene after TP53, APC, and KRAS. FAT3 encodes a cadherin protein involved in cell-cell interactions and adhesion. In the institutional cohort, FAT3 was found to be mutated in 16% of all GI cases. Further analysis of TCGA datasets comprising 1808 GI cancer patients revealed FAT3 mutations in 12% of the cases. To identify additional prognostic biomarkers associated with FAT3, we performed network analysis and identified a 6-gene FAT3-related signature (FAT3, RYK, FAT2, EGFLAM, NTRK3, IGSF9, and HMGA2) that significantly stratified GI cancer patients based on overall survival, progression-free survival, and disease-specific survival. The perturbation profile of the 6-gene signature was associated with 509 patients or 28% of total GI patients. Further, immune deconvolution analysis of stratified patients revealed increased infiltration of immune cells with immunosuppressive phenotypic properties in GI cancer patients with higher expression of the FAT3-related gene signature. In summary, this analysis reveals the distribution of FAT3 mutational profiles in GI patients and identifies a gene signature that stratifies patients based on survival outcomes. This analysis can provide new tools for patient stratification and therapy implementation, leading to better outcomes for cancer patients. Citation Format: Pankaj Kumar Ahluwalia, Tiffanie Leeman, Ashis Mondal, Ashutosh Vashisht, Harmanpreet Singh, Ravindra Kolhe. Comprehensive profiling of the FAT3-associated gene signature in deciphering immunophenotypes of gastrointestinal cancers: Analysis of institutional cohort and TCGA dataset [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1761.

Whole genome sequencing and molecular epidemiology of methicillin-resistant Staphylococcus aureus isolated from patients with bacteraemia in Slovenia.

Data on the molecular epidemiology of methicillin-resistant Staphylococcus aureus isolates from patients with bacteraemia in Slovenia are lacking. The aim of this study was to phenotypically and genotypically investigate 82 MRSA strains isolated from patients with bloodstream infections in central Slovenia between 2019 and 2022. Whole-genome sequencing of selected strains was performed to characterize the strains based on sequence typing, antimicrobial resistance, toxin, and virulence factors genes. Most MRSA carried SCCmec II (63.4%), followed by SCCmec IV (34.1%) and SCCmec V (2.5%). A high proportion of strains belonging to the ST225 lineage (45.1%) was observed, followed by ST97 (18.3%), ST2883 (15.9%), ST22 (9.8%), ST5 (3.7%), and the ST1, ST398 and ST45 lineages (2.4% each). Sixteen different spa types were identified, predominantly ST225-t003 (31.7%), ST97-t359 (15.9%), and ST2883-t4336 (14.6%). None of the strains carried Panton-Valentine leukocidin, exfoliative toxins, or toxic shock toxin. All MRSA strains were susceptible to linezolid, rifampicin, vancomycin, teicoplanin, and trimethoprim-sulfamethoxazole. MRSA strains were resistant to erythromycin, clindamycin, tetracycline and gentamicin, with a frequency of 74.4%, 74.4%, 8.5%, and 1.2%, respectively. This study demonstrates that bacteraemia in central Slovenia is caused by diverse MRSA lineages. Identification of newly emerged lineages should be followed in the future to detect changes in the molecular epidemiology of MRSA in our country.

Abstract 4596: Targeting arginine methylome in 9p21/MTAP-deleted malignant cancers with a next generation PRMT5-specific inhibitor CTS3497

Abstract It remains largely challenging to treat cancer patients with chromosome 9p21-deletion which occurs in approximately 15% of human cancers. Recently, specific targeting of protein arginine methyltransferase 5 (PRMT5), a key member of the type-II PRMT family and a master epigenetic modulator of arginine methylome essential for cancer progression, has emerged as a promising therapy in various hematological and solid tumors with MTAP deletion serving as a precision biomarker. To optimize the potential of targeting PRMT5, we developed CTS3497, a potent, brain-penetrable, orally bioavailable MTA-cooperative PRMT5 inhibitor, using our EpigenPLUSTM platform. The overall excellent drug-like properties of CTS3497 support fast track IND application, and it is anticipated to enter Phase 1 trials in 2024. CTS3497 demonstrated potent PRMT5 inhibitory activity in cellular pharmacodynamics (PD) assays, evidenced by the suppression of symmetric dimethylarginine (SDMA). Additionally, it exhibited profound cell growth inhibition with a low single-digit nM IC50 and 171-fold high selectivity for MTAPnull over isogenic MTAPwt cell lines. Moreover, CTS3497 displayed strong inhibitory effects on cell growth across diverse panels of MTAPnull cancer cell lines and patient-derived xenograft organoids (PDXOs) spanning various lineages. In vivo, oral administration of CTS3497 resulted in significant tumor regression in multiple MTAP-deleted xenograft models, notably MTAPnull brain tumors. The observed anti-tumor activity correlated with PRMT5 modulation, as indicated by reduced SDMA levels. Impressively, sustained tumor growth inhibition led to the complete tumor regression, with no regrowth even after discontinuation of CTS3497 treatment. The brain penetration property of CTS3497 underscores its therapeutic potential for both primary brain tumors, and brain metastases arising from specific cancers like lung cancer. At the molecular level, PRMT5 inhibition by CTS3497 induced substantial and distinctive RNA splicing defects and differentially expressed gene signatures, impacting key cellular activities downstream of arginine and epigenetic reprogramming in MTAPnull cancer cells. Furthermore, a synergistic anti-tumor effect was observed in multiple MTAP-deleted tumors when combining PRMT5 inhibition with CTS2190, an investigational type-I PRMT inhibitor currently in phase I/II clinical trials. In summary, CTS3497, a PRMT5 inhibitor selectively targeting MTAPnull tumor cells, demonstrated robust and durable anti-tumor responses both in vitro and in vivo across diverse lineages of MTAP-deleted tumors. Its superior ADME properties and favorable safety profiles observed in preclinical studies mark a breakthrough in the next generation epigenetic therapy. CTS3497 emerges as a promising therapeutic option for patients with MTAP-deleted malignancies. Citation Format: Yilin Liu, Long Wang, Hui Shi, Yuanyuan Liu, Qiugeng Ouyang, Xingnian Fu, Meng Wang, Jiannan Guo, Yiqin Wang, Youzhen Wang, Guoliang Xu, Yuan Mi, Haiping Wu. Targeting arginine methylome in 9p21/MTAP-deleted malignant cancers with a next generation PRMT5-specific inhibitor CTS3497 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4596.

Higher order polyploids exhibit enhanced desiccation tolerance in the grass Microchloa caffra.

Desiccation tolerance evolved recurrently across diverse plant lineages to enable survival in water limited conditions. Many resurrection plants are polyploid and several groups have hypothesized that polyploidy contributed to the evolution of desiccation tolerance. However, due to the vast phylogenetic distance between resurrection plant lineages, the rarity of desiccation tolerance, and the prevalence of polyploidy in plants, this hypothesis has been difficult to test. Here, we surveyed natural variation in morphological, reproductive, and desiccation tolerance traits across several cytotypes of a single species to test for links between polyploidy and increased resilience. We sampled multiple natural populations of the resurrection grass Microchloa caffra across an environmental gradient ranging from mesic to xeric in South Africa. We describe two distinct ecotypes of M. caffra that occupy different extremes of the environmental gradient and exhibit consistent differences in ploidy, morphological, reproductive, and desiccation tolerance traits in both field and common growth conditions. Interestingly, plants with more polyploid genomes exhibited consistently higher recovery from desiccation, were less reproductive, and larger than plants with smaller genomes and lower ploidy. These data suggest that selective pressures in increasingly xeric sites may play a role in maintaining and increasing desiccation tolerance that are mediated by changes in ploidy.

Contrasting the depths of divergence between gene-tree and coalescent estimates in the North American racers (Colubridae: Coluber constrictor)

Abstract The North American racers (Coluber constrictor) are widely distributed across the Nearctic and numerous studies have demonstrated extensive variation in morphology, ecology, and population genetic structure. Here we take an integrative approach to understand lineage diversification within this taxon by combining genomic sequence capture data, mtDNA sequence data, morphometrics, and ecological niche models. Both the genomic data and mtDNA phylogeographic analyses support five lineages distributed across the range of this species. However, demographic model selection based on these two datasets strongly conflict in both the model of divergence and estimates of timing of lineage divergence. While mtDNA and concatenated genomic data suggest a Miocene origin of these distinct groups, coalescent-based demographic models with the sequence capture data suggest lineage diversification occurred at ~33 kya in allopatry without gene flow. Using linear morphological measurements of head shape we demonstrate that lineages distributed largely east and west of the Mississippi River are distinguishable. Furthermore, ecological niche models demonstrate that lineages distributed in subtropical habitats have environmental niche space that is significantly differentiated from lineages distributed across the continent. Taken together, these results suggest that ecology is an important axis of lineage divergence within this group and that more fine-scale analyses may find even greater differentiation between the populations identified here. This abstract translated to Spanish is avaliable in the Supporting Infromation section (Este resumen traducido al español está disponible en la sección, Supporting Infromation).

RrA, an enzyme from Rhodospirillum rubrum, is a prototype of a new family of short-chain L-asparaginases.

L-Asparaginases (ASNases) catalyze the hydrolysis of L-Asn to L-Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300-400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum, reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N-terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L-Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short-chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild-type RrA, both with and without bound L-Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand-alone subfamily of ASNases that can catalyze the hydrolysis of L-Asn to L-Asp despite the lack of the C-terminal domain that is present in all ASNases described structurally to date.

Combining Molecular, Macroevolutionary, and Macroecological Perspectives on the Generation of Diversity.

Charles Darwin presented a unified process of diversification driven by the gradual accumulation of heritable variation. The growth in DNA databases and the increase in genomic sequencing, combined with advances in molecular phylogenetic analyses, gives us an opportunity to realize Darwin's vision, connecting the generation of variation to the diversification of lineages. The rate of molecular evolution is correlated with the rate of diversification across animals and plants, but the relationship between genome change and speciation is complex: Mutation rates evolve in response to life history and niche; substitution rates are influenced by mutation, selection, and population size; rates of acquisition of reproductive isolation vary between populations; and traits, niches, and distribution can influence diversification rates. The connection between mutation rate and diversification rate is one part of the complex and varied story of speciation, which has theoretical importance for understanding the generation of biodiversity and also practical impacts on the use of DNA to understand the dynamics of speciation over macroevolutionary timescales.

How much should we sequence? An analysis of the Swiss SARS-CoV-2 surveillance effort.

During the SARS-CoV-2 pandemic, many countries directed substantial resources toward genomic surveillance to detect and track viral variants. There is a debate over how much sequencing effort is necessary in national surveillance programs for SARS-CoV-2 and future pandemic threats. We aimed to investigate the effect of reduced sequencing on surveillance outcomes in a large genomic data set from Switzerland, comprising more than 143k sequences. We employed a uniform downsampling strategy using 100 iterations each to investigate the effects of fewer available sequences on the surveillance outcomes: (i) first detection of variants of concern (VOCs), (ii) speed of introduction of VOCs, (iii) diversity of lineages, (iv) first cluster detection of VOCs, (v) density of active clusters, and (vi) geographic spread of clusters. The impact of downsampling on VOC detection is disparate for the three VOC lineages, but many outcomes including introduction and cluster detection could be recapitulated even with only 35% of the original sequencing effort. The effect on the observed speed of introduction and first detection of clusters was more sensitive to reduced sequencing effort for some VOCs, in particular Omicron and Delta, respectively. A genomic surveillance program needs a balance between societal benefits and costs. While the overall national dynamics of the pandemic could be recapitulated by a reduced sequencing effort, the effect is strongly lineage-dependent-something that is unknown at the time of sequencing-and comes at the cost of accuracy, in particular for tracking the emergence of potential VOCs.IMPORTANCESwitzerland had one of the most comprehensive genomic surveillance systems during the COVID-19 pandemic. Such programs need to strike a balance between societal benefits and program costs. Our study aims to answer the question: How would surveillance outcomes have changed had we sequenced less? We find that some outcomes but also certain viral lineages are more affected than others by sequencing less. However, sequencing to around a third of the original effort still captured many important outcomes for the variants of concern such as their first detection but affected more strongly other measures like the detection of first transmission clusters for some lineages. Our work highlights the importance of setting predefined targets for a national genomic surveillance program based on which sequencing effort should be determined. Additionally, the use of a centralized surveillance platform facilitates aggregating data on a national level for rapid public health responses as well as post-analyses.

Turkish isolates of alfalfa mosaic virus belong to a distinct lineage among global population

Complete RNA3 sequences of 15 alfalfa mosaic virus (AMV) isolates sampled from Afyonkarahisar, Ankara, Konya, Tokat, Kayseri, and Nevşehir Provinces of Türkiye were obtained using RT-PCR with two primer pairs designed in this study. Molecular and population analyses were performed on them together with 63 global isolates which full genome sequences available in NCBI GenBank. Phylogenetic analysis based on movement protein (MP) and coat protein (CP) regions showed unrelated isolates clustering between the constructed trees, and the new Turkish isolates were dispersed in different phylogroups. Particularly, three new Turkish potato isolates: B6, K1, and K2 occupied a divergent MP lineage, together with one alfalfa isolate from China: ZM1 NX3. S, η, and k parameters of genetic diversity confirmed that CP were more conserved than MP sequences. The 15 new Turkish isolates shared 95.4–100% nucleotide and 94–100% amino acid, and 92.8–100% nucleotide and 93.3–100% amino acid identities among themselves at MP and CP regions, respectively. All tested populations were found to be under very strict purifying constraints (dN/dS ratios of 0.061–0.219), and had undergone recent expansion, probably due to new mutations in their genome, according to diversity and neutrality tests. Fixation index (FST) values obtained by pairwise comparisons of phylogroups were all above 0.33, suggested that the isolates clustering has been done convincingly.

Geomorphic impacts within Red River Fault and island shifting as witnessed by the phylogeography of the largest water strider

Palaeogeological events and climate oscillations profoundly impact the demographics and distributions of small-range species, increasing the extinction risk. The largest water strider worldwide, Gigantometra gigas (Hemiptera: Gerridae), exhibits restricted distributions in Vietnam and southern China. Herein, we generated three genomic datasets (mitogenomes, 146 nuclear protein-coding genes and single nucleotide polymorphisms) with ecological niche modelling (ENM) to explicitly test whether the present-day distribution of G. gigas actually resulted from geographical and climatic effects. We found that the origin of this largest water strider reached the divergence time of the genus within Gerridae, providing a greater opportunity to explore its response to geographic movements. The right-lateral motion of the Red River Fault facilitated the divergence of two phylogeographic lineages, resulting in the “north–south component” genetic pattern in G. gigas. The Hainan and southeast Vietnam populations of the southern linage were completely separated by the Beibu Gulf but exhibited similar genetic compositions, confirming that Hainan had a continental origin and that Hainan Island joined with the Indo-China Peninsula to promote gene exchange among populations. Additionally, we noticed the low genetic diversity but long demographic history of the northern lineage, which displayed population dynamics opposite to those of other organisms. Integrating the demographic changes and ENM findings revealed that suitable habitat contraction and rapid demographic decline during the Last Glacial Maximum (LGM) triggered the low genetic diversity of the northern lineage. Overall, the demographic history of the largest water strider was mainly shaped by geographical features, and first provided evidence from the phylogeographic perspective of aquatic insects to support the hypothesis of Hainan Island shifting.

Phylogenomics reveals Adeleorina are an ancient and distinct subgroup of Apicomplexa

Apicomplexans are a diverse phylum of unicellular eukaryotes that share obligate relationships with terrestrial and aquatic animal hosts. Many well-studied apicomplexans are responsible for several deadly zoonotic and human diseases, most notably malaria caused by Plasmodium. Interest in the evolutionary origin of apicomplexans has also spurred recent work on other more deeply-branching lineages, especially gregarines and sister groups like squirmids and chrompodellids. But a full picture of apicomplexan evolution is still lacking several lineages, and one major, diverse lineage that is notably absent is the adeleorinids. Adeleorina apicomplexans comprises hundreds of described species that infect invertebrate and vertebrate hosts across the globe. Although historically considered coccidians, phylogenetic trees based on limited data have shown conflicting branch positions for this subgroup, leaving this question unresolved. Phylogenomic trees and large-scale analyses comparing cellular functions and metabolism between major subgroups of apicomplexans have not incorporated Adeleorina because only a handful of molecular markers and a couple organellar genomes are available, ultimately excluding this group from contributing to our understanding of apicomplexan evolution and biology. To address this gap, we have generated complete genomes from mitochondria and plastids, as well as multiple deep-coverage single-cell transcriptomes of nuclear genes from two Adeleorina species, Klossia helicina and Legerella nova, and inferred a 206-protein phylogenomic tree of Apicomplexa. We observed distinct structures reported in species descriptions as remnant host structures surrounding adeleorinid oocysts. Klossia helicina and L. nova branched, as expected, with monoxenous adeleorinids within the Adeleorina and their mitochondrial and plastid genomes exhibited similarity to published organellar adeleorinid genomes. We show with a phylogeneomic tree and subsequent phylogenomic analyses that Adeleorina are not closely related to any of the currently sampled apicomplexan subgroups, and instead fall as a sister to a large clade encompassing Coccidia, Protococcidia, Hematozoa, and Nephromycida, collectively. This resolves Adeleorina as a key independently-branching group, separate from coccidians, on the tree of Apicomplexa, which now has all known major lineages sampled.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.