Evolution of alpha-amylases: architectural features and key residues in the stabilization of the (beta/alpha)(8) scaffold

We provide a comprehensive analysis of the current enzymes with alpha-amylase activity (AAMYs) that belong to family 13 glycoside hydrolase (GH-13; 144 Archaea, Bacteria, and Eukaryota sequences from 87 different species). This study aims to further knowledge of the evolutionary molecular relationships among the sequences of their A and B domains with special emphasis on the correlation between what is observed in the structures and protein evolution. Multialignments for the A domain distinguish two clusters for sequences from Archaea organisms, eight for sequences from Bacteria organisms, and three for sequences from Eukaryota organisms. The clusters for Bacteria do not follow any strict taxonomic pathway; in fact, they are rather scattered. When we compared the A domains of sequences belonging to different kingdoms, we found that various pairs of clusters were significantly similar. Using either sequence similarity with crystallized structures or secondary-structure prediction methods, we identified in all AAMYs the eight putative beta-strands that constitute the beta-sheet in the TIM barrel of the A domain and studied the packing in its interior. We also discovered a "hidden homology" in the TIM barrel, an invariant Gly located upstream in the sequence before the conserved Asp in beta-strand 3. This Gly precedes an alpha-helix and is actively involved in capping its N-terminal end with a capping box. In all cases, a Schellman motif caps the C-terminal end of this helix.     

Distribution of sequences homologous to the impCAB operon of TP110 among bacterial plasmids of different incompatibility groups

Summary Mutagenic DNA repair is a function of many naturally occurring plasmids belonging to several different incompatibility groups. A DNA probe corresponding to the impCAB operon of the IncIl plasmid TP110, which encodes such functions, was used to investigate the distribution of homologous sequences in both related and unrelated plasmids. Southern blotting was used to demonstrate considerable sequence conservation amongst a number of plasmid types, with imp-related sequences being found on plasmids belonging to the I1, I1/B, B and FIV incompatibility groups. However, no homology was detected amongst plasmids of the N and L/M incompatibility groups, many of which carry functionally similar gene clusters. It appears that sequences determining mutagenic repair functions have been largely conserved within any one incompatibility group, but that significant divergent evolution has occurred between groups.     

Unravelling the genetic diversity of ruminal bacteria belonging to the CFB phylum

Molecular biology approaches were employed to examine the genetic diversity of bacteria from the Cytophaga/Flexibacter/Bacteroides (CFB) phylum in the rumen of cattle. By this means we were able to identify cultured strains that represent some of the larger CFB clusters previously identified only by PCR amplification and sequencing. Complete 16S rDNA sequences were obtained for 16 previously isolated rumen strains, including the type strains of Prevotella ruminicola, P. bryantii, P. brevis and P. albensis to represent a wide range of diversity. Phylogenetic analysis of cultured strains revealed the existence of three clusters of ruminal CFB: (i) a cluster of Prevotella strains, which have been found only in the rumen, including the two type strains, P. brevis GA33(T) and P. ruminicola 23(T); (ii) Prevotella spp. that cluster with prevotellas from other ecological niches such as the oral cavity and which include the type strains, P. bryantii B(1)4(T) and P. albensis M384(T); (iii) two Bacteroides spp. strains clustering with B. forsythus of oral origin. In order to establish whether the cultivated isolates cover the whole range of ruminal CFB genetic diversity, 16S rRNA gene sequences were amplified and cloned from DNA extracted from the same rumen samples (one cow in Slovenia, one in Scotland and three in Japan). Sequencing and phylogenetic analysis of 16S rRNA genes confirmed the existence of two superclusters of ruminal Prevotella, one exclusively ruminal and the other including non-ruminal species. In the case of ruminal Bacteroides spp., however, phylogenetic analysis revealed the existence of three new superclusters, one of which has as yet no cultivable counterpart. Interestingly, these Bacteroides clusters were represented almost exclusively by clone libraries from the Japanese cattle and only three sequences were from the European cattle. This study agrees with previous analyses in showing that rumen Prevotella/Bacteroides strains exhibit a remarkable degree of genetic diversity and suggests that different strain groupings may differ greatly in their recovery by cultural methods. The most important conclusion, however, is that cultured strains can be identified that represent some of the larger clusters previously identified only by PCR amplification and sequencing.     

Identification of a Conserved Non-Protein-Coding Genomic Element that Plays an Essential Role in Alphabaculovirus Pathogenesis

Highly homologous sequences 154–157 bp in length grouped under the name of “conserved non-protein-coding element” (CNE) were revealed in all of the sequenced genomes of baculoviruses belonging to the genus Alphabaculovirus. A CNE alignment led to the detection of a set of highly conserved nucleotide clusters that occupy strictly conserved positions in the CNE sequence. The significant length of the CNE and conservation of both its length and cluster architecture were identified as a combination of characteristics that make this CNE different from known viral non-coding functional sequences. The essential role of the CNE in the Alphabaculovirus life cycle was demonstrated through the use of a CNE-knockout Autographa californica multiple nucleopolyhedrovirus (AcMNPV) bacmid. It was shown that the essential function of the CNE was not mediated by the presumed expression activities of the protein- and non-protein-coding genes that overlap the AcMNPV CNE. On the basis of the presented data, the AcMNPV CNE was categorized as a complex-structured, polyfunctional genomic element involved in an essential DNA transaction that is associated with an undefined function of the baculovirus genome.     

节肢动物α-淀粉酶系统进化树构建及分析

目的:构建节肢动物α-淀粉酶的系统进化树,探讨其进化关系,找出进化树中聚类在一起的α-淀粉酶的特异性序列。方法:在美国国立生物技术信息中心(National Center for Biotechnology Information,NCBI)数据库中选取了56个节肢动物的α-淀粉酶氨基酸序列,利用CLUSTALX2.0进行序列比对、MEGA6.0建立进化树,通过BOXSHADE找到聚类的α-淀粉酶特异性序列。结果:56个α-淀粉酶聚类成A、B、C、D四大簇,A簇特异性序列为"VD NHD NQ",B簇特异性序列为"ID NHD NX",C簇特异性序列为"ID NHD NQ",D簇特异性序列为"XGN NHD X"。A、B、C、D四簇都含有保守的NHD(天冬酰胺-组氨酸-天冬氨酸)序列,但序列两端氨基酸种类不同。结论:56个节肢动物α-淀粉酶分为4簇,每簇都有其特异性序列,但都含有保守序列NHD。     

Molecular and cytologic analysis of clusters of moderate repeats of bird genomes

Pigeon genome long sequences containing clusters of moderately repeating elements have been cloned. Molecular analysis has shown a dispersed distribution of the repeats in both pigeon and chicken genomes. Within a single cluster, a scrambled distribution of elements belonging to different families of repeats has been shown. Similar repeated sequences have been revealed within clusters. The analysed clusters of repeats are characterized by a limited structural variability in the genomes. In situ hybridization revealed the localization of sequences complementary to the cloned clusters in pigeon and chicken macrochromosomes. Preferential localization has been demonstrated in telomeric and centromeric chromosome regions as well as in the region of R-bands.     

Features of Active Electrical States of Asymmetrical Dendrites of Brainstem Motoneurons during Generation of Stochastic Implulse Patterns: a Simulation Study

On models of motoneurons of the n. abducens nucleus with reconstructed dendritic arborizations having an active membrane, we investigated features of the relationships between passive transfer properties and dynamics of excitation states of asymmetrical dendrites during generation of complex periodical and stochastic impulse patterns (output neuronal codes). Various patterns were obtained by varying the intensity of tonic synaptic excitation homogeneously distributed over the dendrites. The electrical states of sites belonging to branches of the same dendrite or different dendrites were compared. For this comparison, branches were selected, which, according to the earlier performed cluster analysis, were assigned to the groups (electrotonic clusters) with a high and a low effectiveness of passive transfer of the somatopetal current. The selection took into account features of the dendritic structure of neurons of the exemined type. These were: (i) the presence of groups of the asymmetrical branches differing from each other according to their belonging to different clusters (high or low transfer effectiveness) in different dendrites, and (ii) the presence of branches belonging to different dendrites characterized by significantly different orientations in three-dimensional space of the brainstem within each electrical cluster. Comparative analysis showed that, in a given dendrite during generation of a complex periodical pattern, the asymmetrical branches belonging to high- or low-efficiency clusters were characterized by being in different states (high or low depolarization) in different phases of generation of repeated sequences of action potentials (APs). This relationship was consistent with those previously detected in neurons of other types and in other specimens of neurons of the above-mentioned type. During generation of such periodical spike patterns, the branches of different dendrites belonging to the same electrotonic cluster were in similar states. Similar relationships between the states of the branches of the same dendrite belonging to different clusters were also observed during generation of complex stochastic (non-periodical) impulse patterns. In the latter case, however, the essential feature was that the branches of different dendrites belonging to the same electrotonic cluster were often in opposite states. Thus, the number of combinations of discrete electrical states of asymmetrical parts of the dendritic arborization was much greater. Probably, it is precisely this circumstance that determined the quasi-stochastic nature of the output impulse pattern.     

Phylogenetic Relationships Among JC Virus Strains in Japanese/Koreans and Native Americans Speaking Amerind or Na-Dene

Abstract Many genetic studies using human mtDNA or the Y chromosome have been conducted to elucidate the relationships among the three Native American groups speaking Amerind, Na-Dene, and Eskimo-Aleut. Human polyomavirus JC (JCV) may also help to gain insights into this issue. JCV isolates are classified into more than 10 geographically distinct genotypes (designated subtypes here), which were generated by splits in the three superclusters, Types A, B, and C. A particular subtype of JCV (named MY) belonging to Type B is spread in both Japanese/Koreans and Native Americans speaking Amerind or Na-Dene. In this study, we evaluated the phylogenetic relationships among MY isolates worldwide, using the whole-genome approach, with which a highly reliable phylogeny of JCV isolates can be reconstructed. Thirty-six complete sequences belonging to MY (10 from Japanese/Koreans, 24 from Native Americans, and 2 from others), together with 54 belonging to other subtypes around the world, were aligned and subjected to phylogenetic analysis using the neighbor-joining and maximum-likelihood methods. In the resultant phylogenetic trees, the MY sequences diverged into two Japanese/Korean and five Native American clades with high bootstrap probabilities. Two of the Native American clades contained isolates mainly from Na-Denes and the others contained isolates mainly from Amerinds. The Na-Dene clades were not clustered together, nor were the Amerind clades. In contrast, the two Japanese/Korean clades were clustered at a high bootstrap probability. We concluded that there is no distinction between Amerinds and Na-Denes in terms of indigenous JCVs, although they are linguistically distinguished from each other.     

Phylogenetic relationships among JC virus strains in Japanese/Koreans and Native Americans speaking Amerind or Na-Dene

Many genetic studies using human mtDNA or the Y chromosome have been conducted to elucidate the relationships among the three Native American groups speaking Amerind, Na-Dene, and Eskimo-Aleut. Human polyomavirus JC (JCV) may also help to gain insights into this issue. JCV isolates are classified into more than 10 geographically distinct genotypes (designated subtypes here), which were generated by splits in the three superclusters, Types A, B, and C. A particular subtype of JCV (named MY) belonging to Type B is spread in both Japanese/Koreans and Native Americans speaking Amerind or Na-Dene. In this study, we evaluated the phylogenetic relationships among MY isolates worldwide, using the whole-genome approach, with which a highly reliable phylogeny of JCV isolates can be reconstructed. Thirty-six complete sequences belonging to MY (10 from Japanese/Koreans, 24 from Native Americans, and 2 from others), together with 54 belonging to other subtypes around the world, were aligned and subjected to phylogenetic analysis using the neighbor-joining and maximum-likelihood methods. In the resultant phylogenetic trees, the MY sequences diverged into two Japanese/Korean and five Native American clades with high bootstrap probabilities. Two of the Native American clades contained isolates mainly from Na-Denes and the others contained isolates mainly from Amerinds. The Na-Dene clades were not clustered together, nor were the Amerind clades. In contrast, the two Japanese/Korean clades were clustered at a high bootstrap probability. We concluded that there is no distinction between Amerinds and Na-Denes in terms of indigenous JCVs, although they are linguistically distinguished from each other.     

Metagenomic analyses reveal phylogenetic diversity of carboxypeptidase gene sequences in activated sludge of a wastewater treatment plant in Shanghai,China

Activated sludge of wastewater treatment plants carries a diverse microflora. However, up to 80–90 % of microorganisms in activated sludge cannot be cultured by current laboratory techniques, leaving an enzyme reservoir largely unexplored. In this study, we investigated carboxypeptidase diversity in activated sludge of a wastewater treatment plant in Shanghai, China, by a culture-independent metagenomic approach. Three sets of consensus degenerate hybrid oligonucleotide primers (CODEHOPs) targeting conserved domains of public carboxypeptidases have been designed to amplify carboxypeptidase gene sequences in the metagenomic DNA of activated sludge by PCR. The desired amplicons were evaluated by carboxypeptidase sequence clone libraries and phylogenetic analyses. We uncovered a significant diversity of carboxypeptidases present in the activated sludge. Deduced carboxypeptidase amino acid sequences (127–208 amino acids) were classified into three distinct clusters, α, β, and γ. Sequences belonging to clusters α and β shared 58–97 % identity to known carboxypeptidase sequences from diverse species, whereas sequences in the cluster γ were remarkably less related to public carboxypeptidase homologous in the GenBank database, strongly suggesting that novel carboxypeptidase families or microbial niches exist in the activated sludge. We also observed numerous carboxypeptidase sequences that were much closer to those from representative strains present in industrial and sewage treatment and bioremediation. Thermostable and halotolerant carboxypeptidase sequences were also detected in clusters α and β. Coexistence of various carboxypeptidases is evidence of a diverse microflora in the activated sludge, a feature suggesting a valuable gene resource to be further explored for biotechnology application.     

Discovery and profiling of novel and conserved microRNAs during flower development in Carya cathayensis via deep sequencing

Hickory (Carya cathayensis Sarg.) is an economically important woody plant in China, but its long juvenile phase delays yield. MicroRNAs (miRNAs) are critical regulators of genes and important for normal plant development and physiology, including flower development. We used Solexa technology to sequence two small RNA libraries from two floral differentiation stages in hickory to identify miRNAs related to flower development. We identified 39 conserved miRNA sequences from 114 loci belonging to 23 families as well as two novel and ten potential novel miRNAs belonging to nine families. Moreover, 35 conserved miRNA*s and two novel miRNA*s were detected. Twenty miRNA sequences from 49 loci belonging to 11 families were differentially expressed; all were up-regulated at the later stage of flower development in hickory. Quantitative real-time PCR of 12 conserved miRNA sequences, five novel miRNA families, and two novel miRNA*s validated that all were expressed during hickory flower development, and the expression patterns were similar to those detected with Solexa sequencing. Finally, a total of 146 targets of the novel and conserved miRNAs were predicted. This study identified a diverse set of miRNAs that were closely related to hickory flower development and that could help in plant floral induction.     

Multilocus sequence typing scheme for the characterization of 936-like phages infecting Lactococcus lactis

Lactococcus lactis phage infections are costly for the dairy industry because they can slow down the fermentation process and adversely impact product safety and quality. Although many strategies have been developed to better control phage populations, new virulent phages continue to emerge. Thus, it is beneficial to develop an efficient method for the routine identification of new phages within a dairy plant to rapidly adapt antiphage tactics. Here, we present a multilocus sequence typing (MLST) scheme for the characterization of the 936-like phages, the most prevalent phage group infecting L. lactis strains worldwide. The proposed MLST system targets the internal portion of five highly conserved genomic sequences belonging to the packaging, morphogenesis, and lysis modules. Our MLST scheme was used to analyze 100 phages with different restriction fragment length polymorphism (RFLP) patterns isolated from 11 different countries between 1971 and 2010. PCR products were obtained for all the phages analyzed, and sequence analysis highlighted the high discriminatory power of the MLST system, detecting 93 different sequence types. A conserved locus within the lys gene (coding for endolysin) was the most discriminative, with 65 distinct alleles. The locus within the mcp gene (major capsid protein) was the most conserved (54 distinct alleles). Phylogenetic analyses of the concatenated sequences exhibited a strong concordance of the clusters with the phage host range, indicating the clonal evolution of these phages. A public database has been set up for the proposed MLST system, and it can be accessed at http://pubmlst.org/bacteriophages/.     

The endoxylanases from family 11: computer analysis of protein sequences reveals important structural and phylogenetic relationships

Eighty-two amino acid sequences of the catalytic domains of mature endoxylanases belonging to family 11 have been aligned using the programs MATCHBOX and CLUSTAL. The sequences range in length from 175 to 233 residues. The two glutamates acting as catalytic residues are conserved in all sequences. A very good correlation is found between the presence (at position 100) of an asparagine in the so-called 'alkaline' xylanases, or an aspartic acid in those with a more acidic pH optimum. Four boxes defining segments of highest similarity were detected; they correspond to regions of defined secondary structure: B5, B6, B8 and the carboxyl end of the alpha helix, respectively. Cysteine residues are not common in these sequences (0.7% of all residues), and disulfide bridges are not important in explaining the stability of several thermophilic xylanases. The alignment allows the classification of the enzymes in groups according to sequence similarity. Fungal and bacterial enzymes were found to form mostly separate clusters of higher similarity.     

Human and mouse introns are linked to the same processes and functions through each genome's most frequent non-conserved motifs

We identified the most frequent, variable-length DNA sequence motifs in the human and mouse genomes and sub-selected those with multiple recurrences in the intergenic and intronic regions and at least one additional exonic instance in the corresponding genome. We discovered that these motifs have virtually no overlap with intronic sequences that are conserved between human and mouse, and thus are genome-specific. Moreover, we found that these motifs span a substantial fraction of previously uncharacterized human and mouse intronic space. Surprisingly, we found that these genome-specific motifs are over-represented in the introns of genes belonging to the same biological processes and molecular functions in both the human and mouse genomes even though the underlying sequences are not conserved between the two genomes. In fact, the processes and functions that are linked to these genome-specific sequence-motifs are distinct from the processes and functions which are associated with intronic regions that are conserved between human and mouse. The findings show that intronic regions from different genomes are linked to the same processes and functions in the absence of underlying sequence conservation. We highlight the ramifications of this observation with a concrete example that involves the microsatellite instability gene MLH1.     

Clusters in alpha/beta barrel proteins: implications for protein structure, function, and folding: a graph theoretical approach

The alpha/beta barrel fold is adopted by most enzymes performing a variety of catalytic reactions, but with very low sequence similarity. In order to understand the stabilizing interactions important in maintaining the alpha/beta barrel fold, we have identified residue clusters in a dataset of 36 alpha/beta barrel proteins that have less than 10% sequence identity within themselves. A graph theoretical algorithm is used to identify backbone clusters. This approach uses the global information of the nonbonded interaction in the alpha/beta barrel fold for the clustering procedure. The nonbonded interactions are represented mathematically in the form of an adjacency matrix. On diagonalizing the adjacency matrix, clusters and cluster centers are obtained from the highest eigenvalue and its corresponding vector components. Residue clusters are identified in the strand regions forming the beta barrel and are topologically conserved in all 36 proteins studied. The residues forming the cluster in each of the alpha/beta protein are also conserved among the sequences belonging to the same family. The cluster centers are found to occur in the middle of the strands or in the C-terminal of the strands. In most cases, the residues forming the clusters are part of the active site or are located close to the active site. The folding nucleus of the alpha/beta fold is predicted based on hydrophobicity index evaluation of residues and identification of cluster centers. The predicted nucleation sites are found to occur mostly in the middle of the strands. Proteins 2001;43:103-112.     

Surveying phylogenetic footprints in large gene clusters: applications to Hox cluster duplications

Evolutionarily conserved non-coding genomic sequences represent a potentially rich source for the discovery of gene regulatory regions. Since these elements are subject to stabilizing selection they evolve much more slowly than adjacent non-functional DNA. These so-called phylogenetic footprints can be detected by comparison of the sequences surrounding orthologous genes in different species. Therefore the loss of phylogenetic footprints as well as the acquisition of conserved non-coding sequences in some lineages, but not in others, can provide evidence for the evolutionary modification of cis-regulatory elements. We introduce here a statistical model of footprint evolution that allows us to estimate the loss of sequence conservation that can be attributed to gene loss and other structural reasons. This approach to studying the pattern of cis-regulatory element evolution, however, requires the comparison of relatively long sequences from many species. We have therefore developed an efficient software tool for the identification of corresponding footprints in long sequences from multiple species. We apply this novel method to the published sequences of HoxA clusters of shark, human, and the duplicated zebrafish and Takifugu clusters as well as the published HoxB cluster sequences. We find that there is a massive loss of sequence conservation in the intergenic region of the HoxA clusters, consistent with the finding in [Chiu et al., PNAS 99 (2002) 5492]. The loss of conservation after cluster duplication is more extensive than expected from structural reasons. This suggests that binding site turnover and/or adaptive modification may also contribute to the loss of sequence conservation.     

Analysis and validation of genome-specific DNA variations in 5' flanking conserved sequences of wheat low-molecular-weight glutenin subunit genes

The thirty-three 5' flanking conserved sequences of the known low-molecular-weight subunit (LMW-GS) genes have been divided into eight clusters, which was in agreement with the classification based on the deduced N-terminal protein sequences. The DNA polymorphism between the eight clusters was obtained by sequence alignment, and a total of 34 polymorphic positions were observed in the approximately 200 bp regions, among which 18 polymorphic positions were candidate SNPs. Seven cluster-specific primer sets were designed for seven out of eight clusters containing cluster-specific bases, with which the genomic DNA of the ditelosomic lines of group 1 chromosomes of a wheat variety 'Chinese Spring' was employed to carry out chromosome assignment. The subsequent cloning and DNA sequencing of PCR fragments validated the sequences specificity of the 5' flanking conserved sequences between LMW-GS gene groups in different genomes. These results suggested that the coding and 5' flanking regions of LMW-GS genes are likely to have evolved in a concerted fashion. The seven primer sets developed in this study could be used to isolate the complete ORFs of seven groups of LMW-GS genes, respectively, and therefore possess great value for further research in the contributions of a single LMW-GS gene to wheat quality in the complex genetic background and the efficient selections of quality-related components in breeding programs.     

Nucleotide sequences of Caenorhabditis elegans core histone genes. Genes for different histone classes share common flanking sequence elements

We have determined the nucleotide sequence of core histone genes and flanking regions from two of approximately 11 different genomic histone clusters of the nematode Caenorhabditis elegans. Four histone genes from one cluster (H3, H4, H2B, H2A) and two histone genes from another (H4 and H2A) were analyzed. The predicted amino acid sequences of the two H4 and H2A proteins from the two clusters are identical, whereas the nucleotide sequences of the genes have diverged 9% (H2A) and 12% (H4). Flanking sequences, which are mostly not similar, were compared to identify putative regulatory elements. A conserved sequence of 34 base-pairs is present 19 to 42 nucleotides 3' of the termination codon of all the genes. Within the conserved sequence is a 16-base dyad sequence homologous to the one typically found at the 3' end of histone genes from higher eukaryotes. The C. elegans core histone genes are organized as divergently transcribed pairs of H3-H4 and H2A-H2B and contain 5' conserved sequence elements in the shared spacer regions. One of the sequence elements, 5' CTCCNCCTNCCCACCNCANA 3', is located immediately upstream from the canonical TATA homology of each gene. Another sequence element, 5' CTGCGGGGACACATNT 3', is present in the spacer of each heterotypic pair. These two 5' conserved sequences are not present in the promoter region of histone genes from other organisms, where 5' conserved sequences are usually different for each histone class. They are also not found in non-histone genes of C. elegans. These putative regulatory sequences of C. elegans core histone genes are similar to the regulatory elements of both higher and lower eukaryotes. The coding regions of the genes and the 3' regulatory sequences are similar to those of higher eukaryotes, whereas the presence of common 5' sequence elements upstream from genes of different histone classes is similar to histone promoter elements in yeast.     

The conservation of polyol transporter proteins and their involvement in lichenized Ascomycota

In lichen symbiosis, polyol transfer from green algae is important for acquiring the fungal carbon source. However, the existence of polyol transporter genes and their correlation with lichenization remain unclear. Here, we report candidate polyol transporter genes selected from the genome of the lichen-forming fungus (LFF) Ramalina conduplicans. A phylogenetic analysis using characterized polyol and monosaccharide transporter proteins and hypothetical polyol transporter proteins of R. conduplicans and various ascomycetous fungi suggested that the characterized yeast’ polyol transporters form multiple clades with the polyol transporter-like proteins selected from the diverse ascomycetous taxa. Thus, polyol transporter genes are widely conserved among Ascomycota, regardless of lichen-forming status. In addition, the phylogenetic clusters suggested that LFFs belonging to Lecanoromycetes have duplicated proteins in each cluster. Consequently, the number of sequences similar to characterized yeast’ polyol transporters were evaluated using the genomes of 472 species or strains of Ascomycota. Among these, LFFs belonging to Lecanoromycetes had greater numbers of deduced polyol transporter proteins. Thus, various polyol transporters are conserved in Ascomycota and polyol transporter genes appear to have expanded during the evolution of Lecanoromycetes.