Parsing parallel evolution:Ecological divergence and differential gene expression in the adaptive radiations of thick-lipped Midas cichlid fishes from Nicaragua

The study of parallel evolution facilitates the discovery of common rules of diversification. Here, we examine the repeated evolution of thick lips in Midas cichlid fishes (the Amphilophus citrinellus species complex) - from two Great Lakes and two crater lakes in Nicaragua - to assess whether similar changes in ecology, phenotypic trophic traits and gene expression accompany parallel trait evolution. Using next-generation sequencing technology, we characterize transcriptome-wide differential gene expression in the lips of wild-caught sympatric thick- and thin-lipped cichlids from all four instances of repeated thick-lip evolution. Six genes (apolipoprotein D, myelin-associated glycoprotein precursor, four-and-a-half LIM domain protein 2, calpain-9, GTPase IMAP family member 8-like and one hypothetical protein) are significantly underexpressed in the thick-lipped morph across all four lakes. However, other aspects of lips' gene expression in sympatric morphs differ in a lake-specific pattern, including the magnitude of differentially expressed genes (97-510). Generally, fewer genes are differentially expressed among morphs in the younger crater lakes than in those from the older Great Lakes. Body shape, lower pharyngeal jaw size and shape, and stable isotopes (dC and dN) differ between all sympatric morphs, with the greatest differentiation in the Great Lake Nicaragua. Some ecological traits evolve in parallel (those related to foraging ecology; e.g. lip size, body and head shape) but others, somewhat surprisingly, do not (those related to diet and food processing; e.g. jaw size and shape, stable isotopes). Taken together, this case of parallelism among thick- and thin-lipped cichlids shows a mosaic pattern of parallel and nonparallel evolution. © 2012 Blackwell Publishing Ltd.

Small RNA-mediated chromatin silencing directed to the 3' region of the Arabidopsis gene encoding the developmental regulator, FLC.

Small RNA-mediated chromatin silencing is well characterized for repeated sequences and transposons, but its role in regulating single-copy endogenous genes is unclear. We have identified two small RNAs (30 and 24 nucleotides) corresponding to the reverse strand 3' to the canonical poly(A) site of FLOWERING LOCUS C (FLC), an Arabidopsis gene encoding a repressor of flowering. Genome searches suggest that these RNAs originate from the FLC locus in a genomic region lacking repeats. The 24-nt small RNA, which is most abundant in developing fruits, is absent in mutants defective in RNA polymerase IVa, RNA-DEPENDENT RNA POLYMERASE 2, and DICER-LIKE 3, components required for RNAi-mediated chromatin silencing. The corresponding genomic region shows histone 3 lysine 9 dimethylation, which was reduced in a dcl2,3,4 triple mutant. Investigations into the origins of the small RNAs revealed a polymerase IVa-dependent spliced, antisense transcript covering the 3' FLC region. Mutation of this genomic region by T-DNA insertion led to FLC misexpression and delayed flowering, suggesting that RNAi-mediated chromatin modification is an important component of endogenous pathways that function to suppress FLC expression.

Variants in the catechol-o-methyltransferase (COMT) gene are associated with schizophrenia in Irish high-density families

The enzyme catechol-o-methyltransferase (COMT) transfers a methyl group from adenosylmethionine to catecholamines including the neurotransmitters dopamine, epinephrine and norepinephrine. This methylation results in the degradation of catecholamines. The involvement of the COMT gene in the metabolic pathway of these neurotransmitters has made it an attractive candidate gene for many psychiatric disorders. In this article, we reported our study of association of COMT with schizophrenia in Irish families with a high density of schizophrenia. Three single nucleotide polymorphisms (SNPs) were genotyped for the 274 such families and within-family transmission disequilibrium tests were performed. SNP rs4680, which is the functional Val/Met polymorphism, showed modest association with the disease by the TRANSMIT, FBAT and PDT programs, while the other two SNPs were negative. These SNPs showed lower level of LDs with each other in the Irish subjects than in Ashkenazi Jews. Haplotype analysis indicated that a haplotype, haplotype A-G-A for SNPs rs737865-rs4680-rs165599, was preferentially transmitted to the affected subjects. This was different from the reported G-G-G haplotype found in Ashkenazi Jews, but both haplotypes shared the Val allele. We concluded that COMT gene is associated with schizophrenia and carries a small but significant risk to the susceptibility in the Irish subjects.

Vitamin D receptor gene variants and esophageal adenocarcinoma risk:A population-based case-control study

Purpose: Polymorphisms in the vitamin D receptor (VDR) gene may be of etiological importance in determining cancer risk. The aim of this study was to assess the association between common VDR gene polymorphisms and esophageal adenocarcinoma (EAC) risk in an all-Ireland population-based case-control study. Methods: EAC cases and frequency-matched controls by age and gender recruited between March 2002 and December 2004 throughout Ireland were included. Participants were interviewed, and a blood sample collected for DNA extraction. Twenty-seven single nucleotide polymorphisms in the VDR gene were genotyped using Sequenom or TaqMan assays while the poly(A) microsatellite was genotyped by fluorescent fragment analysis. Unconditional logistic regression was applied to assess the association between VDR polymorphisms and EAC risk. Results: A total of 224 cases of EAC and 256 controls were involved in analyses. After adjustment for potential confounders, TT homozygotes at rs2238139 and rs2107301 had significantly reduced risks of EAC compared with CC homozygotes. In contrast, SS alleles of the poly(A) microsatellite had significantly elevated risks of EAC compared with SL/LL alleles. However, following permutation analyses to adjust for multiple comparisons, no significant associations were observed between any VDR gene polymorphism and EAC risk. Conclusions: VDR gene polymorphisms were not significantly associated with EAC development in this Irish population. Confirmation is required from larger studies. © Springer Science+Business Media, LLC 2011.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Retroviral-mediated gene transfer of the human erythropoietin gene into a factor-dependent cell line, TF-1

The factor-dependent cell line, TF-1, established from a patient with erythroleukaemia, shows characteristics of immature erythroblasts. Addition of granulocyte-macrophage colony stimulating factor (GM-CSF) to the culture medium is required for long-term growth of the cells. Erythropoietin (Epo) can also be used to sustain TF-1 cells but for only limited periods (approximately a week). Low levels of both growth factors can act synergistically to maintain proliferation for a longer period of time than Epo alone. To eliminate the requirement of exogenous Epo for growth, TF-1 cells were co-cultured with a retroviral secreting cell line containing the human erythropoietin (hEpo) gene and a neomycin (neo) selectable marker. TF-1 cells which exhibited neo resistance (indicating infection by the retrovirus) were then grown in low concentrations of GM-CSF without the addition of Epo. Under these conditions growth of normal TF-1 cells was not sustained. The neo-resistant cells survived for more than 14 days indicating synergy between GM-CSF and the Epo synthesised by the co-cultured TF-1 cells. Radioimmunoassays performed on growth media detected concentrations up to 1 mU/ml of Epo, implying that stable integration of the retroviral vector and expression of the hEpo gene have been achieved.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100kb), rare copy number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1,564 cases and 1,748 controls all from Ireland, and further extended the analysis to include an additional 5,196 UK controls. We found association with duplications at chr20p12.2 (P=0.007) and evidence of replication in large independent European schizophrenia (P=0.052) and UK bipolar disorder case-control cohorts (P=0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11,707 cases and 10 carriers in 21,204 controls (meta-analysis CMH P value=2x10(-4) (odds ratio (OR)=11.3, 95% CI=3.7, ∞)). Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68Mb with similar breakpoints across samples. By haplotype analysis and sequencing we identified a tandem ∼149kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P=2.5x10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

Identification of a collagen type I adhesin of Bacteroides fragilis

Bacteroides fragilis is an opportunistic pathogen which can cause life threatening infections in humans and animals. The ability to adhere to components of the extracellular matrix, including collagen, is related to bacterial host colonisation. Collagen Far Western analysis of the B. fragilis outer membrane protein (OMP) fraction revealed the presence two collagen adhesin bands of ∼31 and ∼34 kDa. The collagen adhesins in the OMP fraction were separated and isolated by two-dimensional SDS-PAGE and also purified by collagen affinity chromatography. The collagen binding proteins isolated by both these independent methods were subjected to tandem mass spectroscopy for peptide identification and matched to a single hypothetical protein encoded by B. fragilis NCTC 9343 (BF0586), conserved in YCH46 (BF0662) and 638R (BF0633) and which is designated in this study as cbp1 (collagen binding protein). Functionality of the protein was confirmed by targeted insertional mutagenesis of the cbp1 gene in B. fragilis GSH18 which resulted in the specific loss of both the ∼31 kDa and the ∼34 kDa adhesin bands. Purified his-tagged Cbp1, expressed in a B. fragilis wild-type and a glycosylation deficient mutant, confirmed that the cbp1 gene encoded the observed collagen adhesin, and showed that the 34 kDa band represents a glycosylated version of the ∼31 kDa protein. Glycosylation did not appear to be required for binding collagen. This study is the first to report the presence of collagen type I adhesin proteins in B. fragilis and to functionally identify a gene encoding a collagen binding protein. © 2014 Galvão et al.

Validation and application of a reporter gene assay for the determination of estrogenic endocrine disruptor activity in milk

Endocrine disruptors (EDs) are compounds known to interfere with the endocrine system by disturbing the action or pathways of natural hormones which may lead to infertility or cancer.Our diet is considered to be one of the main exposure routes to EDs. Since milk and dairy products are major components of our diet they should be monitored for ED contamination. Most assays developed to date utilise targeted, chromatography based methods which lack information on the biological activity and mixture effects of the monitored compounds.A biological reporter gene assay (RGA) was developed to assess the total estrogen hormonal load in milk. It has been validated according to EU decision 2002/657/EC. Analytes were extracted by liquid-liquid extraction with acetonitrile followed by clean up on a HLB column which yielded good recovery and small matrix effects. The method has been shown to be estrogen specific, repeatable and reproducible, with covariance values below 20%. In conclusion, this method enables the detection of low levels of estrogen hormonal activity in milk with a detection capability of 36pgg EEQ and has been successfully applied in testing a range of milk samples. © 2014 Elsevier Ltd.

RNAi dynamics in Juvenile Fasciola spp. Liver flukes reveals the persistence of gene silencing in vitro

Background: Fasciola spp. liver fluke cause pernicious disease in humans and animals. Whilst current control is unsustainable due to anthelmintic resistance, gene silencing (RNA interference, RNAi) has the potential to contribute to functional validation of new therapeutic targets. The susceptibility of juvenile Fasciola hepatica to double stranded (ds)RNA-induced RNAi has been reported. To exploit this we probe RNAi dynamics, penetrance and persistence with the aim of building a robust platform for reverse genetics in liver fluke. We describe development of standardised RNAi protocols for a commercially-available liver fluke strain (the US Pacific North West Wild Strain), validated via robust transcriptional silencing of seven virulence genes, with in-depth experimental optimisation of three: cathepsin L (FheCatL) and B (FheCatB) cysteine proteases, and a σ-class glutathione transferase (FheσGST).

Methodology/Principal Findings: Robust transcriptional silencing of targets in both F. hepatica and Fasciola gigantica juveniles is achievable following exposure to long (200–320 nt) dsRNAs or 27 nt short interfering (si)RNAs. Although juveniles are highly RNAi-susceptible, they display slower transcript and protein knockdown dynamics than those reported previously. Knockdown was detectable following as little as 4h exposure to trigger (target-dependent) and in all cases silencing persisted for ≥25 days following long dsRNA exposure. Combinatorial silencing of three targets by mixing multiple long dsRNAs was similarly efficient. Despite profound transcriptional suppression, we found a significant time-lag before the occurrence of protein suppression; FheσGST and FheCatL protein suppression were only detectable after 9 and 21 days, respectively.

Conclusions/Significance: In spite of marked variation in knockdown dynamics, we find that a transient exposure to long dsRNA or siRNA triggers robust RNAi penetrance and persistence in liver fluke NEJs supporting the development of multiple-throughput phenotypic screens for control target validation. RNAi persistence in fluke encourages in vivo studies on gene function using worms exposed to RNAi-triggers prior to infection.

Type-I interferons induce lung protease responses following respiratory syncytial virus infection via RIG-I-like receptors

The role of proteases in viral infection of the lung is poorly understood. Thus, we examined matrix metalloproteinases (MMPs) and cathepsin proteases in respiratory syncytial virus (RSV)-infected mouse lungs. RSV-induced gene expression for MMPs -2, -3, -7, -8, -9, -10, -12, -13, -14, -16, -17, -19, -20, -25, -27, and -28 and cathepsins B, C, E, G, H, K, L1, S, W, and Z in the airways of Friend leukemia virus B sensitive strain mice. Increased proteases were present in the bronchoalveolar lavage fluid (BALF) and lung tissue during infection. Mitochondrial antiviral-signaling protein (MAVS) and TIR-domain-containing adapter-inducing interferon-β-deficient mice were exposed to RSV. Mavs-deficient mice had significantly lower expression of airway MMP-2, -3, -7, -8, -9, -10, -12, -13, and -28 and cathepsins C, G, K, S, W, and Z. In lung epithelial cells, retinoic acid-inducible gene-1 (RIG-I) was identified as the major RIG-I-like receptor required for RSV-induced protease expression via MAVS. Overexpression of RIG-I or treatment with interferon-β in these cells induced MMP and cathepsin gene and protein expression. The significance of RIG-1 protease induction was demonstrated by the fact that inhibiting proteases with batimastat, E64 or ribavirin prevented airway hyperresponsiveness and enhanced viral clearance in RSV-infected mice.

Exendin-4 protects against post-myocardial infarction remodelling via specific actions on inflammation and the extracellular matrix

Glucagon-like peptide-1 (GLP-1) is an insulin-releasing hormone clinically exploited for glycaemic control in diabetes, which also confers acute cardioprotection and benefits in experimental/clinical heart failure. We specifically investigated the role of the GLP-1 mimetic, exendin-4, in post-myocardial infarction (MI) remodelling, which is a key contributor to heart failure. Adult female normoglycaemic mice underwent coronary artery ligation/sham surgery prior to infusion with exendin-4/vehicle for 4 weeks. Metabolic parameters and infarct sizes were comparable between groups. Exendin-4 protected against cardiac dysfunction and chamber dilatation post-MI and improved survival. Furthermore, exendin-4 modestly decreased cardiomyocyte hypertrophy/apoptosis but markedly attenuated interstitial fibrosis and myocardial inflammation post-MI. This was associated with altered extracellular matrix (procollagen IαI/IIIαI, connective tissue growth factor, fibronectin, TGF-β3) and inflammatory (IL-10, IL-1β, IL-6) gene expression in exendin-4-treated mice, together with modulation of both Akt/GSK-3β and Smad2/3 signalling. Exendin-4 also altered macrophage response gene expression in the absence of direct actions on cardiac fibroblast differentiation, suggesting cardioprotective effects occurring secondary to modulation of inflammation. Our findings indicate that exendin-4 protects against post-MI remodelling via preferential actions on inflammation and the extracellular matrix independently of its established actions on glycaemic control, thereby suggesting that selective targeting of GLP-1 signalling may be required to realise its clear therapeutic potential for post-MI heart failure.

Giving structure to the biofilm matrix: an overview of individual strategies and emerging common themes

Biofilms are communities of microbial cells that underpin diverse processes including sewage bioremediation, plant growth promotion, chronic infections and industrial biofouling. The cells resident in the biofilm are encased within a self-produced exopolymeric matrix that commonly comprises lipids, proteins that frequently exhibit amyloid-like properties, eDNA and exopolysaccharides. This matrix fulfils a variety of functions for the community, from providing structural rigidity and protection from the external environment to controlling gene regulation and nutrient adsorption. Critical to the development of novel strategies to control biofilm infections, or the capability to capitalize on the power of biofilm formation for industrial and biotechnological uses, is an in-depth knowledge of the biofilm matrix. This is with respect to the structure of the individual components, the nature of the interactions between the molecules and the three-dimensional spatial organization. We highlight recent advances in the understanding of the structural and functional role that carbohydrates and proteins play within the biofilm matrix to provide three-dimensional architectural integrity and functionality to the biofilm community. We highlight, where relevant, experimental techniques that are allowing the boundaries of our understanding of the biofilm matrix to be extended using Escherichia coli, Staphylococcus aureus, Vibrio cholerae, and Bacillus subtilis as exemplars.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.