Whole-genome sequencing reveals clinically relevant insights into the aetiology of familial breast cancers.

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.

The transcriptomic response of the coral Acropora digitifera to a competent Symbiodinium strain: the symbiosome as an arrested early phagosome.

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; <3% of the transcriptome) 4 h after the exposure of Acropora digitifera planulae to a competent strain of Symbiodinium (a clade B strain). This phenomenon has not previously been detected as a consequence of both the lower sensitivity of the microarray approaches used and the sampling times used. The results indicate that complex changes occur, including transient suppression of mitochondrial metabolism and protein synthesis, but are also consistent with the hypothesis that the symbiosome is a phagosome that has undergone early arrest, raising the possibility of common mechanisms in the symbiotic interactions of corals and symbiotic sea anemones with their endosymbionts.

Annotated genes and nonannotated genomes: cross-species use of Gene Ontology in ecology and evolution research.

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.

Molecular epidemiology of clinical Cryptococcus neoformans strains from India.

Little is known about the molecular epidemiology of the human pathogenic fungus Cryptococcus neoformans in India, a country now in the midst of an epidemic of AIDS-related cryptococcosis. We studied 57 clinical isolates from several regions in India, of which 51 were C. neoformans var. grubii, 1 was C. neoformans var. neoformans, and 5 were C. neoformans var. gattii. This strain set included 18 additional sequential isolates from 14 patients. Strains were characterized phenotypically by measuring the polysaccharide capsule and by determining the MICs of standard antifungals. Molecular typing was performed by a PCR-based method using the minisatellite-specific core sequence (M13), by electrophoretic karyotyping, by restriction fragment length polymorphisms with the C. neoformans transposon 1 (TCN-1), and by URA5 DNA sequence analysis. Overall, Indian isolates were less heterogeneous than isolates from other regions and included a subset that clustered into one group based on URA5 DNA sequence analysis. In summary, our results demonstrate (i) differences in genetic diversity of C. neoformans isolates from India compared to isolates from other regions in the world; (ii) that DNA typing with the TCN-1 probe can adequately distinguish C. neoformans var. grubii strains; (iii) that TCN-1 sequences are absent in many C. neoformans var. gattii strains, supporting previous studies indicating that these strains have a limited geographical dispersal; and (iv) that human cryptococcal infection can be associated with microevolution of the infecting strain and by simultaneous coinfection with two distinct C. neoformans strains.

Detection of lateral gene transfer among microbial genomes.

An increasingly comprehensive assessment is being developed of the extent and potential significance of lateral gene transfer among microbial genomes. Genomic sequences can be identified as being of putatively lateral origin by their unexpected phyletic distribution, atypical sequence composition, differential presence or absence in closely related genomes, or incongruent phylogenetic trees. These complementary approaches sometimes yield inconsistent results. Not only more data but also quantitative models and simulations are needed urgently.

On surrogate methods for detecting lateral gene transfer.

Surrogate methods for detecting lateral gene transfer are those that do not require inference of phylogenetic trees. Herein I apply four such methods to identify open reading frames (ORFs) in the genome of Escherichia coli K12 that may have arisen by lateral gene transfer. Only two of these methods detect the same ORFs more frequently than expected by chance, whereas several intersections contain many fewer ORFs than expected. Each of the four methods detects a different non-random set of ORFs. The methods may detect lateral ORFs of different relative ages; testing this hypothesis will require rigorous inference of trees.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2.

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

Protistan parasite QPX of hard-shell clam Mercenaria mercenaria is a member of Labyrinthulomycota.

Biomass of the protistan parasite QPX (quahaug parasite X) of hard-shell clam Mercenaria mercenaria was enriched from in vitro culture. The nuclear gene encoding the 18S RNA of the small-subunit ribosomal (ssu-rDNA) was recovered using the polymerase chain reaction (PCR) and sequenced. Phylogenetic analysis clearly showed that QPX is a member of phylum Labyrinthulomycota, within which it appears as a specific relative of Thraustochytrium pachydermum. These results confirm the provisional assignment of QPX to the Labyrinthulomycota made previously on the basis of morphological and ultrastructural characters found in some, but not all, geographic isolates.

Ichthyophonus irregularis sp. nov. from the yellowtail flounder Limanda ferruginea from the Nova Scotia shelf.

A previously described unusual form of the protistan parasite Ichthyophonus, differing in morphological and developmental features from I. hoferi sensu Plehn & Mulsow, was recovered from yellowtail flounder Limanda ferruginea Storer from the Brown's Bank area of the Nova Scotia shelf. The nuclear gene encoding the rRNA of the small ribosomal subunit was amplified from this unusual form of Ichthyophonus using the polymerase chain reaction, sequenced and aligned with other eukaryote small subunit (ssu)-rDNAs. Inferred phylogenetic trees clearly show that its ssu-rDNA is distinct from those of 2 isolates of I. hoferi sensu Plehn & Mulsow from different hosts and geographical locations (herring in the North Sea, and yellowtail flounder from the Nova Scotia shelf). We consider the unusual form to be a separate species, I. irregularis. The occurrence of a second, distinct type of Ichthyophonus within a single host species raises the possibility that ichthyophoniasis could be produced by different (although related) pathogens, and in some cases, by concurrent infections of the two.

Gene content and organization of a 281-kbp contig from the genome of the extremely thermophilic archaeon, Sulfolobus solfataricus P2.

The sequence of a 281-kbp contig from the crenarchaeote Sulfolobus solfataricus P2 was determined and analysed. Notable features in this region include 29 ribosomal protein genes, 12 tRNA genes (four of which contain archaeal-type introns), operons encoding enzymes of histidine biosynthesis, pyrimidine biosynthesis, and arginine biosynthesis, an ATPase operon, numerous genes for enzymes of lipopolysaccharide biosynthesis, and six insertion sequences. The content and organization of this contig are compared with sequences from crenarchaeotes, euryarchaeotes, bacteria, and eukaryotes.

Completing the sequence of the Sulfolobus solfataricus P2 genome.

The Sulfolobus solfataricus P2 genome collaborators are poised to sequence the entire 3-Mbp genome of this crenarchaeote archaeon. About 80% of the genome has been sequenced to date, with the rest of the sequence being assembled fast. In this publication we introduce the genomic sequencing and automated analysis strategy and present intial data derived from the sequence analysis. After an overview of the general sequence features, metabolic pathway studies are explained, using sugar metabolism as an example. The paper closes with an overview of repetitive elements in S. solfataricus.

Constructing multigenome views of whole microbial genomes.

We have designed and implemented a system to carry out cross-genome comparisons of open reading frames (ORFs) from multiple genomes. This implementation includes a genome profiling system that allows us to explore pairwise comparisons at different levels of match similarity and ask biologically motivated queries involving number and identity of ORFs, their function, functional category, distribution in genomes or in biological domains, and statistics on their matches and match families. This analysis required precise definition of new classification terms and concepts. We define the terms genomic signature, summary signature, biologic domain signature, domain class, match level, match family, and extended match family, then use these terms to define concepts, including genomically universal proteins and proteins characteristics of sets of genomes. We initiate an analysis based on automated FASTA (Pearson, 1996) comparison of 22,419 conceptually translated protein sequences from nine microbial genomes.

Cloning and characterization of a nuclear gene encoding a starch-branching enzyme from the marine red alga Gracilaria gracilis.

The biosynthesis of starch in red algae occurs in the cytosol, in contrast to green plants where it takes place in the plastid. We have cloned a nuclear gene from the red alga Gracilaria gracilis that encodes a homolog of starch-branching enzymes (SBEs); this gene, which is apparently intron-free, was designated as GgSBE1. A potential TATA box, CAAT boxes, and other potential regulatory elements were observed in its 5' flanking region. The encoded 766-aa peptide shares significant sequence similarity with SBEs from green plants (at least 40%), and with glycogen-branching enzymes (GBEs) from human (46%) and Saccharomyces cerevisiae (45%). Southern-hybridization analysis indicates that the gene is single-copy, although weaker signals suggest that related genes exist in the genome of G. gracilis. Phylogenetic analyses indicate that GgSBE1 groups within the eukaryote branching enzymes (BEs) and not with eubacterial GBEs, suggesting that its gene has not been derived directly from an endosymbiotic cyanobacterium, but instead is ancestrally eukaryotic.

Characterization of a galactose-1-phosphate uridylyltransferase gene from the marine red alga Gracilaria gracilis.

The metabolism of D-galactose is a major feature of red-algal physiology. We have cloned and sequenced a gene from the red alga Gracilaria gracilis that encodes a key enzyme of D-galactose metabolism, galactose-1-phosphate uridylyltransferase (GALT). This gene, designated GgGALT1, is apparently devoid of introns. A potential TATA box, four potential CAAT boxes, and a repeated sequence occur in the 5'-flanking region. The predicted 369-aa peptide shares significant sequence similarity with GALTs from other organisms (human, 47%; Saccharomyces cerevisiae, 49%; Solanum tuberosum, 49%). Southern-hybridization analysis reveals two related, but apparently not identical, GALT genes in the nuclear genome of G. gracilis. Sequence analysis indicates that the GgGALT1 enzyme lacks a rubredoxin "knuckle" motif, which in bacterial and fungal GALTs is involved in binding zinc. An open reading frame encoding a potential peptidyl tRNA hydrolase occurs 179 bp downstream from the GgGALT1 gene.

Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America.

Epidermal changes caused by a chytridiomycete fungus (Chytridiomycota; Chytridiales) were found in sick and dead adult anurans collected from montane rain forests in Queensland (Australia) and Panama during mass mortality events associated with significant population declines. We also have found this new disease associated with morbidity and mortality in wild and captive anurans from additional locations in Australia and Central America. This is the first report of parasitism of a vertebrate by a member of the phylum Chytridiomycota. Experimental data support the conclusion that cutaneous chytridiomycosis is a fatal disease of anurans, and we hypothesize that it is the proximate cause of these recent amphibian declines.

Microbial genescapes: phyletic and functional patterns of ORF distribution among prokaryotes.

We have implemented a statistically based approach to comparative genomics that allows us to define and characterize distributional patterns of conceptually translated open reading frames (ORFs) at different confidence levels based on pairwise FASTA matches. In this report, we apply this methodology to nine microbial genomes, focusing particularly on phyletic and functional patterns of ORF distribution within and between the two prokaryotic domains of life, Bacteria and Archaea. We examine patterns of presence and absence of matches, determine the universal ORF set, analyze features of genome specialization between closely related organisms, and present genomic evidence for the monophyly of Archaea. These analyses illustrate how a quantitative approach to comparative genomics can illuminate questions of fundamental biological significance.

Microbial genescapes: a prokaryotic view of the yeast genome.

We examine the translated open reading frames (ORFs) of the yeast Saccharomyces cerevisiae, focusing on those that have FASTA matches in phyletically defined sets of completely sequenced genomes. On this basis, we identify archaeal yeast, bacterial yeast, universal yeast, and yeast ORFs that do not have a match in any of nine prokaryote genomes. Similarly, we examine the yeast mitochondrial genome and the subset of the yeast nuclear ORFs identified as being involved in mitochondrial biogenesis. For the yeast ORFs that match one or more ORFs in these prokaryote genomes, we examine the phyletic and functional distributions of these matches as a function of match strength. These results provide genome level insights into the origin of the eukaryotic cell and the origin of mitochondria. More generally, they exemplify how the growing database of prokaryote genome sequences can help us understand eukaryote genomes.

An archaebacterial homolog of pelota, a meiotic cell division protein in eukaryotes.

An open reading frame (pelA) specifying a homolog of pelota and DOM34, proteins required for meiotic cell division in Drosophila melanogaster and Saccharomyces cerevisiae, respectively, has been cloned, sequenced and identified from the archaebacterium Sulfolobus solfataricus. The S. solfataricus PelA protein is about 20% identical with pelota, DOM34 and the hypothetical protein R74.6 of Caenorhabditis elegans. The presence of a pelota homolog in archaebacteria implies that the meiotic functions of the eukaryotic protein were co-opted from, or added to, other functions existing before the emergence of eukaryotes. The nuclear localization signal and negatively charged carboxy-terminus characteristic of eukaryotic pelota-like proteins are absent from the S. solfataricus homolog, and hence may be indicative of the acquired eukaryotic function(s).

A novel clade of protistan parasites near the animal-fungal divergence.

Sequences of nuclear-encoded small-subunit rRNA genes have been determined for representatives of the enigmatic genera Dermocystidium, Ichthyophonus, and Psorospermium, protistan parasites of fish and crustaceans. The small-subunit rRNA genes from these parasites and from the "rosette agent" (also a parasite of fish) together form a novel, statistically supported clade. Phylogenetic analyses demonstrate this clade to diverge near the animal-fungal dichotomy, although more precise resolution is problematic. In the most parsimonious and maximally likely phylogenetic frameworks inferred from the most stably aligned sequence regions, the clade constitutes the most basal branch of the metazoa; but within a limited range of model parameters, and in some analyses that incorporate less well-aligned sequence regions, an alternative topology in which it diverges immediately before the animal-fungal dichotomy was recovered. Mitochondrial cristae of Dermocystidium spp. are flat, whereas those of Ichthyophonus hoferi appear tubulovesiculate. These results extend our understanding of the types of organisms from which metazoa and fungi may have evolved.

Organizational characteristics and information content of an archaeal genome: 156 kb of sequence from Sulfolobus solfataricus P2.

We have initiated a project to sequence the 3 Mbp genome of the thermoacidophilic archaebacterium Sulfolobus solfataricus P2. Cosmids were selected from a provisional set of minimally overlapping clones, subcloned in pUC18, and sequenced using a hybrid (random plus directed) strategy to give two blocks of contiguous unique sequence, respectively, 100,389 and 56,105 bp. These two contigs contain a total of 163 open reading frames (ORFs) in 26-29 putative operons; 56 ORFs could be identified with reasonable certainty. Clusters of ORFs potentially encode proteins of glycogen biosynthesis, oxidative decarboxylation of pyruvate, ATP-dependent transport across membranes, isoprenoid biosynthesis, protein synthesis, and ribosomes. Putative promoters occur upstream of most ORFs. Thirty per cent of the predicted strong and medium-strength promoters can initiate transcription at the start codon or within 10 nucleotides upstream, indicating a process of initial mRNA-ribosome contact unlike that of most eubacterial genes. A novel termination motif is proposed to account for 15 additional terminations. The two contigs differ in densities of ORFs, insertion elements and repeated sequences; together they contain two copies of the previously reported insertion sequence ISC 1217, five additional IS elements representing four novel types, four classes of long non-IS repeated sequences, and numerous short, perfect repeats.