Search

Olive flounder (Paralichthys olivaceus) is a commercially important flatfish species cultured in East Asia. Female flounders generally grow more rapidly than males, therefore control of the sex ratio seems to be a proposed way to increase production. However, the sex determination gene and sex determination mechanism have yet been elucidated. The brain is an important organ that is involved in gonadal development. To explore the sex differences of gene expression in the brain before and during the flounder gonadal differentiation, we used messenger RNA (mRNA)-seq technology to investigate transcriptomes of male and female brains. Between female and male brains, 103 genes were differentially expressed before ovarian differentiation, 16 genes were differentially expressed before testicular differentiation, and 64 genes were differentially expressed during gonadal differentiation. According to annotation and Kyoto Encyclopedia of Genes and Genomes information, the differentially expressed genes (DEGs) were involved in circadian rhythm, circadian rhythm-fly, circadian entrainment, dopaminergic synapse, calcium signaling, glutamatergic synapse, taste transduction, herpes simplex infection, long-term depression, retrograde endocannabinoid signaling, and the synaptic vesicle cycle pathways. MicroRNA (miRNA)-seq was performed during the gonadal differentiation and the target genes of miRNAs were predicted. Integrated analysis of mRNA-seq and miRNA-seq showed that 29 of the 64 DEGs were regulated by the differentially expressed miRNAs during the gonadal differentiation. Our study provides a basis for further studies of brain sex differentiation and the molecular mechanism of sex determination in olive flounder., (© 2019 Wiley Periodicals, Inc.)

The discovery of Lactylation may pave the way for novel approaches to investigating sepsis. This study focused on the prognostic and diagnostic significance of lactylated genes in sepsis. RNA sequencing was performed on blood samples from 20 sepsis patients and 10 healthy individuals at Southwest Medical University in Luzhou, Sichuan, China. Genes associated with sepsis were identified through analysis of RNA sequencing data. Afterward, the genes that were expressed differently were compared with the lactylation genes, resulting in the identification of 55 lactylation genes linked to sepsis. The overlapping genes underwent analysis using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Protein-Protein Network Interactions were used to screen for the core genes. The datasets GSE65682, GSE69528, GSE54514, GSE63042, and GSE95233 were obtained from the GEO database to validate core genes. Survival analysis evaluated the predictive significance of central genes in sepsis, while Receiver Operating Characteristic (ROC) Curve analysis was employed to establish the diagnostic value of genes. Additionally, a meta-analysis was conducted to confirm the precision of RNA-seq data. We obtained five peripheral blood samples, including two from healthy individuals, one from a patient with systemic inflammatory response syndrome (SIRS), and two from patients with sepsis. These samples were used to identify the specific location of core genes using 10×single-cell sequencing. High-throughput sequencing and bioinformatics techniques identified two lactylation-related genes (S100A11 and CCNA2) associated with sepsis. Survival analysis indicated that septic patients with reduced levels of S100A11 had a decreased 28-day survival rate compared to those with elevated levels. Conversely, individuals exhibiting decreased CCNA2 levels demonstrated a greater likelihood of surviving for 28 days than those in the high expression category, indicating a favorable association with survival rates among sepsis patients (P < 0.05). Both genes showed high sensitivity and specificity based on the ROC curve, with AUC values of 0.961 for S100A11 and 0.890 for CCNA2. The meta-analysis revealed that S100A11 exhibited high levels of expression in the sepsis survivors, whereas it displayed low levels of expression in the non-survivors; on the other hand, CCNA2 demonstrated low expression in the sepsis survivors and high expression in the non-survivors (P < 0.05). Single-cell RNA sequencing ultimately showed that monocyte macrophages, T cells, and B cells exhibited high expression levels of the crucial genes associated with sepsis-induced lactylation. In conclusion, the lactylation genes S100A11 and CCNA2 are strongly linked to sepsis and could be valuable markers for diagnosing, predicting outcomes, and providing guidance for sepsis. [ABSTRACT FROM AUTHOR]

CPSF30, a key polyadenylation factor, also serves as an m6A reader, playing a crucial role in determining RNA fate post-transcription. While its homologs mammals are known to be vital for viral replication and immune evasion, the full scope of CPSF30 in plant, particular in viral regulation, remains less explored. Our study demonstrates that CPSF30 significantly facilitates the infection of turnip mosaic virus (TuMV) in Arabidopsis thaliana, as evidenced by infection experiments on the engineered cpsf30 mutant. Among the two isoforms, CPSF30-L, which were characterized with m6A binding activity, emerged as the primary isoform responding to TuMV infection. Analysis of m6A components revealed potential involvement of the m6A machinery in regulating TuMV infection. In contrast, CPSF30-S exhibited distinct subcellular localization, coalescing with P-body markers (AtDCP1 and AtDCP2) in cytoplasmic granules, suggesting divergent regulatory mechanisms between the isoforms. Furthermore, comprehensive mRNA-Seq and miRNA-Seq analysis of Col-0 and cpsf30 mutants revealed global transcriptional reprogramming, highlighting CPSF30’s role in selectively modulating gene expression during TuMV infection. In conclusion, this research underscores CPSF30’s critical role in the TuMV lifecycle and sets the stage for further exploration of its function in plant viral regulation.

Introduction: Status epilepticus (SE) is a seizure lasting more than 5 min that can have lethal consequences or lead to various neurological disorders, including epilepsy. Using a pilocarpine-induced SE model in mice we investigated temporal changes in the hippocampal transcriptome. Methods: We performed mRNA-seq and microRNA-seq analyses at various times after drug treatment. Results: At 1 h after the start of seizures, hippocampal cells upregulated transcription of immediate early genes and genes involved in the IGF-1, ERK/MAPK and RNAPolII/transcription pathways. At 8 h, we observed changes in the expression of genes associated with oxidative stress, overall transcription downregulation, particularly for genes related to mitochondrial structure and function, initiation of a stress response through regulation of ribosome and translation/EIF2 signaling, and upregulation of an inflammatory response. During the middle of the latent period, 36 h, we identified upregulation of membrane components, cholesterol synthesis enzymes, channels, and extracellular matrix (ECM), as well as an increased inflammatory response. At the end of the latent period, 120 h, most changes in expression were in genes involved in ion transport, membrane channels, and synapses. Notably, we also elucidated the involvement of novel pathways, such as cholesterol biosynthesis pathways, iron/BMP/ferroptosis pathways, and circadian rhythms signaling in SE and epileptogenesis. Discussion: These temporal changes in metabolic reactions indicate an immediate response to injury followed by recovery and regeneration. CREB was identified as the main upstream regulator. Overall, our data provide new insights into molecular functions and cellular processes involved at different stages of seizures and offer potential avenues for effective therapeutic strategies. [ABSTRACT FROM AUTHOR]

IntroductionStatus epilepticus (SE) is a seizure lasting more than 5 min that can have lethal consequences or lead to various neurological disorders, including epilepsy. Using a pilocarpine-induced SE model in mice we investigated temporal changes in the hippocampal transcriptome.MethodsWe performed mRNA-seq and microRNA-seq analyses at various times after drug treatment.ResultsAt 1 h after the start of seizures, hippocampal cells upregulated transcription of immediate early genes and genes involved in the IGF-1, ERK/MAPK and RNA-PolII/transcription pathways. At 8 h, we observed changes in the expression of genes associated with oxidative stress, overall transcription downregulation, particularly for genes related to mitochondrial structure and function, initiation of a stress response through regulation of ribosome and translation/EIF2 signaling, and upregulation of an inflammatory response. During the middle of the latent period, 36 h, we identified upregulation of membrane components, cholesterol synthesis enzymes, channels, and extracellular matrix (ECM), as well as an increased inflammatory response. At the end of the latent period, 120 h, most changes in expression were in genes involved in ion transport, membrane channels, and synapses. Notably, we also elucidated the involvement of novel pathways, such as cholesterol biosynthesis pathways, iron/BMP/ferroptosis pathways, and circadian rhythms signaling in SE and epileptogenesis.DiscussionThese temporal changes in metabolic reactions indicate an immediate response to injury followed by recovery and regeneration. CREB was identified as the main upstream regulator. Overall, our data provide new insights into molecular functions and cellular processes involved at different stages of seizures and offer potential avenues for effective therapeutic strategies.

The standard evolutionary theory of aging predicts a negative relationship (trade-off) between fecundity and longevity. However, in principle, the fecundity–longevity relationship can become positive in populations in which individuals have unequal resources. Positive fecundity–longevity relationships also occur in queens of eusocial insects such as ants and bees. Developmental diet is likely to be central to determining trade-offs as it affects key fitness traits, but its exact role remains uncertain. For example, in Drosophila melanogaster , changes in adult diet can affect fecundity, longevity, and gene expression throughout life, but it is unknown how changes in developmental (larval) diet affect fecundity–longevity relationships and gene expression in adults. Using D. melanogaster , we tested the hypothesis that varying developmental diets alters the directionality of fecundity–longevity relationships in adults, and characterized associated gene expression changes. We reared larvae on low (20%), medium (100%), and high (120%) yeast diets, and transferred adult females to a common diet. We measured fecundity and longevity of individual adult females and profiled gene expression changes with age. Adult females raised on different larval diets exhibited fecundity–longevity relationships that varied from significantly positive to significantly negative, despite minimal differences in mean lifetime fertility or longevity. Treatments also differed in age-related gene expression, including for aging-related genes. Hence, the sign of fecundity–longevity relationships in adult insects can be altered and even reversed by changes in larval diet quality. By extension, larval diet differences may represent a key mechanistic factor underpinning positive fecundity–longevity relationships observed in species such as eusocial insects. [ABSTRACT FROM AUTHOR]

Recently, a novel poplar mosaic disease caused by bean common mosaic virus (BCMV) was investigated in Populus alba var. pyramidalis in China. Symptom characteristics, physiological performance of the host, histopathology, genome sequences and vectors, and gene regulation at the transcriptional and posttranscriptional levels were analyzed and RT-qPCR (quantitative reverse transcription PCR) validation of expression was performed in our experiments. In this work, the mechanisms by which the BCMV pathogen impacts physiological performance and the molecular mechanisms of the poplar response to viral infection were reported. The results showed that BCMV infection decreased the chlorophyll content, inhibited the net photosynthesis rate (Pn) and stomatal conductance (Gs), and significantly changed chlorophyll fluorescence parameters in diseased leaves. Transcriptome analysis revealed that the expression of the majority of DEGs (differentially expressed genes) involved in the flavonoid biosynthesis pathway was promoted, but the expression of all or almost all DEGs associated with photosynthesis-antenna proteins and the photosynthesis pathway was inhibited in poplar leaves, suggesting that BCMV infection increased the accumulation of flavonoids but decreased photosynthesis in hosts. Gene set enrichment analysis (GSEA) illustrated that viral infection promoted the expression of genes involved in the defense response or plant-pathogen interaction. MicroRNA-seq analysis illustrated that 10 miRNA families were upregulated while 6 families were downregulated in diseased poplar leaves; moreover, miR156, the largest family with the most miRNA members and target genes, was only differentially upregulated in longperiod disease (LD) poplar leaves. Integrated transcriptome and miRNA-seq analyses revealed 29 and 145 candidate miRNA-target gene pairs; however, only 17 and 76 pairs, accounting for 2.2% and 3.2% of all DEGs, were authentically negatively regulated in short-period disease (SD) and LD leaves, respectively. Interestingly, 4 miR156/SPL (squamosa promoter-binding-like protein) miRNA -target gene pairs were identified in LD leaves: the miR156 molecules were upregulated, but SPL genes were downregulated. In conclusion, BCMV infection significantly changed transcriptional and posttranscriptional gene expression in poplar leaves, inhibited photosynthesis, increased the accumulation of flavonoids, induced systematic mosaic symptoms, and decreased physiological performance in diseased poplar leaves. This study elucidated the fine-tuned regulation of poplar gene expression by BCMV; moreover, the results also suggested that miR156/SPL modules played important roles in the virus response and development of viral systematic symptoms in plant virus disease. [ABSTRACT FROM AUTHOR]

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection. [ABSTRACT FROM AUTHOR]

The development of techniques for measuring gene expression globally has greatly expanded our understanding of gene regulatory mechanisms in depth and scale. We can now quantify every intermediate and transition in the canonical pathway of gene expression-from DNA to mRNA to protein-genome-wide. Employing such measurements in parallel can produce rich datasets, but extracting the most information requires careful experimental design and analysis. Here, we argue for the value of genome-wide studies that measure multiple outputs of gene expression over many timepoints during the course of a natural developmental process. We discuss our findings from a highly parallel gene expression dataset of meiotic differentiation, and those of others, to illustrate how leveraging these features can provide new and surprising insight into fundamental mechanisms of gene regulation.

• Polarization of cells dramatically increases cytokine secretion in bronchial epithelial cells. • CFTR F508del mutation itself is sufficient to significantly alter the epigenome and transcriptome of epithelial cells. • The F508del mutation is closely associated with immunity-related pathways in the epithelial cells in the absence of infection. • These data may serve as epithelial biomarkers for CFTR mRNA therapies. Cystic fibrosis (CF) is characterized by chronic inflammation and excessive cytokines secretion in the lung. Isogenic human CF bronchial epithelial (CFBE41o-) cell lines stably expressing wt-CFTR (WTBE) or F508del mutant (CFBE) are widely used tools in understanding responses to stimuli or drugs and CF pathogenesis in vitro. However, the intrinsic cellular differences in culture are unknown. We performed integrative analyses of these isogenic cells at the protein, mRNA, and chromatin levels in the submerged and air-liquid interface (ALI) conditions to determine cell intrinsic effects of mutant versus complemented CFTR expression. CFBE and WTBE cells displayed different cytokine secretion patterns, including IL-6, IL-8, CXCL1, CXCL10, and CCL5. The ALI culture dramatically increased cytokine secretion in both cells. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) result showed different chromatin landscapes upon polarization and CFBE cells, compared to WTBE cells, exhibited higher genome-wide chromatin accessibility under both culture methods. At the transcriptome level, differentially expressed genes identified by mRNA sequencing between two cell lines were highly concentrated in immunity-related pathways. This multilayered study shows that expression of wild-type CFTR has an epithelial cell intrinsic effect on the cell's epigenome and transcriptome particularly in immunity relevant activities. These data will serve as a resource for the CF community and may serve as epithelial biomarkers for CFTR mRNA therapy. [ABSTRACT FROM AUTHOR]

Vibrio parahaemolyticus, as one of the main pathogens of marine vibriosis, has brought huge losses to aquaculture. However, the interaction mechanism between V. parahaemolyticus and Epinephelus coioides remains unclear. Moreover, there is a lack of comprehensive multi-omics analysis of the immune response of grouper spleen to V. parahaemolyticus. Herein, E. coioides was artificially injected with V. parahaemolyticus, and it was found that the mortality was 16.7% in the early stage of infection, and accompanied by obvious histopathological lesions in the spleen. Furthermore, 1586 differentially expressed genes were screened by mRNA-seq. KEGG analysis showed that genes were significantly enriched in immune-related pathways, Acute-phase immune response, Apoptosis, Complement system and Cytokine-cytokine receptor interaction. As for miRNA-seq analysis, a total of 55 significantly different miRNAs were identified. Further functional annotation analysis indicated that the target genes of differentially expressed miRNAs were enriched in three important pathways (Phosphatidylinositol signaling system, Lysosome and Focal adhesions). Through mRNA-miRNA integrated analysis, 1427 significant miRNA-mRNA pairs were obtained and "p53 signaling pathway", "Intestinal immune network for IgA production" were considered as two crucial pathways. Finally, miR-144-y, miR-497-x, novel-m0459-5p, miR-7133-y, miR-378-y, novel-m0440-5p and novel-m0084-3p may be as key miRNAs to regulate immune signaling pathways via the miRNA-mRNA interaction network. The above results suggest that the mRNA-miRNA integrated analysis not only sheds new light on the molecular mechanisms underlying the interaction between host and V. parahaemolyticus but also provides valuable and new insights into resistance to vibrio infection. [ABSTRACT FROM AUTHOR]

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner. • Cactus is a new pipeline for comprehensive ATAC-Seq and mRNA-Seq data analysis. • The pipeline architecture ensures efficient, consistent, and reproducible analysis. • Preprocesses raw reads, conducts differential analyses, and enrichment analyses. • Demonstrates robust multi-modal and multi-species applicability in a case study. • Offers unique capabilities over other pipelines, enhancing biological insights. [ABSTRACT FROM AUTHOR]

Israel represents the southern limit of the distribution of the banded newt (Ommatotriton vittatus). The life cycle of O. vittatus includes several distinct phases: eggs, aquatic larvae, a terrestrial phase and an aquatic reproductive phase. We investigated differences in gene expression during the life cycle and transition of banded newts between terrestrial and aquatic habitats using mRNA-seq. We identified ∼10 k genes that were differentially expressed (DE) in one of the pairwise comparisons between 3 groups: 1 - terrestrial newts (males and females), 2 - aquatic newts (males and females), 3 - aquatic larvae before metamorphosis. The groups were clearly defined by Principal Components Analysis (PCA). The greatest difference was between aquatic newts (males and females) and aquatic larvae: ∼7.4 k DE genes. Of special interest were the ∼2.4 k genes DE between the aquatic and terrestrial phenotypes. These included prominent candidates with known roles in kidney function (uromodulin homologs were strongly associated with aquatic lifestyle), tissue structure (keratins), and the thyroid hormone signaling modulator DUOXA1. Additional developmental and metabolic pathways overrepresented among the identified DE genes included "epidermis development", "nervous system development", "nucleotide-sugar biosynthesis". Overall, both metamorphosis and environmental transition of banded newts involve extensive transcriptomic remodeling involving developmental, metabolic, and cellular pathways. Understanding the roles of these pathways and individual genes is instrumental for studies of transition between habitats, especially those affected by climate change. Furthermore, the phenotypic flexibility of the newt and the underlying regulation of gene expression can shed light on the evolution of terrestrial vertebrates. [Display omitted] • The banded newt has several distinct phases: eggs, aquatic larvae, a terrestrial phase and an aquatic reproductive phase. • These transitions cause extensive transcriptomic remodeling of developmental, metabolic, and cellular pathways. • Identified DE genes have known roles in kidney function, tissue structure, and hormone signaling. [ABSTRACT FROM AUTHOR]

Excessive long-term consumption of dietary carbohydrates, including glucose, sucrose, or fructose, has been shown to have significant impact on genome-wide gene expression, which likely results from changes in metabolic substrate flux. However, there has been no comprehensive study on the acute effects of individual sugars on the genome-wide gene expression that may reveal the genetic changes altering signaling pathways, subsequent metabolic processes, and ultimately physiological/pathological responses. Considering that gene expressions in response to acute carbohydrate ingestion might be different in nutrient sensitive and insensitive mammals, we conducted comparative studies of genome-wide gene expression by deep mRNA sequencing of the liver in nutrient sensitive C57BL/6J and nutrient insensitive BALB/cJ mice. Furthermore, to determine the temporal responses, we compared livers from mice in the fasted state and following ingestion of standard laboratory mouse chow supplemented with plain drinking water or water containing 20% glucose, sucrose, or fructose. Supplementation with these carbohydrates induced unique extents and temporal changes in gene expressions in a strain specific manner. Fructose and sucrose stimulated gene changes peaked at 3 h postprandial, whereas glucose effects peaked at 12 h and 6 h postprandial in C57BL/6J and BABL/cJ mice, respectively. Network analyses revealed that fructose changed genes were primarily involved in lipid metabolism and were more complex in C57BL/6J than in BALB/cJ mice. These data demonstrate that there are qualitative and antitative differences in the normal physiological responses of the liver between these two strains of mice and C57BL/6J is more sensitive to sugar intake than BALB/cJ. [ABSTRACT FROM AUTHOR]

In mammals, DNA methylation patterns are established by various types of DNA methyltransferases and can be stably passed on during cell division, thus creating a paradigm for epigenetic regulation that can mediate long-lasting changes in gene expression even when the initial triggering signal has disappeared. Although functional deficiency of DNMT3A, one of the methyltransferases, leads to abnormal DNA methylation patterns that result in developmental deficits in mammals, the impacts of its overexpression on tissue gene expression and DNA methylation patterns remain unclear. Here, our previously established hDNMT3A transgenic rat model and mRNA sequencing and bisulphite sequencing PCR were used to analyse the impact of hDNMT3A overexpression on tissue transcriptome and methylome, and whether the impact could be inherited intergenerationally was subsequently investigated. Our results revealed that the overexpression of hDNMT3A could induce notable gene expression variations in rat testis and brain. More importantly, 36.02% and 38.89% of these variations could be intergenerationally inherited to offspring without the transmission of the initial endogenic trigger in the brain and testis, respectively. Furthermore, we found that intergenerationally inherited DNA methylation variations in their promoters and exons could be the underlying mechanism. Compared with inheritable variations that were passively induced by environmental factors, these variations were actively induced by endogenous epigenetic modifiers. This study provided evidence for the epigenetic inheritance of endogenous factors that actively induce gene expression and DNA methylation variations; however, more studies are needed to determine the number of generations that these variations can be stably inherited.

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor. [ABSTRACT FROM AUTHOR]

Streptococcus mutans is a gram-positive bacterium that causes tooth decay. The exopolyssacharides, mostly glucans synthesized by the bacterium are responsible for establishing pathogenic bio-films associated with dental caries disease. The regulatory immune and inflammatory reactions implicated by the synthesized glucans are still not clearly understood. In this study, a water-soluble exopolyssacharide (WSP) was extracted from culture of Str. mutans. The structural properties of WSP, [α-(1 → 3, 1 → 6)-D-glucan] were confirmed using Fourier-transform infrared spectroscopy and 13C-nuclear magnetic resonance spectroscopy. Furthermore, the effects of WSP on the global gene expression of the macrophage-like RAW 264.7 cells were analyzed using mRNA-seq analysis. Using Gene Ontology analysis, we compiled a total of 24,421 genes that were upregulated or downregulated by more than 5.0-fold and 0.3-fold, respectively. Most of the transcripts were grouped under immune response and inflammation-related gene categories. Among the 802 immunity-related genes analyzed, chemokine ligand 7 (Ccl7), interleukin-1β (IL-1β), interleukin-1α (IL-1α) and interleukin-6 (IL-6) were upregulated after WSP exposure. In addition, among a total of 344 genes related to inflammation, Ccl7, IL-1α and IL-6 were upregulated. These results suggest that [α-(1 → 3, 1 → 6)-D-glucan] from Str. mutans produces activates macrophages and may contribute to the immune and inflammatory response to periodontal disease. [ABSTRACT FROM AUTHOR]

In mammals, DNA methylation patterns are established by various types of DNA methyltransferases and can be stably passed on during cell division, thus creating a paradigm for epigenetic regulation that can mediate long-lasting changes in gene expression even when the initial triggering signal has disappeared. Although functional deficiency of DNMT3A, one of the methyltransferases, leads to abnormal DNA methylation patterns that result in developmental deficits in mammals, the impacts of its overexpression on tissue gene expression and DNA methylation patterns remain unclear. Here, our previously established hDNMT3A transgenic rat model and mRNA sequencing and bisulphite sequencing PCR were used to analyse the impact of hDNMT3A overexpression on tissue transcriptome and methylome, and whether the impact could be inherited intergenerationally was subsequently investigated. Our results revealed that the overexpression of hDNMT3A could induce notable gene expression variations in rat testis and brain. More importantly, 36.02% and 38.89% of these variations could be intergenerationally inherited to offspring without the transmission of the initial endogenic trigger in the brain and testis, respectively. Furthermore, we found that intergenerationally inherited DNA methylation variations in their promoters and exons could be the underlying mechanism. Compared with inheritable variations that were passively induced by environmental factors, these variations were actively induced by endogenous epigenetic modifiers. This study provided evidence for the epigenetic inheritance of endogenous factors that actively induce gene expression and DNA methylation variations; however, more studies are needed to determine the number of generations that these variations can be stably inherited. [ABSTRACT FROM AUTHOR]

The shapes of plant organs reflect the evolution of each lineage and have been diversified according to lineage-specific adaptations to environment. Research on the molecular pathways responsible for organ shapes has traditionally been focused mainly on leaves or flowers. Thus, little is known about the pathways controlling fruit shapes, despite their diversity in some plant species. In this study, we analyzed oriental persimmon (Diospyros kaki), which exhibits considerable diversity in fruit shapes among cultivars, to elucidate the underlying molecular mechanism using transcriptomic data and quantitative evaluation. First, to filter the candidate genes associated with persimmon fruit shapes, the whole gene expression patterns obtained using mRNA-Seq analysis from 100 individuals, including a segregated population and various cultivars, were assessed to detect correlations with principal component scores for fruit shapes characterized with elliptic Fourier descriptors. Next, a gene co-expression network analysis with weighted gene co-expression network analysis (WGCNA) package revealed that class 1 KNOX family genes and SEEDSTICK function as integrators along with some phytohormone-related genes, to regulate the fruit shape diversity. On the other hand, the OVATE family genes also contribute to fruit shape diversity, of which pathway would be potentially shared with other plant species. Evolutionary aspects suggest that acquisition of a high lineage-specific and variable expression of class 1 KNOX gene, knotted-like homeobox of Arabidopsis thaliana 1 (KNAT1), in young fruit is important for establishing the persimmon-specific mechanism that determines fruit shape diversity. [ABSTRACT FROM AUTHOR]

Bangpungtongseong-san (BT), an oriental herbal medicine, is used to treat obesity in Korea and East Asia and its antiobesity effects have been examined by several researchers. However, the molecular mechanisms of the antihepatic steatosis effects of BT are unclear. In this study, we examined the effects of BT on obesity, particularly nonalcoholic fatty liver disease, by analyzing metabolic and transcriptional responses using mRNA-sequencing profiles. C57BL/6J mice were fed a high-fat diet (HFD) or HFD + BT (1.5%, w/w, BT) for 12 weeks. Phenotype characteristics were estimated, and the antiobesity mechanism was examined using mRNA sequencing transcriptomic profiles in HFD-induced obese mice. BT treatment ameliorated dyslipidemia and hepatic steatosis in HFD-induced obese mice and reduced body weight gain. The levels of hepatic lipotoxicity markers were significantly decreased, while hepatic antioxidant enzyme activities were augmented by BT compared with in the HFD group. BT attenuated HFD-induced fatty liver through transcriptional changes in the liver. BT treatment downregulated mitochondrial oxidative phosphorylation-related genes in the liver, suggesting improved mitochondrial function. BT treatment also decreased the hepatic fibrosis-related transcriptome. Our findings provide insight into the antiobesity effects of BT, an alternative oriental medicine, for treating obesity-related conditions. Metabolic and transcriptional responses to diet-induced obesity with BT treatment improved liver function. [ABSTRACT FROM AUTHOR]

Background: Feline injection-site sarcoma (FISS), an aggressive iatrogenic subcutaneous malignancy, is challenging to manage clinically and little is known about the molecular basis of its pathogenesis. Tumor transcriptome profiling has proved valuable for gaining insights into the molecular basis of cancers and for identifying new therapeutic targets. Here, we report the first study of the FISS transcriptome and the first cross-species comparison of the FISS transcriptome with those of anatomically similar soft-tissue sarcomas in dogs and humans.Methods: Using high-throughput short-read paired-end sequencing, we comparatively profiled FISS tumors vs. normal tissue samples as well as cultured FISS-derived cell lines vs. skin-derived fibroblasts. We analyzed the mRNA-seq data to compare cancer/normal gene expression level, identify biological processes and molecular pathways that are associated with the pathogenesis of FISS, and identify multimegabase genomic regions of potential somatic copy number alteration (SCNA) in FISS. We additionally conducted cross-species analyses to compare the transcriptome of FISS to those of soft-tissue sarcomas in dogs and humans, at the level of cancer/normal gene expression ratios.Results: We found: (1) substantial differential expression biases in feline orthologs of human oncogenes and tumor suppressor genes suggesting conserved functions in FISS; (2) a genomic region with recurrent SCNA in human sarcomas that is syntenic to a feline genomic region of probable SCNA in FISS; and (3) significant overlap of the pattern of transcriptional alterations in FISS with the patterns of transcriptional alterations in soft-tissue sarcomas in humans and in dogs. We demonstrated that a protein, BarH-like homeobox 1 (BARX1), has increased expression in FISS cells at the protein level. We identified 11 drugs and four target proteins as potential new therapies for FISS, and validated that one of them (GSK-1059615) inhibits growth of FISS-derived cells in vitro.Conclusions: (1) Window-based analysis of mRNA-seq data can uncover SCNAs. (2) The transcriptome of FISS-derived cells is highly consistent with that of FISS tumors. (3) FISS is highly similar to soft-tissue sarcomas in dogs and humans, at the level of gene expression. This work underscores the potential utility of comparative oncology in improving understanding and treatment of FISS. [ABSTRACT FROM AUTHOR]

Abstract: Macroalgae (seaweeds) are essential for the functioning of temperate marine ecosystems, but there is increasing evidence to suggest that their survival is under threat from anthropogenic stressors and disease. Nautella italica R11 is recognized as an aetiological agent of bleaching disease in the red alga, Delisea pulchra. Yet, there is a lack of knowledge surrounding the molecular mechanisms involved in this model host–pathogen interaction. Here we report that mutations in the gene encoding for a LuxR‐type quorum sensing transcriptional regulator, RaiR, render N. italica R11 avirulent, suggesting this gene is important for regulating the expression of virulence phenotypes. Using an RNA sequencing approach, we observed a strong transcriptional response of N. italica R11 towards the presence of D. pulchra. In particular, genes involved in oxidative stress resistance, carbohydrate and central metabolism were upregulated in the presence of the host, suggesting a role for these functions in the opportunistic pathogenicity of N. italica R11. Furthermore, we show that RaiR regulates a subset of genes in N. italica R11, including those involved in metabolism and the expression of phage‐related proteins. The outcome of this research reveals new functions important for virulence of N. italica R11 and contributes to our greater understanding of the complex factors mitigating microbial diseases in macroalgae. [ABSTRACT FROM AUTHOR]

Objectives: To analyze the transcriptome of Spodoptera frugiperda 9 (Sf9) cells infected with AcMNPV or AcMNPV- BmK IT. Results: A comprehensive transcriptome profile for Sf9 cells infected with AcMNPV or AcMNPV- BmK IT is shown. 43127572, 46529744 and 47235310 RNA-Seq profiles permitted the quantification of expression levels for control (C), AcMNPV- BmK IT treatment (ABT) and AcMNPV treatment (AT) groups. There were 997 up-regulated or down-regulated candidate genes with significant different expression level in these treatment groups. Conclusion: These results provide a broad spectrum of candidate genes that are critically involved in the molecular regulation mechanism of Sf9 cells infected with AcMNPV- BmK IT. [ABSTRACT FROM AUTHOR]

Lgr5+ supporting cells (SCs) are enriched hair cell (HC) progenitors in the cochlea. Both in vitro and in vivo studies have shown that HC injury can spontaneously activate Lgr5+ progenitors to regenerate HCs in the neonatal mouse cochlea. Promoting HC regeneration requires the understanding of the mechanism of HC regeneration, and this requires knowledge of the key genes involved in HC injury-induced self-repair responses that promote the proliferation and differentiation of Lgr5+ progenitors. Here, as expected, we found that neomycin-treated Lgr5+ progenitors (NLPs) had significantly greater HC regeneration ability, and greater but not significant proliferation ability compared to untreated Lgr5+ progenitors (ULPs) in response to neomycin exposure. Next, we used RNA-seq analysis to determine the differences in the gene-expression profiles between the transcriptomes of NLPs and ULPs from the neonatal mouse cochlea. We first analyzed the genes that were enriched and differentially expressed in NLPs and ULPs and then analyzed the cell cycle genes, the transcription factors, and the signaling pathway genes that might regulate the proliferation and differentiation of Lgr5+ progenitors. We found 9 cell cycle genes, 88 transcription factors, 8 microRNAs, and 16 cell-signaling pathway genes that were significantly upregulated or downregulated after neomycin injury in NLPs. Lastly, we constructed a protein-protein interaction network to show the interaction and connections of genes that are differentially expressed in NLPs and ULPs. This study has identified the genes that might regulate the proliferation and HC regeneration of Lgr5+ progenitors after neomycin injury, and investigations into the roles and mechanisms of these genes in the cochlea should be performed in the future to identify potential therapeutic targets for HC regeneration. [ABSTRACT FROM AUTHOR]

The ban on bisphenol A (BPA) has led to a rapid increase in the use of BPA analogs, and they are increasingly being detected in the natural environment and biological organisms. Studies have pointed out that BPA analogs can lead to adverse health outcomes. However, their interference with ovarian tissue has not been fully elucidated. In this study, seven- to eight-week-old CD-1 mice were exposed to corn oil containing 300 μg/kg/day bisphenol B (BPB) or bisphenol AF (BPAF) through oral gavage, and ovarian tissues were collected at 14 and 28 days of exposure. Ovarian toxicity was evaluated by the ovarian index, ovarian area, and follicle number. mRNA-seq was used to identify differentially expressed genes (DEGs) and infer the association of DEGs with ovarian diseases. BPB or BPAF exposure induced morphological changes in ovarian tissue in CD-1 mice. In addition, Gene Ontology (GO) analysis revealed disturbances in biological processes (BP) associated with steroid biosynthetic process (GO:0006694) and cellular calcium ion homeostasis (GO:0006874). Subsequently, regulatory networks of BPA analogs (BPB or BPAF)-DEGs-ovarian diseases were constructed. Importantly, the expression levels of DEGs and transcription factors (TFs) associated with ovarian disease were altered. BPB or BPAF exposure causes damage to ovarian morphology through the synergistic effects of multiple biological processes and may be associated with altered mRNA expression profiles as a risk factor for ovarian diseases. [Display omitted] • BPB or BPAF exposure causes morphological changes in the ovary of CD-1 mice. • BPB or BPAF exposure alters the transcriptomic profile of ovarian tissue. • Dysregulation of DEGs and TFs is associated with ovarian disease. [ABSTRACT FROM AUTHOR]

Preparation of an Illumina sequencing library for gene expression analysis (mRNA-Seq) requires microgram amounts of starting total RNA or PCR-based amplification. Here we describe a protocol based on T7 linear RNA amplification that does not introduce significant bias, requires only 10 ng total RNA, and generates a directional, fully representative, whole-transcript mRNA-Seq Illumina library that is highly consistent across over three orders of magnitude of input RNA.

Background: Invasive aspergillosis is started after germination of Aspergillus fumigatus conidia that are inhaled by susceptible individuals. Fungal hyphae can grow in the lung through the epithelial tissue and disseminate hematogenously to invade into other organs. Low fungaemia indicates that fungal elements do not reside in the bloodstream for long. Results: We analyzed whether blood represents a hostile environment to which the physiology of A. fumigatus has to adapt. An in vitro model of A. fumigatus infection was established by incubating mycelium in blood. Our model allowed to discern the changes of the gene expression profile of A. fumigatus at various stages of the infection. The majority of described virulence factors that are connected to pulmonary infections appeared not to be activated during the blood phase. Three active processes were identified that presumably help the fungus to survive the blood environment in an advanced phase of the infection: iron homeostasis, secondary metabolism, and the formation of detoxifying enzymes. Conclusions: We propose that A. fumigatus is hardly able to propagate in blood. After an early stage of sensing the environment, virtually all uptake mechanisms and energy-consuming metabolic pathways are shut-down. The fungus appears to adapt by trans-differentiation into a resting mycelial stage. This might reflect the harsh conditions in blood where A. fumigatus cannot take up sufficient nutrients to establish self-defense mechanisms combined with significant growth. [ABSTRACT FROM AUTHOR]

Wood formation is a complex developmental process governed by genetic and environmental stimuli. Populus and Eucalyptus are fast-growing, high-yielding tree genera that represent ecologically and economically important species suitable for generating significant lignocellulosic biomass., Comparative analysis of the developing xylem and leaf transcriptomes of Populus trichocarpa and Eucalyptus grandis together with phylogenetic analyses identified clusters of homologous genes preferentially expressed during xylem formation in both species., A conserved set of 336 single gene pairs showed highly similar xylem preferential expression patterns, as well as evidence of high functional constraint. Individual members of multi-gene orthologous clusters known to be involved in secondary cell wall biosynthesis also showed conserved xylem expression profiles. However, species-specific expression as well as opposite (xylem versus leaf) expression patterns observed for a subset of genes suggest subtle differences in the transcriptional regulation important for xylem development in each species., Using sequence similarity and gene expression status, we identified functional homologs likely to be involved in xylem developmental and biosynthetic processes in Populus and Eucalyptus. Our study suggests that, while genes involved in secondary cell wall biosynthesis show high levels of gene expression conservation, differential regulation of some xylem development genes may give rise to unique xylem properties. [ABSTRACT FROM AUTHOR]

Analyses of gene expressions in single cells are important for understanding detailed biological phenomena. Here, a highly sensitive and accurate method by sequencing (called “bead-seq”) to obtain a whole gene expression profile for a single cell is proposed. A key feature of the method is to use a complementary DNA (cDNA) library on magnetic beads, which enables adding washing steps to remove residual reagents in a sample preparation process. By adding the washing steps, the next steps can be carried out under the optimal conditions without losing cDNAs. Error sources were carefully evaluated to conclude that the first several steps were the key steps. It is demonstrated that bead-seq is superior to the conventional methods for single-cell gene expression analyses in terms of reproducibility, quantitative accuracy, and biases caused during sample preparation and sequencing processes. [ABSTRACT FROM AUTHOR]

Bats are the only mammals capable of true flight. Critical adaptations for flight include a pair of dramatically elongated hands with broad wing membranes. To study the molecular mechanisms of bat wing evolution, we perform genomewide mRNA sequencing and in situ hybridization for embryonic bat limbs. We identify seven key genes that display unique expression patterns in embryonic bat wings and feet, compared with mouse fore- and hindlimbs. The expression of all 5′HoxD genes (Hoxd9-13) and Tbx3, six known crucial transcription factors for limb and digit development, is extremely high and prolonged in the elongating wing area. The expression of Fam5c, a tumour suppressor, in bat limbs is bat-specific and significantly high in all short digit regions (the thumb and foot digits). These results suggest multiple genetic changes occurred independently during the evolution of bat wings to elongate the hand digits, promote membrane growth and keep other digits short. Our findings also indicate that the evolution of limb morphology depends on the complex integration of multiple gene regulatory networks and biological processes that control digit formation and identity, chondrogenesis, and interdigital regression or retention. [ABSTRACT FROM AUTHOR]

Background Spaceflight environment have been shown to generate reactive oxygen species (ROS) and induce oxidative stress in plants, but little is known about the gene expression of the ROS gene network in plants grown in long-term spaceflight. The molecular response and adaptation to the spaceflight environment of Mizuna plants harvested after 27 days of cultivation onboard the International Space Station (ISS) were measured using genome-wide mRNA expression analysis (mRNA-Seq). Results Total reads of transcripts from the Mizuna grown in the ISS as well as on the ground by mRNA-Seq showed 8,258 and 14,170 transcripts up-regulated and down-regulated, respectively, in the space-grown Mizuna when compared with those from the ground-grown Mizuna. A total of 20 in 32 ROS oxidative marker genes were up-regulated, including high expression of four hallmarks, and preferentially expressed genes associated with ROSscavenging including thioredoxin, glutaredoxin, and alternative oxidase genes. In the transcription factors of the ROS gene network, MEKK1-MKK4-MPK3, OXI1-MKK4- MPK3, and OXI1-MPK3 of MAP cascades, induction of WRKY22 by MEKK1-MKK4- MPK3 cascade, induction of WRKY25 and repression of Zat7 by Zat12 were suggested. RbohD and RbohF genes were up-regulated preferentially in NADPH oxidase genes, which produce ROS. Conclusions This large-scale transcriptome analysis revealed that the spaceflight environment induced oxidative stress and the ROS gene network activation in the space-grown Mizuna. Among transcripts altered in expression by space conditions, some were common genes response to abiotic and biotic stress. Furthermore, certain genes were exclusively up-regulated in Mizuna grown on the ISS. Surprisingly, Mizuna grew in space normally, as well as on the ground, demonstrating that plants can acclimate to long-term exposure in the spaceflight environment by reprogramming the expression of the ROS gene network. [ABSTRACT FROM AUTHOR]

Salinity is a threat to crops in many parts of the world, and together with drought, it is predicted to be a serious constraint to food security. However, understanding the impact of this stressor on plants is a major challenge due to the involvement of numerous genes and regulatory pathways. While transcriptomic analyses of barley ( Hordeum vulgare L.) under salt stress have been reported with microarrays, there are no reports as yet of the use of mRNA-Seq. We demonstrate the utility of mRNA-Seq by analysing cDNA libraries derived from acutely salt-stressed and unstressed leaf material of H. vulgare cv. Hindmarsh. The data yielded >50 million sequence tags which aligned to 26,944 sequences in the Unigene reference database. To gain maximum information, we performed de novo assembly of unaligned reads and discovered >3,800 contigs, termed novel tentative consensus sequences, which are either new, or significant improvements on current databases. Differential gene expression screening found 48 significantly up-regulated and 62 significantly down-regulated transcripts. The work provides comprehensive insights into genome-wide effects of salinity and is a new resource for the study of gene regulation in barley and wheat. Further, the bioinformatics workflow may be applicable to other non-model plants to establish their transcriptomes and identify unique sequences. [ABSTRACT FROM AUTHOR]

The earliest stages of Huntington disease are marked by changes in gene expression that are caused in an indirect and poorly understood manner by polyglutamine expansions in the huntingtin (HIT) protein. To explore the hypothesis that DNA methylation may be altered in cells expressing mutated HIT, we use reduced representation bisulfite sequencing (RRBS) to map sites of DNA methylation in cells carrying either wild-type or mutant HIT. We find that a large fraction of the genes that change in expression in the presence of mutant huntingtin demonstrate significant changes in DNA methylation. Regions with low CpG content, which have previously been shown to undergo methylation changes in response to neuronal activity, are disproportionately affected. On the basis of the sequence of regions that change in methylation. we identify AP-1 and SOX2 as transcriptional regulators associated with DNA methylation changes, and we confirm these hypotheses using genome-wide chromatin immunoprecipitation sequencing (ChIP-Seq). Our findings suggest new mechanisms for the effects of polyglutamine-expanded HIT. These results also raise important questions about the potential effects of changes in DNA methylation on neurogenesis and cognitive decline in patients with Huntington disease. [ABSTRACT FROM AUTHOR]

A laticifer is a cell involved in plant defense against biotic stresses such as herbivores and microorganisms; however, its gene expression is poorly understood. We compared protein accumulation and transcriptomes among laticifers of lignified and unlignified organs of mulberry ( Morus alba), which has a non-articulated, branched type of laticifer. LA-a (equivalent to MLX56) and its homolog LA-b (insecticidal chitinase-like proteins containing two chitin-binding domains) were major proteins in laticifers of unlignified organs, and another protein (LA-c) was a major protein in laticifers of lignified organs. Purification, cDNA cloning, and bioassay of LA-c revealed that LA-c was an acidic class I chitinase having antifungal but not insecticidal activity. Comparative mRNA-Seq analysis using a GS-FLX revealed transcripts of other possible defense-related proteins. Jacalin-like lectin, galacturonase-inhibitor, and pathogenesis-related proteins were also abundant; however, the relative amounts differed among laticifers of lignified and unlignified organs. The results suggest a discontinuous laticifer network in planta and adaptation to different potential enemies among these organs. [ABSTRACT FROM AUTHOR]

Abstract: To compare the two RNA-sequencing protocols, ribo-minus RNA-sequencing (rmRNA-seq) and polyA-selected RNA-sequencing (mRNA-seq), we acquired transcriptomic data–52 and 32 million alignable reads of 35 bases in length–from the mouse cerebrum, respectively. We found that a higher proportion, 44% and 25%, of the uniquely alignable rmRNA-seq reads, is in intergenic and intronic regions, respectively, as compared to 23% and 15% from the mRNA-seq dataset. Further analysis made an additional discovery of transcripts of protein-coding genes (such as Histone, Heg1, and Dux), ncRNAs, snoRNAs, snRNAs, and novel ncRNAs as well as repeat elements in rmRNA-seq dataset. This result suggests that rmRNA-seq method should detect more polyA- or bimorphic transcripts. Finally, through comparative analyses of gene expression profiles among multiple datasets, we demonstrated that different RNA sample preparations may result in significant variations in gene expression profiles. [ABSTRACT FROM AUTHOR]

• TBT has toxic effects on the juvenile development of Japanese medaka. • TBT induced bent spine during the development of juvenile medaka. • Low dosage of TBT (1 ng/g bw/d) has adipogenic effect on medaka juveniles. • The mRNA-Seq results suggested adipogenesis was caused by metabolic processes. The widely used antifoulant tributyltin chloride (TBT) is highly toxic to aquatic organisms. In the present study, four-week-old Japanese medaka (Oryzias latipes) juveniles were orally exposed to TBT at 1 and 10 ng/g bw/d for 1, 2, and 4 weeks, respectively. Half of the tested medaka juveniles showed bone morphology alterations in both 1 and 10 ng/g bw/d TBT 4-week exposure groups. Nile Red (NR) staining showed that the juveniles exposed to 1 ng/g bw/d TBT for 2 and 4 weeks had significantly enlarged adipocyte areas. The mRNA-Seq analysis indicated that 1 ng/g bw/d TBT exposure for 2 weeks affected bone morphology through developmental processes. The GO and KEGG analyses suggested that the adipogenic effect of TBT observed in this study may be induced by metabolic processes, oxidative phosphorylation, and fatty acid degradation and metabolism pathways. Therefore, both morphological observation and mRNA-Seq analysis showed obesogenic effects and developmental toxicity of TBT to juvenile Japanese medaka. [ABSTRACT FROM AUTHOR]