Search

The Molecular Pathology Section, Cleveland Clinic (Cleveland, OH), has undergone enhancement of its testing portfolio and processes. An Excel 2013- and paper-based data-management system was replaced with a commercially available laboratory information-management system (LIMS) software application, a separate bioinformatics platform, customized test-interpretation applications, a dedicated sample-accessioning service, and a results-releasing software application. The customized LIMS solution manages complex workflows, large-scale data packets, and process automation. A customized approach was required because, in a survey of commercially available off-the-shelf software products, none met the diverse and complex needs of this molecular diagnostics service. The project utilized the expertise of clinical laboratorians, pathologists, genetics counselors, bioinformaticians, and systems analysts in partnering with software-engineering consultants to design and implement a solution. Concurrently, Agile software-building best practices were formulated, which may be emulated for scalable and cost-effective laboratory-authored software.

Mutations in the viral genome of SARS-CoV-2 can impact the performance of molecular diagnostic assays. In some cases, such as S gene target failure, the impact can serve as a unique indicator of a particular SARS-CoV-2 variant and provide a method for rapid detection. Here we describe partial ORF1ab gene target failure (pOGTF) on the cobas® SARS-CoV-2 assays, defined by a ≥2 thermocycles delay in detection of the ORF1ab gene compared to the E gene. We demonstrate that pOGTF is 97% sensitive and 99% specific for SARS-CoV-2 lineage BA.2.12.1, an emerging variant in the United States with spike L452Q and S704L mutations that may impact transmission, infectivity, and/or immune evasion. Increasing rates of pOGTF closely mirrored rates of BA.2.12.1 sequences uploaded to public databases, and, importantly increasing local rates of pOGTF also mirrored increasing overall test positivity. Use of pOGTF as a proxy for BA.2.12.1 provides faster tracking of the variant than whole-genome sequencing and can benefit laboratories without sequencing capabilities.

(1) Background: EWS fusion genes are associated with Ewing sarcoma and other Ewing family tumors including desmoplastic small round tumor, DSRCT. We utilize a clinical genomics workflow to reveal real-world frequencies of EWS fusion events, cataloging events that are similar, or divergent at the EWS breakpoint. (2) Methods: EWS fusion events from our next-generation sequencing panel (NGS) samples were first sorted by breakpoint or fusion junctions to map out the frequency of breakpoints. Fusion results were illustrated as in-frame fusion peptides involving EWS and a partner gene. (3) Results: From 2471 patient pool samples for fusion analysis at the Cleveland Clinic Molecular Pathology Laboratory, we identified 182 fusion samples evolved with the EWS gene. They are clustered in several breakpoints: chr22:29683123 (65.9%), and chr22:29688595 (2.7%). About 3/4 of Ewing sarcoma and DSRCT tumors have an identical EWS breakpoint motif at Exon 7 (SQQSSSYGQQ-) fused to a specific part of FLI1 (NPSYDSVRRG or-SSLLAYNTSS), ERG (NLPYEPPRRS), FEV (NPVGDGLFKD) or WT1 (SEKPYQCDFK). Our method also worked with Caris transcriptome data, too. Our primary clinical utility is to use this information to identify neoantigens for therapeutic purposes. (4) Conclusions and future perspectives: our method allows interpretation of what peptides result from the in-frame translation of EWS fusion junctions. These sequences, coupled with HLA-peptide binding data, are used to identify potential sequences of cancer-specific immunogenic peptides for Ewing sarcoma or DSRCT patients. This information may also be useful for immune monitoring (e.g., circulating T-cells with fusion-peptide specificity) to detect vaccine candidates, responses, or residual disease.

FLT3-internal tandem duplication occurs in 20–30% of acute myeloid leukemia and confers an adverse prognosis with its allelic ratio being a key risk stratifier. The US Food and Drug Administration recently approved FLT3 inhibitors midostaurin and gilteritinib in FLT3 mutation-positive acute myeloid leukemia. Historically, FLT3 was tested by fragment analysis, which has become the standard method endorsed by international guidelines. However, next generation sequencing is increasingly used at acute myeloid leukemia diagnosis given its ability to simultaneously evaluate multiple clinically informative markers. As FLT3-internal tandem duplication detection was known to be challenging by next generation sequencing and the results carry profound prognostic and therapeutic implications, it is important to thoroughly examine its performance in FLT3-internal tandem duplication detection and allelic ratio classification. In a comparative study with fragment analysis, we retrospectively reviewed our experience using a custom-designed, hybridization capture-based, targeted next generation sequencing panel. Among 7902 cases, FLT3-internal tandem duplication was detected in 335 with variable sizes (3–231 bp) and insertion sites. Fragment analysis was also performed in 402 cases, demonstrating 100% concordance in FLT3-internal tandem duplication detection. In 136 dual-tested, positive cases, 128/136 (94%) exhibited concordant high/low allelic ratio classifications. The remaining 6% showed borderline low allelic ratio by next generation sequencing. The two methods were concordant in FLT3-tyrosine kinase domain mutation detection at the hotspot D835/I836 targeted by fragment analysis. Furthermore, seven mutations which may benefit from FLT3 inhibitor therapy were detected by next generation sequencing, in regions not covered by fragment analysis. Our study demonstrates that using a hybridization capture-based chemistry and optimized bioinformatics pipeline, next generation sequencing can reliably detect FLT3-internal tandem duplication and classify its allelic ratio for acute myeloid leukemia risk stratification. Next generation sequencing also exhibits superior comprehensiveness in FLT3 mutation detection and may further improve personalized, targeted therapy in acute myeloid leukemia.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to mutate during the ongoing COVID-19 pandemic, and some of the nucleotide polymorphisms may result in diagnostic detection failures.….

Context.— Bone and soft tissue tumors are heterogeneous, diagnostically challenging, and often defined by gene fusions. Objective.— To present our experience using a custom 34-gene targeted sequencing fusion panel. Design.— Total nucleic acid extracted from formalin-fixed, paraffin-embedded (FFPE) tumor specimens was subjected to open-ended, nested anchored multiplex polymerase chain reaction and enrichment of 34 gene targets, thus enabling detection of known and novel fusion partners. Results.— During a 12-month period, 147 patients were tested as part of routine clinical care. Tumor percentage ranged from 10% to 100% and turnaround time ranged from 3 to 15 (median, 7.9) days. The most common diagnostic groups were small round blue cell tumors, tumors of uncertain differentiation, fibroblastic/myofibroblastic tumors, and adipocytic tumors. In-frame fusion transcripts were identified in 64 of 142 cases sequenced (45%): in 62 cases, the detection of a disease-defining fusion confirmed the morphologic impression; in 2 cases, a germline TFG-GPR128 polymorphic fusion variant was detected. Several genes in the panel partnered with multiple fusion partners specific for different diagnoses, for example, EWSR1, NR4A3, FUS, NCOA2, and TFE3. Interesting examples are presented to highlight how fusion detection or lack thereof was instrumental in establishing accurate diagnoses. Novel fusion partners were detected for 2 cases of solid aneurysmal bone cysts (PTBP1-USP6, SLC38A2-USP6). Conclusions.— Multiplex detection of fusions in total nucleic acid purified from FFPE specimens facilitates diagnosis of bone and soft tissue tumors. This technology is particularly useful for morphologically challenging entities and in the absence of prior knowledge of fusion partners, and has the potential to discover novel fusion partners.

// Ferga C. Gleeson 1 , Jesse S. Voss 2 , Benjamin R. Kipp 2 , Sarah E. Kerr 2 , John S. Van Arnam 2 , John R. Mills 2 , Cherisse A. Marcou 2 , Amber R. Schneider 2 , Zheng Jin Tu 3 , Michael R. Henry 2 and Michael J. Levy 1 1 Division of Gastroenterology & Hepatology, Mayo Clinic Rochester, MN, USA 2 Department of Laboratory Medicine-Pathology, Mayo Clinic Rochester, MN, USA 3 Divisions of Biomedical Statics & Informatics, Department of Health Sciences Research, Mayo Clinic Rochester, MN, USA Correspondence to: Ferga C. Gleeson, email: gleeson.ferga@mayo.edu Keywords: pancreas neuroendocrine tumor, endoscopic ultrasound fine needle aspiration, targeted next-generation sequencing, prognostic biomarker, predictive biomarker Received: April 14, 2017 Accepted: May 23, 2017 Published: June 28, 2017 ABSTRACT Background: Recent genetic studies have highlighted that alterations in MEN1, chromatin remodeling genes, and mammalian target of rapamycin (mTOR) pathway genes are the most frequent molecular events identified in pancreas neuroendocrine tumors (pNETs). The prognostic or predictive impact of these biomarkers and other less frequently observed aberrations, i.e. PTEN, TSC2 and PIK3CA are relatively unknown. The aims of this targeted next generation sequencing (NGS) study were to assess tumor cytology genotype diversity, to survey for potential adverse prognostic biomarkers and the prevalence of mTOR pathway variants. Methods: Using a custom 15 gene gastroenteropancreatic neuroendocrine tumor panel, targeted NGS of archived (2002-2013) primary pNETs (n=90) and pNET liver metastasis (n=32) cytology smears was performed. Results: The genetic variant landscape revealed that 21% and 28% of primary and metastatic liver pNETs harbored ≥ 2 variants per tumor, respectively. The most prevalent primary tumor variants were in the MEN1 (42%), DAXX (11%), ATRX (10%), and TSC2 (8%) genes. Patients harboring aberrations in TSC2, KRAS or TP53 were more likely to experience disease progression and reduced overall survival, when compared to individuals who were wild-type. The prevalence of these potential prognostic biomarkers in early disease was observed in 3.3% of the primary tumor cohort. mTOR pathway variants including alterations in PTEN, TSC2 and PIK3CA were identified in 10% and 12.5% of primary tumors and pNET liver metastasis, respectively. Conclusion: Cytology based tumor genotyping revealed a broad spectrum of genetic variants including possible adverse prognostic biomarkers, reflective of an aggressive phenotype. It also demonstrated the prevalence of potential predictive biomarkers for mTOR pathway inhibitor sensitivity.

Key Points Question Are SARS-CoV-2 variants, virus clades, or clade groups associated with disease severity and patient outcomes? Findings In this cross-sectional study of 302 SARS-CoV-2 isolates, 6 different Global Initiative on Sharing All Influenza Data clades circulated in the community followed by a rapid reduction in clade diversity. Several variants, including 23403A>G (D614G), were significantly associated with lower hospitalization rates and increased patient survival. Meaning These findings suggest that SARS-CoV-2 clade assignment is an important factor that may aid in estimating patient outcomes., This cross-sectional study examines the association of identified SARS-CoV-2 variants, virus clades, and clade groups with disease severity and patient outcomes., Importance Understanding of SARS-CoV-2 variants that alter disease outcomes are important for clinical risk stratification and may provide important clues to the complex virus-host relationship. Objective To examine the association of identified SARS-CoV-2 variants, virus clades, and clade groups with disease severity and patient outcomes. Design, Setting, and Participants In this cross-sectional study, viral genome analysis of clinical specimens obtained from patients at the Cleveland Clinic infected with SARS-CoV-2 during the initial wave of infection (March 11 to April 22, 2020) was performed. Identified variants were matched with clinical outcomes. Data analysis was performed from April to July 2020. Main Outcomes and Measures Hospitalization, intensive care unit (ICU) admission, mortality, and laboratory outcomes were matched with SARS-CoV-2 variants. Results Specimens sent for viral genome sequencing originated from 302 patients with SARS-CoV-2 infection (median [interquartile range] age, 52.6 [22.8 to 82.5] years), of whom 126 (41.7%) were male, 195 (64.6%) were White, 91 (30.1%) required hospitalization, 35 (11.6%) needed ICU admission, and 17 (5.6%) died. From these specimens, 2531 variants (484 of which were unique) were identified. Six different SARS-CoV-2 clades initially circulated followed by a rapid reduction in clade diversity. Several variants were associated with lower hospitalization rate, and those containing 23403A>G (D614G Spike) were associated with increased survival when the patient was hospitalized (64 of 74 patients [86.5%] vs 10 of 17 patients [58.8%]; χ21 = 6.907; P = .009). Hospitalization and ICU admission were similar regardless of clade. Infection with Clade V variants demonstrated higher creatinine levels (median [interquartile range], 2.6 [−0.4 to 5.5] mg/dL vs 1.0 [0.2 to 2.2] mg/dL; mean creatinine difference, 2.9 mg/dL [95% CI, 0.8 to 5.0 mg/dL]; Kruskal-Wallis P = .005) and higher overall mortality rates (3 of 14 patients [21.4%] vs 17 of 302 patients [5.6%]; χ21 = 5.640; P = .02) compared with other variants. Infection by strains lacking the 23403A>G variant showed higher mortality in multivariable analysis (odds ratio [OR], 22.4; 95% CI, 0.6 to 5.6; P = .01). Increased variants of open reading frame (ORF) 3a were associated with decreased hospitalization frequency (OR, 0.4; 95% CI, 0.2 to 0.96; P = .04), whereas increased variants of Spike (OR, 0.01; 95% CI, G (D614G) specific genotypes occurred. Replaced clades were associated with worse clinical outcomes, including mortality. These findings help explain persistent hospitalization yet decreasing mortality as the pandemic progresses. SARS-CoV-2 clade assignment is an important factor that may aid in estimating patient outcomes.

FLT3-internal tandem duplication occurs in 20–30% of acute myeloid leukemia and confers an adverse prognosis with its allelic ratio being a key risk stratifier. The US Food and Drug Administration recently approved FLT3 inhibitors midostaurin and gilteritinib in FLT3 mutation-positive acute myeloid leukemia. Historically, FLT3 was tested by fragment analysis, which has become the standard method endorsed by international guidelines. However, next generation sequencing is increasingly used at acute myeloid leukemia diagnosis given its ability to simultaneously evaluate multiple clinically informative markers. As FLT3-internal tandem duplication detection was known to be challenging by next generation sequencing and the results carry profound prognostic and therapeutic implications, it is important to thoroughly examine its performance in FLT3-internal tandem duplication detection and allelic ratio classification. In a comparative study with fragment analysis, we retrospectively reviewed our experience using a custom-designed, hybridization capture-based, targeted next generation sequencing panel. Among 7902 cases, FLT3-internal tandem duplication was detected in 335 with variable sizes (3–231 bp) and insertion sites. Fragment analysis was also performed in 402 cases, demonstrating 100% concordance in FLT3-internal tandem duplication detection. In 136 dual-tested, positive cases, 128/136 (94%) exhibited concordant high/low allelic ratio classifications. The remaining 6% showed borderline low allelic ratio by next generation sequencing. The two methods were concordant in FLT3-tyrosine kinase domain mutation detection at the hotspot D835/I836 targeted by fragment analysis. Furthermore, seven mutations which may benefit from FLT3 inhibitor therapy were detected by next generation sequencing, in regions not covered by fragment analysis. Our study demonstrates that using a hybridization capture-based chemistry and optimized bioinformatics pipeline, next generation sequencing can reliably detect FLT3-internal tandem duplication and classify its allelic ratio for acute myeloid leukemia risk stratification. Next generation sequencing also exhibits superior comprehensiveness in FLT3 mutation detection and may further improve personalized, targeted therapy in acute myeloid leukemia.

International audience; The fungus Cochliobolus miyabeanus causes severe leaf spot disease on rice (Oryza sativa) and two North American specialty crops, American wildrice (Zizania palustris) and switchgrass (Panicum virgatum). Despite the importance of C. miyabeanus as a disease-causing agent in wildrice, little is known about either the mechanisms of pathogenicity or host defense responses. To start bridging these gaps, the genome of C. miyabeanus strain TG12bL2 was shotgun sequenced using Illumina technology. The genome assembly consists of 31.79 Mbp in 2,378 scaffolds with an N50 = 74,921. It contains 11,000 predicted genes of which 94.5% were annotated. Approximately 10% of total gene number is expected to be secreted. The C. miyabeanus genome is rich in carbohydrate active enzymes, and harbors 187 small secreted peptides (SSPs) and some fungal effector homologs. Detoxification systems were represented by a variety of enzymes that could offer protection against plant defense compounds. The non-ribosomal peptide synthetases and polyketide synthases (PKS) present were common to other Cochliobolus species. Additionally, the fungal transcriptome was analyzed at 48 hours after inoculation in planta. A total of 10,674 genes were found to be expressed, some of which are known to be involved in pathogenicity or response to host defenses including hydrophobins, cutinase, cell wall degrading enzymes, enzymes related to reactive oxygen species scavenging, PKS, detoxification systems, SSPs, and a known fungal effector. This work will facilitate future research on C. miyabeanus pathogen-associated molecular patterns and effectors, and in the identification of their corresponding wildrice defense mechanisms.

// Ferga C. Gleeson 1 , Sarah E. Kerr 2 , Benjamin R. Kipp 2 , Jesse S. Voss 2 , Douglas M. Minot 2 , Zheng Jin Tu 3 , Michael R. Henry 2 , Rondell P. Graham 2 , George Vasmatzis 4 , John C. Cheville 4 , Konstantinos N. Lazaridis 1, 4 , Michael J. Levy 1 1 Division of Gastroenterology & Hepatology, Mayo Clinic Rochester, MN, USA 2 Department of Laboratory Medicine & Pathology, Mayo Clinic Rochester, MN, USA 3 Division of Biomedical Statics & Informatics, Department of Health Sciences Research, Mayo Clinic Rochester, MN, USA 4 Center for Individualized Medicine, Mayo Clinic, Rochester, MN, USA Correspondence to: Ferga C. Gleeson, email: gleeson.ferga@mayo.edu Keywords: endoscopic ultrasound fine needle aspiration, pancreatic adenocarcinoma, targeted next-generation sequencing, mutation concordance, personalized medicine Received: April 04, 2016 Accepted: April 24, 2016 Published: May 18, 2016 ABSTRACT Background & Aims: Less than 10% of registered drug intervention trials for pancreatic ductal adenocarcinoma (PDAC) include a biomarker stratification strategy. The ability to identify distinct mutation subsets via endoscopic ultrasound fine needle aspiration (EUS FNA) molecular cytology could greatly aid clinical trial patient stratification and offer predictive markers. We identified chemotherapy treatment naive ampullary adenocarcinoma and PDAC patients who underwent EUS FNA to assess multigene mutational frequency and diversity with a surgical resection concordance assessment, where available. Methods: Following strict cytology smear screening criteria, targeted next generation sequencing (NGS) using a 160 cancer gene panel was performed. Results: Complete sequencing was achieved in 29 patients, whereby 83 pathogenic alterations were identified in 21 genes. Cytology genotyping revealed that the majority of mutations were identified in KRAS (93%), TP53 (72%), SMAD4 (31%), and GNAS (10%). There was 100% concordance for the following pathogenic alterations: KRAS, TP53, SMAD4, KMT2D, NOTCH2, MSH2, RB1, SMARCA4, PPP2R1A, PIK3R1, SCL7A8, ATM, and FANCD2. Absolute multigene mutational concordance was 83%. Incremental cytology smear mutations in GRIN2A, GATA3 and KDM6A were identified despite re-examination of raw sequence reads in the corresponding resection specimens. Conclusions: EUS FNA cytology genotyping using a 160 cancer gene NGS panel revealed a broad spectrum of pathogenic alterations. The fidelity of cytology genotyping to that of paired surgical resection specimens suggests that EUS FNA represents a suitable surrogate and may complement the conventional stratification criteria in decision making for therapies and may guide future biomarker driven therapeutic development.

Phosphorus, in its orthophosphate form (Pi), is one of the most limiting macronutrients in soils for plant growth and development. However, the whole-genome molecular mechanisms contributing to plant acclimation to Pi deficiency remain largely unknown. White lupin (Lupinus albus) has evolved unique adaptations for growth in Pi-deficient soils, including the development of cluster roots to increase root surface area. In this study, we utilized RNA-Seq technology to assess global gene expression in white lupin cluster roots, normal roots, and leaves in response to Pi supply. We de novo assembled 277,224,180 Illumina reads from 12 complementary DNA libraries to build what is to our knowledge the first white lupin gene index (LAGI 1.0). This index contains 125,821 unique sequences with an average length of 1,155 bp. Of these sequences, 50,734 were transcriptionally active (reads per kilobase per million reads ≥ 3), representing approximately 7.8% of the white lupin genome, using the predicted genome size of Lupinus angustifolius as a reference. We identified a total of 2,128 sequences differentially expressed in response to Pi deficiency with a 2-fold or greater change and P ≤ 0.05. Twelve sequences were consistently differentially expressed due to Pi deficiency stress in three species, Arabidopsis (Arabidopsis thaliana), potato (Solanum tuberosum), and white lupin, making them ideal candidates to monitor the Pi status of plants. Additionally, classic physiological experiments were coupled with RNA-Seq data to examine the role of cytokinin and gibberellic acid in Pi deficiency-induced cluster root development. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to Pi deficiency.

In sexually reproducing animals, oocytes arrest at diplotene or diakinesis and resume meiosis (meiotic maturation) in response to hormones. In Caenorhabditis elegans, major sperm protein triggers meiotic resumption through a mechanism involving somatic Gαs–adenylate cyclase signaling and soma-to-germline gap-junctional communication. Using genetic mosaic analysis, we show that the major effector of Gαs–adenylate cyclase signaling, protein kinase A (PKA), is required in gonadal sheath cells for oocyte meiotic maturation and dispensable in the germ line. This result rules out a model in which cyclic nucleotides must transit through sheath-oocyte gap junctions to activate PKA in the germ line, as proposed in vertebrate systems. We conducted a genetic screen to identify regulators of oocyte meiotic maturation functioning downstream of Gαs–adenylate cyclase–PKA signaling. We molecularly identified 10 regulatory loci, which include essential and nonessential factors. sacy-1, which encodes a highly conserved DEAD-box helicase, is an essential germline factor that negatively regulates meiotic maturation. SACY-1 is a multifunctional protein that establishes a mechanistic link connecting the somatic control of meiotic maturation to germline sex determination and gamete maintenance. Modulatory factors include multiple subunits of a CoREST-like complex and the TWK-1 two-pore potassium channel. These factors are not absolutely required for meiotic maturation or its negative regulation in the absence of sperm, but function cumulatively to enable somatic control of meiotic maturation. This work provides insights into the genetic control of meiotic maturation signaling in C. elegans, and the conserved factors identified here might inform analysis in other systems through either homology or analogy.

The human salivary proteome is extremely complex, including proteins from salivary glands, serum, and oral microbes. Much has been learned about the host component, but little is known about the microbial component. Here we report a metaproteomic analysis of salivary supernatant pooled from six healthy subjects. For deep interrogation of the salivary proteome, we combined protein dynamic range compression (DRC), multidimensional peptide fractionation, and high-mass accuracy MS/MS with a novel two-step peptide identification method using a database of human proteins plus those translated from oral microbe genomes. Peptides were identified from 124 microbial species as well as uncultured phylotypes such as TM7. Streptococcus, Rothia, Actinomyces, Prevotella, Neisseria, Veilonella, Lactobacillus, Selenomonas, Pseudomonas, Staphylococcus, and Campylobacter were abundant among the 65 genera from 12 phyla represented. Taxonomic diversity in our study was broadly consistent with metagenomic studies of saliva. Proteins mapped to 20 KEGG pathways, with carbohydrate metabolism, amino acid metabolism, energy metabolism, translation, membrane transport, and signal transduction most represented. The communities sampled appear to be actively engaged in glycolysis and protein synthesis. This first deep metaproteomic catalog from human salivary supernatant provides a baseline for future studies of shifts in microbial diversity and protein activities potentially associated with oral disease.

In this study, we compared the genomic integration efficiencies and transposition site preferences of Sleeping Beauty (SB or SB11), Tol2, and piggyBac (PB) transposon systems in primary T cells derived from peripheral blood lymphocytes (PBL) and umbilical cord blood (UCB). We found that PB demonstrated the highest efficiency of stable gene transfer in PBL-derived T cells, whereas SB11 and Tol2 mediated intermediate and lowest efficiencies, respectively. Southern hybridization analysis demonstrated that PB generated the highest number of integrants when compared to SB and Tol2 in both PBL and UCB T cells. Tol2 and PB appeared more likely to promote clonal expansion than SB, which may be in part due to the dysregulated expression of cancer-related genes near the insertion sites. Genome-wide integration analysis demonstrated that SB, Tol2, and PB integrations occurred in all the chromosomes without preference. Additionally, Tol2 and PB integration sites were mainly localized near transcriptional start sites (TSSs), CpG islands and DNaseI hypersensitive sites, whereas SB integrations were randomly distributed. These results suggest that SB may be a preferential choice of the delivery vector in T cells due to its random integration site preference and relatively high efficiency, and support continuing development of SB-mediated T-cell phase I trials.

Connexin (Cx) 26, a major gap junction protein expressed in mammary epithelial cells, has been considered to be a tumor suppressor gene candidate. This study investigated the molecular mechanism of transcriptional up-regulation of Cx26 by phorbol ester (TPA) in human immortalized MCF-10 mammary epithelial cells and MDA-MB-231 mammary cancer cells. Such up-regulation was mediated through the protein kinase C pathway and could be blocked by the PKC inhibitor, calphostin C. Based on the results of the nuclear run-on assay, there was a TPA-induced increase in the rate of transcriptional initiation. We identified a TPA-induced DNase I hypersensitivity (DH) region approximately 1 kb 5′ upstream of the ATG translation starting site. Sequence analysis revealed that this DH region was located in intron 1 and contained two TRE-like TGAT/ATCA elements, two 5′TTCA3′ motifs and a 5′AGGAAG3′ PEA3 motif. Both TRE-like elements were capable of binding AP1. TPA inducibility of this DH region was seen by the CAT reporter assay and appeared to be direction-dependent suggesting a functional cooperation between PEA3/TTCA and TRE.

A yellow-green alga

Antimicrobials have been used extensively as growth promoters (AGPs) in agricultural animal production. However, the specific mechanism of action for AGPs has not yet been determined. The work presented here was to determine and characterize the microbiome of pigs receiving one AGP, tylosin, compared with untreated pigs. We hypothesized that AGPs exerted their growth promoting effect by altering gut microbial population composition. We determined the fecal microbiome of pigs receiving tylosin compared with untreated pigs using pyrosequencing of 16S rRNA gene libraries. The data showed microbial population shifts representing both microbial succession and changes in response to the use of tylosin. Quantitative and qualitative analyses of sequences showed that tylosin caused microbial population shifts in both abundant and less abundant species. Our results established a baseline upon which mechanisms of AGPs in regulation of health and growth of animals can be investigated. Furthermore, the data will aid in the identification of alternative strategies to improve animal health and consequently production.

Background Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling. Results Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes. Conclusions Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

Background The nutritional and economic value of many crops is effectively a function of seed protein and oil content. Insight into the genetic and molecular control mechanisms involved in the deposition of these constituents in the developing seed is needed to guide crop improvement. A quantitative trait locus (QTL) on Linkage Group I (LG I) of soybean (Glycine max (L.) Merrill) has a striking effect on seed protein content. Results A soybean near-isogenic line (NIL) pair contrasting in seed protein and differing in an introgressed genomic segment containing the LG I protein QTL was used as a resource to demarcate the QTL region and to study variation in transcript abundance in developing seed. The LG I QTL region was delineated to less than 8.4 Mbp of genomic sequence on chromosome 20. Using Affymetrix® Soy GeneChip and high-throughput Illumina® whole transcriptome sequencing platforms, 13 genes displaying significant seed transcript accumulation differences between NILs were identified that mapped to the 8.4 Mbp LG I protein QTL region. Conclusions This study identifies gene candidates at the LG I protein QTL for potential involvement in the regulation of protein content in the soybean seed. The results demonstrate the power of complementary approaches to characterize contrasting NILs and provide genome-wide transcriptome insight towards understanding seed biology and the soybean genome.

Previously, we reported a rat S1 protein that is antigenically related to statin, a nonproliferating cell-specific marker; however, it shares high homology with the known human elongation factor-1 alpha (EF-1 alpha). To differentiate S1 from rat EF-1 alpha and to study their respective regulation for expression, a rat EF-1 alpha cDNA clone was isolated and characterized. The nucleotide and deduced amino acid sequences of this partial rat EF-1 alpha cDNA are compared with that of human and mouse as well as with rat S1. Both their messages were detected in rat brain by EF-1 alpha- or S1-specific probes. However, the mRNA encoding EF-1 alpha is more abundant than that encoding S1. S1 and EF-1 alpha expression were investigated in the parotid and submandibular glands of untreated rats and those treated with isoproterenol, a proliferation-inducing catecholamine. Quantitative solution hybridization demonstrated a dramatic reduction (approximately 68 ) in the S1 mRNA following isoproterenol injection in proliferation-responsive parotid glands and a mild reduction (approximately 20 ) of S1 steady-state messages in the proliferation-refractile submandibular glands. A slight increase or no changes of EF-1 alpha levels in both parotid and submandibular glands following isoproterenol treatment are also observed. Therefore, the EF-1 alpha and S1 genes are different genes, both expressed and regulated in vivo, but in differential quantitative and qualitative patterns.

Background The rapidly expanding availability of de novo sequencing technologies can greatly facilitate efforts to monitor the relatively high mutation rates of influenza A viruses and the detection of quasispecies. Both the mutation rates and the lineages of influenza A viruses are likely to play an important role in the natural history of these viruses and the emergence of phenotypically and antigenically distinct strains. Methodology and Principal Findings We evaluated quasispecies and mixed infections by de novo sequencing the whole genomes of 10 virus isolates, including eight avian influenza viruses grown in embryonated chicken eggs (six waterfowl isolates - five H3N2 and one H4N6; an H7N3 turkey isolate; and a bald eagle isolate with H1N1/H2N1 mixed infection), and two tissue cultured H3N2 swine influenza viruses. Two waterfowl cloacal swabs were included in the analysis. Full-length sequences of all segments were obtained with 20 to 787-X coverage for the ten viruses and one cloacal swab. The second cloacal swab yielded 15 influenza reads of ∼230 bases, sufficient for bioinformatic inference of mixed infections or quasispecies. Genomic subpopulations or quasispecies of viruses were identified in four egg grown avian influenza isolates and one cell cultured swine virus. A bald eagle isolate and the second cloacal swab showed evidence of mixed infections with two (H1 and H2) and three (H1, H3, and H4) HA subtypes, respectively. Multiple sequence differences were identified between cloacal swab and the virus recovered using embryonated chicken eggs. Conclusions We describe a new approach to comprehensively identify mixed infections and quasispecies in low passage influenza A isolates and cloacal swabs and add to the understanding of the ecology of influenza A virus populations.

The mRNA and protein expressions of connexin 26 (Cx26) in rat mammary gland and uterus can be up-regulated during pregnancy as well as by the administration of human CG (hCG). In the present study, we found that the time course and magnitude of Cx26 induction by hCG was different in these two tissues. The molecular mechanism underscoring this difference was therefore investigated. We had previously demonstrated that both Sp1 and Sp3 transcription factors play a functional role in Cx26 expression. By the electrophoretic mobility shift assay, nuclear extracts from both virgin mammary gland and uterus were capable of binding to a labeled oligonucleotide probe that contained the proximal GC box and formed three protein-DNA complexes (C1, C2, and C3). In the mammary gland, pregnancy enhanced the intensity of all three complexes, whereas in the uterine tissue there was a decrease in the C2 and C3 complexes and an emergence of a new major component, C4 complex. In the supershift study, the C1 complex could be supershifted only by an antibody against Sp1, whereas C2, C3, and C4 could all be supershifted by an antibody against Sp3, suggesting a potential presence of Sp3 isoforms of various sizes. We therefore conclude that the basal Sp profiles in virgin mammary gland and uterine tissue are similar. However, in response to pregnancy, the changes in Sp profile are tissue specific and may account for the temporal and quantitative differences between these two tissues in Cx26 induction.

Human connexin 26 (Cx26) has been considered to be a candidate suppressor gene in mammary epithelial cells. To gain insight into the transcriptional regulation of this gene, we have cloned and sequenced the 5′ portion of the gene, which extends 4.8 kb upstream from the ATG translation start site. The 3′ end of the non-coding exon 1 (160 bp) is located at 3149 bp upstream from the 5′ end of exon 2. Comparison between the human Cx26 gene and the mouse gene reveals a highly conserved promoter region with 81% homology. In addition to six GC boxes and two GT boxes, a TTAAAA box is located at −24 to −19 bp upstream of the transcription start point. Analogous to the mouse β-casein gene, the promoter region of the human Cx26 gene also contains a YY1-like binding site and a consensus mammary gland factor binding site.

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP–encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence., Author Summary The filamentous ascomycete genus Cochliobolus includes highly aggressive necrotrophic and hemibiotrophic pathogens with particular specificity to their host plants, often associated with production of host selective toxins (HST) that allow necrotrophs to trigger host cell death. Hemibiotrophs must keep their hosts alive during initial infection stages and rely on subverting host defenses by secreting small protein effectors. Many Cochliobolus species have emerged rapidly as devastating pathogens due to HSTs. The genomes of Cochliobolus and related pathogens that differ in host preference, host specificity, and virulence strategies have been sequenced. Our comparative results, at the whole-genome level, and with a spotlight on core genes for secondary metabolism and small secreted proteins, touch on how pathogens develop and hone these tools, according to host or lifestyle. We suggest that, while necrotrophs and hemibiotrophs employ fundamentally contrasting mechanisms of promoting disease, the tools they utilize (HSTs and protein effectors) overlap. The suites of secondary metabolite and SSP genes that each possesses reflect astounding diversity among species, hinting that gene products, particularly those associated with unique genomic regions, are candidates for pathogenic lifestyle differences. Manipulations of strain-unique secondary metabolite genes associated with host-specific virulence provide tangible examples.

Proline-rich proteins (PRPs) are selectively expressed in the acinar cells of the salivary glands and are inducible by beta-agonist isoproterenol and dietary tannins. In the previous studies of rat PRP gene, R15, the 5'-flanking region up to -1.7 kilobase pairs (kb), which was thought to contain the necessary proximal regulatory elements, failed to confer the catecholamine isoproterenol and dietary tannin inducibility to the transgene expression in the salivary glands of transgenic mice. Here we analyzed distal 5'-flanking region of R15 in order to understand the mechanisms of tissue-specific and inducible gene regulation. An upstream regulatory region located between -2.4 and -1.7 kb of the R15 5'-flanking region is demonstrated to be indispensable for the salivary-specific and inducible reporter gene expression in vivo, by transgenic approach. Element(s) within the 0.7-kb (-2.4 to -1.7) region that is able to cis-activate the expression of a heterologous reporter gene expression is further elucidated by transient transfection assays in vitro. Three distinct nuclear orphan receptor NGFI-B regulatory sequences are identified within a 184-base pair (bp) minimal control region extended from -1995 to -1812 nucleotides relative to the transcription start site. When reporter gene containing this 184-bp control region and heterologous promoter was cotransfected with the NGFI-B expression construct, a transactivation that mimics the effect of cAMP is observed in the parotid cells. Finally, mutations on all three identified NGFI-B binding sites and coexpression of a dominant negative mutant construct, pCMV-NGFI-B(Delta25-195), abolish this transactivation mediated by NGFI-B. In summary, these data suggest that the inducible nuclear orphan receptor NGFI-B may participate in the regulation of salivary acinar cell-specific and inducible expression of the rat R15 gene via three distinct distal NGFI-B sites.

Background Next generation sequencing provides new opportunities to explore transcriptomes. However, challenges remain for accurate differentiation of homoeoalleles and paralogs, particularly in polyploid organisms with no supporting genome sequence. In this study, RNA-Seq was employed to generate and characterize the first gene expression atlas for hexaploid oat. Results The software packages Trinity and Oases were used to produce a transcript assembly from nearly 134 million 100-bp paired-end reads from developing oat seeds. Based on the quality-parameters employed, Oases assemblies were superior. The Oases 67-kmer assembly, denoted dnOST (de novo Oat Seed Transcriptome), is over 55 million nucleotides in length and the average transcript length is 1,043 nucleotides. The 74.8× sequencing depth was adequate to differentiate a large proportion of putative homoeoalleles and paralogs. To assess the robustness of dnOST, we successfully identified gene transcripts associated with the biosynthetic pathways of three compounds with health-promoting properties (avenanthramides, tocols, β-glucans), and quantified their expression. Conclusions To our knowledge, this study provides the first direct performance comparison between two major assemblers in a polyploid organism. The workflow we developed provides a useful guide for comparable analyses in other organisms. The transcript assembly developed here is a major advance. It expands the number of oat ESTs 3-fold, and constitutes the first comprehensive transcriptome study in oat. This resource will be a useful new tool both for analysis of genes relevant to nutritional enhancement of oat, and for improvement of this crop in general.

Background The late blight pathogen Phytophthora infestans can attack both potato foliage and tubers. Although interaction transcriptome dynamics between potato foliage and various pathogens have been reported, no transcriptome study has focused specifically upon how potato tubers respond to pathogen infection. When inoculated with P. infestans, tubers of nontransformed ‘Russet Burbank’ (WT) potato develop late blight disease while those of transgenic ‘Russet Burbank’ line SP2211 (+RB), which expresses the potato late blight resistance gene RB (Rpi-blb1), do not. We compared transcriptome responses to P. infestans inoculation in tubers of these two lines. Results We demonstrated the practicality of RNA-seq to study tetraploid potato and present the first RNA-seq study of potato tuber diseases. A total of 483 million paired end Illumina RNA-seq reads were generated, representing the transcription of around 30,000 potato genes. Differentially expressed genes, gene groups and ontology bins that exhibited differences between the WT and +RB lines were identified. P. infestans transcripts, including those of known effectors, were also identified. Conclusion Faster and stronger activation of defense related genes, gene groups and ontology bins correlate with successful tuber resistance against P. infestans. Our results suggest that the hypersensitive response is likely a general form of resistance against the hemibiotrophic P. infestans—even in potato tubers, organs that develop below ground.