Extracellular Vesicles in Head and Neck Cancer: A Potential New Trend in Diagnosis, Prognosis, and Treatment.

Head and neck cancer (HNC) is a fatal and debilitating disease that is characterized by steady, poor survival rates despite advances in treatment. There is an urgent and unmet need to improve our understanding of what drives this insidious cancer and causes poor outcomes. Extracellular vesicles (EVs) are small vesicles that originate from tumor cells, immune cells, and other cell types and are secreted into plasma, saliva, and other bio-fluids. EVs represent dynamic, real-time changes of cells and offer an exciting opportunity to improve our understanding of HNC biology that may translate to improved clinical practice. Considering the amplified interest in EVs, we have sought to provide a contemporary review of the most recent and salient literature that is shaping the field. Herein, we discuss the functionality of EVs in HNCs and their clinical potential with regards to biomarker and therapeutic capabilities.

Hypoxia Causes Transgenerational Impairment of Ovarian Development and Hatching Success in Fish.

Hypoxia is a pressing environmental problem in both marine and freshwater ecosystems globally, and this problem will be further exacerbated by global warming in the coming decades. Recently, we reported that hypoxia can cause transgenerational impairment of sperm quality and quantity in fish (in F0, F1, and F2 generations) through DNA methylome modifications. Here, we provide evidence that female fish ( Oryzias melastigma) exposed to hypoxia exhibit reproductive impairments (follicle atresia and retarded oocyte development), leading to a drastic reduction in hatching success in the F2 generation of the transgenerational group, although they have never been exposed to hypoxia. Further analyses show that the observed transgenerational impairments in ovarian functions are related to changes in the DNA methylation and expression pattern of two gene clusters that are closely associated with stress-induced cell cycle arrest and cell apoptosis. The observed epigenetic and transgenerational alterations suggest that hypoxia may pose a significant threat to the sustainability of natural fish populations.

Tools and data services registry: a community effort to document bioinformatics resources.

Life sciences are yielding huge data sets that underpin scientific discoveries fundamental to improvement in human health, agriculture and the environment. In support of these discoveries, a plethora of databases and tools are deployed, in technically complex and diverse implementations, across a spectrum of scientific disciplines. The corpus of documentation of these resources is fragmented across the Web, with much redundancy, and has lacked a common standard of information. The outcome is that scientists must often struggle to find, understand, compare and use the best resources for the task at hand.Here we present a community-driven curation effort, supported by ELIXIR-the European infrastructure for biological information-that aspires to a comprehensive and consistent registry of information about bioinformatics resources. The sustainable upkeep of this Tools and Data Services Registry is assured by a curation effort driven by and tailored to local needs, and shared amongst a network of engaged partners.As of November 2015, the registry includes 1785 resources, with depositions from 126 individual registrations including 52 institutional providers and 74 individuals. With community support, the registry can become a standard for dissemination of information about bioinformatics resources: we welcome everyone to join us in this common endeavour. The registry is freely available at https://bio.tools.

Pathogenesis of POLR1C-dependent Type 3 Treacher Collins Syndrome revealed by a zebrafish model.

Treacher Collins Syndrome (TCS) is a rare congenital birth disorder (1 in 50,000 live births) characterized by severe craniofacial defects, including the downward slanting palpebral fissures, hypoplasia of the facial bones, and cleft palate (CP). Over 90% of patients with TCS have a mutation in the TCOF1 gene. However, some patients exhibit mutations in two new causative genes, POLR1C and POLR1D, which encode subunits of RNA polymerases I and III, that affect ribosome biogenesis. In this study, we examine the role of POLR1C in TCS using zebrafish as a model system. Our data confirmed that polr1c is highly expressed in the facial region, and dysfunction of this gene by knockdown or knock-out resulted in mis-expression of neural crest cells during early development that leads to TCS phenotype. Next generation sequencing and bioinformatics analysis of the polr1c mutants further demonstrated the up-regulated p53 pathway and predicted skeletal disorders. Lastly, we partially rescued the TCS facial phenotype in the background of p53 mutants, which supported the hypothesis that POLR1C-dependent type 3 TCS is associated with the p53 pathway.

Transcriptomic and Functional Analyses on the Effects of Dioxin on Insulin Secretion of Pancreatic Islets and ß-Cells.

In this study, transcriptomic and Ingenuity Pathway Analysis (IPA) underlined that an ex-vivo TCDD treatment (0.1 nM) stimulated insulin-release in mouse pancreatic islets via the effect on the Akt-mTOR-p70S6K, AMPK and ERK1/2 pathways. Functional studies using both ex-vivo islets and the mouse ß-cell-line (Min-6) validated the stimulatory effects of TCDD (0.1 and 1 nM) on basal-insulin secretion. At 0.1 nM TCDD treatment on Min-6, Western blot analysis showed activation of ERK1/2 and decreased expression of pyruvate dehydrogenase kinase (PDK). A reduction of PDK expression is associated with an increase of pyruvate dehydrogenase flux. This observation was supported by the detection of significantly higher cellular ATP levels, an increase of glucose-stimulated-insulin-secretion (GSIS), and an inhibition of the AMPK pathway. At 1 nM TCDD treatment on Min-6, significant inhibitions of the Akt-mTOR pathway, cellular ATP production, and GSIS were evident. The experimental studies in Min-6 supported the IPA of transcriptomic data in pancreatic islets. Collectively, TCDD treatment caused an elevated basal-insulin release in both islets and ß-cell cultures. Moreover, our data revealed that the modulation of the Akt-mTOR-p70S6K, AMPK and ERK1/2 pathways might be an important component of the mechanism for the TCDD-perturbing effects on ATP production in ß-cells in affecting insulin secretion.

Effects of in Utero PFOS Exposure on Transcriptome, Lipidome, and Function of Mouse Testis.

Transcriptomic and LC-MS/MS-based targeted lipidomic analyses were conducted to identify the effects of in utero PFOS exposure on neonatal testes and its relation to testicular dysfunction in adult offspring. Pregnant mice were orally administered 0.3 and 3 µg PFOS/g body weight until term. Neonatal testes (P1) were collected for the detection of PFOS, and were subjected to omics study. Integrated pathway analyses using DAVID, KEGG, and IPA underlined the effects of PFOS exposure on lipid metabolism, oxidative stress and cell junction signaling in testes. LC-MS/MS analysis showed that the levels of adrenic acid and docosahexaenoic acid (DHA) in testes were significantly reduced in the PFOS treatment groups. A significant linear decreasing trend in eicosapentaenoic acid and DHA with PFOS concentrations was observed. Moreover, LOX-mediated 5-hydroxyeicosatetraenoic acids (HETE) and 15-HETE from arachidonic acid in the testes were significantly elevated and a linear increasing trend of 15-HETE concentrations was detected with doses of PFOS. The perturbations of lipid mediators suggested that PFOS has potential negative impacts on testicular functions. Postnatal analysis of male offspring at P63 showed significant reductions in serum testosterone and epididymal sperm count. This study sheds light into the as yet unrevealed action of PFOS on lipid mediators in affecting testicular functions.

Transcriptomic analysis reveals specific osmoregulatory adaptive responses in gill mitochondria-rich cells and pavement cells of the Japanese eel.

BACKGROUND: Homeostasis of ions and water is important for the maintenance of cellular functions. The regulation of the homeostasis is particularly important in euryhaline fish that migrate between freshwater (FW) and seawater (SW) environments. The fish gill, the major tissue that forms an interface separating the extracellular fluids and external water environment, has an effective transport system to maintain and regulate a constant body osmolality. In fish gills, the two major epithelial cells, pavement cells (PVCs) and mitochondria-rich cells (MRCs), are known to play key and complementary roles in ion transport at the interface. Discovering the robust mechanisms underlying the two cell types' response to osmotic stress would benefit our understanding of the fundamental mechanism allowing PVCs and MRCs to handle osmotic stress. Owing to the limited genomic data available on estuarine species, existing knowledge in this area is slim. In this study, transcriptome analyses were conducted using PVCs and MRCs isolated from Japanese eels adapted to FW or SW environments to provide a genome-wide molecular study to unravel the fundamental processes at work. RESULTS: The study identified more than 12,000 transcripts in the gill cells. Interestingly, remarkable differential expressed genes (DEGs) were identified in PVCs (970 transcripts) instead of MRCs (400 transcripts) in gills of fish adapted to FW or SW. Since PVCs cover more than 90 % of the gill epithelial surface, the greater change in gene expression patterns in PVCs in response to external osmolality is anticipated. In the integrity pathway analysis, 19 common biological functions were identified in PVCs and MRCs. In the enriched signaling pathways analysis, most pathways differed between PVCs and MRCs; 14 enriched pathways were identified in PVCs and 12 in MRCs. The results suggest that the osmoregulatory responses in PVCs and MRCs are cell-type specific, which supports the complementary functions of the cells in osmoregulation. CONCLUSIONS: This is the first study to provide transcriptomic analysis of PVCs and MRCs in gills of eels adapted to FW or SW environments. It describes the cell-type specific transcriptomic network in different tonicity. The findings consolidate the known osmoregulatory pathways and provide molecular insight in osmoregulation. The presented data will be useful for researchers to select their targets for further studies.

Tissue-specific transcriptome assemblies of the marine medaka Oryzias melastigma and comparative analysis with the freshwater medaka Oryzias latipes.

BACKGROUND: The marine medaka Oryzias melastigma has been demonstrated as a novel model for marine ecotoxicological studies. However, the lack of genome and transcriptome reference has largely restricted the use of O. melastigma in the assessment of in vivo molecular responses to environmental stresses and the analysis of biological toxicity in the marine environment. Although O. melastigma is believed to be phylogenetically closely related to Oryzias latipes, the divergence between these two species is still largely unknown. Using Illumina high-throughput RNA sequencing followed by de novo assembly and comprehensive gene annotation, we provided transcriptomic resources for the brain, liver, ovary and testis of O. melastigma. We also investigated the possible extent of divergence between O. melastigma and O. latipes at the transcriptome level. RESULTS: More than 14,000 transcripts across brain, liver, ovary and testis in marine medaka were annotated, of which 5880 transcripts were orthologous between O. melastigma and O. latipes. Tissue-enriched genes were identified in O. melastigma, and Gene Ontology analysis demonstrated the functional specificity of the annotated genes in respective tissue. Lastly, the identification of marine medaka-enriched transcripts suggested the necessity of generating transcriptome dataset of O. melastigma. CONCLUSIONS: Orthologous transcripts between O. melastigma and O. latipes, tissue-enriched genes and O. melastigma-enriched transcripts were identified. Genome-wide expression studies of marine medaka require an assembled transcriptome, and this sequencing effort has generated a valuable resource of coding DNA for a non-model species. This transcriptome resource will aid future studies assessing in vivo molecular responses to environmental stresses and those analyzing biological toxicity in the marine environment.

Integrative identification of Epstein-Barr virus-associated mutations and epigenetic alterations in gastric cancer.

BACKGROUND & AIMS: The mechanisms by which Epstein-Barr virus (EBV) contributes to the development of gastric cancer are unclear. We investigated EBV-associated genomic and epigenomic variations in gastric cancer cells and tumors. METHODS: We performed whole-genome, transcriptome, and epigenome sequence analyses of a gastric adenocarcinoma cell line (AGS cells), before and after EBV infection. We then looked for alterations in gastric tumor samples, with (n = 34) or without (n = 100) EBV infection, collected from patients at the Prince of Wales Hospital, Chinese University of Hong Kong (from 1998 through 2004), or the First Affiliated Hospital of Sun Yat-sen University, Guangzhou, China (from 1999 through 2006). RESULTS: Transcriptome analysis showed that infected cells expressed 9 EBV genes previously detected in EBV-associated gastric tumors and 71 EBV genes not previously reported in gastric tumors. Ten viral genes that had not been reported previously in gastric cancer but were expressed most highly in EBV-infected cells also were expressed in primary EBV-positive gastric tumors. Whole-genome sequence analysis identified 45 EBV-associated nonsynonymous mutations. These mutations, in genes such as AKT2, CCNA1, MAP3K4, and TGFBR1, were associated significantly with EBV-positive gastric tumors, compared with EBV-negative tumors. An activating mutation in AKT2 was associated with reduced survival times of patients with EBV-positive gastric cancer (P = .006); this mutation was found to dysregulate mitogen-activated protein kinase signaling. Integrated epigenome and transcriptome analyses identified 216 genes transcriptionally down-regulated by EBV-associated hypermethylation; methylation of ACSS1, FAM3B, IHH, and TRABD increased significantly in EBV-positive tumors. Overexpression of Indian hedgehog (IHH) and TraB domain containing (TRABD) increased proliferation and colony formation of gastric cancer cells, whereas knockdown of these genes reduced these activities. We found 5 signaling pathways (axon guidance, focal adhesion formation, interactions among cytokines and receptors, mitogen-activated protein kinase signaling, and actin cytoskeleton regulation) to be affected commonly by EBV-associated genomic and epigenomic alterations. CONCLUSIONS: By using genomic, transcriptome, and epigenomic comparisons of EBV infected vs noninfected gastric cancer cells and tumor samples, we identified alterations in genes, gene expression, and methylation that affect different signaling networks. These might be involved in EBV-associated gastric carcinogenesis.

The NGS WikiBook: a dynamic collaborative online training effort with long-term sustainability.

Next-generation sequencing (NGS) is increasingly being adopted as the backbone of biomedical research. With the commercialization of various affordable desktop sequencers, NGS will be reached by increasing numbers of cellular and molecular biologists, necessitating community consensus on bioinformatics protocols to tackle the exponential increase in quantity of sequence data. The current resources for NGS informatics are extremely fragmented. Finding a centralized synthesis is difficult. A multitude of tools exist for NGS data analysis; however, none of these satisfies all possible uses and needs. This gap in functionality could be filled by integrating different methods in customized pipelines, an approach helped by the open-source nature of many NGS programmes. Drawing from community spirit and with the use of the Wikipedia framework, we have initiated a collaborative NGS resource: The NGS WikiBook. We have collected a sufficient amount of text to incentivize a broader community to contribute to it. Users can search, browse, edit and create new content, so as to facilitate self-learning and feedback to the community. The overall structure and style for this dynamic material is designed for the bench biologists and non-bioinformaticians. The flexibility of online material allows the readers to ignore details in a first read, yet have immediate access to the information they need. Each chapter comes with practical exercises so readers may familiarize themselves with each step. The NGS WikiBook aims to create a collective laboratory book and protocol that explains the key concepts and describes best practices in this fast-evolving field.

Transcriptome sequencing of Chinese and Caucasian population identifies ethnic-associated differential transcript abundance of heterogeneous nuclear ribonucleoprotein K (hnRNPK).

Gene expression variations (GEV) among different ethnic groups have been a subject matter for extensive study. Relatively less known is the extent of alternative splicing variations (ASV) in the context of ethnicity. We conducted a transcriptome sequencing study of 20 lymphoblastoid cell lines obtained from Caucasian and Han Chinese, and identified known genes that exhibit differential isoform abundance between the two ethnic groups. Among them hnRNPK, a co-factor of p53 (TP53), could be further replicated in a 39-sample cohort with TaqMan assay. Although within-population novel splice variants are common, inter-population novel splice variants are rare. We further analyzed 5.63 billion sequencing reads retrieved from the NCBI Sequence Read Archive and identified potential ethnic-specific transcribed regions.

ViralFusionSeq: accurately discover viral integration events and reconstruct fusion transcripts at single-base resolution.

SUMMARY: Insertional mutagenesis from virus infection is an important pathogenic risk for the development of cancer. Despite the advent of high-throughput sequencing, discovery of viral integration sites and expressed viral fusion events are still limited. Here, we present ViralFusionSeq (VFS), which combines soft-clipping information, read-pair analysis and targeted de novo assembly to discover and annotate viral-human fusions. VFS was used in an RNA-Seq experiment, simulated DNA-Seq experiment and re-analysis of published DNA-Seq datasets. Our experiments demonstrated that VFS is both sensitive and highly accurate. AVAILABILITY: VFS is distributed under GPL version 3 at http://hkbic.cuhk.edu.hk/software/viralfusionseq

The SEQanswers wiki: a wiki database of tools for high-throughput sequencing analysis.

Recent advances in sequencing technology have created unprecedented opportunities for biological research. However, the increasing throughput of these technologies has created many challenges for data management and analysis. As the demand for sophisticated analyses increases, the development time of software and algorithms is outpacing the speed of traditional publication. As technologies continue to be developed, methods change rapidly, making publications less relevant for users. The SEQanswers wiki (SEQwiki) is a wiki database that is actively edited and updated by the members of the SEQanswers community (http://SEQanswers.com/). The wiki provides an extensive catalogue of tools, technologies and tutorials for high-throughput sequencing (HTS), including information about HTS service providers. It has been implemented in MediaWiki with the Semantic MediaWiki and Semantic Forms extensions to collect structured data, providing powerful navigation and reporting features. Within 2 years, the community has created pages for over 500 tools, with approximately 400 literature references and 600 web links. This collaborative effort has made SEQwiki the most comprehensive database of HTS tools anywhere on the web. The wiki includes task-focused mini-reviews of commonly used tools, and a growing collection of more than 100 HTS service providers. SEQwiki is available at: http://wiki.SEQanswers.com/.

SEQanswers: an open access community for collaboratively decoding genomes.

SUMMARY: The affordability of high-throughput sequencing has created an unprecedented surge in the use of genomic data in basic, translational and clinical research. The rapid evolution of sequencing technology, coupled with its broad adoption across biology and medicine, necessitates fast, collaborative interdisciplinary discussion. SEQanswers provides a real-time knowledge-sharing resource to address this need, covering experimental and computational aspects of sequencing and sequence analysis. Developers of popular analysis tools are among the >4000 active members, and ~40 peer-reviewed publications have referenced SEQanswers. AVAILABILITY: The SEQanswers community is freely accessible at http://SEQanswers.com/

Detection of splicing events and multiread locations from RNA-seq data based on a geometric-tail (GT) distribution of intron length.

BACKGROUND: RNA sequencing (RNA-seq) measures gene expression levels and permits splicing analysis. Many existing aligners are capable of mapping millions of sequencing reads onto a reference genome. For reads that can be mapped to multiple positions along the reference genome (multireads), these aligners may either randomly assign them to a location, or discard them altogether. Either way could bias downstream analyses. Meanwhile, challenges remain in the alignment of reads spanning across splice junctions. Existing splicing-aware aligners that rely on the read-count method in identifying junction sites are inevitably affected by sequencing depths. RESULTS: The distance between aligned positions of paired-end (PE) reads or two parts of a spliced read is dependent on the experiment protocol and gene structures. We here proposed a new method that employs an empirical geometric-tail (GT) distribution of intron lengths to make a rational choice in multireads selection and splice-sites detection, according to the aligned distances from PE and sliced reads. CONCLUSIONS: GT models that combine sequence similarity from alignment, and together with the probability of length distribution, could accurately determine the location of both multireads and spliced reads.

Restoration of the Oral Microbiota After Surgery for Head and Neck Squamous Cell Carcinoma Is Associated With Patient Outcomes.

OBJECTIVE: To evaluate the dynamics of the oral microbiome and associated patient outcomes following treatment of head and neck squamous cell carcinoma (HNSCC). MATERIALS AND METHODS: This was a prospective cohort study at a tertiary academic center in Hong Kong SAR of patients with head and neck squamous cell carcinoma evaluating the oral microbiome in pre- and postsurgery oral rinses (at 1, 3, and 6 months) with 16S rRNA gene V3-V4 amplicon sequencing. RESULTS: In total, 76 HNSCC patients were evaluated. There was a significantly depressed alpha diversities of oral microbial communities observed in HNSCC oral rinse samples within the first 6 months post-surgery when compared to presurgery or healthy controls. Distant clustering between pre- and postsurgery was also observed (p < 0.022). Following treatment, eight oral bacterial genera showed a trend towards the restoration in the relative abundances that approximate healthy persons. In evaluating patient outcomes, the decreased relative abundance of three periodontal bacteria (Capnocytophaga, Prevotella 7, and Leptotrichia) and the increased relative abundance of two commensal bacteria (Streptococcus and Rothia) at 6 months postsurgery compared to presurgery showed a better 3-year disease-specific survival (a cutoff of Kaplan-Meier survival curve test p < 0.3 at 36 months). In particular, the postsurgery restoration of Prevotella 7 was statistically significant in the surveyed patients (survival rate of 84% vs. 56% at 36 months, p = 0.0065). CONCLUSIONS: Oral microbiome dysbiosis associated with HNSCC is dynamic. These dynamics of the oral microbiome postsurgery are also associated with patient treatment and outcomes and may serve as potential biomarkers for patient management in HNSCC.

Transcriptomic analysis reveals the oncogenic role of S6K1 in hepatocellular carcinoma.

The p70 ribosomal protein S6 kinase 1 (S6K1), a serine/threonine kinase, is commonly overexpressed in a variety of cancers. However, its expression level and functional roles in hepatocellular carcinoma (HCC), which ranks as the third leading cause of cancer-related death worldwide, is still largely unknown. In the current report, we show the in vivo and in vitro overexpression of S6K1 in HCC. In the functional analysis, we demonstrate that S6K1 is required for the proliferation and colony formation abilities in HCC. By using comparative transcriptomic analysis followed by gene ontology enrichment analysis and Ingenuity Pathway Analysis, we find that the depletion of S6K1 can elevate the expression of a cluster of apoptotic genes, tumor suppressor genes and immune responsive genes. Moreover, the knockdown of S6K1 is predicted to reduce the tumorigenicity of HCC through the regulation of hubs of genes including STAT1, HDAC4, CEBPA and ONECUT1. In conclusion, we demonstrate the oncogenic role of S6K1 in HCC, suggesting the possible use of S6K1 as a therapeutic target for HCC treatment.

Replication Study for the Association of GWAS-associated Loci With Adolescent Idiopathic Scoliosis Susceptibility and Curve Progression in a Chinese Population.

STUDY DESIGN: A genetic association (replication) study. OBJECTIVE: The aim of this study was to replicate and further evaluate the association among seven genome-wide association studies (GWAS)-identified single nucleotide polymorphisms (SNPs) in Chinese girls with adolescent idiopathic scoliosis (AIS) with disease onset, curve types, and progression. SUMMARY OF BACKGROUND DATA: AIS is the most common pediatric spinal deformity with a strong genetic predisposition. Recent GWAS identified 10 new disease predisposition loci for AIS. METHODS: Three hundred nineteen female AIS patients with Cobb angle ≥ 10 and 201 healthy controls were studied for the association with disease onset. Seven GWAS-identified SNPs (rs11190870 in LBX1, rs12946942 in SOX9/KCNJ2, rs13398147 in PAX3/EPH4, rs241215 in AJAP1, rs3904778 in BNC2, rs6570507 in GPR126, and rs678741 in LBX1-AS1) were analyzed. In subgroup analysis, AIS patients were subdivided by curve types and disease progression to examine for genotype association. RESULTS: We replicated the association with disease onset in four common SNPs rs11190870, rs3904778, rs6570507, and rs678741. In addition, rs1190870 and rs678741 remained significantly associated in the right thoracic curves only subgroup. However, no significant difference was observed with both clinical curve progression or Cobb angle. CONCLUSION: This study replicated the associations of four GWAS-associated SNPs with occurrence of AIS in our Chinese population. However, none of these SNPs was associated with curve severity and progression. The results suggest that curve progression may be determined by environmental (nongenetic) factor, but further study with a larger sample size is required to address this issue. LEVEL OF EVIDENCE: 4.

A Targeted Gene Panel That Covers Coding, Non-coding and Short Tandem Repeat Regions Improves the Diagnosis of Patients With Neurodegenerative Diseases.

Genetic testing for neurodegenerative diseases (NDs) is highly challenging because of genetic heterogeneity and overlapping manifestations. Targeted-gene panels (TGPs), coupled with next-generation sequencing (NGS), can facilitate the profiling of a large repertoire of ND-related genes. Due to the technical limitations inherent in NGS and TGPs, short tandem repeat (STR) variations are often ignored. However, STR expansions are known to cause such NDs as Huntington's disease and spinocerebellar ataxias type 3 (SCA3). Here, we studied the clinical utility of a custom-made TGP that targets 199 NDs and 311 ND-associated genes on 118 undiagnosed patients. At least one known or likely pathogenic variation was found in 54 patients; 27 patients demonstrated clinical profiles that matched the variants; and 16 patients whose original diagnosis were refined. A high concordance of variant calling were observed when comparing the results from TGP and whole-exome sequencing of four patients. Our in-house STR detection algorithm has reached a specificity of 0.88 and a sensitivity of 0.82 in our SCA3 cohort. This study also uncovered a trove of novel and recurrent variants that may enrich the repertoire of ND-related genetic markers. We propose that a combined comprehensive TGPs-bioinformatics pipeline can improve the clinical diagnosis of NDs.

Transcriptomic analysis reveals transgenerational effect of hypoxia on the neural control of testicular functions.

There are over 400 hypoxic zones in the ocean worldwide. Both laboratory and field studies have shown that hypoxia causes endocrine disruption and reproductive impairments in vertebrates. More importantly, our recent study discovered that parental (F0) hypoxia exposure resulted in the transgenerational impairment of sperm quality in the F2 generation through the epigenetic regulation of germ cells. In the present study, we aim to test the hypothesis that the brain, as the major regulator of the brain-pituitary-gonad (BPG) axis, is also involved in the observed transgenerational effect. Using comparative transcriptomic analysis on brain tissues of marine medaka Oryzias melastigma, 45 common differentially expressed genes caused by parental hypoxia exposure were found in the hypoxic group of the F0 and F2 generations, and the transgenerational groups of the F2 generation. The bioinformatic analysis on this deregulated gene cluster further highlighted the possible involvement of the brain in the transgenerational effect of hypoxia on testicular structure, including abnormal morphologies of the epididymis and the seminal vesicle, and degeneration of the seminiferous tubule. This finding is concordant to the result of hematoxylin and eosin staining, which showed the reduction of testicular lobular diameter in the F0 and F2 generations. Our study demonstrated for the first time the involvement of the brain in the transgenerational effect of hypoxia.