The exception that proves the rule: Virulence gene expression at the onset of Plasmodium falciparum blood stage infections.

Controlled human malaria infections (CHMI) are a valuable tool to study parasite gene expression in vivo under defined conditions. In previous studies, virulence gene expression was analyzed in samples from volunteers infected with the Plasmodium falciparum (Pf) NF54 isolate, which is of African origin. Here, we provide an in-depth investigation of parasite virulence gene expression in malaria-naïve European volunteers undergoing CHMI with the genetically distinct Pf 7G8 clone, originating in Brazil. Differential expression of var genes, encoding major virulence factors of Pf, PfEMP1s, was assessed in ex vivo parasite samples as well as in parasites from the in vitro cell bank culture that was used to generate the sporozoites (SPZ) for CHMI (Sanaria PfSPZ Challenge (7G8)). We report broad activation of mainly B-type subtelomeric located var genes at the onset of a 7G8 blood stage infection in naïve volunteers, mirroring the NF54 expression study and suggesting that the expression of virulence-associated genes is generally reset during transmission from the mosquito to the human host. However, in 7G8 parasites, we additionally detected a continuously expressed single C-type variant, Pf7G8_040025600, that was most highly expressed in both pre-mosquito cell bank and volunteer samples, suggesting that 7G8, unlike NF54, maintains expression of some previously expressed var variants during transmission. This suggests that in a new host, the parasite may preferentially express the variants that previously allowed successful infection and transmission. Trial registration: ClinicalTrials.gov - NCT02704533; 2018-004523-36.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients.

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104 cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104 cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104 cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO2. Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean.

Identification of novel inner membrane complex and apical annuli proteins of the malaria parasite Plasmodium falciparum.

The inner membrane complex (IMC) is a defining feature of apicomplexan parasites, which confers stability and shape to the cell, functions as a scaffolding compartment during the formation of daughter cells and plays an important role in motility and invasion during different life cycle stages of these single-celled organisms. To explore the IMC proteome of the malaria parasite Plasmodium falciparum we applied a proximity-dependent biotin identification (BioID)-based proteomics approach, using the established IMC marker protein Photosensitized INA-Labelled protein 1 (PhIL1) as bait in asexual blood-stage parasites. Subsequent mass spectrometry-based peptide identification revealed enrichment of 12 known IMC proteins and several uncharacterized candidate proteins. We validated nine of these previously uncharacterized proteins by endogenous GFP-tagging. Six of these represent new IMC proteins, while three proteins have a distinct apical localization that most likely represents structures described as apical annuli in Toxoplasma gondii. Additionally, various Kelch13 interacting candidates were identified, suggesting an association of the Kelch13 compartment and the IMC in schizont and merozoite stages. This work extends the number of validated IMC proteins in the malaria parasite and reveals for the first time the existence of apical annuli proteins in P. falciparum. Additionally, it provides evidence for a spatial association between the Kelch13 compartment and the IMC in late blood-stage parasites.

Structure of Prototypic Peptide Transporter DtpA from E. coli in Complex with Valganciclovir Provides Insights into Drug Binding of Human PepT1.

Members of the solute carrier 15 family (SLC15) transport di- and tripeptides as well as peptidomimetic drugs across the cell membrane. Structures of bacterial homologues have provided valuable information on the binding and transport of their natural substrates, but many do not transport medically relevant drugs. In contrast, a homologue from Escherichia coli, DtpA (dipeptide and tripeptide permease), shows a high similarity to human PepT1 (SLC15A1) in terms of ligand selectivity and transports a similar set of drugs. Here, we present the crystal structure of DtpA in ligand-free form (at 3.30 Å resolution) and in complex with the antiviral prodrug valganciclovir (at 2.65 Å resolution) supported by biochemical data. We show that valganciclovir unexpectedly binds with the ganciclovir moiety mimicking the N-terminal residue of a canonical peptide substrate. On the basis of a homology model we argue that this binding mode also applies to the human PepT1 transporter. Our results provide new insights into the binding mode of prodrugs and will assist the rational design of drugs with improved absorption rates.

Recurring seasonality exposes dominant species and niche partitioning strategies of open ocean picoeukaryotic algae.

Ocean spring phytoplankton blooms are dynamic periods important to global primary production. We document vertical patterns of a diverse suite of eukaryotic algae, the prasinophytes, in the North Atlantic Subtropical Gyre with monthly sampling over four years at the Bermuda Atlantic Time-series Study site. Water column structure was used to delineate seasonal stability periods more ecologically relevant than seasons defined by calendar dates. During winter mixing, tiny prasinophytes dominated by Class II comprise 46 ± 24% of eukaryotic algal (plastid-derived) 16S rRNA V1-V2 amplicons, specifically Ostreococcus Clade OII, Micromonas commoda, and Bathycoccus calidus. In contrast, Class VII are rare and Classes I and VI peak during warm stratified periods when surface eukaryotic phytoplankton abundances are low. Seasonality underpins a reservoir of genetic diversity from multiple prasinophyte classes during warm periods that harbor ephemeral taxa. Persistent Class II sub-species dominating the winter/spring bloom period retreat to the deep chlorophyll maximum in summer, poised to seed the mixed layer upon winter convection, exposing a mechanism for initiating high abundances at bloom onset. Comparisons to tropical oceans reveal broad distributions of the dominant sub-species herein. This unparalleled window into temporal and spatial niche partitioning of picoeukaryotic primary producers demonstrates how key prasinophytes prevail in warm oceans.

Plasticity of the binding pocket in peptide transporters underpins promiscuous substrate recognition.

Proton-dependent oligopeptide transporters (POTs) are promiscuous transporters of the major facilitator superfamily that constitute the main route of entry for a wide range of dietary peptides and orally administrated peptidomimetic drugs. Given their clinical and pathophysiological relevance, several POT homologs have been studied extensively at the structural and molecular level. However, the molecular basis of recognition and transport of diverse peptide substrates has remained elusive. We present 14 X-ray structures of the bacterial POT DtpB in complex with chemically diverse di- and tripeptides, providing novel insights into the plasticity of the conserved central binding cavity. We analyzed binding affinities for more than 80 peptides and monitored uptake by a fluorescence-based transport assay. To probe whether all 8400 natural di- and tripeptides can bind to DtpB, we employed state-of-the-art molecular docking and machine learning and conclude that peptides with compact hydrophobic residues are the best DtpB binders.

Global analysis of putative phospholipases in Plasmodium falciparum reveals an essential role of the phosphoinositide-specific phospholipase C in parasite maturation.

For its replication within red blood cells, the malaria parasite depends on a highly active and regulated lipid metabolism. Enzymes involved in lipid metabolic processes such as phospholipases are, therefore, potential drug targets. Here, using reverse genetics approaches, we show that only 1 out of the 19 putative phospholipases expressed in asexual blood stages of Plasmodium falciparum is essential for proliferation in vitro, pointing toward a high level of redundancy among members of this enzyme family. Using conditional mislocalization and gene disruption techniques, we show that this essential phosphoinositide-specific phospholipase C (PI-PLC, PF3D7_1013500) has a previously unrecognized essential role during intracellular parasite maturation, long before its previously perceived role in parasite egress and invasion. Subsequent lipidomic analysis suggests that PI-PLC mediates cleavage of phosphatidylinositol bisphosphate (PIP2) in schizont-stage parasites, underlining its critical role in regulating phosphoinositide levels in the parasite. IMPORTANCE The clinical symptoms of malaria arise due to repeated rounds of replication of Plasmodium parasites within red blood cells (RBCs). Central to this is an intense period of membrane biogenesis. Generation of membranes not only requires de novo synthesis and acquisition but also the degradation of phospholipids, a function that is performed by phospholipases. In this study, we investigate the essentiality of the 19 putative phospholipase enzymes that the human malaria parasite Plasmodium falciparum expresses during its replication within RBCs. We not only show that a high level of functional redundancy exists among these enzymes but, at the same time, also identify an essential role for the phosphoinositide-specific phospholipase C in parasite development and cleavage of the phospholipid phosphatidylinositol bisphosphate.

Plastid-localized xanthorhodopsin increases diatom biomass and ecosystem productivity in iron-limited surface oceans.

Microbial rhodopsins are photoreceptor proteins that convert light into biological signals or energy. Proteins of the xanthorhodopsin family are common in eukaryotic photosynthetic plankton including diatoms. However, their biological role in these organisms remains elusive. Here we report on a xanthorhodopsin variant (FcR1) isolated from the polar diatom Fragilariopsis cylindrus. Applying a combination of biophysical, biochemical and reverse genetics approaches, we demonstrate that FcR1 is a plastid-localized proton pump which binds the chromophore retinal and is activated by green light. Enhanced growth of a Thalassiora pseudonana gain-of-function mutant expressing FcR1 under iron limitation shows that the xanthorhodopsin proton pump supports growth when chlorophyll-based photosynthesis is iron-limited. The abundance of xanthorhodopsin transcripts in natural diatom communities of the surface oceans is anticorrelated with the availability of dissolved iron. Thus, we propose that these proton pumps convey a fitness advantage in regions where phytoplankton growth is limited by the availability of dissolved iron.

Cell biological analysis reveals an essential role for Pfcerli2 in erythrocyte invasion by malaria parasites.

Merozoite invasion of host red blood cells (RBCs) is essential for survival of the human malaria parasite Plasmodium falciparum. Proteins involved with RBC binding and invasion are secreted from dual-club shaped organelles at the apical tip of the merozoite called the rhoptries. Here we characterise P. falciparum Cytosolically Exposed Rhoptry Leaflet Interacting protein 2 (PfCERLI2), as a rhoptry bulb protein that is essential for merozoite invasion. Phylogenetic analyses show that cerli2 arose through an ancestral gene duplication of cerli1. We show that PfCERLI2 is essential for blood-stage growth and localises to the cytosolic face of the rhoptry bulb. Inducible knockdown of PfCERLI2 led to a proportion of merozoites failing to invade and was associated with elongation of the rhoptry organelle during merozoite development and inhibition of rhoptry antigen processing. These findings identify PfCERLI2 as a protein that has key roles in rhoptry biology during merozoite invasion.

PMRT1, a Plasmodium-Specific Parasite Plasma Membrane Transporter, Is Essential for Asexual and Sexual Blood Stage Development.

Membrane transport proteins perform crucial roles in cell physiology. The obligate intracellular parasite Plasmodium falciparum, an agent of human malaria, relies on membrane transport proteins for the uptake of nutrients from the host, disposal of metabolic waste, exchange of metabolites between organelles, and generation and maintenance of transmembrane electrochemical gradients for its growth and replication within human erythrocytes. Despite their importance for Plasmodium cellular physiology, the functional roles of a number of membrane transport proteins remain unclear, which is particularly true for orphan membrane transporters that have no or limited sequence homology to transporter proteins in other evolutionary lineages. Therefore, in the current study, we applied endogenous tagging, targeted gene disruption, conditional knockdown, and knockout approaches to investigate the subcellular localization and essentiality of six membrane transporters during intraerythrocytic development of P. falciparum parasites. They are localized at different subcellular structures-the food vacuole, the apicoplast, and the parasite plasma membrane-and four out of the six membrane transporters are essential during asexual development. Additionally, the plasma membrane resident transporter 1 (PMRT1; PF3D7_1135300), a unique Plasmodium-specific plasma membrane transporter, was shown to be essential for gametocytogenesis and functionally conserved within the genus Plasmodium. Overall, we reveal the importance of four orphan transporters to blood stage P. falciparum development, which have diverse intracellular localizations and putative functions. IMPORTANCE Plasmodium falciparum-infected erythrocytes possess multiple compartments with designated membranes. Transporter proteins embedded in these membranes not only facilitate movement of nutrients, metabolites, and other molecules between these compartments, but also are common therapeutic targets and can confer antimalarial drug resistance. Orphan membrane transporters in P. falciparum without sequence homology to transporters in other evolutionary lineages and divergent from host transporters may constitute attractive targets for novel intervention approaches. Here, we localized six of these putative transporters at different subcellular compartments and probed their importance during asexual parasite growth by using reverse genetic approaches. In total, only two candidates turned out to be dispensable for the parasite, highlighting four candidates as putative targets for therapeutic interventions. This study reveals the importance of several orphan transporters to blood stage P. falciparum development.

Characterization of Apicomplexan Amino Acid Transporters (ApiATs) in the Malaria Parasite Plasmodium falciparum.

During the symptomatic human blood phase, malaria parasites replicate within red blood cells. Parasite proliferation relies on the uptake of nutrients, such as amino acids, from the host cell and blood plasma, requiring transport across multiple membranes. Amino acids are delivered to the parasite through the parasite-surrounding vacuolar compartment by specialized nutrient-permeable channels of the erythrocyte membrane and the parasitophorous vacuole membrane (PVM). However, further transport of amino acids across the parasite plasma membrane (PPM) is currently not well characterized. In this study, we focused on a family of Apicomplexan amino acid transporters (ApiATs) that comprises five members in Plasmodium falciparum. First, we localized four of the P. falciparum ApiATs (PfApiATs) at the PPM using endogenous green fluorescent protein (GFP) tagging. Next, we applied reverse genetic approaches to probe into their essentiality during asexual replication and gametocytogenesis. Upon inducible knockdown and targeted gene disruption, a reduced asexual parasite proliferation was detected for PfApiAT2 and PfApiAT4. Functional inactivation of individual PfApiATs targeted in this study had no effect on gametocyte development. Our data suggest that individual PfApiATs are partially redundant during asexual in vitro proliferation and fully redundant during gametocytogenesis of P. falciparum parasites. IMPORTANCE Malaria parasites live and multiply inside cells. To facilitate their extremely fast intracellular proliferation, they hijack and transform their host cells. This also requires the active uptake of nutrients, such as amino acids, from the host cell and the surrounding environment through various membranes that are the consequence of the parasite's intracellular lifestyle. In this paper, we focus on a family of putative amino acid transporters termed ApiAT. We show expression and localization of four transporters in the parasite plasma membrane of Plasmodium falciparum-infected erythrocytes that represent one interface of the pathogen to its host cell. We probed into the impact of functional inactivation of individual transporters on parasite growth in asexual and sexual blood stages of P. falciparum and reveal that only two of them show a modest but significant reduction in parasite proliferation but no impact on gametocytogenesis, pointing toward dispensability within this transporter family.

Common virulence gene expression in adult first-time infected malaria patients and severe cases.

Sequestration of Plasmodium falciparum(P. falciparum)-infected erythrocytes to host endothelium through the parasite-derived P. falciparum erythrocyte membrane protein 1 (PfEMP1) adhesion proteins is central to the development of malaria pathogenesis. PfEMP1 proteins have diversified and expanded to encompass many sequence variants, conferring each parasite a similar array of human endothelial receptor-binding phenotypes. Here, we analyzed RNA-seq profiles of parasites isolated from 32 P. falciparum-infected adult travellers returning to Germany. Patients were categorized into either malaria naive (n = 15) or pre-exposed (n = 17), and into severe (n = 8) or non-severe (n = 24) cases. For differential expression analysis, PfEMP1-encoding var gene transcripts were de novo assembled from RNA-seq data and, in parallel, var-expressed sequence tags were analyzed and used to predict the encoded domain composition of the transcripts. Both approaches showed in concordance that severe malaria was associated with PfEMP1 containing the endothelial protein C receptor (EPCR)-binding CIDRα1 domain, whereas CD36-binding PfEMP1 was linked to non-severe malaria outcomes. First-time infected adults were more likely to develop severe symptoms and tended to be infected for a longer period. Thus, parasites with more pathogenic PfEMP1 variants are more common in patients with a naive immune status, and/or adverse inflammatory host responses to first infections favor the growth of EPCR-binding parasites.

Genome-Scale Metabolic Reconstruction and in Silico Perturbation Analysis of the Polar Diatom Fragilariopsis cylindrus Predicts High Metabolic Robustness.

Diatoms are major primary producers in polar environments where they can actively grow under extremely variable conditions. Integrative modeling using a genome-scale model (GSM) is a powerful approach to decipher the complex interactions between components of diatom metabolism and can provide insights into metabolic mechanisms underlying their evolutionary success in polar ecosystems. We developed the first GSM for a polar diatom, Fragilariopsis cylindrus, which enabled us to study its metabolic robustness using sensitivity analysis. We find that the predicted growth rate was robust to changes in all model parameters (i.e., cell biochemical composition) except the carbon uptake rate. Constraints on total cellular carbon buffer the effect of changes in the input parameters on reaction fluxes and growth rate. We also show that single reaction deletion of 20% to 32% of active (nonzero flux) reactions and single gene deletion of 44% to 55% of genes associated with active reactions affected the growth rate, as well as the production fluxes of total protein, lipid, carbohydrate, DNA, RNA, and pigments by less than 1%, which was due to the activation of compensatory reactions (e.g., analogous enzymes and alternative pathways) with more highly connected metabolites involved in the reactions that were robust to deletion. Interestingly, including highly divergent alleles unique for F. cylindrus increased its metabolic robustness to cellular perturbations even more. Overall, our results underscore the high robustness of metabolism in F. cylindrus, a feature that likely helps to maintain cell homeostasis under polar conditions.

Structure-Based Identification and Functional Characterization of a Lipocalin in the Malaria Parasite Plasmodium falciparum.

Proteins of the lipocalin family are known to bind small hydrophobic ligands and are involved in various physiological processes ranging from lipid transport to oxidative stress responses. The genome of the malaria parasite Plasmodium falciparum contains a single protein PF3D7_0925900 with a lipocalin signature. Using crystallography and small-angle X-ray scattering, we show that the protein has a tetrameric structure of typical lipocalin monomers; hence we name it P. falciparum lipocalin (PfLCN). We show that PfLCN is expressed in the intraerythrocytic stages of the parasite and localizes to the parasitophorous and food vacuoles. Conditional knockdown of PfLCN impairs parasite development, which can be rescued by treatment with the radical scavenger Trolox or by temporal inhibition of hemoglobin digestion. This suggests a key function of PfLCN in counteracting oxidative stress-induced cell damage during multiplication of parasites within erythrocytes.

Structural Insights Into PfARO and Characterization of its Interaction With PfAIP.

Apicomplexan parasites contain rhoptries, which are specialized secretory organelles that coordinate host cell invasion. During the process of invasion, rhoptries secrete their contents to facilitate interaction with, and entry into, the host cell. Here, we report the crystal structure of the rhoptry protein Armadillo Repeats-Only (ARO) from the human malaria parasite, Plasmodium falciparum (PfARO). The structure of PfARO comprises five tandem Armadillo-like (ARM) repeats, with adjacent ARM repeats stacked in a head-to-tail orientation resulting in PfARO adopting an elongated curved shape. Interestingly, the concave face of PfARO contains two distinct patches of highly conserved residues that appear to play an important role in protein-protein interaction. We functionally characterized the P. falciparum homolog of ARO interacting protein (PfAIP) and demonstrate that it localizes to the rhoptries. We show that conditional mislocalization of PfAIP leads to deficient red blood cell invasion. Guided by the structure, we identified mutations of PfARO that lead to mislocalization of PfAIP. Using proximity-based biotinylation we probe into PfAIP interacting proteins.

Dissecting the Gene Expression, Localization, Membrane Topology, and Function of the Plasmodium falciparum STEVOR Protein Family.

During its intraerythrocytic development, the malaria parasite Plasmodium falciparum exposes variant surface antigens (VSAs) on infected erythrocytes to establish and maintain an infection. One family of small VSAs is the polymorphic STEVOR proteins, which are marked for export to the host cell surface through their PEXEL signal peptide. Interestingly, some STEVORs have also been reported to localize to the parasite plasma membrane and apical organelles, pointing toward a putative function in host cell egress or invasion. Using deep RNA sequencing analysis, we characterized P. falciparumstevor gene expression across the intraerythrocytic development cycle, including free merozoites, in detail and used the resulting stevor expression profiles for hierarchical clustering. We found that most stevor genes show biphasic expression oscillation, with maximum expression during trophozoite stages and a second peak in late schizonts. We selected four STEVOR variants, confirmed the expected export of these proteins to the host cell membrane, and tracked them to a secondary location, either to the parasite plasma membrane or the secretory organelles of merozoites in late schizont stages. We investigated the function of a particular STEVOR that showed rhoptry localization and demonstrated its role at the parasite-host interface during host cell invasion by specific antisera and targeted gene disruption. Experimentally determined membrane topology of this STEVOR revealed a single transmembrane domain exposing the semiconserved as well as variable protein regions to the cell surface.IMPORTANCE Malaria claims about half a million lives each year. Plasmodium falciparum, the causative agent of the most severe form of the disease, uses proteins that are translocated to the surface of infected erythrocytes for immune evasion. To circumvent the detection of these gene products by the immune system, the parasite evolved a complex strategy that includes gene duplications and elaborate sequence polymorphism. STEVORs are one family of these variant surface antigens and are encoded by about 40 genes. Using deep RNA sequencing of blood-stage parasites, including free merozoites, we first established stevor expression of the cultured isolate and compared it with published transcriptomes. We reveal a biphasic expression of most stevor genes and confirm this for individual STEVORs at the protein level. The membrane topology of a rhoptry-associated variant was experimentally elucidated and linked to host cell invasion, underlining the importance of this multifunctional protein family for parasite proliferation.

Identifying metabolic pathways for production of extracellular polymeric substances by the diatom Fragilariopsis cylindrus inhabiting sea ice.

Diatoms are significant primary producers in sea ice, an ephemeral habitat with steep vertical gradients of temperature and salinity characterizing the ice matrix environment. To cope with the variable and challenging conditions, sea ice diatoms produce polysaccharide-rich extracellular polymeric substances (EPS) that play important roles in adhesion, cell protection, ligand binding and as organic carbon sources. Significant differences in EPS concentrations and chemical composition corresponding to temperature and salinity gradients were present in sea ice from the Weddell Sea and Eastern Antarctic regions of the Southern Ocean. To reconstruct the first metabolic pathway for EPS production in diatoms, we exposed Fragilariopsis cylindrus, a key bi-polar diatom species, to simulated sea ice formation. Transcriptome profiling under varying conditions of EPS production identified a significant number of genes and divergent alleles. Their complex differential expression patterns under simulated sea ice formation was aligned with physiological and biochemical properties of the cells, and with field measurements of sea ice EPS characteristics. Thus, the molecular complexity of the EPS pathway suggests metabolic plasticity in F. cylindrus is required to cope with the challenging conditions of the highly variable and extreme sea ice habitat.

Building a locally diploid genome and transcriptome of the diatom Fragilariopsis cylindrus.

The genome of the cold-adapted diatom Fragilariopsis cylindrus is characterized by highly diverged haplotypes that intersperse its homozygous genome. Here, we describe how a combination of PacBio DNA and Illumina RNA sequencing can be used to resolve this complex genomic landscape locally into the highly diverged haplotypes, and how to map various environmentally controlled transcripts onto individual haplotypes. We assembled PacBio sequence data with the FALCON assembler and created a haplotype resolved annotation of the assembly using annotations of a Sanger sequenced F. cylindrus genome. RNA-seq datasets from six different growth conditions were used to resolve allele-specifc gene expression in F. cylindrus. This approach enables to study differential expression of alleles in a complex genomic landscape and provides a useful tool to study how diverged haplotypes in diploid organisms are used for adaptation and evolution to highly variable environments.