Microarray results from culture-amplified mixtures of serotypes correlate poorly with microarray results from uncultured mixtures.

Empirical Bayesian models for analysing molecular serotyping microarrays

BACKGROUND Microarrays offer great potential as a platform for molecular diagnostics, testing clinical samples for the presence of numerous biomarkers in highly multiplexed assays. In this study applied to infectious diseases, data from a microarray designed for molecular serotyping of Streptococcus pneumoniae was used, identifying the presence of any one of 91 known pneumococcal serotypes from DNA extracts. This microarray incorporated oligonucleotide probes for all known capsular polysaccharide synthesis genes and required a statistical analysis of the microarray intensity data to determine which serotype, or combination of serotypes, were present within a sample based on the combination of genes detected. RESULTS We propose an empirical Bayesian model for calculating the probabilities of combinations of serotypes from the microarray data. The model takes into consideration the dependencies between serotypes, induced by genes they have in common, and by homologous genes which, although not identical, are similar to each other in sequence. For serotypes which are very similar in capsular gene composition, extra probes are included on the microarray, providing additional information which is integrated into the Bayesian model. For each serotype combination with high probability, a second model, a Bayesian random effects model is applied to determine the relative abundance of each serotype. CONCLUSIONS To assess the accuracy of the proposed analysis we applied our methods to experimental data from samples containing individual serotypes and samples containing combinations of serotypes with known levels of abundance. All but two of the known serotypes of S. pneumoniae that were tested as individual samples could be uniquely determined by the Bayesian model. The model also enabled the presence of combinations of serotypes within samples to be determined. Serotypes with very low abundance within a combination of serotypes can be detected (down to 2% abundance in this study). As well as detecting the presence of serotype combinations, an approximate measure of the percentage abundance of the serotypes within the combination can be obtained.

Frozen robust multiarray analysis (fRMA).

Robust multiarray analysis (RMA) is the most widely used preprocessing algorithm for Affymetrix and Nimblegen gene expression microarrays. RMA performs background correction, normalization, and summarization in a modular way. The last 2 steps require multiple arrays to be analyzed simultaneously. The ability to borrow information across samples provides RMA various advantages. For example, the summarization step fits a parametric model that accounts for probe effects, assumed to be fixed across arrays, and improves outlier detection. Residuals, obtained from the fitted model, permit the creation of useful quality metrics. However, the dependence on multiple arrays has 2 drawbacks: (1) RMA cannot be used in clinical settings where samples must be processed individually or in small batches and (2) data sets preprocessed separately are not comparable. We propose a preprocessing algorithm, frozen RMA (fRMA), which allows one to analyze microarrays individually or in small batches and then combine the data for analysis. This is accomplished by utilizing information from the large publicly available microarray databases. In particular, estimates of probe-specific effects and variances are precomputed and frozen. Then, with new data sets, these are used in concert with information from the new arrays to normalize and summarize the data. We find that fRMA is comparable to RMA when the data are analyzed as a single batch and outperforms RMA when analyzing multiple batches. The methods described here are implemented in the R package fRMA and are currently available for download from the software section of http://rafalab.jhsph.edu.

Data analysis methods for detection of differential protein expression in two-dimensional gel electrophoresis.

The recent development of microarray technology has led statisticians and bioinformaticians to develop new statistical methodologies for comparing different biological samples. The objective is to identify a small number of differentially expressed genes from among thousands. In quantitative proteomics, analysis of protein expression using two-dimensional gel electrophoresis shows some similarities with transcriptomic studies. Thus, the goal of this study was to evaluate different data analysis methodologies widely used in array analysis using different proteomic data sets of hundreds of proteins. Even with few replications, the significance analysis of microarrays method appeared to be more powerful than the Student's t test in truly declaring differentially expressed proteins. This procedure will avoid wasting time due to false positives and losing information with false negatives.

The PneuCarriage Project: A Multi-Centre Comparative Study to Identify the Best Serotyping Methods for Examining Pneumococcal Carriage in Vaccine Evaluation Studies

BACKGROUND The pneumococcus is a diverse pathogen whose primary niche is the nasopharynx. Over 90 different serotypes exist, and nasopharyngeal carriage of multiple serotypes is common. Understanding pneumococcal carriage is essential for evaluating the impact of pneumococcal vaccines. Traditional serotyping methods are cumbersome and insufficient for detecting multiple serotype carriage, and there are few data comparing the new methods that have been developed over the past decade. We established the PneuCarriage project, a large, international multi-centre study dedicated to the identification of the best pneumococcal serotyping methods for carriage studies. METHODS AND FINDINGS Reference sample sets were distributed to 15 research groups for blinded testing. Twenty pneumococcal serotyping methods were used to test 81 laboratory-prepared (spiked) samples. The five top-performing methods were used to test 260 nasopharyngeal (field) samples collected from children in six high-burden countries. Sensitivity and positive predictive value (PPV) were determined for the test methods and the reference method (traditional serotyping of >100 colonies from each sample). For the alternate serotyping methods, the overall sensitivity ranged from 1% to 99% (reference method 98%), and PPV from 8% to 100% (reference method 100%), when testing the spiked samples. Fifteen methods had ≥70% sensitivity to detect the dominant (major) serotype, whilst only eight methods had ≥70% sensitivity to detect minor serotypes. For the field samples, the overall sensitivity ranged from 74.2% to 95.8% (reference method 93.8%), and PPV from 82.2% to 96.4% (reference method 99.6%). The microarray had the highest sensitivity (95.8%) and high PPV (93.7%). The major limitation of this study is that not all of the available alternative serotyping methods were included. CONCLUSIONS Most methods were able to detect the dominant serotype in a sample, but many performed poorly in detecting the minor serotype populations. Microarray with a culture amplification step was the top-performing method. Results from this comprehensive evaluation will inform future vaccine evaluation and impact studies, particularly in low-income settings, where pneumococcal disease burden remains high.

Microarrays for the study of viral gene expression during human cytomegalovirus latent infection.

Human cytomegalovirus (HCMV) is one of the largest known DNA viruses. It is ubiquitous, and following resolution of primary productive infection, it persists in the human host by establishing a lifelong latent infection in myeloid lineage cells such as monocytes and their progenitors. Most adults with HCMV infection are healthy but it can cause neurologic deficits in infants, and remains an important cause of morbidity and mortality in the immunosuppressed patient. Microarray-based studies of HCMV have provided useful information about genes that are transcriptionally active during both productive and latent phases of infection. This chapter describes how to study genes in HCMV using microarrays and two cell types (productively infected human foreskin fibroblasts, and latently infected primary human myeloid progenitor cells).

Development of the Human Infant Intestinal Microbiota

Almost immediately after a human being is born, so too is a new microbial ecosystem, one that resides in that person's gastrointestinal tract. Although it is a universal and integral part of human biology, the temporal progression of this process, the sources of the microbes that make up the ecosystem, how and why it varies from one infant to another, and how the composition of this ecosystem influences human physiology, development, and disease are still poorly understood. As a step toward systematically investigating these questions, we designed a microarray to detect and quantitate the small subunit ribosomal RNA (SSU rRNA) gene sequences of most currently recognized species and taxonomic groups of bacteria. We used this microarray, along with sequencing of cloned libraries of PCR-amplified SSU rDNA, to profile the microbial communities in an average of 26 stool samples each from 14 healthy, full-term human infants, including a pair of dizygotic twins, beginning with the first stool after birth and continuing at defined intervals throughout the first year of life. To investigate possible origins of the infant microbiota, we also profiled vaginal and milk samples from most of the mothers, and stool samples from all of the mothers, most of the fathers, and two siblings. The composition and temporal patterns of the microbial communities varied widely from baby to baby. Despite considerable temporal variation, the distinct features of each baby's microbial community were recognizable for intervals of weeks to months. The strikingly parallel temporal patterns of the twins suggested that incidental environmental exposures play a major role in determining the distinctive characteristics of the microbial community in each baby. By the end of the first year of life, the idiosyncratic microbial ecosystems in each baby, although still distinct, had converged toward a profile characteristic of the adult gastrointestinal tract.

SEROSURVEY FOR FELINE LEUKEMIA VIRUS AND LENTIVIRUSES IN CAPTIVE SMALL NEOTROPIC FELIDS IN SÃO PAULO STATE, BRAZIL

Abstract Feline leukemia virus (FeLV), Gammaretrovirus, and feline immunodeficiency virus, a Lentivirus, are members of the family Retroviridae, and may establish persistent infections in the domestic cat (Felis catus). Cytoproliferative and cytosuppressive disorders may result from infection with these viruses. Morbidity and mortality rates are high in domestic cats worldwide. Infection of endangered neotropic small felids with these viruses could be devastating. To investigate the prevalence of FeLV and feline lentiviruses in neotropic small felids kept in captivity in São Paulo state, Brazil, serum samples from 104 animals belonging to the species Leopardus pardalis, Leopardus tigrinus, Leopardus wiedii, Herpailurus yaguarondi, and Oncifelis geoffroyi were tested for FeLV and feline lentiviruses by commercially available immunoassays. All results were negative, suggesting that retrovirus infection is not an important clinical problem in these populations. Because domestic cats in São Paulo city are naturally infected with these pathogens, and feral cats are commonly found in zoologic facilities in Brazil, preventive measures should be taken to avoid transmission of retroviruses to naive populations of wild and captive neotropic felids in Brazil.

Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise

A major goal of biology is to provide a quantitative description of cellular behaviour. This task, however, has been hampered by the difficulty in measuring protein abundances and their variation. Here we present a strategy that pairs high-throughput flow cytometry and a library of GFP-tagged yeast strains to monitor rapidly and precisely protein levels at single-cell resolution. Bulk protein abundance measurements of >2,500 proteins in rich and minimal media provide a detailed view of the cellular response to these conditions, and capture many changes not observed by DNA microarray analyses. Our single-cell data argue that noise in protein expression is dominated by the stochastic production/destruction of messenger RNAs. Beyond this global trend, there are dramatic protein-specific differences in noise that are strongly correlated with a protein's mode of transcription and its function. For example, proteins that respond to environmental changes are noisy whereas those involved in protein synthesis are quiet. Thus, these studies reveal a remarkable structure to biological noise and suggest that protein noise levels have been selected to reflect the costs and potential benefits of this variation.

Inter-viral conflicts that exploit host CRISPR immune systems of Sulfolobus.

Infection of Sulfolobus islandicus REY15A with mixtures of different Sulfolobus viruses, including STSV2, did not induce spacer acquisition by the host CRISPR immune system. However, coinfection with the tailed fusiform viruses SMV1 and STSV2 generated hyperactive spacer acquisition in both CRISPR loci, exclusively from STSV2, with the resultant loss of STSV2 but not SMV1. SMV1 was shown to activate adaptation while itself being resistant to CRISPR-mediated adaptation and DNA interference. Exceptionally, a single clone S-1 isolated from an SMV1 + STSV2-infected culture, that carried STSV2-specific spacers and had lost STSV2 but not SMV1, acquired spacers from SMV1. This effect was also reproducible on reinfecting wild-type host cells with a variant SMV1 isolated from the S-1 culture. The SMV1 variant lacked a virion protein ORF114 that was shown to bind DNA. This study also provided evidence for: (i) limits on the maximum sizes of CRISPR loci; (ii) spacer uptake strongly retarding growth of infected cultures; (iii) protospacer selection being essentially random and non-directional, and (iv) the reversible uptake of spacers from STSV2 and SMV1. A hypothesis is presented to explain the interactive conflicts between SMV1 and the host CRISPR immune system.

Monte Carlo tests for associations between disease and alleles at highly polymorphic loci.

In an association analysis comparing cases and controls with respect to allele frequencies at a highly polymorphic locus, a potential problem is that the conventional chi-squared test may not be valid for a large, sparse contingency table. However, reliance on statistics with known asymptotic distribution is now unnecessary, as Monte Carlo simulations can be performed to estimate the significance level of any test statistic. We have implemented a Monte Carlo method for four 'chi-squared' test statistics, three of which involved combination of alleles, and evaluated their performance on a real data set. Combining rare alleles to avoid small expected cell counts, and considering each allele in turn against the rest, reduced the power to detect a genuine association when the number of alleles was very large. We should either not combine alleles at all, or combine them in such a way that preserves the evidence for an association.

Chromatin immunoprecipitation and microarray-based analysis of protein location

Genome-wide location analysis, also known as ChIP-Chip, combines chromatin immunoprecipitation and DNA microarray analysis to identify protein-DNA interactions that occur in living cells. Protein-DNA interactions are captured in vivo by chemical crosslinking. Cell lysis, DNA fragmentation and immunoaffinity purification of the desired protein will co-purify DNA fragments that are associated with that protein. The enriched DNA population is then labeled, combined with a differentially labeled reference sample and applied to DNA microarrays to detect enriched signals. Various computational and bioinformatic approaches are then applied to normalize the enriched and reference channels, to connect signals to the portions of the genome that are represented on the DNA microarrays, to provide confidence metrics and to generate maps of protein-genome occupancy. Here, we describe the experimental protocols that we use from crosslinking of cells to hybridization of labeled material, together with insights into the aspects of these protocols that influence the results. These protocols require approximately 1 week to complete once sufficient numbers of cells have been obtained, and have been used to produce robust, high-quality ChIP-chip results in many different cell and tissue types.

Population density profiles of nasopharyngeal carriage of 5 bacterial species in pre-school children measured using quantitative PCR offer potential insights into the dynamics of transmission.

Bacterial vaccines can reduce carriage rates. Colonization is usually a binary endpoint. Real time quantitative PCR (qPCR) can quantify bacterial DNA in mucosal samples over a wide range. Using culture and single-gene species-specific qPCRs for Streptococcus pneumoniae (lytA), Streptococcus pyogenes (ntpC), Moraxella catarrhalis (ompJ), Haemophilus influenzae (hdp) and Staphylococcus aureus (nuc) and standard curves against log-phase reference strain broth cultures we described frequency and peak density distributions of carriage in nasopharyngeal swabs from 161 healthy 2-4 y old children collected into STGG broth. In general, detection by qPCR and culture was consistent. Discordance mostly occurred at lower detection thresholds of both methods, although PCR assays for S. pyogenes and S. aureus were less sensitive. Density varied across 5-7 orders of magnitude for the 5 species with the abundant species skewed toward high values (modes: S. pneumoniae log3-4, M. catarrhalis & H. influenzae log4-5 CFU/ml broth). Wide ranges of bacterial DNA concentrations in healthy children carrying these bacteria could mean that different individuals at different times vary greatly in infectiousness. Understanding the host, microbial and environmental determinants of colonization density will permit more accurate prediction of vaccine effectiveness.

Macrophages induce COX-2 expression in breast cancer cells: role of IL-1β autoamplification.

Tumor-associated macrophages and high levels of cyclooxygenase-2 (COX-2) are associated with poor prognosis in breast cancer patients, but their potential interdependence has not been evaluated. The objective of this study was to determine whether macrophages regulate COX-2 expression in breast cancer cells. For this purpose, THP-1 cells were cocultured with HCC1954 breast cancer cells. Coculture led to increased COX-2 expression in the HCC1954 cells and elevated prostaglandin E(2) levels in conditioned media. Similar results were observed when THP-1 cells were incubated with HCC1937 breast cancer cells or when human monocyte-derived macrophages were cocultured with HCC1954 cells. Coculture triggered production of reactive oxygen species (ROS) in HCC1954 cells. COX-2 induction was blocked in cells preincubated with an reduced nicotinamide adenine dinucleotide phosphate (NADPH) oxidase inhibitor or by silencing p67PHOX, a subunit of NADPH oxidase. ROS production triggered activation of Src and mitogen-activated protein kinases (MAPKs). Blocking Src or MAPK activities or antagonizing the activator protein-1 (AP-1) transcription factor attenuated COX-2 induction in HCC1954 cells. Coculture caused rapid induction of interleukin-1β (IL-1β) in both breast cancer cells and macrophages. Increased IL-1β expression was blocked by an interleukin-1 receptor antagonist (IL-1Ra), suggesting autocrine and paracrine effects. Importantly, macrophage-induced COX-2 expression was blocked in HCC1954 cells preincubated with IL-1Ra or anti-IL-1β IgG. Together, these results indicate that macrophage-mediated induction of COX-2 in breast cancer cells is a consequence of IL-1β-mediated stimulation of ROS→Src→MAPK→AP-1 signaling. IL-1β-dependent induction of COX-2 in breast cancer cells provides a mechanism whereby macrophages contribute to tumor progression and potential therapeutic targets in breast cancer.

Comparison of agar dilution, microdilution, and disk elution methods for measuring the synergy of cefotaxime and its metabolite against anaerobes.

The activities of cefotaxime (CTX) and desacetyl cefotaxime (des-CTX) were tested both singly and in combination against 173 anaerobic clinical isolates. The MIC of CTX for 50% of 60 Bacteroides fragilis isolates was 22.4 micrograms/ml in broth, compared with 47.4 micrograms/ml in agar. This reduced efficacy in agar was seen with all species tested and is in apparent conflict with reported clinical efficacy of the drug. Synergy between CTX and des-CTX was observed with 70 to 100% of the isolates, including 60% of all Bacteroides spp. tested. The susceptibility results in a synergy system correlated well with those noted in a broth-disk elution method incorporating 32 micrograms of CTX and 8 micrograms of des-CTX per ml. The correlation was poorer when the broth-disk method contained 16 micrograms of CTX and 8 micrograms of des-CTX per ml.

Differential distribution of IL28B.rs12979860 single-nucleotide polymorphism among Egyptian healthcare workers with and without a hepatitis C virus-specific cellular immune response

The CC genotype of the interleukin (IL)-28B.rs12979860 gene has been associated with spontaneous hepatitis C virus (HCV) clearance and treatment response. The distribution and correlation of an IL28B.rs12979860 single-nucleotide polymorphism (SNP) with HCV-specific cell-mediated immune (CMI) responses among Egyptian healthcare workers (HCWs) is not known. We determined this relationship in 402 HCWs who serve a patient cohort with ~85 % HCV prevalence. We enrolled 402 HCWs in four groups: group 1 (n = 258), seronegative aviremic subjects; group 2 (n = 25), seronegative viremic subjects; group 3 (n = 41), subjects with spontaneously resolved HCV infection; and group 4 (n = 78), chronic HCV patients. All subjects were tested for an HCV-specific CMI response using an ex-vivo interferon-gamma (IFNγ) ELISpot assay with nine HCV genotype-4a overlapping 15-mer peptide pools corresponding to all of the HCV proteins. All subjects were tested for IL28B.rs12979860 SNP by real-time PCR. An HCV-specific CMI was demonstrated in ~27 % of the seronegative aviremic HCWs (group 1), suggesting clearance of infection after low-level exposure to HCV. The frequency of IL28B.rs12979860 C allele homozygosity in the four groups was 49 %, 48 %, 49 %, and 23 %, while that of the T allele was 14 %, 16 %, 12 and 19 %, respectively, suggesting differential distributions among subjects with different HCV status. As reported, IL28B.rs12979860 predicted the outcome of HCV infection (p < 0.05), but we did not find any relationship between the IL28B genotypes and the outcome of HCV-specific CMI responses in the four groups (p > 0.05). The data show differential IL28B.rs12979860 genotype distribution among Egyptian HCWs with different HCV status and could not predict the outcome of HCV-specific CMI responses.

Distinct and essential roles of transcription factors IRF-3 and IRF-7 in response to viruses for IFN-alpha/beta gene induction.

Induction of the interferon (IFN)-alpha/beta gene transcription in virus-infected cells is an event central to innate immunity. Mice lacking the transcription factor IRF-3 are more vulnerable to virus infection. In embryonic fibroblasts, virus-induced IFN-alpha/beta gene expression levels are reduced and the spectrum of the IFN-alpha mRNA subspecies altered. Furthermore, cells additionally defective in IRF-7 expression totally fail to induce these genes in response to infections by any of the virus types tested. In these cells, a normal profile of IFN-alpha/beta mRNA induction can be achieved by coexpressing both IRF-3 and IRF-7. These results demonstrate the essential and distinct roles of thetwo factors, which together ensure the transcriptional efficiency and diversity of IFN-alpha/beta genes for the antiviral response.

RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome

BACKGROUND RNA-Seq is revolutionizing the way transcript abundances are measured. A key challenge in transcript quantification from RNA-Seq data is the handling of reads that map to multiple genes or isoforms. This issue is particularly important for quantification with de novo transcriptome assemblies in the absence of sequenced genomes, as it is difficult to determine which transcripts are isoforms of the same gene. A second significant issue is the design of RNA-Seq experiments, in terms of the number of reads, read length, and whether reads come from one or both ends of cDNA fragments. RESULTS We present RSEM, an user-friendly software package for quantifying gene and isoform abundances from single-end or paired-end RNA-Seq data. RSEM outputs abundance estimates, 95% credibility intervals, and visualization files and can also simulate RNA-Seq data. In contrast to other existing tools, the software does not require a reference genome. Thus, in combination with a de novo transcriptome assembler, RSEM enables accurate transcript quantification for species without sequenced genomes. On simulated and real data sets, RSEM has superior or comparable performance to quantification methods that rely on a reference genome. Taking advantage of RSEM's ability to effectively use ambiguously-mapping reads, we show that accurate gene-level abundance estimates are best obtained with large numbers of short single-end reads. On the other hand, estimates of the relative frequencies of isoforms within single genes may be improved through the use of paired-end reads, depending on the number of possible splice forms for each gene. CONCLUSIONS RSEM is an accurate and user-friendly software tool for quantifying transcript abundances from RNA-Seq data. As it does not rely on the existence of a reference genome, it is particularly useful for quantification with de novo transcriptome assemblies. In addition, RSEM has enabled valuable guidance for cost-efficient design of quantification experiments with RNA-Seq, which is currently relatively expensive.

Intrinsic genetic characteristics determine tumor-modifying capacity of fibroblasts: matrix metalloproteinase-3 5A/5A genotype enhances breast cancer cell invasion

Stromal fibroblasts can contribute to tumor invasion through the release of matrix metalloproteinases (MMPs). Population studies have suggested that single nucleotide polymorphisms (SNPs) in MMP genes influence levels of expression and may be associated with breast cancer risk and with disease progression. This study directly examined the impact of MMP SNP genotype on the ability of host fibroblasts to promote tumor cell invasion. Primary breast fibroblasts were isolated from patients with (n = 13) or without (n = 19) breast cancer, and their ability to promote breast cancer cell invasion was measured in in vitro invasion assays. Fibroblast invasion-promoting capacity (IPC) was analyzed in relation to donor type (tumor or non-tumor patient), MMP-1, MMP-3, and MMP-9 SNP genotype and MMP activity using independent samples t test and analysis of variance. All statistical tests were two-sided. Tumor-derived fibroblasts promoted higher levels of invasion than normal fibroblasts (p = 0.041). When IPC was related to genotype, higher levels of IPC were generated by tumor fibroblasts with the high-expressing MMP-3 5A/5A genotype compared with the 5A/6A and 6A/6A genotypes (p = 0.05 and 0.07, respectively), and this was associated with enhanced MMP-3 release. The functional importance of MMP-3 was demonstrated by enhanced invasion in the presence of recombinant MMP-3, whereas reduction occurred in the presence of a specific MMP-3 inhibitor. An inverse relationship was demonstrated between fibroblast IPC and the high-expressing MMP-1 genotype (p = 0.031), but no relationship was seen with MMP-9 SNP status. In contrast, normal fibroblasts showed no variation in IPC in relation to MMP genotype, with MMP-3 5A/5A fibroblasts exhibiting significantly lower levels of IPC than their tumor-derived counterparts (p = 0.04). This study has shown that tumor-derived fibroblasts exhibit higher levels of IPC than normal fibroblasts and that the MMP-3 5A/5A genotype contributes to this through enhanced MMP-3 release. Despite a high-expressing genotype, normal fibroblasts do not exhibit higher IPC or enhanced MMP release. This suggests that more complex changes occur in tumor-derived fibroblasts, enabling full expression of the MMP SNP genotype and these possibly are epigenetic in nature. The results do suggest that, in women with breast cancer, a high-expressing MMP-3 genotype may promote tumor progression more effectively.

Androgen regulation of micro-RNAs in prostate cancer.

BACKGROUND Androgens play a critical role in the growth of both androgen dependent and castration-resistant prostate cancer (CRPC). Only a few micro-RNAs (miRNAs) have been suggested to be androgen regulated. We aim to identify androgen regulated miRNAs. METHODS We utilized LNCaP derived model, we have established, and which overexpresses the androgen receptor (AR), the VCaP cell line, and 13 intact-castrated prostate cancer (PC) xenograft pairs, as well as clinical specimens of untreated (PC) and CRPC. The expression of miRNAs was analyzed by microarrays and quantitative RT-PCR (Q-RT-PCR). Transfection of pre-miR-141 and anti-miR-141 was also used. RESULTS Seventeen miRNAs were > 1.5-fold up- or downregulated upon dihydrotestosterone (DHT) treatment in the cell lines, and 42 after castration in the AR-positive xenografts. Only four miRNAs (miR-10a, miR-141, miR-150*, and miR-1225-5p) showed similar androgen regulation in both cell lines and xenografts. Of those, miR-141 was found to be expressed more in PC and CRPC compared to benign prostate hyperplasia. Additionally, the overexpression of miR-141 enhanced growth of parental LNCaP cells while inhibition of miR-141 by anti-miR-141 suppressed the growth of the LNCaP subline overexpressing AR. CONCLUSIONS Only a few miRNAs were found to be androgen-regulated in both cell lines and xenografts models. Of those, the expression of miR-141 was upregulated in cancer. The ectopic overexpression of miR-141 increased growth of LNCaP cell suggesting it may contribute to the progression of PC.

Three common functional polymorphisms in microRNA encoding genes in the susceptibility to hepatocellular carcinoma: a systematic review and meta-analysis.

Emerging evidences have shown that common genetic polymorphisms in microRNAs may be associated with the development of hepatocellular carcinoma (HCC); but individually published studies and previous meta-analyses revealed inconclusive results. The aims of this review and meta-analysis are to assess whether common single-nucleotide polymorphisms (SNPs) in the genes encoding the microRNAs are associated with susceptibility to HCC development and clinicopathologic characteristics of hepatitis B virus (HBV) related HCC. A computerized search was performed in PubMed, Embase, Web of Science and China BioMedicine (CBM) databases to identify relevant articles published before January 1st 2013. Ten case-control studies were assessed with a total of 3437 cases and 3437 healthy controls. Three common functional SNPs in miRNA-encoding genes were found, including miR-146a G>C (rs2910164), miR-196a-2 C>T (rs11614913) and miR-499 T>C (rs3746444). This meta-analysis revealed that the miR-146a C variant was associated with a decrease in HCC risk, especially among Asian and male populations; while the miR-196a-2 T variant was associated with susceptibility to HCC among Caucasian populations. However, we failed to find any significant correlations between the miR-499 C polymorphism and HCC risks. When further stratification on HBV status was conducted, a similar trend of association between the three SNPs and the HBV-related HCC risks was observed, but these results were not statistically significant due to small sample sizes. The current meta-analysis demonstrates that SNPs contained in the genes encoding miR-146a and miR-196a-2 may play a major role in genetic susceptibility to HCC.

Matrix eQTL: Ultra fast eQTL analysis via large matrix operations

Expression quantitative trait loci (eQTL) mapping aims to determine genomic regions that regulate gene transcription. Expression QTL is used to study the regulatory structure of normal tissues and to search for genetic factors in complex diseases such as cancer, diabetes, and cystic fibrosis. A modern eQTL dataset contains millions of SNPs and thousands of transcripts measured for hundreds of samples. This makes the analysis computationally complex as it involves independent testing for association for every transcript-SNP pair. The heavy computational burden makes eQTL analysis less popular, often forces analysts to restrict their attention to just a subset of transcripts and SNPs. As larger genotype and gene expression datasets become available, the demand for fast tools for eQTL analysis increases. We present a new method for fast eQTL analysis via linear models, called Matrix eQTL. Matrix eQTL can model and test for association using both linear regression and ANOVA models. The models can include covariates to account for such factors as population structure, gender, and clinical variables. It also supports testing of heteroscedastic models and models with correlated errors. In our experiment on large datasets Matrix eQTL was thousands of times faster than the existing popular software for QTL/eQTL analysis. Matrix eQTL is implemented as both Matlab and R packages and thus can easily be run on Windows, Mac OS, and Linux systems. The software is freely available at the following address: http://www.bios.unc.edu/research/genomic_software/Matrix_eQTL

Corresponding author:

Over the years, methods of cytogenetic analysis evolved and became part of routine laboratory testing, providing valuable diagnostic and prognostic information in hematologic disorders. Karyotypic aberrations contribute to the understanding of the molecular pathogenesis of disease and thereby to rational application of therapeutic modalities. Most of the progress in this field stems from the application of metaphase cytogenetics (MC), but recently, novel molecular technologies have been introduced that complement MC and overcome many of the limitations of traditional cytogenetics, including a need for cell culture. Whole genome scanning using comparative genomic hybridization and single nucleotide polymorphism arrays (CGH-A; SNP-A) can be used for analysis of somatic or clonal unbalanced chromosomal defects. In SNP-A, the combination of copy number detection and genotyping enables diagnosis of copy-neutral loss of heterozygosity, a lesion that cannot be detected using MC but may have important pathogenetic implications. Overall, whole genome scanning arrays, despite the drawback of an inability to detect balanced translocations, allow for discovery of chromosomal defects in a higher proportion of patients with hematologic malignancies. Newly detected chromosomal aberrations, including somatic uniparental disomy, may lead to more precise prognostic schemes in many diseases.

Shaping the landscape of the Escherichia coli chromosome: replication-transcription encounters in cells with an ectopic replication origin

Each cell division requires the unwinding of millions of DNA base pairs to allow chromosome duplication and gene transcription. As DNA replication and transcription share the same template, conflicts between both processes are unavoidable and head-on collisions are thought to be particularly problematic. Surprisingly, a recent study reported unperturbed cell cycle progression in Escherichia coli cells with an ectopic replication origin in which highly transcribed rrn operons were forced to be replicated opposite to normal. In this study we have re-generated a similar strain and found the doubling time to be twice that of normal cells. Replication profiles of this background revealed significant deviations in comparison to wild-type profiles, particularly in highly transcribed regions and the termination area. These deviations were alleviated by mutations that either inactivate the termination area or destabilise RNA polymerase complexes and allow their easier displacement by replication forks. Our data demonstrate that head-on replication-transcription conflicts are highly problematic. Indeed, analysis of the replication profile of the previously published E. coli construct revealed a chromosomal rearrangement that alleviates replication-transcription conflicts in an intriguingly simple way. Our data support the idea that avoiding head-on collisions has significantly contributed to shaping the distinct architecture of bacterial chromosomes.

Circulating microRNAs as biomarkers for hepatocellular carcinoma.

GOALS We investigated whether measurement of serum levels of the microRNAs (miRNAs) miR-16, miR-195, and miR-199a, alone or in combination with conventional serum markers, can help to differentiate hepatocellular carcinoma (HCC) from chronic liver diseases (CLDs). BACKGROUND Recent reports suggest a link between aberrant expression of miRNA, and HCC. STUDY This retrospective analysis was conducted using sera from 105 HCC patients, 107 CLD patients, and 71 normal control subjects. The miRNAs were measured using real-time reverse transcription-polymerase chain reaction. The conventional HCC markers α-fetoprotein (AFP), lens culinaris agglutinin-reactive AFP (AFP-L3%), and des-γ-carboxyprothrombin (DCP) were measured with commercial kits. RESULTS Serum levels of miR-16 and miR-199a were significantly lower in HCC than in CLD patients or control subjects (P<0.01). As a single marker, miR-16 had the highest sensitivity for HCC, followed by miR-199a, AFP, DCP, AFP-L3%, and miR-195. The combination of miR-16, AFP, AFP-L3%, and DCP yielded the optimal combination of sensitivity (92.4%) and specificity (78.5%) for HCC, overall and when analysis was restricted to patients with tumors size smaller than 3 cm. As a second-line HCC marker, miR-16 yielded positive HCC predictions in 18 of the 26 (69.2%) HCC patients with negative results on all 3 conventional markers, most of whom had tumors size smaller than 3 cm; miR-16 was falsely positive in only 12 of 96 (12.5%) CLD patients. CONCLUSIONS The addition of miR-16 to conventional serum markers improved sensitivity and specificity for HCC. Use of miR-16 for second-line testing in cases considered negative on the basis of conventional HCC markers should be explored in larger, prospective studies.

Rare Variants Create Synthetic Genome-Wide Associations

Genome-wide association studies (GWAS) have now identified at least 2,000 common variants that appear associated with common diseases or related traits (http://www.genome.gov/gwastudies), hundreds of which have been convincingly replicated. It is generally thought that the associated markers reflect the effect of a nearby common (minor allele frequency >0.05) causal site, which is associated with the marker, leading to extensive resequencing efforts to find causal sites. We propose as an alternative explanation that variants much less common than the associated one may create "synthetic associations" by occurring, stochastically, more often in association with one of the alleles at the common site versus the other allele. Although synthetic associations are an obvious theoretical possibility, they have never been systematically explored as a possible explanation for GWAS findings. Here, we use simple computer simulations to show the conditions under which such synthetic associations will arise and how they may be recognized. We show that they are not only possible, but inevitable, and that under simple but reasonable genetic models, they are likely to account for or contribute to many of the recently identified signals reported in genome-wide association studies. We also illustrate the behavior of synthetic associations in real datasets by showing that rare causal mutations responsible for both hearing loss and sickle cell anemia create genome-wide significant synthetic associations, in the latter case extending over a 2.5-Mb interval encompassing scores of "blocks" of associated variants. In conclusion, uncommon or rare genetic variants can easily create synthetic associations that are credited to common variants, and this possibility requires careful consideration in the interpretation and follow up of GWAS signals.

The developmental potential of iPSCs is greatly influenced by reprogramming factor selection.

Induced pluripotent stem cells (iPSCs) are commonly generated by transduction of Oct4, Sox2, Klf4, and Myc (OSKM) into cells. Although iPSCs are pluripotent, they frequently exhibit high variation in terms of quality, as measured in mice by chimera contribution and tetraploid complementation. Reliably high-quality iPSCs will be needed for future therapeutic applications. Here, we show that one major determinant of iPSC quality is the combination of reprogramming factors used. Based on tetraploid complementation, we found that ectopic expression of Sall4, Nanog, Esrrb, and Lin28 (SNEL) in mouse embryonic fibroblasts (MEFs) generated high-quality iPSCs more efficiently than other combinations of factors including OSKM. Although differentially methylated regions, transcript number of master regulators, establishment of specific superenhancers, and global aneuploidy were comparable between high- and low-quality lines, aberrant gene expression, trisomy of chromosome 8, and abnormal H2A.X deposition were distinguishing features that could potentially also be applicable to human.

Cytokinesis failure generating tetraploids promotes tumorigenesis in p53-null cells

A long-standing hypothesis on tumorigenesis is that cell division failure, generating genetically unstable tetraploid cells, facilitates the development of aneuploid malignancies. Here we test this idea by transiently blocking cytokinesis in p53-null (p53-/-) mouse mammary epithelial cells (MMECs), enabling the isolation of diploid and tetraploid cultures. The tetraploid cells had an increase in the frequency of whole-chromosome mis-segregation and chromosomal rearrangements. Only the tetraploid cells were transformed in vitro after exposure to a carcinogen. Furthermore, in the absence of carcinogen, only the tetraploid cells gave rise to malignant mammary epithelial cancers when transplanted subcutaneously into nude mice. These tumours all contained numerous non-reciprocal translocations and an 8–30-fold amplification of a chromosomal region containing a cluster of matrix metalloproteinase (MMP) genes. MMP overexpression is linked to mammary tumours in humans and animal models. Thus, tetraploidy enhances the frequency of chromosomal alterations and promotes tumour development in p53-/- MMECs.

Probable person to person transmission of novel avian influenza A (H7N9) virus in Eastern China, 2013: epidemiological investigation

OBJECTIVE To determine whether the novel avian influenza H7N9 virus can transmit from person to person and its efficiency. DESIGN Epidemiological investigations conducted after a family cluster of two patients with avian H7N9 in March 2013. SETTING Wuxi, Eastern China. PARTICIPANTS Two patients, their close contacts, and relevant environments. Samples from the patients and environments were collected and tested by real time reverse transcriptase-polymerase chain reaction (rRT-PCR), viral culture, and haemagglutination inhibition assay. Any contacts who became ill had samples tested for avian H7N9 by rRT-PCR. Paired serum samples were obtained from contacts for serological testing by haemagglutination inhibition assays. MAIN OUTCOMES MEASURES Clinical data, history of exposure before the onset of illnesses, and results of laboratory testing of pathogens and further analysis of sequences and phylogenetic tree to isolated strains. RESULTS The index patient became ill five to six days after his last exposure to poultry. The second patient, his daughter aged 32, who provided unprotected bedside care in the hospital, had no known exposure to poultry. She developed symptoms six days after her last contact with her father. Two strains were isolated successfully from the two patients. Genome sequence and analyses of phylogenetic trees showed that both viruses were almost genetically identical. Forty three close contacts of both patients were identified. One had mild illness but had negative results for avian H7N9 by rRT-PCR. All 43 close contacts tested negative for haemagglutination inhibition antibodies specific for avian H7N9. CONCLUSIONS The infection of the daughter probably resulted from contact with her father (the index patient) during unprotected exposure, suggesting that in this cluster the virus was able to transmit from person to person. The transmissibility was limited and non-sustainable.

Genome-wide association of IL28B with response to pegylated interferon-α and ribavirin therapy for chronic hepatitis C

The recommended treatment for patients with chronic hepatitis C, pegylated interferon-α (PEG-IFN-α) plus ribavirin (RBV), does not provide sustained virologic response (SVR) in all patients. We report a genome-wide association study (GWAS) to null virological response (NVR) in the treatment of patients with hepatitis C virus (HCV) genotype 1 within a Japanese population. We found two SNPs near the gene IL28B on chromosome 19 to be strongly associated with NVR (rs12980275, P = 1.93 × 10−13, and rs8099917, 3.11 × 10−15). We replicated these associations in an independent cohort (combined P values, 2.84 × 10−27 (OR = 17.7; 95% CI = 10.0–31.3) and 2.68 × 10−32 (OR = 27.1; 95% CI = 14.6–50.3), respectively). Compared to NVR, these SNPs were also associated with SVR (rs12980275, P = 3.99 × 10−24, and rs8099917, P = 1.11 × 10−27). In further fine mapping of the region, seven SNPs (rs8105790, rs11881222, rs8103142, rs28416813, rs4803219, rs8099917 and rs7248668) located in the IL28B region showed the most significant associations (P = 5.52 × 10−28–2.68 × 10−32; OR = 22.3–27.1). Real-time quantitative PCR assays in peripheral blood mononuclear cells showed lower IL28B expression levels in individuals carrying the minor alleles (P = 0.015).

Osteoclast nuclei of myeloma patients show chromosome translocations specific for the myeloma cell clone: a new type of cancer-host partnership?

A major clinical manifestation of bone cancers is bone destruction. It is widely accepted that this destruction is not caused by the malignant cells themselves, but by osteoclasts, multinucleated cells of monocytic origin that are considered to be the only cells able to degrade bone. The present study demonstrates that bone-resorbing osteoclasts from myeloma patients contain nuclei with translocated chromosomes of myeloma B-cell clone origin, in addition to nuclei without these translocations, by using combined FISH and immunohistochemistry on bone sections. These nuclei of malignant origin are transcriptionally active and appear fully integrated amongst the other nuclei. The contribution of malignant nuclei to the osteoclast population analysed in this study was greater than 30%. Osteoclast-myeloma clone hybrids contained more nuclei than normal osteoclasts and their occurrence correlated with the proximity of myeloma cells. Similar hybrid cells were generated in myeloma cell-osteoclast co-cultures, as revealed by tracing myeloma nuclei using translocations, bromo-deoxyuridine, or the Y chromosome of male myeloma cells in female osteoclasts. These observations indicate that hybrid cells can originate through fusion between myeloma cells and osteoclasts. In conclusion, malignant cells contribute significantly to the formation of bone-resorbing osteoclasts in multiple myeloma. Osteoclast-myeloma clone hybrids reflect a previously unrecognized mechanism of bone destruction in which malignant cells participate directly. The possibility that malignant cells corrupt host cells by the transfer of malignant DNA may have been underestimated to date in cancer research.

Human papillomavirus genotype 31 does not express detectable microRNA levels during latent or productive virus replication.

It has recently become clear that several pathogenic DNA viruses express virally encoded microRNAs in infected cells. In particular, numerous microRNAs have been identified in a range of human and animal herpesviruses, and individual microRNAs have also been identified in members of the polyoma- and adenovirus families. Although their functions remain largely unknown, it seems likely that these viral microRNAs play an important role in viral replication in vivo. Here we present an analysis of the microRNAs expressed in human cells during the latent and productive phases of the human papillomavirus genotype 31 (HPV31) replication cycle. Although over 500 cellular microRNAs were cloned and identified, not a single HPV31-specific microRNA was obtained. We therefore concluded that HPV31, and possibly human papillomaviruses in general, does not express viral microRNAs.

Age‐associated microRNA expression in human peripheral blood is associated with all‐cause mortality and age‐related traits

Recent studies provide evidence of correlations of DNA methylation and expression of protein-coding genes with human aging. The relations of microRNA expression with age and age-related clinical outcomes have not been characterized thoroughly. We explored associations of age with whole-blood microRNA expression in 5221 adults and identified 127 microRNAs that were differentially expressed by age at P < 3.3 × 10-4 (Bonferroni-corrected). Most microRNAs were underexpressed in older individuals. Integrative analysis of microRNA and mRNA expression revealed changes in age-associated mRNA expression possibly driven by age-associated microRNAs in pathways that involve RNA processing, translation, and immune function. We fitted a linear model to predict 'microRNA age' that incorporated expression levels of 80 microRNAs. MicroRNA age correlated modestly with predicted age from DNA methylation (r = 0.3) and mRNA expression (r = 0.2), suggesting that microRNA age may complement mRNA and epigenetic age prediction models. We used the difference between microRNA age and chronological age as a biomarker of accelerated aging (Δage) and found that Δage was associated with all-cause mortality (hazards ratio 1.1 per year difference, P = 4.2 × 10-5 adjusted for sex and chronological age). Additionally, Δage was associated with coronary heart disease, hypertension, blood pressure, and glucose levels. In conclusion, we constructed a microRNA age prediction model based on whole-blood microRNA expression profiling. Age-associated microRNAs and their targets have potential utility to detect accelerated aging and to predict risks for age-related diseases.

Isolation of major pancreatic cell types and long-term culture-initiating cells using novel human surface markers.

We have developed a novel panel of cell-surface markers for the isolation and study of all major cell types of the human pancreas. Hybridomas were selected after subtractive immunization of Balb/C mice with intact or dissociated human islets and assessed for cell-type specificity and cell-surface reactivity by immunohistochemistry and flow cytometry. Antibodies were identified by specific binding of surface antigens on islet (panendocrine or alpha-specific) and nonislet pancreatic cell subsets (exocrine and duct). These antibodies were used individually or in combination to isolate populations of alpha, beta, exocrine, or duct cells from primary human pancreas by FACS and to characterize the detailed cell composition of human islet preparations. They were also employed to show that human islet expansion cultures originated from nonendocrine cells and that insulin expression levels could be increased to up to 1% of normal islet cells by subpopulation sorting and overexpression of the transcription factors Pdx-1 and ngn3, an improvement over previous results with this culture system. These methods permit the analysis and isolation of functionally distinct pancreatic cell populations with potential for cell therapy.

Immunogenicity of induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs), reprogrammed from somatic cells with defined factors, hold great promise for regenerative medicine as the renewable source of autologous cells. Whereas it has been generally assumed that these autologous cells should be immune-tolerated by the recipient from whom the iPSCs are derived, their immunogenicity has not been vigorously examined. We show here that, whereas embryonic stem cells (ESCs) derived from inbred C57BL/6 (B6) mice can efficiently form teratomas in B6 mice without any evident immune rejection, the allogeneic ESCs from 129/SvJ mice fail to form teratomas in B6 mice due to rapid rejection by recipients. B6 mouse embryonic fibroblasts (MEFs) were reprogrammed into iPSCs by either retroviral approach (ViPSCs) or a novel episomal approach (EiPSCs) that causes no genomic integration. In contrast to B6 ESCs, teratomas formed by B6 ViPSCs were mostly immune-rejected by B6 recipients. In addition, the majority of teratomas formed by B6 EiPSCs were immunogenic in B6 mice with T cell infiltration, and apparent tissue damage and regression were observed in a small fraction of teratomas. Global gene expression analysis of teratomas formed by B6 ESCs and EiPSCs revealed a number of genes frequently overexpressed in teratomas derived from EiPSCs, and several such gene products were shown to contribute directly to the immunogenicity of the B6 EiPSC-derived cells in B6 mice. These findings indicate that, in contrast to derivatives of ESCs, abnormal gene expression in some cells differentiated from iPSCs can induce T-cell-dependent immune response in syngeneic recipients. Therefore, the immunogenicity of therapeutically valuable cells derived from patient-specific iPSCs should be evaluated before any clinic application of these autologous cells into the patients.

Single-RNA counting reveals alternative modes of gene expression in yeast

Proper execution of transcriptional programs is a key requirement of gene expression regulation, demanding accurate control of timing and amplitude. How precisely the transcription machinery fulfills this task is not known. Using an in situ hybridization approach that detects single mRNA molecules, we measured mRNA abundance and transcriptional activity within single Saccharomyces cerevisiae cells. We found that expression levels for particular genes are higher than initially reported and can vary substantially among cells. However, variability for most constitutively expressed genes is unexpectedly small. Combining single-transcript measurements with computational modeling indicates that low expression variation is achieved by transcribing genes using single transcription-initiation events that are clearly separated in time, rather than by transcriptional bursts. In contrast, PDR5, a gene regulated by the transcription coactivator complex SAGA, is expressed using transcription bursts, resulting in larger variation. These data directly demonstrate the existence of multiple expression modes used to modulate the transcriptome.

Developmental heterogeneity of cardiac fibroblasts does not predict pathological proliferation and activation.

RATIONALE Fibrosis is mediated partly by extracellular matrix-depositing fibroblasts in the heart. Although these mesenchymal cells are reported to have multiple embryonic origins, the functional consequence of this heterogeneity is unknown. OBJECTIVE We sought to validate a panel of surface markers to prospectively identify cardiac fibroblasts. We elucidated the developmental origins of cardiac fibroblasts and characterized their corresponding phenotypes. We also determined proliferation rates of each developmental subset of fibroblasts after pressure overload injury. METHODS AND RESULTS We showed that Thy1(+)CD45(-)CD31(-)CD11b(-)Ter119(-) cells constitute the majority of cardiac fibroblasts. We characterized these cells using flow cytometry, epifluorescence and confocal microscopy, and transcriptional profiling (using reverse transcription polymerase chain reaction and RNA-seq). We used lineage tracing, transplantation studies, and parabiosis to show that most adult cardiac fibroblasts derive from the epicardium, a minority arises from endothelial cells, and a small fraction from Pax3-expressing cells. We did not detect generation of cardiac fibroblasts by bone marrow or circulating cells. Interestingly, proliferation rates of fibroblast subsets on injury were identical, and the relative abundance of each lineage remained the same after injury. The anatomic distribution of fibroblast lineages also remained unchanged after pressure overload. Furthermore, RNA-seq analysis demonstrated that Tie2-derived and Tbx18-derived fibroblasts within each operation group exhibit similar gene expression profiles. CONCLUSIONS The cellular expansion of cardiac fibroblasts after transaortic constriction surgery was not restricted to any single developmental subset. The parallel proliferation and activation of a heterogeneous population of fibroblasts on pressure overload could suggest that common signaling mechanisms stimulate their pathological response.

Seroepidemiologic survey of captive Old-World primates for antibodies to human and simian retroviruses, and isolation of a lentivirus from sooty mangabeys (Cercocebus atys).

Sera from 526 Old-World monkeys and apes, representing 50 species and 20 genera and living in US zoos and vivaria, were screened for antibodies to HTLV-I, HTLV-III/LAV, and simian-AIDS retrovirus, type I (SRV-I). Sera were screened initially by ELISA, and ELISA-positive sera, as well as ELISA-negative sera from cage contacts, were further tested by Western blotting. A large number of false-positive and a small number of false-negative ELISA sera were identified. Although most true positive reactions were directed to a single retrovirus, a number of individuals from 4 species were positive for more than one retrovirus. Specific seroreactivity to HTLV-I was found in 39/526 (7%) animals of 15 species. True positive reactions to SRV-I were found in 21/516 (4%) animals, including talapoins and 2 species of macaques. Specific serologic reactions to HTLV-III/LAV were detected in 23/526 (4%) monkeys. Many of the HTLV-III/LAV seropositive animals were from one mixed-species zoo exhibit, containing sooty mangabeys, mandrills, Kolb's guenons, and talapoins. A type D virus was isolated from the blood of 3/10 SRV-I antibody-positive Tonkeana macaques, but from none of 11 seropositive talapoins. A lentivirus was isolated from the blood of 4/7 HTLV-III/LAV seropositive sooty mangabeys, but not from seropositive talapoins in the same exhibit or from 2 seropositive colobus from another zoo. The sooty mangabey lentivirus produced generalized lymphadenopathy, leukopenia, and decreased levels of T4 lymphocytes in 2 experimentally infected rhesus macaques.

Immune complex relay by subcapsular sinus macrophages and non-cognate B cells drives antibody affinity maturation

Subcapsular sinus (SCS) macrophages capture antigens from lymph and present them intact for B cell encounter and follicular delivery. However, the properties of SCS macrophages are poorly defined. Here we show SCS macrophage development depended on lymphotoxin-alpha1beta2, and the cells had low lysosomal enzyme expression and retained opsonized antigens on their surface. Intravital imaging revealed immune complexes moving along macrophage processes into the follicle. Moreover, noncognate B cells relayed antigen opsonized by newly produced antibodies from the subcapsular region to the germinal center, and affinity maturation was impaired when this transport process was disrupted. Thus, we characterize SCS macrophages as specialized antigen-presenting cells functioning at the apex of an antigen transport chain that promotes humoral immunity.

Human cytomegalovirus encodes an MHC class I-like molecule (UL142) that functions to inhibit NK cell lysis.

Clinical and low passage strains of human CMV (HCMV) encode an additional MHC class I-related molecule UL142, in addition to the previously described UL18. The UL142 open reading frame is encoded within the ULb' region which is missing from a number of common high passage laboratory strains. Cells expressing UL142 following transfection, and fibroblasts infected with a recombinant adenovirus-expressing UL142, were used to screen both polyclonal NK cells and NK cell clones, in a completely autologous system. Analysis of 100 NK cell clones derived from five donors, revealed 23 clones that were inhibited by fibroblasts expressing UL142 alone. Small-interfering RNA-mediated knockdown of UL142 mRNA expression in HCMV-infected cells resulted in increased sensitivity to lysis. From these data we conclude that UL142 is a novel HCMV-encoded MHC class I-related molecule which inhibits NK cell killing in a clonally dependent manner.

Radiation and Dual Checkpoint Blockade Activates Non-Redundant Immune Mechanisms in Cancer

Immune checkpoint inhibitors result in impressive clinical responses, but optimal results will require combination with each other and other therapies. This raises fundamental questions about mechanisms of non-redundancy and resistance. Here we report major tumour regressions in a subset of patients with metastatic melanoma treated with an anti-CTLA4 antibody (anti-CTLA4) and radiation, and reproduced this effect in mouse models. Although combined treatment improved responses in irradiated and unirradiated tumours, resistance was common. Unbiased analyses of mice revealed that resistance was due to upregulation of PD-L1 on melanoma cells and associated with T-cell exhaustion. Accordingly, optimal response in melanoma and other cancer types requires radiation, anti-CTLA4 and anti-PD-L1/PD-1. Anti-CTLA4 predominantly inhibits T-regulatory cells (Treg cells), thereby increasing the CD8 T-cell to Treg (CD8/Treg) ratio. Radiation enhances the diversity of the T-cell receptor (TCR) repertoire of intratumoral T cells. Together, anti-CTLA4 promotes expansion of T cells, while radiation shapes the TCR repertoire of the expanded peripheral clones. Addition of PD-L1 blockade reverses T-cell exhaustion to mitigate depression in the CD8/Treg ratio and further encourages oligoclonal T-cell expansion. Similarly to results from mice, patients on our clinical trial with melanoma showing high PD-L1 did not respond to radiation plus anti-CTLA4, demonstrated persistent T-cell exhaustion, and rapidly progressed. Thus, PD-L1 on melanoma cells allows tumours to escape anti-CTLA4-based therapy, and the combination of radiation, anti-CTLA4 and anti-PD-L1 promotes response and immunity through distinct mechanisms.

Improved variant discovery through local re-alignment of short-read next-generation sequencing data using SRMA

A primary component of next-generation sequencing analysis is to align short reads to a reference genome, with each read aligned independently. However, reads that observe the same non-reference DNA sequence are highly correlated and can be used to better model the true variation in the target genome. A novel short-read micro realigner, SRMA, that leverages this correlation to better resolve a consensus of the underlying DNA sequence of the targeted genome is described here.

Polymorphisms in EGFR and IL28B are associated with spontaneous clearance in an HCV-infected iranian population

Although most hepatitis C virus (HCV)-infected individuals develop chronic infection, about 25% of them are able to clear the virus spontaneously without any therapeutic intervention. The aim of the present study was to identify genes associated with spontaneous HCV clearance in a population of Iranian patients. We genotyped 110 single-nucleotide polymorphisms (SNPs) in 59 selected—candidate—genes in a cohort of 107 HCV-infected participants who spontaneously cleared the infection and 176 participants whose infection persisted. Three out of the 110 SNPs were found to be associated with HCV outcome (P-values<0.03). rs11506105 in EGFR (epidermal growth factor receptor gene), and rs11881222 and rs12979860 in IL28B (interferon-λ3 gene). Multivariate logistic regression of the three markers showed that the A/A genotypes in both rs11506105 (EFGR) and rs11881222 (IL28B), and the C/C genotype in rs12979860 (IL28B) are associated with HCV clearance (recessive model: odds ratio (OR)=2.06, 95% confidence interval (95% CI)=1.09–3.88, P=0.025; OR=2.09, 95% CI=1.23–3.60, P=0.007; and OR=1.95, 95% CI=1.15–3.35, P=0.014 for rs11506105, rs12979860 and rs11881222, respectively). In conclusion, EGFR and IL28B SNPs are strong independent predictive markers of spontaneous viral clearance.

Visualization of cell cycle in mouse embryos with Fucci2 reporter directed by Rosa26 promoter.

Fucci technology makes possible the distinction between live cells in the G(1) and S/G(2)/M phases by dual-color imaging. This technology relies upon ubiquitylation-mediated proteolysis, and transgenic mice expressing Fucci provide a powerful model system with which to study the coordination of the cell cycle and development. The mice were initially generated using the CAG promoter; lines expressing the G(1) and S/G(2)/M phase probes that emitted orange (mKO2) and green (mAG) fluorescence, respectively, were separately constructed. Owing to cell type-biased strength of the CAG promoter as well as the positional effects of random transgenesis, however, we noticed some variability in Fucci expression levels. To control more reliably the expression of cell cycle probes, we used different genetic approaches to create two types of reporter mouse lines with Fucci2 and Rosa26 transcriptional machinery. Fucci2 is a recently developed Fucci derivative, which emits red (mCherry) and green (mVenus) fluorescence and provides better color contrast than Fucci. A new transgenic line, R26p-Fucci2, utilizes the Rosa26 promoter and harbors the G(1) and S/G(2)/M phase probes in a single transgene to preserve their co-inheritance. In the other R26R-Fucci2 approach, the two probes are incorporated into Rosa26 locus conditionally. The Cre-mediated loxP recombination technique thus allows researchers to design cell-type-specific Fucci2 expression. By performing time-lapse imaging experiments using R26p-Fucci2 and R26-Fucci2 in which R26R-Fucci2 had undergone germline loxP recombination, we demonstrated the great promise of these mouse reporters for studying cell cycle behavior in vivo.

Experimental Evolution of an RNA Virus in Wild Birds: Evidence for Host-Dependent Impacts on Population Structure and Competitive Fitness

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.

Apolipoprotein E genotype and concussion in college athletes.

OBJECTIVE To evaluate the association between apolipoprotein E (APOE) polymorphisms (E2, C/T Arg158Cys; E4, T/C Cys112Arg; and promoter, g-219t) and the history of concussion in college athletes. We hypothesized that carrying 1 or more APOE rare (or minor) allele assessed in this study would be associated with having a history of 1 or more concussions. DESIGN Multicenter cross-sectional study. SETTING University athletic facilities. PARTICIPANTS One hundred ninety-six male football (n = 163) and female soccer (n = 33) college athletes volunteered. INTERVENTIONS Written concussion history questionnaire and saliva samples for genotyping. MAIN OUTCOME MEASURES Self-reported history of a documented concussion and rare APOE genotype (E2, E4, promoter). RESULTS There was a significant association (Wald χ² = 3.82; P = 0.05; odds ratio = 9.8) between carrying all APOE rare alleles and the history of a previous concussion. There was also a significant association (Wald χ² = 3.96, P = 0.04, odds ratio = 8.4) between carrying the APOE promoter minor allele and experiencing 2 or more concussions. CONCLUSIONS Carriers of all 3 APOE rare (or minor) alleles assessed in this study were nearly 10 times more likely to report a previous concussion and may be at a greater risk of concussion versus noncarriers. Promoter minor allele carriers were 8.4 times more likely to report multiple concussions and may be at a greater risk of multiple concussions versus noncarriers. Research involving larger samples of individuals with multiple concussions and carriers of multiple APOE rare alleles is warranted.

Generation of large numbers of dendritic cells from mouse bone marrow cultures supplemented with granulocyte/macrophage colony-stimulating factor

Antigen-presenting, major histocompatibility complex (MHC) class II-rich dendritic cells are known to arise from bone marrow. However, marrow lacks mature dendritic cells, and substantial numbers of proliferating less-mature cells have yet to be identified. The methodology for inducing dendritic cell growth that was recently described for mouse blood now has been modified to MHC class II-negative precursors in marrow. A key step is to remove the majority of nonadherent, newly formed granulocytes by gentle washes during the first 2-4 d of culture. This leaves behind proliferating clusters that are loosely attached to a more firmly adherent "stroma. " At days 4-6 the clusters can be dislodged, isolated by 1-g sedimentation, and upon reculture, large numbers of dendritic cells are released. The latter are readily identified on the basis of their distinct cell shape, ultrastructure, and repertoire of antigens, as detected with a panel of monoclonal antibodies. The dendritic cells express high levels of MHC class II products and act as powerful accessory cells for initiating the mixed leukocyte reaction. Neither the clusters nor mature dendritic cells are generated if macrophage colony-stimulating factor rather than granulocyte/macrophage colony-stimulating factor (GM-CSF) is applied. Therefore, GM-CSF generates all three lineages of myeloid cells (granulocytes, macrophages, and dendritic cells). Since > 5 x 10(6) dendritic cells develop in 1 wk from precursors within the large hind limb bones of a single animal, marrow progenitors can act as a major source of dendritic cells. This feature should prove useful for future molecular and clinical studies of this otherwise trace cell type.

Key questions about the checkpoint blockade-are microRNAs an answer?

The introduction of immune-checkpoint blockade in the cancer therapy led to a paradigm change of the management of late stage cancers. There are already multiple FDA approved checkpoint inhibitors and many other agents are undergoing phase 2 and early phase 3 clinical trials. The therapeutic indication of immune checkpoint inhibitors expanded in the last years, but still remains unclear who can benefit. MicroRNAs are small RNAs with no coding potential. By complementary pairing to the 3' untranslated region of messenger RNA, microRNAs exert posttranscriptional control of protein expression. A network of microRNAs directly and indirectly controls the expression of checkpoint receptors and several microRNAs can target multiple checkpoint molecules, mimicking the therapeutic effect of a combined immune checkpoint blockade. In this review, we will describe the microRNAs that control the expression of immune checkpoints and we will present four specific issues of the immune checkpoint therapy in cancer: (1) imprecise therapeutic indication, (2) difficult response evaluation, (3) numerous immunologic adverse-events, and (4) the absence of response to immune therapy. Finally, we propose microRNAs as possible solutions for these pitfalls. We consider that in the near future microRNAs could become important therapeutic partners of the immune checkpoint therapy.

Spatial and Temporal Clustering of Dengue Virus Transmission in Thai Villages

BACKGROUND Transmission of dengue viruses (DENV), the leading cause of arboviral disease worldwide, is known to vary through time and space, likely owing to a combination of factors related to the human host, virus, mosquito vector, and environment. An improved understanding of variation in transmission patterns is fundamental to conducting surveillance and implementing disease prevention strategies. To test the hypothesis that DENV transmission is spatially and temporally focal, we compared geographic and temporal characteristics within Thai villages where DENV are and are not being actively transmitted. METHODS AND FINDINGS Cluster investigations were conducted within 100 m of homes where febrile index children with (positive clusters) and without (negative clusters) acute dengue lived during two seasons of peak DENV transmission. Data on human infection and mosquito infection/density were examined to precisely (1) define the spatial and temporal dimensions of DENV transmission, (2) correlate these factors with variation in DENV transmission, and (3) determine the burden of inapparent and symptomatic infections. Among 556 village children enrolled as neighbors of 12 dengue-positive and 22 dengue-negative index cases, all 27 DENV infections (4.9% of enrollees) occurred in positive clusters (p < 0.01; attributable risk [AR] = 10.4 per 100; 95% confidence interval 1-19.8 per 100]. In positive clusters, 12.4% of enrollees became infected in a 15-d period and DENV infections were aggregated centrally near homes of index cases. As only 1 of 217 pairs of serologic specimens tested in positive clusters revealed a recent DENV infection that occurred prior to cluster initiation, we attribute the observed DENV transmission subsequent to cluster investigation to recent DENV transmission activity. Of the 1,022 female adult Ae. aegypti collected, all eight (0.8%) dengue-infected mosquitoes came from houses in positive clusters; none from control clusters or schools. Distinguishing features between positive and negative clusters were greater availability of piped water in negative clusters (p < 0.01) and greater number of Ae. aegypti pupae per person in positive clusters (p = 0.04). During primarily DENV-4 transmission seasons, the ratio of inapparent to symptomatic infections was nearly 1:1 among child enrollees. Study limitations included inability to sample all children and mosquitoes within each cluster and our reliance on serologic rather than virologic evidence of interval infections in enrollees given restrictions on the frequency of blood collections in children. CONCLUSIONS Our data reveal the remarkably focal nature of DENV transmission within a hyperendemic rural area of Thailand. These data suggest that active school-based dengue case detection prompting local spraying could contain recent virus introductions and reduce the longitudinal risk of virus spread within rural areas. Our results should prompt future cluster studies to explore how host immune and behavioral aspects may impact DENV transmission and prevention strategies. Cluster methodology could serve as a useful research tool for investigation of other temporally and spatially clustered infectious diseases.

Transmission of atherosclerosis susceptibility with gut microbial transplantation.

Recent studies indicate both clinical and mechanistic links between atherosclerotic heart disease and intestinal microbial metabolism of certain dietary nutrients producing trimethylamine N-oxide (TMAO). Here we test the hypothesis that gut microbial transplantation can transmit choline diet-induced TMAO production and atherosclerosis susceptibility. First, a strong association was noted between atherosclerotic plaque and plasma TMAO levels in a mouse diversity panel (n = 22 strains, r = 0.38; p = 0.0001). An atherosclerosis-prone and high TMAO-producing strain, C57BL/6J, and an atherosclerosis-resistant and low TMAO-producing strain, NZW/LacJ, were selected as donors for cecal microbial transplantation into apolipoprotein e null mice in which resident intestinal microbes were first suppressed with antibiotics. Trimethylamine (TMA) and TMAO levels were initially higher in recipients on choline diet that received cecal microbes from C57BL/6J inbred mice; however, durability of choline diet-dependent differences in TMA/TMAO levels was not maintained to the end of the study. Mice receiving C57BL/6J cecal microbes demonstrated choline diet-dependent enhancement in atherosclerotic plaque burden as compared with recipients of NZW/LacJ microbes. Microbial DNA analyses in feces and cecum revealed transplantation of donor microbial community features into recipients with differences in taxa proportions between donor strains that were transmissible to recipients and that tended to show coincident proportions with TMAO levels. Proportions of specific taxa were also identified that correlated with plasma TMAO levels in donors and recipients and with atherosclerotic lesion area in recipients. Atherosclerosis susceptibility may be transmitted via transplantation of gut microbiota. Gut microbes may thus represent a novel therapeutic target for modulating atherosclerosis susceptibility.

Increased density of tumor-associated macrophages is associated with decreased survival in advanced thyroid cancer.

Thyroid cancers are infiltrated with tumor-associated macrophages (TAMs), yet their role in cancer progression is not known. The objectives of this study were to characterize the density of TAMs in well-differentiated (WDTC), poorly differentiated (PDTC), and anaplastic thyroid cancers (ATC) and to correlate TAM density with clinicopathologic parameters. Immunohistochemistry was performed on tissue microarray sections from WDTC (n=33), PDTC (n=37), and ATC (n=20) using macrophage-specific markers. Electronic medical records were used to gather clinical and pathologic data. Follow-up information of PDTC patients was available for 0-12 years. In total, 9 out of 33 WDTC (27%), 20 out of 37 PDTC (54%), and 19 out of 20 ATC (95%) had an increased density of CD68(+) TAMs (> or = 10 per 0.28 mm(2); WDTC versus PDTC, P=0.03; WDTC versus ATC, P<0.0001; PDTC versus ATC, P<0.002). Increased TAMs in PDTC was associated with capsular invasion (P=0.034), extrathyroidal extension (P=0.009), and decreased cancer-related survival (P=0.009) compared with PDTC with a low density of TAMs. In conclusion, the density of TAMs is increased in advanced thyroid cancers. The presence of a high density of TAMs in PDTC correlates with invasion and decreased cancer-related survival. These results suggest that TAMs may facilitate tumor progression. As novel therapies directed against thyroid tumor cell-specific targets are being tested, the potential role of TAMs as potential modulators of the thyroid cancer behavior will need to be considered.

Circulating plasma microRNAs as a screening method for detection of colorectal adenomas.

BACKGROUND MicroRNAs (miRNAs) are small non-coding RNA molecules. Reduced or increased levels of specific miRNAs are observed in colon and other cancers, supporting their role in carcinogenesis. Detection of colorectal polyps is the cornerstone of the Bowel Cancer Screening Programme in the UK. However, uptake of screening nationally remains under 60%. We aimed to see whether circulating plasma miRNAs can be used to screen for patients with colorectal polyps, adenomas, or both. METHODS Blood samples were taken from patients from the Bowel Cancer Screening Programme (asymptomatic but faecal occult blood testing [FOBt] positive). Plasma RNA was extracted, target miRNAs (19a, 98, 146b, 186, 191, 222*, 331-5p, 452, 625, 664, 1247) were identified on pooled case miRNA assay cards, and miRNA fraction was quantified by quantitative RT-PCR assay. Results were compared with endoscopy reports and with histology of any polyps identified and removed. Analysis was done with Excel (2011) and SPSS (version 20) software. FINDINGS 210 patients were included (117 with polyps, 12 with cancer, 81 healthy controls [FOBt positive]). The miRNA panel showed significant differences in expression (on t testing) for patients compared with controls for those with polyps, cancer, or both (miR-19a, p=0·0184; miR-98, p=0·0206; miR-146b, p=0·0029; miR-186, p=0·0006; miR-62,5 p=0·0008), polyps (miR-19a, p=0·0233; miR-98, p=0·0224; miR-146b, p=0·003; miR-186, p=0·0004; miR-625, p=0·001), adenomas (miR-19a, p=0·0339; miR-98, p=0·0266; miR-146b, p=0·0045; miR-186, p=0·0008; miR-625, p=0·0049), multiple adenomas (both sides of colon; miR-146b, p=0·0194; miR-186, p=0·0226; miR-625, p=0·0013), and right-sided adenomas (miR-98, p=0·031; miR-146b, p=0·0076; miR-186, p=0·0041; miR-331-5p, p=0·0142; miR-625, p=0·0049). Receiver operating characteristic analysis showed sensitivity of 60% or more, and specificity of 86% or more for men with polyps, men with adenomas, all patients with haemorrhoids or diverticulosis and polyps, and all patients with haemorrhoids or diverticulosis and adenomas. INTERPRETATION The target miRNAs that we identified showed significant differences in expression levels for patients with polyps and patients with adenomas from controls. Use of this panel has potential as a screening test. FUNDING Bowel Disease Research Foundation.

Integrated molecular genetic profiling of pediatric high-grade gliomas reveals key differences with the adult disease.

PURPOSE To define copy number alterations and gene expression signatures underlying pediatric high-grade glioma (HGG). PATIENTS AND METHODS We conducted a high-resolution analysis of genomic imbalances in 78 de novo pediatric HGGs, including seven diffuse intrinsic pontine gliomas, and 10 HGGs arising in children who received cranial irradiation for a previous cancer using single nucleotide polymorphism microarray analysis. Gene expression was analyzed with gene expression microarrays for 53 tumors. Results were compared with publicly available data from adult tumors. RESULTS Significant differences in copy number alterations distinguish childhood and adult glioblastoma. PDGFRA was the predominant target of focal amplification in childhood HGG, including diffuse intrinsic pontine gliomas, and gene expression analyses supported an important role for deregulated PDGFRalpha signaling in pediatric HGG. No IDH1 hotspot mutations were found in pediatric tumors, highlighting molecular differences with adult secondary glioblastoma. Pediatric and adult glioblastomas were clearly distinguished by frequent gain of chromosome 1q (30% v 9%, respectively) and lower frequency of chromosome 7 gain (13% v 74%, respectively) and 10q loss (35% v 80%, respectively). PDGFRA amplification and 1q gain occurred at significantly higher frequency in irradiation-induced tumors, suggesting that these are initiating events in childhood gliomagenesis. A subset of pediatric HGGs showed minimal copy number changes. CONCLUSION Integrated molecular profiling showed substantial differences in the molecular features underlying pediatric and adult HGG, indicating that findings in adult tumors cannot be simply extrapolated to younger patients. PDGFRalpha may be a useful target for pediatric HGG, including diffuse pontine gliomas.

Haplotype estimation using sequencing reads.

High-throughput sequencing technologies produce short sequence reads that can contain phase information if they span two or more heterozygote genotypes. This information is not routinely used by current methods that infer haplotypes from genotype data. We have extended the SHAPEIT2 method to use phase-informative sequencing reads to improve phasing accuracy. Our model incorporates the read information in a probabilistic model through base quality scores within each read. The method is primarily designed for high-coverage sequence data or data sets that already have genotypes called. One important application is phasing of single samples sequenced at high coverage for use in medical sequencing and studies of rare diseases. Our method can also use existing panels of reference haplotypes. We tested the method by using a mother-father-child trio sequenced at high-coverage by Illumina together with the low-coverage sequence data from the 1000 Genomes Project (1000GP). We found that use of phase-informative reads increases the mean distance between switch errors by 22% from 274.4 kb to 328.6 kb. We also used male chromosome X haplotypes from the 1000GP samples to simulate sequencing reads with varying insert size, read length, and base error rate. When using short 100 bp paired-end reads, we found that using mixtures of insert sizes produced the best results. When using longer reads with high error rates (5-20 kb read with 4%-15% error per base), phasing performance was substantially improved.

JTK_CYCLE: an efficient nonparametric algorithm for detecting rhythmic components in genome-scale data sets.

Circadian rhythms are oscillations of physiology, behavior, and metabolism that have period lengths near 24 hours. In several model organisms and humans, circadian clock genes have been characterized and found to be transcription factors. Because of this, researchers have used microarrays to characterize global regulation of gene expression and algorithmic approaches to detect cycling. This article presents a new algorithm, JTK_CYCLE, designed to efficiently identify and characterize cycling variables in large data sets. Compared with COSOPT and the Fisher's G test, two commonly used methods for detecting cycling transcripts, JTK_CYCLE distinguishes between rhythmic and nonrhythmic transcripts more reliably and efficiently. JTK_CYCLE's increased resistance to outliers results in considerably greater sensitivity and specificity. Moreover, JTK_CYCLE accurately measures the period, phase, and amplitude of cycling transcripts, facilitating downstream analyses. Finally, JTK_CYCLE is several orders of magnitude faster than COSOPT, making it ideal for large-scale data sets. JTK_CYCLE was used to analyze legacy data sets including NIH3T3 cells, which have comparatively low amplitude oscillations. JTK_CYCLE's improved power led to the identification of a novel cluster of RNA-interacting genes whose abundance is under clear circadian regulation. These data suggest that JTK_CYCLE is an ideal tool for identifying and characterizing oscillations in genome-scale data sets.

Banana contains a diverse array of endogenous badnaviruses.

Banana streak disease is caused by several distinct badnavirus species, one of which is Banana streak Obino l'Ewai virus. Banana streak Obino l'Ewai virus has severely hindered international banana (Musa spp.) breeding programmes, as new hybrids are frequently infected with this virus, curtailing any further exploitation. This infection is thought to arise from viral DNA integrated in the nuclear genome of Musa balbisiana (B genome), one of the wild species contributing to many of the banana cultivars currently grown. In order to determine whether the DNA of other badnavirus species is integrated in the Musa genome, PCR-amplified DNA fragments from Musa acuminata, M. balbisiana and Musa schizocarpa, as well as cultivars 'Obino l'Ewai' and 'Klue Tiparot', were cloned. In total, 103 clones were sequenced and all had similarity to open reading frame III in the badnavirus genome, although there was remarkable variation, with 36 distinct sequences being recognized with less than 85 % nucleotide identity to each other. There was no commonality in the sequences amplified from M. acuminata and M. balbisiana, suggesting that integration occurred following the separation of these species. Analysis of rates of non-synonymous and synonymous substitution suggested that the integrated sequences evolved under a high degree of selective constraint as might be expected for a living badnavirus, and that each distinct sequence resulted from an independent integration event.

Gene expression in human skeletal muscle: alternative normalization method and effect of repeated biopsies

The reverse transcriptase-polymerase chain reaction (RT-PCR) method has lately become widely used to determine transcription and mRNA content in rodent and human muscle samples. However, the common use of endogenous controls for correcting for variance in cDNA between samples is not optimal. Specifically, we investigated (1) a new normalization method based on determining the cDNA content by the flourophores PicoGreen and OliGreen, (2) effect of repeated muscle biopsies on mRNA gene expression, and (3) the spatial heterogeneity in mRNA expression across the muscle. Standard curves using oligo standards revealed a high degree of sensitivity and linearity (2.5–45 ng; R 2>0.99) with OliGreen reagent, as was the case for OliGreen analyses with standard curves constructed from serial dilutions of representative RT samples (R 2 >0.99 for a ten times dilution range of a representative reversed transcribed (RT) sample). Likewise, PicoGreen reagent detected the RNA:DNA hybrid content in RT samples with great sensitivity. Standard curves constructed from both double-stranded lambda DNA (1–10 ng) and from serial dilutions of representative RT samples consistently resulted in linearity with R 2 >0.99. The present determination of cDNA content in reversed transcribed human skeletal muscle RNA samples by both PicoGreen and OliGreen analyses suggests that these fluorophores provide a potential alternative normalization procedure for human gene expression studies. In addition, the present study shows that multiple muscle biopsies obtained from the same muscle do not influence the mRNA response induced by an acute exercise bout for any of the genes examined.

Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing.

Knowledge of the rate and nature of spontaneous mutation is fundamental to understanding evolutionary and molecular processes. In this report, we analyze spontaneous mutations accumulated over thousands of generations by wild-type Escherichia coli and a derivative defective in mismatch repair (MMR), the primary pathway for correcting replication errors. The major conclusions are (i) the mutation rate of a wild-type E. coli strain is ~1 × 10(-3) per genome per generation; (ii) mutations in the wild-type strain have the expected mutational bias for G:C > A:T mutations, but the bias changes to A:T > G:C mutations in the absence of MMR; (iii) during replication, A:T > G:C transitions preferentially occur with A templating the lagging strand and T templating the leading strand, whereas G:C > A:T transitions preferentially occur with C templating the lagging strand and G templating the leading strand; (iv) there is a strong bias for transition mutations to occur at 5'ApC3'/3'TpG5' sites (where bases 5'A and 3'T are mutated) and, to a lesser extent, at 5'GpC3'/3'CpG5' sites (where bases 5'G and 3'C are mutated); (v) although the rate of small (≤4 nt) insertions and deletions is high at repeat sequences, these events occur at only 1/10th the genomic rate of base-pair substitutions. MMR activity is genetically regulated, and bacteria isolated from nature often lack MMR capacity, suggesting that modulation of MMR can be adaptive. Thus, comparing results from the wild-type and MMR-defective strains may lead to a deeper understanding of factors that determine mutation rates and spectra, how these factors may differ among organisms, and how they may be shaped by environmental conditions.

siRNA specificity searching incorporating mismatch tolerance data.

UNLABELLED Artificially synthesized short interfering RNAs (siRNAs) are widely used in functional genomics to knock down specific target genes. One ongoing challenge is to guarantee that the siRNA does not elicit off-target effects. Initial reports suggested that siRNAs were highly sequence-specific; however, subsequent data indicates that this is not necessarily the case. It is still uncertain what level of similarity and other rules are required for an off-target effect to be observed, and scoring schemes have not been developed to look beyond simple measures such as the number of mismatches or the number of consecutive matching bases present. We created design rules for predicting the likelihood of a non-specific effect and present a web server that allows the user to check the specificity of a given siRNA in a flexible manner using a combination of methods. The server finds potential off-target matches in the corresponding RefSeq database and ranks them according to a scoring system based on experimental studies of specificity. AVAILABILITY The server is available at http://informatics-eskitis.griffith.edu.au/SpecificityServer.

Integrated genomic profiling of endometrial carcinoma associates aggressive tumors with indicators of PI3 kinase activation.

Although 75% of endometrial cancers are treated at an early stage, 15% to 20% of these recur. We performed an integrated analysis of genome-wide expression and copy-number data for primary endometrial carcinomas with extensive clinical and histopathological data to detect features predictive of recurrent disease. Unsupervised analysis of the expression data distinguished 2 major clusters with strikingly different phenotypes, including significant differences in disease-free survival. To identify possible mechanisms for these differences, we performed a global genomic survey of amplifications, deletions, and loss of heterozygosity, which identified 11 significantly amplified and 13 significantly deleted regions. Amplifications of 3q26.32 harboring the oncogene PIK3CA were associated with poor prognosis and segregated with the aggressive transcriptional cluster. Moreover, samples with PIK3CA amplification carried signatures associated with in vitro activation of PI3 kinase (PI3K), a signature that was shared by aggressive tumors without PIK3CA amplification. Tumors with loss of PTEN expression or PIK3CA overexpression that did not have PIK3CA amplification also shared the PI3K activation signature, high protein expression of the PI3K pathway member STMN1, and an aggressive phenotype in test and validation datasets. However, mutations of PTEN or PIK3CA were not associated with the same expression profile or aggressive phenotype. STMN1 expression had independent prognostic value. The results affirm the utility of systematic characterization of the cancer genome in clinically annotated specimens and suggest the particular importance of the PI3K pathway in patients who have aggressive endometrial cancer.

Overlapping patterns of IGF2 and H19 expression during human development: biallelic IGF2 expression correlates with a lack of H19 expression.

The spatial patterns of IGF2 and H19 gene expression are strikingly similar during parts of human embryonic/fetal and early postnatal development. Notable exceptions were found with the ciliary anlage of the embryonic retina and the choroid plexus/leptomeninges, where transcripts from the IGF2 but not the H19 locus could be detected. Moreover, in contrast to the other tissue samples examined, the choroid plexus/leptomeninges expressed both parental IGF2 alleles. Whilst RNase protection analysis revealed a weak activity of the P1 promoter in the choroid plexus/leptomeninges, the P2, P3 and P4 promoters were all active wherever IGF2 was expressed. We discuss these observations with respect to a hypothesized coordinated control of the reciprocally imprinted and closely linked IGF2 and H19 loci.

Seroepidemiological survey of feline retrovirus infections in cats in Taiwan in 1993 and 1994.

In order to investigate the prevalence of infections with three feline retroviruses feline immunodeficiency virus (FIV), feline leukemia virus (FeLV) and feline syncytial virus (FSV) in Taiwan, we collected a total of 75 blood samples from cats from veterinary hospitals, a breeding cattery and a homeless shelter in 1993 and 1994. We examined the presences of anti-FIV and FSV antibodies and FeLV-p27 antigen in these samples by the indirect immunofluorescence and/or enzyme-linked immunosorbent assays. All of the serum samples positive for FIV were obtained from homeless cats and the overall FIV positive rate was 4%. The overall positive rates of FSV and FeLV were 28% and 1.3%, respectively. From these results, together with previous seroepidemiological surveys by others, it was revealed that the prevalence of FIV and FeLV infections appeared to be lower in Taiwan than in the United States or Japan. In contrast, the prevalence of FSV infection in Taiwan was as high as that in Japan.

Detection of Sugarcane bacilliform virus isolate and its influence on yield and quality of cane in Yunnan

Sugarcane bacilliform virus(SCBV) was detected by PCR from sugarcane showing chlorosis and mottle symptom from Kaiyuan,Yunnan Province. Part sequence of replicase gene of the isolate SCBV-Kaiyuan was determined. Sequence analysis indicated that the 589 bp of SCBV-Kaiyuan shared identities of 73.2%-74.0% and 83.1%-84.1% at nucleotide and amino acid levels with SCBV-Australia respectively,66.7%-68.4% and 65.6%-67.7% with SCBV-Morocco. The quality and yield of the sugarcane infected with SCBV-Kaiyuan was also investigated. The juice extraction,sucrose content,gravity purity and average stalk weight were decreased 1.55%,1.24%,2.22% and 0.26 kg in plants infected with SCBV-Kaiyuan,but reducing sugar was increased by 0.21% in infected plants.

Exploitative and Hierarchical Antagonism in a Cooperative Bacterium

Social organisms that cooperate with some members of their own species, such as close relatives, may fail to cooperate with other genotypes of the same species. Such noncooperation may take the form of outright antagonism or social exploitation. Myxococcus xanthus is a highly social prokaryote that cooperatively develops into spore-bearing, multicellular fruiting bodies in response to starvation. Here we have characterized the nature of social interactions among nine developmentally proficient strains of M. xanthus isolated from spatially distant locations. Strains were competed against one another in all possible pairwise combinations during starvation-induced development. In most pairings, at least one competitor exhibited strong antagonism toward its partner and a majority of mixes showed bidirectional antagonism that decreased total spore production, even to the point of driving whole populations to extinction. Differential response to mixing was the primary determinant of competitive superiority rather than the sporulation efficiencies of unmixed populations. In some competitive pairings, the dominant partner sporulated more efficiently in mixed populations than in clonal isolation. This finding represents a novel form of exploitation in bacteria carried out by socially competent genotypes and is the first documentation of social exploitation among natural bacterial isolates. Patterns of antagonistic superiority among these strains form a highly linear dominance hierarchy. At least some competition pairs construct chimeric, rather than segregated, fruiting bodies. The cooperative prokaryote M. xanthus has diverged into a large number of distinct social types that cooperate with clone-mates but exhibit intense antagonism toward distinct social types of the same species. Most lengthy migration events in nature may thus result in strong antagonism between migratory and resident populations, and this antagonism may have large effects on local population sizes and dynamics. Intense mutual antagonism appears to be more prevalent in this prokaryotic social species than has been observed in the eukaryotic social slime mold Dictyostelium discoideum, which also exhibits multicellular development. The finding of several cases of facultative social exploitation among these natural isolates suggests that such exploitation may occur frequently in nature in many prokaryotes with cooperative traits.

Analysis of the interaction of primate retroviruses with the human RNA interference machinery.

The question of whether retroviruses, including human immunodeficiency virus type 1 (HIV-1), interact with the cellular RNA interference machinery has been controversial. Here, we present data showing that neither HIV-1 nor human T-cell leukemia virus type 1 (HTLV-1) expresses significant levels of either small interfering RNAs or microRNAs in persistently infected T cells. We also demonstrate that the retroviral nuclear transcription factors HIV-1 Tat and HTLV-1 Tax, as well as the Tas transactivator encoded by primate foamy virus, fail to inhibit RNA interference in human cells. Moreover, the stable expression of physiological levels of HIV-1 Tat did not globally inhibit microRNA production or expression in infected human cells. Together, these data argue that HIV-1 and HTLV-1 neither induce the production of viral small interfering RNAs or microRNAs nor repress the cellular RNA interference machinery in infected cells.

Interplay between microtubule bundling and sorting factors ensures acentriolar spindle stability during C. elegans oocyte meiosis

In many species, oocyte meiosis is carried out in the absence of centrioles. As a result, microtubule organization, spindle assembly, and chromosome segregation proceed by unique mechanisms. Here, we report insights into the principles underlying this specialized form of cell division, through studies of C. elegans KLP-15 and KLP-16, two highly homologous members of the kinesin-14 family of minus-end-directed kinesins. These proteins localize to the acentriolar oocyte spindle and promote microtubule bundling during spindle assembly; following KLP-15/16 depletion, microtubule bundles form but then collapse into a disorganized array. Surprisingly, despite this defect we found that during anaphase, microtubules are able to reorganize into a bundled array that facilitates chromosome segregation. This phenotype therefore enabled us to identify factors promoting microtubule organization during anaphase, whose contributions are normally undetectable in wild-type worms; we found that SPD-1 (PRC1) bundles microtubules and KLP-18 (kinesin-12) likely sorts those bundles into a functional orientation capable of mediating chromosome segregation. Therefore, our studies have revealed an interplay between distinct mechanisms that together promote spindle formation and chromosome segregation in the absence of structural cues such as centrioles.

Generation of gene-edited rats by delivery of CRISPR/Cas9 protein and donor DNA into intact zygotes using electroporation

The generation of gene-edited animals using the CRISPRs/Cas9 system is based on microinjection into zygotes which is inefficient, time consuming and demands high technical skills. We report the optimization of an electroporation method for intact rat zygotes using sgRNAs and Cas9 protein in combination or not with ssODNs (~100 nt). This resulted in high frequency of knockouts, between 15 and 50% of analyzed animals. Importantly, using ssODNs as donor template resulted in precise knock-in mutations in 25-100% of analyzed animals, comparable to microinjection. Electroporation of long ssDNA or dsDNA donors successfully used in microinjection in the past did not allow generation of genome-edited animals despite dsDNA visualization within zygotes. Thus, simultaneous electroporation of a large number of intact rat zygotes is a rapid, simple, and efficient method for the generation of a variety of genome-edited rats.

Hoxa9 transforms primary bone marrow cells through specific collaboration with Meis1a but not Pbx1b.

Hoxa9, Meis1 and Pbx1 encode homeodomaincontaining proteins implicated in leukemic transformation in both mice and humans. Hoxa9, Meis1 and Pbx1 proteins have been shown to physically interact with each other, as Hoxa9 cooperatively binds consensus DNA sequences with Meis1 and with Pbx1, while Meis1 and Pbx1 form heterodimers in both the presence and absence of DNA. In this study, we sought to determine if Hoxa9 could transform hemopoietic cells in collaboration with either Pbx1 or Meis1. Primary bone marrow cells, retrovirally engineered to overexpress Hoxa9 and Meis1a simultaneously, induced growth factor-dependent oligoclonal acute myeloid leukemia in <3 months when transplanted into syngenic mice. In contrast, overexpression of Hoxa9, Meis1a or Pbx1b alone, or the combination of Hoxa9 and Pbx1b failed to transform these cells acutely within 6 months post-transplantation. Similar results were obtained when FDC-P1 cells, engineered to overexpress these genes, were transplanted to syngenic recipients. Thus, these studies demonstrate a selective collaboration between a member of the Hox family and one of its DNA-binding partners in transformation of hemopoietic cells.

Prediction of guide strand of microRNAs from its sequence and secondary structure

BACKGROUND MicroRNAs (miRNAs) are produced by the sequential processing of a long hairpin RNA transcript by Drosha and Dicer, an RNase III enzymes, and form transitory small RNA duplexes. One strand of the duplex, which incorporates into RNA-induced silencing complex (RISC) and silences the gene expression is called guide strand, or miRNA; while the other strand of duplex is degraded and called the passenger strand, or miRNA*. Predicting the guide strand of miRNA is important for better understanding the RNA interference pathways. RESULTS This paper describes support vector machine (SVM) models developed for predicting the guide strands of miRNAs. All models were trained and tested on a dataset consisting of 329 miRNA and 329 miRNA* pairs using five fold cross validation technique. Firstly, models were developed using mono-, di-, and tri-nucleotide composition of miRNA strands and achieved the highest accuracies of 0.588, 0.638 and 0.596 respectively. Secondly, models were developed using split nucleotide composition and achieved maximum accuracies of 0.553, 0.641 and 0.602 for mono-, di-, and tri-nucleotide respectively. Thirdly, models were developed using binary pattern and achieved the highest accuracy of 0.708. Furthermore, when integrating the secondary structure features with binary pattern, an accuracy of 0.719 was seen. Finally, hybrid models were developed by combining various features and achieved maximum accuracy of 0.799 with sensitivity 0.781 and specificity 0.818. Moreover, the performance of this model was tested on an independent dataset that achieved an accuracy of 0.80. In addition, we also compared the performance of our method with various siRNA-designing methods on miRNA and siRNA datasets. CONCLUSION In this study, first time a method has been developed to predict guide miRNA strands, of miRNA duplex. This study demonstrates that guide and passenger strand of miRNA precursors can be distinguished using their nucleotide sequence and secondary structure. This method will be useful in understanding microRNA processing and can be implemented in RNA silencing technology to improve the biological and clinical research. A web server has been developed based on SVM models described in this study (http://crdd.osdd.net:8081/RISCbinder/).

Anti-alphav integrin monoclonal antibody intetumumab enhances the efficacy of radiation therapy and reduces metastasis of human cancer xenografts in nude rats.

We previously reported that intetumumab (CNTO 95), a fully human anti-αv integrin monoclonal antibody, is a radiosensitizer in mice with xenograft tumors. Because intetumumab does not cross-react with mouse integrins, but has cross-reactivity with rat integrins, we next studied the potential combined use of radiation therapy and intetumumab in human cancer xenograft models in nude rats to assess effects on both tumor cells and the tumor microenvironment. Nude rats bearing human head and neck cancer and non-small cell lung cancer (NSCLC) xenografts were treated with intetumumab and fractionated local tumor radiotherapy. Effects on tumor growth and metastasis, blood perfusion, oxygenation, and gastrointestinal toxicity were studied. Intetumumab alone had a moderate effect on tumor growth. When combined with fractionated radiation therapy, intetumumab significantly inhibited tumor growth and produced a tumor response rate that was significantly better than with radiation therapy alone. Treatment with intetumumab also significantly reduced lung metastasis in the A549 NSCLC xenograft model. The oxygenation and blood perfusion in xenograft tumors measured by microbubble-enhanced ultrasound imaging were substantially increased after treatment with intetumumab. The combined use of intetumumab and radiation therapy reduced the microvessel density and increased apoptosis in tumor cells and the tumor microenvironment. Toxicity studies showed that treatment with intetumumab did not cause the histopathologic changes in the lungs and did not sensitize the sensitive gastrointestinal epithelium to the effect of radiation therapy. Intetumumab can potentiate the efficacy of fractionated radiation therapy in human cancer xenograft tumors in nude rats without increased toxicity.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing.

SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.

A genome-wide short hairpin RNA screening of jurkat T-cells for human proteins contributing to productive HIV-1 replication.

Short interfering RNAs (siRNAs) have been used to inhibit HIV-1 replication. The durable inhibition of HIV-1 replication by RNA interference has been impeded, however, by a high mutation rate when viral sequences are targeted and by cytotoxicity when cellular genes are knocked down. To identify cellular proteins that contribute to HIV-1 replication that can be chronically silenced without significant cytotoxicity, we employed a shRNA library that targets 54,509 human transcripts. We used this library to select a comprehensive population of Jurkat T-cell clones, each expressing a single discrete shRNA. The Jurkat clones were then infected with HIV-1. Clones that survived viral infection represent moieties silenced for a human mRNA needed for virus replication, but whose chronic knockdown did not cause cytotoxicity. Overall, 252 individual Jurkat mRNAs were identified. Twenty-two of these mRNAs were secondarily verified for their contributions to HIV-1 replication. Five mRNAs, NRF1, STXBP2, NCOA3, PRDM2, and EXOSC5, were studied for their effect on steps of the HIV-1 life cycle. We discuss the similarities and differences between our shRNA findings for HIV-1 using a spreading infection assay in human Jurkat T-cells and results from other investigators who used siRNA-based screenings in HeLa or 293T cells.

COHCAP: an integrative genomic pipeline for single-nucleotide resolution DNA methylation analysis

COHCAP (City of Hope CpG Island Analysis Pipeline) is an algorithm to analyze single-nucleotide resolution DNA methylation data produced by either an Illumina methylation array or targeted bisulfite sequencing. The goal of the COHCAP algorithm is to identify CpG islands that show a consistent pattern of methylation among CpG sites. COHCAP is currently the only DNA methylation package that provides integration with gene expression data to identify a subset of CpG islands that are most likely to regulate downstream gene expression, and it can generate lists of differentially methylated CpG islands with ∼50% concordance with gene expression from both cell line data and heterogeneous patient data. For example, this article describes known breast cancer biomarkers (such as estrogen receptor) with a negative correlation between DNA methylation and gene expression. COHCAP also provides visualization for quality control metrics, regions of differential methylation and correlation between methylation and gene expression. This software is freely available at https://sourceforge.net/projects/cohcap/.

Rapid purification of homodimer and heterodimer HIV-1 reverse transcriptase by metal chelate affinity chromatography.

We have modified an Escherichia coli vector expressing 66-kDa HIV-1 reverse transcriptase (p66) so that it simultaneously expresses this and the pol-coded protease. The twin expression cassette yields high quantities of both reverse transcriptase and protease; however, under these conditions, 50% of the over-expressed p66 reverse transcriptase is processed, resulting in accumulation of large quantities of p66/p51 enzyme. Furthermore, addition of a poly(histidine) affinity label at the amino terminus of the reverse-transcriptase-coding sequence (His-p66) permits a simple, rapid purification of milligram quantities of either p66 or p66/p51 enzyme from a crude lysate by metal chelate affinity chromatography. Purified His-p66 and His-p66/His-p51 reverse transcriptase exhibit both reverse transcriptase and RNase H activity. Purification by metal chelate chromatography of a p66/p51 enzyme wherein only the p66 component is labelled strengthens the argument for the existence of a heterodimer.

MIRA-seq for DNA methylation analysis of CpG islands

AIM To develop a reliable method for whole genome analysis of DNA methylation. MATERIALS & METHODS Genome-scale analysis of DNA methylation includes affinity-based approaches such as enrichment using methyl-CpG-binding proteins. One of these methods, the methylated-CpG island recovery assay (MIRA), is based on the high affinity of the MBD2b-MBD3L1 complex for CpG-methylated DNA. Here we provide a detailed description of MIRA and combine it with next generation sequencing platforms (MIRA-seq). RESULTS We assessed the performance of MIRA-seq and compared the data with whole genome bisulfite sequencing. CONCLUSION MIRA-seq is a reliable, genome-scale DNA methylation analysis platform for scoring DNA methylation differences at CpG-rich genomic regions. The method is not limited by primer or probe design and is cost effective.

Identification of Gene Positioning Factors Using High-Throughput Imaging Mapping

Genomes are arranged non-randomly in the 3D space of the cell nucleus. Here, we have developed HIPMap, a high-precision, high-throughput, automated fluorescent in situ hybridization imaging pipeline, for mapping of the spatial location of genome regions at large scale. High-throughput imaging position mapping (HIPMap) enabled an unbiased siRNA screen for factors involved in genome organization in human cells. We identify 50 cellular factors required for proper positioning of a set of functionally diverse genomic loci. Positioning factors include chromatin remodelers, histone modifiers, and nuclear envelope and pore proteins. Components of the replication and post-replication chromatin re-assembly machinery are prominently represented among positioning factors, and timely progression of cells through replication, but not mitosis, is required for correct gene positioning. Our results establish a method for the large-scale mapping of genome locations and have led to the identification of a compendium of cellular factors involved in spatial genome organization.

Diverse tumorigenic pathways in ovarian serous carcinoma.

This study was undertaken to analyze genetic alterations in 108 sporadic serous ovarian neoplasms to elucidate ovarian serous carcinogenesis. Our results demonstrate that K-ras mutations occur in approximately 50% of serous borderline tumors (SBTs), non-invasive micropapillary serous carcinomas (MPSCs), and invasive micropapillary serous carcinomas, which represent a morphological continuum of tumor progression. Moreover, progressive increase in the degree of allelic imbalance of chromosomes 1p, 5q, 8p, 18q, 22q, and Xp was observed comparing serous borderline tumors to noninvasive and invasive micropapillary serous carcinomas. In contrast, high-grade (conventional serous carcinoma) tumors contained wild-type K-ras in all 23 cases studied and a high frequency of allelic imbalance even in small (early) primary tumors similar to that found in advanced stage tumors. Based on these findings, we propose a dualistic model for ovarian serous carcinogenesis. One pathway involves a stepwise progression from SBT to noninvasive and then invasive MPSC. The other pathway is characterized by rapid progression from the ovarian surface epithelium or inclusion cysts to a conventional (high-grade) serous carcinoma.

[Detection of CRISPR and its relationship to drug resistance in Shigella].

OBJECTIVE To detect clustered regularly interspaced short palindromic repeats (CRISPR) in Shigella, and to analyze its relationship to drug resistance. METHODS Four pairs of primers were used for the detection of convincing CRISPR structures CRISPR-S2 and CRISPR-S4, questionable CRISPR structures CRISPR-S1 and CRISPR-S3 in 60 Shigella strains. All primers were designed using sequences in CRISPR database. CRISPR Finder was used to analyze CRISPR and susceptibilities of Shigella strains were tested by agar diffusion method. Furthermore, we analyzed the relationship between drug resistance and CRISPR-S4. RESULTS The positive rate of convincing CRISPR structures was 95%. The four CRISPR loci formed 12 spectral patterns (A-L), all of which contained convincing CRISPR structures except type K. We found one new repeat and 12 new spacers. The multi-drug resistance rate was 53. 33% . We found no significant difference between CRISPR-S4 and drug resistant. However, the repeat sequence of CRISPR-S4 in multi- or TE-resistance strains was mainly R4.1 with AC deletions in the 3' end, and the spacer sequences of CRISPR-S4 in multi-drug resistance strains were mainly Sp5.1, Sp6.1 and Sp7. CONCLUSION CRISPR was common in Shigella. Variations df repeat sequences and diversities of spacer sequences might be related to drug resistance in Shigella.

Fra-1 protooncogene regulates IL-6 expression in macrophages and promotes the generation of M2d macrophages

The tumor microenvironment (TME) plays a prominent role in the growth of tumor cells. As the major inflammatory component of the TME, M2d macrophages are educated by the TME such that they adopt an immunosuppressive role that promotes tumor metastasis and progression. Fra-1 forms activator protein-1 heterodimers with Jun partners and drives gene transcription. Fra-1 is thought to drastically induce tumorigenesis and progression. However, the functional role of Fra-1 in the generation of M2d macrophages is poorly understood to date. Here, we demonstrate that 4T1 mammary carcinoma cells, when co-cultured with RAW264.7 macrophage cells, skew the RAW264.7 macrophage cell differentiation into M2d macrophages. The 4T1 cells stimulate de novo overexpression of Fra-1 in RAW264.7 cells, and then Fra-1 binds to the interleukin 6 (IL-6) promoter to increase the production of the cytokine IL-6 in RAW264.7 cells. IL-6 acts in an autocrine fashion to skew RAW264.7 macrophage cell differentiation into M2d macrophages. These findings open new insights into how to reverse M2d macrophage-induced immune tolerance to improve the efficacy of immunotherapeutic approaches.

Thymic expression of the golli-myelin basic protein gene in the SJL/J mouse

The T cell antigen-specific repertoire is thought to be shaped by thymic expression of self molecules. Since a myelin basic protein (MBP)-like gene (golli-MBP) has been reported to be expressed by cells of the immune system, the present study was undertaken to determine whether the golli-MBP gene was expressed in the mouse thymus and, if so, to characterize transcripts of this gene in this organ. Using exon-specific primers for MBP and golli-MBP, cDNA from thymus and other tissues was amplified, and the amplified products analyzed by Southern blotting with exon-specific oligonucleotide probes. The amplified products were subcloned, and the inserts characterized by DNA sequencing. The thymic transcripts were found to contain golli-MBP exons 1, 2, 3, 5A, 5B, 5C, 6, 7, 8, and 11.

Human and other faeces as breeding media of the trachoma vector Musca sorbens.

The fly Musca sorbens Wiedemann (Diptera: Muscidae) apparently transmits Chlamydia trachomatis, causing human trachoma. The literature indicates that M. sorbens breeds predominantly in isolated human faeces on the soil surface, but not in covered pit latrines. We sought to identify breeding media of M. sorbens in a rural Gambian village endemic for trachoma. Test breeding media were presented for oviposition on soil-filled buckets and monitored for adult emergence. Musca sorbens emerged from human (6/9 trials), calf (3/9), cow (3/9), dog (2/9) and goat (1/9) faeces, but not from horse faeces, composting kitchen scraps or a soil control (0/9 of each). After adjusting for mass of medium, the greatest number of flies emerged from human faeces (1426 flies/kg). Median time for emergence was 9 (inter quartile range = 8-9.75) days post-oviposition. Of all flies emerging from faeces 81% were M. sorbens. Male and female flies emerging from human faeces were significantly larger than those from other media, suggesting that they would be more fecund and live longer than smaller flies from other sources. Female flies caught from children's eyes were of a similar size to those from human faeces, but significantly larger than those from other media. We consider that human faeces are the best larval medium for M. sorbens, although some breeding also occurs in animal faeces. Removal of human faeces from the environment, through the provision of basic sanitation, is likely to greatly reduce fly density, eye contact and hence trachoma transmission, but if faeces of other animals are present M. sorbens will persist.

Scaffold-based three-dimensional human fibroblast culture provides a structural matrix that supports angiogenesis in infarcted heart tissue.

BACKGROUND We have developed techniques to implant angiogenic patches onto the epicardium over regions of infarcted cardiac tissue to stimulate revascularization of the damaged tissue. These experiments used a scaffold-based 3D human dermal fibroblast culture (3DFC) as an epicardial patch. The 3DFC contains viable cells that secrete angiogenic growth factors and has previously been shown to stimulate angiogenic activity. The hypothesis tested was that a viable 3DFC cardiac patch would stimulate an angiogenic response within an area of infarcted cardiac tissue. METHODS AND RESULTS A coronary occlusion of a branch of the left anterior descending coronary artery was performed by thermal ligation in severe combined immunodeficient mice. 3DFCs with or without viable cells were sized to the damaged area, implanted in replicate mice onto the epicardium at the site of tissue injury, and compared with animals that received infarct surgery but no implant. Fourteen and 30 days after surgery, hearts were exposed and photographed, and tissue samples were prepared for histology and cytochemistry. Fourteen and 30 days after surgery, the damaged myocardium receiving viable 3DFC exhibited a significantly greater angiogenic response (including arterioles, venules, and capillaries) than nonviable and untreated control groups. CONCLUSIONS In this animal model, viable 3DFC stimulates angiogenesis within a region of cardiac infarction and can augment a repair response in damaged tissue. Therefore, a potential use for 3DFC is the repair of myocardial tissue damaged by infarction.

DNA methylation and healthy human aging

The process of aging results in a host of changes at the cellular and molecular levels, which include senescence, telomere shortening, and changes in gene expression. Epigenetic patterns also change over the lifespan, suggesting that epigenetic changes may constitute an important component of the aging process. The epigenetic mark that has been most highly studied is DNA methylation, the presence of methyl groups at CpG dinucleotides. These dinucleotides are often located near gene promoters and associate with gene expression levels. Early studies indicated that global levels of DNA methylation increase over the first few years of life and then decrease beginning in late adulthood. Recently, with the advent of microarray and next-generation sequencing technologies, increases in variability of DNA methylation with age have been observed, and a number of site-specific patterns have been identified. It has also been shown that certain CpG sites are highly associated with age, to the extent that prediction models using a small number of these sites can accurately predict the chronological age of the donor. Together, these observations point to the existence of two phenomena that both contribute to age-related DNA methylation changes: epigenetic drift and the epigenetic clock. In this review, we focus on healthy human aging throughout the lifetime and discuss the dynamics of DNA methylation as well as how interactions between the genome, environment, and the epigenome influence aging rates. We also discuss the impact of determining 'epigenetic age' for human health and outline some important caveats to existing and future studies.

Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag.

Human artificial chromosomes have been used to model requirements for human chromosome segregation and to explore the nature of sequences competent for centromere function. Normal human centromeres require specialized chromatin that consists of alpha satellite DNA complexed with epigenetically modified histones and centromere-specific proteins. While several types of alpha satellite DNA have been used to assemble de novo centromeres in artificial chromosome assays, the extent to which they fully recapitulate normal centromere function has not been explored. Here, we have used two kinds of alpha satellite DNA, DXZ1 (from the X chromosome) and D17Z1 (from chromosome 17), to generate human artificial chromosomes. Although artificial chromosomes are mitotically stable over many months in culture, when we examined their segregation in individual cell divisions using an anaphase assay, artificial chromosomes exhibited more segregation errors than natural human chromosomes (P < 0.001). Naturally occurring, but abnormal small ring chromosomes derived from chromosome 17 and the X chromosome also missegregate more than normal chromosomes, implicating overall chromosome size and/or structure in the fidelity of chromosome segregation. As different artificial chromosomes missegregate over a fivefold range, the data suggest that variable centromeric DNA content and/or epigenetic assembly can influence the mitotic behavior of artificial chromosomes.

Implementing the ESHRE 'poor responder' criteria in research studies: methodological implications.

The Bologna criteria for defining poor ovarian response (POR) during IVF provide a useful template for new research in this field of assisted conception. However, designing studies around the European Society for Human Reproduction and Embryology POR criteria can be methodologically challenging, as the new definition includes various POR subpopulations with diverse baseline characteristics and unknown clinical prognosis. When designing RCTs, potential result bias may be introduced if women from each subpopulation are not evenly allocated between intervention groups. In the case of small or moderate-size RCTs, a single-sequence randomization method may not ensure balanced allocation between groups. Stratified randomization methods provide an alternative methodological approach. Depending on the chosen methodology, patient characteristics and outcomes within each intervention group may be better reported according to relevant subpopulations.

An Extensive MicroRNA-Mediated Network of RNA-RNA Interactions Regulates Established Oncogenic Pathways in Glioblastoma

By analyzing gene expression data in glioblastoma in combination with matched microRNA profiles, we have uncovered a posttranscriptional regulation layer of surprising magnitude, comprising more than 248,000 microRNA (miR)-mediated interactions. These include ∼7,000 genes whose transcripts act as miR "sponges" and 148 genes that act through alternative, nonsponge interactions. Biochemical analyses in cell lines confirmed that this network regulates established drivers of tumor initiation and subtype implementation, including PTEN, PDGFRA, RB1, VEGFA, STAT3, and RUNX1, suggesting that these interactions mediate crosstalk between canonical oncogenic pathways. siRNA silencing of 13 miR-mediated PTEN regulators, whose locus deletions are predictive of PTEN expression variability, was sufficient to downregulate PTEN in a 3'UTR-dependent manner and to increase tumor cell growth rates. Thus, miR-mediated interactions provide a mechanistic, experimentally validated rationale for the loss of PTEN expression in a large number of glioma samples with an intact PTEN locus.

Candidate Causal Regulatory Effects by Integration of Expression QTLs with Complex Trait Genetic Associations

The recent success of genome-wide association studies (GWAS) is now followed by the challenge to determine how the reported susceptibility variants mediate complex traits and diseases. Expression quantitative trait loci (eQTLs) have been implicated in disease associations through overlaps between eQTLs and GWAS signals. However, the abundance of eQTLs and the strong correlation structure (LD) in the genome make it likely that some of these overlaps are coincidental and not driven by the same functional variants. In the present study, we propose an empirical methodology, which we call Regulatory Trait Concordance (RTC) that accounts for local LD structure and integrates eQTLs and GWAS results in order to reveal the subset of association signals that are due to cis eQTLs. We simulate genomic regions of various LD patterns with both a single or two causal variants and show that our score outperforms SNP correlation metrics, be they statistical (r(2)) or historical (D'). Following the observation of a significant abundance of regulatory signals among currently published GWAS loci, we apply our method with the goal to prioritize relevant genes for each of the respective complex traits. We detect several potential disease-causing regulatory effects, with a strong enrichment for immunity-related conditions, consistent with the nature of the cell line tested (LCLs). Furthermore, we present an extension of the method in trans, where interrogating the whole genome for downstream effects of the disease variant can be informative regarding its unknown primary biological effect. We conclude that integrating cellular phenotype associations with organismal complex traits will facilitate the biological interpretation of the genetic effects on these traits.

Single-Construct Polycistronic Doxycycline-Inducible Vectors Improve Direct Cardiac Reprogramming and Can Be Used to Identify the Critical Timing of Transgene Expression

Direct reprogramming is a promising approach in regenerative medicine. Overexpression of the cardiac transcription factors Gata4, Mef2c, and Tbx5 (GMT) or GMT plus Hand2 (GHMT) directly reprogram fibroblasts into cardiomyocyte-like cells (iCMs). However, the critical timing of transgene expression and the molecular mechanisms for cardiac reprogramming remain unclear. The conventional doxycycline (Dox)-inducible temporal transgene expression systems require simultaneous transduction of two vectors (pLVX-rtTA/pLVX-cDNA) harboring the reverse tetracycline transactivator (rtTA) and the tetracycline response element (TRE)-controlled transgene, respectively, leading to inefficient cardiac reprogramming. Herein, we developed a single-construct-based polycistronic Dox-inducible vector (pDox-cDNA) expressing both the rtTA and TRE-controlled transgenes. Fluorescence activated cell sorting (FACS) analyses, quantitative RT-PCR, and immunostaining revealed that pDox-GMT increased cardiac reprogramming three-fold compared to the conventional pLVX-rtTA/pLVX-GMT. After four weeks, pDox-GMT-induced iCMs expressed multiple cardiac genes, produced sarcomeric structures, and beat spontaneously. Co-transduction of pDox-Hand2 with retroviral pMX-GMT increased cardiac reprogramming three-fold compared to pMX-GMT alone. Temporal Dox administration revealed that Hand2 transgene expression is critical during the first two weeks of cardiac reprogramming. Microarray analyses demonstrated that Hand2 represses cell cycle-promoting genes and enhances cardiac reprogramming. Thus, we have developed an efficient temporal transgene expression system, which could be invaluable in the study of cardiac reprogramming.

Predictors of mammographic density: insights gained from a novel regression analysis of a twin study.

Understanding which factors influence mammographically dense and nondense areas is important because percent mammographic density adjusted for age is a strong, continuously distributed risk factor for breast cancer, especially when adjusted for weight or body mass index. Using computer-assisted methods, we measured mammographically dense areas for 571 monozygotic and 380 dizygotic Australian and North American twin pairs ages 40 to 70 years. We used a novel regression modeling approach in which each twin's measure of dense and nondense area was regressed against one or both of the twin's and co-twin's covariates. The nature of changes to regression estimates with the inclusion of the twin and/or co-twin's covariates can be evaluated for consistency with causal and/or other models. By causal, we mean that if it were possible to vary a covariate experimentally then the expected value of the outcome measure would change. After adjusting for the individual's weight, the co-twin associations with weight were attenuated, consistent with a causal effect of weight on mammographic measures, which in absolute log cm(2)/kg was thrice stronger for nondense than dense area. After adjusting for weight, later age at menarche, and greater height were associated with greater dense and lesser nondense areas in a manner inconsistent with causality. The associations of dense and nondense areas with parity are consistent with a causal effect and/or within-person confounding. The associations between mammographic density measures and height are consistent with shared early life environmental factors that predispose to both height and percent mammographic density and possibly breast cancer risk.

Lentiviral vector design and imaging approaches to visualize the early stages of cellular reprogramming.

Induced pluripotent stem cells (iPSCs) can be derived from somatic cells by gene transfer of reprogramming transcription factors. Expression levels of these factors strongly influence the overall efficacy to form iPSC colonies, but additional contribution of stochastic cell-intrinsic factors has been proposed. Here, we present engineered color-coded lentiviral vectors in which codon-optimized reprogramming factors are co-expressed by a strong retroviral promoter that is rapidly silenced in iPSC, and imaged the conversion of fibroblasts to iPSC. We combined fluorescence microscopy with long-term single cell tracking, and used live-cell imaging to analyze the emergence and composition of early iPSC clusters. Applying our engineered lentiviral vectors, we demonstrate that vector silencing typically occurs prior to or simultaneously with the induction of an Oct4-EGFP pluripotency marker. Around 7 days post-transduction (pt), a subfraction of cells in clonal colonies expressed Oct4-EGFP and rapidly expanded. Cell tracking of single cell-derived iPSC colonies supported the concept that stochastic epigenetic changes are necessary for reprogramming. We also found that iPSC colonies may emerge as a genetic mosaic originating from different clusters. Improved vector design with continuous cell tracking thus creates a powerful system to explore the subtle dynamics of biological processes such as early reprogramming events.

Microarray analysis of pneumococcal gene expression during invasive disease.

Streptococcus pneumoniae is a leading cause of invasive bacterial disease. This is the first study to examine the expression of S. pneumoniae genes in vivo by using whole-genome microarrays available from The Institute for Genomic Research. Total RNA was collected from pneumococci isolated from infected blood, infected cerebrospinal fluid, and bacteria attached to a pharyngeal epithelial cell line in vitro. Microarray analysis of pneumococcal genes expressed in these models identified body site-specific patterns of expression for virulence factors, transporters, transcription factors, translation-associated proteins, metabolism, and genes with unknown function. Contributions to virulence predicted for several unknown genes with enhanced expression in vivo were confirmed by insertion duplication mutagenesis and challenge of mice with the mutants. Finally, we cross-referenced our results with previous studies that used signature-tagged mutagenesis and differential fluorescence induction to identify genes that are potentially required by a broad range of pneumococcal strains for invasive disease.

Subpopulations of mouse blood monocytes differ in maturation stage and inflammatory response.

Blood monocytes are well-characterized precursors for macrophages and dendritic cells. Subsets of human monocytes with differential representation in various disease states are well known. In contrast, mouse monocyte subsets have been characterized minimally. In this study we identify three subpopulations of mouse monocytes that can be distinguished by differential expression of Ly-6C, CD43, CD11c, MBR, and CD62L. The subsets share the characteristics of extensive phagocytosis, similar expression of M-CSF receptor (CD115), and development into macrophages upon M-CSF stimulation. By eliminating blood monocytes with dichloromethylene-bisphosphonate-loaded liposomes and monitoring their repopulation, we showed a developmental relationship between the subsets. Monocytes were maximally depleted 18 h after liposome application and subsequently reappeared in the circulation. These cells were exclusively of the Ly-6C(high) subset, resembling bone marrow monocytes. Serial flow cytometric analyses of newly released Ly-6C(high) monocytes showed that Ly-6C expression on these cells was down-regulated while in circulation. Under inflammatory conditions elicited either by acute infection with Listeria monocytogenes or chronic infection with Leishmania major, there was a significant increase in immature Ly-6C(high) monocytes, resembling the inflammatory left shift of granulocytes. In addition, acute peritoneal inflammation recruited preferentially Ly-6C(med-high) monocytes. Taken together, these data identify distinct subpopulations of mouse blood monocytes that differ in maturation stage and capacity to become recruited to inflammatory sites.

A map of human circular RNAs in clinically relevant tissues

Cellular circular RNAs (circRNAs) are generated by head-to-tail splicing and are present in all multicellular organisms studied so far. Recently, circRNAs have emerged as a large class of RNA which can function as post-transcriptional regulators. It has also been shown that many circRNAs are tissue- and stage-specifically expressed. Moreover, the unusual stability and expression specificity make circRNAs important candidates for clinical biomarker research. Here, we present a circRNA expression resource of 20 human tissues highly relevant to disease-related research: vascular smooth muscle cells (VSMCs), human umbilical vein cells (HUVECs), artery endothelial cells (HUAECs), atrium, vena cava, neutrophils, platelets, cerebral cortex, placenta, and samples from mesenchymal stem cell differentiation. In eight different samples from a single donor, we found highly tissue-specific circRNA expression. Circular-to-linear RNA ratios revealed that many circRNAs were expressed higher than their linear host transcripts. Among the 71 validated circRNAs, we noticed potential biomarkers. In adenosine deaminase-deficient, severe combined immunodeficiency (ADA-SCID) patients and in Wiskott-Aldrich-Syndrome (WAS) patients' samples, we found evidence for differential circRNA expression of genes that are involved in the molecular pathogenesis of both phenotypes. Our findings underscore the need to assess circRNAs in mechanisms of human disease. KEY MESSAGES: circRNA resource catalog of 20 clinically relevant tissues. circRNA expression is highly tissue-specific. circRNA transcripts are often more abundant than their linear host RNAs. circRNAs can be differentially expressed in disease-associated genes.

Streptococcus pneumoniae early response genes to human lung epithelial cells

BACKGROUND Streptococcus pneumoniae infection starts from colonization of the host respiratory tract where interaction with host respiratory tract epithelial cells occurs. To investigate pneumococcal genes that are involved in the early stage of interaction with host epithelial cells, transcriptional responses of an encapsulated pathogenic pneumococcal strain TIGR4 upon exposure to human lung epithelial cells A549 for 0.5 h and 1 h time periods were investigated by using TIGR (JCVI) microarray technology. Gene expression changes were validated by quantitative real-time PCR (qRT-PCR) analysis. FINDINGS We observed different transcriptional profiles at two incubation time periods in which most gene expressions were down-regulated at 0.5 h but up-regulated at 1 h. Many genes associated with ribonucleotide biosynthesis were down-regulated at both time points, whereas the genes associated with cell envelope, energy metabolism, transport and protein synthesis were mostly up-regulated at 1 h. Furthermore, these profiles were compared to the transcriptomes of a TIGR4-derived strain in response to human macrophages for the same time periods. We found one set of genes that exhibited similar expression changes upon exposure to both types of host cells, including cell envelope-associated bgaA (SP0648) and nanA (SP1693), and uncharacterized gene clusters such as SP1677-SP1680 and SP1688-SP1690. CONCLUSION These data indicate that at the early stage of interaction with host epithelial cells, a complex gene regulation and expression change occur in bacteria. Some of them might play an essential role during pathogen-host interactions and for the establishment of infection.

Nef induces multiple genes involved in cholesterol synthesis and uptake in human immunodeficiency virus type 1-infected T cells.

Several recent reports indicate that cholesterol might play an important role in human immunodeficiency virus type 1 (HIV-1) replication. We investigated the effects of HIV-1 infection on cholesterol biosynthesis and uptake using microarrays. HIV-1 increased gene expression of cholesterol genes in both transformed T-cell lines and primary CD4(+) T cells. Consistent with our microarray data, (14)C-labeled mevalonate and acetate incorporation was increased in HIV-1-infected cells. Our data also demonstrate that changes in cholesterol biosynthesis and uptake are only observed in the presence of functional Nef, suggesting that increased cholesterol synthesis may contribute to Nef-mediated enhancement of virion infectivity and viral replication.

Smooth muscle lineage diversity in the chick embryo. Two types of aortic smooth muscle cell differ in growth and receptor-mediated transcriptional responses to transforming growth factor-beta.

Lineage analysis studies in the avian embryo have identified two types of smooth muscle cells (SMCs) in the tunica media of large elastic arteries; one that originates within the cardiac neural crest and is ectoderm in origin (Ect) and another that arises from local mesenchyme of mesodermal origin (Mes). To determine if differences in primary embryonic lineage can give rise to SMCs with stable differences in growth and differentiation properties, we isolated Ect and Mes SMCs from the Day 14 chick embryo aorta. We report that despite different primary embryonic origins, Ect and Mes SMCs express nearly identical levels of seven SMC differentiation markers in vitro, consistent with their common smooth muscle developmental fates in vivo. By contrast, Ect SMCs displayed a greater capacity for growth in serum-free medium than Mes SMCs, but only under conditions permitting short-range cell-cell interactions. Most of the peptide growth factors tested that might account for serum-independent growth (PDGF-AA, PDGF-BB, basic FGF, EGF, or activin) stimulated DNA synthesis to similar extents in Ect and Mes SMCs. However, we found dramatic, lineage-dependent differences in SMC responses to transforming growth factor-beta (TGF-beta). Exposure to TGF-beta 1 (0.4 to 400 pmole/liter) consistently increased DNA synthesis in Ect SMCs, whereas in paired cultures of Mes SMCs, TGF-beta 1 was growth inhibitory. In SMC cultures transfected with p3TP-lux, a luciferase reporter controlled by the TGF-beta 1-response elements of the human PAI-1 promoter, TGF-beta 1 (120 pM) produced 12 +/- 2-fold increases in luciferase activity in Ect SMCs and only 3 +/- 1.5-fold increases in Mes SMCs. Analysis of TGF-beta receptor phenotypes by Northern blot, radioligand binding, and crosslinking assays showed that Ect and Mes SMCs expressed similar levels of types I, II, and III TGF-beta receptors. However, using a polyclonal antibody specific for the chick type II TGF-beta receptor subunit, we demonstrate that Mes SMCs produce a fully glycosylated form of this protein while Ect SMCs elaborate only an unglycosylated type II TGF-beta receptor. These results show that Ect and Mes SMCs exhibit lineage-dependent differences in growth and receptor-mediated transcriptional responses to at least one important class of SMC morphogens and growth modifiers, e.g., the TGF-betas. Our findings suggest that different SMC populations within a common vessel wall may respond in lineage-dependent ways to signals that direct formation of the tunica media in the embryo and to factors involved in the progression of vascular disease later in life.

Identification of Senescent Cells in the Bone Microenvironment.

Cellular senescence is a fundamental mechanism by which cells remain metabolically active yet cease dividing and undergo distinct phenotypic alterations, including upregulation of p16Ink4a , profound secretome changes, telomere shortening, and decondensation of pericentromeric satellite DNA. Because senescent cells accumulate in multiple tissues with aging, these cells and the dysfunctional factors they secrete, termed the senescence-associated secretory phenotype (SASP), are increasingly recognized as promising therapeutic targets to prevent age-related degenerative pathologies, including osteoporosis. However, the cell type(s) within the bone microenvironment that undergoes senescence with aging in vivo has remained poorly understood, largely because previous studies have focused on senescence in cultured cells. Thus in young (age 6 months) and old (age 24 months) mice, we measured senescence and SASP markers in vivo in highly enriched cell populations, all rapidly isolated from bone/marrow without in vitro culture. In both females and males, p16Ink4a expression by real-time quantitative polymerase chain reaction (rt-qPCR) was significantly higher with aging in B cells, T cells, myeloid cells, osteoblast progenitors, osteoblasts, and osteocytes. Further, in vivo quantification of senescence-associated distension of satellites (SADS), ie, large-scale unraveling of pericentromeric satellite DNA, revealed significantly more senescent osteocytes in old compared with young bone cortices (11% versus 2%, p < 0.001). In addition, primary osteocytes from old mice had sixfold more (p < 0.001) telomere dysfunction-induced foci (TIFs) than osteocytes from young mice. Corresponding with the age-associated accumulation of senescent osteocytes was significantly higher expression of multiple SASP markers in osteocytes from old versus young mice, several of which also showed dramatic age-associated upregulation in myeloid cells. These data show that with aging, a subset of cells of various lineages within the bone microenvironment become senescent, although senescent myeloid cells and senescent osteocytes predominantly develop the SASP. Given the critical roles of osteocytes in orchestrating bone remodeling, our findings suggest that senescent osteocytes and their SASP may contribute to age-related bone loss. © 2016 American Society for Bone and Mineral Research.

Host–Parasite Interactions and the Evolution of Gene Expression

Interactions between hosts and parasites provide an ongoing source of selection that promotes the evolution of a variety of features in the interacting species. Here, we use a genetically explicit mathematical model to explore how patterns of gene expression evolve at genetic loci responsible for host resistance and parasite infection. Our results reveal the striking yet intuitive conclusion that gene expression should evolve along very different trajectories in the two interacting species. Specifically, host resistance loci should frequently evolve to co-express alleles, whereas parasite infection loci should evolve to express only a single allele. This result arises because hosts that co-express resistance alleles are able to recognize and clear a greater diversity of parasite genotypes. By the same token, parasites that co-express antigen or elicitor alleles are more likely to be recognized and cleared by the host, and this favours the expression of only a single allele. Our model provides testable predictions that can help interpret accumulating data on expression levels for genes relevant to host−parasite interactions.

Enterovirus-induced gene expression profile is critical for human pancreatic islet destruction

Virally induced inflammatory responses, beta cell destruction and release of beta cell autoantigens may lead to autoimmune reactions culminating in type 1 diabetes. Therefore, viral capability to induce beta cell death and the nature of virus-induced immune responses are among key determinants of diabetogenic viruses. We hypothesised that enterovirus infection induces a specific gene expression pattern that results in islet destruction and that such a host response pattern is not shared among all enterovirus infections but varies between virus strains. The changes in global gene expression and secreted cytokine profiles induced by lytic or benign enterovirus infections were studied in primary human pancreatic islet using DNA microarrays and viral strains either isolated at the clinical onset of type 1 diabetes or capable of causing a diabetes-like condition in mice. The expression of pro-inflammatory cytokine genes (IL-1-α, IL-1-β and TNF-α) that also mediate cytokine-induced beta cell dysfunction correlated with the lytic potential of a virus. Temporally increasing gene expression levels of double-stranded RNA recognition receptors, antiviral molecules, cytokines and chemokines were detected for all studied virus strains. Lytic coxsackievirus B5 (CBV-5)-DS infection also downregulated genes involved in glycolysis and insulin secretion. The results suggest a distinct, virus-strain-specific, gene expression pattern leading to pancreatic islet destruction and pro-inflammatory effects after enterovirus infection. However, neither viral replication nor cytotoxic cytokine production alone are sufficient to induce necrotic cell death. More likely the combined effect of these and possibly cellular energy depletion lie behind the enterovirus-induced necrosis of islets.

NDRG1 promotes growth of hepatocellular carcinoma cells by directly interacting with GSK-3β and Nur77 to prevent β-catenin degradation

The N-myc downstream regulated gene 1 (NDRG1) is significantly associated with advanced tumor stages and poor survival of hepatocellular carcinoma (HCC), thereby implicating it as a potential target for HCC treatment. We aim to further understand its biological roles in hepatocarcinogenesis, as a means to exploit it for therapeutic purposes. By screening using the ProtoArray® Human Protein Microarrays, we identified glycogen synthase kinase 3β (GSK-3β) and the orphan nuclear receptor (Nur77) as potential interaction partners of NDRG1. These interactions were confirmed in HCC cell lines in vitro by co-immunoprecipitation; and co-localizations of NDRG1 with GSK-3β and Nur77 were observed by immunofluorescence staining. Additionally, high levels of NDRG1 competitively bind to GSK-3β and Nur77 to allow β-catenin to escape degradation, with consequent elevated levels of downstream oncogenic genes. In vivo, we consistently observed that NDRG1 suppression in HCC xenografts decreased β-catenin levels and its downstream target Cyclin D1, with concomitant tumor growth inhibition. Clinically, the over-expression of NDRG1 in HCC patient samples is positively correlated with GSK-3β-9ser (| R | = 0.28, p = 0.01), Nur77 (| R | = 0.42, p < 0.001), and β-catenin (| R |= 0.32, p = 0.003) expressions. In conclusion, we identified GSK-3β and Nur77 as novel interaction partners of NDRG1. These protein-protein interactions regulate the turnover of β-catenin and subsequent downstream signaling mediated by β-catenin in HCC cells, and provides potential targets for future therapeutic interventions.

Uncoupling of leading- and lagging-strand DNA replication during lesion bypass in vivo.

Numerous agents attack DNA, forming lesions that impair normal replication. Specialized DNA polymerases transiently replace the replicative polymerase and copy past lesions, thus generating mutations, the major initiating cause of cancer. We monitored, in Escherichia coli, the kinetics of replication of both strands of DNA molecules containing a single replication block in either the leading or lagging strand. Despite a block in the leading strand, lagging-strand synthesis proceeded further, implying transient uncoupling of concurrent strand synthesis. Replication through the lesion requires specialized DNA polymerases and is achieved with similar kinetics and efficiencies in both strands.