Microsatellite Content in 397 Nuclear Exons and Their Flanking Regions in the Fern Family Ophioglossaceae

Microsatellites or SSRs are small tandem repeats that are 1-6 bp long. They are usually highly polymorphic and form important portions of genomes. They have been extensively analyzed in humans, animals and model plants; however, information from non-flowering plants is generally lacking. Here, we examined 29 samples of Ophioglossaceae ferns, mainly from the genera Botrychium and Sceptridium. We analyzed the SSR distribution, density and composition in almost 400 nuclear exons and their flanking regions. We detected 45 SSRs in exons and 1475 SSRs in the flanking regions. In the exons, only di-, tri- and tetranucleotides were found, and all of them were 12 bp long. The annotation of the exons containing SSRs showed that they were related to various processes, such as metabolism, catalysis, transportation or plant growth. The flanking regions contained SSRs from all categories, with the most numerous being dinucleotides, followed by tetranucleotides. More than one-third of all the SSRs in the flanking regions were 12 bp long. The SSR densities in the exons were very low, ranging from 0 to 0.07 SSRs/kb, while those in the flanking regions ranged from 0.24 to 0.81 SSRs/kb; and those in the combined dataset ranged from 0.2 to 0.81 SSRs/kb. The majority of the detected SSRs in the flanking regions were polymorphic and present at the same loci across two or more samples but differing in the number of repeats. The SSRs detected here may serve as a basis for further population genetic, phylogenetic or evolutionary genetic studies, as well as for further studies focusing on SSRs in the genomes and their roles in adaptation, evolution and diseases.

Organellar phylogenomics of Ophioglossaceae fern genera

Previous phylogenies showed conflicting relationships among the subfamilies and genera within the fern family Ophioglossaceae. However, their classification remains unsettled where contrasting classifications recognize four to 15 genera. Since these treatments are mostly based on phylogenetic evidence using limited, plastid-only loci, a phylogenomic understanding is actually necessary to provide conclusive insight into the systematics of the genera. In this study, we have therefore compiled datasets with the broadest sampling of Ophioglossaceae genera to date, including all fifteen currently recognized genera, especially for the first time the South African endemic genus Rhizoglossum. Notably, our comprehensive phylogenomic matrix is based on both plastome and mitogenome genes. Inferred from the coding sequences of 83 plastid and 37 mitochondrial genes, a strongly supported topology for these subfamilies is presented, and is established by analyses using different partitioning approaches and substitution models. At the generic level, most relationships are well resolved except for few within the subfamily Ophioglossoideae. With this new phylogenomic scheme, key morphological and genomic changes were further identified along this backbone. In addition, we confirmed numerous horizontally transferred (HGT) genes in the genera Botrypus, Helminthostachys, Mankyua, Sahashia, and Sceptridium. These HGT genes are most likely located in mitogenomes and are predominately donated from angiosperm Santalales or non-Ophioglossaceae ferns. By our in-depth searches of the organellar genomes, we also provided phylogenetic overviews for the plastid and mitochondrial MORFFO genes found in these Ophioglossaceae ferns.

PHASTEST: faster than PHASTER, better than PHAST

PHASTEST (PHAge Search Tool with Enhanced Sequence Translation) is the successor to the PHAST and PHASTER prophage finding web servers. PHASTEST is designed to support the rapid identification, annotation and visualization of prophage sequences within bacterial genomes and plasmids. PHASTEST also supports rapid annotation and interactive visualization of all other genes (protein coding regions, tRNA/tmRNA/rRNA sequences) in bacterial genomes. Given that bacterial genome sequencing has become so routine, the need for fast tools to comprehensively annotate bacterial genomes has become progressively more important. PHASTEST not only offers faster and more accurate prophage annotations than its predecessors, it also provides more complete whole genome annotations and much improved genome visualization capabilities. In standardized tests, we found that PHASTEST is 31% faster and 2-3% more accurate in prophage identification than PHASTER. Specifically, PHASTEST can process a typical bacterial genome in 3.2 min (raw sequence) or in 1.3 min when given a pre-annotated GenBank file. Improvements in PHASTEST's ability to annotate bacterial genomes now make it a particularly powerful tool for whole genome annotation. In addition, PHASTEST now offers a much more modern and responsive visualization interface that allows users to generate, edit, annotate and interactively visualize (via zooming, rotating, dragging, panning, resetting), colourful, publication quality genome maps. PHASTEST continues to offer popular options such as an API for programmatic queries, a Docker image for local installations, support for multiple (metagenomic) queries and the ability to perform automated look-ups against thousands of previously PHAST-annotated bacterial genomes. PHASTEST is available online at https://phastest.ca.

Proksee: in-depth characterization and visualization of bacterial genomes

Proksee (https://proksee.ca) provides users with a powerful, easy-to-use, and feature-rich system for assembling, annotating, analysing, and visualizing bacterial genomes. Proksee accepts Illumina sequence reads as compressed FASTQ files or pre-assembled contigs in raw, FASTA, or GenBank format. Alternatively, users can supply a GenBank accession or a previously generated Proksee map in JSON format. Proksee then performs assembly (for raw sequence data), generates a graphical map, and provides an interface for customizing the map and launching further analysis jobs. Notable features of Proksee include unique and informative assembly metrics provided via a custom reference database of assemblies; a deeply integrated high-performance genome browser for viewing and comparing analysis results at individual base resolution (developed specifically for Proksee); an ever-growing list of embedded analysis tools whose results can be seamlessly added to the map or searched and explored in other formats; and the option to export graphical maps, analysis results, and log files for data sharing and research reproducibility. All these features are provided via a carefully designed multi-server cloud-based system that can easily scale to meet user demand and that ensures the web server is robust and responsive.

PlasMapper 3.0-a web server for generating, editing, annotating and visualizing publication quality plasmid maps

PlasMapper 3.0 is a web server that allows users to generate, edit, annotate and interactively visualize publication quality plasmid maps. Plasmid maps are used to plan, design, share and publish critical information about gene cloning experiments. PlasMapper 3.0 is the successor to PlasMapper 2.0 and offers many features found only in commercial plasmid mapping/editing packages. PlasMapper 3.0 allows users to paste or upload plasmid sequences as input or to upload existing plasmid maps from its large database of >2000 pre-annotated plasmids (PlasMapDB). This database can be searched by plasmid names, sequence features, restriction sites, preferred host organisms, and sequence length. PlasMapper 3.0 also supports the annotation of new or never-before-seen plasmids using its own feature database that contains common promoters, terminators, regulatory sequences, replication origins, selectable markers and other features found in most cloning plasmids. PlasMapper 3.0 has several interactive sequence editors/viewers that allow users to select and view plasmid regions, insert genes, modify restriction sites or perform codon optimization. The graphics for PlasMapper 3.0 have also been substantially upgraded. It now offers an interactive, full-color plasmid viewer/editor that allows users to zoom, rotate, re-color, linearize, circularize, edit annotated features and modify plasmid images or labels to improve the esthetic qualities of their plasmid map and textual displays. All the plasmid images and textual displays are downloadable in multiple formats. PlasMapper 3.0 is available online at https://plasmapper.ca.

Population Genetic Structure and Diversity of Cryptic Species of the Plant Genus Macrocarpaea (Gentianaceae) from the Tropical Andes

The Pleistocene climatic oscillations (PCO) that provoked several cycles of glacial-interglacial periods are thought to have profoundly affected species distribution, richness and diversity around the world. While the effect of the PCO on population dynamics at temperate latitudes is well known, considerable questions remain about its impact on the biodiversity of neotropical mountains. Here, we use amplified fragment length polymorphism molecular markers (AFLPs) to investigate the phylogeography and genetic structure of 13 plant species belonging to the gentian genus Macrocarpaea (Gentianaceae) in the tropical Andes. These woody herbs, shrubs or small trees show complex and potentially reticulated relationships, including cryptic species. We show that populations of M. xerantifulva in the dry system of the Rio Marañón in northern Peru have lower levels of genetic diversity compared to other sampled species. We suggest that this is due to a recent demographic bottleneck resulting from the contraction of the montane wet forests into refugia because of the expansion of the dry system into the valley during the glacial cycles of the PCO. This may imply that the ecosystems of different valleys of the Andes might have responded differently to the PCO.

The Evolution of Trait Disparity during the Radiation of the Plant Genus Macrocarpaea (Gentianaceae) in the Tropical Andes

The evolutionary processes responsible for the extraordinary diversity in the middle elevation montane forests of the Tropical Andes (MMF; 1000-3500 m) remain poorly understood. It is not clear whether adaptive divergence, niche conservatism or geographical processes were the main contributors to the radiation of the respective lineages occurring there. We investigated the evolutionary history of plant lineages in the MMF. We used the vascular plant genus Macrocarpaea (Gentianaceae) as a model, as it consists of 118 morphologically diverse species, a majority of which are endemic to the MMF. We used a time-calibrated molecular phylogeny and morphological and climatic data to compare a set of evolutionary scenarios of various levels of complexity in a phylogenetic comparative framework. In this paper, we show that the hypothesis of adaptive radiation for Macrocarpaea in the MMF is unlikely. The genus remained confined to the upper montane forests (UMF > 1800 m) during more than a half of its evolutionary history, possibly due to evolutionary constraints. Later, coinciding with the beginning of the Pleistocene (around 2.58 Ma), a phylogenetically derived (recently branching) clade, here referred to as the M. micrantha clade (25 species), successfully colonized and radiated in the lower montane forests (LMF < 1800 m). This colonization was accompanied by the evolution of a new leaf phenotype that is unique to the species of the M. micrantha clade that likely represents an adaptation to life in this new environment (adaptive zone). Therefore, our results suggest that niche conservatism and geographical processes have dominated most of the diversification history of Macrocarpaea, but that a rare adaptive divergence event allowed a transition into a new adaptive zone and enabled progressive radiation in this zone through geographical processes.

Biotic and abiotic factors shape arbuscular mycorrhizal fungal communities associated with the roots of the widespread fern Botrychium lunaria (Ophioglossaceae)

Arbuscular mycorrhizal fungi (AMF) play central roles in terrestrial ecosystems by interacting with both above and belowground communities as well as by influencing edaphic properties. The AMF communities associated with the roots of the fern Botrychium lunaria (Ophioglossaceae) were sampled in four transects at 2400 m a.s.l. in the Swiss Alps and analyzed using metabarcoding. Members of five Glomeromycota genera were identified across the 71 samples. Our analyses revealed the existence of a core microbiome composed of four abundant Glomus operational taxonomic units (OTUs), as well as a low OTU turnover between samples. The AMF communities were not spatially structured, which contrasts with most studies on AMF associated with angiosperms. pH, microbial connectivity and humus cover significantly shaped AMF beta diversity but only explained a minor fraction of variation in beta diversity. AMF OTUs associations were found to be significant by both cohesion and co-occurrence analyses, suggesting a role for fungus-fungus interactions in AMF community assembly. In particular, OTU co-occurrences were more frequent between different genera than among the same genus, rising the hypothesis of functional complementarity among the AMF associated to B. lunaria. Altogether, our results provide new insights into the ecology of fern symbionts in alpine grasslands.

High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data

Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.

Visualizing and comparing circular genomes using the CGView family of tools

Graphical genome maps are widely used to assess genome features and sequence characteristics. The CGView (Circular Genome Viewer) software family is a popular collection of tools for generating genome maps for bacteria, organelles and viruses. In this review, we describe the capabilities of the original CGView program along with those of subsequent companion applications, including the CGView Server and the CGView Comparison Tool. We also discuss GView, a graphical user interface-enabled rewrite of CGView, and the GView Server, which offers several integrated analyses for identifying shared or unique genome regions relative to a collection of comparison genomes. We conclude with some remarks about our current development efforts related to CGView aimed at adding new functionality while increasing ease of use.

A large interactive visual database of copy number variants discovered in taurine cattle

BACKGROUND

Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations.

RESULTS

CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle.

CONCLUSIONS

We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.

DrugBank 5.0: a major update to the DrugBank database for 2018

DrugBank (www.drugbank.ca) is a web-enabled database containing comprehensive molecular information about drugs, their mechanisms, their interactions and their targets. First described in 2006, DrugBank has continued to evolve over the past 12 years in response to marked improvements to web standards and changing needs for drug research and development. This year's update, DrugBank 5.0, represents the most significant upgrade to the database in more than 10 years. In many cases, existing data content has grown by 100% or more over the last update. For instance, the total number of investigational drugs in the database has grown by almost 300%, the number of drug-drug interactions has grown by nearly 600% and the number of SNP-associated drug effects has grown more than 3000%. Significant improvements have been made to the quantity, quality and consistency of drug indications, drug binding data as well as drug-drug and drug-food interactions. A great deal of brand new data have also been added to DrugBank 5.0. This includes information on the influence of hundreds of drugs on metabolite levels (pharmacometabolomics), gene expression levels (pharmacotranscriptomics) and protein expression levels (pharmacoprotoemics). New data have also been added on the status of hundreds of new drug clinical trials and existing drug repurposing trials. Many other important improvements in the content, interface and performance of the DrugBank website have been made and these should greatly enhance its ease of use, utility and potential applications in many areas of pharmacological research, pharmaceutical science and drug education.

nmrML: A Community Supported Open Data Standard for the Description, Storage, and Exchange of NMR Data

NMR is a widely used analytical technique with a growing number of repositories available. As a result, demands for a vendor-agnostic, open data format for long-term archiving of NMR data have emerged with the aim to ease and encourage sharing, comparison, and reuse of NMR data. Here we present nmrML, an open XML-based exchange and storage format for NMR spectral data. The nmrML format is intended to be fully compatible with existing NMR data for chemical, biochemical, and metabolomics experiments. nmrML can capture raw NMR data, spectral data acquisition parameters, and where available spectral metadata, such as chemical structures associated with spectral assignments. The nmrML format is compatible with pure-compound NMR data for reference spectral libraries as well as NMR data from complex biomixtures, i.e., metabolomics experiments. To facilitate format conversions, we provide nmrML converters for Bruker, JEOL and Agilent/Varian vendor formats. In addition, easy-to-use Web-based spectral viewing, processing, and spectral assignment tools that read and write nmrML have been developed. Software libraries and Web services for data validation are available for tool developers and end-users. The nmrML format has already been adopted for capturing and disseminating NMR data for small molecules by several open source data processing tools and metabolomics reference spectral libraries, e.g., serving as storage format for the MetaboLights data repository. The nmrML open access data standard has been endorsed by the Metabolomics Standards Initiative (MSI), and we here encourage user participation and feedback to increase usability and make it a successful standard.

Genetic architecture of gene expression underlying variation in host response to porcine reproductive and respiratory syndrome virus infection

It has been shown that inter-individual variation in host response to porcine reproductive and respiratory syndrome (PRRS) has a heritable component, yet little is known about the underlying genetic architecture of gene expression in response to PRRS virus (PRRSV) infection. Here, we integrated genome-wide genotype, gene expression, viremia level, and weight gain data to identify genetic polymorphisms that are associated with variation in inter-individual gene expression and response to PRRSV infection in pigs. RNA-seq analysis of peripheral blood samples collected just prior to experimental challenge (day 0) and at 4, 7, 11 and 14 days post infection from 44 pigs revealed 6,430 differentially expressed genes at one or more time points post infection compared to the day 0 baseline. We mapped genetic polymorphisms that were associated with inter-individual differences in expression at each day and found evidence of cis-acting expression quantitative trait loci (cis-eQTL) for 869 expressed genes (qval < 0.05). Associations between cis-eQTL markers and host response phenotypes using 383 pigs suggest that host genotype-dependent differences in expression of GBP5, GBP6, CCHCR1 and CMPK2 affect viremia levels or weight gain in response to PRRSV infection.

Heatmapper: web-enabled heat mapping for all

Heatmapper is a freely available web server that allows users to interactively visualize their data in the form of heat maps through an easy-to-use graphical interface. Unlike existing non-commercial heat map packages, which either lack graphical interfaces or are specialized for only one or two kinds of heat maps, Heatmapper is a versatile tool that allows users to easily create a wide variety of heat maps for many different data types and applications. More specifically, Heatmapper allows users to generate, cluster and visualize: (i) expression-based heat maps from transcriptomic, proteomic and metabolomic experiments; (ii) pairwise distance maps; (iii) correlation maps; (iv) image overlay heat maps; (v) latitude and longitude heat maps and (vi) geopolitical (choropleth) heat maps. Heatmapper offers a number of simple and intuitive customization options for facile adjustments to each heat map's appearance and plotting parameters. Heatmapper also allows users to interactively explore their numeric data values by hovering their cursor over each heat map cell, or by using a searchable/sortable data table view. Heat map data can be easily uploaded to Heatmapper in text, Excel or tab delimited formatted tables and the resulting heat map images can be easily downloaded in common formats including PNG, JPG and PDF. Heatmapper is designed to appeal to a wide range of users, including molecular biologists, structural biologists, microbiologists, epidemiologists, environmental scientists, agriculture/forestry scientists, fish and wildlife biologists, climatologists, geologists, educators and students. Heatmapper is available at http://www.heatmapper.ca.

PHASTER: a better, faster version of the PHAST phage search tool

PHASTER (PHAge Search Tool – Enhanced Release) is a significant upgrade to the popular PHAST web server for the rapid identification and annotation of prophage sequences within bacterial genomes and plasmids. Although the steps in the phage identification pipeline in PHASTER remain largely the same as in the original PHAST, numerous software improvements and significant hardware enhancements have now made PHASTER faster, more efficient, more visually appealing and much more user friendly. In particular, PHASTER is now 4.3× faster than PHAST when analyzing a typical bacterial genome. More specifically, software optimizations have made the backend of PHASTER 2.7X faster than PHAST, while the addition of 80 CPUs to the PHASTER compute cluster are responsible for the remaining speed-up. PHASTER can now process a typical bacterial genome in 3 min from the raw sequence alone, or in 1.5 min when given a pre-annotated GenBank file. A number of other optimizations have also been implemented, including automated algorithms to reduce the size and redundancy of PHASTER's databases, improvements in handling multiple (metagenomic) queries and higher user traffic, along with the ability to perform automated look-ups against 14 000 previously PHAST/PHASTER annotated bacterial genomes (which can lead to complete phage annotations in seconds as opposed to minutes). PHASTER's web interface has also been entirely rewritten. A new graphical genome browser has been added, gene/genome visualization tools have been improved, and the graphical interface is now more modern, robust and user-friendly. PHASTER is available online at www.phaster.ca.

ECMDB 2.0: A richer resource for understanding the biochemistry of E. coli

ECMDB or the Escherichia coli Metabolome Database (http://www.ecmdb.ca) is a comprehensive database containing detailed information about the genome and metabolome of E. coli (K-12). First released in 2012, the ECMDB has undergone substantial expansion and many modifications over the past 4 years. This manuscript describes the most recent version of ECMDB (ECMDB 2.0). In particular, it provides a comprehensive update of the database that was previously described in the 2013 NAR Database Issue and details many of the additions and improvements made to the ECMDB over that time. Some of the most important or significant enhancements include a 13-fold increase in the number of metabolic pathway diagrams (from 125 to 1650), a 3-fold increase in the number of compounds linked to pathways (from 1058 to 3280), the addition of dozens of operon/metabolite signalling pathways, a 44% increase in the number of compounds in the database (from 2610 to 3760), a 7-fold increase in the number of compounds with NMR or MS spectra (from 412 to 3261) and a massive increase in the number of external links to other E. coli or chemical resources. These additions, along with many other enhancements aimed at improving the ease or speed of querying, searching and viewing the data within ECMDB should greatly facilitate the understanding of not only the metabolism of E. coli, but also allow the in-depth exploration of its extensive metabolic networks, its many signalling pathways and its essential biochemistry.

Accurate, fully-automated NMR spectral profiling for metabolomics

Many diseases cause significant changes to the concentrations of small molecules (a.k.a. metabolites) that appear in a person’s biofluids, which means such diseases can often be readily detected from a person’s “metabolic profile"—i.e., the list of concentrations of those metabolites. This information can be extracted from a biofluids Nuclear Magnetic Resonance (NMR) spectrum. However, due to its complexity, NMR spectral profiling has remained manual, resulting in slow, expensive and error-prone procedures that have hindered clinical and industrial adoption of metabolomics via NMR. This paper presents a system, BAYESIL, which can quickly, accurately, and autonomously produce a person’s metabolic profile. Given a 1D ¹HNMR spectrum of a complex biofluid (specifically serum or cerebrospinal fluid), BAYESIL can automatically determine the metabolic profile. This requires first performing several spectral processing steps, then matching the resulting spectrum against a reference compound library, which contains the “signatures” of each relevant metabolite. BAYESIL views spectral matching as an inference problem within a probabilistic graphical model that rapidly approximates the most probable metabolic profile. Our extensive studies on a diverse set of complex mixtures including real biological samples (serum and CSF), defined mixtures and realistic computer generated spectra; involving > 50 compounds, show that BAYESIL can autonomously find the concentration of NMR-detectable metabolites accurately (~ 90% correct identification and ~ 10% quantification error), in less than 5 minutes on a single CPU. These results demonstrate that BAYESIL is the first fully-automatic publicly-accessible system that provides quantitative NMR spectral profiling effectively—with an accuracy on these biofluids that meets or exceeds the performance of trained experts. We anticipate this tool will usher in high-throughput metabolomics and enable a wealth of new applications of NMR in clinical settings. BAYESIL is accessible at http://www.bayesil.ca.

An integrated platform for bovine DNA methylome analysis suitable for small samples

Background

Oocytes and early embryos contain minute amounts of DNA, RNA and proteins, making the study of early mammalian development highly challenging. The study of the embryo epigenome, in particular the DNA methylome, has been made accessible thanks to the possibility of amplifying specific sequences according to their initial methylation status. This paper describes a novel platform dedicated to the genome-wide study of bovine DNA methylation, including a complete pipeline for data analysis and visualization. The platform allows processing and integrating of DNA methylome and transcriptome data from the same sample. Procedures were optimized for genome-wide analysis of 10 ng of DNA (10 bovine blastocysts). Bovine sperm and blastocysts were compared as a test of platform capability.

Results

The hypermethylation of bovine sperm DNA compared to the embryo genome was confirmed. Differentially methylated regions were distributed across various classes of bovine sperm genomic feature including primarily promoter, intronic and exonic regions, non-CpG-island regions (shore, shelf and open-sea) and CpG islands with low-to-intermediate CpG density. The blastocyst genome bore more methylation marks than sperm DNA only in CpG islands with high CpG density. Long-terminal-repeat retrotransposons (LTR), LINE and SINE were more methylated in sperm DNA, as were low-complexity repetitive elements in blastocysts.

Conclusions

This is the first early embryo compatible genome-wide epigenetics platform for bovine. Such platforms should improve the study of the potential epigenetic risks of assisted reproductive technologies (ART), the establishment sequence of embryonic cell lines and potential deviations in both gene expression and DNA methylation capable of having long-term impact.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-451) contains supplementary material, which is available to authorized users.

Molecular phylogenetics supports widespread cryptic species in moonworts (Botrychium s.s., Ophioglossaceae)

PREMISE OF THE STUDY

Previous phylogenetic studies of moonworts (Botrychium sensu stricto (s.s.)) included few taxa from outside of North America. This low geographical representation limited interpretations of relationships of this group rich in cryptic species. With 18 out of 30 species in the genus being polyploid, understanding their evolutionary history remains a major challenge.

METHODS

A new molecular phylogeny was reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI) analyses based on multiple accessions of the most wide-ranging Arctic taxa of Botrychium in North America and Europe using three noncoding plastid DNA regions (psbA-trnH(GUG), trnL(UAA)-trnF(GAA) intergenic spacer, and rpL16 intron).

KEY RESULTS

The new phylogeny confirms the identity of several recently described species and proposed new taxa. Nine subclades are newly identified within the two major clades in Botrychium s.s.: Lanceolatum and Lunaria. Chloroplast DNA was variable enough to separate morphologically cryptic species in the Lunaria clade. On the contrary, much less variation is seen within the morphologically variable Lanceolatum clade despite sampling over the same broad geographic range. The chloroplast region psbA-trnH(GUG) is identified as an efficient DNA barcode for the identification of cryptic taxa in Botrychium s.s.

CONCLUSIONS

The combined increase in species representation, samples from throughout the geographic range of each species, and sequencing of multiple plastid DNA regions supports morphologically cryptic species in Botrychium s.s.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Development of a porcine (Sus scofa) embryo-specific microarray: array annotation and validation

BACKGROUND

The domestic pig is an important livestock species and there is strong interest in the factors that affect the development of viable embryos and offspring in this species. A limited understanding of the molecular mechanisms involved in early embryonic development has inhibited our ability to fully elucidate these factors. Next generation deep sequencing and microarray technologies are powerful tools for delineation of molecular pathways involved in the developing embryo.

RESULTS

Here we present the development of a porcine-embryo-specific microarray platform created from a large expressed sequence tag (EST) analysis generated by Roche/454 next-generation sequencing of cDNAs constructed from critical stages of in vivo or in vitro porcine preimplantation embryos. Two cDNA libraries constructed from in vitro and in vivo produced preimplantation porcine embryos were normalized and sequenced using 454 Titanium pyrosequencing technology. Over one million high-quality EST sequences were obtained and used to develop the EMbryogene Porcine Version 1 (EMPV1) microarray composed of 43,795 probes. Based on an initial probe sequence annotation, the EMPV1 features 17,409 protein-coding, 473 pseudogenes, 46 retrotransposed, 2,359 non-coding RNA, 4,121 splice variants in 2,862 genes and a total of 12,324 Novel Transcript Regions (NTR). After re-annotation, the total unique genes increased from 11,961 to 16,281 and 1.9% of them belonged to a large olfactory receptor (OR) gene family. Quality control on the EMPV1 was performed and revealed an even distribution of ten clusters of spiked-in control spots and array to array (dye-swap) correlation was 0.97.

CONCLUSIONS

Using next-generation deep sequencing we have produced a large EST dataset to allow for the selection of probe sequences for the development of the EMPV1 microarray platform. The quality of this embryo-specific array was confirmed with a high-level of reproducibility using current Agilent microarray technology. With more than an estimated 20,000 unique genes represented on the EMPV1, this platform will provide the foundation for future research into the in vivo and in vitro factors that affect the viability of porcine embryos, as well as the effects of these factors on the live offspring that result from these embryos.

Comparing thousands of circular genomes using the CGView Comparison Tool

BACKGROUND

Continued sequencing efforts coupled with advances in sequencing technology will lead to the completion of a vast number of small genomes. Whole-genome comparisons represent an important part of the analysis of any new genome sequence, as they can provide a better understanding of the biology and evolution of the source organism. Visualization of the results is important, as it allows information from a variety of sources to be integrated and interpreted. However, existing graphical comparison tools lack features needed for efficiently comparing a new genome to hundreds or thousands of existing sequences. Moreover, existing tools are limited in terms of the types of comparisons that can be performed, the extent to which the output can be customized, and the ease with which the entire process can be automated.

RESULTS

The CGView Comparison Tool (CCT) is a package for visually comparing bacterial, plasmid, chloroplast, or mitochondrial sequences of interest to existing genomes or sequence collections. The comparisons are conducted using BLAST, and the BLAST results are presented in the form of graphical maps that can also show sequence features, gene and protein names, COG (Clusters of Orthologous Groups of proteins) category assignments, and sequence composition characteristics. CCT can generate maps in a variety of sizes, including 400 Megapixel maps suitable for posters. Comparisons can be conducted within a particular species or genus, or all available genomes can be used. The entire map creation process, from downloading sequences to redrawing zoomed maps, can be completed easily using scripts included with the CCT. User-defined features or analysis results can be included on maps, and maps can be extensively customized. To simplify program setup, a CCT virtual machine that includes all dependencies preinstalled is available. Detailed tutorials illustrating the use of CCT are included with the CCT documentation.

CONCLUSION

CCT can be used to visually compare a reference sequence to thousands of existing genomes or sequence collections (next-generation sequencing reads for example) on a standard desktop computer. It provides analysis and visualization functionality not available in any existing circular genome visualization tool. By visually presenting sequence conservation information along with functional classifications and sequence composition characteristics, CCT can be a useful tool for identifying rapidly evolving or novel sequences, horizontally transferred sequences, or unusual functional properties in newly sequenced genomes. CCT is freely available for download at http://stothard.afns.ualberta.ca/downloads/CCT/.

Combining resources to obtain a comprehensive survey of the bovine embryo transcriptome through deep sequencing and microarrays

While most assisted reproductive technologies (ART) are considered routine for the reproduction of species of economical importance, such as the bovine, the impact of these manipulations on the developing embryo remains largely unknown. In an effort to obtain a comprehensive survey of the bovine embryo transcriptome and how it is modified by ART, resources were combined to design an embryo-specific microarray. Close to one million high-quality reads were produced from subtracted bovine embryo libraries using Roche 454 Titanium deep sequencing technology, which enabled the creation of an augmented bovine genome catalog. This catalog was enriched with bovine embryo transcripts, and included newly discovered indel type and 3'UTR variants. Using this augmented bovine genome catalog, the EmbryoGENE Bovine Microarray was designed and is composed of a total of 42,242 probes, including 21,139 known reference genes; 9,322 probes for novel transcribed regions (NTRs); 3,677 alternatively spliced exons; 3,353 3'-tiling probes; and 3,723 controls. A suite of bioinformatics tools was also developed to facilitate microrarray data analysis and database creation; it includes a quality control module, a Laboratory Information Management System (LIMS) and microarray analysis software. Results obtained during this study have already led to the identification of differentially expressed blastocyst targets, NTRs, splice variants of the indel type, and 3'UTR variants. We were able to confirm microarray results by real-time PCR, indicating that the EmbryoGENE bovine microarray has the power to detect physiologically relevant changes in gene expression.

In-depth annotation of SNPs arising from resequencing projects using NGS-SNP

Summary: NGS-SNP is a collection of command-line scripts for providing rich annotations for SNPs identified by the sequencing of whole genomes from any organism with reference sequences in Ensembl. Included among the annotations, several of which are not available from any existing SNP annotation tools, are the results of detailed comparisons with orthologous sequences. These comparisons can, for example, identify SNPs that affect conserved residues, or alter residues or genes linked to phenotypes in another species.

Availability: NGS-SNP is available both as a set of scripts and as a virtual machine. The virtual machine consists of a Linux operating system with all the NGS-SNP dependencies pre-installed. The source code and virtual machine are freely available for download at http://stothard.afns.ualberta.ca/downloads/NGS-SNP/.

Contact:stothard@ualberta.ca

Supplementary information:Supplementary data are available at Bioinformatics online.

Comparative RH maps of the river buffalo and bovine Y chromosomes

Radiation hybrid maps were constructed for river buffalo and cattle Y chromosomes. A total of 41 cattle-derived Y-chromosome molecular markers were selected and tested with 2 previously described 5,000-rad whole-genome radiation hybrid (RH) panels (river buffalo - BBURH(5000) and cattle - BTARH(5000)) for generation of maps. Among the initial 41 selected markers, a subset of 26 markers generated PCR products suitable for scoring with the BBURH(5000) panel. Of these, 19 markers (73%) were distributed in 1 linkage group spanning 341.3 cR. Retention frequencies (RF) for individual markers ranged from 17.8% for SMCY to 56.7% for BTY1, with an average RF of 37.6%. From the selected markers, 37 generated reliable scores using the BTARH(5000) panel. The newly constructed BTAY RH map contains 28 markers distributed within 1 linkage group. Twenty-four of these markers had been previously mapped on BTAY using a 7,000-rad cattle-hamster WG-RH panel and 4 markers were mapped for the first time (ZFY, SeqRep, RepSeqS4 and BTY1). The length of the BTAY RH map was estimated to be 602.4 cR. Retention frequencies for individual mapped markers ranged from 10% (INRA126) to 63.3% (SeqRep), with an average RF of 35.3%. RH marker positions along the Y chromosome were compared between BBUY and BTAY, which revealed differences in the order of some of the markers. The BBUY pseudoautosomal region (PAR) is delineated by 3 BTAY PAR markers (MAF45, TGLA325 and UMN2008). These markers are telomeric in both species but are not found in the same order. Here we have demonstrated the effective use of bovine Y chromosome markers for the development of the first BBUY RH map. Likewise, these set of markers can be used for comparative assessment of Y chromosomes in other members of the Bovidae family.

A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls

The detection and mapping of genetic markers linked to quantitative trait loci (QTL) can be utilized to enhance genetic improvement of livestock populations. With the completion of the bovine genome sequence assembly, single nucleotide polymorphisms (SNP) assays spanning the whole bovine genome and research work on large scale identification, validation and analysis of genotypic variation in cattle has become possible. The objective of the present study was to perform a whole genome scan to identify and map QTL affecting milk production traits and somatic cell scores using linkage disequilibrium (LD) regression and 1536 SNP markers. Three and 18 SNP were found to be associated with only milk yield (MY) at a genome and chromosome wise significance (p < 0.05) level respectively. Among the 21 significant SNP, 16 were in a region reported to have QTL for MY in other dairy cattle populations and while the rest five were new QTL finding. Four SNP out of 21 are significant for the milk production traits (MY, fat yield, protein yield (PY), and milk contents) in the present study. Six and nine SNP were associated with PY at a genome and chromosome wise significant (p < 0.05) level respectively. Three and 17 SNP were found to be associated with FY at a genome and chromosome wise significant (p < 0.05) level. Five and seven SNP were mapped with somatic cell score at a genome and chromosome wise significant (p < 0.05) level respectively. The results of this study have revealed QTL for MY, PY, protein percentage, FY, fat percentage, somatic cell score and persistency of milk in the Canadian dairy cattle population. The chromosome regions identified in this study should be further investigated to potentially identify the causative mutations underlying the QTL.

Characterization of bovine miRNAs by sequencing and bioinformatics analysis

Background

MicroRNAs (miRNAs) are a family of ~22 nucleotide small RNA molecules which regulate gene expression by fully or partially binding to their complementary sequences in mRNAs or promoters. A large number of miRNAs and their expression patterns have been reported in human, mouse and rat. However, miRNAs and their expression patterns in live stock species such as beef cattle are not well studied.

Results

We constructed and sequenced small-RNA libraries to yield a total of 13,541 small-RNA sequences from 11 bovine tissues including brain, subcutaneous fat, muscle, liver, kidney, spleen and thymus. In total, 228 miRNAs including 29 novel miRNA candidates were identified. Of the 199 miRNAs, 101 have been previously reported as bovine miRNAs and the other 98 are bovine orthologs of known miRNAs that have been identified in at least one other mammalian species. Of the 29 novel miRNA candidates, 17 appeared at this point in time to be bovine specific, while the remaining 12 had evidence of evolutionary conservation in other mammalian species. Five miRNAs (miR-23a, -23b, -99a, -125b and -126-5p) were very abundant across the 11 tissues, accounting for 44.3% of all small RNA sequences. The expression analysis of selected miRNAs using qRT-PCR also showed that miR-26a and -99a were highly expressed in all tissues, while miR-122 and miR-133a were predominantly expressed in liver and muscle, respectively.

Conclusion

The miRNA expression patterns among 11 tissues from beef cattle revealed that most miRNAs were ubiquitously expressed in all tissues, while only a few miRNAs were tissue specific. Only 60% miRNAs in this study were found to display strand bias, suggesting that there are some key factors for mature miRNA selection other than internal stability. Most bovine miRNAs are highly conserved in other three mammalian species, indicating that these miRNAs may have a role in different species that are potential molecular markers for evolution.

The genome sequence of taurine cattle: a window to ruminant biology and evolution

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle

Background

The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH₅₀₀₀) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species.

Results

A total of 3990 markers were typed on the BBURH₅₀₀₀ panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD ≥ 8. The estimated total length of the WG-RH map is 36,933 cR₅₀₀₀. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score ≥ 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0).

Conclusion

We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.

A whole-genome scan to map quantitative trait loci for conformation and functional traits in Canadian Holstein bulls

Genetic improvement of livestock populations can be achieved through detection and mapping of genetic markers linked to quantitative trait loci (QTL). With the completion of the bovine genome sequence assembly, single nucleotide polymorphism (SNP) assays spanning the whole bovine genome and research work on large-scale identification, validation, and analysis of genotypic variation in cattle has become possible. A total of 462 Canadian Holstein Bulls were used to test the association between SNP and QTL. Single locus linkage disequilibrium regression model was implemented to perform a whole genome scan to identify and map QTL affecting conformation and functional traits. One thousand five hundred thirty-six SNP markers from introns and exons of potential QTL regions for economically important traits across the bovine genome were selected for association analysis. A total of 45 and 151 SNP were found to be associated with 17 conformation and functional traits at a genome- and chromosome-wise significance level, respectively. Among the 196 significant SNP, 169 of them are newly detected in this study, whereas 27 of them have been reported in previous literature and 161 of these were located in genes and are worth further investigating to potentially identify the causative mutations underlying the QTL. The single locus linkage disequilibrium regression method using SNP marker genotypes has proven to be a successful methodology for detecting and mapping QTL in dairy cattle populations.

The CGView Server: a comparative genomics tool for circular genomes

The CGView Server generates graphical maps of circular genomes that show sequence features, base composition plots, analysis results and sequence similarity plots. Sequences can be supplied in raw, FASTA, GenBank or EMBL format. Additional feature or analysis information can be submitted in the form of GFF (General Feature Format) files. The server uses BLAST to compare the primary sequence to up to three comparison genomes or sequence sets. The BLAST results and feature information are converted to a graphical map showing the entire sequence, or an expanded and more detailed view of a region of interest. Several options are included to control which types of features are displayed and how the features are drawn. The CGView Server can be used to visualize features associated with any bacterial, plasmid, chloroplast or mitochondrial genome, and can aid in the identification of conserved genome segments, instances of horizontal gene transfer, and differences in gene copy number. Because a collection of sequences can be used in place of a comparison genome, maps can also be used to visualize regions of a known genome covered by newly obtained sequence reads. The CGView Server can be accessed at http://stothard.afns.ualberta.ca/cgview_server/

Identification of a novel immunosubversion mechanism mediated by a virologue of the B-lymphocyte receptor TACI

TACI (transmembrane activator and calcium modulator and cyclophilin ligand [CAML] interactor) is a part of a novel network of ligands and receptors involved in B-cell survival and isotype switching. The TACI protein mediates its effects through CAML, an endoplasmic reticulum (ER)-localized protein that controls Ca(2+) efflux. The adenovirus E3-6.7K protein prevents inflammatory responses and also confers resistance from a variety of apoptotic stimuli and maintains ER Ca(2+) homeostasis; however, the mechanism of action is unknown. Here, we provide evidence that E3-6.7K shares sequence homology with TACI and inhibits apoptosis and ER Ca(2+) efflux through an interaction with CAML, a Ca(2+)-modulating protein. We demonstrate a direct interaction between E3-6.7K and CAML and reveal that the two proteins colocalize in an ER-like compartment. Furthermore, the interaction between the two proteins is localized to the N-terminal domain of CAML and to a 22-amino-acid region near the C terminus of E3-6.7K termed the CAML-binding domain (CBD). Mutational analysis of the CBD showed that an interaction with CAML is required for E3-6.7K to inhibit thapsigargin-induced apoptosis and ER Ca(2+) efflux. E3-6.7K appears to be the first virologue of TACI to be identified. It targets CAML in a novel immunosubversive mechanism to alter ER Ca(2+) homeostasis, which consequently inhibits inflammation and protects infected cells from apoptosis.

The adenovirus E3-6.7K protein adopts diverse membrane topologies following posttranslational translocation

The E3 region of adenovirus codes for several membrane proteins, most of which are involved in immune evasion and prevention of host cell apoptosis. We explored the topology and targeting mechanisms of E3-6.7K, the most recently described member of this group, by using an in vitro translation system supplemented with microsomes. Here, we present evidence that E3-6.7K, one of the smallest signal-anchor proteins known, translocates across the membrane of the endoplasmic reticulum in a posttranslational, ribosome-independent, yet ATP-dependent manner, reminiscent of the translocation of tail-anchored proteins. Our analysis also demonstrated that E3-6.7K could achieve several distinct topological fates. In addition to the previously postulated type III orientation (N-luminal/C-cytoplasmic, termed NtmE3-6.7K), we detected a tail-anchored form adopting the opposite orientation (N-cytoplasmic/C-luminal, termed CtmE3-6.7K) as well as the possibility of a fully translocated form (N and C termini are both translocated, termed NCE3-6.7K). Due to the translocation of a positively charged domain, both the CtmE3-6.7K and NCE3-6.7K topologies of E3-6.7K constitute exceptions to the "positive inside" rule. The NtmE3-6.7K and NCE3-6.7K are the first examples of posttranslationally translocated proteins in higher eukaryotes that are not tail anchored. Distinct topological forms were also found in transfected cells, as both N and C termini of E3-6.7K were detected on the extracellular surface of transfected cells. The demonstration of unexpected topological forms and translocation mechanisms for E3-6.7K defies conventional thinking about membrane protein topogenesis and advises that both the mode of targeting and topology of signal-anchor proteins should be determined experimentally.

Adenovirus E3-6.7K maintains calcium homeostasis and prevents apoptosis and arachidonic acid release

E3-6.7K is a small and hydrophobic membrane glycoprotein encoded by the E3 region of subgroup C adenovirus. Recently, E3-6.7K has been shown to be required for the downregulation of tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL) receptors by the adenovirus E3/10.4K and E3/14.5K complex of proteins. We demonstrate here that E3-6.7K has additional protective roles, independent of other virus proteins. In transfected Jurkat T-cell lymphoma cells, E3-6.7K was found to maintain endoplasmic reticulum-Ca(2+) homeostasis and inhibit the induction of apoptosis by thapsigargin. The presence of E3-6.7K also lead to a reduction in the TNF-induced release of arachidonic acid from transfected U937 human histiocytic lymphoma cells. In addition, E3-6.7K protected cells against apoptosis induced through Fas, TNF receptor, and TRAIL receptors. Therefore, E3-6.7K confers a wide range of protective effects against both Ca(2+) flux-induced and death receptor-mediated apoptosis.

Human CARD4 protein is a novel CED-4/Apaf-1 cell death family member that activates NF-kappaB

The nematode CED-4 protein and its human homolog Apaf-1 play a central role in apoptosis by functioning as direct activators of death-inducing caspases. A novel human CED-4/Apaf-1 family member called CARD4 was identified that has a domain structure strikingly similar to the cytoplasmic, receptor-like proteins that mediate disease resistance in plants. CARD4 interacted with the serine-threonine kinase RICK and potently induced NF-kappaB activity through TRAF-6 and NIK signaling molecules. In addition, coexpression of CARD4 augmented caspase-9-induced apoptosis. Thus, CARD4 coordinates downstream NF-kappaB and apoptotic signaling pathways and may be a component of the host innate immune response.

Epitope mapping and immunoneutralization of recombinant human stem-cell factor

The epitope regions of three anti-[stem-cell factor (SCF)]g have been mapped by characterization of immunoreactivities against truncated forms of SCF in immunoblots and against synthetic peptides in solution-phase competition ELISA. Two of the antibodies, mAb 7H6 and mAb 8H7A, were raised against Escherichia coli-derived human SCF-(1-164) while the third, polyclonal antibody (pAb) 1337, was raised against a peptide corresponding to residues 3-31 of human SCF. The epitopes of mAbs 7H6 and 8H7A have been mapped to residues 61-95 and 95-110, respectively. The epitope of pAb 1337 has been mapped to residues 21-31. The ability of the anti-SCF Ig to recognize E. coli-derived human SCF presented in various formats, i.e. partially denatured (fixed in standard ELISA or on a western blot) or native (in solution), was studied, mAb 7H6 recognized its epitope in partially denatured or native SCF with equally high affinity, while mAb 8H7A and pAb 1337 recognized their epitopes only when SCF was at least partially denatured, mAb 7H6 was found to neutralize in vitro SCF-mediated cell proliferation and SCF binding to its receptor, when present in equimolar concentrations relative to the ligand, suggesting that the epitope region is functionally significant. Evidence that the mAb 7H6 epitope is represented by discontinuous regions (residues within sequences 61-65 and 91-95 are critically involved) is presented. The observation that the mAb 7H6 epitope is discontinuous has implications for the structure of SCF.

Relationship of sexual and physical abuse to pain description, coping, psychological distress, and health-care utilization in a chronic pain sample

OBJECTIVE

We investigated the association of physical/sexual abuse to pain description, coping, psychological distress, and health-care utilization in a heterogeneous sample of chronic pain patients.

DESIGN

A cross-sectional, retrospective design was used. Patients were categorized as abused (n = 22) or nonabused (n = 58) based on responses to a valid and reliable sexual/physical abuse questionnaire.

MAIN OUTCOME MEASURES

Pain description (Visual Analog Scale measures of pain intensity and frequency, and the McGill Pain Questionnaire); coping ability and attributional style (Functional Interference Estimate, Self-Control Schedule, Pain Locus of Control Scale); psychological distress (SCL-90-R Global Severity Index); and a Health-care utilization measure.

RESULTS AND CONCLUSIONS

No differences between abused and nonabused groups were found for the pain description or functional interference variables. Compared to the nonabused group, the abused group had significantly lower Self-Control Schedule scores, higher Pain Locus of Control Scale Chance Factor scores, and higher SCL-90-R Global Severity Index scores and was more likely to use the emergency room for pain symptoms. These results replicate the findings of previous studies of the effects of abuse in more specific pain samples and underscore the importance of assessment of abuse in patients with chronic pain. The data suggest that interventions which involve coping-skills training or self-control management of pain may be affected by an abuse history via reduced perceptions of efficacy, resourcefulness, and beliefs that external variables are responsible for pain.

Why do women exercise? Factor analysis and further validation of the Reasons for Exercise Inventory

This investigation examined the factor structure and correlates of the Reasons for Exercise Inventory among 101 exercising women. Subjects completed the 24-item inventory (with one added item), reported their weekly frequency of exercise, and completed two standardized body-image measures. Factor analysis indicated that, with minor modifications, the instrument has an internally consistent structure with four factors of Appearance/Weight Management, Fitness/Health Management, Stress/Mood Management, and Socializing. Appearance/Weight Management was associated with a more negative body image independent of actual body mass and was the only motive related to self-reported frequency of exercise.