A novel method of microsatellite genotyping-by-sequencing using individual combinatorial barcoding

This study examines the potential of next-generation sequencing based ‘genotyping-by-sequencing’ (GBS) of microsatellite loci for rapid and cost-effective genotyping in large-scale population genetic studies. The recovery of individual genotypes from large sequence pools was achieved by PCR-incorporated combinatorial barcoding using universal primers. Three experimental conditions were employed to explore the possibility of using this approach with existing and novel multiplex marker panels and weighted amplicon mixture. The GBS approach was validated against microsatellite data generated by capillary electrophoresis. GBS allows access to the underlying nucleotide sequences that can reveal homoplasy, even in large datasets and facilitates cross laboratory transfer. GBS of microsatellites, using individual combinatorial barcoding, is potentially faster and cheaper than current microsatellite approaches and offers better and more data.

Exogenous RNA interference exposes contrasting roles for sugar exudation in host-finding by plant pathogens

Plant parasitic nematodes (PPN) locate host plants by following concentration gradients of root exudate chemicals in the soil. We present a simple method for RNA interference (RNAi)-induced knockdown of genes in tomato seedling roots, facilitating the study of root exudate composition, and PPN responses. Knockdown of sugar transporter genes, STP1 and STP2, in tomato seedlings triggered corresponding reductions of glucose and fructose, but not xylose, in collected root exudate. This corresponded directly with reduced infectivity and stylet thrusting of the promiscuous PPN Meloidogyne incognita, however we observed no impact on the infectivity or stylet thrusting of the selective Solanaceae PPN Globodera pallida. This approach can underpin future efforts to understand the early stages of plant-pathogen interactions in tomato and potentially other crop plants.

LaGomiCs - Lagomorph Genomics Consortium: An International Collaborative Effort for Sequencing the Genomes of an Entire Mammalian Order

The order Lagomorpha comprises about 90 living species, divided in 2 families: the pikas (Family Ochotonidae), and the rabbits, hares, and jackrabbits (Family Leporidae). Lagomorphs are important economically and scientifically as major human food resources, valued game species, pests of agricultural significance, model laboratory animals, and key elements in food webs. A quarter of the lagomorph species are listed as threatened. They are native to all continents except Antarctica, and occur up to 5000 m above sea level, from the equator to the Arctic, spanning a wide range of environmental conditions. The order has notable taxonomic problems presenting significant difficulties for defining a species due to broad phenotypic variation, overlap of morphological characteristics, and relatively recent speciation events. At present, only the genomes of 2 species, the European rabbit (Oryctolagus cuniculus) and American pika (Ochotona princeps) have been sequenced and assembled. Starting from a paucity of genome information, the main scientific aim of the Lagomorph Genomics Consortium (LaGomiCs), born from a cooperative initiative of the European COST Action “A Collaborative European Network on Rabbit Genome Biology—RGB-Net” and the World Lagomorph Society (WLS), is to provide an international framework for the sequencing of the genome of all extant and selected extinct lagomorphs. Sequencing the genomes of an entire order will provide a large amount of information to address biological problems not only related to lagomorphs but also to all mammals. We present current and planned sequencing programs and outline the final objective of LaGomiCs possible through broad international collaboration.

Characterizing the effects of ocean acidification in larval and juvenile Manila clam, Ruditapes philippinarum, using a transcriptomic approach

Ocean acidification as a result of anthropogenic carbon dioxide (CO2) emissions and global climate change poses a risk to the ecological landscape of intertidal and shallow subtidal communities. The organisms that inhabit these waters will have to cope with changing environmental conditions through the appropriate modulation of physiological processes. Calcifying organisms are particularly at risk, as increased atmospheric levels of CO2 in the atmosphere increase the partial pressure of CO2 (pCO2) in the oceans. Increased pCO2 reduces the saturation of carbonate minerals required to form calcified structures. Being able to cope with the increased energetic demand of maintaining these structures, in addition to other vital physiological processes, will be the key driver that determines which organisms will persist. Assessment of larval and juvenile Manila clam mortality and physiology in this study suggests that this species is capable of coping with elevated pCO2 conditions. The use of high throughput sequencing and RNA sequence analysis in larval clams revealed several physiological processes that play important roles in the Manila clam’s ability to tolerate elevated pCO2 conditions during this life stage. Exposure of juvenile Manila clams, acclimated to elevated pCO2 conditions, to a thermal stress revealed that this species might also be capable of coping with multiple stressors associated with global climate change. Manila clams could therefore represent a model for studying physiological mechanisms associated with successful acclimation of populations to ocean acidification.

RT-Bst: an integrated approach for reverse transcription and enrichment of cDNA from viral RNA

The synthesis of cDNA from RNA is challenging due to the inefficiency of reverse transcription (RT). In order to address this, a method was developed known as RT-Bst for sequential RT of RNA and Bst DNA polymerase amplification for enrichment of cDNA in a single tube reaction. Using genomic RNA from bacteriophage MS2, the yield of cDNA produced by RT alone and RT-Bst were compared by analysis of PCR-amplified products. Using random primers a superior performance was observed when amplifying MS2 RNA following RT-Bst compared to RT alone, indicating that greater quantities of cDNA were present after RT-Bst. RT-Bst was also compared with RT alone for their relative ability to produce sufficient cDNA to amplify 8 target regions spanning the respiratory syncytial virus (RSV) genome. Six out of 8 targets were amplified consistently by PCR subsequent to RT-Bst amplification whereas only 3 out of 8 targets could be amplified after RT alone. RSV sequences were selectively amplified using RSV specific primers from a mixed template containing an excess of MS2 RNA in a RT-Bst reaction without amplifying MS2 sequences. This suggests that RT-Bst can be used to amplify RNA sequences non-specifically using random primers and specifically using sequence specific primers and enhances the yield of cDNA when compared to RT alone.

Next-generation sequencing is highly sensitive for the detection of beta-catenin mutations in desmoid-type fibromatoses

Desmoid-type fibromatoses are locally aggressive and frequently recurrent tumours, and an accurate diagnosis is essential for patient management. The majority of sporadic lesions harbour beta-catenin (CTNNB1) mutations. We used next-generation sequencing to detect CTNNB1 mutations and to compare the sensitivity and specificity of next-generation sequencing with currently employed mutation detection techniques: mutation-specific restriction enzyme digestion and polymerase chain reaction amplification. DNA was extracted from formalin-fixed paraffin-embedded needle biopsy or resection tissue sections from 144 patients with sporadic desmoid-type fibromatoses, four patients with syndrome-related desmoid-type fibromatoses and 11 morphological mimics. Two primer pairs were designed for CTNNB1 mutation hotspots. Using ≥10 ng of DNA, libraries were generated by Fluidigm and sequenced on the Ion Torrent Personal Genome Machine. Next-generation sequencing had a sensitivity of 92.36 % (133/144, 95 % CIs: 86.74 to 96.12 %) and a specificity of 100 % for the detection of CTNNB1 mutations in desmoid-type fibromatoses-like spindle cell lesions. All mutations detected by mutation-specific restriction enzyme digestion were identified by next-generation sequencing. Next-generation sequencing identified additional mutations in 11 tumours that were not detected by mutation-specific restriction enzyme digestion, two of which have not been previously described. Next-generation sequencing is highly sensitive for the detection of CTNNB1 mutations. This multiplex assay has the advantage of detecting additional mutations compared to those detected by mutation-specific restriction enzyme digestion (sensitivity 82.41 %). The technology requires minimal DNA and is time- and cost-efficient.

Ultra-deep sequencing provides insights into the virology of hepatitis C super-infections in a case of three sequential infections with different genotypes

The current epidemic of Hepatitis C infection in HIV-positive men who have sex with men is associated with increasing use of recreational drugs. Multiple HCV infections have been reported in haemophiliacs and intravenous drug users. Using ultra-deep sequencing analysis, we present the case of an HIV-positive MSM with evidence of three sequential HCV infections, each occurring during the acute phase of the preceding infection, following risk exposures. We observed rapid replacement of the original strain by the incoming genotype at subsequent time points. The impact of HCV super-infection remains unclear and UDS may provide new insights.

Post-Transcriptional Regulation of BCL2 mRNA by the RNA-Binding Protein ZFP36L1 in Malignant B Cells

The human ZFP36 zinc finger protein family consists of ZFP36, ZFP36L1, and ZFP36L2. These proteins regulate various cellular processes, including cell apoptosis, by binding to adenine uridine rich elements in the 3′ untranslated regions of sets of target mRNAs to promote their degradation. The pro-apoptotic and other functions of ZFP36 family members have been implicated in the pathogenesis of lymphoid malignancies. To identify candidate mRNAs that are targeted in the pro-apoptotic response by ZFP36L1, we reverse-engineered a gene regulatory network for all three ZFP36 family members using the ‘maximum information coefficient’ (MIC) for target gene inference on a large microarray gene expression dataset representing cells of diverse histological origin. Of the three inferred ZFP36L1 mRNA targets that were identified, we focussed on experimental validation of mRNA for the pro-survival protein, BCL2, as a target for ZFP36L1. RNA electrophoretic mobility shift assay experiments revealed that ZFP36L1 interacted with the BCL2 adenine uridine rich element. In murine BCL1 leukemia cells stably transduced with a ZFP36L1 ShRNA lentiviral construct, BCL2 mRNA degradation was significantly delayed compared to control lentiviral expressing cells and ZFP36L1 knockdown in different cell types (BCL1, ACHN, Ramos), resulted in increased levels of BCL2 mRNA levels compared to control cells. 3′ untranslated region luciferase reporter assays in HEK293T cells showed that wild type but not zinc finger mutant ZFP36L1 protein was able to downregulate a BCL2 construct containing the BCL2 adenine uridine rich element and removal of the adenine uridine rich core from the BCL2 3′ untranslated region in the reporter construct significantly reduced the ability of ZFP36L1 to mediate this effect. Taken together, our data are consistent with ZFP36L1 interacting with and mediating degradation of BCL2 mRNA as an important target through which ZFP36L1 mediates its pro-apoptotic effects in malignant B-cells.

Sequencing cutting patterns with colored interval graphs

The problem addressed here originates in the industry of flat glass cutting and wood panel sawing, where smaller items are cut from larger items accordingly to predefined cutting patterns. In this type of industry the smaller pieces that are cut from the patterns are piled around the machine in stacks according to the size of the pieces, which are moved to the warehouse only when all items of the same size have been cut. If the cutting machine can process only one pattern at a time, and the workspace is limited, it is desirable to set the sequence in which the cutting patterns are processed in a way to minimize the maximum number of open stacks around the machine. This problem is known in literature as the minimization of open stacks (MOSP). To find the best sequence of the cutting patterns, we propose an integer programming model, based on interval graphs, that searches for an appropriate edge completion of the given graph of the problem, while defining a suitable coloring of its vertices.

Spectrum and Frequency of GJB2 Mutations in a Cohort of 264 Portuguese Nonsyndromic Sensorineural Hearing Loss Patients

OBJECTIVE: To assess the spectrum and prevalence of mutations in the GJB2 gene in Portuguese nonsyndromic sensorineural hearing loss (NSSHL) patients. DESIGN: Sequencing of the coding region, basal promoter, exon 1, and donor splice site of the GJB2 gene; screening for the presence of the two common GJB6 deletions. STUDY SAMPLE: A cohort of 264 Portuguese NSSHL patients. RESULTS: At least one out of 21 different GJB2 variants was identified in 80 (30.2%) of the 264 patients analysed. Two mutant alleles were found in 53 (20%) of these probands, of which 83% (44/53) harboured at least one c.35delG allele. Twenty-seven (10.2%) of the probands harboured only one mutant allele. Subsequent analysis revealed that the GJB6 deletion del(GJB6-D13S1854) was present in at least 7.4% (2/27) of the patients carrying only one mutant GJB2 allele. Overall, one in five (55/264) of the patients were diagnosed as having DFNB1-related NSSHL, of which the vast majority (53/55) harboured only GJB2 mutations. CONCLUSIONS: This study provides clear demonstration that mutations in the GJB2 gene are an important cause of NSSHL in Portugal, thus representing a valuable indicator as regards therapeutical and rehabilitation options, as well as genetic counseling of these patients and their families.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Regulatory RNA binding molecules : implication for diabetes pathogenesis

Le glucose est notre principale source d'énergie. Après un repas, le taux de glucose dans le sang (glycémie) augmente, ce qui entraine la sécrétion d'insuline. L'insuline est une hormone synthétisée au niveau du pancréas par des cellules dites bêta. Elle agit sur différents organes tels que les muscles, le foie ou le tissu adipeux, induisant ainsi le stockage du glucose en vue d'une utilisation future.¦Le diabète est une maladie caractérisée par un taux élevé de glucose dans le sang (hyperglycémie), résultant d'une incapacité de notre corps à utiliser ou à produire suffisamment d'insuline. A long terme, cette hyperglycémie entraîne une détérioration du système cardio-vasculaire ainsi que de nombreuses complications. On distingue principalement deux type de diabète : le diabète de type 1 et le diabète de type 2, le plus fréquent (environ 90% des cas). Bien que ces deux maladies diffèrent sur beaucoup de points, elles partagent quelques similitudes. D'une part, on décèle une diminution de la quantité de cellules bêta. Cette diminution est cependant partielle dans le cas d'un diabète de type 2, et totale dans celui d'un diabète de type 1. D'autre part, la présence dans la circulation de médiateurs de l'inflammation nommés cytokines est décelée aussi bien chez les patients de type 1 que de type 2. Les cytokines sont sécrétées lors d'une inflammation. Elles servent de moyen de communication entre les différents acteurs de l'inflammation et ont pour certaines un effet néfaste sur la survie des cellules bêta.¦L'objectif principal de ma thèse a été d'étudier en détail l'effet de petites molécules régulatrices de l'expression génique, appelées microARNs. Basé sur le fait que de nombreuses publications ont démontré que les microARNs étaient impliqués dans différentes maladies telles que le cancer, j'ai émis l'hypothèse qu'ils pouvaient également jouer un rôle important dans le développement du diabète.¦Nous avons commencé par mettre des cellules bêta en culture en présence de cytokines, imitant ainsi un environnement inflammatoire. Nous avons pu de ce fait identifier les microARNs dont les niveaux d'expression étaient modifiés. A l'aide de méthodes biochimiques, nous avons ensuite observé que la modulation de certains microARNs par les cytokines avaient des effets néfastes sur la cellule bêta : sur sa production et sa sécrétion d'insuline, ainsi que sur sa mort (apoptose). Nous avons en conséquence pu démontrer que ces petites molécules avaient un rôle important à jouer dans le dysfonctionnement des cellules bêta induit par les cytokines, aboutissant au développement du diabète.¦-¦La cellule bêta pancréatique est une cellule endocrine présente dans les îlots de Langerhans, dans le pancréas. L'insuline, une hormone sécrétée par ces cellules, joue un rôle essentiel dans la régulation de la glycémie. Le diabète se développe si le taux d'insuline relâché par les cellules bêta n'est pas suffisant pour couvrir les besoins métaboliques corporels. Le diabète de type 1, qui représente environ 5 à 10% des cas, est une maladie auto-immune qui se caractérise par une réaction inflammatoire déclenchée par notre système immunitaire envers les cellules bêta. La conséquence de cette attaque est une disparition progressive des cellules bêta. Le diabète de type 2 est, quant à lui, largement plus répandu puisqu'il représente environ 90% des cas. Des facteurs à la fois génétiques et environnementaux sont responsables d'une diminution de la sensibilité des tissus métabolisant l'insuline, ainsi que d'une réduction de la sécrétion de l'insuline par les cellules bêta, ce qui a pour conséquence le développement de la maladie. Malgré les différences entre ces deux types de diabète, ils ont pour points communs la présence d'infiltrat immunitaire et la diminution de l'état fonctionnel des cellules bêta.¦Une meilleure compréhension des mécanismes aboutissant à l'altération de la cellule bêta est primordiale, avant de pouvoir développer de nouvelles stratégies thérapeutiques capables de guérir cette maladie. Durant ma thèse, j'ai donc étudié l'implication de petites molécules d'ARN, régulatrices de l'expression génique, appelées microARNs, dans les conditions physiopathologiques qui aboutissent au développement du diabète. J'ai débuté mon étude par l'identification de microARNs dont le niveau d'expression était modifié lorsque les cellules bêta étaient exposées à des conditions favorisant à la fois le développement du diabète de type 1 (cytokines) et celui du diabète de type 2 (palmitate). Nous avons découvert qu'une modification de l'expression des miR-21, -34a et -146a était commune aux deux traitements. Ces changements d'expressions ont également été confirmés dans deux modèles animaux : les souris NOD qui développent un diabète s'apparentant au diabète de type 1 et les souris db/db qui développent plutôt un diabète de type 2. Puis, à l'aide de puces à ADN, nous avons comparé l'expression de microARNs chez des souris NOD pré-diabétiques. Nous avons alors retrouvé des changements au niveau de l'expression des mêmes microARNs mais également au niveau d'une famille de microARNs : les miR-29a, -29b et -29c. De manière artificielle, nous avons ensuite surexprimé ou inhibé en conditions physiopathologiques l'expression de tous ces microARNs et nous nous sommes intéressés à l'impact d'un tel changement sur différentes fonctions de la cellule bêta comme la synthèse et la sécrétion d'insulinè ainsi que leur survie. Nous avons ainsi pu démontrer que les miR-21, -34a, -29a, -29b, -29c avaient un effet délétère sur la sécrétion d'insuline et que la surexpression de tous ces microARNs (excepté le miR-21) favorisait la mort. Finalement, nous avons démontré que la plupart de ces microARNs étaient impliqués dans la régulation d'importantes voies de signalisation responsables de l'apoptose des cellules bêta telles que les voies de NFKB, BCL2 ou encore JNK.¦Par conséquent, nos résultats démontrent que les microARNs ont un rôle important à jouer dans le dysfonctionnement des cellules bêta lors de la mise en place du diabète.