Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean.

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans.

Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.

New insights into global biogeography, population structure and natural selection from the genome of the epipelagic copepod Oithona.

In the epipelagic ocean, the genus Oithona is considered as one of the most abundant and widespread copepods and plays an important role in the trophic food web. Despite its ecological importance, little is known about Oithona and cyclopoid copepods genomics. Therefore, we sequenced, assembled and annotated the genome of Oithona nana. The comparative genomic analysis integrating available copepod genomes highlighted the expansions of genes related to stress response, cell differentiation and development, including genes coding Lin12-Notch-repeat (LNR) domain proteins. The Oithona biogeography based on 28S sequences and metagenomic reads from the Tara Oceans expedition showed the presence of O. nana mostly in the Mediterranean Sea (MS) and confirmed the amphitropical distribution of Oithona similis. The population genomics analyses of O. nana in the Northern MS, integrating the Tara Oceans metagenomic data and the O. nana genome, led to the identification of genetic structure between populations from the MS basins. Furthermore, 20 loci were found to be under positive selection including four missense and eight synonymous variants, harbouring soft or hard selective sweep patterns. One of the missense variants was localized in the LNR domain of the coding region of a male-specific gene. The variation in the B-allele frequency with respect to the MS circulation pattern showed the presence of genomic clines between O. nana and another undefined Oithona species possibly imported through Atlantic waters. This study provides new approaches and results in zooplankton population genomics through the integration of metagenomic and oceanographic data.

Loss of function of glucocerebrosidase GBA2 is responsible for motor neuron defects in hereditary spastic paraplegia.

Spastic paraplegia 46 refers to a locus mapped to chromosome 9 that accounts for a complicated autosomal-recessive form of hereditary spastic paraplegia (HSP). With next-generation sequencing in three independent families, we identified four different mutations in GBA2 (three truncating variants and one missense variant), which were found to cosegregate with the disease and were absent in controls. GBA2 encodes a microsomal nonlysosomal glucosylceramidase that catalyzes the conversion of glucosylceramide to free glucose and ceramide and the hydrolysis of bile acid 3-O-glucosides. The missense variant was also found at the homozygous state in a simplex subject in whom no residual glucocerebrosidase activity of GBA2 could be evidenced in blood cells, opening the way to a possible measurement of this enzyme activity in clinical practice. The overall phenotype was a complex HSP with mental impairment, cataract, and hypogonadism in males associated with various degrees of corpus callosum and cerebellar atrophy on brain imaging. Antisense morpholino oligonucleotides targeting the zebrafish GBA2 orthologous gene led to abnormal motor behavior and axonal shortening/branching of motoneurons that were rescued by the human wild-type mRNA but not by applying the same mRNA containing the missense mutation. This study highlights the role of ceramide metabolism in HSP pathology.

High-throughput sequencing of a 4.1 Mb linkage interval reveals FLVCR2 deletions and mutations in lethal cerebral vasculopathy.

Rare lethal disease gene identification remains a challenging issue, but it is amenable to new techniques in high-throughput sequencing (HTS). Cerebral proliferative glomeruloid vasculopathy (PGV), or Fowler syndrome, is a severe autosomal recessive disorder of brain angiogenesis, resulting in abnormally thickened and aberrant perforating vessels leading to hydranencephaly. In three multiplex consanguineous families, genome-wide SNP analysis identified a locus of 14 Mb on chromosome 14. In addition, 280 consecutive SNPs were identical in two Turkish families unknown to be related, suggesting a founder mutation reducing the interval to 4.1 Mb. To identify the causative gene, we then specifically enriched for this region with sequence capture and performed HTS in a proband of seven families. Due to technical constraints related to the disease, the average coverage was only 7×. Nonetheless, iterative bioinformatic analyses of the sequence data identified mutations and a large deletion in the FLVCR2 gene, encoding a 12 transmembrane domain-containing putative transporter. A striking absence of alpha-smooth muscle actin immunostaining in abnormal vessels in fetal PGV brains, suggests a deficit in pericytes, cells essential for capillary stabilization and remodeling during brain angiogenesis. This is the first lethal disease-causing gene to be identified by comprehensive HTS of an entire linkage interval.