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.

The alpha and beta subunits of the metalloprotease meprin are expressed in separate layers of human epidermis, revealing different functions in keratinocyte proliferation and differentiation

The zinc endopeptidase meprin (EC 3.4.24.18) is expressed in brush border membranes of intestine and kidney tubules, intestinal leukocytes, and certain cancer cells, suggesting a role in epithelial differentiation and cell migration. Here we show by RT-PCR and immunoblotting that meprin is also expressed in human skin. As visualized by immunohistochemistry, the two meprin subunits are localized in separate cell layers of the human epidermis. Meprin alpha is expressed in the stratum basale, whereas meprin beta is found in cells of the stratum granulosum just beneath the stratum corneum. In hyperproliferative epidermis such as in psoriasis vulgaris, meprin alpha showed a marked shift of expression from the basal to the uppermost layers of the epidermis. The expression patterns suggest distinct functions for the two subunits in skin. This assumption is supported by diverse effects of recombinant meprin alpha and beta on human adult low-calcium high-temperature keratinocytes. Here, beta induced a dramatic change in cell morphology and reduced the cell number, indicating a function in terminal differentiation, whereas meprin alpha did not affect cell viability, and may play a role in basal keratinocyte proliferation.

Messenger RNA levels and binding sites of muscarinic acetylcholine receptors in gastrointestinal muscle layers from healthy dairy cows

Acetylcholine interacts with muscarinic receptors (M) to mediate gastrointestinal (GI) smooth muscle contractions. We have compared mRNA levels and binding sites of M(1)to M(5) in muscle tissues from fundus abomasi, pylorus, ileum, cecum, proximal loop of the ascending colon (PLAC), and external loop of the spiral colon (ELSC) of healthy dairy cows. The mRNA levels were measured by quantitative RT-PCR. The inhibition of [(3)H]-QNB (1-quinuclidinyl-[phenyl-4-(3)H]-benzilate) binding by M antagonists [atropine (M(1 - 5)), pirenzepine (M(1)), methoctramine (M(2)), 4-DAMP (M(3)), and tropicamide (M(4))] was used to identify receptors at the functional level. Maximal binding (B(max)) was determined through saturation binding with atropine as a competitor. The mRNA levels of M(1), M(2), M(3), and M(5) represented 0.2, 48, 50, and 1.8%, respectively, of the total M population, whereas mRNA of M(4) was undetectable. The mRNA levels of M(2) and of M(3) in the ileum were lower (P < 0.05) than in other GI locations, which were similar among each other. Atropine, pirenzepine, methoctramine, and 4-DAMP inhibited [(3)H]-QNB binding according to an either low- or high-affinity receptor pattern, whereas tropicamide had no effect on [(3)H]-QNB binding. The [(3)H]-QNB binding was dose-dependent and saturable. B(max) in fundus, pylorus, and PLAC was lower (P < 0.05) than in the ELSC, and in the pylorus lower (P < 0.05) than in the ileum. B(max) and mRNA levels were negatively correlated (r = -0.3; P < 0.05). In conclusion, densities of M are different among GI locations, suggesting variable importance of M for digestive functions along the GI tract.

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.

Water transport in complex, non-wetting porous layers with applications to water management in low temperature fuel cells

An experimental setup was designed to visualize water percolation inside the porous transport layer, PTL, of proton exchange membrane, PEM, fuel cells and identify the relevant characterization parameters. In parallel with the observation of the water movement, the injection pressure (pressure required to transport water through the PTL) was measured. A new scaling for the drainage in porous media has been proposed based on the ratio between the input and the dissipated energies during percolation. A proportional dependency was obtained between the energy ratio and a non-dimensional time and this relationship is not dependent on the flow regime; stable displacement or capillary fingering. Experimental results show that for different PTL samples (from different manufacturers) the proportionality is different. The identification of this proportionality allows a unique characterization of PTLs with respect to water transport. This scaling has relevance in porous media flows ranging far beyond fuel cells. In parallel with the experimental analysis, a two-dimensional numerical model was developed in order to simulate the phenomena observed in the experiments. The stochastic nature of the pore size distribution, the role of the PTL wettability and morphology properties on the water transport were analyzed. The effect of a second porous layer placed between the porous transport layer and the catalyst layer called microporous layer, MPL, was also studied. It was found that the presence of the MPL significantly reduced the water content on the PTL by enhancing fingering formation. Moreover, the presence of small defects (cracks) within the MPL was shown to enhance water management. Finally, a corroboration of the numerical simulation was carried out. A threedimensional version of the network model was developed mimicking the experimental conditions. The morphology and wettability of the PTL are tuned to the experiment data by using the new energy scaling of drainage in porous media. Once the fit between numerical and experimental data is obtained, the computational PTL structure can be used in different types of simulations where the conditions are representative of the fuel cell operating conditions.

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.