Diatom modulation of select bacteria through use of two unique secondary metabolites.

Unicellular eukaryotic phytoplankton, such as diatoms, rely on microbial communities for survival despite lacking specialized compartments to house microbiomes (e.g., animal gut). Microbial communities have been widely shown to benefit from diatom excretions that accumulate within the microenvironment surrounding phytoplankton cells, known as the phycosphere. However, mechanisms that enable diatoms and other unicellular eukaryotes to nurture specific microbiomes by fostering beneficial bacteria and repelling harmful ones are mostly unknown. We hypothesized that diatom exudates may tune microbial communities and employed an integrated multiomics approach using the ubiquitous diatom Asterionellopsis glacialis to reveal how it modulates its naturally associated bacteria. We show that A. glacialis reprograms its transcriptional and metabolic profiles in response to bacteria to secrete a suite of central metabolites and two unusual secondary metabolites, rosmarinic acid and azelaic acid. While central metabolites are utilized by potential bacterial symbionts and opportunists alike, rosmarinic acid promotes attachment of beneficial bacteria to the diatom and simultaneously suppresses the attachment of opportunists. Similarly, azelaic acid enhances growth of beneficial bacteria while simultaneously inhibiting growth of opportunistic ones. We further show that the bacterial response to azelaic acid is numerically rare but globally distributed in the world's oceans and taxonomically restricted to a handful of bacterial genera. Our results demonstrate the innate ability of an important unicellular eukaryotic group to modulate select bacteria in their microbial consortia, similar to higher eukaryotes, using unique secondary metabolites that regulate bacterial growth and behavior inversely across different bacterial populations.

ß-actin regulates a heterochromatin landscape essential for optimal induction of neuronal programs during direct reprograming.

During neuronal development, ß-actin serves an important role in growth cone mediated axon guidance. Consistent with this notion, in vivo ablation of the ß-actin gene leads to abnormalities in the nervous system. However, whether ß-actin is involved in the regulation of neuronal gene programs is not known. In this study, we directly reprogramed ß-actin+/+ WT, ß-actin+/- HET and ß-actin-/- KO mouse embryonic fibroblast (MEFs) into chemically induced neurons (CiNeurons). Using RNA-seq analysis, we profiled the transcriptome changes among the CiNeurons. We discovered that induction of neuronal gene programs was impaired in KO CiNeurons in comparison to WT ones, whereas HET CiNeurons showed an intermediate levels of induction. ChIP-seq analysis of heterochromatin markers demonstrated that the impaired expression of neuronal gene programs correlated with the elevated H3K9 and H3K27 methylation levels at gene loci in ß-actin deficient MEFs, which is linked to the loss of chromatin association of the BAF complex ATPase subunit Brg1. Together, our study shows that heterochromatin alteration in ß-actin null MEFs impedes the induction of neuronal gene programs during direct reprograming. These findings are in line with the notion that H3K9Me3-based heterochromatin forms a major epigenetic barrier during cell fate change.

Fibroblast Differentiation and Matrix Remodeling Impaired under Simulated Microgravity in 3D Cell Culture Model.

Exposure to microgravity affects astronauts' health in adverse ways. However, less is known about the extent to which fibroblast differentiation during the wound healing process is affected by the lack of gravity. One of the key steps of this process is the differentiation of fibroblasts into myofibroblasts, which contribute functionally through extracellular matrix production and remodeling. In this work, we utilized collagen-based three-dimensional (3D) matrices to mimic interstitial tissue and studied fibroblast differentiation under simulated microgravity (sµG). Our results demonstrated that alpha-smooth muscle actin (αSMA) expression and translocation of Smad2/3 into the cell nucleus were reduced upon exposure to sµG compared to the 1g control, which suggests the impairment of fibroblast differentiation under sµG. Moreover, matrix remodeling and production were decreased under sµG, which is in line with the impaired fibroblast differentiation. We further investigated changes on a transcriptomic level using RNA sequencing. The results demonstrated that sµG has less effect on fibroblast transcriptomes, while sµG triggers changes in the transcriptome of myofibroblasts. Several genes and biological pathways found through transcriptome analysis have previously been reported to impair fibroblast differentiation. Overall, our data indicated that fibroblast differentiation, as well as matrix production and remodeling, are impaired in 3D culture under sµG conditions.

NASQAR: a web-based platform for high-throughput sequencing data analysis and visualization.

BACKGROUND: As high-throughput sequencing applications continue to evolve, the rapid growth in quantity and variety of sequence-based data calls for the development of new software libraries and tools for data analysis and visualization. Often, effective use of these tools requires computational skills beyond those of many researchers. To ease this computational barrier, we have created a dynamic web-based platform, NASQAR (Nucleic Acid SeQuence Analysis Resource). RESULTS: NASQAR offers a collection of custom and publicly available open-source web applications that make extensive use of a variety of R packages to provide interactive data analysis and visualization. The platform is publicly accessible at http://nasqar.abudhabi.nyu.edu/ . Open-source code is on GitHub at https://github.com/nasqar/NASQAR , and the system is also available as a Docker image at https://hub.docker.com/r/aymanm/nasqarall . NASQAR is a collaboration between the core bioinformatics teams of the NYU Abu Dhabi and NYU New York Centers for Genomics and Systems Biology. CONCLUSIONS: NASQAR empowers non-programming experts with a versatile and intuitive toolbox to easily and efficiently explore, analyze, and visualize their Transcriptomics data interactively. Popular tools for a variety of applications are currently available, including Transcriptome Data Preprocessing, RNA-seq Analysis (including Single-cell RNA-seq), Metagenomics, and Gene Enrichment.

ß-Actin-dependent global chromatin organization and gene expression programs control cellular identity.

During differentiation and development, cell fate and identity are established by waves of genetic reprogramming. Although the mechanisms are largely unknown, during these events, dynamic chromatin reorganization is likely to ensure that multiple genes involved in the same cellular functions are coregulated, depending on the nuclear environment. In this study, using high-content screening of embryonic fibroblasts from a ß-actin knockout (KO) mouse, we found major chromatin rearrangements and changes in histone modifications, such as methylated histone (H)3-lysine-(K)9. Genome-wide H3K9 trimethylation-(Me)3 landscape changes correlate with gene up- and down-regulation in ß-actin KO cells. Mechanistically, we found loss of chromatin association by the Brahma-related gene ( Brg)/Brahma-associated factor (BAF) chromatin remodeling complex subunit Brg1 in the absence of ß-actin. This actin-dependent chromatin reorganization was concomitant with the up-regulation of sets of genes involved in angiogenesis, cytoskeletal organization, and myofibroblast features in ß-actin KO cells. Some of these genes and phenotypes were gained in a ß-actin dose-dependent manner. Moreover, reintroducing a nuclear localization signal-containing ß-actin in the knockout cells affected nuclear features and gene expression. Our results suggest that, by affecting the genome-wide organization of heterochromatin through the chromatin-binding activity of the BAF complex, ß-actin plays an essential role in the determination of gene expression programs and cellular identity.-Xie, X., Almuzzaini, B., Drou, N., Kremb, S., Yousif, A., Östlund Farrants, A.-K., Gunsalus, K., Percipalle, P. ß-Actin-dependent global chromatin organization and gene expression programs control cellular identity.

Analyses of pig genomes provide insight into porcine demography and evolution.

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model.

RAMPART: a workflow management system for de novo genome assembly.

MOTIVATION: The de novo assembly of genomes from whole- genome shotgun sequence data is a computationally intensive, multi-stage task and it is not known a priori which methods and parameter settings will produce optimal results. In current de novo assembly projects, a popular strategy involves trying many approaches, using different tools and settings, and then comparing and contrasting the results in order to select a final assembly for publication. RESULTS: Herein, we present RAMPART, a configurable workflow management system for de novo genome assembly, which helps the user identify combinations of third-party tools and settings that provide good results for their particular genome and sequenced reads. RAMPART is designed to exploit High performance computing environments, such as clusters and shared memory systems, where available. AVAILABILITY AND IMPLEMENTATION: RAMPART is available under the GPLv3 license at: https://github.com/TGAC/RAMPART.

Tumor suppressor Par-4 activates autophagy-dependent ferroptosis.

Ferroptosis is a unique iron-dependent form of non-apoptotic cell death characterized by devastating lipid peroxidation. Whilst growing evidence suggests that ferroptosis is a type of autophagy-dependent cell death, the underlying molecular mechanisms regulating ferroptosis are largely unknown. In this study, through an unbiased RNA-sequencing screening, we demonstrate the activation of a multi-faceted tumor-suppressor protein Par-4/PAWR during ferroptosis. Functional studies reveal that genetic depletion of Par-4 effectively blocks ferroptosis, whereas Par-4 overexpression sensitizes cells to undergo ferroptosis. More importantly, we have determined that Par-4-triggered ferroptosis is mechanistically driven by the autophagic machinery. Upregulation of Par-4 promotes activation of ferritinophagy (autophagic degradation of ferritin) via the nuclear receptor co-activator 4 (NCOA4), resulting in excessive release of free labile iron and, hence, enhanced lipid peroxidation and ferroptosis. Inhibition of Par-4 dramatically suppresses the NCOA4-mediated ferritinophagy signaling axis. Our results also establish that Par-4 activation positively correlates with reactive oxygen species (ROS) production, which is critical for ferritinophagy-mediated ferroptosis. Furthermore, Par-4 knockdown effectively blocked ferroptosis-mediated tumor suppression in the mouse xenograft models. Collectively, these findings reveal that Par-4 has a crucial role in ferroptosis, which could be further exploited for cancer therapy.

Targeting oxeiptosis-mediated tumor suppression: a novel approach to treat colorectal cancers by sanguinarine.

Oxeiptosis is a recently identified reactive oxygen species (ROS)-sensitive, caspase independent, non-inflammatory regulated cell death pathway. The activation of Kelch-like ECH-associated protein 1-Phosphoglycerate mutase 5-Apoptosis inducing factor mitochondria associated 1 (KEAP1-PGAM5-AIFM1) pathway is the key signaling event in the execution of oxeiptosis. In the present study, we demonstrate that sanguinarine (SNG), a quaternary benzophenanthridine alkaloid, induces oxeiptosis in human colorectal cancer (CRC) cells via ROS, specifically hydrogen peroxide (H2O2)-dependent activation of KEAP1-PGAM5-AIFM1 signaling axis. Whilst, knockdown of KEAP1, PGAM5, and AIFM1 largely abolishes SNG-induced oxeiptosis, hence reinforcing the importance of the role of this pathway in the SNG-mediated cytotoxicity. Moreover, extracellular addition of H2O2 sensitizes SNG-induced oxeiptosis in CRC cells, while removal of intracellular ROS by ROS scavengers, not only alleviated the overproduction of ROS caused by SNG, but also reversed the biochemical events associated with oxeiptosis. Finally, in vivo study demonstrates that SNG effectively reduces the tumor growth in HT-29 xenograft mouse model through features associated with oxeiptosis. This study highlights oxeiptosis as a novel tumor suppressive mechanism and further investigation of the role of oxeiptosis in cancer treatment is warranted.

High-Throughput Sequencing Identified Distinct Bipartite and Monopartite Begomovirus Variants Associated with DNA-Satellites from Tomato and Muskmelon Plants in Saudi Arabia.

The studies on the prevalence and genetic diversity of begomoviruses in Saudi Arabia are minimal. In this study, field-grown symptomatic tomato and muskmelon plants were collected, and initially, begomovirus infection was confirmed by the core coat protein sequences. Four tomato and two muskmelon plants with viral infections were further evaluated for Illumina MiSeq sequencing, and twelve sequences (2.7-2.8 kb) equivalent to the full-length DNA-A or DNA-B components of begomoviruses were obtained along with eight sequences (~1.3-1.4 kb) equivalent to the begomovirus-associated DNA-satellite components. Four begomovirus sequences obtained from tomato plants were variants of tomato yellow leaf curl virus (TYLCV) with nt sequence identities of 95.3-100%. Additionally, two tomato plants showed a mixed infection of TYLCV and cotton leaf curl Gezira virus (CLCuGeV), okra yellow crinkle Cameroon alphasatellite (OYCrCMA), and okra leaf curl Oman betasatellite (OLCuOMB). Meanwhile, from muskmelon plants, two sequences were closely related (99-99.6%) to the tomato leaf curl Palampur virus (ToLCPalV) DNA-A, whereas two other sequences showed 97.9-100% sequence identities to DNA-B of ToLCPalV, respectively. Complete genome sequences of CLCuGeV and associated DNA-satellites were also obtained from these muskmelon plants. The nt sequence identities of the CLCuGeV, OYCrCMA, and OLCuOMB isolates obtained were 98.3-100%, 99.5-100%, and 95.6-99.7% with their respective available variants. The recombination was only detected in TYLCV and OLCuOMB isolates. To our knowledge, this is the first identification of a mixed infection of bipartite and monopartite begomoviruses associated with DNA-satellites from tomato and muskmelon in Saudi Arabia. The begomovirus variants reported in this study were clustered with Iranian isolates of respective begomovirus components in the phylogenetic dendrogram. Thus, the Iranian agroecological route can be a possible introduction of these begomoviruses and/or their associated DNA-satellites into Saudi Arabia.

Draft genome sequence of Streptomyces strain S4, a symbiont of the leaf-cutting ant Acromyrmex octospinosus.

Streptomyces spp. are common symbionts of the leaf-cutting ant species Acromyrmex octospinosus, which feeds on basidiomycete fungus leaf matter and harvests the lipid- and carbohydrate-rich gongylidia as a food source. A. octospinosus and other ant genera use antifungal compounds produced by Streptomyces spp. and other actinomycetes in order to help defend their fungal gardens from parasitic fungi. Herein, we report the draft genome sequence of Streptomyces strain S4, an antifungal-producing symbiont of A. octospinosus.

Complete genome sequence of the proteolytic Clostridium botulinum type A5 (B3') strain H04402 065.

H04402 065 is one of a very small group of strains of proteolytic Clostridium botulinum that form type A5 neurotoxin. Here, we report the complete 3.9-Mb genome sequence and annotation of strain H04402 065, which was isolated from a botulism patient in the United Kingdom in 2004.

A mixed community of actinomycetes produce multiple antibiotics for the fungus farming ant Acromyrmex octospinosus.

BACKGROUND: Attine ants live in an intensely studied tripartite mutualism with the fungus Leucoagaricus gongylophorus, which provides food to the ants, and with antibiotic-producing actinomycete bacteria. One hypothesis suggests that bacteria from the genus Pseudonocardia are the sole, co-evolved mutualists of attine ants and are transmitted vertically by the queens. A recent study identified a Pseudonocardia-produced antifungal, named dentigerumycin, associated with the lower attine Apterostigma dentigerum consistent with the idea that co-evolved Pseudonocardia make novel antibiotics. An alternative possibility is that attine ants sample actinomycete bacteria from the soil, selecting and maintaining those species that make useful antibiotics. Consistent with this idea, a Streptomyces species associated with the higher attine Acromyrmex octospinosus was recently shown to produce the well-known antifungal candicidin. Candicidin production is widespread in environmental isolates of Streptomyces, so this could either be an environmental contaminant or evidence of recruitment of useful actinomycetes from the environment. It should be noted that the two possibilities for actinomycete acquisition are not necessarily mutually exclusive. RESULTS: In order to test these possibilities we isolated bacteria from a geographically distinct population of A. octospinosus and identified a candicidin-producing Streptomyces species, which suggests that they are common mutualists of attine ants, most probably recruited from the environment. We also identified a Pseudonocardia species in the same ant colony that produces an unusual polyene antifungal, providing evidence for co-evolution of Pseudonocardia with A. octospinosus. CONCLUSIONS: Our results show that a combination of co-evolution and environmental sampling results in the diversity of actinomycete symbionts and antibiotics associated with attine ants.

BRCA1 and BRCA2 genes mutations among high risk breast cancer patients in Jordan.

Familial breast cancer is estimated to account for 15-20% of all cases of breast cancer. Surveillance for familial breast cancer is well-established world-wide. However, this service does not exist in Jordan, due to the scarcity of information with regard to the genetic profiling of these patients, and therefore lack of recommendations for policy-makers. As such, patients with very strong family history of breast or ovarian cancers are not screened routinely; leading to preventable delay in diagnosis. Whole coding sequencing for BCRA1/BCRA2 using next-generation sequencing (NGS)/Ion PGM System was performed. Sanger sequencing were then used to confirm the pathogenic variants detected by NGS. In this study, 192 breast cancer patients (and 8 ovarian cancer cases) were included. The prevalence of recurrent pathogenic mutations was 14.5%, while the prevalence of newly detected mutations was 3.5%. Two novel pathogenic mutations were identified in BRCA2 genes. The common mutations in the Ashkenazi population used for screening may not apply in the Jordanian population, as previously reported mutations were not prevalent, and other new mutations were identified. These data will aid to establish a specific screening test for BRCA 1/BRCA2 in the Jordanian population.

Integrative genomic analysis reveals mechanisms of immune evasion in P. falciparum malaria.

The mechanisms behind the ability of Plasmodium falciparum to evade host immune system are poorly understood and are a major roadblock in achieving malaria elimination. Here, we use integrative genomic profiling and a longitudinal pediatric cohort in Burkina Faso to demonstrate the role of post-transcriptional regulation in host immune response in malaria. We report a strong signature of miRNA expression differentiation associated with P. falciparum infection (127 out of 320 miRNAs, B-H FDR 5%) and parasitemia (72 miRNAs, B-H FDR 5%). Integrative miRNA-mRNA analysis implicates several infection-responsive miRNAs (e.g., miR-16-5p, miR-15a-5p and miR-181c-5p) promoting lymphocyte cell death. miRNA cis-eQTL analysis using whole-genome sequencing data identified 1,376 genetic variants associated with the expression of 34 miRNAs (B-H FDR 5%). We report a protective effect of rs114136945 minor allele on parasitemia mediated through miR-598-3p expression. These results highlight the impact of post-transcriptional regulation, immune cell death processes and host genetic regulatory control in malaria.

Complexity of genome evolution by segmental rearrangement in Brassica rapa revealed by sequence-level analysis.

BACKGROUND: The Brassica species, related to Arabidopsis thaliana, include an important group of crops and represent an excellent system for studying the evolutionary consequences of polyploidy. Previous studies have led to a proposed structure for an ancestral karyotype and models for the evolution of the B. rapa genome by triplication and segmental rearrangement, but these have not been validated at the sequence level. RESULTS: We developed computational tools to analyse the public collection of B. rapa BAC end sequence, in order to identify candidates for representing collinearity discontinuities between the genomes of B. rapa and A. thaliana. For each putative discontinuity, one of the BACs was sequenced and analysed for collinearity with the genome of A. thaliana. Additional BAC clones were identified and sequenced as part of ongoing efforts to sequence four chromosomes of B. rapa. Strikingly few of the 19 inter-chromosomal rearrangements corresponded to the set of collinearity discontinuities anticipated on the basis of previous studies. Our analyses revealed numerous instances of newly detected collinearity blocks. For B. rapa linkage group A8, we were able to develop a model for the derivation of the chromosome from the ancestral karyotype. We were also able to identify a rearrangement event in the ancestor of B. rapa that was not shared with the ancestor of A. thaliana, and is represented in triplicate in the B. rapa genome. In addition to inter-chromosomal rearrangements, we identified and analysed 32 BACs containing the end points of segmental inversion events. CONCLUSION: Our results show that previous studies of segmental collinearity between the A. thaliana, Brassica and ancestral karyotype genomes, although very useful, represent over-simplifications of their true relationships. The presence of numerous cryptic collinear genome segments and the frequent occurrence of segmental inversions mean that inference of the positions of genes in B. rapa based on the locations of orthologues in A. thaliana can be misleading. Our results will be of relevance to a wide range of plants that have polyploid genomes, many of which are being considered according to a paradigm of comprising conserved synteny blocks with respect to sequenced, related genomes.

A newly-developed community microarray resource for transcriptome profiling in Brassica species enables the confirmation of Brassica-specific expressed sequences.

BACKGROUND: The Brassica species include an important group of crops and provide opportunities for studying the evolutionary consequences of polyploidy. They are related to Arabidopsis thaliana, for which the first complete plant genome sequence was obtained and their genomes show extensive, although imperfect, conserved synteny with that of A. thaliana. A large number of EST sequences, derived from a range of different Brassica species, are available in the public database, but no public microarray resource has so far been developed for these species. RESULTS: We assembled unigenes using approximately 800,000 EST sequences, mainly from three species: B. napus, B. rapa and B. oleracea. The assembly was conducted with the aim of co-assembling ESTs of orthologous genes (including homoeologous pairs of genes in B. napus from each of the A and C genomes), but resolving assemblies of paralogous, or paleo-homoeologous, genes (i.e. the genes related by the ancestral genome triplication observed in diploid Brassica species). 90,864 unique sequence assemblies were developed. These were incorporated into the BAC sequence annotation for the Brassica rapa Genome Sequencing Project, enabling the identification of cognate genomic sequences for a proportion of them. A 60-mer oligo microarray comprising 94,558 probes was developed using the unigene sequences. Gene expression was analysed in reciprocal resynthesised B. napus lines and the B. oleracea and B. rapa lines used to produce them. The analysis showed that significant expression could consistently be detected in leaf tissue for 35,386 unigenes. Expression was detected across all four genotypes for 27,355 unigenes, genome-specific expression patterns were observed for 7,851 unigenes and 180 unigenes displayed other classes of expression pattern. Principal component analysis (PCA) clearly resolved the individual microarray datasets for B. rapa, B. oleracea and resynthesised B. napus. Quantitative differences in expression were observed between the resynthesised B. napus lines for 98 unigenes, most of which could be classified into non-additive expression patterns, including 17 that showed cytoplasm-specific patterns. We further characterized the unigenes for which A genome-specific expression was observed and cognate genomic sequences could be identified. Ten of these unigenes were found to be Brassica-specific sequences, including two that originate from complex loci comprising gene clusters. CONCLUSION: We succeeded in developing a Brassica community microarray resource. Although expression can be measured for the majority of unigenes across species, there were numerous probes that reported in a genome-specific manner. We anticipate that some proportion of these will represent species-specific transcripts and the remainder will be the consequence of variation of sequences within the regions represented by the array probes. Our studies demonstrated that the datasets obtained from the arrays can be used for typical analyses, including PCA and the analysis of differential expression. We have also demonstrated that Brassica-specific transcripts identified in silico in the sequence assembly of public EST database accessions are indeed reported by the array. These would not be detectable using arrays designed using A. thaliana sequences.

In Mitochondria ?-Actin Regulates mtDNA Transcription and Is Required for Mitochondrial Quality Control.

In eukaryotic cells, actin regulates both cytoplasmic and nuclear functions. However, whether actin-based structures are present in the mitochondria and are involved in mitochondrial functions has not been investigated. Here, using wild-type ?-actin +/+ and knockout (KO) ?-actin ?/? mouse embryonic fibroblasts we show evidence for the defect in maintaining mitochondrial membrane potential (MMP) in ?-actin-null cells. MMP defects were associated with impaired mitochondrial DNA (mtDNA) transcription and nuclear oxidative phosphorylation (OXPHOS) gene expression. Using super-resolution microscopy we provided direct evidence on the presence of ?-actin-containing structures inside mitochondria. Large aggregates of TFAM-stained nucleoids were observed in bulb-shaped mitochondria in KO cells, suggesting defects in mitochondrial nucleoid segregation without ?-actin. The observation that mitochondria-targeted ?-actin rescued mtDNA transcription and MMP suggests an indispensable functional role of a mitochondrial ?-actin pool necessary for mitochondrial quality control.

The genome and phenome of the green alga Chloroidium sp. UTEX 3007 reveal adaptive traits for desert acclimatization.

To investigate the phenomic and genomic traits that allow green algae to survive in deserts, we characterized a ubiquitous species, Chloroidium sp. UTEX 3007, which we isolated from multiple locations in the United Arab Emirates (UAE). Metabolomic analyses of Chloroidium sp. UTEX 3007 indicated that the alga accumulates a broad range of carbon sources, including several desiccation tolerance-promoting sugars and unusually large stores of palmitate. Growth assays revealed capacities to grow in salinities from zero to 60 g/L and to grow heterotrophically on >40 distinct carbon sources. Assembly and annotation of genomic reads yielded a 52.5 Mbp genome with 8153 functionally annotated genes. Comparison with other sequenced green algae revealed unique protein families involved in osmotic stress tolerance and saccharide metabolism that support phenomic studies. Our results reveal the robust and flexible biology utilized by a green alga to successfully inhabit a desert coastline.

The genomic and phenotypic diversity of Schizosaccharomyces pombe.

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.