ABSTRACT
Narrow-leafed lupin (NLL; Lupinus angustifolius) is a key rotational crop for sustainable farming systems, whose grain is high in protein content. It is a gluten-free, non-genetically modified, alternative protein source to soybean (Glycine max) and as such has gained interest as a human food ingredient. Here, we present a chromosome-length reference genome for the species and a pan-genome assembly comprising 55 NLL lines, including Australian and European cultivars, breeding lines and wild accessions. We present the core and variable genes for the species and report on the absence of essential mycorrhizal associated genes. The genome and pan-genomes of NLL and its close relative white lupin (Lupinus albus) are compared. Furthermore, we provide additional evidence supporting LaRAP2-7 as the key alkaloid regulatory gene for NLL and demonstrate the NLL genome is underrepresented in classical NLR disease resistance genes compared to other sequenced legume species. The NLL genomic resources generated here coupled with previously generated RNA sequencing datasets provide new opportunities to fast-track lupin crop improvement.
Subject(s)
Lupinus , Australia , Chromosomes , Genomics , Humans , Lupinus/genetics , Plant BreedingABSTRACT
KEY MESSAGE: A plant-specific Trimethylguanosine Synthase1-like homologue was identified as a candidate gene for the efl mutation in narrow-leafed lupin, which alters phenology by reducing vernalisation requirement. The vernalisation pathway is a key component of flowering time control in plants from temperate regions but is not well understood in the legume family. Here we examined vernalisation control in the temperate grain legume species, narrow-leafed lupin (Lupinus angustifolius L.), and discovered a candidate gene for an ethylene imine mutation (efl). The efl mutation changes phenology from late to mid-season flowering and additionally causes transformation from obligate to facultative vernalisation requirement. The efl locus was mapped to pseudochromosome NLL-10 in a recombinant inbred line (RIL) mapping population developed by accelerated single seed descent. Candidate genes were identified in the reference genome, and a diverse panel of narrow-leafed lupins was screened to validate mutations specific to accessions with efl. A non-synonymous SNP mutation within an S-adenosyl-L-methionine-dependent methyltransferase protein domain of a Trimethylguanosine Synthase1-like (TGS1) orthologue was identified as the candidate mutation giving rise to efl. This mutation caused substitution of an amino acid within an established motif at a position that is otherwise highly conserved in several plant families and was perfectly correlated with the efl phenotype in F2 and F6 genetic population and a panel of diverse accessions, including the original efl mutant. Expression of the TGS1 homologue did not differ between wild-type and efl genotypes, supporting altered functional activity of the gene product. This is the first time a TGS1 orthologue has been associated with vernalisation response and flowering time control in any plant species.
Subject(s)
Flowers/growth & development , Gene Expression Regulation, Plant , Genetics, Population , Lupinus/growth & development , Methyltransferases/metabolism , Plant Leaves/growth & development , Plant Proteins/metabolism , Flowers/genetics , Lupinus/genetics , Methyltransferases/genetics , Mutation , Phenotype , Phylogeny , Plant Leaves/genetics , Plant Proteins/geneticsABSTRACT
BACKGROUND: Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown. RESULTS: By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse. CONCLUSION: This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.
Subject(s)
Gene Expression Regulation, Developmental , Gene Expression Regulation, Plant , Lupinus/genetics , MicroRNAs/genetics , Seeds/growth & development , Seeds/genetics , Argonaute Proteins/genetics , Biological Variation, Population , Computational Biology , Databases, Genetic , Gene Duplication , Gene Regulatory Networks , Gene Silencing , Genome, Plant , Germination/genetics , High-Throughput Nucleotide Sequencing , Lupinus/growth & development , Phylogeny , Plant Leaves/genetics , Plant Roots/genetics , Transcription Factors/genetics , TranscriptomeABSTRACT
Rhizoctonia solani AG1-IA is a soil-borne necrotrophic pathogen that causes devastating rice sheath blight disease in rice-growing regions worldwide. Sclerotia play an important role in the life cycle of R. solani AG1-IA. In this study, RNA sequencing was used to investigate the transcriptomic dynamics of sclerotial development (SD) of R. solani AG1-IA. Gene ontology and pathway enrichment analyses using the Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to investigate the functions and pathways of differentially expressed genes (DEGs). Six cDNA libraries were generated, and more than 300 million clean reads were obtained and assembled into 15,100 unigenes. In total, 12,575 differentially expressed genes were identified and 34.62% (4353) were significantly differentially expressed with a FDR ≤ 0.01 and |log2Ratio| ≥ 1, which were enriched into eight profiles using Short Time-series Expression Miner. Furthermore, KEGG and gene ontology analyses suggest the DEGs were significantly enriched in several biological processes and pathways, including binding and catalytic functions, biosynthesis of ribosomes, and other biological functions. Further annotation of the DEGs using the Clusters of Orthologous Groups (COG) database found most DEGs were involved in amino acid transport and metabolism, as well as energy production and conversion. Furthermore, DEGs relevant to SD of R. solani AG1-IA were involved in secondary metabolite biosynthesis, melanin biosynthesis, ubiquitin processes, autophagy, and reactive oxygen species metabolism. The gene expression profiles of 10 randomly selected DEGs were validated by quantitative real-time reverse transcription PCR and were consistent with the dynamics in transcript abundance identified by RNA sequencing. The data provide a high-resolution map of gene expression during SD, a key process contributing to the pathogenicity of this devastating pathogen. In addition, this study provides a useful resource for further studies on the genomics of R. solani AG1-IA and other Rhizoctonia species.
Subject(s)
Fungal Proteins/genetics , Gene Expression Regulation, Fungal , Morphogenesis/genetics , Oryza/microbiology , Plant Diseases/genetics , Rhizoctonia/growth & development , Transcriptome , Fungal Proteins/metabolism , Gene Expression Profiling , Genome, Fungal , Plant Diseases/microbiology , Rhizoctonia/genetics , Rhizoctonia/pathogenicityABSTRACT
Narrow-leafed lupin (Lupinus angustifolius L.) cultivation was transformed by 2 dominant vernalization-insensitive, early flowering time loci known as Ku and Julius (Jul), which allowed expansion into shorter season environments. However, reliance on these loci has limited genetic and phenotypic diversity for environmental adaptation in cultivated lupin. We recently predicted that a 1,423-bp deletion in the cis-regulatory region of LanFTc1, a FLOWERING LOCUS T (FT) homologue, derepressed expression of LanFTc1 and was the underlying cause of the Ku phenotype. Here, we surveyed diverse germplasm for LanFTc1 cis-regulatory variation and identified 2 further deletions of 1,208 and 5,162Ā bp in the 5' regulatory region, which overlap the 1,423-bp deletion. Additionally, we confirmed that no other polymorphisms were perfectly associated with vernalization responsiveness. Phenotyping and gene expression analyses revealed that Jul accessions possessed the 5,162-bp deletion and that the Jul and Ku deletions were equally capable of removing vernalization requirement and up-regulating gene expression. The 1,208-bp deletion was associated with intermediate phenology, vernalization responsiveness, and gene expression and therefore may be useful for expanding agronomic adaptation of lupin. This insertion/deletion series may also help resolve how the vernalization response is mediated at the molecular level in legumes.
Subject(s)
Flowers/growth & development , Genes, Plant/genetics , INDEL Mutation/genetics , Lupinus/genetics , Flowers/genetics , Gene Expression Regulation, Plant/genetics , Gene Expression Regulation, Plant/physiology , Genes, Plant/physiology , Genetic Variation/genetics , INDEL Mutation/physiology , Linkage Disequilibrium/genetics , Lupinus/growth & development , Polymorphism, Single Nucleotide/genetics , Promoter Regions, Genetic/genetics , SeasonsABSTRACT
Quinolizidine alkaloids (QAs) are toxic secondary metabolites that complicate the end use of narrow-leafed lupin (NLL; Lupinus angustifolius L.) grain, as levels sometimes exceed the industry limit for its use as a food and feed source. The genotypic and environmental influences on QA production in NLL are poorly understood. Here, the expression of QA biosynthetic genes was analysed in vegetative and reproductive tissues of bitter (high QA) and sweet (low QA) accessions. It was demonstrated that sweet accessions are characterized by lower QA biosynthetic gene expression exclusively in leaf and stem tissues than bitter NLL, consistent with the hypothesis that QAs are predominantly produced in aerial tissues and transported to seeds, rather than synthesized within the seed itself. This analysis informed our identification of additional candidate genes involved in QA biosynthesis. Drought and temperature stress are two major abiotic stresses that often occur during NLL pod set. Hence, we assessed the effect of drought, increased temperature, and their combination, on QA production in three sweet NLL cultivars. A cultivar-specific response to drought and temperature in grain QA levels was observed, including the identification of a cultivar where alkaloid levels did not change with these stress treatments.
Subject(s)
Alkaloids/biosynthesis , Lupinus/genetics , Plant Proteins/genetics , Quinolizidines/metabolism , Alkaloids/metabolism , Droughts , Gene Expression Regulation, Plant , Lupinus/physiology , Plant Leaves/genetics , Plant Leaves/metabolism , Seeds/genetics , Seeds/growth & development , Seeds/metabolism , Stress, Physiological , TemperatureABSTRACT
The root-infecting necrotrophic fungal pathogen Rhizoctoniasolani causes significant disease to all the world's major food crops. As a model for pathogenesis of legumes, we have examined the interaction of R. solani AG8 with Medicago truncatula. RNAseq analysis of the moderately resistant M. truncatula accession A17 and highly susceptible sickle (skl) mutant (defective in ethylene sensing) identified major early transcriptional reprogramming in A17. Responses specific to A17 included components of ethylene signaling, reactive oxygen species metabolism, and consistent upregulation of the isoflavonoid biosynthesis pathway. Mass spectrometry revealed accumulation of the isoflavonoid-related compounds liquiritigenin, formononetin, medicarpin, and biochanin A in A17. Overexpression of an isoflavone synthase in M. truncatula roots increased isoflavonoid accumulation and resistance to R. solani. Addition of exogenous medicarpin suggested this phytoalexin may be one of several isoflavonoids required to contribute to resistance to R. solani. Together, these results provide evidence for the role of ethylene-mediated accumulation of isoflavonoids during defense against root pathogens in legumes. The involvement of ethylene signaling and isoflavonoids in the regulation of both symbiont-legume and pathogen-legume interactions in the same tissue may suggest tight regulation of these responses are required in the root tissue.
Subject(s)
Disease Resistance , Ethylenes/metabolism , Isoflavones/metabolism , Medicago truncatula/microbiology , Plant Diseases/microbiology , Plant Roots/microbiology , Rhizoctonia/physiology , Signal Transduction , Biosynthetic Pathways/genetics , Disease Resistance/genetics , Gene Expression Profiling , Gene Expression Regulation, Plant , Medicago truncatula/genetics , Medicago truncatula/immunology , Medicago truncatula/metabolism , Metabolome/genetics , Mutation/genetics , Phenotype , Plant Roots/genetics , Plant Roots/metabolism , Rhizoctonia/growth & development , Transcription, GeneticABSTRACT
Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609Ā Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.
Subject(s)
Genome, Plant/genetics , Lupinus/genetics , Lupinus/microbiology , Plant Proteins/genetics , Disease Resistance/genetics , Disease Resistance/physiology , Plant Proteins/physiology , Polyploidy , Synteny/geneticsABSTRACT
BACKGROUND: Pathogenic members of the Fusarium oxysporum species complex are responsible for vascular wilt disease on many important crops including legumes, where they can be one of the most destructive disease causing necrotrophic fungi. We previously developed a model legume-infecting pathosystem based on the reference legume Medicago truncatula and a pathogenic F. oxysporum forma specialis (f. sp.) medicaginis (Fom). To dissect the molecular pathogenicity arsenal used by this root-infecting pathogen, we sequenced its transcriptome during infection of a susceptible and resistant host accession. RESULTS: High coverage RNA-Seq of Fom infected root samples harvested from susceptible (DZA315) or resistant (A17) M. truncatula seedlings at early or later stages of infection (2 or 7Ā days post infection (dpi)) and from vegetative (in vitro) samples facilitated the identification of unique and overlapping sets of in planta differentially expressed genes. This included enrichment, particularly in DZA315 in planta up-regulated datasets, for proteins associated with sugar, protein and plant cell wall metabolism, membrane transport, nutrient uptake and oxidative processes. Genes encoding effector-like proteins were identified, including homologues of the F. oxysporum f. sp. lycopersici Secreted In Xylem (SIX) proteins, and several novel candidate effectors based on predicted secretion, small protein size and high in-planta induced expression. The majority of the effector candidates contain no known protein domains but do share high similarity to predicted proteins predominantly from other F. oxysporum ff. spp. as well as other Fusaria (F. solani, F. fujikori, F. verticilloides, F. graminearum and F. pseudograminearum), and from another wilt pathogen of the same class, a Verticillium species. Overall, this suggests these novel effector candidates may play important roles in Fusaria and wilt pathogen virulence. CONCLUSION: Combining high coverage in planta RNA-Seq with knowledge of fungal pathogenicity protein features facilitated the identification of differentially expressed pathogenicity associated genes and novel effector candidates expressed during infection of a resistant or susceptible M. truncatula host. The knowledge from this first in depth in planta transcriptome sequencing of any F. oxysporum ff. spp. pathogenic on legumes will facilitate the dissection of Fusarium wilt pathogenicity mechanisms on many important legume crops.
Subject(s)
Fusarium/genetics , Gene Expression Profiling , Host-Pathogen Interactions , Medicago truncatula/microbiology , Plant Diseases/microbiology , Transcriptome , Fusarium/pathogenicity , Gene Expression Regulation, Fungal , High-Throughput Nucleotide Sequencing , Phenotype , Virulence/geneticsABSTRACT
BACKGROUND: Soil-borne fungi of the Fusarium oxysporum species complex cause devastating wilt disease on many crops including legumes that supply human dietary protein needs across many parts of the globe. We present and compare draft genome assemblies for three legume-infecting formae speciales (ff. spp.): F. oxysporum f. sp. ciceris (Foc-38-1) and f. sp. pisi (Fop-37622), significant pathogens of chickpea and pea respectively, the world's second and third most important grain legumes, and lastly f. sp. medicaginis (Fom-5190a) for which we developed a model legume pathosystem utilising Medicago truncatula. RESULTS: Focusing on the identification of pathogenicity gene content, we leveraged the reference genomes of Fusarium pathogens F. oxysporum f. sp. lycopersici (tomato-infecting) and F. solani (pea-infecting) and their well-characterised core and dispensable chromosomes to predict genomic organisation in the newly sequenced legume-infecting isolates. Dispensable chromosomes are not essential for growth and in Fusarium species are known to be enriched in host-specificity and pathogenicity-associated genes. Comparative genomics of the publicly available Fusarium species revealed differential patterns of sequence conservation across F. oxysporum formae speciales, with legume-pathogenic formae speciales not exhibiting greater sequence conservation between them relative to non-legume-infecting formae speciales, possibly indicating the lack of a common ancestral source for legume pathogenicity. Combining predicted dispensable gene content with in planta expression in the model legume-infecting isolate, we identified small conserved regions and candidate effectors, four of which shared greatest similarity to proteins from another legume-infecting ff. spp. CONCLUSIONS: We demonstrate that distinction of core and potential dispensable genomic regions of novel F. oxysporum genomes is an effective tool to facilitate effector discovery and the identification of gene content possibly linked to host specificity. While the legume-infecting isolates didn't share large genomic regions of pathogenicity-related content, smaller regions and candidate effector proteins were highly conserved, suggesting that they may play specific roles in inducing disease on legume hosts.
Subject(s)
Fabaceae/microbiology , Fusarium/genetics , Genome, Fungal , Comparative Genomic Hybridization , Conserved Sequence , DNA, Fungal/genetics , Fungal Proteins/genetics , Fusarium/classification , Host Specificity , Molecular Sequence Annotation , Phylogeny , Plant Diseases/microbiology , Sequence Analysis, DNAABSTRACT
UNLABELLED: Liver-targeted gene therapy based on recombinant adeno-associated viral vectors (rAAV) shows promising therapeutic efficacy in animal models and adult-focused clinical trials. This promise, however, is not directly translatable to the growing liver, where high rates of hepatocellular proliferation are accompanied by loss of episomal rAAV genomes and subsequently a loss in therapeutic efficacy. We have developed a hybrid rAAV/piggyBac transposon vector system combining the highly efficient liver-targeting properties of rAAV with stable piggyBac-mediated transposition of the transgene into the hepatocyte genome. Transposition efficiency was first tested using an enhanced green fluorescent protein expression cassette following delivery to newborn wild-type mice, with a 20-fold increase in stably gene-modified hepatocytes observed 4 weeks posttreatment compared to traditional rAAV gene delivery. We next modeled the therapeutic potential of the system in the context of severe urea cycle defects. A single treatment in the perinatal period was sufficient to confer robust and stable phenotype correction in the ornithine transcarbamylase-deficient Spf(ash) mouse and the neonatal lethal argininosuccinate synthetase knockout mouse. Finally, transposon integration patterns were analyzed, revealing 127,386 unique integration sites which conformed to previously published piggyBac data. CONCLUSION: Using a hybrid rAAV/piggyBac transposon vector system, we achieved stable therapeutic protection in two urea cycle defect mouse models; a clinically conceivable early application of this technology in the management of severe urea cycle defects could be as a bridging therapy while awaiting liver transplantation; further improvement of the system will result from the development of highly human liver-tropic capsids, the use of alternative strategies to achieve transient transposase expression, and engineered refinements in the safety profile of piggyBac transposase-mediated integration.
Subject(s)
Adenoviridae/genetics , Genetic Therapy/methods , Genetic Vectors/pharmacology , Hyperammonemia/therapy , Urea/metabolism , Animals , Animals, Newborn , Disease Models, Animal , Gene Transfer Techniques , Humans , Hyperammonemia/diagnosis , Liver Diseases/therapy , Mice , Mice, Transgenic , Severity of Illness Index , Statistics, NonparametricABSTRACT
The chromosome-centric human proteome project (C-HPP) aims to define the complete set of proteins encoded in each human chromosome. The neXtProt database (September 2013) lists 20,128 proteins for the human proteome, of which 3831 human proteins (Ć¢ĀĀ¼19%) are considered "missing" according to the standard metrics table (released September 27, 2013). In support of the C-HPP initiative, we have extended the annotation strategy developed for human chromosome 7 "missing" proteins into a semiautomated pipeline to functionally annotate the "missing" human proteome. This pipeline integrates a suite of bioinformatics analysis and annotation software tools to identify homologues and map putative functional signatures, gene ontology, and biochemical pathways. From sequential BLAST searches, we have primarily identified homologues from reviewed nonhuman mammalian proteins with protein evidence for 1271 (33.2%) "missing" proteins, followed by 703 (18.4%) homologues from reviewed nonhuman mammalian proteins and subsequently 564 (14.7%) homologues from reviewed human proteins. Functional annotations for 1945 (50.8%) "missing" proteins were also determined. To accelerate the identification of "missing" proteins from proteomics studies, we generated proteotypic peptides in silico. Matching these proteotypic peptides to ENCODE proteogenomic data resulted in proteomic evidence for 107 (2.8%) of the 3831 "missing proteins, while evidence from a recent membrane proteomic study supported the existence for another 15 "missing" proteins. The chromosome-wise functional annotation of all "missing" proteins is freely available to the scientific community through our web server (http://biolinfo.org/protannotator).
Subject(s)
Automation , Chromosomes, Human , Proteome , Databases, Protein , Humans , SoftwareABSTRACT
The chromosome-centric human proteome project aims to systematically map all human proteins, chromosome by chromosome, in a gene-centric manner through dedicated efforts from national and international teams. This mapping will lead to a knowledge-based resource defining the full set of proteins encoded in each chromosome and laying the foundation for the development of a standardized approach to analyze the massive proteomic data sets currently being generated. The neXtProt database lists 946 proteins as the human proteome of chromosome 7. However, 170 (18%) proteins of human chromosome 7 have no evidence at the proteomic, antibody, or structural levels and are considered "missing" in this study as they lack experimental support. We have developed a protocol for the functional annotation of these "missing" proteins by integrating several bioinformatics analysis and annotation tools, sequential BLAST homology searches, protein domain/motif and gene ontology (GO) mapping, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Using the BLAST search strategy, homologues for reviewed non-human mammalian proteins with protein evidence were identified for 90 "missing" proteins while another 38 had reviewed non-human mammalian homologues. Putative functional annotations were assigned to 27 of the remaining 43 novel proteins. Proteotypic peptides have been computationally generated to facilitate rapid identification of these proteins. Four of the "missing" chromosome 7 proteins have been substantiated by the ENCODE proteogenomic peptide data.
Subject(s)
Chromosomes, Human, Pair 7 , Genome, Human , Human Genome Project , Mammals/genetics , Molecular Sequence Annotation , Proteome/genetics , Algorithms , Animals , Chromosome Mapping , Databases, Protein , Humans , Mammals/metabolism , Proteome/metabolism , Sequence Homology, Amino AcidABSTRACT
The black PĆ©rigord truffle (Tuber melanosporum Vittad.) is a highly prized food today, with its unique scent (i.e., perfume) and texture. Despite these attributes, it remains relatively poorly studied, lacking "omics" information to characterize its biology and biochemistry, especially changes associated with freshness and the proteins/metabolites responsible for its organoleptic properties. In this study, we have functionally annotated the truffle proteome from the 2010 T. melanosporum genome comprising 12,771 putative nonredundant proteins. Using sequential BLAST search strategies, we identified homologues for 2587 proteins with 2486 (96.0%) fungal homologues (available from http://biolinfo.org/protannotator/blacktruffle.php). A combined 1D PAGE and high-accuracy LC-MS/MS proteomic study was employed to validate the results of the functional annotation and identified 836 (6.5%) proteins, of which 47.5% (i.e., 397) were present in our bioinformatics studies. Our study, functionally annotating 6487 black PĆ©rigord truffle proteins and confirming 836 by proteomic experiments, is by far the most comprehensive study to date contributing significantly to the scientific community. This study has resulted in the functional characterization of novel proteins to increase our biological understanding of this organism and to uncover potential biomarkers of authenticity, freshness, and perfume maturation.
Subject(s)
Fungal Proteins/genetics , Genome, Fungal , Proteome , Saccharomycetales/genetics , Software , Fungal Proteins/metabolism , Gene Expression , Molecular Sequence Annotation , Odorants/analysis , Saccharomyces cerevisiae/genetics , Saccharomyces cerevisiae/metabolism , Saccharomycetales/metabolism , Sequence Alignment , Sequence Homology, Amino AcidABSTRACT
BACKGROUND: The transcriptome of an organism can be studied with the analysis of expressed sequence tag (EST) data sets that offers a rapid and cost effective approach with several new and updated bioinformatics approaches and tools for assembly and annotation. The comprehensive analyses comprehend an organism along with the genome and proteome analysis. With the advent of large-scale sequencing projects and generation of sequence data at protein and cDNA levels, automated analysis pipeline is necessary to store, organize and annotate ESTs. RESULTS: TranSeqAnnotator is a workflow for large-scale analysis of transcriptomic data with the most appropriate bioinformatics tools for data management and analysis. The pipeline automatically cleans, clusters, assembles and generates consensus sequences, conceptually translates these into possible protein products and assigns putative function based on various DNA and protein similarity searches. Excretory/secretory (ES) proteins inferred from ESTs/short reads are also identified. The TranSeqAnnotator accepts FASTA format raw and quality ESTs along with protein and short read sequences and are analysed with user selected programs. After pre-processing and assembly, the dataset is annotated at the nucleotide, protein and ES protein levels. CONCLUSION: TranSeqAnnotator has been developed in a Linux cluster, to perform an exhaustive and reliable analysis and provide detailed annotation. TranSeqAnnotator outputs gene ontologies, protein functional identifications in terms of mapping to protein domains and metabolic pathways. The pipeline is applied to annotate large EST datasets to identify several novel and known genes with therapeutic experimental validations and could serve as potential targets for parasite intervention. TransSeqAnnotator is freely available for the scientific community at http://estexplorer.biolinfo.org/TranSeqAnnotator/.
Subject(s)
Expressed Sequence Tags , Gene Expression Profiling/statistics & numerical data , Molecular Sequence Annotation/methods , Software , TranscriptomeABSTRACT
BACKGROUND: Helminths are important socio-economic organisms, responsible for causing major parasitic infections in humans, other animals and plants. These infections impose a significant public health and economic burden globally. Exceptionally, some helminth organisms like Caenorhabditis elegans are free-living in nature and serve as model organisms for studying parasitic infections. Excretory/secretory proteins play an important role in parasitic helminth infections which make these proteins attractive targets for therapeutic use. In the case of helminths, large volume of expressed sequence tags (ESTs) has been generated to understand parasitism at molecular level and for predicting excretory/secretory proteins for developing novel strategies to tackle parasitic infections. However, mostly predicted ES proteins are not available for further analysis and there is no repository available for such predicted ES proteins. Furthermore, predictions have, in the main, focussed on classical secretory pathways while it is well established that helminth parasites also utilise non-classical secretory pathways. RESULTS: We developed a free Helminth Secretome Database (HSD), which serves as a repository for ES proteins predicted using classical and non-classical secretory pathways, from EST data for 78 helminth species (64 nematodes, 7 trematodes and 7 cestodes) ranging from parasitic to free-living organisms. Approximately 0.9 million ESTs compiled from the largest EST database, dbEST were cleaned, assembled and analysed by different computational tools in our bioinformatics pipeline and predicted ES proteins were submitted to HSD. CONCLUSION: We report the large-scale prediction and analysis of classically and non-classically secreted ES proteins from diverse helminth organisms. All the Unigenes (contigs and singletons) and excretory/secretory protein datasets generated from this analysis are freely available. A BLAST server is available at http://estexplorer.biolinfo.org/hsd, for checking the sequence similarity of new protein sequences against predicted helminth ES proteins.
Subject(s)
Databases, Factual , Expressed Sequence Tags , Helminth Proteins/metabolism , Helminths/metabolism , Animals , Computational Biology , Helminth Proteins/chemistry , Internet , User-Computer InterfaceABSTRACT
Background: COVID-19 vaccines have been rolled out recently in several parts of the world. Although the protective efficacy is frequently discussed, little is known about the factors associated with COVID-19 vaccine adverse effects. The study was conducted with the aim to evaluate the occurrence of adverse events following immunization (AEFI) with two doses of covishield and covaxin and to assess factors associated with these adverse effects. Methods: A longitudinal study was conducted for a period of three months in the adults above 18 years of age attending rural health training center (RHTC) either to receive their first or second dose of covishield or covaxin. After vaccination, the participants were observed at the health facility for 30 min for any AEFI and also followed up telephonically on seventh day from vaccination. Data was collected on predesigned and pretested questionnaire and appropriate statistical tests were applied. Results: Out of 532 participants, 250 (47%) came for their first dose while 282 (53%) came for second dose. In both the groups maximum participation was seen by males and those belonging to age group 18-30 years. Majority of the participants reported local tenderness (39.3%) after first dose of covaxin and fever (30.5%) after first dose of covishield. Mainly significant association was observed after vaccination in participants with comorbidities. Conclusion: The short-term adverse events with both the vaccines were observed, but these were mild and short lived. In this context, our study becomes more relevant in disseminating short-term safety data post vaccination. This will help individuals in their decision to accept vaccination.
ABSTRACT
BACKGROUND: Excretory/secretory proteins (ESPs) play a major role in parasitic infection as they are present at the host-parasite interface and regulate host immune system. In case of parasitic helminths, transcriptomics has been used extensively to understand the molecular basis of parasitism and for developing novel therapeutic strategies against parasitic infections. However, none of transcriptomic studies have extensively covered ES protein prediction for identifying novel therapeutic targets, especially as parasites adopt non-classical secretion pathways. RESULTS: We developed a semi-automated computational approach for prediction and annotation of ES proteins using transcriptomic data from next generation sequencing platforms. For the prediction of non-classically secreted proteins, we have used an improved computational strategy, together with homology matching to a dataset of experimentally determined parasitic helminth ES proteins. We applied this protocol to analyse 454 short reads of parasitic nematode, Strongyloides ratti. From 296231 reads, we derived 28901 contigs, which were translated into 20877 proteins. Based on our improved ES protein prediction pipeline, we identified 2572 ES proteins, of which 407 (1.9%) proteins have classical N-terminal signal peptides, 923 (4.4%) were computationally identified as non-classically secreted while 1516 (7.26%) were identified by homology to experimentally identified parasitic helminth ES proteins. Out of 2572 ES proteins, 2310 (89.8%) ES proteins had homologues in the free-living nematode Caenorhabditis elegans and 2220 (86.3%) in parasitic nematodes. We could functionally annotate 1591 (61.8%) ES proteins with protein families and domains and establish pathway associations for 691 (26.8%) proteins. In addition, we have identified 19 representative ES proteins, which have no homologues in the host organism but homologous to lethal RNAi phenotypes in C. elegans, as potential therapeutic targets. CONCLUSION: We report a comprehensive approach using freely available computational tools for the secretome analysis of NGS data. This approach has been applied to S. ratti 454 transcriptomic data for in silico excretory/secretory proteins prediction and analysis, providing a foundation for developing new therapeutic solutions for parasitic infections.
Subject(s)
Sequence Analysis, DNA , Transcriptome , Animals , Caenorhabditis elegans/metabolism , Contig Mapping , Databases, Factual , Expressed Sequence Tags , Helminth Proteins/genetics , Helminth Proteins/metabolism , Rats , Search Engine , Strongyloides ratti/metabolismABSTRACT
Drought stress requires plants to adjust their water balance to maintain tissue water levels. Isohydric plants ('water-savers') typically achieve this through stomatal closure, while anisohydric plants ('water-wasters') use osmotic adjustment and maintain stomatal conductance. Isohydry or anisohydry allows plant species to adapt to different environments. In this paper we show that both mechanisms occur in bread wheat (Triticum aestivum L.). Wheat lines with reproductive drought-tolerance delay stomatal closure and are temporarily anisohydric, before closing stomata and become isohydric at higher threshold levels of drought stress. Drought-sensitive wheat is isohydric from the start of the drought treatment. The capacity of the drought-tolerant line to maintain stomatal conductance correlates with repression of ABA synthesis in spikes and flag leaves. Gene expression profiling revealed major differences in the drought response in spikes and flag leaves of both wheat lines. While the isohydric drought-sensitive line enters a passive growth mode (arrest of photosynthesis, protein translation), the tolerant line mounts a stronger stress defence response (ROS protection, LEA proteins, cuticle synthesis). The drought response of the tolerant line is characterised by a strong response in the spike, displaying enrichment of genes involved in auxin, cytokinin and ethylene metabolism/signalling. While isohydry may offer advantages for longer term drought stress, anisohydry may be more beneficial when drought stress occurs during the critical stages of wheat spike development, ultimately improving grain yield.
Subject(s)
Plant Growth Regulators/physiology , Plant Proteins/genetics , Plant Proteins/metabolism , Plant Stomata/physiology , Stress, Physiological , Transcriptome , Triticum/physiology , Droughts , Gene Expression Profiling , Gene Expression Regulation, Plant , Osmotic Pressure , Plant Physiological PhenomenaABSTRACT
Aphids are virus-spreading insect pests affecting crops worldwide and their fast population build-up and insecticide resistance make them problematic to control. Here, we aim to understand the molecular basis of spotted alfalfa aphid (SAA) or Therioaphis trifolii f. maculata resistance in Medicago truncatula, a model organism for legume species. We compared susceptible and resistant near isogenic Medicago lines upon SAA feeding via transcriptome sequencing. Expression of genes involved in defense and stress responses, protein kinase activity and DNA binding were enriched in the resistant line. Potentially underlying some of these changes in gene expression was the finding that members of the MYB, NAC, AP2 domain and ERF transcription factor gene families were differentially expressed in the resistant versus susceptible lines. A TILLING population created in the resistant cultivar was screened using exome capture sequencing and served as a reverse genetics tool to functionally characterise genes involved in the aphid resistance response. This screening revealed three transcription factors (a NAC, AP2 domain and ERF) as important regulators in the defence response, as a premature stop-codon in the resistant background led to a delay in aphid mortality and enhanced plant susceptibility. This combined functional genomics approach will facilitate the future development of pest resistant crops by uncovering candidate target genes that can convey enhanced aphid resistance.