Global mesozooplankton communities show lower connectivity in deep oceanic layers

Mesozooplankton is a key component of the ocean, regulating global processes such as the carbon pump, and ensuring energy transfer from lower to higher trophic levels. Yet, knowledge on mesozooplankton diversity, distribution and connectivity at global scale is still fragmented. To fill this gap, we applied DNA metabarcoding to mesozooplankton samples collected during the Malaspina-2010 circumnavigation expedition across the Atlantic, Indian and Pacific oceans from the surface to bathypelagic depths. We highlight the still scarce knowledge on global mesozooplankton diversity and identify the Indian Ocean and the deep sea as the oceanic regions with the highest proportion of hidden diversity. We report no consistent alpha-diversity patterns for mesozooplankton at a global scale, neither across vertical nor horizontal gradients. However, beta-diversity analysis suggests horizontal and vertical structuring of mesozooplankton communities mostly attributed to turnover and reveals an increase in mesozooplankton beta-diversity with depth, indicating reduced connectivity at deeper layers. Additionally, we identify a water mass type-mediated structuring of mesozooplankton bathypelagic communities instead of an oceanic basin-mediated as observed at upper layers. This suggests limited dispersal at deep ocean layers, most likely due to weaker currents and lower mixing of water mass types, thus reinforcing the importance of oceanic currents and barriers to dispersal in shaping global plankton communities.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Otolith chemical fingerprints of skipjack tuna (Katsuwonus pelamis) in the Indian Ocean: First insights into stock structure delineation

The chemical composition of otoliths (earbones) can provide valuable information about stock structure and connectivity patterns among marine fish. For that, chemical signatures must be sufficiently distinct to allow accurate classification of an unknown fish to their area of origin. Here we have examined the suitability of otolith microchemistry as a tool to better understand the spatial dynamics of skipjack tuna (Katsuwonus pelamis), a highly valuable commercial species for which uncertainties remain regarding its stock structure in the Indian Ocean. For this aim, we have compared the early life otolith chemical composition of young-of-the-year (<6 months) skipjack tuna captured from the three main nursery areas of the equatorial Indian Ocean (West, Central and East). Elemental (Li:Ca, Sr:Ca, Ba:Ca, Mg:Ca and Mn:Ca) and stable isotopic (δ13C, δ18O) signatures were used, from individuals captured in 2018 and 2019. Otolith Sr:Ca, Ba:Ca, Mg:Ca and δ18O significantly differed among fish from different nurseries, but, in general, the chemical signatures of the three nursery areas largely overlapped. Multivariate analyses of otolith chemical signatures revealed low geographic separation among Central and Eastern nurseries, achieving a maximum overall random forest cross validated classification success of 51%. Cohort effect on otolith trace element signatures was also detected, indicating that variations in chemical signatures associated with seasonal changes in oceanographic conditions must be well understood, particularly for species with several reproductive peaks throughout the year. Otolith microchemistry in conjunction with other techniques (e.g., genetics, particle tracking) should be further investigated to resolve skipjack stock structure, which will ultimately contribute to the sustainable management of this stock in the Indian Ocean.

Pan-regional marine benthic cryptobiome biodiversity patterns revealed by metabarcoding Autonomous Reef Monitoring Structures

Autonomous Reef Monitoring Structures (ARMS) have been applied worldwide to characterize the critical yet frequently overlooked biodiversity patterns of marine benthic organisms. In order to disentangle the relevance of environmental factors in benthic patterns, here, through standardized metabarcoding protocols, we analyse sessile and mobile (<2 mm) organisms collected using ARMS deployed across six regions with different environmental conditions (3 sites × 3 replicates per region): Baltic, Western Mediterranean, Adriatic, Black and Red Seas, and the Bay of Biscay. A total of 27,473 Amplicon Sequence Variants (ASVs) were observed ranging from 1,404 in the Black Sea to 9,958 in the Red Sea. No ASVs were shared among all regions. The highest number of shared ASVs was between the Western Mediterranean and the Adriatic Sea (116) and Bay of Biscay (115). Relatively high numbers of ASVs (103), mostly associated with the genus Amphibalanus, were also shared between the lower salinity seas (Baltic and Black Seas). We found that compositional differences in spatial patterns of rocky-shore benthos are determined slightly more by dispersal limitation than environmental filtering. Dispersal limitation was similar between sessile and mobile groups, while the sessile group had a larger environmental niche breadth than the mobile group. Further, our study can provide a foundation for future evaluations of biodiversity patterns in the cryptobiome, which can contribute up to 70% of the local biodiversity.

Combining genetic markers with stable isotopes in otoliths reveals complexity in the stock structure of Atlantic bluefin tuna (Thunnus thynnus)

Atlantic bluefin tuna (Thunnus thynnus) from the two main spawning populations in the Mediterranean and Gulf of Mexico occur together in the western, central and eastern Atlantic. Stock composition of catches from mixing areas is uncertain, presenting a major challenge to the sustainable management of the fisheries. This study combines genetic and chemical markers to develop an integrated method of population assignment. Stable isotope signatures (δ¹³C and δ¹⁸O) in the otolith core of adults from the two main spawning populations (adult baselines) showed less overlap than those of yearlings (12-18 months old) from western and eastern nursery areas suggesting that some exchange occurs towards the end of the yearling phase. The integrated model combined δ¹⁸O with four genetic markers (SNPs) to distinguish the adult baselines with greater accuracy than chemical or genetic markers alone. When used to assign individuals from the mixing areas to their population of origin, the integrated model resolved some (but not all) discrepancies between the chemistry and genetic methods. Some individuals in the mixing area had otolith δ¹⁸O values and genetic profiles which when taken together, were not representative of either population. These fish may originate from another Atlantic spawning area or may represent population contingents that move away from the main spawning areas during the first year of life. This complexity in stock structure is not captured by the current two-stock model.

Environmental DNA Metabarcoding: A Promising Tool for Ballast Water Monitoring

Nonindigenous species are introduced worldwide with ballast water (BW). To prevent further introductions, oceanic BW exchange and BW treatment systems are utilized, but their performance needs to be evaluated. To that aim, characterizing BW communities is essential but usually relies on exhaustive sampling and morphological taxonomic identification, which does not always allow fine-scale taxonomic resolution. Through the analysis of BW samples from 11 vessels arriving to the Chesapeake Bay (USA), we evaluated the potential of environmental DNA (eDNA) metabarcoding for BW monitoring by assessing whether the impact of BW management type could be identified, analyzing the influence of BW sampling access locations on communities, and comparing the accuracy of eDNA for taxonomic assignment and identification of nonindigenous taxa. We found that (1) different sampling access locations of the same tank resulted in different communities, (2) communities from treated and exchanged BW differ, (3) signals of source port and of ocean exchange are observed, (4) eDNA metabarcoding results in more diversity than morphological taxonomy, and (5) the nonindigenous copepod Oithona davisae, not reported before in the Chesapeake Bay, is detected. Overall, this study highlights the potential of eDNA metabarcoding for BW monitoring, but more comprehensive sampling will be needed to optimize the approach.

Selecting RAD-Seq Data Analysis Parameters for Population Genetics: The More the Better?

Restriction site-associated DNA sequencing (RAD-seq) has become a powerful and widely used tool in molecular ecology studies as it allows to cost-effectively recover thousands of polymorphic sites across individuals of non-model organisms. However, its successful implementation in population genetics relies on correct data processing that would minimize potential loci-assembly biases and consequent genotyping error rates. RAD-seq data processing when no reference genome is available involves the assembly of hundreds of thousands high-throughput sequencing reads into orthologous loci, for which various key parameter values need to be selected by the researcher. Previous studies exploring the effect of these parameter values found or assumed that a larger number of recovered polymorphic loci is associated with a better assembly. Here, using three RAD-seq datasets from different species, we explore the effect of read filtering, loci assembly and polymorphic site selection on number of markers obtained and genetic differentiation inferred using the Stacks software. We find (i) that recovery of higher numbers of polymorphic loci is not necessarily associated with higher genetic differentiation, (ii) that the presence of PCR duplicates, selected loci assembly parameters and selected SNP filtering parameters affect the number of recovered polymorphic loci and degree of genetic differentiation, and (iii) that this effect is different in each dataset, meaning that defining a systematic universal protocol for RAD-seq data analysis may lead to missing relevant information about population differentiation.

A bacterial community-based index to assess the ecological status of estuarine and coastal environments

Biotic indices for monitoring marine ecosystems are mostly based on the analysis of benthic macroinvertebrate communities. Due to their high sensitivity to pollution and fast response to environmental changes, bacterial assemblages could complement the information provided by benthic metazoan communities as indicators of human-induced impacts, but so far, this biological component has not been well explored for this purpose. Here we performed 16S rRNA gene amplicon sequencing to analyze the bacterial assemblage composition of 51 estuarine and coastal stations characterized by different environmental conditions and human-derived pressures. Using the relative abundance of putative indicator bacterial taxa, we developed a biotic index that is significantly correlated with a sediment quality index calculated on the basis of organic and inorganic compound concentrations. This new index based on bacterial assemblage composition can be a sensitive tool for providing a fast environmental assessment and allow a more comprehensive integrative ecosystem approach for environmental management.

Environmental status assessment using DNA metabarcoding: towards a genetics based Marine Biotic Index (gAMBI)

Marine ecosystem protection and conservation initiatives rely on the assessment of ecological integrity and health status of marine environments. The AZTI's Marine Biotic Index (AMBI), which consists on using macroinvertebrate diversity as indicator of ecosystem health, is used worldwide for this purpose. Yet, this index requires taxonomic assignment of specimens, which typically involves a time and resource consuming visual identification of each sample. DNA barcoding or metabarcoding are potential harmonized, faster and cheaper alternatives for species identification, although the suitability of these methods for easing the implementation of the AMBI is yet to be evaluated. Here, we analyze the requirements for the implementation of a genetics based AMBI (gAMBI), and show, using available sequence data, that information about presence/absence of the most frequently occurring species provides accurate AMBI values. Our results set the basics for the implementation of the gAMBI, which has direct implications for a faster and cheaper marine monitoring and health status assessment.

Genomics in marine monitoring: new opportunities for assessing marine health status

This viewpoint paper explores the potential of genomics technology to provide accurate, rapid, and cost efficient observations of the marine environment. The use of such approaches in next generation marine monitoring programs will help achieve the goals of marine legislation implemented world-wide. Genomic methods can yield faster results from monitoring, easier and more reliable taxonomic identification, as well as quicker and better assessment of the environmental status of marine waters. A summary of genomic methods that are ready or show high potential for integration into existing monitoring programs is provided (e.g. qPCR, SNP based methods, DNA barcoding, microarrays, metagenetics, metagenomics, transcriptomics). These approaches are mapped to existing indicators and descriptors and a series of case studies is presented to assess the cost and added value of these molecular techniques in comparison with traditional monitoring systems. Finally, guidelines and recommendations are suggested for how such methods can enter marine monitoring programs in a standardized manner.

Whole transcriptome analysis of Acinetobacter baumannii assessed by RNA-sequencing reveals different mRNA expression profiles in biofilm compared to planktonic cells

Acinetobacterbaumannii has emerged as a dangerous opportunistic pathogen, with many strains able to form biofilms and thus cause persistent infections. The aim of the present study was to use high-throughput sequencing techniques to establish complete transcriptome profiles of planktonic (free-living) and sessile (biofilm) forms of A. baumannii ATCC 17978 and thereby identify differences in their gene expression patterns. Collections of mRNA from planktonic (both exponential and stationary phase cultures) and sessile (biofilm) cells were sequenced. Six mRNA libraries were prepared following the mRNA-Seq protocols from Illumina. Reads were obtained in a HiScanSQ platform and mapped against the complete genome to describe the complete mRNA transcriptomes of planktonic and sessile cells. The results showed that the gene expression pattern of A. baumannii biofilm cells was distinct from that of planktonic cells, including 1621 genes over-expressed in biofilms relative to stationary phase cells and 55 genes expressed only in biofilms. These differences suggested important changes in amino acid and fatty acid metabolism, motility, active transport, DNA-methylation, iron acquisition, transcriptional regulation, and quorum sensing, among other processes. Disruption or deletion of five of these genes caused a significant decrease in biofilm formation ability in the corresponding mutant strains. Among the genes over-expressed in biofilm cells were those in an operon involved in quorum sensing. One of them, encoding an acyl carrier protein, was shown to be involved in biofilm formation as demonstrated by the significant decrease in biofilm formation by the corresponding knockout strain. The present work serves as a basis for future studies examining the complex network systems that regulate bacterial biofilm formation and maintenance.

Fine mapping of a major histocompatibility complex in ankylosing spondylitis: association of the HLA-DPA1 and HLA-DPB1 regions

OBJECTIVE

To investigate the potential association of major histocompatibility complex (MHC) markers other than HLA-B27 with ankylosing spondylitis (AS).

METHODS

A total of 603 patients with AS and 542 healthy control subjects, all of whom were HLA-B27 positive, were selected for this study based on clinical criteria. First, high-density genotyping across the MHC region (2,360 single-nucleotide polymorphisms [SNPs]) was performed in a cohort of 191 patients and 241 control subjects. After a fine-mapping study, 5 SNPs from the HLA-DPA1/DPB1 region were validated in a second cohort of 412 patients with AS and 301 healthy control subjects.

RESULTS

Seventeen SNPs located within or near the HLA-DPA1 and HLA-DPB1 loci showed association with AS (P = 1.38 × 10⁻⁵ to 0.05). In addition, multimarker tests, both linkage disequilibrium and sliding windows, showed association of some groups of adjacent SNPs within the HLA-DPA1/DPB1 region with AS (P = 1.0 × 10⁻⁴ to 3.96 × 10⁻⁷). We validated the association by genotyping 5 SNPs from the DPA1/DPB1 region in an additional cohort and obtained P values from 6.42 × 10⁻⁵ to 0.01 in the analysis of the combined cohorts. Subtyping analysis of HLA-DPA1 and HLA-DPB1 showed that HLA-DPA1*01:03, A1*02:01, and B1*13:01 were the subtypes most susceptible to AS.

CONCLUSION

HLA markers and linkage disequilibrium blocks near HLA-DPA1 and HLA-DPB1 are statistically associated with AS. We identified a region located around the HLA-DPA1 and HLA-DPB1 loci associated with AS, another region within the MHC that is different from HLA-B27.

Whole transcriptome analysis of a reversible neurodegenerative process in Drosophila reveals potential neuroprotective genes

Background

Neurodegenerative diseases are progressive and irreversible and they can be initiated by mutations in specific genes. Spalt-like genes (Sall) encode transcription factors expressed in the central nervous system. In humans, SALL mutations are associated with hereditary syndromes characterized by mental retardation, sensorineural deafness and motoneuron problems, among others. Drosophila sall mutants exhibit severe neurodegeneration of the central nervous system at embryonic stage 16, which surprisingly reverts later in development at embryonic stage 17, suggesting a potential to recover from neurodegeneration. We hypothesize that this recovery is mediated by a reorganization of the transcriptome counteracting SALL lost. To identify genes associated to neurodegeneration and neuroprotection, we used mRNA-Seq to compare the transcriptome of Drosophila sall mutant and wild type embryos from neurodegeneration and reversal stages.

Results

Neurodegeneration stage is associated with transcriptional changes in 220 genes, of which only 5% were already described as relevant for neurodegeneration. Genes related to the groups of Redox, Lifespan/Aging and Mitochondrial diseases are significantly represented at this stage. By contrast, neurodegeneration reversal stage is associated with significant changes in 480 genes, including 424 not previously associated with neuroprotection. Immune response and Salt stress are the most represented groups at this stage.

Conclusions

We identify new genes associated to neurodegeneration and neuroprotection by using an mRNA-Seq approach. The strong homology between Drosophila and human genes raises the possibility to unveil novel genes involved in neurodegeneration and neuroprotection also in humans.

Comparison of methods to detect copy number alterations in cancer using simulated and real genotyping data

Background

The detection of genomic copy number alterations (CNA) in cancer based on SNP arrays requires methods that take into account tumour specific factors such as normal cell contamination and tumour heterogeneity. A number of tools have been recently developed but their performance needs yet to be thoroughly assessed. To this aim, a comprehensive model that integrates the factors of normal cell contamination and intra-tumour heterogeneity and that can be translated to synthetic data on which to perform benchmarks is indispensable.

Results

We propose such model and implement it in an R package called CnaGen to synthetically generate a wide range of alterations under different normal cell contamination levels. Six recently published methods for CNA and loss of heterozygosity (LOH) detection on tumour samples were assessed on this synthetic data and on a dilution series of a breast cancer cell-line: ASCAT, GAP, GenoCNA, GPHMM, MixHMM and OncoSNP. We report the recall rates in terms of normal cell contamination levels and alteration characteristics: length, copy number and LOH state, as well as the false discovery rate distribution for each copy number under different normal cell contamination levels.

Assessed methods are in general better at detecting alterations with low copy number and under a little normal cell contamination levels. All methods except GPHMM, which failed to recognize the alteration pattern in the cell-line samples, provided similar results for the synthetic and cell-line sample sets. MixHMM and GenoCNA are the poorliest performing methods, while GAP generally performed better. This supports the viability of approaches other than the common hidden Markov model (HMM)-based.

Conclusions

We devised and implemented a comprehensive model to generate data that simulate tumoural samples genotyped using SNP arrays. The validity of the model is supported by the similarity of the results obtained with synthetic and real data. Based on these results and on the software implementation of the methods, we recommend GAP for advanced users and GPHMM for a fully driven analysis.

A high density SNP genotyping approach within the 19q13 chromosome region identifies an association of a CNOT3 polymorphism with ankylosing spondylitis

OBJECTIVE

To identify genomic variants in the 19q13 chromosome region associated with ankylosing spondylitis (AS) in human leucocyte antigen (HLA)-B27-positive populations.

METHODS

High-throughput genotyping of 1536 haplotype-tag single nucleotide polymorphisms (SNPs) was performed in 249 patients with AS and 302 healthy controls. Some of the identified associations were validated by genotyping four SNPs in two additional cohorts consisting of 412 cases/301 controls and 144 cases/203 controls. All individuals selected (both cases and controls) were HLA-B27-positive.

RESULTS

Two markers in two different genes (CNOT3 and LAIR2) showed significant association (p<10(-3)) with AS. In addition, sliding windows analysis showed association of groups of adjacent SNPs in regions located around CNOT3 (Chr19: 59347459-59356564, p=2.43 × 10(-4) to 6.54 × 10(-4)). The associations were validated by genotyping four SNPs from regions located near LAIR2 and CNOT3 genes (rs1055234, rs8111398, rs2287828 and rs4591276) in two additional cohorts. The CNOT3 polymorphism (rs1055234) remained associated with AS (combined p=9.73 × 10(-6)). One SNP, located downstream of KIR3DL1, was detected which, tested in combination with HLA-Bw4I80, was associated with AS.

CONCLUSION

A novel significant association was detected between SNP rs1055234 and AS susceptibility.

miRanalyzer: an update on the detection and analysis of microRNAs in high-throughput sequencing experiments

We present a new version of miRanalyzer, a web server and stand-alone tool for the detection of known and prediction of new microRNAs in high-throughput sequencing experiments. The new version has been notably improved regarding speed, scope and available features. Alignments are now based on the ultrafast short-read aligner Bowtie (granting also colour space support, allowing mismatches and improving speed) and 31 genomes, including 6 plant genomes, can now be analysed (previous version contained only 7). Differences between plant and animal microRNAs have been taken into account for the prediction models and differential expression of both, known and predicted microRNAs, between two conditions can be calculated. Additionally, consensus sequences of predicted mature and precursor microRNAs can be obtained from multiple samples, which increases the reliability of the predicted microRNAs. Finally, a stand-alone version of the miRanalyzer that is based on a local and easily customized database is also available; this allows the user to have more control on certain parameters as well as to use specific data such as unpublished assemblies or other libraries that are not available in the web server. miRanalyzer is available at http://bioinfo2.ugr.es/miRanalyzer/miRanalyzer.php.

Phylogenomic evidence for separate acquisition of plastids in cryptophytes, haptophytes, and stramenopiles

According to the chromalveolate hypothesis (Cavalier-Smith T. 1999. Principles of protein and lipid targeting in secondary symbiogenesis: euglenoid, dinoflagellate, and sporozoan plastid origins and the eukaryote family tree. J Eukaryot Microbiol 46:347-366), the four eukaryotic groups with chlorophyll c-containing plastids originate from a single photosynthetic ancestor, which acquired its plastids by secondary endosymbiosis with a red alga. So far, molecular phylogenies have failed to either support or disprove this view. Here, we devise a phylogenomic falsification of the chromalveolate hypothesis that estimates signal strength across the three genomic compartments: If the four chlorophyll c-containing lineages indeed derive from a single photosynthetic ancestor, then similar amounts of plastid, mitochondrial, and nuclear sequences should allow to recover their monophyly. Our results refute this prediction, with statistical support levels too different to be explained by evolutionary rate variation, phylogenetic artifacts, or endosymbiotic gene transfer. Therefore, we reject the chromalveolate hypothesis as falsified in favor of more complex evolutionary scenarios involving multiple higher order eukaryote-eukaryote endosymbioses.

Detecting and overcoming systematic errors in genome-scale phylogenies

Genome-scale data sets result in an enhanced resolution of the phylogenetic inference by reducing stochastic errors. However, there is also an increase of systematic errors due to model violations, which can lead to erroneous phylogenies. Here, we explore the impact of systematic errors on the resolution of the eukaryotic phylogeny using a data set of 143 nuclear-encoded proteins from 37 species. The initial observation was that, despite the impressive amount of data, some branches had no significant statistical support. To demonstrate that this lack of resolution is due to a mutual annihilation of phylogenetic and nonphylogenetic signals, we created a series of data sets with slightly different taxon sampling. As expected, these data sets yielded strongly supported but mutually exclusive trees, thus confirming the presence of conflicting phylogenetic and nonphylogenetic signals in the original data set. To decide on the correct tree, we applied several methods expected to reduce the impact of some kinds of systematic error. Briefly, we show that (i) removing fast-evolving positions, (ii) recoding amino acids into functional categories, and (iii) using a site-heterogeneous mixture model (CAT) are three effective means of increasing the ratio of phylogenetic to nonphylogenetic signal. Finally, our results allow us to formulate guidelines for detecting and overcoming phylogenetic artefacts in genome-scale phylogenetic analyses.

Toward Resolving the Eukaryotic Tree: The Phylogenetic Positions of Jakobids and Cercozoans

Resolving the global phylogeny of eukaryotes has proven to be challenging. Among the eukaryotic groups of uncertain phylogenetic position are jakobids, a group of bacterivorous flagellates that possess the most bacteria-like mitochondrial genomes known. Jakobids share several ultrastructural features with malawimonads and an assemblage of anaerobic protists (e.g., diplomonads and oxymonads). These lineages together with Euglenozoa and Heterolobosea have collectively been designated "excavates". However, published molecular phylogenies based on the sequences of nuclear rRNAs and up to six nucleus-encoded proteins do not provide convincing support for the monophyly of excavates, nor do they uncover their relationship to other major eukaryotic groups. Here, we report the first large-scale eukaryotic phylogeny, inferred from 143 nucleus-encoded proteins comprising 31,604 amino acid positions, that includes jakobids, malawimonads and cercozoans. We obtain compelling support for the monophyly of jakobids, Euglenozoa plus Heterolobosea (JEH group), and for the association of cercozoans with stramenopiles plus alveolates. Furthermore, we observe a sister-group relationship between the JEH group and malawimonads after removing fast-evolving species from the dataset. We discuss the implications of these results for the concept of "excavates" and for the elucidation of eukaryotic phylogeny in general.

Phylogenetic analyses of nuclear, mitochondrial, and plastid multigene data sets support the placement of Mesostigma in the Streptophyta

All extant green plants belong to 1 of 2 major lineages, commonly known as the Chlorophyta (most of the green algae) and the Streptophyta (land plants and their closest green algal relatives). The scaly green flagellate Mesostigma viride has an important place in the debate on the origin of green plants. However, there have been conflicting results from molecular systematics as to whether Mesostigma diverges before the Chlorophyta/Streptophyta split or is an early diverging flagellate member of the Streptophyta. Previous studies employed either a limited taxon sampling (plastid and mitochondrial genomes) or a small number of phylogenetically informative sites (single nuclear genes). Here, we use large data sets from the nuclear (125 proteins; 29,319 positions), mitochondrial (33 proteins; 6,622 positions), and plastid (50 proteins; 10,137 positions) genomes with an expanded taxon sampling (21, 13, and 28 species, respectively) to reevaluate the phylogenetic position of Mesostigma. Our study supports the placement of Mesostigma in the Streptophyta (as an early diverging lineage) and provides evidence that systematic biases have played a role in generating some of the previous conflicting results. Importantly, we demonstrate that using an increased taxon sampling as well as more realistic models of evolution allows increasing congruence among the nuclear, mitochondrial, and plastid data sets.

Monophyly of Primary Photosynthetic Eukaryotes: Green Plants, Red Algae, and Glaucophytes

Between 1 and 1.5 billion years ago, eukaryotic organisms acquired the ability to convert light into chemical energy through endosymbiosis with a Cyanobacterium (e.g.,). This event gave rise to "primary" plastids, which are present in green plants, red algae, and glaucophytes ("Plantae" sensu Cavalier-Smith). The widely accepted view that primary plastids arose only once implies two predictions: (1) all plastids form a monophyletic group, as do (2) primary photosynthetic eukaryotes. Nonetheless, unequivocal support for both predictions is lacking (e.g.,). In this report, we present two phylogenomic analyses, with 50 genes from 16 plastid and 15 cyanobacterial genomes and with 143 nuclear genes from 34 eukaryotic species, respectively. The nuclear dataset includes new sequences from glaucophytes, the less-studied group of primary photosynthetic eukaryotes. We find significant support for both predictions. Taken together, our analyses provide the first strong support for a single endosymbiotic event that gave rise to primary photosynthetic eukaryotes, the Plantae. Because our dataset does not cover the entire eukaryotic diversity (but only four of six major groups in), further testing of the monophyly of Plantae should include representatives from eukaryotic lineages for which currently insufficient sequence information is available.