Theoretical study of ArcB and its dimerization, interaction with anaerobic metabolites, and activation of ArcA.

The complex metabolism of Escherichia coli has been extensively studied, including its response to oxygen availability. The ArcA/B two-component system (TCS) is the key regulator for the transition between these two environmental conditions and has been thoroughly characterized using genetic and biochemical approaches. Still, to date, limited structural data is available. The breakthrough provided by AlphaFold2 in 2021 has brought a reliable tool to the scientific community for assessing the structural features of complex proteins. In this report, we analyzed the structural aspects of the ArcA/B TCS using AlphaFold2 models. The models are consistent with the experimentally determined structures of ArcB kinase. The predicted structure of the dimeric form of ArcB is consistent with the extensive genetic and biochemical data available regarding mechanistic signal perception and regulation. The predicted interaction of the dimeric form of ArcB with its cognate response regulator (ArcA) is also consistent with both the forward and reverse phosphotransfer mechanisms. The ArcB model was used to detect putative binding cavities to anaerobic metabolites, encouraging testing of these predictions experimentally. Finally, the highly accurate models of other ArcB homologs suggest that different experimental approaches are needed to determine signal perception in kinases lacking the PAS domain. Overall, ArcB is a kinase with features that need further testing, especially in determining its crystal structure under different conditions.

Two-Component System Sensor Kinases from Asgardian Archaea May Be Witnesses to Eukaryotic Cell Evolution.

The signal transduction paradigm in bacteria involves two-component systems (TCSs). Asgardarchaeota are archaea that may have originated the current eukaryotic lifeforms. Most research on these archaea has focused on eukaryotic-like features, such as genes involved in phagocytosis, cytoskeleton structure, and vesicle trafficking. However, little attention has been given to specific prokaryotic features. Here, the sequence and predicted structural features of TCS sensor kinases analyzed from two metagenome assemblies and a genomic assembly from cultured Asgardian archaea are presented. The homology of the sensor kinases suggests the grouping of Lokiarchaeum closer to bacterial homologs. In contrast, one group from a Lokiarchaeum and a meta-genome assembly from Candidatus Heimdallarchaeum suggest the presence of a set of kinases separated from the typical bacterial TCS sensor kinases. AtoS and ArcB homologs were found in meta-genome assemblies along with defined domains for other well-characterized sensor kinases, suggesting the close link between these organisms and bacteria that may have resulted in the metabolic link to the establishment of symbiosis. Several kinases are predicted to be cytoplasmic; some contain several PAS domains. The data shown here suggest that TCS kinases in Asgardian bacteria are witnesses to the transition from bacteria to eukaryotic organisms.

Catalases in the pathogenesis of Sporothrix schenckii research.

Pathogenic fungal infection success depends on the ability to escape the immune response. Most strategies for fungal infection control are focused on the inhibition of virulence factors and increasing the effectiveness of antifungal drugs. Nevertheless, little attention has been focused on their physiological resistance to the host immune system. Hints may be found in pathogenic fungi that also inhabit the soil. In nature, the saprophyte lifestyle of fungi is also associated with predators that can induce oxidative stress upon cell damage. The natural sources of nutrients for fungi are linked to cellulose degradation, which in turn generates reactive oxygen species (ROS). Overall, the antioxidant arsenal needed to thrive both in free-living and pathogenic lifestyles in fungi is fundamental for success. In this review, we present recent findings regarding catalases and oxidative stress in fungi and how these can be in close relationship with pathogenesis. Additionally, special focus is placed on catalases of Sporothrix schenckii as a pathogenic model with a dual lifestyle. It is assumed that catalase expression is activated upon exposure to H2O2, but there are reports where this is not always the case. Additionally, it may be relevant to consider the role of catalases in S. schenckii survival in the saprophytic lifestyle and why their study can assess their involvement in the survival and therefore, in the virulence phenotype of different species of Sporothrix and when each of the three catalases are required. Also, studying antioxidant mechanisms in other isolates of pathogenic and free-living fungi may be linked to the virulence phenotype and be potential therapeutic and diagnostic targets. Thus, the rationale for this review to place focus on fungal catalases and their role in pathogenesis in addition to counteracting the effect of immune system reactive oxygen species. Fungi that thrive in soil and have mammal hosts could shed light on the importance of these enzymes in the two types of lifestyles. We look forward to encouraging more research in a myriad of areas on catalase biology with a focus on basic and applied objectives and placing these enzymes as virulence determinants.

Escherichia coli transcription factors of unknown function: sequence features and possible evolutionary relationships.

Organisms need mechanisms to perceive the environment and respond accordingly to environmental changes or the presence of hazards. Transcription factors (TFs) are required for cells to respond to the environment by controlling the expression of genes needed. Escherichia coli has been the model bacterium for many decades, and still, there are features embedded in its genome that remain unstudied. To date, 58 TFs remain poorly characterized, although their binding sites have been experimentally determined. This study showed that these TFs have sequence variation at the third codon position G+C content but maintain the same Codon Adaptation Index (CAI) trend as annotated functional transcription factors. Most of these transcription factors are in areas of the genome where abundant repetitive and mobile elements are present. Sequence divergence points to groups with distinctive sequence signatures but maintaining the same type of DNA binding domain. Finally, the analysis of the promoter sequences of the 58 TFs showed A+T rich regions that agree with the features of horizontally transferred genes. The findings reported here pave the way for future research of these TFs that may uncover their role as spare factors in case of lose-of-function mutations in core TFs and trace back their evolutionary history.

Gene Regulatory Network Inference and Gene Module Regulating Virulence in Fusarium oxysporum.

In this work, we inferred the gene regulatory network (GRN) of the fungus Fusarium oxysporum by using the regulatory networks of Aspergillus nidulans FGSC A4, Neurospora crassa OR74A, Saccharomyces cerevisiae S288c, and Fusarium graminearum PH-1 as templates for sequence comparisons. Topological properties to infer the role of transcription factors (TFs) and to identify functional modules were calculated in the GRN. From these analyzes, five TFs were identified as hubs, including FOXG_04688 and FOXG_05432, which regulate 2,404 and 1,864 target genes, respectively. In addition, 16 communities were identified in the GRN, where the largest contains 1,923 genes and the smallest contains 227 genes. Finally, the genes associated with virulence were extracted from the GRN and exhaustively analyzed, and we identified a giant module with ten TFs and 273 target genes, where the most highly connected node corresponds to the transcription factor FOXG_05265, homologous to the putative bZip transcription factor CPTF1 of Claviceps purpurea, which is involved in ergotism disease that affects cereal crops and grasses. The results described in this work can be used for the study of gene regulation in this organism and open the possibility to explore putative genes associated with virulence against their host.

The Search for Cryptic L-Rhamnosyltransferases on the Sporothrix schenckii Genome.

The fungal cell wall is an attractive structure to look for new antifungal drug targets and for understanding the host-fungus interaction. Sporothrix schenckii is one of the main causative agents of both human and animal sporotrichosis and currently is the species most studied of the Sporothrix genus. The cell wall of this organism has been previously analyzed, and rhamnoconjugates are signature molecules found on the surface of both mycelia and yeast-like cells. Similar to other reactions where sugars are covalently linked to other sugars, lipids, or proteins, the rhamnosylation process in this organism is expected to involve glycosyltransferases with the ability to transfer rhamnose from a sugar donor to the acceptor molecule, i.e., rhamnosyltransferases. However, no obvious rhamnosyltransferase has thus far been identified within the S. schenckii proteome or genome. Here, using a Hidden Markov Model profile strategy, we found within the S. schenckii genome five putative genes encoding for rhamnosyltransferases. Expression analyses indicated that only two of them, named RHT1 and RHT2, were significantly expressed in yeast-like cells and during interaction with the host. These two genes were heterologously expressed in Escherichia coli, and the purified recombinant proteins showed rhamnosyltransferase activity, dependent on the presence of UDP-rhamnose as a sugar donor. To the best of our knowledge, this is the first report about rhamnosyltransferases in S. schenckii.

Universal Primers for Detection of Novel Plant Capsid-Less Viruses: Papaya Umbra-like Viruses as Example.

For diagnosis of positive-sense single-stranded RNA viruses, primers are usually raised against the sequence encoding capsid proteins, since structural proteins are more conserved. This chapter focuses on the design of primers for a group of novel viruses lacking a capsid, known as papaya Umbra-like viruses (unassigned genus) associated with Papaya Sticky Disease, which represent a threat to papaya production. Based on sequence alignments of a region encoding the RNA-dependent RNA Polymerase, universal primers to detect all the known viruses from four countries are proposed. The Forward universal primer can be used in combination with clade- and subclade-specific primers for rapid virus identification. We walk the reader through downloading sequences from nucleotide databases, doing sequence alignments and phylogenetic tree construction to identify conserved and variable regions as valid primer targets; we also show how to design and analyze the primers.

Re-establishment of the genus Pseudalbizzia (Leguminosae, Caesalpinioideae, mimosoid clade): the New World species formerly placed in Albizia.

Following recent mimosoid phylogenetic and phylogenomic studies demonstrating the non-monophyly of the genus Albizia, we present a new molecular phylogeny focused on the neotropical species in the genus, with much denser taxon sampling than previous studies. Our aims were to test the monophyly of the neotropical section Arthrosamanea, resolve species relationships, and gain insights into the evolution of fruit morphology. We perform a Bayesian phylogenetic analysis of sequences of nuclear internal and external transcribed spacer regions and trace the evolution of fruit dehiscence and lomentiform pods. Our results find further support for the non-monophyly of the genus Albizia, and confirm the previously proposed segregation of Hesperalbizia, Hydrochorea, Balizia and Pseudosamanea. All species that were sampled from section Arthrosamanea form a clade that is sister to a clade composed of Jupunba, Punjuba, Balizia and Hydrochorea. We find that lomentiform fruits are independently derived from indehiscent septate fruits in both Hydrochorea and section Arthrosamanea. Our results show that morphological adaptations to hydrochory, associated with shifts into seasonally flooded habitats, have occurred several times independently in different geographic areas and different lineages within the ingoid clade. This suggests that environmental conditions have likely played a key role in the evolution of fruit types in Albizia and related genera. We resurrect the name Pseudalbizzia to accommodate the species of section Arthrosamanea, except for two species that were not sampled here but have been shown in other studies to be more closely related to other ingoid genera and we restrict the name Albizia s.s. to the species from Africa, Madagascar, Asia, Australia, and the Pacific. Twenty-one new nomenclatural combinations in Pseudalbizzia are proposed, including 16 species and 5 infraspecific varietal names. In addition to the type species Pseudalbizziaberteroana, the genus has 17 species distributed across tropical regions of the Americas, including the Caribbean. Finally, a new infrageneric classification into five sections is proposed and a distribution map of the species of Pseudalbizzia is presented.

Auxin perception in Agave is dependent on the species' Auxin Response Factors.

Auxins are one of the most important and studied phytohormones in nature. Auxin signaling and perception take place in the cytosol, where the auxin is sensed. Then, in the nucleus, the auxin response factors (ARF) promote the expression of early-response genes. It is well known that not all plants respond to the same amount and type of auxins and that the response can be very different even among plants of the same species, as we present here. Here we investigate the behavior of ARF in response to various auxins in Agave angustifolia Haw., A. fourcroydes Lem. and A. tequilana Weber var. Azul. By screening the available database of A. tequilana genes, we have identified 32 ARF genes with high sequence identity in the conserved domains, grouped into three main clades. A phylogenetic tree was inferred from alignments of the 32 Agave ARF protein sequences and the evolutionary relationship with other species was analyzed. AteqARF 4, 15, 21, and 29 were selected as a representative diverse sample coming from each of the different subclades that comprise the two main clades of the inferred phylogenetic reconstruction. These ARFs showed differential species-specific expression patterns in the presence of indole-3-acetic acid (IAA) and 2,4-dichlorophenoxyacetic acid (2,4-D). Interestingly, A. angustifolia showed different phenotypes in the presence and absence of auxins. In the absence of auxin, A. angustifolia produces roots, while shoots are developed in the presence of IAA. However, in the presence of 2,4-D, the plant meristem converts into callus. According to our results, it is likely that AteqARF15 participates in this outcome.

Data mining of metagenomes to find novel enzymes: a non-computationally intensive method.

Currently, there is a need of non-computationally-intensive bioinformatics tools to cope with the increase of large datasets produced by Next Generation Sequencing technologies. We present a simple and robust bioinformatics pipeline to search for novel enzymes in metagenomic sequences. The strategy is based on pattern searching using as reference conserved motifs coded as regular expressions. As a case study, we applied this scheme to search for novel proteases S8A in a publicly available metagenome. Briefly, (1) the metagenome was assembled and translated into amino acids; (2) patterns were matched using regular expressions; (3) retrieved sequences were annotated; and (4) diversity analyses were conducted. Following this pipeline, we were able to identify nine sequences containing an S8 catalytic triad, starting from a metagenome containing 9,921,136 Illumina reads. Identity of these nine sequences was confirmed by BLASTp against databases at NCBI and MEROPS. Identities ranged from 62 to 89% to their respective nearest ortholog, which belonged to phyla Proteobacteria, Actinobacteria, Planctomycetes, Bacterioidetes, and Cyanobacteria, consistent with the most abundant phyla reported for this metagenome. All these results support the idea that they all are novel S8 sequences and strongly suggest that our methodology is robust and suitable to detect novel enzymes.

A Lineage of Begomoviruses Encode Rep and AC4 Proteins of Enigmatic Ancestry: Hints on the Evolution of Geminiviruses in the New World.

The begomoviruses (BGVs) are plant pathogens that evolved in the Old World during the Cretaceous and arrived to the New World (NW) in the Cenozoic era. A subgroup of NW BGVs, the "Squash leaf curl virus (SLCV) lineage" (S-Lin), includes viruses with unique characteristics. To get clues on the evolutionary origin of this lineage, a search for divergent members was undertaken. Four novel BGVs were characterized, including one that is basal to the group. Comparative analyses led to discover a ~670 bp genome module that is nearly exclusive of this lineage, encompassing the replication origin, the AC4 gene, and 480 bp of the Rep gene. A similar DNA module was found in two curtoviruses, hence suggesting that the S-Lin ancestor acquired its distinctive genomic segment by recombination with a curtovirus. This hypothesis was definitely disproved by an in-depth sequence analysis. The search for homologs of S-Lin Rep uncover the common origin of Rep proteins encoded by diverse Geminiviridae genera and viral "fossils" integrated at plant genomes. In contrast, no homolog of S-Lin Rep was found in public databases. Consequently, it was concluded that the SLCV clade ancestor evolved by a recombination event between a primitive NW BGV and a virus from a hitherto unknown lineage.

Identification of functional signatures in the metabolism of the three cellular domains of life.

In order to identify common and specific enzymatic activities associated with the metabolism of the three cellular domains of life, the conservation and variations between the enzyme contents of Bacteria, Archaea, and Eukarya organisms were evaluated. To this end, the content of enzymes belonging to a particular pathway and their abundance and distribution in 1507 organisms that have been annotated and deposited in the KEGG database were assessed. In addition, we evaluated the consecutive enzymatic reaction pairs obtained from metabolic pathway reactions and transformed into sequences of enzymatic reactions, with catalytic activities encoded in the Enzyme Commission numbers, which are linked by a substrate. Both analyses are complementary: the first considers individual reactions associated with each organism and metabolic map, and the second evaluates the functional associations between pairs of consecutive reactions. From these comparisons, we found a set of five enzymatic reactions that were widely distributed in all the organisms and considered here as universal to Bacteria, Archaea, and Eukarya; whereas 132 pairs out of 3151 reactions were identified as significant, only 5 of them were found to be widely distributed in all the taxonomic divisions. However, these universal reactions are not widely distributed along the metabolic maps, suggesting their dispensability to all metabolic processes. Finally, we found that universal reactions are also associated with ancestral domains, such as those related to phosphorus-containing groups with a phosphate group as acceptor or those related to the ribulose-phosphate binding barrel, triosephosphate isomerase, and D-ribose-5-phosphate isomerase (RpiA) lid domain, among others. Therefore, we consider that this analysis provides clues about the functional constraints associated with the repertoire of enzymatic functions per organism.

Molecular characterization of laccase genes from the basidiomycete Trametes hirsuta Bm-2 and analysis of the 5' untranslated region (5'UTR).

The aim of this study was to identify and characterize laccase genes produced by Trametes hirsuta Bm-2 in a liquid medium, both with and without induction. The amplification of 5'and 3'regions of laccase sequences was obtained by the RACE-PCR method, and these were assembled to obtain a cDNA of total length. Two new laccase genes were isolated from basal medium (lac-B) and lignocellulosic grapefruit substrate (lac-T), both encoding open reading frames of 2566 bp. Both laccase-predicted proteins consisted of 521 amino acids, four copper-binding regions, a signal peptide, and five potential glycosilation sites (Asn-Xaa-Ser/Tre). Moreover, the deduced amino acid sequences share about 76-85% identity with other laccases of WRF. Sequence comparison showed 47 synonymous point mutations between lac-B and lac-T. In addition, 5' untranslated regions (UTR) of laccase genes lac-B and lac-T showed differences in length and number of regulatory elements that may affect transcriptional or translational expression of these genes.

Ploidy Determination in the Pathogenic Fungus Sporothrix spp.

The pathogenic clade of the Sporothrix genus comprises the etiological agents of sporotrichosis, a worldwide emergent disease. Despite the growing understanding of their successful pathogen traits, there is little information on genome sizes and ploidy within the genus. Therefore, in this work, we evaluated the ploidy of four species of the Sporothrix genus, specifically Sporothrix brasiliensis, Sporothrix schenckii, Sporothrix globosa, and Sporothrix pallida. Through cell cycle analysis of the yeast-phase cells, we showed that the DNA content of G0/G1 cells was similar to the genome size determined by whole genome sequencing. Moreover, ploidy of S. schenckii, S. brasiliensis, and S. pallida that was determined by allele composition using next-generation sequencing (NGS) data is consistent with monomorphic positions at each allele. These data show that the analyzed strains of Sporothrix are haploid, or at least aneuploid, thereby laying the foundation for the development of a molecular toolbox for Sporothrix spp.

Temporal distribution and genetic variants in influenza A(H1N1)pdm09 virus circulating in Mexico, seasons 2012 and 2013.

The 2012 and 2013 annual influenza epidemics in Mexico were characterized by presenting different seasonal patterns. In 2012 the A(H1N1)pdm09 virus caused a high incidence of influenza infections after a two-year period of low circulation; whereas the 2013 epidemic presented circulation of the A(H1N1)pdm09 virus throughout the year. We have characterized the molecular composition of the Hemagglutinin (HA) and Neuraminidase (NA) genes of the A(H1N1)pdm09 virus from both epidemic seasons, emphasizing the genetic characteristics of viruses isolated from Yucatan in Southern Mexico. The molecular analysis of viruses from the 2012 revealed that all viruses from Mexico were predominantly grouped in clade 7. Strikingly, the molecular characterization of viruses from 2013 revealed that viruses circulating in Yucatan were genetically different to viruses from other regions of Mexico. In fact, we identified the occurrence of two genetic variants containing relevant mutations at both the HA and NA surface antigens. There was a difference on the temporal circulation of each genetic variant, viruses containing the mutations HA-A141T / NA-N341S were detected in May, June and July; whereas viruses containing the mutations HA-S162I / NA-L206S circulated in August and September. We discuss the significance of these novel genetic changes.

The high content of ß-carotene present in orange-pulp fruits of Carica papaya L. is not correlated with a high expression of the CpLCY-ß2 gene.

We investigated the transcriptional regulation of six genes involved in carotenoid biosynthesis, together with the carotenoid accumulation during postharvest ripening of three different papaya genotypes of contrasting pulp color. Red-pulp genotype (RPG) showed the lowest content of yellow pigments (YP), such as ß-cryptoxanthin, zeaxanthin, and violaxanthin, together with the lowest relative expression levels (REL) of CpLCY-ß2 and CpCHX-ß genes. On the contrary, the yellow-pulp genotype (YPG) showed the highest content of YP and the highest REL of CpLCY-ß2 and CpCHX-ß genes. Interestingly, the orange-pulp genotype (OPG) showed intermediate content of YP and intermediate REL of CpLCY-ß2 and CpCHX-ß genes. The highest content of ß-carotene shown by OPG despite having an intermediate REL of the CpLCY-ß2 genes, suggests a post-transcriptional regulation. Thus, the transcriptional level of the genes, directing the carotenoid biosynthesis pathway, can partially explain the accumulation of carotenoids during the postharvest ripening in C. papaya genotypes of contrasting pulp color.

Genetic Variability in E6 and E7 Oncogenes from Human Papillomavirus Type 58 in Mexican Women.

BACKGROUND/AIMS: The study aimed to describe human papillomavirus (HPV) 58 genetic variability in E6 and E7 oncogenes from women in southeast Mexico and their phylogenetic relationships with the sequences from other geographical regions. METHODS: The E6-E7 region was amplified by nested PCR, and sequenced for identification of polymorphisms, phylogenetic trees construction, and haplotype and fixation tests. RESULTS: HPV58 positive samples were obtained from a repository, 54 were amplified, 47 sequences for the E6 gene, and 51 sequences for the E7 gene were obtained. Fifteen new E6 mutations were found; the most frequent were G279T (G57V; 29.78%), T249G (F47C; 34.04%), and A270G (Y54C; 34.04%), and previously reported c307t (63.82%). For E7, 17 known mutations were found, the most frequent were C632T (T20I), 35.30%, G760A (G63S), 35.30%, and t744g 74.50%. No significant association with the severity of the lesions was found. The polytomy in the E6 tree did not allow proposing phylogenetic relationship, and E7 tree presented defined branches. All sequences were presumably A lineage, most closely related to A1 and/or A3 sublineage. HPV58 variants are not specific for a geographical area. Population and fixation analyses suggest a possible Asian origin of HPV58 from Yucatan. The most frequent E7 haplotype in Yucatan groups with other populations of the world. CONCLUSION: The genetic variability of HPV58 from Mexico was described for the first time. E7 was more conserved than E6. New mutants present exclusively in Yucatan were identified.

Characterization of N-Acetylglucosamine Biosynthesis in Pneumocystis species. A New Potential Target for Therapy.

N-acetylglucosamine (GlcNAc) serves as an essential structural sugar on the cell surface of organisms. For example, GlcNAc is a major component of bacterial peptidoglycan, it is an important building block of fungal cell walls, including a major constituent of chitin and mannoproteins, and it is also required for extracellular matrix generation by animal cells. Herein, we provide evidence for a uridine diphospho (UDP)-GlcNAc pathway in Pneumocystis species. Using an in silico search of the Pneumocystis jirovecii and P. murina (Pm) genomic databases, we determined the presence of at least four proteins implicated in the Saccharomyces cerevisiae UDP-GlcNAc biosynthetic pathway. These genes, termed GFA1, GNA1, AGM1, and UDP-GlcNAc pyrophosphorylase (UAP1), were either confirmed to be present in the Pneumocystis genomes by PCR, or, in the case of Pm uap1 (Pmuap1), functionally confirmed by direct enzymatic activity assay. Expression analysis using quantitative PCR of Pneumocystis pneumonia in mice demonstrated abundant expression of the Pm uap1 transcript. A GlcNAc-binding recombinant protein and a novel GlcNAc-binding immune detection method both verified the presence of GlcNAc in P. carinii (Pc) lysates. Studies of Pc cell wall fractions using high-performance gas chromatography/mass spectrometry documented the presence of GlcNAc glycosyl residues. Pc was shown to synthesize GlcNAc in vitro. The competitive UDP-GlcNAc substrate synthetic inhibitor, nikkomycin Z, suppressed incorporation of GlcNAc by Pc preparations. Finally, treatment of rats with Pneumocystis pneumonia using nikkomycin Z significantly reduced organism burdens. Taken together, these data support an important role for GlcNAc generation in the cell surface of Pneumocystis organisms.

A computationally simplistic poly-phasic approach to explore microbial communities from the Yucatan aquifer as a potential sources of novel natural products.

The need for new antibiotics has sparked a search for the microbes that might potentially produce them. Current sequencing technologies allow us to explore the biotechnological potential of microbial communities in diverse environments without the need for cultivation, benefitting natural product discovery in diverse ways. A relatively recent method to search for the possible production of novel compounds includes studying the diverse genes belonging to polyketide synthase pathways (PKS), as these complex enzymes are an important source of novel therapeutics. In order to explore the biotechnological potential of the microbial community from the largest underground aquifer in the world located in the Yucatan, we used a polyphasic approach in which a simple, non-computationally intensive method was coupled with direct amplification of environmental DNA to assess the diversity and novelty of PKS type I ketosynthase (KS) domains. Our results suggest that the bioinformatic method proposed can indeed be used to assess the novelty of KS enzymes; nevertheless, this in silico study did not identify some of the KS diversity due to primer bias and stringency criteria outlined by the metagenomics pipeline. Therefore, additionally implementing a method involving the direct cloning of KS domains enhanced our results. Compared to other freshwater environments, the aquifer was characterized by considerably less diversity in relation to known ketosynthase domains; however, the metagenome included a family of KS type I domains phylogenetically related, but not identical, to those found in the curamycin pathway, as well as an outstanding number of thiolases. Over all, this first look into the microbial community found in this large Yucatan aquifer and other fresh water free living microbial communities highlights the potential of these previously overlooked environments as a source of novel natural products.