Complete mitochondrial genome of the aluminum-tolerant fungus Rhodotorula taiwanensis RS1 and comparative analysis of Basidiomycota mitochondrial genomes

The transporter PHO84/NtPT1 is a target of aluminum to affect phosphorus absorption in Saccharomyces cerevisiae and Nicotiana tabacum L

The molecular mechanism of aluminum toxicity in biological systems is not completely understood. Saccharomyces cerevisiae is one of the most used model organisms in the study of environmental metal toxicity. Using an unbiased metallomic approach in yeast, we found that aluminum treatment caused phosphorus deprivation, and the lack of phosphorus increased as the pH of the environment decreased compared to the control strain. By screening the phosphate signaling and response pathway (PHO pathway) in yeast with the synthetic lethality of a new phosphorus-restricted aluminum-sensitive gene, we observed that pho84Δ mutation conferred severe growth defect to aluminum under low-phosphorus conditions, and the addition of phosphate alleviated this sensitivity. Subsequently, the data showed that PHO84 determined the intracellular aluminum-induced phosphorus deficiency, and the expression of PHO84 was positively correlated with aluminum stress, which was mediated by phosphorus through the coordinated regulation of PHO4/PHO2. Moreover, aluminum reduced phosphorus absorption and inhibited tobacco plant growth in acidic media. In addition, the high-affinity phosphate transporter NtPT1 in tobacco exhibited similar effects to PHO84, and overexpression of NtPT1 conferred aluminum resistance in yeast cells. Taken together, positive feedback regulation of the PHO pathway centered on the high-affinity phosphate transporters is a highly conservative mechanism in response to aluminum toxicity. The results may provide a basis for aluminum-resistant microorganisms or plant engineering and acidic soil treatment.

Co-evolution of large inverted repeats and G-quadruplex DNA in fungal mitochondria may facilitate mitogenome stability: the case of Malassezia

Mitogenomes are essential due to their contribution to cell respiration. Recently they have also been implicated in fungal pathogenicity mechanisms. Members of the basidiomycetous yeast genus Malassezia are an important fungal component of the human skin microbiome, linked to various skin diseases, bloodstream infections, and they are increasingly implicated in gut diseases and certain cancers. In this study, the comparative analysis of Malassezia mitogenomes contributed to phylogenetic tree construction for all species. The mitogenomes presented significant size and gene order diversity which correlates to their phylogeny. Most importantly, they showed the inclusion of large inverted repeats (LIRs) and G-quadruplex (G4) DNA elements, rendering Malassezia mitogenomes a valuable test case for elucidating the evolutionary mechanisms responsible for this genome diversity. Both LIRs and G4s coexist and convergently evolved to provide genome stability through recombination. This mechanism is common in chloroplasts but, hitherto, rarely found in mitogenomes.

The complete mitochondrial genome of the lipid-producing yeast Rhodotorula toruloides

Mitochondria are semi-autonomous organelles with their own genome and crucial to cellular material and energy metabolism. Here, we report the complete mitochondrial genome of a lipid-producing basidiomycetous yeast Rhodotorula toruloides NP11. The mitochondrial genome of R. toruloides NP11 was assembled into a circular DNA molecule of 125937bp, encoding 15 proteins, 28 transfer RNAs, 2 ribosomal RNA subunits and 10 open reading frames with unknown function. The G + C content (41%) of the mitochondrial genome is substantially lower than that of the nuclear genome (62%) of R. toruloides NP11. Further reanalysis of the transcriptome data confirmed the transcription of four mitochondrial genes. The comparison of the mitochondrial genomes of R. toruloides NP11 and NBRC0880 revealed a significant genetic divergence. These data can complement our understanding of the genetic background of R. toruloides and provide fundamental information for further genetic engineering of this strain.

Aluminum-Nitrogen Interactions in the Soil-Plant System

Aluminum (Al) is the most abundant metal in the Earth's crust and is not an essential element for plant growth. In contrast, nitrogen (N) is the most important mineral element for plant growth, but this non-metal is often present at low levels in soils, and plants are often N deficient. Aluminum toxicity is dominant in acid soils, and so plants growing in acid soils have to overcome both Al toxicity and N limitation. Because of low N-use efficiency, large amounts of N fertilizers are applied to crop fields to achieve high yields, leading to soil acidification and potential Al toxicity. Aluminum lowers plant N uptake and N-use efficiency because Al inhibits root growth. Although numerous studies have investigated the interactions between Al and N, a complete review of these studies was lacking. This review describes: (1) the link between plant Al tolerance and ammonium/nitrate (NH4+/NO3-) preference; (2) the effects of NH4+/NO3- and pH on Al toxicity; (3) the effects of Al on soil N transformations; and (4) the effects of Al on NH4+/NO3- uptake and assimilation by plants. Acid soils are characterized chemically by a relatively high ratio of NH4+ to NO3- and high concentrations of toxic Al. Aluminum-tolerant plants generally prefer NH4+ as an N source, while Al-sensitive plants prefer NO3-. Compared with NO3-, NH4+ increases the solubilization of toxic Al into soil solutions, but NH4+ generally alleviates Al phytotoxicity under solution culture because the protons from NH4+ compete with Al3+ for adsorption sites on the root surface. Plant NO3- uptake and nitrate reductase activity are both inhibited by Al, while plant NH4+ uptake is inhibited to a smaller degree than NO3-. Together, the results of numerous studies indicate that there is a synergistic interaction between plant Al tolerance and NH4+ nutrition. This has important implications for the adaptation of plants to acid soils that are dominated chemically by toxic Al as well as NH4+. Finally, we discuss how this knowledge can be used to increase plant Al tolerance and N-use efficiency in acid soils.

Prospects for Fungal Bioremediation of Acidic Radioactive Waste Sites: Characterization and Genome Sequence of Rhodotorula taiwanensis MD1149

Highly concentrated radionuclide waste produced during the Cold War era is stored at US Department of Energy (DOE) production sites. This radioactive waste was often highly acidic and mixed with heavy metals, and has been leaking into the environment since the 1950s. Because of the danger and expense of cleanup of such radioactive sites by physicochemical processes, in situ bioremediation methods are being developed for cleanup of contaminated ground and groundwater. To date, the most developed microbial treatment proposed for high-level radioactive sites employs the radiation-resistant bacterium Deinococcus radiodurans. However, the use of Deinococcus spp. and other bacteria is limited by their sensitivity to low pH. We report the characterization of 27 diverse environmental yeasts for their resistance to ionizing radiation (chronic and acute), heavy metals, pH minima, temperature maxima and optima, and their ability to form biofilms. Remarkably, many yeasts are extremely resistant to ionizing radiation and heavy metals. They also excrete carboxylic acids and are exceptionally tolerant to low pH. A special focus is placed on Rhodotorula taiwanensis MD1149, which was the most resistant to acid and gamma radiation. MD1149 is capable of growing under 66 Gy/h at pH 2.3 and in the presence of high concentrations of mercury and chromium compounds, and forming biofilms under high-level chronic radiation and low pH. We present the whole genome sequence and annotation of R. taiwanensis strain MD1149, with a comparison to other Rhodotorula species. This survey elevates yeasts to the frontier of biology's most radiation-resistant representatives, presenting a strong rationale for a role of fungi in bioremediation of acidic radioactive waste sites.

Whole genome sequencing of Rhodotorula mucilaginosa isolated from the chewing stick (Distemonanthus benthamianus): insights into Rhodotorula phylogeny, mitogenome dynamics and carotenoid biosynthesis

In industry, the yeast Rhodotorula mucilaginosa is commonly used for the production of carotenoids. The production of carotenoids is important because they are used as natural colorants in food and some carotenoids are precursors of retinol (vitamin A). However, the identification and molecular characterization of the carotenoid pathway/s in species belonging to the genus Rhodotorula is scarce due to the lack of genomic information thus potentially impeding effective metabolic engineering of these yeast strains for improved carotenoid production. In this study, we report the isolation, identification, characterization and the whole nuclear genome and mitogenome sequence of the endophyte R. mucilaginosa RIT389 isolated from Distemonanthus benthamianus, a plant known for its anti-fungal and antibacterial properties and commonly used as chewing sticks. The assembled genome of R. mucilaginosa RIT389 is 19 Mbp in length with an estimated genomic heterozygosity of 9.29%. Whole genome phylogeny supports the species designation of strain RIT389 within the genus in addition to supporting the monophyly of the currently sequenced Rhodotorula species. Further, we report for the first time, the recovery of the complete mitochondrial genome of R. mucilaginosa using the genome skimming approach. The assembled mitogenome is at least 7,000 bases larger than that of Rhodotorula taiwanensis which is largely attributed to the presence of large intronic regions containing open reading frames coding for homing endonuclease from the LAGLIDADG and GIY-YIG families. Furthermore, genomic regions containing the key genes for carotenoid production were identified in R. mucilaginosa RIT389, revealing differences in gene synteny that may play a role in the regulation of the biotechnologically important carotenoid synthesis pathways in yeasts.

Hydroxy-Al and cell-surface negativity are responsible for the enhanced sensitivity of Rhodotorula taiwanensis to aluminum by increased medium pH

Aluminum (Al) is ubiquitous and toxic to microbes. High Al³⁺ concentration and low pH are two key factors responsible for Al toxicity, but our present results contradict this idea. Here, an Al-tolerant yeast strain Rhodotorula taiwanensis RS1 was incubated in glucose media containing Al with a continuous pH gradient from pH 3.1-4.2. The cells became more sensitive to Al and accumulated more Al when pH increased. Calculations using an electrostatic model Speciation Gouy Chapman Stern indicated that, the increased Al sensitivity of cells was associated with AlOH²⁺ and Al(OH) ₂⁺ rather than Al³⁺. The alcian blue (a positively charged dye) adsorption and zeta potential determination of cell surface indicated that, higher pH than 3.1 increased the negative charge and Al adsorption at the cell surface. Taken together, the enhanced sensitivity of R. taiwanensis RS1 to Al from pH 3.1-4.2 was associated with increased hydroxy-Al and cell-surface negativity.

A potential contribution of the less negatively charged cell wall to the high aluminum tolerance of Rhodotorula taiwanensis RS1

Rhodotorula taiwanensis RS1 (Rt) is a high-aluminum (Al)-tolerant yeast that can survive Al at concentrations up to 200 mM. In this study, we compared Rt with an Al-sensitive congeneric strain, R. mucilaginosa AKU 4812 (Rm) and Al sensitive mutant 1 (alsm1) of Rt, to explore the Al tolerance mechanisms of Rt. The growth of Rm was completely inhibited by 1 mM Al, but that of Rt was not inhibited until Al concentration was more than 70 mM. The growth of alsm1 was inhibited much more by 70 mM and 100 mM Al than that of Rt. Compared with Rm cells, Rt cells accumulated less Al in the cell wall and cytoplasm. A time-course analysis showed that Al was absorbed by Rm cells much more rapidly than by Rt cells when exposed to the same Al concentration. Meanwhile, the Al content of alsm1 was higher than that of Rt. Although the cell wall of Rt was thicker than that of alsm1 and Rm under control and 0.1 mM Al, that of Rt was thinner than that of alsm1 under 70 mM Al despite that their cell walls were thickened. The alcian blue adsorption was lower and cell wall zeta-potential was higher in Rt and alsm1 than in Rm, indicating a less negative charge of cell wall of Rt and alsm1 than that of Rm. Taken together, the less negatively charged cell wall of Rt may restrict the adsorption of cationic Al in cells, potentially contributing to its high Al tolerance. Copyright © 2016 John Wiley & Sons, Ltd.

Complete mitochondrial genome of the Verticillium-wilt causing plant pathogen Verticillium nonalfalfae

Verticillium nonalfalfae is a fungal plant pathogen that causes wilt disease by colonizing the vascular tissues of host plants. The disease induced by hop isolates of V. nonalfalfae manifests in two different forms, ranging from mild symptoms to complete plant dieback, caused by mild and lethal pathotypes, respectively. Pathogenicity variations between the causal strains have been attributed to differences in genomic sequences and perhaps also to differences in their mitochondrial genomes. We used data from our recent Illumina NGS-based project of genome sequencing V. nonalfalfae to study the mitochondrial genomes of its different strains. The aim of the research was to prepare a V. nonalfalfae reference mitochondrial genome and to determine its phylogenetic placement in the fungal kingdom. The resulting 26,139 bp circular DNA molecule contains a full complement of the 14 "standard" fungal mitochondrial protein-coding genes of the electron transport chain and ATP synthase subunits, together with a small rRNA subunit, a large rRNA subunit, which contains ribosomal protein S3 encoded within a type IA-intron and 26 tRNAs. Phylogenetic analysis of this mitochondrial genome placed it in the Verticillium spp. lineage in the Glomerellales group, which is also supported by previous phylogenetic studies based on nuclear markers. The clustering with the closely related Verticillium dahliae mitochondrial genome showed a very conserved synteny and a high sequence similarity. Two distinguishing mitochondrial genome features were also found-a potential long non-coding RNA (orf414) contained only in the Verticillium spp. of the fungal kingdom, and a specific fragment length polymorphism observed only in V. dahliae and V. nubilum of all the Verticillium spp., thus showing potential as a species specific biomarker.

Intronic and plasmid-derived regions contribute to the large mitochondrial genome sizes of Agaricomycetes

Sizes of mitochondrial genomes vary extensively between fungal species although they typically contain a conserved set of core genes. We have characterised the mitochondrial genome of the conifer root rot pathogen Heterobasidion irregulare and compared the size, gene content and structure of 20 Basidiomycete mitochondrial genomes. The mitochondrial genome of H. irregulare was 114, 193 bp and contained a core set of 15 protein coding genes, two rRNA genes and 26 tRNA genes. In addition, we found six non-conserved open reading frames (ORFs) and four putative plasmid genes clustered in three separate regions together with 24 introns and 14 intronic homing endonuclease genes, unequally spread across seven of the core genes. The size differences among the 20 Basidiomycetes can largely be explained by length variation of intergenic regions and introns. The Agaricomycetes contained the nine largest mitochondrial genomes in the Basidiomycete group and Agaricomycete genomes are significantly (p < 0.001) larger than the other Basidiomycetes. A feature of the Agaricomycete mitochondrial genomes in this study was the simultaneous occurrence of putative plasmid genes and non-conserved ORFs, with Cantharellus cibarius as only exception, where no non-conserved ORF was identified. This indicates a mitochondrial plasmid origin of the non-conserved ORFs or increased mitochondrial genome dynamics of species harbouring mitochondrial plasmids. We hypothesise that two independent factors are the driving forces for large mitochondrial genomes: the homing endonuclease genes in introns and integration of plasmid DNA.

High Rhodotorula sequences in skin transcriptome of patients with diffuse systemic sclerosis

Previous studies have suggested a role for pathogens as a trigger of systemic sclerosis (SSc), although neither a pathogen nor a mechanism of pathogenesis is known. Here we show enrichment of Rhodotorula sequences in the skin of patients with early, diffuse SSc compared with that in normal controls. RNA-seq was performed on four SSc patients and four controls, to a depth of 200 million reads per patient. Data were analyzed to quantify the nonhuman sequence reads in each sample. We found little difference between bacterial microbiome and viral read counts, but found a significant difference between the read counts for a mycobiome component, R. glutinis. Normal samples contained almost no detected R. glutinis or other Rhodotorula sequence reads (mean score 0.021 for R. glutinis, 0.024 for all Rhodotorula). In contrast, SSc samples had a mean score of 5.039 for R. glutinis (5.232 for Rhodotorula). We were able to assemble the D1-D2 hypervariable region of the 28S ribosomal RNA (rRNA) of R. glutinis from each of the SSc samples. Taken together, these results suggest that R. glutinis may be present in the skin of early SSc patients at higher levels than in normal skin, raising the possibility that it may be triggering the inflammatory response found in SSc.

Mitochondrial genome of the basidiomycetous yeast Jaminaea angkorensis

Jaminaea angkorensis is an anamorphic basidiomycetous yeast species originally isolated from decaying leaves in Cambodia. Taxonomically, J. angkorensis is affiliated with Microstromatales (Exobasidiomycetes, Ustilaginomycotina, Basidiomycota) and represents a basal phylogenetic lineage of this fungal order. To perform a comparative analysis of J. angkorensis with other basidiomycetes, we determined and analyzed its complete mitochondrial DNA sequence. The mitochondrial genome is represented by 29,999 base pairs long, circular DNA containing 32 % guanine and cytosine residues. Its genetic organization is relatively compact and comprises typical genes for 15 conserved proteins involved in oxidative phosphorylation (atp6, 8, and 9; cob; cox1, 2, and 3; and nad1, 2, 3, 4, 4L, 5, and 6) and translation (rps3), two ribosomal RNAs (rnl and rns) and twenty-two transfer RNAs (trnA-Y). Although the gene content is similar to other basidiomycetes, the gene orders in the examined species exhibit only a limited synteny, reflecting their phylogenetic distances and extensive genome rearrangements. In addition, a comparative analysis of basidiomycete mitochondrial genomes indicates that stop-to-tryptophan reassignment of the UGA codon was accompanied by structural alterations of tRNA-Trp(CCA). These results provide an insight into the evolution of the genetic code in fungal mitochondria.