Functional Genomic Strategies Research Articles

Cellular microRNA networks regulate host dependency of hepatitis C virus infection

Cellular microRNAs (miRNAs) have been shown to regulate hepatitis C virus (HCV) replication, yet a systematic interrogation of the repertoire of miRNAs impacting HCV life cycle is lacking. Here we apply integrative functional genomics strategies to elucidate global HCV–miRNA interactions. Through genome-wide miRNA mimic and hairpin inhibitor phenotypic screens, and miRNA–mRNA transcriptomics analyses, we identify three proviral and nine antiviral miRNAs that interact with HCV. These miRNAs are functionally linked to particular steps of HCV life cycle and related viral host dependencies. Further mechanistic studies demonstrate that miR-25, let-7, and miR-130 families repress essential HCV co-factors, thus restricting viral infection at multiple stages. HCV subverts the antiviral actions of these miRNAs by dampening their expression in cell culture models and HCV-infected human livers. This comprehensive HCV–miRNA interaction map provides fundamental insights into HCV-mediated pathogenesis and unveils molecular pathways linking RNA biology to viral infections.

In Vivo Regulation of the Zebrafish Endoderm Progenitor Niche by T-Box Transcription Factors

SummaryT-box transcription factors T/Brachyury homolog A (Ta) and Tbx16 are essential for correct mesoderm development in zebrafish. The downstream transcriptional networks guiding their functional activities are poorly understood. Additionally, important contributions elsewhere are likely masked due to redundancy. Here, we exploit functional genomic strategies to identify Ta and Tbx16 targets in early embryogenesis. Surprisingly, we discovered they not only activate mesodermal gene expression but also redundantly regulate key endodermal determinants, leading to substantial loss of endoderm in double mutants. To further explore the gene regulatory networks (GRNs) governing endoderm formation, we identified targets of Ta/Tbx16-regulated homeodomain transcription factor Mixl1, which is absolutely required in zebrafish for endoderm formation. Interestingly, we find many endodermal determinants coordinately regulated through common genomic occupancy by Mixl1, Eomesa, Smad2, Nanog, Mxtx2, and Pou5f3. Collectively, these findings augment the endoderm GRN and reveal a panel of target genes underlying the Ta, Tbx16, and Mixl1 mutant phenotypes.

Ascribing Functions to Genes: Journey Towards Genetic Improvement of Rice Via Functional Genomics.

Rice, one of the most important cereal crops for mankind, feeds more than half the world population. Rice has been heralded as a model cereal owing to its small genome size, amenability to easy transformation, high synteny to other cereal crops and availability of complete genome sequence. Moreover, sequence wealth in rice is getting more refined and precise due to resequencing efforts. This humungous resource of sequence data has confronted research fraternity with a herculean challenge as well as an excellent opportunity to functionally validate expressed as well as regulatory portions of the genome. This will not only help us in understanding the genetic basis of plant architecture and physiology but would also steer us towards developing improved cultivars. No single technique can achieve such a mammoth task. Functional genomics through its diverse tools viz. loss and gain of function mutants, multifarious omics strategies like transcriptomics, proteomics, metabolomics and phenomics provide us with the necessary handle. A paradigm shift in technological advances in functional genomics strategies has been instrumental in generating considerable amount of information w.r.t functionality of rice genome. We now have several databases and online resources for functionally validated genes but despite that we are far from reaching the desired milestone of functionally characterizing each and every rice gene. There is an urgent need for a common platform, for information already available in rice, and collaborative efforts between researchers in a concerted manner as well as healthy public-private partnership, for genetic improvement of rice crop better able to handle the pressures of climate change and exponentially increasing population.

Functional Genome Screening to Elucidate the Colistin Resistance Mechanism.

Antibiogram profile of 1590 clinical bacterial isolates based on thirteen different antimicrobial compounds showed that 1.6% of the bacterial isolates are multidrug resistant. Distribution pattern based on 16S rRNA sequence analysis showed that Pseudomonas aeruginosa constituted the largest group (83.6%) followed by Burkholderia pseudomallei sp. A191 (5.17%), Staphylococcus sp. A261 (3.45%). Among the various antibiotics used, colistin appeared to be the most effective against the Gram negative bacteria. Burkholderia pseudomallei sp. A191 and Pseudomonas aeruginosa sp. A111 showed resistance to 1500 μg/ml and 750 μg/ml of colistin respectively which constitutes 7.7% of the bacterial population. A functional genomics strategy was employed to discover the molecular support for colistin resistance in Burkholderia pseudomallei sp. A191. A pUC plasmid-based genomic expression library was constructed with an estimated library size of 2.1 × 107bp. Five colistin resistant clones were obtained after functional screening of the library. Analysis of DNA sequence of five colistin resistant clones showed homology to two component regularity systems (TCRS) encoding for a histidine kinase (mrgS) and its regulatory component (mrgR). Cross complementation assay showed that mutations in mrgS were sufficient enough to confer colistin resistant phenotype in a sensitive strain.

Drought Stress in Grain Legumes during Reproduction and Grain Filling

AbstractWater scarcity is a major constraint limiting grain legume production particularly in the arid and semi‐arid tropics. Different climate models have predicted changes in rainfall distribution and frequent drought spells for the future. Although drought impedes the productivity of grain legumes at all growth stages, its occurrence during reproductive and grain development stages (terminal drought) is more critical and usually results in significant loss in grain yield. However, the extent of yield loss depends on the duration and intensity of the stress. A reduction in the rate of net photosynthesis, and poor grain set and grain development are the principal reasons for terminal drought‐induced loss in grain yield. Insight into the impact and resistance mechanism of terminal drought is required for effective crop improvement programmes aiming to improve resistance to terminal drought in grain legumes. In this article, the impact of terminal drought on leaf development and senescence, light harvesting and carbon fixation, and grain development and grain composition is discussed. The mechanisms of resistance, management options, and innovative breeding and functional genomics strategies to improve resistance to terminal drought in grain legumes are also discussed.

Functional Genomic Strategies for Elucidating Human-Virus Interactions: Will CRISPR Knockout RNAi and Haploid Cells?

Over the last several years a wealth of transformative human–virus interaction discoveries have been produced using loss-of-function functional genomics. These insights have greatly expanded our understanding of how human pathogenic viruses exploit our cells to replicate. Two technologies have been at the forefront of this genetic revolution, RNA interference (RNAi) and random retroviral insertional mutagenesis using haploid cell lines (haploid cell screening), with the former technology largely predominating. Now the cutting edge gene editing of the CRISPR/Cas9 system has also been harnessed for large-scale functional genomics and is poised to possibly displace these earlier methods. Here we compare and contrast these three screening approaches for elucidating host–virus interactions, outline their key strengths and weaknesses including a comparison of an arrayed multiple orthologous RNAi reagent screen to a pooled CRISPR/Cas9 human rhinovirus 14–human cell interaction screen, and recount some notable insights made possible by each. We conclude with a brief perspective on what might lie ahead for the fast evolving field of human–virus functional genomics.

HDAC6 activity is a non-oncogene addiction hub for inflammatory breast cancers.

IntroductionInflammatory breast cancer (IBC) is the most lethal form of breast cancers with a 5-year survival rate of only 40 %. Despite its lethality, IBC remains poorly understood which has greatly limited its therapeutic management. We thus decided to utilize an integrative functional genomic strategy to identify the Achilles’ heel of IBC cells.MethodsWe have pioneered the development of genetic tools as well as experimental and analytical strategies to perform RNAi-based loss-of-function studies at a genome-wide level. Importantly, we and others have demonstrated that these functional screens are able to identify essential functions linked to certain cancer phenotypes. Thus, we decided to use this approach to identify IBC specific sensitivities.ResultsWe identified and validated HDAC6 as a functionally necessary gene to maintain IBC cell viability, while being non-essential for other breast cancer subtypes. Importantly, small molecule inhibitors for HDAC6 already exist and are in clinical trials for other tumor types. We thus demonstrated that Ricolinostat (ACY1215), a leading HDAC6 inhibitor, efficiently controls IBC cell proliferation both in vitro and in vivo. Critically, functional HDAC6 dependency is not associated with genomic alterations at its locus and thus represents a non-oncogene addiction. Despite HDAC6 not being overexpressed, we found that its activity is significantly higher in IBC compared to non-IBC cells, suggesting a possible rationale supporting the observed dependency.ConclusionOur finding that IBC cells are sensitive to HDAC6 inhibition provides a foundation to rapidly develop novel, efficient, and well-tolerated targeted therapy strategies for IBC patients.Electronic supplementary materialThe online version of this article (doi:10.1186/s13058-015-0658-0) contains supplementary material, which is available to authorized users.

Functional genomics to discover antibiotic resistance genes: The paradigm of resistance to colistin mediated by ethanolamine phosphotransferase in Shewanella algae MARS 14

Shewanella algae MARS 14 is a colistin-resistant clinical isolate retrieved from bronchoalveolar lavage of a hospitalised patient. A functional genomics strategy was employed to discover the molecular support for colistin resistance in S. algae MARS 14. A pZE21 MCS-1 plasmid-based genomic expression library was constructed in Escherichia coli TOP10. The estimated library size was 1.30×108bp. Functional screening of colistin-resistant clones was carried out on Luria–Bertani agar containing 8mg/L colistin. Five colistin-resistant clones were obtained after complete screening of the genomic expression library. Analysis of DNA sequencing results found a unique gene in all selected clones. Amino acid sequence analysis of this unique gene using the Integrated Microbial Genomes (IMG) and KEGG databases revealed that this gene encodes ethanolamine phosphotransferase (EptA, or so-called PmrC). Reverse transcription PCR analysis indicated that resistance to colistin in S. algae MARS 14 was associated with overexpression of EptA (27-fold increase), which plays a crucial role in the arrangement of outer membrane lipopolysaccharide.

A Proteomic Analysis of the Upper and Lower Flanks of the Base of Rice Shoot in the Gravitropism

Due to gravitational stimulation, the lower part of a shoot base grows faster than the upper part, leading the shoot to curve upward. Though much research has been done on the mechanism of plant gravitropism, it still requires extensive elucidation. Recently, functional genomic strategies have been applied to study this mechanism in plants. The present study carried out a proteomic analysis to gain a better understanding of gravity stimulation in rice. Three-week-old rice seedlings were gravitropically stimulated and samples were harvested at 4 different time points: 0.5, 3, 6, and 9 h. Then, the total crude proteins were extracted from the lower and upper parts of the shoot base, separated by 2-DE, and silver stained. At each time point, proteins in the lower and upper parts were compared, and the differently expressed proteins were identified using MALDI TOF or ESI-MS/MS. After gravity stimulation, proteins involved in nine different functional categories were either up-regulated or down-regulated. Sugar metabolism, glycolysis, the tricarboxylic acid (TCA/citric) cycle, pyruvate metabolism, and transcription regulation-related proteins were regulated. Although the initiation of defense reactions mainly occurred in roots, some different defense mechanisms were also evoked in the aerial tissues. Interestingly, the abundance of some proteins changed drastically at only 0.5 h after reorientation: inosine monophosphate dehydrogenase (up to 6.49-fold higher in lower flanks at 0.5 h), ATP synthase D (4.25-fold), and ribulose-1,5 -bisphosphate carboxylase oxygenase (3.62-fold). These findings may aid in understanding the mechanism of the gravitropism.

Rare platelet GPCR variants: what can we learn?

Platelet-expressed GPCRs are critical regulators of platelet function. Pharmacological blockade of these receptors forms a powerful therapeutic tool in the treatment and prevention of arterial thrombosis associated with coronary atherosclerosis and ischaemic stroke. However, anti-thrombotic drug therapy is associated with high inter-patient variability in therapeutic response and adverse bleeding side effects. In order to optimize the use of existing anti-platelet drugs and to develop new therapies, more detailed knowledge is required relating to the molecular mechanisms that regulate GPCR and therefore platelet function. One approach has been to identify rare, function-disrupting mutations within key platelet proteins in patients with bleeding disorders. In this review, we describe how an integrated functional genomics strategy has contributed important structure-function information about platelet GPCRs with specific emphasis upon purinergic and thromboxane A2 receptors. We also discuss the potential implications these findings have for pharmacotherapy and for understanding the molecular basis of mild bleeding disorders.

Mapping the Extracellular and Membrane Proteome Associated with the Vasculature and the Stroma in the Embryo

In order to map the extracellular or membrane proteome associated with the vasculature and the stroma in an embryonic organism in vivo, we developed a biotinylation technique for chicken embryo and combined it with mass spectrometry and bioinformatic analysis. We also applied this procedure to implanted tumors growing on the chorioallantoic membrane or after the induction of granulation tissue. Membrane and extracellular matrix proteins were the most abundant components identified. Relative quantitative analysis revealed differential protein expression patterns in several tissues. Through a bioinformatic approach, we determined endothelial cell protein expression signatures, which allowed us to identify several proteins not yet reported to be associated with endothelial cells or the vasculature. This is the first study reported so far that applies in vivo biotinylation, in combination with robust label-free quantitative proteomics approaches and bioinformatic analysis, to an embryonic organism. It also provides the first description of the vascular and matrix proteome of the embryo that might constitute the starting point for further developments.

From promises to practical strategies in epigenetic epidemiology

The epigenome has been heralded as a key 'missing piece' of the aetiological puzzle for complex phenotypes across the biomedical sciences. The standard research approaches developed for genetic epidemiology, however, are not necessarily appropriate for epigenetic studies of common disease. Here, we discuss the optimal execution of population-based studies of epigenetic variation, which will contribute to the emerging field of 'epigenetic epidemiology' and emphasize the importance of establishing a causal role in pathology for disease-associated epigenetic changes. We propose that improved understanding of the molecular mechanisms underlying human health and disease are best achieved through carrying out studies of epigenetics in populations as a part of an integrated functional genomics strategy.

Progress and pitfalls in the functional analysis of OA susceptibility alleles

Better knowledge of variation in the human genome and the availability of novel technologies allowing the simultaneous assessment of millions of single nucleotide polymorphisms (SNP) using DNA microarrays, have dramatically changed the landscape of genetic research. Genome wide association studies (GWAS) for complex diseases have generated a wave of enthusiasm in the community and lead to substantial progress in understanding many chronic diseases. However, analyzing GWAS studies in osteoarthritis is still a major challenge. Large patient and control cohorts are necessary and phenotype definitions, joints of interest as well as the chronic progressive nature of the diseases have raised enormous challenges for genetic researchers in this field. Nevertheless, different genetic regions have been associated with OA. Identification of regions or even of specific SNPs does not directly demonstrate the involvement of a given gene. Therefore, GWAS studies only provide a starting point to better understand the pathophysiology of disease and the identification of therapeutic targets and strategies. The functional genomics approaches includes different in vitro and in vivo models. The challenges following GWAS are clearly illustrated by the identification of an OA locus on Chromosome 9 that harbors different genes which all could involved in OA. Current evidence suggest that GPR22 and COG5 could play a role in OA but other effects from the C9 cluster cannot be excluded. Moreover, associated regions do not only harbor genes but also other regulators of biological processes such as miRNAs. GWAS studies also will only identify common variants often with limited impact on the heritability of disease, whereas rare variants in genes may have a more profound effect on disease development in an affected individual. Therefore, novel technologies such as next generation sequencing likely will provide additional tools to understand processes that define OA. Again, functional genomic strategies identifying the roles of specific genes and their variants in joint and joint-associated studies will be necessary to provide a full picture.

Discovering Thiamine Transporters as Targets of Chloroquine Using a Novel Functional Genomics Strategy

Chloroquine (CQ) and other quinoline-containing antimalarials are important drugs with many therapeutic benefits as well as adverse effects. However, the molecular targets underlying most such effects are largely unknown. By taking a novel functional genomics strategy, which employs a unique combination of genome-wide drug-gene synthetic lethality (DGSL), gene-gene synthetic lethality (GGSL), and dosage suppression (DS) screens in the model organism Saccharomyces cerevisiae and is thus termed SL/DS for simplicity, we found that CQ inhibits the thiamine transporters Thi7, Nrt1, and Thi72 in yeast. We first discovered a thi3Δ mutant as hypersensitive to CQ using a genome-wide DGSL analysis. Using genome-wide GGSL and DS screens, we then found that a thi7Δ mutation confers severe growth defect in the thi3Δ mutant and that THI7 overexpression suppresses CQ-hypersensitivity of this mutant. We subsequently showed that CQ inhibits the functions of Thi7 and its homologues Nrt1 and Thi72. In particular, the transporter activity of wild-type Thi7 but not a CQ-resistant mutant (Thi7T287N) was completely inhibited by the drug. Similar effects were also observed with other quinoline-containing antimalarials. In addition, CQ completely inhibited a human thiamine transporter (SLC19A3) expressed in yeast and significantly inhibited thiamine uptake in cultured human cell lines. Therefore, inhibition of thiamine uptake is a conserved mechanism of action of CQ. This study also demonstrated SL/DS as a uniquely effective methodology for discovering drug targets.

Osteoporosis genetics: year 2011 in review

Increased rates of osteoporotic fractures represent a worldwide phenomenon, which result from a progressing aging in the population around the world and creating socioeconomic problems. This review will focus mostly on human genetic studies identifying genomic regions, genes and mutations associated with osteoporosis (bone mineral density (BMD) and bone loss) and related fractures, which were published during 2011. Although multiple genome-wide association studies (GWAS) were performed to date, the genetic cause of osteoporosis and fractures has not yet been found, and only a small fraction of high heritability of bone mass was successfully explained. GWAS is a successful tool to initially define and prioritize specific chromosomal regions showing associations with the desired traits or diseases. Following the initial discovery and replication, targeted sequencing is needed in order to detect those rare variants which GWAS does not reveal by design. Recent GWAS findings for BMD included WNT16 and MEF2C. The role of bone morphogenetic proteins in fracture healing has been explored by several groups, and new single-nucleotide polymorphisms present in genes such as NOGGIN and SMAD6 were found to be associated with a greater risk of fracture non-union. Finding new candidate genes, and mutations associated with BMD and fractures, also provided new biological connections. Thus, candidates for molecular link between bone metabolism and lactation (for example, RAP1A gene), as well as possible pleiotropic effects for bone and muscle (ACTN3 gene) were suggested. The focus of contemporary studies seems to move toward whole-genome sequencing, epigenetic and functional genomics strategies to find causal variants for osteoporosis.

Integrated approaches to functionally characterize novel factors in lipoprotein metabolism

To discuss if and how the combined analysis of large-scale datasets from multiple independent sources benefits the mapping of novel genetic elements with relevance to lipoprotein metabolism and allows for conclusions on underlying molecular mechanisms. Genome-wide association studies (GWAS) have identified numerous genomic loci associated with plasma lipid levels and cardiovascular disease. Yet, despite being highly successful in mapping novel loci the GWAS approach falls short to systematically extract functional information from genomic data. With the aim to complement GWAS for a better insight into disease mechanisms and identification of the most promising targets for drug development, a number of high-throughput functional genomics strategies have now been applied. These include computational approaches, consideration of gene-gene and gene-environment interactions, as well as unbiased gene-expression analyses in relevant tissues. For a limited number of loci, mechanistic insight has been gained through in-vitro and in-vivo studies by knockdown and overexpression of candidate genes. The integration of GWAS data with existing functional genomics strategies has contributed to ascertain the relevance of a number of novel factors for lipoprotein biology and disease. However, technologies are warranted that provide a more systematic insight into the molecular function and pathogenic relevance of promising candidate genes.

Exploring Sterol and Lipid Biology in the Functional Genomic Era

The growing prevalence of disorders of lipid and sterol homeostasis such as obesity and atherosclerosis in Western societies has increased demand for the development of therapies for their prevention and treatment. Crucial for this development is an understanding of the underlying cellular mechanisms which are involved in both healthy and diseased states. However, gaps remain in our knowledge of fundamental processes, such as intracellular sterol transport, lipid storage and mobilisation. Functional genomic strategies, such as genome-wide gene knockdown screens, are increasingly being exploited as powerful strategies to fill these gaps. This review discusses experimental approaches and findings in recent functional genomics studies of sterol and lipid biology, both in model organisms and human cells.

Foxp2 regulates gene networks implicated in neurite outgrowth in the developing brain.

Forkhead-box protein P2 is a transcription factor that has been associated with intriguing aspects of cognitive function in humans, non-human mammals, and song-learning birds. Heterozygous mutations of the human FOXP2 gene cause a monogenic speech and language disorder. Reduced functional dosage of the mouse version (Foxp2) causes deficient cortico-striatal synaptic plasticity and impairs motor-skill learning. Moreover, the songbird orthologue appears critically important for vocal learning. Across diverse vertebrate species, this well-conserved transcription factor is highly expressed in the developing and adult central nervous system. Very little is known about the mechanisms regulated by Foxp2 during brain development. We used an integrated functional genomics strategy to robustly define Foxp2-dependent pathways, both direct and indirect targets, in the embryonic brain. Specifically, we performed genome-wide in vivo ChIP–chip screens for Foxp2-binding and thereby identified a set of 264 high-confidence neural targets under strict, empirically derived significance thresholds. The findings, coupled to expression profiling and in situ hybridization of brain tissue from wild-type and mutant mouse embryos, strongly highlighted gene networks linked to neurite development. We followed up our genomics data with functional experiments, showing that Foxp2 impacts on neurite outgrowth in primary neurons and in neuronal cell models. Our data indicate that Foxp2 modulates neuronal network formation, by directly and indirectly regulating mRNAs involved in the development and plasticity of neuronal connections.

A genomic-based approach identifies FXYD domain containing ion transport regulator 2 (FXYD2)γa as a pancreatic beta cell-specific biomarker

Non-invasive imaging of the pancreatic beta cell mass (BCM) requires the identification of novel and specific beta cell biomarkers. We have developed a systems biology approach to the identification of promising beta cell markers. We followed a functional genomics strategy based on massive parallel signal sequencing (MPSS) and microarray data obtained in human islets, purified primary rat beta cells, non-beta cells and INS-1E cells to identify promising beta cell markers. Candidate biomarkers were validated and screened using established human and macaque (Macacus cynomolgus) tissue microarrays. After a series of filtering steps, 12 beta cell-specific membrane proteins were identified. For four of the proteins we selected or produced antibodies targeting specifically the human proteins and their splice variants; all four candidates were confirmed as islet-specific in human pancreas. Two splice variants of FXYD domain containing ion transport regulator 2 (FXYD2), a regulating subunit of the Na(+)-K(+)-ATPase, were identified as preferentially present in human pancreatic islets. The presence of FXYD2gammaa was restricted to pancreatic islets and selectively detected in pancreatic beta cells. Analysis of human fetal pancreas samples showed the presence of FXYD2gammaa at an early stage (15 weeks). Histological examination of pancreatic sections from individuals with type 1 diabetes or sections from pancreases of streptozotocin-treated Macacus cynomolgus monkeys indicated a close correlation between loss of FXYD2gammaa and loss of insulin-positive cells. We propose human FXYD2gammaa as a novel beta cell-specific biomarker.

Gata3 regulates trophoblast development downstream of Tead4 and in parallel to Cdx2

The mouse blastocyst and stem cells derived from its tissue lineages provide a unique genetic system for examining the establishment and loss of pluripotency. The transcription factor Cdx2 plays a central role by repressing pluripotency genes, such as Oct4, and promoting extraembryonic trophoblast fate at the blastocyst stage. However, genetic evidence has suggested that Cdx2 does not work alone in the trophoblast lineage. We have used bioinformatic and functional genomic strategies to identify the transcription factor Gata3 as a trophoblast factor. We show Gata3 to be capable of inducing trophoblast fate in embryonic stem cells and driving trophoblast differentiation in trophoblast stem cells. In addition, Cdx2 is not required for Gata3-induced expression of a subset of trophoblast genes in embryonic stem cells. We show that Gata3 is coexpressed with Cdx2 in the blastocyst, but this does not depend on Cdx2. In the embryo, expression of Gata3, like that of Cdx2, depends on Tead4, and the expression of both factors becomes restricted to trophoblast by a mechanism that does not initially rely on Oct4. These observations suggest that Gata3 and Cdx2 can act in parallel pathways downstream of Tead4 to induce the expression of common and independent targets in the trophoblast lineage, whereas Oct4 is required for continued repression of trophoblast fate in the embryonic lineage.