15000 unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis

Gene Expression and Phenotypic Assessment of Egg Quality across Developmental Stages of Atlantic Cod throughout the Spawning Season

Egg quality in fishes is commonly determined by fertilisation success and cleavage patterns as a phenotypic outcome of underlying regulatory mechanisms. Although these phenotypic estimators of egg quality are useful in farming conditions, these "good quality" egg batches do not always translate to good larval growth and survival. The identification of genes involved in embryonic development may help find links between genetic factors of maternal origin and egg quality. Herein, the relative expression of seven stage-specific developmental genes of Atlantic cod was analysed using quantitative PCR to understand the function during embryogenesis and its relationship with egg quality. Genes ccnb2 and pvalb1 showed significant differential expression between developmental stages and significant upregulation from blastula and somite stages, respectively. The comparison of spawning batches showed that the relative gene expression of genes ccnb2, acta, tnnt3 and pvalb1 was significantly higher from the middle of the spawning season where phenotypic quality estimators establish the best egg quality. Moreover, a positive significant correlation was observed between quality estimators based on egg morphology and the genetic expression of genes acta and acta1 during somitogenesis. This study suggests that the combination of quality estimators, genetics and batch timing could help optimise reproductive protocols for commercial stocks of Atlantic cod.

Stability and function of RCL1 are dependent on the interaction with BMS1

During ribosome biogenesis, the small subunit (SSU) processome is responsible for 40S assembly. The BMS1/RCL1 complex is a core component of the SSU processome that plays an important role in 18S rRNA processing and maturation. Genetic studies using zebrafish mutants indicate that both Bms1-like (Bms1l) and Rcl1 are essential for digestive organ development. In spite of vital functions of this complex, the mutual dependence of these two nucleolar proteins for the stability and function remains elusive. In this study, we identified an RCL1-interacting domain in BMS1, which is conserved in zebrafish and humans. Moreover, both the protein stability and nucleolar entry of RCL1 depend on its interaction with BMS1, otherwise RCL1 degraded through the ubiquitination-proteasome pathway. Functional studies revealed that overexpression of RCL1 in BMS1-knockdown cells can partially rescue the defects in 18S rRNA processing and cell proliferation, and hepatocyte-specific overexpression of Rcl1 can resume zebrafish liver development in the bms1l substitution mutant bms1lsq163/sq163but not in the knockout mutant bms1lzju1/zju1, which is attributed to the nucleolar entry of Rcl1 in the former mutant. Our data demonstrate that BMS1 and RCL1 interaction is essential for not only pre-rRNA processing but also the communication between ribosome biogenesis and cell cycle regulation.

Hhex and Prox1a synergistically dictate the hepatoblast to hepatocyte differentiation in zebrafish

The specification of endoderm cells to prospective hepatoblasts is the starting point for hepatogenesis. However, how a prospective hepatoblast gains the hepatic fate remains elusive. Previous studies have shown that loss-of-function of either hhex or prox1a alone causes a small liver phenotype but without abolishing the hepatocyte differentiation, suggesting that absence of either Hhex or Prox1a alone is not sufficient to block the hepatoblast differentiation. Here, via genetic studies of the zebrafish two single (hhex-/- and prox1a-/-) and one double (hhex-/-prox1a-/-) mutants, we show that simultaneous loss-of-function of the hhex and prox1a two genes does not block the endoderm cells to gain the hepatoblast potency but abolishes the hepatic differentiation from the prospective hepatoblast. Consequently, the hhex-/-prox1a-/- double mutant displays a liverless phenotype that cannot be rescued by the injection of bmp2a mRNA. Taken together, we provide strong evidences showing that Hhex teams with Prox1a to act as a master control of the differentiation of the prospective hepatoblasts towards hepatocytes.

Cdx1b protects intestinal cell fate by repressing signaling networks for liver specification

In mammals, the expression of the homeobox family member Cdx2/CDX2 is restricted within the intestine. Conditional ablation of the mouse Cdx2 in the endodermal cells causes a homeotic transformation of the intestine towards the esophagus or gastric fate. In this report, we show that null mutants of zebrafish cdx1b, encoding the counterpart of mammalian CDX2, could survive more than 10 days post fertilization, a stage when the zebrafish digestive system has been well developed. Through RNA sequencing (RNA-seq) and single-cell sequencing (scRNA-seq) of the dissected intestine from the mutant embryos, we demonstrate that the loss-of-function of the zebrafish cdx1b yields hepatocyte-like intestinal cells, a phenotype never observed in the mouse model. Further RNA-seq data analysis, and genetic double mutants and signaling inhibitor studies reveal that Cdx1b functions to guard the intestinal fate by repressing, directly or indirectly, a range of transcriptional factors and signaling pathways for liver specification. Finally, we demonstrate that heat shock-induced overexpression of cdx1b in a transgenic fish abolishes the liver formation. Therefore, we demonstrate that Cdx1b is a key repressor of hepatic fate during the intestine specification in zebrafish.

Mycn regulates intestinal development through ribosomal biogenesis in a zebrafish model of Feingold syndrome 1

Feingold syndrome type 1, caused by loss-of-function of MYCN, is characterized by varied phenotypes including esophageal and duodenal atresia. However, no adequate model exists for studying the syndrome's pathological or molecular mechanisms, nor is there a treatment strategy. Here, we developed a zebrafish Feingold syndrome type 1 model with nonfunctional mycn, which had severe intestinal atresia. Single-cell RNA-seq identified a subcluster of intestinal cells that were highly sensitive to Mycn, and impaired cell proliferation decreased the overall number of intestinal cells in the mycn mutant fish. Bulk RNA-seq and metabolomic analysis showed that expression of ribosomal genes was down-regulated and that amino acid metabolism was abnormal. Northern blot and ribosomal profiling analysis showed abnormal rRNA processing and decreases in free 40S, 60S, and 80S ribosome particles, which led to impaired translation in the mutant. Besides, both Ribo-seq and western blot analysis showed that mTOR pathway was impaired in mycn mutant, and blocking mTOR pathway by rapamycin treatment can mimic the intestinal defect, and both L-leucine and Rheb, which can elevate translation via activating TOR pathway, could rescue the intestinal phenotype of mycn mutant. In summary, by this zebrafish Feingold syndrome type 1 model, we found that disturbance of ribosomal biogenesis and blockage of protein synthesis during development are primary causes of the intestinal defect in Feingold syndrome type 1. Importantly, our work suggests that leucine supplementation may be a feasible and easy treatment option for this disease.

Unraveling Differential Transcriptomes and Cell Types in Zebrafish Larvae Intestine and Liver

The zebrafish intestine and liver, as in other vertebrates, are derived from the endoderm. Great effort has been devoted to deciphering the molecular mechanisms controlling the specification and development of the zebrafish intestine and liver; however, genome-wide comparison of the transcriptomes between these two organs at the larval stage remains unexplored. There is a lack of extensive identification of feature genes marking specific cell types in the zebrafish intestine and liver at 5 days post-fertilization, when the larval fish starts food intake. In this report, through RNA sequencing and single-cell RNA sequencing of intestines and livers separately dissected from wild-type zebrafish larvae at 5 days post-fertilization, together with the experimental validation of 47 genes through RNA whole-mount in situ hybridization, we identified not only distinctive transcriptomes for the larval intestine and liver, but also a considerable number of feature genes for marking the intestinal bulb, mid-intestine and hindgut, and for marking hepatocytes and cholangiocytes. Meanwhile, we identified 135 intestine- and 97 liver-enriched transcription factor genes in zebrafish larvae at 5 days post-fertilization. Our findings provide rich molecular and cellular resources for studying cell patterning and specification during the early development of the zebrafish intestine and liver.

Calcium Signaling in the Cerebellar Radial Glia and Its Association with Morphological Changes during Zebrafish Development

Radial glial cells are a distinct non-neuronal cell type that, during development, span the entire width of the brain walls of the ventricular system. They play a central role in the origin and placement of neurons, since their processes form structural scaffolds that guide and facilitate neuronal migration. Furthermore, glutamatergic signaling in the radial glia of the adult cerebellum (i.e., Bergmann glia), is crucial for precise motor coordination. Radial glial cells exhibit spontaneous calcium activity and functional coupling spread calcium waves. However, the origin of calcium activity in relation to the ontogeny of cerebellar radial glia has not been widely explored, and many questions remain unanswered regarding the role of radial glia in brain development in health and disease. In this study we used a combination of whole mount immunofluorescence and calcium imaging in transgenic (gfap-GCaMP6s) zebrafish to determine how development of calcium activity is related to morphological changes of the cerebellum. We found that the morphological changes in cerebellar radial glia are quite dynamic; the cells are remarkably larger and more elaborate in their soma size, process length and numbers after 7 days post fertilization. Spontaneous calcium events were scarce during the first 3 days of development and calcium waves appeared on day 5, which is associated with the onset of more complex morphologies of radial glia. Blockage of gap junction coupling inhibited the propagation of calcium waves, but not basal local calcium activity. This work establishes crucial clues in radial glia organization, morphology and calcium signaling during development and provides insight into its role in complex behavioral paradigms.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Rcl1 depletion impairs 18S pre-rRNA processing at the A1-site and up-regulates a cohort of ribosome biogenesis genes in zebrafish

Yeast Rcl1 is a potential endonuclease that mediates pre-RNA cleavage at the A₂-site to separate 18S rRNA from 5.8S and 25S rRNAs. However, the biological function of Rcl1 in opisthokonta is poorly defined. Moreover, there is no information regarding the exact positions of 18S pre-rRNA processing in zebrafish. Here, we report that zebrafish pre-rRNA harbours three major cleavage sites in the 5′ETS, namely –477nt (A′-site), –97nt (A₀-site) and the 5′ETS and 18S rRNA link (A₁-site), as well as two major cleavage regions within the ITS1, namely 208–218nt (site 2) and 20–33nt (site E). We also demonstrate that depletion of zebrafish Rcl1 mainly impairs cleavage at the A₁-site. Phenotypically, rcl1^–/– mutants exhibit a small liver and exocrine pancreas and die before 15 days post-fertilization. RNA-seq analysis revealed that the most significant event in rcl1^–/– mutants is the up-regulated expression of a cohort of genes related to ribosome biogenesis and tRNA production. Our data demonstrate that Rcl1 is essential for 18S rRNA maturation at the A₁-site and for digestive organogenesis in zebrafish. Rcl1 deficiency, similar to deficiencies in other ribosome biogenesis factors, might trigger a common mechanism to upregulate the expression of genes responsible for ribosome biogenesis.

Preferential transcription of the mutated allele in NPM1 mutated acute myeloid leukaemia

Nucleophosmin is commonly both over-expressed and mutated in acute myeloid leukemia (AML). NPM1 mutations are always heterozygous. In addition, NPM1 has a number of different splice variants with the major variant encoded by exons 1–9 and 11–12 (NPM1.1). Further variants include NPM1.2 which lacks exons 8 and 10 and NPM1.3 which comprises exons 1–10 (and so lacks the region of sequence mutated in AML). In this study we quantified the expression of these three variants in 108 AML patient samples with and without NPM1 mutations and also assessed the level of expression from the wild-type and mutant alleles in variants NPM1.1 and NPM1.2. The results show that NPM1.1 is the most commonly expressed variant, however transcripts from wild-type and mutated alleles do not occur at equal levels, with a significant bias toward the mutated allele. Considering the involvement of mutant nucleophosmin in the progression and maintenance of AML, a bias towards mutated transcripts could have a significant impact on disease maintenance.

Differential impacts of carp and salmon pituitary extracts on induced oogenesis, egg quality, molecular ontogeny and embryonic developmental competence in European eel

Low egg quality and embryonic survival are critical challenges in aquaculture, where assisted reproduction procedures and other factors may impact egg quality. This includes European eel (Anguilla anguilla), where pituitary extract from carp (CPE) or salmon (SPE) is applied to override a dopaminergic inhibition of the neuroendocrine system, preventing gonadotropin secretion and gonadal development. The present study used either CPE or SPE to induce vitellogenesis in female European eel and compared impacts on egg quality and offspring developmental competence with emphasis on the maternal-to-zygotic transition (MZT). Females treated with SPE produced significantly higher proportions of floating eggs with fewer cleavage abnormalities and higher embryonic survival. These findings related successful embryogenesis to higher abundance of mRNA transcripts of genes involved in cell adhesion, activation of MZT, and immune response (dcbld1, epcam, oct4, igm) throughout embryonic development. The abundance of mRNA transcripts of cldnd, foxr1, cea, ccna1, ccnb1, ccnb2, zar1, oct4, and npm2 was relatively stable during the first eight hours, followed by a drop during MZT and low levels thereafter, indicating transfer and subsequent clearance of maternal mRNA. mRNA abundance of zar1, epcam, and dicer1 was associated with cleavage abnormalities, while mRNA abundance of zar1, sox2, foxr1, cldnd, phb2, neurod4, and neurog1 (before MZT) was associated with subsequent embryonic survival. In a second pattern, low initial mRNA abundance with an increase during MZT and higher levels persisting thereafter indicating the activation of zygotic transcription. mRNA abundance of ccna1, npm2, oct4, neurod4, and neurog1 during later embryonic development was associated with hatch success. A deviating pattern was observed for dcbld1, which mRNA levels followed the maternal-effect gene pattern but only for embryos from SPE treated females. Together, the differences in offspring production and performance reported in this study show that PE composition impacts egg quality and embryogenesis and in particular, the transition from initial maternal transcripts to zygotic transcription.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

rDNA subtypes and their transcriptional expression in zebrafish at different developmental stages

Eukaryotic 18S, 5.8S and 28S rRNAs are processed from a single transcript transcribed from the 45S rDNA gene, which is normally tandemly arrayed over hundred copies in a genome. Recently, a maternal (M) subtype and a somatic (S) subtype of rDNA were identified in zebrafish, with M-subtype on chromosome 4 and S-subtype on chromosome 5. It appears that the M-subtype is only expressed in eggs whilst the expression of the S-subtype is coupled with the initiation of zygotic gene expression. In this report, we identified three novel but transcriptionally inactive 18S variants in zebrafish genome with chromosome location different from the M- and S-subtype, suggesting translocation of 18S rDNA fragment during zebrafish evolution. Furthermore, we confirmed that the unfertilized eggs only have the M-subtype transcripts while brain, heart and liver have only the S-subtype transcripts. Both the M- and S-subtype transcripts were detected in female gonad. Our results support that the expression of different subtypes of rDNA is differentially regulated to meet the requirement for 'specialized ribosomes' during different developmental stages.

Capn3 depletion causes Chk1 and Wee1 accumulation and disrupts synchronization of cell cycle reentry during liver regeneration after partial hepatectomy

Recovery of liver mass to a healthy liver donor by compensatory regeneration after partial hepatectomy (PH) is a prerequisite for liver transplantation. Synchronized cell cycle reentry of the existing hepatocytes after PH is seemingly a hallmark of liver compensatory regeneration. Although the molecular control of the PH-triggered cell cycle reentry has been extensively studied, little is known about how the synchronization is achieved after PH. The nucleolus-localized protein cleavage complex formed by the nucleolar protein Digestive-organ expansion factor (Def) and cysteine proteinase Calpain 3 (Capn3) has been implicated to control wounding healing during liver regeneration through selectively cleaving the tumor suppressor p53 in the nucleolus. However, whether the Def-Capn3 complex participates in regulating the synchronization of cell cycle reentry after PH is unknown. In this report, we generated a zebrafish capn3b null mutant (capn3b^∆19∆14). The homozygous mutant was viable and fertile, but suffered from a delayed liver regeneration after PH. Delayed liver regeneration in capn3b^∆19∆14 was due to disruption of synchronized cell proliferation after PH. Mass spectrometry (MS) analysis of nuclear proteins revealed that a number of negative regulators of cell cycle are accumulated in the capn3b^∆19∆14 liver after PH. Moreover, we demonstrated that Check-point kinase 1 (Chk1) and Wee1, two key negative regulators of G2 to M transition, are substrates of Capn3. We also demonstrated that Chk1 and Wee1 were abnormally accumulated in the nucleoli of amputated capn3b^∆19∆14 liver. In conclusion, our findings suggest that the nucleolar-localized Def-Capn3 complex acts as a novel regulatory pathway for the synchronization of cell cycle reentry, at least partially, through inactivating Chk1 and Wee1 during liver regeneration after PH.

Retinoic Acid-Regulated Target Genes During Development: Integrative Genomics Analysis

Retinoic acid (RA), a major natural active metabolite of vitamin A (VA) is well known to play critical roles in embryonic development. The effects of RA are mediated by nuclear receptors (RARs), which regulate the expression of gene batteries involved in cell growth and differentiation. Since the early 1990s several laboratories have focused on understanding how RA-regulated genes and RAR binding sites operate by studying the differentiation of embryonal carcinoma cells and embryonic stem cells. The development of hybridization-based microarray technology and high performance software analysis programs has allowed the characterization of thousands of RA-regulated genes. During the two last decades, publication of the genome sequence of various organisms has allowed advances in massive parallel sequencing and bioinformatics analysis of genome-wide data sets. These new generation sequencing (NGS) technologies have revolutionized the field by providing a global integrated picture of RA-regulated gene networks and the regulatory programs involved in cell fate decisions during embryonal carcinoma and embryonic stem cells differentiation. Now the challenge is to reconstruct the RA-regulated gene networks at the single cell level during the development of specialized embryonic tissues.

Gene redundancy and gene compensation: An updated view

Gene knockdown approaches using antisense oligo nucleotides or analogs such as siRNAs and morpholinos have been widely adopted to study gene functions although the off-target issue has been always a concern in these studies. On the other hand, classic genetic analysis relies on the availability of loss-of-function or gain-of-function mutants. The fast development of genome editing technologies such as TALEN and CRISPR/Cas9 has greatly facilitated the generation of null mutants for the functional studies of target genes in a variety of organisms such as zebrafish. Surprisingly, an unexpected discrepancy was observed between morphant phenotype and mutant phenotype for many genes in zebrafish, i.e., while the morphant often displays an obvious phenotype, the corresponding null mutant appears relatively normal or only exhibits a mild phenotype due to gene compensation. Two recent reports have partially answered this intriguing question by showing that a pre-mature termination codon and homologous sequence are required to elicit the gene compensation and the histone modifying complex COMPASS is involved in activating the expression of the compensatory genes. Here, I summarize these exciting new progress and try to redefine the concept of genetic compensation and gene compensation.

iSeq 2.0: A Modular and Interchangeable Toolkit for Interaction Screening in Yeast

We developed a flexible toolkit for combinatorial screening in Saccharomyces cerevisiae, which generates large libraries of cells, each uniquely barcoded to mark a combination of DNA elements. This interaction sequencing platform (iSeq 2.0) includes genomic landing pads that assemble combinations through sequential integration of plasmids or yeast mating, 15 barcoded plasmid libraries containing split selectable markers (URA3AI, KanMXAI, HphMXAI, and NatMXAI), and an array of ∼24,000 "double-barcoder" strains that can make existing yeast libraries iSeq compatible. Various DNA elements are compatible with iSeq: DNA introduced on integrating plasmids, engineered genomic modifications, or entire genetic backgrounds. DNA element libraries are modular and interchangeable, and any two libraries can be combined, making iSeq capable of performing many new combinatorial screens by short-read sequencing.

Analysis of public RNA-sequencing data reveals biological consequences of genetic heterogeneity in cell line populations

Meta-analysis of datasets available in public repositories are used to gather and summarise experiments performed across laboratories, as well as to explore consistency of scientific findings. As data quality and biological equivalency across samples may obscure such analyses and consequently their conclusions, we investigated the comparability of 85 public RNA-seq cell line datasets. Thousands of pairwise comparisons of single nucleotide variants in 139 samples revealed variable genetic heterogeneity of the eight cell line populations analysed as well as variable data quality. The H9 and HCT116 cell lines were found to be remarkably stable across laboratories (with median concordances of 99.2% and 98.5%, respectively), in contrast to the highly variable HeLa cells (89.3%). We show that the genetic heterogeneity encountered greatly affects gene expression between same-cell comparisons, highlighting the importance of interrogating the biological equivalency of samples when comparing experimental datasets. Both the number of differentially expressed genes and the expression levels negatively correlate with the genetic heterogeneity. Finally, we demonstrate how comparing genetically heterogeneous datasets affect gene expression analyses and that high dissimilarity between same-cell datasets alters the expression of more than 300 cancer-related genes, which are often the focus of studies using cell lines.

Allelic imbalance of tissue-type plasminogen activator (t-PA) gene expression in human brain tissue

We have identified a single-nucleotide polymorphism (SNP) in the t-PA enhancer (-7351C>T), which is associated with endothelial t-PA release in vivo. In vitro studies demonstrated that this SNP is functional at the level of transcription. In the brain, t-PA has been implicated in both physiologic and pathophysiologic processes. The aim of the present study was to examine the effect of the t-PA –7351C>T SNP on t-PA gene expression in human brain tissue. Allelic mRNA expression was measured in heterozygous post-mortem brain tissues using quantitative TaqMan genotyping assay. Protein-DNA interactions were assessed using electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP). Significantly higher levels of t-PA mRNA were generated from chromosomes that harboured the wild-type –7351C allele, as compared to those generated from the mutant T allele (for the hippocampus, C to T allelic ratio of ~1.3, p=0.010, n=12; and for the cortex, C to T allelic ratio of ~1.2, p=0.017, n=12). EMSA showed reduced neuronal and astrocytic nuclear protein binding affinity to the T allele, and identified Sp1 and Sp3 as the major transcription factors that bound to the –7351 site. ChIP analyses confirmed that Sp1 recognises this site in intact cells. In conclusion, the t-PA –7351C>T SNP affects t-PA gene expression in human brain tissue. This finding might have clinical implications for neurological conditions associated with enhanced t-PA levels, such as in the acute phase of cerebral ischaemia, and also for stroke recovery.

Heritable gene expression differences between lake and stream stickleback include both parallel and antiparallel components

The repeated phenotypic patterns that characterize populations undergoing parallel evolution provide support for a deterministic role of adaptation by natural selection. Determining the level of parallelism also at the genetic level is thus central to our understanding of how natural selection works. Many studies have looked for repeated genomic patterns in natural populations, but work on gene expression is less common. The studies that have examined gene expression have found some support for parallelism, but those studies almost always used samples collected from the wild that potentially confounds the effects of plasticity with heritable differences. Here we use two independent pairs of lake and stream threespine stickleback (Gasterosteus aculeatus) raised in common garden conditions to assess both parallel and antiparallel (that is, similar versus different directions of lake–stream expression divergence in the two watersheds) heritable gene expression differences as measured by total RNA sequencing. We find that more genes than expected by chance show either parallel (22 genes, 0.18% of expressed genes) or antiparallel (24 genes, 0.20% of expressed genes) lake–stream expression differences. These results correspond well with previous genomic studies in stickleback ecotype pairs that found similar levels of parallelism. We suggest that parallelism might be similarly constrained at the genomic and transcriptomic levels.

The IL-8 Gene Polymorphisms in Behçet's Disease Observed in Denizli Province of Turkey

Behçet's disease is a multisystemic inflammatory disorder as a triad of symptoms including recurrent oral and genital aphthous ulceration and uveitis with unknown pathogenesis. IL-8, a proinflammatory cytokine, has been found increased in the active stage of BD. DNA samples were obtained from 88 patients with BD and 112 healthy control subjects in Denizli province of Turkey. All genotyping experiments of SNPs in IL-8 gene were performed using polymerase chain reaction-restriction fragment polymorphism. We found that IL-8 -845 T > C and -738 T > A sites are non-polymorphic. There were no differences in the polymorphisms of IL-8 +396 G/T, +781 C/T, and +1633 C/T sites except IL-8 -251 T > A in between patients and healthy controls. Analysis of IL-8 polymorphisms indicates that the distribution of frequencies seems to be associated with -251 T > A and gender, -251 T > A and erythema nodosum, -251 T > A and ocular involvement, +781 C > T and erythema nodosum, +396 G > T and pathergy positivity, and +1633 C > T and papulopustular lesion. We demonstrated that the frequencies of IL-8 haplotypes were significantly different with BD patients than control group. We found that the distribution of IL-8 haplotypes was significantly different with genital ulcers, ocular involvement, positive pathergy test, erythema nodosum, papulopustular lesions, and arthritis with BD patients than healthy control individuals. Our study suggests that IL-8 gene polymorphisms may affect susceptibility to BD and increase the risk of developing disease. In order to confirm and assess the association of IL-8 and other cytokine gene polymorphisms in the pathophysiology of BD, large cohort studies are needed.

A uniform survey of allele-specific binding and expression over 1000-Genomes-Project individuals

Large-scale sequencing in the 1000 Genomes Project has revealed multitudes of single nucleotide variants (SNVs). Here, we provide insights into the functional effect of these variants using allele-specific behaviour. This can be assessed for an individual by mapping ChIP-seq and RNA-seq reads to a personal genome, and then measuring 'allelic imbalances' between the numbers of reads mapped to the paternal and maternal chromosomes. We annotate variants associated with allele-specific binding and expression in 382 individuals by uniformly processing 1,263 functional genomics data sets, developing approaches to reduce the heterogeneity between data sets due to overdispersion and mapping bias. Since many allelic variants are rare, aggregation across multiple individuals is necessary to identify broadly applicable 'allelic elements'. We also found SNVs for which we can anticipate allelic imbalance from the disruption of a binding motif. Our results serve as an allele-specific annotation for the 1000 Genomes variant catalogue and are distributed as an online resource (alleledb.gersteinlab.org).

Differentially methylated obligatory epialleles modulate context-dependent LAM gene expression in the honeybee Apis mellifera

Differential intragenic methylation in social insects has been hailed as a prime mover of environmentally driven organismal plasticity and even as evidence for genomic imprinting. However, very little experimental work has been done to test these ideas and to prove the validity of such claims. Here we analyze in detail differentially methylated obligatory epialleles of a conserved gene encoding lysosomal α-mannosidase (AmLAM) in the honeybee. We combined genotyping of progenies derived from colonies founded by single drone inseminated queens, ultra-deep allele-specific bisulfite DNA sequencing, and gene expression to reveal how sequence variants, DNA methylation, and transcription interrelate. We show that both methylated and non-methylated states of AmLAM follow Mendelian inheritance patterns and are strongly influenced by polymorphic changes in DNA. Increased methylation of a given allele correlates with higher levels of context-dependent AmLAM expression and appears to affect the transcription of an antisense long noncoding RNA. No evidence of allelic imbalance or imprinting involved in this process has been found. Our data suggest that by generating alternate methylation states that affect gene expression, sequence variants provide organisms with a high level of epigenetic flexibility that can be used to select appropriate responses in various contexts. This study represents the first effort to integrate DNA sequence variants, gene expression, and methylation in a social insect to advance our understanding of their relationships in the context of causality.

ASEQ: fast allele-specific studies from next-generation sequencing data

Background

Single base level information from next-generation sequencing (NGS) allows for the quantitative assessment of biological phenomena such as mosaicism or allele-specific features in healthy and diseased cells. Such studies often present with computationally challenging burdens that hinder genome-wide investigations across large datasets that are now becoming available through the 1,000 Genomes Project and The Cancer Genome Atlas (TCGA) initiatives.

Results

We present ASEQ, a tool to perform gene-level allele-specific expression (ASE) analysis from paired genomic and transcriptomic NGS data without requiring paternal and maternal genome data. ASEQ offers an easy-to-use set of modes that transparently to the user takes full advantage of a built-in fast computational engine. We report its performances on a set of 20 individuals from the 1,000 Genomes Project and show its detection power on imprinted genes. Next we demonstrate high level of ASE calls concordance when comparing it to AlleleSeq and MBASED tools. Finally, using a prostate cancer dataset we report on a higher fraction of ASE genes with respect to healthy individuals and show allele-specific events nominated by ASEQ in genes that are implicated in the disease.

Conclusions

ASEQ can be used to rapidly and reliably screen large NGS datasets for the identification of allele specific features. It can be integrated in any NGS pipeline and runs on computer systems with multiple CPUs, CPUs with multiple cores or across clusters of machines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0084-2) contains supplementary material, which is available to authorized users.

The molecular genetic architecture of attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common childhood behavioral condition which affects 2-10% of school age children worldwide. Although the underlying molecular mechanism for the disorder is poorly understood, familial, twin and adoption studies suggest a strong genetic component. Here we provide a state-of-the-art review of the molecular genetics of ADHD incorporating evidence from candidate gene and linkage designs, as well as genome-wide association (GWA) studies of common single-nucleotide polymorphisms (SNPs) and rare copy number variations (CNVs). Bioinformatic methods such as functional enrichment analysis and protein-protein network analysis are used to highlight biological processes of likely relevance to the aetiology of ADHD. Candidate gene associations of minor effect size have been replicated across a number of genes including SLC6A3, DRD5, DRD4, SLC6A4, LPHN3, SNAP-25, HTR1B, NOS1 and GIT1. Although case-control SNP-GWAS have had limited success in identifying common genetic variants for ADHD that surpass critical significance thresholds, quantitative trait designs suggest promising associations with Cadherin13 and glucose-fructose oxidoreductase domain 1 genes. Further, CNVs mapped to glutamate receptor genes (GRM1, GRM5, GRM7 and GRM8) have been implicated in the aetiology of the disorder and overlap with bioinformatic predictions based on ADHD GWAS SNP data regarding enriched pathways. Although increases in sample size across multi-center cohorts will likely yield important new results, we advocate that this must occur in parallel with a shift away from categorical case-control approaches that view ADHD as a unitary construct, towards dimensional approaches that incorporate endophenotypes and statistical classification methods.

Opportunities and Limitations of Modelling Alzheimer's Disease with Induced Pluripotent Stem Cells

Reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) has opened the way for patient-specific disease modelling. Following their differentiation into neuronal cell types, iPSC have enabled the investigation of human neurodegenerative diseases, such as Alzheimer's disease (AD). While human iPSCs certainly provide great opportunities to repeatedly interrogate specific human brain cell types of individuals with familial and sporadic forms of the disease, the complex aetiology and timescale over which AD develops in humans poses particular challenges to iPSC-based AD models. Here, we discuss the current state-of-play in the context of these and other iPSC model-related challenges and elaborate on likely future developments in this field of research.

Global gene expression analysis in the zebrafish: the challenge and the promise

Microarray technology allows the simultaneous changes in multiple gene expression to be determined as part of a drug discovery platform, or to understand the biological mechanisms of disease. In terms of both of these objectives, those studying the zebrafish as a complement to mouse-based research are only now able to entertain the construction of microarrays due to recent progress in sequencing the zebrafish genome. The advantage of the late entry of zebrafish to the microarray field offers the means of immediately adopting best practice, while also appreciating that what is currently on offer can be improved.:

The Atlantic Bonito (Sarda sarda, Bloch 1793) transcriptome and detection of differential expression during larvae development

The Atlantic bonito (Sarda sarda, Bloch 1793) belongs to the important marine fish species with a wide geographical distribution covering the Atlantic Ocean, the Mediterranean and its bordering seas. Aquaculture practices for this species are still in their infancies and scientific studies are seldom undertaken, mainly because of difficulties in sampling. Thus for small tuna species like the Atlantic bonito only little is known about its biology and regarding the molecular background even less information is available. In the production of marine fish it is known that the most critical period is the larval stages, as high growth rates as well as significant developmental changes take place. In this study we have investigated the transcriptome of the Atlantic bonito of five larvae stages applying Illumina sequencing technology. For non-model species like aquaculture species, transcriptome analysis of RNA samples from individuals using Illumina sequencing technology is technically efficient and cost effective. In the present study a total number of 169,326,711 paired-end reads with a read length of 100 base pairs were generated resulting in a reference transcriptome of 68,220 contigs with an average length of 2054 base pairs. For differential expression analyses single end reads were obtained from different developmental stages and mapped to the constructed reference transcriptome. Differential expression analyses revealed in total 18,657 differentially expressed transcripts and were assigned to five distinguished groups. Each of the five clusters shows stage specific gene expression. We present for the first time in the Atlantic bonito an extensive RNA-Seq based characterization of its transcriptome as well as significant information on differential expression among five developmental larvae stages. The generated transcripts, including SNP and microsatellite information for candidate molecular markers and gene expression information will be a valuable resource for future genetic and molecular studies.

Transcriptional signature of progesterone in the fathead minnow ovary (Pimephales promelas)

A growing number of studies have examined transcriptional responses to sex steroids along the hypothalamic-pituitary-gonadal axis in teleost fishes. However, data are lacking on the molecular cascades that underlie progesterone signaling. The objective of this study was to characterize the transcriptional response in the ovary of fathead minnows (Pimephales promelas) in response to progesterone (P4). Fathead minnow ovaries were exposed in vitro to 500 ng P4/L. Germinal vesicle migration and breakdown (GVBD) was observed and microarrays were used to identify gene cascades affected by P4. Microarray analysis identified 1702 differentially expressed transcripts after P4 treatment. Functional enrichment analysis revealed that transcripts involved in the molecular functions of protein serine/threonine kinase activity, ATP binding, and activity of calcium channels were increased after P4 treatment. There was an overwhelming decrease in levels of transcripts of genes that are structural constituents of ribosomes with P4 treatment. There was also evidence for gene expression changes in steroid and maturation-related transcripts. Pathway analyses identified cell cycle regulation, insulin action, hedgehog, and B cell activation as pathways containing an over-representation of highly regulated transcripts. Significant regulatory sub-networks of P4-mediated transcripts included genes regulated by tumor protein p53 and E2F transcription factor 1. These data provide novel insight into the molecular signaling cascades that underlie P4-signaling in the ovary and identify genes and processes that may indicate premature GVBD due to environmental pollutants that mimic progestins.

Expression of developmental-stage-specific genes in the gilthead sea bream Sparus aurata L

The mechanism of early fish development as well as the control of egg quality is of great importance for the ability of the oocyte to develop after fertilization. Embryonic development is initially regulated by maternally provided mRNAs and later by the zygotic genome. Maternal mRNAs have an important role in initiating processes crucial to patterning the developing fish embryo. Furthermore, it has been shown that maternal RNA plays an important role in egg quality. The identification and characterization of candidate maternal genes in non-model fish species with important aquaculture interest like the gilthead sea bream Sparus aurata L. is of importance for future studies related to egg quality. The broodstock of the gilthead sea bream produces large quantities of eggs with a high and non-controllable quality variation. In the present study, we have studied the gene expression of 16 genes (gapdh 1 and 2, cathepsin D, L, S and Z, erk1, jnk1, p38 alpha and p38 delta, ppar alpha, beta and gamma, tubulin beta, ferritin M, cyclinA2) of different functional categories in seven developmental stages. The 16 genes were chosen based on their putative involvement in egg quality and regulation of early development. In total, 11 showed a characteristic gene expression pattern pinpointing to the possible function as maternal genes and thus may function as molecular biomarker for egg quality.

Deep mRNA sequencing analysis to capture the transcriptome landscape of zebrafish embryos and larvae

Transcriptome analysis is a powerful tool to obtain large amount genome-scale gene expression profiles. Despite its extensive usage to diverse biological problems in the last decade, transcriptomic researches approaching the zebrafish embryonic development have been very limited. Several recent studies have made great progress in this direction, yet the large gap still exists, especially regarding to the transcriptome dynamics of embryonic stages from early gastrulation onwards. Here, we present a comprehensive analysis about the transcriptomes of 9 different stages covering 7 major periods (cleavage, blastula, gastrula, segmentation, pharyngula, hatching and early larval stage) in zebrafish development, by recruiting the RNA-sequencing technology. We detected the expression for at least 24,065 genes in at least one of the 9 stages. We identified 16,130 genes that were significantly differentially expressed between stages and were subsequently classified into six clusters. Each revealed gene cluster had distinct expression patterns and characteristic functional pathways, providing a framework for the understanding of the developmental transcriptome dynamics. Over 4000 genes were identified as preferentially expressed in one of the stages, which could be of high relevance to stage-specific developmental and molecular events. Among the 68 transcription factor families active during development, most had enhanced average expression levels and thus might be crucial for embryogenesis, whereas the inactivation of the other families was likely required by the activation of the zygotic genome. We discussed our RNA-seq data together with previous findings about the Wnt signaling pathway and some other genes with known functions, to show how our data could be used to advance our understanding about these developmental functional elements. Our study provides ample information for further study about the molecular and cellular mechanisms underlying vertebrate development.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.

Analysis of allele-specific gene expression using a target-oriented tiling microarray assay

Advances in molecular and computational biology in recent years have led to the development or the improvement of methods for analyzing global gene expression. In most of these efforts, it is assumed that alleles of different origins contribute equally to the transcript pool. However, accumulating evidence suggests that many genes are not equally expressed from the paternal and maternal chromosomes. In addition to imprinting, the phenomenon of imbalanced allelic expression is widespread in heterozygous individuals. To distinguish transcript pools derived from different alleles present in the same organism, a number of methods have been developed. Here, we describe an oligonucleotide tiling microarray-based assay for analyzing allele-specific gene expression. Specifically targeting single-nucleotide polymorphisms, this two-step assay offers a high-throughput and multiplex method for detecting and quantifying unequal allelic expression that is readily applicable to many experimental systems.

Impact of experimental noise and annotation imprecision on data quality in microarray experiments

Data quality is intrinsically influenced by design, technical, and analytical parameters. Quality parameters have not yet been well defined for gene expression analysis by microarrays, though ad interim, following recommended good experimental practice guidelines should ensure generation of reliable and reproducible data. Here we summarize essential practical recommendations for experimental design, technical considerations, feature annotation issues, and standardization efforts.

iASeq: integrative analysis of allele-specificity of protein-DNA interactions in multiple ChIP-seq datasets

BACKGROUND

ChIP-seq provides new opportunities to study allele-specific protein-DNA binding (ASB). However, detecting allelic imbalance from a single ChIP-seq dataset often has low statistical power since only sequence reads mapped to heterozygote SNPs are informative for discriminating two alleles.

RESULTS

We develop a new method iASeq to address this issue by jointly analyzing multiple ChIP-seq datasets. iASeq uses a Bayesian hierarchical mixture model to learn correlation patterns of allele-specificity among multiple proteins. Using the discovered correlation patterns, the model allows one to borrow information across datasets to improve detection of allelic imbalance. Application of iASeq to 77 ChIP-seq samples from 40 ENCODE datasets and 1 genomic DNA sample in GM12878 cells reveals that allele-specificity of multiple proteins are highly correlated, and demonstrates the ability of iASeq to improve allelic inference compared to analyzing each individual dataset separately.

CONCLUSIONS

iASeq illustrates the value of integrating multiple datasets in the allele-specificity inference and offers a new tool to better analyze ASB.

dsPIG: a tool to predict imprinted genes from the deep sequencing of whole transcriptomes

Background

Dysregulation of imprinted genes, which are expressed in a parent-of-origin-specific manner, plays an important role in various human diseases, such as cancer and behavioral disorder. To date, however, fewer than 100 imprinted genes have been identified in the human genome. The recent availability of high-throughput technology makes it possible to have large-scale prediction of imprinted genes. Here we propose a Bayesian model (dsPIG) to predict imprinted genes on the basis of allelic expression observed in mRNA-Seq data of independent human tissues.

Results

Our model (dsPIG) was capable of identifying imprinted genes with high sensitivity and specificity and a low false discovery rate when the number of sequenced tissue samples was fairly large, according to simulations. By applying dsPIG to the mRNA-Seq data, we predicted 94 imprinted genes in 20 cerebellum samples and 57 imprinted genes in 9 diverse tissue samples with expected low false discovery rates. We also assessed dsPIG using previously validated imprinted and non-imprinted genes. With simulations, we further analyzed how imbalanced allelic expression of non-imprinted genes or different minor allele frequencies affected the predictions of dsPIG. Interestingly, we found that, among biallelically expressed genes, at least 18 genes expressed significantly more transcripts from one allele than the other among different individuals and tissues.

Conclusion

With the prevalence of the mRNA-Seq technology, dsPIG has become a useful tool for analysis of allelic expression and large-scale prediction of imprinted genes. For ease of use, we have set up a web service and also provided an R package for dsPIG at http://www.shoudanliang.com/dsPIG/.

PYK2: a calcium-sensitive protein tyrosine kinase activated in response to fertilization of the zebrafish oocyte

Fertilization begins with binding and fusion of a sperm with the oocyte, a process that triggers a high amplitude calcium transient which propagates through the oocyte and stimulates a series of preprogrammed signal transduction events critical for zygote development. Identification of the pathways downstream of this calcium transient remains an important step in understanding the basis of zygote quality. The present study demonstrates that the calcium-calmodulin sensitive protein tyrosine kinase PYK2 is a target of the fertilization-induced calcium transient in the zebrafish oocyte and that it plays an important role in actin-mediated events critical for sperm incorporation. At fertilization, PYK2 was activated initially at the site of sperm-oocyte interaction and was closely associated with actin filaments forming the fertilization cone. Later PYK2 activation was evident throughout the entire oocyte cortex, however activation was most intense over the animal hemisphere. Fertilization-induced PYK2 activation could be blocked by suppressing calcium transients in the ooplasm via injection of BAPTA as a calcium chelator. PYK2 activation could be artificially induced in unfertilized oocytes by injection of IP3 at concentrations sufficient to induce calcium release. Functionally, suppression of PYK2 activity by chemical inhibition or by injection of a dominant-negative construct encoding the N-terminal ERM domain of PKY2 inhibited formation of an organized fertilization cone and reduced the frequency of successful sperm incorporation. Together, the above findings support a model in which PYK2 responds to the fertilization-induced calcium transient by promoting reorganization of the cortical actin cytoskeleton to form the fertilization cone.

Gene expression profiles of spleen, liver, and head kidney in turbot (Scophthalmus maximus) along the infection process with Philasterides dicentrarchi using an immune-enriched oligo-microarray

We evaluated the expression profiles of turbot in spleen, liver, and head kidney across five temporal points of the Philasterides dicentrarchi infection process using an 8x15K Agilent oligo-microarray. The microarray included 2,176 different fivefold replicated gene probes designed from a turbot 3' sequenced EST database. We were able to identify 221 differentially expressed (DE) genes (8.1% of the whole microarray), 113 in spleen, 83 in liver, and 90 in head kidney, in at least 1 of the 5 temporal points sampled for each organ. Most of these genes could be annotated (83.0%) and functionally categorized using GO terms (69.1%) after the additional sequencing of DE genes from the 5' end. Many DE genes were related to innate and acquired immune functions. A high proportion of DE genes were organ-specific (70.6%), although their associated GO functions showed notable similarities in the three organs. The most striking difference in functional distribution was observed between the up- and downregulated gene groups. Upregulated genes were mostly associated to immune functions, while downregulated ones mainly involved metabolism-related genes. Genetic response appeared clustered in a few groups of genes with similar expression profiles along the temporal series. The information obtained will aid to understand the turbot immune response and will specifically be valuable to develop strategies of defense to P. dicentrarchi to achieve more resistant broodstocks for turbot industry.

Transcriptome analysis of head kidney in grass carp and discovery of immune-related genes

BACKGROUND

Grass carp (Ctenopharyngodon idella) is one of the most economically important freshwater fish, but its production is often affected by diseases that cause serious economic losses. To date, no good breeding varieties have been obtained using the oriented cultivation technique. The ability to identify disease resistance genes in grass carp is important to cultivate disease-resistant varieties of grass carp.

RESULTS

In this study, we constructed a non-normalized cDNA library of head kidney in grass carp, and, after clustering and assembly, we obtained 3,027 high-quality unigenes. Solexa sequencing was used to generate sequence tags from the transcriptomes of the head kidney in grass carp before and after grass carp reovirus (GCRV) infection. After processing, we obtained 22,144 tags that were differentially expressed by more than 2-fold between the uninfected and infected groups. 679 of the differentially expressed tags (3.1%) mapped to 483 of the unigenes (16.0%). The up-regulated and down-regulated unigenes were annotated using gene ontology terms; 16 were annotated as immune-related and 42 were of unknown function having no matches to any of the sequences in the databases that were used in the similarity searches. Semi-quantitative RT-PCR revealed four unknown unigenes that showed significant responses to the viral infection. Based on domain structure predictions, one of these sequences was found to encode a protein that contained two transmembrane domains and, therefore, may be a transmembrane protein. Here, we proposed that this novel unigene may encode a virus receptor or a protein that mediates the immune signalling pathway at the cell surface.

CONCLUSION

This study enriches the molecular basis data of grass carp and further confirms that, based on fish tissue-specific EST databases, transcriptome analysis is an effective route to discover novel functional genes.

Hereditary mixed polyposis syndrome is caused by a 40-kb upstream duplication that leads to increased and ectopic expression of the BMP antagonist GREM1

The hereditary mixed polyposis syndrome (HMPS) was first described about 50 years ago in a large Ashkenazi Jewish family from St Mark’s Hospital, London. The family showed apparent autosomal dominant inheritance of multiple types of colorectal polyp, with colorectal carcinoma in a high proportion of individuals. In the last 15 years, we have mapped the HMPS gene to chromosome 15q13.3 and identified an ancestral haplotype common to all the known HMPS families. Here, we have used genetic mapping, copy number analysis, exclusion of mutations by high-throughput sequencing, gene expression analysis and functional assays to show that HMPS is caused by a large duplication spanning the 3′ end of the SCG5 gene and a region upstream of the GREM1 locus. This mutation has no effect on SCG5 expression, but is associated with greatly increased, allele-specific GREM1 expression. Whilst GREM1 is expressed in intestinal sub-epithelial myofibroblasts in controls, HMPS patients predominantly express GREM1 in the epithelium of the large bowel. The HMPS duplication contains predicted transcriptional enhancer elements; we have shown that some of these interact with the GREM1 promoter and are capable of driving gene expression in vitro. Increased GREM1 expression is predicted to lead to reduced bone morphogenetic protein pathway activity, a mechanism that also underlies tumorigenesis in juvenile polyposis of the large bowel. The pathogenic mechanism in HMPS is extremely unusual in Mendelian cancer syndromes and highlights ectopic gene expression as a mechanism of tumorigenesis.

Effects of BRCA2 cis-regulation in normal breast and cancer risk amongst BRCA2 mutation carriers

INTRODUCTION

Cis-acting regulatory single nucleotide polymorphisms (SNPs) at specific loci may modulate penetrance of germline mutations at the same loci by introducing different levels of expression of the wild-type allele. We have previously reported that BRCA2 shows differential allelic expression and we hypothesize that the known variable penetrance of BRCA2 mutations might be associated with this mechanism.

METHODS

We combined haplotype analysis and differential allelic expression of BRCA2 in breast tissue to identify expression haplotypes and candidate cis-regulatory variants. These candidate variants underwent selection based on in silico predictions for regulatory potential and disruption of transcription factor binding, and were functionally analyzed in vitro and in vivo in normal and breast cancer cell lines. SNPs tagging the expression haplotypes were correlated with the total expression of several genes in breast tissue measured by Taqman and microarray technologies. The effect of the expression haplotypes on breast cancer risk in BRCA2 mutation carriers was investigated in 2,754 carriers.

RESULTS

We identified common haplotypes associated with differences in the levels of BRCA2 expression in human breast cells. We characterized three cis-regulatory SNPs located at the promoter and two intronic regulatory elements which affect the binding of the transcription factors C/EBPα, HMGA1, D-binding protein (DBP) and ZF5. We showed that the expression haplotypes also correlated with changes in the expression of other genes in normal breast. Furthermore, there was suggestive evidence that the minor allele of SNP rs4942440, which is associated with higher BRCA2 expression, is also associated with a reduced risk of breast cancer (per-allele hazard ratio (HR) = 0.85, 95% confidence interval (CI) = 0.72 to 1.00, P-trend = 0.048).

CONCLUSIONS

Our work provides further insights into the role of cis-regulatory variation in the penetrance of disease-causing mutations. We identified small-effect genetic variants associated with allelic expression differences in BRCA2 which could possibly affect the risk in mutation carriers through altering expression levels of the wild-type allele.

Common genetic variants at the 11q13.3 renal cancer susceptibility locus influence binding of HIF to an enhancer of cyclin D1 expression

Although genome-wide association studies (GWAS) have identified the existence of numerous population-based cancer susceptibility loci, mechanistic insights remain limited, particularly for intergenic polymorphisms. Here, we show that polymorphism at a remote intergenic region on chromosome 11q13.3, recently identified as a susceptibility locus for renal cell carcinoma, modulates the binding and function of hypoxia-inducible factor (HIF) at a previously unrecognized transcriptional enhancer of CCND1 (encoding cyclin D1) that is specific for renal cancers characterized by inactivation of the von Hippel-Lindau tumor suppressor (pVHL). The protective haplotype impairs binding of HIF-2, resulting in an allelic imbalance in cyclin D1 expression, thus affecting a link between hypoxia pathways and cell cycle control.

An ovary transcriptome for all maturational stages of the striped bass (Morone saxatilis), a highly advanced perciform fish

Background

The striped bass and its relatives (genus Morone) are important fisheries and aquaculture species native to estuaries and rivers of the Atlantic coast and Gulf of Mexico in North America. To open avenues of gene expression research on reproduction and breeding of striped bass, we generated a collection of expressed sequence tags (ESTs) from a complementary DNA (cDNA) library representative of their ovarian transcriptome.

Results

Sequences of a total of 230,151 ESTs (51,259,448 bp) were acquired by Roche 454 pyrosequencing of cDNA pooled from ovarian tissues obtained at all stages of oocyte growth, at ovulation (eggs), and during preovulatory atresia. Quality filtering of ESTs allowed assembly of 11,208 high-quality contigs ≥ 100 bp, including 2,984 contigs 500 bp or longer (average length 895 bp). Blastx comparisons revealed 5,482 gene orthologues (E-value < 10^-3), of which 4,120 (36.7% of total contigs) were annotated with Gene Ontology terms (E-value < 10^-6). There were 5,726 remaining unknown unique sequences (51.1% of total contigs). All of the high-quality EST sequences are available in the National Center for Biotechnology Information (NCBI) Short Read Archive (GenBank: SRX007394). Informative contigs were considered to be abundant if they were assembled from groups of ESTs comprising ≥ 0.15% of the total short read sequences (≥ 345 reads/contig). Approximately 52.5% of these abundant contigs were predicted to have predominant ovary expression through digital differential display in silico comparisons to zebrafish (Danio rerio) UniGene orthologues. Over 1,300 Gene Ontology terms from Biological Process classes of Reproduction, Reproductive process, and Developmental process were assigned to this collection of annotated contigs.

Conclusions

This first large reference sequence database available for the ecologically and economically important temperate basses (genus Morone) provides a foundation for gene expression studies in these species. The predicted predominance of ovary gene expression and assignment of directly relevant Gene Ontology classes suggests a powerful utility of this dataset for analysis of ovarian gene expression related to fundamental questions of oogenesis. Additionally, a high definition Agilent 60-mer oligo ovary 'UniClone' microarray with 8 × 15,000 probe format has been designed based on this striped bass transcriptome (eArray Group: Striper Group, Design ID: 029004).

Targeted mutation of the talpid3 gene in zebrafish reveals its conserved requirement for ciliogenesis and Hedgehog signalling across the vertebrates

Using zinc-finger nuclease-mediated mutagenesis, we have generated mutant alleles of the zebrafish orthologue of the chicken talpid3 (ta3) gene, which encodes a centrosomal protein that is essential for ciliogenesis. Animals homozygous for these mutant alleles complete embryogenesis normally, but manifest a cystic kidney phenotype during the early larval stages and die within a month of hatching. Elimination of maternally derived Ta3 activity by germline replacement resulted in embryonic lethality of ta3 homozygotes. The phenotype of such maternal and zygotic (MZta3) mutant zebrafish showed strong similarities to that of chick ta3 mutants: absence of primary and motile cilia as well as aberrant Hedgehog (Hh) signalling, the latter manifest by the expanded domains of engrailed and ptc1 expression in the somites, reduction of nkx2.2 expression in the neural tube, symmetric pectoral fins, cyclopic eyes and an ectopic lens. GFP-tagged Gli2a localised to the basal bodies in the absence of the primary cilia and western blot analysis showed that Gli2a protein is aberrantly processed in MZta3 embryos. Zygotic expression of ta3 largely rescued the effects of maternal depletion, but the motile cilia of Kupffer's vesicle remained aberrant, resulting in laterality defects. Our findings underline the importance of the primary cilium for Hh signaling in zebrafish and reveal the conservation of Ta3 function during vertebrate evolution.

The 14q22.2 colorectal cancer variant rs4444235 shows cis-acting regulation of BMP4

Common genetic variation at human 14q22.2 tagged by rs4444235 is significantly associated with colorectal cancer (CRC) risk. Re-sequencing was used to comprehensively annotate the 17kb region of strong linkage disequilibrium encompassing rs4444235. Through bioinformatic analyses using H3K4Me1, H3K4Me3, and DNase-I hypersensitivity chromatin signatures and evolutionary conservation we identified seven candidate disease-causing single-nucleotide polymorphisms mapping to six regions within the 17-kb region predicted to have regulatory potential. Reporter gene studies of these regions demonstrated that the element to which rs4444235 maps acts as an allele-specific transcriptional enhancer. Allele-specific expression studies in CRC cell lines heterozygous for rs4444235 showed significantly increased expression of bone morphogenetic protein-4 (BMP4) associated with the risk allele (P<0.001). These data provide evidence for a functional basis for the non-coding risk variant rs4444235 at 14q22.2 and emphasizes the importance of genetic variation in the BMP pathway genes as determinants of CRC risk.

Gene expression profiles of the spleen, liver, and head kidney in turbot (Scophthalmus maximus) along the infection process with Aeromonas salmonicida using an immune-enriched oligo-microarray

We evaluated the expression profiles of turbot in the spleen, liver, and head kidney across five temporal points of the Aeromonas salmonicida infection process using an 8 × 15 K Agilent oligo-microarray. The microarray included 2,176 different fivefold replicated gene probes designed from a turbot 3' sequenced EST database. We were able to identify 471 differentially expressed (DE) genes (17.3% of the whole microarray), 223 in the spleen, 246 in the liver, and 125 in the head kidney, in at least one of the five temporal points sampled for each organ. Most of these genes could be annotated (83.0%) and functionally categorized using Gene Ontology terms (69.1%) after the additional sequencing of DE genes from the 5' end. Many DE genes were related to innate and acquired immune functions in accordance to previous studies with this pathogen in other fish species. A high proportion of DE genes were organ specific (77.1%), but their associated GO functions were rather similar in the three organs. The most striking difference in functional distribution was observed between the up- and down-regulated gene groups. Up-regulated genes were mostly associated to key immune functions while down-regulated ones mainly involved metabolism- and transport-related genes. Genetic response appeared clustered in groups of genes with similar expression profiles along the temporal series. The spleen showed the most clustering while the liver and head kidney displayed a higher diversification. The information obtained will aid to understand the turbot immune response and will specifically be valuable to develop strategies of defense to A. salmonicida to achieve more resistant broodstocks for turbot industry.

Pharmacogenetics and individualized therapy in children: immunosuppressants, antidepressants, anticancer and anti-inflammatory drugs

Pharmacogenetic polymorphisms that change the amino acid sequences in coding regions only account for part of the interindividual differences in disease susceptibility and drug response. Additional pharmacogenomic and epigenetic factors are also involved. In children, pharmacogenetic studies are limited, although it has been clear for many years that the interactions between developmental patterns of drug-metabolizing enzymes and transporters have a major impact on dose exposure with age-specific dosage requirements. This article will analyze the factors affecting variability in drug response in children and focus on the pharmacogenetic polymorphisms of immunosuppressants, antidepressants, anticancer and anti-inflammatory drugs. Additional pharmacogenetic and epigenetic studies should be performed to allow the individualization of therapy in children.