Description of Naviculavanseea sp. nov. (Naviculales, Naviculaceae), a new species of diatom from the highly alkaline Lake Van (Republic of Türkiye) with complete characterisation of its organellar genomes and multigene phylogeny

The current article describes Naviculavanseeasp. nov., a new species of diatom from Lake Van, a highly alkaline lake in Eastern Anatolia (Türkiye). The description is based on light and scanning electron microscopy performed on two monoclonal cultures. The complete nuclear rRNA clusters and plastid genomes have been sequenced for these two strains and the complete mitogenome for one of them. The plastome of both strains shows the probable loss of a functional ycf35 gene. They also exhibit two IB4 group I introns in their rrl, each encoding for a putative LAGLIDADG homing endonuclease, with the first L1917 IB4 intron reported amongst diatoms. The Maximum Likelihood phylogeny inferred from a concatenated alignment of 18S, rbcL and psbC distinguishes N.vanseea sp. nov. from the morphologically similar species Naviculacincta and Naviculamicrodigitoradiata.

The invasive land flatworm Arthurdendyus triangulatus has repeated sequences in the mitogenome, extra-long cox2 gene and paralogous nuclear rRNA clusters

Using a combination of short- and long-reads sequencing, we were able to sequence the complete mitochondrial genome of the invasive 'New Zealand flatworm' Arthurdendyus triangulatus (Geoplanidae, Rhynchodeminae, Caenoplanini) and its two complete paralogous nuclear rRNA gene clusters. The mitogenome has a total length of 20,309 bp and contains repetitions that includes two types of tandem-repeats that could not be solved by short-reads sequencing. We also sequenced for the first time the mitogenomes of four species of Caenoplana (Caenoplanini). A maximum likelihood phylogeny associated A. triangulatus with the other Caenoplanini but Parakontikia ventrolineata and Australopacifica atrata were rejected from the Caenoplanini and associated instead with the Rhynchodemini, with Platydemus manokwari. It was found that the mitogenomes of all species of the subfamily Rhynchodeminae share several unusual structural features, including a very long cox2 gene. This is the first time that the complete paralogous rRNA clusters, which differ in length, sequence and seemingly number of copies, were obtained for a Geoplanidae.

Mitochondrial genome sequence of the protist Ancyromonas sigmoides Kent, 1881 (Ancyromonadida) from the Sugluk Inlet, Hudson Strait, Nunavik, Québec

Introduction

There is little information on evolutionarily ancient eukaryotes, which are often referred to as basal eukaryotes, in Arctic waters. Despite earlier studies being conducted in the Russian White Sea, only few have been reported.

Methods

Following a shotgun sequence survey of diatom cultures from Sugluk Inlet off the Hudson Strait in Northern Québec, we obtained the complete mitochondrial genome and the operon of nuclear ribosomal RNA genes from a strain that matches that of Ancyromonas sigmoides (Kent, 1881).

Results

The sequence of the mitogenome retrieved was 41,889 bp in length and encoded 38 protein-coding genes, 5 non-conserved open-reading frames, and 2 rRNA and 24 tRNA genes. The mitogenome has retained sdh2 and sdh3, two genes of the succinate dehydrogenase complex, which are sometimes found among basal eukaryotes but seemingly missing among the Malawimonadidae, a lineage sister to Ancyromonadida in some phylogenies. The phylogeny inferred from the 18S rRNA gene associated A. sigmoides from Sugluk Inlet with several other strains originating from the Arctic. The study also unveiled the presence of a metagenomic sequence ascribed to bacteria in GenBank, but it was clearly a mitochondrial genome with a gene content highly similar to that of A. sigmoides, including the non-conserved open-reading frames.

Discussion

After re-annotation, a phylogeny was inferred from mitochondrial protein sequences, and it strongly associated A. sigmoides with the misidentified organism, with the two being possibly conspecific or sibling species as they are more similar to one another than to species of the genus Malawimonas. Overall our phylogeny showed that the ice associated ancryomonads were clearly distinct from more southerly strains.

The mitochondrial genome of the bioluminescent fish Malacosteus niger Ayres, 1848 (Stomiidae, Actinopterygii) is large and complex, and contains an inverted-repeat structure

We determined the complete mitogenome sequence of the bioluminescent fish Malacosteus niger using long-read sequencing technologies. The 21,263 bp mitogenome features a complex structure with two copies of a 1198-bp inverted-repeat and a region of 2616-bp containing alternating copies of 16 and 26 bp repeat elements. Whole mitogenome phylogenies inferred from both nucleotide and amino-acid datasets place M. niger among Melanostomiinae. The need for additional complete mitogenome sequences from the subfamily Malacosteinae is discussed.

Voriconazole Therapeutic Drug Monitoring Among Lung Transplant Recipients Receiving Targeted Therapy for Invasive Aspergillosis

BACKGROUND

Voriconazole is the first line treatment for invasive aspergillosis (IA) Current guidelines suggest performing regular voriconazole therapeutic drug monitoring (TDM) to optimize treatment efficacy. We aimed to determine if TDM was predictive of clinical outcome in LTRs.

METHODS

Retrospective chart review was performed for all LTRs with probable or proven IA, treated with voriconazole monotherapy and who underwent TDM during therapy. Clinical outcome and toxicity were measured at 12 weeks. Classification and regression tree (CART) analysis was used to determine the most predictive voriconazole level thresholds for successful outcome.

RESULTS

118 TDM samples from 30 LTRs with IA were analyzed. Three LTRs were excluded due to early treatment discontinuation. The median TDM level was 1.2 μg/mL (range 0.06-7.3). At 12 weeks, 62% (17/27) of patients had a successful outcome, while 37% (10/27) of patients failed therapy. CART analysis determined that the best predictor for successful outcome was a median TDM level > 0.72 μg/mL. Seventy percent (14/20) of patients with median TDM above 0.72 μg/mL had a successful outcome, compared to 42.9% (3/7) of patients with a median TDM below 0.72 μg/mL (OR 3.11; 95% CI: 0.53-20.4; P = 0.21). CART analysis determined that a TDM level greater than 2.13 μg/mL was predictive of hepatotoxicity.

CONCLUSIONS

Our data suggests that a voriconazole TDM range between 0.72 μg/mL and 2.13 μg/mL may be associated with improved outcomes. Our study is in line with current recommendations on the use of voriconazole TDM in improving outcome and minimizing toxicity in LTR with IA. This article is protected by copyright. All rights reserved.

Nitzschia anatoliensis sp. nov., a cryptic diatom species from the highly alkaline Van Lake (Turkey)

In this article we describe Nitzschia anatoliensis Górecka, Gastineau & Solak sp. nov., an example of a diatom species inhabiting extreme habitats. The new species has been isolated and successfully grown from the highly alkaline Van Lake in East Turkey. The description is based on morphology (light and scanning electron microscopy), the sequencing of its organellar genomes and several molecular phylogenies. This species could easily be overlooked because of its extreme similarity to Nitzschia aurariae but molecular phylogenies indicate that they are only distantly related. Furthermore, molecular data suggest that N. anatoliensis may occur in several alkaline lakes of Asia Minor and Siberia, but was previously misidentified as Nitzschia communis. It also revealed the very close genetic proximity between N. anatoliensis and the endosymbiont of the dinotom Kryptoperidinium foliaceum, providing additional clues on what might have been the original species of diatoms to enter symbiosis.

Mitochondrial and Plastid Genomes of the Monoraphid Diatom Schizostauron trachyderma

We provide for the first time the complete plastid and mitochondrial genomes of a monoraphid diatom: Schizostauron trachyderma. The mitogenome is 41,957 bp in size and displays two group II introns in the cox1 gene. The 187,029 bp plastid genome features the typical quadripartite architecture of diatom genomes. It contains a group II intron in the petB gene that overlaps the large single-copy and the inverted repeat region. There is also a group IB4 intron encoding a putative LAGLIDADG homing endonuclease in the rnl gene. The multigene phylogenies conducted provide more evidence of the proximity between S. trachyderma and fistula-bearing species of biraphid diatoms.

Extreme Enlargement of the Inverted Repeat Region in the Plastid Genomes of Diatoms from the Genus Climaconeis

We sequenced the plastid genomes of three diatoms from the genus Climaconeis, including two strains formerly designated as Climaconeis scalaris. At 208,097 and 216,580 bp, the plastid genomes of the latter strains are the largest ever sequenced among diatoms and their increased size is explained by the massive expansion of the inverted repeat region. Important rearrangements of gene order were identified among the two populations of Climaconeis cf. scalaris. The other sequenced Climaconeis chloroplast genome is 1.5 times smaller compared with those of the Climaconeis cf. scalaris strains and it features an usual quadripartite structure. The extensive structural changes reported here for the genus Climaconeis are compared with those previously observed for other algae and plants displaying large plastid genomes.

Haslea silbo, A Novel Cosmopolitan Species of Blue Diatoms

Specimens of a new species of blue diatoms from the genus Haslea Simonsen were discovered in geographically distant sampling sites, first in the Canary Archipelago, then North Carolina, Gulf of Naples, the Croatian South Adriatic Sea, and Turkish coast of the Eastern Mediterranean Sea. An exhaustive characterization of these specimens, using a combined morphological and genomic approach led to the conclusion that they belong to a single new to science cosmopolitan species, Haslea silbo sp. nov. A preliminary characterization of its blue pigment shows similarities to marennine produced by Haslea ostrearia, as evidenced by UV-visible spectrophotometry and Raman spectrometry. Life cycle stages including auxosporulation were also observed, providing data on the cardinal points of this species. For the two most geographically distant populations (North Carolina and East Mediterranean), complete mitochondrial and plastid genomes were sequenced. The mitogenomes of both strains share a rare atp6 pseudogene, but the number, nature, and positions of the group II introns inside its cox1 gene differ between the two populations. There are also two pairs of genes fused in single ORFs. The plastid genomes are characterized by large regions of recombination with plasmid DNA, which are in both cases located between the ycf35 and psbA genes, but whose content differs between the strains. The two sequenced strains hosts three plasmids coding for putative serine recombinase protein whose sequences are compared, and four out of six of these plasmids were highly conserved.

A gene-rich and compact chloroplast genome of the green alga Nephroselmis pyriformis (N.Carter) Ettl 1982 from the shores of Mersin (Eastern Mediterranean Sea)

We report the complete chloroplast genome of the MED1 strain of Nephroselmis pyriformis from the Eastern Mediterranean Sea. At 111,026 bp, this genome is smaller and more compact than those of Nephroselmis olivacea and Nephroselmis astigmatica, and in contrast to the latter taxa, its inverted repeat contains no complete protein-coding genes. It encodes 3 rRNAs, 33 tRNAs and 94 proteins. Maximum likelihood analysis of a concatenated set of chloroplast genes from green algae belonging to deep-diverging lineages positioned the three Nephroselmis species in a strongly supported clade in which N. pyriformis is sister to N. astigmatica.

Anisakiasis in a Canadian patient with incarcerated epigastric hernia

Anisakiasis is human zoonotic parasitic infection caused by a nematode parasite called Anisakis. This infection is usually reported in Asian countries where consumption of raw seafood is common. Very few cases have been reported in North America. We present the case of a female Canadian patient with an Anisakis larvae in an incarcerated ventral hernia. Cases of Anisakis infections are exceedingly rare in western countries, with very few previous reports describing extra-gastrointestinal cases. Diagnosis is often difficult since the symptoms of anisakiasis are not pathognomonic. As the larvae cannot survive in the body, conservative treatment might be effective in intestinal anisakiasis and surgery is usually performed when complications are encountered. Preventive measures are crucial and include educating the public about the risks of raw fish consumption and the importance of visually inspecting consumed fish and freezing it before ingestion to kill the larvae and prevent the infection.

Complete mitogenomes of the chlorophycean green algae Bulbochaete rectangularis var. hiloensis (Oedogoniales) and Stigeoclonium helveticum (Chaetophorales) provide insight into the sequence of events that led to the acquisition of a reduced-derived pattern of evolution in the Chlamydomonadales and Sphaeropleales

Mitogenome evolution in the Chlorophyceae is characterized by the acquisition of a reduced-derived pattern by the Chlamydomonadales + Sphaeropleales clade. Because no mitogenomes are available for the sister clade Oedogoniales + Chaetophorales + Chaetopeltidales, it remains unclear whether the common ancestor of chlorophycean green algae harbored a reduced-derived or ancestral-type mitogenome. The 70,191 and 46,765-bp mitogenomes reported here for Bulbochaete rectangularis var. hiloensis (Oedogoniales) and Stigeoclonium helveticum (Chaetophorales), respectively, shed light on this question. Both contain the same set of 41 conserved genes, a repertoire lacking numerous protein-coding genes but featuring all 27 tRNA genes typically found in ancestral-type mitogenomes.

Complete mitogenomes of the chlorophyte green algae Scherffelia dubia and Tetraselmis sp. CCMP 881 (Chlorodendrophyceae)

We report here the first mitogenome sequences for the chlorophyte class Chlorodendrophyceae. The mitogenomes of Tetraselmis sp. CCMP 881 and Scherffelia dubia (SAG 17.86) are 46,904 bp and 78,958 bp long, respectively, but their gene repertoires are almost identical. Each genome harbors an inverted repeat (IR). The 14,105-bp IR of S. dubia encodes seven genes in addition to a part of rps19, whereas the 2445-bp IR of Tetraselmis sp. CCMP 881 contains a single gene. Considering that an IR has also been found in the mitogenomes of certain earlier-diverging chlorophytes, the IRs of chlorodendrophycean algae probably represent ancestral features.

Complete mitogenome of the chlorophyte green alga Marsupiomonas sp. NIES 1824 (Pedinophyceae)

The 25,137-bp mitogenome of the green alga Pedinomonas minor (Pedinomonadales, Pedinophyceae), which belongs to a basal class of the core Chlorophyta, is unusual in displaying a reduced gene content as well as other derived traits. Here, we present the mitogenome of Marsupiomonas sp. NIES 1824 (Marsupiomonadales, Pedinophyceae). Despite its smaller size, this 24,252-bp genome encodes twice as many genes (39) as its P. minor homolog. Besides gradual gene erosion, our comparative analyses revealed that major changes in GC content and codon usage led to the gain of distinct, noncanonical genetic codes during evolution of the mitogenome in the Pedinophyceae.

Complete mitogenomes of the marine picoplanktonic green algae Prasinoderma sp. MBIC 10622 and Prasinococcus capsulatus CCMP 1194 (Palmophyllophyceae)

Marine picoalgae from the Prasinococcales order occupy the deepest branch of the Chlorophyta (Palmophyllophyceae). Here, we describe the mitogenomes of Prasinoderma sp. MBIC 10622 and Prasinococcus capsulatus CCMP 1194. At 37,590 and 41,006 bp, respectively, they are smaller than their Prasinoderma coloniale homolog and unlike the latter, lack an inverted repeat. The intronless Prasinoderma sp. mitogenome possesses the largest gene repertoire (68) among all chlorophytes examined to date. At the gene order level, it displays more ancestral traits than its prasinococcalean homologs, closely resembling the mitogenomes of Mamiellophyceae. Remarkably, the P. capsulatus mitogenome features a trans-spliced group II intron.

Complete mitogenome of the streptophyte green alga Coleochaete scutata (Coleochaetophyceae)

We have sequenced the mitogenome of Coleochaete scutata strain SAG 110.80M. This mitogenome is the largest among the streptophyte green algae examined to date. At 242,024 bp, it is 4.3-fold larger than the mitogenome of Chaetosphaeridium globosum, the only other mitogenome available for the Coleochaetophyceae. This size difference is mainly explained by differences in the abundance of introns and in the length of intergenic regions containing vestiges of coding sequences thought to be of foreign origin. With 31 group I and 26 group II introns, the C. scutata mitogenome is the most intron-rich organelle genome known among streptophyte algae.

Complete chloroplast genome of the tiny marine diatom Nanofrustulum shiloi (Bacillariophyta) from the Adriatic Sea

We report the chloroplast genome sequence of Nanofrustulum shiloi, a tiny araphid pennate diatom collected from the Adriatic Sea. The 160,994-bp N. shiloi genome displays a quadripartite structure and its gene repertoire resembles those of other diatom chloroplast genomes. Besides the genes located in the inverted repeat, psbY is duplicated. A gene-poor region in the large single-copy region contains multiple ORFs sharing sequence similarities with plasmids and chloroplast ORFs found in other diatom species. The genome features a single intron, a group II intron in petB. Phylogenomic analysis identified N. shiloi at a basal position within the araphid 2 clade.

Complete mitogenome of the invasive bivalve Rangia cuneata

We sequenced the complete mitogenome of a Baltic Sea specimen of the invasive bivalve Rangia cuneata. The mitogenome is 18,993 bp long and encodes 13 proteins, 2 rRNAs, and 22 tRNAs. A large intergenic region between ND6 and ND2 contains an open reading frame that may originate from duplication of ND2. The R. cuneata mitogenome is rearranged in gene order relative to previously sequenced mitogenomes of Mactroidea.

Tracing the Evolution of the Plastome and Mitogenome in the Chloropicophyceae Uncovered Convergent tRNA Gene Losses and a Variant Plastid Genetic Code

The tiny green algae belonging to the Chloropicophyceae play a key role in marine phytoplankton communities; this newly erected class of prasinophytes comprises two genera (Chloropicon and Chloroparvula) containing each several species. We sequenced the plastomes and mitogenomes of eight Chloropicon and five Chloroparvula species to better delineate the phylogenetic affinities of these taxa and to infer the suite of changes that their organelle genomes sustained during evolution. The relationships resolved in organelle-based phylogenomic trees were essentially congruent with previously reported rRNA trees, and similar evolutionary trends but distinct dynamics were identified for the plastome and mitogenome. Although the plastome sustained considerable changes in gene content and order at the time the two genera split, subsequently it remained stable and maintained a very small size. The mitogenome, however, was remodeled more gradually and showed more fluctuation in size, mainly as a result of expansions/contractions of intergenic regions. Remarkably, the plastome and mitogenome lost a common set of three tRNA genes, with the trnI(cau) and trnL(uaa) losses being accompanied with important variations in codon usage. Unexpectedly, despite the disappearance of trnI(cau) from the plastome in the Chloroparvula lineage, AUA codons (the codons recognized by this gene product) were detected in certain plastid genes. By comparing the sequences of plastid protein-coding genes from chloropicophycean and phylogenetically diverse chlorophyte algae with those of the corresponding predicted proteins, we discovered that the AUA codon was reassigned from isoleucine to methionine in Chloroparvula. This noncanonical genetic code has not previously been uncovered in plastids.

Screening for Exotic Forest Pathogens to Increase Survey Capacity Using Metagenomics

Anthropogenic activities have a major impact on the global environment. Canada's natural resources are threatened by the spread of fungal pathogens, which is facilitated by agricultural practices and international trade. Fungi are introduced to new environments and sometimes become established, in which case they can cause disease outbreaks resulting in extensive forest decline. Here, we describe how a nationwide sample collection strategy coupled to next-generation sequencing (NGS) (i.e., metagenomics) can achieve fast and comprehensive screening for exotic invasive species. This methodology can help provide guidance to phytopathology stakeholders such as regulatory agencies. Several regulated invasive species were monitored by processing field samples collected over 3 years (2013 to 2015) near high-risk areas across Canada. Fifteen sequencing runs were required on the Ion Torrent platform to process 398 samples that yielded 45 million reads. High-throughput screening of fungal and oomycete operational taxonomic units using customized fungi-specific ribosomal internal transcribed spacer 1 barcoded primers was performed. Likewise, Phytophthora-specific barcoded primers were used to amplify the adenosine triphosphate synthase subunit 9-nicotinamide adenine dinucleotide dehydrogenase subunit 9 spacer. Several Phytophthora spp. were detected by NGS and confirmed by species-specific quantitative polymerase chain reaction (qPCR) assays. The target species Heterobasidion annosum sensu stricto could be detected only through metagenomics. We demonstrated that screening target species using a variety of sampling techniques and NGS-the results of which were validated by qPCR-has the potential to increase survey capacity and detection sensitivity, reduce hands-on time and costs, and assist regulatory agencies to identify ports of entry. Considering that early detection and prevention are the keys in mitigating invasive species damage, our method represents a substantial asset in plant pathology management.

Complete mitochondrial genome of a rare diatom (Bacillariophyta) Proschkinia and its phylogenetic and taxonomic implications

We obtained the complete mitogenome of Proschkinia sp. strain SZCZR1824, a strain belonging to a poorly known diatom genus with no previous molecular data. This genome is 48,863 bp long, with two group I introns in rnl and three group II introns in cox1. Using mitogenomic data, Proschkinia sp. was recovered with Fistulifera solaris, far distant from Navicula and Nitzschia, two genera with which Proschkinia has sometimes been associated based on morphology.

Complete mitogenome of Cerithidea obtusa, the red chut-chut snail from the Cần Giờ Mangrove in Vietnam

We sequenced the complete mitogenome of the red chut-chut snail Cerithidea obtusa, from the Cần Giờ mangrove in Vietnam. The mitogenome is 15,708 bp long. It is colinear with the mitogenomes of other members of the superfamily Cerithioidea, and the maximum-likelihood phylogeny obtained with the cox1, cox2 and cox3 genes of several Caenogastropoda associated all Cerithioidea together inside a strongly supported clade.

Mitogenome sequence of a Black Sea isolate of the kinetoplastid Bodo saltans

We obtained the mitogenome sequence of a Black Sea isolate of the kinetoplastid Bodo saltans. This sequence consists of two contigs totaling 24,925 bp and encodes ten protein-coding genes, one conserved ORF and one rRNA gene. Alignment of the Black Sea mitogenome with the limited sequence data currently available in public databases for another strain of B. saltans revealed significant genetic divergence between the two isolates. Maximum likelihood phylogenetic inference clearly resolved the Bodonidae from the Trypanosomatidae.

Genome-wide organellar analyses from the hornwort Leiosporoceros dussii show low frequency of RNA editing

Because hornworts occupy a pivotal position in early land colonization as sister to other bryophytes, sister to tracheophytes, or sister to all other land plants, a renewed interest has arisen in their phylogenetic diversity, morphology, and genomes. To date, only five organellar genome sequences are available for hornworts. We sequenced the plastome (155,956 bp) and mitogenome (212,153 bp) of the hornwort Leiosporoceros dussii, the sister taxon to all hornworts. The Leiosporoceros organellar genomes show conserved gene structure and order with respect to the other hornworts and other bryophytes. Additionally, using RNA-seq data we quantified the frequency of RNA-editing events (the canonical C-to-U and the reverse editing U-to-C) in both organellar genomes. In total, 109 sites were found in the plastome and 108 in the mitogenome, respectively. The proportion of edited sites corresponds to 0.06% of the plastome and 0.05% of the mitogenome (in reference to the total genome size), in contrast to 0.58% of edited sites in the plastome of Anthoceros angustus (161,162 bp). All edited sites in the plastome and 88 of 108 sites in the mitogenome are C-to-U conversions. Twenty reverse edited sites (U-to-C conversions) were found in the mitogenome (17.8%) and none in the plastome. The low frequency of RNA editing in Leiosporoceros, which is nearly 88% less than in the plastome of Anthoceros and the mitogenome of Nothoceros, indicates that the frequency of RNA editing has fluctuated during hornwort diversification. Hornworts are a pivotal land plant group to unravel the genomic implications of RNA editing and its maintenance despite the evident evolutionary disadvantages.

The complete mitochondrial DNA of the tropical oyster Crassostrea belcheri from the Cần Giò' mangrove in Vietnam

The complete mitochondrial genome of the oyster Crassostrea belcheri from the Cần Giò' mangrove in Vietnam has been sequenced. It consists of a circular DNA molecule of 21020 base pairs (bp), coding for 12 proteins, 20 transfer RNAs, and two ribosomal RNAs. Like the mitogenomes of Crassostrea iredalei and Crassostrea sp. DB1, it contains a non-coding region and two ORFs. The C. belcheri mitogenome provides information that could improve the molecular phylogeny of Asian oysters and be useful to the development of oyster aquaculture in South East Asia.

Implementing a web-based introductory bioinformatics course for non-bioinformaticians that incorporates practical exercises

A recent scientific discipline, bioinformatics, defined as using informatics for the study of biological problems, is now a requirement for the study of biological sciences. Bioinformatics has become such a powerful and popular discipline that several academic institutions have created programs in this field, allowing students to become specialized. However, biology students who are not involved in a bioinformatics program also need a solid toolbox of bioinformatics software and skills. Therefore, we have developed a completely online bioinformatics course for non-bioinformaticians, entitled "BIF-1901 Introduction à la bio-informatique et à ses outils (Introduction to bioinformatics and bioinformatics tools)," given by the Department of Biochemistry, Microbiology, and Bioinformatics of Université Laval (Quebec City, Canada). This course requires neither a bioinformatics background nor specific skills in informatics. The underlying main goal was to produce a completely online up-to-date bioinformatics course, including practical exercises, with an intuitive pedagogical framework. The course, BIF-1901, was conceived to cover the three fundamental aspects of bioinformatics: (1) informatics, (2) biological sequence analysis, and (3) structural bioinformatics. This article discusses the content of the modules, the evaluations, the pedagogical framework, and the challenges inherent to a multidisciplinary, fully online course. © 2017 by The International Union of Biochemistry and Molecular Biology, 46(1):31-38, 2018.

Divergent copies of the large inverted repeat in the chloroplast genomes of ulvophycean green algae

The chloroplast genomes of many algae and almost all land plants carry two identical copies of a large inverted repeat (IR) sequence that can pair for flip-flop recombination and undergo expansion/contraction. Although the IR has been lost multiple times during the evolution of the green algae, the underlying mechanisms are still largely unknown. A recent comparison of IR-lacking and IR-containing chloroplast genomes of chlorophytes from the Ulvophyceae (Ulotrichales) suggested that differential elimination of genes from the IR copies might lead to IR loss. To gain deeper insights into the evolutionary history of the chloroplast genome in the Ulvophyceae, we analyzed the genomes of Ignatius tetrasporus and Pseudocharacium americanum (Ignatiales, an order not previously sampled), Dangemannia microcystis (Oltmannsiellopsidales), Pseudoneochloris marina (Ulvales) and also Chamaetrichon capsulatum and Trichosarcina mucosa (Ulotrichales). Our comparison of these six chloroplast genomes with those previously reported for nine ulvophyceans revealed unsuspected variability. All newly examined genomes feature an IR, but remarkably, the copies of the IR present in the Ignatiales, Pseudoneochloris, and Chamaetrichon diverge in sequence, with the tRNA genes from the rRNA operon missing in one IR copy. The implications of this unprecedented finding for the mechanism of IR loss and flip-flop recombination are discussed.

Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae)

Background

The chloroplast genome sustained extensive changes in architecture during the evolution of the Chlorophyceae, a morphologically and ecologically diverse class of green algae belonging to the Chlorophyta; however, the forces driving these changes are poorly understood. The five orders recognized in the Chlorophyceae form two major clades: the CS clade consisting of the Chlamydomonadales and Sphaeropleales, and the OCC clade consisting of the Oedogoniales, Chaetophorales, and Chaetopeltidales. In the OCC clade, considerable variations in chloroplast DNA (cpDNA) structure, size, gene order, and intron content have been observed. The large inverted repeat (IR), an ancestral feature characteristic of most green plants, is present in Oedogonium cardiacum (Oedogoniales) but is lacking in the examined members of the Chaetophorales and Chaetopeltidales. Remarkably, the Oedogonium 35.5-kb IR houses genes that were putatively acquired through horizontal DNA transfer. To better understand the dynamics of chloroplast genome evolution in the Oedogoniales, we analyzed the cpDNA of a second representative of this order, Oedocladium carolinianum.

Methods

The Oedocladium cpDNA was sequenced and annotated. The evolutionary distances separating Oedocladium and Oedogonium cpDNAs and two other pairs of chlorophycean cpDNAs were estimated using a 61-gene data set. Phylogenetic analysis of an alignment of group IIA introns from members of the OCC clade was performed. Secondary structures and insertion sites of oedogonialean group IIA introns were analyzed.

Results

The 204,438-bp Oedocladium genome is 7.9 kb larger than the Oedogonium genome, but its repertoire of conserved genes is remarkably similar and gene order differs by only one reversal. Although the 23.7-kb IR is missing the putative foreign genes found in Oedogonium, it contains sequences coding for a putative phage or bacterial DNA primase and a hypothetical protein. Intergenic sequences are 1.5-fold longer and dispersed repeats are more abundant, but a smaller fraction of the Oedocladium genome is occupied by introns. Six additional group II introns are present, five of which lack ORFs and carry highly similar sequences to that of the ORF-less IIA intron shared with Oedogonium. Secondary structure analysis of the group IIA introns disclosed marked differences in the exon-binding sites; however, each intron showed perfect or nearly perfect base pairing interactions with its target site.

Discussion

Our results suggest that chloroplast genes rearrange more slowly in the Oedogoniales than in the Chaetophorales and raise questions as to what was the nature of the foreign coding sequences in the IR of the common ancestor of the Oedogoniales. They provide the first evidence for intragenomic proliferation of group IIA introns in the Viridiplantae, revealing that intron spread in the Oedocladium lineage likely occurred by retrohoming after sequence divergence of the exon-binding sites.

Mitochondrion-to-Chloroplast DNA Transfers and Intragenomic Proliferation of Chloroplast Group II Introns in Gloeotilopsis Green Algae (Ulotrichales, Ulvophyceae)

To probe organelle genome evolution in the Ulvales/Ulotrichales clade, the newly sequenced chloroplast and mitochondrial genomes of Gloeotilopsis planctonica and Gloeotilopsis sarcinoidea (Ulotrichales) were compared with those of Pseudendoclonium akinetum (Ulotrichales) and of the few other green algae previously sampled in the Ulvophyceae. At 105,236 bp, the G planctonica mitochondrial DNA (mtDNA) is the largest mitochondrial genome reported so far among chlorophytes, whereas the 221,431-bp G planctonica and 262,888-bp G sarcinoidea chloroplast DNAs (cpDNAs) are the largest chloroplast genomes analyzed among the Ulvophyceae. Gains of non-coding sequences largely account for the expansion of these genomes. Both Gloeotilopsis cpDNAs lack the inverted repeat (IR) typically found in green plants, indicating that two independent IR losses occurred in the Ulvales/Ulotrichales. Our comparison of the Pseudendoclonium and Gloeotilopsis cpDNAs offered clues regarding the mechanism of IR loss in the Ulotrichales, suggesting that internal sequences from the rDNA operon were differentially lost from the two original IR copies during this process. Our analyses also unveiled a number of genetic novelties. Short mtDNA fragments were discovered in two distinct regions of the G sarcinoidea cpDNA, providing the first evidence for intracellular inter-organelle gene migration in green algae. We identified for the first time in green algal organelles, group II introns with LAGLIDADG ORFs as well as group II introns inserted into untranslated gene regions. We discovered many group II introns occupying sites not previously documented for the chloroplast genome and demonstrated that a number of them arose by intragenomic proliferation, most likely through retrohoming.

Comparative Chloroplast Genome Analyses of Streptophyte Green Algae Uncover Major Structural Alterations in the Klebsormidiophyceae, Coleochaetophyceae and Zygnematophyceae

The Streptophyta comprises all land plants and six main lineages of freshwater green algae: Mesostigmatophyceae, Chlorokybophyceae, Klebsormidiophyceae, Charophyceae, Coleochaetophyceae and Zygnematophyceae. Previous comparisons of the chloroplast genome from nine streptophyte algae (including four zygnematophyceans) revealed that, although land plant chloroplast DNAs (cpDNAs) inherited most of their highly conserved structural features from green algal ancestors, considerable cpDNA changes took place during the evolution of the Zygnematophyceae, the sister group of land plants. To gain deeper insights into the evolutionary dynamics of the chloroplast genome in streptophyte algae, we sequenced the cpDNAs of nine additional taxa: two klebsormidiophyceans (Entransia fimbriata and Klebsormidium sp. SAG 51.86), one coleocheatophycean (Coleochaete scutata) and six zygnematophyceans (Cylindrocystis brebissonii, Netrium digitus, Roya obtusa, Spirogyra maxima, Cosmarium botrytis and Closterium baillyanum). Our comparative analyses of these genomes with their streptophyte algal counterparts indicate that the large inverted repeat (IR) encoding the rDNA operon experienced loss or expansion/contraction in all three sampled classes and that genes were extensively shuffled in both the Klebsormidiophyceae and Zygnematophyceae. The klebsormidiophycean genomes boast greatly expanded IRs, with the Entransia 60,590-bp IR being the largest known among green algae. The 206,025-bp Entransia cpDNA, which is one of the largest genome among streptophytes, encodes 118 standard genes, i.e., four additional genes compared to its Klebsormidium flaccidum homolog. We inferred that seven of the 21 group II introns usually found in land plants were already present in the common ancestor of the Klebsormidiophyceae and its sister lineages. At 107,236 bp and with 117 standard genes, the Coleochaete IR-less genome is both the smallest and most compact among the streptophyte algal cpDNAs analyzed thus far; it lacks eight genes relative to its Chaetosphaeridium globosum homolog, four of which represent unique events in the evolutionary scenario of gene losses we reconstructed for streptophyte algae. The 10 compared zygnematophycean cpDNAs display tremendous variations at all levels, except gene content. During zygnematophycean evolution, the IR disappeared a minimum of five times, the rDNA operon was broken at four distinct sites, group II introns were lost on at least 43 occasions, and putative foreign genes, mainly of phage/viral origin, were gained.

Distinctive Architecture of the Chloroplast Genome in the Chlorodendrophycean Green Algae Scherffelia dubia and Tetraselmis sp. CCMP 881

The Chlorodendrophyceae is a small class of green algae belonging to the core Chlorophyta, an assemblage that also comprises the Pedinophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Here we describe for the first time the chloroplast genomes of chlorodendrophycean algae (Scherffelia dubia, 137,161 bp; Tetraselmis sp. CCMP 881, 100,264 bp). Characterized by a very small single-copy (SSC) region devoid of any gene and an unusually large inverted repeat (IR), the quadripartite structures of the Scherffelia and Tetraselmis genomes are unique among all core chlorophytes examined thus far. The lack of genes in the SSC region is offset by the rich and atypical gene complement of the IR, which includes genes from the SSC and large single-copy regions of prasinophyte and streptophyte chloroplast genomes having retained an ancestral quadripartite structure. Remarkably, seven of the atypical IR-encoded genes have also been observed in the IRs of pedinophycean and trebouxiophycean chloroplast genomes, suggesting that they were already present in the IR of the common ancestor of all core chlorophytes. Considering that the relationships among the main lineages of the core Chlorophyta are still unresolved, we evaluated the impact of including the Chlorodendrophyceae in chloroplast phylogenomic analyses. The trees we inferred using data sets of 79 and 108 genes from 71 chlorophytes indicate that the Chlorodendrophyceae is a deep-diverging lineage of the core Chlorophyta, although the placement of this class relative to the Pedinophyceae remains ambiguous. Interestingly, some of our phylogenomic trees together with our comparative analysis of gene order data support the monophyly of the Trebouxiophyceae, thus offering further evidence that the previously observed affiliation between the Chlorellales and Pedinophyceae is the result of systematic errors in phylogenetic reconstruction.

Dynamic Evolution of the Chloroplast Genome in the Green Algal Classes Pedinophyceae and Trebouxiophyceae

Previous studies of trebouxiophycean chloroplast genomes revealed little information regarding the evolutionary dynamics of this genome because taxon sampling was too sparse and the relationships between the sampled taxa were unknown. We recently sequenced the chloroplast genomes of 27 trebouxiophycean and 2 pedinophycean green algae to resolve the relationships among the main lineages recognized for the Trebouxiophyceae. These taxa and the previously sampled members of the Pedinophyceae and Trebouxiophyceae are included in the comparative chloroplast genome analysis we report here. The 38 genomes examined display considerable variability at all levels, except gene content. Our results highlight the high propensity of the rDNA-containing large inverted repeat (IR) to vary in size, gene content and gene order as well as the repeated losses it experienced during trebouxiophycean evolution. Of the seven predicted IR losses, one event demarcates a superclade of 11 taxa representing 5 late-diverging lineages. IR expansions/contractions account not only for changes in gene content in this region but also for changes in gene order and gene duplications. Inversions also led to gene rearrangements within the IR, including the reversal or disruption of the rDNA operon in some lineages. Most of the 20 IR-less genomes are more rearranged compared with their IR-containing homologs and tend to show an accelerated rate of sequence evolution. In the IR-less superclade, several ancestral operons were disrupted, a few genes were fragmented, and a subgroup of taxa features a G+C-biased nucleotide composition. Our analyses also unveiled putative cases of gene acquisitions through horizontal transfer.

Six newly sequenced chloroplast genomes from prasinophyte green algae provide insights into the relationships among prasinophyte lineages and the diversity of streamlined genome architecture in picoplanktonic species

Background

Because they represent the earliest divergences of the Chlorophyta, the morphologically diverse unicellular green algae making up the prasinophytes hold the key to understanding the nature of the first viridiplants and the evolutionary patterns that accompanied the radiation of chlorophytes. Nuclear-encoded 18S rDNA phylogenies unveiled nine prasinophyte clades (clades I through IX) but their branching order is still uncertain. We present here the newly sequenced chloroplast genomes of Nephroselmis astigmatica (clade III) and of five picoplanktonic species from clade VI (Prasinococcus sp. CCMP 1194, Prasinophyceae sp. MBIC 106222 and Prasinoderma coloniale) and clade VII (Picocystis salinarum and Prasinophyceae sp. CCMP 1205). These chloroplast DNAs (cpDNAs) were compared with those of the six previously sampled prasinophytes (clades I, II, III and V) in order to gain information both on the relationships among prasinophyte lineages and on chloroplast genome evolution.

Results

Varying from 64.3 to 85.6 kb in size and encoding 100 to 115 conserved genes, the cpDNAs of the newly investigated picoplanktonic species are substantially smaller than those observed for larger-size prasinophytes, are economically packed and contain a reduced gene content. Although the Nephroselmis and Picocystis cpDNAs feature a large inverted repeat encoding the rRNA operon, gene partitioning among the single copy regions is remarkably different. Unexpectedly, we found that all three species from clade VI (Prasinococcales) harbor chloroplast genes not previously documented for chlorophytes (ndhJ, rbcR, rpl21, rps15, rps16 and ycf66) and that Picocystis contains a trans-spliced group II intron. The phylogenies inferred from cpDNA-encoded proteins are essentially congruent with 18S rDNA trees, resolving with robust support all six examined prasinophyte lineages, with the exception of the Pycnococcaceae.

Conclusions

Our results underscore the high variability in genome architecture among prasinophyte lineages, highlighting the strong pressure to maintain a small and compact chloroplast genome in picoplanktonic species. The unique set of six chloroplast genes found in the Prasinococcales supports the ancestral status of this lineage within the prasinophytes. The widely diverging traits uncovered for the clade-VII members (Picocystis and Prasinophyceae sp. CCMP 1205) are consistent with their resolution as separate lineages in the chloroplast phylogeny.

Chloroplast phylogenomic analysis resolves deep-level relationships within the green algal class Trebouxiophyceae

Background

The green algae represent one of the most successful groups of photosynthetic eukaryotes, but compared to their land plant relatives, surprisingly little is known about their evolutionary history. This is in great part due to the difficulty of recognizing species diversity behind morphologically similar organisms. The Trebouxiophyceae is a species-rich class of the Chlorophyta that includes symbionts (e.g. lichenized algae) as well as free-living green algae. Members of this group display remarkable ecological variation, occurring in aquatic, terrestrial and aeroterrestrial environments. Because a reliable backbone phylogeny is essential to understand the evolutionary history of the Trebouxiophyceae, we sought to identify the relationships among the major trebouxiophycean lineages that have been previously recognized in nuclear-encoded 18S rRNA phylogenies. To this end, we used a chloroplast phylogenomic approach.

Results

We determined the sequences of 29 chlorophyte chloroplast genomes and assembled amino acid and nucleotide data sets derived from 79 chloroplast genes of 61 chlorophytes, including 35 trebouxiophyceans. The amino acid- and nucleotide-based phylogenies inferred using maximum likelihood and Bayesian methods and various models of sequence evolution revealed essentially the same relationships for the trebouxiophyceans. Two major groups were identified: a strongly supported clade of 29 taxa (core trebouxiophyceans) that is sister to the Chlorophyceae + Ulvophyceae and a clade comprising the Chlorellales and Pedinophyceae that represents a basal divergence relative to the former group. The core trebouxiophyceans form a grade of strongly supported clades that include a novel lineage represented by the desert crust alga Pleurastrosarcina brevispinosa. The assemblage composed of the Oocystis and Geminella clades is the deepest divergence of the core trebouxiophyceans. Like most of the chlorellaleans, early-diverging core trebouxiophyceans are predominantly planktonic species, whereas core trebouxiophyceans occupying more derived lineages are mostly terrestrial or aeroterrestrial algae.

Conclusions

Our phylogenomic study provides a solid foundation for addressing fundamental questions related to the biology and ecology of the Trebouxiophyceae. The inferred trees reveal that this class is not monophyletic; they offer new insights not only into the internal structure of the class but also into the lifestyle of its founding members and subsequent adaptations to changing environments.

Tracing the evolution of streptophyte algae and their mitochondrial genome

Six monophyletic groups of charophycean green algae are recognized within the Streptophyta. Although incongruent with earlier studies based on genes from three cellular compartments, chloroplast and nuclear phylogenomic analyses have resolved identical relationships among these groups, placing the Zygnematales or the Zygnematales + Coleochaetales as sister to land plants. The present investigation aimed at determining whether this consensus view is supported by the mitochondrial genome and at gaining insight into mitochondrial DNA (mtDNA) evolution within and across streptophyte algal lineages and during the transition toward the first land plants. We present here the newly sequenced mtDNAs of representatives of the Klebsormidiales (Entransia fimbriata and Klebsormidium spec.) and Zygnematales (Closterium baillyanum and Roya obtusa) and compare them with their homologs in other charophycean lineages as well as in selected embryophyte and chlorophyte lineages. Our results indicate that important changes occurred at the levels of genome size, gene order, and intron content within the Zygnematales. Although the representatives of the Klebsormidiales display more similarity in genome size and intron content, gene order seems more fluid and gene losses more frequent than in other charophycean lineages. In contrast, the two members of the Charales display an extremely conservative pattern of mtDNA evolution. Collectively, our analyses of gene order and gene content and the phylogenies we inferred from 40 mtDNA-encoded proteins failed to resolve the relationships among the Zygnematales, Coleochaetales, and Charales; however, they are consistent with previous phylogenomic studies in favoring that the morphologically complex Charales are not sister to land plants.

The mitochondrial genome of the prasinophyte Prasinoderma coloniale reveals two trans-spliced group I introns in the large subunit rRNA gene

Organelle genes are often interrupted by group I and or group II introns. Splicing of these mobile genetic occurs at the RNA level via serial transesterification steps catalyzed by the introns'own tertiary structures and, sometimes, with the help of external factors. These catalytic ribozymes can be found in cis or trans configuration, and although trans-arrayed group II introns have been known for decades, trans-spliced group I introns have been reported only recently. In the course of sequencing the complete mitochondrial genome of the prasinophyte picoplanktonic green alga Prasinoderma coloniale CCMP 1220 (Prasinococcales, clade VI), we uncovered two additional cases of trans-spliced group I introns. Here, we describe these introns and compare the 54,546 bp-long mitochondrial genome of Prasinoderma with those of four other prasinophytes (clades II, III and V). This comparison underscores the highly variable mitochondrial genome architecture in these ancient chlorophyte lineages. Both Prasinoderma trans-spliced introns reside within the large subunit rRNA gene (rnl) at positions where cis-spliced relatives, often containing homing endonuclease genes, have been found in other organelles. In contrast, all previously reported trans-spliced group I introns occur in different mitochondrial genes (rns or coxI). Each Prasinoderma intron is fragmented into two pieces, forming at the RNA level a secondary structure that resembles those of its cis-spliced counterparts. As observed for other trans-spliced group I introns, the breakpoint of the first intron maps to the variable loop L8, whereas that of the second is uniquely located downstream of P9.1. The breakpoint In each Prasinoderma intron corresponds to the same region where the open reading frame (ORF) occurs when present in cis-spliced orthologs. This correlation between the intron breakpoint and the ORF location in cis-spliced orthologs also holds for other trans-spliced introns; we discuss the possible implications of this interesting observation for trans-splicing of group I introns.