Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals.

Cycads are known to host symbiotic cyanobacteria, including Nostocales species, 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 to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were 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 includes Aulosira spp. 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.

Genetic analysis and QTL mapping of domestication-related traits in chili pepper (Capsicum annuum L.).

Chili pepper (Capsicum annuum L.) is one of the oldest and most phenotypically diverse pre-Columbian crops of the Americas. Despite the abundance of genetic resources, the use of wild germplasm and landraces in chili pepper breeding is limited. A better understanding of the evolutionary history in chili peppers, particularly in the context of traits of agronomic interest, can contribute to future improvement and conservation of genetic resources. In this study, an F2:3 mapping population derived from a cross between a C. annuum wild accession (Chiltepin) and a cultivated variety (Puya) was used to identify genomic regions associated with 19 domestication and agronomic traits. A genetic map was constructed consisting of 1023 single nucleotide polymorphism (SNP) markers clustered into 12 linkage groups and spanning a total of 1,263.87 cM. A reciprocal translocation that differentiates the domesticated genome from its wild ancestor and other related species was identified between chromosomes 1 and 8. Quantitative trait locus (QTL) analysis detected 20 marker-trait associations for 13 phenotypes, from which 14 corresponded to previously identified loci, and six were novel genomic regions related to previously unexplored domestication-syndrome traits, including form of unripe fruit, seedlessness, deciduous fruit, and growth habit. Our results revealed that the genetic architecture of Capsicum domestication is similar to other domesticated species with few loci with large effects, the presence of QTLs clusters in different genomic regions, and the predominance of domesticated recessive alleles. Our analysis indicates the domestication process in chili pepper has also had an effect on traits not directly related to the domestication syndrome. The information obtained in this study provides a more complete understanding of the genetic basis of Capsicum domestication that can potentially guide strategies for the exploitation of wild alleles.

The role of non-additive gene action on gene expression variation in plant domestication.

BACKGROUND: Plant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs. RESULTS: We compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene's expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action. CONCLUSIONS: Our study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Microbial Diversity in Cultivated and Feral Vanilla Vanilla planifolia Orchids Affected by Stem and Rot Disease.

The worldwide production of vanilla, a native orchid from Mexico, is greatly affected by stem and root rot disease (SRD), typically associated with Fusarium oxysporum fungi. We hypothesized that the presence of Fusarium species in vanilla is not sufficient for the plant to express symptoms of the disease. We described the taxonomic composition of endophytic microbiomes in symptomatic and asymptomatic vanilla plants using 16S and ITS rDNA metabarcoding, and ITS Sanger sequences generated from fungal isolates. We compared the bacterial and fungal diversity in vanilla plants from a long-term plantation, and from feral plants found near abandoned plantations that did not present SRD symptoms. No significant differences were found in the species richness of the bacterial and fungal microbiome among feral, or asymptomatic and symptomatic cultivated vanilla. However, significant differences were detected in both fungal and bacterial diversity from different organs in the same plant, with roots being more diverse than stems. We found that Proteobacteria and Actinobacteria, as well as the fungal families Nectriaceae and Xylariaceae, constitute the core of the vanilla microbiome that inhabits the root and stem of both cultivated and feral plants. Our work provides information on the microbial diversity associated to root and stem rot in vanilla and lays the groundwork for a better understanding of the role of the microbiome in vanilla fungal diseases.

Phylogenomics and population genomics of SARS-CoV-2 in Mexico during the pre-vaccination stage reveals variants of interest B.1.1.28.4 and B.1.1.222 or B.1.1.519 and the nucleocapsid mutation S194L associated with symptoms.

Understanding the evolution of the SARS-CoV-2 virus in various regions of the world during the Covid-19 pandemic is essential to help mitigate the effects of this devastating disease. We describe the phylogenomic and population genetic patterns of the virus in Mexico during the pre-vaccination stage, including asymptomatic carriers. A real-time quantitative PCR screening and phylogenomic reconstructions directed at sequence/structure analysis of the spike glycoprotein revealed mutation of concern E484K in genomes from central Mexico, in addition to the nationwide prevalence of the imported variant 20C/S:452R (B.1.427/9). Overall, the detected variants in Mexico show spike protein mutations in the N-terminal domain (i.e. R190M), in the receptor-binding motif (i.e. T478K, E484K), within the S1-S2 subdomains (i.e. P681R/H, T732A), and at the basis of the protein, V1176F, raising concerns about the lack of phenotypic and clinical data available for the variants of interest we postulate: 20B/478K.V1 (B.1.1.222 or B.1.1.519) and 20B/P.4 (B.1.1.28.4). Moreover, the population patterns of single nucleotide variants from symptomatic and asymptomatic carriers obtained with a self-sampling scheme confirmed the presence of several fixed variants, and differences in allelic frequencies among localities. We identified the mutation N:S194L of the nucleocapsid protein associated with symptomatic patients. Phylogenetically, this mutation is frequent in Mexican sub-clades. Our results highlight the dual and complementary role of spike and nucleocapsid proteins in adaptive evolution of SARS-CoV-2 to their hosts and provide a baseline for specific follow-up of mutations of concern during the vaccination stage.

Phylogenomics of Salvia L. subgenus Calosphace (Lamiaceae).

The evolutionary relationships of Salvia have been difficult to estimate. In this study, we used the Next Generation Sequencing method Hyb-Seq to evaluate relationships among 90 Lamiaceae samples, including representatives of Mentheae, Ocimeae, Salvia subgenera Audibertia, Leonia, Salvia, and 69 species of subgenus Calosphace, representing 32 of Epling's sections. A bait set was designed in MarkerMiner using available transcriptome data to enrich 119 variable nuclear loci. Nuclear and chloroplast loci were assembled with hybphylomaker (HPM), followed by coalescent approach analyses for nuclear data (ASTRAL, BEAST) and a concatenated Maximum Likelihood analysis of chloroplast loci. The HPM assembly had an average of 1,314,368 mapped reads for the sample and 527 putative exons. Phylogenetic inferences resolved strongly supported relationships for the deep-level nodes, agreeing with previous hypotheses which assumed that subgenus Audibertia is sister to subgenus Calosphace. Within subgenus Calosphace, we recovered eight monophyletic sections sensu Epling, Cardinalis, Hastatae, Incarnatae, and Uricae in all the analyses (nDNA and cpDNA), Biflorae, Lavanduloideae, and Sigmoideae in nuclear analyses (ASTRAL, BEAST) and Curtiflorae in ASTRAL trees. Network analysis supports deep node relationships, some of the main clades, and recovers reticulation within the core Calosphace. The chloroplast phylogeny resolved deep nodes and four monophyletic Calosphace sections. Placement of S. axillaris is distinct in nuclear evidence and chloroplast, as sister to the rest of the S. subg. Calosphace in chloroplast and a clade with "Hastatae clade" sister to the rest of the subgenus in nuclear evidence. We also tested the monophyly of S. hispanica, S. polystachia, S. purpurea, and S. tiliifolia, including two samples of each, and found that S. hispanica and S. purpurea are monophyletic. Our baits can be used in future studies of Lamiaceae phylogeny to estimate relationships between genera and among species. In this study, we presented a Hyb-Seq phylogeny for complex, recently diverged Salvia, which could be implemented in other Lamiaceae.

Comparative transcriptome analysis of cultivated and wild seeds of Salvia hispanica (chia).

Salvia hispanica (chia) constituted an important crop for pre-Columbian civilizations and is considered a superfood for its rich content of essential fatty acids and proteins. In this study, we performed the first comprehensive comparative transcriptome analysis between seeds from cultivated varieties and from accessions collected from native wild populations in Mexico. From the 69,873 annotated transcripts assembled de novo, enriched functional categories and pathways revealed that the lipid metabolism was one of the most activated processes. Expression changes were detected among wild and cultivated groups and among growth conditions in transcripts responsible for triacylglycerol and fatty acid synthesis and degradation. We also quantified storage protein fractions that revealed variation concerning nutraceutical proteins such as albumin and glutelin. Genetic diversity estimated with 23,641 single nucleotide polymorphisms (SNPs) revealed that most of the variation remains in the wild populations, and that a wild-type cultivated variety is genetically related to wild accessions. Additionally, we reported 202 simple sequence repeat (SSRs) markers useful for population genetic studies. Overall, we provided transcript variation that can be used for breeding programs to further develop chia varieties with enhanced nutraceutical traits and tools to explore the genetic diversity and history of this rediscovered plant.

Unlocking a high bacterial diversity in the coralloid root microbiome from the cycad genus Dioon.

Cycads are among the few plants that have developed specialized roots to host nitrogen-fixing bacteria. We describe the bacterial diversity of the coralloid roots from seven Dioon species and their surrounding rhizosphere and soil. Using 16S rRNA gene amplicon sequencing, we found that all coralloid roots are inhabited by a broad diversity of bacterial groups, including cyanobacteria and Rhizobiales among the most abundant groups. The diversity and composition of the endophytes are similar in the six Mexican species of Dioon that we evaluated, suggesting a recent divergence of Dioon populations and/or similar plant-driven restrictions in maintaining the coralloid root microbiome. Botanical garden samples and natural populations have a similar taxonomic composition, although the beta diversity differed between these populations. The rhizosphere surrounding the coralloid root serves as a reservoir and source of mostly diazotroph and plant growth-promoting groups that colonize the coralloid endosphere. In the case of cyanobacteria, the endosphere is enriched with Nostoc spp and Calothrix spp that are closely related to previously reported symbiont genera in cycads and other early divergent plants. The data reported here provide an in-depth taxonomic characterization of the bacterial community associated with coralloid root microbiome. The functional aspects of the endophytes, their biological interactions, and their evolutionary history are the next research step in this recently discovered diversity within the cycad coralloid root microbiome.

Ethnobotany of Mexican and northern Central American cycads (Zamiaceae).

BACKGROUND: This study documents cycad-human relationships in Mexico, Belize, Guatemala, El Salvador, and Honduras over the last 6000 years. The impetus was acute need for a better understanding of previously undocumented uses of cycads in this region, and the need to improve cycad conservation strategies using ethnobotanical data. We hypothesized that cycads are significant dietary items with no long-term neurological effects, are important to religious practice, and contribute to cultural identity and sense of place, but that traditional knowledge and uses are rapidly eroding. Guiding questions focused on nomenclature, food and toxicity, relationships to palms and maize, land management issues, roles in religious ceremony, and medicinal uses, among others, and contributions of these to preservation of cycads. METHODS: From 2000 to 2017, the authors conducted 411 semi-structured ethnographic interviews, engaged in participant-observation in Mexican and Honduran communities, and carried out archival research and literature surveys. RESULTS: We documented 235 terms and associated uses that 28 ethnic groups have for 57 species in 19 languages across 21 Mexican states and 4 Central American nations. Carbohydrate-rich cycads have been both famine foods and staples for at least six millennia across the region and are still consumed in Mexico and Honduras. Certain parts are eaten without removing toxins, while seed and stem starches are detoxified via several complex processes. Leaves are incorporated into syncretic Roman Catholic-Mesoamerican religious ceremonies such as pilgrimages, Easter Week, and Day of the Dead. Cycads are often perceived as ancestors and protectors of maize, revealing a close relationship between both groups. Certain beliefs and practices give cycads prominent roles in conceptions of sense of place and cultural heritage. CONCLUSIONS: Cycads are still used as foods in many places. Though they do not appear to cause long-term neurological damage, their health effects are not fully understood. They are often important to religion and contribute to cultural identity and sense of place. However, because most traditional knowledge and uses are rapidly eroding, new community-based biocultural conservation efforts are needed. These should incorporate tradition where possible and seek inspiration from existing successful cases in Honduras and Mexico.

Increasing Metagenomic Resolution of Microbiome Interactions Through Functional Phylogenomics and Bacterial Sub-Communities.

The genomic composition of the microbiome and its relationship with the environment is an exciting open question in biology. Metagenomics is a useful tool in the discovery of previously unknown taxa, but its use to understand the functional and ecological capacities of the microbiome is limited until taxonomy and function are understood in the context of the community. We suggest that this can be achieved using a combined functional phylogenomics and co-culture-based experimental strategy that can increase our capacity to measure sub-community interactions. Functional phylogenomics can identify and partition the genome such that hidden gene functions and gene clusters with unique evolutionary signals are revealed. We can test these phylogenomic predictions using an experimental model based on sub-community populations that represent a subset of the diversity directly obtained from environmental samples. These populations increase the detection of mechanisms that drive functional forces in the assembly of the microbiome, in particular the role of metabolites from key taxa in community interactions. Our combined approach leverages the potential of metagenomics to address biological questions from ecological systems.

Genetic variation in avocado stem weevils Copturus aguacatae (Coleoptera: Curculionidae) in Mexico.

BACKGROUND AND AIM: The avocado stem weevil Copturus aguacatae is an important pest in avocado plantations. Its presence hinders the production and marketing of avocado in Mexico, the largest avocado producer worldwide. Biological control through pheromone synthesis, a strategy favored over chemical control in crops, is currently limited by difficult field identification of this species. MATERIALS AND METHODS: Using DNA barcoding, we examine the patterns of genetic variation of C. aguacatae in avocado trees in Mexico to help facilitate its identification and biological control. RESULTS: We show that there is one single species of avocado stem weevil throughout the sampled sites in Mexico. Overall, haplotype diversity is high, with Oaxaca forming one distinct group and all other sampled populations are admixed irrespective of geographic origin. CONCLUSION: The results suggest that high gene flow is maintained in this species and that a global strategy for biocontrol can be designed and implemented throughout the sampled range.

Population genetics of the understory fishtail palm Chamaedorea ernesti-augusti in Belize: high genetic connectivity with local differentiation.

BACKGROUND: Developing a greater understanding of population genetic structure in lowland tropical plant species is highly relevant to our knowledge of increasingly fragmented forests and to the conservation of threatened species. Specific studies are particularly needed for taxa whose population dynamics are further impacted by human harvesting practices. One such case is the fishtail or xaté palm (Chamaedorea ernesti-augusti) of Central America, whose wild-collected leaves are becoming progressively more important to the global ornamental industry. We use microsatellite markers to describe the population genetics of this species in Belize and test the effects of climate change and deforestation on its recent and historical effective population size. RESULTS: We found high levels of inbreeding coupled with moderate or high allelic diversity within populations. Overall high gene flow was observed, with a north and south gradient and ongoing differentiation at smaller spatial scales. Immigration rates among populations were more difficult to discern, with minimal evidence for isolation by distance. We infer a tenfold reduction in effective population size ca. 10,000 years ago, but fail to detect changes attributable to Mayan or contemporary deforestation. CONCLUSION: Populations of C. ernesti-augusti are genetically heterogeneous demes at a local spatial scale, but are widely connected at a regional level in Belize. We suggest that the inferred patterns in population genetic structure are the result of the colonization of this species into Belize following expansion of humid forests in combination with demographic and mating patterns. Within populations, we hypothesize that low aggregated population density over large areas, short distance pollen dispersal via thrips, low adult survival, and low fruiting combined with early flowering may contribute towards local inbreeding via genetic drift. Relatively high levels of regional connectivity are likely the result of animal-mediated long-distance seed dispersal. The greatest present threat to the species is the potential onset of inbreeding depression as the result of increased human harvesting activities. Future genetic studies in understory palms should focus on both fine-scale and landscape-level genetic structure.

PERMANENT GENETIC RESOURCES: Development of microsatellite markers of the Mexican understorey palm Chamaedorea elegans, cross-species genotyping, and amplification in congeners.

With striking morphological diversity and adaptability, Chamaedorea palms constitute an ecologically and economically important understorey component of Neotropical forests. Nine loci developed for Chamaedorea elegans evaluated in three Veracruz populations resulted in a large number of alleles (8-18), and high expected heterozygosity (0.49-0.92), but low observed (0.27-0.65) heterozygosity. Deviations from Hardy-Weinberg and high inbreeding suggest a lack of panmixia. Eight loci optimized for Chamaedorea ernesti-augustii showed similar patterns of variation. All nine multiplexing loci amplified in other five congeneric species, which will facilitate comparisons within the genus and contribute to the conservation of its genetic resources.