Branching Lineages Research Articles

Exploring the contribution of Zfp521/ZNF521 on primary hematopoietic stem/progenitor cells and leukemia progression.

Hematopoietic stem cells (HSCs) drive cellular turnover in the hematopoietic system by balancing self-renewal and differentiation. In the adult bone marrow (BM), these cells are regulated by a complex cellular microenvironment known as "niche," which involves dynamic interactions between diverse cellular and non-cellular elements. During blood cell maturation, lineage branching is guided by clusters of genes that interact or counteract each other, forming complex networks of lineage-specific transcription factors. Disruptions in these networks can lead to obstacles in differentiation, lineage reprogramming, and ultimately malignant transformation, including acute myeloid leukemia (AML). Zinc Finger Protein 521 (Znf521/Zfp521), a conserved transcription factor enriched in HSCs in both human and murine hematopoiesis, plays a pivotal role in regulating HSC self-renewal and differentiation. Its enforced expression preserves progenitor cell activity, while inhibition promotes differentiation toward the lymphoid and myeloid lineages. Transcriptomic analysis of human AML patient samples has revealed upregulation of ZNF521 in AMLs with the t(9;11) fusion gene MLL-AF9. In vitro studies have shown that ZNF521 collaborates with MLL-AF9 to enhance the growth of transformed leukemic cells, increase colony formation, and activate MLL target genes. Conversely, inhibition of ZNF521 using short-hairpin RNA (shRNA) results in decreased leukemia proliferation, reduced colony formation, and induction of cell cycle arrest in MLL-rearranged AML cell lines. In vivo experiments have demonstrated that mZFP521-deficient mice transduced with MLL-AF9 experience a delay in leukemia development. This review provides an overview of the regulatory network involving ZNF521, which plays a crucial role in controlling both HSC self-renewal and differentiation pathways. Furthermore, we examine the impact of ZNF521 on the leukemic phenotype and consider it a potential marker for MLL-AF9+ AML.

Leaf functional trait evolution and its putative climatic drivers in African Coffea species.

Leaf traits are known to be strong predictors of plant performance and can be expected to (co)vary along environmental gradients. We investigated the variation, integration, environmental relationships, and evolutionary history of leaf functional traits in the genus Coffea L., typically a rainforest understory shrub, across Africa. A better understanding of the adaptive processes involved in leaf trait evolution can inform the use and conservation of coffee genetic resources in a changing climate. We used phylogenetic comparative methods to investigate the evolution of six leaf traits measured from herbarium specimens of 58 African Coffea species. We added environmental data and data on maximum plant height for each species to test trait-environment correlations in various (sub)clades, and we compared continuous trait evolution models to identify variables driving trait diversification. A substantial leaf trait variation was detected across the genus Coffea in Africa, which was mostly interspecific. Of these traits, stomatal size and stomatal density exhibited a clear trade-off. We observed low densities of large stomata in early branching lineages and higher densities of smaller stomata in more recent taxa, which we hypothesise to be related to declining CO2 levels since the mid-Miocene. Brownian Motion evolution was rejected in favour of White Noise or Ornstein-Uhlenbeck models for all traits, implying these traits are adaptively significant rather than driven by pure drift. The evolution of leaf area was likely driven by precipitation, with smaller leaves in dryer climates across the genus. Generally, Coffea leaf traits appear to be evolutionarily labile and governed by stabilising selection, though evolutionary patterns and correlations differ depending on the traits and clades considered. Our study highlights the importance of a phylogenetic perspective when studying trait relationships across related taxa, as well as the consideration of various taxonomic ranges.

Phylogeny and evolution of streptophyte algae.

The Streptophyta emerged about a billion years ago. Nowadays, this branch of the green lineage is most famous for one of its clades, the land plants (Embryophyta). Although Embryophyta make up the major share of species numbers in Streptophyta, there is a diversity of probably >5000 species of streptophyte algae that form a paraphyletic grade next to land plants. Here, we focus on the deep divergences that gave rise to the diversity of streptophytes, hence particularly on the streptophyte algae. Phylogenomic efforts have not only clarified the position of streptophyte algae relative to land plants, but recent efforts have also begun to unravel the relationships and major radiations within streptophyte algal diversity. We illustrate how new phylogenomic perspectives have changed our view on the evolutionary emergence of key traits, such as intricate signalling networks that are intertwined with multicellular growth and the chemodiverse hotbed from which they emerged. These traits are key for the biology of land plants but were bequeathed from their algal progenitors.

Floral morphology and development of Pteridophyllum racemosum Siebold & Zucc. (Papaveraceae)

ABSTRACT Papaveraceae are known for their often showy flowers and diverse morphologies. Pteridophyllum racemosum (Siebold & Zucc.) is the only member in the Pteridophylloideae within the Papaveraceae, and its phylogenetic position has long been controversial, and a comprehensive analysis of its floral morphology was lacking. Our study focuses on the floral morphogenesis of P. racemosum. Histological sections complemented with scanning electron microscopy allowed a detailed characterization of P. racemosum floral development including description of the vegetative morphology, showing (1) the fine structure of the unusual pinnate leaves; (2) the well-defined floral structure and disymmetry; (3) a series of landmarks in floral development stages, with emphasis on the development of the gynoecium. Interestingly, the P. racemosum floral architecture is similar to that of the proposed ancestral Papaveraceae flower supporting the phylogenetic placement of P. racemosum as an early branching lineage within the Papaveraceae. Our comprehensive overview adds to the colorful library of floral diversity in Papaveraceae, providing a solid base for comparative analyses.

Mitochondrial genomes provide insights into the Euholognatha (Insecta: Plecoptera)

BackgroundEuholognatha is a monophyletic group within stoneflies comprised by a superfamily Nemouroidea and a family Scopuridae. Based on morphological data, the family-level phylogenetic relationships within Euholognatha are widely accepted, but there is still controversy among different molecular studies. To better understand the phylogeny of all six extant euholognathan families, we sequenced and analyzed seven euholognathan mitogenomes.ResultsThe sequence heterogeneity analysis observed a low degree of compositional heterogeneity in euholognathan mitogenomes. Meanwhile, leuctrid mitogenomes were more heterogeneous than other euholognathan families, which may affect the phylogenetic reconstruction. Phylogenetic analyses with various datasets generated three topologies. The Leuctridae was recovered as the earliest branching lineage, and the sister relationship of Capniidae and Taeniopterygidae was supported by most tree topologies and FcLM analyses. When separately excluding sparsely sampled Scopuridae or high heterogeneity leuctrid taxa, phylogenetic analyses under the same methods generated more stable and consistent tree topologies. Finally, based on the results of this study, we reconstructed the relationships within Euholognatha as: Leuctridae + (Scopuridae + ((Taeniopterygidae + Capniidae) + (Nemouridae + Notonemouridae))).ConclusionOur research shows the potential of data optimizing strategies in reconstructing phylogeny within Euholognatha and provides new insight into the phylogeny of this group.

Multiple Cell Lineages Give Rise to a Cell Type.

It has long been assumed that a specific cell type arises following stepwise specification of cells corresponding to the branching of cell lineages. However, accumulating evidence indicates that multiple and even remote cell lineages can lead to the development of the same cells. Four examples giving different yet new insights will be discussed: skeletal muscle development from precursors with distinct initial histories of transcriptional regulation, lens cell development from remote lineages yet sharing basic transcription factors, blood cell development under intersectional pathways, and neural tissue development from cardiac precursors through the manipulation of just one component of epigenetic regulation. These examples provide flexible and nondogmatic perspectives on developmental cell regulation, fundamentally revising the old model relying on cell lineages.

Clinical, epidemiological, and virological features of Mpox in Beijing, China – May 31–June 21, 2023

ABSTRACT The first locally acquired case in the Chinese mainland was reported on May 31, 2023, lagging behind other countries. In this study, we aimed to examine the early clinical and epidemiological characteristics of the earliest cases of Mpox in Beijing, China. Additionally, we investigated the sequence and transmission patterns of the Mpox virus (MPXV). We analyzed 37 reported cases of Mpox in Beijing from May 31, 2023 to June 21, 2023. The age range of the subjects was 24–51 years. Thirty-six cases (97.3%) were identified in men who have sex with men (MSM), and 19 cases (51.4%) tested positive for the human immunodeficiency virus. Thirty-three cases were symptomatic, while four were asymptomatic. Skin lesions were observed in 32 cases (97.0%), fever in 26 (78.8%), and swollen lymph nodes in 17 (51.5%). Rash typically appeared in the genital or perianal area 1–3 days before fever onset, with a minimum incubation period of 2 days. For individuals with skin rashes, the skin lesion samples showed 100% positivity and low Ct values. There were high oropharyngeal swab (75.8%) and blood (84.6%) positivity rates. All MPXV strains belonged to the B.1.3 branch of the West African lineage. These strains carried 76–86 nucleotide substitutions compared with the reference human MPXV genome, and genetic diversity was observed. Our findings provide the first insights into the landscape of early transmission of Mpox in Beijing and help inform policy formulation in the Chinese mainland.

On Intrinsic Mechanism and Practice Approach for Integrating Red Culture of Ta-pieh Mountains into College Students’ Ideological and Political Education

The red culture of Ta-pieh Mountains represents a major branch of the spiritual lineage of the Chinese Communists, an important support for the cultural confidence in the new era, a valuable spiritual legacy as well as a visual textbook for ideological and political education. Integrating the red culture of Ta-pieh mountains into college students’ ideological and political education can not only provide resources for patriotic education of college students, enhance their political literacy, and ameliorate their excellent ideological quality, but also boost college students' scientific understanding, reinforce the college students’ understanding of the true significance of upholding and strengthening the Party's leadership over higher education in that it is also an important means to heighten the ideology and value orientation of college students. The red culture of Ta-pieh Mountains will be deeply integrated into college students' ideological and political education by implementing internal mechanisms, say, operational mechanism, incentive mechanism, and teacher development mechanism, etc., in combination with the specific practices, say, on-campus ideological and political teaching, off-campus practical education and online learning, etc.

A complete species-level phylogeny of the Erythrura parrotfinches (Aves: Estrildidae)

Identifying species boundaries and phylogenetic relationships among groups of closely related species provides a necessary framework for understanding how biodiversity evolves in natural systems. Here we present a complete phylogeny of the avian genus Erythrura (family Estrildidae) commonly known as parrotfinches, which includes species threatened by habitat loss and the pet trade. Using both mitogenome and reduced-representation genome-wide nuclear DNA sequence data, we reconstructed the evolutionary history of the group by sampling all 12 recognized species, four of which had not previously been studied in a phylogenetic context. We included intra-species geographic sampling that allowed us to comment on species limits in some taxa. We recovered the Gouldian Finch (Chloebia gouldiae) of Australia which has often been placed in the monotypic genus Chloebia, as being sister to a clade comprising all Erythrura species. In addition, we recovered a well-supported clade comprising eight species distributed throughout the Pacific Island eco-region, whereas those species occurring in continental southeast Asia, the Greater Sundas, and the Philippines, were recovered as earlier branching lineages. Of note was the early branching of the Fiji-endemic E. kleinschmidti which corroborates its unique phenotype. We also found a deep phylogenetic split (8.59% corrected, 7.89% uncorrected divergence in the mitochondrial gene ND2) between the Java and Philippine populations of E. hyperythra, indicating unrecognized species-level diversity within this taxon. In contrast, genome-wide nuclear data suggested that the New Guinea endemic species E. papuana is embedded within the widespread species E. trichroa in all phylogenetic reconstructions, corroborating previously published mitochondrial data that suggested a similar pattern. By generating a phylogenetic hypothesis for the relationships among all species of Erythrura parrotfinches, we provide a framework for better understanding the extant diversity and evolutionary history of this group.

The first local case of mpox caused by an imported case in the Chinese mainland

Monkeypox (mpox) is a zoonotic disease caused by the mpox virus (MPXV) that has been primarily limited to Central and West African nations since its discovery. The recent spread of the West African lineage of MPXV in historically unaffected countries has raised concerns for global public health. Despite a significant decrease in global mpox cases, there is still a risk of a global resurgence. This study reports the first local case of mpox caused by an imported case in the Chinese mainland. Polymerase chain reaction (PCR) diagnosed the two cases, and the viral genomes were obtained by next-generation sequencing. Genomic analysis revealed that the two strains shared an identical genome sequence and belonged to the B.1.3 branch of the West African lineage, which is the first local case of mpox caused by an imported case in the Chinese mainland, highlighting the potential threat of mpox in China and the immediate need for adequate surveillance measures.

Systematic reappraisal of the diminutive Amphora thumensis and the transfer of Halamphora parathumensis to the genus Amphora

Since its description the amphoroid diatom Amphora thumensis has confounded simple taxonomic classification having been assigned to Amphora, Cymbella, and now Halamphora. Taxonomic ambiguity around this taxon stems from its small size and reduced, often contradictory, morphological features. In this investigation, we confirm the taxonomic position of A. thumensis by presenting a four-marker molecular phylogeny for a population collected from a small spring-fed lake in the Rocky Mountains of Colorado, USA. Our phylogenetic hypothesis places A. thumensis not within the genus Halamphora, but within the deep branching lineages of Amphora. We discuss this phylogenetic hypothesis in the context of the morphological and ecological diversity of the genus and evaluate the utility of putative diagnostic characters within Amphora. In light of this new evidence, the valid name Amphora thumensis (Mayer) Krieger is identified as the preferred designation for this taxon. Consequently, we transfer the morphologically similar Halamphora parathumensis to Amphora.

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets

Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

Whole genome sequencing of Ethiopian Brucella abortus isolates expands the known diversity of an early branching sub-Saharan African lineage.

Brucellosis remains one of the most significant zoonotic diseases globally, responsible for both considerable human morbidity and economic losses due to its impacts on livestock productivity. Despite this, there remain significant evidence gaps in many low- and middle-income countries, including those of sub-Saharan Africa. Here we report the first molecular characterisation of Brucella sp. from Ethiopia. Fifteen Brucella sp. isolates from an outbreak in cattle from a herd in central Ethiopia were identified as Brucella abortus, using bacterial culture and molecular methods. Sequencing of the Ethiopian B. abortus isolates allowed their phylogenetic comparison with 411 B. abortus strains of diverse geographical origins, using whole genome single nucleotide polymorphisms (wgSNP). The Ethiopian isolates belonged to an early-branching lineage (Lineage A) previously only represented by data from two strains, both of sub-Saharan African origin (Kenya and Mozambique). A second B. abortus lineage (Lineage B), also comprised solely of strains originating from sub-Saharan Africa, was identified. The majority of strains belonged to one of two lineages of strains originating from a much broader geographical range. Further analyses based on multi-locus sequence typing (MLST) and multi-locus variable-number tandem repeat analysis (MLVA) expanded the number of B. abortus strains available for comparison with the Ethiopian isolates and were consistent with the findings from wgSNP analysis. MLST profiles of the Ethiopian isolates expanded the sequence type (ST) diversity of the early branching lineage of B. abortus, equivalent to wgSNP Lineage A. A more diverse cluster of STs, equivalent to wgSNP Lineage B, was comprised solely of strains originating from sub-Saharan Africa. Similarly, analysis of B. abortus MLVA profiles (n = 1891) confirmed that the Ethiopian isolates formed a unique cluster, similar to only two existing strains, and distinct from the majority of other strains of sub-Saharan African origin. These findings expand the known diversity of an under-represented lineage of B. abortus and suggest a potential evolutionary origin for the species in East Africa. In addition to providing information concerning Brucella species extant within Ethiopia this work serves as the basis for further studies on the global population structure and evolutionary history of a major zoonotic pathogen.

Predictive nonlinear modeling of malignant myelopoiesis and tyrosine kinase inhibitor therapy.

Chronic myeloid leukemia (CML) is a blood cancer characterized by dysregulated production of maturing myeloid cells driven by the product of the Philadelphia chromosome, the BCR-ABL1 tyrosine kinase. Tyrosine kinase inhibitors (TKIs) have proved effective in treating CML, but there is still a cohort of patients who do not respond to TKI therapy even in the absence of mutations in the BCR-ABL1 kinase domain that mediate drug resistance. To discover novel strategies to improve TKI therapy in CML, we developed a nonlinear mathematical model of CML hematopoiesis that incorporates feedback control and lineage branching. Cell-cell interactions were constrained using an automated model selection method together with previous observations and new in vivo data from a chimeric BCR-ABL1 transgenic mouse model of CML. The resulting quantitative model captures the dynamics of normal and CML cells at various stages of the disease and exhibits variable responses to TKI treatment, consistent with those of CML patients. The model predicts that an increase in the proportion of CML stem cells in the bone marrow would decrease the tendency of the disease to respond to TKI therapy, in concordance with clinical data and confirmed experimentally in mice. The model further suggests that, under our assumed similarities between normal and leukemic cells, a key predictor of refractory response to TKI treatment is an increased maximum probability of self-renewal of normal hematopoietic stem cells. We use these insights to develop a clinical prognostic criterion to predict the efficacy of TKI treatment and design strategies to improve treatment response. The model predicts that stimulating the differentiation of leukemic stem cells while applying TKI therapy can significantly improve treatment outcomes.

Splitting one species into 22: an unusual tripling of molecular, morphological, and geographical differentiation in the fern family Didymochlaenaceae (Polypodiales).

The pantropical fern genus Didymochlaena (Didymochlaenaceae) has long been considered to contain one species only. Recent studies have resolved this genus/family as either sister to the rest of eupolypods I or as the second branching lineage of eupolypods I, and have shown that this genus is not monospecific, but the exact species diversity is unknown. In this study, a new phylogeny is reconstructed based on an expanded taxon sampling and six molecular markers. Our major results include: (i) Didymochlaena is moderately or weakly supported as sister to the rest of eupolypods I, highlighting the difficulty in resolving the relationships of this important fern lineage in the polypods; (ii) species in Didymochlaena are resolved into a New World clade and an Old World clade, and the latter further into an African clade and an Asian-Pacific clade; (iii) an unusual tripling of molecular, morphological and geographical differentiation in Didymochlaena is detected, suggesting single vicariance or dispersal events in individual regions and no evidence for reversals at all, followed by allopatric speciation at more or less homogeneous rates; (iv) evolution of 18 morphological characters is inferred and two morphological synapomorphies defining the family are recognized-the elliptical sori and fewer than 10 sori per pinnule, the latter never having been suggested before; (v) based on morphological and molecular variation, 22 species in the genus are recognized contrasting with earlier estimates of between one and a few; and (vi) our biogeographical analysis suggests an origin for Didymochlaena in the latest Jurassic-earliest Cretaceous and the initial diversification of the extant lineages in the Miocene-all but one species diverged from their sisters within the last 27 Myr, in most cases associated with allopatric speciation owing to geologic and climatic events, or dispersal.

Oct1 cooperates with the Smad family of transcription factors to promote mesodermal lineage specification.

The transition between pluripotent and tissue-specific states is a key aspect of development. Understanding the pathways driving these transitions will facilitate the engineering of properly differentiated cells for experimental and therapeutic uses. Here, we showed that during mesoderm differentiation, the transcription factor Oct1 activated developmental lineage-appropriate genes that were silent in pluripotent cells. Using mouse embryonic stem cells (ESCs) with an inducible knockout of Oct1, we showed that Oct1 deficiency resulted in poor induction of mesoderm-specific genes, leading to impaired mesodermal and terminal muscle differentiation. Oct1-deficient cells exhibited poor temporal coordination of the induction of lineage-specific genes and showed inappropriate developmental lineage branching, resulting in poorly differentiated cell states retaining epithelial characteristics. In ESCs, Oct1 localized with the pluripotency factor Oct4 at mesoderm-associated genes and remained bound to those loci during differentiation after the dissociation of Oct4. Binding events for Oct1 overlapped with those for the histone lysine demethylase Utx, and an interaction between Oct1 and Utx suggested that these two proteins cooperate to activate gene expression. The specificity of the ubiquitous Oct1 for the induction of mesodermal genes could be partially explained by the frequent coexistence of Smad and Oct binding sites at mesoderm-specific genes and the cooperative stimulation of mesodermal gene transcription by Oct1 and Smad3. Together, these results identify Oct1 as a key mediator of mesoderm lineage-specific gene induction.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data.

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes.

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

Molecular phylogeny and description of Hygrobates cyrnusensis sp. nov. reveals multiple colonization of Corsica by water mites of the H. fluviatilis-complex (Acariformes, Hydrachnidia, Hygrobatidae)

Water mites of the Hygrobates fluviatilis-complex are common and often the most abundant water mite group in running waters in the Palaearctic. This study provides the first time-calibrated phylogeny of this complex to date. Fossil-calibrated modelling reveals that the most recent common ancestor of the H. fluviatilis-complex originated in the Oligocene. Phylogenetic and dating analyses show the presence of three monophyletic lineages within the H. fluviatilis-complex, with the H. turcicus-clade as the earliest branching lineage. The diversification and speciation processes mainly started from the Oligocene-Miocene around 24 million years ago and continued during the Miocene. In this study, we tried to link time frames of boundary diversification and appearance of certain species with particular paleogeological events. We found evidence of multiple colonization of Corsica in different phases of its geological history, resulting in the presence of two endemic lineages, H. corsicus and H. cyrnusensis sp. nov., the latter species is described by an integrative approach in this study.

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species’ divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species’ pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.