The Montage Image Mosaic Service: Custom Image Mosaics On-Demand

The Montage software suite has proven extremely useful as a general engine for reprojecting, background matching, and mosaicking astronomical image data from a wide variety of sources. The processing algorithms support all common World Coordinate System (WCS) projections and have been shown to be both astrometrically accurate and flux conserving. The background ‘matching’ algorithm does not remove background flux but rather finds the best compromise background based on all the input and matches the individual images to that. The Infrared Science Archive (IRSA), part of the Infrared Processing and Analysis Center (IPAC) at Caltech, has now wrapped the Montage software as a CGI service and provided a compute and request management infrastructure capable of producing approximately 2 TBytes / day of image mosaic output (e.g. from 2MASS and SDSS data). Besides the basic Montage engine, this service makes use of a 16-node LINUX cluster (dual processor, dual core) and the ROME request management software developed by the National Virtual Observatory (NVO). ROME uses EJB/database technology to manage user requests, queue processing and load balance between users, and managing job monitoring and user notification. The Montage service will be extended to process userdefined data collections, including private data uploads

A Cross-Match of 2MASS and SDSS: Newly-Found L and T Dwarfs and an Estimate of the Space Densitfy of T Dwarfs

We report new L and T dwarfs found in a cross-match of the SDSS Data Release 1 and 2MASS. Our simultaneous search of the two databases effectively allows us to relax the criteria for object detection in either survey and to explore the combined databases to a greater completeness level. We find two new T dwarfs in addition to the 13 already known in the SDSS DR1 footprint. We also identify 22 new candidate and bona-fide L dwarfs, including a new young L2 dwarf and a peculiar L2 dwarf with unusually blue near-IR colors: potentially the result of mildly sub-solar metallicity. These discoveries underscore the utility of simultaneous database cross-correlation in searching for rare objects. Our cross-match completes the census of T dwarfs within the joint SDSS and 2MASS flux limits to the 97% level. Hence, we are able to accurately infer the space density of T dwarfs. We employ Monte Carlo tools to simulate the observed population of SDSS DR1 T dwarfs with 2MASS counterparts and find that the space density of T0-T8 dwarf systems is 0.0070 (-0.0030; +0.0032) per cubic parsec (95% confidence interval), i.e., about one per 140 cubic parsecs. Compared to predictions for the T dwarf space density that depend on various assumptions for the sub-stellar mass function, this result is most consistent with models that assume a flat sub-stellar mass function dN/dM ~ M^0. No >T8 dwarfs were discovered in the present cross-match, though less than one was expected in the limited area (2099 sq. degrees) of SDSS DR1.Comment: To appear in ApJ, Feb 10, 2008 issue. 37 pages, including 12 figures and 14 table

The First Data Release of the KODIAQ Survey

We present and make publicly available the first data release (DR1) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. The KODIAQ survey is aimed at studying galactic and circumgalactic gas in absorption at high-redshift, with a focus on highly-ionized gas traced by OVI, using the HIRES spectrograph on the Keck-I telescope. KODIAQ DR1 consists of a fully-reduced sample of 170 quasars at 0.29 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000) between 2004 and 2012. DR1 contains 247 spectra available in continuum normalized form, representing a sum total exposure time of ~1.6 megaseconds. These co-added spectra arise from a total of 567 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a HIRES data reduction package made available through the XIDL distribution. DR1 is publicly available to the community, housed as a higher level science product at the KOA. We will provide future data releases that make further QSOs, including those with pre-2004 observations taken with the previous-generation HIRES detectors.Comment: 14 pages, 9 figures, Submitted to AJ. All data products available at the Keck Observatory Archive beginning May 15, 2015. URL: https://koa.ipac.caltech.edu/applications/KODIA

The impact of mutation and gene conversion on the local diversification of antigen genes in African trypanosomes

Patterns of genetic diversity in parasite antigen gene families hold important information about their potential to generate antigenic variation within and between hosts. The evolution of such gene families is typically driven by gene duplication, followed by point mutation and gene conversion. There is great interest in estimating the rates of these processes from molecular sequences for understanding the evolution of the pathogen and its significance for infection processes. In this study, a series of models are constructed to investigate hypotheses about the nucleotide diversity patterns between closely related gene sequences from the antigen gene archive of the African trypanosome, the protozoan parasite causative of human sleeping sickness in Equatorial Africa. We use a hidden Markov model approach to identify two scales of diversification: clustering of sequence mismatches, a putative indicator of gene conversion events with other lower-identity donor genes in the archive, and at a sparser scale, isolated mismatches, likely arising from independent point mutations. In addition to quantifying the respective probabilities of occurrence of these two processes, our approach yields estimates for the gene conversion tract length distribution and the average diversity contributed locally by conversion events. Model fitting is conducted using a Bayesian framework. We find that diversifying gene conversion events with lower-identity partners occur at least five times less frequently than point mutations on variant surface glycoprotein (VSG) pairs, and the average imported conversion tract is between 14 and 25 nucleotides long. However, because of the high diversity introduced by gene conversion, the two processes have almost equal impact on the per-nucleotide rate of sequence diversification between VSG subfamily members. We are able to disentangle the most likely locations of point mutations and conversions on each aligned gene pair

The Design and Operation of The Keck Observatory Archive

The Infrared Processing and Analysis Center (IPAC) and the W. M. Keck Observatory (WMKO) operate an archive for the Keck Observatory. At the end of 2013, KOA completed the ingestion of data from all eight active observatory instruments. KOA will continue to ingest all newly obtained observations, at an anticipated volume of 4 TB per year. The data are transmitted electronically from WMKO to IPAC for storage and curation. Access to data is governed by a data use policy, and approximately two-thirds of the data in the archive are public.Comment: 12 pages, 4 figs, 4 tables. Presented at Software and Cyberinfrastructure for Astronomy III, SPIE Astronomical Telescopes + Instrumentation 2014. June 2014, Montreal, Canad

Shifts in Attack Behavior of an Important Kelp Forest Predator Within Marine Reserves

Marine reserves have become increasingly valuable tools with which to manage ecosystems. These reserves consistently restore populations of top predators, often reducing availability of their favored prey. We hypothesized that such prey reduction in reserves causes protected predators to alter their attack behavior to include less palatable prey, potentially amplifying top-down effects on community structure. To test this hypothesis, we presented the relatively unpalatable sea hare Aplysia californica to freely foraging spiny lobsters Panulirus interruptus in 4 marine no-take reserves, each paired with an adjacent fished area. We found that lobsters only attacked sea hares inside reserves, where lobster density was significantly greater than that of the adjacent fished areas. Attacks on otherwise unpalatable prey exclusively in no-take reserves was likely caused by increased hunger, since in the laboratory only food-deprived lobsters attacked sea hares. These findings are the first to suggest that management involving no-take reserves may have unintended consequences on community structure that result from behavioral changes in key predators in the face of increased competition for food. We suspect that these effects may become more widely detected as reserves across the globe grow older and are researched further

Montage: a grid portal and software toolkit for science-grade astronomical image mosaicking

Montage is a portable software toolkit for constructing custom, science-grade mosaics by composing multiple astronomical images. The mosaics constructed by Montage preserve the astrometry (position) and photometry (intensity) of the sources in the input images. The mosaic to be constructed is specified by the user in terms of a set of parameters, including dataset and wavelength to be used, location and size on the sky, coordinate system and projection, and spatial sampling rate. Many astronomical datasets are massive, and are stored in distributed archives that are, in most cases, remote with respect to the available computational resources. Montage can be run on both single- and multi-processor computers, including clusters and grids. Standard grid tools are used to run Montage in the case where the data or computers used to construct a mosaic are located remotely on the Internet. This paper describes the architecture, algorithms, and usage of Montage as both a software toolkit and as a grid portal. Timing results are provided to show how Montage performance scales with number of processors on a cluster computer. In addition, we compare the performance of two methods of running Montage in parallel on a grid.Comment: 16 pages, 11 figure

Gene expression in Leishmania is regulated predominantly by gene dosage

ABSTRACT Leishmania tropica, a unicellular eukaryotic parasite present in North and East Africa, the Middle East, and the Indian subcontinent, has been linked to large outbreaks of cutaneous leishmaniasis in displaced populations in Iraq, Jordan, and Syria. Here, we report the genome sequence of this pathogen and 7,863 identified protein-coding genes, and we show that the majority of clinical isolates possess high levels of allelic diversity, genetic admixture, heterozygosity, and extensive aneuploidy. By utilizing paired genome-wide high-throughput DNA sequencing (DNA-seq) with RNA-seq, we found that gene dosage, at the level of individual genes or chromosomal “somy” (a general term covering disomy, trisomy, tetrasomy, etc.), accounted for greater than 85% of total gene expression variation in genes with a 2-fold or greater change in expression. High gene copy number variation (CNV) among membrane-bound transporters, a class of proteins previously implicated in drug resistance, was found for the most highly differentially expressed genes. Our results suggest that gene dosage is an adaptive trait that confers phenotypic plasticity among natural Leishmania populations by rapid down- or upregulation of transporter proteins to limit the effects of environmental stresses, such as drug selection. IMPORTANCE Leishmania is a genus of unicellular eukaryotic parasites that is responsible for a spectrum of human diseases that range from cutaneous leishmaniasis (CL) and mucocutaneous leishmaniasis (MCL) to life-threatening visceral leishmaniasis (VL). Developmental and strain-specific gene expression is largely thought to be due to mRNA message stability or posttranscriptional regulatory networks for this species, whose genome is organized into polycistronic gene clusters in the absence of promoter-mediated regulation of transcription initiation of nuclear genes. Genetic hybridization has been demonstrated to yield dramatic structural genomic variation, but whether such changes in gene dosage impact gene expression has not been formally investigated. Here we show that the predominant mechanism determining transcript abundance differences (>85%) in Leishmania tropica is that of gene dosage at the level of individual genes or chromosomal somy

Montage: An Astronomical Image Mosaic Service for the NVO

Montage is a software system for generating astronomical image mosaics according to user-specified size, rotation, WCS-compliant projection and coordinate system, with background modeling and rectification capabilities. Its architecture has been described in the proceedings of ADASS XII and XIII (Berriman et al. 2003, 2004). It has been designed as a toolkit, with independent modules for image reprojection, background rectification and co-addition, and will run on workstations, clusters and grids. The primary limitation of Montage thus far has been in the projection algorithm. It uses a spherical trigonometry approach that is general at the expense of speed. The reprojection algorithm has now been made 30 times faster for commonly used tangent plane to tangent plane reprojections that cover up to several square degrees, through modification of a custom algorithm first derived for the Spitzer Space Telescope. This focus session will describe this algorithm, demonstrate the generation of mosaics in real time, and describe applications of the software. In particular, we will highlight one case study which shows how Montage is supporting the generation of science-grade mosaics of images measured with the Infrared Array Camera aboard the Spitzer Space Telescope

Data reduction pipelines for the Keck Observatory Archive

The Keck Observatory Archive (KOA) currently serves ~ 42 TB of data spanning over 20 years from all ten past and current facility instruments at Keck. Although most of the available data are in the raw form, for four instruments (HIRES, NIRC2, OSIRIS, LWS), quick-look, browse products generated by automated pipelines are also offered to facilitate assessment of the scientific content and quality of the data. KOA underwrote the update of the MAKEE package to support reduction of the CCD upgrade to HIRES, developed scripts for reduction of NIRC2 data and automated the existing OSIRIS and LWS data reduction packages. We describe in some detail the recently completed automated pipeline for NIRSPEC, which will be used to create browse products in KOA and made available for quicklook of the data by the observers at the telescope. We review the currently available data reduction tools for Keck data, and present our plans and anticipated priorities for the development of automated pipelines and release of reduced data products for the rest of the current and future instruments. We also anticipate that Keck's newest instrument, NIRES, which will be delivered with a fully automated pipeline, will be the first to have both raw and level-1 data ingested at commissioning