Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phylogenetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

GUN4 appeared early in cyanobacterial evolution

Photosynthesis relies on chlorophylls, which are synthesized via a common tetrapyrrole trunk pathway also leading to heme, vitamin B12, and other pigmented cofactors. The first committed step for chlorophyll biosynthesis is insertion of magnesium into protoporphyrin IX by magnesium chelatase. Magnesium chelatase is composed of H-, I-, and D-subunits, with the tetrapyrrole substrate binding to the H-subunit. This subunit is rapidly inactivated in the presence of substrate, light, and oxygen, so oxygenic photosynthetic organisms require mechanisms to protect magnesium chelatase from similar loss of function. An additional protein, GUN4, binds to the H-subunit and to tetrapyrroles. GUN4 has been proposed to serve this protective role via its ability to bind linear tetrapyrroles (bilins). In the current work, we probe the origins of bilin binding by GUN4 via comparative phylogenetic analysis and biochemical validation of a conserved bilin-binding motif. Based on our results, we propose that bilin-binding GUN4 proteins arose early in cyanobacterial evolution and that this early acquisition represents an ancient adaptation for maintaining chlorophyll biosynthesis in the presence of light and oxygen.

59 BASIC BIOLOGICAL SCIENCES↗

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations.

59 BASIC BIOLOGICAL SCIENCES↗

Unisexual flower initiation in the monoecious Quercus suber L.: a molecular approach

Abstract Several plant species display a temporal separation of the male and female flower organ development to enhance outbreeding; however, little is known regarding the genetic mechanisms controlling this temporal separation. Quercus suber is a monoecious oak tree with accentuated protandry: in late winter, unisexual male flowers emerge adjacent to the swollen buds, whereas unisexual female flowers emerge in the axils of newly formed leaves formed during spring (4–8 weeks after male flowering). Here, a phylogenetic profiling has led to the identification of cork oak homologs of key floral regulatory genes. The role of these cork oak homologs during flower development was identified with functional studies in Arabidopsis thaliana. The expression profile throughout the year of flower regulators (inducers and repressors), in leaves and buds, suggests that the development of male and female flowers may be preceded by separated induction events. Female flowers are most likely induced during the vegetative flush occurring in spring, whereas male flowers may be induced in early summer. Male flowers stay enclosed within the pre-dormant buds, but complete their development before the vegetative flush of the following year, displaying a long period of anthesis that spans the dormant period. Our results portray a genetic mechanism that may explain similar reproductive habits in other monoecious tree species.

Sobral, Rómulo↗

The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

60 APPLIED LIFE SCIENCES↗

Reekeekee- and roodoodooviruses, two different Microviridae clades constituted by the smallest DNA phages

Small circular single-stranded DNA viruses of the Microviridae family are both prevalent and diverse in all ecosystems. They usually harbor a genome between 4.3 and 6.3 kb, with a microvirus recently isolated from a marine Alphaproteobacteria being the smallest known genome of a DNA phage (4.248 kb). A subfamily, Amoyvirinae, has been proposed to classify this virus and other related small Alphaproteobacteria-infecting phages. Here, we report the discovery, in meta-omics data sets from various aquatic ecosystems, of sixteen complete microvirus genomes significantly smaller (2.991–3.692 kb) than known ones. Phylogenetic analysis reveals that these sixteen genomes represent two related, yet distinct and diverse, novel groups of microviruses—amoyviruses being their closest known relatives. We propose that these small microviruses are members of two tentatively named subfamilies Reekeekeevirinae and Roodoodoovirinae. As known microvirus genomes encode many overlapping and overprinted genes that are not identified by gene prediction software, we developed a new methodology to identify all genes based on protein conservation, amino acid composition, and selection pressure estimations. Surprisingly, only four to five genes could be identified per genome, with the number of overprinted genes lower than that in phiX174. These small genomes thus tend to have both a lower number of genes and a shorter length for each gene, leaving no place for variable gene regions that could harbor overprinted genes. Even more surprisingly, these two Microviridae groups had specific and different gene content, and major differences in their conserved protein sequences, highlighting that these two related groups of small genome microviruses use very different strategies to fulfill their lifecycle with such a small number of genes. The discovery of these genomes and the detailed prediction and annotation of their genome content expand our understanding of ssDNA phages in nature and are further evidence that these viruses have explored a wide range of possibilities during their long evolution.

59 BASIC BIOLOGICAL SCIENCES↗

Intra- and inter-subtype HIV diversity between 1994 and 2018 in southern Uganda: a longitudinal population-based study

There is limited data on human immunodeficiency virus (HIV) evolutionary trends in African populations. We evaluated changes in HIV viral diversity and genetic divergence in southern Uganda over a 24-year period spanning the introduction and scale-up of HIV prevention and treatment programs using HIV sequence and survey data from the Rakai Community Cohort Study, an open longitudinal population-based HIV surveillance cohort. Gag (p24) and env (gp41) HIV data were generated from people living with HIV (PLHIV) in 31 inland semi-urban trading and agrarian communities (1994–2018) and four hyperendemic Lake Victoria fishing communities (2011–2018) under continuous surveillance. HIV subtype was assigned using the Recombination Identification Program with phylogenetic confirmation. Inter-subtype diversity was evaluated using the Shannon diversity index, and intra-subtype diversity with the nucleotide diversity and pairwise TN93 genetic distance. Genetic divergence was measured using root-to-tip distance and pairwise TN93 genetic distance analyses. Demographic history of HIV was inferred using a coalescent-based Bayesian Skygrid model. Evolutionary dynamics were assessed among demographic and behavioral population subgroups, including by migration status. 9931 HIV sequences were available from 4999 PLHIV, including 3060 and 1939 persons residing in inland and fishing communities, respectively. In inland communities, subtype A1 viruses proportionately increased from 14.3% in 1995 to 25.9% in 2017 (P < .001), while those of subtype D declined from 73.2% in 1995 to 28.2% in 2017 (P < .001). The proportion of viruses classified as recombinants significantly increased by nearly four-fold from 12.2% in 1995 to 44.8% in 2017. Inter-subtype HIV diversity has generally increased. While intra-subtype p24 genetic diversity and divergence leveled off after 2014, intra-subtype gp41 diversity, effective population size, and divergence increased through 2017. Intra- and inter-subtype viral diversity increased across all demographic and behavioral population subgroups, including among individuals with no recent migration history or extra-community sexual partners. This study provides insights into population-level HIV evolutionary dynamics following the scale-up of HIV prevention and treatment programs. Continued molecular surveillance may provide a better understanding of the dynamics driving population HIV evolution and yield important insights for epidemic control and vaccine development.

60 APPLIED LIFE SCIENCES↗

Seasonal dynamics of core fungi in the switchgrass phyllosphere, and co-occurrence with leaf bacteria

Plant leaves harbor complex microbial communities that influence plant health and productivity. Nevertheless, a detailed understanding of phyllosphere community assembly and drivers is needed, particularly for phyllosphere fungi. Here, we investigated seasonal dynamics of epiphytic phyllosphere fungal communities in switchgrass ( Panicum virgatum L.), a focal bioenergy crop. We also leverage previously published data on switchgrass phyllosphere bacterial communities from the same experimental plants, allowing us to compare fungal and bacterial dynamics and explore inter-Domain network associations in the switchgrass phyllosphere. Overall, we found a strong impact of sampling date on fungal community composition, with multiple taxonomic levels exhibiting clear temporal patterns in relative abundance. In addition, leaf nitrogen concentration, leaf dry matter content, plant height, and minimum daily air temperature explained significant variation in phyllosphere fungal communities, likely due to their correlation with sampling date. Finally, among the core taxa, fungal-bacterial network associations were much more common than bacteria-bacteria associations, suggesting the importance of inter-Domain phylogenetic diversity in microbiome assembly. Although our findings highlight the complexity of phyllosphere microbiome assembly, the clear temporal patterns in lineage-specific fungal abundances give promise to the potential for accurately predicting shifts in fungal phyllosphere communities throughout the growing season, a key research priority for sustainable agriculture.

16S rRNA↗

FUSARIUM-ID v.3.0: An Updated, Downloadable Resource for Fusarium Species Identification

Species within Fusarium are of global agricultural, medical, and food/feed safety concern and have been extensively characterized. However, accurate identification of species is challenging and usually requires DNA sequence data. FUSARIUM-ID ( http://isolate.fusariumdb.org/blast.php ) is a publicly available database designed to support the identification of Fusarium species using sequences of multiple phylogenetically informative loci, especially the highly informative ~680-bp 5' portion of the translation elongation factor 1-alpha (TEF1) gene that has been adopted as the primary barcoding locus in the genus. However, FUSARIUM-ID v.1.0 and 2.0 had several limitations, including inconsistent metadata annotation for the archived sequences and poor representation of some species complexes and marker loci. Here, we present FUSARIUM-ID v.3.0, which provides the following improvements: (i) additional and updated annotation of metadata for isolates associated with each sequence, (ii) expanded taxon representation in the TEF1 sequence database, (iii) availability of the sequence database as a downloadable file to enable local BLAST queries, and (iv) a tutorial file for users to perform local BLAST searches using either freely available software, such as SequenceServer, BLAST+ executable in the command line, and Galaxy, or the proprietary Geneious software. FUSARIUM-ID will be updated on a regular basis by archiving sequences of TEF1 and other loci from newly identified species and greater in-depth sampling of currently recognized species.

Plant Sciences↗

Biomechanical drivers of the evolution of butterflies and moths with a coilable proboscis

Current biomechanical models suggest that butterflies and moths use their proboscis as a drinking straw pulling nectar as a continuous liquid column. Our analyses revealed an alternative mode for fluid uptake: drinking bubble trains that help defeat drag. We combined X-ray phase-contrast imaging, optical video microscopy, micro-computed tomography, phylogenetic models of evolution and fluid mechanics models of bubble-train formation to understand the biomechanics of butterfly and moth feeding. Our models suggest that the bubble-train mechanism appeared in the early evolution of butterflies and moths with a proboscis long enough to coil. We propose that, in addition to the ability to drink a continuous column of fluid from pools, the ability to exploit fluid films by capitalizing on bubble trains would have expanded the range of available food sources, facilitating diversification of Lepidoptera.

Palaoro, Alexandre V.↗

Evolutionary trajectory of transcription factors and selection of targets for metabolic engineering

Transcription factors (TFs) provide potentially powerful tools for plant metabolic engineering as they often control multiple genes in a metabolic pathway. However, selecting the best TF for a particular pathway has been challenging, and the selection often relies significantly on phylogenetic relationships. Here, we offer examples where evolutionary relationships have facilitated the selection of the suitable TFs, alongside situations where such relationships are misleading from the perspective of metabolic engineering. We argue that the evolutionary trajectory of a particular TF might be a better indicator than protein sequence homology alone in helping decide the best targets for plant metabolic engineering efforts. This article is part of the theme issue ‘The evolution of plant metabolism’.

Life Sciences & Biomedicine - Other Topics↗

Turicibacter bilis sp. nov., a novel bacterium isolated from the chicken eggshell and swine ileum

Three novel, anaerobic, Gram-positive bacteria were isolated from the eggshell of two separate white leghorn chicken flocks and the ileum of a healthy pig, and designated MMM721 T , ISU324 and PIG517 respectively. Cells were pleomorphic and capable of forming long chains of rods or coccoid clusters. Phylogenetic analysis of the 16S rRNA gene sequences identified these strains to be within the genus Turicibacter, of which only one species, Turicibacter sanguinis, has been formally described. However, whole genome sequencing of novel isolates returned a digital DNA–DNA hybridization value of 22.5% and average nucleotide identity (ANI) values of 76.4% (ANIb) and 86.0% (ANIm), indicating divergence between the type strain MMM721 T and T. sanguinis, suggesting the strains represented a novel species. The major fatty acid methyl esters of strain MMM721 T were C 16:0 , C 18:1 ω7c and C 18:0 . The strains mainly produced the volatile fatty acid lactate, along with smaller amounts of acetate and butyrate. Together, these data indicate that MMM721 T , along with ISU324 and PIG517, represent a novel species within the genus Turicibacter. We propose the name Turicibacter bilis sp. nov. for the species. The type strain is MMM721 T (=ATCC TSD-238 T =CCUG 74757 T ).

59 BASIC BIOLOGICAL SCIENCES↗

Characterization and rank assignment criteria for the anaerobic fungi (Neocallimastigomycota)

Establishing a solid taxonomic framework is crucial for enabling discovery and documentation efforts. This ensures effective communication between scientists as well as reproducibility of results between laboratories, and facilitates the exchange and preservation of biological material. Such framework can only be achieved by establishing clear criteria for taxa characterization and rank assignment. Within the anaerobic fungi (phylum Neocallimastigomycota), the need for such criteria is especially vital. Difficulties associated with their isolation, maintenance and long-term storage often result in limited availability and loss of previously described taxa. To this end, we provide here a list of morphological, microscopic, phylogenetic and phenotypic criteria for assessment and documentation when characterizing newly obtained Neocallimastigomycota isolates. We also recommend a polyphasic rank-assignment scheme for novel genus-, species- and strain-level designations for newly obtained Neocallimastigomycota isolates.

59 BASIC BIOLOGICAL SCIENCES↗

Proposed minimal standards for description of methanogenic archaea

Methanogenic archaea are a diverse, polyphyletic group of strictly anaerobic prokaryotes capable of producing methane as their primary metabolic product. It has been over three decades since minimal standards for their taxonomic description have been proposed. In light of advancements in technology and amendments in systematic microbiology, revision of the older criteria for taxonomic description is essential. Most of the previously recommended minimum standards regarding phenotypic characterization of pure cultures are maintained. Electron microscopy and chemotaxonomic methods like whole-cell protein and lipid analysis are desirable but not required. Because of advancements in DNA sequencing technologies, obtaining a complete or draft whole genome sequence for type strains and its deposition in a public database are now mandatory. Genomic data should be used for rigorous comparison to close relatives using overall genome related indices such as average nucleotide identity and digital DNA–DNA hybridization. Phylogenetic analysis of the 16S rRNA gene is also required and can be supplemented by phylogenies of the mcrA gene and phylogenomic analysis using multiple conserved, single-copy marker genes. Additionally, it is now established that culture purity is not essential for studying prokaryotes, and description of Candidatus methanogenic taxa using single-cell or metagenomics along with other appropriate criteria is a viable alternative. The revisions to the minimal criteria proposed here by the members of the Subcommittee on the Taxonomy of Methanogenic Archaea of the International Committee on Systematics of Prokaryotes should allow for rigorous yet practical taxonomic description of these important and diverse microbes.

Microbiology↗

Sulfurospirillum diekertiae sp. nov., a tetrachloroethene-respiring bacterium isolated from contaminated soil

Two anaerobic, tetrachloroethene- (PCE-) respiring bacterial isolates, designated strain ACS DCE T and strain ACS TCE , were characterized using a polyphasic approach. Cells were Gram-stain-negative, motile, non-spore-forming and shared a vibrioid- to spirillum-shaped morphology. Optimum growth occurred at 30°C and 0.1–0.4% salinity. The pH range for growth was pH 5.5–7.5, with an optimum at pH 7.2. Hydrogen, formate, pyruvate and lactate as electron donors supported respiratory reductive dechlorination of PCE to cis-1,2-dichloroethene (cDCE) in strain ACS DCE T and of PCE to trichloroethene (TCE) in strain ACS TCE . Both strains were able to grow with pyruvate under microaerobic conditions. Nitrate, elemental sulphur, and thiosulphate were alternative electron acceptors. Autotrophic growth was not observed and acetate served as carbon source for both strains. The major cellular fatty acids were C 16:1 ω7c, C 16:0 , C 14:0 and C 18:1 ω7c. Both genomes feature a circular plasmid. Strains ACS DCE T and ACS TCE were previously assigned to the candidate species 'Sulfurospirillum acididehalogenans'. Here, based on key genomic features and pairwise comparisons of whole-genome sequences, including average nucleotide identity, digital DNA–DNA hybridization and average amino acid identity, strains ACS DCE T and ACS TCE , 'Ca. Sulfurospirillum diekertiae' strains SL2-1 and SL2-2, and the unclassified Sulfurospirillum sp. strain SPD-1 are grouped into one distinct species separate from previously described Sulfurospirillum species. Compared to Sulfurospirillum multivorans and Sulfurospirillum halorespirans, which dechlorinate PCE to cDCE without substantial TCE accumulation, these five strains produce TCE or cDCE as the end product. In addition, some cellular fatty acids (e.g., C 16:0 3OH, C 17:0 iso 3OH, C 17:0 2OH) were detected in strains ACS DCE T and ACS TCE but not in other Sulfurospirillum species. On the basis of phylogenetic, physiological and phenotypic characteristics, 'Ca. Sulfurospirillum acididehalogenans' and 'Ca. Sulfurospirillum diekertiae' are proposed to be merged into one novel species within the genus Sulfurospirillum, for which the name Sulfurospirillum diekertiae sp. nov. is proposed. Finally, the type strain is ACS DCE T (=JCM 33349 T = KCTC 15819 T =CGMCC 1.5292 T ).

59 BASIC BIOLOGICAL SCIENCES↗

Saccharomycopsis praedatoria sp. nov., a predacious yeast isolated from soil and rotten wood in an Amazonian rainforest biome

Three yeast isolates were obtained from soil and rotting wood samples collected in an Amazonian rainforest biome in Brazil. Comparison of the intergenic spacer 5.8S region and the D1/D2 domains of the large subunit rRNA gene showed that the isolates represent a novel species of the genus Saccharomycopsis. A tree inferred from the D1/D2 sequences placed the novel species near a subclade containing Saccharomycopsis lassenensis, Saccharomycopsis fermentans, Saccharomycopsis javanensis, Saccharomycopsis babjevae, Saccharomycopsis schoenii and Saccharomycopsis oosterbeekiorum, but with low bootstrap support. In terms of sequence divergence, the novel species had the highest identity in the D1/D2 domains with Saccharomycopsis capsularis, from which it differed by 36 substitutions. In contrast, a phylogenomic analysis based on 1061 single-copy orthologs for a smaller set of Saccharomycopsis species whose whole genome sequences are available indicated that the novel species represented by strain UFMG-CM-Y6991 is phylogenetically closer to Saccharomycopsis fodiens and Saccharomycopsis sp. TF2021a (=Saccharomycopsis phalluae). The novel yeast is homothallic and produces asci with one spheroidal ascospore with an equatorial or subequatorial ledge. The name Saccharomycopsis praedatoria sp. nov. is proposed to accommodate the novel species. The holotype of Saccharomycopsis praedatoria is CBS 16589 T . The MycoBank number is MB849369. S. praedatoria was able to kill cells of Saccharomyces cerevisiae by means of penetration with infection pegs, a trait common to most species of Saccharomycopsis.

60 APPLIED LIFE SCIENCES↗

Spirochaete genome identified in red abalone sample represents a novel genus Candidatus Haliotispira gen. nov. within the order Spirochaetales

A fully assembled spirochaete genome was identified as a contaminating scaffold in our red abalone ( Haliotis rufescens ) genome assembly. In this paper, we describe the analysis of this bacterial genome. The assembled spirochaete genome is 3.25 Mb in size with 48.5 mol% G+C content. The proteomes of 38 species were compared with the spirochaete genome and it was discovered to form an independent branch within the family Spirochaetaceae on the phylogenetic tree. The comparison of 16S rRNA sequences and average nucleotide identity scores between the spirochaete genome with known species of different families in Spirochaetia indicate that it is an unknown species. Further, the percentage of conserved proteins compared to neighbouring taxa confirm that it does not belong to a known genus within Spirochaetaceae . We propose the name Candidatus Haliotispira prima gen. nov., sp. nov. based on its taxonomic placement and origin. We also tested for the presence of this species in different species of abalone and found that it is also present in white abalone ( Haliotis sorenseni ). In addition, we highlight the need for better classification of taxa within the class Spirochaetia .

Microbiology↗

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

59 BASIC BIOLOGICAL SCIENCES↗

Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals

Cycads are known to host symbiotic cyanobacteria, includingNostocalesspecies, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition toNostocspp., symbiotic cyanobacteria of the genusAulosirawere identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includesAulosiraspp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Genetics & Heredity↗