Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline.

Summary: DNA sequencing is becoming more affordable and faster through advances in high-throughput technologies. This rise in data availability has contributed to the development of novel algorithms to elucidate previously obscure features and led to an increased reliance on complex workflows to integrate such tools into analyses pipelines. To facilitate the analysis of DNA sequencing data, we created metapipeline-DNA, a highly configurable and extensible pipeline. It encompasses a broad range of processing including raw sequencing read alignment and recalibration, variant calling, quality control and subclonal reconstruction. Metapipeline-DNA also contains configuration options to select and tune analyses while being robust to failures. This standardizes and simplifies the ability to analyze large DNA sequencing in both clinical and research settings. Availability: Metapipeline-DNA is an open-source Nextflow pipeline under the GPLv2 license and is freely available at https://github.com/uclahs-cds/metapipeline-DNA.

Moxibustion pre-treatment attenuates seizure severity during status epilepticus and counteracts the proconvulsant function of the purinergic P2X7 receptor.

Moxibustion, traditional Chinese medicine treatment, involves the warming of specific acupuncture points of the body using ignited herbal materials. Evidence suggests beneficial effects of moxibustion in several brain diseases including epilepsy, however, whether moxibustion pretreatment impacts on seizures and what are the underlying mechanisms remains to be established. Evidence has suggested the purinergic ATP-gated P2X7 receptor (P2X7R) to be involved in the actions of moxibustion. Moreover, P2X7R signalling is now well established to contribute to long-lasting brain hyperexcitability underlying epilepsy development. Whether P2X7R signalling is involved in the seizure-reducing actions of moxibustion has not been investigated to date. For our studies we used C57BL/6 male mice that received moxibustion pre-treatments at the acupoints Zusanli (ST36) and Dazhui (GV14) once daily for either 7, 14, or 21 days. This was followed by an intraperitoneal injection of kainic acid (KA, 30 mg/kg) to induce status epilepticus. Behavioral changes during KA-induced status epilepticus were analyzed according to the Racine scale. Changes in electrographic seizures were analyzed via cortical implanted electroencephalogram (EEG) electrodes. While no effect on seizure severity was observed following 7 days of moxibustion pre-treatment, moxibustion pre-treatment at both ST36 and GV14 for 14 or 21 days significantly reduced KA-induced behavior seizures at a similar rate. Cortical EEG recordings showed that 14 days of moxibustion pre-treatments also reduced electrographic seizures, confirming the anticonvulsant actions of moxibustion pre-treatment. To determine whether moxibustion impacts the pro-convulsant actions of P2X7R signaling, mice were treated with the P2X7R agonist BzATP or P2X7R antagonist A438079. While treatment with the P2X7R agonist BzATP exacerbated seizure severity, treatment with the P2X7R antagonist reduced seizure severity. We further found that moxibustion pre-treatment attenuated epileptic seizures by counteracting the effects of BzATP. These results suggest that moxibustion pre-treatment at the acupoints ST36 and GV14 for 14 days has anti-epileptic effects, which may counteract the proconvulsant functions of the P2X7R.

Upgrading of Grade Group 1 Prostate Cancer at Prostatectomy: Germline Risk Factors in a Prospective Cohort.

BACKGROUND: Localized prostate tumors show significant spatial heterogeneity, with regions of high-grade disease adjacent to lower-grade disease. Consequently, prostate cancer biopsies are prone to sampling bias, potentially leading to underestimation of tumor grade. To study the clinical, epidemiologic and molecular hallmarks of this phenomenon, we conducted a prospective study of grade upgrading: differences in detected prostate cancer grade between biopsy and surgery. METHODS: We established a prospective, multi-institutional cohort of men with Grade Group 1 (GG1) prostate cancer on biopsy who underwent radical prostatectomy. Upgrading was defined as detection of GG2+ in the resected tumor. Germline DNA from 192 subjects was subjected to whole-genome sequencing to quantify ancestry, pathogenic variants in DNA damage response genes and polygenic risk. RESULTS: Of 285 men, 67% upgraded at surgery. PSA density and percent of cancer in pre-prostatectomy positive biopsy cores were significantly associated with upgrading. No assessed genetic risk factor was predictive of upgrading, including polygenic risk scores for prostate cancer diagnosis. CONCLUSIONS: In a cohort of low-grade prostate cancer patients, a majority upgraded at radical prostatectomy. PSA density and percent of cancer in pre-prostatectomy positive biopsy cores portended the presence of higher-grade disease, while germline genetics was not informative in this setting. Patients with low-risk prostate cancer, but elevated PSA density or percent cancer in positive biopsy cores, may benefit from repeat biopsy, additional imaging or other approaches to complement active surveillance. IMPACT: Further risk stratification of patients with low-risk prostate cancer may provide useful context for active surveillance decision-making.

Habitat changes due to glacial freezing and melting reshape microbial networks.

The phenomenon of glacial freezing and thawing involves microbial sequestration, release, and colonization, which has the potential to impact ecosystem functioning through changes in microbial diversity and interactions. In this study, we examined the structural features of microbial communities of the Dongkemadi glacier, including bacteria, fungi, and archaea, in four distinct glacial environments (snow, ice, meltwater, and frontier soil). The sequestration, release, and colonization of glacial microbes have been found to significantly impact the diversity and structure of glacial microbial communities, as well as the complexity of microbial networks. Specifically, the complexity of bacterial networks has been observed to increase in a sequential manner during these processes. Utilizing the Inter-Domain Ecological Network approach, researchers have further explored the cross-trophic interactions among bacteria, fungi, and archaea. The complexity of the bacteria-fungi-archaea network exhibited a sequential increase due to the processes of sequestration, release, and colonization of glacial microbes. The release and colonization of glacial microbes led to a shift in the role of archaea as key species within the network. Additionally, our findings suggest that the hierarchical interactions among various microorganisms contributed to the heightened complexity of the bacteria-fungi-archaea network. The primary constituents of the glacial microbial ecosystem are unclassified species associated with the Polaromonas. It is noteworthy that various key species in glacial ecosystems are influenced by the distinct environmental factors. Moreover, our findings suggest that key species are not significantly depleted in response to abrupt alterations in individual environmental factors, shedding light on the dynamics of microbial cross-trophic interactions within glacial ecosystems.

Bayesian analysis of joint quantile regression for multi-response longitudinal data with application to primary biliary cirrhosis sequential cohort study.

This article proposes a Bayesian approach for jointly estimating marginal conditional quantiles of multi-response longitudinal data with multivariate mixed effects model. The multivariate asymmetric Laplace distribution is employed to construct the working likelihood of the considered model. Penalization priors on regression parameters are incorporated into the working likelihood to conduct Bayesian high-dimensional inference. Markov chain Monte Carlo algorithm is used to obtain the fully conditional posterior distributions of all parameters and latent variables. Monte Carlo simulations are conducted to evaluate the sample performance of the proposed joint quantile regression approach. Finally, we analyze a longitudinal medical dataset of the primary biliary cirrhosis sequential cohort study to illustrate the real application of the proposed modeling method.

Saxibacter everestensis gen. nov., sp. nov., A Novel Member of the Family Brevibacteriaceae, Isolated from the North Slope of Mount Everest.

We isolated and analyzed a novel, Gram-stain-positive, aerobic, rod-shaped, non-motile actinobacterium, designated as strain ZFBP1038T, from rock sampled on the north slope of Mount Everest. The growth requirements of this strain were 10-37 °C, pH 4-10, and 0-6% (w/v) NaCl. The sole respiratory quinone was MK-9, and the major fatty acids were anteiso-C15:0 and iso-C17:0. Peptidoglycan containing meso-diaminopimelic acid, ribose, and glucose were the major cell wall sugars, while polar lipids included diphosphatidyl glycerol, phosphatidyl glycerol, an unidentified phospholipid, and an unidentified glycolipid. A phylogenetic analysis based on 16S rRNA gene sequences showed that strain ZFBP1038T has the highest similarity with Spelaeicoccus albus DSM 26341 T (96.02%). ZFBP1038T formed a distinct monophyletic clade within the family Brevibacteriaceae and was distantly related to the genus Spelaeicoccus. The G + C content of strain ZFBP1038T was 63.65 mol% and the genome size was 4.05 Mb. Digital DNA-DNA hybridization, average nucleotide identity, and average amino acid identity values between the genomes of strain ZFBP1038T and representative reference strains were 19.3-25.2, 68.0-71.0, and 52.8-60.1%, respectively. Phylogenetic, phenotypic, and chemotaxonomic characteristics as well as comparative genome analyses suggested that strain ZFBP1038T represents a novel species of a new genus, for which the name Saxibacter gen. nov., sp. nov. was assigned with the type strain Saxibacter everestensis ZFBP1038T (= EE 014 T = GDMCC 1.3024 T = JCM 35335 T).

The genomic and evolutionary landscapes of anaplastic thyroid carcinoma.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

NFTest: automated testing of Nextflow pipelines.

MOTIVATION: The ongoing expansion in the volume of biomedical data has contributed to a growing complexity in the tools and technologies used in research with an increased reliance on complex workflows written in orchestration languages such as Nextflow to integrate algorithms into processing pipelines. The growing use of workflows involving various tools and algorithms has led to increased scrutiny of software development practices to avoid errors in individual tools and in the connections between them. RESULTS: To facilitate test-driven development of Nextflow pipelines, we created NFTest, a framework for automated pipeline testing and validation with customizability options for Nextflow features. It is open-source, easy to initialize and use, and customizable to allow for testing of complex workflows with test success configurable through a broad range of assertions. NFTest simplifies the testing burden on developers by automating tests once defined and providing a flexible interface for running tests to validate workflows. This reduces the barrier to rigorous biomedical workflow testing and paves the way toward reducing computational errors in biomedicine. AVAILABILITY AND IMPLEMENTATION: NFTest is an open-source Python framework under the GPLv2 license and is freely available at https://github.com/uclahs-cds/tool-NFTest. The call-sSNV Nextflow pipeline is available at: https://github.com/uclahs-cds/pipeline-call-sSNV.

Role of TRAK1 variants in epilepsy: genotype-phenotype analysis in a pediatric case of epilepsy with developmental disorder.

Purpose: The TRAK1 gene is mapped to chromosome 3p22.1 and encodes trafficking protein kinesin binding 1. The aim of this study was to investigate the genotype-phenotype of TRAK1-associated epilepsy. Methods: Trio-based whole-exome sequencing was performed on a cohort of 98 patients with epilepsy of unknown etiologies. Protein modeling and the VarCards database were used to predict the damaging effects of the variants. Detailed neurological phenotypes of all patients with epilepsy having TRAK1 variants were analyzed to assess the genotype-phenotype correlations. Results: A novel TRAK1 compound heterozygous variant comprising variant c.835C > T, p.Arg279Cys and variant c.2560A > C, p.Lys854Gln was identified in one pediatric patient. Protein modeling and VarCards database analyses revealed that the variants were damaging. The patient received a diagnosis of early infantile epileptic spasms with a developmental disorder; he became seizure-free through valproate and adrenocorticotropic hormone treatment. Further results for six variants in 12 patients with epilepsy indicated that biallelic TRAK1 variants (including homozygous or compound heterozygous variants) were associated with epilepsy with developmental disorders. Among these patients, eight (67%) had epileptic spasms and seven (58%) were intractable to anti-seizure medicines. Moreover, eight patients experienced refractory status epilepticus, of which seven (88%) died in early life. To our knowledge, this is the first reported case of epilepsy caused by TRAK1 compound heterozygous variants. Conclusion: Biallelic TRAK1 variants can cause epilepsy and developmental disorders. In these patients, seizures progress to status epilepticus, suggesting a high risk for poor outcomes and the requirement of early treatment.

Investigation of FRMPD4 variants associated with X-linked epilepsy.

BACKGROUND: The etiology of unexplained epilepsy in most patients remains unclear. Variants of FRMPD4 are suggested to be associated with neurodevelopmental disorders. Therefore, we screened for disease-causing FRMPD4 variants in patients with epilepsy. METHODS: Trios-based whole-exome sequencing was conducted on a cohort of 85 patients with unexplained epilepsy, their parents, and extended family members. Additional cases with FRMPD4 variants were identified from the China Epilepsy Gene Matching Platform V.1.0. The frequency of variants was analyzed, and their subregional effects were predicted using in silico tools. The genotype-phenotype correlation of the newly defined causative genes and protein stability were analyzed using I-Mutant V.3.0 and Grantham scores. RESULTS: Two novel missense variants of FRMPD4 were identified in two families. Using the gene matching platform, we identified three additional novel missense variants. These variants presented at low or no allele frequencies in the gnomAD database. All the variants were located outside the three FRMPD4 main domains (WW, PDZ, and FERM). In silico analyses revealed that the variants were damaging and were predicted to be the least stable. All patients eventually became seizure-free. Eight of the 21 patients with FRMPD4 variants had epilepsy, of which five (63%) had missense variants located outside the domains, two had deletions involving exon 2, and one had a frameshift variant located outside the domains. Patients with epilepsy caused by missense variants were often free of intellectual disabilities (4/5), whereas patients with epilepsy caused by truncated variants had intellectual disabilities and structural brain abnormalities (3/3). CONCLUSIONS: The FRMPD4 gene is potentially associated with epilepsy. The genotype-phenotype correlation of FRMPD4 variants indicated that differences in variant types and locations of FRMPD4 may explain their phenotypic variation.

Prostate lineage-specific metabolism governs luminal differentiation and response to antiandrogen treatment.

Lineage transitions are a central feature of prostate development, tumourigenesis and treatment resistance. While epigenetic changes are well known to drive prostate lineage transitions, it remains unclear how upstream metabolic signalling contributes to the regulation of prostate epithelial identity. To fill this gap, we developed an approach to perform metabolomics on primary prostate epithelial cells. Using this approach, we discovered that the basal and luminal cells of the prostate exhibit distinct metabolomes and nutrient utilization patterns. Furthermore, basal-to-luminal differentiation is accompanied by increased pyruvate oxidation. We establish the mitochondrial pyruvate carrier and subsequent lactate accumulation as regulators of prostate luminal identity. Inhibition of the mitochondrial pyruvate carrier or supplementation with exogenous lactate results in large-scale chromatin remodelling, influencing both lineage-specific transcription factors and response to antiandrogen treatment. These results establish reciprocal regulation of metabolism and prostate epithelial lineage identity.

A novel UV-resistant bacterium Sphingomonas endolithica sp. nov., and genomic analysis, isolated from the north slope of Mount Everest.

Gram-stain-negative, aerobic, rod-shaped, non-motile bacterium strain ZFBP2030T was isolated from a rock on the North slope of Mount Everest. This strain contained a unique ubiquinone-10 (Q-10) as a predominant respiratory quinone. Among the tested fatty acids, the strain contained summed feature 8, C14:0 2OH, and C16:0, as major cellular fatty acids. The polar lipid profile contained phosphatidyl glycerol, phosphatidyl ethanolamine, three unidentified phospholipids, two unidentified aminolipids, and six unidentified lipids. The cell-wall peptidoglycan was a meso-diaminopimelic acid, and cell-wall sugars were ribose and galactose. Phylogenetic analyses based on 16S rRNA gene sequence revealed that strain ZFBP2030T was a member of the genus Sphingomonas, exhibiting high sequence similarity to the 16S rRNA gene sequences of Sphingomonas aliaeris DH-S5T (97.9%), Sphingomonas alpina DSM 22537T (97.3%) and Sphingomonas hylomeconis CCTCC AB 2013304T (97.0%). The 16S rRNA gene sequence similarity between ZFBP2030T and other typical strains was less than 97.0%. The average amino acid identity values, average nucleotide identity, and digital DNA-DNA hybridization values between strain ZFBP2030T and its highest sequence similarity strains were 56.9-79.9%, 65.1-82.2%, and 19.3-25.8%, respectively. The whole-genome size of the novel strain ZFBP2030T was 4.1 Mbp, annotated with 3838 protein-coding genes and 54 RNA genes. Moreover, DNA G + C content was 64.7 mol%. Stress-related functions predicted in the subsystem classification of the strain ZFBP2030T genome included osmotic, oxidative, cold/heat shock, detoxification, and periplasmic stress responses. The overall results of this study clearly showed that strain ZFBP2030T is a novel species of the genus Sphingomonas, for which the name Sphingomonas endolithica sp. nov. is proposed. The type of strain is ZFBP2030T (= EE 013T = GDMCC 1.3123T = JCM 35386T).

Spatio-temporal variation and environmental drivers of chlorophyll a concentration and diatom community in four small urban lakes of Kunming, China.

Chlorophyll a (Chla) and diatom community structure are two indicators of lake water quality. In this study, we investigated the environmental parameters, chlorophyll a, and diatom community of four small urban lakes in Kunming (Beitan, Beihu, Nanhu and Longtan lakes in the campus of Yunnan Normal University) between March 2017 and December 2019. The results showed that the concentrations of total nitrogen (TN), total phosphorus (TP), and Chla in the four lakes showed significant seasonal fluctuation. The Chla concentration in Nanhu Lake, which had the highest nutrient level among the four lakes, was significantly higher than that in the other three lakes and largely affected by TN. In comparison, water temperature significantly contributed to the increases in Chla concentration in the other three lakes. Water temperature and TN were significantly correlated with Chla across the four lakes. Diatom assemblages in Beitan, Nanhu, and Longtan lakes were dominated by planktonic diatoms, and benthic diatoms were dominant in the shallowest lake Beihu, suggesting that water depth significantly affected the proportion of planktonic diatoms and dominant taxa. Water depth, TN, TP, transparency, and water temperature affected the spatio-temporal changes of diatom community structure, with water temperature as the major factor in causing the seasonal variation in diatom community, and TN and TP as the major drivers for community variation among lakes within the same season.

STING is a cell-intrinsic metabolic checkpoint restricting aerobic glycolysis by targeting HK2.

Evasion of antitumour immunity is a hallmark of cancer. STING, a putative innate immune signalling adaptor, has a pivotal role in mounting antitumour immunity by coordinating innate sensing and adaptive immune surveillance in myeloid cells. STING is markedly silenced in various human malignancies and acts as a cell-intrinsic tumour suppressor. How STING exerts intrinsic antitumour activity remains unclear. Here, we report that STING restricts aerobic glycolysis independent of its innate immune function. Mechanistically, STING targets hexokinase II (HK2) to block its hexokinase activity. As such, STING inhibits HK2 to restrict tumour aerobic glycolysis and promote antitumour immunity in vivo. In human colorectal carcinoma samples, lactate, which can be used as a surrogate for aerobic glycolysis, is negatively correlated with STING expression level and antitumour immunity. Taken together, this study reveals that STING functions as a cell-intrinsic metabolic checkpoint that restricts aerobic glycolysis to promote antitumour immunity. These findings have important implications for the development of STING-based therapeutic modalities to improve antitumour immunotherapy.

CD40LG-associated X-linked Hyper-IgM Syndrome (XHIGM) with pulmonary alveolar proteinosis: a case report.

BACKGROUND: D40LG-associated X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis has rarely been reported, and its genotype-phenotypic correlation remains elusive. CASE PRESENTATION: We describe a five-month-old boy with CD40LG mutation (c.516T > A, p.Tyr172Ter) X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis as the first manifestation. The patient completely recovered after immunotherapy and allogeneic hematopoietic stem cell transplantation. In addition, four previously reported patients with CD40LG mutation with pulmonary alveolar proteinosis were also analyzed. All of these patients presented with early onset of pulmonary infections and a good response to immunotherapy. The structural model of CD40LG indicated that all mutations caused the X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis to be located within the tumor necrosis factor homology domain. CONCLUSIONS: A case was presented, and the characteristics of four cases of CD40LG-associated X-linked hyper-IgM syndrome with pulmonary alveolar proteinosis were summarized. The variant locations may explain the phenotypic heterogeneity of patients with the CD40LG mutation.

The increase of Nrf2 m6A modification induced by FTO downregulation promotes hippocampal neuron injury and aggravates the progression of epilepsy in a rat model.

Epilepsy is a common chronic neurological disorder characterized by widespread neuronal death. The purpose of this study was to investigate the role of nuclear factor erythroid 2-related factor 2 (Nrf2) m6A methylation in epilepsy. To create epileptic models, the rats were given Lithium chloride and pilocarpine, and isolated primary rat hippocampal neurons were cultured in an Mg2+ -free medium. The frequency of seizures was recorded in the epilepsy group of rats. The functional tests included TUNEL, MTT, and flow cytometry. Mechanistically, RNA degradation assay, RNA immunoprecipitation, and methylated RNA immunoprecipitation were performed. In epileptic models, Nrf2 and fat mass and obesity-associated (FTO) levels were downregulated, whereas YT521-B homology (YTH) domain family protein 2 (YTHDF2) was upregulated. Additionally, in epileptic models, there was a rise in the m6A methylation level of Nrf2 mRNA. Overexpressing FTO increased cell viability and reduced apoptosis, but Nrf2 interference reversed these effects. Meanwhile, FTO overexpression decreased the m6A methylation of Nrf2 mRNA. Moreover, YTHDF2 bound to Nrf2 mRNA and decreased its stability. Furthermore, FTO overexpression reduced seizure frequency in rats and inhibited hippocampal neuron apoptosis via lowering the m6A methylation level of Nrf2 mRNA. Overexpressing FTO reduced m6A methylation of Nrf2 mRNA, increased cell viability, suppressed apoptosis, and slowed the progression of epileptic diseases, which is linked to YTHDF2 binding to m6A-modified Nrf2 and promoting its degradation, as well as downregulating Nrf2 expression in hippocampal neurons.

Anti-MDA5 antibody dermatomyositis-associated rapidly progressive interstitial lung disease patient complicated with mixed connective tissue disease: A case report.

Anti-MDA5 antibody dermatomyositis (DM) is a special type of myositis, which can potentially cause rapidly progressive interstitial lung disease (RP-ILD). Mixed connective tissue disease (MCTD) is a complex disease with different characteristics of autoimmune connective tissue disease, associated with ILD. Both are rare diseases, and few patients with both diseases have been reported. A 71-year-old woman complained of palpitations, with a 2 months history of rash around her hands, extensor surface of right elbow, and the nape of her neck. Subsequently, the patient had acute exacerbation of dyspnea and tachypnea. Anti-Ro52, U1 RNP and MDA5 antibodies were positive; the presenting evidence was suggestive of anti-MDA5+ DM-RP-ILD complicated with MCTD. Our patient deteriorated rapidly and had a fatal outcome, despite "triple therapy" for RP-ILD. This case illustrates that patients with coexisting anti-MDA5+ DM and MCTD have the former's typical clinical manifestations, and may develop ILD quickly rather than slowly as in MCTD, especially with the coexistence of anti-Ro52 antibodies.

Detection and Quantification of Candidatus Methanoperedens-Like Archaea in Freshwater Wetland Soils.

Candidatus Methanoperedens-like archaea, which can use multiple electron acceptors (nitrate, iron, manganese, and sulfate) for anaerobic methane oxidation, could play an important role in reducing methane emissions from freshwater wetlands. Currently, very little is known about the distribution and community composition of Methanoperedens-like archaea in freshwater wetlands, particularly based on their alpha subunit of methyl-coenzyme M reductase (mcrA) genes. Here, the community composition, diversity, and abundance of Methanoperedens-like archaea were investigated in a freshwater wetland through high-throughput sequencing and quantitative PCR on their mcrA genes. A large number of Methanoperedens-like mcrA gene sequences (119,250) were recovered, and a total of 31 operational taxonomic units (OTUs) were generated based on 95% sequence similarity cut-off. The majority of Methanoperedens-like sequences can be grouped into three distinct clusters that were closely associated with the known Methanoperedens species which can couple anaerobic methane oxidation to nitrate or iron reduction. The community composition of Methanoperedens-like archaea differed significantly among different sampling sites, and their mcrA gene abundance was 1.49 × 106 ~ 4.62 × 106 copies g-1 dry soil in the examined wetland. In addition, the community composition of Methanoperedens-like archaea was significantly affected by the soil water content, and the archaeal abundance was significantly positively correlated with the water content. Our results suggest that the mcrA gene is a good biomarker for detection and quantification of Methanoperedens-like archaea, and provide new insights into the distribution and environmental regulation of these archaea in freshwater wetlands.

Effect of gradual increase of atmospheric CO2 concentration on nitrification potential and communities of ammonia oxidizers in paddy fields.

Currently, the influence of elevated atmospheric CO2 concentration (eCO2) on ammonia oxidation to nitrite, the rate-limiting step of nitrification in paddy soil, is poorly known. Previous studies that simulate the effect of eCO2 on nitrification are primarily based on an abrupt increase of atmospheric CO2 concentration. However, paddy ecosystems are experiencing a gradual increase of CO2 concentration. To better understand how the nitrification potential, abundance and communities of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) respond to eCO2 in paddy ecosystems, a field experiment was conducted using the following two treatments: a gradual increase of CO2 (EC, increase of 40 ppm per year until 200 ppm above ambient) and ambient CO2 (CK). The results demonstrated that the EC treatment significantly (P < 0.05) stimulated the soil potential nitrification rate (PNR) at the jointing and milky stages, which increased by 127.83% and 27.35%, respectively, compared with CK. Furthermore, the EC treatment significantly (P < 0.05) stimulated the AOA and AOB abundance by 56.60% and 133.84%, respectively, at the jointing stage. Correlation analysis showed that the PNR correlated well with the abundance of AOB (R2 = 0.7389, P < 0.001). In addition, the EC treatment significantly (P < 0.05) altered the community structure of AOB, while it had little effect on that of AOA. A significant difference in the proportion of Nitrosospira was observed between CO2 treatments. In conclusion, the gradual increase of CO2 positively influenced the PNR and abundance of ammonia oxidizers, and AOB could be more important than AOA in nitrification under eCO2.