Shared and Compartment-Specific Processes in Nucleus Pulposus and Annulus Fibrosus During Intervertebral Disc Degeneration

Elucidating how cell populations promote onset and progression of intervertebral disc degeneration (IDD) has the potential to enable more precise therapeutic targeting of cells and mechanisms. Single-cell RNA-sequencing (scRNA-seq) is performed on surgically separated annulus fibrosus (AF) (19,978; 26,983 cells) and nucleus pulposus (NP) (20,884; 24,489 cells) from healthy and diseased human intervertebral discs (IVD). In both tissue types, depletion of cell subsets involved in maintenance of healthy IVD is observed, specifically the immature cell subsets - fibroblast progenitors and stem cells - indicative of an impairment of normal tissue self-renewal. Tissue-specific changes are also identified. In NP, several fibrotic populations are increased in degenerated IVD, indicating tissue-remodeling. In degenerated AF, a novel disease-associated subset is identified, which expresses disease-promoting genes. It is associated with pathogenic biological processes and the main gene regulatory networks include thrombospondin signaling and FOXO1 transcription factor. In NP and AF cells thrombospondin protein promoted expression of genes associated with TGFβ/fibrosis signaling, angiogenesis, and nervous system development. The data reveal new insights of both shared and tissue-specific changes in specific cell populations in AF and NP during IVD degeneration. These identified mechanisms and molecules are novel and more precise targets for IDD prevention and treatment.

Human-specific duplicate CHRFAM7A gene is associated with more severe osteoarthritis and amplifies pain behaviours

OBJECTIVES

CHRFAM7A is a uniquely human fusion gene that functions as a dominant negative regulator of alpha 7 acetylcholine nicotinic receptor (α7nAChR) in vitro. This study determined the impact of CHRFAM7A on α7nAChR agonist responses, osteoarthritis (OA) severity and pain behaviours and investigated mechanisms.

METHODS

Transgenic CHRFAM7A (TgCHRFAM7A) mice were used to determine the impact of CHRFAM7A on knee OA histology, pain severity in OA and other pain models, response to nAchR agonist and IL-1β. Mouse and human cells were used for mechanistic studies.

RESULTS

Transgenic (Tg) TgCHRFAM7A mice developed more severe structural damage and increased mechanical allodynia than wild type (WT) mice in the destabilisation of medial meniscus model of OA. This was associated with a decreased suppression of inflammation by α7nAchR agonist. TgCHRFAM7A mice displayed a higher basal sensitivity to pain stimuli and increased pain behaviour in the monoiodoacetate and formalin models. Dorsal root ganglia of TgCHRFAM7A mice showed increased macrophage infiltration and expression of the chemokine fractalkine and also had a compromised antinociceptive response to the α7nAchR agonist nicotine. Both native CHRNA7 and CHRFAM7A subunits were expressed in human joint tissues and the CHRFAM7A/CHRNA7 ratio was increased in OA cartilage. Human chondrocytes with two copies of CHRFAM7A had reduced anti-inflammatory responses to nicotine.

CONCLUSION

CHRFAM7A is an aggravating factor for OA-associated inflammation and tissue damage and a novel genetic risk factor and therapeutic target for pain.

Senescent cell population with ZEB1 transcription factor as its main regulator promotes osteoarthritis in cartilage and meniscus

OBJECTIVES

Single-cell level analysis of articular cartilage and meniscus tissues from human healthy and osteoarthritis (OA) knees.

METHODS

Single-cell RNA sequencing (scRNA-seq) analyses were performed on articular cartilage and meniscus tissues from healthy (n=6, n=7) and OA (n=6, n=6) knees. Expression of genes of interest was validated using immunohistochemistry and RNA-seq and function was analysed by gene overexpression and depletion.

RESULTS

scRNA-seq analyses of human knee articular cartilage (70 972 cells) and meniscus (78 017 cells) identified a pathogenic subset that is shared between both tissues. This cell population is expanded in OA and has strong OA and senescence gene signatures. Further, this subset has critical roles in extracellular matrix (ECM) and tenascin signalling and is the dominant sender of signals to all other cartilage and meniscus clusters and a receiver of TGFβ signalling. Fibroblast activating protein (FAP) is also a dysregulated gene in this cluster and promotes ECM degradation. Regulons that are controlled by transcription factor ZEB1 are shared between the pathogenic subset in articular cartilage and meniscus. In meniscus and cartilage cells, FAP and ZEB1 promote expression of genes that contribute to OA pathogenesis, including senescence.

CONCLUSIONS

These single-cell studies identified a senescent pathogenic cell cluster that is present in cartilage and meniscus and has FAP and ZEB1 as main regulators which are novel and promising therapeutic targets for OA-associated pathways in both tissues.

Transcriptomic analyses of joint tissues during osteoarthritis development in a rat model reveal dysregulated mechanotransduction and extracellular matrix pathways

OBJECTIVE

Transcriptomic changes in joint tissues during the development of osteoarthritis (OA) are of interest for the discovery of biomarkers and mechanisms of disease. The objective of this study was to use the rat medial meniscus transection (MMT) model to discover stage and tissue-specific transcriptomic changes.

DESIGN

Sham or MMT surgeries were performed in mature rats. Cartilage, menisci and synovium were scored for histopathological changes at 2, 4 and 6 weeks post-surgery and processed for RNA-sequencing. Differentially expressed genes (DEG) were used to identify pathways and mechanisms. Published transcriptomic datasets from animal models and human OA were used to confirm and extend present findings.

RESULTS

The total number of DEGs was already high at 2 weeks (723 in meniscus), followed by cartilage (259) and synovium (42) and declined to varying degrees in meniscus and synovium but increased in cartilage at 6 weeks. The most upregulated genes included tenascins. The 'response to mechanical stimulus' and extracellular matrix-related pathways were enriched in both cartilage and meniscus. Pathways that were enriched in synovium at 4 weeks indicate processes related to synovial hyperplasia and fibrosis. Synovium also showed upregulation of IL-11 and several MMPs. The mechanical stimulus pathway included upregulation of the mechanoreceptors PIEZO1, PIEZO2 and TRPV4 and nerve growth factor. Analysis of data from prior RNA-sequencing studies of animal models and human OA support these findings.

CONCLUSION

These results indicate several shared pathways that are affected during OA in cartilage and meniscus and support the role of mechanotransduction and other pathways in OA pathogenesis.

Low-Cost Peptide Microarrays for Mapping Continuous Antibody Epitopes

Understanding antibody specificity and defining response profiles to antigens continue to be essential to both vaccine research and therapeutic antibody development. Peptide scanning assays enable mapping of continuous epitopes in order to delineate antibody-antigen interactions beyond traditional immunoassay formats. We have developed a relatively low-cost method to generate peptide microarray slides for antibody binding studies that allow for interrogation of up to 1536 overlapping peptides derived from the target antigens on a single microslide. Using an IntavisAG MultiPep RS peptide synthesizer and a Digilab MicroGrid II 600 microarray printer robot, each peptide is tagged with a polyethylene glycol aminooxy terminus to improve peptide solubility, orientation, and conjugation efficiency to the slide surface. Interrogation of the surface can then be performed using polyclonal immune sera or monoclonal antibodies, and sensitive detection using an InnoScan 1100 AL scanner with fluorescent-conjugated secondary reagents maximizes conservation of reagents.

Krüppel-like factor-4 and Krüppel-like factor-2 are important regulators of joint tissue cells and protect against tissue destruction and inflammation in osteoarthritis

OBJECTIVES

Analysing expression patterns of Krüppel-like factor (KLF) transcription factors in normal and osteoarthritis (OA) human cartilage, and determining functions and mechanisms of KLF4 and KLF2 in joint homoeostasis and OA pathogenesis.

METHODS

Experimental approaches included human joint tissues cells, transgenic mice and mouse OA model with viral KLF4 gene delivery to demonstrate therapeutic benefit in structure and pain improvement. Mechanistic studies applied global gene expression analysis and chromatin immunoprecipitation sequencing (ChIP-seq).

RESULTS

Several KLF genes were significantly decreased in OA cartilage. Among them, KLF4 and KLF2 were strong inducers of cartilage collagen genes and Proteoglycan-4. Cartilage-specific deletion of Klf2 in mature mice aggravated severity of experimental OA. Transduction of human chondrocytes with Adenovirus (Ad) expressing KLF4 or KLF2 enhanced expression of major cartilage extracellular matrix (ECM) genes and SRY-box transcription factor-9, and suppressed mediators of inflammation and ECM-degrading enzymes. Ad-KLF4 and Ad-KLF2 enhanced similar protective functions in meniscus cells and synoviocytes, and promoted chondrocytic differentiation of human mesenchymal stem cells. Viral KLF4 delivery into mouse knees reduced severity of OA-associated changes in cartilage, meniscus and synovium, and improved pain behaviours. ChIP-seq analysis suggested that KLF4 directly bound cartilage signature genes. Ras-related protein-1 signalling was the most enriched pathway in KLF4-transduced cells, and its signalling axis was involved in upregulating cartilage ECM genes by KLF4 and KLF2.

CONCLUSIONS

KLF4 and KLF2 may be central transcription factors that increase protective and regenerative functions in joint tissue cells, suggesting that KLF gene transfer or molecules upregulating KLFs are therapeutic candidates for OA.

Real-time digital polymerase chain reaction (PCR) as a novel technology improves limit of detection for rare allele assays

Background

Tumor heterogeneity may lead to false negative test results for tissue biopsy-based companion diagnostic tests. Real-time polymerase chain reaction (PCR) and digital PCR assays are used to detect rare alleles in cell-free circulating DNA for liquid biopsies; however, those tests lack strong sensitivity at low allele frequencies. We show here a novel real-time digital PCR instrument that utilizes cycle-based amplification curves to further improve the sensitivity and quantification accuracy of digital PCR.

Methods

The novel real-time digital PCR instrument was compared to an endpoint digital PCR system to determine the sensitivity and quantification accuracy of both instruments. Samples were all thermal cycled on the real-time digital PCR instrument but were analyzed on both endpoint and real-time digital PCR instruments to compare the performance without introducing other variables. Contrived samples for epidermal growth factor receptor (EGFR) exon 19 deletion, T790M, and L858R point mutations as well as human epidermal growth factor receptor 2 (HER2) amplification were tested. Different mutant allele frequencies and wildtype to mutant gene copy number ratios were tested for EGFR and HER2, respectively.

Results

By removing false positive datapoints using real-time amplification curves, real-time digital PCR improved sensitivity by lowering the baseline for wildtype samples. For EGFR 19del assay, samples with 2 or more fluorescein amidite (FAM) labeled positive wells are determined positive by real-time digital PCR, while a minimum of 5 FAM positive datapoints is needed by endpoint digital PCR. Improved limit of detection for EGFR 19del mutation was also observed. Real-time digital PCR also had better quantification accuracy and sensitivity, resulting in the mutant allele frequencies being closer to the expected values for all EGFR mutations, especially at very low allele frequencies. However, at high allele frequencies or for gene amplification assays, real-time digital PCR is comparable with endpoint digital PCR.

Conclusions

This novel technology with improved sensitivity is important and needed because it addresses current issues with liquid biopsy tests. Due to limited amounts of circulating tumor DNA (ctDNA) obtained for liquid biopsy tests, few copies of mutant alleles are expected. With the lower baseline of real-time digital PCR, false negative test results from tissue biopsy would be more effectively reduced, leading to more patients receiving the targeted therapy they need for better survival.

Extended survival versus accelerated rejection of nonhuman primate islet allografts: Effect of mesenchymal stem cell source and timing

Mesenchymal stem cells (MSC) have been shown to be immunomodulatory, tissue regenerative, and graft promoting; however, several questions remain with regard to ideal MSC source and timing of administration. In this study, we utilized a rigorous preclinical model of allogeneic islet cell transplantation, incorporating reduced immune suppression and near to complete mismatch of major histocompatibility antigens between the diabetic cynomolgus monkey recipient and the islet donor, to evaluate both the graft promoting impact of MSC source, that is, derived from the islet recipient, the islet donor or an unrelated third party as well as the impact of timing. Co-transplant of MSC and islets on post-operative day 0, followed by additional IV MSC infusions in the first posttransplant month, resulted in prolongation of rejection free and overall islet survival and superior metabolic control for animals treated with recipient as compared to donor or third-party MSC. Immunological analyses demonstrated that infusion of MSC from either source did not prevent alloantibody formation to the islet or MSC donor; however, treatment with recipient MSC resulted in significant downregulation of memory T cells, decreased anti-donor T cell proliferation, and a trend toward increased Tregulatory:Tconventional ratios.

Increased sensitivity using real-time dPCR for detection of SARS-CoV-2

A real-time dPCR system was developed to improve the sensitivity, specificity and quantification accuracy of end point dPCR. We compared three technologies - real-time qPCR, end point dPCR and real-time dPCR - in the context of SARS-CoV-2. Some improvement in limit of detection was obtained with end point dPCR compared with real-time qPCR, and the limit of detection was further improved with the newly developed real-time dPCR technology through removal of false-positive signals. Real-time dPCR showed increased linear dynamic range compared with end point dPCR based on quantitation from amplification curves. Real-time dPCR can improve the performance of TaqMan assays beyond real-time qPCR and end point dPCR with better sensitivity and specificity, absolute quantification and a wider linear range of detection.

Feasibility and Comparison Study of Fecal Sample Collection Methods in Healthy Volunteers and Solid Organ Transplant Recipients Using 16S rRNA and Metagenomics Approaches

The human microbiome encompasses a variety of microorganisms that change dynamically and are in close contact with the body. The microbiome influences health and homeostasis, as well as the immune system, and any significant change in this equilibrium (dysbiosis) triggers both acute and chronic health conditions. Microbiome research has surged, in part, due to advanced sequencing technologies enabling rapid, accurate, and cost-effective identification of the microbiome. A major prerequisite for stool sample collection to study the gut microbiome in longitudinal prospective studies requires standardized protocols that can be easily replicated. However, there are still significant bottlenecks to stool specimen collection that contribute to low patient retention rates in microbiome studies. These barriers are further exacerbated in solid organ transplant recipients where diarrhea is estimated to occur in up to half the patient population. We sought to test two relatively easy sample collection methods (fecal swab and wipes) and compare them to the more cumbersome "gold" standard collection method (scoop) using two different sequencing technologies (16S ribosomal RNA sequencing and shotgun metagenomics). Our comparison of the collection methods shows that both the swabs and the wipes are comparable to the scoop method in terms of bacterial abundance and diversity. The swabs, however, were closer in representation to the scoop and were easier to collect and process compared to the wipes. Potential contamination of the swab and the wipe samples by abundant skin commensals was low in our analysis. Comparison of the two sequencing technologies showed that they were complementary, and that 16S sequencing provided enough coverage to detect and differentiate between bacterial species identified in the collected samples. Our pilot study demonstrates that alternative collection methods for stool sampling are a viable option in clinical applications, such as organ transplant studies. The use of these methods may result in better patient retention recruitment rates in serial microbiome studies.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Seasonal dynamics of DNA and RNA viral bioaerosol communities in a daycare center

BACKGROUND

Viruses play an important role in ecosystems, including the built environment (BE). While numerous studies have characterized bacterial and fungal microbiomes in the BE, few have focused on the viral microbiome (virome). Longitudinal microbiome studies provide insight into the stability and dynamics of microbial communities; however, few such studies exist for the microbiome of the BE, and most have focused on bacteria. Here, we present a longitudinal, metagenomic-based analysis of the airborne DNA and RNA virome of a children's daycare center. Specifically, we investigate how the airborne virome varies as a function of season and human occupancy, and we identify possible sources of the viruses and their hosts, mainly humans, animals, plants, and insects.

RESULTS

Season strongly influenced the airborne viral community composition, and a single sample collected when the daycare center was unoccupied suggested that occupancy also influenced the community. The pattern of influence differed between DNA and RNA viromes. Human-associated viruses were much more diverse and dominant in the winter, while the summertime virome contained a high relative proportion and diversity of plant-associated viruses.

CONCLUSIONS

This airborne microbiome in this building exhibited seasonality in its viral community but not its bacterial community. Human occupancy influenced both types of communities. By adding new data about the viral microbiome to complement burgeoning information about the bacterial and fungal microbiomes, this study contributes to a more complete understanding of the airborne microbiome.

Loss of GCNT2/I-branched glycans enhances melanoma growth and survival

Cancer cells often display altered cell-surface glycans compared to their nontransformed counterparts. However, functional contributions of glycans to cancer initiation and progression remain poorly understood. Here, from expression-based analyses across cancer lineages, we found that melanomas exhibit significant transcriptional changes in glycosylation-related genes. This gene signature revealed that, compared to normal melanocytes, melanomas downregulate I-branching glycosyltransferase, GCNT2, leading to a loss of cell-surface I-branched glycans. We found that GCNT2 inversely correlated with clinical progression and that loss of GCNT2 increased melanoma xenograft growth, promoted colony formation, and enhanced cell survival. Conversely, overexpression of GCNT2 decreased melanoma xenograft growth, inhibited colony formation, and increased cell death. More focused analyses revealed reduced signaling responses of two representative glycoprotein families modified by GCNT2, insulin-like growth factor receptor and integrins. Overall, these studies reveal how subtle changes in glycan structure can regulate several malignancy-associated pathways and alter melanoma signaling, growth, and survival.

A genome-wide survey of mutations in the Jurkat cell line

Background

The Jurkat cell line has an extensive history as a model of T cell signaling. But at the turn of the 21st century, some expression irregularities were observed, raising doubts about how closely the cell line paralleled normal human T cells. While numerous expression deficiencies have been described in Jurkat, genetic explanations have only been provided for a handful of defects.

Results

Here, we report a comprehensive catolog of genomic variation in the Jurkat cell line based on whole-genome sequencing. With this list of all detectable, non-reference sequences, we prioritize potentially damaging mutations by mining public databases for functional effects. We confirm documented mutations in Jurkat and propose links from detrimental gene variants to observed expression abnormalities in the cell line.

Conclusions

The Jurkat cell line harbors many mutations that are associated with cancer and contribute to Jurkat’s unique characteristics. Genes with damaging mutations in the Jurkat cell line are involved in T-cell receptor signaling (PTEN, INPP5D, CTLA4, and SYK), maintenance of genome stability (TP53, BAX, and MSH2), and O-linked glycosylation (C1GALT1C1). This work ties together decades of molecular experiments and serves as a resource that will streamline both the interpretation of past research and the design of future Jurkat studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4718-6) contains supplementary material, which is available to authorized users.

Microdroplet PCR for Highly Multiplexed Targeted Bisulfite Sequencing

Many methods exist for examining CpG DNA methylation. However, many of these are qualitative, laborious to apply to a large number of genes simultaneously, or are not easy to target to specific regions of interest. Microdroplet PCR-based bisulfite sequencing allows for quantitative single base resolution analysis of investigator selected regions of interest. Following bisulfite conversion of genomic DNA, targeted microdroplet PCR is conducted with custom primer libraries. Samples are then fragmented, concatenated, and sequenced by high-throughput sequencing. The most recent technology allows for this method to be conducted with as little as 250 ng of bisulfite-converted DNA. The primary advantage of this method is the ability to hand-select the targeted regions covered by up to 10,000 amplicons of 500-600 bp. Moreover, the nature of microdroplet PCR virtually eliminates PCR bias and allows for the amplification of all targets simultaneously in a single tube.

Allele-specific wild-type TP53 expression in the unaffected carrier parent of children with Li-Fraumeni syndrome

Li-Fraumeni syndrome (LFS) is an autosomal dominant disorder where an oncogenic TP53 germline mutation is passed from parent to child. Tumor protein p53 is a key tumor suppressor regulating cell cycle arrest in response to DNA damage. Paradoxically, some mutant TP53 carriers remain unaffected, while their children develop cancer within the first few years of life. To address this paradox, response to UV stress was compared in dermal fibroblasts (dFb) from an affected LFS patient vs. their unaffected carrier parent. UV induction of CDKN1A/p21, a regulatory target of p53, in LFS patient dFb was significantly reduced compared to the unaffected parent. UV exposure also induced significantly greater p53[Ser15]-phosphorylation in LFS patient dFb, a reported property of some mutant p53 variants. Taken together, these results suggested that unaffected parental dFb may express an increased proportion of wild-type vs. mutant p53. Indeed, a significantly increased ratio of wild-type to mutant TP53 allele-specific expression in the unaffected parent dFb was confirmed by RT-PCR-RFLP and RNA-seq analysis. Hence, allele-specific expression of wild-type TP53 may allow an unaffected parent to mount a response to genotoxic stress more characteristic of homozygous wild-type TP53 individuals than their affected offspring, providing protection from the oncogenesis associated with LFS.

Differential Sensitivity of Target Genes to Translational Repression by miR-17~92

MicroRNAs (miRNAs) are thought to exert their functions by modulating the expression of hundreds of target genes and each to a small degree, but it remains unclear how small changes in hundreds of target genes are translated into the specific function of a miRNA. Here, we conducted an integrated analysis of transcriptome and translatome of primary B cells from mutant mice expressing miR-17~92 at three different levels to address this issue. We found that target genes exhibit differential sensitivity to miRNA suppression and that only a small fraction of target genes are actually suppressed by a given concentration of miRNA under physiological conditions. Transgenic expression and deletion of the same miRNA gene regulate largely distinct sets of target genes. miR-17~92 controls target gene expression mainly through translational repression and 5’UTR plays an important role in regulating target gene sensitivity to miRNA suppression. These findings provide molecular insights into a model in which miRNAs exert their specific functions through a small number of key target genes.

MicroRNAs (miRNAs) are small RNAs encoded by our genome. Each miRNA binds hundreds of target mRNAs and performs specific functions. It is thought that miRNAs exert their function by reducing the expression of all these target genes and each to a small degree. However, these target genes often have very diverse functions. It has been unclear how small changes in hundreds of target genes with diverse functions are translated into the specific function of a miRNA. Here we take advantage of recent technical advances to globally examine the mRNA and protein levels of 868 target genes regulated by miR-17~92, the first oncogenic miRNA, in mutant mice with transgenic overexpression or deletion of this miRNA gene. We show that miR-17~92 regulates target gene expression mainly at the protein level, with little effect on mRNA. Surprisingly, only a small fraction of target genes respond to miR-17~92 expression changes. Further studies show that the sensitivity of target genes to miR-17~92 is determined by a non-coding region of target mRNA. Our findings demonstrate that not every target gene is equal, and suggest that the function of a miRNA is mediated by a small number of key target genes.

“I”-branched Carbohydrates Negatively Regulate Galectin-3-binding and Melanoma Malignancy

Malignant transformation is often associated with aberrant glycosylation of cell surface proteins. Commonly observed changes in glycan structures include increased presence of sialic acid, altered expression of mucins, and abnormal branching of N-glycans. Increasing evidence supports the idea that the presence of certain glycans correlate with cancer progression by affecting tumor cell invasiveness, ability to bind pro-tumorigenic galectins, and by promoting metastasis to distant organs. Taking this relationship into consideration, we conducted MALDI-TOF mass spectrometry on cell surface glycans as well as glycomic gene profiling from normal human epidermal melanocytes (NHEM) and human metastatic melanoma (MM) cells. We found that a defining feature between NHEMs and MMs was the relative level of “I”-branched and “i”-linear poly-N-acetyllactosamines on N-glycans. While NHEM predominantly expressed “I”-branched structures, MMs expressed “i”-linear modified N-glycans. Furthermore, glycomic gene profiling and confirmatory RT-qPCR data showed that the “I”-branching β1,6 N-acetylglucosaminyltransferase 2, GCNT2, is decreased in MM compared with NHEMs. With regards to its role in malignancy, we found that GCNT2 functioned as a negative regulator of galectin-3 (Gal-3) binding. GCNT2 overexpression in human MM cells inhibited Gal-3 mediated ERK1/2 phosphorylation and significantly reduced primary tumor growth. While GCNT2 knockdown in human MM cells promoted Gal-3 mediated ERK1/2 phosphorylation and significantly enhanced primary tumor growth. These findings highlight new glycobiological insights into melanoma development, implicating GCNT2 as a negative regulator of Gal-3 binding and melanoma malignancy.

The Activation-Induced Assembly of an RNA/Protein Interactome Centered on the Splicing Factor U2AF2 Regulates Gene Expression in Human CD4 T Cells

Activation of CD4 T cells is a reaction to challenges such as microbial pathogens, cancer and toxins that defines adaptive immune responses. The roles of T cell receptor crosslinking, intracellular signaling, and transcription factor activation are well described, but the importance of post-transcriptional regulation by RNA-binding proteins (RBPs) has not been considered in depth. We describe a new model expanding and activating primary human CD4 T cells and applied this to characterizing activation-induced assembly of splicing factors centered on U2AF2. We immunoprecipitated U2AF2 to identify what mRNA transcripts were bound as a function of activation by TCR crosslinking and costimulation. In parallel, mass spectrometry revealed the proteins incorporated into the U2AF2-centered RNA/protein interactome. Molecules that retained interaction with the U2AF2 complex after RNAse treatment were designated as "central" interactome members (CIMs). Mass spectrometry also identified a second class of activation-induced proteins, "peripheral" interactome members (PIMs), that bound to the same transcripts but were not in physical association with U2AF2 or its partners. siRNA knockdown of two CIMs and two PIMs caused changes in activation marker expression, cytokine secretion, and gene expression that were unique to each protein and mapped to pathways associated with key aspects of T cell activation. While knocking down the PIM, SYNCRIP, impacts a limited but immunologically important set of U2AF2-bound transcripts, knockdown of U2AF1 significantly impairs assembly of the majority of protein and mRNA components in the activation-induced interactome. These results demonstrated that CIMs and PIMs, either directly or indirectly through RNA, assembled into activation-induced U2AF2 complexes and play roles in post-transcriptional regulation of genes related to cytokine secretion. These data suggest an additional layer of regulation mediated by the activation-induced assembly of RNA splicing interactomes that is important for understanding T cell activation.

Transfection of microRNA Mimics Should Be Used with Caution

Transient transfection of chemically synthesized microRNA (miRNA) mimics is being used extensively to study the functions and mechanisms of endogenous miRNAs. However, it remains unclear whether transfected miRNAs behave similarly to endogenous miRNAs. Here we show that transient transfection of miRNA mimics into HeLa cells by a commonly used method led to the accumulation of high molecular weight RNA species and a few hundred fold increase in mature miRNA levels. In contrast, expression of the same miRNAs through lentiviral infection or plasmid transfection of HeLa cells, transgenic expression in primary lymphocytes, and endogenous overexpression in lymphoma and leukemia cell lines did not lead to the appearance of high molecular weight RNA species. The increase of mature miRNA levels in these cells was below 10-fold, which was sufficient to suppress target gene expression and to drive lymphoma development in mice. Moreover, transient transfection of miRNA mimics at high concentrations caused non-specific alterations in gene expression, while at low concentrations achieved expression levels comparable to other methods but failed to efficiently suppress target gene expression. Small RNA deep sequencing analysis revealed that the guide strands of miRNA mimics were frequently mutated, while unnatural passenger strands of some miRNA mimics accumulated to high levels. The high molecular weight RNA species were a heterogeneous mixture of several classes of RNA species generated by concatemerization, 5'- and 3'-end tailing of miRNA mimics. We speculate that the supraphysiological levels of mature miRNAs and these artifactual RNA species led to non-specific changes in gene expression. Our results have important implications for the design and interpretation of experiments primarily employing transient transfection of miRNA mimics.

ClickSeq: Fragmentation-Free Next-Generation Sequencing via Click Ligation of Adaptors to Stochastically Terminated 3'-Azido cDNAs

We present a simple method called "ClickSeq" for NGS (next-generation sequencing) library synthesis that uses click chemistry rather than enzymatic reactions for the ligation of Illumina sequencing adaptors. In ClickSeq, randomly primed reverse transcription reactions are supplemented with azido-2',3'-dideoxynucleotides that randomly terminate DNA synthesis and release 3'-azido-blocked cDNA fragments in a process akin to dideoxy-Sanger sequencing. Purified fragments are "click ligated" via copper-catalyzed alkyne-azide cycloaddition to DNA oligos modified with a 5'-alkyne group. This generates ssDNA molecules containing an unnatural triazole-linked DNA backbone that is sufficiently biocompatible for PCR amplification to generate a cDNA library for RNAseq. Here, we analyze viral RNAs and mRNA to demonstrate that ClickSeq produces unbiased NGS libraries with low error rates comparable to standard methods. Importantly, ClickSeq is robust against common artifacts of NGS such as chimera formation and artifactual recombination with fewer than 3 aberrant events detected per million reads.

RNA purification and expression analysis using microarrays and RNA deep sequencing

Transcriptome analysis or global gene expression profiling is a powerful tool for discovery as well as -understanding biological mechanisms in health and disease. We present in this chapter a description of methods used to isolate mRNA from cells and tissues that has been optimized for preservation of RNA quality using clinical materials and implemented successfully in several large, multicenter studies by the authors. In addition, two methods, gene expression microarrays and RNAseq, are described for mRNA profiling of cells and tissues from clinical or laboratory sources.

An evolutionarily conserved long noncoding RNA TUNA controls pluripotency and neural lineage commitment

Here, we generated a genome-scale shRNA library targeting long intergenic noncoding RNAs (lincRNAs) in the mouse. We performed an unbiased loss-of-function study in mouse embryonic stem cells (mESCs) and identified 20 lincRNAs involved in the maintenance of pluripotency. Among these, TUNA (Tcl1 Upstream Neuron-Associated lincRNA, or megamind) was required for pluripotency and formed a complex with three RNA-binding proteins (RBPs). The TUNA-RBP complex was detected at the promoters of Nanog, Sox2, and Fgf4, and knockdown of TUNA or the individual RBPs inhibited neural differentiation of mESCs. TUNA showed striking evolutionary conservation of both sequence- and CNS-restricted expression in vertebrates. Accordingly, knockdown of tuna in zebrafish caused impaired locomotor function, and TUNA expression in the brains of Huntington's disease patients was significantly associated with disease grade. Our results suggest that the lincRNA TUNA plays a vital role in pluripotency and neural differentiation of ESCs and is associated with neurological function of adult vertebrates.

Library construction for next-generation sequencing: overviews and challenges

High-throughput sequencing, also known as next-generation sequencing (NGS), has revolutionized genomic research. In recent years, NGS technology has steadily improved, with costs dropping and the number and range of sequencing applications increasing exponentially. Here, we examine the critical role of sequencing library quality and consider important challenges when preparing NGS libraries from DNA and RNA sources. Factors such as the quantity and physical characteristics of the RNA or DNA source material as well as the desired application (i.e., genome sequencing, targeted sequencing, RNA-seq, ChIP-seq, RIP-seq, and methylation) are addressed in the context of preparing high quality sequencing libraries. In addition, the current methods for preparing NGS libraries from single cells are also discussed.

Technical variations in low-input RNA-seq methodologies

Recent advances in RNA-seq methodologies from limiting amounts of mRNA have facilitated the characterization of rare cell-types in various biological systems. So far, however, technical variations in these methods have not been adequately characterized, vis-à-vis sensitivity, starting with reduced levels of mRNA. Here, we generated sequencing libraries from limiting amounts of mRNA using three amplification-based methods, viz. Smart-seq, DP-seq and CEL-seq, and demonstrated significant technical variations in these libraries. Reduction in mRNA levels led to inefficient amplification of the majority of low to moderately expressed transcripts. Furthermore, noise in primer hybridization and/or enzyme incorporation was magnified during the amplification step resulting in significant distortions in fold changes of the transcripts. Consequently, the majority of the differentially expressed transcripts identified were either high-expressed and/or exhibited high fold changes. High technical variations ultimately masked subtle biological differences mandating the development of improved amplification-based strategies for quantitative transcriptomics from limiting amounts of mRNA.

Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua

BACKGROUND

The medicinal plant Artemisia annua is covered with filamentous trichomes and glandular, artemisinin producing trichomes. A high artemisinin supply is needed at a reduced cost for treating malaria. Artemisinin production in bioreactors can be facilitated if a better insight is obtained in the biosynthesis of artemisinin and other metabolites. Therefore, metabolic activities of glandular and filamentous trichomes were investigated at the transcriptome level.

RESULTS

By laser pressure catapulting, glandular and filamentous trichomes as well as apical and sub-apical cells from glandular trichomes were collected and their transcriptome was sequenced using Illumina RNA-Seq. A de novo transcriptome was assembled (Trinity) and studied with a differential expression analysis (edgeR).A comparison of the transcriptome from glandular and filamentous trichomes shows that MEP, MVA, most terpene and lipid biosynthesis pathways are significantly upregulated in glandular trichomes. Conversely, some transcripts coding for specific sesquiterpenoid and triterpenoid enzymes such as 8-epi-cedrol synthase and an uncharacterized oxidosqualene cyclase were significantly upregulated in filamentous trichomes. All known artemisinin biosynthesis genes are upregulated in glandular trichomes and were detected in both the apical and sub-apical cells of the glandular trichomes. No significant differential expression could be observed between the apical and sub-apical cells.

CONCLUSIONS

Our results underscore the vast metabolic capacities of A. annua glandular trichomes but nonetheless point to the existence of specific terpene metabolic pathways in the filamentous trichomes. Candidate genes that might be involved in artemisinin biosynthesis are proposed based on their putative function and their differential expression level.

Optimization of peptide arrays for studying antibodies to hepatitis C virus continuous epitopes

Accurate and in-depth mapping of antibody responses is of great value in vaccine and antibody research. Using hepatitis C virus (HCV) as a model, we developed an affordable and high-throughput microarray-based assay for mapping antibody specificities to continuous antibody epitopes of HCV at high resolution. Important parameters in the chemistry for conjugating peptides/antigens to the array surface, the array layout, fluorophore choice and the methods for data analysis were investigated. Microscopic glass slide pre-coated with N-Hydroxysuccinimide (NHS)-ester (Slide H) was the preferred surface for conjugation of aminooxy-tagged peptides. This combination provides a simple chemical means to orient the peptides to the conjugation surface via an orthogonal covalent linkage at the N- or C-terminus of each peptide. The addition of polyvinyl alcohol to printing buffer gave uniform spot morphology and improved sensitivity and specificity of binding signals. Libraries of overlapping peptides covering the HCV E1 and E2 glycoprotein polypeptides (15-mer, 10 amino acids overlap) of 6 major HCV genotypes and the entire polypeptide sequence of the prototypic strain H77 were synthesized and printed in quadruplets in the assays. The utility of the peptide arrays was confirmed using HCV monoclonal antibodies (mAbs) specific to known continuous epitopes and immune sera of rabbits immunized with HCV antigens. The methods developed here can be easily adapted to studying antibody responses to antigens relevant in vaccine and autoimmune research.

VHL-dependent regulation of a β-dystroglycan glycoform and glycogene expression in renal cancer

Identification of novel biomarkers and targets in renal cell carcinoma (RCC) remains a priority and one cellular compartment that is a rich potential source of such molecules is the plasma membrane. A shotgun proteomic analysis of cell surface proteins enriched by cell surface biotinylation and avidin affinity chromatography was explored using the UMRC2- renal cancer cell line, which lacks von Hippel-Lindau (VHL) tumour suppressor gene function, to determine whether proteins of interest could be detected. Of the 814 proteins identified ~22% were plasma membrane or membrane-associated, including several with known associations with cancer. This included β-dystroglycan, the transmembrane subunit of the DAG1 gene product. VHL-dependent changes in the form of β-dystroglycan were detected in UMRC2-/+VHL transfectants. Deglycosylation experiments showed that this was due to differential sialylation. Analysis of normal kidney cortex and conventional RCC tissues showed that a similar change also occurred in vivo. Investigation of the expression of genes involved in glycosylation in UMRC2-/+VHL cells using a focussed microarray highlighted a number of enzymes involved in sialylation; upregulation of bifunctional UDP-N-acetylglucosamine 2-epimerase/N-acetylmannosamine kinase (GNE) was validated in UMRC2- cells compared with their +VHL counterparts and also found in conventional RCC tissue. These results implicate VHL in the regulation of glycosylation and raise interesting questions regarding the extent and importance of such changes in RCC.

Genomic meta-analysis of growth factor and integrin pathways in chronic kidney transplant injury

BACKGROUND

Chronic Allograft Nephropathy (CAN) is a clinical entity of progressive kidney transplant injury. The defining histology is tubular atrophy with interstitial fibrosis (IFTA). Using a meta-analysis of microarrays from 84 kidney transplant biopsies, we revealed growth factor and integrin adhesion molecule pathways differentially expressed and correlated with histological progression. A bioinformatics approach mining independent datasets leverages new and existing data to identify correlative changes in integrin and growth factor signaling pathways.

RESULTS

Analysis of CAN/IFTA Banff grades showed that hepatocyte growth factor (HGF), and epidermal growth factor (EGF) pathways are significantly differentially expressed in all classes of CAN/IFTA. MAPK-dependent pathways were also significant. However, the TGFβ pathways, albeit present, failed to differentiate CAN/IFTA progression. The integrin subunits β8, αv, αμ and β5 are differentially expressed, but β1, β6 and α6 specifically correlate with progression of chronic injury. Results were validated using our published proteomic profiling of CAN/IFTA.

CONCLUSIONS

CAN/IFTA with chronic kidney injury is characterized by expression of distinct growth factors and specific integrin adhesion molecules as well as their canonical signaling pathways. Drug target mapping suggests several novel candidates for the next generation of therapeutics to prevent or treat progressive transplant dysfunction with interstitial fibrosis.

Noncoding transcription within the Igh distal V(H) region at PAIR elements affects the 3D structure of the Igh locus in pro-B cells

Noncoding sense and antisense germ-line transcription within the Ig heavy chain locus precedes V(D)J recombination and has been proposed to be associated with Igh locus accessibility, although its precise role remains elusive. However, no global analysis of germ-line transcription throughout the Igh locus has been done. Therefore, we performed directional RNA-seq, demonstrating the locations and extent of both sense and antisense transcription throughout the Igh locus. Surprisingly, the majority of antisense transcripts are localized around two Pax5-activated intergenic repeat (PAIR) elements in the distal IghV region. Importantly, long-distance loops measured by chromosome conformation capture (3C) are observed between these two active PAIR promoters and Eμ, the start site of Iμ germ-line transcription, in a lineage- and stage-specific manner, even though this antisense transcription is Eμ-independent. YY1(-/-) pro-B cells are greatly impaired in distal V(H) gene rearrangement and Igh locus compaction, and we demonstrate that YY1 deficiency greatly reduces antisense transcription and PAIR-Eμ interactions. ChIP-seq shows high level YY1 binding only at Eμ, but low levels near some antisense promoters. PAIR-Eμ interactions are not disrupted by DRB, which blocks transcription elongation without disrupting transcription factories once they are established, but the looping is reduced after heat-shock treatment, which disrupts transcription factories. We propose that transcription-mediated interactions, most likely at transcription factories, initially compact the Igh locus, bringing distal V(H) genes close to the DJ(H) rearrangement which is adjacent to Eμ. Therefore, we hypothesize that one key role of noncoding germ-line transcription is to facilitate locus compaction, allowing distal V(H) genes to undergo efficient rearrangement.

Expression of glycogenes in differentiating human NT2N neurons. Downregulation of fucosyltransferase 9 leads to decreased Lewis(x) levels and impaired neurite outgrowth

BACKGROUND

Several glycan structures are functionally relevant in biological events associated with differentiation and regeneration which occur in the central nervous system. Here we have analysed the glycogene expression and glycosylation patterns during human NT2N neuron differentiation. We have further studied the impact of downregulating fucosyltransferase 9 (FUT9) on neurite outgrowth.

METHODS

The expression of glycogenes in human NT2N neurons differentiating from teratocarcinoma NTERA-2/cl.D1 cells has been analysed using the GlycoV4 GeneChip expression microarray. Changes in glycosylation have been monitored by immunoblot, immunofluorescence microscopy, HPLC and MALDI-TOF MS. Peptide mass fingerprinting and immunoprecipitation have been used for protein identification. FUT9 was downregulated using silencing RNA.

RESULTS AND CONCLUSIONS

One hundred twelve mRNA transcripts showed statistically significant up-regulation, including the genes coding for proteins involved in the synthesis of the Lewis(x) motif (FUT9), polysialic acid (ST8SIA2 and ST8SIA4) and HNK-1 (B3GAT2). Accordingly, increased levels of the corresponding carbohydrate epitopes have been observed. The Lewis(x) structure was found in a carrier glycoprotein that was identified as the CRA-a isoform of human neural cell adhesion molecule 1. Downregulation of FUT9 caused significant decreases in the levels of Lewis(x), as well as GAP-43, a marker of neurite outgrowth. Concomitantly, a reduction in neurite formation and outgrowth has been observed that was reversed by FUT9 overexpression.

GENERAL SIGNIFICANCE

These results provided information about the regulation of glycogenes during neuron differentiation and they showed that the Lewis(x) motif plays a functional role in neurite outgrowth from human neurons.

Glycophenotypic alterations induced by Pteridium aquilinum in mice gastric mucosa: synergistic effect with Helicobacter pylori infection

The bracken fern Pteridium aquilinum is a plant known to be carcinogenic to animals. Epidemiological studies have shown an association between bracken fern exposure and gastric cancer development in humans. The biological effects of exposure to this plant within the gastric carcinogenesis process are not fully understood. In the present work, effects in the gastric mucosa of mice treated with Pteridium aquilinum were evaluated, as well as molecular mechanisms underlying the synergistic role with Helicobacter pylori infection. Our results showed that exposure to Pteridium aquilinum induces histomorphological modifications including increased expression of acidic glycoconjugates in the gastric mucosa. The transcriptome analysis of gastric mucosa showed that upon exposure to Pteridium aquilinum several glycosyltransferase genes were differently expressed, including Galntl4, C1galt1 and St3gal2, that are mainly involved in the biosynthesis of simple mucin-type carbohydrate antigens. Concomitant treatment with Pteridium aquilinum and infection with Helicobacter pylori also resulted in differently expressed glycosyltransferase genes underlying the biosynthesis of terminal sialylated Lewis antigens, including Sialyl-Lewis^x. These results disclose the molecular basis for the altered pattern of glycan structures observed in the mice gastric mucosa. The gene transcription alterations and the induced glycophenotypic changes observed in the gastric mucosa contribute for the understanding of the molecular mechanisms underlying the role of Pteridium aquilinum in the gastric carcinogenesis process.

Illumina mate-paired DNA sequencing-library preparation using Cre-Lox recombination

Standard Illumina mate-paired libraries are constructed from 3- to 5-kb DNA fragments by a blunt-end circularization. Sequencing reads that pass through the junction of the two joined ends of a 3–5-kb DNA fragment are not easy to identify and pose problems during mapping and de novo assembly. Longer read lengths increase the possibility that a read will cross the junction. To solve this problem, we developed a mate-paired protocol for use with Illumina sequencing technology that uses Cre-Lox recombination instead of blunt end circularization. In this method, a LoxP sequence is incorporated at the junction site. This sequence allows screening reads for junctions without using a reference genome. Junction reads can be trimmed or split at the junction. Moreover, the location of the LoxP sequence in the reads distinguishes mate-paired reads from spurious paired-end reads. We tested this new method by preparing and sequencing a mate-paired library with an insert size of 3 kb from Saccharomyces cerevisiae. We present an analysis of the library quality statistics and a new bio-informatics tool called DeLoxer that can be used to analyze an IlluminaCre-Lox mate-paired data set. We also demonstrate how the resulting data significantly improves a de novo assembly of the S. cerevisiae genome.

Quantitative bias in Illumina TruSeq and a novel post amplification barcoding strategy for multiplexed DNA and small RNA deep sequencing

Here we demonstrate a method for unbiased multiplexed deep sequencing of RNA and DNA libraries using a novel, efficient and adaptable barcoding strategy called Post Amplification Ligation-Mediated (PALM). PALM barcoding is performed as the very last step of library preparation, eliminating a potential barcode-induced bias and allowing the flexibility to synthesize as many barcodes as needed. We sequenced PALM barcoded micro RNA (miRNA) and DNA reference samples and evaluated the quantitative barcode-induced bias in comparison to the same reference samples prepared using the Illumina TruSeq barcoding strategy. The Illumina TruSeq small RNA strategy introduces the barcode during the PCR step using differentially barcoded primers, while the TruSeq DNA strategy introduces the barcode before the PCR step by ligation of differentially barcoded adaptors. Results show virtually no bias between the differentially barcoded miRNA and DNA samples, both for the PALM and the TruSeq sample preparation methods. We also multiplexed miRNA reference samples using a pre-PCR barcode ligation. This barcoding strategy results in significant bias.

Genome-wide identification of Bcl11b gene targets reveals role in brain-derived neurotrophic factor signaling

B-cell leukemia/lymphoma 11B (Bcl11b) is a transcription factor showing predominant expression in the striatum. To date, there are no known gene targets of Bcl11b in the nervous system. Here, we define targets for Bcl11b in striatal cells by performing chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) in combination with genome-wide expression profiling. Transcriptome-wide analysis revealed that 694 genes were significantly altered in striatal cells over-expressing Bcl11b, including genes showing striatal-enriched expression similar to Bcl11b. ChIP-seq analysis demonstrated that Bcl11b bound a mixture of coding and non-coding sequences that were within 10 kb of the transcription start site of an annotated gene. Integrating all ChIP-seq hits with the microarray expression data, 248 direct targets of Bcl11b were identified. Functional analysis on the integrated gene target list identified several zinc-finger encoding genes as Bcl11b targets, and further revealed a significant association of Bcl11b to brain-derived neurotrophic factor/neurotrophin signaling. Analysis of ChIP-seq binding regions revealed significant consensus DNA binding motifs for Bcl11b. These data implicate Bcl11b as a novel regulator of the BDNF signaling pathway, which is disrupted in many neurological disorders. Specific targeting of the Bcl11b-DNA interaction could represent a novel therapeutic approach to lowering BDNF signaling specifically in striatal cells.

MicroRNA regulation of molecular networks mapped by global microRNA, mRNA, and protein expression in activated T lymphocytes

MicroRNAs (miRNAs) regulate specific immune mechanisms, but their genome-wide regulation of T lymphocyte activation is largely unknown. We performed a multidimensional functional genomics analysis to integrate genome-wide differential mRNA, miRNA, and protein expression as a function of human T lymphocyte activation and time. We surveyed expression of 420 human miRNAs in parallel with genome-wide mRNA expression. We identified a unique signature of 71 differentially expressed miRNAs, 57 of which were previously not known as regulators of immune activation. The majority of miRNAs are upregulated, mRNA expression of these target genes is downregulated, and this is a function of binding multiple miRNAs (combinatorial targeting). Our data reveal that consideration of this complex signature, rather than single miRNAs, is necessary to construct a full picture of miRNA-mediated regulation. Molecular network mapping of miRNA targets revealed the regulation of activation-induced immune signaling. In contrast, pathways populated by genes that are not miRNA targets are enriched for metabolism and biosynthesis. Finally, we specifically validated miR-155 (known) and miR-221 (novel in T lymphocytes) using locked nucleic acid inhibitors. Inhibition of these two highly upregulated miRNAs in CD4(+) T cells was shown to increase proliferation by removing suppression of four target genes linked to proliferation and survival. Thus, multiple lines of evidence link top functional networks directly to T lymphocyte immunity, underlining the value of mapping global gene, protein, and miRNA expression.

Application of microdroplet PCR for large-scale targeted bisulfite sequencing

Cytosine methylation of DNA CpG dinucleotides in gene promoters is an epigenetic modification that regulates gene transcription. While many methods exist to interrogate methylation states, few current methods offer large-scale, targeted, single CpG resolution. We report an approach combining bisulfite treatment followed by microdroplet PCR with next-generation sequencing to assay the methylation state of 50 genes in the regions 1 kb upstream of and downstream from their transcription start sites. This method yielded 96% coverage of the targeted CpGs and demonstrated high correlation between CpG island (CGI) DNA methylation and transcriptional regulation. The method was scaled to interrogate the methylation status of 77,674 CpGs in the promoter regions of 2100 genes in primary CD4 T cells. The 2100 gene library yielded 97% coverage of all targeted CpGs and 99% of the target amplicons.

Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation

Background

Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth.

Results

Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism.

Conclusions

Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation.

Identification and molecular analysis of glycosaminoglycans in cutaneous lupus erythematosus and dermatomyositis

Glycosaminoglycans (GAGs), also known histologically as dermal mucin, accumulate in several inflammatory skin conditions. Because different GAG species have distinct immunologic effects, the authors examined two GAGs, hyaluronan (HA) and chondroitin sulfate (CS), using specific stains in cutaneous lupus erythematosus (CLE) and dermatomyositis (DM). In the dermis of one CLE subtype, tumid LE (TLE), they found only increased HA, but both HA and CS were significantly elevated in another CLE subtype, discoid LE (DLE). DM lesional dermis accumulated mainly CS but not HA. The authors then used glycomic gene expression microarrays to assess the expression of HA- and CS-related genes in CLE skin. Real-time quantitative PCR confirmed significantly increased expression of HAS2, CHSY1, and C4ST1 in the combined groups of CLE lesions (n = 8) compared to healthy controls (n = 4). Thus, the increase in HA in CLE presumably results from upregulation of HAS2, whereas CHSY1 and C4ST1 appear to contribute to increased CS. Based on their known immunomodulatory effects in other systems, HA and CS may thus participate in the pathophysiology of these inflammatory skin conditions.

Rad3 decorates critical chromosomal domains with gammaH2A to protect genome integrity during S-Phase in fission yeast

Schizosaccharomyces pombe Rad3 checkpoint kinase and its human ortholog ATR are essential for maintaining genome integrity in cells treated with genotoxins that damage DNA or arrest replication forks. Rad3 and ATR also function during unperturbed growth, although the events triggering their activation and their critical functions are largely unknown. Here, we use ChIP-on-chip analysis to map genomic loci decorated by phosphorylated histone H2A (γH2A), a Rad3 substrate that establishes a chromatin-based recruitment platform for Crb2 and Brc1 DNA repair/checkpoint proteins. Unexpectedly, γH2A marks a diverse array of genomic features during S-phase, including natural replication fork barriers and a fork breakage site, retrotransposons, heterochromatin in the centromeres and telomeres, and ribosomal RNA (rDNA) repeats. γH2A formation at the centromeres and telomeres is associated with heterochromatin establishment by Clr4 histone methyltransferase. We show that γH2A domains recruit Brc1, a factor involved in repair of damaged replication forks. Brc1 C-terminal BRCT domain binding to γH2A is crucial in the absence of Rqh1^Sgs1, a RecQ DNA helicase required for rDNA maintenance whose human homologs are mutated in patients with Werner, Bloom, and Rothmund–Thomson syndromes that are characterized by cancer-predisposition or accelerated aging. We conclude that Rad3 phosphorylates histone H2A to mobilize Brc1 to critical genomic domains during S-phase, and this pathway functions in parallel with Rqh1 DNA helicase in maintaining genome integrity.

Eukaryotic genomes, which range in size from ∼10⁷ to ∼10¹¹ base pairs, are replicated with nearly absolute fidelity every cell cycle. This amazing feat happens despite the frequent stalling or collapse of replication forks. The checkpoint kinase ATR is activated by replication fork stalling and phosphorylates histone H2A in nucleosomes surrounding damaged DNA. As the genomic regions triggering ATR activation are largely unknown, we used a whole-genome microarray to map chromosomal domains enriched with phospho-H2A during DNA replication in fission yeast. This analysis identified specific sites, including natural replication fork barriers in ribosomal DNA repeats, retrotransposon elements, and most surprisingly, all heterochromatin regions. Phospho-H2A binds the genome maintenance protein Brc1, and our genetic studies reveal that this molecular pathway becomes crucial in the absence of Rqh1, a conserved DNA helicase that is linked to cancer predisposition. As the fission yeast and human genomes share many similarities, our study reveals genomic landmarks that could similarly trigger ATR activation in human cells and shows that phospho-H2A and Brc1 are a critical part of the network that maintains genome integrity during DNA replication.

Uncovering the interplay between CD8, CD4 and antibody responses to complex pathogens

Vaccinia virus (VACV) was used as the vaccine strain to eradicate smallpox. VACV is still administered to healthcare workers or researchers who are at risk of contracting the virus, and to military personnel. Thus, VACV represents a weapon against outbreaks, both natural (e.g., monkeypox) or man-made (bioterror). This virus is also used as a vector for experimental vaccine development (cancer/infectious disease). As a prototypic poxvirus, VACV is a model system for studying host-pathogen interactions. Until recently, little was known about the targets of host immune responses, which was likely owing to VACVs large genome (>200 open reading frames). However, the last few years have witnessed an explosion of data, and VACV has quickly become a useful model to study adaptive immune responses. This review summarizes and highlights key findings based on identification of VACV antigens targeted by the immune system (CD4, CD8 and antibodies) and the complex interplay between responses.

Production and application of glycan microarrays

Glycans are vital elements of living organisms and are involved in recognition, communication, cell growth and development, motility, and other significant processes. The interactions of glycans with the proteins that bind them provide valuable information about protein interaction and specificity. By printing glycans on microarrays, investigators are able to effectively determine the binding specificity of certain proteins with an extremely efficient and precise result. Such chips are performed by standard robotic microarray printing. Incubating the slides with various GBP-containing substances not only reveals clear receptor preferences of the proteins, but also detects minute differences in structure specificity.

Whole genome transcript profiling from fingerstick blood samples: a comparison and feasibility study

Background

Whole genome gene expression profiling has revolutionized research in the past decade especially with the advent of microarrays. Recently, there have been significant improvements in whole blood RNA isolation techniques which, through stabilization of RNA at the time of sample collection, avoid bias and artifacts introduced during sample handling. Despite these improvements, current human whole blood RNA stabilization/isolation kits are limited by the requirement of a venous blood sample of at least 2.5 mL. While fingerstick blood collection has been used for many different assays, there has yet to be a kit developed to isolate high quality RNA for use in gene expression studies from such small human samples. The clinical and field testing advantages of obtaining reliable and reproducible gene expression data from a fingerstick are many; it is less invasive, time saving, more mobile, and eliminates the need of a trained phlebotomist. Furthermore, this method could also be employed in small animal studies, i.e. mice, where larger sample collections often require sacrificing the animal. In this study, we offer a rapid and simple method to extract sufficient amounts of high quality total RNA from approximately 70 μl of whole blood collected via a fingerstick using a modified protocol of the commercially available Qiagen PAXgene RNA Blood Kit.

Results

From two sets of fingerstick collections, about 70 uL whole blood collected via finger lancet and capillary tube, we recovered an average of 252.6 ng total RNA with an average RIN of 9.3. The post-amplification yields for 50 ng of total RNA averaged at 7.0 ug cDNA. The cDNA hybridized to Affymetrix HG-U133 Plus 2.0 GeneChips had an average % Present call of 52.5%. Both fingerstick collections were highly correlated with r² values ranging from 0.94 to 0.97. Similarly both fingerstick collections were highly correlated to the venous collection with r² values ranging from 0.88 to 0.96 for fingerstick collection 1 and 0.94 to 0.96 for fingerstick collection 2.

Conclusions

Our comparisons of RNA quality and gene expression data of the fingerstick method with traditionally processed sample workflows demonstrate excellent RNA quality from the capillary collection as well as very high correlations of gene expression data.

Detection of differentially expressed wound-healing-related glycogenes in galectin-3-deficient mice

PURPOSE

A prior study showed that exogenous galectin-3 (Gal-3) stimulates re-epithelialization of corneal wounds in wild-type (Gal-3(+/+)) mice but, surprisingly, not in galectin-3-deficient (Gal-3(-/-)) mice. In an effort to understand why the injured corneas of Gal-3(-/-) mice are unresponsive to exogenous Gal-3, the present study was designed to determine whether genes encoding the enzymes that regulate the synthesis of glycan ligands of Gal-3 are differentially expressed in Gal-3(-/-) corneas compared with the Gal-3(+/+) corneas.

METHODS

Glycogene microarray technology was used to identify differentially expressed glycosyltransferases in healing Gal-3(+/+) and Gal-3(-/-) corneas.

RESULTS

Of approximately 2000 glycogenes on the array, the expression of 8 was upregulated and that of 14 was downregulated more than 1.3-fold in healing Gal-3(-/-) corneas. A galactosyltransferase, beta3GalT5, which has the ability to synthesize Gal-3 ligands was markedly downregulated in healing Gal-3(-/-) corneas. The genes for polypeptide galactosaminyltransferases (ppGalNAcT-3 and -7) that are known to initiate O-linked glycosylation and N-aspartyl-beta-glucosaminidase, which participates in the removal of N-glycans, were found to be upregulated in healing Gal-3(-/-) corneas. Microarray data were validated by qRT-PCR.

CONCLUSIONS

Based on the known functions of the differentially expressed glycogenes, it appears that the glycan structures on glycoproteins and glycolipids, synthesized as a result of the differential glycogene expression pattern in healing Gal-3(-/-) corneas may lead to the downregulation of specific counterreceptors for Gal-3. This may explain, at least in part, why, unlike healing Gal-3(+/+) corneas, the healing Gal-3(-/-) corneas are unresponsive to the stimulatory effect of exogenous Gal-3 on re-epithelialization of corneal wounds.

Genetics of platelet reactivity in normal, healthy individuals

BACKGROUND

The Platelet Function Analyzer-100 (PFA-100) is widely used to measure platelet reactivity in whole blood under high shear.

OBJECTIVE

To characterize the genetic component of platelet reactivity among normal individuals, using the PFA-100.

METHODS

We compared baseline platelet reactivity with sex, age, platelet count, hematocrit, plasma von Willebrand factor antigen (VWF:Ag), and alleles of seven candidate genes: integrin subunits alpha2 (ITGA2) and beta3 (ITGB3), platelet glycoproteins GPIbalpha (GP1BA) and GPVI (GP6), purinogenic receptors (P2RY1 and P2RY12) and cyclooxygenase-1 (COX1).

RESULTS

Based on linear and logistic regression models, we report an inverse correlation between baseline closure time (CT) initiated by collagen plus epinephrine (CEPI) and plasma VWF:Ag level, ITGA2 807T and P2RY1 893C, and an inverse correlation between baseline CT initiated by collagen plus adenosine diphosphate (CADP) and P2RY1 893C or GP1BA -5C.

CONCLUSIONS

These results indicate that genetic polymorphisms in ITGA2 and P2RY1 combine with plasma VWF:Ag levels to modulate baseline platelet reactivity in response to collagen plus EPI, while genetic differences in P2RY1 and GP1BA significantly effect platelet responses to collagen plus ADP. Our results demonstrate that the PFA-100 can be used to evaluate the effects of genetic predictors of platelet function.

Analysis of differential expression of glycosyltransferases in healing corneas by glycogene microarrays

It is generally accepted that the glycans on the cell surface and extracellular matrix proteins play a pivotal role in the events that mediate re-epithelialization of wounds. Yet, the global alteration in the structure and composition of glycans, specifically occurring during corneal wound closure remains unknown. In this study, GLYCOv2 glycogene microarray technology was used for the first time to identify the differentially expressed glycosylation-related genes in healing mouse corneas. Of approximately 2000 glycogenes on the array, the expression of 11 glycosytransferase and glycosidase enzymes was upregulated and that of 19 was downregulated more than 1.5-fold in healing corneas compared with the normal, uninjured corneas. Among them, notably, glycosyltransferases, beta3GalT5, T-synthase, and GnTIVb, were all found to be induced in the corneas in response to injury, whereas, GnTIII and many sialyltransferases were downregulated. Interestingly, it appears that the glycan structures on glycoproteins and glycolipids, expressed in healing corneas as a result of differential regulation of these glycosyltransferases, may serve as specific counter-receptors for galectin-3, a carbohydrate-binding protein, known to play a key role in re-epithelialization of corneal wounds. Additionally, many glycogenes including a proteoglycan, glypican-3, cell adhesion proteins dectin-1 and -2, and mincle, and mucin 1 were identified for the first time to be differentially regulated during corneal wound healing. Results of glycogene microarray data were confirmed by qRT-PCR and lectin blot analyses. The differentially expressed glycogenes identified in the present study have not previously been investigated in the context of wound healing and represent novel factors for investigating the role of carbohydrate-mediated recognition in corneal wound healing.