Anaerobic methane oxidation is quantitatively important in deeper peat layers of boreal peatlands: Evidence from anaerobic incubations, in situ stable isotopes depth profiles, and microbial communities.

Boreal peatlands store most of their carbon in layers deeper than 0.5 m under anaerobic conditions, where carbon dioxide and methane are produced as terminal products of organic matter degradation. Since the global warming potential of methane is much greater than that of carbon dioxide, the balance between the production rates of these gases is important for future climate predictions. Herein, we aimed to understand whether anaerobic methane oxidation (AMO) could explain the high CO2/CH4 anaerobic production ratios that are widely observed for the deeper peat layers of boreal peatlands. Furthermore, we quantified the metabolic pathways of methanogenesis to examine whether hydrogenotrophic methanogenesis is a dominant methane production pathway for the presumably recalcitrant deeper peat. To assess the CH4 cycling in deeper peat, we combined laboratory anaerobic incubations with a pathway-specific inhibitor, in situ depth patterns of stable isotopes in CH4, and 16S rRNA gene amplicon sequencing for three representative boreal peatlands in Western Siberia. We found up to a 69 % reduction in CH4 production due to AMO, which largely explained the high CO2/CH4 anaerobic production ratios and the in situ depth-related patterns of δ13C and δD in methane. The absence of acetate accumulation after inhibiting acetotrophic methanogenesis and the presence of sulfate- and nitrate-reducing anaerobic acetate oxidizers in the deeper peat indicated that these microorganisms use SO42- and NO3- as electron acceptors. Acetotrophic methanogenesis dominated net CH4 production in the deeper peat, accounting for 81 ± 13 %. Overall, anaerobic oxidation is quantitatively important for the methane cycle in the deeper layers of boreal peatlands, affecting both methane and its main precursor concentrations.

Spatial Changes in Microbial Communities along Different Functional Zones of a Free-Water Surface Wetland.

Constructed wetlands (CWs) are complicated ecosystems that include vegetation, sediments, and the associated microbiome mediating numerous processes in wastewater treatment. CWs have various functional zones where contrasting biochemical processes occur. Since these zones are characterized by different particle-size composition, physicochemical conditions, and vegetation, one can expect the presence of distinct microbiomes across different CW zones. Here, we investigated spatial changes in microbiomes along different functional zones of a free-water surface wetland located in Moscow, Russia. The microbiome structure was analyzed using Illumina MiSeq amplicon sequencing. We also determined particle diameter and surface area of sediments, as well as chemical composition of organic pollutants in different CW zones. Specific organic particle aggregates similar to activated sludge flocs were identified in the sediments. The highest accumulation of hydrocarbons was found in the zones with predominant sedimentation of fine fractions. Phytofilters had the highest rate of organic pollutants decomposition and predominance of Smithella, Ignavibacterium, and Methanothrix. The sedimentation tank had lower microbial diversity, and higher relative abundances of Parcubacteria, Proteiniclasticum, and Macellibacteroides, as well as higher predicted abundances of genes related to methanogenesis and methanotrophy. Thus, spatial changes in microbiomes of constructed wetlands can be associated with different types of wastewater treatment processes.

Special vulnerability of somatic niche cells to transposable element activation in Drosophila larval ovaries.

In the Drosophila ovary, somatic escort cells (ECs) form a niche that promotes differentiation of germline stem cell (GSC) progeny. The piRNA (Piwi-interacting RNA) pathway, which represses transposable elements (TEs), is required in ECs to prevent the accumulation of undifferentiated germ cells (germline tumor phenotype). The soma-specific piRNA cluster flamenco (flam) produces a substantial part of somatic piRNAs. Here, we characterized the biological effects of somatic TE activation on germ cell differentiation in flam mutants. We revealed that the choice between normal and tumorous phenotypes of flam mutant ovaries depends on the number of persisting ECs, which is determined at the larval stage. Accordingly, we found much more frequent DNA breaks in somatic cells of flam larval ovaries than in adult ECs. The absence of Chk2 or ATM checkpoint kinases dramatically enhanced oogenesis defects of flam mutants, in contrast to the germline TE-induced defects that are known to be mostly suppressed by сhk2 mutation. These results demonstrate a crucial role of checkpoint kinases in protecting niche cells against deleterious TE activation and suggest substantial differences between DNA damage responses in ovarian somatic and germ cells.

Highly photostable fluorescent labeling of proteins in live cells using exchangeable coiled coils heterodimerization.

Fluorescent proteins are commonly used to label target proteins in live cells. However, the conventional approach based on covalent fusion of targeted proteins with fluorescent protein probes is limited by the slow rate of fluorophore maturation and irretrievable loss of fluorescence due to photobleaching. Here, we report a genetically encoded protein labeling system utilizing transient interactions of small, 21-28 residues-long helical protein tags (K/E coils, KEC). In this system, a protein of interest, covalently tagged with a single coil, is visualized through binding to a cytoplasmic fluorescent protein carrying a complementary coil. The reversible heterodimerization of KECs, whose affinity can be tuned in a broad concentration range from nanomolar to micromolar, allows continuous exchange and replenishment of the tag bound to a targeted protein with the entire cytosolic pool of soluble fluorescent coils. We found that, under conditions of partial illumination of living cells, the photostability of labeling with KECs exceeds that of covalently fused fluorescent probes by approximately one order of magnitude. Similarly, single-molecule localization microscopy with KECs provided higher labeling density and allowed a much longer duration of imaging than with conventional fusion to fluorescent proteins. We also demonstrated that this method is well suited for imaging newly synthesized proteins, because the labeling efficiency by KECs is not dependent on the rate of fluorescent protein maturation. In conclusion, KECs can be used to visualize various target proteins which are directly exposed to the cytosol, thereby enabling their advanced characterization in time and space.

Novel potential causative genes in carotid paragangliomas.

BACKGROUND: Carotid paragangliomas (CPGLs) are rare neuroendocrine tumors that arise from the paraganglion at the bifurcation of the carotid artery and are responsible for approximately 65% of all head and neck paragangliomas. CPGLs can occur sporadically or along with different hereditary tumor syndromes. Approximately 30 genes are known to be associated with CPGLs. However, the genetic basis behind the development of these tumors is not fully elucidated, and the molecular mechanisms underlying CPGL pathogenesis remain unclear. METHODS: Whole exome and transcriptome high-throughput sequencing of CPGLs was performed on an Illumina platform. Exome libraries were prepared using a Nextera Rapid Capture Exome Kit (Illumina) and were sequenced under 75 bp paired-end model. For cDNA library preparation, a TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Gold (Illumina) was used; transcriptome sequencing was carried out with 100 bp paired-end read length. Obtained data were analyzed using xseq which estimates the influence of mutations on gene expression profiles allowing to identify potential causative genes. RESULTS: We identified a total of 16 candidate genes (MYH15, CSP1, MYH3, PTGES3L, CSGALNACT2, NMD3, IFI44, GMCL1, LSP1, PPFIBP2, RBL2, MAGED1, CNIH3, STRA6, SLC6A13, and ATM) whose variants potentially influence their expression (cis-effect). The strongest cis-effect of loss-of-function variants was found in MYH15, CSP1, and MYH3, and several likely pathogenic variants in these genes associated with CPGLs were predicted. CONCLUSIONS: Using the xseq probabilistic model, three novel potential causative genes, namely MYH15, CSP1, and MYH3, were identified in carotid paragangliomas.

The CIMP-high phenotype is associated with energy metabolism alterations in colon adenocarcinoma.

BACKGROUND: CpG island methylator phenotype (CIMP) is found in 15-20% of malignant colorectal tumors and is characterized by strong CpG hypermethylation over the genome. The molecular mechanisms of this phenomenon are not still fully understood. The development of CIMP is followed by global gene expression alterations and metabolic changes. In particular, CIMP-low colon adenocarcinoma (COAD), predominantly corresponded to consensus molecular subtype 3 (CMS3, "Metabolic") subgroup according to COAD molecular classification, is associated with elevated expression of genes participating in metabolic pathways. METHODS: We performed bioinformatics analysis of RNA-Seq data from The Cancer Genome Atlas (TCGA) project for CIMP-high and non-CIMP COAD samples with DESeq2, clusterProfiler, and topGO R packages. Obtained results were validated on a set of fourteen COAD samples with matched morphologically normal tissues using quantitative PCR (qPCR). RESULTS: Upregulation of multiple genes involved in glycolysis and related processes (ENO2, PFKP, HK3, PKM, ENO1, HK2, PGAM1, GAPDH, ALDOA, GPI, TPI1, and HK1) was revealed in CIMP-high tumors compared to non-CIMP ones. Most remarkably, the expression of the PKLR gene, encoding for pyruvate kinase participating in gluconeogenesis, was decreased approximately 20-fold. Up to 8-fold decrease in the expression of OGDHL gene involved in tricarboxylic acid (TCA) cycle was observed in CIMP-high tumors. Using qPCR, we confirmed the increase (4-fold) in the ENO2 expression and decrease (2-fold) in the OGDHL mRNA level on a set of COAD samples. CONCLUSIONS: We demonstrated the association between CIMP-high status and the energy metabolism changes at the transcriptomic level in colorectal adenocarcinoma against the background of immune pathway activation. Differential methylation of at least nine CpG sites in OGDHL promoter region as well as decreased OGDHL mRNA level can potentially serve as an additional biomarker of the CIMP-high status in COAD.

Bioinformatic identification of differentially expressed genes associated with prognosis of locally advanced lymph node-positive prostate cancer.

Prostate cancer (PCa) is one of the primary causes of cancer-related mortality in men worldwide. Patients with locally advanced PCa with metastases in regional lymph nodes are usually marked as a high-risk group. One of the chief concerns for this group is to make an informed decision about the necessity of conducting adjuvant androgen deprivation therapy after radical surgical treatment. During the oncogenic transformation and progression of the disease, the expression of many genes is altered. Some of these genes can serve as markers for diagnosis, predicting the prognosis or effectiveness of drug therapy, as well as possible therapeutic targets. We undertook bioinformatic analysis of the RNA-seq data deposited in The Cancer Genome Atlas consortium database to identify possible prognostic markers. We compared the groups with favorable and unfavorable prognosis for the cohort of patients with PCa showing lymph node metastasis (pT2N1M0, pT3N1M0, and pT4N1M0) and for the most common molecular type carrying the fusion transcript TMPRSS2-ERG. For the entire cohort, we revealed at least six potential markers (IDO1, UGT2B15, IFNG, MUC6, CXCL11, and GBP1). Most of these genes are involved in the positive regulation of immune response. For the TMPRSS2-ERG subtype, we also identified six genes, the expression of which may be associated with prognosis: TOB1, GALNT7, INAFM1, APELA, RAC3, and NNMT. The identified genes, after additional studies and validation in the extended cohort, could serve as a prognostic marker of locally advanced lymph node-positive PCa.

Mutational load in carotid body tumor.

BACKGROUND: Carotid body tumor (CBT) is a rare neoplasm arising from paraganglion located near the bifurcation of the carotid artery. There is great intra-tumor heterogeneity, and CBT development could be associated with both germline and somatic allelic variants. Studies on the molecular genetics of CBT are limited, and the molecular mechanisms of its pathogenesis are not fully understood. This work is focused on the estimation of mutational load (ML) in CBT. METHODS: Using the NextSeq 500 platform, we performed exome sequencing of tumors with matched lymph node tissues and peripheral blood obtained from six patients with CBT. To obtain reliable results in tumors with low ML, we developed and successfully applied a complex approach for the analysis of sequencing data. ML was evaluated as the number of somatic variants per megabase (Mb) of the target regions covered by the Illumina TruSeq Exome Library Prep Kit. RESULTS: The ML in CBT varied in the range of 0.09-0.28/Mb. Additionally, we identified several pathogenic/likely pathogenic somatic and germline allelic variants across six patients studied (including TP53 variants). CONCLUSIONS: Using the developed approach, we estimated the ML in CBT, which is much lower than in common malignant tumors. Identified variants in known paraganglioma/pheochromocytoma-causative genes and novel genes could be associated with the pathogenesis of CBT. The obtained results expand our knowledge of the mutation process in CBT as well as the biology of tumor development.

HK3 overexpression associated with epithelial-mesenchymal transition in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a common cancer worldwide. The main cause of death in CRC includes tumor progression and metastasis. At molecular level, these processes may be triggered by epithelial-mesenchymal transition (EMT) and necessitates specific alterations in cell metabolism. Although several EMT-related metabolic changes have been described in CRC, the mechanism is still poorly understood. RESULTS: Using CrossHub software, we analyzed RNA-Seq expression profile data of CRC derived from The Cancer Genome Atlas (TCGA) project. Correlation analysis between the change in the expression of genes involved in glycolysis and EMT was performed. We obtained the set of genes with significant correlation coefficients, which included 21 EMT-related genes and a single glycolytic gene, HK3. The mRNA level of these genes was measured in 78 paired colorectal cancer samples by quantitative polymerase chain reaction (qPCR). Upregulation of HK3 and deregulation of 11 genes (COL1A1, TWIST1, NFATC1, GLIPR2, SFPR1, FLNA, GREM1, SFRP2, ZEB2, SPP1, and RARRES1) involved in EMT were found. The results of correlation study showed that the expression of HK3 demonstrated a strong correlation with 7 of the 21 examined genes (ZEB2, GREM1, TGFB3, TGFB1, SNAI2, TWIST1, and COL1A1) in CRC. CONCLUSIONS: Upregulation of HK3 is associated with EMT in CRC and may be a crucial metabolic adaptation for rapid proliferation, survival, and metastases of CRC cells.

Exome analysis of carotid body tumor.

BACKGROUND: Carotid body tumor (CBT) is a form of head and neck paragangliomas (HNPGLs) arising at the bifurcation of carotid arteries. Paragangliomas are commonly associated with germline and somatic mutations involving at least one of more than thirty causative genes. However, the specific functionality of a number of these genes involved in the formation of paragangliomas has not yet been fully investigated. METHODS: Exome library preparation was carried out using Nextera® Rapid Capture Exome Kit (Illumina, USA). Sequencing was performed on NextSeq 500 System (Illumina). RESULTS: Exome analysis of 52 CBTs revealed potential driver mutations (PDMs) in 21 genes: ARNT, BAP1, BRAF, BRCA1, BRCA2, CDKN2A, CSDE1, FGFR3, IDH1, KIF1B, KMT2D, MEN1, RET, SDHA, SDHB, SDHC, SDHD, SETD2, TP53BP1, TP53BP2, and TP53I13. In many samples, more than one PDM was identified. There are also 41% of samples in which we did not identify any PDM; in these cases, the formation of CBT was probably caused by the cumulative effect of several not highly pathogenic mutations. Estimation of average mutation load demonstrated 6-8 mutations per megabase (Mb). Genes with the highest mutation rate were identified. CONCLUSIONS: Exome analysis of 52 CBTs for the first time revealed the average mutation load for these tumors and also identified potential driver mutations as well as their frequencies and co-occurrence with the other PDMs.

Upregulation of NETO2 gene in colorectal cancer.

BACKGROUND: Neuropilin and tolloid-like 2 (NETO2) is a single-pass transmembrane protein that has been shown primarily implicated in neuron-specific processes. Upregulation of NETO2 gene was also detected in several cancer types. In colorectal cancer (CRC), it was associated with tumor progression, invasion, and metastasis, and seems to be involved in epithelial-mesenchymal transition (EMT). However, the mechanism of NETO2 action is still poorly understood. RESULTS: We have revealed significant increase in the expression of NETO2 gene and deregulation of eight EMT-related genes in CRC. Four of them were upregulated (TWIST1, SNAIL1, LEF1, and FOXA2); the mRNA levels of other genes (FOXA1, BMP2, BMP5, and SMAD7) were decreased. Expression of NETO2 gene was weakly correlated with that of genes involved in the EMT process. CONCLUSIONS: We found considerable NETO2 upregulation, but no significant correlation between the expression of NETO2 and EMT-related genes in CRC. Thus, NETO2 may be involved in CRC progression, but is not directly associated with EMT.

Microbial Community Structure of Activated Sludge in Treatment Plants with Different Wastewater Compositions.

Activated sludge (AS) plays a crucial role in the treatment of domestic and industrial wastewater. AS is a biocenosis of microorganisms capable of degrading various pollutants, including organic compounds, toxicants, and xenobiotics. We performed 16S rRNA gene sequencing of AS and incoming sewage in three wastewater treatment plants (WWTPs) responsible for processing sewage with different origins: municipal wastewater, slaughterhouse wastewater, and refinery sewage. In contrast to incoming wastewater, the taxonomic structure of AS biocenosis was found to become stable in time, and each WWTP demonstrated a unique taxonomic pattern. Most pathogenic microorganisms (Streptococcus, Trichococcus, etc.), which are abundantly represented in incoming sewage, were significantly decreased in AS of all WWTPs, except for the slaughterhouse wastewater. Additional load of bioreactors with influent rich in petroleum products and organic matter was associated with the increase of bacteria responsible for AS bulking and foaming. Here, we present a novel approach enabling the prediction of the metabolic potential of bacterial communities based on their taxonomic structures and MetaCyc database data. We developed a software application, XeDetect, to implement this approach. Using XeDetect, we found that the metabolic potential of the three bacterial communities clearly reflected the substrate composition. We revealed that the microorganisms responsible for AS bulking and foaming (most abundant in AS of slaughterhouse wastewater) played a leading role in the degradation of substrates such as fatty acids, amino acids, and other bioorganic compounds. Moreover, we discovered that the chemical, rather than the bacterial composition of the incoming wastewater was the main factor in AS structure formation. XeDetect (freely available: https://sourceforge.net/projects/xedetect) represents a novel powerful tool for the analysis of the metabolic capacity of bacterial communities. The tool will help to optimize bioreactor performance and avoid some most common technical problems.