Structure of the Maturing 90S Pre-ribosome in Association with the RNA Exosome.

Ribosome assembly is catalyzed by numerous trans-acting factors and coupled with irreversible pre-rRNA processing, driving the pathway toward mature ribosomal subunits. One decisive step early in this progression is removal of the 5' external transcribed spacer (5'-ETS), an RNA extension at the 18S rRNA that is integrated into the huge 90S pre-ribosome structure. Upon endo-nucleolytic cleavage at an internal site, A1, the 5'-ETS is separated from the 18S rRNA and degraded. Here we present biochemical and cryo-electron microscopy analyses that depict the RNA exosome, a major 3'-5' exoribonuclease complex, in a super-complex with the 90S pre-ribosome. The exosome is docked to the 90S through its co-factor Mtr4 helicase, a processive RNA duplex-dismantling helicase, which strategically positions the exosome at the base of 5'-ETS helices H9-H9', which are dislodged in our 90S-exosome structures. These findings suggest a direct role of the exosome in structural remodeling of the 90S pre-ribosome to drive eukaryotic ribosome synthesis.

ZAKα Recognizes Stalled Ribosomes through Partially Redundant Sensor Domains.

Impairment of ribosome function activates the MAPKKK ZAK, leading to activation of mitogen-activated protein (MAP) kinases p38 and JNK and inflammatory signaling. The mechanistic basis for activation of this ribotoxic stress response (RSR) remains completely obscure. We show that the long isoform of ZAK (ZAKα) directly associates with ribosomes by inserting its flexible C terminus into the ribosomal intersubunit space. Here, ZAKα binds helix 14 of 18S ribosomal RNA (rRNA). An adjacent domain in ZAKα also probes the ribosome, and together, these sensor domains are critically required for RSR activation after inhibition of both the E-site, the peptidyl transferase center (PTC), and ribotoxin action. Finally, we show that ablation of the RSR response leads to organismal phenotypes and decreased lifespan in the nematode Caenorhabditis elegans (C. elegans). Our findings yield mechanistic insight into how cells detect ribotoxic stress and provide experimental in vivo evidence for its physiological importance.

The kinase Rio1 and a ribosome collision-dependent decay pathway survey the integrity of 18S rRNA cleavage.

The 18S rRNA sequence is highly conserved, particularly at its 3'-end, which is formed by the endonuclease Nob1. How Nob1 identifies its target sequence is not known, and in vitro experiments have shown Nob1 to be error-prone. Moreover, the sequence around the 3'-end is degenerate with similar sites nearby. Here, we used yeast genetics, biochemistry, and next-generation sequencing to investigate a role for the ATPase Rio1 in monitoring the accuracy of the 18S rRNA 3'-end. We demonstrate that Nob1 can miscleave its rRNA substrate and that miscleaved rRNA accumulates upon bypassing the Rio1-mediated quality control (QC) step, but not in healthy cells with intact QC mechanisms. Mechanistically, we show that Rio1 binding to miscleaved rRNA is weaker than its binding to accurately processed 18S rRNA. Accordingly, excess Rio1 results in accumulation of miscleaved rRNA. Ribosomes containing miscleaved rRNA can translate, albeit more slowly, thereby inviting collisions with trailing ribosomes. These collisions result in degradation of the defective ribosomes utilizing parts of the machinery for mRNA QC. Altogether, the data support a model in which Rio1 inspects the 3'-end of the nascent 18S rRNA to prevent miscleaved 18S rRNA-containing ribosomes from erroneously engaging in translation, where they induce ribosome collisions. The data also demonstrate how ribosome collisions purify cells of altered ribosomes with different functionalities, with important implications for the concept of ribosome heterogeneity.

Thermophile 90S Pre-ribosome Structures Reveal the Reverse Order of Co-transcriptional 18S rRNA Subdomain Integration.

Eukaryotic ribosome biogenesis involves RNA folding and processing that depend on assembly factors and small nucleolar RNAs (snoRNAs). The 90S (SSU-processome) is the earliest pre-ribosome structurally analyzed, which was suggested to assemble stepwise along the growing pre-rRNA from 5' > 3', but this directionality may not be accurate. Here, by analyzing the structure of a series of 90S assembly intermediates from Chaetomium thermophilum, we discover a reverse order of 18S rRNA subdomain incorporation. Large parts of the 18S rRNA 3' and central domains assemble first into the 90S before the 5' domain is integrated. This final incorporation depends on a contact between a heterotrimer Enp2-Bfr2-Lcp5 recruited to the flexible 5' domain and Kre33, which reconstitutes the Kre33-Enp-Brf2-Lcp5 module on the compacted 90S. Keeping the 5' domain temporarily segregated from the 90S scaffold could provide extra time to complete the multifaceted 5' domain folding, which depends on a distinct set of snoRNAs and processing factors.

Enhanced ac4C detection in RNA via chemical reduction and cDNA synthesis with modified dNTPs.

The functional analysis of epitranscriptomic modifications in RNA is constrained by a lack of methods that accurately capture their locations and levels. We previously demonstrated that the RNA modification N4-acetylcytidine (ac4C) can be mapped at base resolution through sodium borohydride reduction to tetrahydroacetylcytidine (tetrahydro-ac4C), followed by cDNA synthesis to misincorporate adenosine opposite reduced ac4C sites, culminating in C:T mismatches at acetylated cytidines (RedaC:T). However, this process is relatively inefficient, resulting in <20% C:T mismatches at a fully modified ac4C site in 18S rRNA. Considering that ac4C locations in other substrates including mRNA are unlikely to reach full penetrance, this method is not ideal for comprehensive mapping. Here, we introduce "RetraC:T" (reduction to tetrahydro-ac4C and reverse transcription with amino-dATP to induce C:T mismatches) as a method with enhanced ability to detect ac4C in cellular RNA. In brief, RNA is reduced through NaBH4 or the closely related reagent sodium cyanoborohydride (NaCNBH3) followed by cDNA synthesis in the presence of a modified DNA nucleotide, 2-amino-dATP, that preferentially binds to tetrahydro-ac4C. Incorporation of the modified dNTP substantially improved C:T mismatch rates, reaching stoichiometric detection of ac4C in 18S rRNA. Importantly, 2-amino-dATP did not result in truncated cDNA products nor increase mismatches at other locations. Thus, modified dNTPs are introduced as a new addition to the toolbox for detecting ac4C at base resolution.

DNA-PKcs has KU-dependent function in rRNA processing and haematopoiesis.

The DNA-dependent protein kinase (DNA-PK), which comprises the KU heterodimer and a catalytic subunit (DNA-PKcs), is a classical non-homologous end-joining (cNHEJ) factor1. KU binds to DNA ends, initiates cNHEJ, and recruits and activates DNA-PKcs. KU also binds to RNA, but the relevance of this interaction in mammals is unclear. Here we use mouse models to show that DNA-PK has an unexpected role in the biogenesis of ribosomal RNA (rRNA) and in haematopoiesis. The expression of kinase-dead DNA-PKcs abrogates cNHEJ2. However, most mice that both expressed kinase-dead DNA-PKcs and lacked the tumour suppressor TP53 developed myeloid disease, whereas all other previously characterized mice deficient in both cNHEJ and TP53 expression succumbed to pro-B cell lymphoma3. DNA-PK autophosphorylates DNA-PKcs, which is its best characterized substrate. Blocking the phosphorylation of DNA-PKcs at the T2609 cluster, but not the S2056 cluster, led to KU-dependent defects in 18S rRNA processing, compromised global protein synthesis in haematopoietic cells and caused bone marrow failure in mice. KU drives the assembly of DNA-PKcs on a wide range of cellular RNAs, including the U3 small nucleolar RNA, which is essential for processing of 18S rRNA4. U3 activates purified DNA-PK and triggers phosphorylation of DNA-PKcs at T2609. DNA-PK, but not other cNHEJ factors, resides in nucleoli in an rRNA-dependent manner and is co-purified with the small subunit processome. Together our data show that DNA-PK has RNA-dependent, cNHEJ-independent functions during ribosome biogenesis that require the kinase activity of DNA-PKcs and its phosphorylation at the T2609 cluster.

The Efg1-Bud22 dimer associates with the U14 snoRNP contacting the 5' rRNA domain of an early 90S pre-ribosomal particle.

The DEAD-box helicase Dbp4 plays an essential role during the early assembly of the 40S ribosome, which is only poorly understood to date. By applying the yeast two-hybrid method and biochemical approaches, we discovered that Dbp4 interacts with the Efg1-Bud22 dimer. Both factors associate with early pre-90S particles and smaller complexes, each characterized by a high presence of the U14 snoRNA. A crosslink analysis of Bud22 revealed its contact to the U14 snoRNA and the 5' domain of the nascent 18S rRNA, close to its U14 snoRNA hybridization site. Moreover, depletion of Bud22 or Efg1 specifically affects U14 snoRNA association with pre-ribosomal complexes. Accordingly, we concluded that the role of the Efg1-Bud22 dimer is linked to the U14 snoRNA function on early 90S ribosome intermediates chaperoning the 5' domain of the nascent 18S rRNA. The successful rRNA folding of the 5' domain and the release of Efg1, Bud22, Dpb4, U14 snoRNA and associated snoRNP factors allows the subsequent recruitment of the Kre33-Bfr2-Enp2-Lcp5 module towards the 90S pre-ribosome.

High resolution landscape of ribosomal RNA processing and surveillance.

Ribosomal RNAs are processed in a complex pathway. We profiled rRNA processing intermediates in yeast at single-molecule and single-nucleotide levels with circularization, targeted amplification and deep sequencing (CircTA-seq), gaining significant mechanistic insights into rRNA processing and surveillance. The long form of the 5' end of 5.8S rRNA is converted to the short form and represents an intermediate of a unified processing pathway. The initial 3' end processing of 5.8S rRNA involves trimming by Rex1 and Rex2 and Trf4-mediated polyadenylation. The 3' end of 25S rRNA is formed by sequential digestion by four Rex proteins. Intermediates with an extended A1 site are generated during 5' degradation of aberrant 18S rRNA precursors. We determined precise polyadenylation profiles for pre-rRNAs and show that the degradation efficiency of polyadenylated 20S pre-rRNA critically depends on poly(A) lengths and degradation intermediates released from the exosome are often extensively re-polyadenylated.

Arabidopsis HOT3/eIF5B1 constrains rRNA RNAi by facilitating 18S rRNA maturation.

Ribosome biogenesis is essential for protein synthesis in gene expression. Yeast eIF5B has been shown biochemically to facilitate 18S ribosomal RNA (rRNA) 3' end maturation during late-stage 40S ribosomal subunit assembly and gate the transition from translation initiation to elongation. But the genome-wide effects of eIF5B have not been studied at the single-nucleotide resolution in any organism, and 18S rRNA 3' end maturation is poorly understood in plants. Arabidopsis HOT3/eIF5B1 was found to promote development and heat stress acclimation by translational regulation, but its molecular function remained unknown. Here, we show that HOT3 is a late-stage ribosome biogenesis factor that facilitates 18S rRNA 3' end processing and is a translation initiation factor that globally impacts the transition from initiation to elongation. By developing and implementing 18S-ENDseq, we revealed previously unknown events in 18S rRNA 3' end maturation or metabolism. We quantitatively defined processing hotspots and identified adenylation as the prevalent nontemplated RNA addition at the 3' ends of pre-18S rRNAs. Aberrant 18S rRNA maturation in hot3 further activated RNA interference to generate RDR1- and DCL2/4-dependent risiRNAs mainly from a 3' portion of 18S rRNA. We further showed that risiRNAs in hot3 were predominantly localized in ribosome-free fractions and were not responsible for the 18S rRNA maturation or translation initiation defects in hot3. Our study uncovered the molecular function of HOT3/eIF5B1 in 18S rRNA maturation at the late 40S assembly stage and revealed the regulatory crosstalk among ribosome biogenesis, messenger RNA (mRNA) translation initiation, and siRNA biogenesis in plants.

Maize RNA 3'-terminal phosphate cyclase-like protein promotes 18S pre-rRNA cleavage and is important for kernel development.

Plant ribosomes contain four specialized ribonucleic acids, the 5S, 5.8S, 18S, and 25S ribosomal RNAs (rRNAs). Maturation of the latter three rRNAs requires cooperative processing of a single transcript by several endonucleases and exonucleases at specific sites. In maize (Zea mays), the exact nucleases and components required for rRNA processing remain poorly understood. Here, we characterized a conserved RNA 3'-terminal phosphate cyclase (RCL)-like protein, RCL1, that functions in 18S rRNA maturation. RCL1 is highly expressed in the embryo and endosperm during early seed development. Loss of RCL1 function resulted in lethality due to aborted embryo cell differentiation. We also observed pleiotropic defects in the rcl1 endosperm, including abnormal basal transfer cell layer growth and aleurone cell identity, and reduced storage reserve accumulation. The rcl1 seeds had lower levels of mature 18S rRNA and the related precursors were altered in abundance compared with wild type. Analysis of transcript levels and protein accumulation in rcl1 revealed that the observed lower levels of zein and starch synthesis enzymes mainly resulted from effects at the transcriptional and translational levels, respectively. These results demonstrate that RCL1-mediated 18S pre-rRNA processing is essential for ribosome function and messenger RNA translation during maize seed development.

Dissecting Trypanosoma brucei RRP44 function in the maturation of segmented ribosomal RNA using a regulated genetic complementation system.

Trypanosoma brucei belongs to a group of protozoans presenting fragmented large subunit rRNA. Its LSU rRNA equivalent to the 25S/28S rRNA of other eukaryotes is split into six fragments, requiring additional processing for removal of the extra spacer sequences. We have used a genetic complementation strategy to further investigate the T. brucei RRP44 nuclease in pre-rRNA maturation. TbRRP44 contains both a PIN and a RNB domain whose homologues are found in association with the exosome complex. We found that the exonucleolytic activity of the RNB domain as well as the physical presence of the PIN domain are essential for TbRRP44 function, while a catalytic site mutation in the PIN domain has no detectable effect on cell growth. A new endonucleolytic cleavage site in ITS1 was identified. In addition to the 5.8S rRNA 3'-end maturation, TbRRP44 is required for degradation of the excised 5'-ETS and for removal of part of ITS1 during maturation of the 18S rRNA 3'-end. TbRRP44 deficiency leads to accumulation of many LSU intermediate precursors, most of them not detected in control cells. TbRRP44 is also required for U3 snoRNA and spliced leader processing, indicating that TbRRP44 may have a wide role in RNA processing in T. brucei.

Tafenoquine for Relapsing Babesiosis: A Case Series.

BACKGROUND: Relapsing babesiosis often occurs in highly immunocompromised patients and has been attributed to the acquisition of resistance against drugs commonly used for treatment such as atovaquone, azithromycin, and clindamycin. Tafenoquine, which is approved for malaria prophylaxis and presumptive antirelapse treatment of Plasmodium vivax malaria, has shown activity against Babesia microti in several animal models of acute infection and in a single human case of relapsing babesiosis. Here, we report 5 cases of relapsing babesiosis treated with tafenoquine, including the previous case, and begin to define the conditions for optimal use of tafenoquine in relapsing babesiosis. METHODS: A definitive diagnosis of babesiosis was made by microscopic examination of Giemsa-stained thin blood smears or a real-time polymerase chain reaction (PCR) that targets the parasite 18S rRNA gene. Clearance of B. microti infection was ascertained by use of blood smear and real-time PCR. RESULTS: Tafenoquine was initiated with a loading dose of 600 mg. A weekly maintenance dose consisted of 200 mg or 300 mg; the lower dose was associated with a delayed clearance of B. microti. In 2 cases, all antimicrobial agents but tafenoquine were discontinued prior to clearance of infection. In 2 other cases, clearance was achieved while tafenoquine was administered along with other antimicrobial agents. In 3 of these 4 cases, tafenoquine was used in combination with atovaquone-proguanil. Other agents included atovaquone, azithromycin, and/or clindamycin. In 1 case, tafenoquine was administered alone and failed to prevent relapse. CONCLUSIONS: Tafenoquine can be a useful adjunct for the treatment of highly immunocompromised patients experiencing relapsing babesiosis caused by B. microti.

Autochthonous Human Babesia divergens Infection, England.

We describe a case of autochthonous human Babesia divergens infection in an immunocompetent woman in England. The patient had fever, hemolysis, multiorgan failure, and 18% parasitemia. We confirmed B. divergens by 18S rDNA PCR and sequencing. Clinicians should consider babesiosis as a differential diagnosis in patients with unexplained hemolysis.

Confirmed Case of Autochthonous Human Babesiosis, Hungary.

We report a case of autochthonous human babesiosis in Hungary, confirmed by PCR and partial sequencing of the Babesia spp. 18S rRNA gene. Babesiosis should be considered during the differential diagnosis of febrile illnesses, and peripheral blood smears to detect Babesia spp. should be part of the routine clinical workup.

Macroalgal eukaryotic microbiome composition indicates novel phylogenetic diversity and broad host spectrum of oomycete pathogens.

Seaweeds are important components of marine ecosystems with emerging potential in aquaculture and as sources of biofuel, food products and pharmacological compounds. However, an increasingly recognised threat to natural and industrial seaweed populations is infection with parasitic single-celled eukaryotes from the relatively understudied oomycete lineage. Here we examine the eukaryomes of diverse brown, red and green marine macroalgae collected from polar (Baffin Island), cold-temperate (Falkland Islands) and tropical (Ascension Island) locations, with a focus on oomycete and closely related diatom taxa. Using 18S rRNA gene amplicon sequencing, we show unexpected genetic and taxonomic diversity of the eukaryomes, a strong broad-brush association between eukaryome composition and geographic location, and some evidence of association between eukaryome structure and macroalgal phylogenetic relationships (phylosymbiosis). However, the oomycete fraction of the eukaryome showed disparate patterns of diversity and structure, highlighting much weaker association with geography and no evidence of phylosymbiosis. We present several novel haplotypes of the most common oomycete Eurychasma dicksonii and report for the first time a cosmopolitan distribution and absence of host specificity of this important pathogen. This indicates rich diversity in macroalgal oomycete pathogens and highlights that these pathogens may be generalist and highly adaptable to diverse environmental conditions.

Dnttip2 is essential for 18S rRNA processing and digestive organ development in zebrafish.

Dnttip2 is one of the components of the small subunit (SSU) processome. In yeast, depletion of dnttip2 leads to an inefficient processing of pre-rRNA and a decrease in synthesis of the mature 18S rRNA. However, the biological roles of Dnttip2 in higher organisms are poorly defined. In this study, we demonstrate that dnttip2 is a maternal gene in zebrafish. Depletion of Dnttip2 leads to embryonic lethal with severe digestive organs hypoplasia. The loss of function of Dnttip2 also leads to partial defects in cleavage at the A0-site and E-site during 18S rRNA processing. In conclusion, Dnttip2 is essential for 18S rRNA processing and digestive organ development in zebrafish.

First Karyotypic Insights into Potamotrygon schroederi Fernández-Yépez, 1958: Association of Different Classes of Repetitive DNA.

INTRODUCTION: Currently, there are 38 valid species of freshwater stingrays, and these belong to the subfamily Potamotrygoninae. However, cytogenetic information about this group is limited, with studies mainly using classical techniques, Giemsa, and C-banding. METHODS: In this study, we used classical and molecular cytogenetic techniques - mapping of 18S and 5S rDNA and simple sequence repeats (SSRs) - in order to investigate the karyotypic composition of Potamotrygon schroederi and reveal the karyoevolutionary trends of this group. RESULTS: The species presented 2n = 66 chromosomes with 18m + 12sm + 16st + 20a, heterochromatic blocks distributed in the centromeric regions of all the chromosomes, and terminal blocks in the q arm of pairs 2 and 3. Mapping of 18S rDNA regions revealed multiple clusters on pairs 2 and 7 and a homolog of pair 24. The 5S rDNA region was found in the pericentromeric portion of the subtelocentric pair 16. Furthermore, dinucleotide SSRs sequences were found in the centromeric and terminal regions of different chromosomal pairs, with preferential accumulation in pair 17. In addition, we identified conspicuous blocks of (GATA)n and (GACA)n sequences colocalized with the 5S rDNA (pair 16). CONCLUSION: In general, this study corroborates the general trend of a reduction in 2n in the species of Potamotrygoninae subfamily. Moreover, we found that the location of rDNA regions is very similar among Potamotrygon species, and the SSRs accumulation in the second subtelocentric pair (17) seems to be a common trait in this genus.

Comparative Cytogenetics in Tyrannidae (Aves, Passeriformes): High Genetic Diversity despite Conserved Karyotype Organization.

INTRODUCTION: Passeriformes has the greatest species diversity among Neoaves, and the Tyrannidae is the richest in this order with about 600 valid species. The diploid number of this family remains constant, ranging from 2n = 76 to 84, but the chromosomal morphology varies, indicating the occurrence of different chromosomal rearrangements. Cytogenetic studies of the Tyrannidae remain limited, with approximately 20 species having been karyotyped thus far. This study aimed to describe the karyotypes of two species from this family, Myiopagis viridicata and Sirystes sibilator. METHODS: Skin biopsies were taken from each individual to establish fibroblast cell cultures and to obtain chromosomal preparations using the standard methodology. The chromosomal distribution of constitutive heterochromatin was investigated by C-banding, while the location of simple repetitive sequences (SSRs), 18S rDNA, and telomeric sequences was found through fluorescence in situ hybridization. RESULTS: The karyotypes of both species are composed of 2n = 80. The 18S rDNA probes hybridized into two pairs of microchromosomes in M. viridicata, but only a single pair in S. sibilator. Only the telomeric portions of each chromosome in both species were hybridized by the telomere sequence probes. Most of the SSRs were found accumulated in the centromeric and telomeric regions of several macro- and microchromosomes in both species, which likely correspond to the heterochromatin-rich regions. CONCLUSION: Although both species analyzed showed a conserved karyotype organization (2n = 80), our study revealed significant differences in their chromosomal architecture, rDNA distribution, and SSR accumulation. These findings were discussed in the context of the evolution of Tyrannidae karyotypes.

Antisense pairing and SNORD13 structure guide RNA cytidine acetylation.

N4-acetylcytidine (ac4C) is an RNA nucleobase found in all domains of life. The establishment of ac4C in helix 45 (h45) of human 18S ribosomal RNA (rRNA) requires the combined activity of the acetyltransferase NAT10 and the box C/D snoRNA SNORD13. However, the molecular mechanisms governing RNA-guided nucleobase acetylation in humans remain unexplored. After applying comparative sequence analysis and site-directed mutagenesis to provide evidence that SNORD13 folds into three main RNA helices, we report two assays that enable the study of SNORD13-dependent RNA acetylation in human cells. First, we demonstrate that ectopic expression of SNORD13 rescues h45 in a SNORD13 knockout cell line. Next, we show that mutant snoRNAs can be used in combination with nucleotide resolution ac4C sequencing to define structure and sequence elements critical for SNORD13 function. Finally, we develop a second method that reports on the substrate specificity of endogenous NAT10-SNORD13 via mutational analysis of an ectopically expressed pre-rRNA substrate. By combining mutational analysis of these reconstituted systems with nucleotide resolution ac4C sequencing, our studies reveal plasticity in the molecular determinants underlying RNA-guided cytidine acetylation that is distinct from deposition of other well-studied rRNA modifications (e.g., pseudouridine). Overall, our studies provide a new approach to reconstitute RNA-guided cytidine acetylation in human cells as well as nucleotide resolution insights into the mechanisms governing this process.

The enhanced neutral process with decreasing cell size: a study on phytoplankton metacommunities from the glacier-fed river of Qinghai-Xizang Plateau.

The cell size of phytoplankton is an important defining functional trait that can serve as a driver and sentinel of phytoplankton community structure and function. However, the study of the assembly patterns and drivers of phytoplankton metacommunities with different cell sizes has not been widely carried out. In this study, we systematically investigated the biodiversity patterns, drivers, and assembly processes of the three phytoplankton cell sizes (micro: 20-200 µm; nano: 2-20 µm; pico: 0.2-2 µm) in the Za'gya Zangbo River from the source to the estuary using 18S rDNA amplicon sequencing. The results demonstrated that the alpha diversity and co-occurrence network complexity for all three sizes of phytoplankton increased to a peak downstream of the glacier sources and then decreased to the estuary. The nanophytoplankton subcommunity consistently had the highest alpha diversity and co-occurrence network complexity. On the other hand, total beta diversity followed a unimodal trend of decreasing and then increasing from source to estuary, and was dominated by species replacement components. In addition, deterministic processes driven mainly by physiochemical indices (PCIs) and biogenic elements (BGEs) dominated the assembly of micro- and nanophytoplankton subcommunities, whereas stochastic processes driven by geographical factors (GGFs) dominated the assembly of picophytoplankton subcommunities. The results explained the contradictions in previous studies of phytoplankton community assembly processes in highland aquatic ecosystems, elucidating the different contributions of deterministic and stochastic processes, and the complexity of compositional mechanisms in shaping the assembly of micro-, nano-, and picophytoplankton in this highland glacial river. IMPORTANCE: The cell size of phytoplankton is a key life-history trait and key determinant, and phytoplankton of different cell sizes are differentially affected by ecological processes. However, the study of the assembly patterns and drivers of phytoplankton metacommunities with different cell sizes has not been widely carried out. We provide an in-depth analysis of phytoplankton community diversity across three cell sizes in the glacier-fed river, describing how the pattern of phytoplankton communities differs across cell sizes in response to geochemical gradients. The results show that the smaller phytoplankton (picophytoplankton) are relatively more influenced by dispersal-based stochastic processes, whereas larger ones (microphytoplankton and nanophytoplankton) are more structured by selection-based deterministic processes.