Search

Abstract Background Though genomic-level data are becoming widely available, many of the metazoan species sequenced are laboratory systems whose natural history is not well documented. In contrast, the wide array of species with very well-characterized natural history have, until recently, lacked genomics tools. It is now possible to address significant evolutionary genomics questions by applying high-throughput sequencing to discover the majority of genes for ecologically tractable species, and by subsequently developing microarray platforms from which to investigate gene regulatory networks that function in natural systems. We used GS-FLX Titanium Sequencing (Roche/454-Sequencing) of two normalized libraries of pooled RNA samples to characterize a transcriptome of the dark-eyed junco (Junco hyemalis), a North American sparrow that is a classically studied species in the fields of photoperiodism, speciation, and hormone-mediated behavior. Results From a broad pool of RNA sampled from tissues throughout the body of a male and a female junco, we sequenced a total of 434 million nucleotides from 1.17 million reads that were assembled de novo into 31,379 putative transcripts representing 22,765 gene sets covering 35.8 million nucleotides with 12-fold average depth of coverage. Annotation of roughly half of the putative genes was accomplished using sequence similarity, and expression was confirmed for the majority with a preliminary microarray analysis. Of 716 core bilaterian genes, 646 (90 %) were recovered within our characterized gene set. Gene Ontology, orthoDB orthology groups, and KEGG Pathway annotation provide further functional information about the sequences, and 25,781 potential SNPs were identified. Conclusions The extensive sequence information returned by this effort adds to the growing store of genomic data on diverse species. The extent of coverage and annotation achieved and confirmation of expression, show that transcriptome sequencing provides useful information for ecological model systems that have historically lacked genomic tools. The junco-specific microarray developed here is allowing investigations of gene expression responses to environmental and hormonal manipulations – extending the historic work on natural history and hormone-mediated phenotypes in this system.

Vibrio coralliilyticus is a pathogen of coral and shellfish, leading to devastating economic and ecological consequences worldwide. Although rising ocean temperatures correlate with increased V. coralliilyticus pathogenicity, the specific molecular mechanisms and determinants contributing to virulence remain poorly understood. Here, we systematically analyzed the type VI secretion system (T6SS), a contact-dependent toxin delivery apparatus, in V. coralliilyticus. We identified 2 omnipresent T6SSs that are activated at temperatures in which V. coralliilyticus becomes virulent; T6SS1 is an antibacterial system mediating interbacterial competition, whereas T6SS2 mediates anti-eukaryotic toxicity and contributes to mortality during infection of an aquatic model organism, Artemia salina. Using comparative proteomics, we identified the T6SS1 and T6SS2 toxin arsenals of 3 V. coralliilyticus strains with distinct disease etiologies. Remarkably, T6SS2 secretes at least 9 novel anti-eukaryotic toxins comprising core and accessory repertoires. We propose that T6SSs differently contribute to V. coralliilyticus's virulence: T6SS2 plays a direct role by targeting the host, while T6SS1 plays an indirect role by eliminating competitors., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Mass et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The type III secretion system (T3SS) is a needle-like complex used by numerous bacterial pathogens in host infection by directly injecting exotoxins into the host cell cytoplasm, leading to cell death. The T3SS is a known virulence factor in the shrimp pathogen Vibrio campbellii . The ∼40 genes comprising the V. campbellii T3SS are regulated by a network of transcription factors in response to changes in the cell's environment: cell density (quorum sensing; QS), temperature, calcium, and host cell contact. Under positive environmental stimuli, the master T3SS transcription factor ExsA activates expression of the four structural T3SS operons required for needle formation. Previous studies identified a key role of the master QS transcription factor LuxR: repression of exsA transcription via DNA binding at the exsBA promoter. Here we uncovered a new regulatory role of LuxR: indirect post-translational repression of ExsA activity via direct transcriptional repression of the gene encoding the anti-anti-activator ExsC. In V. campbellii , ExsC is a positive regulator of T3SS transcription: high ExsC expression leads to full ExsA transcription activation of the T3SS structural promoters. LuxR binding at the exsC promoter represses transcription of exsC through disruption of ExsA binding. Our findings collectively show that V. campbellii responds to high cell density signals to shut down the expression of the T3SS. We postulate that this dual regulatory mechanism by LuxR enables both the rapid inactivation of existing ExsA protein and blocks its further synthesis, leading to a rapid shutdown of T3SS activity at high cell density., Importance: Vibrio campbellii utilizes the type III secretion system (T3SS) as a mechanism of pathogenesis, which is a highly studied 'injectisome' complex that delivers exotoxins into host cells during infection. The T3SS pathogenicity island in V. campbellii comprises ∼40 genes that are organized into four structural operons. In this study, we determined that quorum sensing - a method of bacterial communication - regulates T3SS genes both at the transcriptional and post-translational levels to shut down T3SS gene expression at high population densities.

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., (© 2023. The Author(s).)

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., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Genetics Society of America.)

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., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

Ovarian cancer is a deadly disease attributed to late-stage detection as well as recurrence and the development of chemoresistance. Ovarian cancer stem cells (OCSCs) are hypothesized to be largely responsible for the emergence of chemoresistant tumors. Although chemotherapy may initially succeed at decreasing the size and number of tumors, it leaves behind residual malignant OCSCs. In this study, we demonstrate that aldehyde dehydrogenase 1A1 (ALDH1A1) is essential for the survival of OCSCs. We identified a first-in-class ALDH1A1 inhibitor, compound 974 , and used 974 as a tool to decipher the mechanism of stemness regulation by ALDH1A1. The treatment of OCSCs with 974 significantly inhibited ALDH activity, the expression of stemness genes, and spheroid and colony formation. An in vivo limiting dilution assay demonstrated that 974 significantly inhibited CSC frequency. A transcriptomic sequencing of cells treated with 974 revealed a significant downregulation of genes related to stemness and chemoresistance as well as senescence and the senescence-associated secretory phenotype (SASP). We confirmed that 974 inhibited the senescence and stemness induced by platinum-based chemotherapy in functional assays. Overall, these data establish that ALDH1A1 is essential for OCSC survival and that ALDH1A1 inhibition suppresses chemotherapy-induced senescence and stemness. Targeting ALDH1A1 using small-molecule inhibitors in combination with chemotherapy therefore presents a promising strategy to prevent ovarian cancer recurrence and has the potential for clinical translation.

Microbiota studies have reported changes in the microbial composition of the breast upon cancer development. However, results are inconsistent and limited to the later phases of cancer development (after diagnosis). We analyzed and compared the resident bacterial taxa of histologically normal breast tissue (healthy, H, n  = 49) with those of tissues donated prior to (prediagnostic, PD, n  = 15) and after (adjacent normal, AN, n  = 49, and tumor, T, n  = 46) breast cancer diagnosis ( n total = 159). DNA was isolated from tissue samples and submitted for Illumina MiSeq paired-end sequencing of the V3-V4 region of the 16S gene. To infer bacterial function in breast cancer, we predicted the functional bacteriome from the 16S sequencing data using PICRUSt2. Bacterial compositional analysis revealed an intermediary taxonomic signature in the PD tissue relative to that of the H tissue, represented by shifts in Bacillaceae , Burkholderiaceae , Corynebacteriaceae , Streptococcaceae , and Staphylococcaceae . This compositional signature was enhanced in the AN and T tissues. We also identified significant metabolic reprogramming of the microbiota of the PD, AN, and T tissue compared with the H tissue. Further, preliminary correlation analysis between host transcriptome profiling and microbial taxa and genes in H and PD tissues identified altered associations between the human host and mammary microbiota in PD tissue compared with H tissue. These findings suggest that compositional shifts in bacterial abundance and metabolic reprogramming of the breast tissue microbiota are early events in breast cancer development that are potentially linked with cancer susceptibility. IMPORTANCE The goal of this study was to determine the role of resident breast tissue bacteria in breast cancer development. We analyzed breast tissue bacteria in healthy breast tissue and breast tissue donated prior to (precancerous) and after (postcancerous) breast cancer diagnosis. Compared to healthy tissue, the precancerous and postcancerous breast tissues demonstrated differences in the amounts of breast tissue bacteria. In addition, breast tissue bacteria exhibit different functions in pre-cancerous and post-cancerous breast tissues relative to healthy tissue. These differences in function are further emphasized by altered associations of the breast tissue bacteria with gene expression in the human host prior to cancer development. Collectively, these analyses identified shifts in bacterial abundance and metabolic function (dysbiosis) prior to breast tumor diagnosis. This dysbiosis may serve as a therapeutic target in breast cancer prevention.

Extensive portions of the human genome have unknown function, including those derived from transposable elements. One such element, the DNA transposon Hsmar1, entered the primate lineage approximately 50 million years ago leaving behind terminal inverted repeat (TIR) sequences and a single intact copy of the Hsmar1 transposase, which retains its ancestral TIR-DNA-binding activity, and is fused with a lysine methyltransferase SET domain to constitute the chimeric SETMAR gene. Here, we provide a structural basis for recognition of TIRs by SETMAR and investigate the function of SETMAR through genome-wide approaches. As elucidated in our 2.37 Å crystal structure, SETMAR forms a dimeric complex with each DNA-binding domain bound specifically to TIR-DNA through the formation of 32 hydrogen bonds. We found that SETMAR recognizes primarily TIR sequences (∼5000 sites) within the human genome as assessed by chromatin immunoprecipitation sequencing analysis. In two SETMAR KO cell lines, we identified 163 shared differentially expressed genes and 233 shared alternative splicing events. Among these genes are several pre-mRNA-splicing factors, transcription factors, and genes associated with neuronal function, and one alternatively spliced primate-specific gene, TMEM14B, which has been identified as a marker for neocortex expansion associated with brain evolution. Taken together, our results suggest a model in which SETMAR impacts differential expression and alternative splicing of genes associated with transcription and neuronal function, potentially through both its TIR-specific DNA-binding and lysine methyltransferase activities, consistent with a role for SETMAR in simian primate development., Competing Interests: Conflict of interest R. C. W. has received grant support from Eli Lilly and Company, and R. C. W. serves as a scientific advisor to HiberCell; all other authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Family with sequence similarity 83 member A (FAM83A) presents oncogenic properties in several cancers including breast cancer. Recently, we reported FAM83A overexpression in normal breast tissues from women at high risk of breast cancer. We now hypothesize that FAM83A is a key factor in breast cancer initiation., Methods: Immunohistochemical staining was used to evaluate FAM83A protein levels in both a normal breast tissue microarray (TMA, N = 411) and a breast tumor TMA (N = 349). EGFR staining and its correlation with FAM83A expression were also assessed. Lentivirus-mediated manipulation of FAM83A expression in primary and hTERT-immortalized breast epithelial cells was employed. Biological and molecular alterations upon FAM83A overexpression/downregulation and FAM83A's interaction partners were investigated., Results: TMA analysis revealed a 1.5-fold increase in FAM83A expression level in breast cancer cases as compared with normal breast tissues (p < 0.0001). FAM83A protein expression was directly correlated with EGFR level in both normal and breast cancer tissues. In in vitro assays, exogenous expression of FAM83A in either primary or immortalized breast epithelial cells promoted cell viability and proliferation. Additionally, Ingenuity Pathway Analysis (IPA) revealed that FAM83A overexpression in primary cells affected the expression of genes involved in cellular morphology and metabolism. Mass spectrometry analysis identified DDX3X and LAMB3 as potential FAM83A interaction partners in primary cells, while we detected FAM83A interaction with cytoskeleton reorganization factors, including LIMA1, MYH10, PLEC, MYL6 in the immortalized cells., Conclusions: This study shows that FAM83A promotes metabolic activation in primary breast epithelial cells and cell proliferation in both primary and immortalized cells. These findings support its role in early breast oncogenesis., (© 2022. The Author(s).)

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., (© 2022. The Author(s).)

Swimming motility is a critical virulence factor in pathogenesis for numerous Vibrio species. Vibrio campbellii DS40M4 is a wild-type isolate that has been recently established as a highly tractable model strain for bacterial genetics studies. We sought to exploit the tractability and relevance of this strain for characterization of flagellar gene regulation in V. campbellii . Using comparative genomics, we identified homologs of V. campbellii flagellar and chemotaxis genes conserved in other members of the Vibrionaceae and determined the transcriptional profile of these loci using differential RNA-seq. We systematically deleted all 63 predicted flagellar and chemotaxis genes in V. campbellii and examined their effects on motility and flagellum production. We specifically focused on the core regulators of the flagellar hierarchy established in other vibrios: RpoN (σ 54 ), FlrA, FlrC, and FliA. Our results show that V. campbellii transcription of flagellar and chemotaxis genes is governed by a multitiered regulatory hierarchy similar to other motile Vibrio species. However, there are several critical differences in V. campbellii : (i) the σ 54 -dependent regulator FlrA is dispensable for motility; (ii) the flgA , fliEFGHIJ , flrA , and flrBC operons do not require σ 54 for expression; and (iii) FlrA and FlrC coregulate class II genes. Our model proposes that the V. campbellii flagellar transcriptional hierarchy has three classes of genes, in contrast to the four-class hierarchy in Vibrio cholerae. Our genetic and phenotypic dissection of the V. campbellii flagellar regulatory network highlights the differences that have evolved in flagellar regulation across the Vibrionaceae . IMPORTANCE Vibrio campbellii is a Gram-negative bacterium that is free-living and ubiquitous in marine environments and is an important global pathogen of fish and shellfish. Disruption of the flagellar motor significantly decreases host mortality of V. campbellii , suggesting that motility is a key factor in pathogenesis. Using this model organism, we identified >60 genes that encode proteins with predicted structural, mechanical, or regulatory roles in function of the single polar flagellum in V. campbellii . We systematically tested strains containing single deletions of each gene to determine the impact on motility and flagellum production. Our studies have uncovered differences in the regulatory network and function of several genes in V. campbellii compared to established systems in Vibrio cholerae and Vibrio parahaemolyticus.

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.

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., (© 2021 Society for Applied Microbiology and John Wiley & Sons Ltd.)

High-grade serous ovarian cancer (HGSOC) is characterized by chromosomal instability, DNA damage, oxidative stress, and high metabolic demand that exacerbate misfolded, unfolded, and damaged protein burden resulting in increased proteotoxicity. However, the underlying mechanisms that maintain protein homeostasis to promote HGSOC growth remain poorly understood. This study reports that the neuronal deubiquitinating enzyme, ubiquitin carboxyl-terminal hydrolase L1 (UCHL1), is overexpressed in HGSOC and maintains protein homeostasis. UCHL1 expression was markedly increased in HGSOC patient tumors and serous tubal intraepithelial carcinoma (HGSOC precursor lesions). High UCHL1 levels correlated with higher tumor grade and poor patient survival. UCHL1 inhibition reduced HGSOC cell proliferation and invasion, as well as significantly decreased the in vivo metastatic growth of ovarian cancer xenografts. Transcriptional profiling of UCHL1-silenced HGSOC cells revealed downregulation of genes implicated with proteasome activity along with upregulation of endoplasmic reticulum stress-induced genes. Reduced expression of proteasome subunit alpha 7 (PSMA7) and acylaminoacyl peptide hydrolase (APEH), upon silencing of UCHL1, resulted in a significant decrease in proteasome activity, impaired protein degradation, and abrogated HGSOC growth. Furthermore, the accumulation of polyubiquitinated proteins in the UCHL1-silenced cells led to attenuation of mTORC1 activity and protein synthesis, and induction of terminal unfolded protein response. Collectively, these results indicate that UCHL1 promotes HGSOC growth by mediating protein homeostasis through the PSMA7-APEH-proteasome axis. IMPLICATIONS: This study identifies the novel links in the proteostasis network to target protein homeostasis in HGSOC and recognizes the potential of inhibiting UCHL1 and APEH to sensitize cancer cells to proteotoxic stress in solid tumors., (©2021 American Association for Cancer Research.)

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 (Min ApcΔ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 Cre Min (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 Cre Min 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 ., (©2021 American Association for Cancer Research.)

Periods of social instability can elicit adaptive phenotypic plasticity to promote success in future competition. However, the underlying molecular mechanisms have primarily been studied in captive and laboratory-reared animals, leaving uncertainty as to how natural competition among free-living animals affects gene activity. Here, we experimentally generated social competition among wild, cavity-nesting female birds (tree swallows, Tachycineta bicolor ). After territorial settlement, we reduced the availability of key breeding resources (i.e., nest boxes), generating heightened competition; within 24 h we reversed the manipulation, causing aggressive interactions to subside. We sampled females during the peak of competition and 48 h after it ended, along with date-matched controls. We measured transcriptomic and epigenomic responses to competition in two socially relevant brain regions (hypothalamus and ventromedial telencephalon). Gene network analyses suggest that processes related to energy mobilization and aggression (e.g., dopamine synthesis) were up-regulated during competition, the latter of which persisted 2 d after competition had ended. Cellular maintenance processes were also down-regulated after competition. Competition additionally altered methylation patterns, particularly in pathways related to hormonal signaling, suggesting those genes were transcriptionally poised to respond to future competition. Thus, experimental competition among free-living animals shifts gene expression in ways that may facilitate the demands of competition at the expense of self-maintenance. Further, some of these effects persisted after competition ended, demonstrating the potential for epigenetic biological embedding of the social environment in ways that may prime individuals for success in future social instability., Competing Interests: The authors declare no competing interest.

Sequence-specific DNA-binding domains (DBDs) are conserved in all domains of life. These proteins carry out a variety of cellular functions, and there are a number of distinct structural domains already described that allow for sequence-specific DNA binding, including the ubiquitous helix-turn-helix (HTH) domain. In the facultative pathogen Vibrio cholerae , the chitin sensor ChiS is a transcriptional regulator that is critical for the survival of this organism in its marine reservoir. We recently showed that ChiS contains a cryptic DBD in its C terminus. This domain is not homologous to any known DBD, but it is a conserved domain present in other bacterial proteins. Here, we present the crystal structure of the ChiS DBD at a resolution of 1.28 Å. We find that the ChiS DBD contains an HTH domain that is structurally similar to those found in other DNA-binding proteins, like the LacI repressor. However, one striking difference observed in the ChiS DBD is that the canonical tight turn of the HTH is replaced with an insertion containing a β-sheet, a variant which we term the helix-sheet-helix. Through systematic mutagenesis of all positively charged residues within the ChiS DBD, we show that residues within and proximal to the ChiS helix-sheet-helix are critical for DNA binding. Finally, through phylogenetic analyses we show that the ChiS DBD is found in diverse proteobacterial proteins that exhibit distinct domain architectures. Together, these results suggest that the structure described here represents the prototypical member of the ChiS-family of DBDs. IMPORTANCE Regulating gene expression is essential in all domains of life. This process is commonly facilitated by the activity of DNA-binding transcription factors. There are diverse structural domains that allow proteins to bind to specific DNA sequences. The structural basis underlying how some proteins bind to DNA, however, remains unclear. Previously, we showed that in the major human pathogen Vibrio cholerae , the transcription factor ChiS directly regulates gene expression through a cryptic DNA-binding domain. This domain lacked homology to any known DNA-binding protein. In the current study, we determined the structure of the ChiS DNA-binding domain (DBD) and found that the ChiS-family DBD is a cryptic variant of the ubiquitous helix-turn-helix (HTH) domain. We further demonstrate that this domain is conserved in diverse proteins that may represent a novel group of transcriptional regulators., (Copyright © 2021 Klancher et al.)

In Vibrio species, chitin-induced natural transformation enables bacteria to take up DNA from the external environment and integrate it into their genome. Expression of the master competence regulator TfoX bypasses the need for chitin induction and drives expression of the genes required for competence in several Vibrio species. Here, we show that TfoX expression in Vibrio campbellii strains DS40M4 and NBRC 15631 enables high natural transformation frequencies. Conversely, transformation was not achieved in the model quorum-sensing strain V. campbellii BB120 (previously classified as Vibrio harveyi ). Surprisingly, we find that quorum sensing is not required for transformation in V. campbellii DS40M4 or Vibrio parahaemolyticus in contrast to the established regulatory pathway in Vibrio cholerae in which quorum sensing is required to activate the competence regulator QstR. Similar to V. cholerae , expression of both QstR and TfoX is necessary for transformation in DS40M4. There is a wide disparity in transformation frequencies among even closely related Vibrio strains, with V. vulnificus having the lowest functional transformation frequency. Ectopic expression of both TfoX and QstR is sufficient to produce a significant increase in transformation frequency in Vibrio vulnificus To explore differences in competence regulation, we used previously studied V. cholerae competence genes to inform a comparative genomics analysis coupled with transcriptomics. We find that transformation capability cannot necessarily be predicted by the level of gene conservation but rather correlates with competence gene expression following TfoX induction. Thus, we have uncovered notable species- and strain-level variations in the competence gene regulation pathway across the Vibrio genus. IMPORTANCE Naturally transformable, or competent, bacteria are able to take up DNA from their environment, a key method of horizontal gene transfer for acquisition of new DNA sequences. Our research shows that Vibrio species that inhabit marine environments exhibit a wide diversity in natural transformation capability ranging from nontransformability to high transformation rates in which 10% of cells measurably incorporate new DNA. We show that the role of regulatory systems controlling the expression of competence genes (e.g., quorum sensing) differs throughout both the species and strain levels. We explore natural transformation capabilities of Vibrio campbellii species which have been thus far uncharacterized and find novel regulation of competence. Expression of two key transcription factors, TfoX and QstR, is necessary to stimulate high levels of transformation in Vibrio campbellii and recover low rates of transformation in Vibrio vulnificus ., (Copyright © 2019 Simpson et al.)

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.

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.

Small RNAs (sRNAs) that are 21 to 24 nucleotides (nt) in length are found in most eukaryotic organisms and regulate numerous biological functions, including transposon silencing, development, reproduction, and stress responses, typically via control of the stability and/or translation of target mRNAs. Major classes of sRNAs in plants include microRNAs (miRNAs) and small interfering RNAs (siRNAs); sRNAs are known to travel as a silencing signal from cell to cell, root to shoot, and even between host and pathogen. In mammals, sRNAs are transported inside extracellular vesicles (EVs), which are mobile membrane-bound compartments that participate in intercellular communication. In addition to sRNAs, EVs carry proteins, lipids, metabolites, and potentially other types of nucleic acids. Here we report that Arabidopsis ( Arabidopsis thaliana ) EVs also contain diverse species of sRNA. We found that specific miRNAs and siRNAs are preferentially loaded into plant EVs. We also report a previously overlooked class of "tiny RNAs" (10 to 17 nt) that are highly enriched in EVs. This RNA category of unknown function has a broad and very diverse genome origin and might correspond to degradation products., (© 2019 American Society of Plant Biologists. All rights reserved.)

Noncoding RNAs perform diverse regulatory functions in living cells. In plants, two RNA polymerase II-related enzymes, RNA polymerases IV and V (Pol IV and V), specialize in the synthesis of noncoding RNAs that silence a subset of transposable elements and genes via RNA-directed DNA methylation (RdDM). In this process, Pol IV partners with RNA-dependent RNA polymerase 2 (RDR2) to produce double-stranded RNAs that are then cut by an RNase III enzyme, Dicer-like 3 (DCL3), into 24 nt small interfering RNAs (siRNAs). The siRNAs are loaded into an Argonaute family protein, primarily AGO4, and guide the complex to complementary DNA target sequences where RdDM and repressive chromatin modifications ensue. The dependence of 24 nt siRNA biogenesis on Pol IV and RDR2 has been known for more than a decade, but the elusive pre-siRNA transcripts synthesized by Pol IV and RDR2 have only recently been identified. This chapter describes the approaches that enabled our identification of Pol IV/RDR2-dependent RNAs (P4R2 RNAs) in Arabidopsis thaliana. These included the use of a triple Dicer mutant (dcl2 dcl3 dcl4) to cause P4R2 RNAs to accumulate, genome-wide identification and mapping of P4R2 RNAs using a modified Illumina small RNA-Seq protocol, and multiple bioinformatic pipelines for data analysis and displaying results.

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.

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., (© 2017 John Wiley & Sons Ltd.)

Defense responses of peanut ( Arachis hypogaea ) to biotic and abiotic stresses include the synthesis of prenylated stilbenoids. Members of this compound class show several protective activities in human disease studies, and the list of potential therapeutic targets continues to expand. Despite their medical and biological importance, the biosynthetic pathways of prenylated stilbenoids remain to be elucidated, and the genes encoding stilbenoid-specific prenyltransferases have yet to be identified in any plant species. In this study, we combined targeted transcriptomic and metabolomic analyses to discover prenyltransferase genes in elicitor-treated peanut hairy root cultures. Transcripts encoding five enzymes were identified, and two of these were functionally characterized in a transient expression system consisting of Agrobacterium -infiltrated leaves of Nicotiana benthamiana We observed that one of these prenyltransferases, AhR4DT-1, catalyzes a key reaction in the biosynthesis of prenylated stilbenoids, in which resveratrol is prenylated at its C-4 position to form arachidin-2, whereas another, AhR3'DT-1, added the prenyl group to C-3' of resveratrol. Each of these prenyltransferases was highly specific for stilbenoid substrates, and we confirmed their subcellular location in the plastid by fluorescence microscopy. Structural analysis of the prenylated stilbenoids suggested that these two prenyltransferase activities represent the first committed steps in the biosynthesis of a large number of prenylated stilbenoids and their derivatives in peanut. In summary, we have identified five candidate prenyltransferases in peanut and confirmed that two of them are stilbenoid-specific, advancing our understanding of this specialized enzyme family and shedding critical light onto the biosynthesis of bioactive stilbenoids., (© 2018 by The American Society for Biochemistry and Molecular Biology, Inc.)

In eukaryotes, transcriptionally inactive loci are enriched within highly condensed heterochromatin. In plants, as in mammals, the DNA of heterochromatin is densely methylated and wrapped by histones displaying a characteristic subset of post-translational modifications. Growing evidence indicates that these chromatin modifications are not sufficient for silencing. Instead, they are prerequisites for further assembly of higher-order chromatin structures that are refractory to transcription but not fully understood. We show that silencing of transposons in the pericentromeric heterochromatin of Arabidopsis thaliana requires SMC4, a core subunit of condensins I and II, acting in conjunction with CG methylation by MET1 (DNA METHYLTRANSFERASE 1), CHG methylation by CMT3 (CHROMOMETHYLASE 3), the chromatin remodeler DDM1 (DECREASE IN DNA METHYLATION 1), and histone modifications, including histone H3 Lys 27 monomethylation (H3K27me1), imparted by ATXR5 and ATXR6. SMC4/condensin also acts within the mostly euchromatic chromosome arms to suppress conditionally expressed genes involved in flowering or DNA repair, including the DNA glycosylase ROS1, which facilitates DNA demethylation. Collectively, our genome-wide analyses implicate condensin in the suppression of hundreds of loci, acting in both DNA methylation-dependent and methylation-independent pathways., (© 2017 Wang et al.; Published by Cold Spring Harbor Laboratory Press.)

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 ., (©2017 American Association for Cancer Research.)

The full exploration of gene-environment interactions requires model organisms with well-characterized ecological interactions in their natural environment, manipulability in the laboratory and genomic tools. The waterflea Daphnia magna is an established ecological and toxicological model species, central to the food webs of freshwater lentic habitats and sentinel for water quality. Its tractability and cyclic parthenogenetic life-cycle are ideal to investigate links between genes and the environment. Capitalizing on this unique model system, the STRESSFLEA consortium generated a comprehensive RNA-Seq data set by exposing two inbred genotypes of D. magna and a recombinant cross of these genotypes to a range of environmental perturbations. Gene models were constructed from the transcriptome data and mapped onto the draft genome of D. magna using EvidentialGene. The transcriptome data generated here, together with the available draft genome sequence of D. magna and a high-density genetic map will be a key asset for future investigations in environmental genomics.

Correct gene expression requires tight RNA quality control both at transcriptional and post-transcriptional levels. Using a splicing-defective allele of PASTICCINO2 (PAS2), a gene essential for plant development, we isolated suppressor mutations modifying pas2-1 mRNA profiles and restoring wild-type growth. Three suppressor of pas2 (sop) mutations modified the degradation of mis-spliced pas2-1 mRNA species, allowing the synthesis of a functional protein. Cloning of the suppressor mutations identified the core subunit of the exosome SOP2/RRP4, the exosome nucleoplasmic cofactor SOP3/HEN2 and a novel zinc-finger protein SOP1 that colocalizes with HEN2 in nucleoplasmic foci. The three SOP proteins counteract post-transcriptional (trans)gene silencing (PTGS), which suggests that they all act in RNA quality control. In addition, sop1 mutants accumulate some, but not all of the misprocessed mRNAs and other types of RNAs that are observed in exosome mutants. Taken together, our data show that SOP1 is a new component of nuclear RNA surveillance that is required for the degradation of a specific subset of nuclear exosome targets.

Eukaryotic origins of DNA replication are bound by the origin recognition complex (ORC), which scaffolds assembly of a pre-replicative complex (pre-RC) that is then activated to initiate replication. Both pre-RC assembly and activation are strongly influenced by developmental changes to the epigenome, but molecular mechanisms remain incompletely defined. We have been examining the activation of origins responsible for developmental gene amplification in Drosophila. At a specific time in oogenesis, somatic follicle cells transition from genomic replication to a locus-specific replication from six amplicon origins. Previous evidence indicated that these amplicon origins are activated by nucleosome acetylation, but how this affects origin chromatin is unknown. Here, we examine nucleosome position in follicle cells using micrococcal nuclease digestion with Ilumina sequencing. The results indicate that ORC binding sites and other essential origin sequences are nucleosome-depleted regions (NDRs). Nucleosome position at the amplicons was highly similar among developmental stages during which ORC is or is not bound, indicating that being an NDR is not sufficient to specify ORC binding. Importantly, the data suggest that nucleosomes and ORC have opposite preferences for DNA sequence and structure. We propose that nucleosome hyperacetylation promotes pre-RC assembly onto adjacent DNA sequences that are disfavored by nucleosomes but favored by ORC., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

In Arabidopsis thaliana, abundant 24 nucleotide small interfering RNAs (24 nt siRNA) guide the cytosine methylation and silencing of transposons and a subset of genes. 24 nt siRNA biogenesis requires nuclear RNA polymerase IV (Pol IV), RNA-dependent RNA polymerase 2 (RDR2) and DICER-like 3 (DCL3). However, siRNA precursors are mostly undefined. We identified Pol IV and RDR2-dependent RNAs (P4R2 RNAs) that accumulate in dcl3 mutants and are diced into 24 nt RNAs by DCL3 in vitro. P4R2 RNAs are mostly 26-45 nt and initiate with a purine adjacent to a pyrimidine, characteristics shared by Pol IV transcripts generated in vitro. RDR2 terminal transferase activity, also demonstrated in vitro, may account for occasional non-templated nucleotides at P4R2 RNA 3' termini. The 24 nt siRNAs primarily correspond to the 5' or 3' ends of P4R2 RNAs, suggesting a model whereby siRNAs are generated from either end of P4R2 duplexes by single dicing events.

Background: The mechanism by which plants synthesize and store high amounts of triacylglycerols (TAG) in tissues other than seeds is not well understood. The comprehension of controls for carbon partitioning and oil accumulation in nonseed tissues is essential to generate oil-rich biomass in perennial bioenergy crops. Persea americana (avocado), a basal angiosperm with unique features that are ancestral to most flowering plants, stores ~ 70 % TAG per dry weight in its mesocarp, a nonseed tissue. Transcriptome analyses of select pathways, from generation of pyruvate and leading up to TAG accumulation, in mesocarp tissues of avocado was conducted and compared with that of oil-rich monocot (oil palm) and dicot (rapeseed and castor) tissues to identify tissue- and species-specific regulation and biosynthesis of TAG in plants., Results: RNA-Seq analyses of select lipid metabolic pathways of avocado mesocarp revealed patterns similar to that of other oil-rich species. However, only some predominant orthologs of the fatty acid biosynthetic pathway genes in this basal angiosperm were similar to those of monocots and dicots. The accumulation of TAG, rich in oleic acid, was associated with higher transcript levels for a putative stearoyl-ACP desaturase and endoplasmic reticulum (ER)-associated acyl-CoA synthetases, during fruit development. Gene expression levels for enzymes involved in terminal steps to TAG biosynthesis in the ER further indicated that both acyl-CoA-dependent and -independent mechanisms might play a role in TAG assembly, depending on the developmental stage of the fruit. Furthermore, in addition to the expression of an ortholog of WRINKLED1 (WRI1), a regulator of fatty acid biosynthesis, high transcript levels for WRI2-like and WRI3-like suggest a role for additional transcription factors in nonseed oil accumulation. Plastid pyruvate necessary for fatty acid synthesis is likely driven by the upregulation of genes involved in glycolysis and transport of its intermediates. Together, a comparative transcriptome analyses for storage oil biosynthesis in diverse plants and tissues suggested that several distinct and conserved features in this basal angiosperm species might contribute towards its rich TAG content., Conclusions: Our work represents a comprehensive transcriptome resource for a basal angiosperm species and provides insight into their lipid metabolism in mesocarp tissues. Furthermore, comparison of the transcriptome of oil-rich mesocarp of avocado, with oil-rich seed and nonseed tissues of monocot and dicot species, revealed lipid gene orthologs that are highly conserved during evolution. The orthologs that are distinctively expressed in oil-rich mesocarp tissues of this basal angiosperm, such as WRI2, ER-associated acyl-CoA synthetases, and lipid-droplet associated proteins were also identified. This study provides a foundation for future investigations to increase oil-content and has implications for metabolic engineering to enhance storage oil content in nonseed tissues of diverse species.

As approaches are sought for more efficient and democratized uses of non-model and expanded model genomics references, ease of integration of genomic feature datasets is especially desirable in multidisciplinary research communities. Valuable conclusions are often missed or slowed when researchers refer experimental results to a single reference sequence that lacks integrated pan-genomic and multi-experiment data in accessible formats. Association of genomic positional information, such as results from an expansive variety of next-generation sequencing experiments, with annotated reference features such as genes or predicted protein binding sites, provides the context essential for conclusions and ongoing research. When the experimental system includes polymorphic genomic inputs, rapid calculation of gene structural and protein translational effects of sequence variation from the reference can be invaluable. Here we present FEATnotator, a lightweight, fast and easy to use open source software program that integrates and reports overlap and proximity in genomic information from any user-defined datasets including those from next generation sequencing applications. We illustrate use of the tool by summarizing whole genome sequence variation of a widely used natural isolate of Arabidopsis thaliana in the context of gene models of the reference accession. Previous discovery of a protein coding deletion influencing root development is replicated rapidly. Appropriate even in investigations of a single gene or genic regions such as QTL, comprehensive reports provided by FEATnotator better prepare researchers for interpretation of their experimental results. The tool is available for download at http://featnotator.sourceforge.net., (Copyright © 2015 Elsevier Inc. All rights reserved.)

In Arabidopsis (Arabidopsis thaliana), the Pseudomonas syringae effector proteins AvrB and AvrRpm1 are both detected by the RESISTANCE TO PSEUDOMONAS MACULICOLA1 (RPM1) disease resistance (R) protein. By contrast, soybean (Glycine max) can distinguish between these effectors, with AvrB and AvrRpm1 being detected by the Resistance to Pseudomonas glycinea 1b (Rpg1b) and Rpg1r R proteins, respectively. We have been using these genes to investigate the evolution of R gene specificity and have previously identified RPM1 and Rpg1b. Here, we report the cloning of Rpg1r, which, like RPM1 and Rpg1b, encodes a coiled-coil (CC)-nucleotide-binding (NB)-leucine-rich repeat (LRR) protein. As previously found for Rpg1b, we determined that Rpg1r is not orthologous with RPM1, indicating that the ability to detect both AvrB and AvrRpm1 evolved independently in soybean and Arabidopsis. The tightly linked soybean Rpg1b and Rpg1r genes share a close evolutionary relationship, with Rpg1b containing a recombination event that combined a NB domain closely related to Rpg1r with CC and LRR domains from a more distantly related CC-NB-LRR gene. Using structural modeling, we mapped polymorphisms between Rpg1b and Rpg1r onto the predicted tertiary structure of Rpg1b, which revealed highly polymorphic surfaces within both the CC and LRR domains. Assessment of chimeras between Rpg1b and Rpg1r using a transient expression system revealed that AvrB versus AvrRpm1 specificity is determined by the C-terminal portion of the LRR domain. The P. syringae effector AvrRpt2, which targets RPM1 INTERACTOR4 (RIN4) proteins in both Arabidopsis and soybean, partially blocked recognition of both AvrB and AvrRpm1 in soybean, suggesting that both Rpg1b and Rpg1r may detect these effectors via modification of a RIN4 homolog., (© 2014 American Society of Plant Biologists. All Rights Reserved.)

We have adapted the methodology of CLIP-seq (Crosslinking-Immunoprecipitation and DNA Sequencing) to map the segments of encapsidated RNAs that contact the protein shells of virions. Results from the protocol report on the RNA sequences that contact the viral capsid.

In Arabidopsis, multisubunit RNA polymerases IV and V orchestrate RNA-directed DNA methylation (RdDM) and transcriptional silencing, but what identifies the loci to be silenced is unclear. We show that heritable silent locus identity at a specific subset of RdDM targets requires HISTONE DEACETYLASE 6 (HDA6) acting upstream of Pol IV recruitment and siRNA biogenesis. At these loci, epigenetic memory conferring silent locus identity is erased in hda6 mutants such that restoration of HDA6 activity cannot restore siRNA biogenesis or silencing. Silent locus identity is similarly lost in mutants for the cytosine maintenance methyltransferase, MET1. By contrast, pol IV or pol V mutants disrupt silencing without erasing silent locus identity, allowing restoration of Pol IV or Pol V function to restore silencing. Collectively, these observations indicate that silent locus specification and silencing are separable steps that together account for epigenetic inheritance of the silenced state., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Camelina (Camelina sativa), a Brassicaceae oilseed, has received recent interest as a biofuel crop and production platform for industrial oils. Limiting wider production of camelina for these uses is the need to improve the quality and content of the seed protein-rich meal and oil, which is enriched in oxidatively unstable polyunsaturated fatty acids that are deleterious for biodiesel. To identify candidate genes for meal and oil quality improvement, a transcriptome reference was built from 2047 Sanger ESTs and more than 2 million 454-derived sequence reads, representing genes expressed in developing camelina seeds. The transcriptome of approximately 60K transcripts from 22 597 putative genes includes camelina homologues of nearly all known seed-expressed genes, suggesting a high level of completeness and usefulness of the reference. These sequences included candidates for 12S (cruciferins) and 2S (napins) seed storage proteins (SSPs) and nearly all known lipid genes, which have been compiled into an accessible database. To demonstrate the utility of the transcriptome for seed quality modification, seed-specific RNAi lines deficient in napins were generated by targeting 2S SSP genes, and high oleic acid oil lines were obtained by targeting FATTY ACID DESATURASE 2 (FAD2) and FATTY ACID ELONGASE 1 (FAE1). The high sequence identity between Arabidopsis thaliana and camelina genes was also exploited to engineer high oleic lines by RNAi with Arabidopsis FAD2 and FAE1 sequences. It is expected that these transcriptomic data will be useful for breeding and engineering of additional camelina seed traits and for translating findings from the model Arabidopsis to an oilseed crop., (© 2013 Society for Experimental Biology, Association of Applied Biologists and John Wiley & Sons Ltd.)

Background: Theobroma cacao L. cultivar Matina 1-6 belongs to the most cultivated cacao type. The availability of its genome sequence and methods for identifying genes responsible for important cacao traits will aid cacao researchers and breeders., Results: We describe the sequencing and assembly of the genome of Theobroma cacao L. cultivar Matina 1-6. The genome of the Matina 1-6 cultivar is 445 Mbp, which is significantly larger than a sequenced Criollo cultivar, and more typical of other cultivars. The chromosome-scale assembly, version 1.1, contains 711 scaffolds covering 346.0 Mbp, with a contig N50 of 84.4 kbp, a scaffold N50 of 34.4 Mbp, and an evidence-based gene set of 29,408 loci. Version 1.1 has 10x the scaffold N50 and 4x the contig N50 as Criollo, and includes 111 Mb more anchored sequence. The version 1.1 assembly has 4.4% gap sequence, while Criollo has 10.9%. Through a combination of haplotype, association mapping and gene expression analyses, we leverage this robust reference genome to identify a promising candidate gene responsible for pod color variation. We demonstrate that green/red pod color in cacao is likely regulated by the R2R3 MYB transcription factor TcMYB113, homologs of which determine pigmentation in Rosaceae, Solanaceae, and Brassicaceae. One SNP within the target site for a highly conserved trans-acting siRNA in dicots, found within TcMYB113, seems to affect transcript levels of this gene and therefore pod color variation., Conclusions: We report a high-quality sequence and annotation of Theobroma cacao L. and demonstrate its utility in identifying candidate genes regulating traits.

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

GBrowse is widely used by biologists to visualize genome annotation. We have developed WebGBrowse to facilitate hosting multiple GBrowse instances specific for different users on the same Web server. WebGBrowse automatically sets up each user-specific GBrowse instance by extracting information from user-supplied data to produce the appropriate configuration. This unit describes installation and administration of WebGBrowse for bioinformaticians who plan to manage local WebGBrowse servers in their institutions., (© 2011 by John Wiley & Sons, Inc.)

Summary: Ergatis is a flexible workflow management system for designing and executing complex bioinformatics pipelines. However, its complexity restricts its usage to only highly skilled bioinformaticians. We have developed a web-based prokaryotic genome annotation server, Integrative Services for Genomics Analysis (ISGA), which builds upon the Ergatis workflow system, integrates other dynamic analysis tools and provides intuitive web interfaces for biologists to customize and execute their own annotation pipelines. ISGA is designed to be installed at genomics core facilities and be used directly by biologists., Availability: ISGA is accessible at http://isga.cgb.indiana.edu/ and the system is also freely available for local installation.

Summary: The Generic Genome Browser (GBrowse) is one of the most widely used tools for visualizing genomic features along a reference sequence. However, the installation and configuration of GBrowse is not trivial for biologists. We have developed a web server, WebGBrowse that allows users to upload genome annotation in the GFF3 format, configure the display of each genomic feature by simply using a web browser and visualize the configured genomic features with the integrated GBrowse software., Availability: WebGBrowse is accessible via http://webgbrowse.cgb.indiana.edu/ and the system is also freely available for local installations.

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.