Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays

Genomic aberrations recurrent in a particular cancer type can be important prognostic markers for tumor progression. Typically in early tumorigenesis, cells incur a breakdown of the DNA replication machinery that results in an accumulation of genomic aberrations in the form of duplications, deletions, translocations, and other genomic alterations. Microarray methods allow for finer mapping of these aberrations than has previously been possible; however, data processing and analysis methods have not taken full advantage of this higher resolution. Attention has primarily been given to analysis on the single sample level, where multiple adjacent probes are necessarily used as replicates for the local region containing their target sequences. However, regions of concordant aberration can be short enough to be detected by only one, or very few, array elements. We describe a method called Multiple Sample Analysis for assessing the significance of concordant genomic aberrations across multiple experiments that does not require a-priori definition of aberration calls for each sample. If there are multiple samples, representing a class, then by exploiting the replication across samples our method can detect concordant aberrations at much higher resolution than can be derived from current single sample approaches. Additionally, this method provides a meaningful approach to addressing population-based questions such as determining important regions for a cancer subtype of interest or determining regions of copy number variation in a population. Multiple Sample Analysis also provides single sample aberration calls in the locations of significant concordance, producing high resolution calls per sample, in concordant regions. The approach is demonstrated on a dataset representing a challenging but important resource: breast tumors that have been formalin-fixed, paraffin-embedded, archived, and subsequently UV-laser capture microdissected and hybridized to two-channel BAC arrays using an amplification protocol. We demonstrate the accurate detection on simulated data, and on real datasets involving known regions of aberration within subtypes of breast cancer at a resolution consistent with that of the array. Similarly, we apply our method to previously published datasets, including a 250K SNP array, and verify known results as well as detect novel regions of concordant aberration. The algorithm has been fully implemented and tested and is freely available as a Java application at http://www.cbil.upenn.edu/MSA.     

Temporal organization of multi-whisker contact in rats

Previous work has established that during exploration and discrimination, rats move their whiskers at frequencies between 6 and 12 Hz and that whisking frequency changes during contact. One critical component of any tactile system is contact. In the rat whisker system, such contacts may involve one or more vibrissa in the whisker array and contact duration of each whisker may vary over a considerable range, depending upon the behavioral context. However, little is known about the variables controlling contact duration or about the temporal relationships among contacts by adjacent whiskers. To address these issues head fixed rats were trained to touch a piezo-contact-sensor with the shaft of their whiskers (Bermejo and Zeigler, Somatosens Mot Res 17: 373-377, 2000). During the task, whisker movements and contacts were monitored with a high-speed camera at 500 frames/s and stored on videotape. To facilitate analysis, animals had their whiskers selectively trimmed. Data are reported from animals with C1 & C2, D1 & D2, or Arc2 (E2, D2, C2, B2) whiskers intact. For both row and arc animals, when just a single whisker touched the sensor the duration of contact was significantly shorter than when multiple whiskers made contact. When multiple whiskers made contact, onset was rarely simultaneous. Furthermore, in row-intact animals, contact progressed in an orderly fashion such that the rostral whisker in a row made contact first followed 24 ms (SE = 1.9 ms) later by the caudal whisker. When contact reversed the caudal whisker lifted off first, followed by the rostral whisker. Thus, the order in which whiskers touch an object regulates contact duration: the first whisker to touch the sensor stays in contact longer than any other whisker. The temporal discharge properties of neurons in the trigeminal system are expected to reflect position of whiskers on the nose.     

Binding of simple carbohydrates and some N-acetyllactosamine-containing oligosaccharides to Erythrina cristagalli agglutinin as followed with a fluorescent indicator ligand

Erythrina cristagalli agglutinin, a dimeric lectin [J. L. Iglesias, et al. (1982) Eur. J. Biochem.123, 247–252] was shown by equilibrium dialysis to be bivalent for 4-methylumbelliferyl-β-d-galactoside. Upon binding to the lectin, this ligand showed a difference absorption spectrum with two maxima (at 322 and 336 nm) of equal intensity (Δ? = 1.2 × 10³m^?1 cm^?1). A similar spectrum with a comparable value of Δ? was obtained with 4-methylumbelliferyl-N-acetyl-β-d-galactosaminide. Binding of methyl-α-d-galactoside, lactose, and N-acetyllactosamine all produced small but equally intense protein difference spectra with a maximum (Δ? = 2.8 × 10² M^?1 cm^?1) at 291.6 nm. Upon binding of N-dansyl-d-galactosamine to the lectin, there was a fivefold increase in fluorescence intensity of this ligand. The association constant for N-dansyl-d-galactosamine was caused by a very favorable ΔS° of the dansyl group without affecting the strictly carbohydrate-specific character of binding. N-Dansyl-d-galactosamine was employed as a fluorescent indicator ligand in substitution titrations. This involved the use of simple carbohydrates, N-acetyllactosamine, and oligosaccharides which occur in the carbohydrate units of N-glycoproteins; the latter were Gal(β → 4)GlcNAc(β1 → 2)Man, Gal(β1 → 4)GlcNAc(β1 → 6)Man, and Gal(β1 → 4)GlcNAc(β1 → 6)[Gal(β1 → 4)GlcNAc(β1 → 2)]Man. The titrations were performed at two temperatures to determine the thermodynamic parameters. In the series N-acetyl-d-galactosamine, methyl-α-d-galactoside, and lactose, ?ΔH° increased from 24 to 41 kJ mol^?1; it increased further for N-acetyllactosamine and then remained unchanged for the N-acetyllactosamine-containing oligosaccharides (55 ± 1 kJ mol^?1). This indicated that the site specifically accommodated the disaccharide structure with an important contribution of the 2-acetamido group in the penultimate sugar. Beyond this, no additional contacts seemed to be formed. This conclusion also followed from considerations of ΔS° values which became more unfavorable in the above series (?23 to ?101 ± 4 J mol^?1 K^?1); the most negative value of ΔS° was observed with N-acetyllactosamine and the three N-acetyllactosamine-containing oligosaccharides.     

The N-terminal domain of the regulatory subunit is sufficient for complete activation of acetohydroxyacid synthase III from Escherichia coli

We have previously proposed a model for the fold of the N-terminal domain of the small, regulatory subunit (SSU) of acetohydroxyacid synthase isozyme III. The fold is an alpha-beta sandwich with betaalphabetabetaalphabeta topology, structurally homologous to the C-terminal regulatory domain of 3-phosphoglycerate dehydrogenase. We suggested that the N-terminal domains of a pair of SSUs interact in the holoenzyme to form two binding sites for the feedback inhibitor valine in the interface between them. The model was supported by mutational analysis and other evidence. We have now examined the role of the C-terminal portion of the SSU by construction of truncated polypeptides (lacking 35, 48, 80, 95, or 112 amino acid residues from the C terminus) and examining the properties of holoenzymes reconstituted using these constructs. The Delta35, Delta48, and Delta80 constructs all lead to essentially complete activation of the catalytic subunits. The Delta80 construct, corresponding to the putative N-terminal domain, has the highest level of affinity for the catalytic subunits and leads to a reconstituted enzyme with k(cat)/K(M) about twice that of the wild-type enzyme. On the other hand, none of these constructs binds valine or leads to a valine-sensitive enzyme on reconstitution. The enzyme reconstituted with the Delta80 construct does not bind valine, either. The N-terminal portion (about 80 amino acid residues) of the SSU is thus necessary and sufficient for recognition and activation of the catalytic subunits, but the C-terminal half of the SSU is required for valine binding and response. We suggest that the C-terminal region of the SSU contributes to monomer-monomer interactions, and provide additional experimental evidence for this suggestion.     

Listeria monocytogenes sigma B regulates stress response and virulence functions

While the stress-responsive alternative sigma factor sigma(B) has been identified in different species of Bacillus, Listeria, and Staphylococcus, the sigma(B) regulon has been extensively characterized only in B. subtilis. We combined biocomputing and microarray-based strategies to identify sigma(B)-dependent genes in the facultative intracellular pathogen Listeria monocytogenes. Hidden Markov model (HMM)-based searches identified 170 candidate sigma(B)-dependent promoter sequences in the strain EGD-e genome sequence. These data were used to develop a specialized, 208-gene microarray, which included 166 genes downstream of HMM-predicted sigma(B)-dependent promoters as well as selected virulence and stress response genes. RNA for the microarray experiments was isolated from both wild-type and Delta sigB null mutant L. monocytogenes cells grown to stationary phase or exposed to osmotic stress (0.5 M KCl). Microarray analyses identified a total of 55 genes with statistically significant sigma(B)-dependent expression under the conditions used in these experiments, with at least 1.5-fold-higher expression in the wild type over the sigB mutant under either stress condition (51 genes showed at least 2.0-fold-higher expression in the wild type). Of the 55 genes exhibiting sigma(B)-dependent expression, 54 were preceded by a sequence resembling the sigma(B) promoter consensus sequence. Rapid amplification of cDNA ends-PCR was used to confirm the sigma(B)-dependent nature of a subset of eight selected promoter regions. Notably, the sigma(B)-dependent L. monocytogenes genes identified through this HMM/microarray strategy included both stress response genes (e.g., gadB, ctc, and the glutathione reductase gene lmo1433) and virulence genes (e.g., inlA, inlB, and bsh). Our data demonstrate that, in addition to regulating expression of genes important for survival under environmental stress conditions, sigma(B) also contributes to regulation of virulence gene expression in L. monocytogenes. These findings strongly suggest that sigma(B) contributes to L. monocytogenes gene expression during infection.     

Patterns of female dominance in Propithecus diadema edwardsi of Ranomafana National Park, Madagascar

Many lemur species are characterized by some form of female dominance, ranging from female feeding priority to complete female dominance, although this is a rare trait in primates and other mammals. The status of the Milne-Edwards' sifaka (Propithecus diadema edwardsi), a diurnal lemur, is ambiguous. Some short-term studies have found little or no aggression. The aim of the current, long-term study was to quantify the intersexual-dominance patterns of this sifaka. The distribution, outcome, and context of aggressive interactions were studied in four groups of wild sifakas. The majority of intersexual aggressive interactions were decided, with the loser expressing submissive behavior. Intersexual aggressive interactions occurred in all social contexts, and within all social contexts the females won the vast majority (92.7-96.0%) of aggressive interactions. While aggression rates were low (0.22/hr), this evidence suggests female dominance. We propose that female dominance exists because it provides a fitness advantage to both males and females.     

Structural studies of [2',6'-dimethyl-L-tyrosine1]endomorphin-2 analogues: enhanced activity and cis orientation of the Dmt-Pro amide bond

Analogues of endomorphin-2 (EM-2: Tyr-Pro-Phe-Phe-NH(2)) (1) were designed to examine the importance of each residue on mu-opioid receptor interaction. Replacement of Tyr(1) by 2',6'-dimethyl-L-tyrosine (Dmt) (9-12) exerted profound effects: [Dmt(1)]EM-2 (9) elevated mu-opioid affinity 4.6-fold (K(i mu=0.15 nM) yet selectivity fell 330-fold as delta-affinity rose (K(i)delta=28.2 nM). This simultaneous increased mu- and delta-receptor bioactivities resulted in dual agonism (IC(50)=0.07 and 1.87 nM, respectively). While substitution of Phe(4) by a phenethyl group (4) decreased mu affinity (K(i)mu=13.3 nM), the same derivative containing Dmt (12) was comparable to EM-2 but also acquired weak delta antagonism (pA(2)=7.05). 1H NMR spectroscopy revealed a trans configuration (1:2 to 1:3, cis/trans) in the Tyr-Pro amide bond, but a cis configuration (5:3 to 13:7, cis/trans) with Dmt-Pro analogues.     

The NFP locus of Medicago truncatula controls an early step of Nod factor signal transduction upstream of a rapid calcium flux and root hair deformation

Establishment of the Rhizobium-legume symbiosis depends on a molecular dialogue, in which rhizobial nodulation (Nod) factors act as symbiotic signals, playing a key role in the control of specificity of infection and nodule formation. Using nodulation-defective (Nod-) mutants of Medicago truncatula to study the mechanisms controlling Nod factor perception and signalling, we have previously identified five genes that control components of a Nod factor-activated signal transduction pathway. Characterisation of a new M. truncatula Nod- mutant led to the identification of the Nod Factor Perception (NFP) locus. The nfp mutant has a novel phenotype among Nod- mutants of M. truncatula, as it does not respond to Nod factors by any of the responses tested. The nfp mutant thus shows no rapid calcium flux, the earliest detectable Nod factor response of wild-type plants, and no root hair deformation. The nfp mutant is also deficient in Nod factor-induced calcium spiking and early nodulin gene expression. While certain genes controlling Nod factor signal transduction also control the establishment of an arbuscular mycorrhizal symbiosis, the nfp mutant shows a wild-type mycorrhizal phenotype. These data indicate that the NFP locus controls an early step of Nod factor signal transduction, upstream of previously identified genes and specific to nodulation.