A generalized expression for accelerating beamlet decomposition simulations.

Paraxial diffraction modeling based on the Fourier transform has seen widespread implementation for simulating the response of a diffraction-limited optical system. For systems where the paraxial assumption is not sufficient, a class of algorithms has been developed that employs hybrid propagation physics to compute the propagation of an elementary beamlet along geometric ray paths. These "beamlet decomposition" algorithms include the well-known Gaussian beamlet decomposition (GBD) algorithm, of which several variants have been created. To increase the computational efficiency of the GBD algorithm, we derive an alternative expression of the technique that utilizes the analytical propagation of beamlets to tilted planes. We then use this accelerated algorithm to conduct a parameter-space search to find the optimal combination of free parameters in GBD to construct the analytical Airy function. The experiment is conducted on a consumer-grade CPU, and a high-performance GPU, where the new algorithm is 34 times faster than the previously published algorithm on CPUs, and 67,513 times faster on GPUs.

Exascale integrated modeling of low-order wavefront sensing and control for the Roman Coronagraph instrument.

Astronomical instruments to detect exoplanets require extreme wavefront stability. For these missions to succeed, comprehensive and precise modeling is required to design and analyze suitable coronagraphs and their wavefront control systems. In this paper, we describe techniques for integrated modeling at scale that is, to the best of our knowledge, 1000 times faster than previously published works. We show how this capability has been used to validate performance and perform uncertainty quantification for the Roman Coronagraph instrument. Finally, we show how this modeling capacity may be necessary to design and build the next generation of space-based coronagraph instruments.

PAP-LMPCR for improved, allele-specific footprinting and automated chromatin fine structure analysis.

The analysis of chromatin fine structure and transcription factor occupancy of differentially expressed genes by in vivo footprinting and ligation-mediated-PCR (LMPCR) is a powerful tool to understand the impact of chromatin on gene expression. However, as with all PCR-based techniques, the accuracy of the experiments has often been reduced by sequence similarities and the presence of GC-rich or repeat sequences, and some sequences are completely refractory to analysis. Here we describe a novel method, pyrophosphorolysis activated polymerization LMPCR or PAP-LMPCR, which is capable of generating accurate and reproducible footprints specific for individual alleles and can read through sequences previously not accessible for analysis. In addition, we have adapted this technique for automation, thus enabling the simultaneous and rapid analysis of chromatin structure at many different genes.

Multipotent progenitor cells isolated from adult human pancreatic tissue.

The supply of islet cells is a limiting factor for the widespread application of islet transplantation of type-1 diabetes. Islets constitute 1% to 2% of pancreatic tissue, leaving approximately 98% as discard after islet isolation and purification. In this report we present our data on the isolation of multipotent progenitor cells from discarded adult human pancreatic tissue. The collected cells from discarded nonislet fractions, after enzymatic digestion and gradient purification of islets, were dissociated for suspension culture in a serum-free medium. The cell clusters grown to a size of 100 to 150 mum contained cells staining for stage-specific embryonic antigens, but not insulin or C-peptide. To direct cell differentiation toward islets, clusters were recultured in a pancreatic differentiation medium. Insulin and C-peptide-positive cells by immunocytochemistry appeared within a week, reaching over 10% of the cell population. Glucagon and somatostatin-positive cells were also detected. The cell clusters were found to secrete insulin in response to glucose stimulation. Cells from the same clusters also had the capacity for differentiation into neural cells, as documented by staining for neural and glial cell markers when cultured as monolayers in media containing neurotrophic factors. These data suggest that multipotent pancreatic progenitor cells exist within the human pancreatic tissue that is typically discarded during islet isolation procedures. These adult progenitor cells can be successfully differentiated into insulin-producing cells, and thus they have the potential for treatment of type-1 diabetes mellitus.

Inhibition of p38 mitogen-activated protein kinase protects human islets from cryoinjury and improves the yield, viability, and quality of frozen-thawed islets.

The development of an optimal islet cryopreservation method will permit transplantation of islets from multiple donors in a single procedure and contribute to alleviation of the islet shortage. In this study, we have improved human islet cryopreservation methods under serum-free conditions using an intracellular-based islet cryopreservation solution (ICS), especially supplemented with a p38 pathway inhibitor (p38IH) to suppress p38 mitogen-activated protein kinase (MAPK) activation. Three different solutions were compared for freezing and thawing of human islets (1) conventional RPMI1640 medium, (2) ICS, and (3) ICS supplemented with a p38IH, SD-282 (ICS-p38IH). Islet cryopreservation with ICS-p38IH significantly improved islet recovery, viability, and quality after thawing of cryopreserved islets. This improvement may allow the use of cryopreserved islets in clinical islet transplantation.

Mice conditionally lacking the Wolfram gene in pancreatic islet beta cells exhibit diabetes as a result of enhanced endoplasmic reticulum stress and apoptosis.

AIMS/HYPOTHESIS: Wolfram syndrome is an autosomal recessive disorder characterised by childhood diabetes mellitus, optic atrophy and severe neurodegeneration, resulting in premature death. The aim of this study was to investigate the mechanisms responsible for the phenotype of carbohydrate intolerance and loss of pancreatic beta cells in this disorder. MATERIALS AND METHODS: To study the role of the Wolfram gene (Wfs1) in beta cells, we developed a mouse model with conditional deletion of Wfs1 in beta cells by crossing floxed Wfs1 exon 8 animals with mice expressing Cre recombinase under the control of a rat insulin promoter (RIP2-Cre). Complementary experiments using RNA interference of Wfs1 expression were performed in mouse insulinoma (MIN6) cell lines (WfsKD). RESULTS: Male knockout mice (betaWfs(-/-)) began developing variable and progressive glucose intolerance and concomitant insulin deficiency, compared with littermate controls, by 12 weeks of age. Analysis of islets from betaWfs(-/-) mice revealed a reduction in beta cell mass, enhanced apoptosis, elevation of a marker of endoplasmic reticulum stress (immunoglobulin heavy chain-binding protein [BiP]), and dilated endoplasmic reticulum with decreased secretory granules by electron microscopy. WfsKD cell lines had significantly increased apoptosis and elevated expression of the genes encoding BiP and C/EBP-homologous protein (CHOP), two markers of endoplasmic reticulum stress. CONCLUSIONS/INTERPRETATION: These results indicate that (1) lack of expression of Wfs1 in beta cells was sufficient to result in the diabetes mellitus phenotype; (2) beta cell death occurred by an accelerated process of apoptosis; and (3) lack of Wfs1 was associated with dilated endoplasmic reticulum and increased markers of endoplasmic reticulum stress, which appears to be a significant contributor to the reduction in beta cell survival.

Epigenetic changes and repositioning determine the evolutionary fate of duplicated genes.

Consideration of epigenetic silencing, perhaps by DNA methylation, led to an epigenetic complementation (EC) model for evolution by gene duplication (Rodin and Riggs (2003) J. Mol. Evol., 56, 718-729). This and subsequent work on genome-wide analyses of gene duplicates in several eukaryotic species pointed to a fundamental link between localization in the genome, epigenetic regulation of expression, and the evolutionary fate of new redundant gene copies, which can be either non- or neo-functionalization. Our main message in this report is that repositioning of a new duplicate to an ectopic site epigenetically alters its expression pattern, and concomitantly the rate and direction of mutations. Furthermore, comparison of syntenic vs. non-syntenic pairs of gene duplicates of different age unambiguously indicates that repositioning saves redundant gene duplicates from pseudogenization and hastens their evolution towards a new development-time and tissue-specific pattern of function.

Real-time PCR assay for quantitative mismatch detection.

We describe here a quantitative real-time PCR assay for the detection of single-base-pair differences that does not require fluorescently labeled gene-specific probes or complicated primer combinations. Following PCR or RT-PCR of a gene segment that may contain allele-specific differences, 100 pg amplified product are used for a real-time PCR with allele-specific primers and SYBR Green. The use of HEPES buffer at a pH of 6.95 together with AmpliTaq DNA polymerase results in a threshold difference between the correct template and the mismatched template of as many as 20 cycles, depending on the mismatch. Correct matches can be detected in an excess of mismatched template at least at the 0.01 level for the six primer-template matches versus mismatches tested: GC vs. A.C, AT vs. G.T, GC vs. C.C, GC vs. G.G, AT vs. C.T, and GC vs. G.A. Because the initial amplification is separate from real-time detection, conditions can be independently optimized for each step, making the assay particularly suitable for the detection of allele-specific expression in single cells.

Ligation-mediated PCR: robotic liquid handling for DNA damage and repair.

The investigation of in vivo DNA repair in mammalian cells at nucleotide resolution requires the quantification of break frequencies less than one per kilobase. By optimizing several parameters of the ligation-mediated PCR technique, we find that the required sensitivity can be achieved. We also report details of a one-day procedure that can be performed either with or without a robotic liquid handling workstation. The use of near-infrared fluorescent-labeled primers with detection by a LI-COR DNA sequencer provides for safe, nonradioactive detection, similar in sensitivity to the use of 32P-labeled primers but with the additional advantage that high-quality digitized data are obtained directly. Multiplexing can be performed; that is, more than one sequence can be analyzed in a single reaction, and multiple reactions can be processed robotically. Primer sets for exons 5-8 of the tumor suppressor gene, p53, were designed for simultaneous thermal cycling. The improved procedure with infrared detection was used to monitor low-frequency damage (<1 break/kb) and/or repair of UVB, UVC, and chemical methylation. Quantitative data on the linearity of response and reproducibility are described. The coefficient of variation for technical replicates was typically 10%. The methods described here will permit high sample throughput for the detection of DNA damage and repair as well as in vivo protein footprints.

X chromosome inactivation, differentiation, and DNA methylation revisited, with a tribute to Susumu Ohno.

X chromosome inactivation and DNA methylation are reviewed, with emphasis on the contributions of Susumu Ohno and the predictions made in my 1975 paper (Riggs, 1975), in which I proposed the "maintenance methylase" model for somatic inheritance of methylation patterns and suggested that DNA methylation would be involved in mammalian X chromosome inactivation and development. The maintenance methylase model is discussed and updated to consider methylation patterns in cell populations that have occasional, stochastic methylation changes by de novo methylation or demethylation, either active or passive. The "way station" model for the spread of X inactivation by LINE-1 elements is also considered, and some recent results from my laboratory are briefly reviewed.

Terminal transferase-dependent PCR (TDPCR) for in vivo UV photofootprinting of vertebrate cells.

Terminal transferase-dependent PCR (TDPCR) is a versatile, sensitive method for detecting DNA lesions such as those generated by the footprinting agents commonly used to detect in vivo protein-DNA interactions. Data similar to those obtained by ligation-mediated PCR (LMPCR) are obtained, but one advantage of TDPCR is that no special enzymes are needed other than terminal deoxynucleotide transferase, T4 DNA ligase, and thermostable DNA polymerases. A detailed TDPCR protocol is given for using UV photofootprinting to detect in vivo footprints and chromatin fine structure in vertebrate cells. One version of the protocol makes use of nonradioactive labeling by near-infrared fluorochromes and detection by a LI-COR DNA sequencing instrument. Sensitivity similar to that of (32)P-labeling is obtained, but with superior band resolution and quantitation.

Detection of antisense and ribozyme accessible sites on native mRNAs: application to NCOA3 mRNA.

The efficacies of antisense oligonucleotides and ribozymes are greatly dependent on the accessibility of their mRNA targets. Target site accessibility is affected by both RNA structure and the proteins associated along the length of the RNA. To mimic the native state of mRNA for site identification, we have previously used endogenous mRNAs in cellular extracts as targets for defined sequence oligodeoxynucleotides (ODNs) designed to identify both antisense pairing and potential ribozyme cleavage sites. The rationale for this approach is that the specific pairing of an ODN with a mRNA forms a DNA:RNA hybrid that is cleaved by the endogenous RNaseH in the cell extract. To extend the usefulness of this basic approach, we report here the use of semi-random ODN libraries to identify hammerhead ribozyme cleavage sites. Thus, the most accessible sites for antisense and ribozyme base pairing are selected by this approach. A novel feature of the approach described here is the use of terminal transferase-dependent PCR (TDPCR) after reverse transcription to estimate the cleavage efficiency and to precisely determine the RNaseH and ribozyme cleavage sites on mRNAs in cell extracts following treatment with ODN or ribozyme libraries. As a model system, we have targeted the NCOA3 (also known as AIB-1) mRNA in cell extracts. The NCOA3 mRNA encodes a nuclear receptor co-activator that is amplified and over-expressed in a high proportion of breast and ovarian cancers. A highly accessible site on this mRNA was identified, and a ribozyme targeted to this site was demonstrated to effectively downregulate NCOA3 function in cells.

Crystallization and preliminary X-ray analysis of neural haemoglobin from the nemertean worm Cerebratulus lacteus.

The nemertean worm Cerebratulus lacteus neural tissue haemoglobin (109 amino acids, the shortest known haemoglobin) has been overexpressed in Escherichia coli, purified and crystallized. A highly redundant native data set has been collected at the Cu K(alpha) wavelength to 2.05 A resolution. The crystals belong to the orthorhombic P2(1)2(1)2(1) space group, with unit-cell parameters a = 42.5, b = 43.1, c = 60.2 A and one molecule per asymmetric unit. The anomalous difference Patterson map clearly reveals the position of the haem Fe atom, thus paving the way for MAD/SAD structure determination.

A complex duplication created by gene targeting at the imprinted H19 locus results in two classes of methylation and correlated Igf2 expression phenotypes.

Imprinting of the mouse H19 and Igf2 genes is dependent on the presence of an intervening imprinting control region (ICR) situated 2 kb upstream of H19 and approximately 70 kb downstream of Igf2. Several recent studies have provided substantial evidence that the unmethylated maternal ICR acts as an insulator that prevents activation of Igf2 by a suite of enhancers downstream of the H19 gene. The methylated paternal ICR and H19 promoter have no activity, allowing sole activation of Igf2 expression. We have produced mice in which a duplication of the H19/Igf2 ICR produces, in each generation, two classes of methylation levels that correlated with two Igf2 imprinting phenotypes. One hypermethylated class also shows activation of the normally silent Igf2 gene, whereas the other hypomethylated class shows only slight activation of Igf2, in agreement with methylation's role in ICR function. This study describes a rare, possibly unique type of mutation that induces two distinct phenotypes in each generation.

Mapping psoralen cross-links at the nucleotide level in mammalian cells: suppression of cross-linking at transcription factor- or nucleosome-binding sites.

We have developed a new genomic sequencing method for detecting, with resolution at the nucleotide level, the interstrand DNA cross-links induced by 4,5',8-trimethylpsoralen along single-copy genes in mammalian cells. The cross-links (diadducts) initially formed are converted into monoadducts by alkali reversal prior to the use of terminal transferase-dependent PCR (TD-PCR). After alkali reversal, but not before, the DNA strands can be separated and used as templates for gene-specific primer extension, which is the first step in the TD-PCR procedure. The converted psoralen adducts block primer extension, and the prematurely terminated single-stranded products are then amplified by TD-PCR and visualized on a sequencing gel. Adducts formed by angelicin, a psoralen derivative that forms only monoadducts, were also investigated by use of TD-PCR. Comparison of the adduct distribution patterns of in vivo-treated DNA with those of in vitro-treated DNA revealed that the binding of transcription factors inhibited both psoralen cross-linking and angelicin monoadduct formation in the c-JUN and c-FOS promoters in living human cells. Adduct formation was also inhibited in the region of a putative positioned nucleosome in the c-FOS promoter. These methods should be of general use for study of in vivo protein-DNA interactions and DNA repair.

Ligation-mediated PCR for quantitative in vivo footprinting.

Ligation-mediated polymerase chain reaction (LM-PCR) is a genomic analysis technique for determination of (1) primary DNA nucleotide sequences (2) cytosine methylation patterns (3) DNA lesion formation and repair, and (4) in vivo protein-DNA footprints. However, LM-PCR can be limited by the multiple steps required and the relatively short stretch of sequence (usually <200 bp) that can be analyzed per reaction. We report here a simplified, one-day LM-PCR protocol in which all pipetting steps can be performed by a robotic workstation and which, moreover, provides longer reads (>350 bp) and enhanced signal quality by use of nonradioactive detection and a LI-COR DNA sequencing instrument. Sensitivity comparable to radiolabeling is achieved using oligonucleotide primers that are 5'-end labeled with infrared fluorochromes. We showed that the technique could be used for sensitive and reproducible in vivo photofootprinting of the human phosphoglycerate kinase 1 (PGK1) promoter, as well as providing good Maxam-Gilbert sequence information. The methods described here should allow high-throughput, high-resolution analysis of transcription factor binding and chromatin structure, and also may be useful for sequencing gaps that are refractory to cloning.

Examination of the DNA substrate selectivity of DNA cytosine methyltransferases using mass tagging.

The biological significance of cytosine methylation is as yet incompletely understood, but substantial and growing evidence strongly suggests that perturbation of methylation patterns, resulting from the infidelity of DNA cytosine methyltransferase, is an important component of the development of human cancer. We have developed a novel in vitro assay that allows us to quantitatively determine the DNA substrate preferences of cytosine methylases. This approach, which we call mass tagging, involves the labeling of target cytosine residues in synthetic DNA duplexes with stable isotopes, such as (15)N. Methylation is then measured by the formation of 5-methylcytosine (5mC) by gas chromatography/mass spectrometry. The DNA substrate selectivity is determined from the mass spectrum of the product 5mC. With the non-symmetrical duplex DNA substrate examined in this study we find that the bacterial methyltransferase HPA:II (duplex DNA recognition sequence CCGG) methylates the one methylatable cytosine of each strand similarly. Introduction of an A-C mispair at the methylation site shifts methylation exclusively to the mispaired cytosine residue. In direct competition assays with HPA:II methylase we observe that the mispaired substrate is methylated more extensively than the fully complementary, normal substrate, although both have one HPA:II methylation site. Through the use of this approach we will be able to learn more about the mechanisms by which methylation patterns can become altered.

Chromatin fine structure profiles for a developmentally regulated gene: reorganization of the lysozyme locus before trans-activator binding and gene expression.

The chicken lysozyme locus is activated in a stepwise fashion during myeloid differentiation. We have used this locus as a model to study at high resolution changes in chromatin structure both in chicken cell lines representing various stages of macrophage differentiation and in primary cells from transgenic mice. In this study we have addressed the question of whether chromatin rearrangements can be detected in myeloid precursor cells at a stage well before overt transcription of the lysozyme gene begins. In addition to restriction enzyme accessibility assays and DMS footprinting, we have applied new, very sensitive techniques to assay for chromatin changes. Particularly informative was UV photofootprinting, using terminal transferase-dependent PCR and nonradioactive detection. We find that the basic chromatin structure in lysozyme nonexpressing hematopoietic precursor cells is highly similar to the pattern found in fully differentiated lysozyme-expressing cells. In addition, we find that only in nonexpressing cells are dimethylsulfate footprints and UV photofootprints affected by trichostatin, an inhibitor of histone deacetylation. These results are interpreted to mean that most chromatin pattern formation is complete before the binding of end-stage trans-activators, supporting the notion that heritable chromatin structure is central to the stable epigenetic programs that guide development.

Nitric-oxide dioxygenase activity and function of flavohemoglobins. sensitivity to nitric oxide and carbon monoxide inhibition.

Widely distributed flavohemoglobins (flavoHbs) function as NO dioxygenases and confer upon cells a resistance to NO toxicity. FlavoHbs from Saccharomyces cerevisiae, Alcaligenes eutrophus, and Escherichia coli share similar spectra, O(2), NO, and CO binding kinetics, and steady-state NO dioxygenation kinetics. Turnover numbers (V(max)) for S. cerevisiae, A. eutrophus, and E. coli flavoHbs are 112, 290, and 365 NO heme(-1) s(-1), respectively, at 37 degrees C with 200 microm O(2). The K(M) values for NO are low and range from 0.1 to 0.25 microm. V(max)/K(M)(NO) ratios of 900-2900 microm(-1) s(-1) indicate an extremely efficient dioxygenation mechanism. Approximate K(M) values for O(2) range from 60 to 90 microm. NO inhibits the dioxygenases at NO:O(2) ratios of > or =1:100 and makes true K(M)(O(2)) values difficult to determine. High and roughly equal second order rate constants for O(2) and NO association with the reduced flavoHbs (17-50 microm(-1) s(-1)) and small NO dissociation rate constants suggest that NO inhibits the dioxygenase reaction by forming inactive flavoHbNO complexes. Carbon monoxide also binds reduced flavoHbs with high affinity and competitively inhibits NO dioxygenases with respect to O(2) (K(I)(CO) = approximately 1 microm). These results suggest that flavoHbs and related hemoglobins evolved as NO detoxifying components of nitrogen metabolism capable of discriminating O(2) from inhibitory NO and CO.