Identification and quantification of a new family of peptide endocannabinoids (Pepcans) showing negative allosteric modulation at CB1 receptors

The α-hemoglobin-derived dodecapeptide RVD-hemopressin (RVDPVNFKLLSH) has been proposed to be an endogenous agonist for the cannabinoid receptor type 1 (CB(1)). To study this peptide, we have raised mAbs against its C-terminal part. Using an immunoaffinity mass spectrometry approach, a whole family of N-terminally extended peptides in addition to RVD-Hpα were identified in rodent brain extracts and human and mouse plasma. We designated these peptides Pepcan-12 (RVDPVNFKLLSH) to Pepcan-23 (SALSDLHAHKLRVDPVNFKLLSH), referring to peptide length. The most abundant Pepcans found in the brain were tested for CB(1) receptor binding. In the classical radioligand displacement assay, Pepcan-12 was the most efficacious ligand but only partially displaced both [(3)H]CP55,940 and [(3)H]WIN55,212-2. The data were fitted with the allosteric ternary complex model, revealing a cooperativity factor value α < 1, thus indicating a negative allosteric modulation. Dissociation kinetic studies of [(3)H]CP55,940 in the absence and presence of Pepcan-12 confirmed these results by showing increased dissociation rate constants induced by Pepcan-12. A fluorescently labeled Pepcan-12 analog was synthesized to investigate the binding to CB(1) receptors. Competition binding studies revealed K(i) values of several Pepcans in the nanomolar range. Accordingly, using competitive ELISA, we found low nanomolar concentrations of Pepcans in human plasma and ∼100 pmol/g in mouse brain. Surprisingly, Pepcan-12 exhibited potent negative allosteric modulation of the orthosteric agonist-induced cAMP accumulation, [(35)S]GTPγS binding, and CB(1) receptor internalization. Pepcans are the first endogenous allosteric modulators identified for CB(1) receptors. Given their abundance in the brain, Pepcans could play an important physiological role in modulating endocannabinoid signaling.

Radiation metabolomics. 5. Identification of urinary biomarkers of ionizing radiation exposure in nonhuman primates by mass spectrometry-based metabolomics

Mass spectrometry-based metabolomics has previously demonstrated utility for identifying biomarkers of ionizing radiation exposure in cellular, mouse and rat in vivo radiation models. To provide a valuable link from small laboratory rodents to humans, γ-radiation-induced urinary biomarkers were investigated using a nonhuman primate total-body-irradiation model. Mass spectrometry-based metabolomics approaches were applied to determine whether biomarkers could be identified, as well as the previously discovered rodent biomarkers of γ radiation. Ultra-performance liquid chromatography-electrospray ionization quadrupole time-of-flight mass spectrometry analysis was carried out on a time course of clean-catch urine samples collected from nonhuman primates (n = 6 per cohort) exposed to sham, 1.0, 3.5, 6.5 or 8.5 Gy doses of (60)Co γ ray (∼0.55 Gy/min) ionizing radiation. By multivariate data analysis, 13 biomarkers of radiation were discovered: N-acetyltaurine, isethionic acid, taurine, xanthine, hypoxanthine, uric acid, creatine, creatinine, tyrosol sulfate, 3-hydroxytyrosol sulfate, tyramine sulfate, N-acetylserotonin sulfate, and adipic acid. N-Acetyltaurine, isethionic acid, and taurine had previously been identified in rats, and taurine and xanthine in mice after ionizing radiation exposure. Mass spectrometry-based metabolomics has thus successfully revealed and verified urinary biomarkers of ionizing radiation exposure in the nonhuman primate for the first time, which indicates possible mechanisms for ionizing radiation injury.

Identification of SNPs in the cystic fibrosis interactome influencing pulmonary progression in cystic fibrosis

There is growing evidence that the great phenotypic variability in patients with cystic fibrosis (CF) not only depends on the genotype, but apart from a combination of environmental and stochastic factors predominantly also on modifier gene effects. It has been proposed that genes interacting with CF transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) are potential modifiers. Therefore, we assessed the impact of single-nucleotide polymorphisms (SNPs) of several of these interacters on CF disease outcome. SNPs that potentially alter gene function were genotyped in 95 well-characterized p.Phe508del homozygous CF patients. Linear mixed-effect model analysis was used to assess the relationship between sequence variants and the repeated measurements of lung function parameters. In total, we genotyped 72 SNPs in 10 genes. Twenty-five SNPs were used for statistical analysis, where we found strong associations for one SNP in PPP2R4 with the lung clearance index (P ≤ 0.01), the specific effective airway resistance (P ≤ 0.005) and the forced expiratory volume in 1 s (P ≤ 0.005). In addition, we identified one SNP in SNAP23 to be significantly associated with three lung function parameters as well as one SNP in PPP2R1A and three in KRT19 to show a significant influence on one lung function parameter each. Our findings indicate that direct interacters with CFTR, such as SNAP23, PPP2R4 and PPP2R1A, may modify the residual function of p.Phe508del-CFTR while variants in KRT19 may modulate the amount of p.Phe508del-CFTR at the apical membrane and consequently modify CF disease.

Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace I: carcass traits

We report the identification of quantitative trait loci (QTL) affecting carcass composition, carcass length, fat deposition and lean meat content using a genome scan across 462 animals from a combined intercross and backcross between Hampshire and Landrace pigs. Data were analysed using multiple linear regression fitting additive and dominance effects. This model was compared with a model including a parent-of-origin effect to spot evidence of imprinting. Several precisely defined muscle phenotypes were measured in order to dissect body composition in more detail. Three significant QTL were detected in the study at the 1% genome-wide level, and twelve significant QTL were detected at the 5% genome-wide level. These QTL comprise loci affecting fat deposition and lean meat content on SSC1, 4, 9, 10, 13 and 16, a locus on SSC2 affecting the ratio between weight of meat and bone in back and weight of meat and bone in ham and two loci affecting carcass length on SSC12 and 17. The well-defined phenotypes in this study enabled us to detect QTL for sizes of individual muscles and to obtain information of relevance for the description of the complexity underlying other carcass traits.

Characterization of retinal damage in the episcleral vein cauterization rat glaucoma model

Episcleral vein cauterization (EVC) is used in rats to generate a glaucoma model with high intraocular pressure (IOP). The long-term retinal damage in this glaucoma model, however, has not been accurately quantified. We report the location and amount of retinal ganglion cell (RGC) damage caused by (EVC) induced IOP elevation in two rat strains. IOP was raised in one eye of Wistar (N = 5) and Brown-Norway(B-N)(N = 7) rats by EVC and monitored monthly until IOP in contralateral eyes equalized at 5 months post-surgery. Animals were maintained for 3.5-4.5 additional months. B-N rats (N = 7) that had no EVC served as controls for this strain. Scotopic flash ERGs were recorded at baseline and just prior to euthanasia. Automated counts of all retrogradely labeled RGCs in retinal flat-mounts were determined and compared between contralateral eyes. RGC density maps were constructed and RGC size distribution was determined. Oscillatory potentials in the group of eyes which had elevated IOP were decreased at the time of euthanasia, when IOP had returned to normal. The group of normal B-N rats had similar RGC counts between contralateral eyes. In the experimental group the mean number of RGCs was not significantly different between control and experimental eyes, but 1 of 5 Wistar and 2 of 7 B-N experimental eyes had at least 30% fewer RGCs than contralateral control eyes. Total retinal area in B-N experimental eyes was higher compared to contralateral eyes. Cumulative IOP exposure of the experimental eyes was modestly correlated with RGC loss while oscillatory potentials appeared to be inversely related to RGC loss. In retinas with extensive (> 30% RGC loss) but not complete damage, smaller cells were preserved better than larger ones. The above results indicate that RGC loss in both Wistar and B-N strains is variable after a prolonged elevation of IOP via EVC. Such variability despite equivalent IOP levels and ERG abnormalities, suggests unknown factors that can protect IOP-stressed RGCs. Identification and enhancement of such factors could prove useful for glaucoma therapy.

Identification and in silico analyses of novel TGFBR1 and TGFBR2 mutations in Marfan syndrome-related disorders

Very recently, heterozygous mutations in the genes encoding transforming growth factor beta receptors I (TGFBR1) and II (TGFBR2) have been reported in Loeys-Dietz aortic aneurysm syndrome (LDS). In addition, dominant TGFBR2 mutations have been identified in Marfan syndrome type 2 (MFS2) and familial thoracic aortic aneurysms and dissections (TAAD). In the past, mutations of these genes were associated with atherosclerosis and several human cancers. Here, we report a total of nine novel and one known heterozygous sequence variants in the TGFBR1 and TGFBR2 genes in nine of 70 unrelated individuals with MFS-like phenotypes who previously tested negative for mutations in the gene encoding the extracellular matrix protein fibrillin-1 (FBN1). To assess the pathogenic impact of these sequence variants, in silico analyses were performed by the PolyPhen, SIFT, and Fold-X algorithms and by means of a 3D homology model of the TGFBR2 kinase domain. Our results showed that in all but one of the patients the pathogenic effect of at least one sequence variant is highly probable (c.722C > T, c.799A > C, and c.1460G > A in TGFBR1 and c.773T > G, c.1106G > T, c.1159G > A, c.1181G > A, and c.1561T > C in TGFBR2). These deleterious alleles occurred de novo or segregated with the disease in the families, indicating a causative association between the sequence variants and clinical phenotypes. Since TGFBR2 mutations found in patients with MFS-related disorders cannot be distinguished from heterozygous TGFBR2 mutations reported in tumor samples, we emphasize the importance of segregation analysis in affected families. In order to be able to find the mutation that is indeed responsible for a MFS-related phenotype, we also propose that genetic testing for sequence alterations in TGFBR1 and TGFBR2 should be complemented by mutation screening of the FBN1 gene.

Identification of markers to characterize and sort human articular chondrocytes with enhanced in vitro chondrogenic capacity

OBJECTIVE: To identify markers associated with the chondrogenic capacity of expanded human articular chondrocytes and to use these markers for sorting of more highly chondrogenic subpopulations. METHODS: The chondrogenic capacity of chondrocyte populations derived from different donors (n = 21) or different clonal strains from the same cartilage biopsy specimen (n = 21) was defined based on the glycosaminoglycan (GAG) content of tissues generated using a pellet culture model. Selected cell populations were analyzed by microarray and flow cytometry. In some experiments, cells were sorted using antibodies against molecules found to be associated with differential chondrogenic capacity and again assessed in pellet cultures. RESULTS: Significance Analysis of Microarrays indicated that chondrocytes with low chondrogenic capacity expressed higher levels of insulin-like growth factor 1 and of catabolic genes (e.g., matrix metalloproteinase 2, aggrecanase 2), while chondrocytes with high chondrogenic capacity expressed higher levels of genes involved in cell-cell or cell-matrix interactions (e.g., CD49c, CD49f). Flow cytometry analysis showed that CD44, CD151, and CD49c were expressed at significantly higher levels in chondrocytes with higher chondrogenic capacity. Flow cytometry analysis of clonal chondrocyte strains indicated that CD44 and CD151 could also identify more chondrogenic clones. Chondrocytes sorted for brighter CD49c or CD44 signal expression produced tissues with higher levels of GAG per DNA (up to 1.4-fold) and type II collagen messenger RNA (up to 3.4-fold) than did unsorted cells. CONCLUSION: We identified markers that allow characterization of the capacity of monolayer-expanded chondrocytes to form in vitro cartilaginous tissue and enable enrichment for subpopulations with higher chondrogenic capacity. These markers might be used as a means to predict and possibly improve the outcome of cell-based cartilage repair techniques.

Observation and identification of metabolites emerging during postmortem decomposition of brain tissue by means of in situ 1H-magnetic resonance spectroscopy

Postmortem decomposition of brain tissue was investigated by (1)H-magnetic resonance spectroscopy (MRS) in a sheep head model and selected human cases. Aiming at the eventual estimation of postmortem intervals in forensic medicine, this study focuses on the characterization and identification of newly observed metabolites. In situ single-voxel (1)H-MRS at 1.5 T was complemented by multidimensional homo- and heteronuclear high-resolution NMR spectroscopy of an extract of sheep brain tissue. The inclusion of spectra of model solutions in the program LC Model confirmed the assignments in situ. The first postmortem phase was characterized mainly by changes in the concentrations of metabolites usually observed in vivo and by the appearance of previously reported decay products. About 3 days postmortem, new metabolites, including free trimethylammonium, propionate, butyrate, and iso-butyrate, started to appear in situ. Since the observed metabolites and the time course is comparable in sheep and human brain tissue, the model system seems to be appropriate.

Multi-antenna non-line-of-sight identification techniques for target localization in mobile ad-hoc networks

Target localization has a wide range of military and civilian applications in wireless mobile networks. Examples include battle-field surveillance, emergency 911 (E911), traffc alert, habitat monitoring, resource allocation, routing, and disaster mitigation. Basic localization techniques include time-of-arrival (TOA), direction-of-arrival (DOA) and received-signal strength (RSS) estimation. Techniques that are proposed based on TOA and DOA are very sensitive to the availability of Line-of-sight (LOS) which is the direct path between the transmitter and the receiver. If LOS is not available, TOA and DOA estimation errors create a large localization error. In order to reduce NLOS localization error, NLOS identifcation, mitigation, and localization techniques have been proposed. This research investigates NLOS identifcation for multiple antennas radio systems. The techniques proposed in the literature mainly use one antenna element to enable NLOS identifcation. When a single antenna is utilized, limited features of the wireless channel can be exploited to identify NLOS situations. However, in DOA-based wireless localization systems, multiple antenna elements are available. In addition, multiple antenna technology has been adopted in many widely used wireless systems such as wireless LAN 802.11n and WiMAX 802.16e which are good candidates for localization based services. In this work, the potential of spatial channel information for high performance NLOS identifcation is investigated. Considering narrowband multiple antenna wireless systems, two xvNLOS identifcation techniques are proposed. Here, the implementation of spatial correlation of channel coeffcients across antenna elements as a metric for NLOS identifcation is proposed. In order to obtain the spatial correlation, a new multi-input multi-output (MIMO) channel model based on rough surface theory is proposed. This model can be used to compute the spatial correlation between the antenna pair separated by any distance. In addition, a new NLOS identifcation technique that exploits the statistics of phase difference across two antenna elements is proposed. This technique assumes the phases received across two antenna elements are uncorrelated. This assumption is validated based on the well-known circular and elliptic scattering models. Next, it is proved that the channel Rician K-factor is a function of the phase difference variance. Exploiting Rician K-factor, techniques to identify NLOS scenarios are proposed. Considering wideband multiple antenna wireless systems which use MIMO-orthogonal frequency division multiplexing (OFDM) signaling, space-time-frequency channel correlation is exploited to attain NLOS identifcation in time-varying, frequency-selective and spaceselective radio channels. Novel NLOS identi?cation measures based on space, time and frequency channel correlation are proposed and their performances are evaluated. These measures represent a better NLOS identifcation performance compared to those that only use space, time or frequency.

Development of an optimization model for biofuel facility size and location and a simulation model for design of a biofuel supply chain

To mitigate greenhouse gas (GHG) emissions and reduce U.S. dependence on imported oil, the United States (U.S.) is pursuing several options to create biofuels from renewable woody biomass (hereafter referred to as “biomass”). Because of the distributed nature of biomass feedstock, the cost and complexity of biomass recovery operations has significant challenges that hinder increased biomass utilization for energy production. To facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization and tapping unused forest residues, it is proposed to develop biofuel supply chain models based on optimization and simulation approaches. The biofuel supply chain is structured around four components: biofuel facility locations and sizes, biomass harvesting/forwarding, transportation, and storage. A Geographic Information System (GIS) based approach is proposed as a first step for selecting potential facility locations for biofuel production from forest biomass based on a set of evaluation criteria, such as accessibility to biomass, railway/road transportation network, water body and workforce. The development of optimization and simulation models is also proposed. The results of the models will be used to determine (1) the number, location, and size of the biofuel facilities, and (2) the amounts of biomass to be transported between the harvesting areas and the biofuel facilities over a 20-year timeframe. The multi-criteria objective is to minimize the weighted sum of the delivered feedstock cost, energy consumption, and GHG emissions simultaneously. Finally, a series of sensitivity analyses will be conducted to identify the sensitivity of the decisions, such as the optimal site selected for the biofuel facility, to changes in influential parameters, such as biomass availability and transportation fuel price. Intellectual Merit The proposed research will facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization in the renewable biofuel industry. The GIS-based facility location analysis considers a series of factors which have not been considered simultaneously in previous research. Location analysis is critical to the financial success of producing biofuel. The modeling of woody biomass supply chains using both optimization and simulation, combing with the GIS-based approach as a precursor, have not been done to date. The optimization and simulation models can help to ensure the economic and environmental viability and sustainability of the entire biofuel supply chain at both the strategic design level and the operational planning level. Broader Impacts The proposed models for biorefineries can be applied to other types of manufacturing or processing operations using biomass. This is because the biomass feedstock supply chain is similar, if not the same, for biorefineries, biomass fired or co-fired power plants, or torrefaction/pelletization operations. Additionally, the research results of this research will continue to be disseminated internationally through publications in journals, such as Biomass and Bioenergy, and Renewable Energy, and presentations at conferences, such as the 2011 Industrial Engineering Research Conference. For example, part of the research work related to biofuel facility identification has been published: Zhang, Johnson and Sutherland [2011] (see Appendix A). There will also be opportunities for the Michigan Tech campus community to learn about the research through the Sustainable Future Institute.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Radiation metabolomics. 1. Identification of minimally invasive urine biomarkers for gamma-radiation exposure in mice

Gamma-radiation exposure has both short- and long-term adverse health effects. The threat of modern terrorism places human populations at risk for radiological exposures, yet current medical countermeasures to radiation exposure are limited. Here we describe metabolomics for gamma-radiation biodosimetry in a mouse model. Mice were gamma-irradiated at doses of 0, 3 and 8 Gy (2.57 Gy/min), and urine samples collected over the first 24 h after exposure were analyzed by ultra-performance liquid chromatography-time-of-flight mass spectrometry (UPLC-TOFMS). Multivariate data were analyzed by orthogonal partial least squares (OPLS). Both 3- and 8-Gy exposures yielded distinct urine metabolomic phenotypes. The top 22 ions for 3 and 8 Gy were analyzed further, including tandem mass spectrometric comparison with authentic standards, revealing that N-hexanoylglycine and beta-thymidine are urinary biomarkers of exposure to 3 and 8 Gy, 3-hydroxy-2-methylbenzoic acid 3-O-sulfate is elevated in urine of mice exposed to 3 but not 8 Gy, and taurine is elevated after 8 but not 3 Gy. Gene Expression Dynamics Inspector (GEDI) self-organizing maps showed clear dose-response relationships for subsets of the urine metabolome. This approach is useful for identifying mice exposed to gamma radiation and for developing metabolomic strategies for noninvasive radiation biodosimetry in humans.

Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.

A disease model for wheezing disorders in preschool children based on clinicians' perceptions

BACKGROUND: Wheezing disorders in childhood vary widely in clinical presentation and disease course. During the last years, several ways to classify wheezing children into different disease phenotypes have been proposed and are increasingly used for clinical guidance, but validation of these hypothetical entities is difficult. METHODOLOGY/PRINCIPAL FINDINGS: The aim of this study was to develop a testable disease model which reflects the full spectrum of wheezing illness in preschool children. We performed a qualitative study among a panel of 7 experienced clinicians from 4 European countries working in primary, secondary and tertiary paediatric care. In a series of questionnaire surveys and structured discussions, we found a general consensus that preschool wheezing disorders consist of several phenotypes, with a great heterogeneity of specific disease concepts between clinicians. Initially, 24 disease entities were described among the 7 physicians. In structured discussions, these could be narrowed down to three entities which were linked to proposed mechanisms: a) allergic wheeze, b) non-allergic wheeze due to structural airway narrowing and c) non-allergic wheeze due to increased immune response to viral infections. This disease model will serve to create an artificial dataset that allows the validation of data-driven multidimensional methods, such as cluster analysis, which have been proposed for identification of wheezing phenotypes in children. CONCLUSIONS/SIGNIFICANCE: While there appears to be wide agreement among clinicians that wheezing disorders consist of several diseases, there is less agreement regarding their number and nature. A great diversity of disease concepts exist but a unified phenotype classification reflecting underlying disease mechanisms is lacking. We propose a disease model which may help guide future research so that proposed mechanisms are measured at the right time and their role in disease heterogeneity can be studied.

Tropopause level Rossby wave breaking in the Northern Hemisphere: a feature-based validation of the ECHAM5-HAM climate model

Breaking synoptic-scale Rossby waves (RWB) at the tropopause level are central to the daily weather evolution in the extratropics and the subtropics. RWB leads to pronounced meridional transport of heat, moisture, momentum, and chemical constituents. RWB events are manifest as elongated and narrow structures in the tropopause-level potential vorticity (PV) field. A feature-based validation approach is used to assess the representation of Northern Hemisphere RWB in present-day climate simulations carried out with the ECHAM5-HAM climate model at three different resolutions (T42L19, T63L31, and T106L31) against the ERA-40 reanalysis data set. An objective identification algorithm extracts RWB events from the isentropic PV field and allows quantifying the frequency of occurrence of RWB. The biases in the frequency of RWB are then compared to biases in the time mean tropopause-level jet wind speeds. The ECHAM5-HAM model captures the location of the RWB frequency maxima in the Northern Hemisphere at all three resolutions. However, at coarse resolution (T42L19) the overall frequency of RWB, i.e. the frequency averaged over all seasons and the entire hemisphere, is underestimated by 28%.The higher-resolution simulations capture the overall frequency of RWB much better, with a minor difference between T63L31 and T106L31 (frequency errors of −3.5 and 6%, respectively). The number of large-size RWB events is significantly underestimated by the T42L19 experiment and well represented in the T106L31 simulation. On the local scale, however, significant differences to ERA-40 are found in the higher-resolution simulations. These differences are regionally confined and vary with the season. The most striking difference between T106L31 and ERA-40 is that ECHAM5-HAM overestimates the frequency of RWB in the subtropical Atlantic in all seasons except for spring. This bias maximum is accompanied by an equatorward extension of the subtropical westerlies.