Species delimitation of the liquorice tribe (Leguminosae: Glycyrrhizeae) based on phylogenomic and machine learning analyses

The liquorice tribe Glycyrrhizeae is a leguminous herbaceous group of plants comprised of the genera Glycyrrhiza and Glycyrrhizopsis. Some Glycyrrhiza taxa contain glycyrrhizin, a pharmacologically significant sweet substance that also has applications in crafting industrial materials. Here, we utilized an expanded taxon sampling of Glycyrrhizeae to reconstruct the phylogenetic relationships in the tribe based on genome skimming data, including whole chloroplast genomes, nuclear ribosomal DNA, and low-copy nuclear DNA. We also launched machine learning analysis (MLA) for one species pair with controversial taxonomic boundary. The integrated results indicated Glycyrrhizopsis should be split from Glycyrrhiza, while the former genus Meristotropis should be treated as part of Glycyrrhiza. Glycyrrhizopsis includes two species, Glycyrrhizopsis asymmetrica and Glycyrrhizopsis flavescens, and we recognize 13 species in Glycyrrhiza: Glycyrrhiza acanthocarpa, Glycyrrhiza astragalina, Glycyrrhiza bucharica, Glycyrrhiza echinata, Glycyrrhiza foetida, Glycyrrhiza glabra, Glycyrrhiza gontscharovii, Glycyrrhiza lepidota, Glycyrrhiza macedonica, Glycyrrhiza pallidiflora, Glycyrrhiza squamulosa, Glycyrrhiza triphylla, and Glycyrrhiza yunnanensis. We propose a broader G. glabra that includes former Glycyrrhiza aspera, G. glabra s.s., Glycyrrhiza inflata, and Glycyrrhiza uralensis, and represents the glycyrrhizin-contained medicinal group. Our ancestral state inferences show the ancestor of Glycyrrhiza lacked glycyrrhizin, and the presence of glycyrrhizin evolved twice within Glycyrrhiza during the last one million years. Our integrative phylogenomics-MLA study not only provides new insights into long-standing taxonomic controversies of Glycyrrhizeae, but also represents a useful approach for future taxonomic studies on other plant taxa.     

Nucleotide sequence analysis of the lemur beta-globin gene family: evidence for major rate fluctuations in globin polypeptide evolution

Lemur beta-related globin genes have been isolated and sequenced. Orthology of prosimian and human epsilon-, gamma-, and beta-related globin genes was established by dot-matrix analysis. All of these lemur globin genes potentially encode functional beta-related globin polypeptides, though precisely when the gamma-globin gene is expressed remains unknown. The organization of the 18-kb brown lemur beta-globin gene cluster (5' epsilon-gamma-[psi eta-delta]-beta 3') is consistent with its evolution by contraction via unequal crossing-over from the putative ancestral mammalian beta-globin gene cluster (5' epsilon-gamma- eta-delta-beta 3'). The dwarf lemur nonadult globin genes are arranged as in the brown lemur. Similar levels of synonymous (silent) nucleotide substitutions and noncoding DNA sequence differences have accumulated between species in all of these genes, suggesting a uniform rate of noncoding DNA divergence throughout primate beta-globin gene clusters. These differences are comparable with those observed in the nonfunctional psi eta pseudogene and have therefore accumulated at the presumably maximal neutral rate. In contrast, nonsynonymous (replacement) nucleotide substitutions show a significant heterogeneity in distribution for both the same gene in different lineages and different genes in the same lineage. These major fluctuations in replacement but not silent substitution rates cannot be attributed to changes in mutation rate, suggesting that changes in the rate of globin polypeptide evolution in primates is not governed solely by variable mutation rates.      

The utility of the morphological variation of pollen for resolving the evolutionary history of Billia (subfam. Hippocastanoideae,Sapindaceae)

In this study, we examined the utility of pollen morphology for resolving questions about the evolutionary history of Billia, which is a poorly known genus of Neotropical trees. Billia has been traditionally circumscribed with two species and treated as sister to Aesculus L. However, the number of species in Billia is uncertain, because the genus exhibits abundant morphological diversity but little discontinuous variation. Therefore, Billia may be monotypic and highly polymorphic, or it may have two species with blurred boundaries due to incipient speciation and/or hybridization. Moreover, one recent molecular phylogenetic study shows Billia nested withinAesculus. Our work sought to address the following questions: (i) Are there discontinuities in the pollen of Billia that may suggest species boundaries? (ii) Does the pollen of Billia show evidence for inter-specific hybridization? (iii) Do the exine morphology and size of pollen in Billia differ from those in Aesculus? Our results from scanning electron microscopy showed that pollen exine morphology is not taxonomically informative in Billia but that there are significant differences in pollen size between red- and white-flowered individuals. Thus, our pollen data support the utility of flower color in Billia for species delimitation. Our assessments of pollen viability do not support hybridization in the genus, but cannot be used to rule it out. Finally, pollen exine morphology may lend some support to an evolutionary origin ofBillia within eastern North American Aesculus. In contrast, data on pollen size suggest that Billia may belong in a topological position outside of Aesculus.     

A network-based approach to classify the three domains of life

Background

Chromosomal orthologs can reveal the shared ancestral gene set and their evolutionary trends. Additionally, physico-chemical properties of encoded proteins could provide information about functional adaptation and ecological niche requirements.

Results

We analyzed 7080 genes (five groups of 1416 orthologs each) from Rhizobiales species (S. meliloti, R. etli, and M. loti, plant symbionts; A. tumefaciens, a plant pathogen; and B. melitensis, an animal pathogen). We evaluated their phylogenetic relationships and observed three main topologies. The first, with closer association of R. etli to A. tumefaciens; the second with R. etli closer to S. meliloti; and the third with A. tumefaciens and S. meliloti as the closest pair. This was not unusual, given the close relatedness of these three species. We calculated the synonymous (dS) and nonsynonymous (dN) substitution rates of these orthologs, and found that informational and metabolic functions showed relatively low dN rates; in contrast, genes from hypothetical functions and cellular processes showed high dN rates. An alternative measure of sequence variability, percentage of changes by species, was used to evaluate the most specific proportion of amino acid residues from alignments. When dN was compared with that measure a high correlation was obtained, revealing that much of evolutive information was extracted with the percentage of changes by species at the amino acid level. By analyzing the sequence variability of orthologs with a set of five properties (polarity, electrostatic charge, formation of secondary structures, molecular volume, and amino acid composition), we found that physico-chemical characteristics of proteins correlated with specific functional roles, and association of species did not follow their typical phylogeny, probably reflecting more adaptation to their life styles and niche preferences. In addition, orthologs with low dN rates had residues with more positive values of polarity, volume and electrostatic charge.

Conclusions

These findings revealed that even when orthologs perform the same function in each genomic background, their sequences reveal important evolutionary tendencies and differences related to adaptation. This article was reviewed by: Dr. Purificación López-García, Prof. Jeffrey Townsend (nominated by Dr. J. Peter Gogarten), and Ms. Olga Kamneva.     

Stoichiometry and compartmentation of NADH metabolism in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, reduction of NAD(+) to NADH occurs in dissimilatory as well as in assimilatory reactions. This review discusses mechanisms for reoxidation of NADH in this yeast, with special emphasis on the metabolic compartmentation that occurs as a consequence of the impermeability of the mitochondrial inner membrane for NADH and NAD(+). At least five mechanisms of NADH reoxidation exist in S. cerevisiae. These are: (1) alcoholic fermentation; (2) glycerol production; (3) respiration of cytosolic NADH via external mitochondrial NADH dehydrogenases; (4) respiration of cytosolic NADH via the glycerol-3-phosphate shuttle; and (5) oxidation of intramitochondrial NADH via a mitochondrial 'internal' NADH dehydrogenase. Furthermore, in vivo evidence indicates that NADH redox equivalents can be shuttled across the mitochondrial inner membrane by an ethanol-acetaldehyde shuttle. Several other redox-shuttle mechanisms might occur in S. cerevisiae, including a malate-oxaloacetate shuttle, a malate-aspartate shuttle and a malate-pyruvate shuttle. Although key enzymes and transporters for these shuttles are present, there is as yet no consistent evidence for their in vivo activity. Activity of several other shuttles, including the malate-citrate and fatty acid shuttles, can be ruled out based on the absence of key enzymes or transporters. Quantitative physiological analysis of defined mutants has been important in identifying several parallel pathways for reoxidation of cytosolic and intramitochondrial NADH. The major challenge that lies ahead is to elucidate the physiological function of parallel pathways for NADH oxidation in wild-type cells, both under steady-state and transient-state conditions. This requires the development of techniques for accurate measurement of intracellular metabolite concentrations in separate metabolic compartments.