Sequencing refractory regions in bird genomes are hotspots for accelerated protein evolution.

BACKGROUND: Approximately 1000 protein encoding genes common for vertebrates are still unannotated in avian genomes. Are these genes evolutionary lost or are they not yet found for technical reasons? Using genome landscapes as a tool to visualize large-scale regional effects of genome evolution, we reexamined this question. RESULTS: On basis of gene annotation in non-avian vertebrate genomes, we established a list of 15,135 common vertebrate genes. Of these, 1026 were not found in any of eight examined bird genomes. Visualizing regional genome effects by our sliding window approach showed that the majority of these "missing" genes can be clustered to 14 regions of the human reference genome. In these clusters, an additional 1517 genes (often gene fragments) were underrepresented in bird genomes. The clusters of "missing" genes coincided with regions of very high GC content, particularly in avian genomes, making them "hidden" because of incomplete sequencing. Moreover, proteins encoded by genes in these sequencing refractory regions showed signs of accelerated protein evolution. As a proof of principle for this idea we experimentally characterized the mRNA and protein products of four "hidden" bird genes that are crucial for energy homeostasis in skeletal muscle: ALDOA, ENO3, PYGM and SLC2A4. CONCLUSIONS: A least part of the "missing" genes in bird genomes can be attributed to an artifact caused by the difficulty to sequence regions with extreme GC% ("hidden" genes). Biologically, these "hidden" genes are of interest as they encode proteins that evolve more rapidly than the genome wide average. Finally we show that four of these "hidden" genes encode key proteins for energy metabolism in flight muscle.

Functional evaluation of prevascularization in one-stage versus two-stage tissue engineering approach of human bio-artificial muscle.

A common shortcoming of current tissue engineered constructs is the lack of a functional vasculature, limiting their size and functionality. Prevascularization is a possible strategy to introduce vascular networks in these constructs. It includes among others co-culturing target cells with endothelial (precursor) cells that are able to form endothelial networks through vasculogenesis. In this paper, we compared two different prevascularization approaches of bio-artificial skeletal muscle tissue (BAM) in vitro and in vivo. In a one-stage approach, human muscle cells were directly co-cultured with endothelial cells in 3D. In a two-stage approach, a one week old BAM containing differentiated myotubes was coated with a fibrin hydrogel containing endothelial cells. The obtained endothelial networks were longer and better interconnected with the two-stage approach. We evaluated whether prevascularization had a beneficial effect on in vivo perfusion of the BAM and improved myotube survival by implantation on the fascia of the latissimus dorsi muscle of NOD/SCID mice for 5 or 14 d. Also in vivo, the two-stage approach displayed the highest vascular density. At day 14, anastomosis of implanted endothelial networks with the host vasculature was apparent. BAMs without endothelial networks contained longer and thicker myotubes in vitro, but their morphology degraded in vivo. In contrast, maintenance of myotube morphology was well supported in the two-stage prevascularized BAMs. To conclude, a two-stage prevascularization approach for muscle engineering improved the vascular density in the construct and supported myotube maintenance in vivo.

Vascularization of tissue-engineered skeletal muscle constructs.

Skeletal muscle tissue can be created in vitro by tissue engineering approaches, based on differentiation of muscle stem cells. Several approaches exist and generally result in three dimensional constructs composed of multinucleated myofibers to which we refer as myooids. Engineering methods date back to 3 decades ago and meanwhile a wide range of cell types and scaffold types have been evaluated. Nevertheless, in most approaches, myooids remain very small to allow for diffusion-mediated nutrient supply and waste product removal, typically less than 1 mm thick. One of the shortcomings of current in vitro skeletal muscle organoid development is the lack of a functional vascular structure, thus limiting the size of myooids. This is a challenge which is nowadays applicable to almost all organoid systems. Several approaches to obtain a vascular structure within myooids have been proposed. The purpose of this review is to give a concise overview of these approaches.

GC content of vertebrate exome landscapes reveal areas of accelerated protein evolution.

BACKGROUND: Rapid accumulation of vertebrate genome sequences render comparative genomics a powerful approach to study macro-evolutionary events. The assessment of phylogenic relationships between species routinely depends on the analysis of sequence homology at the nucleotide or protein level. RESULTS: We analyzed mRNA GC content, codon usage and divergence of orthologous proteins in 55 vertebrate genomes. Data were visualized in genome-wide landscapes using a sliding window approach. Landscapes of GC content reveal both evolutionary conservation of clustered genes, and lineage-specific changes, so that it was possible to construct a phylogenetic tree that closely matched the classic "tree of life". Landscapes of GC content also strongly correlated to landscapes of amino acid usage: positive correlation with glycine, alanine, arginine and proline and negative correlation with phenylalanine, tyrosine, methionine, isoleucine, asparagine and lysine. Peaks of GC content correlated strongly with increased protein divergence. CONCLUSIONS: Landscapes of base- and amino acid composition of the coding genome opens a new approach in comparative genomics, allowing identification of discrete regions in which protein evolution accelerated over deep evolutionary time. Insight in the evolution of genome structure may spur novel studies assessing the evolutionary benefit of genes in particular genomic regions.

Human tissue-engineered skeletal muscle: a novel 3D in vitro model for drug disposition and toxicity after intramuscular injection.

The development of laboratory-grown tissues, referred to as organoids, bio-artificial tissue or tissue-engineered constructs, is clearly expanding. We describe for the first time how engineered human muscles can be applied as a pre- or non-clinical model for intramuscular drug injection to further decrease and complement the use of in vivo animal studies. The human bio-artificial muscle (BAM) is formed in a seven day tissue engineering procedure during which human myoblasts fuse and differentiate to aligned myofibers in an extracellular matrix. The dimensions of the BAM constructs allow for injection and follow-up during several days after injection. A stereotactic setup allows controllable injection at multiple sites in the BAM. We injected several compounds; a dye, a hydrolysable compound, a reducible substrate and a wasp venom toxin. Afterwards, direct reflux, release and metabolism were assessed in the BAM constructs in comparison to 2D cell culture and isolated human muscle strips. Spectrophotometry and luminescence allowed to measure the release of the injected compounds and their metabolites over time. A release profile over 40 hours was observed in the BAM model in contrast to 2D cell culture, showing the capacity of the BAM model to function as a drug depot. We also determined compound toxicity on the BAMs by measuring creatine kinase release in the medium, which increased with increasing toxic insult. Taken together, we show that the BAM is an injectable human 3D cell culture model that can be used to measure release and metabolism of injected compounds in vitro.

Endoplasmic reticulum-associated degradation of the mouse PC1/3-N222D hypomorph and human PCSK1 mutations contributes to obesity.

BACKGROUND: The proprotein convertase 1/3 (PC1/3), encoded by proprotein convertase subtilisin/kexin type 1 (PCSK1), cleaves and hence activates several orexigenic and anorexigenic proproteins. Congenital inactivation of PCSK1 leads to obesity in human but not in mice. However, a mouse model harboring the hypomorphic mutation N222D is obese. It is not clear why the mouse models differ in phenotype. METHODS: Gene expression analysis was performed with pancreatic islets from Pcsk1(N222D/N222D) mice. Subsequently, biosynthesis, maturation, degradation and activity were studied in islets, pituitary, hypothalamus and cell lines. Coimmunoprecipitation of PC1/3-N222D and human PC1/3 variants associated with obesity with the endoplasmic reticulum (ER) chaperone BiP was studied in cell lines. RESULTS: Gene expression analysis of islets of Pcsk1(N222D/N222D) mice showed enrichment of gene sets related to the proteasome and the unfolded protein response. Steady-state levels of PC1/3-N222D and in particular the carboxy-terminally processed form were strongly reduced in islets, pituitary and hypothalamus. However, impairment of substrate cleavage was tissue dependent. Proinsulin processing was drastically reduced, while processing of proopiomelanocortin (POMC) to adrenocorticotropic hormone (ACTH) in pituitary was only mildly impaired. Growth hormone expression and IGF-1 levels were normal, indicating near-normal processing of hypothalamic proGHRH. PC1/3-N222D binds to BiP and is rapidly degraded by the proteasome. Analysis of human PC1/3 obesity-associated mutations showed increased binding to BiP and prolonged intracellular retention for all investigated mutations, in particular for PC1/3-T175M, PC1/3-G226R and PC1/3-G593R. CONCLUSIONS: This study demonstrates that the hypomorphic mutation in Pcsk1(N222D) mice has an effect on catalytic activity in pancreatic islets, pituitary and hypothalamus. Reduced substrate processing activity in Pcsk1(N222D/N222D) mice is due to enhanced degradation in addition to reduced catalytic activity of the mutant. PC1/3-N222D binds to BiP, suggesting impaired folding and reduced stability. Enhanced BiP binding is also observed in several human obesity-associated PC1/3 variants, suggesting a common mechanism.

Tissue clearing for confocal imaging of native and bio-artificial skeletal muscle.

Novel clearing techniques have revolutionized three-dimensional confocal imaging of the brain without the need for physical tissue sectioning. We evaluated three clearing methods, ScaleA2, Clear(T2), and 3DISCO for visualizing native and tissue engineered muscle by confocal microscopy. We found that Clear(T2) treatment improved the depth of visualization of immunohistochemical staining slightly, but did not improve depth of visualization of endogenous green fluorescent protein (GFP). ScaleA2 preserved endogenous GFP signal better and permitted significantly deeper GFP imaging, but it was incompatible with tropomyosin immunohistochemical staining. 3DISCO treatment preserved both endogenous GFP and immunohistochemical staining, and permitted significantly deeper imaging. Clearing time for the 3DISCO procedure is short compared to ScaleA2 and Clear(T2). We suggest that 3DISCO is the preferable clearing method for native and tissue engineered skeletal muscle tissue.

mRNA expression analysis of cell cycle genes in islets of pregnant mice.

AIMS/HYPOTHESIS: Pregnancy requires an increase in the functional beta cell mass to match metabolic needs for insulin. To understand this adaptation at the molecular level, we undertook a time course analysis of mRNA expression in mice. METHODS: Total RNA extracted from C57Bl6/J mouse islets every 3 days during pregnancy was hybridised on commercially available expression arrays. Gene network analysis was performed and changes in functional clusters over time visualised. The function of putative novel cell cycle genes was assessed via silencing in replicating mouse insulinoma 6 (MIN6) cells. RESULTS: Gene network analysis identified a large gene cluster associated with cell cycle control (67 genes, all upregulated by ≥ 1.5-fold, p < 0.001). The number of upregulated cell cycle genes and the mRNA expression levels of individual genes peaked at pregnancy day (P)9.5. Filtering of poorly annotated genes with enhanced expression in islets at P9.5, and in MIN6 cells and thymus resulted in further studies with G7e (also known as D17H6S56E-5) and Fignl1. Gene knock-down experiments in MIN6 cells suggested that these genes are indeed involved in adequate cell cycle accomplishment. CONCLUSIONS/INTERPRETATION: A sharp peak of cell cycle-related mRNA expression in islets occurs around P9.5, after which beta cell replication is increased. As illustrated by the identification of G7e and Fignl1 in islets of pregnant mice, further study of this distinct transcriptional peak should help to unravel the complex process of beta cell replication.

Testing the hypothesis of tissue selectivity: the intersection-union test and a Bayesian approach.

MOTIVATION: Finding genes that are preferentially expressed in a particular tissue or condition is a problem that cannot be solved by standard statistical testing procedures. A relatively unknown procedure that can be used is the intersection-union test (IUT). However, two disadvantages of the IUT are that it is conservative and it conveys only the information of the least differing target tissue-other tissue pair. RESULTS: We propose a Bayesian procedure that quantifies how much evidence there is in the overall expression profile for selective over-expression. In a small simulation study, it is shown that the proposed method outperforms the IUT when it comes to finding selectively expressed genes. An application to publicly available data consisting of 22 tissues shows that the Bayesian method indeed selects genes with functions that reflect the specific tissue functions. The proposed method can also be used to find genes that are underexpressed in a particular tissue. AVAILABILITY: Both MATLAB and R code that implement the IUT and the Bayesian procedure in an efficient way, can be downloaded at http://ppw.kuleuven.be/okp/software/BayesianIUT/.

Mucosal gene signatures to predict response to infliximab in patients with ulcerative colitis.

BACKGROUND AND AIMS: Infliximab is an effective treatment for ulcerative colitis with over 60% of patients responding to treatment and up to 30% reaching remission. The mechanism of resistance to anti-tumour necrosis factor alpha (anti-TNFalpha) is unknown. This study used colonic mucosal gene expression to provide a predictive response signature for infliximab treatment in ulcerative colitis. METHODS: Two cohorts of patients who received their first treatment with infliximab for refractory ulcerative colitis were studied. Response to infliximab was defined as endoscopic and histological healing. Total RNA from pre-treatment colonic mucosal biopsies was analysed with Affymetrix Human Genome U133 Plus 2.0 Arrays. Quantitative RT-PCR was used to confirm microarray data. RESULTS: For predicting response to infliximab treatment, pre-treatment colonic mucosal expression profiles were compared for responders and non-responders. Comparative analysis identified 179 differentially expressed probe sets in cohort A and 361 in cohort B with an overlap of 74 probe sets, representing 53 known genes, between both analyses. Comparative analysis of both cohorts combined, yielded 212 differentially expressed probe sets. The top five differentially expressed genes in a combined analysis of both cohorts were osteoprotegerin, stanniocalcin-1, prostaglandin-endoperoxide synthase 2, interleukin 13 receptor alpha 2 and interleukin 11. All proteins encoded by these genes are involved in the adaptive immune response. These markers separated responders from non-responders with 95% sensitivity and 85% specificity. CONCLUSION: Gene array studies of ulcerative colitis mucosal biopsies identified predictive panels of genes for (non-)response to infliximab. Further study of the pathways involved should allow a better understanding of the mechanisms of resistance to infliximab therapy in ulcerative colitis. ClinicalTrials.gov number, NCT00639821.

Efficacy and safety of adeno-associated viral vectors based on serotype 8 and 9 vs. lentiviral vectors for hemophilia B gene therapy.

BACKGROUND: Adeno-associated viral (AAV) and lentiviral vectors are promising vectors for gene therapy for hemophilia because they are devoid of viral genes and have the potential for long-term gene expression. OBJECTIVES: To compare the performance of different AAV serotypes (AAV8 and AAV9) vs. lentiviral vectors expressing factor (F) IX. METHODS AND RESULTS: AAV-based and lentiviral vectors were generated that express FIX from the same hepatocyte-specific expression cassette. AAV9 transduced the liver as efficiently as AAV8 and resulted in supra-physiological FIX levels (3000-6000% of normal) stably correcting the bleeding diathesis. Surprisingly, AAV9 resulted in unprecedented and widespread cardiac gene transfer, which was more efficient than with AAV8. AAV8 and AAV9 were not associated with any proinflammatory cytokine induction, in accordance with their minimal interactions with innate immune effectors. In contrast, lentiviral transduction resulted in modest and stable FIX levels near the therapeutic threshold (1%) and triggered a rapid self-limiting proinflammatory response (interleukin-6), which probably reflected their ability to efficiently interact with the innate immune system. CONCLUSIONS: AAV8 and 9 result in significantly higher FIX expression levels and have a reduced proinflammatory risk in comparison with lentiviral vectors. The unexpected cardiotropic properties of AAV9 have implications for gene therapy for heart disease.