Changes in Mutations of Cell-Free DNA and Liver Tumor Tissue in Patients with Advanced Hepatocellular Carcinoma before and after Introduction of Lenvatinib.

INTRODUCTION: Cell-free DNA (cfDNA) is expected to contribute to the decision for treatment and prediction of effects with minimally invasion. We investigated the correlation between gene mutations before and after lenvatinib (LEN) treatment and its effectiveness, in order to find advanced hepatocellular carcinoma (HCC) patients who would benefit greatly from the therapy. METHODS: We analyzed cfDNA before and 6-8 weeks after the start of treatment in 20 advanced HCC patients who started LEN. A next-generation sequencer was used for CTNNB1 and TP53. Concerning TERT promoter, -124C>T and -146C>T mutations are researched using digital PCR. In addition, we examined liver tumor biopsy tissues by the same method. Computerized tomography evaluation was performed at 6-8 weeks and 3-4 months to assess the efficacy. RESULTS: Frequencies of TERT promoter, CTNNB1, and TP53 mutations in pretreatment cfDNA were 45%, 65%, and 65%, but 53%, 41%, and 47% in HCC tissues, respectively. There were no clear correlations between these gene mutations and the disease-suppressing effect or progression-free survival. Overall, there were many cases showing a decrease in mutations after LEN treatment. Integrating the reduction of CTNNB1 and TP53 genetic mutations increased the potential for disease suppression. CONCLUSION: This study suggests that analysis of cfDNA in advanced HCC patients may be useful for identifying LEN responders and determining therapeutic efficacy. Furthermore, it has potential for selecting responders for other molecular-targeted drugs.

Simple tape-stripping method for highly reliable and quantitative analysis of skin microbiome.

The composition of human skin microbiome profoundly impacts host skin health and disease. However, the relationship between skin homeostasis or the development of skin diseases and daily changes in skin microbial composition is poorly understood. Longitudinal samplings at more frequent intervals would address this issue, while conventional sampling methods have technical difficulties, leading to limitations in sampling opportunities. Here, we developed a simple and stable tape-stripping method regardless of the operator's skill. Our method enables skin microbial sampling within 30 seconds and taking multiple skin microbial samples from the same body site. The amount of microbial DNA among multiple sampling sites could be measured within 13.5%. The sequencing results of multiple sampling showed high consistency, Pearson's correlation coefficient between multiple samples of 0.98. Furthermore, these results were comparable to those collected by the conventional swabbing method. These results demonstrate that our tape-stripping method enables simple microbiome collection and highly reliable quantitative skin microbiome analysis. These features of our method would lead to a further understanding of skin disease development or diagnosis of skin conditions in clinical research by increasing the opportunities for microbial sampling.

Analyses of tertiary lymphoid structures observed in cases of Merkel cell carcinoma showing spontaneous regression.

Merkel cell carcinoma (MCC) is a high-grade skin cancer, but spontaneous regression is observed at a markedly higher frequency than in other carcinomas. Although spontaneous regression is a phenomenon that greatly impacts treatment planning, we still cannot predict it. We previously reported on the prognostic impact of the presence or absence of tertiary lymphoid structures (TLS) and of Merkel cell polyomavirus (MCPyV) infection. To learn more about the spontaneous regression of MCC, detailed analyses were performed focusing on spontaneous regression cases. We collected 71 Japanese patients with MCC including 6 cases of spontaneous regression. Samples were analysed by immunostaining, spatial single-cell analysis using PhenoCycler, and RNA sequencing using the next-generation sequencer (NGS). All 6 cases of spontaneous regression were positive for MCPyV. TLS was positive in all 5 cases analysed. Spatial single-cell analyses revealed that PD-L1-positive tumour cells were in close proximity to CD20-positive B cell and CD3-, 4-positive T cells. Gene set enrichment analysis between MCPyV-positive and TLS-positive samples and other samples showed significantly high enrichment of "B-cell-mediated immunity" gene sets in the MCPyV-positive and TLS-positive groups. In conclusion, TLS may play an important role in the spontaneous regression of MCC.

Noninvasive graft monitoring using donor-derived cell-free DNA in Japanese liver transplantation.

AIM: To evaluate the use of donor-derived cell-free DNA (dd-cfDNA) in diagnosing graft injuries in Japanese liver transplantation (LTx), including family-related living donors. METHODS: A total of 321 samples from 10 newly operated LTx recipients were collected to monitor the early dynamics of dd-cfDNA levels after LTx. Fifty-five samples from 55 recipients were collected during protocol biopsies (PB), whereas 36 samples from 27 recipients were collected during event biopsies, consisting of 11 biopsy-proven acute rejection (AR), 20 acute dysfunctions without rejection (ADWR), and 5 chronic rejections. The levels of dd-cfDNA were quantified using a next-generation sequencer based on single nucleotide polymorphisms. RESULTS: The dd-cfDNA levels were elevated significantly after LTx, followed by a rapid decline to the baseline in patients without graft injury within 30 days post-LTx. The dd-cfDNA levels were significantly higher in the 11 samples obtained during AR than those obtained during PB (p < 0.0001), which decreased promptly after treatment. The receiver operator characteristic curve analysis of diagnostic ability yielded areas under the curve of 0.975 and 0.897 for AR (rejection activity index [RAI] ≥3) versus PB and versus non-AR (ADWR + PB). The dd-cfDNA levels during AR were elevated earlier and correlated more strongly with the RAI (r = 0.740) than aspartate aminotransferase/alanine aminotransferase. The dd-cfDNA levels were neither associated with graft fibrosis based on histology nor the status of donor-specific antibodies in PB samples. CONCLUSIONS: Donor-derived cell-free DNA serves as a sensitive biomarker for detecting graft injuries in LTx. Further large-scale cohort studies are warranted to optimize its use in differentiating various post-LTx etiologies.

The DNA integrity number and concentration are useful parameters for successful comprehensive genomic profiling test for cancer using formalin-fixed paraffin embedded tissue.

The acquisition of high-quality biospecimens and the appropriate handling of these materials are indispensable for successful clinical sequencing. We developed a cancer clinical sequencing system targeting 160 cancer genes: PleSSision-Rapid. Through the PleSSision-Rapid system, we have analyzed DNA quality evaluated by DIN (DNA integrity number) with 1329 formalin-fixed paraffin embedded (FFPE) samples including 477 prospectively collected tissues for genomic test (P) and 852 archival samples after routine pathological diagnosis (A1/A2). As a result, the samples with more than DIN 2.1 was 92.0% (439/477) in prospectively collected sample (P), while it was 85.6% (332/388) and 76.7% (356/464) in two types of archival samples (A1/A2). We performed the PleSSision-Rapid sequence using the samples with over DIN 2.1 and DNA concentration >10 ng/µL with which we were able to construct a DNA library, and the probability of sequence success was almost equivalent during all types of specimen processing, at 90.7% (398/439) in (P), 92.5% (307/332) in (A1) and 90.2% (321/356) in (A2), respectively. Our result indicated the clinical benefit to prepare the prospective collection of FFPE materials for indisputable clinical sequence, and that DIN ≥ 2.1 would be a solid parameter for sample preparation of comprehensive genomic profiling tests.

Emerging RUNX2-Mediated Gene Regulatory Mechanisms Consisting of Multi-Layered Regulatory Networks in Skeletal Development.

Skeletal development is tightly coordinated by chondrocytes and osteoblasts, which are derived from skeletal progenitors, and distinct cell-type gene regulatory programs underlie the specification and differentiation of cells. Runt-related transcription factor 2 (Runx2) is essential to chondrocyte hypertrophy and osteoblast differentiation. Genetic studies have revealed the biological functions of Runx2 and its involvement in skeletal genetic diseases. Meanwhile, molecular biology has provided a framework for our understanding of RUNX2-mediated transactivation at a limited number of cis-regulatory elements. Furthermore, studies using next-generation sequencing (NGS) have provided information on RUNX2-mediated gene regulation at the genome level and novel insights into the multiple layers of gene regulatory mechanisms, including the modes of action of RUNX2, chromatin accessibility, the concept of pioneer factors and phase separation, and three-dimensional chromatin organization. In this review, I summarize the emerging RUNX2-mediated regulatory mechanism from a multi-layer perspective and discuss future perspectives for applications in the treatment of skeletal diseases.

Appropriate use of cancer comprehensive genome profiling assay using circulating tumor DNA.

Comprehensive genomic profiling (CGP) is being increasingly used for the routine clinical management of solid cancers. In July 2018, the use of tumor tissue-based CGP assays became available for all solid cancers under the universal health insurance system in Japan. Several restrictions presently exist, such as patient eligibility and limitations on the opportunities to perform such assays. The clinical implementation of CGP based on plasma circulating tumor DNA (ctDNA) is also expected to raise issues regarding the selection and use of tissue DNA and ctDNA CGP. A Joint Task Force for the Promotion of Cancer Genome Medicine comprised of three Japanese cancer-related societies has formulated a policy proposal for the appropriate use of plasma CGP (in Japanese), available at https://www.jca.gr.jp/researcher/topics/2021/files/20210120.pdf, http://www.jsco.or.jp/jpn/user_data/upload/File/20210120.pdf, and https://www.jsmo.or.jp/file/dl/newsj/2765.pdf. Based on these recommendations, the working group has summarized the respective advantages and cautions regarding the use of tissue DNA CGP and ctDNA CGP with reference to the advice of a multidisciplinary expert panel, the preferred use of plasma specimens over tissue, and multiple ctDNA testing. These recommendations have been prepared to maximize the benefits of performing CGP assays and might be applicable in other countries and regions.

The application of shotgun metagenomics to the diagnosis of granulomatous amoebic encephalitis due to Balamuthia mandrillaris: a case report.

BACKGROUND: Granulomatous amoebic encephalitis (GAE) is an infrequent and fatal infectious disease worldwide. Antemortem diagnosis in this condition is very difficult because clinical manifestations and neuroimaging are nonspecific. CASE PRESENTATION: A 60-year-old Japanese woman was admitted with a chief complaint of left homonymous hemianopsia. Brain-MRI showed extensive necrotizing lesions enhanced by gadolinium, in the right frontal lobe, right occipital lobe, and left parietal lobe. Epithelioid granulomas of unknown etiology were found in the biopsied brain specimens. Shotgun metagenomic sequencing using a next-generation sequencer detected DNA fragments of Balamuthia mandrillaris in the tissue specimens. The diagnosis of granulomatous amoebic encephalitis was confirmed using an amoeba-specific polymerase chain reaction and immunostaining on the biopsied tissues. CONCLUSIONS: Shotgun metagenomics is useful for the diagnosis of central nervous system infections such as GAE wherein the pathogens are difficult to identify.

Whole-genome resequencing using genomic DNA extracted from horsehair roots for gene-doping control in horse sports.

Gene doping is prohibited in horseracing and equestrian sports. In previous studies, we developed non-targeted transgene and genome editing detection methods based on whole genome resequencing (WGR) using genomic DNA extracted from whole blood. In this study, we aimed to develop a WGR method using DNA extracts from hair roots. Hair roots are a preferred substrate because their collection is less invasive than blood collection. Hair is also easier to store for long periods of time. Although almost all genomic DNA extracted from hair root samples stored for years at room temperature was degraded, the quality of genomic DNA from samples stored for years at refrigerated temperatures (4-8°C) was maintained. High-molecular-weight genomic DNA was isolated from hair roots using a magnetic silica beads method of extraction, enabling WGR from horsehair root extracts. Nucleotide sequencing results and numbers of single-nucleotide polymorphisms and insertions/deletions concurred with those previously reported for WGR of DNA extracted from whole blood. Therefore, we consider that storing hair samples at refrigerated temperatures prevents degradation of DNA, allowing the detection of gene doping in these samples based on WGR. It is likely this finding will also have a deterrent effect, as it is now possible to test horses with archived samples even if they or their parents are deceased. To our knowledge, this is the first report employing WGR on horsehair roots stored for a long term.

Performance of Oncomine Fusion Transcript kit for formalin-fixed, paraffin-embedded lung cancer specimens.

Gene fusions play an important role in the carcinogenesis of lung adenocarcinoma. The recent association of four oncogenic driver genes, ALK, ROS1, RET, and NTRK1, as lung tumor predictive biomarkers has increased the need for precision medicine. We used formalin-fixed, paraffin-embedded tissue samples of non-small cell lung cancer from 150 EGFR mutation-negative cases and 10 fusion status-known cases and compared the performance of the Oncomine Dx Fusion Transcript Test (ODxFT) with FISH break-apart for the detection of ALK, RET, and ROS1 fusion genes. RNA was extracted from the paraffin-embedded tissue samples with or without macrodissection under hematoxylin and eosin staining, and the ALK fusion gene was independently determined using these assays. Fusion detection analyses were successfully carried out using ODxFT in 150 cases, with only one invalid case. ALK fusion genes were detected at a frequency of 7.3% (11/150) in the lung cancer specimens. Concordance rate between the ODxFT and ALK-FISH analyses was 99.3% (148/149). Sensitivity and specificity were 91.7% and 99.3%, respectively. All the samples with a known fusion status were accurately matched between the two assays. Our results show a high concordance rate between the ODxFT and ALK-FISH analyses. ODxFT was thus validated as an effective method for detecting clinically significant ALK fusion genes in paraffin-embedded tissue samples.

Urinary tract infection due to anaerobic bacteria in a two-month-old infant.

The significance of anaerobic bacteria as a pathogen in urinary tract infection (UTI) in children is unclear. A two-month-old infant presenting with poor feeding received a diagnosis of polymicrobial anaerobic UTI by next-generation sequencing and was found to have obstructive uropathy. Anaerobic bacteria may be a cause of UTI in children with urinary tract obstruction.

FGFR Signaling as a Candidate Therapeutic Target for Cancers Resistant to Carbon Ion Radiotherapy.

Radiotherapy is an essential component of cancer therapy. Carbon ion radiotherapy (CIRT) promises to improve outcomes compared with standard of care in many cancers. Nevertheless, clinicians often observe in-field recurrence after CIRT. This indicates the presence of a subset of cancers that harbor intrinsic resistance to CIRT. Thus, the development of methods to identify and sensitize CIRT-resistant cancers is needed. To address this issue, we analyzed a unique donor-matched pair of clinical specimens: a treatment-naïve tumor, and the tumor that recurred locally after CIRT in the same patient. Exon sequencing of 409 cancer-related genes identified enrichment of somatic mutations in FGFR3 and FGFR4 in the recurrent tumor compared with the treatment-naïve tumor, indicating a pivotal role for FGFR signaling in cancer cell survival through CIRT. Inhibition of FGFR using the clinically available pan-FGFR inhibitor LY2874455 sensitized multiple cancer cell lines to carbon ions at 3 Gy (RBE: relative biological effectiveness), the daily dose prescribed to the patient. The sensitizer enhancement ratio was 1.66 ± 0.17, 1.27 ± 0.09, and 1.20 ± 0.18 in A549, H1299, and H1703 cells, respectively. Our data indicate the potential usefulness of the analytical pipeline employed in this pilot study to identify targetable mutations associated with resistance to CIRT, and of LY21874455 as a sensitizer for CIRT-resistant cancers. The results warrant validation in larger cohorts.

Deep-Sequence Identification and Role in Virus Replication of a JC Virus Quasispecies in Patients with Progressive Multifocal Leukoencephalopathy.

JC virus (JCV) is a DNA virus causing progressive multifocal leukoencephalopathy (PML) in immunodeficient patients. In the present study, 22 genetic quasispecies with more than 1.5% variant frequency were detected in JCV genomes from six clinical samples of PML by next-generation sequencing. A mutation from A to C at nucleotide (nt) 3495 in JCV Mad1 resulting in a V-to-G amino acid substitution at amino acid (aa) position 392 of the large T antigen (TAg) was identified in all six cases of PML at 3% to 19% variant frequencies. Transfection of JCV Mad1 DNA possessing the V392G substitution in TAg into IMR-32 and human embryonic kidney 293 (HEK293) cells resulted in dramatically decreased production of JCV-encoded proteins. The virus DNA copy number was also reduced in supernatants of the mutant virus-transfected cells. Transfection of the IMR-32 and HEK293 cells with a virus genome containing a revertant mutation recovered viral production and protein expression. Cotransfection with equal amounts of wild-type genome and mutated JCV genome did not reduce the expression of viral proteins or viral replication, suggesting that the mutation did not have any dominant-negative function. Finally, immunohistochemistry demonstrated that TAg was expressed in all six pathological samples in which the quasispecies were detected. In conclusion, the V392G amino acid substitution in TAg identified frequently in PML lesions has a function in suppressing JCV replication, but the frequency of the mutation was restricted and its role in PML lesions was limited. IMPORTANCE: DNA viruses generally have lower mutation frequency than RNA viruses, and the detection of quasispecies in JCV has rarely been reported. In the present study, a next-generation sequencer identified a JCV quasispecies with an amino acid substitution in the T antigen in patients with PML. In vitro studies showed that the mutation strongly repressed the expression of JC viral proteins and reduced the viral replication. However, because the frequency of the mutation was low in each case, the total expression of virus proteins was sustained in vivo. Thus, JC virus replicates in PML lesions in the presence of a mutant virus which is able to repress virus replication.

Bacterial community structures in air conditioners installed in Japanese residential buildings.

The bacterial community structures in four Japanese split-type air conditioners were analyzed using a next-generation sequencer. A variety of bacteria were detected in the air filter of an air conditioner installed on the first floor. In the evaporator of this air conditioner, bacteria belonging to the genus Methylobacterium, or the family of Sphingomonadaceae, were predominantly detected. On the other hand, the majority of bacteria detected in the air filters and evaporators of air conditioners installed on the fifth and twelfth floors belonged to the family Enterobacteriaceae. The source of bacteria belonging to the family Enterobacteriaceae may have been aerosols generated by toilet flushing in the buildings. Our results suggested the possibility that the bacterial contamination in the air conditioners was affected by the floor level on which they were installed. The air conditioner installed on the lower floor, near the ground, may have been contaminated by a variety of outdoor bacteria, whereas the air conditioners installed on floors more distant from the ground may have been less contaminated by outdoor bacteria. However, these suppositions may apply only to the specific split-type air conditioners that we analyzed, because our sample size was small.

A case of pulmonary thromboembolism due to coagulation factor V Leiden in Japan ~ usefulness of next generation sequencing~.

BACKGROUND: Because the venous thromboembolisms (VTEs) due to the coagulation factor V R506Q (FV Leiden) mutation is often seen in Caucasians, the VTE onset in Japan has not been reported. CASE PRESENTATION: A 34-year-old man from north Africa experiencing sudden dyspnea went to a hospital for advice. The patient had pain in his right leg and a high plasma D-dimer level. A contrast-enhanced computed tomography scan revealed a contrast deficit in the bilateral pulmonary artery and in the right lower extremity. The patient was diagnosed with VTE, and anticoagulation therapy was initiated. Our targeted gene panel sequencing revealed that the occurrence of VTE was attributed to a presence of the FV Leiden mutation. CONCLUSIONS: This is the first report demonstrating VTE caused by the FV Leiden mutation in Japan.

SOX2 nonsense mutation in a patient clinically diagnosed with non-syndromic hypogonadotropic hypogonadism.

Hypogonadotropic hypogonadism (HH) is a genetically heterogeneous condition that occurs either as an isolated disorder or as a component of congenital malformation syndromes. SOX2 is a causative gene of syndromic HH characterized by anophthalmia, microphthalmia, or coloboma and other neurological defects such as epilepsy. To date, the causal relationship between SOX2 abnormalities and non-syndromic HH remains speculative. Here, we identified a nonsense mutation of SOX2 in a male patient clinically diagnosed with non-syndromic HH. The patient had epilepsy but no additional clinical features. Ophthalmological examination revealed no abnormalities except for decreased thickness of the retinal nerve fiber layer. Audiometry showed mild sensorineural hearing impairment of both ears. Hormonal evaluation suggested isolated gonadotropin deficiency. Next-generation sequencing-based mutation screening of 13 major causative genes for HH identified a p.Lys35∗ mutation in SOX2 and excluded pathogenic mutations in other tested genes. The p.Lys35∗ mutation appeared to encode a non-functioning SOX2 protein that lacks 283 of 317 amino acids. The SOX2 mutation was absent in the maternal DNA sample, while a paternal sample was unavailable for sequence analysis. These results expand the clinical consequences of SOX2 haploinsufficiency to include non-syndromic HH. Systematic mutation screening using a next-generation sequencer and detailed evaluation of nonspecific ocular/neurological features may help identify SOX2 mutation-positive individuals among HH patients.

Distribution of norovirus genotypes and subtypes in river water by ultra-deep sequencing-based analysis.

To determine the distribution of Norovirus (NoV) genotypes in natural river water in Thailand, we conducted a genome analysis using a next-generation sequencer. Twenty-five river water samples were collected at five different sites of the Khlong Klon River in the suburbs of Bangkok between August 2013 and December 2014. The partial genome of NoV was detected in 15 of the 25 samples (60·0%). Seven of these 15 samples (46·7%) contained multiple NoV GII genotypes: GII.4, GII.6, and GII.17. Our data showed that GII.17 had already emerged in August 2013 as a minor population, and it became a major genotype in December 2014. Our findings indicate that the virus was likely to have been circulating in the community before it appeared in the river water. SIGNIFICANCE AND IMPACT OF THE STUDY: Our study was to investigate the frequencies of multiple genogroups and genotypes of norovirus in the river water near Bangkok, Thailand, by ultra-deep sequencing-based analysis. This study revealed that the epidemic strain was likely to have been circulating in the community before it appeared in the river water. Monitoring of the Norovirus (NoV) genomes in the natural environment may contribute to an understanding of the emergence of new epidemic NoV strains in human populations.

Periodic Pattern of Genetic and Fitness Diversity during Evolution of an Artificial Cell-Like System.

Genetic and phenotypic diversity are the basis of evolution. Despite their importance, however, little is known about how they change over the course of evolution. In this study, we analyzed the dynamics of the adaptive evolution of a simple evolvable artificial cell-like system using single-molecule real-time sequencing technology that reads an entire single artificial genome. We found that the genomic RNA population increases in fitness intermittently, correlating with a periodic pattern of genetic and fitness diversity produced by repeated diversification and domination. In the diversification phase, a genomic RNA population spreads within a genetic space by accumulating mutations until mutants with higher fitness are generated, resulting in an increase in fitness diversity. In the domination phase, the mutants with higher fitness dominate, decreasing both the fitness and genetic diversity. This study reveals the dynamic nature of genetic and fitness diversity during adaptive evolution and demonstrates the utility of a simplified artificial cell-like system to study evolution at an unprecedented resolution.

A robust phylogeny among major lineages of the East African cichlids.

The huge monophyletic group of the East African cichlid radiations (EAR) consists of thousands of species belonging to 12-14 tribes; the number of tribes differs among studies. Many studies have inferred phylogenies of EAR tribes using various genetic markers. However, these phylogenies partly contradict one another and can have weak statistic support. In this study, we conducted maximum-likelihood (ML) phylogenetic analyses using restriction site-associated DNA (RAD) sequences and propose a new robust phylogenetic hypothesis among Lake Tanganyika cichlid fishes, which cover most EAR tribes. Data matrices can vary in size and contents depending on the strategies used to process RAD sequences. Therefore, we prepared 23 data matrices with various processing strategies. The ML phylogenies inferred from 15 large matrices (2.0×10(6) to 1.1×10(7) base pairs) resolved every tribe as a monophyletic group with 100% bootstrap support and shared the same topology regarding relationships among the tribes. Most nodes among the tribes were supported by 100% bootstrap values, and the bootstrap support for the other node varied among the 15 ML trees from 70% to 100%. These robust ML trees differ partly in topology from those in earlier studies, and these phylogenetic relationships have important implications for the tribal classification of EAR.

Impact of different antibiotics on methane production using waste-activated sludge: mechanisms and microbial community dynamics.

Anaerobic digestion is an effective method for reducing the by-product of waste-activated sludge (WAS) from wastewater treatment plants and for producing bioenergy from WAS. However, only a limited number of studies have attempted to improve anaerobic digestion by targeting the microbial interactions in WAS. In this study, we examined whether different antibiotics positively, negatively, or neutrally influence methane fermentation by evaluating changes in the microbial community and functions in WAS. Addition of azithromycin promoted the microbial communities related to the acidogenic and acetogenic stages, and a high concentration of soluble proteins and a high activity of methanogens were detected. Chloramphenicol inhibited methane production but did not affect the bacteria that contribute to the hydrolysis, acidogenesis, and acetogenesis digestion stages. The addition of kanamycin, which exhibits the same methane productivity as a control (antibiotic-free WAS), did not affect all of the microbial communities during anaerobic digestion. This study demonstrates the simultaneous functions and interactions of diverse bacteria and methanogenic Archaea in different stages of the anaerobic digestion of WAS. The ratio of Caldilinea, Methanosarcina, and Clostridium may correspond closely to the trend of methane production in each antibiotic. The changes in microbial activities and function by antibiotics facilitate a better understanding of bioenergy production.