Search

Background Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1–5 Gbp) and idiosyncratic genome features. Results Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species. Conclusions Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Reef-building corals live in a mutualistic relationship with photosynthetic algae (family Symbiodiniaceae) that usually provide most of the energy required by the coral host. This relationship is sensitive to temperature stress; as little as a 1°C increase often leads to the collapse of the association. This sensitivity has led to an interest in the potential of more stress-tolerant algae to supplement or substitute for the normal Symbiodiniaceae mutualists. In this respect, the apicomplexan-like microalga Chromera is of particular interest due to its greater temperature tolerance. We generated a de novo transcriptome for a Chromera strain isolated from a GBR coral ('GBR Chromera') and compared with those of the reference strain of Chromera ('Sydney Chromera'), and to those of Symbiodiniaceae (Fugacium kawagutii, Cladocopium goreaui and Breviolum minutum), as well as the apicomplexan parasite, Plasmodium falciparum. In contrast to the high sequence divergence amongst representatives of different genera within the family Symbiodiniaceae, the two Chromera strains featured low sequence divergence at orthologous genes, implying that they are likely to be conspecifics. Although KEGG categories provide few criteria by which true coral mutualists might be identified, they do supply a molecular rationalization that explains the ecological dominance of Cladocopium spp. amongst Indo-Pacific reef corals. The presence of HSP20 genes may contribute to the high thermal tolerance of Chromera.

Corals and the reef ecosystems that they support are in global decline due to increasing anthropogenic pressures such as climate change1. However, effective reef conservation strategies are hampered by a limited mechanistic understanding of coral biology and the functional roles of the diverse microbial communities that underpin coral health2,3. Here, we present an integrated genomic characterization of the coral species Porites lutea and its microbial partners. High-quality genomes were recovered from P. lutea, as well as a metagenome-assembled Cladocopium C15 (the dinoflagellate symbiont) and 52 bacterial and archaeal populations. Comparative genomic analysis revealed that many of the bacterial and archaeal genomes encode motifs that may be involved in maintaining association with the coral host and in supplying fixed carbon, B-vitamins and amino acids to their eukaryotic partners. Furthermore, mechanisms for ammonia, urea, nitrate, dimethylsulfoniopropionate and taurine transformation were identified that interlink members of the holobiont and may be important for nutrient acquisition and retention in oligotrophic waters. Our findings demonstrate the critical and diverse roles that microorganisms play within the coral holobiont and underscore the need to consider all of the components of the holobiont if we are to effectively inform reef conservation strategies.

Coral reefs are sustained by symbioses between corals and symbiodiniacean dinoflagellates. These symbioses vary in the extent of their permanence in and specificity to the host. Although dinoflagellates are primarily free-living, Symbiodiniaceae diversified mainly as symbiotic lineages. Their genomes reveal conserved symbiosis-related gene functions and high sequence divergence. However, the evolutionary mechanisms that underpin the transition from the free-living lifestyle to symbiosis remain poorly understood. Here, we discuss the genome evolution of Symbiodiniaceae in diverse ecological niches across the broad spectrum of symbiotic associations, from free-living to putative obligate symbionts. We pose key questions regarding genome evolution vis-à-vis the transition of dinoflagellates from free-living to symbiotic and propose strategies for future research to better understand coral-dinoflagellate and other eukaryote-eukaryote symbioses.

Modern microbial taxonomy generally relies on the use of single marker genes or sets of concatenated genes to generate a framework for the delineation and classification of organisms at different taxonomic levels. However, given that DNA is the 'blueprint of life', and hence the ultimate arbiter of taxonomy, classification systems should attempt to use as much of the blueprint as possible to capture a comprehensive phylogenetic signal. Recent analysis of whole-genome sequences from coral reef symbionts (dinoflagellates of the family Symbiodiniaceae) and other microalgal groups has uncovered extensive divergence not recognised by current algal taxonomic approaches. In the era of 'sequence everything', we argue that whole-genome data are pivotal to guide informed taxonomic inference, particularly for microbial eukaryotes.

Coral reefs are the epitome of species diversity, yet the number of described scleractinian coral species, the framework-builders of coral reefs, remains moderate by comparison. DNA sequencing studies are rapidly challenging this notion by exposing a wealth of undescribed diversity, but the evolutionary and ecological significance of this diversity remains largely unclear. Here, we present an annotated genome for one of the most ubiquitous corals in the Indo-Pacific (Pachyseris speciosa), and uncover through a comprehensive genomic and phenotypic assessment that it comprises morphologically indistinguishable, but ecologically divergent cryptic lineages. Demographic modelling based on whole-genome resequencing disproved that morphological crypsis was due to recent divergence, and instead indicated ancient morphological stasis. Although the lineages occur sympatrically across shallow and mesophotic habitats, extensive genotyping using a rapid diagnostic assay revealed differentiation of their ecological distributions. Leveraging “common garden” conditions facilitated by the overlapping distributions, we assessed physiological and quantitative skeletal traits and demonstrated concurrent phenotypic differentiation. Lastly, spawning observations of genotyped colonies highlighted the potential role of temporal reproductive isolation in the limited admixture, with consistent genomic signatures in genes related to morphogenesis and reproduction. Overall, our findings demonstrate how ecologically and phenotypically divergent coral species can evolve despite morphological stasis, and provide new leads into the potential mechanisms facilitating such divergence in sympatry. More broadly, they indicate that our current taxonomic framework for reef-building corals may be scratching the surface of the ecologically relevant diversity on coral reefs, consequently limiting our ability to protect or restore this diversity effectively.

Genomes of 150 coral colonies reveal evolutionary processes related to past climatic change on the Great Barrier Reef., Genetic signatures caused by demographic and adaptive processes during past climatic shifts can inform predictions of species’ responses to anthropogenic climate change. To identify these signatures in Acropora tenuis, a reef-building coral threatened by global warming, we first assembled the genome from long reads and then used shallow whole-genome resequencing of 150 colonies from the central inshore Great Barrier Reef to inform population genomic analyses. We identify population structure in the host that reflects a Pleistocene split, whereas photosymbiont differences between reefs most likely reflect contemporary (Holocene) conditions. Signatures of selection in the host were associated with genes linked to diverse processes including osmotic regulation, skeletal development, and the establishment and maintenance of symbiosis. Our results suggest that adaptation to post-glacial climate change in A. tenuis has involved selection on many genes, while differences in symbiont specificity between reefs appear to be unrelated to host population structure.

Approximately 10-15% of breast cancers are associated with a strong family history of disease. Pathogenic variants in BRCA1, BRCA2 or other moderate to highly penetrant susceptibility genes (e.g. TP53, ATM, CHEK2, PALB2 and PTEN) account for a number of breast cancer families. However, for over 50% of families the underlying genetic contribution to their risk remains unknown (termed here as non-BRCA1/2). This has a profound impact for how individuals and their families are managed in the clinic. We applied whole genome sequencing (WGS) to determine whether somatic mutation analysis can reveal insight into the aetiology of familial breast cancer. The full repertoire of somatic mutations was evaluated in 26 BRCA1, 22 BRCA2 and 32 non-BRCA1/2 tumours; including SNPs, indels, copy number changes and structural rearrangements, and mutational signatures. Genomes were also analysed using the HRD Index and HRDetect, as predictors of homologous recombination deficiency. BRCA1, BRCA2 and non-BRCA1/2 tumours exhibited a different burden of mutations, a different spectrum of mutational signatures and different telomere length. Based on collective patterns of mutation signatures, tumours were classified as 'BRCA1-like', 'BRCA2-like' or 'non-BRCA1/2-like' with a 15% rate of tumour re-classification from their original clinical BRCA status. The results demonstrate the power of WGS to differentiate between BRCA1 and BRCA2 driven tumours; in the identification of double-pathogenic germline mutation carriers based on the resulting somatic mutation signature; and in the interpretation of BRCA unclassified variants. WGS of tumour genomes reveals fascinating insights into tumour aetiology and could compliment current genetic testing of breast cancer families. Citation Format: Simpson P, Nones K, Johnson J, Newell F, Patch A-M, Thorne H, Kazakoff S, De Luca X, Parsons M, Ferguson K, Reid L, McCart Reed A, Srihari S, Lakis V, Davidson A, Mukhopadhyay P, Holmes O, Xu Q, Wood S, Leonard C, Beasley J, Degasperi A, Nik-Zainal S, Ragan M, Spurdle A, Khanna KK, Lakhani S, Pearson J, Chenevix-Trench G, Waddell N. Using whole genome sequencing and somatic mutation signatures to unravel insight into familial breast cancer aetiology [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P5-10-01.

Symbiosis between dinoflagellates of the genus Symbiodinium and reef-building corals forms the trophic foundation of the world’s coral reef ecosystems. Here we present the first draft genome of Symbiodinium goreaui (Clade C, type C1: 1.03 Gbp), one of the most ubiquitous endosymbionts associated with corals, and an improved draft genome of Symbiodinium kawagutii (Clade F, strain CS-156: 1.05 Gbp) to further elucidate genomic signatures of this symbiosis. Comparative analysis of four available Symbiodinium genomes against other dinoflagellate genomes led to the identification of 2460 nuclear gene families (containing 5% of Symbiodinium genes) that show evidence of positive selection, including genes involved in photosynthesis, transmembrane ion transport, synthesis and modification of amino acids and glycoproteins, and stress response. Further, we identify extensive sets of genes for meiosis and response to light stress. These draft genomes provide a foundational resource for advancing our understanding of Symbiodinium biology and the coral-algal symbiosis.

Since the discovery of Chromera velia as a novel coral-associated microalga, this organism has attracted interest because of its unique evolutionary position between the photosynthetic dinoflagellates and the parasitic apicomplexans. The nature of the relationship between Chromera and its coral host is controversial. Is it a mutualism, from which both participants benefit, a parasitic relationship, or a chance association? To better understand the interaction, larvae of the common Indo-Pacific reef-building coral Acropora digitifera were experimentally infected with Chromera, and the impact on the host transcriptome was assessed at 4, 12, and 48 h post-infection using Illumina RNA-Seq technology. The transcriptomic response of the coral to Chromera was complex and implies that host immunity is strongly suppressed, and both phagosome maturation and the apoptotic machinery is modified. These responses differ markedly from those described for infection with a competent strain of the coral mutualist Symbiodinium, instead resembling those of vertebrate hosts to parasites and/or pathogens such as Mycobacterium tuberculosis. Consistent with ecological studies suggesting that the association may be accidental, the transcriptional response of A. digitifera larvae leads us to conclude that Chromera could be a coral parasite, commensal, or accidental bystander, but certainly not a beneficial mutualist., 論文

Symbiodinium is best-known as the photosynthetic symbiont of corals, but some clades are symbiotic in other organisms or include free-living forms. Identifying similarities and differences among these clades can help us understand their relationship with corals, and thereby inform on measures to manage coral reefs in a changing environment. Here, using sequences from 24 publicly available transcriptomes and genomes of Symbiodinium, we assessed 78,389 gene families in Symbiodinium clades and the immediate outgroup Polarella glacialis, and identified putative overrepresented functions in gene families that (1) distinguish Symbiodinium from other members of Order Suessiales, (2) are shared by all of the Symbiodinium clades for which we have data, and (3) based on available information, are specific to each clade. Our findings indicate that transmembrane transport, mechanisms of response to reactive oxygen species, and protection against UV radiation are functions enriched in all Symbiodinium clades but not in P. glacialis. Enrichment of these functions indicates the capability of Symbiodinium to establish and maintain symbiosis, and to respond and adapt to its environment. The observed differences in lineage-specific gene families imply extensive genetic divergence among clades. Our results provide a platform for future investigation of lineage- or clade-specific adaptation of Symbiodinium to their environment.

Climate change at the Pleistocene/Holocene boundary reshaped many coastal landscapes, and provides an opportunity to study recent adaptive processes in marine species and ecosystems including coral reefs. On the Great Barrier Reef (GBR) sea level rise flooded a vast shelf creating a distinct inshore region which now harbours extensive coral assemblages despite being subject to relatively high turbidity, freshwater input and thermal fluctuations. To investigate how the coral holobiont has adapted to these conditions we first generated a highly contiguous genome assembly forAcropora tenuisbased on long-read sequencing, and then used shallow whole-genome resequencing of 148Acropora tenuiscolonies from five inshore locations to model demographic history, identify signatures of selection and profile symbiont communities. We show that corals from Magnetic Island, located in the central inshore region of the GBR, are genetically distinct from those 50-500km further north, reflecting a Pleistocene (250-600Kya) split, whereas photosymbiont genotypes differ between reefs in a pattern more likely to reflect contemporary (Holocene) conditions. We also identified loci in the coral host genome with signatures of positive selection in the northern population and used coalescent simulations to show that these are unlikely to be accounted for by demographic history. Genes at these loci have roles in a diverse range of processes that includes heterotrophic nutrition, osmotic regulation, skeletal development and the establishment and maintenance of symbiosis. Our results show that, in the case ofA. tenuisholobionts from the inshore GBR, the genomes of both the coral host and the primary photosymbiont of have been significantly shaped by their environment and illustrate the complexity of adaptations that have occurred in response to past climate change.

SummaryReef-building corals live in a mutualistic relationship with photosynthetic algae (family Symbiodiniaceae) that usually provide most of the energy required by the coral host. This relationship is sensitive to temperature stress; as little as a 1°C increase often leading to collapse of the association. This sensitivity has led to interest in the potential of more stress tolerant algae to supplement or substitute for the normal Symbiodiniaceae mutualists. In this respect, the apicomplexan-like microalga Chromera is of particular interest due to its greater temperature tolerance. We generated a de novo transcriptome for a Chromera strain isolated from a GBR coral (“GBR Chromera”) and compared to those of the reference strain of Chromera (“Sydney Chromera”), and to those of Symbiodiniaceae (Fugacium, Cladocopium and Breviolum), as well as the apicomplexan parasite, Plasmodium falciparum. By contrast with the Symbiodiniaceae, the two Chromera strains had a high level of sequence similarity evident by very low levels of divergence in orthologous genes. Although KEGG categories provide few criteria by which true coral mutualists might be identified, they do supply a molecular rationalization for the ubiquitous association of Cladocopium strains with Indo-Pacific reef corals. The presence of HSP20 genes may underlie the higher thermal tolerance of Chromera.

Background Dinoflagellates are taxonomically diverse and ecologically important phytoplankton that are ubiquitously present in marine and freshwater environments. Mostly photosynthetic, dinoflagellates provide the basis of aquatic primary production; most taxa are free-living, while some can form symbiotic and parasitic associations with other organisms. However, knowledge of the molecular mechanisms that underpin the adaptation of these organisms to diverse ecological niches is limited by the scarce availability of genomic data, partly due to their large genome sizes estimated up to 250 Gbp. Currently available dinoflagellate genome data are restricted to Symbiodiniaceae (particularly symbionts of reef-building corals) and parasitic lineages, from taxa that have smaller genome size ranges, while genomic information from more diverse free-living species is still lacking. Results Here, we present two draft diploid genome assemblies of the free-living dinoflagellate Polarella glacialis, isolated from the Arctic and Antarctica. We found that about 68% of the genomes are composed of repetitive sequence, with long terminal repeats likely contributing to intra-species structural divergence and distinct genome sizes (3.0 and 2.7 Gbp). For each genome, guided using full-length transcriptome data, we predicted > 50,000 high-quality protein-coding genes, of which ~40% are in unidirectional gene clusters and ~25% comprise single exons. Multi-genome comparison unveiled genes specific to P. glacialis and a common, putatively bacterial origin of ice-binding domains in cold-adapted dinoflagellates. Conclusions Our results elucidate how selection acts within the context of a complex genome structure to facilitate local adaptation. Because most dinoflagellate genes are constitutively expressed, Polarella glacialis has enhanced transcriptional responses via unidirectional, tandem duplication of single-exon genes that encode functions critical to survival in cold, low-light polar environments. These genomes provide a foundational reference for future research on dinoflagellate evolution.

Dinoflagellates of the family Symbiodiniaceae (Order Suessiales) are predominantly symbiotic, and many are known for their association with corals. The genetic and functional diversity among Symbiodiniaceae is well acknowledged, but the genome-wide sequence divergence among these lineages remains little known. Here, we present de novo genome assemblies of five isolates from the basal genus Symbiodinium, encompassing distinct ecological niches. Incorporating existing data from Symbiodiniaceae and other Suessiales (15 genome datasets in total), we investigated genome features that are common or unique to these Symbiodiniaceae, to genus Symbiodinium, and to the individual species S. microadriaticum and S. tridacnidorum. Our whole-genome comparisons reveal extensive sequence divergence, with no sequence regions common to all 15. Based on similarity of k-mers from whole-genome sequences, the distances among Symbiodinium isolates are similar to those between isolates of distinct genera. We observed extensive structural rearrangements among symbiodiniacean genomes; those from two distinct Symbiodinium species share the most (853) syntenic gene blocks. Functions enriched in genes core to Symbiodiniaceae are also enriched in those core to Symbiodinium. Gene functions related to symbiosis and stress response exhibit similar relative abundance in all analysed genomes. Our results suggest that structural rearrangements contribute to genome sequence divergence in Symbiodiniaceae even within a same species, but the gene functions have remained largely conserved in Suessiales. This is the first comprehensive comparison of Symbiodiniaceae based on whole-genome sequence data, including comparisons at the intra-genus and intra-species levels.

Symbiodiniaceae are predominantly symbiotic dinoflagellates critical to corals and other reef organisms. Symbiodinium is a basal symbiodiniacean lineage and includes symbiotic and free-living taxa. However, the molecular mechanisms underpinning these distinct lifestyles remain little known. Here, we present high-quality de novo genome assemblies for the symbiotic Symbiodinium tridacnidorum CCMP2592 (genome size 1.3 Gbp) and the free-living Symbiodinium natans CCMP2548 (genome size 0.74 Gbp). These genomes display extensive sequence divergence, sharing only ~1.5% conserved regions (≥90% identity). We predicted 45,474 and 35,270 genes for S. tridacnidorum and S. natans, respectively; of the 58,541 homologous gene families, 28.5% are common to both genomes. We recovered a greater extent of gene duplication and higher abundance of repeats, transposable elements and pseudogenes in the genome of S. tridacnidorum than in that of S. natans. These findings demonstrate that genome structural rearrangements are pertinent to distinct lifestyles in Symbiodinium, and may contribute to the vast genetic diversity within the genus, and more broadly in Symbiodiniaceae. Moreover, the results from our whole-genome comparisons against a free-living outgroup support the notion that the symbiotic lifestyle is a derived trait in, and that the free-living lifestyle is ancestral to, Symbiodinium.

Background Whole-genome sequencing (WGS) is a powerful method for revealing the diversity and complexity of the somatic mutation burden of tumours. Here, we investigated the utility of tumour and matched germline WGS for understanding aetiology and treatment opportunities for high-risk individuals with familial breast cancer. Patients and methods We carried out WGS on 78 paired germline and tumour DNA samples from individuals carrying pathogenic variants in BRCA1 (n = 26) or BRCA2 (n = 22) or from non-carriers (non-BRCA1/2; n = 30). Results Matched germline/tumour WGS and somatic mutational signature analysis revealed patients with unreported, dual pathogenic germline variants in cancer risk genes (BRCA1/BRCA2; BRCA1/MUTYH). The strategy identified that 100% of tumours from BRCA1 carriers and 91% of tumours from BRCA2 carriers exhibited biallelic inactivation of the respective gene, together with somatic mutational signatures suggestive of a functional deficiency in homologous recombination. A set of non-BRCA1/2 tumours also had somatic signatures indicative of BRCA-deficiency, including tumours with BRCA1 promoter methylation, and tumours from carriers of a PALB2 pathogenic germline variant and a BRCA2 variant of uncertain significance. A subset of 13 non-BRCA1/2 tumours from early onset cases were BRCA-proficient, yet displayed complex clustered structural rearrangements associated with the amplification of oncogenes and pathogenic germline variants in TP53, ATM and CHEK2. Conclusions Our study highlights the role that WGS of matched germline/tumour DNA and the somatic mutational signatures can play in the discovery of pathogenic germline variants and for providing supporting evidence for variant pathogenicity. WGS-derived signatures were more robust than germline status and other genomic predictors of homologous recombination deficiency, thus impacting the selection of platinum-based or PARP inhibitor therapy. In this first examination of non-BRCA1/2 tumours by WGS, we illustrate the considerable heterogeneity of these tumour genomes and highlight that complex genomic rearrangements may drive tumourigenesis in a subset of cases.

C.H. Waddington introduced the epigenetic landscape as a metaphor to represent cellular decision-making during development. Like a population of balls rolling down a rough hillside, developing cells follow specific trajectories (valleys) and eventually come to rest in one or another low-energy state that represents a mature cell type. Waddington depicted the topography of this landscape as determined by interactions among gene products, thereby connecting genotype to phenotype. In modern terms, each point on the landscape represents a state of the underlying genetic regulatory network, which in turn is described by a gene expression profile. In this chapter we demonstrate how the mathematical formalism of Hopfield networks can be used to model this epigenetic landscape. Hopfield networks are auto-associative artificial neural networks; input patterns are stored as attractors of the network and can be recalled from noisy or incomplete inputs. The resulting models capture the temporal dynamics of a gene regulatory network, yielding quantitative insight into cellular development and phenotype.

Despite the ecological significance of the relationship between reef-building corals and intracellular photosynthetic dinoflagellates of the genus Symbiodinium, very little is known about the molecular mechanisms involved in its establishment. Indeed, microarray-based analyses point to the conclusion that host gene expression is largely or completely unresponsive during the establishment of symbiosis with a competent strain of Symbiodinium. In this study, the use of Illumina RNA-Seq technology allowed detection of a transient period of differential expression involving a small number of genes (1073 transcripts

C.H. Waddington introduced the epigenetic landscape as a metaphor to represent cellular decision-making during development. Like a population of balls rolling down a rough hillside, developing cells follow specific trajectories (valleys) and eventually come to rest in one or another low-energy state that represents a mature cell type. Waddington depicted the topography of this landscape as determined by interactions among gene products, thereby connecting genotype to phenotype. In modern terms, each point on the landscape represents a state of the underlying genetic regulatory network, which in turn is described by a gene expression profile. In this chapter we demonstrate how the mathematical formalism of Hopfield networks can be used to model this epigenetic landscape. Hopfield networks are auto-associative artificial neural networks; input patterns are stored as attractors of the network and can be recalled from noisy or incomplete inputs. The resulting models capture the temporal dynamics of a gene regulatory network, yielding quantitative insight into cellular development and phenotype.

16 páginas, 1 tabla, 3 figuras., Mesenteric infection by the parasitic blood fluke Schistosoma bovis is a common veterinary problem in Africa and the Middle East and occasionally in the Mediterranean Region. The species also has the ability to form interspecific hybrids with the human parasite S. haematobium with natural hybridisation observed in West Africa, presenting possible zoonotic transmission. Additionally, this exchange of alleles between species may dramatically influence disease dynamics and parasite evolution. We have generated a 374 Mb assembly of the S. bovis genome using Illumina and PacBio-based technologies. Despite infecting different hosts and organs, the genome sequences of S. bovis and S. haematobium appeared strikingly similar with 97% sequence identity. The two species share 98% of protein-coding genes, with an average sequence identity of 97.3% at the amino acid level. Genome comparison identified large continuous parts of the genome (up to several 100 kb) showing almost 100% sequence identity between S. bovis and S. haematobium. It is unlikely that this is a result of genome conservation and provides further evidence of natural interspecific hybridization between S. bovis and S. haematobium. Our results suggest that foreign DNA obtained by interspecific hybridization was maintained in the population through multiple meiosis cycles and that hybrids were sexually reproductive, producing viable offspring. The S. bovis genome assembly forms a highly valuable resource for studying schistosome evolution and exploring genetic regions that are associated with species-specific phenotypic traits., This work was supported by grants from the QIMR Berghofer Medical Research Institute (Chenhall Estate) and the Australian Infectious Diseases Research Centre. NDY was supported by a National Health and Medical Research Council (NHMRC) CDF1, PKK by an NHMRC ECRF and DPM by an NHMRC SPRF (APP1102926). The work has been supported by NHMRC project grant APP1102322.

Background Foodborne infections caused by lung flukes of the genus Paragonimus are a significant and widespread public health problem in tropical areas. Approximately 50 Paragonimus species have been reported to infect animals and humans, but Paragonimus westermani is responsible for the bulk of human disease. Despite their medical and economic importance, no genome sequence for any Paragonimus species is available. Results We sequenced and assembled the genome of P. westermani, which is among the largest of the known pathogen genomes with an estimated size of 1.1 Gb. A 922.8 Mb genome assembly was generated from Illumina and Pacific Biosciences (PacBio) sequence data, covering 84% of the estimated genome size. The genome has a high proportion (45%) of repeat-derived DNA, particularly of the long interspersed element and long terminal repeat subtypes, and the expansion of these elements may explain some of the large size. We predicted 12,852 protein coding genes, showing a high level of conservation with related trematode species. The majority of proteins (80%) had homologs in the human liver fluke Opisthorchis viverrini, with an average sequence identity of 64.1%. Assembly of the P. westermani mitochondrial genome from long PacBio reads resulted in a single high-quality circularized 20.6 kb contig. The contig harbored a 6.9 kb region of non-coding repetitive DNA comprised of three distinct repeat units. Our results suggest that the region is highly polymorphic in P. westermani, possibly even within single worm isolates. Conclusions The generated assembly represents the first Paragonimus genome sequence and will facilitate future molecular studies of this important, but neglected, parasite group.

Dinoflagellates are a diverse group of unicellular primary producers and grazers that exhibit some of the most remarkable features known among eukaryotes. These include gigabase-sized nuclear genomes, permanently condensed chromosomes and highly reduced organelle DNA. However, the genetic inventory that allows dinoflagellates to thrive in diverse ecological niches is poorly characterised. Here we systematically assess the functional capacity of 3,368,684 predicted proteins from 47 transcriptome datasets spanning eight dinoflagellate orders. We find that 1,232,023 proteins do not share significant sequence similarity to known sequences, i.e. are “dark”. Of these, we consider 441,006 (13.1% of overall proteins) that are found in multiple taxa, or occur as alternative splice variants, to comprise the high-confidence dark proteins. Even with unknown function, 43.3% of these dark proteins can be annotated with conserved structural features using an exhaustive search against available data, validating their existence and importance. Furthermore, these dark proteins and their putative homologs are largely lineage-specific and recovered in multiple taxa. We also identified conserved functions in all dinoflagellates, and those specific to toxin-producing, symbiotic, and cold-adapted lineages. Our results demonstrate the remarkable divergence of gene functions in dinoflagellates, and provide a platform for investigations into the diversification of these ecologically important organisms.

Corals and the reef ecosystems that they support are in global decline due to increasing anthropogenic pressures such as climate change

[This corrects the article DOI: 10.1038/s42003-018-0098-3.].

The Homologous Recombination (HR) pathway is crucial for the repair of DNA double-strand breaks (DSBs) generated during DNA replication. Defects in HR repair have been linked to the initiation and development of a wide variety of human malignancies, and exploited in chemical, radiological and targeted therapies. In this study, we performed a personalised pathway analysis independently for four large sporadic breast cancer cohorts to investigate the status of HR pathway dysregulation in individual sporadic breast tumours, its association with HR repair deficiency and its impact on tumour characteristics. Specifically, we first manually curated a list of HR genes according to our recent review on this pathway (Liu et al., 2014), and then applied a personalised pathway analysis method named Pathifier (Drier et al., 2013) on the expression levels of the curated genes to obtain an HR score quantifying HR pathway dysregulation in individual tumours. Based on the score, we observed a great diversity in HR dysregulation between and within gene expression-based breast cancer subtypes, and by using two published HR-defect signatures, we found HR pathway dysregulation reflects HR repair deficiency. Furthermore, we identified a novel association between HR pathway dysregulation and chromosomal instability (CIN) in sporadic breast cancer. Although CIN has long been considered as a hallmark of most solid tumours, with recent extensive studies highlighting its importance in tumour evolution and drug resistance, the molecular basis of CIN in sporadic cancers remains poorly understood. Our results imply that HR pathway dysregulation might contribute to CIN in sporadic breast cancer.

Diazotrophic bacteria potentially supply substantial amounts of biologically fixed nitrogen to crops, but their occurrence may be suppressed by high nitrogen fertilizer application. Here, we explored the impact of high nitrogen fertilizer rates on the presence of diazotrophs in field-grown sugarcane with industry-standard or reduced nitrogen fertilizer application. Despite large differences in soil microbial communities between test sites, a core sugarcane root microbiome was identified. The sugarcane root-enriched core taxa overlap with those of Arabidopsis thaliana raising the possibility that certain bacterial families have had long association with plants. Reduced nitrogen fertilizer application had remarkably little effect on the core root microbiome and did not increase the relative abundance of root-associated diazotrophs or nif gene counts. Correspondingly, low nitrogen fertilizer crops had lower biomass and nitrogen content, reflecting a lack of major input of biologically fixed nitrogen, indicating that manipulating nitrogen fertilizer rates does not improve sugarcane yields by enriching diazotrophic populations under the test conditions. Standard nitrogen fertilizer crops had improved biomass and nitrogen content, and corresponding soils had higher abundances of nitrification and denitrification genes. These findings highlight that achieving a balance in maximizing crop yields and minimizing nutrient pollution associated with nitrogen fertilizer application requires understanding of how microbial communities respond to fertilizer use.

Bioinformatics has emerged as new discipline at the interface of molecular bioscience with mathematics, computer science and information technology. Bioinformatics is driven by data arising from high-throughput technologies in molecular bioscience. To enable biological discovery, bioinformatics draws on and extends technologies for data capture, management, integration and mining, computing, and communication technology. The rise of genomics has been a key driver for bioinformatics. Genomics, however, was never an end unto itself, but rather was intended to enable the understanding of complex biological systems. Bioinformatics continues to evolve in support of its constituent domains and, increasingly, their integration into genome-scale molecular systems biology. This article presents bioinformatics first from the perspective of computer science and information technology, then from the perspective of bioscience. In practice these perspectives often merge, making bioinformatics a rich, vibrant area of multidisciplinary research and application.

Symbiosis between dinoflagellates of the genusSymbiodiniumand reef-building corals forms the trophic foundation of the world’s coral reef ecosystems. Here we present the first draft genome ofSymbiodinium goreaui(Clade C, type C1: 1.03 Gbp), one of the most ubiquitous endosymbionts associated with corals, and an improved draft genome ofSymbiodinium kawagutii(Clade F, strain CS-156: 1.05 Gbp), previously sequenced as strain CCMP2468, to further elucidate genomic signatures of this symbiosis. Comparative analysis of four availableSymbiodiniumgenomes against other dinoflagellate genomes led to the identification of 2460 nuclear gene families that show evidence of positive selection, including genes involved in photosynthesis, transmembrane ion transport, synthesis and modification of amino acids and glycoproteins, and stress response. Further, we identified extensive sets of genes for meiosis and response to light stress. These draft genomes provide a foundational resource for advancing our understandingSymbiodiniumbiology and the coral-algal symbiosis.

Since the discovery of Chromera velia as a novel coral-associated microalga, this organism has attracted interest because of its unique evolutionary position between the photosynthetic dinoflagellates and the parasitic apicomplexans. The nature of the relationship between Chromera and its coral host is controversial. Is it a mutualism, from which both participants benefit, a parasitic relationship, or a chance association? To better understand the interaction, larvae of the common Indo-Pacific reef-building coral Acropora digitifera were experimentally infected with Chromera, and the impact on the host transcriptome was assessed at 4, 12, and 48 h post-infection using Illumina RNA-Seq technology. The transcriptomic response of the coral to Chromera was complex and implies that host immunity is strongly suppressed, and both phagosome maturation and the apoptotic machinery is modified. These responses differ markedly from those described for infection with a competent strain of the coral mutualist Symbiodinium, instead resembling those of vertebrate hosts to parasites and/or pathogens such as Mycobacterium tuberculosis. Consistent with ecological studies suggesting that the association may be accidental, the transcriptional response of A. digitifera larvae leads us to conclude that Chromera could be a coral parasite, commensal, or accidental bystander, but certainly not a beneficial mutualist.

Genome-wide regulatory networks enable cells to function, develop, and survive. Perturbation of these networks can lead to appearance of a disease phenotype. Inspired by Conrad Waddington's epigenetic landscape of cell development, we use a Hopfield network formalism to construct an attractor landscape model of disease progression based on protein- or gene-correlation networks of Parkinson's disease, glioma, and colorectal cancer. Attractors in this landscape correspond to normal and disease states of the cell. We introduce approaches to estimate the size and robustness of these attractors, and take a network-based approach to study their biological features such as the key genes and their functions associated with the attractors. Our results show that the attractor of cancer cells is wider than the attractor of normal cells, suggesting a heterogeneous nature of cancer. Perturbation analysis shows that robustness depends on characteristics of the input data (number of samples per time-point, and the fraction which converge to an attractor). We identify unique gene interactions at each stage, which reflect the temporal rewiring of the gene regulatory network (GRN) with disease progression. Our model of the attractor landscape, constructed from large-scale gene expression profiles of individual patients, captures snapshots of disease progression and identifies gene interactions specific to different stages, opening the way for development of stage-specific therapeutic strategies.

Genome evolution of microbes involves parent-to-offspring descent, and lateral genetic transfer that convolutes the phylogenomic signal. This study investigated phylogenomic signals among thousands of microbial genomes based on short subsequences without using multiple-sequence alignment. The signal from ribosomal RNAs is strong across all taxa, and the signal of plasmids is strong only in closely related groups, particularly Proteobacteria. However, the signal from other chromosomal regions (∼99% of the genomes) is remarkably restricted in breadth. The similarity of subsequences is found to correlate with taxonomic rank and informs on conserved and differential core functions relative to niche specialization and evolutionary diversification of microbes. These results provide a comprehensive, alignment-free view of microbial genome evolution as a network, beyond a tree-like structure., Microbial genomes have been shaped by parent-to-offspring (vertical) descent and lateral genetic transfer. These processes can be distinguished by alignment-based inference and comparison of phylogenetic trees for individual gene families, but this approach is not scalable to whole-genome sequences, and a tree-like structure does not adequately capture how these processes impact microbial physiology. Here we adopted alignment-free approaches based on k-mer statistics to infer phylogenomic networks involving 2,783 completely sequenced bacterial and archaeal genomes and compared the contributions of rRNA, protein-coding, and plasmid sequences to these networks. Our results show that the phylogenomic signal arising from ribosomal RNAs is strong and extends broadly across all taxa, whereas that from plasmids is strong but restricted to closely related groups, particularly Proteobacteria. However, the signal from the other chromosomal regions is restricted in breadth. We show that mean k-mer similarity can correlate with taxonomic rank. We also link the implicated k-mers to genome annotation (thus, functions) and define core k-mers (thus, core functions) in specific phyletic groups. Highly conserved functions in most phyla include amino acid metabolism and transport as well as energy production and conversion. Intracellular trafficking and secretion are the most prominent core functions among Spirochaetes, whereas energy production and conversion are not highly conserved among the largely parasitic or commensal Tenericutes. These observations suggest that differential conservation of functions relates to niche specialization and evolutionary diversification of microbes. Our results demonstrate that k-mer approaches can be used to efficiently identify phylogenomic signals and conserved core functions at the multigenome scale. IMPORTANCE Genome evolution of microbes involves parent-to-offspring descent, and lateral genetic transfer that convolutes the phylogenomic signal. This study investigated phylogenomic signals among thousands of microbial genomes based on short subsequences without using multiple-sequence alignment. The signal from ribosomal RNAs is strong across all taxa, and the signal of plasmids is strong only in closely related groups, particularly Proteobacteria. However, the signal from other chromosomal regions (∼99% of the genomes) is remarkably restricted in breadth. The similarity of subsequences is found to correlate with taxonomic rank and informs on conserved and differential core functions relative to niche specialization and evolutionary diversification of microbes. These results provide a comprehensive, alignment-free view of microbial genome evolution as a network, beyond a tree-like structure.

On November 18-19, 2016, the Human Frontier Science Program Organization (HFSPO) hosted a meeting of senior managers of key data resources and leaders of several major funding organizations to discuss the challenges associated with sustaining biological and biomedical (i.e., life sciences) data resources and associated infrastructure. A strong consensus emerged from the group that core data resources for the life sciences should be supported through coordinated international efforts that better ensure long-term sustainability and that appropriately align funding with scientific impact. Ideally, funding for such data resources should allow for access at no charge, as is presently the usual (and preferred) mechanism. Designing and implementing a plan for long-term accessibility is the vision. Below, the rationale for this vision is described, and some important considerations for developing a new international funding model to support core data resources for the life sciences are presented.

Bacteria and archaea can exchange genetic material across lineages through processes of lateral genetic transfer (LGT). Collectively, these exchange relationships can be modeled as a network and analyzed using concepts from graph theory. In particular, densely connected regions within an LGT network have been defined as genetic exchange communities (GECs). However, it has been problematic to construct networks in which edges solely represent LGT. Here we apply term frequency-inverse document frequency (TF-IDF), an alignment-free method originating from document analysis, to infer regions of lateral origin in bacterial genomes. We examine four empirical datasets of different size (number of genomes) and phyletic breadth, varying a key parameter (word length k) within bounds established in previous work. We map the inferred lateral regions to genes in recipient genomes, and construct networks in which the nodes are groups of genomes, and the edges natively represent LGT. We then extract maximum and maximal cliques (i.e., GECs) from these graphs, and identify nodes that belong to GECs across a wide range of k. Most surviving lateral transfer has happened within these GECs. Using Gene Ontology enrichment tests we demonstrate that biological processes associated with metabolism, regulation and transport are often over-represented among the genes affected by LGT within these communities. These enrichments are largely robust to change of k.

In breast cancer metastasis, the dynamic continuum involving pro- and anti-inflammatory regulators can become compromised. Over 600 genes have been implicated in metastasis to bone, lung or brain but how these genes might contribute to perturbation of immune function is poorly understood. To gain insight, we adopted a gene co-expression network approach that draws on the functional parallels between naturally occurring bone marrow-derived mesenchymal stem cells (BM-MSCs) and cancer stem cells (CSCs). Our network analyses indicate a key role for metastasis suppressor RARRES3, including potential to regulate the immunoproteasome (IP), a specialized proteasome induced under inflammatory conditions. Knockdown of RARRES3 in near-normal mammary epithelial and breast cancer cell lines increases overall transcript and protein levels of the IP subunits, but not of their constitutively expressed counterparts. RARRES3 mRNA expression is controlled by interferon regulatory factor IRF1, an inducer of the IP, and is sensitive to depletion of the retinoid-related receptor RORA that regulates various physiological processes including immunity through modulation of gene expression. Collectively, these findings identify a novel regulatory role for RARRES3 as an endogenous inhibitor of IP expression, and contribute to our evolving understanding of potential pathways underlying breast cancer driven immune modulation.

Lateral genetic transfer (LGT) is the process by which genetic material moves between organisms (and viruses) in the biosphere. Among the many approaches developed for the inference of LGT events from DNA sequence data, methods based on the comparison of phylogenetic trees remain the gold standard for many types of problem. Identifying LGT events from sequenced genomes typically involves a series of steps in which homologous sequences are identified and aligned, phylogenetic trees are inferred, and their topologies are compared to identify unexpected or conflicting relationships. These types of approach have been used to elucidate the nature and extent of LGT and its physiological and ecological consequences throughout the Tree of Life. Advances in DNA sequencing technology have led to enormous increases in the number of sequenced genomes, including ultra-deep sampling of specific taxonomic groups and single cell-based sequencing of unculturable "microbial dark matter." Environmental shotgun sequencing enables the study of LGT among organisms that share the same habitat.This abundance of genomic data offers new opportunities for scientific discovery, but poses two key problems. As ever more genomes are generated, the assembly and annotation of each individual genome receives less scrutiny; and with so many genomes available it is tempting to include them all in a single analysis, but thousands of genomes and millions of genes can overwhelm key algorithms in the analysis pipeline. Identifying LGT events of interest therefore depends on choosing the right dataset, and on algorithms that appropriately balance speed and accuracy given the size and composition of the chosen set of genomes.

Lateral genetic transfer (LGT) is the process by which genetic material moves between organisms (and viruses) in the biosphere. Among the many approaches developed for the inference of LGT events from DNA sequence data, methods based on the comparison of phylogenetic trees remain the gold standard for many types of problem. Identifying LGT events from sequenced genomes typically involves a series of steps in which homologous sequences are identified and aligned, phylogenetic trees are inferred, and their topologies are compared to identify unexpected or conflicting relationships. These types of approach have been used to elucidate the nature and extent of LGT and its physiological and ecological consequences throughout the Tree of Life. Advances in DNA sequencing technology have led to enormous increases in the number of sequenced genomes, including ultra-deep sampling of specific taxonomic groups and single cell-based sequencing of unculturable “microbial dark matter.” Environmental shotgun sequencing enables the study of LGT among organisms that share the same habitat. This abundance of genomic data offers new opportunities for scientific discovery, but poses two key problems. As ever more genomes are generated, the assembly and annotation of each individual genome receives less scrutiny; and with so many genomes available it is tempting to include them all in a single analysis, but thousands of genomes and millions of genes can overwhelm key algorithms in the analysis pipeline. Identifying LGT events of interest therefore depends on choosing the right dataset, and on algorithms that appropriately balance speed and accuracy given the size and composition of the chosen set of genomes.

Ernst Haeckel based his landmark Tree of Life on the supposed ontogenic recapitulation of phylogeny, i.e. that successive embryonic stages during the development of an organism re-trace the morphological forms of its ancestors over the course of evolution. Much of this idea has since been discredited. Today, phylogenies are often based on families of molecular sequences. The standard approach starts with a multiple sequence alignment, in which the sequences are arranged relative to each other in a way that maximises a measure of similarity position-by-position along their entire length. A tree (or sometimes a network) is then inferred. Rigorous multiple sequence alignment is computationally demanding, and evolutionary processes that shape the genomes of many microbes (bacteria, archaea and some morphologically simple eukaryotes) can add further complications. In particular, recombination, genome rearrangement and lateral genetic transfer undermine the assumptions that underlie multiple sequence alignment, and imply that a tree-like structure may be too simplistic. Here, using genome sequences of 143 bacterial and archaeal genomes, we construct a network of phylogenetic relatedness based on the number of shared k-mers (subsequences at fixed length k). Our findings suggest that the network captures not only key aspects of microbial genome evolution as inferred from a tree, but also features that are not treelike. The method is highly scalable, allowing for investigation of genome evolution across a large number of genomes. Instead of using specific regions or sequences from genome sequences, or indeed Haeckel’s idea of ontogeny, we argue that genome phylogenies can be inferred using k-mers from whole-genome sequences. Representing these networks dynamically allows biological questions of interest to be formulated and addressed quickly and in a visually intuitive manner.

Bacterial species in the plant-beneficial-environmental clade of Burkholderia represent a substantial component of rhizosphere microbes in many plant species. To better understand the molecular mechanisms of the interaction, we combined functional studies with high-resolution dual transcriptome analysis of sugarcane and root-associated diazotrophic Burkholderia strain Q208. We show that Burkholderia Q208 forms a biofilm at the root surface and suppresses the virulence factors that typically trigger immune response in plants. Up-regulation of bd-type cytochromes in Burkholderia Q208 suggests an increased energy production and creates the microaerobic conditions suitable for BNF. In this environment, a series of metabolic pathways are activated in Burkholderia Q208 implicated in oxalotrophy, microaerobic respiration, and formation of PHB granules, enabling energy production under microaerobic conditions. In the plant, genes involved in hypoxia survival are up-regulated and through increased ethylene production, larger aerenchyma is produced in roots which in turn facilitates diffusion of oxygen within the cortex. The detected changes in gene expression, physiology and morphology in the partnership are evidence of a sophisticated interplay between sugarcane and a plant-growth promoting Burkholderia species that advance our understanding of the mutually beneficial processes occurring in the rhizosphere.

Protein S-acylation is a reversible post-translational lipid modification that involves linkage of a fatty acid chain predominantly to a cysteine amino acid via a thioester bond. The fatty acid molecule is primarily palmitate, thus the term ‘palmitoylation’ is more commonly used. Palmitoylation has been found to modulate all stages of protein function including maturational processing, trafficking, membrane anchoring, signaling range and efficacy, and degradation. In breast cancer, palmitoylation has been shown to control the function of commonly dysregulated genes including estrogen receptors, the epidermal growth factor (EGF) family of receptors, and cancer stem cell markers. Importantly, palmitoylation is a critical factor controlling the formation of complexes at the plasma membrane involving tetraspanins, integrins, and gene products that are key to cell–cell communication. During metastasis, cancer cells enhance their metastatic capacity by interacting with stroma and immune cells. Although aberrant palmitoylation could contribute to tumor initiation and growth, its potential role in these cell–cell interactions is of particular interest, as it may provide mechanistic insight into metastasis, including cancer cell-driven immune modulation. Compelling evidence for a role for aberrant palmitoylation in breast cancer remains to be established. To this end, in this review we summarize emerging evidence and highlight pertinent knowledge gaps, suggesting directions for future research.