Variant Calling from RNA-seq Data Using the GATK Joint Genotyping Workflow.

The Genome Analysis Toolkit (GATK) developed at the Broad Institute provides state-of-the-art pipelines for germline and somatic variant discovery and genotyping. Unfortunately, the fully validated GATK pipeline for calling variant on RNAseq data is a Per-sample workflow that does not include the recent improvements seen in modern workflows, especially the possibility to perform joint genotyping analysis. Here, we describe how modern GATK commands from distinct workflows can be combined to call variants on RNAseq samples. We provide a detailed tutorial that starts with raw RNAseq reads and ends with filtered variants, of which some were shown to be associated with bovine paratuberculosis.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages.

BACKGROUND: Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne's disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (-); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip. RESULTS: The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/-) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10- 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids. CONCLUSION: In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Resistance determinants and their genetic context in enterobacteria from a longitudinal study of pigs reared under various husbandry conditions.

Pigs are major reservoirs of resistant Enterobacteriaceae that can reach humans through consumption of contaminated meat or vegetables grown in manure-fertilized soil. Samples were collected from sows during lactation and their piglets at five time points spanning the production cycle. Cefotaxime-resistant bacteria were quantified and isolated from feed, feces, manures and carcasses of pigs reared with penicillin-using or antibiotic-free husbandries. The isolates were characterized by antibiotic susceptibility testing, whole genome sequencing and conjugation assays. The extended spectrum ß-lactamase (ESBL) phenotype was more frequent in isolates originating from antibiotic-free animals, while the bacteria isolated from penicillin-using animals were on average resistant to a greater number of antibiotics. The ESBL-encoding genes identified were bla CTX-M-1, bla CTX-M-15 and bla CMY-2 and they co-localised on plasmids with various genes encoding resistance to ß-lactams, co-trimoxazole, phenicols and tetracycline, all antibiotics used in pig production. Groups of genes conferring the observed resistance and the mobile elements disseminating multidrug resistance were determined. The observed resistance to ß-lactams was mainly due to the complementary actions of penicillin-binding proteins, an efflux pump and ß-lactamases. Most resistance determinants were shared by animals raised with or without antimicrobials. This suggests a key contribution of indigenous enterobacteria maternally transmitted along the sow lineage, regardless of antimicrobial use. It is unclear if the antimicrobial resistance observed in the enterobacteria populations of the commercial pig herds studied were present before the use of antibiotics, or the extent to which historical antimicrobial use exerted a selective pressure defining the resistant bacterial populations in farms using penicillin prophylaxis.Importance: Antimicrobial resistance is a global threat that needs to be fought on numerous fronts along the One Health continuum. Vast quantities of antimicrobials are used in agriculture to ensure animal welfare and productivity, and are arguably a driving force for the persistence of environmental and food-borne resistant bacteria. This study evaluated the impact of conventional, organic and other antibiotic-free husbandry practices on the frequency and nature of antimicrobial resistance genes and multidrug resistant enterobacteria. It provides knowledge about the relative contribution of specific resistance determinants to observed antibiotic resistance. It also showed the clear co-selection of genes coding for extended-spectrum beta-lactamases and genes coding for the resistance to antibiotics commonly used for prophylaxis or in curative treatments in pig operations.

The GATK joint genotyping workflow is appropriate for calling variants in RNA-seq experiments.

The Genome Analysis Toolkit (GATK) is a popular set of programs for discovering and genotyping variants from next-generation sequencing data. The current GATK recommendation for RNA sequencing (RNA-seq) is to perform variant calling from individual samples, with the drawback that only variable positions are reported. Versions 3.0 and above of GATK offer the possibility of calling DNA variants on cohorts of samples using the HaplotypeCaller algorithm in Genomic Variant Call Format (GVCF) mode. Using this approach, variants are called individually on each sample, generating one GVCF file per sample that lists genotype likelihoods and their genome annotations. In a second step, variants are called from the GVCF files through a joint genotyping analysis. This strategy is more flexible and reduces computational challenges in comparison to the traditional joint discovery workflow. Using a GVCF workflow for mining SNP in RNA-seq data provides substantial advantages, including reporting homozygous genotypes for the reference allele as well as missing data. Taking advantage of RNA-seq data derived from primary macrophages isolated from 50 cows, the GATK joint genotyping method for calling variants on RNA-seq data was validated by comparing this approach to a so-called "per-sample" method. In addition, pair-wise comparisons of the two methods were performed to evaluate their respective sensitivity, precision and accuracy using DNA genotypes from a companion study including the same 50 cows genotyped using either genotyping-by-sequencing or with the Bovine SNP50 Beadchip (imputed to the Bovine high density). Results indicate that both approaches are very close in their capacity of detecting reference variants and that the joint genotyping method is more sensitive than the per-sample method. Given that the joint genotyping method is more flexible and technically easier, we recommend this approach for variant calling in RNA-seq experiments.

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.

Low-depth genotyping-by-sequencing (GBS) in a bovine population: strategies to maximize the selection of high quality genotypes and the accuracy of imputation.

BACKGROUND: Genotyping-by-sequencing (GBS) has emerged as a powerful and cost-effective approach for discovering and genotyping single-nucleotide polymorphisms. The GBS technique was largely used in crop species where its low sequence coverage is not a drawback for calling genotypes because inbred lines are almost homozygous. In contrast, only a few studies used the GBS technique in animal populations (with sizeable heterozygosity rates) and many of those that have been published did not consider the quality of the genotypes produced by the bioinformatic pipelines. To improve the sequence coverage of the fragments, an alternative GBS preparation protocol that includes selective primers during the PCR amplification step has been recently proposed. In this study, we compared this modified protocol with the conventional two-enzyme GBS protocol. We also described various procedures to maximize the selection of high quality genotypes and to increase the accuracy of imputation. RESULTS: The in silico digestions of the bovine genome showed that the combination of PstI and MspI is more suitable for sequencing bovine GBS libraries than the use of single digestions with PstI or ApeKI. The sequencing output of the GBS libraries generated a total of 123,666 variants with the selective-primer approach and 272,103 variants with the conventional approach. Validating our data with genotypes obtained from mass spectrometry and Illumina's bovine SNP50 array, we found that the genotypes produced by the conventional GBS method were concordant with those produced by these alternative genotyping methods, whereas the selective-primer method failed to call heterozygotes with confidence. Our results indicate that high accuracy in genotype calling (>97%) can be obtained using low read-depth thresholds (3 to 5 reads) provided that markers are simultaneously filtered for genotype quality scores. We also show that factors such as the minimum call rate and the minor allele frequency positively influence the accuracy of imputation of missing GBS data. The highest accuracies (around 85%) of imputed GBS markers were obtained with the FIMPUTE program when GBS and SNP50 array genotypes were combined (80,190 to 100,297 markers) before imputation. CONCLUSIONS: We discovered that the conventional two-enzyme GBS protocol could produce a large number of high-quality genotypes provided that appropriate filtration criteria were used. In contrast, the selective-primer approach resulted in a substantial proportion of miscalled genotypes and should be avoided for livestock genotyping studies. Overall, our study demonstrates that carefully adjusting the different filtering parameters applied to the GBS data is critical to maximize the selection of high quality genotypes and to increase the accuracy of imputation of missing data. The strategies and results presented here provide a framework to maximize the output of the GBS technique in animal populations and qualified the PstI/MspI GBS assay as a low-cost high-density genotyping platform. The conclusions reported here regarding read-depth and genotype quality filtering could benefit many GBS applications, notably genome-wide association studies, where there is a need to increase the density of markers genotyped across the target population while preserving the quality of genotypes.

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.

A 3' external transcribed spacer in a tRNA transcript acts as a sponge for small RNAs to prevent transcriptional noise.

During ribosomal and transfer RNA maturation, external transcribed spacer (ETS) and internal transcribed spacer (ITS) sequences are excised and, as non-functional by-products, are rapidly degraded. However, we report that the 3'ETS of the glyW-cysT-leuZ polycistronic tRNA precursor is highly and specifically enriched by co-purification with at least two different small regulatory RNAs (sRNAs), RyhB and RybB. Both sRNAs are shown to base pair with the same region in the 3'ETS of leuZ (3'ETS(leuZ)). Disrupting the pairing by mutating 3'ETS(leuZ) strongly increased the activity of sRNAs, even under non-inducing conditions. Our results indicate that 3'ETS(leuZ) prevents sRNA-dependent remodeling of tricarboxylic acid (TCA) cycle fluxes and decreases antibiotic sensitivity when sRNAs are transcriptionally repressed. This suggests that 3'ETS(leuZ) functions as a sponge to absorb transcriptional noise from repressed sRNAs. Additional data showing RybB and MicF sRNAs are co-purified with ITS(metZ-metW) and ITS(metW-metV) strongly suggest a wide distribution of this phenomenon.

The chloroplast genome of the green alga Schizomeris leibleinii (Chlorophyceae) provides evidence for bidirectional DNA replication from a single origin in the chaetophorales.

In the Chlorophyceae, the chloroplast genome is extraordinarily fluid in architecture and displays unique features relative to other groups of green algae. For the Chaetophorales, 1 of the 5 major lineages of the Chlorophyceae, it has been shown that the distinctive architecture of the 223,902-bp genome of Stigeoclonium helveticum is consistent with bidirectional DNA replication from a single origin. Here, we report the 182,759-bp chloroplast genome sequence of Schizomeris leibleinii, a member of the earliest diverging lineage of the Chaetophorales. Like its Stigeoclonium homolog, the Schizomeris genome lacks a large inverted repeat encoding the rRNA operon and displays a striking bias in coding regions that is associated with a bias in base composition along each strand. Our results support the notion that these two chaetophoralean genomes replicate bidirectionally from a putative origin located in the vicinity of the small subunit ribosomal RNA gene. Their shared structural characteristics were most probably inherited from the common ancestor of all chaetophoralean algae. Short dispersed repeats account for most of the 41-kb size variation between the Schizomeris and Stigeoclonium genomes, and there is no indication that homologous recombination between these repeated elements led to the observed gene rearrangements. A comparison of the extent of variation sustained by the Stigeoclonium and Schizomeris chloroplast DNAs (cpDNAs) with that observed for the cpDNAs of the chlamydomonadalean Chlamydomonas and Volvox suggests that gene rearrangements as well as changes in the abundance of intergenic and intron sequences occurred at a slower pace in the Chaetophorales than in the Chlamydomonadales.

The exceptionally large chloroplast genome of the green alga Floydiella terrestris illuminates the evolutionary history of the Chlorophyceae.

The Chlorophyceae, an advanced class of chlorophyte green algae, comprises five lineages that form two major clades (Chlamydomonadales + Sphaeropleales and Oedogoniales + Chaetopeltidales + Chaetophorales). The four complete chloroplast DNA (cpDNA) sequences currently available for chlorophyceans uncovered an extraordinarily fluid genome architecture as well as many structural features distinguishing this group from other green algae. We report here the 521,168-bp cpDNA sequence from a member of the Chaetopeltidales (Floydiella terrestris), the sole chlorophycean lineage not previously sampled for chloroplast genome analysis. This genome, which contains 97 conserved genes and 26 introns (19 group I and 7 group II introns), is the largest chloroplast genome ever sequenced. Intergenic regions account for 77.8% of the genome size and are populated by short repeats. Numerous genomic features are shared with the cpDNA of the chaetophoralean Stigeoclonium helveticum, notably the absence of a large inverted repeat and the presence of unique gene clusters and trans-spliced group II introns. Although only one of the Floydiella group I introns encodes a homing endonuclease gene, our finding of five free-standing reading frames having similarity with such genes suggests that chloroplast group I introns endowed with mobility were once more abundant in the Floydiella lineage. Parsimony analysis of structural genomic features and phylogenetic analysis of chloroplast sequence data unambiguously resolved the Oedogoniales as sister to the Chaetopeltidales and Chaetophorales. An evolutionary scenario of the molecular events that shaped the chloroplast genome in the Chlorophyceae is presented.

Chloroplast DNA sequence of the green alga Oedogonium cardiacum (Chlorophyceae): unique genome architecture, derived characters shared with the Chaetophorales and novel genes acquired through horizontal transfer.

BACKGROUND: To gain insight into the branching order of the five main lineages currently recognized in the green algal class Chlorophyceae and to expand our understanding of chloroplast genome evolution, we have undertaken the sequencing of chloroplast DNA (cpDNA) from representative taxa. The complete cpDNA sequences previously reported for Chlamydomonas (Chlamydomonadales), Scenedesmus (Sphaeropleales), and Stigeoclonium (Chaetophorales) revealed tremendous variability in their architecture, the retention of only few ancestral gene clusters, and derived clusters shared by Chlamydomonas and Scenedesmus. Unexpectedly, our recent phylogenies inferred from these cpDNAs and the partial sequences of three other chlorophycean cpDNAs disclosed two major clades, one uniting the Chlamydomonadales and Sphaeropleales (CS clade) and the other uniting the Oedogoniales, Chaetophorales and Chaetopeltidales (OCC clade). Although molecular signatures provided strong support for this dichotomy and for the branching of the Oedogoniales as the earliest-diverging lineage of the OCC clade, more data are required to validate these phylogenies. We describe here the complete cpDNA sequence of Oedogonium cardiacum (Oedogoniales). RESULTS: Like its three chlorophycean homologues, the 196,547-bp Oedogonium chloroplast genome displays a distinctive architecture. This genome is one of the most compact among photosynthetic chlorophytes. It has an atypical quadripartite structure, is intron-rich (17 group I and 4 group II introns), and displays 99 different conserved genes and four long open reading frames (ORFs), three of which are clustered in the spacious inverted repeat of 35,493 bp. Intriguingly, two of these ORFs (int and dpoB) revealed high similarities to genes not usually found in cpDNA. At the gene content and gene order levels, the Oedogonium genome most closely resembles its Stigeoclonium counterpart. Characters shared by these chlorophyceans but missing in members of the CS clade include the retention of psaM, rpl32 and trnL(caa), the loss of petA, the disruption of three ancestral clusters and the presence of five derived gene clusters. CONCLUSION: The Oedogonium chloroplast genome disclosed additional characters that bolster the evidence for a close alliance between the Oedogoniales and Chaetophorales. Our unprecedented finding of int and dpoB in this cpDNA provides a clear example that novel genes were acquired by the chloroplast genome through horizontal transfers, possibly from a mitochondrial genome donor.

DEEP DIVISION IN THE CHLOROPHYCEAE (CHLOROPHYTA) REVEALED BY CHLOROPLAST PHYLOGENOMIC ANALYSES(1).

The Chlorophyceae (sensu Mattox and Stewart) is a morphologically diverse class of the Chlorophyta displaying biflagellate and quadriflagellate motile cells with varying configurations of the flagellar apparatus. Phylogenetic analyses of 18S rDNA data and combined 18S and 26S rDNA data from a broad range of chlorophycean taxa uncovered five major monophyletic groups (Chlamydomonadales, Sphaeropleales, Oedogoniales, Chaetophorales, and Chaetopeltidales) but could not resolve their branching order. To gain insight into the interrelationships of these groups, we analyzed multiple genes encoded by the chloroplast genomes of Chlamydomonas reinhardtii P. A. Dang. and Chlamydomonas moewusii Gerloff (Chlamydomonadales), Scenedesmus obliquus (Turpin) Kütz. (Sphaeropleales), Oedogonium cardiacum Wittr. (Oedogoniales), Stigeoclonium helveticum Vischer (Chaetophorales), and Floydiella terrestris (Groover et Hofstetter) Friedl et O'Kelly (Chaetopeltidales). The C. moewusii, Oedogonium, and Floydiella chloroplast DNAs were partly sequenced using a random strategy. Trees were reconstructed from nucleotide and amino acid data sets derived from 44 protein-coding genes of 11 chlorophytes and nine streptophytes as well as from 57 protein-coding genes of the six chlorophycean taxa. All best trees identified two robustly supported major lineages within the Chlorophyceae: a clade uniting the Chlamydomonadales and Sphaeropleales, and a clade uniting the Oedogoniales, Chaetophorales, and Chaetopeltidales (OCC clade). This dichotomy is independently supported by molecular signatures in chloroplast genes, such as insertions/deletions and the distribution of trans-spliced group II introns. Within the OCC clade, the sister relationship observed for the Chaetophorales and Chaetopeltidales is also strengthened by independent data. Character state reconstruction of basal body orientation allowed us to refine hypotheses regarding the evolution of the flagellar apparatus.

Keratin 8 modulation of desmoplakin deposition at desmosomes in hepatocytes.

Keratins, the intermediate filament proteins of epithelial cells, connect to desmosomes, the cell-cell adhesion structures at the surface membrane. The building elements of desmosomes include desmoglein and desmocollin, which provide the actual cell adhesive properties, and desmoplakins, which anchor the keratin intermediate filaments to desmosomes. In the work reported here, we address the role of keratin 8 in modulating desmoplakin deposition at surface membrane in mouse hepatocytes. The experimental approach is based on the use of keratin 8- and keratin 18-null mouse hepatocytes as cell models. In wild-type mouse hepatocytes, desmoplakin is aligned with desmoglein and keratin 8 at the surface membrane. In keratin 8-null hepatocytes, the intermediate filament loss leads to alterations in desmoplakin distribution at the surface membrane, but not of desmoglein. Intriguingly, a significant proportion of keratin 18-null hepatocytes express keratin 8 at the surface membrane, associated with a proper desmoplakin alignment with desmoglein at desmosomes. A Triton treatment of the monolayer reveals that most of the desmoplakin present in either wild-type, keratin 8- or keratin 18-null hepatocytes is insoluble. Deletion analysis of keratin 8 further suggests that the recovery of desmoplakin alignment requires the keratin 8 rod domain. In addition, similarly to other works revealing a key role of desmoplakin phosphorylation on its interaction with intermediate filaments, we find that the phosphorylation status of the keratin 8 head domain affects desmoplakin distribution at desmosomes. Together, the data indicate that a proper alignment/deposition of desmoplakin with keratins and desmoglein in hepatocytes requires keratin 8, through a reciprocal phosphoserine-dependent process.

Distinctive architecture of the chloroplast genome in the chlorophycean green alga Stigeoclonium helveticum.

The chloroplast genome has experienced many architectural changes during the evolution of chlorophyte green algae, with the class Chlorophyceae displaying the lowest degree of ancestral traits. We have previously shown that the completely sequenced chloroplast DNAs (cpDNAs) of Chamydomonas reinhardtii (Chlamydomonadales) and Scenedesmus obliquus (Sphaeropleales) are highly scrambled in gene order relative to one another. Here, we report the complete cpDNA sequence of Stigeoclonium helveticum (Chaetophorales), a member of a third chlorophycean lineage. This genome, which encodes 97 genes and contains 21 introns (including four putatively trans-spliced group II introns inserted at novel sites), is remarkably rich in derived features and extremely rearranged relative to its chlorophycean counterparts. At 223,902 bp, Stigeoclonium cpDNA is the largest chloroplast genome sequenced thus far, and in contrast to those of Chlamydomonas and Scenedesmus, features no large inverted repeat. Interestingly, the pattern of gene distribution between the DNA strands and the bias in base composition along each strand suggest that the Stigeoclonium genome replicates bidirectionally from a single origin. Unlike most known trans-spliced group II introns, those of Stigeoclonium exhibit breaks in domains I and II. By placing our comparative genome analyses in a phylogenetic framework, we inferred an evolutionary scenario of the mutational events that led to changes in genome architecture in the Chlorophyceae.