Entire Set Of Genes Research Articles

Shorter Androgen Receptor PolyQ Alleles Protect Against Life-Threatening COVID-19 Disease in Males

Background: COVID-19 presentation ranges from asymptomatic to fatal. The variability in severity is due in part to specific mutations in the host genome. GWAS effectively identifies genetic variability due to common biallelic polymorphisms. Efforts in genetic research are trying to identify significant associations in patients infected by SARS-CoV-2. Methods: We developed a synthetic approach to genetic data representation using machine learning methods to investigate complementary genetic variability in COVID-19 infected patients that might explain disease severity due to rare variants and poly-amino acids repeat polymorphisms. Using host whole-exome sequencing data, we compared extreme phenotypic presentations of an Italian cohort of 939 subjects infected with SARS-CoV-2. We then applied the LASSO Logistic Regression model on Boolean gene-based representation of the entire set of human genes. Findings: Polymorphisms/rare variants in certain genes, including short polyQ (≤22) of the androgen receptor ( AR ), conferred protection against severe forms of COVID-19. We then demonstrated that testosterone was higher in males with AR long-polyQ (≥23), confirming receptor resistance (p=0.004 Mann-Whitney U test). Finally, long-polyQ (≥23) correlates with increased inflammation markers (p=0.021) and fibrinogen consumption (p=0.039), confirming the anti-inflammatory role of testosterone. Interpretation: Our results contribute to designing reliable clinical and public health measures and provide a rationale to test testosterone supplementation as adjuvant treatment in symptomatic COVID-19 men expressing AR polyQ longer than 23. Funding: MIUR project “Dipartimenti di Eccellenza 2018-2020” to Department of Medical Biotechnologies University of Siena, Italy; Private donors for COVID research (Italian D.L. n.18 March 17, 2020).Declaration of Interests: The authors declare no competing interests.Ethics Approval Statement: The GEN-COVID study was approved by the University Hospital of Siena Ethics Review Board (Protocol n. 16917, dated March 16, 2020).

USP32 confers cancer cell resistance to YM155 via promoting ER-associated degradation of solute carrier protein SLC35F2.

Background: The most commonly preferred chemotherapeutic agents to treat cancers are small-molecule drugs. However, the differential sensitivity of various cancer cells to small molecules and untargeted delivery narrow the range of potential therapeutic applications. The mechanisms responsible for drug resistance in a variety of cancer cells are also largely unknown. Several deubiquitinating enzymes (DUBs) are the main determinants of drug resistance in cancer cells.Methods: We used CRISPR-Cas9 to perform genome-scale knockout of the entire set of genes encoding ubiquitin-specific proteases (USPs) and systematically screened for DUBs resistant to the clinically evaluated anticancer compound YM155. A series of in vitro and in vivo experiments were conducted to reveal the relationship between USP32 and SLC35F2 on YM155-mediated DNA damage in cancer cells.Results: CRISPR-based dual-screening method identified USP32 as a novel DUB that governs resistance for uptake of YM155 by destabilizing protein levels of SLC35F2, a solute-carrier protein essential for the uptake of YM155. The expression of USP32 and SLC35F2 was negatively correlated across a panel of tested cancer cell lines. YM155-resistant cancer cells in particular exhibited elevated expression of USP32 and low expression of SLC35F2.Conclusion: Collectively, our DUB-screening strategy revealed a resistance mechanism governed by USP32 associated with YM155 resistance in breast cancers, one that presents an attractive molecular target for anti-cancer therapies. Targeted genome knockout verified that USP32 is the main determinant of SLC35F2 protein stability in vitro and in vivo, suggesting a novel way to treat tumors resistant to small-molecule drugs.

Prometheus, an omics portal for interkingdom comparative genomic analyses.

Functional analyses of genes are crucial for unveiling biological responses, genetic engineering, and developing new medicines. However, functional analyses have largely been restricted to model organisms, representing a major hurdle for functional studies and industrial applications. To resolve this, comparative genome analyses can be used to provide clues to gene functions as well as their evolutionary history. To this end, we present Prometheus, a web-based omics portal that contains more than 17,215 sequences from prokaryotic and eukaryotic genomes. This portal supports interkingdom comparative analyses via a domain architecture-based gene identification system and Gene Search, and users can easily and rapidly identify single or entire gene sets in specific pathways. Bioinformatics tools for further analyses are provided in Prometheus or through Bio-Express, a cloud-based bioinformatics analysis platform. Prometheus is a new paradigm for comparative analyses of large amounts of genomic information.

Petal abscission in fragrant roses is associated with large scale differential regulation of the abscission zone transcriptome

Flowers of fragrant roses such as Rosa bourboniana are ethylene-sensitive and undergo rapid petal abscission while hybrid roses show reduced ethylene sensitivity and delayed abscission. To understand the molecular mechanism underlying these differences, a comparative transcriptome of petal abscission zones (AZ) of 0 h and 8 h ethylene-treated flowers from R. bourboniana was performed. Differential regulation of 3700 genes (1518 up, 2182 down) representing 8.5% of the AZ transcriptome was observed between 0 and 8 h ethylene-treated R. bourboniana petal AZ. Abscission was associated with large scale up-regulation of the ethylene pathway but prominent suppression of the JA, auxin and light-regulated pathways. Regulatory genes encoding kinases/phosphatases/F-box proteins and transcription factors formed the major group undergoing differential regulation besides genes for transporters, wall modification, defense and phenylpropanoid pathways. Further comparisons with ethylene-treated petals of R. bourboniana and 8 h ethylene-treated AZ (R. hybrida) identified a core set of 255 genes uniquely regulated by ethylene in R. bourboniana AZ. Almost 23% of these encoded regulatory proteins largely conserved with Arabidopsis AZ components. Most of these were up-regulated while an entire set of photosystem genes was prominently down-regulated. The studies provide important information on regulation of petal abscission in roses.

Spatial Metagenomics of Three Geothermal Sites in Pisciarelli Hot Spring Focusing on the Biochemical Resources of the Microbial Consortia.

Terrestrial hot springs are of great interest to the general public and to scientists alike due to their unique and extreme conditions. These have been sought out by geochemists, astrobiologists, and microbiologists around the globe who are interested in their chemical properties, which provide a strong selective pressure on local microorganisms. Drivers of microbial community composition in these springs include temperature, pH, in-situ chemistry, and biogeography. Microbes in these communities have evolved strategies to thrive in these conditions by converting hot spring chemicals and organic matter into cellular energy. Following our previous metagenomic analysis of Pisciarelli hot springs (Naples, Italy), we report here the comparative metagenomic study of three novel sites, formed in Pisciarelli as result of recent geothermal activity. This study adds comprehensive information about phylogenetic diversity within Pisciarelli hot springs by peeking into possible mechanisms of adaptation to biogeochemical cycles, and high applicative potential of the entire set of genes involved in the carbohydrate metabolism in this environment (CAZome). This site is an excellent model for the study of biodiversity on Earth and biosignature identification, and for the study of the origin and limits of life.

Sequencing and Analysis of the Complete Organellar Genomes of Prototheca wickerhamii.

Of the Prototheca genus, Prototheca wickerhamii has the highest clinical significance in humans. However, neither nuclear nor organellar genomes of this species were sequenced until now. The hitherto determined and analyzed mitochondrial and plastid genomes of the alleged P. wickerhamii species belong in fact to another species, recently named Prototheca xanthoriae. This study provides a first insight into the organellar genomes of a true P. wickerhamii (type strain ATCC 16529). The P. wickerhamii mitochondrion had a 53.8-kb genome, which was considerably larger than that of Prototheca ciferrii (formerly Prototheca zopfii gen. 1) and Prototheca bovis (formerly Prototheca zopfii gen. 2), yet similarly functional, with the differences in size attributable to a higher number of introns and the presence of extra unique putative genes. The 48-kb plastid genome of P. wickerhamii, compared to autotrophic Trebouxiophyceae, was highly reduced due to the elimination of the photosynthesis-related genes. The gene content of the plastid genome of P. wickerhamii was, however, very similar to other colorless Prototheca species. Plastid genome-based phylogeny reinforced the polyphyly of the genus Prototheca, with Helicosporidium and Auxenochlorella branching within clades of Prototheca species. Phylogenetic reconstruction also confirmed the close relationship of P. wickerhamii and P. xanthoriae, which is reflected in the synteny of their organellar genomes. Interestingly, the entire set of atp genes was lost in P. wickerhamii plastid genome while being preserved in P. xanthoriae.

A GlycoGene CRISPR-Cas9 lentiviral library to study lectin binding and human glycan biosynthesis pathways.

Glycan biosynthesis on cell surface proteins and lipids is orchestrated by different classes of enzymes and proteins including the following: i. glycosyltransferases that add saccharides; ii. glycosidases that trim glycans; iii. conserved oligomeric golgi complex members that regulate intracellular transport; iv. enzymes aiding the biosynthesis of sugar-nucleotides; and v. sulfotransferases. This manuscript describes a pooled "glycoGene CRISPR" lentiviral library that targets 347 human genes involved in the above processes. Approximately 10 single-guide RNA (sgRNA) are included against each glycogene, with the putative editing site spanning the length of the target. A data analysis scheme is presented in order to determine glycosylation pathways regulating biological processes. As proof of principle, forward genetic screen results are presented to identify penetrating glycogenes that regulate the binding of P-/E-selectin, anti-sialyl Lewis-X monoclonal antibody HECA-452 and selected lectins (phaseolus vulgaris leucoagglutinin, vicia villosa lectin, peanut agglutinin) to HL-60 promyelocytic cells. Besides validating previously established biology, the study identifies three enzymes, PAPSS1, SLC35B2 and TPST2, as key molecules regulating sulfation of the major P-selectin glycoprotein ligand-1 in leukocytes. Approximately 80-90% of the sgRNA used in this study displayed high editing efficiency, and the CRISPR library picked up entire gene sets regulating specific biosynthetic pathways rather than only isolated genes. These data suggest that the glycoGene CRISPR library contains high-efficiency sgRNA. Further, this resource could be useful for the rapid screening of glycosylation-related genes and pathways that control lectin recognition in a variety of contexts.

Genome-scale screening of deubiquitinase subfamily identifies USP3 as a stabilizer of Cdc25A regulating cell cycle in cancer.

Conventional screening methods for deubiquitinating enzymes (DUBs) have important limitations. A loss-of-function study based on the knockout of DUB genes in mammalian cells can provide an excellent model for exploring DUB function. Here, we used CRISPR-Cas9 to perform genome-scale knockout of the entire set of genes encoding ubiquitin-specific proteases (USPs), a DUB subfamily, and then systematically screened for DUBs that stabilize the Cdc25A oncoprotein. USP3 was identified as a deubiquitinase of Cdc25A. USP3 depletion reduces the Cdc25A protein level, resulting in a significant delay in cell-cycle progression, and reduces the growth of cervical tumor xenografts in nude mice. Clinically, USP3 expression is positively correlated with Cdc25A protein expression and the poorest survival in breast cancer. We envision that our DUB knockout library kit will facilitate genome-scale screening of functional DUBs for target proteins of interest in a wide range of biomedical fields.

Cancer subtype classification and modeling by pathway attention and propagation.

Biological pathway is an important curated knowledge of biological processes. Thus, cancer subtype classification based on pathways will be very useful to understand differences in biological mechanisms among cancer subtypes. However, pathways include only a fraction of the entire gene set, only one-third of human genes in KEGG, and pathways are fragmented. For this reason, there are few computational methods to use pathways for cancer subtype classification. We present an explainable deep-learning model with attention mechanism and network propagation for cancer subtype classification. Each pathway is modeled by a graph convolutional network. Then, a multi-attention-based ensemble model combines several hundreds of pathways in an explainable manner. Lastly, network propagation on pathway-gene network explains why gene expression profiles in subtypes are different. In experiments with five TCGA cancer datasets, our method achieved very good classification accuracies and, additionally, identified subtype-specific pathways and biological functions. The source code is available at http://biohealth.snu.ac.kr/software/GCN_MAE. Supplementary data are available at Bioinformatics online.

The evolutionarily conserved HtrA is associated with stress tolerance and protein homeostasis in the halotolerant cyanobacterium Halothece sp. PCC7418.

The HtrA protein family represents an important class of serine proteases that are widely distributed across taxa. These evolutionarily conserved proteins are crucial for survival and function as monitors of protein synthesis during various stresses. Here, we performed gene expression analysis of the entire set of putative serine protease genes in Halothece sp. PCC7418 under salt stress conditions. The gene-encoding HtrA2 (H3553) was highly upregulated. This gene was cloned and functionally characterized, and its sub-cellular localization was determined. The recombinant H3553 protein (rH3553) displayed a pH optimum of 8.0, remained stable at 45°C, and its proteolytic activity was not affected by salts. H3553 completely degraded the unfolded model protein, β-casein. In contrast, the folded model substrates (lysozyme or BSA) were not degraded by rH3553. Denaturation of BSA at a high temperature significantly increased its degradation by rH3553. H3553 was detected in the soluble protein fraction as well as the plasma membrane and thylakoid membrane fractions. Interestingly, the majority of H3553 was present in the plasma membrane under salt and heat stress conditions. Thus, H3553 resides in multiple sub-cellular locations and its localization drastically changes after exposure to stresses. Taken together, H3553 underpins protein quality-control process and is involved in the response and adaptation to salinity and heat stresses.

Flexibility and Adaptability of Quorum Sensing in Nature.

Quorum sensing (QS), a type of chemical communication, allows bacteria to sense and coordinate activities in natural biofilm communities using N-acyl homoserine lactones (AHLs) as one type of signaling molecule. For AHL-based communication to occur, bacteria must produce and recognize the same signals, which activate similar genes in different species. Our current understanding of AHL-QS suggests that signaling between species would arise randomly, which is not probable. We propose that AHL-QS signaling is a mutable and adaptable process, within limits. AHLs are highly-conserved signals, however, their corresponding receptor proteins (LuxR) are highly variable. We suggest that both flexibility and adaptation occur among receptor proteins, allowing for complex signaling networks to develop in biofilms over time.

Freshwater Perkinsea: diversity, ecology and genomic information

AbstractStudies on freshwater Perkinsea are scarce compared to their marine counterparts; they are therefore not well ecologically characterized. In this study, we investigated the diversity, distribution and ecological role of Perkinsea in freshwater ecosystems. Our approach included (1) the phylogenetic analyses of near full-length SSU and LSU sequences of freshwater Perkinsea, (2) a meta-analysis of public Perkinsea 18S ribosomal RNA gene sequences available from the freshwater environments (25 lakes, 4 rivers), (3) microscopic observations of Perkinsea associated with planktonic communities and (4) single amplified genome analysis. Whereas Perkinsea appear to be rare in river ecosystems (85 reads), they are found in almost all of the lakes studied. However, their diversity does vary considerably between lakes (from 0 to 2 463 Operational Taxonomic Units (OTUs)). Phylogenetic analysis showed that the Parvilucifera/Dinovorax/Snorkelia and Perkinsus/Xcellia/Gadixcellia clades resulted from an initial speciation event. This second clade is further split into well-supported, monophyletic groups, including a clade dominated by freshwater representatives, which is further structured into three distinct subclades: freshwater clade 1, freshwater clade 2 and a freshwater and brackish clade. The Perkinsea Single Amplified Genome (SAG) as well as most of the abundant Operational Taxonomic Units (OTUs) fall into freshwater clade 2. The tyramide signal amplification-fluorescent in situ hybridization method showed an internal association between Perkinsea and the colonial phytoplankton Sphaerocystis. The Single Amplified Genome (SAG) annotation contained 698 genes and gene ontology terms could be assigned to 486 protein-coding genes. Although the number of genes appears to be low (10.6% of the entire gene set assessed by BUSCO), the analysis of the proteome revealed some putative secreted virulence factors. This study showed a large distribution of Perkinsea across lake ecosystems and potential parasitic association with phytoplankton. However, further investigations are needed for a better knowledge on the role of these microorganisms in freshwater ecosystems.

Normalization Methods for the Analysis of Unbalanced Transcriptome Data: A Review.

Dozens of normalization methods for correcting experimental variation and bias in high-throughput expression data have been developed during the last two decades. Up to 23 methods among them consider the skewness of expression data between sample states, which are even more than the conventional methods, such as loess and quantile. From the perspective of reference selection, we classified the normalization methods for skewed expression data into three categories, data-driven reference, foreign reference, and entire gene set. We separately introduced and summarized these normalization methods designed for gene expression data with global shift between compared conditions, including both microarray and RNA-seq, based on the reference selection strategies. To our best knowledge, this is the most comprehensive review of available preprocessing algorithms for the unbalanced transcriptome data. The anatomy and summarization of these methods shed light on the understanding and appropriate application of preprocessing methods.

Collaborative Annotation Redefines Gene Sets for Crucial Phytopathogens.

Accurate and comprehensive annotation of genomic sequences underpins advances in managing plant disease. However, important plant pathogens still have incomplete and inconsistent gene sets and lack dedicated funding or teams to improve this annotation. This paper describes a collaborative approach to gene curation to address this shortcoming. In the first instance, over 40 members of the Botrytis cinerea community from eight countries, with training and infrastructural support from Ensembl Fungi, used the gene editing tool Apollo to systematically review the entire gene set (11,707 protein coding genes) in 6–7 months. This has subsequently been checked and disseminated. Following this, a similar project for another pathogen, Blumeria graminis f. sp. hordei, also led to a completely redefined gene set. Currently, we are working with the Zymoseptoria tritici community to enable them to achieve the same. While the tangible outcome of these projects is improved gene sets, it is apparent that the inherent agreement and ownership of a single gene set by research teams as they undergo this curation process are consequential to the acceleration of research in the field. With the generation of large data sets increasingly affordable, there is value in unifying both the divergent data sets and their associated research teams, pooling time, expertise, and resources. Community-driven annotation efforts can pave the way for a new kind of collaboration among pathogen research communities to generate well-annotated reference data sets, beneficial not just for the genome being examined but for related species and the refinement of automatic gene prediction tools.

Abstract 218: Pathophysiological Significance of Browning of Perivascular Adipose Tissue in the Development of Atherosclerosis

Background: Perivascular adipose tissue (PVAT) possibly plays a pivotal role in the development of atherosclerosis through its direct local action to the vessels. However, its physiological and molecular mechanisms are less understood. Visceral and subcutaneous white adipose tissues exert a brown adipose tissue-like phenotype under specific conditions such as cold exposure, i.e. browning. We here investigated the pathophysiological role of browning of PVAT against the development of atherosclerosis after endovascular injury. Methods and Results: Endovascular injury was generated by wire insertion into the femoral artery of C57BL/6 female mice. Transcriptome analysis revealed robust upregulation of brown adipose tissue markers such as Ucp1 , Elovl3 , Cox8b and Cidea in injured arteries. Notably, administration of an atheroprotective agent, 17-beta estradiol, significantly inhibited these changes; Ucp1 was the most down-regulated gene by 17-beta estradiol administration in the entire gene set (log2-fold change = -6.58, false discovery rate = 0.003). Consistently, the upregulation of browning markers after endovascular injury and these inhibitions by 17-beta estradiol were confirmed in PVAT by quantitative real-time PCR, western blot and immunohistochemical staining. The present study also demonstrated that vascular injury promotes macrophage infiltration in PVAT accompanied with upregulation of inflammatory cytokine production. Furthermore, we confirmed spatiotemporal synchronicity between browning and inflammation in PVAT by immunohistochemical staining. Conclusions: We observed that endovascular injury elicits browning in PVAT accompanied with exacerbated inflammation, and that atheroprotective 17-beta estradiol strongly inhibits this phenomenon. These findings may provide novel insights into pathogenesis of the atherosclerosis.

Expression of socially sensitive genes: The multi-ethnic study of atherosclerosis.

BackgroundGene expression may be an important biological mediator in associations between social factors and health. However, previous studies were limited by small sample sizes and use of differing cell types with heterogeneous expression patterns. We use a large population-based cohort with gene expression measured solely in monocytes to investigate associations between seven social factors and expression of genes previously found to be sensitive to social factors.MethodsWe employ three methodological approaches: 1) omnibus test for the entire gene set (Global ANCOVA), 2) assessment of each association individually (linear regression), and 3) machine learning method that performs variable selection with correlated predictors (elastic net).ResultsIn global analyses, significant associations with the a priori defined socially sensitive gene set were detected for major or lifetime discrimination and chronic burden (p = 0.019 and p = 0.047, respectively). Marginally significant associations were detected for loneliness and adult socioeconomic status (p = 0.066, p = 0.093, respectively). No associations were significant in linear regression analyses after accounting for multiple testing. However, a small percentage of gene expressions (up to 11%) were associated with at least one social factor using elastic net.ConclusionThe Global ANCOVA and elastic net findings suggest that a small percentage of genes may be “socially sensitive,” (i.e. demonstrate differential expression by social factor), yet single gene approaches such as linear regression may be ill powered to capture this relationship. Future research should further investigate the biological mechanisms through which social factors act to influence gene expression and how systemic changes in gene expression affect overall health.

Abstract P3-08-10: A unique coding and non-coding benign breast transcriptome in post-menopausal ER+ breast cancer

Abstract Background: Differences in ER+ and ER- breast cancer tumor biology are well-documented, but little is known about the field of background benign breast in which those cancers arise. We evaluated the transcriptome of benign breast tissues from women with concurrent ipsilateral ER+ and ER- breast cancers (BC) to characterize coding and non-coding RNA profiles. This pilot study provides insight into the transcriptomic landscapes of benign breast tissues in patients with BC, and the microenvironment of the at-risk benign tissue. Methods: With institutional approval, cryobanked breast tissues from patients with concurrent ipsilateral ER+ BC (benign ER+BC, N=14) or ER- BC (benign ER-BC, N=10), were selected for benign tissues with similar epithelial:stromal ratios and grouped into pre (preM) and post-menopausal (PM) groups. Following RNA sequencing (Illumina TruSeq Stranded mRNA kit & Illumina HiSeq 4000), reads were processed (MAP-RSeq v3.0.0) and aligned (STAR aligner; hg38). Differential expression (DE) analysis (edgeR 2.6.2) identified DE genes from normalized RPKM reads (absolute log2 fold change (FC) > 1 and false discovery rate (FDR) < 0.10), corrected for intra-group biases. Over-representation analysis [Ingenuity pathway analysis (IPA), Ingenuity® Systems] and gene set enrichment analysis [(GSEA), GeneTrail 2.0] identified significantly-enriched pathways. Results: In the PM group, there were 144 DE transcripts between benign ER+BC and benign ER-BC, including coding RNAs (40%), antisense RNAs (35%) and lncRNAs (7%). In contrast, the preM group had no significantly DE genes between benign ER+BC and benign ER-BC. In the PM DE coding gene set, the top DE transcripts in benign ER+BC included many genes implicated in BC development# or ER+ BC progression (*) (e.g. up-regulated: KAAG1*, DNAJB7/HSP40*, TMEM151B, ZBTB32# p < 0.001 ; down-regulated: CPB1*, FOS, TPPP3#, CLEC3B#p < 0.001). Top canonical pathways altered in benign ER+BC included MAPK, PI3K, and acute phase response pathways (p<0.05). GSEA of the entire gene set (N=15,223; ranked in order of 144 DE genes) identified 72 altered pathways (P < 0.005); those with the highest normalization enrichment scores (NES) (> 0.4) functionally grouped as immune function-related (T cell function and antigen presentation). Depleted pathways with NES > 0.4, (N=6) functionally grouped into proteasome-related, fatty acid biosynthesis and mitochondrial energy metabolism. Among the non-coding DE gene set, notably, the entire DE antisense RNA gene set (N=51 transcripts) was up-regulated in benign ER+ BC compared to benign ER-BC (P< 0.001) with a subset (N=11) showing marked up-regulation (> 4 log2FC). Among the DE antisense RNAs, 70% have reported roles in carcinogenesis or BC progression (e.g. KRT7-AS, NAV2-AS2, CCDC144NL-AS1, RP11-66B24.4, HSA-MIR4454). Conclusion: The benign breast transcriptome differs between postmenopausal women with ER+ vs ER- BC, with distinctive coding and non-coding RNA signatures. In postmenopausal women with ER+ BC, benign breast expresses a unique antisense RNA set and is enriched in genes implicated in BC development or progression. These data provide insight into at-risk benign breast and facilitate identification of potential biomarkers of carcinogenesis. Citation Format: Carter JM, Nair AA, Davila JI, Heinzen EP, Hoskin TL, Winham SJ, Radisky DC, Visscher DW, Degnim AC. A unique coding and non-coding benign breast transcriptome in post-menopausal ER+ breast cancer [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-08-10.

Abstract P5-06-01: Distinctive coding and non-coding RNA profiles of pre-menopausal and post-menopausal benign breast

Abstract Background: Despite a paucity of data on the benign breast transcriptome, reduction mammoplasty tissue is commonly used as a control for RNA biomarker discovery in breast cancer (BC). We evaluated the transcriptome of benign breast to characterize pre-menopausal and post-menopausal coding and non-coding RNA profiles. These profiles will provide insight into the metabolic landscape of benign breast and inform transcriptomic studies of BC and peri-tumoral microenvironments. Methods: With institutional approval, cryobanked fresh breast tissues from reduction mammoplasties (age 15-62y; follow-up: 4-10 y with no cancer events), with histologically-confirmed benign tissues with similar epithelial: stromal tissues were grouped into 4 age sets (< 30y, N=10; 30-39y, N=11; 40-49y, N=10; ≥50y (clinically post-menopausal), N=10). Following RNA extraction, library preparation (Illumina TruSeq Stranded mRNA kit), and sequencing (Illumina HiSeq 4000), reads were processed (MAP-RSeq v3.0.0) and aligned (STAR aligner; hg38). Differential expression (DE) analysis (edgeR 2.6.2) identified DE genes from normalized RPKM reads (absolute log2 fold change (FC) > 1 and false discovery rate (FDR) < 0.10), corrected for intra-group biases using medians (absolute FC cut-off of >1.5). Over-representation analysis [Ingenuity pathway analysis (IPA), Ingenuity® Systems] and gene set enrichment analysis [(GSEA), GeneTrail 2.0] identified significantly-enriched pathways. Results: Across the 4 age sets, 561 DE genes were identified. Compared to the post-menopausal (PM) set, the number of DE genes was highest in <30 y set (N= 372) and decreased with increasing age (N= 170, age 30-39 set and N=20, age 40-49 set), generating up-regulated (PMup) and down-regulated (PMdown) PM transcriptomic profiles. The top PMdown DE genes included RANKL, WNT4, MKI67, extracellular matrix, and lactation-related genes (lactoferrin ,MUC4, MUC16, p < 0.01). Top PMdown canonical pathways were cell cycle-related (CDK1, CCNA2, CCNB2, ESPL1, TOP2A)(p< 0.001). Top PMup genes included those involved in adipogenesis, NPY1R, NPY2R, unsaturated fatty acid synthesis and eicosanoid signaling (P< 0.001). Top PMup canonical pathways included acyl co-A hydrolysis, stearate biosynthesis, acute phase response and RXR signaling (P < 0.001). GSEA of the entire gene set (N=15,466; ranked in order of 561 DE genes) identified 16 significant pathways, functionally grouped as mitochondrial energy metabolism, proteasome, and unsaturated fatty acid biosynthesis and signaling (PPAR pathway). While coding RNAs comprised 85% of DE transcripts, lncRNA comprised ˜5%. Top DE lncRNAs and precursor miRNAs included MEG3, a tumor suppressor, LIN00092, miR22 and miR1182 (p < 0.0001). Conclusions: Pre and post-menopausal benign breast tissues have distinctive transcriptomic profiles. PM benign breast has down-regulated proliferation/cell-cycle-related pathways, and up-regulated genes in mitochondrial energy pathways, lipogenesis and inflammatory pathways. Notably, many of the top DE coding and non-coding RNAs in PM benign breast have been implicated in BC progression, highlighting their value in better understanding breast carcinogenesis and the need to characterize their functional roles in normal aging and menopause. Citation Format: Carter JM, Nair AA, Davila JI, Heinzen EP, Hoskin TL, Winham SJ, Radisky DC, Degnim AC, Visscher DW. Distinctive coding and non-coding RNA profiles of pre-menopausal and post-menopausal benign breast [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P5-06-01.

The Core- and Pan-Genomic Analyses of the Genus Comamonas: From Environmental Adaptation to Potential Virulence.

Comamonas is often reported to be one of the major members of microbial communities in various natural and engineered environments. Versatile catabolic capabilities of Comamonas have been studied extensively in the last decade. In contrast, little is known about the ecological roles and adaptation of Comamonas to different environments as well as the virulence of potentially pathogenic Comamonas strains. In this study, we provide genomic insights into the potential ecological roles and virulence of Comamonas by analysing the entire gene set (pangenome) and the genes present in all genomes (core genome) using 34 genomes of 11 different Comamonas species. The analyses revealed that the metabolic pathways enabling Comamonas to acquire energy from various nutrient sources are well conserved. Genes for denitrification and ammonification are abundant in Comamonas, suggesting that Comamonas plays an important role in the nitrogen biogeochemical cycle. They also encode sophisticated redox sensory systems and diverse c-di-GMP controlling systems, allowing them to be able to effectively adjust their biofilm lifestyle to changing environments. The virulence factors in Comamonas were found to be highly species-specific. The conserved strategies used by potentially pathogenic Comamonas for surface adherence, motility control, nutrient acquisition and stress tolerance were also revealed.

Whole Genome Characterization of Orthopoxvirus (OPV) Abatino, a Zoonotic Virus Representing a Putative Novel Clade of Old World Orthopoxviruses.

Orthopoxviruses (OPVs) are diffused over the complete Eurasian continent, but previously described strains are mostly from northern Europe, and few infections have been reported from Italy. Here we present the extended genomic characterization of OPV Abatino, a novel OPV isolated in Italy from an infected Tonkean macaque, with zoonotic potential. Phylogenetic analysis based on 102 conserved OPV genes (core gene set) showed that OPV Abatino is most closely related to the Ectromelia virus species (ECTV), although placed on a separate branch of the phylogenetic tree, bringing substantial support to the hypothesis that this strain may be part of a novel OPV clade. Extending the analysis to the entire set of genes (coding sequences, CDS) further substantiated this hypothesis. In fact the genome of OPV Abatino included more CDS than ECTV; most of the extra genes (mainly located in the terminal genome regions), showed the highest similarity with cowpox virus (CPXV); however vaccinia virus (VACV) and monkeypox virus (MPXV) were the closest OPV for certain CDS. These findings suggest that OPV Abatino could be the result of complex evolutionary events, diverging from any other previously described OPV, and may indicate that previously reported cases in Italy could represent the tip of the iceberg yet to be explored.