Kisspeptin-10 binding to Gpr54 in osteoclasts prevents bone loss by activating Dusp18-mediated dephosphorylation of Src.

Osteoclasts are over-activated as we age, which results in bone loss. Src deficiency in mice leads to severe osteopetrosis due to a functional defect in osteoclasts, indicating that Src function is essential in osteoclasts. G-protein-coupled receptors (GPCRs) are the targets for ∼35% of approved drugs but it is still unclear how GPCRs regulate Src kinase activity. Here, we reveal that GPR54 activation by its natural ligand Kisspeptin-10 (Kp-10) causes Dusp18 to dephosphorylate Src at Tyr 416. Mechanistically, Gpr54 recruits both active Src and the Dusp18 phosphatase at its proline/arginine-rich motif in its C terminus. We show that Kp-10 binding to Gpr54 leads to the up-regulation of Dusp18. Kiss1, Gpr54 and Dusp18 knockout mice all exhibit osteoclast hyperactivation and bone loss, and Kp-10 abrogated bone loss by suppressing osteoclast activity in vivo. Therefore, Kp-10/Gpr54 is a promising therapeutic target to abrogate bone resorption by Dusp18-mediated Src dephosphorylation.

Association between CYP3A4, CYP3A5 and ABCB1 genotype and tacrolimus treatment outcomes among allogeneic HSCT patients.

Aim: Successful treatment with tacrolimus to prevent graft versus host disease (GVHD) and minimize tacrolimus-related toxicities among allogeneic hematopoietic cell transplantation (alloHCT) recipients is contingent upon quickly achieving and maintaining concentrations within a narrow therapeutic range. The primary objective was to investigate associations between CYP3A4, CYP3A5 or ABCB1 genotype and the proportion of patients that attained an initial tacrolimus goal concentration following initiation of intravenous (iv.) and conversion to oral administration. Materials & methods: We retrospectively evaluated 86 patients who underwent HLA-matched (8/8) related donor alloHCT and were prescribed a tacrolimus-based regimen for GVHD prophylaxis. Results & conclusion: The findings of the present study suggests that CYP3A5 genotype may impact attainment of initial therapeutic tacrolimus concentrations with oral administration in alloHCT recipients.

"Alterations in the Skin Microbiota Are Associated With Symptom Severity in Mycosis Fungoides".

Cutaneous T cell lymphoma (CTCL), a non-Hodgkin lymphoma, is thought to arise from mature tissue-resident memory T cells. The most common subtypes include Mycosis Fungoides and Sezary Syndrome. The role of skin microbiota remains unclear in the symptom manifestation of MF. Among 39 patients with MF, we analyzed bacteria colonizing MF lesions and non-lesional skin in the contralateral side and characterized regional changes in the skin microbiota related to MF involvement using the difference in relative abundance of each genus between lesional and contralateral non-lesional skin. We investigated the relationship between these skin microbiota alterations and symptom severity. No statistically significant difference was found in bacterial diversity and richness between lesional and non-lesional skin. Different skin microbiota signatures were associated with different symptoms. More pronounced erythema in the lesions was associated with an increase in Staphylococcus. Pain and thick skin in the lesions were associated with a decrease in Propionibacterium. The results of this pilot study suggest that the skin microbiota plays an important role in changing skin phenotypes among patients with MF. Larger skin microbiota studies are needed to confirm these findings and support the use of antibiotic treatment to mitigate CTCL symptoms.

Tumor-infiltrating lymphocyte treatment for anti-PD-1-resistant metastatic lung cancer: a phase 1 trial.

Adoptive cell therapy using tumor-infiltrating lymphocytes (TILs) has shown activity in melanoma, but has not been previously evaluated in metastatic non-small cell lung cancer. We conducted a single-arm open-label phase 1 trial ( NCT03215810 ) of TILs administered with nivolumab in 20 patients with advanced non-small cell lung cancer following initial progression on nivolumab monotherapy. The primary end point was safety and secondary end points included objective response rate, duration of response and T cell persistence. Autologous TILs were expanded ex vivo from minced tumors cultured with interleukin-2. Patients received cyclophosphamide and fludarabine lymphodepletion, TIL infusion and interleukin-2, followed by maintenance nivolumab. The end point of safety was met according to the prespecified criteria of ≤17% rate of severe toxicity (95% confidence interval, 3-29%). Of 13 evaluable patients, 3 had confirmed responses and 11 had reduction in tumor burden, with a median best change of 35%. Two patients achieved complete responses that were ongoing 1.5 years later. In exploratory analyses, we found T cells recognizing multiple types of cancer mutations were detected after TIL treatment and were enriched in responding patients. Neoantigen-reactive T cell clonotypes increased and persisted in peripheral blood after treatment. Cell therapy with autologous TILs is generally safe and clinically active and may constitute a new treatment strategy in metastatic lung cancer.

Normal IgH Repertoire Diversity in an Infant with ADA Deficiency After Gene Therapy.

PURPOSE: Adenosine deaminase (ADA) deficiency causes severe combined immunodeficiency (SCID) through an accumulation of toxic metabolites within lymphocytes. Recently, ADA deficiency has been successfully treated using lentiviral-transduced autologous CD34+ cells carrying the ADA gene. T and B cell function appears to be fully restored, but in many patients' B cell numbers remain low, and assessments of the immunoglobulin heavy (IgHV) repertoire following gene therapy are lacking. METHODS: We performed deep sequencing of IgHV repertoire in peripheral blood lymphocytes from a child following lentivirus-based gene therapy for ADA deficiency and compared to the IgHV repertoire in healthy infants and adults. RESULTS: After gene therapy, Ig diversity increased over time as evidenced by V, D, and J gene usage, N-additions, CDR3 length, extent of somatic hypermutation, and Ig class switching. There was the emergence of predominant IgHM, IgHG, and IgHA CDR3 lengths after gene therapy indicating successful oligoclonal expansion in response to antigens. This provides proof of concept for the feasibility and utility of molecular monitoring in following B cell reconstitution following gene therapy for ADA deficiency. CONCLUSION: Based on deep sequencing, gene therapy resulted in an IgHV repertoire with molecular diversity similar to healthy infants.

Fucosylated Proteome Profiling Identifies a Fucosylated, Non-Ribosomal, Stress-Responsive Species of Ribosomal Protein S3.

Alterations in genes encoding for proteins that control fucosylation are known to play causative roles in several developmental disorders, such as Dowling-Degos disease 2 and congenital disorder of glycosylation type IIc (CDGIIc). Recent studies have provided evidence that changes in fucosylation can contribute to the development and progression of several different types of cancers. It is therefore important to gain a detailed understanding of how fucosylation is altered in disease states so that interventions may be developed for therapeutic purposes. In this report, we find that fucosylation occurs on many intracellular proteins. This is an interesting finding, as the fucosylation machinery is restricted to the secretory pathway and is thought to predominately affect cell-membrane-bound and secreted proteins. We find that Ribosomal protein S3 (RPS3) is fucosylated in normal tissues and in cancer cells, and that the extent of its fucosylation appears to respond to stress, including MAPK inhibitors, suggesting a new role in posttranslational protein function. Our data identify a new ribosome-independent species of fucosylated RPS3 that interacts with proteins involved in posttranscriptional regulation of RNA, such as Heterogeneous nuclear ribonucleoprotein U (HNRNPU), as well as with a predominance of non-coding RNAs. These data highlight a novel role for RPS3, which, given previously reported oncogenic roles for RPS3, might represent functions that are perturbed in pathologies such as cancer. Together, our findings suggest a previously unrecognized role for fucosylation in directly influencing intracellular protein functions.

A Mutational Survey of Acral Nevi.

IMPORTANCE: Acral skin may develop nevi, but their mutational status and association with acral melanoma is unclear. OBJECTIVE: To perform targeted next-generation sequencing on a cohort of acral nevi to determine their mutational spectrum. DESIGN, SETTING, AND PARTICIPANTS: Acral nevi specimens (n = 50) that had been obtained for diagnostic purposes were identified from the pathology archives of a tertiary care academic cancer center and a university dermatology clinic. Next-generation sequencing was performed on DNA extracted from the specimens, and mutations called. A subset of samples was stained immunohistochemically for the BRAF V600E mutation. RESULTS: A total of 50 nevi from 49 patients (19 males and 30 females; median [range] age, 48 [13-85] years) were examined. Analysis of the sequencing data revealed a high prevalence of BRAF mutations (n = 43), with a lower frequency of NRAS mutations (n = 5). Mutations in BRAF and NRAS were mutually exclusive. CONCLUSIONS AND RELEVANCE: In this cohort study, nevi arising on mostly sun-protected acral skin showed a rate of BRAF mutation similar to that of acquired nevi on sun-exposed skin but far higher than that of acral melanoma. These findings are in contrast to the well-characterized mutational landscape of acral melanoma.

Eprenetapopt (APR-246) and Azacitidine in TP53-Mutant Myelodysplastic Syndromes.

PURPOSE: Approximately 20% of patients with TP53-mutant myelodysplastic syndromes (MDS) achieve complete remission (CR) with hypomethylating agents. Eprenetapopt (APR-246) is a novel, first-in-class, small molecule that restores wild-type p53 functions in TP53-mutant cells. METHODS: This was a phase Ib/II study to determine the safety, recommended phase II dose, and efficacy of eprenetapopt administered in combination with azacitidine in patients with TP53-mutant MDS or acute myeloid leukemia (AML) with 20%-30% marrow blasts (ClinicalTrials.gov identifier: NCT03072043). RESULTS: Fifty-five patients (40 MDS, 11 AML, and four MDS/myeloproliferative neoplasms) with at least one TP53 mutation were treated. The overall response rate was 71% with 44% achieving CR. Of patients with MDS, 73% (n = 29) responded with 50% (n = 20) achieving CR and 58% (23/40) a cytogenetic response. The overall response rate and CR rate for patients with AML was 64% (n = 7) and 36% (n = 4), respectively. Patients with only TP53 mutations by next-generation sequencing had higher rates of CR (69% v 25%; P = .006). Responding patients had significant reductions in TP53 variant allele frequency and p53 expression by immunohistochemistry, with 21 (38%) achieving complete molecular remission (variant allele frequency < 5%). Median overall survival was 10.8 months with significant improvement in responding versus nonresponding patients by landmark analysis (14.6 v 7.5 months; P = .0005). Overall, 19/55 (35%) patients underwent allogeneic hematopoietic stem-cell transplant, with a median overall survival of 14.7 months. Adverse events were similar to those reported for azacitidine or eprenetapopt monotherapy, with the most common grade ≥ 3 adverse events being febrile neutropenia (33%), leukopenia (29%), and neutropenia (29%). CONCLUSION: Combination treatment with eprenetapopt and azacitidine is well-tolerated yielding high rates of clinical response and molecular remissions in patients with TP53-mutant MDS and oligoblastic AML.

Artificial selection for host resistance to tumour growth and subsequent cancer cell adaptations: an evolutionary arms race.

BACKGROUND: Cancer progression is governed by evolutionary dynamics in both the tumour population and its host. Since cancers die with the host, each new population of cancer cells must reinvent strategies to overcome the host's heritable defences. In contrast, host species evolve defence strategies over generations if tumour development limits procreation. METHODS: We investigate this "evolutionary arms race" through intentional breeding of immunodeficient SCID and immunocompetent Black/6 mice to evolve increased tumour suppression. Over 10 generations, we injected Lewis lung mouse carcinoma cells [LL/2-Luc-M38] and selectively bred the two individuals with the slowest tumour growth at day 11. Their male progeny were hosts in the subsequent round. RESULTS: The evolved SCID mice suppressed tumour growth through biomechanical restriction from increased mesenchymal proliferation, and the evolved Black/6 mice suppressed tumour growth by increasing immune-mediated killing of cancer cells. However, transcriptomic changes of multicellular tissue organisation and function genes allowed LL/2-Luc-M38 cells to adapt through increased matrix remodelling in SCID mice, and reduced angiogenesis, increased energy utilisation and accelerated proliferation in Black/6 mice. CONCLUSION: Host species can rapidly evolve both immunologic and non-immunologic tumour defences. However, cancer cell plasticity allows effective phenotypic and population-based counter strategies.

Sirt2 Inhibition Enhances Metabolic Fitness and Effector Functions of Tumor-Reactive T Cells.

Dysregulated metabolism is a key driver of maladaptive tumor-reactive T lymphocytes within the tumor microenvironment. Actionable targets that rescue the effector activity of antitumor T cells remain elusive. Here, we report that the Sirtuin-2 (Sirt2) NAD+-dependent deacetylase inhibits T cell metabolism and impairs T cell effector functions. Remarkably, upregulation of Sirt2 in human tumor-infiltrating lymphocytes (TILs) negatively correlates with response to TIL therapy in advanced non-small-cell lung cancer. Mechanistically, Sirt2 suppresses T cell metabolism by targeting key enzymes involved in glycolysis, tricarboxylic acid-cycle, fatty acid oxidation, and glutaminolysis. Accordingly, Sirt2-deficient murine T cells exhibit increased glycolysis and oxidative phosphorylation, resulting in enhanced proliferation and effector functions and subsequently exhibiting superior antitumor activity. Importantly, pharmacologic inhibition of Sirt2 endows human TILs with these superior metabolic fitness and effector functions. Our findings unveil Sirt2 as an unexpected actionable target for reprogramming T cell metabolism to augment a broad spectrum of cancer immunotherapies.

The Homeobox gene, HOXB13, Regulates a Mitotic Protein-Kinase Interaction Network in Metastatic Prostate Cancers.

HOXB13, a homeodomain transcription factor, is linked to recurrence following radical prostatectomy. While HOXB13 regulates Androgen Receptor (AR) functions in a context dependent manner, its critical effectors in prostate cancer (PC) metastasis remain largely unknown. To identify HOXB13 transcriptional targets in metastatic PCs, we performed integrative bioinformatics analysis of differentially expressed genes (DEGs) in the proximity of the human prostate tumor-specific AR binding sites. Unsupervised Principal Component Analysis (PCA) led to a focused core HOXB13 target gene-set referred to as HOTPAM9 (HOXB13 Targets separating Primary And Metastatic PCs). HOTPAM9 comprised 7 mitotic kinase genes overexpressed in metastatic PCs, TRPM8, and the heat shock protein HSPB8, whose levels were significantly lower in metastatic PCs compared to the primary disease. The expression of a two-gene set, CIT and HSPB8 with an overall balanced accuracy of 98.8% and a threshold value of 0.2347, was sufficient to classify metastasis. HSPB8 mRNA expression was significantly increased following HOXB13 depletion in multiple metastatic CRPC models. Increased expression of HSPB8 by the microtubule inhibitor Colchicine or by exogenous means suppressed migration of mCRPC cells. Collectively, our results indicate that HOXB13 promotes metastasis of PCs by coordinated regulation of mitotic kinases and blockade of a putative tumor suppressor gene.

Targeting the BRD4-HOXB13 Coregulated Transcriptional Networks with Bromodomain-Kinase Inhibitors to Suppress Metastatic Castration-Resistant Prostate Cancer.

Resistance to androgen receptor (AR) antagonists is a significant problem in the treatment of castration-resistant prostate cancers (CRPC). Identification of the mechanisms by which CRPCs evade androgen deprivation therapies (ADT) is critical to develop novel therapeutics. We uncovered that CRPCs rely on BRD4-HOXB13 epigenetic reprogramming for androgen-independent cell proliferation. Mechanistically, BRD4, a member of the BET bromodomain family, epigenetically promotes HOXB13 expression. Consistently, genetic disruption of HOXB13 or pharmacological suppression of its mRNA and protein expression by the novel dual-activity BET bromodomain-kinase inhibitors directly correlates with rapid induction of apoptosis, potent inhibition of tumor cell proliferation and cell migration, and suppression of CRPC growth. Integrative analysis revealed that the BRD4-HOXB13 transcriptome comprises a proliferative gene network implicated in cell-cycle progression, nucleotide metabolism, and chromatin assembly. Notably, although the core HOXB13 target genes responsive to BET inhibitors (HOTBIN10) are overexpressed in metastatic cases, in ADT-treated CRPC cell lines and patient-derived circulating tumor cells (CTC) they are insensitive to AR depletion or blockade. Among the HOTBIN10 genes, AURKB and MELK expression correlates with HOXB13 expression in CTCs of mCRPC patients who did not respond to abiraterone (ABR), suggesting that AURKB inhibitors could be used additionally against high-risk HOXB13-positive metastatic prostate cancers. Combined, our study demonstrates that BRD4-HOXB13-HOTBIN10 regulatory circuit maintains the malignant state of CRPCs and identifies a core proproliferative network driving ADT resistance that is targetable with potent dual-activity bromodomain-kinase inhibitors.

Reprogramming of a defense signaling pathway in rough lemon and sweet orange is a critical element of the early response to 'Candidatus Liberibacter asiaticus'.

Huanglongbing (HLB) in citrus infected by Candidatus Liberibacter asiaticus (CLas) has caused tremendous losses to the citrus industry. No resistant genotypes have been identified in citrus species or close relatives. Among citrus varieties, rough lemon (Citrus jambhiri) has been considered tolerant due to its ability to produce a healthy flush of new growth after infection. The difference between tolerance and susceptibility is often defined by the speed and intensity of a plant's response to a pathogen, especially early defense responses. RNA-seq data were collected from three biological replicates of CLas- and mock-inoculated rough lemon and sweet orange at week 0 and 7 following infection. Functional analysis of the differentially expressed genes (DEGs) indicated that genes involved in the mitogen activated protein kinase (MAPK) signaling pathway were highly upregulated in rough lemon. MAPK induces the transcription of WRKY and other transcription factors which potentially turn on multiple defense-related genes. A Subnetwork Enrichment Analysis further revealed different patterns of regulation of several functional categories, suggesting DEGs with different functions were subjected to reprogramming. In general, the amplitude of the expression of defense-related genes is much greater in rough lemon than in sweet orange. A quantitative disease resistance response may contribute to the durable tolerance level to HLB observed in rough lemon.

ARG1 Functions in the Physiological Adaptation of Undifferentiated Plant Cells to Spaceflight.

Scientific access to spaceflight and especially the International Space Station has revealed that physiological adaptation to spaceflight is accompanied or enabled by changes in gene expression that significantly alter the transcriptome of cells in spaceflight. A wide range of experiments have shown that plant physiological adaptation to spaceflight involves gene expression changes that alter cell wall and other metabolisms. However, while transcriptome profiling aptly illuminates changes in gene expression that accompany spaceflight adaptation, mutation analysis is required to illuminate key elements required for that adaptation. Here we report how transcriptome profiling was used to gain insight into the spaceflight adaptation role of Altered response to gravity 1 (Arg1), a gene known to affect gravity responses in plants on Earth. The study compared expression profiles of cultured lines of Arabidopsis thaliana derived from wild-type (WT) cultivar Col-0 to profiles from a knock-out line deficient in the gene encoding ARG1 (ARG1 KO), both on the ground and in space. The cell lines were launched on SpaceX CRS-2 as part of the Cellular Expression Logic (CEL) experiment of the BRIC-17 spaceflight mission. The cultured cell lines were grown within 60 mm Petri plates in Petri Dish Fixation Units (PDFUs) that were housed within the Biological Research In Canisters (BRIC) hardware. Spaceflight samples were fixed on orbit. Differentially expressed genes were identified between the two environments (spaceflight and comparable ground controls) and the two genotypes (WT and ARG1 KO). Each genotype engaged unique genes during physiological adaptation to the spaceflight environment, with little overlap. Most of the genes altered in expression in spaceflight in WT cells were found to be Arg1-dependent, suggesting a major role for that gene in the physiological adaptation of undifferentiated cells to spaceflight. Key Words: ARG1-Spaceflight-Gene expression-Physiological adaptation-BRIC. Astrobiology 17, 1077-1111.

Tumor exome sequencing and copy number alterations reveal potential predictors of intrinsic resistance to multi-targeted tyrosine kinase inhibitors.

Multi-targeted tyrosine kinase inhibitors (TKIs) have broad efficacy and similar FDA-approved indications, suggesting shared molecular drug targets across cancer types. Irrespective of tumor type, 20-30% of patients treated with multi-targeted TKIs demonstrate intrinsic resistance, with progressive disease as a best response. We conducted a retrospective cohort study to identify tumor (somatic) point mutations, insertion/deletions, and copy number alterations (CNA) associated with intrinsic resistance to multi-targeted TKIs. Using a candidate gene approach (n=243), tumor next-generation sequencing and CNA data was associated with resistant and non-resistant outcomes. Resistant individuals (n=11) more commonly harbored somatic point mutations in NTRK1, KDR, TGFBR2, and PTPN11 and CNA in CDK4, CDKN2B, and ERBB2 compared to non-resistant (n=26, p<0.01). Using a random forest classification model for variable reduction and a decision tree classification model, we were able to differentiate intrinsically resistant from non-resistant patients. CNA in CDK4 and CDKN2B were the most important analytical features, implicating the cyclin D pathway as a potentially important factor in resistance to multi-targeted TKIs. Replication of these results in a larger, independent patient cohort has potential to inform personalized prescribing of these widely utilized agents.

Elongator subunit 3 positively regulates plant immunity through its histone acetyltransferase and radical S-adenosylmethionine domains.

BACKGROUND: Pathogen infection triggers a large-scale transcriptional reprogramming in plants, and the speed of this reprogramming affects the outcome of the infection. Our understanding of this process has significantly benefited from mutants that display either delayed or accelerated defense gene induction. In our previous work we demonstrated that the Arabidopsis Elongator complex subunit 2 (AtELP2) plays an important role in both basal immunity and effector-triggered immunity (ETI), and more recently showed that AtELP2 is involved in dynamic changes in histone acetylation and DNA methylation at several defense genes. However, the function of other Elongator subunits in plant immunity has not been characterized. RESULTS: In the same genetic screen used to identify Atelp2, we found another Elongator mutant, Atelp3-10, which mimics Atelp2 in that it exhibits a delay in defense gene induction following salicylic acid treatment or pathogen infection. Similarly to AtELP2, AtELP3 is required for basal immunity and ETI, but not for systemic acquired resistance (SAR). Furthermore, we demonstrate that both the histone acetyltransferase and radical S-adenosylmethionine domains of AtELP3 are essential for its function in plant immunity. CONCLUSION: Our results indicate that the entire Elongator complex is involved in basal immunity and ETI, but not in SAR, and support that Elongator may play a role in facilitating the transcriptional induction of defense genes through alterations to their chromatin.

The Arabidopsis elongator complex subunit2 epigenetically regulates plant immune responses.

The Arabidopsis thaliana Elongator complex subunit2 (ELP2) genetically interacts with NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 (NPR1), a key transcription coactivator of plant immunity, and regulates the induction kinetics of defense genes. However, the mechanistic relationship between ELP2 and NPR1 and how ELP2 regulates the kinetics of defense gene induction are unclear. Here, we demonstrate that ELP2 is an epigenetic regulator required for pathogen-induced rapid transcriptome reprogramming. We show that ELP2 functions in a transcriptional feed-forward loop regulating both NPR1 and its target genes. An elp2 mutation increases the total methylcytosine number, reduces the average methylation levels of methylcytosines, and alters (increases or decreases) methylation levels of specific methylcytosines. Interestingly, infection of plants with the avirulent bacterial pathogen Pseudomonas syringae pv tomato DC3000/avrRpt2 induces biphasic changes in DNA methylation levels of NPR1 and PHYTOALEXIN DEFICIENT4 (PAD4), which encodes another key regulator of plant immunity. These dynamic changes are blocked by the elp2 mutation, which is correlated with delayed induction of NPR1 and PAD4. The elp2 mutation also reduces basal histone acetylation levels in the coding regions of several defense genes. Together, our data demonstrate a new role for Elongator in somatic DNA demethylation/methylation and suggest a function for Elongator-mediated chromatin regulation in pathogen-induced transcriptome reprogramming.

Characterization of Capsicum annuum genetic diversity and population structure based on parallel polymorphism discovery with a 30K unigene Pepper GeneChip.

The widely cultivated pepper, Capsicum spp., important as a vegetable and spice crop world-wide, is one of the most diverse crops. To enhance breeding programs, a detailed characterization of Capsicum diversity including morphological, geographical and molecular data is required. Currently, molecular data characterizing Capsicum genetic diversity is limited. The development and application of high-throughput genome-wide markers in Capsicum will facilitate more detailed molecular characterization of germplasm collections, genetic relationships, and the generation of ultra-high density maps. We have developed the Pepper GeneChip® array from Affymetrix for polymorphism detection and expression analysis in Capsicum. Probes on the array were designed from 30,815 unigenes assembled from expressed sequence tags (ESTs). Our array design provides a maximum redundancy of 13 probes per base pair position allowing integration of multiple hybridization values per position to detect single position polymorphism (SPP). Hybridization of genomic DNA from 40 diverse C. annuum lines, used in breeding and research programs, and a representative from three additional cultivated species (C. frutescens, C. chinense and C. pubescens) detected 33,401 SPP markers within 13,323 unigenes. Among the C. annuum lines, 6,426 SPPs covering 3,818 unigenes were identified. An estimated three-fold reduction in diversity was detected in non-pungent compared with pungent lines, however, we were able to detect 251 highly informative markers across these C. annuum lines. In addition, an 8.7 cM region without polymorphism was detected around Pun1 in non-pungent C. annuum. An analysis of genetic relatedness and diversity using the software Structure revealed clustering of the germplasm which was confirmed with statistical support by principle components analysis (PCA) and phylogenetic analysis. This research demonstrates the effectiveness of parallel high-throughput discovery and application of genome-wide transcript-based markers to assess genetic and genomic features among Capsicum annuum.

Changes in the transcriptome of morula-stage bovine embryos caused by heat shock: relationship to developmental acquisition of thermotolerance.

BACKGROUND: While initially sensitive to heat shock, the bovine embryo gains thermal resistance as it progresses through development so that physiological heat shock has little effect on development to the blastocyst stage by Day 5 after insemination. Here, experiments using 3' tag digital gene expression (3'DGE) and real-time PCR were conducted to determine changes in the transcriptome of morula-stage bovine embryos in response to heat shock (40 degrees C for 8 h) that could be associated with thermotolerance. RESULTS: Using 3'DGE, expression of 173 genes were modified by heat shock, with 94 genes upregulated by heat shock and 79 genes downregulated by heat shock. A total of 38 differentially-regulated genes were associated with the ubiquitin protein, UBC. Heat shock increased expression of one heat shock protein gene, HSPB11, and one heat shock protein binding protein, HSPBP1, tended to increase expression of HSPA1A and HSPB1, but did not affect expression of 64 other genes encoding heat shock proteins, heat shock transcription factors or proteins interacting with heat shock proteins. Moreover, heat shock increased expression of five genes associated with oxidative stress (AKR7A2, CBR1, GGH, GSTA4, and MAP2K5), decreased expression of HIF3A, but did not affect expression of 42 other genes related to free radical metabolism. Heat shock also had little effect on genes involved in embryonic development. Effects of heat shock for 2, 4 and 8 h on selected heat shock protein and antioxidant genes were also evaluated by real-time PCR. Heat shock increased steady-state amounts of mRNA for HSPA1A (P<0.05) and tended to increase expression of HSP90AA1 (P<0.07) but had no effect on expression of SOD1 or CAT. CONCLUSIONS: Changes in the transcriptome of the heat-shocked bovine morula indicate that the embryo is largely resistant to effects of heat shock. As a result, transcription of genes involved in thermal protection is muted and there is little disruption of gene networks involved in embryonic development. It is likely that the increased resistance of morula-stage embryos to heat shock as compared to embryos at earlier stages of development is due in part to developmental acquisition of mechanisms to prevent accumulation of denatured proteins and free radical damage.

De novo assembly of the pepper transcriptome (Capsicum annuum): a benchmark for in silico discovery of SNPs, SSRs and candidate genes.

BACKGROUND: Molecular breeding of pepper (Capsicum spp.) can be accelerated by developing DNA markers associated with transcriptomes in breeding germplasm. Before the advent of next generation sequencing (NGS) technologies, the majority of sequencing data were generated by the Sanger sequencing method. By leveraging Sanger EST data, we have generated a wealth of genetic information for pepper including thousands of SNPs and Single Position Polymorphic (SPP) markers. To complement and enhance these resources, we applied NGS to three pepper genotypes: Maor, Early Jalapeño and Criollo de Morelos-334 (CM334) to identify SNPs and SSRs in the assembly of these three genotypes. RESULTS: Two pepper transcriptome assemblies were developed with different purposes. The first reference sequence, assembled by CAP3 software, comprises 31,196 contigs from >125,000 Sanger-EST sequences that were mainly derived from a Korean F1-hybrid line, Bukang. Overlapping probes were designed for 30,815 unigenes to construct a pepper Affymetrix GeneChip® microarray for whole genome analyses. In addition, custom Python scripts were used to identify 4,236 SNPs in contigs of the assembly. A total of 2,489 simple sequence repeats (SSRs) were identified from the assembly, and primers were designed for the SSRs. Annotation of contigs using Blast2GO software resulted in information for 60% of the unigenes in the assembly. The second transcriptome assembly was constructed from more than 200 million Illumina Genome Analyzer II reads (80-120 nt) using a combination of Velvet, CLC workbench and CAP3 software packages. BWA, SAMtools and in-house Perl scripts were used to identify SNPs among three pepper genotypes. The SNPs were filtered to be at least 50 bp from any intron-exon junctions as well as flanking SNPs. More than 22,000 high-quality putative SNPs were identified. Using the MISA software, 10,398 SSR markers were also identified within the Illumina transcriptome assembly and primers were designed for the identified markers. The assembly was annotated by Blast2GO and 14,740 (12%) of annotated contigs were associated with functional proteins. CONCLUSIONS: Before availability of pepper genome sequence, assembling transcriptomes of this economically important crop was required to generate thousands of high-quality molecular markers that could be used in breeding programs. In order to have a better understanding of the assembled sequences and to identify candidate genes underlying QTLs, we annotated the contigs of Sanger-EST and Illumina transcriptome assemblies. These and other information have been curated in a database that we have dedicated for pepper project.