Association of HERPUD1 genetic variant rs2217332 with age-related macular degeneration and polypoidal choroidal vasculopathy in an Indian cohort.

PURPOSE: Age-related macular degeneration (AMD) and polypoidal choroidal vasculopathy (PCV) are sister diseases and have several similar clinical features and still have few genetic differences. The association of HERPUD1 (homocysteine inducible ER protein with ubiquitin like domain 1) gene variant rs2217332 with PCV is known; however, such association with AMD has not been reported in the Indian population. We analyzed the association of rs2217332 with PCV and AMD to identify the preferential association of this variant with these diseases. METHODS: This is a population-based case-control study consisting of 422 patients (129 AMD cases; 101 PCV cases, 192 healthy controls) recruited from the vitreoretinal clinic Sankara Nethralaya. The sample size for the study was calculated using appropriate power calculation methods. Genotype was determined using PCR-based Sanger sequencing. The SPSS V23.0 statistical package tool was used to calculate chi-square and ROC to determine the association of rs2217332 with control, AMD, and PCV. RESULTS: Here, we report for the first time the association of this genetic variant (rs2217332) with AMD and PCV in the Indian population. The case-control study shows a significant association of this SNP with PCV (P value = 0.002); however, this variant is not significantly associated with AMD (P value = 0.602). Comparison between AMD (as control) and PCV (as case) also showed significant association of the SNP with PCV (P value = 0.02). Minor allele A conferred to increase the risk of PCV. CONCLUSIONS: The study concludes that the genetic variant rs2217332 in HERPUD1 gene is highly significantly associated with PCV and not with AMD in Indian populations.

Cxcl9l and Cxcr3.2 regulate recruitment of osteoclast progenitors to bone matrix in a medaka osteoporosis model.

Bone homeostasis requires continuous remodeling of bone matrix to maintain structural integrity. This involves extensive communication between bone-forming osteoblasts and bone-resorbing osteoclasts to orchestrate balanced progenitor cell recruitment and activation. Only a few mediators controlling progenitor activation are known to date and have been targeted for intervention of bone disorders such as osteoporosis. To identify druggable pathways, we generated a medaka (Oryzias latipes) osteoporosis model, where inducible expression of receptor-activator of nuclear factor kappa-Β ligand (Rankl) leads to ectopic formation of osteoclasts and excessive bone resorption, which can be assessed by live imaging. Here we show that upon Rankl induction, osteoblast progenitors up-regulate expression of the chemokine ligand Cxcl9l. Ectopic expression of Cxcl9l recruits mpeg1-positive macrophages to bone matrix and triggers their differentiation into osteoclasts. We also demonstrate that the chemokine receptor Cxcr3.2 is expressed in a distinct subset of macrophages in the aorta-gonad-mesonephros (AGM). Live imaging revealed that upon Rankl induction, Cxcr3.2-positive macrophages get activated, migrate to bone matrix, and differentiate into osteoclasts. Importantly, mutations in cxcr3.2 prevent macrophage recruitment and osteoclast differentiation. Furthermore, Cxcr3.2 inhibition by the chemical antagonists AMG487 and NBI-74330 also reduced osteoclast recruitment and protected bone integrity against osteoporotic insult. Our data identify a mechanism for progenitor recruitment to bone resorption sites and Cxcl9l and Cxcr3.2 as potential druggable regulators of bone homeostasis and osteoporosis.

Cytoplasmic polyadenylation-mediated translational control of maternal mRNAs directs maternal-to-zygotic transition.

In the earliest stages of animal development following fertilization, maternally deposited mRNAs direct biological processes to the point of zygotic genome activation (ZGA). These maternal mRNAs undergo cytoplasmic polyadenylation (CPA), suggesting translational control of their activation. To elucidate the biological role of CPA during embryogenesis, we performed genome-wide polysome profiling at several stages of zebrafish development. Our analysis revealed a correlation between CPA and polysome-association dynamics, demonstrating a coupling of translation to the CPA of maternal mRNAs. Pan-embryonic CPA inhibition disrupted the maternal-to-zygotic transition (MZT), causing a failure of developmental progression beyond the mid-blastula transition and changes in global gene expression that indicated a failure of ZGA and maternal mRNA clearance. Among the genes that were differentially expressed were those encoding chromatin modifiers and key transcription factors involved in ZGA, including nanog, pou5f3 and sox19b, which have distinct CPA dynamics. Our results establish the necessity of CPA for ensuring progression of the MZT. The RNA-seq data generated in this study represent a valuable zebrafish resource for the discovery of novel elements of the early embryonic transcriptome.

Prevalence of primary mutations in Leber hereditary optic neuropathy: A five-year report from a tertiary eye care center in India.

PURPOSE: Genetic testing for primary mutations m.3460G>A, m.11778G>A, and m.14484T>C in ND1, ND4, and ND6 genes of mitochondrial DNA is the recommended assay for Leber hereditary optic neuropathy (LHON; OMIM 535000). This report discusses the outcome of molecular genetic screening for these three primary mutations in suspected LHON cases in India. METHODS: Two hundred and seventy-eight unrelated presumed LHON patients who were seen at the neuro-ophthalmology clinic of a tertiary eye care center from 2014-2018 were analyzed. They were genotyped for the three common variants by polymerase chain reaction-based direct sequencing, and their plasmy status was also determined by restriction enzyme digestion. RESULTS: Eighty two of 278 patients were positive for one of the 3 common mutations with m.11778G>A in ND4 gene more frequently distributed (N=72) in homoplasmic state (N=59/82). The mean onset age of visual loss was 21.1years (SD, 9.8 years; range, 5-58 years) in patients harboring the primary mutation. The most common clinical presentation was bilateral sequential painless vision loss with central and cecocentral scotomas in the visual field due to optic disc atrophy. CONCLUSIONS: The study subjects are a sample of a much larger number of suspected LHON cases tested for primary mutations in India. (N= 278) and 29.4% (82/278) of patients harbour one of the 3 common mutations. Screening the entire mitochondrial genome and the other nuclear genes encoding mitochondrial protein, would probably aid in identifying the other less common mtDNA mutations causing LHON in Indian population.

Transcriptomic analyses of oncogenic hepatocytes reveal common and different molecular pathways of hepatocarcinogenesis in different developmental stages and genders in krasG12V transgenic zebrafish.

Hepatocellular carcinoma (HCC), the most common type of primary liver cancer, is mainly due to genetic changes in hepatocytes. However, molecular expression in hepatocytes during hepatocarcinogenesis has not been characterized. In this study, using an inducible kras transgenic zebrafish models for HCC, transcriptomic profiles of oncogenic hepatocytes from larvae, male and female adult fish following a brief induction of oncogenic kras were investigated. We found that oncogenic hepatocytes from all the three sources possess most of the cancer hallmarks at molecular level, including Sustaining proliferative signaling, Evading growth suppressors, Resisting cell death, Avoiding immune destruction, Inflammation, Reprogramming of energy metabolism, Angiogenesis, and Activating invasion and metastasis, suggesting the malignant transformation at molecular level could occur at the early stage of hepatocarcinogensis and can be captured in hepatocytes. However, each group of oncogenic hepatocytes also had their own characteristics. Larval oncogenic hepatocytes have cancer stem cell features. Female oncogenic hepatocytes showed resemblance to a mild human HCC subtype while male oncogenic hepatocytes resembled a severe HCC subtype, consistent with the observed sex disparity of HCC in both zebrafish and human. Finally, the two adult groups were more similar to each other than to the larval group, indicating an overwhelming effect of development over the gender.

IFT88 mutations identified in individuals with non-syndromic recessive retinal degeneration result in abnormal ciliogenesis.

Whole genome sequencing (WGS) was performed to identify the variants responsible for inherited retinal degeneration (IRD) in a Caucasian family. Segregation analysis of selected rare variants with pathogenic potential identified a set of compound heterozygous changes p.Arg266*:c.796C>T and p.Ala568Thr:c.1702G>A in the intraflagellar transport protein-88 (IFT88) gene segregating with IRD. Expression of IFT88 with the p.Arg266* and p.Ala568Thr mutations in mIMDC3 cells by transient transfection and in HeLa cells by introducing the mutations using CRISPR-cas9 system suggested that both mutations result in the formation of abnormal ciliary structures. The introduction of the IFT88 p.Arg266* variant in the homozygous state in HeLa cells by CRISPR-Cas9 genome-editing revealed that the mutant transcript undergoes nonsense-mediated decay leading to a significant depletion of IFT88 transcript. Additionally, abnormal ciliogenesis was observed in these cells. These observations suggest that the rare and unique combination of IFT88 alleles observed in this study provide insight into the physiological role of IFT88 in humans and the likely mechanism underlying retinal pathology in the pedigree with IRD.

SMN deficiency alters Nrxn2 expression and splicing in zebrafish and mouse models of spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a progressive neurodegenerative disease affecting lower motor neurons. SMA is caused by mutations in the Survival Motor Neuron 1 (SMN1) gene, which result in reduced levels of functional SMN protein. Biochemical studies have linked the ubiquitously expressed SMN protein to the assembly of pre-mRNA processing U snRNPs, raising the possibility that aberrant splicing is a major defect in SMA. Accordingly, several transcripts affected upon SMN deficiency have been reported. A second function for SMN in axonal mRNA transport has also been proposed that may likewise contribute to the SMA phenotype. The underlying etiology of SMA, however, is still not fully understood. Here, we have used a combination of genomics and live Ca(2+) imaging to investigate the consequences of SMN deficiency in a zebrafish model of SMA. In a transcriptome analyses of SMN-deficient zebrafish, we identified neurexin2a (nrxn2a) as strongly down-regulated and displaying changes in alternative splicing patterns. Importantly, the knock-down of two distinct nrxn2a isoforms phenocopies SMN-deficient fish and results in a significant reduction of motor axon excitability. Interestingly, we observed altered expression and splicing of Nrxn2 also in motor neurons from the Smn(-/-);SMN2(+/+) mouse model of SMA, suggesting conservation of nrxn2 regulation by SMN in mammals. We propose that SMN deficiency affects splicing and abundance of nrxn2a. This may explain the pre-synaptic defects at neuromuscular endplates in SMA pathophysiology.

Genome wide analysis reveals Zic3 interaction with distal regulatory elements of stage specific developmental genes in zebrafish.

Zic3 regulates early embryonic patterning in vertebrates. Loss of Zic3 function is known to disrupt gastrulation, left-right patterning, and neurogenesis. However, molecular events downstream of this transcription factor are poorly characterized. Here we use the zebrafish as a model to study the developmental role of Zic3 in vivo, by applying a combination of two powerful genomics approaches--ChIP-seq and microarray. Besides confirming direct regulation of previously implicated Zic3 targets of the Nodal and canonical Wnt pathways, analysis of gastrula stage embryos uncovered a number of novel candidate target genes, among which were members of the non-canonical Wnt pathway and the neural pre-pattern genes. A similar analysis in zic3-expressing cells obtained by FACS at segmentation stage revealed a dramatic shift in Zic3 binding site locations and identified an entirely distinct set of target genes associated with later developmental functions such as neural development. We demonstrate cis-regulation of several of these target genes by Zic3 using in vivo enhancer assay. Analysis of Zic3 binding sites revealed a distribution biased towards distal intergenic regions, indicative of a long distance regulatory mechanism; some of these binding sites are highly conserved during evolution and act as functional enhancers. This demonstrated that Zic3 regulation of developmental genes is achieved predominantly through long distance regulatory mechanism and revealed that developmental transitions could be accompanied by dramatic changes in regulatory landscape.

Zebrafish Rnf111 is encoded by multiple transcripts and is required for epiboly progression and prechordal plate development.

Arkadia (also known as RING finger 111) encodes a nuclear E3 ubiquitin ligase that targets intracellular effectors and modulators of TGFß/Nodal-related signaling for polyubiquitination and proteasome-dependent degradation. In the mouse, loss of Arkadia results in early embryonic lethality, with defects attributed to compromised Nodal signaling. Here, we report the isolation of zebrafish arkadia/rnf111, which is represented by 5 transcript variants. arkadia/rnf111 is broadly expressed during the blastula and gastrula stages, with eventual enrichment in the anterior mesendoderm, including the prechordal plate. Morpholino knockdown experiments reveal an unexpected role for Arkadia/Rnf111 in both early blastula organization and epiboly progression. Using a splice junction morpholino, we present additional evidence that arkadia/rnf111 transcript variants containing a 3' alternative exon are specifically required for epiboly progression in the late gastrula. This result suggests that arkadia/rnf111 transcript variants encode functionally relevant protein isoforms that provide additional intracellular flexibility and regulation to the Nodal signaling pathway.

Genome-wide mapping of Hif-1α binding sites in zebrafish.

BACKGROUND: Hypoxia Inducible Factor (HIF) regulates a cascade of transcriptional events in response to decreased oxygenation, acting from the cellular to the physiological level. This response is evolutionarily conserved, allowing the use of zebrafish (Danio rerio) as a model for studying the hypoxic response. Activation of the hypoxic response can be achieved in zebrafish by homozygous null mutation of the von Hippel-Lindau (vhl) tumour suppressor gene. Previous work from our lab has focused on the phenotypic characterisation of this mutant, establishing the links between vhl mutation, the hypoxic response and cancer. To further develop fish as a model for studying hypoxic signalling, we examine the transcriptional profile of the vhl mutant with respect to Hif-1α. As our approach uses embryos consisting of many cell types, it has the potential to uncover additional HIF regulated genes that have escaped detection in analogous mammalian cell culture studies. RESULTS: We performed high-density oligonucleotide microarray analysis of the gene expression changes in von Hippel-Lindau mutant zebrafish, which identified up-regulation of well-known hypoxia response genes and down-regulation of genes primarily involved in lipid processing. To identify the dependency of these transcriptional changes on HIF, we undertook Chromatin Immunoprecipitation linked next generation sequencing (ChIP-seq) for the transcription factor Hypoxia Inducible Factor 1α (HIF-1α). We identified HIF-1α binding sites across the genome, with binding sites showing enrichment for an RCGTG motif, showing conservation with the mammalian hypoxia response element. CONCLUSIONS: Transcriptome analysis of vhl mutant embryos detected activation of key hypoxia response genes seen in human cell models of hypoxia, but also suppression of many genes primarily involved in lipid processing. ChIP-seq analysis of Hif-1α binding sites unveiled an unprecedented number of loci, with a high proportion containing a canonical hypoxia response element. Whether these sites are functional remains unknown, nevertheless their frequent location near transcriptional start sites suggests functionality, and will allow for investigation into the potential hypoxic regulation of genes in their vicinity. We expect that our data will be an excellent starting point for analysis of both fish and mammalian gene regulation by HIF.

Deep sequencing of small RNA facilitates tissue and sex associated microRNA discovery in zebrafish.

BACKGROUND: The role of microRNAs in gene regulation has been well established. The extent of miRNA regulation also increases with increasing genome complexity. Though the number of genes appear to be equal between human and zebrafish, substantially less microRNAs have been discovered in zebrafish compared to human (miRBase Release 19). It appears that most of the miRNAs in zebrafish are yet to be discovered. RESULTS: We sequenced small RNAs from brain, gut, liver, ovary, testis, eye, heart and embryo of zebrafish. In brain, gut and liver sequencing was done sex specifically. Majority of the sequenced reads (16-62 %) mapped to known miRNAs, with the exception of ovary (5.7 %) and testis (7.8 %). Using the miRNA discovery tool (miRDeep2), we discovered novel miRNAs from the unannotated reads that ranged from 7.6 to 23.0 %, with exceptions of ovary (51.4 %) and testis (55.2 %). The prediction tool identified a total of 459 novel pre-miRNAs. We compared expression of miRNAs between different tissues and between males and females to identify tissue associated and sex associated miRNAs respectively. These miRNAs could serve as putative biomarkers for these tissues. The brain and liver had highest number of tissue associated (22) and sex associated (34) miRNAs, respectively. CONCLUSIONS: This study comprehensively identifies tissue and sex associated miRNAs in zebrafish. Further, we have discovered 459 novel pre-miRNAs (~30 % seed homology to human miRNA) as a genomic resource which can facilitate further investigations to understand miRNA-mRNA gene regulatory networks in zebrafish which will have implications in understanding the function of human homologs.

Dorsal activity of maternal squint is mediated by a non-coding function of the RNA.

Despite extensive study, the earliest steps of vertebrate axis formation are only beginning to be elucidated. We previously showed that asymmetric localization of maternal transcripts of the conserved zebrafish TGFß factor Squint (Sqt) in 4-cell stage embryos predicts dorsal, preceding nuclear accumulation of ß-catenin. Cell ablations and antisense oligonucleotides that deplete Sqt lead to dorsal deficiencies, suggesting that localized maternal sqt functions in dorsal specification. However, based upon analysis of sqt and Nodal signaling mutants, the function and mechanism of maternal sqt was debated. Here, we show that sqt RNA may function independently of Sqt protein in dorsal specification. sqt insertion mutants express localized maternal sqt RNA. Overexpression of mutant/non-coding sqt RNA and, particularly, the sqt 3'UTR, leads to ectopic nuclear ß-catenin accumulation and expands dorsal gene expression. Dorsal activity of sqt RNA requires Wnt/ß-catenin but not Oep-dependent Nodal signaling. Unexpectedly, sqt ATG morpholinos block both sqt RNA localization and translation and abolish nuclear ß-catenin, providing a mechanism for the loss of dorsal identity in sqt morphants and placing maternal sqt RNA upstream of ß-catenin. The loss of early dorsal gene expression can be rescued by the sqt 3'UTR. Our findings identify new non-coding functions for the Nodal genes and support a model wherein sqt RNA acts as a scaffold to bind and deliver/sequester maternal factors to future embryonic dorsal.

Impaired development of neural-crest cell-derived organs and intellectual disability caused by MED13L haploinsufficiency.

MED13L is a component subunit of the Mediator complex, an important regulator of transcription that is highly conserved across eukaryotes. Here, we report MED13L disruption in a translocation t(12;19) breakpoint of a patient with Pierre-Robin syndrome, moderate intellectual disability, craniofacial anomalies, and muscular defects. The phenotype is similar to previously described patients with MED13L haploinsufficiency. Knockdown of MED13L orthologue in zebrafish, med13b, showed early defective migration of cranial neural crest cells (NCCs) that contributed to cartilage structure deformities in the later stage, recapitulating craniofacial anomalies seen in human patients. Notably, we observed abnormal distribution of developing neurons in different brain regions of med13b morphant embryos, which could be rescued upon introduction of full-length human MED13L mRNA. To compare with mammalian system, we suppressed MED13L expression by short-hairpin RNA in ES-derived human neural progenitors, and differentiated them into neurons. Transcriptome analysis revealed differential expression of components of Wnt and FGF signaling pathways in MED13L-deficient neurons. Our finding provides a novel insight into the mechanism of overlapping phenotypic outcome targeting NCCs derivatives organs in patients with MED13L haploinsufficiency, and emphasizes a clinically recognizable syndromic phenotype in these patients.

Functional mapping of the zebrafish early embryo proteome and transcriptome.

Zebrafish is a popular system for studying vertebrate development and disease that shows high genetic conservation with humans. Molecular level studies at different stages of development are essential for understanding the processes deployed during ontogeny. Here, we performed comparative analysis of the whole proteome and transcriptome of the early stage (24 h post-fertilization) zebrafish embryo. We identified 8363 proteins with their approximate cellular abundances (the largest number of zebrafish embryo proteins quantified thus far), through a combination of thorough deyolking and extensive fractionation procedures, before resolving the peptides by mass spectrometry. We performed deep sequencing of the transcripts and found that the expressed proteome and transcriptome displayed a moderate correlation for the majority of cellular processes. Integrative functional mapping of the quantified genes demonstrated that embryonic developmental systems differentially exploit transcriptional and post-transcriptional regulatory mechanisms to modulate protein abundance. Using network mapping of the low-abundance proteins, we identified various signal transduction pathways important in embryonic development and also revealed genes that may be regulated at the post-transcriptional level. Our data set represents a deep coverage of the functional proteome and transcriptome of the developing zebrafish, and our findings unveil molecular regulatory mechanisms that underlie embryonic development.

Molecular insights into 4-nitrophenol-induced hepatotoxicity in zebrafish: transcriptomic, histological and targeted gene expression analyses.

BACKGROUND: 4-Nitrophenol (4-NP) is a prioritized environmental pollutant and its toxicity has been investigated using zebrafish, advocated as an alternative toxicological model. However, molecular information of 4-NP induced hepatotoxicity is still limited. This study aimed to obtain molecular insights into 4-NP-induced hepatotoxicity using zebrafish as a model. METHODS: Adult male zebrafish were exposed to 4-NP for 8, 24, 48 and 96h. Livers were sampled for microarray experiment, qRT-PCR and various histological analyses. RESULTS: Transcriptomic analysis revealed that genes associated with oxidative phosphorylation and electron transport chain were significantly up-regulated throughout early and late stages of 4-NP exposure due to oxidative phosphorylation uncoupling by 4-NP. This in turn induced oxidative stress damage and up-regulated pathways associated with tumor suppressors Rb and p53, cell cycle, DNA damage, proteasome degradation and apoptosis. Pathways associated with cell adhesion and morphology were deregulated. Carbohydrate and lipid metabolisms were down-regulated while methionine and aromatic amino acid metabolisms as well as NFKB pathway associated with chronic liver conditions were up-regulated. Up-regulation of NFKB, NFAT and interleukin pathways suggested hepatitis. Histological analyses with specific staining methods and qRT-PCR analysis of selected genes corroborated with the transcriptomic analysis suggesting 4-NP induced liver injury. CONCLUSION: Our findings allowed us to propose a plausible model and provide a broader understanding of the molecular events leading to 4-NP induced acute hepatotoxicity for future studies involving other nitrophenol derivatives. GENERAL SIGNIFICANCE: This is the first transcriptomic report on 4-NP induced hepatotoxicity.

Differential transcriptomic analyses revealed genes and signaling pathways involved in iono-osmoregulation and cellular remodeling in the gills of euryhaline Mozambique tilapia, Oreochromis mossambicus.

BACKGROUND: The Mozambique tilapia Oreochromis mossambicus has the ability to adapt to a broad range of environmental salinities and has long been used for investigating iono-osmoregulation. However, to date most studies have focused mainly on several key molecules or parameters hence yielding a limited perspective of the versatile iono-osmoregulation in the euryhaline fish. This study aimed to capture transcriptome-wide differences between the freshwater- and seawater-acclimated gills of the Mozambique tilapia. RESULTS: We have identified over 5000 annotated gene transcripts with high homology (E-value <1.0E-50) to human genes that were differentially expressed in freshwater- and seawater-acclimated gills of the Mozambique tilapia. These putative human homologs were found to be significantly associated with over 50 canonical signaling pathways that are operating in at least 23 biological processes in relation to branchial iono-osmoregulation and cellular remodeling. The analysis revealed multiple signaling pathways in freshwater-acclimated gills acting in concert to maintain cellular homeostasis under hypo-osmotic environment while seawater-acclimated gills abounded with molecular signals to cope with the higher cellular turn-over rate, energetics and iono-regulatory demands under hyper-osmostic stress. Additionally, over 100 transcripts encoding putative inorganic ion transporters/channels were identified, of which several are well established in gill iono-regulation while the remainder are lesser known. We have also validated the expression profiles of 47 representative genes in freshwater- and seawater-acclimated gills, as well as in hypersaline-acclimated (two-fold salinity of seawater) gills. The findings confirmed that many of these responsive genes retained their expression profiles in hypersaline-acclimated gills as in seawater-acclimated gills, although several genes had changed significantly in their expression level/direction in hypersaline-acclimated gills. CONCLUSIONS: This is the first study that has provided an unprecedented transcriptomic-wide perspective of gill iono-osmoregulation since such studies were initiated more than 80 years ago. It has expanded our molecular perspective from a relatively few well-studied molecules to a plethora of gene transcripts and a myriad of canonical signaling pathways driving various biological processes that are operating in gills under hypo-osmotic and hyper-osmotic stresses. These findings would provide insights and resources to fuel future studies on gill iono-osmoregulation and cellular remodeling in response to salinity challenge and acclimation.

Transcriptomic analysis of a transgenic zebrafish hepatocellular carcinoma model reveals a prominent role of immune responses in tumour progression and regression.

Using our previously established xmrk transgenic zebrafish, hepatocellular carcinoma (HCC) was generated by induced expression of xmrk, which encoded a hyperactive epidermal growth factor receptor (EGFR) homolog, and regressed by suppression of xmrk expression. To investigate molecular changes in liver tumour progression and regression, RNA-Seq was performed for induced HCC and early and late stages of liver tissues during tumour regression. We found that Xmrk-induced zebrafish HCC shared strong molecular characteristics with a human HCC subtype (S2), which shows activated Myc signalling, upregulated phosphor-S6 and epithelial cell adhesion molecule. In the HCC stage, there were enhanced proteasome, antigen processing and presentation, aminosugars metabolisms, p53 and cell cycle pathways. During tumour regression, the transcriptomic profile showed a reversed trend of molecular changes compared with human HCC progression. Interestingly, distinct immune responses in tumour progression and regression were observed, including increased major histocompatibility complex class I (MHCI) at the HCC stage, enriched immune cell trafficking signals and inflammation in early regression and enhanced MHCII in late regression. Both neutrophils and macrophages were enriched during tumour progression and regression; however, the distribution of neutrophils and macrophages in HCC was relatively uniform, whereas both types of immune cells were regionally clustered during tumour regression, especially with dominant blood vessel association of macrophage in late regression, suggesting differential functions of these immune cells in tumour progression and regression. As tumour regression in our model resembles the targeted inhibition of EGFR in cancer therapy, our observations may provide molecular insights into the targeted inhibition and highlight the importance of immune response in tumour regression.

Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition.

Maternally deposited mRNAs direct early development before the initiation of zygotic transcription during mid-blastula transition (MBT). To study mechanisms regulating this developmental event in zebrafish, we applied mRNA deep sequencing technology and generated comprehensive information and valuable resources on transcriptome dynamics during early embryonic (egg to early gastrulation) stages. Genome-wide transcriptome analysis documented at least 8000 maternal genes and identified the earliest cohort of zygotic transcripts. We determined expression levels of maternal and zygotic transcripts with the highest resolution possible using mRNA-seq and clustered them based on their expression pattern. We unravel delayed polyadenylation in a large cohort of maternal transcripts prior to the MBT for the first time in zebrafish. Blocking polyadenylation of these transcripts confirms their role in regulating development from the MBT onward. Our study also identified a large number of novel transcribed regions in annotated and unannotated regions of the genome, which will facilitate reannotation of the zebrafish genome. We also identified splice variants with an estimated frequency of 50%-60%. Taken together, our data constitute a useful genomic information and valuable transcriptome resource for gene discovery and for understanding the mechanisms of early embryogenesis in zebrafish.

Liver X receptor agonist T0901317 induced liver perturbation in zebrafish: histological, gene set enrichment and expression analyses.

BACKGROUND: Liver X receptor (LXR), a ligand-activated transcription factor, regulates important biological processes. It has been associated with pathology and proposed as a therapeutic target. The zebrafish is a new vertebrate model for disease modeling, drug and toxicity screening and will be interesting to test for its potential for LXR-related studies. METHODS: Adult male fish were exposed to LXR agonist T0901317 at 20, 200 and 2000nM for 96h and the livers were sampled for histological, microarray and qRT-PCR analyses. RESULTS: Histological analysis suggests dose-dependent perturbation of carbohydrate and lipid metabolisms by T0901317 in the liver, which lead to hepatocyte swelling and cell death. Microarray data revealed several conserved effects of T0901317 with mammalian models, including up-regulation of LXR-targeted genes, modulation of biological pathways associated with proteasome, cell death, extracellular matrix and adhesions, maturity onset diabetes of the young and lipid beta oxidation. Interestingly, this study identified the complement and coagulation systems as down-regulated by T0901317 for the first time, potentially via transcriptional repression by LXR activation. qRT-PCR validated the expression of 16 representative genes, confirming activation of LXR signaling and down-regulation of these biological pathways by T0901317 which could be linked to the anti-thrombogenic, anti-atherogenic and anti-inflammatory actions, as well as metabolic disruptions via LXR activation. CONCLUSION AND GENERAL SIGNIFICANCE: Our study underscores the potential of using zebrafish model coupled with transcriptomic analysis to capture pharmacological and toxicological or pathological events induced by LXR modulators.

Differential transcript isoform usage pre- and post-zygotic genome activation in zebrafish.

BACKGROUND: Zebrafish embryos are transcriptionally silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here, we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA). RESULTS: We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3' untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3' terminal exon when transcribed by the zygote compared to their maternal transcript counterparts. ChIP-sequencing data mapped onto skipped exon events reveal a correlation between histone H3K36 trimethylation peaks and skipped exons, suggesting epigenetic marks being part of alternative splicing regulation. CONCLUSIONS: The novel isoforms and isoform switches reported here include regulators of transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.