Towards accurate indel calling for oncopanel sequencing through an international pipeline competition at precisionFDA

Accurately calling indels with next-generation sequencing (NGS) data is critical for clinical application. The precisionFDA team collaborated with the U.S. Food and Drug Administration's (FDA's) National Center for Toxicological Research (NCTR) and successfully completed the NCTR Indel Calling from Oncopanel Sequencing Data Challenge, to evaluate the performance of indel calling pipelines. Top performers were selected based on precision, recall, and F1-score. The performance of many other pipelines was close to the top performers, which produced a top cluster of performers. The performance was significantly higher in high confidence regions and coding regions, and significantly lower in low complexity regions. Oncopanel capture and other issues may have occurred that affected the recall rate. Indels with higher variant allele frequency (VAF) may generally be called with higher confidence. Many of the indel calling pipelines had good performance. Some of them performed generally well across all three oncopanels, while others were better for a specific oncopanel. The performance of indel calling can further be improved by restricting the calls within high confidence intervals (HCIs) and coding regions, and by excluding low complexity regions (LCR) regions. Certain VAF cut-offs could be applied according to the applications.

Genomic epidemiology of the SARS-CoV-2 epidemic in Brazil

The high numbers of COVID-19 cases and deaths in Brazil have made Latin America an epicentre of the pandemic. SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, but important gaps remain in our understanding of virus transmission dynamics at a national scale. We use 17,135 near-complete genomes sampled from 27 Brazilian states and bordering country Paraguay. From March to November 2020, we detected co-circulation of multiple viral lineages that were linked to multiple importations (predominantly from Europe). After November 2020, we detected large, local transmission clusters within the country. In the absence of effective restriction measures, the epidemic progressed, and in January 2021 there was emergence and onward spread, both within and abroad, of variants of concern and variants under monitoring, including Gamma (P.1) and Zeta (P.2). We also characterized a genomic overview of the epidemic in Paraguay and detected evidence of importation of SARS-CoV-2 ancestor lineages and variants of concern from Brazil. Our findings show that genomic surveillance in Brazil enabled assessment of the real-time spread of emerging SARS-CoV-2 variants.

Genomic epidemiology reveals the impact of national and international restrictions measures on the SARS-CoV-2 epidemic in Brazil

Brazil has experienced some of the highest numbers of COVID-19 cases and deaths globally and from May 2021 made Latin America a pandemic epicenter. Although SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, important gaps remain in our understanding of virus transmission dynamics at the national scale. Here, we describe the genomic epidemiology of SARS-CoV-2 using near-full genomes sampled from 27 Brazilian states and a bordering country - Paraguay. We show that the early stage of the pandemic in Brazil was characterised by the co-circulation of multiple viral lineages, linked to multiple importations predominantly from Europe, and subsequently characterized by large local transmission clusters. As the epidemic progressed under an absence of effective restriction measures, there was a local emergence and onward international spread of Variants of Concern (VOC) and Variants Under Monitoring (VUM), including Gamma (P.1) and Zeta (P.2). In addition, we provide a preliminary genomic overview of the epidemic in Paraguay, showing evidence of importation from Brazil. These data reinforce the usefulness and need for the implementation of widespread genomic surveillance in South America as a toolkit for pandemic monitoring that provides a means to follow the real-time spread of emerging SARS-CoV-2 variants with possible implications for public health and immunization strategies.

Genomic monitoring of the SARS-CoV-2 B1.1.7 (WHO VOC Alpha) in the Sao Paulo state, Brazil

The SARS-CoV-2 alpha VOC (also known as lineage B.1.1.7) initially described in the autumn, 2020 in UK, rapidly became the dominant lineage across much of Europe. Despite multiple studies reporting molecular evidence suggestive of its circulation in Brazil, much is still unknown about its genomic diversity in the state of São Paulo, the main Brazilian economic and transportation hub. To get more insight regarding its transmission dynamics into the State we performed phylogenetic analysis on all alpha VOC strains obtained between February and August 2021 from the Sao Paulo state Network for Pandemic Alert of Emerging SARS-CoV-2 variants. The performed phylogenetic analysis showed that most of the alpha VOC genomes were interspersed with viral strains sampled from different Brazilian states and other countries suggesting that multiple independent Alpha VOC introductions from Brazil and overseas have occurred in the São Paulo State over time. Nevertheless, large monophyletic clusters were also observed especially from the Central-West part of the São Paulo State (the city of Bauru) and the metropolitan region of the São Paulo city. Our results highlight the Alpha VOC molecular epidemiology in the São Paulo state and reinforce the need for continued genomic surveillance strategies for the real-time monitoring of potential emerging SARS-CoV-2 variants during the ever-growing vaccination process.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

A Novel Multisystem Proteinopathy Caused by a Missense ANXA11 Variant

OBJECTIVE

Protein misfolding plays a central role not only in amyotrophic lateral sclerosis (ALS), but also in other conditions, such as frontotemporal dementia (FTD), inclusion body myopathy (hIBM) or Paget's disease of bone. The concept of multisystem proteinopathies (MSP) was created to account for those rare families that segregate at least 2 out of these 4 conditions in the same pedigree. The calcium-dependent phospholipid-binding protein annexin A11 was recently associated to ALS in European pedigrees. Herein, we describe in detail 3 Brazilian families presenting hIBM (isolated or in combination with ALS/FTD) caused by the novel p.D40Y change in the gene encoding annexin A11 (ANXA11).

METHODS

We collected clinical, genetic, pathological and skeletal muscle imaging from 11 affected subjects. Neuroimaging was also obtained from 8 patients and 8 matched controls.

RESULTS

Clinico-radiological phenotype of this novel hIBM reveals a slowly progressive predominant limb-girdle syndrome, but with frequent axial (ptosis/dropped head) and distal (medial gastrocnemius) involvement as well. Muscle pathology identified numerous rimmed vacuoles with positive annexin A11, TDP-43 and p62 inclusions, but no inflammation. Central nervous system was also involved: two patients had FTD, but diffusion tensor imaging uncovered multiple areas of cerebral white matter damage in the whole group (including the corticospinal tracts and frontal subcortical regions).

INTERPRETATION

These findings expand the phenotypic spectrum related to ANXA11. This gene should be considered the cause of a novel multisystem proteinopathy (MSP type 6), rather than just ALS. ANN NEUROL 2021;90:239-252.

Complete Chromosome-Scale Genome Sequence Resource for Sporisorium panici-leucophaei, the Causal Agent of Sourgrass Smut Disease

Here, we present the first complete chromosome-level genome assembly of the smut fungus strain Sporisorium panici-leucophaei SPL10A, the causal agent of the sourgrass (Digitaria insularis) smut disease. Combining Illumina paired-end and Nanopore long reads, we generated a final assembly composed of 23 chromosomes (22 nuclear and one mitochondrial) with 18,915,934 bp. Gene prediction accomplished using extrinsic evidence from the sugarcane smut fungus Sporisorium scitamineum originated a total of 6,402 protein-encoding genes. The secretome (388 proteins) and the effectorome repertoires (68 candidates) were also predicted, given their crucial roles in plant-pathogen interactions. The complete telomere-to-telomere chromosome sequences of this poorly studied fungus will provide a valuable resource for future comparative genomic studies among smuts to unravel their underlying pathogenicity mechanisms.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

High-Quality Draft Genome Sequence Resources of Eight Xylella fastidiosa Strains Isolated from Citrus, Coffee, Plum, and Hibiscus in South America

Xylella fastidiosa subsp. pauca, once confined to South America and infecting mainly citrus and coffee plants, has been found to be associated with other hosts and in other geographic regions. We present high-quality draft genome sequences of X. fastidiosa subsp. pauca strains J1a12, B111, U24D, and XRB isolated from citrus plants in Brazil, strain Fb7 isolated from a citrus plant in Argentina and strains 3124, Pr8x, and Hib4 isolated, respectively, from coffee, plum, and hibiscus plants in Brazil. Sequencing was performed using Roche 454-GS FLX, MiSeq-Illumina or Pacific Biosciences platforms. These high-quality genome assemblies will be useful for further studies about the genomic diversity, evolution, and biology of X. fastidiosa.

Additional observation of a de novo pathogenic variant in KCNT2 leading to epileptic encephalopathy with clinical features of frontal lobe epilepsy

INTRODUCTION

KCNT2 was recently recognized as a gene associated with neurodevelopmental disorder and epilepsy.

CASE REPORT

We present an additional observation of a 16-year-old male patient with a novel de novo KCNT2 likely pathogenic variant and review the five previously reported cases of de novo variants in this gene.

DISCUSSION

Whole exome sequencing identified the missense variant c.725C > A p.(Thr242Asn), which was confirmed by Sanger sequencing. Our patient has a refractory stereotyped and monomorphic type of hyperkinetic focal motor seizure, similar to what is seen in frontal lobe epilepsy, occurring only during sleep. This type of seizure is not usually seen in epileptic encephalopathies.

Intragenic variants in the SMN1 gene determine the clinical phenotype in 5q spinal muscular atrophy

Objective

The aim of the study was to report the proportion of homozygous and compound heterozygous variants in the survival motor neuron 1 (SMN1) gene in a large population of patients with spinal muscular atrophy (SMA) and to correlate the severity of the disease with the presence of specific intragenic variants in SMN1 and with the SMN2 copy number.

Methods

Four hundred fifty Brazilian patients with SMA were included in a retrospective study, and clinical data were analyzed compared with genetic data; the SMN2 copy number was obtained by multiplex ligation-dependent probe amplification and pathogenic variants in SMN1 by next-generation sequencing.

Results

Four hundred two patients (89.3%) presented homozygous exon 7-SMN1 deletion, and 48 (10.7%) were compound heterozygous for the common deletion in one allele and a point mutation in the other allele. Recurrent variants in exons 3 and 6 (c.460C>T, c.770_780dup and c.734_735insC) accounted for almost 80% of compound heterozygous patients. Another recurrent pathogenic variant was c.5C>G at exon 1. Patients with c.770_780dup and c.734_735insC had a clinical phenotype correlated with SMN2 copy number, whereas the variants c.460C>T and c.5C>G determined a milder phenotype independently of the SMN2 copies.

Conclusions

Patients with specific pathogenic variants (c.460C>T and c.5C>G) presented a milder phenotype, and the SMN2 copy number did not correlate with disease severity in this group.

Neurodevelopmental disorder associated with de novo SCN3A pathogenic variants: two new cases and review of the literature

SCN3A was recently recognized as a gene associated with neurodevelopmental disorder and epilepsy. We present two additional patients with a novel de novo SCN3A pathogenic variant, and a review of all published cases of de novo variants. In one of our patients brain magnetic resonance imaging (MRI) disclosed a severe polymicrogyria and in the other it was normal. The clinical phenotype was characterized by a severe developmental delay and refractory epilepsy in the patient with polymicrogyria and intellectual disability with autistic features and pharmacoresponsive epilepsy in the subject with normal MRI. Polymicrogyria, a disorder of progenitor cells proliferation and migration, is an unanticipated finding for an ion channel dysfunction.

Analysis of Three Sugarcane Homo/Homeologous Regions Suggests Independent Polyploidization Events of Saccharum officinarum and Saccharum spontaneum

Whole genome duplication has played an important role in plant evolution and diversification. Sugarcane is an important crop with a complex hybrid polyploid genome, for which the process of adaptation to polyploidy is still poorly understood. In order to improve our knowledge about sugarcane genome evolution and the homo/homeologous gene expression balance, we sequenced and analyzed 27 BACs (Bacterial Artificial Chromosome) of sugarcane R570 cultivar, containing the putative single-copy genes LFY (seven haplotypes), PHYC (four haplotypes), and TOR (seven haplotypes). Comparative genomic approaches showed that these sugarcane loci presented a high degree of conservation of gene content and collinearity (synteny) with sorghum and rice orthologous regions, but were invaded by transposable elements (TE). All the homo/homeologous haplotypes of LFY, PHYC, and TOR are likely to be functional, because they are all under purifying selection (dN/dS ≪ 1). However, they were found to participate in a nonequivalently manner to the overall expression of the corresponding gene. SNPs, indels, and amino acid substitutions allowed inferring the S. officinarum or S. spontaneum origin of the TOR haplotypes, which further led to the estimation that these two sugarcane ancestral species diverged between 2.5 and 3.5 Ma. In addition, analysis of shared TE insertions in TOR haplotypes suggested that two autopolyploidization may have occurred in the lineage that gave rise to S. officinarum, after its divergence from S. spontaneum.

Integrative Variation Analysis Reveals that a Complex Genotype May Specify Phenotype in Siblings with Syndromic Autism Spectrum Disorder

It has been proposed that copy number variations (CNVs) are associated with increased risk of autism spectrum disorder (ASD) and, in conjunction with other genetic changes, contribute to the heterogeneity of ASD phenotypes. Array comparative genomic hybridization (aCGH) and exome sequencing, together with systems genetics and network analyses, are being used as tools for the study of complex disorders of unknown etiology, especially those characterized by significant genetic and phenotypic heterogeneity. Therefore, to characterize the complex genotype-phenotype relationship, we performed aCGH and sequenced the exomes of two affected siblings with ASD symptoms, dysmorphic features, and intellectual disability, searching for de novo CNVs, as well as for de novo and rare inherited point variations—single nucleotide variants (SNVs) or small insertions and deletions (indels)—with probable functional impacts. With aCGH, we identified, in both siblings, a duplication in the 4p16.3 region and a deletion at 8p23.3, inherited by a paternal balanced translocation, t(4, 8) (p16; p23). Exome variant analysis found a total of 316 variants, of which 102 were shared by both siblings, 128 were in the male sibling exome data, and 86 were in the female exome data. Our integrative network analysis showed that the siblings’ shared translocation could explain their similar syndromic phenotype, including overgrowth, macrocephaly, and intellectual disability. However, exome data aggregate genes to those already connected from their translocation, which are important to the robustness of the network and contribute to the understanding of the broader spectrum of psychiatric symptoms. This study shows the importance of using an integrative approach to explore genotype-phenotype variability.

The complete mitochondrial genome of Bothrops jararaca (Reptilia, Serpentes, Viperidae)

The complete mitochondrial genome, containing 17,526 bp, was determined from the pitviper Bothrops jararaca. It is the first mitogenome for the most medically important genus of snake in Latin America. This mitogenome has common snake mitochondrial features such as a duplicated control region that has nearly identical sequences at two different locations of the mitogenome and a translocation of tRNA-Leu (UUR). Besides, we found a translocation of the tRNA-Pro compared to Colubridae snakes. Finally, an unusual possible duplication containing a tRNA-Phe was observed for the first time and may represent a marker of the genus.

A homozygous loss-of-function mutation in inositol monophosphatase 1 (IMPA1) causes severe intellectual disability

The genetic basis of intellectual disability (ID) is extremely heterogeneous and relatively little is known about the role of autosomal recessive traits. In a field study performed in a highly inbred area of Northeastern Brazil, we identified and investigated a large consanguineous family with nine adult members affected by severe ID associated with disruptive behavior. The Genome-Wide Human SNP Array 6.0 microarray was used to determine regions of homozygosity by descent from three affected and one normal family member. Whole-exome sequencing (WES) was performed in one affected patient using the Nextera Rapid-Capture Exome kit and Illumina HiSeq2500 system to identify the causative mutation. Potentially deleterious variants detected in regions of homozygosity by descent and not present in either 59 723 unrelated individuals from the Exome Aggregation Consortium (Browser) or 1484 Brazilians were subject to further scrutiny and segregation analysis by Sanger sequencing. Homozygosity-by-descent analysis disclosed a 20.7-Mb candidate region at 8q12.3-q21.2 (lod score: 3.11). WES identified a homozygous deleterious variant in inositol monophosphatase 1 (IMPA1) (NM_005536), consisting of a 5-bp duplication (c.489_493dupGGGCT; chr8: 82,583,247; GRCh37/hg19) leading to a frameshift and a premature stop codon (p.Ser165Trpfs*10) that cosegregated with the disease in 26 genotyped family members. The IMPA1 gene product is responsible for the final step of biotransformation of inositol triphosphate and diacylglycerol, two second messengers. Despite its many physiological functions, no clinical phenotype has been assigned to this gene dysfunction to date. Additionally, IMPA1 is the main target of lithium, a drug that is at the forefront of treatment for bipolar disorder.

Homozygous missense mutation in MED25 segregates with syndromic intellectual disability in a large consanguineous family

BACKGROUND

Intellectual disability (ID) is a highly heterogeneous condition affecting 2% of the population worldwide. In a field study conducted in a highly inbred area of Northeastern Brazil, we investigated a consanguineous family in which seven adults presented syndromic ID.

METHODS

Genome-Wide Human SNP Array 6.0 (Affymetrix) microarray was used to determine regions of homozygosity-by-descent and whole exome sequencing (WES) was performed in one affected individual using Extended Nextera Rapid-Capture Exome and Illumina HiSeq2500.

RESULTS

We found two regions with an logarithm of the odds (LOD) score of 3.234: a region spanning 4.0 Mb in 19q13.32-q13.33 and a pericentromeric 20 Mb area in chromosome 2 (2p12-q11.2). WES disclosed in the critical region of chromosome 19 a homozygous variant (c.418C>T, p.Arg140Trp) in Mediator complex subunit 25 (MED25), predicted as deleterious by PolyPhen-2, Provean, Mutation Taster and Sorting Intolerant From Tolerant (SIFT). MED25 is a component of the Mediator complex, involved in regulation of transcription of nearly all RNA polymerase II-dependent genes. Deleterious mutations in MED12, MED17 and MED23 have already been associated with ID.

CONCLUSIONS

These findings demonstrate that the combination of field investigation of families in highly inbred regions with modern technologies is an effective way for identifying new genes associated with ID.

Exome sequencing for mucolipidosis III: Detection of a novel GNPTAB gene mutation in a patient with a very mild phenotype

Mucolipidosis II and III alpha/beta (ML II/III alpha/beta) are rare autosomal recessive lysosomal storage diseases that are caused by a deficiency of UDP-GlcNAc:lysosomal enzyme N-acetylglucosamine-1-phosphotransferase, the enzyme responsible for the synthesis of the mannose 6-phosphate targeting signal on lysosomal hydrolases. A Brazilian patient suspected of having a very mild ML III was investigated using whole next-generation sequencing (NGS). Two mutations in the GNPTAB gene were detected and confirmed to be in trans status by parental analysis: c.1208T>C (p.Ile403Thr), previously reported as being pathogenic, and the novel mutation c.1723G>A (p.Gly575Arg). This study demonstrates the effectiveness of using whole NGS for the molecular diagnosis of very mild ML III alpha/beta patients.

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.

A defect in the RNA-processing protein HNRPDL causes limb-girdle muscular dystrophy 1G (LGMD1G)

Limb-girdle muscular dystrophies (LGMD) are a heterogeneous group of genetically determined muscle disorders with a primary or predominant involvement of the pelvic or shoulder girdle musculature. More than 20 genes with autosomal recessive (LGMD2A to LGMD2Q) and autosomal dominant inheritance (LGMD1A to LGMD1H) have been mapped/identified to date. Mutations are known for six among the eight mapped autosomal dominant forms: LGMD1A (myotilin), LGMD1B (lamin A/C), LGMD1C (caveolin-3), LGMD1D (desmin), LGMD1E (DNAJB6), and more recently for LGMD1F (transportin-3). Our group previously mapped the LGMD1G gene at 4q21 in a Caucasian-Brazilian family. We now mapped a Uruguayan family with patients displaying a similar LGMD1G phenotype at the same locus. Whole genome sequencing identified, in both families, mutations in the HNRPDL gene. HNRPDL is a heterogeneous ribonucleoprotein family member, which participates in mRNA biogenesis and metabolism. Functional studies performed in S. cerevisiae showed that the loss of HRP1 (yeast orthologue) had pronounced effects on both protein levels and cell localizations, and yeast proteome revealed dramatic reorganization of proteins involved in RNA-processing pathways. In vivo analysis showed that hnrpdl is important for muscle development in zebrafish, causing a myopathic phenotype when knocked down. The present study presents a novel association between a muscular disorder and a RNA-related gene and reinforces the importance of RNA binding/processing proteins in muscle development and muscle disease. Understanding the role of these proteins in muscle might open new therapeutic approaches for muscular dystrophies.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15×-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome.

The Complete Nucleotide Sequence and Genomic Organization of Citrus Leprosis Associated Virus, Cytoplasmatic type (CiLV-C)

The Citrus leprosis disease (CiL) is associated to a virus (CiLV) transmitted by Brevipalpus spp. mites (Acari: Tenuipalpidae). CiL is endemic in Brazil and its recently spreading to Central America represents a threat to citrus industry in the USA. Electron microscopy images show two forms of CiLV: a rare nuclear form, characterized by rod-shaped naked particle (CiLV-N) and a common cytoplasmic form (CiLV-C) associated with bacilliform-enveloped particle and cytoplasmic viroplasm. Due to this morphological feature, CiLV-C has been treated as Rhabdovirus-like. In this paper we present the complete nucleotide sequence and genomic organization of CiLV-C. It is a bipartite virus with sequence similarity to ssRNA positive plant virus. RNA1 encodes a putative replicase polyprotein and an ORF with no known function. RNA2 encodes 4 ORFs. pl5, p24 and p61 have no significant similarity to any known proteins and p32 encodes a protein with similarity to a viral movement protein. The CiLV-C sequences are associated with typical symptoms of CiL by RT-PCR. Phylogenetic analysis suggests that CiLV-C is probably a member of a new family of plant virus evolutionarily related to Tobamovirus.

Brazilian coffee genome project: an EST-based genomic resource

Coffee is one of the most valuable agricultural commodities and ranks second on international trade exchanges. The genus Coffea belongs to the Rubiaceae family which includes other important plants. The genus contains about 100 species but commercial production is based only on two species, Coffea arabica and Coffea canephora that represent about 70 % and 30 % of the total coffee market, respectively. The Brazilian Coffee Genome Project was designed with the objective of making modern genomics resources available to the coffee scientific community, working on different aspects of the coffee production chain. We have single-pass sequenced a total of 214,964 randomly picked clones from 37 cDNA libraries of C. arabica, C. canephora and C. racemosa, representing specific stages of cells and plant development that after trimming resulted in 130,792, 12,381 and 10,566 sequences for each species, respectively. The ESTs clustered into 17,982 clusters and 32,155 singletons. Blast analysis of these sequences revealed that 22 % had no significant matches to sequences in the National Center for Biotechnology Information database (of known or unknown function). The generated coffee EST database resulted in the identification of close to 33,000 different unigenes. Annotated sequencing results have been stored in an online database at <A HREF="http://www.lge.ibi.unicamp.br/cafe">http://www.lge.ibi.unicamp.br/cafe</A>. Resources developed in this project provide genetic and genomic tools that may hold the key to the sustainability, competitiveness and future viability of the coffee industry in local and international markets.

Identification and genomic characterization of a new virus (Tymoviridae family) associated with citrus sudden death disease

Citrus sudden death (CSD) is a new disease that has killed approximately 1 million orange trees in Brazil. Here we report the identification of a new virus associated with the disease. RNAs isolated from CSD-affected and nonaffected trees were used to construct cDNA libraries. A set of viral sequences present exclusively in libraries of CSD-affected trees was used to obtain the complete genome sequence of the new virus. Phylogenetic analysis revealed that this virus is a new member of the genus Marafivirus. Antibodies raised against the putative viral coat proteins allowed detection of viral antigens of expected sizes in affected plants. Electron microscopy of purified virus confirmed the presence of typical isometric Marafivirus particles. The screening of 773 affected and nonaffected citrus trees for the presence of the virus showed a 99.7% correlation between disease symptoms and the presence of the virus. We also detected the virus in aphids feeding on affected trees. These results suggest that this virus is likely to be the causative agent of CSD. The virus was named Citrus sudden death-associated virus.

Structural features and transcript-editing analysis of sugarcane (Saccharum officinarum L.) chloroplast genome

The complete nucleotide sequence of the chloroplast genome of sugarcane (Saccharum officinarum) was determined. It consists of 141,182 base-pairs (bp), containing a pair of inverted repeat regions (IR(A), IR(B)) of 22,794 bp each. The IR(A) and IR(B) sequences separate a small single copy region (12,546 bp) and a large single copy (83,048 bp) region. The gene content and relative arrangement of the 116 identified genes (82 peptide-encoding genes, four ribosomal RNA genes, 30 tRNA genes), with the 16 ycf genes, are highly similar to maize. Editing events, defined as C-to-U transitions in the mRNA sequences, were comparable with those observed in maize, rice and wheat. The conservation of gene organization and mRNA editing suggests a common ancestor for the sugarcane and maize plastomes. These data provide the basis for functional analysis of plastid genes and plastid metabolism within the Poaceae. The sugarcane chloroplast DNA sequence is available at GenBank under accession NC005878.

Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis

Leptospira species colonize a significant proportion of rodent populations worldwide and produce life-threatening infections in accidental hosts, including humans. Complete genome sequencing of Leptospira interrogans serovar Copenhageni and comparative analysis with the available Leptospira interrogans serovar Lai genome reveal that despite overall genetic similarity there are significant structural differences, including a large chromosomal inversion and extensive variation in the number and distribution of insertion sequence elements. Genome sequence analysis elucidates many of the novel aspects of leptospiral physiology relating to energy metabolism, oxygen tolerance, two-component signal transduction systems, and mechanisms of pathogenesis. A broad array of transcriptional regulation proteins and two new families of afimbrial adhesins which contribute to host tissue colonization in the early steps of infection were identified. Differences in genes involved in the biosynthesis of lipopolysaccharide O side chains between the Copenhageni and Lai serovars were identified, offering an important starting point for the elucidation of the organism's complex polysaccharide surface antigens. Differences in adhesins and in lipopolysaccharide might be associated with the adaptation of serovars Copenhageni and Lai to different animal hosts. Hundreds of genes encoding surface-exposed lipoproteins and transmembrane outer membrane proteins were identified as candidates for development of vaccines for the prevention of leptospirosis.

ORESTES are enriched in rare exon usage variants affecting the encoded proteins

A significant fraction of the variability found in the human transcriptome is due to alternative splicing, including alternative exon usage (AEU), intron retention and use of cryptic splice sites. We present a comparison of a large-scale analysis of AEU in the human transcriptome through genome mapping of Open Reading Frame ESTs (ORESTES) and conventional ESTs. It is shown here that ORESTES probe low abundant messages more efficiently. In addition, most of the variants detected by ORESTES affect the structure of the corresponding proteins.

Comparative analyses of the complete genome sequences of Pierce's disease and citrus variegated chlorosis strains of Xylella fastidiosa

Xylella fastidiosa is a xylem-dwelling, insect-transmitted, gamma-proteobacterium that causes diseases in many plants, including grapevine, citrus, periwinkle, almond, oleander, and coffee. X. fastidiosa has an unusually broad host range, has an extensive geographical distribution throughout the American continent, and induces diverse disease phenotypes. Previous molecular analyses indicated three distinct groups of X. fastidiosa isolates that were expected to be genetically divergent. Here we report the genome sequence of X. fastidiosa (Temecula strain), isolated from a naturally infected grapevine with Pierce's disease (PD) in a wine-grape-growing region of California. Comparative analyses with a previously sequenced X. fastidiosa strain responsible for citrus variegated chlorosis (CVC) revealed that 98% of the PD X. fastidiosa Temecula genes are shared with the CVC X. fastidiosa strain 9a5c genes. Furthermore, the average amino acid identity of the open reading frames in the strains is 95.7%. Genomic differences are limited to phage-associated chromosomal rearrangements and deletions that also account for the strain-specific genes present in each genome. Genomic islands, one in each genome, were identified, and their presence in other X. fastidiosa strains was analyzed. We conclude that these two organisms have identical metabolic functions and are likely to use a common set of genes in plant colonization and pathogenesis, permitting convergence of functional genomic strategies.

Comparison of the genomes of two Xanthomonas pathogens with differing host specificities

The genus Xanthomonas is a diverse and economically important group of bacterial phytopathogens, belonging to the gamma-subdivision of the Proteobacteria. Xanthomonas axonopodis pv. citri (Xac) causes citrus canker, which affects most commercial citrus cultivars, resulting in significant losses worldwide. Symptoms include canker lesions, leading to abscission of fruit and leaves and general tree decline. Xanthomonas campestris pv. campestris (Xcc) causes black rot, which affects crucifers such as Brassica and Arabidopsis. Symptoms include marginal leaf chlorosis and darkening of vascular tissue, accompanied by extensive wilting and necrosis. Xanthomonas campestris pv. campestris is grown commercially to produce the exopolysaccharide xanthan gum, which is used as a viscosifying and stabilizing agent in many industries. Here we report and compare the complete genome sequences of Xac and Xcc. Their distinct disease phenotypes and host ranges belie a high degree of similarity at the genomic level. More than 80% of genes are shared, and gene order is conserved along most of their respective chromosomes. We identified several groups of strain-specific genes, and on the basis of these groups we propose mechanisms that may explain the differing host specificities and pathogenic processes.

Comparative genomic analysis of plant-associated bacteria

This review deals with a comparative analysis of seven genome sequences from plant-associated bacteria. These are the genomes of Agrobacterium tumefaciens, Mesorhizobium loti, Sinorhizobium meliloti, Xanthomonas campestris pv campestris, Xanthomonas axonopodis pv citri, Xylella fastidiosa, and Ralstonia solanacearum. Genome structure and the metabolism pathways available highlight the compromise between the genome size and lifestyle. Despite the recognized importance of the type III secretion system in controlling host compatibility, its presence is not universal in all necrogenic pathogens. Hemolysins, hemagglutinins, and some adhesins, previously reported only for mammalian pathogens, are present in most organisms discussed. Different numbers and combinations of cell wall degrading enzymes and genes to overcome the oxidative burst generally induced by the plant host are characterized in these genomes. A total of 19 genes not involved in housekeeping functions were found common to all these bacteria.

The genome of the natural genetic engineer Agrobacterium tumefaciens C58

The 5.67-megabase genome of the plant pathogen Agrobacterium tumefaciens C58 consists of a circular chromosome, a linear chromosome, and two plasmids. Extensive orthology and nucleotide colinearity between the genomes of A. tumefaciens and the plant symbiont Sinorhizobium meliloti suggest a recent evolutionary divergence. Their similarities include metabolic, transport, and regulatory systems that promote survival in the highly competitive rhizosphere; differences are apparent in their genome structure and virulence gene complement. Availability of the A. tumefaciens sequence will facilitate investigations into the molecular basis of pathogenesis and the evolutionary divergence of pathogenic and symbiotic lifestyles.

The genome sequence of the plant pathogen Xylella fastidiosa. The Xylella fastidiosa Consortium of the Organization for Nucleotide Sequencing and Analysis

Xylella fastidiosa is a fastidious, xylem-limited bacterium that causes a range of economically important plant diseases. Here we report the complete genome sequence of X. fastidiosa clone 9a5c, which causes citrus variegated chlorosis--a serious disease of orange trees. The genome comprises a 52.7% GC-rich 2,679,305-base-pair (bp) circular chromosome and two plasmids of 51,158 bp and 1,285 bp. We can assign putative functions to 47% of the 2,904 predicted coding regions. Efficient metabolic functions are predicted, with sugars as the principal energy and carbon source, supporting existence in the nutrient-poor xylem sap. The mechanisms associated with pathogenicity and virulence involve toxins, antibiotics and ion sequestration systems, as well as bacterium-bacterium and bacterium-host interactions mediated by a range of proteins. Orthologues of some of these proteins have only been identified in animal and human pathogens; their presence in X. fastidiosa indicates that the molecular basis for bacterial pathogenicity is both conserved and independent of host. At least 83 genes are bacteriophage-derived and include virulence-associated genes from other bacteria, providing direct evidence of phage-mediated horizontal gene transfer.