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.

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.

Characterization of the polyubiquitin gene in the marine red alga Gracilaria verrucosa.

We have cloned a nuclear gene (UBI6R) and corresponding cDNAs that encode polyubiquitin in the florideophycidean red alga Gracilaria verrucosa. The gene encodes a polyubiquitin composed of six tandem ubiquitin units, followed by a single glutamine residue. The deduced amino acid sequences are identical among all six units, and identical to the ubiquitin of the florideophyte Aglaothamnion neglectum. There is high sequence similarity among the red algal ubiquitins and those of animals, green plants, fungi and several protists. Only one polyubiquitin gene was found by Southern hybridization analysis of G. verrucosa nuclear DNA. The upstream region of the gene is rich in putative cis-acting transcription-regulatory elements, including a putative heat-responsive element. Poly(A) addition to UBI6R mRNA was observed in cDNAs at four different sites, implicating the sequences AATAAA and (or) AGTAAA as poly(A) addition signals. The polyubiquitin genes of red algae show features of concerted evolution, but appear to be subject to less sequence homogenization than those of animals.

The Sulfolobus solfataricus P2 genome project.

Over 800 kbp of the 3-Mbp genome of Sulfolobus solfataricus have been sequenced to date. Our approach is to sequence subclones of mapped cosmids, followed by sequencing directly on cosmid templates with custom primers. Using a prototype automated system for genome-scale analysis, known as MAGPIE, along with other tools, we have discovered one open reading frame of at least 100 amino acids per kbp of sequence, and have been able to associate 50% of these with known genes through database searches. An examination of completely sequenced cosmids suggests a clustering of genes by function in the S. solfataricus genome.

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.

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).

Ribosomal RNA and the major lines of evolution: a perspective.

Does the "universal tree" based on small-subunit ribosomal RNA sequences show the phylogenetic relationship of all modern organisms? The answer is "yes" only if all these rRNAs are orthologous. Herein I argue that the major rRNA lineages (e.g. eubacterial, one or more archaebacterial and eukaryotic nucleocytoplasmic) probably arose from a divergent population of rRNAs in the progenote, antedating the universal common ancestral organism. Thus the major lineages of rRNA are probably not orthologous, but paralogous. The extrapolated date for the origin of the common ancestral small-subunit rRNA (3.6-4.7 x 10(9) years ago) is consistent with major rRNA lineages being paralogous. This perspective on the early evolution of genes and organisms rationalizes the presence of unexpected ribosomal characters in microsporidia, and bears on xenogenous and endogenous theories of the origin of the organelles in eukaryotes.

Matrix representation in reconstructing phylogenetic relationships among the eukaryotes.

The application of the method of matrix representation (MRP) addresses several problems associated with reconstructing phylogenetic relationships among composite organisms such as the eukaryotes. By allowing multiple molecular-sequence and non-molecular data sets to be combined into a single, relatively compact matrix for joint analysis, MRP addresses problems associated with the depth, diversification and specialization of eukaryotic lineages. Characters whose distributions may have been affected by endosymbiotic lateral transfer can be identified by compatibility analysis of MRP-generated hybrid matrices. In conjunction with variant-filtering techniques, MRP can be used to map patterns of lateral transfer in a given phylogenetic tree. These applications are illustrated with molecular-sequence and non-molecular data sets for eukaryotes.

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.

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.

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.

Phylogenetic inference based on matrix representation of trees.

Rooted phylogenetic trees can be represented as matrices in which the rows correspond to termini, and columns correspond to internal nodes (elements of the n-tree). Parsimony analysis of such a matrix will fully recover the topology of the original tree. The maximum size of the represented matrix depends only on the number of termini in the tree; for a tree derived from molecular sequences, the represented matrix may be orders of magnitude smaller than the original data matrix. Representations of multiple trees (which may or may not have identical termini) can readily be combined into a single matrix; columns of discrete-character-state data can be added and, if desired, weighted differentially. Parsimony analysis of the resulting composite matrix yields a hybrid supertree which typically provides greater resolution than conventional consensus trees. Use of this method is illustrated with examples involving multiple tRNA genes in organelles and multiple protein-coding genes in eukaryotes.

Variant forms of a group I intron in nuclear small-subunit rRNA genes of the marine red alga Porphyra spiralis var. amplifolia.

A group IC1 intron occurs in nuclear small-subunit (18S) ribosomal RNA (SSU rRNA) genes of the marine red alga Porphyra spiralis var. amplifolia. This intron occurs at the same position as the self-splicing group IC1 introns in nuclear SSU rDNAs of the fungus Pneumocystis carinii and in the green alga Chlorella ellipsoidea and shares sequence identity with the Pneumocystis carinii intron in domains L1, P1, P2, and L2, outside the conserved core. Three size variants, differing in amount of sequence in L1, exist and are differentially distributed in geographically distinct populations. Preliminary data suggest that the largest variant can self-splice in vitro. Short open reading frames are present but do not correspond to known genes. Repeated nucleotide motifs, reminiscent of duplicated target sites of transposons or Alu elements, are associated with the intron and with one of the variant forms of L1. Insertions are present in nuclear SSU rDNAs of several other Porphyra species and of the red alga Bangia atropurpurea; insertionless rDNA variants also occur in several Porphyra species. Our observations are most readily explained by intron mobility, although it remains unclear how transfer could have been mediated between genomes of organisms as ecologically diverse as marine red algae, freshwater green algae, and a mammalian-pathogenic fungus.

Characterization of the nuclear gene encoding mitochondrial aconitase in the marine red alga Gracilaria verrucosa.

We have cloned a nuclear gene from the marine red alga Gracilaria verrucosa that encodes the complete 779 amino-acid mitochondrial aconitase (m-ACN), the first characterized from a photosynthetic organism. The N-terminal 28 deduced amino acids are predicted to constitute the mitochondrial transit peptide, the first described from a red alga. Putative transcriptional cis-acting elements were identified in the upstream untranslated region. The G. verrucosa m-ACN gene (m-ACN) is present in a single copy and is located ca. 1.5 kb upstream from the single-copy polyubiquitin gene. The single spliceosomal intron is located near the 5' end of the region encoding the mature m-ACN in precisely the same location and phase as intron 2 in Caenorhabditis elegans m-ACN; sequences at its 3' and 5' splice junctions and at the predicted lariat branch point conform well to the eukaryote consensus sequences. Multiple protein-sequence alignment of m-ACN, bacterial aconitase (b-ACN) and iron-responsive element-binding protein (IRE-BP), and phylogenetic analyses, revealed that m-ACN does not share a recent common ancestry with either b-ACN or IRE-BP.

cDNA cloning and characterization of the nuclear gene encoding chloroplast glyceraldehyde-3-phosphate dehydrogenase from the marine red alga Gracilaria verrucosa.

Using a PCR-generated homologous probe, we have recovered a cDNA (GapA cDNA) encoding the complete 338 amino-acid chloroplast GAPDH of the marine red alga Gracilaria verrucosa, together with its 78 amino-acid transit peptide. This cDNA was readily aligned with chloroplast-localized GAPDH genes (GapA and GapB) of green plants. The proline residue which contributes to the specificity of NAD+ binding to cytosolic GAPDHs is absent from the deduced polypeptide chain of G. verrucosa GapA as is also the case in the chloroplast GAPDHs of plants. The transit peptide shows a high proportion of random coil, an amino-terminal Met-Ala dipeptide, a high content of hydroxylamino acids, and a net positive charge. The polyadenylation signal appears to be AGTAAA. Genomic Southern-hybridization data indicate that only one chloroplast-GAPDH gene may occur in G. verrucosa. Bootstrapped parsimony trees indicate that the G. verrucosa GapA gene is a sister group to plant chloroplast-GAPDH genes, and are most readily interpreted as showing that red algal and plant chloroplast-localized GAPDHs arose in a single endosymbiotic event.

Cloning and characterization of the nuclear gene and cDNAs for triosephosphate isomerase of the marine red alga Gracilaria verrucosa.

cDNAs and an intronless single-copy nuclear gene (TPI1) encoding triosephosphate isomerase have been cloned and sequenced from the marine red alga Gracilaria verrucosa. The predicted amino-acid sequence of TPI1 is readily alignable with those of other known TPIs; 26 of 27 active-site residues and 19 of 26 intersubunit-contact residues are identical between TPIs of G. verrucosa and/or animals and green plants. A partial cDNA sequence of a second TPI gene (TPI2), presumably encoding plastid-localized TPI, was recovered by PCR and demonstrated by phylogenetic analysis to be red algal; no TP12 cDNA or genomic clones could be recovered. Genomic Southern analysis demonstrated that at least two TPI-like genes are present in the nuclear DNA of G. verrucosa.

The nuclear gene and cDNAs encoding cytosolic glyceraldehyde-3-phosphate dehydrogenase from the marine red alga Gracilaria verrucosa: cloning, characterization and phylogenetic analysis.

We have cloned and sequenced the single-copy nuclear gene (GapC) encoding the complete 335-amino acid cytosolic glyceraldehyde-3-phosphate dehydrogenase (GAPC) from the red alga Gracilaria verrucosa. The proline residue which contributes to the specificity of NAD+ binding in other GAPC-like proteins is present. Putative regulatory regions, including GC-rich regions, a GATA element, and 11-base T- and T/G-clusters, but excluding TATA- and CCAAT-boxes, were identified upstream. Two types of GapC cDNAs differing in polyadenylation site were characterized. An 80-bp phase-two spliceosomal intron was identified in a novel position interrupting the highly conserved cofactor-coding region I. The G. verrucosa GAPC was easily aligned with other known GAPC-type sequences. Inferred phylogenetic trees place red algae among the eukaryote crown taxa, although with modest bootstrap support and without stable resolution among related GAPC lineages.

Paired-ion reversed-phase high-performance liquid chromatography of phenol sulfates in synthetic mixtures, algal extracts and urine.

Phenol sulfate esters have been analyzed by paired-ion reversed-phase high-performance liquid chromatography. The method provided direct, rapid chromatography of phenol sulfates in crude extracts of the red alga (Polysiphonia lanosa (2,3-dibromo-4,5-dihydroxybenzyl alcohol 1',4-disulfate), of the brown alga Ascophyllum nodosum (1,2,3,5-tetrahydroxybenzene 2,5-disulfate), and in rat urine (resorcinol mono- and disulfates). Detector response (254 nm) was linear within the approximate range from 30--125 ng to 5--10 microgram. Semipreparative scale chromatography provided sufficient amounts of purified phenol sulfates for further analysis by paper electrophoresis.