Divergent Folding-Mediated Epistasis Among Unstable Membrane Protein Variants

Many membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for the mammalian gonadotropin-releasing hormone receptor GPCRs (GnRHR). We recently demonstrated that evolutionary GnRHR modifications appear to have coincided with adaptive changes in cotranslational folding efficiency. Though protein stability is known to shape evolution, it is unclear how cotranslational folding constraints modulate the synergistic, epistatic interactions between mutations. We therefore compared the pairwise interactions formed by mutations that disrupt the membrane topology (V276T) or tertiary structure (W107A) of GnRHR. Using deep mutational scanning, we evaluated how the plasma membrane expression of these variants is modified by hundreds of secondary mutations. An analysis of 251 mutants in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the severity and the mechanism of destabilization. V276T forms predominantly negative epistatic interactions with destabilizing mutations in soluble loops. In contrast, W107A forms positive interactions with mutations in both loops and transmembrane domains that reflect the diminishing impacts of the destabilizing mutations in variants that are already unstable. These findings reveal how epistasis is remodeled by conformational defects in membrane proteins and in unstable proteins more generally.

Transcriptomics analysis of allergen-induced inflammatory gene expression in the Four-Core Genotype mouse model

Sex differences in allergic inflammation have been reported, but the mechanisms underlying these differences remain unknown. Contributions of both sex hormones and sex-related genes to these mechanisms have been previously suggested in clinical and animal studies. Here, Four-Core Genotypes (FCG) mouse model was used to study the inflammatory response to house dust mite (HDM) challenge and identify differentially expressed genes (DEGs) and regulatory pathways in lung tissue. Briefly, adult mice (8-10 wk old) of the FCG (XXM, XXF, XYM, XYF) were challenged intranasally with 25 μg of HDM or vehicle (PBS-control group) 5 days/wk for 5 wk (n = 3/10 group). At 72 h after the last exposure, we analyzed the eosinophils and neutrophils in the bronchoalveolar lavage (BAL) of FCG mice. We extracted lung tissue and determined DEGs using Templated Oligo-Sequencing (TempO-Seq). DEG analysis was performed using the DESeq2 package and gene enrichment analysis was done using Ingenuity Pathway Analysis. A total of 2,863 DEGs were identified in the FCG. Results revealed increased eosinophilia and neutrophilia in the HDM-treated group with the most significantly expressed genes in XYF phenotype and a predominant effect of female hormones vs. chromosomes. Regardless of the sex hormones, mice with female chromosomes had more downregulated genes in the HDM group but this was reversed in the control group. Interestingly, genes associated with inflammatory responses were overrepresented in the XXM and XYF genotypes treated with HDM. Sex hormones and chromosomes contribute to inflammatory responses to HDM challenge, with female hormones exerting a predominant effect mediated by inflammatory DEGs.NEW & NOTEWORTHY Gene expression profiling helps to provide deep insight into the global view of disease-related mechanisms and responses to therapy. Using the Four-Core Genotype mouse model, our findings revealed the influence of sex hormones and sex chromosomes in the gene expression of lungs exposed to an aeroallergen (House Dust Mite) and identified sex-specific pathways to better understand sex disparities associated with allergic airway inflammation.

Lung proinflammatory microRNA and cytokine expression in a mouse model of allergic inflammation: role of sex chromosome complement and gonadal hormones

Epigenetic alterations such as dysregulation of miRNAs have been reported to play important roles in interactions between genetic and environmental factors. In this study, we tested the hypothesis that induction of lung inflammation by inhaled allergens triggers a sex-specific miRNA regulation that is dependent on chromosome complement and hormonal milieu. We challenged the four core genotypes (FCGs) model through intranasal sensitization with a house dust mite (HDM) solution (or PBS as a control) for 5 wk. The FCG model allows four combinations of gonads and sex chromosomes: 1) XX mice with ovaries (XXF), 2) XY mice with testes (XYM), 3) XX mice with testes (XXM), and 4) XY mice with ovaries (XYF). Following the challenge (n = 5-7/group), we assessed the expression of 84 inflammatory miRNAs in lung tissue using a PCR array and cytokine levels in bronchoalveolar lavage fluid (BAL) by a multiplex protein assay (n = 4-7 animals/group). Our results showed higher levels of the chemokine KC (an Il-8 homolog) and IL-7 in BAL from XYF mice challenged with HDM. In addition, IL-17A was significantly higher in BAL from both XXF and XYF mice. A three-way interaction among treatment, gonads, and sex chromosome revealed 60 of 64 miRNAs that differed in expression depending on genotype; XXF, XXM, XYF, and XYM mice had 45, 32, 4, and 52 differentially expressed miRNAs, respectively. Regulatory networks of miRNAs identified in this study were implicated in pathways associated with asthma. Female gonadal hormonal effects may alter miRNA expression and contribute to the higher susceptibility of females to asthma.NEW & NOTEWORTHY miRNAs play important roles in regulating gene and environmental interactions. However, their role in mediating sex differences in allergic responses and lung diseases has not been elucidated. Our study used a targeted omics approach to characterize the contributions of gonadal hormones and chromosomal components to lung responses to an allergen challenge. Our results point to the influence of sex hormones in miRNA expression and proinflammatory markers in allergic airway inflammation.

Ravidomycin Analogs from Streptomyces sp. Exhibit Altered Antimicrobial and Cytotoxic Selectivity

Six new ravidomycin analogs (1-4, 6, and 7) were isolated from Streptomyces sp. Am59 using UV- and LCMS-guided separation based on Global Natural Products Social (GNPS) molecular networking analysis. Furthermore, we isolated fucomycin V (9), which possesses the same chromophore as ravidomycin but features a d-fucopyranose instead of d-ravidosamine. This is the first report of 9 as a natural product. Four new analogs (10-13) of 9 were also isolated. The structures were elucidated by combined spectroscopic and computational methods. We also found an inconsistency with the published [α]D25 of deacetylravidomycin, which is reported to have a (-) sign. Instead, we observed a (+) specific rotation for the reported absolute configuration of deacetylravidomycin (containing d-ravidosamine). We confirmed the positive sign by reisolating deacetylravidomycin from S. ravidus and by deacetylating ravidomycin. Finally, antibacterial, antifungal, and cytotoxicity activities were determined for the compounds. Compared to deacetylravidomycin, the compounds 4-6, 9, 11, and 12 exhibited greater antibacterial selectivity.

Transcriptomic analyses of ovarian clear-cell carcinoma with concurrent endometriosis

Introduction

Endometriosis, a benign inflammatory disease whereby endometrial-like tissue grows outside the uterus, is a risk factor for endometriosis-associated ovarian cancers. In particular, ovarian endometriomas, cystic lesions of deeply invasive endometriosis, are considered the precursor lesion for ovarian clear-cell carcinoma (OCCC).

Methods

To explore this transcriptomic landscape, OCCC from women with pathology-proven concurrent endometriosis (n = 4) were compared to benign endometriomas (n = 4) by bulk RNA and small-RNA sequencing.

Results

Analysis of protein-coding genes identified 2449 upregulated and 3131 downregulated protein-coding genes (DESeq2, P< 0.05, log2 fold-change > |1|) in OCCC with concurrent endometriosis compared to endometriomas. Gene set enrichment analysis showed upregulation of pathways involved in cell cycle regulation and DNA replication and downregulation of pathways involved in cytokine receptor signaling and matrisome. Comparison of pathway activation scores between the clinical samples and publicly-available datasets for OCCC cell lines revealed significant molecular similarities between OCCC with concurrent endometriosis and OVTOKO, OVISE, RMG1, OVMANA, TOV21G, IGROV1, and JHOC5 cell lines. Analysis of miRNAs revealed 64 upregulated and 61 downregulated mature miRNA molecules (DESeq2, P< 0.05, log2 fold-change > |1|). MiR-10a-5p represented over 21% of the miRNA molecules in OCCC with endometriosis and was significantly upregulated (NGS: log2fold change = 4.37, P = 2.43e-18; QPCR: 8.1-fold change, P< 0.05). Correlation between miR-10a expression level in OCCC cell lines and IC50 (50% inhibitory concentration) of carboplatin in vitro revealed a positive correlation (R2 = 0.93). MiR-10a overexpression in vitro resulted in a significant decrease in proliferation (n = 6; P< 0.05) compared to transfection with a non-targeting control miRNA. Similarly, the cell-cycle analysis revealed a significant shift in cells from S and G2 to G1 (n = 6; P< 0.0001). Bioinformatic analysis predicted that miR-10a-5p target genes that were downregulated in OCCC with endometriosis were involved in receptor signaling pathways, proliferation, and cell cycle progression. MiR-10a overexpression in vitro was correlated with decreased expression of predicted miR-10a target genes critical for proliferation, cell-cycle regulation, and cell survival including [SERPINE1 (3-fold downregulated; P< 0.05), CDK6 (2.4-fold downregulated; P< 0.05), and RAP2A (2-3-fold downregulated; P< 0.05)].

Discussion

These studies in OCCC suggest that miR-10a-5p is an impactful, potentially oncogenic molecule, which warrants further studies.

Host Derivation of Sindbis Virus Influences Mammalian Type I Interferon Response to Infection

Arboviruses are defined by their ability to replicate in both mosquito vectors and mammalian hosts. There is good evidence that arboviruses "prime" their progeny for infection of the next host, such as via differential glycosylation of their outer glycoproteins or packaging of host ribosomal subunits. We and others have previously shown that mosquito-derived viruses more efficiently infect mammalian cells than mammalian-derived viruses. These observations are consistent with arboviruses acquiring host-specific adaptations, and we hypothesized that a virus derived from either the mammalian host or mosquito vector elicits different responses when infecting the mammalian host. Here, we perform an RNA-sequencing analysis of the transcriptional response of Human Embryonic Kidney 293 (HEK-293) cells to infection with either mosquito (Aedes albopictus, C7/10)- or mammalian (Baby Hamster Kidney, BHK-21)-derived Sindbis virus (SINV). We show that the C7/10-derived virus infection leads to a more robust transcriptional response in HEK-293s compared to infection with the BHK-derived virus. Surprisingly, despite more efficient infection, we found an increase in interferon-β (IFN-β) and interferon-stimulated gene (ISG) transcripts in response to the C7/10-derived virus infection versus the BHK-derived virus infection. However, translation of interferon-stimulated genes was lower in HEK-293s infected with the C7/10-derived virus, starkly contrasting with the transcriptional response. This inhibition of ISG translation is reflective of a more rapid overall shut-off of host cell translation following infection with the C7/10-derived virus. Finally, we show that the C7/10-derived virus infection of HEK-293 cells leads to elevated levels of phosphorylated eukaryotic translation elongation factor-2 (eEF2), identifying a potential mechanism leading to the more rapid shut-off of host translation. We postulate that the rapid shut-off of host translation in mammalian cells infected with the mosquito-derived virus acts to counter the IFN-β-stimulated transcriptional response.

Exploring breast tissue microbial composition and the association with breast cancer risk factors

BACKGROUND

Microbial dysbiosis has emerged as an important element in the development and progression of various cancers, including breast cancer. However, the microbial composition of the breast from healthy individuals, even relative to risk of developing breast cancer, remains unclear. Here, we performed a comprehensive analysis of the microbiota of the normal breast tissue, which was analyzed in relation to the microbial composition of the tumor and adjacent normal tissue.

METHODS

The study cohorts included 403 cancer-free women (who donated normal breast tissue cores) and 76 breast cancer patients (who donated tumor and/or adjacent normal tissue samples). Microbiome profiling was obtained by sequencing the nine hypervariable regions of the 16S rRNA gene (V1V2, V2V3, V3V4, V4V5, V5V7, and V7V9). Transcriptome analysis was also performed on 190 normal breast tissue samples. Breast cancer risk score was assessed using the Tyrer-Cuzick risk model.

RESULTS

The V1V2 amplicon sequencing resulted more suitable for the analysis of the normal breast microbiome and identified Lactobacillaceae (Firmicutes phylum), Acetobacterraceae, and Xanthomonadaceae (both Proteobacteria phylum) as the most abundant families in the normal breast. However, Ralstonia (Proteobacteria phylum) was more abundant in both breast tumors and histologically normal tissues adjacent to malignant tumors. We also conducted a correlation analysis between the microbiome and known breast cancer risk factors. Abundances of the bacterial taxa Acetotobacter aceti, Lactobacillus vini, Lactobacillus paracasei, and Xanthonomas sp. were associated with age (p < 0.0001), racial background (p < 0.0001), and parity (p < 0.0001). Finally, transcriptome analysis of normal breast tissues showed an enrichment in metabolism- and immune-related genes in the tissues with abundant Acetotobacter aceti, Lactobacillus vini, Lactobacillus paracasei, and Xanthonomas sp., whereas the presence of Ralstonia in the normal tissue was linked to dysregulation of genes involved in the carbohydrate metabolic pathway.

CONCLUSIONS

This study defines the microbial features of normal breast tissue, thus providing a basis to understand cancer-related dysbiosis. Moreover, the findings reveal that lifestyle factors can significantly affect the normal breast microbial composition.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Transcriptome-wide analysis of pseudouridylation in Drosophila melanogaster

Pseudouridine (Psi) is one of the most frequent post-transcriptional modification of RNA. Enzymatic Psi modification occurs on rRNA, snRNA, snoRNA, tRNA, and non-coding RNA and has recently been discovered on mRNA. Transcriptome-wide detection of Psi (Psi-seq) has yet to be performed for the widely studied model organism Drosophila melanogaster. Here, we optimized Psi-seq analysis for this species and have identified thousands of Psi modifications throughout the female fly head transcriptome. We find that Psi is widespread on both cellular and mitochondrial rRNAs. In addition, more than a thousand Psi sites were found on mRNAs. When pseudouridylated, mRNAs frequently had many Psi sites. Many mRNA Psi sites are present in genes encoding for ribosomal proteins, and many are found in mitochondrial encoded RNAs, further implicating the importance of pseudouridylation for ribosome and mitochondrial function. The 7SLRNA of the signal recognition particle is the non-coding RNA most enriched for Psi. The 3 mRNAs most enriched for Psi encode highly expressed yolk proteins (Yp1, Yp2, and Yp3). By comparing the pseudouridine profiles in the RluA-2 mutant and the w1118 control genotype, we identified Psi sites that were missing in the mutant RNA as potential RluA-2 targets. Finally, differential gene expression analysis of the mutant transcriptome indicates a major impact of loss of RluA-2 on the ribosome and translational machinery.

Single cell evaluation of endocardial Hand2 gene regulatory networks reveals HAND2-dependent pathways that impact cardiac morphogenesis

The transcription factor HAND2 plays essential roles during cardiogenesis. Hand2 endocardial deletion (H2CKO) results in tricuspid atresia or double inlet left ventricle with accompanying intraventricular septum defects, hypo-trabeculated ventricles and an increased density of coronary lumens. To understand the regulatory mechanisms of these phenotypes, single cell transcriptome analysis of mouse E11.5 H2CKO hearts was performed revealing a number of disrupted endocardial regulatory pathways. Using HAND2 DNA occupancy data, we identify several HAND2-dependent enhancers, including two endothelial enhancers for the shear-stress master regulator KLF2. A 1.8 kb enhancer located 50 kb upstream of the Klf2 TSS imparts specific endothelial/endocardial expression within the vasculature and endocardium. This enhancer is HAND2-dependent for ventricular endocardium expression but HAND2-independent for Klf2 vascular and valve expression. Deletion of this Klf2 enhancer results in reduced Klf2 expression within ventricular endocardium. These data reveal that HAND2 functions within endocardial gene regulatory networks including shear-stress response.

High-resolution phylogenetic and population genetic analysis of microbial communities with RoC-ITS

Microbial communities are inter-connected systems of incredible complexity and dynamism that play crucial roles in health, energy, and the environment. To better understand microbial communities and how they respond to change, it is important to know which microbes are present and their relative abundances at the greatest taxonomic resolution possible. Here, we describe a novel protocol (RoC-ITS) that uses the single-molecule Nanopore sequencing platform to assay the composition of microbial communities at the subspecies designation. Using rolling-circle amplification, this methodology produces long-read sequences from a circular construct containing the complete 16S ribosomal gene and the neighboring internally transcribed spacer (ITS). These long reads can be used to generate a high-fidelity circular consensus sequence. Generally, the ribosomal 16S gene provides phylogenetic information down to the species-level, while the much less conserved ITS region contains strain-level information. When linked together, this combination of markers allows for the identification of individual ribosomal units within a specific organism and the assessment of their relative stoichiometry, as well as the ability to monitor subtle shifts in microbial community composition with a single generic assay. We applied RoC-ITS to an artificial microbial community that was also sequenced using the Illumina platform, to assess its accuracy in quantifying the relative abundance and identity of each species.

How experimental competition changes ovarian gene activity in free-living birds: Implications for steroidogenesis, maternal effects, and beyond

The ovary plays an important role in mediating both a female's response to her social environment and communicating it to her developing offspring via maternal effects. Past work has focused on how ovarian hormones respond to competition, but we know little about how the broader ovarian transcriptomic landscape changes, either during or after competition, giving us a narrow perspective on how socially induced phenotypes arise. Here, we experimentally generated social competition among wild, cavity-nesting female birds (tree swallows, Tachycineta bicolor), a species in which females lack a socially induced rise in circulating testosterone but they nevertheless increase allocation to eggs. After territory settlement, we reduced availability of nesting cavities, generating heightened competition; within 24 h we reversed the manipulation, causing aggressive interactions to subside. We measured ovarian transcriptomic responses at the peak of competition and 48 h later, along with date-matched controls. Network analyses indicated that competing females experienced an immediate and temporary decrease in the expression of genes involved in the early stages of steroidogenesis, and this was moderately correlated with plasma testosterone; however, two days after competition had ended, there was a marked increase in the expression of genes involved in the final stages of steroidogenesis, including HSD17B1. Gene networks related to the cell cycle, muscle performance, and extracellular matrix organization also displayed altered activity. Although the functional consequences of these findings are unclear, they shed light on socially responsive ovarian genomic mechanisms that could potentially exert lasting effects on behavior, reproduction, and maternal effects.

How thermal challenges change gene regulation in the songbird brain and gonad: implications for sexual selection in our changing world

In a rapidly warming world, exposure to high temperatures may impact fitness, but the gene regulatory mechanisms that link sublethal heat to sexually selected traits are not well understood, particularly in endothermic animals. Our experiment used zebra finches (Taeniopygia guttata), songbirds that experience extreme temperature fluctuations in their native Australia. We exposed captive males to an acute thermal challenge (43°C) compared with thermoneutral (35°C) and lower (27°C) temperatures. We found significantly more heat dissipation behaviors at 43°C, a temperature previously shown to reduce song production and fertility, and more heat retention behaviors at 27°C. Next, we characterized transcriptomic responses in tissues important for mating effort - the posterior telencephalon, for its role in song production, and the testis, for its role in fertility and hormone production. Differential expression of hundreds of genes in the testes, but few in the brain, suggest the brain is less responsive to extreme temperatures. Nevertheless, gene network analyses revealed that expression related to dopaminergic signaling in the brain co-varied with heat dissipation behaviors, providing a mechanism by which temporary thermal challenges may alter motivational circuits for song production. In both brain and testis, we observed correlations between thermally sensitive gene networks and individual differences in thermoregulatory behavior. Although we cannot directly relate these gene regulatory changes to mating success, our results suggest that individual variation in response to thermal challenges could impact sexually selected traits in a warming world.

Sulfur cycling and host-virus interactions in Aquificales-dominated biofilms from Yellowstone's hottest ecosystems

Modern linkages among magmatic, geochemical, and geobiological processes provide clues about the importance of thermophiles in the origin of biogeochemical cycles. The aim of this study was to identify the primary chemoautotrophs and host-virus interactions involved in microbial colonization and biogeochemical cycling at sublacustrine, vapor-dominated vents that represent the hottest measured ecosystems in Yellowstone National Park (~140 °C). Filamentous microbial communities exposed to extreme thermal and geochemical gradients were sampled using a remotely operated vehicle and subjected to random metagenome sequencing and microscopic analyses. Sulfurihydrogenibium (phylum Aquificae) was the predominant lineage (up to 84% relative abundance) detected at vents that discharged high levels of dissolved H₂, H₂S, and CO₂. Metabolic analyses indicated carbon fixation by Sulfurihydrogenibium spp. was powered by the oxidation of reduced sulfur and H₂, which provides organic carbon for heterotrophic community members. Highly variable Sulfurihydrogenibium genomes suggested the importance of intra-population diversity under extreme environmental and viral pressures. Numerous lytic viruses (primarily unclassified taxa) were associated with diverse archaea and bacteria in the vent community. Five circular dsDNA uncultivated virus genomes (UViGs) of ~40 kbp length were linked to the Sulfurihydrogenibium metagenome-assembled genome (MAG) by CRISPR spacer matches. Four UViGs contained consistent genome architecture and formed a monophyletic cluster with the recently proposed Pyrovirus genus within the Caudovirales. Sulfurihydrogenibium spp. also contained CRISPR arrays linked to plasmid DNA with genes for a novel type IV filament system and a highly expressed β-barrel porin. A diverse suite of transcribed secretion systems was consistent with direct microscopic analyses, which revealed an extensive extracellular matrix likely critical to community structure and function. We hypothesize these attributes are fundamental to the establishment and survival of microbial communities in highly turbulent, extreme-gradient environments.

Aberrant epigenetic and transcriptional events associated with breast cancer risk

BACKGROUND

Genome-wide association studies have identified several breast cancer susceptibility loci. However, biomarkers for risk assessment are still missing. Here, we investigated cancer-related molecular changes detected in tissues from women at high risk for breast cancer prior to disease manifestation. Disease-free breast tissue cores donated by healthy women (N = 146, median age = 39 years) were processed for both methylome (MethylCap) and transcriptome (Illumina's HiSeq4000) sequencing. Analysis of tissue microarray and primary breast epithelial cells was used to confirm gene expression dysregulation.

RESULTS

Transcriptomic analysis identified 69 differentially expressed genes between women at high and those at average risk of breast cancer (Tyrer-Cuzick model) at FDR < 0.05 and fold change ≥ 2. Majority of the identified genes were involved in DNA damage checkpoint, cell cycle, and cell adhesion. Two genes, FAM83A and NEK2, were overexpressed in tissue sections (FDR < 0.01) and primary epithelial cells (p < 0.05) from high-risk breasts. Moreover, 1698 DNA methylation changes were identified in high-risk breast tissues (FDR < 0.05), partially overlapped with cancer-related signatures, and correlated with transcriptional changes (p < 0.05, r ≤ 0.5). Finally, among the participants, 35 women donated breast biopsies at two time points, and age-related molecular alterations enhanced in high-risk subjects were identified.

CONCLUSIONS

Normal breast tissue from women at high risk of breast cancer bears molecular aberrations that may contribute to breast cancer susceptibility. This study is the first molecular characterization of the true normal breast tissues, and provides an opportunity to investigate molecular markers of breast cancer risk, which may lead to new preventive approaches.

A novel transposable element-based authentication protocol for Drosophila cell lines

Drosophila cell lines are used by researchers to investigate various cell biological phenomena. It is crucial to exercise good cell culture practice. Poor handling can lead to both inter- and intra-species cross-contamination. Prolonged culturing can lead to introduction of large- and small-scale genomic changes. These factors, therefore, make it imperative that methods to authenticate Drosophila cell lines are developed to ensure reproducibility. Mammalian cell line authentication is reliant on short tandem repeat (STR) profiling; however, the relatively low STR mutation rate in Drosophila melanogaster at the individual level is likely to preclude the value of this technique. In contrast, transposable elements (TEs) are highly polymorphic among individual flies and abundant in Drosophila cell lines. Therefore, we investigated the utility of TE insertions as markers to discriminate Drosophila cell lines derived from the same or different donor genotypes, divergent sub-lines of the same cell line, and from other insect cell lines. We developed a PCR-based next-generation sequencing protocol to cluster cell lines based on the genome-wide distribution of a limited number of diagnostic TE families. We determined the distribution of five TE families in S2R+, S2-DRSC, S2-DGRC, Kc167, ML-DmBG3-c2, mbn2, CME W1 Cl.8+, and ovarian somatic sheath Drosophila cell lines. Two independent downstream analyses of the next-generation sequencing data yielded similar clustering of these cell lines. Double-blind testing of the protocol reliably identified various Drosophila cell lines. In addition, our data indicate minimal changes with respect to the genome-wide distribution of these five TE families when cells are passaged for at least 50 times. The protocol developed can accurately identify and distinguish the numerous Drosophila cell lines available to the research community, thereby aiding reproducible Drosophila cell culture research.

The Potential for a Released Autosomal X-Shredder Becoming a Driving-Y Chromosome and Invasively Suppressing Wild Populations of Malaria Mosquitoes

Sex-ratio distorters based on X-chromosome shredding are more efficient than sterile male releases for population suppression. X-shredding is a form of sex distortion that skews spermatogenesis of XY males towards the preferential transmission of Y-bearing gametes, resulting in a higher fraction of sons than daughters. Strains harboring X-shredders on autosomes were first developed in the malaria mosquito Anopheles gambiae, resulting in strong sex-ratio distortion. Since autosomal X-shredders are transmitted in a Mendelian fashion and can be selected against, their frequency in the population declines once releases are halted. However, unintended transfer of X-shredders to the Y-chromosome could produce an invasive meiotic drive element, that benefits from its biased transmission to the predominant male-biased offspring and its effective shielding from female negative selection. Indeed, linkage to the Y-chromosome of an active X-shredder instigated the development of the nuclease-based X-shredding system. Here, we analyze mechanisms whereby an autosomal X-shredder could become unintentionally Y-linked after release by evaluating the stability of an established X-shredder strain that is being considered for release, exploring its potential for remobilization in laboratory and wild-type genomes of An. gambiae and provide data regarding expression on the mosquito Y-chromosome. Our data suggest that an invasive X-shredder resulting from a post-release movement of such autosomal transgenes onto the Y-chromosome is unlikely.

S-RNase Alleles Associated With Self-Compatibility in the Tomato Clade: Structure, Origins, and Expression Plasticity

The self-incompatibility (SI) system in the Solanaceae is comprised of cytotoxic pistil S-RNases which are countered by S-locus F-box (SLF) resistance factors found in pollen. Under this barrier-resistance architecture, mating system transitions from SI to self-compatibility (SC) typically result from loss-of-function mutations in genes encoding pistil SI factors such as S-RNase. However, the nature of these mutations is often not well characterized. Here we use a combination of S-RNase sequence analysis, transcript profiling, protein expression and reproductive phenotyping to better understand different mechanisms that result in loss of S-RNase function. Our analysis focuses on 12 S-RNase alleles identified in SC species and populations across the tomato clade. In six cases, the reason for gene dysfunction due to mutations is evident. The six other alleles potentially encode functional S-RNase proteins but are typically transcriptionally silenced. We identified three S-RNase alleles which are transcriptionally silenced under some conditions but actively expressed in others. In one case, expression of the S-RNase is associated with SI. In another case, S-RNase expression does not lead to SI, but instead confers a reproductive barrier against pollen tubes from other tomato species. In the third case, expression of S-RNase does not affect self, interspecific or inter-population reproductive barriers. Our results indicate that S-RNase expression is more dynamic than previously thought, and that changes in expression can impact different reproductive barriers within or between natural populations.

A Hybrid Nanofiber/Paper Cell Culture Platform for Building a 3D Blood-brain Barrier Model

The blood brain barrier (BBB) protects the central nervous system from toxins and pathogens in the blood by regulating permeation of molecules through the barrier interface. In vitro BBB models described to date reproduce some aspects of BBB functionality, but also suffer from incomplete phenotypic expression of brain endothelial traits, difficulty in reproducibility and fabrication, or overall cost. To address these limitations, we describe a three-dimensional (3D) BBB model based on a hybrid paper/nanofiber scaffold. The cell culture platform utilizes lens paper as a framework to accommodate 3D culture of astrocytes. An electrospun nanofiber layer is coated onto one face of the paper to mimic the basement membrane and support growth of an organized two-dimensional layer of endothelial cells (ECs). Human induced pluripotent stem cell-derived ECs and astrocytes are co-cultured to develop a human BBB model. Morphological and spatial organization of model are validated with confocal microscopy. Measurements of transendothelial resistance and permeability demonstrate the BBB model develops a high-quality barrier and responds to hyperosmolar treatments. RNA-sequencing shows introduction of astrocytes both regulates EC tight junction proteins and improves endothelial phenotypes related to vasculogenesis. This model shows promise as a model platform for future in vitro studies of the BBB.

The quorum-sensing systems of Vibrio campbellii DS40M4 and BB120 are genetically and functionally distinct

Vibrio campbellii BB120 (previously classified as Vibrio harveyi) is a fundamental model strain for studying quorum sensing in vibrios. A phylogenetic evaluation of sequenced Vibrio strains in Genbank revealed that BB120 is closely related to the environmental isolate V. campbellii DS40M4. We exploited DS40M4's competence for exogenous DNA uptake to rapidly generate greater than 30 isogenic strains with deletions of genes encoding BB120 quorum-sensing system homologues. Our results show that the quorum-sensing circuit of DS40M4 is distinct from BB120 in three ways: (i) DS40M4 does not produce an acyl homoserine lactone autoinducer but encodes an active orphan LuxN receptor, (ii) the quorum regulatory small RNAs (Qrrs) are not solely regulated by autoinducer signalling through the response regulator LuxO and (iii) the DS40M4 quorum-sensing regulon is much smaller than BB120 (~100 genes vs. ~400 genes, respectively). Using comparative genomics to expand our understanding of quorum-sensing circuit diversity, we observe that conservation of LuxM/LuxN proteins differs widely both between and within Vibrio species. These strains are also phenotypically distinct: DS40M4 exhibits stronger interbacterial cell killing, whereas BB120 forms more robust biofilms and is bioluminescent. These results underscore the need to examine wild isolates for a broader view of bacterial diversity in the marine ecosystem.

ZEB2 regulates endocrine therapy sensitivity and metastasis in luminal a breast cancer cells through a non-canonical mechanism

PURPOSE

The transcription factors ZEB1 and ZEB2 mediate epithelial-to-mesenchymal transition (EMT) and metastatic progression in numerous malignancies including breast cancer. ZEB1 and ZEB2 drive EMT through transcriptional repression of cell-cell junction proteins and members of the tumor suppressive miR200 family. However, in estrogen receptor positive (ER +) breast cancer, the role of ZEB2 as an independent driver of metastasis has not been fully investigated.

METHODS

In the current study, we induced exogenous expression of ZEB2 in ER + MCF-7 and ZR-75-1 breast cancer cell lines and examined EMT gene expression and metastasis using dose-response qRT-PCR, transwell migration assays, proliferation assays with immunofluorescence of Ki-67 staining. We used RNA sequencing to identify pathways and genes affected by ZEB2 overexpression. Finally, we treated ZEB2-overexpressing cells with 17β-estradiol (E2) or ICI 182,780 to evaluate how ZEB2 affects estrogen response.

RESULTS

Contrary to expectation, we found that ZEB2 did not increase canonical epithelial nor decrease mesenchymal gene expressions. Furthermore, ZEB2 overexpression did not promote a mesenchymal cell morphology. However, ZEB1 and ZEB2 protein expression induced significant migration of MCF-7 and ZR-75-1 breast cancer cells in vitro and MCF-7 xenograft metastasis in vivo. Transcriptomic (RNA sequencing) pathway analysis revealed alterations in estrogen signaling regulators and pathways, suggesting a role for ZEB2 in endocrine sensitivity in luminal A breast cancer. Expression of ZEB2 was negatively correlated with estrogen receptor complex genes in luminal A patient tumors. Furthermore, treatment with 17β-estradiol (E2) or the estrogen receptor antagonist ICI 182,780 had no effect on growth of ZEB2-overexpressing cells.

CONCLUSION

ZEB2 is a multi-functional regulator of drug sensitivity, cell migration, and metastasis in ER + breast cancer and functions through non-canonical mechanisms.

Bacterial-Driven Inflammation and Mutant BRAF Expression Combine to Promote Murine Colon Tumorigenesis That Is Sensitive to Immune Checkpoint Therapy

Colorectal cancer is multifaceted, with subtypes defined by genetic, histologic, and immunologic features that are potentially influenced by inflammation, mutagens, and/or microbiota. Colorectal cancers with activating mutations in BRAF are associated with distinct clinical characteristics, although the pathogenesis is not well understood. The Wnt-driven multiple intestinal neoplasia (MinApcΔ716/+) enterotoxigenic Bacteroides fragilis (ETBF) murine model is characterized by IL17-dependent, distal colon adenomas. Herein, we report that the addition of the BRAF V600E mutation to this model results in the emergence of a distinct locus of midcolon tumors. In ETBF-colonized BRAF V600E Lgr5 CreMin (BLM) mice, tumors have similarities to human BRAF V600E tumors, including histology, CpG island DNA hypermethylation, and immune signatures. In comparison to Min ETBF tumors, BLM ETBF tumors are infiltrated by CD8+ T cells, express IFNγ signatures, and are sensitive to anti-PD-L1 treatment. These results provide direct evidence for critical roles of host genetic and microbiota interactions in colorectal cancer pathogenesis and sensitivity to immunotherapy. SIGNIFICANCE: Colorectal cancers with BRAF mutations have distinct characteristics. We present evidence of specific colorectal cancer gene-microbial interactions in which colonization with toxigenic bacteria drives tumorigenesis in BRAF V600E Lgr5 CreMin mice, wherein tumors phenocopy aspects of human BRAF-mutated tumors and have a distinct IFNγ-dominant immune microenvironment uniquely responsive to immune checkpoint blockade.This article is highlighted in the In This Issue feature, p. 1601.

Constitutive activation of MEK5 promotes a mesenchymal and migratory cell phenotype in triple negative breast cancer

Triple negative breast cancer (TNBC) is an aggressive subtype of breast cancer with limited targeted therapeutic options. A defining feature of TNBC is the propensity to metastasize and acquire resistance to cytotoxic agents. Mitogen activated protein kinase (MAPK) and extracellular regulated kinase (ERK) signaling pathways have integral roles in cancer development and progression. While MEK5/ERK5 signaling drives mesenchymal and migratory cell phenotypes in breast cancer, the specific mechanisms underlying these actions remain under-characterized. To elucidate the mechanisms through which MEK5 regulates the mesenchymal and migratory phenotype, we generated stably transfected constitutively active MEK5 (MEK5-ca) TNBC cells. Downstream signaling pathways and candidate targets of MEK5-ca cells were based on RNA sequencing and confirmed using qPCR and Western blot analyses. MEK5 activation drove a mesenchymal cell phenotype independent of cell proliferation effects. Transwell migration assays demonstrated MEK5 activation significantly increased breast cancer cell migration. In this study, we provide supporting evidence that MEK5 functions through FRA-1 to regulate the mesenchymal and migratory phenotype in TNBC.

Dual inhibition of MEK1/2 and MEK5 suppresses the EMT/migration axis in triple-negative breast cancer through FRA-1 regulation

Triple-negative breast cancer (TNBC) presents a clinical challenge due to the aggressive nature of the disease and a lack of targeted therapies. Constitutive activation of the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) pathway has been linked to chemoresistance and metastatic progression through distinct mechanisms, including activation of epithelial-to-mesenchymal transition (EMT) when cells adopt a motile and invasive phenotype through loss of epithelial markers (CDH1), and acquisition of mesenchymal markers (VIM, CDH2). Although MAPK/ERK1/2 kinase inhibitors (MEKi) are useful antitumor agents in a clinical setting, including the Food and Drug Administration (FDA)-approved MEK1,2 dual inhibitors cobimetinib and trametinib, there are limitations to their clinical utility, primarily adaptation of the BRAF pathway and ocular toxicities. The MEK5 (HGNC: MAP2K5) pathway has important roles in metastatic progression of various cancer types, including those of the prostate, colon, bone and breast, and elevated levels of ERK5 expression in breast carcinomas are linked to a worse prognoses in TNBC patients. The purpose of this study is to explore MEK5 regulation of the EMT axis and to evaluate a novel pan-MEK inhibitor on clinically aggressive TNBC cells. Our results show a distinction between the MEK1/2 and MEK5 cascades in maintenance of the mesenchymal phenotype, suggesting that the MEK5 pathway may be necessary and sufficient in EMT regulation while MEK1/2 signaling further sustains the mesenchymal state of TNBC cells. Furthermore, additive effects on MET induction are evident through the inhibition of both MEK1/2 and MEK5. Taken together, these data demonstrate the need for a better understanding of the individual roles of MEK1/2 and MEK5 signaling in breast cancer and provide a rationale for the combined targeting of these pathways to circumvent compensatory signaling and subsequent therapeutic resistance.

Nitrogen cycling microbiomes are structured by plant mycorrhizal associations with consequences for nitrogen oxide fluxes in forests

Volatile nitrogen oxides (N2 O, NO, NO2 , HONO, …) can negatively impact climate, air quality, and human health. Using soils collected from temperate forests across the eastern United States, we show microbial communities involved in nitrogen (N) cycling are structured, in large part, by the composition of overstory trees, leading to predictable N-cycling syndromes, with consequences for emissions of volatile nitrogen oxides to air. Trees associating with arbuscular mycorrhizal (AM) fungi promote soil microbial communities with higher N-cycle potential and activity, relative to microbial communities in soils dominated by trees associating with ectomycorrhizal (ECM) fungi. Metagenomic analysis and gene expression studies reveal a 5 and 3.5 times greater estimated N-cycle gene and transcript copy numbers, respectively, in AM relative to ECM soil. Furthermore, we observe a 60% linear decrease in volatile reactive nitrogen gas flux (NOy  ≡ NO, NO2 , HONO) as ECM tree abundance increases. Compared to oxic conditions, gas flux potential of N2 O and NO increase significantly under anoxic conditions for AM soil (30- and 120-fold increase), but not ECM soil-likely owing to small concentrations of available substrate ( NO 3 - ) in ECM soil. Linear mixed effects modeling shows that ECM tree abundance, microbial process rates, and geographic location are primarily responsible for variation in peak potential NOy flux. Given that nearly all tree species associate with either AM or ECM fungi, our results indicate that the consequences of tree species shifts associated with global change may have predictable consequences for soil N cycling.

Pten and Dicer1 loss in the mouse uterus causes poorly differentiated endometrial adenocarcinoma

Endometrial cancer remains the most common gynecological malignancy in the United States. While the loss of the tumor suppressor, PTEN (phosphatase and tensin homolog), is well studied in endometrial cancer, recent studies suggest that DICER1, the endoribonuclease responsible for miRNA genesis, also plays a significant role in endometrial adenocarcinoma. Conditional uterine deletion of Dicer1 and Pten in mice resulted in poorly differentiated endometrial adenocarcinomas, which expressed Napsin A and HNF1B (hepatocyte nuclear factor 1 homeobox B), markers of clear-cell adenocarcinoma. Adenocarcinomas were hormone-independent. Treatment with progesterone did not mitigate poorly differentiated adenocarcinoma, nor did it affect adnexal metastasis. Transcriptomic analyses of DICER1 deleted uteri or Ishikawa cells revealed unique transcriptomic profiles and global miRNA downregulation. Computational integration of miRNA with mRNA targets revealed deregulated let-7 and miR-16 target genes, similar to published human DICER1-mutant endometrial cancers from TCGA (The Cancer Genome Atlas). Similar to human endometrial cancers, tumors exhibited dysregulation of ephrin-receptor signaling and transforming growth factor-beta signaling pathways. LIM kinase 2 (LIMK2), an essential molecule in p21 signal transduction, was significantly upregulated and represents a novel mechanism for hormone-independent pathogenesis of endometrial adenocarcinoma. This preclinical mouse model represents the first genetically engineered mouse model of poorly differentiated endometrial adenocarcinoma.

Epigenetic control of breast cancer susceptibility

1560

Background: Epigenetic mechanisms such as DNA methylation are important regulators of gene expression and are frequently dysregulated early in breast carcinogenesis. The relationship between DNA methylation aberrations in normal breast tissue and breast cancer risk remains unclear. Methods: Disease-free breast tissue cores donated by 71 high-risk (Tyrer-Cuzick lifetime risk ≥20%) and 79 average-risk women were obtained from the Komen Tissue Bank and processed for whole methylome (Diagenode's MethylCap Library and single-end 75-bp sequencing on Illumina Nextseq) and whole transcriptome (Illumina Nextseq) profiling. Reads from RNA-seq data were aligned to the human genome reference, GRCh38.p12 using STAR v.2.5.2b and tested for differential gene expression using DESeq2 ver. 1.24.0. For DNA methylation data, difference of variation in deduplicated read coverage among 250-bp fixed sized bins spanning CpG islands between high- and average-risk libraries was computed as z-ratios to identify differentially methylated regions. Pathway analysis was performed using IPA v06_01. Results: We identified 1355 CpGs that were differentially methylated between high- and average-risk breast tissues (ΔZ > 0.5, FDR < 0.05). Hypomethylated CpGs were overrepresented in high-risk tissue and were found predominantly (68%) in non-coding regions. Hypermethylated CpG sites were found equally in the gene body and non-coding regions. Transcriptomic analysis identified 112 differentially expressed genes (fold change≥2, FDR < 0.05), involved in chemokines signaling, metabolism and estrogen biosynthesis. Among those, FAM83A (logfc = 2.3, FDR = 0.004) was previously described as epigenetically dysregulated in multiple cancers and transforms breast epithelial cell in vitro. Methylation-expression correlations revealed 11 epigenetically regulated genes including cellular transformation-associated BMPR1B. Two hypomethylated/upregulated long non-coding RNAs were also identified in high-risk breasts. Conclusions: This is the first gene expression/DNA methylation analysis of normal breasts from women at either high or average risk of breast cancer. Our discovery of epigenetically regulated genes associated with breast cancer risk provides an opportunity to mechanistically dissect breast cancer susceptibility and risk-associated molecular alterations. Unlike the current focus of identifying germline mutations or single nucleotide polymorphisms responsible for higher risk, our studies reveal an epigenetic mechanism, which is not discernable through simple genomic sequencing.

Global H-NS counter-silencing by LuxR activates quorum sensing gene expression

Bacteria coordinate cellular behaviors using a cell-cell communication system termed quorum sensing. In Vibrio harveyi, the master quorum sensing transcription factor LuxR directly regulates >100 genes in response to changes in population density. Here, we show that LuxR derepresses quorum sensing loci by competing with H-NS, a global transcriptional repressor that oligomerizes on DNA to form filaments and bridges. We first identified H-NS as a repressor of bioluminescence gene expression, for which LuxR is a required activator. In an hns deletion strain, LuxR is no longer necessary for transcription activation of the bioluminescence genes, suggesting that the primary role of LuxR is to displace H-NS to derepress gene expression. Using RNA-seq and ChIP-seq, we determined that H-NS and LuxR co-regulate and co-occupy 28 promoters driving expression of 63 genes across the genome. ChIP-PCR assays show that as autoinducer concentration increases, LuxR protein accumulates at co-occupied promoters while H-NS protein disperses. LuxR is sufficient to evict H-NS from promoter DNA in vitro, which is dependent on LuxR DNA binding activity. From these findings, we propose a model in which LuxR serves as a counter-silencer at H-NS-repressed quorum sensing loci by disrupting H-NS nucleoprotein complexes that block transcription.

Probing biophysical sequence constraints within the transmembrane domains of rhodopsin by deep mutational scanning

Membrane proteins must balance the sequence constraints associated with folding and function against the hydrophobicity required for solvation within the bilayer. We recently found the expression and maturation of rhodopsin are limited by the hydrophobicity of its seventh transmembrane domain (TM7), which contains polar residues that are essential for function. On the basis of these observations, we hypothesized that rhodopsin's expression should be less tolerant of mutations in TM7 relative to those within hydrophobic TM domains. To test this hypothesis, we used deep mutational scanning to compare the effects of 808 missense mutations on the plasma membrane expression of rhodopsin in HEK293T cells. Our results confirm that a higher proportion of mutations within TM7 (37%) decrease rhodopsin's plasma membrane expression relative to those within a hydrophobic TM domain (TM2, 25%). These results in conjunction with an evolutionary analysis suggest solvation energetics likely restricts the evolutionary sequence space of polar TM domains.

The neurogenomic transition from territory establishment to parenting in a territorial female songbird

Background

The brain plays a critical role in upstream regulation of processes central to mating effort, parental effort, and self-maintenance. For seasonally breeding animals, the brain is likely mediating trade-offs among these processes within a short breeding season, yet research thus far has only explored neurogenomic changes from non-breeding to breeding states or select pathways (e.g., steroids) in male and/or lab-reared animals. Here, we use RNA-seq to explore neural plasticity in three behaviorally relevant neural tissues (ventromedial telencephalon [VmT], hypothalamus [HYPO], and hindbrain [HB]), comparing free-living female tree swallows (Tachycineta bicolor) as they shift from territory establishment to incubation. We additionally highlight changes in aggression-related genes to explore the potential for a neurogenomic shift in the mechanisms regulating aggression, a critical behavior both in establishing and maintaining a territory and in defense of offspring.

Results

HB had few differentially expressed genes, but VmT and HYPO had hundreds. In particular, VmT had higher expression of genes related to neuroplasticity and processes beneficial for competition during territory establishment, but down-regulated immune processes. HYPO showed signs of high neuroplasticity during incubation, and a decreased potential for glucocorticoid signaling. Expression of aggression-related genes also shifted from steroidal to non-steroidal pathways across the breeding season.

Conclusions

These patterns suggest trade-offs between enhanced activity and immunity in the VmT and between stress responsiveness and parental care in the HYPO, along with a potential shift in the mechanisms regulating aggression. Collectively, these data highlight important gene regulatory pathways that may underlie behavioral plasticity in females.

Tissue-specific expression profiles and positive selection analysis in the tree swallow (Tachycineta bicolor) using a de novo transcriptome assembly

Tree swallows (Tachycineta bicolor) are one of the most commonly studied wild birds in North America. They have advanced numerous research areas, including life history, physiology, and organismal responses to global change; however, transcriptomic resources are scarce. To further advance the utility of this system for biologists across disciplines, we generated a transcriptome for the tree swallow using six tissues (brain, blood, ovary, spleen, liver, and muscle) collected from breeding females. We de novo assembled 207,739 transcripts, which we aligned to 14,717 high confidence protein-coding genes. We then characterized each tissue with regard to its unique genes and processes and applied this transcriptome to two fundamental questions in evolutionary biology and endocrinology. First, we analyzed 3,015 single-copy orthologs and identified 46 genes under positive selection in the tree swallow lineage, including those with putative links to adaptations in this species. Second, we analyzed tissue-specific expression patterns of genes involved in sex steroidogenesis and processing. Enzymes capable of synthesizing these behaviorally relevant hormones were largely limited to the ovary, whereas steroid binding genes were found in nearly all other tissues, highlighting the potential for local regulation of sex steroid-mediated traits. These analyses provide new insights into potential sources of phenotypic variation in a free-living female bird and advance our understanding of fundamental questions in evolutionary and organismal biology.

Transcriptome Profiling Reveals Matrisome Alteration as a Key Feature of Ovarian Cancer Progression

BACKGROUND

Ovarian cancer is the most lethal gynecologic malignancy. There is a lack of comprehensive investigation of disease initiation and progression, including gene expression changes during early metastatic colonization.

METHODS

RNA-sequencing (RNA-seq) was done with matched primary tumors and fallopian tubes (n = 8 pairs) as well as matched metastatic and primary tumors (n = 11 pairs) from ovarian cancer patients. Since these are end point analyses, it was combined with RNA-seq using high-grade serous ovarian cancer cells seeded on an organotypic three-dimensional (3D) culture model of the omentum, mimicking early metastasis. This comprehensive approach revealed key changes in gene expression occurring in ovarian cancer initiation and metastasis, including early metastatic colonization.

RESULTS

2987 genes were significantly deregulated in primary tumors compared to fallopian tubes, 845 genes were differentially expressed in metastasis compared to primary tumors and 304 genes were common to both. An assessment of patient metastasis and 3D omental culture model of early metastatic colonization revealed 144 common genes that were altered during early colonization and remain deregulated even in the fully developed metastasis. Deregulation of the matrisome was a key process in early and late metastasis.

CONCLUSION

These findings will help in understanding the key pathways involved in ovarian cancer progression and eventually targeting those pathways for therapeutic interventions.

DNA methyltransferase inhibition reduces inflammation-induced colon tumorigenesis

Chronic inflammation is strongly associated with an increased risk of developing colorectal cancer. DNA hypermethylation of CpG islands alters the expression of genes in cancer cells and plays an important role in carcinogenesis. Chronic inflammation is also associated with DNA methylation alterations and in a mouse model of inflammation-induced colon tumorigenesis, we previously demonstrated that inflammation-induced tumours have 203 unique regions with DNA hypermethylation compared to uninflamed epithelium. To determine if altering inflammation-induced DNA hypermethylation reduces tumorigenesis, we used the same mouse model and treated mice with the DNA methyltransferase (DNMT) inhibitor decitabine (DAC) throughout the tumorigenesis time frame. DAC treatment caused a significant reduction in colon tumorigenesis. The tumours that did form after DAC treatment had reduced inflammation-specific DNA hypermethylation and alteration of expression of associated candidate genes. When compared, inflammation-induced tumours from control (PBS-treated) mice were enriched for cell proliferation associated gene expression pathways whereas inflammation-induced tumours from DAC-treated mice were enriched for interferon gene signatures. To further understand the altered tumorigenesis, we derived tumoroids from the different tumour types. Interestingly, tumoroids derived from inflammation-induced tumours from control mice maintained many of the inflammation-induced DNA hypermethylation alterations and had higher levels of DNA hypermethylation at these regions than tumoroids from DAC-treated mice. Importantly, tumoroids derived from inflammation-induced tumours from the DAC-treated mice proliferated more slowly than those derived from the inflammation-induced tumours from control mice. These studies suggest that inhibition of inflammation-induced DNA hypermethylation may be an effective strategy to reduce inflammation-induced tumorigenesis.

Co-occurring genomic capacity for anaerobic methane and dissimilatory sulfur metabolisms discovered in the Korarchaeota

Phylogenetic and geological evidence supports the hypothesis that life on Earth originated in thermal environments and conserved energy through methanogenesis or sulfur reduction. Here we describe two populations of the deeply rooted archaeal phylum Korarchaeota, which were retrieved from the metagenome of a circumneutral, suboxic hot spring that contains high levels of sulfate, sulfide, methane, hydrogen and carbon dioxide. One population is closely related to 'Candidatus Korarchaeum cryptofilum OPF8', while the more abundant korarchaeote, 'Candidatus Methanodesulfokores washburnensis', contains genes that are necessary for anaerobic methane and dissimilatory sulfur metabolisms. Phylogenetic and ancestral reconstruction analyses suggest that methane metabolism originated in the Korarchaeota, whereas genes for dissimilatory sulfite reduction were horizontally transferred to the Korarchaeota from the Firmicutes. Interactions among enzymes involved in both metabolisms could facilitate exergonic, sulfite-dependent, anaerobic oxidation of methane to methanol; alternatively, 'Ca. M. washburnensis' could conduct methanogenesis and sulfur reduction independently. Metabolic reconstruction suggests that 'Ca. M. washburnensis' is a mixotroph, capable of amino acid uptake, assimilation of methane-derived carbon and/or CO₂ fixation by archaeal type III-b RuBisCO for scavenging ribose carbon. Our findings link anaerobic methane metabolism and dissimilatory sulfur reduction within a single deeply rooted archaeal population and have implications for the evolution of these traits throughout the Archaea.

Marsarchaeota are an aerobic archaeal lineage abundant in geothermal iron oxide microbial mats

The discovery of archaeal lineages is critical to our understanding of the universal tree of life and evolutionary history of the Earth. Geochemically diverse thermal environments in Yellowstone National Park provide unprecedented opportunities for studying archaea in habitats that may represent analogues of early Earth. Here, we report the discovery and characterization of a phylum-level archaeal lineage proposed and herein referred to as the 'Marsarchaeota', after the red planet. The Marsarchaeota contains at least two major subgroups prevalent in acidic, microaerobic geothermal Fe(III) oxide microbial mats across a temperature range from ~50-80 °C. Metagenomics, single-cell sequencing, enrichment culturing and in situ transcriptional analyses reveal their biogeochemical role as facultative aerobic chemoorganotrophs that may also mediate the reduction of Fe(III). Phylogenomic analyses of replicate assemblies corresponding to two groups of Marsarchaeota indicate that they branch between the Crenarchaeota and all other major archaeal lineages. Transcriptomic analyses of several Fe(III) oxide mat communities reveal that these organisms were actively transcribing two different terminal oxidase complexes in situ and genes comprising an F₄₂₀-dependent butanal catabolism. The broad distribution of Marsarchaeota in geothermal, microaerobic Fe(III) oxide mats suggests that similar habitat types probably played an important role in the evolution of archaea.

In children, the microbiota of the nasopharynx and bronchoalveolar lavage fluid are both similar and different

RATIONALE

Sputum and bronchoalveolar lavage fluid (BALF) are often obtained to elucidate the lower airway microbiota in adults. Acquiring sputum samples from children is difficult and obtaining samples via bronchoscopy in children proves challenging due to the need for anesthesia and specialized procedural expertise; therefore nasopharyngeal (NP) swabs are often used as surrogates when investigating the pediatric airway microbiota. In adults, the airway microbiota differs significantly between NP and BALF samples however, minimal data exist in children.

OBJECTIVES

To compare NP and BALF samples in children undergoing clinically indicated bronchoscopy.

METHODS

NP and BALF samples were collected during clinically indicated bronchoscopy. Bacterial DNA was extracted from 72 samples (36 NP/BALF pairs); the bacterial V1-V3 region of the 16S rRNA gene was amplified and sequenced on the Illumina Miseq platform. Analysis was performed using mothur software.

RESULTS

Compared to NP samples, BALF had increased richness and diversity. Similarity between paired NP and BALF (intra-subject) samples was greater than inter-subject samples (P = 0.0006). NP samples contained more Actinobacteria (2.2% vs 21%; adjusted P = 1.4 × 10^-6 ), while BALF contained more Bacteroidetes (29.5% vs 3.2%; adjusted P = 1.2 × 10^-9 ). At the genus level several differences existed, however Streptococcus abundance was similar in both sample types (NP 37.3% vs BAL 36.1%; adjusted P = 0.8).

CONCLUSION

Our results provide evidence that NP samples can be used to distinguish differences between children, but the relative abundance of organisms may differ between the nasopharynx and lower airway in pediatric patients. Studies utilizing NP samples as surrogates for the lower airway should be interpreted with caution.

RNA polymerases IV and V influence the 3' boundaries of Polymerase II transcription units in Arabidopsis

Nuclear multisubunit RNA polymerases IV and V (Pol IV and Pol V) evolved in plants as specialized forms of Pol II. Their functions are best understood in the context of RNA-directed DNA methylation (RdDM), a process in which Pol IV-dependent 24 nt siRNAs direct the de novo cytosine methylation of regions transcribed by Pol V. Pol V has additional functions, independent of Pol IV and 24 nt siRNA biogenesis, in maintaining the repression of transposons and genomic repeats whose silencing depends on maintenance cytosine methylation. Here we report that Pol IV and Pol V play unexpected roles in defining the 3' boundaries of Pol II transcription units. Nuclear run-on assays reveal that in the absence of Pol IV or Pol V, Pol II occupancy downstream of poly A sites increases for approximately 12% of protein-coding genes. This effect is most pronounced for convergently transcribed gene pairs. Although Pols IV and V are detected near transcript ends of the affected Pol II - transcribed genes, their role in limiting Pol II read-through is independent of siRNA biogenesis or cytosine methylation for the majority of these genes. Interestingly, we observed that splicing was less efficient in pol IV or pol V mutant plants, compared to wild-type plants, suggesting that Pol IV or Pol V might affect pre-mRNA processing. We speculate that Pols IV and V (and/or their associated factors) play roles in Pol II transcription termination and pre-mRNA splicing by influencing polymerase elongation rates and/or release at collision sites for convergent genes.

Genes linked to species diversity in a sexually dimorphic communication signal in electric fish

Sexually dimorphic behaviors are often regulated by androgens and estrogens. Steroid receptors and metabolism are control points for evolutionary changes in sexual dimorphism. Electric communication signals of South American knifefishes are a model for understanding the evolution and physiology of sexually dimorphic behavior. These signals are regulated by gonadal steroids and controlled by a simple neural circuit. Sexual dimorphism of the signals varies across species. We used transcriptomics to examine mechanisms for sex differences in electric organ discharges (EODs) of two closely related species, Apteronotus leptorhynchus and Apteronotus albifrons, with reversed sexual dimorphism in their EODs. The pacemaker nucleus (Pn), which controls EOD frequency (EODf), expressed transcripts for steroid receptors and metabolizing enzymes, including androgen receptors, estrogen receptors, aromatase, and 5α-reductase. The Pn expressed mRNA for ion channels likely to regulate the high-frequency activity of Pn neurons and for neuromodulator and neurotransmitter receptors that may regulate EOD modulations used in aggression and courtship. Expression of several ion channel genes, including those for Kir3.1 inward-rectifying potassium channels and sodium channel β1 subunits, was sex-biased or correlated with EODf in ways consistent with EODf sex differences. Our findings provide a basis for future studies to characterize neurogenomic mechanisms by which sex differences evolve.

Early transcriptional response pathways in Daphnia magna are coordinated in networks of crustacean-specific genes

Natural habitats are exposed to an increasing number of environmental stressors that cause important ecological consequences. However, the multifarious nature of environmental change, the strength and the relative timing of each stressor largely limit our understanding of biological responses to environmental change. In particular, early response to unpredictable environmental change, critical to survival and fitness in later life stages, is largely uncharacterized. Here, we characterize the early transcriptional response of the keystone species Daphnia magna to twelve environmental perturbations, including biotic and abiotic stressors. We first perform a differential expression analysis aimed at identifying differential regulation of individual genes in response to stress. This preliminary analysis revealed that a few individual genes were responsive to environmental perturbations and they were modulated in a stressor and genotype-specific manner. Given the limited number of differentially regulated genes, we were unable to identify pathways involved in stress response. Hence, to gain a better understanding of the genetic and functional foundation of tolerance to multiple environmental stressors, we leveraged the correlative nature of networks and performed a weighted gene co-expression network analysis. We discovered that approximately one-third of the Daphnia genes, enriched for metabolism, cell signalling and general stress response, drives transcriptional early response to environmental stress and it is shared among genetic backgrounds. This initial response is followed by a genotype- and/or condition-specific transcriptional response with a strong genotype-by-environment interaction. Intriguingly, genotype- and condition-specific transcriptional response is found in genes not conserved beyond crustaceans, suggesting niche-specific adaptation.

Mismatch Repair Proteins Initiate Epigenetic Alterations during Inflammation-Driven Tumorigenesis

Aberrant silencing of genes by DNA methylation contributes to cancer, yet how this process is initiated remains unclear. Using a murine model of inflammation-induced tumorigenesis, we tested the hypothesis that inflammation promotes recruitment of epigenetic proteins to chromatin, initiating methylation and gene silencing in tumors. Compared with normal epithelium and noninflammation-induced tumors, inflammation-induced tumors gained DNA methylation at CpG islands, some of which are associated with putative tumor suppressor genes. Hypermethylated genes exhibited enrichment of repressive chromatin marks and reduced expression prior to tumorigenesis, at a time point coinciding with peak levels of inflammation-associated DNA damage. Loss of MutS homolog 2 (MSH2), a mismatch repair (MMR) protein, abrogated early inflammation-induced epigenetic alterations and DNA hypermethylation alterations observed in inflammation-induced tumors. These results indicate that early epigenetic alterations initiated by inflammation and MMR proteins lead to gene silencing during tumorigenesis, revealing a novel mechanism of epigenetic alterations in inflammation-driven cancer. Understanding such mechanisms will inform development of pharmacotherapies to reduce carcinogenesis. Cancer Res; 77(13); 3467-78. ©2017 AACR.

MicroRNA-335-5p and -3p synergize to inhibit estrogen receptor alpha expression and promote tamoxifen resistance

microRNAs (miRNAs) are small noncoding RNA molecules involved in the regulation of gene expression and play critical roles in human malignancies. Next-generation sequencing analysis of the MCF-7 breast cancer cell line overexpressing miR-335-5p and miR-335-3p demonstrated that the miRNA duplex repressed genes involved in the ERα signaling pathway, and enhanced resistance of MCF-7 cells to the growth inhibitory effects of tamoxifen. These data suggest that despite its conventional role in tumor suppression, the miR-335 transcript can also play an oncogenic role in promoting agonistic estrogen signaling in a cancerous setting.

Complete Transcriptome RNA-Seq

RNA-Seq is the leading technology for analyzing gene expression on a global scale across a broad spectrum of sample types. However, due to chemical modifications by fixation or degradation due to collection methods, samples often contain an abundance of RNA that is no longer intact, and the capability of current RNA-Seq protocols to accurately quantify such samples is often limited. We have developed an RNA-Seq protocol to address these key issues as well as quantify gene expression from the whole transcriptome. Furthermore, for compatibility with improved sequencing platforms, we use restructured adapter sequences to generate libraries for Illumina HiSeq, MiSeq, and NextSeq platforms. Our protocol utilizes duplex-specific nuclease (DSN) to remove abundant ribosomal RNA sequences while retaining other types of RNA for superior transcriptome profiling from low quantity input. We employ the Illumina sequencing platform, but this method is described in sufficient detail to adapt to other platforms.

The distribution, diversity and function of predominant Thermoproteales in high-temperature environments of Yellowstone National Park

High-temperature environments (> 70°C) contain diverse and abundant members of the crenarchaeal order Thermoproteales. However, a comprehensive study of the distribution and function of diverse members of this group across different habitat types has not been conducted. Consequently, the goals of this study were to determine the distribution of different Thermoproteales genera across geochemically distinct geothermal habitats of Yellowstone National Park, and to identify key functional attributes of major genera that correlate with environmental parameters. Curated sequence assemblies belonging to five genera were characterized in replicate samples of 11 high-temperature communities ranging in pH from 3 to 9. Thermocladium, Vulcanisaeta and Caldivirga spp. were the primary Thermoproteales populations present in low pH (pH < 5) habitats, whereas Thermoproteus populations were found in mildly-acidic (pH 5-6) sulfur sediments, and Pyrobaculum populations were confined to higher pH (pH > 6) sulfur sediments and/or filamentous 'streamer' communities. Metabolic reconstruction and comparative genomics among assemblies show that these populations are primarily chemoorganotrophs that utilize different electron acceptors depending on geochemical conditions. The presence of potential CO₂ fixation pathways in some Thermoproteales populations appears to be linked with NiFe hydrogenases, which combined with high levels of H₂ in many sulfidic systems, may provide the energy required to fix inorganic C.

Dual regulation by microRNA-200b-3p and microRNA-200b-5p in the inhibition of epithelial-to-mesenchymal transition in triple-negative breast cancer

Epithelial to mesenchymal transition (EMT) involves loss of an epithelial phenotype and activation of a mesenchymal one. Enhanced expression of genes associated with a mesenchymal transition includes ZEB1/2, TWIST, and FOXC1. miRNAs are known regulators of gene expression and altered miRNA expression is known to enhance EMT in breast cancer. Here we demonstrate that the tumor suppressive miRNA family, miR-200, is not expressed in triple negative breast cancer (TNBC) cell lines and that miR-200b-3p over-expression represses EMT, which is evident through decreased migration and increased CDH1 expression. Despite the loss of migratory capacity following re-expression of miR-200b-3p, no subsequent loss of the conventional miR-200 family targets and EMT markers ZEB1/2 was observed. Next generation RNA-sequencing analysis showed that enhanced expression of pri-miR-200b lead to ectopic expression of both miR-200b-3p and miR-200b-5p with multiple isomiRs expressed for each of these miRNAs. Furthermore, miR-200b-5p was expressed in the receptor positive, epithelial breast cancer cell lines but not in the TNBC (mesenchymal) cell lines. In addition, a compensatory mechanism for miR-200b-3p/200b-5p targeting, where both miRNAs target the RHOGDI pathway leading to non-canonical repression of EMT, was demonstrated. Collectively, these data are the first to demonstrate dual targeting by miR-200b-3p and miR-200b-5p and a previously undescribed role for microRNA processing and strand expression in EMT and TNBC, the most aggressive breast cancer subtype.

Novel, Deep-Branching Heterotrophic Bacterial Populations Recovered from Thermal Spring Metagenomes

Thermal spring ecosystems are a valuable resource for the discovery of novel hyperthermophilic Bacteria and Archaea, and harbor deeply-branching lineages that provide insight regarding the nature of early microbial life. We characterized bacterial populations in two circumneutral (pH ~8) Yellowstone National Park thermal (T ~80°C) spring filamentous “streamer” communities using random metagenomic DNA sequence to investigate the metabolic potential of these novel populations. Four de novo assemblies representing three abundant, deeply-branching bacterial phylotypes were recovered. Analysis of conserved phylogenetic marker genes indicated that two of the phylotypes represent separate groups of an uncharacterized phylum (for which we propose the candidate phylum name “Pyropristinus”). The third new phylotype falls within the proposed Calescamantes phylum. Metabolic reconstructions of the “Pyropristinus” and Calescamantes populations showed that these organisms appear to be chemoorganoheterotrophs and have the genomic potential for aerobic respiration and oxidative phosphorylation via archaeal-like V-type, and bacterial F-type ATPases, respectively. A survey of similar phylotypes (>97% nt identity) within 16S rRNA gene datasets suggest that the newly described organisms are restricted to terrestrial thermal springs ranging from 70 to 90°C and pH values of ~7–9. The characterization of these lineages is important for understanding the diversity of deeply-branching bacterial phyla, and their functional role in high-temperature circumneutral “streamer” communities.

Recombination Does Not Hinder Formation or Detection of Ecological Species of Synechococcus Inhabiting a Hot Spring Cyanobacterial Mat

Recent studies of bacterial speciation have claimed to support the biological species concept—that reduced recombination is required for bacterial populations to diverge into species. This conclusion has been reached from the discovery that ecologically distinct clades show lower rates of recombination than that which occurs among closest relatives. However, these previous studies did not attempt to determine whether the more-rapidly recombining close relatives within the clades studied may also have diversified ecologically, without benefit of sexual isolation. Here we have measured the impact of recombination on ecological diversification within and between two ecologically distinct clades (A and B') of Synechococcus in a hot spring microbial mat in Yellowstone National Park, using a cultivation-free, multi-locus approach. Bacterial artificial chromosome (BAC) libraries were constructed from mat samples collected at 60°C and 65°C. Analysis of multiple linked loci near Synechococcus 16S rRNA genes showed little evidence of recombination between the A and B' lineages, but a record of recombination was apparent within each lineage. Recombination and mutation rates within each lineage were of similar magnitude, but recombination had a somewhat greater impact on sequence diversity than mutation, as also seen in many other bacteria and archaea. Despite recombination within the A and B' lineages, there was evidence of ecological diversification within each lineage. The algorithm Ecotype Simulation identified sequence clusters consistent with ecologically distinct populations (ecotypes), and several hypothesized ecotypes were distinct in their habitat associations and in their adaptations to different microenvironments. We conclude that sexual isolation is more likely to follow ecological divergence than to precede it. Thus, an ecology-based model of speciation appears more appropriate than the biological species concept for bacterial and archaeal diversification.

Identification of novel immunomodulatory tumor biology through comprehensive characterization of a metastases-specific epigenome in patients with metachronous primary and metastatic urothelial carcinoma (UC) tumor pairs

452

Background: The Cancer Genome Atlas project identified a distinct cluster of hypermethylated muscle-invasive UC tumors in one third of patients. It is unknown if this epigenetic biology persists or changes in UC metastases. The present study aimed to identify uniquely hypermethylated regions in UC metastases compared to their matched primary tumors. Methods: UC patients with paired formalin fixed paraffin embedded tissue samples ( > 75% tumor) from primary and metachronous metastatic tumors were identified. DNA methylation was analyzed on the Illumina HiSeq platform by double-enzyme reduced representation bisulfite sequencing (dRRBS). Results were analyzed by BS-seq within BSmooth open source software. Alignment and methylation estimation was done using Bismark open source software. CpG regions with coverage < 2 in more than 66% of samples were removed to reduce false positive results. Differentially methylated regions (DMRs) in metastases compared to primary tumors were determined by paired t-test. Gene set enrichment analysis (GSEA) of the top 5% DMRs was performed utilizing multiple gene set collections including c7IMMUNO. A false discovery rate of < 10% defined significant DMRs. Results: 15 UC primary/metachronous metastases pairs were analyzed. After filtering for low coverage, 1,781,762 loci remained for analysis. After merging loci within 2500 bp of each other, the top 5% DMRs resulted in 18,452 DMRs. GSEA including the c7IMMUNO gene set identified multiple previously undescribed hypermethylated genes modulating immune function through TGF-beta and Treg signaling including: TGFBR2, TGFBR3, SMAD1, SMAD3, SMAD4, BACH1, BACH2, and VDR all q < 0.05. Conclusions: We identified numerous genes with immunomodulatory functions significantly hypermethylated in UC metastases compared to their matched primary tumors. Our findings provide rationale to examine epigenetic approaches as a means to improve clinical outcomes of UC patients treated with immunotherapy. Our small sample size limits definitive conclusions and warrants validation in independent data sets.

Ecophysiology of an uncultivated lineage of Aigarchaeota from an oxic, hot spring filamentous 'streamer' community

The candidate archaeal phylum 'Aigarchaeota' contains microorganisms from terrestrial and subsurface geothermal ecosystems. The phylogeny and metabolic potential of Aigarchaeota has been deduced from several recent single-cell amplified genomes; however, a detailed description of their metabolic potential and in situ transcriptional activity is absent. Here, we report a comprehensive metatranscriptome-based reconstruction of the in situ metabolism of Aigarchaeota in an oxic, hot spring filamentous 'streamer' community. Fluorescence in situ hybridization showed that these newly discovered Aigarchaeota are filamentous, which is consistent with the presence and transcription of an actin-encoding gene. Aigarchaeota filaments are intricately associated with other community members, which include both bacteria (for example, filamentous Thermocrinis spp.) and archaea. Metabolic reconstruction of genomic and metatranscriptomic data suggests that this aigarchaeon is an aerobic, chemoorganoheterotroph with autotrophic potential. A heme copper oxidase complex was identified in the environmental genome assembly and highly transcribed in situ. Potential electron donors include acetate, fatty acids, amino acids, sugars and aromatic compounds, which may originate from extracellular polymeric substances produced by other microorganisms shown to exist in close proximity and/or autochthonous dissolved organic carbon (OC). Transcripts related to genes specific to each of these potential electron donors were identified, indicating that this aigarchaeon likely utilizes several OC substrates. Characterized members of this lineage cannot synthesize heme, and other cofactors and vitamins de novo, which suggests auxotrophy. We propose the name Candidatus 'Calditenuis aerorheumensis' for this aigarchaeon, which describes its filamentous morphology and its primary electron acceptor, oxygen.