Approaches and Strategy for Cancer Research and Surveillance Data: Integration, Information Pipeline, Data Models, and Informatics Opportunities

Thesis (Ph.D.)--University of Washington, 2016-04

Partitioned likelihood support and the evaluation of data set conflict

In simultaneous analyses of multiple data partitions, the trees relevant when measuring support for a clade are the optimal tree, and the best tree lacking the clade (i.e., the most reasonable alternative). The parsimony-based method of partitioned branch support (PBS) forces each data set to arbitrate between the two relevant trees. This value is the amount each data set contributes to clade support in the combined analysis, and can be very different to support apparent in separate analyses. The approach used in PBS can also be employed in likelihood: a simultaneous analysis of all data retrieves the maximum likelihood tree, and the best tree without the clade of interest is also found. Each data set is fitted to the two trees and the log-likelihood difference calculated, giving partitioned likelihood support (PLS) for each data set. These calculations can be performed regardless of the complexity of the ML model adopted. The significance of PLS can be evaluated using a variety of resampling methods, such as the Kishino-Hasegawa test, the Shimodiara-Hasegawa test, or likelihood weights, although the appropriateness and assumptions of these tests remains debated.

Visualization of molecular fingerprints

A visualization plot of a data set of molecular data is a useful tool for gaining insight into a set of molecules. In chemoinformatics, most visualization plots are of molecular descriptors, and the statistical model most often used to produce a visualization is principal component analysis (PCA). This paper takes PCA, together with four other statistical models (NeuroScale, GTM, LTM, and LTM-LIN), and evaluates their ability to produce clustering in visualizations not of molecular descriptors but of molecular fingerprints. Two different tasks are addressed: understanding structural information (particularly combinatorial libraries) and relating structure to activity. The quality of the visualizations is compared both subjectively (by visual inspection) and objectively (with global distance comparisons and local k-nearest-neighbor predictors). On the data sets used to evaluate clustering by structure, LTM is found to perform significantly better than the other models. In particular, the clusters in LTM visualization space are consistent with the relationships between the core scaffolds that define the combinatorial sublibraries. On the data sets used to evaluate clustering by activity, LTM again gives the best performance but by a smaller margin. The results of this paper demonstrate the value of using both a nonlinear projection map and a Bernoulli noise model for modeling binary data.

Genomic variation of HIV-1 and its effects on drug resistance and molecular epidemiology analyses

In the first part of this study human immunodeficiency virus type 1 (HIV-1) proviral DNA sequences derived from 201 clones of the C2-V3 env region and the first exon of the tat gene were obtained from six MV-1 infected heterosexual couples. These molecular data were used to confirm the epidemiological relationships. The ability of the molecular data to draw such conclusions was also tested with multiple phylogenetic analyses. The tat region was much more useful in establishing epidemiological relationships than the commonly used C2-V3.^ Subsequently, using nucleotide sequences from the first exon of the Tat gene, we tested the hypothesis that a Florida dentist (a common source) infected five of his patients in the course of dental procedures, against the null hypothesis that the dentist and each individual of the dental group independently acquired the virus within the local community. Multiple phylogenetic analyses demonstrated that the sequences of the five patients were significantly more related to each other than to sequences of the controls. Our results using Tat sequences, combined with envelope sequence data, strongly support a common phylogenetic epidemiological relationship among these five patients.^ A third study is presented, which deals with the effects of genomic variations in drug resistance. HIV-1 reverse transcriptase (RT) mutations were detected in DNA from peripheral blood mononuclear cells from 11 of 12 HIV-infected children after 11-20 months of zidovudine monotherapy. The codon 41/215 mutant combination was associated with general decline in health status. Patients developing the codon 70 mutation tended to have a better health status. ^

Molecular characterization of complete and incomplete immunoglobulin heavy chain gene rearrangements in hairy cell leukemia.

PURPOSE: We analyzed patients with hairy cell leukemia (HCL) to achieve a better understanding of the differentiation stage reached by HCL cells and to define the key role of the diversification of cell surface makers, especially CD25 expression. PATIENTS AND METHODS: We analyzed 38 previously untreated patients with HCL to characterize their complete (VDJ(H)) and incomplete (DJ(H)) immunoglobulin (Ig) heavy chain (IgH) rearrangements, including somatic hypermutation pattern and gene segment use. RESULTS: A correlation between immunophenotypic profile and molecular data was seen. All 38 cases showed monoclonal amplifications: VDJ(H) in 97%, DJ(H) in 42%, and both in 39%. Segments from the D(H)3 family were used more in complete compared with incomplete rearrangements (45% vs. 12%; P <.005). Furthermore, comparison between molecular and immunophenotypic characteristics disclosed differences in the expression of CD25 antigen; CD25(-) cases, a phenotype associated with HCL variant, showed complete homology to the germline in 3 of 5 cases (60%), whereas this characteristic was never observed in CD25(+) cases (P <.005). Moreover, V(H)4-34, V(H)1-08, and J(H)3 segments appeared in 2, 1, and 2 CD25(-) cases, respectively, whereas they were absent in all CD25(+) cases. CONCLUSION: These results support that HCL is a heterogeneous entity including subgroups with different molecular characteristics, which reinforces the need for additional studies with a larger number of patients to clarify the real role of gene rearrangements in HCL.

Distal 22q11.2 Microduplication Combined With Typical 22q11.2 Proximal Deletion: A Case Report.

The 22q11 chromosomal region contains low copy repeats (LCRs) sequences that mediate non-allelic homologous recombination, which predisposes to copy number variations (CNVs) at this locus. Hemizygous deletions of the proximal 22q11.2 region result in the 22q11.2 deletion syndrome (22q11.2 DS). In addition, 22q11.2 duplications involving the distal LCR22s have been reported. This article describes a patient presenting a 2.5-Mb de novo deletion at proximal 22q11.21 region (between LCRs A-D), combined with a 1.3-Mb maternally inherited duplication at distal 22q11.23 region (between LCRs F-H). The presence of concomitant chromosomal imbalances found in this patient has not been reported previously. Clinical and molecular data were compared with literature, in order to contribute to genotype-phenotype correlation. These findings exemplify the complexity and genetic heterogeneity observed in 22q11.2 deletion syndrome and highlights the difficulty to make genetic counseling and predict phenotypic consequences in these situations.

The Sperm Of Hexagenia (pseudeatonica) Albivitta Walker (ephemeroptera: Fossoriae: Ephemeridae)

This study describes the sperm morphology of the mayfly Hexagenia (Pseudeatonica) albivitta (Ephemeroptera). Its spermatozoon measures approximately 30 μm of which 9 μm corresponds to the head. The head is composed of an approximately round acrosomal vesicle and a cylindrical nucleus. The nucleus has two concavities, one in the anterior tip, where the acrosomal vesicle is inserted and a deeper one at its base, where the flagellum components are inserted. The flagellum is composed of an axoneme, a mitochondrion and a dense rod adjacent to the mitochondrion. A centriolar adjunct is also observed surrounding the axoneme in the initial portion of the flagellum and extends along the flagellum for at least 2 μm, surrounding the axoneme in a half-moon shape. The axoneme is the longest component of the flagellum, and it follows the 9+9+0 pattern, with no central pair of microtubules. At the posterior region of the flagellum, the mitochondrion has a dumb-bell shape in cross sections that, together with the rectangular mitochondrial-associated rod, is responsible for the flattened shape of the flagellum. An internal membrane is observed surrounding both mitochondrion and its associated structure.

Resurrection of Anopheles goeldii from synonymy with Anopheles nuneztovari (Diptera, Culicidae) and a new record for Anopheles dunhami in the Brazilian Amazon

Nucleotide sequences of the internal transcribed spacer 2 (ITS2) rDNA and partial sequences of the cytochrome coxidase subunit I (COI) mtDNA and white gene nDNA were obtained from specimens of Anopheles nuneztovari A collected in Macapá (state of Amapá), Óbidos, Prainha and Almeirim (state of Pará), Itacoatiara and Parintins (state of Amazonas), Brazil, and compared with previously published sequences of A. nuneztovari s.l. Results of the Bayesian phylogenetic analyses performed using either COI or combined ITS2, COI and white gene sequences suggest that An. nuneztovari B/C is distinct from specimens obtained in the Amazonas/Solimões River basin. Anopheles goeldii, currently in synonymy with An. nuneztovari, was described from individuals collected in Belterra (= Fordlândia) in the Tapajós River, state of Pará, Southern Amazonas River. Morphological comparisons of the characteristics of the male genitalia indicated that An. nuneztovari A and An. goeldii are similar but distinct from An. nuneztovariB/C by the apex of the aedeagus. In considering the results of the phylogenetic analyses and morphological comparisons, An. goeldii is resurrected from synonymy with An. nuneztovari. Additionally, Anopheles dunhamiis reported for the first time in Parintins. This species can be distinguished from An. goeldiiby characters of the male genitalia and molecular data

Phylogenetic relationships within the speciose family Characidae (Teleostei: Ostariophysi: Characiformes) based on multilocus analysis and extensive ingroup sampling

Background: With nearly 1,100 species, the fish family Characidae represents more than half of the species of Characiformes, and is a key component of Neotropical freshwater ecosystems. The composition, phylogeny, and classification of Characidae is currently uncertain, despite significant efforts based on analysis of morphological and molecular data. No consensus about the monophyly of this group or its position within the order Characiformes has been reached, challenged by the fact that many key studies to date have non-overlapping taxonomic representation and focus only on subsets of this diversity. Results: In the present study we propose a new definition of the family Characidae and a hypothesis of relationships for the Characiformes based on phylogenetic analysis of DNA sequences of two mitochondrial and three nuclear genes (4,680 base pairs). The sequences were obtained from 211 samples representing 166 genera distributed among all 18 recognized families in the order Characiformes, all 14 recognized subfamilies in the Characidae, plus 56 of the genera so far considered incertae sedis in the Characidae. The phylogeny obtained is robust, with most lineages significantly supported by posterior probabilities in Bayesian analysis, and high bootstrap values from maximum likelihood and parsimony analyses. Conclusion: A monophyletic assemblage strongly supported in all our phylogenetic analysis is herein defined as the Characidae and includes the characiform species lacking a supraorbital bone and with a derived position of the emergence of the hyoid artery from the anterior ceratohyal. To recognize this and several other monophyletic groups within characiforms we propose changes in the limits of several families to facilitate future studies in the Characiformes and particularly the Characidae. This work presents a new phylogenetic framework for a speciose and morphologically diverse group of freshwater fishes of significant ecological and evolutionary importance across the Neotropics and portions of Africa.

Genetic variability of the Amazon River prawn Macrobrachium amazonicum (Decapoda, Caridea, Palaemonidae)

The freshwater prawn Macrobrachium amazonicum is widely distributed in South America, and occupies habitats with a wide range of salinities. Several investigations have revealed the existence of wide intraspecific variability among different populations, although the understanding of this variability is still fragmentary and incomplete. We compared and characterized inland and coastal populations of M. amazonicum from Brazil, using molecular data (16S and COI mtDNA) to describe the degree of variability, structure, and relationships among them. Genetic divergence rates among populations showed variability at the intraspecific level. All the analyses evidenced significant genetic divergence among populations, structuring them in three groups: I-inland waters of the Amazonian Hydrographic Region (HR); II-Parana/Paraguay HR; and III-coastal systems of northern and northeastern Brazil. Phylogenetic reconstructions revealed that the populations form a single monophyletic clade, which supports their characterization as a single species. Clade I was a sister clade of that formed by clades II and III, which were themselves sister clades. Populations from Sertaozinho/Miguelopolis and Avare, introduced into the state of Sao Paulo, may have originated from natural populations in the states of Mato Grosso do Sul and Para, respectively. Geographical isolation probably contributed to the observed variation, and if this isolation continues. M. amazonicum may undergo speciation within its broad geographical distribution. The sequences obtained here can be used as name-tags for population identification, and the DNA barcodes are useful to identify the origin of specimens used in different freshwater-prawn cultures or introduced populations of unknown origin.

Description and Phylogenetic Relationships of a New Genus and Species of Lizard (Squamata, Gymnophthalmidae) from the Amazonian Rainforest of Northern Brazil

Marinussaurus curupira, a new genus and species of Gymnophthalmidae lizard is described from Iranduba, state of Amazonas, Brazil. The genus is characterized by an elongate body; short and stout pentadactyl limbs; all digits clawed; single frontonasal; two prefrontals; absence of frontoparietals; interparietal and parietals forming a straight posterior margin, with interparietal shorter than parietals; distinctive ear opening and eyelid; few temporals; three pairs of chin shields; nasal divided; a distinct collar; smooth, mainly hexagonal, dorsal scales; smooth quadrangular ventral scales; two precloacal and three femoral pores on each side in males; pores between three or four scales. Parsimony (PAR) and partitioned Bayesian (BA) phylogenetic analyses with morphological and molecular data recovered the new genus as a member of the Ecpleopodini radiation of the Cercosaurinae. A close relationship of the new genus with Arthrosaura is postulated.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Molecules Clarify a Cnidarian Life Cycle - The ""Hydrozoan'' Microhydrula limopsicola Is an Early Life Stage of the Staurozoan Haliclystus antarcticus

Background: Life cycles of medusozoan cnidarians vary widely, and have been difficult to document, especially in the most recently proposed class Staurozoa. However, molecular data can be a useful tool to elucidate medusozoan life cycles by tying together different life history stages. Methodology/Principal Findings: Genetic data from fast-evolving molecular markers (mitochondrial 16S, nuclear ITS1, and nuclear ITS2) show that animals that were presumed to be a hydrozoan, Microhydrula limopsicola (Limnomedusae, Microhydrulidae), are actually an early stage of the life cycle of the staurozoan Haliclystus antarcticus (Stauromedusae, Lucernariidae). Conclusions/Significance: Similarity between the haplotypes of three markers of Microhydrula limopsicola and Haliclystus antarcticus settles the identity of these taxa, expanding our understanding of the staurozoan life cycle, which was thought to be more straightforward and simple. A synthetic discussion of prior observations makes sense of the morphological, histological and behavioral similarities/congruence between Microhydrula and Haliclystus. The consequences are likely to be replicated in other medusozoan groups. For instance we hypothesize that other species of Microhydrulidae are likely to represent life stages of other species of Staurozoa.45

Anatomical development of the pericarp and seed of Oncidium flexuosum Sims (ORCHIDACEAE)

Interpretation of the anatomical structure of the ovary and fruit of the Orchidaceae family is still controversial, which makes it difficult to understand the development and dehiscence of the fruit. The genus Oncidium is polyphyletic and is currently the subject of taxonomic studies. In this study, we have investigated the anatomical development of the pericarp and seed of Oncidium flexuosum Sims to determine important diagnostic characters that, along with molecular data, can assist in defining this group. We have found a new anatomical characteristic of the family: the presence of precursor cells for fruit dehiscence, which were visible from the beginning of development and located on the outer walls of the sterile valves. In contrast with what has been observed by different authors with other species, in the mature fruit of O. flexuosum, only the endocarp of the fertile valves and a few cells near the exocarp and the vascular bundle in the sterile valves show parietal thickening, while the rest remains parenchymatous. During the development of the ovule and embryo, we have shown that the embryonic sac of this species has eight nuclei and that the embryo has a long and elaborate suspensor. (C) 2011 Elsevier GmbH. All rights reserved.

Eudesmanolides and Methyl Ester Derivatives from Dimerostemma arnottii

The phytochemical investigation of Dimerostemma arnottii (Asteraceae) afforded, in addition to a known eudesmanolide, two unusual eudesmane methyl ester derivatives and a new eudesmanolide. Structural elucidation of the compounds was based on their ID and 2D NMR spectroscopic as well as HR-ESI-MS data. There is a remarkable similarity between the structures of the eudesmanes from D. arnottii and those previously encountered in other Dimerostemma species, which is in agreement with the results of a previous phylogenetic study based on molecular data. The chemotaxonomic relevance of the isolated compounds is briefly discussed.