Patterns of tandem repetition in plant whole genome assemblies

Tandem repeats often confound large genome assemblies. A survey of tandemly arrayed repetitive sequences was carried out in whole genome sequences of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the monocots rice and sorghum, and the dicots Arabidopsis thaliana, poplar, grapevine, and papaya, in order to test how these assemblies deal with this fraction of DNA. Our results suggest that plant genome assemblies preferentially include tandem repeats composed of shorter monomeric units (especially dinucleotide and 9–30-bp repeats), while higher repetitive units pose more difficulties to assemble. Nevertheless, notwithstanding that currently available sequencing technologies struggle with higher arrays of repeated DNA, major well-known repetitive elements including centromeric and telomeric repeats as well as high copy-number genes, were found to be reasonably well represented. A database including all tandem repeat sequences characterized here was created to benefit future comparative genomic analyses.     

Conserved gene clusters in the highly rearranged chloroplast genomes of Chlamydomonas moewusii and Chlamydomonas reinhardtii

We have extended to about 75 the number of genes mapped on the Chlamydomonas moewusii and Chlamydomonas reinhardtii chloroplast DNAs (cpDNAs) by partial sequencing of the very closely related C. eugametos and C. moewusii cpDNAs and by hybridizations with Chlamydomonas chloroplast gene-specific sequences. Only four of these genes (tscA and three reading frames) have not been identified in any other algal cpDNAs and thus may be specific to Chlamydomonas. Although the C. moewusii and C. reinhardtii cpDNAs differ by complex sequence rearrangements, 38 genes scattered throughout the genome define 12 conserved clusters of closely linked loci. Aside from the rRNA operon, four of these gene clusters share similarity to evolutionarily primitive operons found in other cpDNAs, representing in fact remnants of these operons. Our results thus indicate that most of the ancestral bacterial operons that characterize the chloroplast genome organization of land plants and early-diverging photosynthetic eukaryotes have been disrupted before the emergence of the polyphyletic genus Chlamydomonas. All gene rearrangements between the C. moewusii and C. reinhardtii cpDNAs, with the exception of those accounting for the relocations of atpA, psbI and rbcL, occurred within corresponding regions of the genome. One of these rearrangements seems to have led to disruption of the ancestral region containing rpl23, rpl2, rps19, rpl16, rpl14, rpl5, rps8 and the psaA exon 1. This gene cluster, which bears striking similarity to the Escherichia coli S10 and spc operons, spans a continuous DNA segment in C. reinhardtii, while it maps to two separate fragments in C. moewusii.     

Regulatory Sequences of Orthologous petD Chloroplast mRNAs are Highly Specific among Chlamydomonas Species

The 5′ untranslated regions (UTR) of chloroplast mRNAs often contain regulatory sequences that control RNA stability and/or translation. The petD chloroplast mRNA in Chlamydomonas reinhardtii has three such essential regulatory elements in its 362-nt long 5′ UTR. To further analyze these elements, we compared 5′ UTR sequences from four Chlamydomonas species (C. reinhardtii, C. incerta, C. moewusii and C. eugametos) and five independent strains of C. reinhardtii. Overall, these petD 5′ UTRs have relatively low sequence conservation across these species. In contrast, sequences of the three regulatory elements and their relative positions appear partially conserved. Functionality of the 5′ UTRs was tested in C. reinhardtii chloroplasts using β-glucuronidase reporter genes, and the nearly identical C. incerta petD functioned for mRNA stability and translation in C. reinhardtii chloroplasts while the more divergent C. eugametos petD did not. This identified what may be key features in these elements. We conclude that these petD regulatory elements, and possibly the corresponding trans-acting factors, function via mechanisms highly specific and surprisingly sensitive to minor sequence changes. This provides a new and broader perspective of these important regulatory sequences that affect photosynthesis in these algae.     

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10⁻¹⁰, and is highly variable between MA lines, ranging from μ = 0.35 × 10⁻¹⁰ to μ = 131.7 × 10⁻¹⁰. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.     

A marine Chlamydomonas sp. emerging as an algal model

The freshwater microalga Chlamydomonas reinhardtii, which lives in wet soil, has served for decades as a model for numerous biological processes, and many tools have been introduced for this organism. Here, we have established a stable nuclear transformation for its marine counterpart, Chlamydomonas sp. SAG25.89, by fusing specific cis‐acting elements from its Actin gene with the gene providing hygromycin resistance and using an elaborated electroporation protocol. Like C. reinhardtii, Chlamydomonas sp. has a high GC content, allowing reporter genes and selection markers to be applicable in both organisms. Chlamydomonas sp. grows purely photoautotrophically and requires ammonia as a nitrogen source because its nuclear genome lacks some of the genes required for nitrogen metabolism. Interestingly, it can grow well under both low and very high salinities (up to 50 g · L^‐1) rendering it as a model for osmotolerance. We further show that Chlamydomonas sp. grows well from 15 to 28°C, but halts its growth at 32°C. The genome of Chlamydomonas sp. contains some gene homologs the expression of which is regulated according to the ambient temperatures and/or confer thermal acclimation in C. reinhardtii. Thus, knowledge of temperature acclimation can now be compared to the marine species. Furthermore, Chlamydomonas sp. can serve as a model for studying marine microbial interactions and for comparing mechanisms in freshwater and marine environments. Chlamydomonas sp. was previously shown to be immobilized rapidly by a cyclic lipopeptide secreted from the antagonistic bacterium Pseudomonas protegens PF‐5, which deflagellates C. reinhardtii.     

Functional Characterization of the Chlamydomonas reinhardtii ERG3 Ortholog,a Gene Involved in the Biosynthesis of Ergosterol

Background

The predominant sterol in the membranes of the alga Chlamydomonas reinhardtii is ergosterol, which is commonly found in the membranes of fungi, but is rarely found in higher plants. Higher plants and fungi synthesize sterols by different pathways, with plants producing cycloartenol as a precursor to end-product sterols, while non-photosynthesizing organisms like yeast and humans produce lanosterol as a precursor. Analysis of the C. reinhardtii genome sequence reveals that this algae is also likely to synthesize sterols using a pathway resembling the higher plant pathway, indicating that its sterols are synthesized somewhat differently than in fungi. The work presented here seeks to establish experimental evidence to support the annotated molecular function of one of the sterol biosynthetic genes in the Chlamydomonas genome.

Methodology/Principal Findings

A gene with homology to the yeast sterol C-5 desaturase, ERG3, is present in the Chlamydomonas genome. To test whether the ERG3 ortholog of C. reinhardtii encodes a sterol C-5 desaturase, Saccharomyces cerevisiae ERG3 knockout strains were created and complemented with a plasmid expressing the Chlamydomonas ERG3. Expression of C. reinhardtii ERG3 cDNA in erg3 null yeast was able to restore ergosterol biosynthesis and reverse phenotypes associated with lack of ERG3 function.

Conclusions/Significance

Complementation of the yeast erg3 null phenotypes strongly suggests that the gene annotated as ERG3 in C. reinhardtii functions as a sterol C-5 desaturase.     

A comparative, BAC end sequence enabled map of the genome of the American mink (Neovison vison)

In this report we present the results of the analysis of approximately 2.7 Mb of genomic information for the American mink (Neovison vison) derived through BAC end sequencing. Our study, which encompasses approximately 1/1000^th of the mink genome, suggests that simple sequence repeats (SSRs) are less common in the mink than in the human genome, whereas the average GC content of the mink genome is slightly higher than that of its human counterpart. The 2.7 Mb mink genomic dataset also contained 2,416 repeat elements (retroids and DNA transposons) occupying almost 31% of the sequence space. Among repeat elements, LINEs were over-represented and endogenous viruses (aka LTRs) under-represented in comparison to the human genome. Finally, we present a virtual map of the mink genome constructed with reference to the human and canine genome assemblies using a comparative genomics approach and incorporating over 200 mink BESs with unique hits to the human genome.     

Identification and analysis of smORFs in Chlamydomonas reinhardtii

Small open reading frames (smORFs) have been acknowledged as an important partner in organism functions ranging from bacteria to higher eukaryotes. However, there is a lack of investigation of smORFs in green algae, despite their importance in ecology and evolution. We applied bioinformatic analysis, ribosome profiling, and small peptide proteomics to provide a genome-wide and high-confident smORF database in the model green alga Chlamydomonas reinhardtii. The whole genome was screened first to mine potential coding smORFs. Then conservative analysis, ribosome profiling, and proteomics data were processed to identify conserved smORFs and generate translation evidence. The combination of procedures resulted in 2014 smORFs that might exist in the C. reinhardtii genome. The expression of smORFs in Cd treatment suggested that two smORFs might participate in redox reaction, three in inorganic phosphate transport, and one in DNA repair under stress. Our study built a genome-widely database in C. reinhardtii, providing target smORFs for further research.     

Variation in assembly stoichiometry in non‐metazoan homologs of the hub domain of Ca2+/calmodulin‐dependent protein kinase II

The multi‐subunit Ca²⁺/calmodulin‐dependent protein kinase II (CaMKII) holoenzyme plays a critical role in animal learning and memory. The kinase domain of CaMKII is connected by a flexible linker to a C‐terminal hub domain that assembles into a 12‐ or 14‐subunit scaffold that displays the kinase domains around it. Studies on CaMKII suggest that the stoichiometry and dynamic assembly/disassembly of hub oligomers may be important for CaMKII regulation. Although CaMKII is a metazoan protein, genes encoding predicted CaMKII‐like hub domains, without associated kinase domains, are found in the genomes of some green plants and bacteria. We show that the hub domains encoded by three related green algae, Chlamydomonas reinhardtii, Volvox carteri f. nagarensis, and Gonium pectoral, assemble into 16‐, 18‐, and 20‐subunit oligomers, as assayed by native protein mass spectrometry. These are the largest known CaMKII hub domain assemblies. A crystal structure of the hub domain from C. reinhardtii reveals an 18‐subunit organization. We identified four intra‐subunit hydrogen bonds in the core of the fold that are present in the Chlamydomonas hub domain, but not in metazoan hubs. When six point mutations designed to recapitulate these hydrogen bonds were introduced into the human CaMKII‐α hub domain, the mutant protein formed assemblies with 14 and 16 subunits, instead of the normal 12‐ and 14‐subunit assemblies. Our results show that the stoichiometric balance of CaMKII hub assemblies can be shifted readily by small changes in sequence.     

Hydrogen evolution as a consumption mode of reducing equivalents in green algal fermentation

Dark anaerobic fermentation in the green algae Chlamydomonas MGA 161, Chlamydomonas reinhardtii, Chlorella pyrenoidosa, and Chlorococcum minutum was studied. Our isolate, Chlamydomonas MGA 161, was unusual in having high H₂ but almost no formate. The fermentation pattern in Chlamydomonas MGA 161 was altered by changes in the NaCl or NH₄Cl concentration. Glycerol formation increased at low (0.1%) and high (7%) NaCl concentrations; starch degradation, and formation of ethanol, H₂, and CO₂ increased with the addition of NH₄Cl to above 5 millimolar in N-deficient cells. C. reinhardtii and C. pyrenoidosa exhibited a very similar anaerobic metabolism, forming formate, acetate and ethanol in a ratio of about 2:2:1. C. minutum was also unusual in forming acetate, glycerol, and CO₂ as its main products, with H₂, formate, and ethanol being formed in negligible amounts. In the presence of CO, ethanol formation increased twofold in Chlamydomonas MGA 161 and C. reinhardtii, but the fermentation pattern in C. minutum did not change. An experiment with hypophosphite addition showed that dark H₂ evolution of the Escherichia coli type could be ruled out in Chlamydomonas MGA 161 and C. reinhardtii. Among the green algae investigated, three fermentation types were identified by the distribution pattern of the end products, which reflected the consumption mode of reducing equivalents in the cells.     

Detection and characterization of phosphatidylcholine in various strains of the genus Chlamydomonas (Volvocales,Chlorophyceae)

The laboratory strains of Chlamydomonas reinhardtii have been reported to contain no phosphatidylcholine (PC), which is considered to be replaced by another zwitterionic lipid, diacylglyceryl-N,N,N-trimethylhomoserine (DGTS). According to the recently published classification, the strains belonged to the subgroup Reinhardtinia. Screening for PC in 13 selected strains of Chlamydomonas in the NIES Algal Collection, which are different in habitats and belong to different phylogenetic subgroups in the genus, revealed the presence of PC in four strains: a strain in the subgroup Polytominia, and three strains in Reinhardtinia. PC was not detected in three other strains of Reinhardtinia analyzed. The presence/absence of PC was not related to the phylogenetic relationship based on 18S rRNA. DGTS was detected in all the strains analyzed. The rare isomer of linolenic acid, 18:3(5,9,12), which has been found in the DGTS of C. reinhardtii, was found in the PC of the two strains and in the DGTS of the five strains. The occurrence of this fatty acid seems limited to a branch of Reinhardtinia. Acquisition and loss of PC in various strains of Chlamydomonas are discussed from the viewpoint of evolution of PC biosynthetic pathway.     

Alternative Acetate Production Pathways in Chlamydomonas reinhardtii during Dark Anoxia and the Dominant Role of Chloroplasts in Fermentative Acetate Production

Chlamydomonas reinhardtii insertion mutants disrupted for genes encoding acetate kinases (EC 2.7.2.1) (ACK1 and ACK2) and a phosphate acetyltransferase (EC 2.3.1.8) (PAT2, but not PAT1) were isolated to characterize fermentative acetate production. ACK1 and PAT2 were localized to chloroplasts, while ACK2 and PAT1 were shown to be in mitochondria. Characterization of the mutants showed that PAT2 and ACK1 activity in chloroplasts plays a dominant role (relative to ACK2 and PAT1 in mitochondria) in producing acetate under dark, anoxic conditions and, surprisingly, also suggested that Chlamydomonas has other pathways that generate acetate in the absence of ACK activity. We identified a number of proteins associated with alternative pathways for acetate production that are encoded on the Chlamydomonas genome. Furthermore, we observed that only modest alterations in the accumulation of fermentative products occurred in the ack1, ack2, and ack1 ack2 mutants, which contrasts with the substantial metabolite alterations described in strains devoid of other key fermentation enzymes.     

Target of rapamycin and LST8 proteins associate with membranes from the endoplasmic reticulum in the unicellular green alga Chlamydomonas reinhardtii

The highly conserved target of rapamycin (TOR) kinase is a central controller of cell growth in all eukaryotes. TOR exists in two functionally and structurally distinct complexes, termed TOR complex 1 (TORC1) and TORC2. LST8 is a TOR-interacting protein that is present in both TORC1 and TORC2. Here we report the identification and characterization of TOR and LST8 in large protein complexes in the model photosynthetic green alga Chlamydomonas reinhardtii. We demonstrate that Chlamydomonas LST8 is part of a rapamycin-sensitive TOR complex in this green alga. Biochemical fractionation and indirect immunofluorescence microscopy studies indicate that TOR and LST8 exist in high-molecular-mass complexes that associate with microsomal membranes and are particularly abundant in the peri-basal body region in Chlamydomonas cells. A Saccharomyces cerevisiae complementation assay demonstrates that Chlamydomonas LST8 is able to functionally and structurally replace endogenous yeast LST8 and allows us to propose that binding of LST8 to TOR is essential for cell growth.     

Mitochondrial and Plastid Genomes of the Colonial Green Alga Gonium pectorale Give Insights into the Origins of Organelle DNA Architecture within the Volvocales

Volvocalean green algae have among the most diverse mitochondrial and plastid DNAs (mtDNAs and ptDNAs) from the eukaryotic domain. However, nearly all of the organelle genome data from this group are restricted to unicellular species, like Chlamydomonas reinhardtii, and presently only one multicellular species, the ∼4,000-celled Volvox carteri, has had its organelle DNAs sequenced. The V. carteri organelle genomes are repeat rich, and the ptDNA is the largest plastome ever sequenced. Here, we present the complete mtDNA and ptDNA of the colonial volvocalean Gonium pectorale, which is comprised of ∼16 cells and occupies a phylogenetic position closer to that of V. carteri than C. reinhardtii within the volvocine line. The mtDNA and ptDNA of G. pectorale are circular-mapping AT-rich molecules with respective lengths and coding densities of 16 and 222.6 kilobases and 73 and 44%. They share some features with the organelle DNAs of V. carteri, including palindromic repeats within the plastid compartment, but show more similarities with those of C. reinhardtii, such as a compact mtDNA architecture and relatively low organelle DNA intron contents. Overall, the G. pectorale organelle genomes raise several interesting questions about the origin of linear mitochondrial chromosomes within the Volvocales and the relationship between multicellularity and organelle genome expansion.     

The plastid aldolase gene fromChlamydomonas reinhardtii: Intron/exon organization,evolution, and promoter structure

Genomic clones encoding the plastidic fructose- 1,6-bisphosphate aldolase ofChlamydomonas reinhardtii were isolated and sequenced. The gene contains three introns which are located within the coding sequence for the mature protein. No introns are located within or near the sequence encoding the transit-peptide, in contrast to the genes for plastidic aldolases of higher plants. Neither the number nor the positions of the three introns of theC. reinhardtii aldolase gene are conserved in the plastidic or cytosolic aldolase genes of higher plants and animals. The 5′ border sequences of introns in the aldolase gene ofC. reinhardtii exhibit the conserved plant consensus sequence. The 3′ acceptor splice sites for introns 1 and 3 show much less similarity to the eukaryotic consensus sequences than do those of intron 2. The plastidic aldolase gene has two tandemly repeated CAAT box motifs in the promoter region. Genomic Southern blots indicate that the gene is encoded by a single locus in theC. reinhardtii genome.     

Rapid Induction of Lipid Droplets in Chlamydomonas reinhardtii and Chlorella vulgaris by Brefeldin A

Algal lipids are the focus of intensive research because they are potential sources of biodiesel. However, most algae produce neutral lipids only under stress conditions. Here, we report that treatment with Brefeldin A (BFA), a chemical inducer of ER stress, rapidly triggers lipid droplet (LD) formation in two different microalgal species, Chlamydomonas reinhardtii and Chlorella vulgaris. LD staining using Nile red revealed that BFA-treated algal cells exhibited many more fluorescent bodies than control cells. Lipid analyses based on thin layer chromatography and gas chromatography revealed that the additional lipids formed upon BFA treatment were mainly triacylglycerols (TAGs). The increase in TAG accumulation was accompanied by a decrease in the betaine lipid diacylglyceryl N,N,N-trimethylhomoserine (DGTS), a major component of the extraplastidic membrane lipids in Chlamydomonas, suggesting that at least some of the TAGs were assembled from the degradation products of membrane lipids. Interestingly, BFA induced TAG accumulation in the Chlamydomonas cells regardless of the presence or absence of an acetate or nitrogen source in the medium. This effect of BFA in Chlamydomonas cells seems to be due to BFA-induced ER stress, as supported by the induction of three homologs of ER stress marker genes by the drug. Together, these results suggest that ER stress rapidly triggers TAG accumulation in two green microalgae, C. reinhardtii and C. vulgaris. A further investigation of the link between ER stress and TAG synthesis may yield an efficient means of producing biofuel from algae.     

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.     

Cardiolipin of Chlamydomonas reinhardtii 137

The fatty acids of cardiolipin from the phototrophic green alga Chlamydomonas reinhardtii 137⁺ have been quantitatively analysed. Comparison is made at the molecular level between the cardiolipin of Chlamydomonas and that of higher plant tissue.