USE OF HIDDEN MARKOV MODELS FOR QTL MAPPING

An important aspect of the QTL mapping problem is the treatment of missing genotype data. If complete genotype data were available, QTL mapping would reduce to the problem of model selection in linear regression. However, in the consideration of loci in the intervals between the available genetic markers, genotype data is inherently missing. Even at the typed genetic markers, genotype data is seldom complete, as a result of failures in the genotyping assays or for the sake of economy (for example, in the case of selective genotyping, where only individuals with extreme phenotypes are genotyped). We discuss the use of algorithms developed for hidden Markov models (HMMs) to deal with the missing genotype data problem.

Bayesian Hidden Markov Modeling of Array CGH Data

Genomic alterations have been linked to the development and progression of cancer. The technique of Comparative Genomic Hybridization (CGH) yields data consisting of fluorescence intensity ratios of test and reference DNA samples. The intensity ratios provide information about the number of copies in DNA. Practical issues such as the contamination of tumor cells in tissue specimens and normalization errors necessitate the use of statistics for learning about the genomic alterations from array-CGH data. As increasing amounts of array CGH data become available, there is a growing need for automated algorithms for characterizing genomic profiles. Specifically, there is a need for algorithms that can identify gains and losses in the number of copies based on statistical considerations, rather than merely detect trends in the data. We adopt a Bayesian approach, relying on the hidden Markov model to account for the inherent dependence in the intensity ratios. Posterior inferences are made about gains and losses in copy number. Localized amplifications (associated with oncogene mutations) and deletions (associated with mutations of tumor suppressors) are identified using posterior probabilities. Global trends such as extended regions of altered copy number are detected. Since the posterior distribution is analytically intractable, we implement a Metropolis-within-Gibbs algorithm for efficient simulation-based inference. Publicly available data on pancreatic adenocarcinoma, glioblastoma multiforme and breast cancer are analyzed, and comparisons are made with some widely-used algorithms to illustrate the reliability and success of the technique.

A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

A rare anomaly of the anterior communicating artery complex hidden by a large broad-neck aneurysm and disclosed by three-dimensional rotational angiography

Double fenestration of the anterior communicating artery (ACoA) complex associated with an aneurysm is a very rare finding and is usually caused by ACoA duplication and the presence of a median artery of the corpus callosum (MACC). We present a patient in whom double fenestration was not associated with ACoA duplication or even with MACC, representing therefore, a previously unreported anatomic variation. A 43 year old woman experienced sudden headache and the CT scans showed subarachnoid haemorrhage (SAH). On admission, her clinical condition was consistent with Hunt and Hess grade II. Conventional digital subtraction angiography (DSA) was performed and revealed multiple intracranial aneurysms arising from both middle cerebral arteries (MCA) and from the ACoA. Three-dimensional rotational angiography (3D-RA) disclosed a double fenestration of the ACoA complex which was missed by DSA. The patient underwent a classic pterional approach in order to achieve occlusion of both left MCA and ACoA aneurysms by surgical clipping. The post-operative period was uneventful. A rare anatomical variation characterised by a double fenestration not associated with ACoA duplication or MACC is described. The DSA images missed the double fenestration which was disclosed by 3D-RA, indicating the importance of 3D-RA in the diagnosis and surgical planning of intracranial aneurysms.

Light Transport in PT-invariant Photonic Structures with Hidden Symmetries

We introduce a recursive bosonic quantization technique for generating classical PT photonic structures that possess hidden symmetries and higher order exceptional points. We study light transport in these geometries and we demonstrate that perfect state transfer is possible only for certain initial conditions. Moreover, we show that for the same propagation direction, left and right coherent transports are not symmetric with field amplitudes following two different trajectories. A general scheme for identifying the conservation laws in such PT-symmetric photonic networks is also presented.

Early and Middle Holocene Hunter-Gatherer Occupations in Western Amazonia: The Hidden Shell Middens

We report on previously unknown early archaeological sites in the Bolivian lowlands, demonstrating for the first time early and middle Holocene human presence in western Amazonia. Multidisciplinary research in forest islands situated in seasonally-inundated savannahs has revealed stratified shell middens produced by human foragers as early as 10,000 years ago, making them the oldest archaeological sites in the region. The absence of stone resources and partial burial by recent alluvial sediments has meant that these kinds of deposits have, until now, remained unidentified. We conducted core sampling, archaeological excavations and an interdisciplinary study of the stratigraphy and recovered materials from three shell midden mounds. Based on multiple lines of evidence, including radiocarbon dating, sedimentary proxies (elements, steroids and black carbon), micromorphology and faunal analysis, we demonstrate the anthropogenic origin and antiquity of these sites. In a tropical and geomorphologically active landscape often considered challenging both for early human occupation and for the preservation of hunter-gatherer sites, the newly discovered shell middens provide evidence for early to middle Holocene occupation and illustrate the potential for identifying and interpreting early open-air archaeological sites in western Amazonia. The existence of early hunter-gatherer sites in the Bolivian lowlands sheds new light on the region’s past and offers a new context within which the late Holocene “Earthmovers” of the Llanos de Moxos could have emerged.