Evolution of genome sequencing techniques

The quality and the speed for genome sequencing has advanced at the same time that technology boundaries are stretched. This advancement has been divided so far in three generations. The first-generation methods enabled sequencing of clonal DNA populations. The second-generation massively increased throughput by parallelizing many reactions while the third-generation methods allow direct sequencing of single DNA molecules. The first techniques to sequence DNA were not developed until the mid-1970s, when two distinct sequencing methods were developed almost simultaneously, one by Alan Maxam and Walter Gilbert, and the other one by Frederick Sanger. The first one is a chemical method to cleave DNA at specific points and the second one uses ddNTPs, which synthesizes a copy from the DNA chain template. Nevertheless, both methods generate fragments of varying lengths that are further electrophoresed. Moreover, it is important to say that until the 1990s, the sequencing of DNA was relatively expensive and it was seen as a long process. Besides, using radiolabeled nucleotides also compounded the problem through safety concerns and prevented the automation. Some advancements within the first generation include the replacement of radioactive labels by fluorescent labeled ddNTPs and cycle sequencing with thermostable DNA polymerase, which allows automation and signal amplification, making the process cheaper, safer and faster. Another method is Pyrosequencing, which is based on the “sequencing by synthesis” principle. It differs from Sanger sequencing, in that it relies on the detection of pyrophosphate release on nucleotide incorporation. By the end of the last millennia, parallelization of this method started the Next Generation Sequencing (NGS) with 454 as the first of many methods that can process multiple samples, calling it the 2º generation sequencing. Here electrophoresis was completely eliminated. One of the methods that is sometimes used is SOLiD, based on sequencing by ligation of fluorescently dye-labeled di-base probes which competes to ligate to the sequencing primer. Specificity of the di-base probe is achieved by interrogating every 1st and 2nd base in each ligation reaction. The widely used Solexa/Illumina method uses modified dNTPs containing so called “reversible terminators” which blocks further polymerization. The terminator also contains a fluorescent label, which can be detected by a camera. Now, the previous step towards the third generation was in charge of Ion Torrent, who developed a technique that is based in a method of “sequencing-by-synthesis”. Its main feature is the detection of hydrogen ions that are released during base incorporation. Likewise, the third generation takes into account nanotechnology advancements for the processing of unique DNA molecules to a real time synthesis sequencing system like PacBio; and finally, the NANOPORE, projected since 1995, also uses Nano-sensors forming channels obtained from bacteria that conducts the sample to a sensor that allows the detection of each nucleotide residue in the DNA strand. The advancements in terms of technology that we have nowadays have been so quick, that it makes wonder: ¿How do we imagine the next generation?

Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era

The class Kinetoplastea encompasses both free-living and parasitic species from a wide range of hosts. Several representatives of this group are responsible for severe human diseases and for economic losses in agriculture and livestock. While this group encompasses over 30 genera, most of the available information has been derived from the vertebrate pathogenic genera Leishmania and Trypanosoma. Recent studies of the previously neglected groups of Kinetoplastea indicated that the actual diversity is much higher than previously thought. This article discusses the known segment of kinetoplastid diversity and how gene-directed Sanger sequencing and next-generation sequencing methods can help to deepen our knowledge of these interesting protists.

Iis--integrated Interactome System: A Web-based Platform For The Annotation, Analysis And Visualization Of Protein-metabolite-gene-drug Interactions By Integrating A Variety Of Data Sources And Tools.

High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

Algorithms for sequence analysis via mutagenesis

Despite many successes of conventional DNA sequencing methods, some DNAs remain difficult or impossible to sequence. Unsequenceable regions occur in the genomes of many biologically important organisms, including the human genome. Such regions range in length from tens to millions of bases, and may contain valuable information such as the sequences of important genes. The authors have recently developed a technique that renders a wide range of problematic DNAs amenable to sequencing. The technique is known as sequence analysis via mutagenesis (SAM). This paper presents a number of algorithms for analysing and interpreting data generated by this technique.

Unlocking hidden genomic sequence

Despite the success of conventional Sanger sequencing, significant regions of many genomes still present major obstacles to sequencing. Here we propose a novel approach with the potential to alleviate a wide range of sequencing difficulties. The technique involves extracting target DNA sequence from variants generated by introduction of random mutations. The introduction of mutations does not destroy original sequence information, but distributes it amongst multiple variants. Some of these variants lack problematic features of the target and are more amenable to conventional sequencing. The technique has been successfully demonstrated with mutation levels up to an average 18% base substitution and has been used to read previously intractable poly(A), AT-rich and GC-rich motifs.

The Genetic Landscape of the Childhood Cancer Medulloblastoma

Medulloblastoma (MB) is the most common malignant brain tumor of children. To identify the genetic alterations in this tumor type, we searched for copy number alterations using high-density microarrays and sequenced all known protein-coding genes and microRNA genes using Sanger sequencing in a set of 22 MBs. We found that, on average, each tumor had 11 gene alterations, fewer by a factor of 5 to 10 than in the adult solid tumors that have been sequenced to date. In addition to alterations in the Hedgehog and Wnt pathways, our analysis led to the discovery of genes not previously known to be altered in MBs. Most notably, inactivating mutations of the histone-lysine N-methyltransferase genes MLL2 or MLL3 were identified in 16% of MB patients. These results demonstrate key differences between the genetic landscapes of adult and childhood cancers, highlight dysregulation of developmental pathways as an important mechanism underlying MBs, and identify a role for a specific type of histone methylation in human tumorigenesis.

Diagnóstico molecular de cancros hereditários por sequenciação de nova geração: cancro da mama e cancro colorretal

O cancro da mama e o cancro colorretal constituem duas das principais causas de morte a nível mundial. Entre 5 a 10% destes casos estão associados a variantes germinais/hereditárias em genes de suscetibilidade para cancro. O objetivo deste trabalho consistiu em validar a utilização da sequenciação de nova geração (NGS) para identificar variantes previamente detetadas pelo método de Sanger em diversos genes de suscetibilidade para cancro da mama e colorretal. Foram sequenciadas por NGS 64 amostras de DNA de utentes com suspeita clínica de predisposição hereditária para cancro da mama ou colorretal, utilizando o painel de sequenciação TruSight Cancer e a plataforma MiSeq (Illumina). Estas amostras tinham sido previamente sequenciadas pelo método de Sanger para os genes BRCA1, BRCA2, TP53, APC, MUTYH, MLH1, MSH2 e STK11. A análise bioinformática dos resultados foi realizada com os softwares MiSeq Reporter, VariantStudio, Isaac Enrichment (Illumina) e Integrative Genomics Viewer (Broad Institute). A NGS demonstrou elevada sensibilidade e especificidade analíticas para a deteção de variantes de sequência em 8 genes de suscetibilidade para cancro colorretal e da mama, uma vez que permitiu identificar a totalidade das 412 variantes (93 únicas, incluindo 27 variantes patogénicas) previamente detetadas pelo método de Sanger. A utilização de painéis de sequenciação de genes de predisposição para cancro por NGS vem possibilitar um diagnóstico molecular mais abrangente, rápido e custo-eficiente, relativamente às metodologias convencionais.

Base molecular da hemocromatose hereditária não-clássica em Portugal

A Hemocromatose Hereditária (HH) é uma doença autossómica recessiva caracterizada pela absorção excessiva de ferro a nível intestinal e sua acumulação em órgãos vitais, podendo originar cardiomiopatia, cirrose e carcinoma hepatocelular. O correspondente diagnóstico molecular é obtido pela associação com genótipos específicos no gene HFE (homozigotia para p.Cys282Tyr ou heterozigotia composta p.Cys282Tyr/p.His63Asp). Contudo, nos países do sul da Europa, cerca de um terço dos doentes com diagnóstico clínico de HH não apresenta os referidos genótipos. Para identificar a base molecular da HH não-clássica em Portugal usaram-se metodologias de pesquisa geral de variantes genéticas (SSCP e dHPLC), Next-Generation Sequencing (NGS) e sequenciação de Sanger, cobrindo seis genes relacionados com o metabolismo do ferro em 303 doentes. Identificaram-se 69 variantes diferentes e de vários tipos, por ex. missense, nonsense, de splicing, que perturbam a transcrição do gene ou a regulação da tradução do mRNA. Seguidamente, realizaram-se estudos in silico e in vitro para esclarecer o significado etiológico de algumas das novas variantes. Concluiu-se que a base molecular desta patologia é bastante heterogénea e que a NGS é uma ferramenta adequada para efetuar a análise simultânea dos vários genes num grande número de amostras. Contudo, o estabelecimento da relevância clínica de algumas variantes requer a realização de estudos funcionais.

Estudo de dislipidemias familiares monogénicas raras

A dislipidemia é um distúrbio do perfil lipídico, seja por elevação ou diminuição de partículas lipídicas. O objetivo deste trabalho é fazer uma revisão dos casos com dislipidemia rara em estudo no Instituto Nacional de Saúde Doutor Ricardo Jorge, apresentando os dados clínicos e moleculares mais relevantes. O perfil lipídico foi determinado para cada caso índex e familiares e o estudo molecular dos genes envolvidos foi realizado por amplificação por PCR e sequenciação de Sanger. Foram estudados, ou está em curso o estudo, de 14 casos índex com os seguintes diagnósticos clínicos: Deficiência familiar em lipoproteína lípase (3), Lipodistrofia familiar parcial de Dunningan Tipo 2 (1), Deficiência em lípase ácida lisossomal (3), Abeta/hipobetalipoproteinemia (2), Deficiência em HDL (1), Hipertrigliceridemia autossómica recessiva (3), Sitosterolemia (1). O fenótipo clínico de cada caso índex é variável dependendo de cada condição. Foi encontrada a causa genética da doença em 8/14 doentes, estando os restantes ainda em estudo. Doentes com as várias dislipidemias raras apresentadas têm um risco acrescido de ter outras doenças graves como pancreatite, doença cardiovascular ou complicações neurológicas e devem, por esta razão, ser identificados o mais precocemente possível, de forma a minimizar ou prevenir os efeitos nefastos destas condições.

Development and validation of gold nanoprobes for human SNP detection towards commercial application

Conventional molecular techniques for detection and characterization of relevant nucleic acid (i.e. DNA) sequences are, nowadays, cumbersome, expensive and with reduced portability. The main objective of this dissertation consisted in the optimization and validation of a fast and low-cost colorimetric nanodiagnostic methodology for the detection of single nucleotide polymorphisms (SNPs). This was done considering SNPs associated to obesity of commercial interest for STAB VIDA, and subsequent evaluation of other clinically relevant targets. Also, integration of this methodology into a microfluidic platform envisaging portability and application on points-of-care (POC) was achieved. To warrant success in pursuing these objectives, the experimental work was divided in four sections: i) genetic association of SNPs to obesity in the Portuguese population; ii) optimization and validation of the non-cross-linking approach for complete genotype characterization of these SNPs; iii) incorporation into a microfluidic platform; and iv) translation to other relevant commercial targets. FTO dbSNP rs#:9939609 carriers had higher body mass index (BMI), total body fat mass, waist perimeter and 2.5 times higher risk to obesity. AuNPs functionalized with thiolated oligonucleotides (Au-nanoprobes) were used via the non-cross-linking to validate a diagnostics approach against the gold standard technique - Sanger Sequencing - with high levels of sensitivity (87.50%) and specificity (91.67%). A proof-of-concept POC microfluidic device was assembled towards incorporation of the molecular detection strategy. In conclusion a successful framework was developed and validated for the detection of SNPs with commercial interest for STAB VIDA, towards future translation into a POC device.

Genes associated with retinitis pigmentosa and allied diseases are frequently mutated in the general population.

BACKGROUND: Retinitis pigmentosa and other hereditary retinal degenerations (HRD) are rare genetic diseases leading to progressive blindness. Recessive HRD are caused by mutations in more than 100 different genes. Laws of population genetics predict that, on a purely theoretical ground, such a high number of genes should translate into an extremely elevated frequency of unaffected carriers of mutations. In this study we estimate the proportion of these individuals within the general population, via the analyses of data from whole-genome sequencing. METHODOLOGY/PRINCIPAL FINDINGS: We screened complete and high-quality genome sequences from 46 control individuals from various world populations for HRD mutations, using bioinformatic tools developed in-house. All mutations detected in silico were validated by Sanger sequencing. We identified clear-cut, null recessive HRD mutations in 10 out of the 46 unaffected individuals analyzed (∼22%). CONCLUSIONS/SIGNIFICANCE: Based on our data, approximately one in 4-5 individuals from the general population may be a carrier of null mutations that are responsible for HRD. This would be the highest mutation carrier frequency so far measured for a class of Mendelian disorders, especially considering that missenses and other forms of pathogenic changes were not included in our assessment. Among other things, our results indicate that the risk for a consanguineous couple of generating a child with a blinding disease is particularly high, compared to other genetic conditions.

Novel ADAM9 homozygous mutation in a consanguineous Egyptian family with severe cone-rod dystrophy and cataract.

OBJECTIVE: To genetically and phenotypically describe a new ADAM9 homozygous mutation in a consanguineous family from Egypt with autosomal recessive cone-rod dystrophy (arCRD), anterior polar and posterior subcapsular cataract. DESIGN, SETTING AND PARTICIPANTS: The parents and their six children were included. They underwent a complete ophthalmic examination with fundus photography and optical coherence tomography (OCT). INTERVENTION: DNA was extracted from peripheral blood from all family members. Screening for mutations in genes known to be implicated in retinal disorders was done with the IROme, an in-solution enrichment array, followed by high-throughput sequencing. Validation of the results was done by bidirectional Sanger sequencing of ADAM9 exon 14, including exon-intron junctions. Screening of normal controls was done by denaturing high-performance liquid chromatography. RESULTS: arCRD was diagnosed in the mother and two of her children. Bilateral anterior polar and posterior subcapsular cataract was observed in the mother and bilateral dot cataract was diagnosed in three of the four children not affected with arCRD, one of whom also had glaucoma. The characteristics of the arCRD were childhood-onset visual impairment, reorganisation of the retinal pigment epithelium with mid-periphery greyish-white discolouration, attenuated retinal vasculatur and optic disc pallor. A coloboma-like macular lesion was observed in one of the arCRD-affected children. IROme analysis identified a c.1396-2A>G homozygous mutation in the splice acceptor site of intron 13 of ADAM9. This mutation was homozygous in the two children affected by arCRD and in their affected mother. This mutation was heterozygous in the unaffected father and the four unaffected children. CONCLUSIONS AND RELEVANCE: We identified a novel autosomal recessive ADAM9 mutation causing arCRD in a consanguineous Egyptian family. The percentage of arCRD cases caused by mutation in ADAM9 remains to be determined. Few families are reported in the literature to date; hence extensive clinical descriptions of families with ADAM9 mutations are of significant importance.

Glucose transporter-1 deficiency syndrome: the expanding clinical and genetic spectrum of a treatable disorder.

Glucose transporter-1 deficiency syndrome is caused by mutations in the SLC2A1 gene in the majority of patients and results in impaired glucose transport into the brain. From 2004-2008, 132 requests for mutational analysis of the SLC2A1 gene were studied by automated Sanger sequencing and multiplex ligation-dependent probe amplification. Mutations in the SLC2A1 gene were detected in 54 patients (41%) and subsequently in three clinically affected family members. In these 57 patients we identified 49 different mutations, including six multiple exon deletions, six known mutations and 37 novel mutations (13 missense, five nonsense, 13 frame shift, four splice site and two translation initiation mutations). Clinical data were retrospectively collected from referring physicians by means of a questionnaire. Three different phenotypes were recognized: (i) the classical phenotype (84%), subdivided into early-onset (<2 years) (65%) and late-onset (18%); (ii) a non-classical phenotype, with mental retardation and movement disorder, without epilepsy (15%); and (iii) one adult case of glucose transporter-1 deficiency syndrome with minimal symptoms. Recognizing glucose transporter-1 deficiency syndrome is important, since a ketogenic diet was effective in most of the patients with epilepsy (86%) and also reduced movement disorders in 48% of the patients with a classical phenotype and 71% of the patients with a non-classical phenotype. The average delay in diagnosing classical glucose transporter-1 deficiency syndrome was 6.6 years (range 1 month-16 years). Cerebrospinal fluid glucose was below 2.5 mmol/l (range 0.9-2.4 mmol/l) in all patients and cerebrospinal fluid : blood glucose ratio was below 0.50 in all but one patient (range 0.19-0.52). Cerebrospinal fluid lactate was low to normal in all patients. Our relatively large series of 57 patients with glucose transporter-1 deficiency syndrome allowed us to identify correlations between genotype, phenotype and biochemical data. Type of mutation was related to the severity of mental retardation and the presence of complex movement disorders. Cerebrospinal fluid : blood glucose ratio was related to type of mutation and phenotype. In conclusion, a substantial number of the patients with glucose transporter-1 deficiency syndrome do not have epilepsy. Our study demonstrates that a lumbar puncture provides the diagnostic clue to glucose transporter-1 deficiency syndrome and can thereby dramatically reduce diagnostic delay to allow early start of the ketogenic diet.

Mutation screening of SEMA3A and SEMA7A in patients with congenital hypogonadotropic hypogonadism.

Background:Congenital hypogonadotropic hypogonadism (HH), a rare disorder characterized by absent, partial, or delayed puberty, can be caused by the lack or deficient number of hypothalamic gonadotropin-releasing hormone (GnRH) neurons. SEMA3A was recently implicated in the etiology of the disorder, and Sema7A-deficient mice have a reduced number of GnRH neurons in their brains.Methods:SEMA3A and SEMA7A were screened by Sanger sequencing in altogether 50 Finnish HH patients (34 with Kallmann syndrome (KS; HH with hyposmia/anosmia) and 16 with normosmic HH (nHH)). In 20 patients, mutation(s) had already been found in genes known to be implicated in congenital HH.Results:Three heterozygous variants (c.458A>G (p.Asn153Ser), c.1253A>G (p.Asn418Ser), and c.1303G>A (p.Val435Ile)) were found in SEMA3A in three KS patients, two of which also had a mutation in FGFR1. Two rare heterozygous variants (c.442C>T (p.Arg148Trp) and c.1421G>A (p.Arg474Gln)) in SEMA7A were found in one male nHH patient with a previously identified KISS1R nonsense variant and one male KS patient with a previously identified mutation in KAL1, respectively.Conclusion:Our results suggest that heterozygous missense variants in SEMA3A and SEMA7A may modify the phenotype of KS but most likely are not alone sufficient to cause the disorder.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection.

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.