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.

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.

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

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.

In the original HTML version of the paper, the following affiliation was missing for author Cheong Xin Chan: School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD 4072, Australia. This affiliation was incorrectly assigned to author Pim Bongaerts as “Present Address”. All affiliations were published correctly in the PDF version of the paper and have now been corrected in the HTML.

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.

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.

Culture-independent molecular surveys of plant root microbiomes indicate that soil type generally has a stronger influence on microbial communities than host phylogeny. However, these studies have mostly focussed on model plants and crops. Here, we examine the root microbiomes of multiple plant phyla including lycopods, ferns, gymnosperms, and angiosperms across a soil chronosequence using 16S rRNA gene amplicon profiling. We confirm that soil type is the primary determinant of root-associated bacterial community composition, but also observe a significant correlation with plant phylogeny. A total of 47 bacterial genera are associated with roots relative to bulk soil microbial communities, including well-recognized plant-associated genera such as Bradyrhizobium, Rhizobium, and Burkholderia, and major uncharacterized lineages such as WPS-2, Ellin329, and FW68. We suggest that these taxa collectively constitute an evolutionarily conserved core root microbiome at this site. This lends support to the inference that a core root microbiome has evolved with terrestrial plants over their 400 million year history., Yeoh et al. study root microbiomes of different plant phyla across a tropical soil chronosequence. They confirm that soil type is the primary determinant of root-associated bacterial communities, but also observe a clear correlation with plant phylogeny and define a core root microbiome at this site.

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.

Accurate subtyping or classification of breast cancer is important for ensuring proper treatment of patients and also for understanding the molecular mechanisms driving this disease. While there have been several gene signatures proposed in the literature to classify breast tumours, these signatures show very low overlaps, different classification performance, and not much relevance to the underlying biology of these tumours. Here we evaluate DNA-damage response (DDR) and cell cycle pathways, which are critical pathways implicated in a considerable proportion of breast tumours, for their usefulness and ability in breast tumour subtyping. We think that subtyping breast tumours based on these two pathways could lead to vital insights into molecular mechanisms driving these tumours. Here, we performed a systematic evaluation of DDR and cell-cycle pathways for subtyping of breast tumours into the five known intrinsic subtypes. Homologous Recombination (HR) pathway showed the best performance in subtyping breast tumours, indicating that HR genes are strongly involved in all breast tumours. Comparisons of pathway based signatures and two standard gene signatures supported the use of known pathways for breast tumour subtyping. Further, the evaluation of these standard gene signatures showed that breast tumour subtyping, prognosis and survival estimation are all closely related. Finally, we constructed an all-inclusive super-signature by combining (union of) all genes and performing a stringent feature selection, and found it to be reasonably accurate and robust in classification as well as prognostic value. Adopting DDR and cell cycle pathways for breast tumour subtyping achieved robust and accurate breast tumour subtyping, and constructing a super-signature which contains feature selected mix of genes from these molecular pathways as well as clinical aspects is valuable in clinical practice.

It has been previously demonstrated that dicotyledonous plants perform rhizophagy, a process in which live microbial cells are engulfed by root cells and digested to acquire the nutrients from the microbes. Here we tested the hypothesis that rhizophagy is a mechanism of nutrient acquisition that is not restricted to dicotyledonous plants. We report that the monocotyledonous species sugarcane (Saccharum officinarum x spontaneum), grown in controlled axenic conditions, incorporated yeast cells into root cells. This suggests that rhizophagy is an evolutionarily conserved trait that predates the divergence of dicot and monocot species. To explore the potential relevance and practical application of rhizophagy, we investigated brewers' yeast (Saccharomyces cerevisiae), a waste product of the brewing industry, for its role as biofertilizer. The addition of live or dead yeast to fertilized soil substantially increased the nitrogen (N) and phosphorus (P) content of roots and shoots of tomato (Solanum lycopersicum) and young sugarcane plants. Yeast addition to soil also increased the root-to-shoot ratio in both species and induced species-specific morphological changes that included increased tillering in sugarcane and greater shoot biomass in tomato plants. These findings support the notion that brewers' yeast is a cost-effective biofertilizer that improves not only plant nutrition but also plant vigor during the early growth phase. It remains to be established which yeast-derived substances trigger the observed plant growth effects, and how rhizophagy contributes to plant nutrient acquisition.

From 1971 to 1985, Carl Woese and colleagues generated oligonucleotide catalogs of 16S/18S rRNAs from more than 400 organisms. Using these incomplete and imperfect data, Carl and his colleagues developed unprecedented insights into the structure, function, and evolution of the large RNA components of the translational apparatus. They recognized a third domain of life, revealed the phylogenetic backbone of bacteria (and its limitations), delineated taxa, and explored the tempo and mode of microbial evolution. For these discoveries to have stood the test of time, oligonucleotide catalogs must carry significant phylogenetic signal; they thus bear re-examination in view of the current interest in alignment-free phylogenetics based on k-mers. Here we consider the aims, successes, and limitations of this early phase of molecular phylogenetics. We computationally generate oligonucleotide sets (e-catalogs) from 16S/18S rRNA sequences, calculate pairwise distances between them based on D 2 statistics, compute distance trees, and compare their performance against alignment-based and k-mer trees. Although the catalogs themselves were superseded by full-length sequences, this stage in the development of computational molecular biology remains instructive for us today.

Sugarcane is a globally important food, biofuel and biomaterials crop. High nitrogen (N) fertilizer rates aimed at increasing yield often result in environmental damage because of excess and inefficient application. Inoculation with diazotrophic bacteria is an attractive option for reducing N fertilizer needs. However, the efficacy of bacterial inoculants is variable, and their effective formulation remains a knowledge frontier. Here, we take a new approach to investigating diazotrophic bacteria associated with roots using culture-independent microbial community profiling of a commercial sugarcane variety (Q208(A) ) in a field setting. We first identified bacteria that were markedly enriched in the rhizosphere to guide isolation and then tested putative diazotrophs for the ability to colonize axenic sugarcane plantlets (Q208(A) ) and promote growth in suboptimal N supply. One isolate readily colonized roots, fixed N2 and stimulated growth of plantlets, and was classified as a new species, Burkholderia australis sp. nov. Draft genome sequencing of the isolate confirmed the presence of nitrogen fixation. We propose that culture-independent identification and isolation of bacteria that are enriched in rhizosphere and roots, followed by systematic testing and confirming their growth-promoting capacity, is a necessary step towards designing effective microbial inoculants.

Recent advances in molecular technologies have opened up unprecedented opportunities for molecular ecologists to better understand the molecular basis of traits of ecological and evolutionary importance in almost any organism. Nevertheless, reliable and systematic inference of functionally relevant information from these masses of data remains challenging. The aim of this review is to highlight how the Gene Ontology (GO) database can be of use in resolving this challenge. The GO provides a largely species-neutral source of information on the molecular function, biological role and cellular location of tens of thousands of gene products. As it is designed to be species-neutral, the GO is well suited for cross-species use, meaning that, functional annotation derived from model organisms can be transferred to inferred orthologues in newly sequenced species. In other words, the GO can provide gene annotation information for species with nonannotated genomes. In this review, we describe the GO database, how functional information is linked with genes/gene products in model organisms, and how molecular ecologists can utilize this information to annotate their own data. Then, we outline various applications of GO for enhancing the understanding of molecular basis of traits in ecologically relevant species. We also highlight potential pitfalls, provide step-by-step recommendations for conducting a sound study in nonmodel organisms, suggest avenues for future research and outline a strategy for maximizing the benefits of a more ecological and evolutionary genomics-oriented ontology by ensuring its compatibility with the GO.

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.

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.

Lateral genetic transfer (LGT) plays an important role in the evolution of microbes. Existing computational methods for detecting genomic regions of putative lateral origin scale poorly to large data. Here, we propose a novel method based on TF-IDF (Term Frequency-Inverse Document Frequency) statistics to detect not only regions of lateral origin, but also their origin and direction of transfer, in sets of hierarchically structured nucleotide or protein sequences. This approach is based on the frequency distributions of k-mers in the sequences. If a set of contiguous k-mers appears sufficiently more frequently in another phyletic group than in its own, we infer that they have been transferred from the first group to the second. We performed rigorous tests of TF-IDF using simulated and empirical datasets. With the simulated data, we tested our method under different parameter settings for sequence length, substitution rate between and within groups and post-LGT, deletion rate, length of transferred region and k size and found that we can detect LGT events with high precision and recall. Our method performs better than an established method, ALFY, which has high recall but low precision. Our method is efficient, with runtime increasing approximately linearly with sequence length.

Many microbes can acquire genetic material from their environment and incorporate it into their genome, a process known as lateral genetic transfer (LGT). Computational approaches have been developed to detect genomic regions of lateral origin, but typically lack sensitivity, ability to distinguish donor from recipient and scalability to very large datasets. To address these issues we have introduced an alignment-free method based on ideas from document analysis, term frequency-inverse document frequency (TF-IDF). Here we examine the performance of TF-IDF on three empirical datasets: 27 genomes of Escherichia coli and Shigella, 110 genomes of enteric bacteria and 143 genomes across 12 bacterial and three archaeal phyla. We investigate the effect of k-mer size, gap size and delineation of groups on the inference of genomic regions of lateral origin, finding an interplay among these parameters and sequence divergence. Because TF-IDF identifies donor groups and delineates regions of lateral origin within recipient genomes, aggregating these regions by gene enables us to explore, for the first time, the mosaic nature of lateral genes including the multiplicity of biological sources, ancestry of transfer and over-writing by subsequent transfers. We carry out Gene Ontology enrichment tests to investigate which biological processes are potentially affected by LGT.

Alignment-free (AF) approaches have recently been highlighted as alternatives to methods based on multiple sequence alignment in phylogenetic inference. However, the sensitivity of AF methods to genome-scale evolutionary scenarios is little known. Here, using simulated microbial genome data we systematically assess the sensitivity of nine AF methods to three important evolutionary scenarios: sequence divergence, lateral genetic transfer (LGT) and genome rearrangement. Among these, AF methods are most sensitive to the extent of sequence divergence, less sensitive to low and moderate frequencies of LGT and most robust against genome rearrangement. We describe the application of AF methods to three well-studied empirical genome datasets and introduce a new application of the jackknife to assess node support. Our results demonstrate that AF phylogenomics is computationally scalable to multi-genome data and can generate biologically meaningful phylogenies and insights into microbial evolution.

Small RNAs (sRNAs) are widespread in bacteria and play critical roles in regulating physiological processes. They are best characterized in Escherichia coli K-12 MG1655, where 83 sRNAs constitute nearly 2% of the gene complement. Most sRNAs act by base pairing with a target mRNA, modulating its translation and/or stability; many of these RNAs share only limited complementarity to their mRNA target, and require the chaperone Hfq to facilitate base pairing. Little is known about the evolutionary dynamics of bacterial sRNAs. Here, we apply phylogenetic and network analyses to investigate the evolutionary processes and principles that govern sRNA gene distribution in 27 E. coli and Shigella genomes. We identify core (encoded in all 27 genomes) and variable sRNAs; more than two-thirds of the E. coli K-12 MG1655 sRNAs are core, whereas the others show patterns of presence and absence that are principally due to genetic loss, not duplication or lateral genetic transfer. We present evidence that variable sRNAs are less tightly integrated into cellular genetic regulatory networks than are the core sRNAs, and that Hfq facilitates posttranscriptional cross talk between the E. coli-Shigella core and variable genomes. Finally, we present evidence that more than 80% of genes targeted by Hfq-associated core sRNAs have been transferred within the E. coli-Shigella clade, and that most of these genes have been transferred intact. These results suggest that Hfq and sRNAs help integrate laterally acquired genes into established regulatory networks.

The widespread presence of antibiotic resistance and virulence among Staphylococcus isolates has been attributed in part to lateral genetic transfer (LGT), but little is known about the broader extent of LGT within this genus. Here we report the first systematic study of the modularity of genetic transfer among 13 Staphylococcus genomes covering four distinct named species. Using a topology-based phylogenetic approach, we found, among 1,354 sets of homologous genes examined, strong evidence of LGT in 368 (27.1%) gene sets, and weaker evidence in another 259 (19.1%). Within-gene and whole-gene transfer contribute almost equally to the topological discordance of these gene sets against a reference phylogeny. Comparing genetic transfer in single-copy and in multicopy gene sets, we observed a higher frequency of LGT in the latter, and a substantial functional bias in cases of whole-gene transfer (little such bias was observed in cases of fragmentary genetic transfer). We found evidence that lateral transfer, particularly of entire genes, impacts not only functions related to antibiotic, drug, and heavy-metal resistance, as well as membrane transport, but also core informational and metabolic functions not associated with mobile elements. Although patterns of sequence similarity support the cohesion of recognized species, LGT within S. aureus appears frequently to disrupt clonal complexes. Our results demonstrate that LGT and gene duplication play important parts in functional innovation in staphylococcal genomes.

On Dec 13, 2010, the phycological community lost one of its most distinguished members, with the passing of Jack L. McLachlan. McLachlan was born April 1, 1930, in Huron, SD, and lived in that state until his family moved to Anacortes, WA, in 1938. His postsecondary education was undertaken first at Western Washington College in Bellingham and then at Oregon State College, Corvallis. There he earned his undergraduate and graduate degrees, completing his PhD in 1957 on culture and autecology of Aphanizomenon flos-aquae, under the direction of Harry K. Phinney. Teaching was a large part of McLachlan's career and, like most National Research Council of Canada scientists, he had cross-appointments to local universities as a lecturer and student research supervisor, the latter being supported in part by his external research grants. Among other major honors McLachlan received were election to the Royal Society of Canada (1979) and the George Lawson Medal from the Canadian Botanical Association, for cumulative lifetime contributions to Canadian botany (1992).

Protein-protein interactions (PPIs) underpin key biological processes in cells. Describing and interpreting this network of interactions is a key focus of computational systems biology. While mouse is commonly used as a model system for mammalian biology, the description of mouse PPIs available in public interaction databases is remarkably poor. Collectively, public resources such as BIND, IntACT and MINT contain only a few thousand mouse PPIs, far behind the many tens or hundreds of thousands likely to exist. To supplement this lack and to take advantage of other high-throughput omic data sets in mouse, here we identify that portion of the human interactome with orthologs in mouse, and from that infer a mouse interolog network. By inferring interactions in mouse based on only the most closely related species with abundant PPI data (human), we create a view of mouse interactions enriched for shared mammalian biological processes. We also demonstrate that available methods for determining orthologs between even closely related species produce distinctly different results, and we propose an integrated view of mouse-human orthology from which to infer a broader interolog network.

Lateral genetic transfer (LGT) is an important adaptive force in evolution, contributing to metabolic, physiological and ecological innovation in most prokaryotes and some eukaryotes. Genomic sequences and other data have begun to illuminate the processes, mechanisms, quantitative extent and impact of LGT in diverse organisms, populations, taxa and environments; deep questions are being posed, and the provisional answers sometimes challenge existing paradigms. At the same time, there is an enhanced appreciation of the imperfections, biases and blind spots in the data and in analytical approaches. Here we identify and consider significant open questions concerning the role of LGT in genome evolution.

The significance of whole-genome duplications (WGD) for vertebrate evolution remains controversial, in part because the mechanisms by which WGD contributed to functional evolution or speciation are still incompletely characterized. Fish genomes provide an ideal context in which to examine the consequences of WGD, because the teleost lineage experienced an additional WGD soon after divergence from tetrapods and because five teleost genomes are available for comparative analysis. Here we present an integrated approach to characterize these post-duplication genomes based on genome-scale synteny, phylogenetic, temporal, and spatial gene expression and on protein sequence data. A minimum of 3%–4% of protein-coding loci have been retained in two copies in each of the five fish genomes, and many of these duplicates are key developmental genes that function as transcription factors or signaling molecules. Almost all duplicate gene pairs we examined have diverged in spatial and/or temporal expression during embryogenesis. A quarter of duplicate pairs have diverged in function via the acquisition of novel protein domains or via changes in the subcellular localization of their encoded proteins. We compared the spatial expression and protein domain architecture of zebrafish WGD-duplicates to those of their single mouse ortholog and found many examples supporting a model of neofunctionalization. WGD-duplicates have acquired novel protein domains more often than have single-copy genes. Post-WGD changes at the gene regulatory level were more common than changes at the protein level. We conclude that the most significant consequence of WGD for vertebrate evolution has been to enable more-specialized regulatory control of development via the acquisition of novel spatiotemporal expression domains. We find limited evidence that reciprocal gene loss led to reproductive isolation and speciation in this lineage.

Pairwise local sequence alignment methods have been the prevailing technique to identify homologous nucleotides between related species. However, existing methods that identify and align all homologous nucleotides in one or more genomes have suffered from poor scalability and limited accuracy. We propose a novel method that couples a gapped extension heuristic with an efficient filtration method for identifying interspersed repeats in genome sequences. During gapped extension, we use the MUSCLE implementation of progressive global multiple alignment with iterative refinement. The resulting gapped extensions potentially contain alignments of unrelated sequence. We detect and remove such undesirable alignments using a hidden Markov model (HMM) to predict the posterior probability of homology. The HMM emission frequencies for nucleotide substitutions can be derived from any time-reversible nucleotide substitution matrix. We evaluate the performance of our method and previous approaches on a hybrid data set of real genomic DNA with simulated interspersed repeats. Our method outperforms a related method in terms of sensitivity, positive predictive value, and localizing boundaries of homology. The described methods have been implemented in freely available software, Repeatoire, available from: http://wwwabi.snv.jussieu.fr/public/Repeatoire. © 2006 IEEE.

Lateral genetic transfer (LGT) involves the movement of genetic material from one lineage into another and its subsequent incorporation into the new host genome via genetic recombination. Studies in individual taxa have indicated lateral origins for stretches of DNA of greatly varying length, from a few nucleotides to chromosome size. Here we analyze 1,462 sets of single-copy, putatively orthologous genes from 144 fully sequenced prokaryote genomes, asking to what extent complete genes and fragments of genes have been transferred and recombined in LGT. Using a rigorous phylogenetic approach, we find evidence for LGT in at least 476 (32.6%) of these 1,462 gene sets: 286 (19.6%) clearly show one or more "observable recombination breakpoints" within the boundaries of the open reading frame, while a further 190 (13.0%) yield trees that are topologically incongruent with the reference tree but do not contain a recombination breakpoint within the open reading frame. We refer to these gene sets as observable recombination breakpoint positive (ORB(+)) and negative (ORB(-)) respectively. The latter are prima facie instances of lateral transfer of an entire gene or beyond. We observe little functional bias between ORB(+) and ORB(-) gene sets, but find that incorporation of entire genes is potentially more frequent in pathogens than in nonpathogens. As ORB(+) gene sets are about 50% more common than ORB(-) sets in our data, the transfer of gene fragments has been relatively frequent, and the frequency of LGT may have been systematically underestimated in phylogenetic studies.

The genome sequence is an icon of early twenty-first century biology. Genomes of nearly 2000 cellular organisms, and from many thousands of organelles and viruses, are now in the public domain. For biological research in individual species, the genome sequence increasingly provides the common

The epigenetic landscape was introduced by Conrad Waddington as a metaphor of cellular development. Like a ball rolling down a hillside is channelled through a succession of valleys until it reaches the bottom, cells follow specific trajectories from a pluripotent state to a committed state. Transcription factors (TFs) interacting as a network (the gene regulatory network (GRN)) orchestrate this developmental process within each cell. Here, we quantitatively model the epigenetic landscape using a kind of artificial neural network called the Hopfield network (HN). An HN is composed of nodes (genes/TFs) and weighted undirected edges, resulting in a weight matrix (W) that stores interactions among the nodes over the entire network. We used gene co-expression to compute the edge weights. Through W, we then associate an energy score (E) to each input pattern (pattern of co-expression for a specific developmental stage) such that each pattern has a specific E. We propose that, based on the co-expression values stored in W, HN associates lower E values to stable phenotypic states and higher E to transient states. We validate our model using time course gene-expression data sets representing stages of development across 12 biological processes including differentiation of human embryonic stem cells into specialized cells, differentiation of THP1 monocytes to macrophages during immune response and trans-differentiation of epithelial to mesenchymal cells in cancer. We observe that transient states have higher energy than the stable phenotypic states, yielding an arc-shaped trajectory. This relationship was confirmed by perturbation analysis. HNs offer an attractive framework for quantitative modelling of cell differentiation (as a landscape) from empirical data. Using HNs, we identify genes and TFs that drive cell-fate transitions, and gain insight into the global dynamics of GRNs.

Copy number alterations (CNAs) are thought to account for 85% of the variation in gene expression observed among breast tumours. The expression of cis-associated genes is impacted by CNAs occurring at proximal loci of these genes, whereas the expression of trans-associated genes is impacted by CNAs occurring at distal loci. While a majority of these CNA-driven genes responsible for breast tumourigenesis are cis-associated, trans-associated genes are thought to further abet the development of cancer and influence disease outcomes in patients. Here we present a network-based approach that integrates copy-number and expression profiles to identify putative cis- and trans-associated genes in breast cancer pathogenesis. We validate these cis- and trans-associated genes by employing them to subtype a large cohort of breast tumours obtained from the METABRIC consortium, and demonstrate that these genes accurately reconstruct the ten subtypes of breast cancer. We observe that individual breast cancer subtypes are driven by distinct sets of cis- and trans-associated genes. Among the cis-associated genes, we recover several known drivers of breast cancer (e.g. CCND1, ERRB2, MDM2 and ZNF703) and some novel putative drivers (e.g. BRF2 and SF3B3). siRNA-mediated knockdown of BRF2 across a panel of breast cancer cell lines showed significant reduction specifically in cell proliferation in HER2+ lines, thereby indicating that BRF2 could be a context-dependent oncogene and potentially targetable in these lines. Among the trans-associated genes, we identify modules of immune-response (CD2, CD19, CD38 and CD79B), mitotic/cell-cycle kinases (e.g. AURKB, MELK, PLK1 and TTK), and DNA-damage response genes (e.g. RFC4 and FEN1)., Comment: 23 pages, 2 tables, 7 figures

The process of inferring phylogenetic trees from molecular sequences almost always starts with a multiple alignment of these sequences but can also be based on methods that do not involve multiple sequence alignment. Very little is known about the accuracy with which such alignment-free methods recover the correct phylogeny or about the potential for increasing their accuracy. We conducted a large-scale comparison of ten alignment-free methods, among them one new approach that does not calculate distances and a faster variant of our pattern-based approach; all distance-based alignment-free methods are freely available from http://www.bioinformatics.org.au (as Python package decaf+py). We show that most methods exhibit a higher overall reconstruction accuracy in the presence of high among-site rate variation. Under all conditions that we considered, variants of the pattern-based approach were significantly better than the other alignment-free methods. The new pattern-based variant achieved a speed-up of an order of magnitude in the distance calculation step, accompanied by a small loss of tree reconstruction accuracy. A method of Bayesian inference from k-mers did not improve on classical alignment-free (and distance-based) methods but may still offer other advantages due to its Bayesian nature. We found the optimal word length k of word-based methods to be stable across various data sets, and we provide parameter ranges for two different alphabets. The influence of these alphabets was analyzed to reveal a trade-off in reconstruction accuracy between long and short branches. We have mapped the phylogenetic accuracy for many alignment-free methods, among them several recently introduced ones, and increased our understanding of their behavior in response to biologically important parameters. In all experiments, the pattern-based approach emerged as superior, at the expense of higher resource consumption. Nonetheless, no alignment-free method that we examined recovers the correct phylogeny as accurately as does an approach based on maximum-likelihood distance estimates of multiply aligned sequences.